Publications

Male sex is a major risk factor for SARS-CoV-2 infection severity. To understand the basis for this sex difference, we studied SARS-CoV-2 infection in a young adult cohort of United States Marine recruits. Among 2,641 male and 244 female unvaccinated and seronegative recruits studied longitudinally, SARS-CoV-2 infections occurred in 1,033 males and 137 females. We identified sex differences in symptoms, viral load, blood transcriptome, RNA splicing, and proteomic signatures. Females had higher pre-infection expression of antiviral interferon-stimulated gene (ISG) programs. Causal mediation analysis implicated ISG differences in number of symptoms, levels of ISGs, and differential splicing of CD45 lymphocyte phosphatase during infection. Our results indicate that the antiviral innate immunity set point causally contributes to sex differences in response to SARS-CoV-2 infection. A record of this paper’s transparent peer review process is included in the supplemental information.

Molecular docking, which aims to find the most stable interacting configuration of a set of molecules, is of critical importance to drug discovery. Although a considerable number of classical algorithms have been developed to carry out molecular docking, most focus on the limiting case of docking two molecules. Since the number of possible configurations of N molecules is exponential in N, those exceptions which permit docking of more than two molecules scale poorly, requiring exponential resources to find high-quality solutions. Here, we introduce a one-hot encoded quadratic unconstrained binary optimization formulation (QUBO) of the multibody molecular docking problem, which is suitable for solution by quantum annealer. Our approach involves a classical pre-computation of pairwise interactions, which scales only quadratically in the number of bodies while permitting well-vetted scoring functions like the Rosetta REF2015 energy function to be used. In a second step, we use the quantum annealer to sample low-energy docked configurations efficiently, considering all possible docked configurations simultaneously through quantum superposition. We show that we are able to minimize the time needed to find diverse low-energy docked configurations by tuning the strength of the penalty used to enforce the one-hot encoding, demonstrating a 3-4 fold improvement in solution quality and diversity over performance achieved with conventional penalty strengths. By mapping the configurational search to a form compatible with current- and future-generation quantum annealers, this work provides an alternative means of solving multibody docking problems that may prove to have performance advantages for large problems, potentially circumventing the exponential scaling of classical approaches and permitting a much more efficient solution to a problem central to drug discovery and validation pipelines.

The cell nucleus is enveloped by a complex membrane, whose wrinkling has been implicated in disease and cellular aging. The biophysical dynamics and spectral evolution of nuclear wrinkling during multicellular development remain poorly understood due to a lack of direct quantitative measurements. Here, we combine live-imaging experiments, theory, and simulations to characterize the onset and dynamics of nuclear wrinkling during egg development in the fruit fly, Drosophila melanogaster, when nurse cell nuclei increase in size and display stereotypical wrinkling behavior. A spectral analysis of three-dimensional high-resolution data from several hundred nuclei reveals a robust asymptotic power-law scaling of angular fluctuations consistent with renormalization and scaling predictions from a nonlinear elastic shell model. We further demonstrate that nuclear wrinkling can be reversed through osmotic shock and suppressed by microtubule disruption, providing tunable physical and biological control parameters for probing mechanical properties of the nuclear envelope. Our findings advance the biophysical understanding of nuclear membrane fluctuations during early multicellular development.

How does breaking the symmetry of an equation alter the symmetry of its solutions? Here, we systematically examine how reducing underlying symmetries from spherical to axisymmetric influences the dynamics of an archetypal model of cell polarization, a key process of biological spatial self-organization. Cell polarization is characterized by nonlinear and non-local dynamics, but we overcome the theory challenges these traits pose by introducing a broadly applicable numerical scheme allowing us to efficiently study continuum models in a wide range of geometries. Guided by numerical results, we discover a dynamical hierarchy of timescales that allows us to reduce relaxation to a purely geometric problem of area-preserving geodesic curvature flow. Through application of variational results, we analytically construct steady states on a number of biologically relevant shapes. In doing so, we reveal non-trivial solutions for symmetry breaking.

We use computational design coupled with experimental characterization to systematically investigate the design principles for macrocycle membrane permeability and oral bioavailability. We designed 184 6–12 residue macrocycles with a wide range of predicted structures containing noncanonical backbone modifications and experimentally determined structures of 35; 29 are very close to the computational models. With such control, we show that membrane permeability can be systematically achieved by ensuring all amide (NH) groups are engaged in internal hydrogen bonding interactions. 84 designs over the 6–12 residue size range cross membranes with an apparent permeability greater than 1 × 10−6 cm/s. Designs with exposed NH groups can be made membrane permeable through the design of an alternative isoenergetic fully hydrogen-bonded state favored in the lipid membrane. The ability to robustly design membrane-permeable and orally bioavailable peptides with high structural accuracy should contribute to the next generation of designed macrocycle therapeutic

Numerous engineered and natural systems form through reinforcement and stabilization of a deformed configuration that was generated by a transient force. An important class of such structures arises during gametogenesis, when a dividing cell undergoes incomplete cytokinesis, giving rise to daughter cells that remain connected through a stabilized intercellular bridge (ICB). ICBs can form through arrest of the contractile cytokinetic furrow and its subsequent stabilization. Despite knowledge of the molecular components, the mechanics underlying robust ICB assembly and the interplay between ring contractility and stiffening are poorly understood. Here, we report joint experimental and theoretical work that explores the physics underlying robust ICB assembly. We develop a continuum mechanics model that reveals the minimal requirements for the formation of stable ICBs, and validate the model’s equilibrium predictions through a tabletop experimental analog. With insight into the equilibrium states, we turn to the dynamics: we demonstrate that contractility and stiffening are in dynamic competition and that the time intervals of their action must overlap to ensure assembly of ICBs of biologically observed proportions. Our results highlight a mechanism in which deformation and remodeling are tightly coordinated—one that is applicable to several mechanics-based applications and is a common theme in biological systems spanning several length scales.

We re-examine the model proposed by Alexander et al. (Phys. Fluids, vol. 18, 2006, 062103) for the closing of a circular hole in a molecularly thin incompressible Langmuir film situated on a Stokesian subfluid. For simplicity their model assumes that the surface phase is inviscid which leads to the result that the cavity area decreases at a constant rate determined by the ratio of edge tension to subfluid viscosity. We reformulate the problem, allowing for a regularising monolayer viscosity. The viscosity-dependent corrections to the hole dynamics are analysed and found to be non-trivial, even when the monolayer viscosity is small; these corrections may explain the departure of experimental data from the theoretical prediction when the hole radius becomes comparable to the Saffman–Delbrück length. Through fitting, under these relaxed assumptions, we find the edge tension could be as much as six times larger ( ∼
4.0 pN) than reported previously.

Small cell clusters exhibit numerous phenomena typically associated with complex systems, such as division of labour and programmed cell death. A conserved class of such clusters occurs during oogenesis in the form of germline cysts that give rise to oocytes. Germline cysts form through cell divisions with incomplete cytokinesis, leaving cells intimately connected through intercellular bridges that facilitate cyst generation, cell fate determination and collective growth dynamics. Using the well-characterized Drosophila melanogaster female germline cyst as a foundation, we present mathematical models rooted in the dynamics of cell cycle proteins and their interactions to explain the generation of germline cell lineage trees (CLTs) and highlight the diversity of observed CLT sizes and topologies across species. We analyse competing models of symmetry breaking in CLTs to rationalize the observed dynamics and robustness of oocyte fate specification, and highlight remaining gaps in knowledge. We also explore how CLT topology affects cell cycle dynamics and synchronization and highlight mechanisms of intercellular coupling that underlie the observed collective growth patterns during oogenesis. Throughout, we point to similarities across organisms that warrant further investigation and comment on the extent to which experimental and theoretical findings made in model systems extend to other species.

To demonstrate a clear link between predicted blood shear forces during valve closure and thrombogenicity that explains the thrombogenic difference between tissue and mechanical valves and provides a practical metric to develop and refine prosthetic valve designs for reduced thrombogenicity.

Cytosolic Ca2+ is a highly dynamic, tightly regulated and broadly conserved cellular signal. Ca2+ dynamics have been studied widely in cellular monocultures, yet organs in vivo comprise heterogeneous populations of stem and differentiated cells. Here, we examine Ca2+ dynamics in the adult Drosophila intestine, a self-renewing epithelial organ in which stem cells continuously produce daughters that differentiate into either enteroendocrine cells or enterocytes. Live imaging of whole organs ex vivo reveals that stem-cell daughters adopt strikingly distinct patterns of Ca2+ oscillations after differentiation: enteroendocrine cells exhibit single-cell Ca2+ oscillations, whereas enterocytes exhibit rhythmic, long-range Ca2+ waves. These multicellular waves do not propagate through immature progenitors (stem cells and enteroblasts), of which the oscillation frequency is approximately half that of enteroendocrine cells. Organ-scale inhibition of gap junctions eliminates Ca2+ oscillations in all cell types – even, intriguingly, in progenitor and enteroendocrine cells that are surrounded only by enterocytes. Our findings establish that cells adopt fate-specific modes of Ca2+ dynamics as they terminally differentiate and reveal that the oscillatory dynamics of different cell types in a single, coherent epithelium are paced independently.