Computational Biology seminars

Michelle M. Li examines how computational analyses can preserve the biological context of a proposed treatment, including cell type, disease state, genetic background and patient characteristics. She introduces Medea, an AI system that combines biological software, predictive models and literature retrieval while checking intermediate steps and reconciling evidence. The seminar presents evaluations involving cell-specific target selection, cancer-cell synthetic lethality and immunotherapy response. A separate yeast experiment tests predictions against previously unpublished measurements of gene-pair interactions under DNA-damaging treatments. The research addresses whether an agent can transfer useful evidence between contexts while recognizing when that transfer is unsupported. Reported comparisons cover predictive performance, computational failures and the ability to abstain. This recording retains the original seminar date.

Jiaqi Zhang examines how computational models can learn the effects of genetic interventions from single-cell experiments. Such experiments reveal causal relationships, but their high-dimensional measurements are costly to collect and difficult to interpret. The seminar connects identifiable causal representations with a predictive method for previously unseen perturbations. The approach incorporates prior biological knowledge and changes in data distributions to estimate responses at individual-cell resolution. It also uses predictions to guide subsequent experiments. An application identifies and experimentally validates previously unknown T-cell regulators with potential relevance to cancer immunotherapy. The recording follows the original July seminar; the series lists Zhang at MIT, while the recording biography describes her incoming Columbia appointment.

David Baker, Demis Hassabis and John Jumper examine two complementary computational problems: designing proteins with desired properties and predicting the structures adopted by natural amino acid sequences. Baker describes the development of methods for creating new protein structures and functions. Hassabis discusses the use of artificial intelligence to accelerate scientific discovery, while Jumper explains the integration of chemical and biological reasoning into protein-structure prediction. The programme connects sequence, three-dimensional structure and function, and shows how computational approaches changed what researchers can attempt in protein science. It preserves the distinction between de novo design and prediction while explaining why advances in both can support biological understanding and the creation of useful molecules.

Fate decisions in developing tissues involve cells transitioning between a set of discrete cell states. Geometric models, often referred to as Waddington landscapes, are an appealing way to describe differentiation dynamics and developmental decisions. We consider the differentiation of neural and mesodermal cells from pluripotent mouse embryonic stem cells exposed to different combinations and durations of signalling factors. We developed a principled statistical approach using flow cytometry data to quantify differentiating cell states. Then, using a framework based on Catastrophe Theory and approximate Bayesian computation, we constructed the corresponding dynamical landscape. The result was a quantitative model that accurately predicted the proportions of neural and mesodermal cells differentiating in response to specific signalling regimes. Taken together, the approach we describe is broadly applicable for the quantitative analysis of differentiation dynamics and for determining the logic of developmental cell fate decisions.

Recent advances of single cell techniques catalyzed quantitative studies on the dynamics of cell phenotypic transitions (CPT) emerging as a new field. However, fixed cell-based approaches have fundamental limits on revealing temporal information, and fluorescence-based live cell imaging approaches are technically challenging for multiplex long-term imaging. To tackle the challenges, we developed an integrated experimental/computational platform for reconstructing single cell phenotypic transition dynamics. Experimentally, we developed a live-cell imaging platform to record the phenotypic transition path of A549 VIM-RFP reporter cell line and unveil parallel paths of epithelial-to-mesenchymal transition (EMT). Computationally, we modified a finite temperature string method to reconstruct the reaction coordinate from the paths, and reconstruct a corresponding quasi-potential, which reveals that the EMT process resembles a barrier-less relaxation process. Our work demonstrates the necessity of extracting dynamical information of phenotypic transitions and the existence of a unified theoretical framework describing transition and relaxation dynamics in systems with and without detailed balance.

Dictyostelium Discoideum (DD) are a fascinating single-cellular organism. When nutrients are plentiful, the DD cells act as autonomous individuals foraging their local vicinity. At the onset of starvation, a few (<0.1%) cells begin communicating with others by emitting a spike in the chemoattractant protein cyclic-AMP. Nearby cells sense the chemical gradient and respond by moving toward it and emitting a cyclic-AMP spike of their own. Cyclic-AMP activity increases over time, and eventually a spiral wave emerges, attracting hundreds of thousands of cells to an aggregation center. How DD cells go from autonomous individuals to a collective entity remains an open question for more than 60 years--a question whose answer would shed light on the emergence of multi-cellular life. Recently, trans-scale imaging has allowed the ability to sense the cyclic-AMP activity at both cell and colony levels. Using both the images as well as toy simulation models, this research aims to clarify whether the activity at the colony level is in fact initiated by a few cells, which may be deemed "leader" or "pacemaker" cells. In this talk, I will demonstrate the use of information-theoretic techniques to classify leaders and followers based on trajectory data, as well as to infer the domain of interaction of leader cells. We validate the techniques on toy models where leaders and followers are known, and then try to answer the question in real data--do leader cells drive collective behavior in DD colonies?

In vivo, the human genome folds into a characteristic ensemble of 3D structures. The mechanism driving the folding process remains unknown. A theoretical model for chromatin (the minimal chromatin model) explains the folding of interphase chromosomes and generates chromosome conformations consistent with experimental data is presented. The energy landscape of the model was derived by using the maximum entropy principle and relies on two experimentally derived inputs: a classification of loci into chromatin types and a catalog of the positions of chromatin loops. This model was generalized by utilizing a neural network to infer these chromatin types using epigenetic marks present at a locus, as assayed by ChIP-Seq. The ensemble of structures resulting from these simulations completely agree with HI-C data and exhibits unknotted chromosomes, phase separation of chromatin types, and a tendency for open chromatin to lie at the periphery of chromosome territories. Although this theoretical methodology was trained in one cell line, the human GM12878 lymphoblastoid cells, it has successfully predicted the structural ensembles of multiple human cell lines. Finally, going beyond Hi-C, our predicted structures are also consistent with microscopy measurements. Analysis of both structures from simulation and microscopy reveals that short segments of chromatin make two-state transitions between closed conformations and open dumbbell conformations. For gene active segments, the vast majority of genes appear clustered in the linker region of the chromatin segment, allowing us to speculate possible mechanisms by which chromatin structure and dynamics may be involved in controlling gene expression. * Supported by the NSF

Many microorganisms and cells function in complex (non-Newtonian) fluids, which are mixtures of different materials and exhibit both viscous and elastic stresses. For example, mammalian sperm swim through cervical mucus on their journey through the female reproductive tract, and they must penetrate the viscoelastic gel outside the ovum to fertilize. In micro-scale swimming the dynamics emerge from the coupled interactions between the complex rheology of the surrounding media and the passive and active body dynamics of the swimmer. We use computational models of swimmers in viscoelastic fluids to investigate and provide mechanistic explanations for emergent swimming behaviors. I will discuss how flexible filaments (such as flagella) can store energy from a viscoelastic fluid to gain stroke boosts due to fluid elasticity. I will also describe 3D simulations of model organisms such as C. Reinhardtii and mammalian sperm, where we use experimentally measured stroke data to separate naturally coupled stroke and fluid effects. We explore why strokes that are adapted to Newtonian fluid environments might not do well in viscoelastic environments.

Bacteria exhibit a bewildering diversity of morphologies, but despite their impact on nearly all aspects of life, they are frequently classified into a few general categories, usually just “spheres” and “rods.” Curved-rod bacteria are one simple variation observed in many environments, particularly the ocean. However, why so many species have evolved this shape is unknown. We used a regularized Stokeslet Boundary Element Method to model the motility of flagellated, curved bacteria. We show that curvature can increase swimming efficiency, revealing a widely applicable selective advantage. Furthermore, we show that the distribution of cell lengths and curvatures observed across bacteria in nature is predicted by evolutionary trade-offs between three tasks influenced by shape: efficient swimming, the ability to detect chemical gradients, and reduced cost of cell construction. We therefore reveal shape as an important component of microbial fitness.

The precise spatial localization of molecular signals within tissues richly informs the mechanisms of tissue formation and function. Here, we’ll introduce Slide-seq, a technology which enables transcriptome-wide measurements with near-single cell spatial resolution. We’ll describe recent experimental and computational advances to enable Slide-seq in biological contexts in biological contexts where high detection sensitivity is important. More broadly, we’ll discuss the promise and challenges of spatial transcriptomics for tissue genomics. Lastly, we’ll touch upon novel molecular recording technologies, which allows recording of the absolute time dynamics of gene expression in live systems into DNA sequences.

Human genes associated with brain-related diseases are being discovered at an accelerating pace. A major challenge is an identification of the mechanisms through which these genes act, and of potential therapeutic strategies. To elucidate such mechanisms in human cells, we established a CRISPR-based platform for genetic screening in human iPSC-derived neurons, astrocytes and microglia. Our approach relies on CRISPR interference (CRISPRi) and CRISPR activation (CRISPRa), in which a catalytically dead version of the bacterial Cas9 protein recruits transcriptional repressors or activators, respectively, to endogenous genes to control their expression, as directed by a small guide RNA (sgRNA). Complex libraries of sgRNAs enable us to conduct genome-wide or focused loss-of-function and gain-of-function screens. Such screens uncover molecular players for phenotypes based on survival, stress resistance, fluorescent phenotypes, high-content imaging and single-cell RNA-Seq. To uncover disease mechanisms and therapeutic targets, we are conducting genetic modifier screens for disease-relevant cellular phenotypes in patient-derived neurons and glia with familial mutations and isogenic controls. In a genome-wide screen, we have uncovered genes that modulate the formation of disease-associated aggregates of tau in neurons with a tauopathy-linked mutation (MAPT V337M). CRISPRi/a can also be used to model and functionally evaluate disease-associated changes in gene expression, such as those caused by eQTLs, haploinsufficiency, or disease states of brain cells. We will discuss an application to Alzheimer’s Disease-associated genes in microglia.