Molecular Dynamics Simulations Suggest a Non-Doublet Decoding Model of -1 Frameshifting by tRNASer3.

In-frame decoding in the ribosome occurs through canonical or wobble Watson-Crick pairing of three mRNA codon bases (a triplet) with a triplet of anticodon bases in tRNA. Departures from the triplet-triplet interaction can result in frameshifting, meaning downstream mRNA codons are then read in a different register. There are many mechanisms to induce frameshifting, and most are insufficiently understood. One previously proposed mechanism is doublet decoding, in which only codon bases 1 and 2 are read by anticodon bases 34 and 35, which would lead to -1 frameshifting. In E. coli, tRNASer3GCU can induce -1 frameshifting at alanine (GCA) codons. The logic of the doublet decoding model is that the Ala codon's GC could pair with the tRNASer3's GC, leaving the third anticodon residue U36 making no interactions with mRNA. Under that model, a U36C mutation would still induce -1 frameshifting, but experiments refute this. We perform all-atom simulations of wild-type tRNASer3, as well as a U36C mutant. Our simulations revealed a hydrogen bond between U36 of the anticodon and G1 of the codon. The U36C mutant cannot make this interaction, as it lacks the hydrogen-bond-donating H3. The simulation thus suggests a novel, non-doublet decoding mechanism for -1 frameshifting by tRNASer3 at Ala codons.

MMB-GUI: a fast morphing method demonstrates a possible ribosomal tRNA translocation trajectory.

Easy-to-use macromolecular viewers, such as UCSF Chimera, are a standard tool in structural biology. They allow rendering and performing geometric operations on large complexes, such as viruses and ribosomes. Dynamical simulation codes enable modeling of conformational changes, but may require considerable time and many CPUs. There is an unmet demand from structural and molecular biologists for software in the middle ground, which would allow visualization combined with quick and interactive modeling of conformational changes, even of large complexes. This motivates MMB-GUI. MMB uses an internal-coordinate, multiscale approach, yielding as much as a 2000-fold speedup over conventional simulation methods. We use Chimera as an interactive graphical interface to control MMB. We show how this can be used for morphing of macromolecules that can be heterogeneous in biopolymer type, sequence, and chain count, accurately recapitulating structural intermediates. We use MMB-GUI to create a possible trajectory of EF-G mediated gate-passing translocation in the ribosome, with all-atom structures. This shows that the GUI makes modeling of large macromolecules accessible to a wide audience. The morph highlights similarities in tRNA conformational changes as tRNA translocates from A to P and from P to E sites and suggests that tRNA flexibility is critical for translocation completion.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems.

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Formation of raft-like assemblies within clusters of influenza hemagglutinin observed by MD simulations.

The association of hemagglutinin (HA) with lipid rafts in the plasma membrane is an important feature of the assembly process of influenza virus A. Lipid rafts are thought to be small, fluctuating patches of membrane enriched in saturated phospholipids, sphingolipids, cholesterol and certain types of protein. However, raft-associating transmembrane (TM) proteins generally partition into Ld domains in model membranes, which are enriched in unsaturated lipids and depleted in saturated lipids and cholesterol. The reason for this apparent disparity in behavior is unclear, but model membranes differ from the plasma membrane in a number of ways. In particular, the higher protein concentration in the plasma membrane may influence the partitioning of membrane proteins for rafts. To investigate the effect of high local protein concentration, we have conducted coarse-grained molecular dynamics (CG MD) simulations of HA clusters in domain-forming bilayers. During the simulations, we observed a continuous increase in the proportion of raft-type lipids (saturated phospholipids and cholesterol) within the area of membrane spanned by the protein cluster. Lateral diffusion of unsaturated lipids was significantly attenuated within the cluster, while saturated lipids were relatively unaffected. On this basis, we suggest a possible explanation for the change in lipid distribution, namely that steric crowding by the slow-diffusing proteins increases the chemical potential for unsaturated lipids within the cluster region. We therefore suggest that a local aggregation of HA can be sufficient to drive association of the protein with raft-type lipids. This may also represent a general mechanism for the targeting of TM proteins to rafts in the plasma membrane, which is of functional importance in a wide range of cellular processes.

Game on, science - how video game technology may help biologists tackle visualization challenges.

The video games industry develops ever more advanced technologies to improve rendering, image quality, ergonomics and user experience of their creations providing very simple to use tools to design new games. In the molecular sciences, only a small number of experts with specialized know-how are able to design interactive visualization applications, typically static computer programs that cannot easily be modified. Are there lessons to be learned from video games? Could their technology help us explore new molecular graphics ideas and render graphics developments accessible to non-specialists? This approach points to an extension of open computer programs, not only providing access to the source code, but also delivering an easily modifiable and extensible scientific research tool. In this work, we will explore these questions using the Unity3D game engine to develop and prototype a biological network and molecular visualization application for subsequent use in research or education. We have compared several routines to represent spheres and links between them, using either built-in Unity3D features or our own implementation. These developments resulted in a stand-alone viewer capable of displaying molecular structures, surfaces, animated electrostatic field lines and biological networks with powerful, artistic and illustrative rendering methods. We consider this work as a proof of principle demonstrating that the functionalities of classical viewers and more advanced novel features could be implemented in substantially less time and with less development effort. Our prototype is easily modifiable and extensible and may serve others as starting point and platform for their developments. A webserver example, standalone versions for MacOS X, Linux and Windows, source code, screen shots, videos and documentation are available at the address: http://unitymol.sourceforge.net/.

Mixing Atomistic and Coarse Grain Solvation Models for MD Simulations: Let WT4 Handle the Bulk.

Accurate simulation of biomolecular systems requires the consideration of solvation effects. The arrangement and dynamics of water close to a solute are strongly influenced by the solute itself. However, as the solute-solvent distance increases, the water properties tend to those of the bulk liquid. This suggests that bulk regions can be treated at a coarse grained (CG) level, while keeping the atomistic details around the solute. Since water represents about 80% of any biological system, this approach may offer a significant reduction in the computational cost of simulations without compromising atomistic details. We show here that mixing the popular SPC water model with a CG model for solvation (called WatFour) can effectively mimic the hydration, structure, and dynamics of molecular systems composed of pure water, simple electrolyte solutions, and solvated macromolecules. As a nontrivial example, we present simulations of the SNARE membrane fusion complex, a trimeric protein-protein complex embedded in a double phospholipid bilayer. Comparison with a fully atomistic reference simulation illustrates the equivalence between both approaches.

GPU-accelerated atom and dynamic bond visualization using hyperballs: a unified algorithm for balls, sticks, and hyperboloids.

Ray casting on graphics processing units (GPUs) opens new possibilities for molecular visualization. We describe the implementation and calculation of diverse molecular representations such as licorice, ball-and-stick, space-filling van der Waals spheres, and approximated solvent-accessible surfaces using GPUs. We introduce HyperBalls, an improved ball-and-stick representation replacing tubes, linking the atom spheres by hyperboloids that can smoothly connect them. This type of depiction is particularly useful to represent dynamic phenomena, such as the evolution of noncovalent bonds. It is furthermore well suited to represent coarse-grained models and spring networks. All these representations can be defined by a single general algebraic equation that is adapted for the ray-casting technique and is well suited for execution on the GPU. Using GPU capabilities, this implementation can routinely, accurately, and interactively render molecules ranging from a few atoms up to huge macromolecular assemblies with more than 500,000 particles. In simple cases, based only on spheres, we have been able to display up to two million atoms smoothly.