Hybrid computational methods for combining structural data from different sources and resolutions are becoming an essential part of structural biology, especially as the field moves toward the study of large macromolecular assemblies. We have developed the molecular dynamics flexible fitting (MDFF) method for combining high-resolution atomic structures with cryo-electron microscopy (cryo-EM) maps, that results in atomic models representing the conformational state captured by cryo-EM. The method has been applied successfully to the ribosome, a ribonucleoprotein complex responsible for protein synthesis. MDFF involves a molecular dynamics simulation in which a guiding potential, based on the cryo-EM map, is added to the standard force field. Forces proportional to the gradient of the density map guide an atomic structure, available from X-ray crystallography, into high-density regions of a cryo-EM map. In this paper we describe the necessary steps to set up, run, and analyze MDFF simulations and the software packages that implement the corresponding functionalities.

The ribosomal L1 stalk is a mobile structure implicated in directing tRNA movement during translocation through the ribosome. This article investigates three aspects of L1 stalk:tRNA interaction. First, by combining through the molecular dynamics flexible fitting method data from cryo-electron microscopy, X-ray crystallography, and molecular dynamics simulations, atomic models of different tRNAs occupying the hybrid P/E state interacting with the L1 stalk are obtained. These models confirm the assignment of FRET states from previous single-molecule investigations of L1 stalk dynamics. Second, the models reconcile how initiator tRNAfMet interacts less strongly with the L1 stalk than elongator tRNAPhe, as seen in previous single-molecule experiments. Third, results from a simulation of the entire ribosome in which the L1 stalk is moved from a half-closed to its open conformation are found to support the hypothesis that L1 stalk opening is involved in tRNA release from the ribosome.

Ab initio molecular dynamics simulations were performed in order to study chemisorption, electronic properties, and desorption of glycine at wet pyrite surfaces focusing on the role of surface point defects. The change in the electronic structure and its influence on the chemical reactivity of the free FeS(2)(100) surface due to sulfur vacancies was studied in detail yielding several adsorption modes of glycine and water molecules. Energetically preferred adsorption modes were furthermore investigated in the presence of hot pressurized water mimicking "Iron Sulfur World" prebiotic conditions. The metadynamics Car-Parrinello technique was employed to map the free energy landscape including paths and barriers for desorption of glycine from such wet defective surfaces. The ubiquitous sulfur vacancies are found to increase the retention time of the adsorbed amino acid by many orders of magnitudes in comparison to the ideal pyrite-water interface. The importance of these findings in terms of a possible two-dimensional primordial chemistry on mineral surfaces is discussed.

No abstract

A joint experimental/theoretical investigation of the elastin-like octapeptide GVG(VPGVG) was carried out. In this article a comprehensive molecular-dynamics study of the temperature-dependent folding and unfolding of the octapeptide is presented. The current study, as well as its experimental counterpart (see companion article in this issue) find that this peptide undergoes an inverse temperature transition (ITT), leading to a folding at approximately 40-60 degrees C. In addition, an unfolding transition is identified at unusually high temperatures approaching the normal boiling point of water. Due to the small size of the system, two broad temperature regimes are found: the ITT regime at approximately 10-60 degrees C and the unfolding regime at approximately T > 60 degrees C, where the peptide has a maximum probability of being folded at T approximately 60 degrees C. A detailed molecular picture involving a thermodynamic order parameter, or reaction coordinate, for this process is presented along with a time-correlation function analysis of the hydrogen-bond dynamics within the peptide as well as between the peptide and solvating water molecules. Correlation with experimental evidence and ramifications on the properties of elastin are discussed.