Conformational pathway for the kissing complex-->extended dimer transition of the SL1 stem-loop from genomic HIV-1 RNA as monitored by targeted molecular dynamics techniques

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Electric-Field-Induced Protein Translocation via a Conformational Transition in SecDF: An MD Study

SecDF is an important component of the Sec protein translocation machinery embedded in the bacterial membrane, which is associated with many functions, such as stabilizing other Sec translocon components within the membrane, maintaining the transmembrane (TM) potential, and facilitating the ATP-independent stage of the translocation mechanism. Related studies suggest that SecDF undergoes functionally important conformational changes that involve mainly its P1-head domain and that these changes are coupled with the proton motive force (Δp). However, there still is not a clear understanding of how SecDF functions, its exact role in the translocation machinery, and how its function is related to Δp. Here, using all-atom molecular dynamics simulations combined with umbrella sampling, we study the P1-head conformational change and how it is coupled to the proton motive force. We report potentials of mean force along a root-mean-square-distance-based reaction coordinate obtained in the presence and absence of the TM electrical potential. Our results show that the interaction of the P1 domain dipole moment with the TM electrical field considerably lowers the free-energy barrier in the direction of F-form to I-form transition.

The study of interactions between DNA and PcrA DNA helicase by using targeted molecular dynamic simulations

DNA helicases are important enzymes involved in all aspects of nucleic acid metabolism, ranging from DNA replication and repair to recombination, rescue of stalled replication and translation. DNA helicases are molecular motors. Through conformational changes caused by ATP hydrolysis and binding, they move along the template double helix, break the hydrogen bonds between the two strands and separate the template chains, so that the genetic information can be accessed. In this paper, targeted molecular dynamic simulations were performed to study the important interactions between DNA and PcrA DNA helicase, which can not be observed from the crystal structures. The key residues on PcrA DNA helicase that have strong interactions with both double stranded DNA (ds-DNA) and single stranded DNA (ss-DNA) have been identified, and it was found that such interactions mostly exist between the protein and DNA backbone, which indicates that the translocation of PcrA is independent of the DNA sequence. The simulations indicate that the ds-DNA is separated upon ATP rebinding, rather than ATP hydrolysis, which suggests that the two strokes in the mechanism have two different major roles. Firstly, in the power stroke (ATP hydrolysis), most of the translocations of the bases from one pocket to the next occur. In the relaxation stroke (ATP binding), most of the 'work' is being done to 'melt' the DNA at the separation fork. Therefore, we propose a mechanism whereby the translocation of the ss-DNA is powered by ATP hydrolysis and the separation of the ds-DNA is powered by ATP binding.

The inhibitory effects of vinylphosphonate-linked thymidine dimers on the unidirectional translocation of PcrA helicase along DNA: a molecular modelling study

The PcrA DNA helicases are important bacterial enzymes and quintessential examples of molecular motors. Through conformational changes caused by ATP hydrolysis, they move along the template double helix, breaking the hydrogen bonds holding the two strands together, and separating the template chains so that the genetic information can be accessed. The flexibility of the DNA backbone is essential for the unidirectional translocation of PcrA. A modified DNA substrate with reduced backbone rotational flexibility (via an incorporated vinylphosphonate linkage) has previously been designed and tested as a helicase substrate. The results show that a single modification on the backbone is sufficient to inhibit the activity of PcrA. In this paper a range of molecular simulation methods have been applied to examine the structural origins of this inhibitory effect, as it tests our theories of the mechanism of action of this motor. We observe that the chemical modification has different effects on the energetics of DNA translocation through the protein as it reaches different sub-sites.

The Amber ff99 Force Field Predicts Relative Free Energy Changes for RNA Helix Formation

The ability of the Amber ff99 force field to predict relative free energies of RNA helix formation was investigated. The test systems were three hexaloop RNA hairpins with identical loops and varying stems. The potential of mean force of stretching the hairpins from the native state to an extended conformation was calculated with umbrella sampling. Because the hairpins have identical loop sequence, the differences in free energy changes are only from the stem composition. The Amber ff99 force field was able to correctly predict the order of stabilities of the hairpins, although the magnitude of the free energy change is larger than that determined by optical melting experiments. The two measurements cannot be compared directly because the unfolded state in the optical melting experiments is a random coil, while the end state in the umbrella sampling simulations was an elongated chain. The calculations can be compared to reference data by using a thermodynamic cycle. By applying the thermodynamic cycle to the transitions between the hairpins using simulations and nearest neighbor data, agreement was found to be within the sampling error of simulations, thus demonstrating that ff99 force field is able to accurately predict relative free energies of RNA helix formation.

Insight into the intrinsic flexibility of the SL1 stem-loop from genomic RNA of HIV-1 as probed by molecular dynamics simulation

The SL1 stem-loop is the dimerization initiation site for linking the two copies of the RNA forming the HIV-1 genome. The 26 nucleotides stem contains a defect consisting on a highly conserved G-rich 1-3 asymmetrical internal loop, which is a major site for nucleocapsid protein binding. Several NMR attempts were undertaken to determine the internal loop structure in the SL1 monomer. However, the RNA constructs used in the different studies were largely mutated, in particular with replacement of the nine nucleotides apical loop by a tetraloop, and divergent results were obtained ranging from a rigid structure to a particularly large flexibility. To investigate the reasons for such discrepancies, we used molecular dynamics simulation of the SL1 monomer to probe the effect of mutations on the conformational stability of the internal loop and of the whole stem. It is found that in the wild-type sequence, the internal loop displays conformational variability originating mainly from the nine nucleotides apical loop flexibility that causes large conformational fluctuations (without changing the average structure) in the 7 bp duplex linking the apical and internal loops. The large amplitude atomic motions in the duplex are transmitted to the internal loop in which they induce conformational changes characterized by a labile hydrogen bond network such as G5 successively H-bonded to A29 and G30. On the contrary, with a four nucleotides apical loop, conformational fluctuations in the duplex are reduced by a factor of 2 and are not sufficiently energizing for promoting changes in the internal loop structure at the time scale of the simulations.

Recent advances in the study of nucleic acid flexibility by molecular dynamics

The recent use of molecular dynamics (MD) simulations to study flexibility of nucleic acids has been reviewed from an analysis of the publications appearing in the past two years (from 2005 till date). Despite the existence of some unsolved problems in the methodologies, these years have been witness to major advances in the field. Based on a critical review of the most recent contributions, excitement exists on the expected evolution of the field in the next years.

Understanding the isomerization of the HIV-1 dimerization initiation domain by the nucleocapsid protein

The specific binding of HIV-1 nucleocapsid (NC) to the hinge region of the kissing-loop (KL) dimer formed by stemloop 1 (SL1) can have significant consequences on its ability to isomerize into the corresponding extended duplex (ED) form. The binding determinants and the effects on the isomerization process were investigated in vitro by a concerted strategy involving ad hoc RNA mutants and electrospray ionization-Fourier transform ion cyclotron resonance (ESI-FTICR) mass spectrometry, which enabled us to characterize the stoichiometry and conformational state of all possible protein-RNA and RNA-RNA assemblies present simultaneously in solution. For the first time, NC-hinge interactions were observed in constructs including at least one unpaired guanine at the 5'-end of the loop-loop duplex, whereas no interactions were detected when the unpaired guanine was placed at its 3'-end. This binding mode is supported by the presence of a grip-like motif described by recent crystal structures, which is formed by the 5'-purines of both hairpins held together by mutual stacking interactions. Using tandem mass spectrometry, hinge interactions were clearly shown to reduce the efficiency of KL/ED isomerization without inducing its complete block. This outcome is consistent with the partial stabilization of the extra-helical grip by the bound protein, which can hamper the purine components from parting ways and initiate the strand exchange process. These findings confirm that the broad binding and chaperone activities of NC induce unique effects that are clearly dependent on the structural context of the cognate nucleic acid substrate. For this reason, the presence of multiple binding sites on the different forms assumed by SL1 can produce seemingly contrasting effects that contribute to a fine modulation of the two-step process of RNA dimerization and isomerization.

Resolving fast and slow motions in the internal loop containing stem-loop 1 of HIV-1 that are modulated by Mg2+ binding: role in the kissing-duplex structural transition

Stem loop 1 (SL1) is a highly conserved hairpin in the 5'-leader of the human immunodeficiency virus type I that forms a metastable kissing dimer that is converted during viral maturation into a stable duplex with the aid of the nucleocapsid (NC) protein. SL1 contains a highly conserved internal loop that promotes the kissing-duplex transition by a mechanism that remains poorly understood. Using NMR, we characterized internal motions induced by the internal loop in an SL1 monomer that may promote the kissing-duplex transition. This includes micro-to-millisecond secondary structural transitions that cause partial melting of three base-pairs above the internal loop making them key nucleation sites for exchanging strands and nanosecond rigid-body stem motions that can help bring strands into spatial register. We show that while Mg2+ binds to the internal loop and arrests these internal motions, it preserves and/or activates local mobility at internal loop residues G272 and G273 which are implicated in NC binding. By stabilizing SL1 without compromising the accessibility of G272 and G273 for NC binding, Mg2+ may increase the dependence of the kissing-duplex transition on NC binding thus preventing spontaneous transitions from taking place and ensuring that viral RNA and protein maturation occur in concert.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

Modeling H3 histone N-terminal tail and linker DNA interactions

Molecular dynamics computer simulations were performed for the 25-residue N-terminal tail of the H3 histone protein in the proximity of a DNA segment of 10 base pairs (bp), representing a model for the linker DNA in chromatin. Several least biased configurations were used as initial configurations. The secondary structure content of the protein was increased by the presence of DNA close to it, but the locations of the secondary motifs were different for different initial orientations of the DNA grooves with respect to the protein. As a common feature to all simulations, the electrostatic attraction between negatively charged DNA and positively charged protein was screened by the water solvent and counterbalanced by the intrinsic compaction of the protein due to hydrophobic effects. The protein secondary structure limited the covering of DNA by the protein to 4-5 bp. The degree of compaction and charge density of the bound protein suggests a possible role of H3 tail in a nonspecific bending and plasticity of the linker DNA when the protein is located in the crowded dense chromatin.

Dissecting the protein-RNA and RNA-RNA interactions in the nucleocapsid-mediated dimerization and isomerization of HIV-1 stemloop 1

The specific binding of HIV-1 nucleocapsid protein (NC) to the different forms assumed in vitro by the stemloop 1 (Lai variant) of the genome's packaging signal has been investigated using electrospray ionization-Fourier transform mass spectrometry (ESI-FTMS). The simultaneous observation of protein-RNA and RNA-RNA interactions in solution has provided direct information about the role of NC in the two-step model of RNA dimerization and isomerization. In particular, two distinct binding sites have been identified on the monomeric stemloop structure, corresponding to the apical loop and stem-bulge motifs. These sites share similar binding affinities that are intermediate between those of stemloop 3 (SL3) and the putative stemloop 4 (SL4) of the packaging signal. Binding to the apical loop, which contains the dimerization initiation site (DIS), competes directly with the annealing of self-complementary sequences to form a metastable kissing-loop (KL) dimer. In contrast, binding to the stem-bulge affects indirectly the monomer-dimer equilibrium by promoting the rearrangement of KL into the more stable extended duplex (ED) conformer. This process is mediated by the duplex-melting activity of NC, which destabilizes the intramolecular base-pairs surrounding the KL stem-bulges and enables their exchange to form the inter-strand pairs that define the ED structure. In this conformer, high-affinity binding takes place at stem-bulge sites that are identical to those present in the monomeric and KL forms. In this case, however, the NC-induced "breathing" does not result in dissociation of the double-stranded structure because of the large number of intermolecular base-pairs. The different binding modes manifested by conformer-specific mutants have shown that NC can also provide low affinity interactions with the bulged-out adenine bases flanking the DIS region of the ED conformer, thus supporting the hypothesis that these exposed nucleotides may constitute "base-grips" for protein contacts during the late stages of the viral lifecycle.