Publications

The first cell fate bifurcation in mammalian development directs cells toward either the trophectoderm (TE) or inner cell mass (ICM) compartments in pre-implantation embryos. This decision is regulated by the subcellular localization of a transcriptional co-activator YAP and takes place over several progressively asynchronous cleavage divisions. As a result of this asynchrony and variable arrangement of blastomeres, reconstructing the dynamics of the TE/ICM cell specification from fixed embryos is extremely challenging. To address this, we developed a live-imaging approach and applied it to measure pairwise dynamics of nuclear YAP and its direct target genes, CDX2 and SOX2, which are key transcription factors of the TE and ICM, respectively. Using these datasets, we constructed a generative model of the first cell fate bifurcation, which reveals the time-dependent statistics of the TE and ICM cell allocation. In addition to making testable predictions for the joint dynamics of the full YAP/CDX2/SOX2 motif, the model revealed the stochastic nature of the induction timing of the key cell fate determinants and identified the features of YAP dynamics that are necessary or sufficient for this induction. Notably, temporal heterogeneity was particularly prominent for SOX2 expression among ICM cells. As heterogeneities within the ICM have been linked to the initiation of the second cell fate decision in the embryo, understanding the origins of this variability is of key significance. The presented approach reveals the dynamics of the first cell fate choice and lays the groundwork for dissecting the next cell fate decisions in mouse development.

A fundamental problem in computer system performance, as well as in the natural sciences, concerns inferring from observations an understanding of the behavior of stochastic processes of interacting system components whose dynamics are driven by an unknown underlying stochastic differential equation (SDE). The objective in solving this problem is to infer the underlying equations of the dynamics of the system from sets of system measurements, indexed over time. Given the stochastic nature of such systems, together with a lack of information on stochastic trajectories in many cases [1, 3], this represents a very challenging problem in general.

A fundamental problem in computer system performance, as well as in the natural sciences, concerns inferring from observations an understanding of the behavior of stochastic processes of interacting system components whose dynamics are driven by an unknown underlying stochastic differential equation (SDE). The objective in solving this problem is to infer the underlying equations of the dynamics of the system from sets of system measurements, indexed over time. Given the stochastic nature of such systems, together with a lack of information on stochastic trajectories in many cases [1, 3], this represents a very challenging problem in general.

During the first cell fate decision in mammalian embryos, the inner cell mass cells, which will give rise to the embryo proper and other extraembryonic tissues, segregate from the trophectoderm cells, the precursors of the placenta. Cell fate segregation proceeds in a gradual manner encompassing two rounds of cell division, as well as cell positional and morphological changes. While it is known that the activity of the Hippo signaling pathway and the subcellular localization of its downstream effector YAP dictate lineage specific gene expression, the response of YAP to these dynamic cellular changes remains incompletely understood. Here we address these questions by quantitative live imaging of endogenously tagged YAP while simultaneously monitoring geometric cellular features and cell cycle progression throughout cell fate segregation. We apply a probabilistic model to our dynamic data, providing a quantitative characterization of the mutual effects of YAP and cellular relative exposed area, which has previously been shown to correlate with subcellular YAP localization in fixed samples. Additionally, we study how nuclear YAP levels are influenced by other factors, such as the decreasing pool of maternally provided YAP that is partitioned to daughter cells through cleavage divisions, cell cycle-associated nuclear volume changes, and a delay after divisions in adjusting YAP levels to new cell positions. Interestingly, we find that establishing low nuclear YAP levels required for the inner cell mass fate is largely achieved by passive cell cycle-associated mechanisms. Moreover, contrary to expectations, we find that mechanical perturbations that result in cell and nuclear shape changes do not influence YAP localization in the embryo. Together our work identifies how various inputs are integrated over a dynamic developmental time course to shape the levels of a key molecular determinant of the first cell fate choice.

The translation and shape deformations of a passive viscous Newtonian droplet immersed in an active nematic liquid crystal under circular confinement are analyzed using a linear stability analysis. We focus on the case of a sharply aligned active nematic in the limit of strong elastic relaxation in two dimensions. Using an active liquid crystal model, we employ the Lorentz reciprocal theorem for Stokes flow to study the growth of interfacial perturbations as a result of both active and elastic stresses. Instabilities are uncovered in both extensile and contractile systems, for which growth rates are calculated and presented in terms of the dimensionless ratios of active, elastic, and capillary stresses, as well as the viscosity ratio between the two fluids. We also extend our theory to analyze the inverse scenario, namely, the stability of an active nematic droplet surrounded by a passive viscous layer. Our results highlight the subtle interplay of capillary, active, elastic, and viscous stresses in governing droplet stability. The instabilities uncovered here may be relevant to a plethora of biological active systems, from the dynamics of passive droplets in bacterial suspensions to the organization of subcellular compartments inside the cell and the cell nucleus.

Cystic fibrosis (CF) is due to loss-of-function variants of the CF transmembrane conductance regulator (CFTR) channel. The most effective treatment for people with CF carrying the F508del mutation is the triple combination of elexacaftor–tezacaftor–ivacaftor (ETI). ETI can correct the underlying defect(s) in other CFTR mutants. The use of disease-relevant predictive models such as patient-derived human nasal epithelial cells allow to investigate the response to CFTR modulators of specific genotypes, possibly supporting patients' access to treatment.

Accurate prediction of drug-target binding and unbinding kinetics and thermodynamics is essential for guiding drug discovery and lead optimization. However, traditional atomistic simulations are often too computationally expensive to capture rare events that govern ligand (un)binding. Several enhanced sampling methods exist to overcome these limitations, but they require extensive manual intervention and introduce variability and artifacts in free energy and kinetic estimates that limit high-throughput scalability. The present work introduces seekrflow, an automated multiscale milestoning simulation pipeline that streamlines the entire workflow from a single receptor-ligand input structure to kinetic and thermodynamic predictions in a single step. This integrated approach minimizes manual intervention, reduces computational overhead, and enhances the reproducibility and accuracy of kinetic and thermodynamic predictions. The accuracy and efficiency of the pipeline is demonstrated on multiple receptor-ligand complexes, including inhibitors of heat shock protein 90, threonine-tyrosine kinase, and the trypsin protein, with predicted kinetic parameters closely matching experimental estimates. seekrflow establishes a new benchmark for automated and high-throughput physics-based predictions of kinetics and thermodynamics.

Accurate prediction of drug-target binding and unbinding kinetics and thermodynamics is essential for guiding drug discovery and lead optimization. However, traditional atomistic simulations are often too computationally expensive to capture rare events that govern ligand (un)binding. Several enhanced sampling methods exist to overcome these limitations, but they require extensive manual intervention and introduce variability and artifacts in free energy and kinetic estimates that limit high-throughput scalability. The present work introduces seekrflow, an automated multiscale milestoning simulation pipeline that streamlines the entire workflow from a single receptor-ligand input structure to kinetic and thermodynamic predictions in a single step. This integrated approach minimizes manual intervention, reduces computational overhead, and enhances the reproducibility and accuracy of kinetic and thermodynamic predictions. The accuracy and efficiency of the pipeline is demonstrated on multiple receptor-ligand complexes, including inhibitors of heat shock protein 90, threonine-tyrosine kinase, and the trypsin protein, with predicted kinetic parameters closely matching experimental estimates. seekrflow establishes a new benchmark for automated and high-throughput physics-based predictions of kinetics and thermodynamics.

Septins are a conserved family of cytoskeletal proteins known for sensing micron-scale membrane curvature via amphipathic helix (AH) domains. While cooperative interactions in septin assembly have been suggested, the molecular mechanisms governing membrane binding and assembly remain unclear. Building on prior findings, we use all-atom molecular dynamics simulations to examine how single and paired extended AH domains, derived from Cdc12, interact with lipid bilayers. A single membrane-bound AH adopts a curved conformation. In solution, a second AH peptide preferentially interacts with the bound peptide through conserved salt bridges, favoring an antiparallel arrangement. Simulations of covalently linked AH tandems confirm this configuration. Dual membrane-bound peptides induce lipid packing defects, reduce tail order, and exhibit slight membrane displacement, suggesting curved membranes may better accommodate multiple AH domains. Our findings advance the mechanistic understanding of septin-membrane interactions and highlight the role of cooperative AH domain binding in stabilizing higher-order structures.

Septins are a conserved family of cytoskeletal proteins known for sensing micron-scale membrane curvature via amphipathic helix (AH) domains. While cooperative interactions in septin assembly have been suggested, the molecular mechanisms governing membrane binding and assembly remain unclear. Building on prior findings, we use all-atom molecular dynamics simulations to examine how single and paired extended AH domains, derived from Cdc12, interact with lipid bilayers. A single membrane-bound AH adopts a curved conformation. In solution, a second AH peptide preferentially interacts with the bound peptide through conserved salt bridges, favoring an antiparallel arrangement. Simulations of covalently linked AH tandems confirm this configuration. Dual membrane-bound peptides induce lipid packing defects, reduce tail order, and exhibit slight membrane displacement, suggesting curved membranes may better accommodate multiple AH domains. Our findings advance the mechanistic understanding of septin-membrane interactions and highlight the role of cooperative AH domain binding in stabilizing higher-order structures.