Moltemplate: A Tool for Coarse-Grained Modeling of Complex Biological Matter and Soft Condensed Matter Physics

Coarse-grained models have long been considered indispensable tools in the investigation of biomolecular dynamics and assembly. However, the process of simulating such models is arduous because unconventional force fields and particle attributes are often needed, and some systems are not in thermal equilibrium. Although modern molecular dynamics programs are highly adaptable, software designed for preparing all-atom simulations typically makes restrictive assumptions about the nature of the particles and the forces acting on them. Consequently, the use of coarse-grained models has remained challenging. Moltemplate is a file format for storing coarse-grained molecular models and the forces that act on them, as well as a program that converts moltemplate files into input files for LAMMPS, a popular molecular dynamics engine. Moltemplate has broad scope and an emphasis on generality. It accommodates new kinds of forces as they are developed for LAMMPS, making moltemplate a popular tool with thousands of users in computational chemistry, materials science, and structural biology. To demonstrate its wide functionality, we provide examples of using moltemplate to prepare simulations of fluids using many-body forces, coarse-grained organic semiconductors, and the motor-driven supercoiling and condensation of an entire bacterial chromosome

Uncharacterized bacterial structures revealed by electron cryotomography

Electron cryotomography (ECT) can reveal the native structure and arrangement of macromolecular complexes inside intact cells. This technique has greatly advanced our understanding of the ultrastructure of bacterial cells. We now view bacteria as structurally complex assemblies of macromolecular machines rather than as undifferentiated bags of enzymes. To date, our group has applied ECT to nearly 90 different bacterial species, collecting more than 15,000 cryotomograms. In addition to known structures, we have observed, to our knowledge, several uncharacterized features in these tomograms. Some are completely novel structures; others expand the features or species range of known structure types. Here, we present a survey of these uncharacterized bacterial structures in the hopes of accelerating their identification and study, and furthering our understanding of the structural complexity of bacterial cells

Chromatin extrusion explains key features of loop and domain formation in wild-type and engineered genomes

We recently used in situ Hi-C to create kilobase-resolution 3D maps of mammalian genomes. Here, we combine these maps with new Hi-C, microscopy, and genome-editing experiments to study the physical structure of chromatin fibers, domains, and loops. We find that the observed contact domains are inconsistent with the equilibrium state for an ordinary condensed polymer. Combining Hi-C data and novel mathematical theorems, we show that contact domains are also not consistent with a fractal globule. Instead, we use physical simulations to study two models of genome folding. In one, intermonomer attraction during polymer condensation leads to formation of an anisotropic “tension globule.” In the other, CCCTC-binding factor (CTCF) and cohesin act together to extrude unknotted loops during interphase. Both models are consistent with the observed contact domains and with the observation that contact domains tend to form inside loops. However, the extrusion model explains a far wider array of observations, such as why loops tend not to overlap and why the CTCF-binding motifs at pairs of loop anchors lie in the convergent orientation. Finally, we perform 13 genome-editing experiments examining the effect of altering CTCF-binding sites on chromatin folding. The convergent rule correctly predicts the affected loops in every case. Moreover, the extrusion model accurately predicts in silico the 3D maps resulting from each experiment using only the location of CTCF-binding sites in the WT. Thus, we show that it is possible to disrupt, restore, and move loops and domains using targeted mutations as small as a single base pair.National Science Foundation (U.S.) (Grant PHY-1427654)National Institutes of Health (U.S.) (New Innovator Award 1DP2OD008540-01)Cancer Prevention and Research Institute of Texas (Scholar Award R1304)Baylor College of Medicine (McNair Medical Institute Scholar Award)Presidential Early Career Award for Scientists and Engineer

De novo design of protein logic gates

The design of modular protein logic for regulating protein function at the posttranscriptional level is a challenge for synthetic biology. Here, we describe the design of two-input AND, OR, NAND, NOR, XNOR, and NOT gates built from de novo–designed proteins. These gates regulate the association of arbitrary protein units ranging from split enzymes to transcriptional machinery in vitro, in yeast and in primary human T cells, where they control the expression of the TIM3 gene related to T cell exhaustion. Designed binding interaction cooperativity, confirmed by native mass spectrometry, makes the gates largely insensitive to stoichiometric imbalances in the inputs, and the modularity of the approach enables ready extension to three-input OR, AND, and disjunctive normal form gates. The modularity and cooperativity of the control elements, coupled with the ability to de novo design an essentially unlimited number of protein components, should enable the design of sophisticated posttranslational control logic over a wide range of biological functions

The Circadian Response of Intrinsically Photosensitive Retinal Ganglion Cells

Intrinsically photosensitive retinal ganglion cells (ipRGC) signal environmental light level to the central circadian clock and contribute to the pupil light reflex. It is unknown if ipRGC activity is subject to extrinsic (central) or intrinsic (retinal) network-mediated circadian modulation during light entrainment and phase shifting. Eleven younger persons (18–30 years) with no ophthalmological, medical or sleep disorders participated. The activity of the inner (ipRGC) and outer retina (cone photoreceptors) was assessed hourly using the pupil light reflex during a 24 h period of constant environmental illumination (10 lux). Exogenous circadian cues of activity, sleep, posture, caffeine, ambient temperature, caloric intake and ambient illumination were controlled. Dim-light melatonin onset (DLMO) was determined from salivary melatonin assay at hourly intervals, and participant melatonin onset values were set to 14 h to adjust clock time to circadian time. Here we demonstrate in humans that the ipRGC controlled post-illumination pupil response has a circadian rhythm independent of external light cues. This circadian variation precedes melatonin onset and the minimum ipRGC driven pupil response occurs post melatonin onset. Outer retinal photoreceptor contributions to the inner retinal ipRGC driven post-illumination pupil response also show circadian variation whereas direct outer retinal cone inputs to the pupil light reflex do not, indicating that intrinsically photosensitive (melanopsin) retinal ganglion cells mediate this circadian variation

A MISSING-LINK IN THE SUPERNOVA-GRB CONNECTION: THE CASE OF SN 2012ap

Gamma Ray Bursts (GRBs) are characterized by ultra-relativistic outflows, while supernovae are generally characterized by non-relativistic ejecta. GRB afterglows decelerate rapidly usually within days, because their low-mass ejecta rapidly sweep up a comparatively larger mass of circumstellar material. However supernovae, with heavy ejecta, can be in nearly free expansion for centuries. Supernovae were thought to have non-relativistic outflows except for few relativistic ones accompanied by GRBs. This clear division was blurred by SN 2009bb, the first supernova with a relativistic outflow without an observed GRB. Yet the ejecta from SN 2009bb was baryon loaded, and in nearly-free expansion for a year, unlike GRBs. We report the first supernova discovered without a GRB, but with rapidly decelerating mildly relativistic ejecta, SN 2012ap. We discovered a bright and rapidly evolving radio counterpart driven by the circumstellar interaction of the relativistic ejecta. However, we did not find any coincident GRB with an isotropic fluence of more than a sixth of the fluence from GRB 980425. This shows for the first time that central engines in type Ic supernovae, even without an observed GRB, can produce both relativistic and rapidly decelerating outflows like GRBs.Comment: 8 pages, 5 figures, 1 table, accepted for publication in Ap

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead

Short FtsZ filaments can drive asymmetric cell envelope constriction at the onset of bacterial cytokinesis

FtsZ, the bacterial homologue of eukaryotic tubulin, plays a central role in cell division in nearly all bacteria and many archaea. It forms filaments under the cytoplasmic membrane at the division site where, together with other proteins it recruits, it drives peptidoglycan synthesis and constricts the cell. Despite extensive study, the arrangement of FtsZ filaments and their role in division continue to be debated. Here, we apply electron cryotomography to image the native structure of intact dividing cells and show that constriction in a variety of Gram‐negative bacterial cells, including Proteus mirabilis and Caulobacter crescentus, initiates asymmetrically, accompanied by asymmetric peptidoglycan incorporation and short FtsZ‐like filament formation. These results show that a complete ring of FtsZ is not required for constriction and lead us to propose a model for FtsZ‐driven division in which short dynamic FtsZ filaments can drive initial peptidoglycan synthesis and envelope constriction at the onset of cytokinesis, later increasing in length and number to encircle the division plane and complete constriction

Integrated Tools

this paper is on new tools for structural and sequence analysis and visualization. AlignPlot, written in C and Python, provides a graphical representation of the RMSD values for each alignment in the set, allowing the user to quickly identify the regions of two structures that are most similar. Particularly important, it provides a user-friendly way to display specific alignments on the screen and navigate among them. MSFviewer, written in Python, provides an integrated link to sequence space, displaying multiple alignments of related sequences on the screen and providing for interactive highlighting of a selected structural align- ment and the associated multiple sequence alignmen