Publications

To prepare gametes with the appropriate number of chromosomes, mammalian oocytes undergo two sequential cell divisions. During each division, a large, long-lived, microtubule-based organelle called the meiotic spindle assembles around condensed chromosomes. Although meiotic spindles have been intensively studied for several decades, as force-generating mechanical objects, they remain very poorly understood. In materials physics, coarse-grained theories have been essential in understanding the large-scale behavior of systems composed of many interacting particles. It is unclear, however, if this approach can succeed in capturing the properties of active, biochemically complex, living materials like the spindle. Here, we show that a class of models based on nematic liquid crystal theory can describe important aspects of the organelle-scale structure and dynamics of spindles in living mouse oocytes. Using our models to interpret quantitative polarization microscopy data, we measure for the first time material properties relating to stress propagation in living oocytes, including the nematic diffusivities corresponding to splay and bend deformations. Unlike the reconstituted amphibian spindles that were previously studied in vitro, nematic elastic stress is exponentially screened in the microtubule network of living mammalian oocytes, with a screening length of order one micron. This observation can be explained by the relatively high volume fraction of embedded chromosomes in mammalian meiotic spindles, which cause long voids in the microtubule network and so disrupt orientational stress propagation.

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

Well-conditioned boundary integral methods for the solution of elliptic boundary value problems (BVPs) are powerful tools for static and dynamic physical simulations. When there are many close-to-touching boundaries (e.g., in complex fluids) or when the solution is needed in the bulk, nearly singular integrals must be evaluated at many targets. We show that precomputing a linear map from surface density to an effective source representation renders this task highly efficient, in the common case where each object is “simple”, i.e., its smooth boundary needs only moderately many nodes. We present a kernel-independent method needing only an upsampled smooth surface quadrature, and one dense factorization, for each distinct shape. No (near-)singular quadrature rules are needed. The resulting effective sources are drop-in compatible with fast algorithms, with no local corrections nor bookkeeping. Our extensive numerical tests include 2D FMM-based Helmholtz and Stokes BVPs with up to 1000 objects (281000 unknowns), and a 3D Laplace BVP with 10 ellipsoids separated by 1/30 of a diameter. We include a rigorous analysis for analytic data in 2D and 3D.

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.