Publications

Tumor-initiating cells (TICs) share features and regulatory pathways with normal stem cells, yet how the stem cell niche contributes to tumorigenesis remains unclear. Here, we identify CXCR4+ macrophages as a niche population enriched in normal mammary ducts, where they promote the regenerative activity of basal cells in response to luminal cell-derived CXCL12. CXCL12 triggers AKT-mediated stabilization of β-catenin, which induces Wnt ligands and pro-migratory genes, enabling intraductal macrophage infiltration and supporting regenerative activity of basal cells. Notably, these same CXCR4+ niche macrophages regulate the tumor-initiating activity of various breast cancer subtypes by enhancing TIC survival and tumor-forming capacity, while promoting early immune evasion through regulatory T cell induction. Furthermore, a CXCR4+ niche macrophage gene signature correlates with poor prognosis in human breast cancer. These findings highlight the pivotal role of the CXCL12-CXCR4 axis in orchestrating interactions between niche macrophages, mammary epithelial cells, and immune cells, thereby establishing a supportive niche for both normal tissue regeneration and mammary tumor initiation.

Suspensions of swimming particles exhibit complex collective behaviors driven by hydrodynamic interactions, showing persistent large-scale flows and long-range correlations. While heavily studied, it remains unclear how such structures depend on the system size and swimmer concentration. To address these issues, we simulate very large systems of suspended swimmers across a range of system sizes and volume fractions. For this we use high-performance simulation tools that build on slender body theory and implicit resolution of steric interactions. At low volume fractions and long times, the particle simulations reveal dynamic flow structures and correlation functions that scale with the system size. These results are consistent with a mean-field limit and agree well with a corresponding kinetic theory. At higher concentrations, the system departs from mean-field behavior. Flow structures become cellular, and correlation lengths scale with the particle size. Here, translational motion is suppressed, while rotational dynamics dominate. These findings highlight the limitations of dilute mean-field models and reveal new behaviors in dense active suspensions.

We examine theoretically the flow interactions and forward flight dynamics of tandem or in-line flapping wings. Two wings are driven vertically with prescribed heaving-and-plunging motions, and the horizontal propulsion speeds and positions are dynamically selected through aero- or hydro-dynamic interactions. Our simulations employ an improved vortex sheet method to solve for the locomotion of the pair within the collective flow field, and we identify 'schooling states' in which the wings travel together with nearly constant separation. Multiple terminal configurations are achieved by varying the initial conditions, and the emergent separations are approximately integer multiples of the wavelength traced out by each wing. We explain the stability of these states by perturbing the follower and mapping out an effective potential for its position in the leader's wake. Each equilibrium position is stabilized since smaller separations are associated with in-phase follower-wake motions that constructively reinforce the flow but lead to decreased thrust on the follower; larger separations are associated with antagonistic follower-wake motions, increased thrust, and a weakened collective wake. The equilibria and their stability are also corroborated by a linearized theory for the motion of the leader, the wake it produces, and its effect on the follower. We also consider a weakly-flapping follower driven with lower heaving amplitude than the leader. We identify 'keep-up' conditions for which the wings may still 'school' together despite their dissimilar kinematics, with the 'freeloading' follower passively assuming a favorable position within the wake that permits it to travel significantly faster than it would in isolation.

Time-resolved single-cell omics data offers high-throughput, genome-wide measurements of cellular states, which are instrumental to reverse-engineer the processes underpinning cell fate. Such technologies are inherently destructive, allowing only cross-sectional measurements of the underlying stochastic dynamical system. Furthermore, cells may divide or die in addition to changing their molecular state. Collectively these present a major challenge to inferring realistic biophysical models. We present a novel approach, \emph{unbalanced} probability flow inference, that addresses this challenge for biological processes modelled as stochastic dynamics with growth. By leveraging a Lagrangian formulation of the Fokker-Planck equation, our method accurately disentangles drift from intrinsic noise and growth. We showcase the applicability of our approach through evaluation on a range of simulated and real single-cell RNA-seq datasets. Comparing to several existing methods, we find our method achieves higher accuracy while enjoying a simple two-step training scheme.

Cells must limit RNA–RNA interactions to avoid irreversible RNA entanglement. Cells may prevent deleterious RNA-RNA interactions by genome organization to avoid complementarity however, RNA viruses generate long, perfectly complementary antisense RNA during replication. How do viral RNAs avoid irreversible entanglement? One possibility is RNA sequestration into biomolecular condensates. To test this, we reconstituted critical SARS-CoV-2 RNA–RNA interactions in Nucleocapsid condensates. We observed that RNAs with low propensity RNA–RNA interactions resulted in more round, liquid-like condensates while those with high sequence complementarity resulted in more heterogeneous networked morphology independent of RNA structure stability. Residue-resolution molecular simulations and direct sequencing-based detection of RNA–RNA interactions support that these properties arise from degree of trans RNA contacts. We propose that extensive RNA–RNA interactions in cell and viral replication are controlled via a combination of genome organization, timing, RNA sequence content, RNA production ratios, and emergent biomolecular condensate material properties.

This work probes the role of cell geometry in orienting self-organized fluid flows in the late-stage Drosophila oocyte. Recent theoretical work has shown that a model, which relies only on hydrodynamic interactions of flexible, cortically anchored microtubules and the mechanical loads from molecular motors moving upon them, is sufficient to generate observed flows. While the emergence of flows has been studied in spheres, oocytes change shape during streaming, and it was unclear how robust these flows are to the geometry of the cell. Here we use biophysical theory and computational analysis to investigate the role of geometry and find that the axis of rotation is set by the shape of the domain and that the flow is robust to biologically relevant perturbations of the domain shape. Using live imaging and three-dimensional flow reconstruction, we test the predictions of the theory/simulation, finding consistency between the model and live experiments, further demonstrating a geometric dependence on flow direction in late-stage Drosophila oocytes.

Huntington's disease (HD) is an inherited neurodegenerative disorder associated with motor and cognitive decline, caused by a mutation in the poly-glutamine (polyQ) region near the N-terminus of the huntingtin (htt) protein. Expansion of the polyQ region results in the disease that is characterized by oligomeric and fibrillar aggregates of mutated protein. The first 17 amino acids (Nt17) of htt, which are adjacent to the polyQ tract, function as a lipid-binding domain, facilitated by the formation of an amphipathic α-helix. There is increasing evidence that lipid interactions may play a role in the toxic gain of function associated with the htt polyQ expansion, as membrane-related changes, including structural abnormalities of several organelles, are observed in HD. Given the uneven and curved shapes of organelles, it is important to examine the mechanistic preferences that drive the preferential partitioning of Nt17 to curved membranes. To better understand the role of the cell membrane environment in the interaction and aggregation of htt, circular dichroism, fluorescence microscopy, and coarse-grained molecular dynamics were employed to measure the association of Nt17 with phospholipid vesicles and subsequent effects throughout time. In zwitterionic curved membranes, sensing was driven by the bulky sidechains of phenylalanine residues, which are able to sense lipid packing defects in the curved regions of the membrane. However, in a mixture of zwitterionic and anionic lipids, curvature sensing is affected by the anionic lipid content, implying the surface charge of membranes affects the curvature sensing process. Salt screening experiments suggest a balance between the electrostatic and hydrophobic interactions that governs the extent to which Nt17 can sense physiologically relevant regions of curvature.

Rational computational design is crucial to the pursuit of novel drugs and therapeutic agents. Meso-scale cyclic peptides, which consist of 7-40 amino acid residues, are of particular interest due to their conformational rigidity, binding specificity, degradation resistance, and potential cell permeability. Because there are few natural cyclic peptides, de novo design involving non-canonical amino acids is a potentially useful goal. Here, we develop an efficient pipeline (CyclicChamp) for cyclic peptide design. After converting the cyclic constraint into an error function, we employ a variant of simulated annealing to search for low-energy peptide backbones while maintaining peptide closure. Compared to the previous random sampling approach, which was capable of sampling conformations of cyclic peptides of up to 14 residues, our method both greatly accelerates the computation speed for sampling conformations of small macrocycles (ca. 7 residues), and addresses the high-dimensionality challenge that large macrocycle designs often encounter. As a result, CyclicChamp makes conformational sampling tractable for 15- to 24-residue cyclic peptides, thus permitting the design of macrocycles in this size range. Microsecond-length molecular dynamics simulations on the resulting 15, 20, and 24 amino acid cyclic designs identify designs with kinetic stability. To test their thermodynamic stability, we perform additional replica exchange molecular dynamics simulations and generate free energy surfaces. Three 15-residue designs, one 20-residue and one 24-residue design emerge as promising candidates.

Amphipathic helices (AHs) are secondary structures that can facilitate binding of proteins to the membrane by folding into a helix with hydrophobic and hydrophilic faces that interact with the same surfaces in the lipid membrane. Septins are cytoskeletal proteins that preferentially bind to domains of micron-scale curvature on the cell membrane. Studies have shown that AH domains in septin are essential for curvature sensing. We present the first computational study of septin AH interactions with lipid bilayers. Using all-atom simulations and metadynamics-enhanced sampling, we study the effect of charge distribution at the flanking ends of septin AH on the energy for helical folding and its consequences on the binding configuration and affinity to the membrane. This is relevant to septins, since the net positive charge on the flanking C-terminal amino acids is a conserved property across several organisms. Simulations revealed that the energy barrier for folding in the neutral-capped AH is much larger than the charge-capped AH, leading to a small fraction of AH folding and integration to the membrane compared to a significantly folded configuration in the bound charge-capped AH. These observations are consistent with the binding measurements of synthetic AH constructs with variable helicity to lipid vesicles. Additionally, we examined an extended AH sequence including eight amino acids upstream and downstream of the AH to mimic the native protein. Again, simulations and experiments show that the extended peptide, with a net positive charge at C-terminus, adopts a strong helical configuration in solution, giving rise to a higher membrane affinity. Altogether, these results identify the energy cost for folding of AHs as a regulator of AH binding configuration and affinity and provide a basic template for parameterizing AH-membrane interactions as a starting point for the future multiscale simulations for septin-membrane interactions.

Germline-soma segregation is crucial for fertility. Primordial germ cells (PGCs) arise early in development and are the very first cells to form in the Drosophila embryo. At the time of PGC formation, the embryo is a syncytium where nuclei divide within a common cytoplasm. Whereas invaginating plasma membrane furrows enclose nuclei to form somatic lineages during the 14th nuclear division cycle, PGCs emerge from the syncytium during the 9th division cycle in a mechanistically distinct process. PGC formation depends on maternally deposited germ granules localized at the embryo’s posterior pole. Germ granules trigger protrusion of membrane buds that enlarge to surround several nuclei that reach the posterior pole. Buds are remodeled to cells through mitotic division and constriction of the bud neck. Previous studies implicated F-actin,1 actin regulators,2,3 and contractile ring components4 in mitotic furrow formation, but what drives bud emergence and how germ granules provoke reshaping of the plasma membrane remain unknown. Here, we investigate the mechanism of germ-granule-induced bud formation. Treating the embryo as a pressurized elastic shell, we used mathematical modeling to examine possible mechanical mechanisms for local membrane protrusion. One mechanism, outward buckling produced by polymerization of a branched F-actin network, is supported by experimental data. Further, we show that germ granules modify membrane lipid composition, promoting local branched F-actin polymerization that initiates PGC formation. We propose that a mechanism for membrane lipid regulation of F-actin dynamics in migrating cells has been adapted for PGC formation in response to spatial cues provided by germ granules.