Artificial intelligence algorithms are extremely powerful at fitting massive and complex datasets. But their interpretation, rationalizing how the machine performs specific predictions when presented a given input, is notoriously hard. This black box behavior hampers wide acceptation of AI in medical diagnostics, where justifications matter, and restrain their utility in natural sciences where understanding mechanisms is the goal.
Now, an interdisciplinary team of biologists and computational researchers from the Technical University of Munich, the Stowers Institute for Medical Research and the Stanford University has shown that applying neural networks, such as those used for facial recognition, together with newly developed model interpretation techniques can be used to decipher complex instructions encoded in DNA.
One of the big unsolved problems in biology is the genome’s second code, its regulatory code. The DNA bases encode not only the instructions for how to build proteins, but also when and where to make these proteins in an organism.
The regulatory code is read by proteins called transcription factors that bind to short stretches of DNA called motifs. However, how particular combinations and arrangements of motifs specify regulatory activity is an extremely complex problem that has been hard to pin down.
DNA binding experiments and computational modeling going hand in hand
The key was to perform transcription factor-DNA binding experiments and computational modeling at the highest possible resolution, down to the level of individual DNA bases. The increased resolution allowed the team not only to train highly accurate neural network models, but also to extract the key elements and patterns from the models, including transcription factor binding motifs and the combinatorial rules by which they function together as code.
“Neural networks are black boxes, but they can be interrogated digitally. So, with a large number of virtual experiments we figured out the rules the neural net learned” says first author Dr. Žiga Avsec, member of the group of Julien Gagneur, professor of computational molecular medicine at the Technical University in Munich. Together with Anshul Kundaje, professor at the Stanford University, he created the first version of the model when he visited Stanford as a guest scientist.
Applied to master regulators of stem cell differentiation and confirmed experimentally by CRISPR genomic edition, the approach revealed complex rules involving a precise positioning along the DNA double helix and specific ordering of events.
“This was extremely satisfying,” says project leader Julia Zeitlinger, investigator at the Stowers Institute and professor at the University of Kansas Medical Center, “as the results fit beautifully with existing experimental results, and also revealed novel insights that surprised us.”
A pattern becomes visible: how Nanog binds to DNA
For example, the researchers found that a well-studied transcription factor called Nanog binds cooperatively to DNA when multiples of its motif are present in a periodic fashion such that they appear on the same side of the spiraling DNA helix.
“There has been a long trail of experimental evidence that such motif periodicity sometimes exists in the regulatory code,” Zeitlinger says. However, the exact circumstances were elusive, and Nanog had not been a suspect. Discovering that Nanog has such a pattern, and seeing additional details of its interactions, was surprising because we did not specifically search for this pattern.”
“This is the key advantage of using neural networks for this task. A classic computational model is built on hand-crafted, rigid rules to ensure that it can be interpreted”, says Avsec. “However, biology is extremely rich and complicated. By abandoning the need to interpret individual parameters, we can train much more flexible and nuanced models that capture any biological phenomena, including those yet unknown.”
A powerful bottom-up approach
This neural net model – named BPNet for Base Pair Network – is a powerful bottom-up approach similar to facial recognition in images, where a neural network first detects edges in the pixels, then learns how edges form facial elements like the eye, nose or mouth, and finally how facial elements together form a face.
Instead of learning from pixels, BPNet learns from the raw DNA sequence and learns to detect sequence motifs and eventually the higher-order rules by which the elements predict the base-resolution binding data.
Both the Zeitlinger Lab and the Kundaje Lab are already using BPNet to reliably identify binding motifs for other cell types, relate motifs to biophysical parameters and learn other structural features in the genome such as those associated with DNA packaging. To enable other scientists to use BPNet and adapt it for their own needs, the researchers have made the entire software framework available with documentation and tutorials.
“This work is a technological tour-de-force,” says Julien Gagneur. “It combines deep learning modeling of genome-wide assays down to single-nucleotide resolutions, together with advanced explainable AI techniques allowing to interpret what “the black box” has learned. The methodology will help biologist studying the full regulatory grammar.”
Žiga Avsec, Melanie Weilert, Avanti Shrikumar, Sabrina Krueger, Amr Alexandari, Khyati Dalal, Robin Fropf, Charles McAnany, Julien Gagneur, Anshul Kundaje, and Julia Zeitlinger
nature genetics, Feb. 18, 2021 – DOI: 10.1038/s41588-021-00782-6