Structural changes in 16S RNA from Escherichia coli upon unfolding by urea

Modeling the 3D structure and conformational dynamics of very large RNAs using coarse-grained molecular simulations

We describe a computational approach to building and simulating realistic 3D models of very large RNA molecules (>1000 nucleotides) at a resolution of one "bead" per nucleotide. The method starts with a predicted secondary structure and uses several stages of energy minimization and Brownian dynamics (BD) simulation to build 3D models. A key step in the protocol is the temporary addition of a 4 ^th spatial dimension that allows all predicted helical elements to become disentangled from each other in an effectively automated way. We then use the resulting 3D models as input to Brownian dynamics simulations that include hydrodynamic interactions (HIs) that allow the diffusive properties of the RNA to be modelled as well as enabling its conformational dynamics to be simulated. To validate the dynamics part of the method, we first show that when applied to small RNAs with known 3D structures the BD-HI simulation models accurately reproduce their experimental hydrodynamic radii (Rh). We then apply the modelling and simulation protocol to a variety of RNAs for which experimental Rh values have been reported ranging in size from 85 to 3569 nucleotides. We show that the 3D models, when used in BD-HI simulations, produce hydrodynamic radii that are usually in good agreement with experimental estimates for RNAs that do not contain tertiary contacts that persist even under very low salt conditions. Finally, we show that sampling of the conformational dynamics of large RNAs on timescales of 100 µs is computationally feasible with BD-HI simulations.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

Multiple endonucleases improve MALDI-MS signature digestion product detection of bacterial transfer RNAs

Individual transfer ribonucleic acids (tRNAs) in a complex mixture can be identified by the matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) detection of their signature digestion products. Signature digestion products are endonuclease digestion products whose mass-to-charge value is unique thus corresponding to only a single tRNA. To improve the effectiveness of this approach, we have expanded the applicable endonucleases and examined the use of multiple endonucleases for tRNA identification. The purine specific endonucleases RNase T1 and RNase TA generate the largest number of predicted signature digestion products. Experimentally, MALDI-MS analysis of endonuclease digests from Escherichia coli and Bacillus subtilis finds that any two endonucleases used in combination increases tRNA identification by about 25% over the number identified with a single endonuclease. Using three endonucleases, RNase T1, RNase A, and RNase TA, further improves the number of tRNAs identified by 10-15% over those found with two endonucleases. Limitations in the MALDI-MS approach for complex mixtures were revealed in this study, suggesting that the direct MALDI-MS analysis of signature digestion products is more effective for organisms having 30 or less unique tRNAs.

Evaluation of chemical labeling strategies for monitoring HCV RNA using vibrational microscopy

Raman and coherent anti-Stokes Raman scattering (CARS) microscopies have the potential to aid in detailed longitudinal studies of RNA localization. Here, we evaluate the use of carbon-deuterium and benzonitrile functional group labels as contrast agents for vibrational imaging of hepatitis C virus (HCV) replicon RNA. Dynamic light scattering and atomic force microscopy were used to evaluate the structural consequences of altering HCV subgenomic replicon RNA. Modification with benzonitrile labels caused the replicon RNA tertiary structure to partially unfold. Conversely, deuterium-modified replicon RNA was structurally similar to unmodified replicon RNA. Furthermore, the deuterated replicon RNA provided promising vibrational contrast in Raman imaging experiments. The functional effect of modifying subgenomic HCV replicon RNA was evaluated using the luciferase gene as a genetic reporter of translation. Benzonitrile labeling of the replicon RNA prevented translation in cell-based luciferase assays, while the deuterated replicon RNA retained both translation and replication competency. Thus, while the scattering cross-section for benzonitrile labels was higher, only carbon-deuterium labels proved to be non-perturbative to the function of HCV replicon RNA.

Geometrical factor and physical reasons for its influence on the kinetic and thermodynamic properties of RNA-like heteropolymers

BACKGROUND

It has been shown that not only the primary sequence but also the geometry of the native conformation influences the kinetic and thermodynamic properties of a protein-like model. The purpose of this paper is to elucidate the geometrical factor, which affects the kinetic and thermodynamic properties of RNA-like heteropolymers, and physical reasons for this phenomenon.

RESULTS

It is shown that increasing the strength of long-range contacts can accelerate the finding of the native secondary structure for RNA-like heteropolymers. Decreasing the strength of long-range contacts results in deceleration of this process. The physical reason for this phenomenon is the increase in the average energy of non-native structures and, as a consequence, the increase (in absolute value) of the relative value of native energy from strengthening long-range contacts under constant energy of the native state. Statistical analysis of natural RNAs has shown that the mean stability of helices formed by the pairing of regions that are remote along the chain prevails over the mean stability of helices formed by regions that are close along the chain.

CONCLUSIONS

The kinetics simulations suggest that the folding of the RNA secondary structure depends on the structural factors of the native state. One of these factors accelerating large RNA folding is the existence of strong long-range helices in the native secondary structure.