Search

Conformational changes are important in RNA for binding and catalysis and understanding these changes is important for understanding how RNA functions. Computational techniques using all-atom molecular models can be used to characterize conformational changes in RNA. These techniques are applied to an RNA conformational change involving a single base pair within a nine base pair RNA duplex. The Adenine-Adenine (AA) non-canonical pair in the sequence 5'GGUGAAGGCU3' paired with 3'PCCGAAGCCG5', where P is Purine, undergoes conformational exchange between two conformations on the timescale of tens of microseconds, as demonstrated in a previous NMR solution structure [Chen, G., et al., Biochemistry, 2006. 45: 6889-903]. The more populated, major, conformation was estimated to be 0.5 to 1.3 kcal/mol more stable at 30 °C than the less populated, minor, conformation. Both conformations are trans-Hoogsteen/sugar edge pairs, where the interacting edges on the adenines change with the conformational change. Targeted Molecular Dynamics (TMD) and Nudged Elastic Band (NEB) were used to model the pathway between the major and minor conformations using the AMBER software package. The adenines were predicted to change conformation via intermediates in which they are stacked as opposed to hydrogen-bonded. The predicted pathways can be described by an improper dihedral angle reaction coordinate. Umbrella sampling along the reaction coordinate was performed to model the free energy profile for the conformational change using a total of 1800 ns of sampling. Although the barrier height between the major and minor conformations was reasonable, the free energy difference between the major and minor conformations was the opposite of that expected based on the NMR experiments. Variations in the force field applied did not improve the misrepresentation of the free energies of the major and minor conformations. As an alternative, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approximation was applied to predict free energy differences between the two conformations using a total of 800 ns of sampling. MM-PBSA also incorrectly predicted the major conformation to be higher in free energy than the minor conformation.

Human immunodeficiency virus-1 employs strand transfer for recombination between two viral genomes. We have previously provided evidence that strand transfer proceeds by an invasion-mediated mechanism in which a DNA segment on the original RNA template is invaded by a second RNA template at a gap site. The initial RNA-DNA hybrid then expands until the DNA is fully transferred. Ribonuclease H (RNase H) cleavages and nucleocapsid protein (NC) were required for long-distance propagation of the hybrid. Evaluation was performed on a unique substrate, with a short gap serving as a precreated invasion site. In our current work, this substrate provided an opportunity for us to test what factors influence a specific invasion site to support transfer, and to distinguish factors that influence invasion site creation from those that impact later steps. RNase H can act in a polymerization-dependent or polymerization-independent mode. Polymerization-dependent and polymerization-independent RNase H were found to be important in creating efficiently used invasion sites in the primer-donor complex, with or without NC. Propagation and terminus transfer steps, emanating from a precreated invasion site in the presence of NC, were stimulated by polymerization-dependent, but not polymerization-independent, RNase H. RNase H can carry out primary and secondary cleavages during synthesis. While both modes of cleavage promoted invasion, only primary cleavage promoted propagation in the presence of NC in our system. These observations suggest that once invasion is initiated at a short gap, it can propagate through an adjacent region interrupted only by nicks, with help by NC. We considered the possibility that propagation solely by strand exchange was a significant contributor to transfers. However, it did not promote transfer even if synthetic progress of reverse transcriptase was intentionally slowed, consistent with strand exchange by random walk in which rate declines precipitously with distance.