Publications

In the fruit fly Drosophila melanogaster, two cells in a cyst of 16 interconnected cells have the potential to become the oocyte, but only one of these will assume an oocyte fate as the cysts transition through regions 2a and 2b of the germarium. The mechanism of specification depends on a polarized microtubule network, a dynein dependent Egl:BicD mRNA cargo complex, a special membranous structure called the fusome and its associated proteins, and the translational regulator orb. In this work, we have investigated the role of orb and the fusome in oocyte specification. We show here that specification is a stepwise process. Initially, orb mRNAs accumulate in the two pro-oocytes in close association with the fusome. This association is accompanied by the activation of the orb autoregulatory loop, generating high levels of Orb. Subsequently, orb mRNAs become enriched in only one of the pro-oocytes, the presumptive oocyte, and this is followed, with a delay, by Orb localization to the oocyte. We find that fusome association of orb mRNAs is essential for oocyte specification in the germarium, is mediated by the orb 3′ UTR, and requires Orb protein. We also show that the microtubule minus end binding protein Patronin functions downstream of orb in oocyte specification. Finally, in contrast to a previously proposed model for oocyte selection, we find that the choice of which pro-oocyte becomes the oocyte does not seem to be predetermined by the amount of fusome material in these two cells, but instead depends upon a competition for orb gene products.

The Kidney Precision Medicine Project (KPMP) aims to create a kidney tissue atlas, define disease subgroups, and identify critical cells, pathways, and targets for novel therapies through molecular investigation of human kidney biopsies obtained from participants with AKI or CKD. We present the case of a 66-year-old woman with diabetic kidney disease who underwent a protocol KPMP kidney biopsy. Her clinical history included diabetes mellitus complicated by neuropathy and eye disease, increased insulin resistance, hypertension, albuminuria, and relatively preserved glomerular filtration rate (early CKD stage 3a). The patient's histopathology was consistent with diabetic nephropathy and arterial and arteriolar sclerosis. Three-dimensional, immunofluorescence imaging of the kidney biopsy specimen revealed extensive periglomerular neovascularization that was underestimated by standard histopathologic approaches. Spatial transcriptomics was performed to obtain gene expression signatures at discrete areas of the kidney biopsy. Gene expression in the areas of glomerular neovascularization revealed increased expression of genes involved in angiogenic signaling, proliferation, and survival of endothelial cells, as well as new vessel maturation and stability. This molecular correlation provides additional insights into the development of kidney disease in patients with diabetes and spotlights how novel molecular techniques used by the KPMP can supplement and enrich the histopathologic diagnosis obtained from a kidney biopsy.

Cytokinesis of animal and fungi cells depends crucially on the anillin scaffold proteins. Fission yeast anillin-related Mid1 anchors cytokinetic ring precursor nodes to the membrane. However, it is unclear if both of its Pleckstrin Homology (PH) and C2 C-terminal domains bind to the membrane as monomers or dimers, and if one domain plays a dominant role. We studied Mid1 membrane binding with all-atom molecular dynamics near a membrane with yeast-like lipid composition. In simulations with the full C terminal region started away from the membrane, Mid1 binds through the disordered L3 loop of C2 in a vertical orientation, with the PH away from the membrane. However, a configuration with both C2 and PH initially bound to the membrane remains associated with the membrane. Simulations of C2-PH dimers show extensive asymmetric membrane contacts. These multiple modes of binding may reflect Mid1’s multiple interactions with membranes, node proteins, and ability to sustain mechanical forces.

Life in complex systems, such as cities and organisms, comes to a standstill when global coordination of mass, energy and information flows is disrupted. Global coordination is no less important in single cells, especially in large oocytes and newly formed embryos, which commonly use fast fluid flows for dynamic reorganization of their cytoplasm. These cytoplasmic streaming flows have been proposed to spontaneously arise from hydrodynamic interactions among cortically anchored microtubules loaded with cargo-carrying molecular motors. Here, we combine modelling and simulation with live imaging to investigate such flows in the Drosophila oocyte. Using a fast, accurate and scalable numerical approach to investigate fluid–structure interactions of thousands of flexible fibres, we demonstrate the robust emergence and evolution of cell-spanning vortices—or twisters—in three-dimensional cellular geometries. These twister flows, dominated by a near-rigid-body rotation with secondary toroidal components, reproduce the variety of experimental observations. In cells, these flows are probably involved in rapid mixing and transport of ooplasmic components.

Single same cell RNAseq/ATACseq multiome data provide unparalleled potential to develop high resolution maps of the cell-type specific transcriptional regulatory circuitry underlying gene expression. We present CREMA, a framework that recovers the full cis-regulatory circuitry by modeling gene expression and chromatin activity in individual cells without peak-calling or cell type labeling constraints. We demonstrate that CREMA overcomes the limitations of existing methods that fail to identify about half of functional regulatory elements which are outside the called chromatin ‘peaks’. These circuit sites outside called peaks are shown to be important cell type specific functional regulatory loci, sufficient to distinguish individual cell types. Analysis of mouse pituitary data identifies a Gata2-circuit for the gonadotrope-enriched disease-associated Pcsk1 gene, which is experimentally validated by reduced gonadotrope expression in a gonadotrope conditional Gata2-knockout model. We present a web accessible human immune cell regulatory circuit resource, and provide CREMA as an R package.

Fluid-structure interactions between active and passive components are important for many biological systems to function. A particular example is chromatin in the cell nucleus, where ATP-powered processes drive coherent motions of the chromatin fiber over micron lengths. Motivated by this system, we develop a multiscale model of a long flexible polymer immersed in a suspension of active force dipoles as an analog to a chromatin fiber in an active fluid—the nucleoplasm. Linear analysis identifies an orientational instability driven by hydrodynamic and alignment interactions between the fiber and the suspension, and numerical simulations show activity can drive coherent motions and structured conformations. These results demonstrate how active and passive components, connected through fluid-structure interactions, can generate coherent structures and self-organize on large scales.

Allostery is the mechanism by which information and control are propagated in biomolecules. It regulates ligand binding, chemical reactions, and conformational changes. An increasing level of experimental resolution and control over allosteric mechanisms promises a deeper understanding of the molecular basis for life and powerful new therapeutics. In this review, we survey the literature for an up-to-date biological and theoretical understanding of protein allostery. By delineating five ways in which the energy landscape or the kinetics of a system may change to give rise to allostery, we aim to help the reader grasp its physical origins. To illustrate this framework, we examine three systems that display these forms of allostery: allosteric inhibitors of beta-lactamases, thermosensation of TRP channels, and the role of kinetic allostery in the function of kinases. Finally, we summarize the growing power of computational tools available to investigate the different forms of allostery presented in this review.

In this work, we propose and compare three numerical methods to handle the one-phase Hele-Shaw problem with surface tension in dimension two by using three variational approaches in the spirit of the seminal works \cite{Otto, Gia_Otto}.

There is a need to define regions of gene activation or repression that control human kidney cells in states of health, injury, and repair to understand the molecular pathogenesis of kidney disease and design therapeutic strategies. Comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measure dual single nucleus RNA expression and chromatin accessibility, DNA methylation, and H3K27ac, H3K4me1, H3K4me3, and H3K27me3 histone modifications to decipher the chromatin landscape and gene regulation of the kidney in reference and adaptive injury states. We establish a spatially-anchored epigenomic atlas to define the kidney’s active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we note distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3, KLF6, and KLF10 regulates the transition between health and injury, while in thick ascending limb cells this transition is regulated by NR2F1. Further, combined perturbation of ELF3, KLF6, and KLF10 distinguishes two adaptive proximal tubular cell subtypes, one of which manifested a repair trajectory after knockout. This atlas will serve as a foundation to facilitate targeted cell-specific therapeutics by reprogramming gene regulatory networks.

The self-assembly of amyloid-beta (Aβ) peptides into fibrillar structures in the brain is a signature of Alzheimer's disease. Recent studies have reported correlations between Alzheimer's disease and type-2 diabetes. Structurally, hyperglycemia induces covalent protein crosslinkings by advanced glycation end products (AGE), which can affect the stability of Aβ oligomers. In this work, we leverage physics-based coarse-grained molecular simulations to probe alternate thermodynamic pathways that affect peptide aggregation propensities at varying concentrations of glucose molecules. Similar to previous experimental reports, our simulations show a glucose concentration-dependent increase in Aβ aggregation rates, without changes in the overall secondary structure content. We discovered that glucose molecules prefer partitioning onto the aggregate–water interface at a specific orientation, resulting in a loss of molecular rotational entropy. This effectively hastens the aggregation rates, as peptide self-assembly can reduce the available surface area for peptide–glucose interactions. This work introduces a new thermodynamic-driven pathway, beyond chemical cross-linking, that can modulate Aβ aggregation.