Conformational ensembles of an RNA hairpin using molecular dynamics and sparse NMR data.

Solution nuclear magnetic resonance (NMR) experiments allow RNA dynamics to be determined in an aqueous environment. However, when a limited number of peaks are assigned, it is difficult to obtain structural information. We here show a protocol based on the combination of experimental data (Nuclear Overhauser Effect, NOE) and molecular dynamics simulations with enhanced sampling methods. This protocol allows to (a) obtain a maximum entropy ensemble compatible with NMR restraints and (b) obtain a minimal set of metastable conformations compatible with the experimental data (maximum parsimony). The method is applied to a hairpin of 29 nt from an inverted SINEB2, which is part of the SINEUP family and has been shown to enhance protein translation. A clustering procedure is introduced where the annotation of base-base interactions and glycosidic bond angles is used as a metric. By reweighting the contributions of the clusters, minimal sets of four conformations could be found which are compatible with the experimental data. A motif search on the structural database showed that some identified low-population states are present in experimental structures of other RNA transcripts. The introduced method can be applied to characterize RNA dynamics in systems where a limited amount of NMR information is available.

Barnaba: software for analysis of nucleic acid structures and trajectories.

RNA molecules are highly dynamic systems characterized by a complex interplay between sequence, structure, dynamics, and function. Molecular simulations can potentially provide powerful insights into the nature of these relationships. The analysis of structures and molecular trajectories of nucleic acids can be nontrivial because it requires processing very high-dimensional data that are not easy to visualize and interpret. Here we introduce Barnaba, a Python library aimed at facilitating the analysis of nucleic acid structures and molecular simulations. The software consists of a variety of analysis tools that allow the user to (i) calculate distances between three-dimensional structures using different metrics, (ii) back-calculate experimental data from three-dimensional structures, (iii) perform cluster analysis and dimensionality reductions, (iv) search three-dimensional motifs in PDB structures and trajectories, and (v) construct elastic network models for nucleic acids and nucleic acids-protein complexes. In addition, Barnaba makes it possible to calculate torsion angles, pucker conformations, and to detect base-pairing/base-stacking interactions. Barnaba produces graphics that conveniently visualize both extended secondary structure and dynamics for a set of molecular conformations. The software is available as a command-line tool as well as a library, and supports a variety of file formats such as PDB, dcd, and xtc files. Source code, documentation, and examples are freely available at https://github.com/srnas/barnaba under GNU GPLv3 license.

2H-NMR and MD Simulations Reveal Membrane-Bound Conformation of Magainin 2 and Its Synergy with PGLa.

Magainin 2 (MAG2) and PGLa are two α-helical antimicrobial peptides found in the skin of the African frog Xenopus laevis. They act by permeabilizing bacterial membranes and exhibit an exemplary synergism. Here, we determined the detailed molecular alignment and dynamical behavior of MAG2 in oriented lipid bilayers by using 2H-NMR on Ala-d3-labeled peptides, which yielded orientation-dependent quadrupolar splittings of the labels. The amphiphilic MAG2 helix was found to lie flat on the membrane surface in 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC)/1,2-dimyristoyl-sn-glycero-3-phosphatidylglycerol (DMPG) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC)/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG), as expected, with a tilt angle close to 90°. This orientation fits well with all-atom molecular-dynamics simulations of MAG2 performed in DMPC and DMPC/DMPG. In the presence of an equimolar amount of PGLa, the NMR analysis showed that MAG2 becames tilted at an angle of 120°, and its azimuthal rotation angle also changes. Since this interaction was found to occur in a concentration range where the peptides per se do not interact with their own type, we propose that MAG2 forms a stable heterodimer with PGLa. Given that the PGLa molecules in the complex are known to be flipped into a fully upright orientation, with a helix tilt close to 180°, they must make up the actual transmembrane pore. We thus suggest that the two negative charges on the C-terminus of the obliquely tilted MAG2 peptides neutralize some of the cationic groups on the upright PGLa helices. This would stabilize the assembly of PGLa into a toroidal pore with an overall reduced charge density, which could explain the mechanism of synergy.

Structure of the membrane anchor of pestivirus glycoprotein E(rns), a long tilted amphipathic helix.

E(rns) is an essential virion glycoprotein with RNase activity that suppresses host cellular innate immune responses upon being partially secreted from the infected cells. Its unusual C-terminus plays multiple roles, as the amphiphilic helix acts as a membrane anchor, as a signal peptidase cleavage site, and as a retention/secretion signal. We analyzed the structure and membrane binding properties of this sequence to gain a better understanding of the underlying mechanisms. CD spectroscopy in different setups, as well as Monte Carlo and molecular dynamics simulations confirmed the helical folding and showed that the helix is accommodated in the amphiphilic region of the lipid bilayer with a slight tilt rather than lying parallel to the surface. This model was confirmed by NMR analyses that also identified a central stretch of 15 residues within the helix that is fully shielded from the aqueous layer, which is C-terminally followed by a putative hairpin structure. These findings explain the strong membrane binding of the protein and provide clues to establishing the E(rns) membrane contact, processing and secretion.

3D hydrophobic moment vectors as a tool to characterize the surface polarity of amphiphilic peptides.

The interaction of membranes with peptides and proteins is largely determined by their amphiphilic character. Hydrophobic moments of helical segments are commonly derived from their two-dimensional helical wheel projections, and the same is true for ß-sheets. However, to the best of our knowledge, there exists no method to describe structures in three dimensions or molecules with irregular shape. Here, we define the hydrophobic moment of a molecule as a vector in three dimensions by evaluating the surface distribution of all hydrophilic and lipophilic regions over any given shape. The electrostatic potential on the molecular surface is calculated based on the atomic point charges. The resulting hydrophobic moment vector is specific for the instantaneous conformation, and it takes into account all structural characteristics of the molecule, e.g., partial unfolding, bending, and side-chain torsion angles. Extended all-atom molecular dynamics simulations are then used to calculate the equilibrium hydrophobic moments for two antimicrobial peptides, gramicidin S and PGLa, under different conditions. We show that their effective hydrophobic moment vectors reflect the distribution of polar and nonpolar patches on the molecular surface and the calculated electrostatic surface potential. A comparison of simulations in solution and in lipid membranes shows how the peptides undergo internal conformational rearrangement upon binding to the bilayer surface. A good correlation with solid-state NMR data indicates that the hydrophobic moment vector can be used to predict the membrane binding geometry of peptides. This method is available as a web application on http://www.ibg.kit.edu/HM/.

Controlling biological activity with light: diarylethene-containing cyclic peptidomimetics.

Photobiological processes in nature are usually triggered by nonpeptidic chromophores or by modified side chains. A system is presented in which the polypeptide backbone itself can be conformationally switched by light. An amino acid analogue was designed and synthesized based on a reversibly photoisomerizable diarylethene scaffold. This analogue was incorporated into the cyclic backbone of the antimicrobial peptide gramicidin S at several sites. The biological activity of the resulting peptidomimetics could then be effectively controlled by ultraviolet/visible light within strictly defined spatial and temporal limits.

Protein ORIGAMI: A program for the creation of 3D paper models of folded peptides.

Protein ORIGAMI (http://ibg.kit.edu/protein_origami) is a browser-based web application that allows the user to create straightforward 3D paper models of folded peptides for research, teaching and presentations. An amino acid sequence can be turned into α-helices, ß-strands and random coils that can be printed out and folded into properly scaled models, with a color code denoting the biophysical characteristics of each amino acid residue (hydrophobicity, charge, etc.). These models provide an intuitive visual and tactile understanding of peptide interactions with other partners, such as helix-helix assembly, oligomerization, membrane binding, or pore formation. Helices can also be displayed as a helical wheel or helical mesh in 2D graphics, to be used in publications or presentations. The highly versatile programme Protein ORIGAMI is also suited to create less conventional helices with arbitrary pitch (e.g., 310 -helix, π-helix, or left-handed helices). Noncanonical amino acids, labels and different terminal modifications can be defined and displayed at will, and different protonation states can be shown. In addition to the web application, the program source code can be downloaded and installed locally on a PC. The printed paper models can be readily used for daily research and discussions, just as for educational purposes and teaching. © 2018 by The International Union of Biochemistry and Molecular Biology, 46:403-409, 2018.

Best of Two Worlds? How MD Simulations of Amphiphilic Helical Peptides in Membranes Can Complement Data from Oriented Solid-State NMR.

The membrane alignment of helical amphiphilic peptides in oriented phospholipid bilayers can be obtained as ensemble and time averages from solid state 2H NMR by fitting the quadrupolar splittings to ideal α-helices. At the same time, molecular dynamics (MD) simulations can provide atomistic insight into peptide-membrane systems. Here, we evaluate the potential of MD simulations to complement the experimental NMR data that is available on three exemplary systems: the natural antimicrobial peptide PGLa and the two designer-made peptides MSI-103 and KIA14, whose sequences were derived from PGLa. Each peptide was simulated for 1 µs in a DMPC lipid bilayer. We calculated from the MD simulations the local angles which define the side chain geometry with respect to the peptide helix. The peptide orientation was then calculated (i) directly from the simulation, (ii) from back-calculated MD-derived NMR splittings, and (iii) from experimental 2H NMR splittings. Our findings are that (1) the membrane orientation and secondary structure of the peptides found in the NMR analysis are generally well reproduced by the simulations; (2) the geometry of the side chains with respect to the helix backbone can deviate significantly from the ideal structure depending on the specific residue, but on average all side chains have the same orientation; and (3) for all of our peptides, the azimuthal rotation angle found from the MD-derived splittings is about 15° smaller than the experimental value.

Structural determinants of the SINE B2 element embedded in the long non-coding RNA activator of translation AS Uchl1.

Pervasive transcription of mammalian genomes leads to a previously underestimated level of complexity in gene regulatory networks. Recently, we have identified a new functional class of natural and synthetic antisense long non-coding RNAs (lncRNA) that increases translation of partially overlapping sense mRNAs. These molecules were named SINEUPs, as they require an embedded inverted SINE B2 element for their UP-regulation of translation. Mouse AS Uchl1 is the representative member of natural SINEUPs. It was originally discovered for its role in increasing translation of Uchl1 mRNA, a gene associated with neurodegenerative diseases. Here we present the secondary structure of the SINE B2 Transposable Element (TE) embedded in AS Uchl1. We find that specific structural regions, containing a short hairpin, are required for the ability of AS Uchl1 RNA to increase translation of its target mRNA. We also provide a high-resolution structure of the relevant hairpin, based on NMR observables. Our results highlight the importance of structural determinants in embedded TEs for their activity as functional domains in lncRNAs.

Real Cost of Speed: The Effect of a Time-Saving Multiple-Time-Stepping Algorithm on the Accuracy of Molecular Dynamics Simulations.

To enhance efficiency in molecular dynamics simulations, forces that vary slowly are often evaluated less often than those that vary rapidly. We show that the multiple-time-step algorithm implemented in recent versions of GROMACS results in significant differences in the collective properties of a system under conditions where the system was previously stable. The implications of changing the simulation algorithm without assessment of potential artifacts on the parametrization and transferability of effective force fields are discussed.