Most cited article - PubMed ID 30199638
QM/MM Calculations on Protein-RNA Complexes: Understanding Limitations of Classical MD Simulations and Search for Reliable Cost-Effective QM Methods
Pseudouridine is the most frequently naturally occurring RNA modification, found in all classes of biologically functional RNAs. Compared to uridine, pseudouridine contains an additional hydrogen bond donor group and is therefore widely regarded as a structure stabilizing modification. However, the effects of pseudouridine modifications on the structure and dynamics of RNAs have so far only been investigated in a limited number of different structural contexts. Here, we introduced pseudouridine modifications into the U-turn motif and the adjacent U:U closing base pair of the neomycin-sensing riboswitch (NSR)-an extensively characterized model system for RNA structure, ligand binding, and dynamics. We show that the effects of replacing specific uridines with pseudouridines on RNA dynamics crucially depend on the exact location of the replacement site and can range from destabilizing to locally or even globally stabilizing. By using a combination of NMR spectroscopy, MD simulations and QM calculations, we rationalize the observed effects on a structural and dynamical level. Our results will help to better understand and predict the consequences of pseudouridine modifications on the structure and function of biologically important RNAs.
- Keywords
- MD simulations, NMR, RNA structure, U-turn, U:U base pair, pseudouridine,
- MeSH
- Nucleic Acid Conformation MeSH
- Base Pairing MeSH
- Pseudouridine * genetics MeSH
- RNA * genetics chemistry MeSH
- Uridine MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
- Names of Substances
- Pseudouridine * MeSH
- RNA * MeSH
- Uridine MeSH
RNA-peptide/protein interactions have been of utmost importance to life since its earliest forms, reaching even before the last universal common ancestor (LUCA). However, the ancient molecular mechanisms behind this key biological interaction remain enigmatic because extant RNA-protein interactions rely heavily on positively charged and aromatic amino acids that were absent (or heavily under-represented) in the early pre-LUCA evolutionary period. Here, an RNA-binding variant of the ribosomal uL11 C-terminal domain was selected from an approximately 1010 library of partially randomized sequences, all composed of ten prebiotically plausible canonical amino acids. The selected variant binds to the cognate RNA with a similar overall affinity although it is less structured in the unbound form than the wild-type protein domain. The variant complex association and dissociation are both slower than for the wild-type, implying different mechanistic processes involved. The profile of the wild-type and mutant complex stabilities along with molecular dynamics simulations uncovers qualitative differences in the interaction modes. In the absence of positively charged and aromatic residues, the mutant uL11 domain uses ion bridging (K+/Mg2+) interactions between the RNA sugar-phosphate backbone and glutamic acid residues as an alternative source of stabilization. This study presents experimental support to provide a new perspective on how early protein-RNA interactions evolved, where the lack of aromatic/basic residues may have been compensated by acidic residues plus metal ions.
- Keywords
- RNA–protein interaction, genetic code evolution, mRNA-display, protein evolution,
- MeSH
- Amino Acids * genetics MeSH
- Ions MeSH
- RNA * genetics MeSH
- Molecular Dynamics Simulation MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
- Names of Substances
- Amino Acids * MeSH
- Ions MeSH
- RNA * MeSH
Molecular dynamics (MD) simulations became a leading tool for investigation of structural dynamics of nucleic acids. Despite recent efforts to improve the empirical potentials (force fields, ffs), RNA ffs have persisting deficiencies, which hamper their utilization in quantitatively accurate simulations. Previous studies have shown that at least two salient problems contribute to difficulties in the description of free-energy landscapes of small RNA motifs: (i) excessive stabilization of the unfolded single-stranded RNA ensemble by intramolecular base-phosphate and sugar-phosphate interactions and (ii) destabilization of the native folded state by underestimation of stability of base pairing. Here, we introduce a general ff term (gHBfix) that can selectively fine-tune nonbonding interaction terms in RNA ffs, in particular, the H bonds. The gHBfix potential affects the pairwise interactions between all possible pairs of the specific atom types, while all other interactions remain intact; i.e., it is not a structure-based model. In order to probe the ability of the gHBfix potential to refine the ff nonbonded terms, we performed an extensive set of folding simulations of RNA tetranucleotides and tetraloops. On the basis of these data, we propose particular gHBfix parameters to modify the AMBER RNA ff. The suggested parametrization significantly improves the agreement between experimental data and the simulation conformational ensembles, although our current ff version still remains far from being flawless. While attempts to tune the RNA ffs by conventional reparametrizations of dihedral potentials or nonbonded terms can lead to major undesired side effects, as we demonstrate for some recently published ffs, gHBfix has a clear promising potential to improve the ff performance while avoiding introduction of major new imbalances.
- MeSH
- RNA chemistry MeSH
- Molecular Dynamics Simulation * MeSH
- Hydrogen Bonding MeSH
- Publication type
- Journal Article MeSH
- Names of Substances
- RNA MeSH