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.

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.

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.IMPORTANCE Bacteria are more structurally complex than is commonly appreciated. Here we present a survey of previously uncharacterized structures that we observed in bacterial cells by electron cryotomography, structures that will initiate new lines of research investigating their identities and roles.

Bdellovibrio predation cycle characterized at nanometre-scale resolution with cryo-electron tomography.

Bdellovibrio bacteriovorus is a microbial predator that offers promise as a living antibiotic for its ability to kill Gram-negative bacteria, including human pathogens. Even after six decades of study, fundamental details of its predation cycle remain mysterious. Here we used cryo-electron tomography to comprehensively image the lifecycle of B. bacteriovorus at nanometre-scale resolution. With high-resolution images of predation in a native (hydrated, unstained) state, we discover several surprising features of the process, including macromolecular complexes involved in prey attachment/invasion and a flexible portal structure lining a hole in the prey peptidoglycan that tightly seals the prey outer membrane around the predator during entry. Unexpectedly, we find that B. bacteriovorus does not shed its flagellum during invasion, but rather resorbs it into its periplasm for degradation. Finally, following growth and division in the bdelloplast, we observe a transient and extensive ribosomal lattice on the condensed B. bacteriovorus nucleoid.

Integrative structural modelling and visualisation of a cellular organelle.

Models of insulin secretory vesicles from pancreatic beta cells have been created using the cellPACK suite of tools to research, curate, construct and visualise the current state of knowledge. The model integrates experimental information from proteomics, structural biology, cryoelectron microscopy and X-ray tomography, and is used to generate models of mature and immature vesicles. A new method was developed to generate a confidence score that reconciles inconsistencies between three available proteomes using expert annotations of cellular localisation. The models are used to simulate soft X-ray tomograms, allowing quantification of features that are observed in experimental tomograms, and in turn, allowing interpretation of X-ray tomograms at the molecular level.

Assisted peptide folding by surface pattern recognition.

Natively disordered proteins belong to a unique class of biomolecules whose function is related to their flexibility and their ability to adopt desired conformations upon binding to substrates. In some cases these proteins can bind multiple partners, which can lead to distinct structures and promiscuity in functions. In other words, the capacity to recognize molecular patterns on the substrate is often essential for the folding and function of intrinsically disordered proteins. Biomolecular pattern recognition is extremely relevant both in vivo (e.g., for oligomerization, immune response, induced folding, substrate binding, and molecular switches) and in vitro (e.g., for biosensing, catalysis, chromatography, and implantation). Here, we use a minimalist computational model system to investigate how polar/nonpolar patterns on a surface can induce the folding of an otherwise unstructured peptide. We show that a model peptide that exists in the bulk as a molten globular state consisting of many interconverting structures can fold into either a helix-coil-helix or an extended helix structure in the presence of a complementary designed patterned surface at low hydrophobicity (3.7%) or a uniform surface at high hydrophobicity (50%). However, we find that a carefully chosen surface pattern can bind to and catalyze the folding of a natively unfolded protein much more readily or effectively than a surface with a noncomplementary or uniform distribution of hydrophobic residues.

Reconciling theories of chaperonin accelerated folding with experimental evidence.

For the last 20 years, a large volume of experimental and theoretical work has been undertaken to understand how chaperones like GroEL can assist protein folding in the cell. The most accepted explanation appears to be the simplest: GroEL, like most other chaperones, helps proteins fold by preventing aggregation. However, evidence suggests that, under some conditions, GroEL can play a more active role by accelerating protein folding. A large number of models have been proposed to explain how this could occur. Focused experiments have been designed and carried out using different protein substrates with conclusions that support many different mechanisms. In the current article, we attempt to see the forest through the trees. We review all suggested mechanisms for chaperonin-mediated folding and weigh the plausibility of each in light of what we now know about the most stringent, essential, GroEL-dependent protein substrates.

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.

Structural diversity of dimers of the Alzheimer amyloid-beta(25-35) peptide and polymorphism of the resulting fibrils.

The 25-35 fragment of the Alzheimer amyloid beta (Abeta) peptide is a naturally occurring proteolytic by-product that retains the toxicity of its larger, better-known counterpart, Abeta (1-40). Soluble oligomers of the amyloid-beta peptide have been implicated in the pathogenesis of Alzheimer's disease as a primary source of neurotoxicity. These oligomers are difficult to characterize experimentally due to their transient nature. As a result, a detailed knowledge of oligomeric structures at the atomic level is lacking. Using replica exchange molecular dynamics simulations, we investigated the conformations adopted by dimers, the smallest soluble oligomers of Abeta(25-35). Our simulations, which total 4 mus in length, reveal a diverse ensemble of well-organized dimers with high beta-sheet content coexisting with unstructured dimer complexes. The structured dimers comprise parallel and antiparallel extended beta-strand, beta-hairpin, and V-shaped beta-strand conformations. Protofibril models constructed from the extended and V-shaped dimers lead to stable structures consistent with experimentally available data from H/D exchange NMR and AFM spectroscopy. Our simulations suggest that fibril polymorphism may be encoded in the early stages of aggregation for the Abeta(25-35) peptide.

The effect of surface tethering on the folding of the src-SH3 protein domain.

The effect of surface tethering on the folding mechanism of the src-SH3 protein domain was investigated using a coarse-grained Go-type protein model. The protein was tethered at various locations along the protein chain and the thermodynamics and kinetics of folding were studied using replica exchange and constant temperature Langevin dynamics. Our simulations reveal that tethering in a structured part of the transition state can dramatically alter the folding mechanism, while tethering in an unstructured part leaves the folding mechanism unaltered as compared to bulk folding. Interestingly, there is only modest correlation between the tethering effect on the folding mechanism and its effect on thermodynamic stability and folding rates. We suggest locations on the protein at which tethering could be performed in single-molecule experiments so as to leave the folding mechanism unaltered from the bulk.

Cooperativity, smooth energy landscapes and the origins of topology-dependent protein folding rates.

The relative folding rates of simple, single-domain proteins, proteins whose folding energy landscapes are smooth, are highly dispersed and strongly correlated with native-state topology. In contrast, the relative folding rates of small, Go-potential lattice polymers, which also exhibit smooth energy landscapes, are poorly dispersed and insignificantly correlated with native-state topology. Here, we investigate this discrepancy in light of a recent, quantitative theory of two-state folding kinetics, the topomer search model. This model stipulates that the topology-dependence of two-state folding rates is a direct consequence of the extraordinarily cooperative equilibrium folding of simple proteins. We demonstrate that traditional Go polymers lack the extreme cooperativity that characterizes the folding of naturally occurring, two-state proteins and confirm that the folding rates of a diverse set of Go 27-mers are poorly dispersed and effectively uncorrelated with native state topology. Upon modestly increasing the cooperativity of the Go-potential, however, significantly increased dispersion and strongly topology-dependent kinetics are observed. These results support previous arguments that the cooperative folding of simple, single-domain proteins gives rise to their topology-dependent folding rates. We speculate that this cooperativity, and thus, indirectly, the topology-rate relationship, may have arisen in order to generate the smooth energetic landscapes upon which rapid folding can occur.

Relative stability of de novo four-helix bundle proteins: insights from coarse grained molecular simulations.

We use a recently developed coarse-grained computational model to investigate the relative stability of two different sets of de novo designed four-helix bundle proteins. Our simulations suggest a possible explanation for the experimentally observed increase in stability of the four-helix bundles with increasing sequence length. In details, we show that both short subsequences composed only by polar residues and additional nonpolar residues inserted, via different point mutations in ad hoc positions, seem to play a significant role in stabilizing the four-helix bundle conformation in the longer sequences. Finally, we propose an additional mutation that rescues a short amino acid sequence that would otherwise adopt a compact misfolded state. Our work suggests that simple computational models can be used as a complementary tool in the design process of de novo proteins.

MINRMS: an efficient algorithm for determining protein structure similarity using root-mean-squared-distance.

MOTIVATION: Existing algorithms for automated protein structure alignment generate contradictory results and are difficult to interpret. An algorithm which can provide a context for interpreting the alignment and uses a simple method to characterize protein structure similarity is needed. RESULTS: We describe a heuristic for limiting the search space for structure alignment comparisons between two proteins, and an algorithm for finding minimal root-mean-squared-distance (RMSD) alignments as a function of the number of matching residue pairs within this limited search space. Our alignment algorithm uses coordinates of alpha-carbon atoms to represent each amino acid residue and requires a total computation time of O(m(3) n(2)), where m and n denote the lengths of the protein sequences. This makes our method fast enough for comparisons of moderate-size proteins (fewer than approximately 800 residues) on current workstation-class computers and therefore addresses the need for a systematic analysis of multiple plausible shape similarities between two proteins using a widely accepted comparison metric.