- Open Access
Kissing loop interaction in adenine riboswitch: insights from umbrella sampling simulations
© Palma et al.; licensee BioMed Central Ltd. 2015
- Published: 1 June 2015
Riboswitches are cis-acting regulatory RNA elements prevalently located in the leader sequences of bacterial mRNA. An adenine sensing riboswitch cis-regulates adeninosine deaminase gene (add) in Vibrio vulnificus. The structural mechanism regulating its conformational changes upon ligand binding mostly remains to be elucidated. In this open framework it has been suggested that the ligand stabilizes the interaction of the distal "kissing loop" complex. Using accurate full-atom molecular dynamics with explicit solvent in combination with enhanced sampling techniques and advanced analysis methods it could be possible to provide a more detailed perspective on the formation of these tertiary contacts.
In this work, we used umbrella sampling simulations to study the thermodynamics of the kissing loop complex in the presence and in the absence of the cognate ligand. We enforced the breaking/formation of the loop-loop interaction restraining the distance between the two loops. We also assessed the convergence of the results by using two alternative initialization protocols. A structural analysis was performed using a novel approach to analyze base contacts.
Contacts between the two loops were progressively lost when larger inter-loop distances were enforced. Inter-loop Watson-Crick contacts survived at larger separation when compared with non-canonical pairing and stacking interactions. Intra-loop stacking contacts remained formed upon loop undocking. Our simulations qualitatively indicated that the ligand could stabilize the kissing loop complex. We also compared with previously published simulation studies.
Discussion and Conclusions
Kissing complex stabilization given by the ligand was compatible with available experimental data. However, the dependence of its value on the initialization protocol of the umbrella sampling simulations posed some questions on the quantitative interpretation of the results and called for better converged enhanced sampling simulations.
- adenine riboswitch
- molecular dynamics
- umbrella sampling
- loop-loop interaction
- kissing complex
In this paper we use atomistic MD with the latest variant of the Amber force field  in combination with enhanced sampling techniques  to provide a more detailed perspective on the formation of the kissing loop complex. The combined approach allows this contribution to be dissected from the other ligand-aptamer interactions and the impact of the ligand on the stability of the loop-loop interaction to be quantified. We reproduce exactly the same protocol that has been used by Allner et al.  in order to perform a fair comparison between the two force-fields on this particular system. Effects of the initialization protocol on the results of umbrella sampling simulations are also discussed in detail.
In this work, we used umbrella sampling (US) simulations  to study the thermodynamics of the kissing loop in the presence and in the absence of the cognate ligand. We enforced the breaking/formation of the loop-loop interaction steering the distance between the two loops and then used the resulting structures as starting conformations for US with multiple restraints . Simulations were carried out with the Gromacs 4.6.3 program package  combined with the PLUMED 2.0 plug-in . All the simulation parameters are discussed in detail in the following subsections.
System set up
All simulations reported hereafter were performed on two systems: the add aptamer domain complexed with adenine (Holo form) and without adenine bound (Apo form). In both cases we used the X-ray structure solved by Serganov et al. [PDB:1Y26] . The ligand was removed to simulate the unbound state. MD simulations were performed using the Amber99 force field  refined with the parmbsc0 α/γ corrections  and the latest χ torsional parameters . The general Amber force field  was used to parametrize the ligand. Partial atomic charges were assigned using the restricted electrostatic potential fit method  based on an electronic structure calculation at the HF/6-31G* level of theory performed with Gaussian03 . The electrostatic interactions were calculated using the particle-mesh Ewald method  and bond-lengths were constrained with LINCS . The systems were set-up following exactly the protocol described in Allnér et al. : aptamers were solvated in a rhombic dodecahedron having 8 nm as box vector lenghts, with a Mg2+-H2O solution using approximately 11000 TIP3 molecules , and a recent parametrization for divalent cations . The 5 crystallographic Mg2+ were initially kept at their respective position, whereas the additional 30 ions added to neutralize the system ([Mg2+] = 0.18 M) were randomly placed. A steepest descent minimization (150 steps) was performed followed by 200 ps of MD at constant temperature (298 K, using stochastic velocity rescaling ) and pressure (1 atm, using the Berendsen barostat ) with positional restraints on both RNA and ions so as to equilibrate water. This procedure was repeated first removing the constraints on the ions and then removing all the remaining constraints. Finally, 12 ns unrestrained simulations at constant volume were performed for each system.
In order to compute the thermodynamic stability of the loop-loop interaction we employed US simulations with multiple harmonic restraints. The distance between the center of mass (CoM) of the backbone atoms of the two loops (L2: bases 20-26; L3: bases 48-54; Figure 2B) was used as a collective variable (CV). We will refer hereafter to this distance as L. 44 uniformly spaced reference values were taken in the range spanning from 12.5 to 34 Å, and restraints with stiffness k = 20 (kcal/mol)/Å2 were employed. In the production phase of the US simulations each of the 44 windows was run for 5 ns. A very important issue in US simulation is the generation of the starting conformations. We here performed two independent US sampling simulations, using starting conformations generated with two different protocols (hereafter referred to as forward and backward). To generate the starting points for the forward US simulations we employed the same protocol used by Allnér et al. . Namely, we ran a series of 44 short (0.25 ns) simulations with a stiffer restraint (k = 40 (kcal/mol)/Å2) keeping the CV at the 44 reference values, where each simulation was initialized from the last frame of the previous one. In this way, before each US window starting structure was sampled, we let the system equilibrate. The reference values were iterated allowing an increasing distance between the loops. The backward US simulation was initialized with an equivalent procedure but iterating the restraints in the opposite order, i.e. starting from the structure with undocked loops. In principle, if US simulations are converged, the result should be independent of the initialization procedure.
The data were analyzed using the last 4 ns of each window. The potential of mean force (PMF) profiles were constructed using the weighted histogram analysis method (WHAM)  implemented by Grossfield  taking the CV values distribution resulting from the US simulations. This implementation of WHAM allows to compute errors with a bootstrapping procedure that assumes uncorrelated samples. To avoid artifacts due to possible correlations we instead adopted a blocking procedure. Namely, we split the final 4 ns of each trajectory in four blocks of 1 ns each and performed the WHAM calculation using only a single block from each simulations. The four resulting profiles are aligned at their CV starting value (12.1 Å for the forward profiles, 34.4 Å for the backward profiles) and error at each point is computed as the standard deviation among the four profiles divided by .
To define the number of stacking interactions and the number of base-pair contacts a local coordinate system was constructed in the centre of each six-membered rings, with the x axis pointing towards the C2 atom and the z axis orthogonal to the ring plane. The pairing and stacking relationship between two bases j and k is based on the vector r jk , i.e. the position of ring center k relative to the coordinate system constructed on base j. The criteria for determining the canonical WC base pairs are: 1) the base pair must be AU or GC; 2) The relative position of the bases is compatible with the geometry of a WC interaction. The latter condition is considered satisfied when the product of the Gaussian function . Mean µ and covariance σ were obtained from the empirical distribution of WC pairs in the crystal structure of the large ribosomal subunit . The criteria for determining the non-canonical base pairs are: 1) the ellipsoidal distance , and ) and ) it is not a WC pair. The criteria for determining the stacking base pairs are: 1) the ellipsoidal distance , and ) and ) and . This procedure yields similar results compared to the MC-annotate software  and was shown to be useful for characterizing both structural and dynamical properties of RNA molecules . The software used to perform this structural analysis is available online (http://github.com/srnas/barnaba).
loops is preserved during undocking. It can be observed that in the Holo simulation the number of intra-strand stacking slightly decreases (≈ 5) for because of the distortion in the structure induced by the one of the two inter-loop WC pairings. After this residual interaction is lost all the intra-loop stacking contacts are recovered.
In order to better assess the convergence of the free energy landscape for kissing complex formation, we also reconstructed the PMF profiles of the Holo and Apo structure from US simulations initialized with the backward process (Figure 2). The forward and backward profiles were aligned at , since at that distance the starting structure of the backward process is equal to the final structure of the forward one. The free-energy change upon docking is estimated taking the difference between the minimum value of the PMF (Holo:; Apo: ) and its value for the undocked structure (): for the Holo kcal/mol, for the Apo kcal/mol. Albeit negative, these numbers are too small and not compatible with the ones found in the forward process. This is a clear signature of hysteresis in the pulling procedure that strongly biases the initial starting points of the US simulation. The reason for this discrepancy can be better understood by performing a structural analysis of the interactions on the different US windows. As it can be seen in Figure 4, in the backward process the native WC base pairs are not reformed. In general, a few contacts are formed between the two loops but they are not enough to stabilize the kissing complex. To be sure that this is a systematic effect we also tried a few alternative settings for backward simulations. Results are presented in the Appendix.
Discussion and Conclusions
Our calculations provide quantitative and atomistic details on the mechanism of kissing loop breaking and formation in the add riboswitch aptamer domain. The results can be compared with those recently obtained by Allnér et al.  on the same system using the CHARMM36 force field [22, 23]. In particular the free-energy computed with the forward process has been obtained with an identical protocol so as to allow a fair comparison between the force fields. In our work the estimated stability of the kissing loop complex is kcal/mol (Holo) and kcal/mol (Apo), so that upon ligand binding kcal/mol > 0. On the contrary, Allnér et al. reported kcal/mol < 0. The sign of ΔΔG indicates whether the ligand binding and the formation of the kissing loop complex are cooperative (positive) or anticooperative (negative). Results obtained with the two force fields thus interestingly lead to two opposite pictures.
Recent experiments probed the differential ligand affinity in aptamers with mutations hindering the formation of the kissing complex . The change in affinity indicates a cooperativity between ligand binding and kissing complex formation. This stabilization has been estimated to be kcal/mol. This number should be interpreted with caution since it is based on the assumption that the mutated aptamer mimics a ligand-bound state that is accessible to the wild type aptamer . Results obtained with Amber force field are in qualitative agreement with this picture.
Recently, the thermodynamics of other stand-alone kissing complexes have also been studied using different biophysical techniques [45, 46]. In these two experimental works the stability of the loop-loop interactions were found to be in the range kcal/mol. Stability depends on the exact sequence and set of intra-loop interactions, but is always on the order of ten kcal/mol. The estimated stability of the kissing loop complex in our calculation, namely kcal/mol and kcal/mol for Holo and Apo respectively, is thus much larger than expected. Results obtained with CHARMM indicate a lower ΔG for both the systems , in better agreement with experimental results, even if still overestimated. Our result could be affected by the known overestimation of stacking interactions in the Amber force field . Additionally, we would like to point out that the overestimation found in our calculations could also be a consequence of difficult convergence in the US simulations. To test if the US simulation are effectively converged, we tried to recover the profiles from simulations initialized with the backward process, with a procedure inspired by two directional pulling in steered MD [48–50]. The discrepancy between forward and backward process is an index of high dependence of the PMF on the initialization procedure and poses some questions on the actual convergence of the US simulations. Similarly to steered MD, one can expect that simulations in the forward and backward process are respectively overestimating and underestimating the kissing complex stability. Optimal results can be obtained in steered MD by combining simulations performed with both protocols [49, 50]. We stress here that even if the forward simulations apparently recover the qualitative behavior of the general accepted model, they cannot be trusted for a quantitative estimation of the free-energy change. The fact that the backward process cannot reach the native docked state is a signature of a barrier in an orthogonal degree of freedom that is not properly sampled. A possible candidate is the barrier related to the desolvation of the loops, required to form the correct interstrand interactions. Additionally, we observe that pulling on the distance between the two loops does not necessarily induce the entropic reduction required upon docking. These issues are expected to affect both forward and backward pulling. Our simulation could not give an estimate of the additional barriers, but we can assume that these issues equally affects the Holo and the Apo systems. Thus, the converged ΔΔG upon ligand binding should be somewhere in between results from the forward and backward simulations. Thus we can expect the ΔΔG to be in a wide range between -2.7 kcal/mol and 17 kcal/mol, which is in qualitative agreement with already mentioned experiments .
The convergence problem is not related to the US method itself but to the difficulty of describing such a complex docking event using a single distance as a CV. This variable is not sufficient to drive the system through the appropriate transition states. This is likely due to the existence of additional barriers on hidden degrees of freedom (e.g. solvation). We believe that in order to reliably quantify the ΔΔG for this system with US or other biased sampling methods one should employ more complex CVs which are closer to the actual reaction coordinate.
In conclusion, in this work we addressed the formation of the kissing loop complex in the A-riboswitch aptamer by means of accurate molecular dynamics simulations in explicit solvent combined with enhanced sampling techniques in presence or absence of the cognate ligand. Results are compatible with experiments and suggest that the ligand stabilizes the kissing loop formation. However, our results also spot some weakness of the umbrella sampling method and call for calculations performed with more advanced techniques, which will be the subject of future investigations.
In order to assess the backward procedure for the US method, we repeated it for the Apo form using a softer restraint (k = 20 (kcal/mol)/Å2). We performed two additional simulations:
A Starting from the final snapshot of the forward procedure explained above we perform a backward procedure with the softer restraint.
B We repeated both the forward and the backward procedures using the softer restraint.
Alessandra Villa is acknowledged for reading the manuscript and providing useful suggestions.
The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013) / ERC Grant Agreement n. 306662, S-RNA-S. Publication of this article has been funded by the same ERC grant.
We acknowledge the CINECA award no. HP10BI6HS2, 2013 under the ISCRA initiative for the availability of high performance computing resources.
This article has been published as part of BMC Bioinformatics Volume 16 Supplement 9, 2015: Proceedings of the Italian Society of Bioinformatics (BITS): Annual Meeting 2014: Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/16/S9.
- Breaker RR: Riboswitches and the RNA world. Cold Spring Harbor Perspectives in Biology. 2012, 4 (2): a003566-PubMed CentralView ArticlePubMedGoogle Scholar
- Serganov A, Nudler E: A decade of riboswitches. Cell. 2013, 152 (1-2): 17-24. 10.1016/j.cell.2012.12.024.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrick JE, Breaker RR: The distributions, mechanisms, and structures of metabolite-binding riboswitches. Genome Biology. 2007, 8 (11): R239-10.1186/gb-2007-8-11-r239.PubMed CentralView ArticlePubMedGoogle Scholar
- Garst AD, Batey RT: A switch in time: detailing the life of a riboswitch. Biochim Biophys Acta. 2009, 1789 (9-10): 584-591. 10.1016/j.bbagrm.2009.06.004.PubMed CentralView ArticlePubMedGoogle Scholar
- Porter EB, Marcano-Velázquez JG, Batey RT: The purine riboswitch as a model system for exploring RNA biology and chemistry. Biochim Biophys Acta. 2014, 1839 (10): 919-930. 10.1016/j.bbagrm.2014.02.014.PubMed CentralView ArticlePubMedGoogle Scholar
- Mandal M, Breaker RR: Adenine riboswitches and gene activation by disruption of a transcription terminator. Nat Struct Mol Biol. 2004, 11 (1): 29-35. 10.1038/nsmb710.View ArticlePubMedGoogle Scholar
- Serganov A, Yuan Y-R, Pikovskaya O, Polonskaia A, Malinina L, Phan AT, et al: Structural basis for discriminative regulation of gene expression by adenine-and guanine-sensing mRNAs. Chem Biol. 2004, 11 (12): 1729-1741. 10.1016/j.chembiol.2004.11.018.View ArticlePubMedGoogle Scholar
- Mandal M, Breaker RR: Gene regulation by riboswitches. Nature Reviews Molecular Cell Biology. 2004, 5 (6): 451-463. 10.1038/nrm1403.View ArticlePubMedGoogle Scholar
- Forsdyke DR: A stem-loop "kissing" model for the initiation of recombination and the origin of introns. Mol Biol Evol. 1995, 12 (5): 949-958.PubMedGoogle Scholar
- Nowakowski J, Tinoco I: RNA structure and stability. Seminars in Virology. 1997, 8 (3): 153-165. 10.1006/smvy.1997.0118.View ArticleGoogle Scholar
- Rieder R, Lang K, Graber D, Micura R: Ligand-induced folding of the adenosine deaminase A-riboswitch and implications on riboswitch translational control. Chembiochem. 2007, 8 (8): 896-902. 10.1002/cbic.200700057.View ArticlePubMedGoogle Scholar
- Lee MK, Gal M, Frydman L, Varani G: Real-time multidimensional NMR follows RNA folding with second resolution. Proc Natl Acad Sci U S A. 2010, 107 (20): 9192-9197. 10.1073/pnas.1001195107.PubMed CentralView ArticlePubMedGoogle Scholar
- Neupane K, Yu H, Foster DA, Wang F, Woodside MT: Single-molecule force spectroscopy of the add adenine riboswitch relates folding to regulatory mechanism. Nucleic Acids Res. 2011, 39 (17): 7677-7687. 10.1093/nar/gkr305.PubMed CentralView ArticlePubMedGoogle Scholar
- Leipply D, Draper DE: Effects of Mg2+ on the free energy landscape for folding a purine riboswitch RNA. Biochemistry. 2011, 50 (14): 2790-2799. 10.1021/bi101948k.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharma M, Bulusu G, Mitra A: MD simulations of ligand-bound and ligand-free aptamer: Molecular level insights into the binding and switching mechanism of the add A-riboswitch. RNA. 2009, 15 (9): 1673-1692. 10.1261/rna.1675809.PubMed CentralView ArticlePubMedGoogle Scholar
- Priyakumar UD, MacKerell AD: Role of the adenine ligand on the stabilization of the secondary and tertiary interactions in the adenine riboswitch. Journal of Molecular Biology. 2010, 396 (5): 1422-1438. 10.1016/j.jmb.2009.12.024.PubMed CentralView ArticlePubMedGoogle Scholar
- Gong Z, Zhao Y, Chen C, Xiao Y: Role of ligand binding in structural organization of add A-riboswitch aptamer: A molecular dynamics simulation. Journal of Biomolecular Structure and Dynamics. 2011, 29 (2): 403-416. 10.1080/07391102.2011.10507394.View ArticlePubMedGoogle Scholar
- Lin JC, Hyeon C, Thirumalai D: Sequence-dependent folding landscapes of adenine riboswitch aptamers. Phys Chem Chem Phys. 2014, 16 (14): 6376-6382. 10.1039/C3CP53932F.View ArticlePubMedGoogle Scholar
- Lin J-C, Thirumalai D: Relative stability of helices determines the folding landscape of adenine riboswitch aptamers. Journal of the American Chemical Society. 2008, 130 (43): 14080-14081. 10.1021/ja8063638.View ArticlePubMedGoogle Scholar
- Allnér O, Nilsson L, Villa A: Loop-loop interaction in an adenine-sensing riboswitch: a molecular dynamics study. RNA. 2013, 19 (7): 916-926. 10.1261/rna.037549.112.PubMed CentralView ArticlePubMedGoogle Scholar
- Di Palma F, Colizzi F, Bussi G: Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch. RNA. 2013, 19 (11): 1517-1524. 10.1261/rna.040493.113.PubMed CentralView ArticlePubMedGoogle Scholar
- Foloppe N, MacKerell AD: All-atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data. Journal of Computational Chemistry. 2000, 21 (2): 86-104. 10.1002/(SICI)1096-987X(20000130)21:2<86::AID-JCC2>3.0.CO;2-G.View ArticleGoogle Scholar
- MacKerell AD, Banavali NK: All-atom empirical force field for nucleic acids: II. Application to molecular dynamics simulations of DNA and RNA in solution. Journal of Computational Chemistry. 2000, 21 (2): 105-120. 10.1002/(SICI)1096-987X(20000130)21:2<105::AID-JCC3>3.0.CO;2-P.View ArticleGoogle Scholar
- Zgarbová M, Otyepka M, Mládek A, Banáš P, Cheatham TE, Jurečka P: Refinement of the Cornell et al. nucleic acids force field based on reference quantum chemical calculations of glycosidic torsion profiles. Journal of Chemical Theory and Computation. 2011, 7 (9): 2886-2902. 10.1021/ct200162x.PubMed CentralView ArticlePubMedGoogle Scholar
- Abrams C, Bussi G: Enhanced sampling in molecular dynamics using metadynamics, replica-exchange, and temperature-acceleration. Entropy. 2014, 16 (1): 163-199.View ArticleGoogle Scholar
- Torrie GM, Valleau JP: Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. Journal of Computational Physics. 1977, 23 (2): 187-199. 10.1016/0021-9991(77)90121-8.View ArticleGoogle Scholar
- Kumar S, Bouzida D, Swendsen RH, Kollman PA, Rosenberg JM: The weighted histogram analysis method for free-energy calculations on biomolecules. I: the method. Journal of Computational Chemistry. 1992, 13 (8): 1011-1021. 10.1002/jcc.540130812.View ArticleGoogle Scholar
- Pronk S, P´áll S, Schulz R, Larsson P, Bjelkmar P, Apostolov R, et al: GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics. 2013, 29 (7): 845-854. 10.1093/bioinformatics/btt055.PubMed CentralView ArticlePubMedGoogle Scholar
- Tribello GA, Bonomi M, Branduardi D, Camilloni C, Bussi G: PLUMED 2: New feathers for an old bird. Computer Physics Communications. 2014, 185 (2): 604-613. 10.1016/j.cpc.2013.09.018.View ArticleGoogle Scholar
- Wang J, Cieplak P, Kollman PA: How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules?. Journal of Computational Chemistry. 2000, 21 (12): 1049-1074. 10.1002/1096-987X(200009)21:12<1049::AID-JCC3>3.0.CO;2-F.View ArticleGoogle Scholar
- Pérez A, Marchán I, Svozil D, Sponer J, Cheatham TE, Laughton CA, Orozco M: Refinement of the AMBER force field for nucleic acids: Improving the description of alpha/gamma conformers. Biophysical Journal. 2007, 92 (11): 3817-3829. 10.1529/biophysj.106.097782.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA: Development and testing of a general amber force field. Journal of Computational Chemistry. 2004, 25 (9): 1157-1174. 10.1002/jcc.20035.View ArticlePubMedGoogle Scholar
- Bayly CI, Cieplak P, Cornell W, Kollman PA: A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP model. Journal of Physical Chemistry. 1993, 97 (40): 10269-10280. 10.1021/j100142a004.View ArticleGoogle Scholar
- Frisch M, Trucks G, Schlegel H, Scuseria G, Robb M, Cheeseman J, et al: Gaussian 03, Revision C.02. 2004, Gaussian, Inc., Wallingford, CTGoogle Scholar
- Darden T, York D, Pedersen L: Particle mesh Ewald: An N·log(N) method for Ewald sums in large systems. Journal of Chemical Physics. 1993, 98: 10089-10.1063/1.464397.View ArticleGoogle Scholar
- Hess B, Bekker H, Berendsen HJC, Fraaije JGEM: LINCS: a linear constraint solver for molecular simulations. Journal of Computational Chemistry. 1997, 18 (12): 1463-1472. 10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-H.View ArticleGoogle Scholar
- Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: Comparison of simple potential functions for simulating liquid water. Journal of Chemical Physics. 1983, 79: 926-10.1063/1.445869.View ArticleGoogle Scholar
- Allnér O, Nilsson L, Villa A: Magnesium ion-water coordination and exchange in biomolecular simulations. J Chem Theory Comput. 2012, 8 (4): 1493-1502. 10.1021/ct3000734.View ArticleGoogle Scholar
- Bussi G, Donadio D, Parrinello M: Canonical sampling through velocity rescaling. Journal of Chemical Physics. 2007, 126 (1): 014101-10.1063/1.2408420.View ArticlePubMedGoogle Scholar
- Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR: Molecular dynamics with coupling to an external bath. Journal of Chemical Physics. 1984, 81: 3684-10.1063/1.448118.View ArticleGoogle Scholar
- Grossfield A: WHAM: an Implementation of the Weighted Histogram Analysis Method. Version 2.0.9, [http://membrane.urmc.rochester.edu/content/wham]
- Klein DJ, Moore PB, Steitz TA: The roles of ribosomal proteins in the structure assembly, and evolution of the large ribosomal subunit. Journal of Molecular Biology. 2004, 340 (1): 141-177. 10.1016/j.jmb.2004.03.076.View ArticlePubMedGoogle Scholar
- Gendron P, Lemieux S, Major F: Quantitative analysis of nucleic acid three-dimensional structures. Journal of Molecular Biology. 2001, 308 (5): 919-936. 10.1006/jmbi.2001.4626.View ArticlePubMedGoogle Scholar
- Bottaro S, Di Palma F, Bussi G: The role of nucleobase interactions in RNA structure and dynamics. Nucleic Acids Research. 2014, 42 (21): 13306-13314. 10.1093/nar/gku972.PubMed CentralView ArticlePubMedGoogle Scholar
- Salim N, Lamichhane R, Zhao R, Banerjee T, Philip J, Rueda D, Feig AL: Thermodynamic and kinetic analysis of an rna kissing interaction and its resolution into an extended duplex. Biophys J. 2012, 102 (5): 1097-1107. 10.1016/j.bpj.2011.12.052.PubMed CentralView ArticlePubMedGoogle Scholar
- Stephenson W, Asare-Okai PN, Chen AA, Keller S, Santiago R, Tenenbaum SA, et al: The essential role of stacking adenines in a two-base-pair rna kissing complex. Journal of the American Chemical Society. 2013, 135 (15): 5602-5611. 10.1021/ja310820h.PubMed CentralView ArticlePubMedGoogle Scholar
- Banáš P, Mládek A, Otyepka M, Zgarbová M, Jurečka P, Svozil D, et al: Can we accurately describe the structure of adenine tracts in B-DNA? reference quantum-chemical computations reveal overstabilization of stacking by molecular mechanics. Journal of Chemical Theory and Computation. 2012, 8 (7): 2448-2460. 10.1021/ct3001238.View ArticlePubMedGoogle Scholar
- Grubmüller H, Heymann B, Tavan P: Ligand-binding-molecular mechanics calculation of the streptavidin biotin rupture force. Science. 1996, 271 (5251): 997-999. 10.1126/science.271.5251.997.View ArticlePubMedGoogle Scholar
- Minh DD, Adib AB: Optimized free energies from bidirectional single-molecule force spectroscopy. Physical Review Letters. 2008, 100 (18): 180602-PubMed CentralView ArticlePubMedGoogle Scholar
- Do TN, Carloni P, Varani G, Bussi G: RNA/peptide binding driven by electrostatics - Insight from bidirectional pulling simulations. J Chem Theory Comput. 2013, 9 (3): 1720-1730. 10.1021/ct3009914.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.