Simulation and analysis of in vitro DNA evolution

We study theoretically the in vitro evolution of a DNA sequence by binding to a transcription factor. Using a simple model of protein-DNA binding and available binding constants for the Mnt protein, we perform large-scale, realistic simulations of evolution starting from a single DNA sequence. We identify different parameter regimes characterized by distinct evolutionary behaviors. For each regime we find analytical estimates which agree well with simulation results. For small population sizes, the DNA evolutional path is a random walk on a smooth landscape. While for large population sizes, the evolution dynamics can be well described by a mean-field theory. We also study how the details of the DNA-protein interaction affect the evolution.Comment: 11 pages, 11 figures. Submitted to PNA

Specific Recognition of p53 Tetramers by Peptides Derived from p53 Interacting Proteins

Oligomerization plays a major role in regulating the activity of many proteins, and in modulating their interactions. p53 is a homotetrameric transcription factor that has a pivotal role in tumor suppression. Its tetramerization domain is contained within its C-terminal domain, which is a site for numerous protein-protein interactions. Those can either depend on or regulate p53 oligomerization. Here we screened an array of peptides derived from proteins known to bind the tetrameric p53 C-terminal domain (p53CTD) and identified ten binding peptides. We quantitatively characterized their binding to p53CTD using fluorescence anisotropy. The peptides bound tetrameric p53CTD with micromolar affinities. Despite the high charge of the binding peptides, electrostatics contributed only mildly to the interactions. NMR studies indicated that the peptides bound p53CTD at defined sites. The most significant chemical shift deviations were observed for the peptides WS100B(81–92), which bound directly to the p53 tetramerization domain, and PKCα(281–295), which stabilized p53CTD in circular dichroism thermal denaturation studies. Using analytical ultracentrifugation, we found that several of the peptides bound preferentially to p53 tetramers. Our results indicate that the protein-protein interactions of p53 are dependent on the oligomerization state of p53. We conclude that peptides may be used to regulate the oligomerization of p53

A Coarse-Grained Biophysical Model of E. coli and Its Application to Perturbation of the rRNA Operon Copy Number

We propose a biophysical model of Escherichia coli that predicts growth rate and an effective cellular composition from an effective, coarse-grained representation of its genome. We assume that E. coli is in a state of balanced exponential steadystate growth, growing in a temporally and spatially constant environment, rich in resources. We apply this model to a series of past measurements, where the growth rate and rRNA-to-protein ratio have been measured for seven E. coli strains with an rRNA operon copy number ranging from one to seven (the wild-type copy number). These experiments show that growth rate markedly decreases for strains with fewer than six copies. Using the model, we were able to reproduce these measurements. We show that the model that best fits these data suggests that the volume fraction of macromolecules inside E. coli is not fixed when the rRNA operon copy number is varied. Moreover, the model predicts that increasing the copy number beyond seven results in a cytoplasm densely packed with ribosomes and proteins. Assuming that under such overcrowded conditions prolonged diffusion times tend to weaken binding affinities, the model predicts that growth rate will not increase substantially beyond the wild-type growth rate, as indicated by other experiments. Our model therefore suggests that changing the rRNA operon copy number of wild-type E. coli cells growing in a constant rich environment does not substantially increase their growth rate. Other observations regarding strains with an altered rRNA operon copy number, such as nucleoid compaction and the rRNA operon feedback response, appear to be qualitatively consistent with this model. In addition, we discuss possible design principles suggested by the model and propose further experiments to test its validity

Extreme disorder in an ultrahigh-affinity protein complex

Molecular communication in biology is mediated by protein interactions. According to the current paradigm, the specificity and affinity required for these interactions are encoded in the precise complementarity of binding interfaces. Even proteins that are disordered under physiological conditions or that contain large unstructured regions commonly interact with well-structured binding sites on other biomolecules. Here we demonstrate the existence of an unexpected interaction mechanism: the two intrinsically disordered human proteins histone H1 and its nuclear chaperone prothymosin-α associate in a complex with picomolar affinity, but fully retain their structural disorder, long-range flexibility and highly dynamic character. On the basis of closely integrated experiments and molecular simulations, we show that the interaction can be explained by the large opposite net charge of the two proteins, without requiring defined binding sites or interactions between specific individual residues. Proteome-wide sequence analysis suggests that this interaction mechanism may be abundant in eukaryotes

Immunoglobulin, glucocorticoid, or combination therapy for multisystem inflammatory syndrome in children: a propensity-weighted cohort study

Background Multisystem inflammatory syndrome in children (MIS-C), a hyperinflammatory condition associated with SARS-CoV-2 infection, has emerged as a serious illness in children worldwide. Immunoglobulin or glucocorticoids, or both, are currently recommended treatments. Methods The Best Available Treatment Study evaluated immunomodulatory treatments for MIS-C in an international observational cohort. Analysis of the first 614 patients was previously reported. In this propensity-weighted cohort study, clinical and outcome data from children with suspected or proven MIS-C were collected onto a web-based Research Electronic Data Capture database. After excluding neonates and incomplete or duplicate records, inverse probability weighting was used to compare primary treatments with intravenous immunoglobulin, intravenous immunoglobulin plus glucocorticoids, or glucocorticoids alone, using intravenous immunoglobulin as the reference treatment. Primary outcomes were a composite of inotropic or ventilator support from the second day after treatment initiation, or death, and time to improvement on an ordinal clinical severity scale. Secondary outcomes included treatment escalation, clinical deterioration, fever, and coronary artery aneurysm occurrence and resolution. This study is registered with the ISRCTN registry, ISRCTN69546370. Findings We enrolled 2101 children (aged 0 months to 19 years) with clinically diagnosed MIS-C from 39 countries between June 14, 2020, and April 25, 2022, and, following exclusions, 2009 patients were included for analysis (median age 8·0 years [IQR 4·2–11·4], 1191 [59·3%] male and 818 [40·7%] female, and 825 [41·1%] White). 680 (33·8%) patients received primary treatment with intravenous immunoglobulin, 698 (34·7%) with intravenous immunoglobulin plus glucocorticoids, 487 (24·2%) with glucocorticoids alone; 59 (2·9%) patients received other combinations, including biologicals, and 85 (4·2%) patients received no immunomodulators. There were no significant differences between treatments for primary outcomes for the 1586 patients with complete baseline and outcome data that were considered for primary analysis. Adjusted odds ratios for ventilation, inotropic support, or death were 1·09 (95% CI 0·75–1·58; corrected p value=1·00) for intravenous immunoglobulin plus glucocorticoids and 0·93 (0·58–1·47; corrected p value=1·00) for glucocorticoids alone, versus intravenous immunoglobulin alone. Adjusted average hazard ratios for time to improvement were 1·04 (95% CI 0·91–1·20; corrected p value=1·00) for intravenous immunoglobulin plus glucocorticoids, and 0·84 (0·70–1·00; corrected p value=0·22) for glucocorticoids alone, versus intravenous immunoglobulin alone. Treatment escalation was less frequent for intravenous immunoglobulin plus glucocorticoids (OR 0·15 [95% CI 0·11–0·20]; p<0·0001) and glucocorticoids alone (0·68 [0·50–0·93]; p=0·014) versus intravenous immunoglobulin alone. Persistent fever (from day 2 onward) was less common with intravenous immunoglobulin plus glucocorticoids compared with either intravenous immunoglobulin alone (OR 0·50 [95% CI 0·38–0·67]; p<0·0001) or glucocorticoids alone (0·63 [0·45–0·88]; p=0·0058). Coronary artery aneurysm occurrence and resolution did not differ significantly between treatment groups. Interpretation Recovery rates, including occurrence and resolution of coronary artery aneurysms, were similar for primary treatment with intravenous immunoglobulin when compared to glucocorticoids or intravenous immunoglobulin plus glucocorticoids. Initial treatment with glucocorticoids appears to be a safe alternative to immunoglobulin or combined therapy, and might be advantageous in view of the cost and limited availability of intravenous immunoglobulin in many countries. Funding Imperial College London, the European Union's Horizon 2020, Wellcome Trust, the Medical Research Foundation, UK National Institute for Health and Care Research, and National Institutes of Health