Publications

Self-diffusion in dense rod suspensions are subject to strong geometric constraints because of steric interactions. This topological effect is essentially anisotropic when rods are nematically aligned with their neighbors, raising considerable challenges in understanding and analyzing their impacts on the bulk physical properties. Via a classical Doi-Onsager kinetic model with the Maier-Saupe potential, we characterize the long-time rotational Brownian diffusivity for dense suspensions of hard rods of finite aspect ratios, based on quadratic orientation autocorrelation functions. Furthermore, we show that the computed nonmonotonic scalings of the diffusivity as a function of volume fraction can be accurately predicted by an alternative tube model in the nematic phase.

The famous Arrhenius equation is well motivated to describe the temperature dependence of chemical reactions but has also been used for complicated biological processes. Here, we evaluate how well the simple Arrhenius equation predicts complex multistep biological processes, using frog and fruit fly embryogenesis as two canonical models. We find the Arrhenius equation provides a good approximation for the temperature dependence of embryogenesis, even though individual developmental stages scale differently with temperature. At low and high temperatures, however, we observed significant departures from idealized Arrhenius Law behavior. When we model multistep reactions of idealized chemical networks we are unable to generate comparable deviations from linearity. In contrast, we find the single enzyme GAPDH shows non-linearity in the Arrhenius plot similar to our observations of embryonic development. Thus, we find that complex embryonic development can be well approximated by the simple Arrhenius Law and propose that the observed departure from this law results primarily from non-idealized individual steps rather than the complexity of the system.

The double-helical topology of DNA molecules observed at room temperature in the absence of any external loads can be disrupted by increasing the bath temperature or by applying tensile forces, leading to spontaneous strand separation known as DNA melting. Here, continuum mechanics of a 2D birod is combined with statistical mechanics to formulate a unified framework for studying both thermal melting and tensile force induced melting of double-stranded molecules: it predicts the variation of melting temperature with tensile load, provides a mechanics-based understanding of the cooperativity observed in melting transitions, and reveals an interplay between solution electrostatics and micromechanical deformations of DNA which manifests itself as an increase in the melting temperature with increasing ion concentration. This novel predictive framework sheds light on the micromechanical aspects of DNA melting and predicts trends that were observed experimentally or extracted phenomenologically using the Clayperon equation.

Current estimates of COVID-19 prevalence are largely based on symptomatic, clinically diagnosed cases. The existence of a large number of undiagnosed infections hampers population-wide investigation of viral circulation. Here, we use longitudinal wastewater analysis to track SARS-CoV-2 dynamics in wastewater at a major urban wastewater treatment facility in Massachusetts, between early January and May 2020. SARS-CoV-2 was first detected in wastewater on March 3. Viral titers in wastewater increased exponentially from mid-March to mid-April, after which they began to decline. Viral titers in wastewater correlated with clinically diagnosed new COVID-19 cases, with the trends appearing 4-10 days earlier in wastewater than in clinical data. We inferred viral shedding dynamics by modeling wastewater viral titers as a convolution of back-dated new clinical cases with the viral shedding function of an individual. The inferred viral shedding function showed an early peak, likely before symptom onset and clinical diagnosis, consistent with emerging clinical and experimental evidence. Finally, we found that wastewater viral titers at the neighborhood level correlate better with demographic variables than with population size. This work suggests that longitudinal wastewater analysis can be used to identify trends in disease transmission in advance of clinical case reporting, and may shed light on infection characteristics that are difficult to capture in clinical investigations, such as early viral shedding dynamics.

We demonstrate that an achiral stretching force transforms disk-shaped colloidal membranes composed of chiral rods into twisted ribbons with handedness opposite the preferred twist of the rods. Using an experimental technique that enforces torque-free boundary conditions we simultaneously measure the force-extension curve and the ribbon shape. An effective theory that accounts for the membrane bending energy and uses geometric properties of the edge to model the internal liquid crystalline degrees of freedom explains both the measured force-extension curve and the force-induced twisted shape.

A genetic etiology is identified for one-third of patients with congenital heart disease (CHD), with 8% of cases attributable to coding de novo variants (DNVs). To assess the contribution of noncoding DNVs to CHD, we compared genome sequences from 749 CHD probands and their parents with those from 1,611 unaffected trios. Neural network prediction of noncoding DNV transcriptional impact identified a burden of DNVs in individuals with CHD (n = 2,238 DNVs) compared to controls (n = 4,177; P = 8.7 × 10−4). Independent analyses of enhancers showed an excess of DNVs in associated genes (27 genes versus 3.7 expected, P = 1 × 10−5). We observed significant overlap between these transcription-based approaches (odds ratio (OR) = 2.5, 95% confidence interval (CI) 1.1–5.0, P = 5.4 × 10−3). CHD DNVs altered transcription levels in 5 of 31 enhancers assayed. Finally, we observed a DNV burden in RNA-binding-protein regulatory sites (OR = 1.13, 95% CI 1.1–1.2, P = 8.8 × 10−5). Our findings demonstrate an enrichment of potentially disruptive regulatory noncoding DNVs in a fraction of CHD at least as high as that observed for damaging coding DNVs.

Microbial diversity in the cystic fibrosis (CF) lung decreases over decades as pathogenic bacteria such as Pseudomonas aeruginosa take over. The dynamics of the CF microbiome and metabolome over shorter time frames, however, remain poorly studied. Here, we analyze paired microbiome and metabolome data from 594 sputum samples collected over 401 days from six adult CF subjects (subject mean = 179 days) through periods of clinical stability and 11 CF pulmonary exacerbations (CFPE). While microbiome profiles were personalized (permutational multivariate analysis of variance [PERMANOVA] r2 = 0.79, P < 0.001), we observed significant intraindividual temporal variation that was highest during clinical stability (linear mixed-effects [LME] model, P = 0.002). This included periods where the microbiomes of different subjects became highly similar (UniFrac distance, <0.05). There was a linear increase in the microbiome alpha-diversity and in the log ratio of anaerobes to pathogens with time (n = 14 days) during the development of a CFPE (LME P = 0.0045 and P = 0.029, respectively). Collectively, comparing samples across disease states showed there was a reduction of these two measures during antibiotic treatment (LME P = 0.0096 and P = 0.014, respectively), but the stability data and CFPE data were not significantly different from each other. Metabolome alpha-diversity was higher during CFPE than during stability (LME P = 0.0085), but no consistent metabolite signatures of CFPE across subjects were identified. Virulence-associated metabolites from P. aeruginosa were temporally dynamic but were not associated with any disease state. One subject died during the collection period, enabling a detailed look at changes in the 194 days prior to death. This subject had over 90% Pseudomonas in the microbiome at the beginning of sampling, and that level gradually increased to over 99% prior to death. This study revealed that the CF microbiome and metabolome of some subjects are dynamic through time. Future work is needed to understand what drives these temporal dynamics and if reduction of anaerobes correlate to clinical response to CFPE therapy.

How bulk cytoplasm generates forces to separate post-anaphase microtubule (MT) asters in Xenopus laevis and other large eggs remains unclear. Previous models proposed dynein-based organelle transport generates length-dependent forces on astral MTs that pull centrosomes through the cytoplasm, away from the midplane. In Xenopus egg extracts, we co-imaged MTs, endoplasmic reticulum (ER), mitochondria, acidic organelles, F-actin, keratin, and fluorescein in moving and stationary asters. In asters that were moving in response to dynein and actomyosin forces, we observed that all cytoplasmic components moved together, i.e., as a continuum. Dynein-mediated organelle transport was restricted by interior MTs and F-actin. Organelles exhibited a burst of dynein-dependent inward movement at the growing aster surface, then mostly halted inside the aster. Dynein-coated beads were slowed by F-actin, but in contrast to organelles, beads did not halt inside asters. These observations call for new models of aster positioning based on surface forces and internal stresses.

A mechanistic understanding of the SARS-CoV-2 viral replication cycle is essential to develop new therapies for the COVID-19 global health crisis. In this study, we show that the SARS-CoV-2 nucleocapsid protein (N-protein) undergoes liquid-liquid phase separation (LLPS) with the viral genome, and propose a model of viral packaging through LLPS. N-protein condenses with specific RNA sequences in the first 1000 nts (5’-End) under physiological conditions and is enhanced at human upper airway temperatures. N-protein condensates exclude non-packaged RNA sequences. We comprehensively map sites bound by N-protein in the 5’-End and find preferences for single-stranded RNA flanked by stable structured elements. Liquid-like N-protein condensates form in mammalian cells in a concentration-dependent manner and can be altered by small molecules. Condensation of N-protein is sequence and structure specific, sensitive to human body temperature, and manipulatable with small molecules thus presenting screenable processes for identifying antiviral compounds effective against SARS-CoV-2.

The large number of available sequences and the diversity of protein functions challenge current experimental and computational approaches to determining and predicting protein function. We present a deep learning Graph Convolutional Network (GCN) for predicting protein functions and concurrently identifying functionally important residues. This model is initially trained using experimentally determined structures from the Protein Data Bank (PDB) but has significant de-noising capability, with only a minor drop in performance observed when structure predictions are used. We take advantage of this denoising property to train the model on > 200,000 protein structures, including many homology-predicted structures, greatly expanding the reach and applications of the method. Our model learns general structure-function relationships by robustly predicting functions of proteins with ≤ 40% sequence identity to the training set. We show that our GCN architecture predicts functions more accurately than Convolutional Neural Networks trained on sequence data alone and previous competing methods. Using class activation mapping, we automatically identify structural regions at the residue-level that lead to each function prediction for every confidently predicted protein, advancing site-specific function prediction. We use our method to annotate PDB and SWISS-MODEL proteins, making several new confident function predictions spanning both fold and function classifications.