Large-amplitude bending motions in phenylalanine transfer RNA

Prediction of solution properties and dynamics of RNAs by means of Brownian dynamics simulation of coarse-grained models: Ribosomal 5S RNA and phenylalanine transfer RNA

Background

The possibility of validating biological macromolecules with locally disordered domains like RNA against solution properties is helpful to understand their function. In this work, we present a computational scheme for predicting global properties and mimicking the internal dynamics of RNA molecules in solution. A simple coarse-grained model with one bead per nucleotide and two types of intra-molecular interactions (elastic interactions and excluded volume interactions) is used to represent the RNA chain. The elastic interactions are modeled by a set of Hooke springs that form a minimalist elastic network. The Brownian dynamics technique is employed to simulate the time evolution of the RNA conformations.

Results

That scheme is applied to the 5S ribosomal RNA of E. Coli and the yeast phenylalanine transfer RNA. From the Brownian trajectory, several solution properties (radius of gyration, translational diffusion coefficient, and a rotational relaxation time) are calculated. For the case of yeast phenylalanine transfer RNA, the time evolution and the probability distribution of the inter-arm angle is also computed.

Conclusions

The general good agreement between our results and some experimental data indicates that the model is able to capture the tertiary structure of RNA in solution. Our simulation results also compare quite well with other numerical data. An advantage of the scheme described here is the possibility of visualizing the real time macromolecular dynamics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13628-015-0025-7) contains supplementary material, which is available to authorized users.

Conformational energy and structure in canonical and noncanonical forms of tRNA determined by temperature analysis of the rate of s(4)U8-C13 photocrosslinking

Bacterial tRNAs frequently have 4-thiouridine (s(4)U) modification at position 8, which is adjacent to the C13-G22-m(7)G46 base triple in the elbow region of the tRNA tertiary structure. Irradiation with light in the UVA range induces an efficient photocrosslink between s(4)U8 and C13. The temperature dependence of the rate constants for photocrosslinking between the s(4)U8 and C13 has been used to investigate the tRNA conformational energy and structure in Escherichia coli tRNA(Val), tRNA(Phe), and tRNA(fMet) under different conditions. Corrections have been made in the measured rate constants to compensate for differences in the excited state lifetimes due to tRNA identity, buffer conditions, and temperature. The resulting rate constants are related to the rate at which the s(4)U8 and C13 come into the alignment needed for photoreaction; this depends on an activation energy, attributable to the conformational potential energy that occurs during the photoreaction, and on the extent of the structural change. Different photocrosslinking rate constants and temperature dependencies occur in the three tRNAs, and these differences are due both to modest differences in the activation energies and in the apparent s(4)U8-C13 geometries. Analysis of tRNA(Val) in buffers without Mg(2+) indicate a smaller activation energy (~13 kJ mol(-1)) and a larger apparent s(4)U8-C13 distance (~12 A) compared to values for the same parameters in buffers with Mg(2+) (~26 kJ mol(-1) and 0.36 A, respectively). These measurements are a quantitative indication of the strong constraint that Mg(2+) imposes on the tRNA flexibility and structure.

Simulating Protein Motions with Rigidity Analysis

Protein motions, ranging from molecular flexibility to large-scale conformational change, play an essential role in many biochemical processes. Despite the explosion in our knowledge of structural and functional data, our understanding of protein movement is still very limited. In previous work, we developed and validated a motion planning based method for mapping protein folding pathways from unstructured conformations to the native state. In this paper, we propose a novel method based on rigidity theory to sample conformation space more effectively, and we describe extensions of our framework to automate the process and to map transitions between specified conformations. Our results show that these additions both improve the accuracy of our maps and enable us to study a broader range of motions for larger proteins. For example, we show that rigidity-based sampling results in maps that capture subtle folding differences between protein G and its mutants, NuG1 and NuG2, and we illustrate how our technique can be used to study large-scale conformational changes in calmodulin, a 148 residue signaling protein known to undergo conformational changes when binding to Ca(2+). Finally, we announce our web-based protein folding server which includes a publicly available archive of protein motions: (http://parasol.tamu.edu/foldingserver/).

Probing non-selective cation binding in the hairpin ribozyme with Tb(III)

Catalysis by the hairpin ribozyme is stimulated by a wide range of both simple and complex metallic and organic cations. This independence from divalent metal ion binding unequivocally excludes inner-sphere coordination to RNA as an obligatory role for metal ions in catalysis. Hence, the hairpin ribozyme is a unique model to study the role of outer-sphere coordinated cations in folding of a catalytically functional RNA structure. Here, we demonstrate that micromolar concentrations of a deprotonated aqueous complex of the lanthanide metal ion terbium(III), Tb(OH)(aq)(2+), reversibly inhibit the ribozyme by competing for a crucial, yet non-selective cation binding site. Tb(OH)(aq)(2+) also reports a likely location of this binding site through backbone hydrolysis, and permits the analysis of metal binding through sensitized luminescence. We propose that the critical cation-binding site is located at a position within the catalytic core that displays an appropriately-sized pocket and a high negative charge density. We show that cationic occupancy of this site is required for tertiary folding and catalysis, yet the site can be productively occupied by a wide variety of cations. It is striking that micromolar Tb(OH)(aq)(2+) concentrations are compatible with tertiary folding, yet interfere with catalysis. The motif implicated here in cation-binding has also been found to organize the structure of multi-helix loops in evolutionary ancient ribosomal RNAs. Our findings, therefore, illuminate general principles of non-selective outer-sphere cation binding in RNA structure and function that may have prevailed in primitive ribozymes of an early "RNA world".

The morph server: a standardized system for analyzing and visualizing macromolecular motions in a database framework

The number of solved structures of macromolecules that have the same fold and thus exhibit some degree of conformational variability is rapidly increasing. It is consequently advantageous to develop a standardized terminology for describing this variability and automated systems for processing protein structures in different conformations. We have developed such a system as a 'front-end' server to our database of macromolecular motions. Our system attempts to describe a protein motion as a rigid-body rotation of a small 'core' relative to a larger one, using a set of hinges. The motion is placed in a standardized coordinate system so that all statistics between any two motions are directly comparable. We find that while this model can accommodate most protein motions, it cannot accommodate all; the degree to which a motion can be accommodated provides an aid in classifying it. Furthermore, we perform an adiabatic mapping (a restrained interpolation) between every two conformations. This gives some indication of the extent of the energetic barriers that need to be surmounted in the motion, and as a by-product results in a 'morph movie'. We make these movies available over the Web to aid in visualization. Many instances of conformational variability occur between proteins with somewhat different sequences. We can accommodate these differences in a rough fashion, generating an 'evolutionary morph'. Users have already submitted hundreds of examples of protein motions to our server, producing a comprehensive set of statistics. So far the statistics show that the median submitted motion has a rotation of approximately 10 degrees and a maximum Calpha displacement of 17 A. Almost all involve at least one large torsion angle change of >140 degrees. The server is accessible at http://bioinfo.mbb.yale. edu/MolMovDB

Vibrational dynamics of transfer RNAs: comparison of the free and synthetase-bound forms

The vibrational dynamics of transfer RNAs, both free, and complexed with the cognate synthetase, are analyzed using a model (Gaussian network model) which recently proved to satisfactorily describe the collective motions of folded proteins. The approach is similar to a normal mode analysis, with the major simplification that no residue specificity is taken into consideration, which permits us (i) to cast the problem into an analytical form applicable to biomolecular systems including about 10(3 )residues, and (ii) to acquire information on the essential dynamics of such large systems within computational times at least two orders of magnitude shorter than conventional simulations. On a local scale, the fluctuations calculated for yeast tRNAPhe and tRNAAsp in the free state, and for tRNAGln complexed with glutaminyl-tRNA synthetase (GlnRS) are in good agreement with the corresponding crystallographic B factors. On a global scale, a hinge-bending region comprising nucleotides U8 to C12 in the D arm, G20 to G22 in the D loop, and m7G46 to C48 in the variable loop (for tRNAPhe), is identified in the free tRNA, conforming with previous observations. The two regions subject to the largest amplitude anticorrelated fluctuations in the free form, i.e. the anticodon region and the acceptor arm are, at the same time, the regions that experience the most severe suppression in their flexibilities upon binding to synthetase, suggesting that their sampling of the conformational space facilitates their recognition by the synthetase. Likewise, examination of the global mode of motion of GlnRS in the complex indicates that residues 40 to 45, 260 to 270, 306 to 314, 320 to 327 and 478 to 485, all of which cluster near the ATP binding site, form a hinge-bending region controlling the cooperative motion, and thereby the catalytic function, of the enzyme. The distal beta-barrel and the tRNA acceptor binding domain, on the other hand, are distinguished by their high mobilities in the global modes of motion, a feature typical of recognition sites, also observed for other proteins. Most of the conserved bases and residues of tRNA and GlnRS are severely constrained in the global motions of the molecules, suggesting their having a role in stabilizing and modulating the global motion.

Global flexibility of tertiary structure in RNA: yeast tRNAPhe as a model system

The study of RNA structure using x-ray crystallography or NMR has yielded a wealth of detailed structural information; however, such approaches do not generally yield quantitative information regarding long-range flexibility in solution. To address this issue, we describe a solution-based method that is capable of characterizing the global flexibilities of nonhelix elements in RNA, provided that such elements are flanked by helix (e.g., bulges, internal loops, or branches). The "phased tau ratio" method is based on the principle that, for RNA molecules possessing two variably phased bends, the relative birefringence decay times depend on the flexibility of each bend, not simply the mean bend angles. The method is used to examine the overall flexibility of the yeast tRNAPhe core (as unmodified transcript). In the presence of magnesium ions, the tRNA core is not significantly more flexible than an equivalent length of RNA helix. In the absence of divalent ions, the tRNA core gains flexibility under conditions where its secondary structure is likely to be largely preserved. The phased tau ratio approach should be broadly applicable to nonhelix elements in both RNA and DNA and to protein-nucleic acid interactions.

A 5-nanosecond molecular dynamics trajectory for B-DNA: analysis of structure, motions, and solvation

We report the results of four new molecular dynamics (MD) simulations on the DNA duplex of sequence d(CGCGAATTCGCG)2, including explicit consideration of solvent water, and a sufficient number of Na+ counterions to provide electroneutrality to the system. Our simulations are configured particularly to characterize the latest MD models of DNA, and to provide a basis for examining the sensitivity of MD results to the treatment of boundary conditions, electrostatics, initial placement of solvent, and run lengths. The trajectories employ the AMBER 4.1 force field. The simulations use particle mesh Ewald summation for boundary conditions, and range in length from 500 ps to 5.0 ns. Analysis of the results is carried out by means of time series for conformationalm, helicoidal parameters, newly developed indices of DNA axis bending, and groove widths. The results support a dynamically stable model of B-DNA for d(CGCGAATTCGCG)2 over the entire length of the trajectory. The MD results are compared with corresponding crystallographic and NMR studies on the d(CGCGAATTCGCG)2 duplex, and placed in the context of observed behavior of B-DNA by comparisons with the complete crystallographic data base of B-form structures. The calculated distributions of mobile solvent molecules, both water and counterions, are displayed. The calculated solvent structure of the primary solvation shell is compared with the location of ordered solvent positions in the corresponding crystal structure. The results indicate that ordered solvent positions in crystals are roughly twice as structured as bulk water. Detailed analysis of the solvent dynamics reveals evidence of the incorporation of ions in the primary solvation of the minor groove B-form DNA. The idea of localized complexation of otherwise mobile counterions in electronegative pockets in the grooves of DNA helices introduces an additional source of sequence-dependent effects on local conformational, helicoidal, and morphological structure, and may have important implications for understanding the functional energetics and specificity of the interactions of DNA and RNA with regulatory proteins, pharmaceutical agents, and other ligands.

Two structurally different RNA molecules are bound by the spliceosomal protein U1A using the same recognition strategy

BACKGROUND

Human U1A protein binds to hairpin II of U1 small nuclear RNA (snRNA) and, together with other proteins, forms the U1 snRNP essential in pre-mRNA splicing. U1A protein also binds to the 3' untranslated region (3'UTR) of its own pre-mRNA, inhibiting polyadenylation of the 3'end and thereby downregulating its own expression. The 3'UTR folds into an evolutionarily conserved secondary structure with two internal loops; one loop contains the sequence AUUGCAC and the other its variant AUUGUAC. The sequence AUUGCAC is also found in hairpin II of U1 snRNA; hence, U1A protein recognizes the same heptanucleotide sequence in two different structural contexts. In order to better understand the control mechanism of the polyadenylation process, we have built a model of the U1A protein-3'UTR complex based on the crystal structure of the U1A protein-hairpin II RNA complex which we determined previously.

RESULTS

In the crystal structure of the U1A protein-hairpin II RNA complex the AUUGCAC sequence fits tightly into a groove on the surface of U1A protein. The conservation of the heptanucleotide in the 3'UTR strongly suggests that U1A protein forms identical sequence-specific contacts with the heptanucleotide sequence when complexed with the 3'UTR. The crystal structure of the hairpin II complex and the twofold symmetry in the 3'UTR RNA provide sufficient information to restrict the conformation of the 3'UTR RNA and have enabled us to build a model of the 3'UTR complex.

CONCLUSIONS

In the U1A-3'UTR complex, sequence-specific interactions are made entirely by the conserved heptanucleotide and the last base pair (C:G) of the stem. The structure is stabilized by protein-protein contacts and by electrostatic interactions between basic amino acids of the protein and the phosphate backbone of the RNA stem regions. The formation of a protein dimer necessary for the inhibition of poly(A) polymerase requires a conformational change of the C termini of the proteins upon RNA binding. This mechanism could prevent the inhibition of poly(A) polymerase by free U1A protein. The model is consistent with biochemical data, and the protein-protein interactions within the 3'UTR complex account for the cooperativity of U1A protein binding to the 3'UTR. The model also serves as an important structural guide for designing further experiments to understand the interaction between the U1A-3'UTR complex and poly(A) polymerase.

Study of global motions in proteins by weighted masses molecular dynamics: adenylate kinase as a test case

The weighted masses molecular dynamics (WMMD) technique is applied to the protein adenylate kinase. A novel set of restraints has been developed to allow the use of this technique with proteins. The WMMD simulation is successful in predicting the flexibility of the two mobile domains of the protein. The end product of the simulation is similar to the known open and AMP bound forms of the enzyme. The biological relevance of the restraints used and potential methods of improving the technique are discussed.

Dynamics of transfer RNAs analyzed by normal mode calculation

Normal mode calculation is applied to tRNAPhe and tRNAAsp, and their structural and vibrational aspects are analyzed. Dihedral angles along the phosphate-ribose backbone (alpha, beta, gamma, epsilon, zeta) and dihedral angles of glycosyl bonds (chi) are selected as movable parameters. The calculated displacement of each atom agrees with experimental data. In modes with frequencies higher than 130 cm-1, the motions are localized around each stem and the elbow region of the L-shape. On the other hand, collective motions such as bending or twisting of arms are seen in modes with lower frequencies. Hinge axes and bend angles are calculated without prior knowledge. Movements in modes with very low frequencies are combinations of hinge bending motions with various hinge axes and bend angles. The thermal fluctuations of dihedral angles well reflect the structural characters of transfer RNAs. There are some dihedral angles of nucleotides located around the elbow region of L-shape, which fluctuate about five to six times more than the average value. Nucleotides in the position seem to be influential in the dynamics of the entire structure. The normal mode calculation seems to provide much information for the study of conformational changes of transfer RNAs induced by aminoacyl-tRNA synthetase or codon during molecular recognition.

Conformational transitions using molecular dynamics with minimum biasing

The molecular dynamics algorithm (MD), which simulates intramolecular motions on the subnanosecond timescale, has been modified to allow the investigation of slow conformational transitions that do not necessarily occur spontaneously in MD simulations. The method is designated CONTRA MD (CONformational TRAnsitions by Molecular Dynamics with minimum biasing). The method requires the prior definition of a single conformational variable that is required to vary monotonically from an initial conformation to a final target conformation. The simulation is broken up into a series of short free MD segments, and we determine, after each segment of MD, whether or not the system has evolved toward the final conformation. Those segments that do not move the system in that direction are deleted. Those that do move it toward the final conformation are patched together sequentially to generate a single representative trajectory along the transition pathway. The CONTRA MD method is demonstrated first by application to the simultaneous C2'-endo to C3'-endo repucker and anti to syn N-glycosidic torsion transitions in 2'-deoxyadenosine and then to the large-scale bending in phenylalanine transfer RNA.

Models for two tRNAs bound to successive codons on mRNA on the ribosome

We have investigated the structural changes necessary to build a model complex of two molecules of phenylalanine transfer RNA (tRNA(Phe) bound to successive codons in a short segment of a model messenger RNA (mRNA), consisting of U6. We keep the mRNA in an ideal helical conformation, deforming the tRNAs as necessary to eliminate steric overlaps while bringing the two 3' termini together. The resulting model has the two tRNAs oriented relative to one another in a manner that is very similar to a model developed by McDonald and Rein (1) in which the tRNAs maintain their ideal crystallographic conformations and all of the deformations are introduced into the mRNA. Consequently, regardless of how one divides the deformations between the tRNAs and the mRNA it is clear that, on the ribosome, the tRNA in the P site has its "front" side (that side with the variable loop) close to the "back" side of the tRNA in the A site (that side with the D loop). The space between the two molecules must be left free on the ribosome, in order to facilitate the transition from the A site to the P site. A detailed pathway is also proposed for changing the anticodon loop structure from that of the A site to that of the P site. The anticodon loop is always kept in a 3'-stacked conformation, since we find that the shift between the 3'-stacked and 5'-stacked structures proposed by Woese (2) is not feasible.

Mobility of individual 5-fluorouridine residues in 5-fluorouracil-substituted Escherichia coli valine transfer RNA. A 19F nuclear magnetic resonance relaxation study

19F nuclear magnetic resonance (n.m.r.) relaxation parameters of 5-fluorouracil-substituted Escherichia coli tRNA(Val)1 were measured and used to characterize the internal mobility of individual 5-fluorouridine (FUrd) residues in terms of several models of molecular motion. Measured relaxation parameters include the spin-lattice (T1) relaxation time at 282 MHz, the 19F(1H) NOE at 282 MHz, and the spin-spin (T2) relaxation time, estimated from linewidth data at 338 MHz, 282 MHz and 84 MHz. Dipolar and chemical shift anisotropy contributions to the 19F relaxation parameters were determined from the field-dependence of T2. The results demonstrate a large chemical shift anisotropy contribution to the 19F linewidths at 282 and 338 MHz. Analysis of chemical shift anisotropy relaxation data shows that, relative to overall tumbling of the macromolecule, negligible torsional motion occurs about the glycosidic bond of FUrd residues in 19F-labeled tRNA(Val)1, consistent with the maintenance of base-base hydrogen-bond and/or stacking interactions at all fluorouracil residues in the molecule. The dipolar relaxation data are analyzed by using the "two-state jump" and "diffusion in a cone" formalisms. Motional amplitudes (theta) are interpreted as being due to pseudorotational fluctuations within the ribose ring of the fluorinated nucleoside. These amplitudes range from approximately 30 degrees to 60 degrees, assuming a correlation time (tau i,2) of 1.6 ns. By using available 19F n.m.r. assignment data for the 14 FUrd residues in 5-fluorouracil-substituted tRNA(Val)1, these motional amplitudes can be correlated directly with the environmental domain of the residue. Residues located in tertiary and helical structural domains show markedly less motion (theta approximately equal to 30 to 35 degrees) than residues located in loops (theta approximately equal to 45 to 60 degrees). A correlation between residue mobility and solvent exposure is also demonstrated. The amplitudes of internal motion for specific residues agree quite well with those derived from X-ray diffraction and molecular dynamics data for yeast tRNA(Phe).

Conformational dynamics of the anticodon loop in yeast tRNAPhe as sensed by the fluorescence of wybutine

Conformational and dynamic properties of the anticodon loop of yeast tRNAPhe were investigated by analyzing the time resolved fluorescence of wybutine serving as a local structural probe adjacent to the anticodon GmAA on its 3' side. The influence of Mg2+, important for stabilizing the tertiary structure of tRNA, and of the complementary anticodon s2UUC of E. coli tRNA2Glu were investigated. Fluorescence lifetimes and anisotropies were measured with ps time resolution using time correlated single photon counting and a mode locked synchronously pumped and frequency doubled dye laser as excitation source. From the analysis of lifetimes (tau) and rotational relaxation times (tau R) we conclude that wybutine occurs in various structural states: one stacked conformation where the base has no free mobility and the only rotational motion reflects the mobility of the whole tRNA molecule (tau = 6 ns, tau R = 19 ns), an unstacked conformation where the base can freely rotate (tau = 100 ps, tau R = 370 ps) and an intermediary state (tau = 2 ns, tau R = 1.6 ns). Under biological conditions, i.e. in the presence of Mg2+ and neutral salts, wybutine is found in a stacked and immobile state which is consistent with the crystallographic picture. In the presence of the complementary codon however, as exemplified by the E. coli-tRNA2Glu anticodon, our analysis indicates that the codon-anticodon complex exists in an equilibrium of structural states with different rotational mobility of wybutine. The conformation with wybutine freely mobile is the predominant one and suggests that this conformation of the codon-anticodon structure differs from the canonical 3'-5' stack.