Phenylalanine transfer RNA: molecular dynamics simulation

Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods

Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.

Understanding the mechanistic basis of non-coding RNA through molecular dynamics simulations

Noncoding RNA (ncRNA) has a key role in regulating gene expression, mediating fundamental processes and diseases via a variety of yet unknown mechanisms. Here, we review recent applications of conventional and enhanced Molecular Dynamics (MD) simulations methods to address the mechanistic function of large biomolecular systems that are tightly involved in the ncRNA function and that are of key importance in life sciences. This compendium focuses of three biomolecular systems, namely the CRISPR-Cas9 genome editing machinery, group II intron ribozyme and the ribonucleoprotein complex of the spliceosome, which edit and process ncRNA. We show how the application of a novel accelerated MD simulations method has been key in disclosing the conformational transitions underlying RNA binding in the CRISPR-Cas9 complex, suggesting a mechanism for RNA recruitment and clarifying the conformational changes required for attaining genome editing. As well, we discuss the use of mixed quantum-classical MD simulations in deciphering the catalytic mechanism of RNA splicing as operated by group II intron ribozyme, one of the largest ncRNA structures crystallized so far. Finally, we debate the future challenges and opportunities in the field, discussing the recent application of MD simulations for unraveling the functional biophysics of the spliceosome, a multi-mega Dalton complex of proteins and small nuclear RNAs that performs RNA splicing in humans. This showcase of applications highlights the current talent of MD simulations to dissect atomic-level details of complex biomolecular systems instrumental for the design of finely engineered genome editing machines. As well, this review aims at inspiring future investigations of several other ncRNA regulatory systems, such as micro and small interfering RNAs, which achieve their function and specificity using RNA-based recognition and targeting strategies.

Computer-based techniques for lead identification and optimization I: Basics

This chapter focuses on computational techniques for identifying and optimizing lead molecules, with a special emphasis on natural compounds. A number of case studies have been specifically discussed, such as the case of the naphthyridine scaffold, discovered through a structure-based virtual screening (SBVS) and proposed as the starting point for further lead optimization process, to enhance its telomeric RNA selectivity. Another example is the case of Liphagal, a tetracyclic meroterpenoid extracted fromAka coralliphaga, known as PI3Kα inhibitor, provide an evidence for the design of new active congeners against PI3Kα using molecular dynamics (MD) simulations. These are only two of the numerous examples of the computational techniques’ powerful in drug design and drug discovery fields. Finally, the design of drugs that can simultaneously interact with multiple targets as a promising approach for treating complicated diseases has been reported. An example of polypharmacological agents are the compounds extracted from mushrooms identified by means of molecular docking experiments. This chapter may be a useful manual of molecular modeling techniques used in the lead-optimization and lead identification processes.

RNA folding pathways in stop motion

We introduce a method for predicting RNA folding pathways, with an application to the most important RNA tetraloops. The method is based on the idea that ensembles of three-dimensional fragments extracted from high-resolution crystal structures are heterogeneous enough to describe metastable as well as intermediate states. These ensembles are first validated by performing a quantitative comparison against available solution nuclear magnetic resonance (NMR) data of a set of RNA tetranucleotides. Notably, the agreement is better with respect to the one obtained by comparing NMR with extensive all-atom molecular dynamics simulations. We then propose a procedure based on diffusion maps and Markov models that makes it possible to obtain reaction pathways and their relative probabilities from fragment ensembles. This approach is applied to study the helix-to-loop folding pathway of all the tetraloops from the GNRA and UNCG families. The results give detailed insights into the folding mechanism that are compatible with available experimental data and clarify the role of intermediate states observed in previous simulation studies. The method is computationally inexpensive and can be used to study arbitrary conformational transitions.

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.

Stacking in RNA: NMR of Four Tetramers Benchmark Molecular Dynamics

Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against 1H-1H NOESY distances and 3J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.

Structural Insights into tRNA Dynamics on the Ribosome

High-resolution structures at different stages, as well as biochemical, single molecule and computational approaches have highlighted the elasticity of tRNA molecules when bound to the ribosome. It is well acknowledged that the inherent structural flexibility of the tRNA lies at the heart of the protein synthesis process. Here, we review the recent advances and describe considerations that the conformational changes of the tRNA molecules offer about the mechanisms grounded in translation.

AWE-WQ: fast-forwarding molecular dynamics using the accelerated weighted ensemble

A limitation of traditional molecular dynamics (MD) is that reaction rates are difficult to compute. This is due to the rarity of observing transitions between metastable states since high energy barriers trap the system in these states. Recently the weighted ensemble (WE) family of methods have emerged which can flexibly and efficiently sample conformational space without being trapped and allow calculation of unbiased rates. However, while WE can sample correctly and efficiently, a scalable implementation applicable to interesting biomolecular systems is not available. We provide here a GPLv2 implementation called AWE-WQ of a WE algorithm using the master/worker distributed computing WorkQueue (WQ) framework. AWE-WQ is scalable to thousands of nodes and supports dynamic allocation of computer resources, heterogeneous resource usage (such as central processing units (CPU) and graphical processing units (GPUs) concurrently), seamless heterogeneous cluster usage (i.e., campus grids and cloud providers), and support for arbitrary MD codes such as GROMACS, while ensuring that all statistics are unbiased. We applied AWE-WQ to a 34 residue protein which simulated 1.5 ms over 8 months with peak aggregate performance of 1000 ns/h. Comparison was done with a 200 μs simulation collected on a GPU over a similar timespan. The folding and unfolded rates were of comparable accuracy.

Influence of Sequence and Covalent Modifications on Yeast tRNA Dynamics

Modified nucleotides are prevalent in tRNA. Experimental studies reveal that these covalent modifications play an important role in tuning tRNA function. In this study, molecular dynamics (MD) simulations were used to investigate how modifications alter tRNA dynamics. The X-ray crystal structures of tRNA(Asp), tRNA(Phe), and tRNA(iMet), both with and without modifications, were used as initial structures for 333 ns explicit solvent MD simulations with AMBER. For each tRNA molecule, three independent trajectory calculations were performed, giving an aggregate of 6 μs of total MD across six molecules. The global root-mean-square deviations (RMSD) of atomic positions show that modifications only introduce significant rigidity to the global structure of tRNA(Phe). Interestingly, RMSDs of the anticodon stem-loop (ASL) suggest that modified tRNA has a more rigid structure compared to the unmodified tRNA in this domain. The anticodon RMSDs of the modified tRNAs, however, are higher than those of corresponding unmodified tRNAs. These findings suggest that the rigidity of the anticodon stem-loop is finely tuned by modifications, where rigidity in the anticodon arm is essential for tRNA translocation in the ribosome, and flexibility of the anticodon is important for codon recognition. Sugar pucker and water residence time of pseudouridines in modified tRNAs and corresponding uridines in unmodified tRNAs were assessed, and the results reinforce that pseudouridine favors the 3′-endo conformation and has a higher tendency to interact with water. Principal component analysis (PCA) was used to examine correlated motions in tRNA. Additionally, covariance overlaps of PCAs were compared for trajectories of the same molecule and between trajectories of modified and unmodified tRNAs. The comparison suggests that modifications alter the correlated motions. For the anticodon bases, the extent of stacking was compared between modified and unmodified molecules, and only unmodified tRNA(Asp) has significantly higher percentage of stacking time. Overall, the simulations reveal that the effect of covalent modification on tRNA dynamics is not simple, with modifications increasing flexibility in some regions of the structure and increasing rigidity in other regions.

Examinations of tRNA Range of Motion Using Simulations of Cryo-EM Microscopy and X-Ray Data

We examined tRNA flexibility using a combination of steered and unbiased molecular dynamics simulations. Using Maxwell's demon algorithm, molecular dynamics was used to steer X-ray structure data toward that from an alternative state obtained from cryogenic-electron microscopy density maps. Thus, we were able to fit X-ray structures of tRNA onto cryogenic-electron microscopy density maps for hybrid states of tRNA. Additionally, we employed both Maxwell's demon molecular dynamics simulations and unbiased simulation methods to identify possible ribosome-tRNA contact areas where the ribosome may discriminate tRNAs during translation. Herein, we collected >500 ns of simulation data to assess the global range of motion for tRNAs. Biased simulations can be used to steer between known conformational stop points, while unbiased simulations allow for a general testing of conformational space previously unexplored. The unbiased molecular dynamics data describes the global conformational changes of tRNA on a sub-microsecond time scale for comparison with steered data. Additionally, the unbiased molecular dynamics data was used to identify putative contacts between tRNA and the ribosome during the accommodation step of translation. We found that the primary contact regions were H71 and H92 of the 50S subunit and ribosomal proteins L14 and L16.

RNA stability under different combinations of amber force fields and solvation models

The proper matching of force field and solvent is critical to obtain correct result in molecular dynamics simulation of bio-molecules. This problem has been intensively investigated for protein but not for RNA yet. In this paper, we use standard molecular dynamics and replica exchange molecular dynamics to take a series of tests on the RNA stability under different combinations of Amber force field parameters (ff98, ff99 and ff99bsc0) and the general Born implicit solvent models (igb1, igb2 and igb5). It is found that only ff98 and ff99bsc0 with igb1 can keep the native conformations of RNA hairpin and duplex. Our results suggest that ff98 plus igb1 may be reasonable choice for molecular dynamics simulation of RNA dynamics.

Experimental and computational determination of tRNA dynamics

As the molecular representation of the genetic code, tRNA plays a central role in the translational machinery where it interacts with several proteins and other RNAs during the course of protein synthesis. These interactions exploit the dynamic flexibility of tRNA. In this minireview, we discuss the effects of modified bases, ions, and proteins on tRNA structure and dynamics and the challenges of observing its motions over the cycle of translation.

Flexibility of Short-Strand RNA in Aqueous Solution as Revealed by Molecular Dynamics Simulation: Are A-RNA and A´-RNA Distinct Conformational Structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A′-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA – namely major groove width, inclination and the number of base pairs in a helical twist – over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A′-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors).

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. .

Spinach on the ceiling: a theoretical chemist's return to biology

I was born in Vienna and came to the United States as a refugee in October 1938. This experience played an important role in my view of the world and my approach to science: It contributed to my realization that it was safe to stop working in fields that I felt I understood and to focus on different areas of research by asking questions that would teach me and others something new. I describe my experiences that led me from chemistry and physics back to my first love, biology, and outline some of the contributions I have made as part of my ongoing learning experience.

Incorporation of aminoacyl-tRNA into the ribosome as seen by cryo-electron microscopy

Aminoacyl-tRNAs (aa-tRNAs) are delivered to the ribosome as part of the ternary complex of aa-tRNA, elongation factor Tu (EF-Tu) and GTP. Here, we present a cryo-electron microscopy (cryo-EM) study, at a resolution of approximately 9 A, showing that during the incorporation of the aa-tRNA into the 70S ribosome of Escherichia coli, the flexibility of aa-tRNA allows the initial codon recognition and its accommodation into the ribosomal A site. In addition, a conformational change observed in the GTPase-associated center (GAC) of the ribosomal 50S subunit may provide the mechanism by which the ribosome promotes a relative movement of the aa-tRNA with respect to EF-Tu. This relative rearrangement seems to facilitate codon recognition by the incoming aa-tRNA, and to provide the codon-anticodon recognition-dependent signal for the GTPase activity of EF-Tu. From these new findings we propose a mechanism that can explain the sequence of events during the decoding of mRNA on the ribosome.

Molecular dynamics of biological macromolecules: a brief history and perspective

A description of the origin of my interest in and the development of molecular dynamics simulations of biomolecules is presented with a historical overview, including the role of my interactions with Shneior Lifson and his group in Israel. Some early applications of the methodology by members of my group are summarized, followed by a description of examples of recent applications and some discussion of possible future directions.

Molecular dynamics of the anticodon domain of yeast tRNA(Phe): codon-anticodon interaction

We have studied the effect of codon-anticodon interaction on the structure and dynamics of transfer RNAs using molecular dynamics simulations over a nanosecond time scale. From our molecular dynamical investigations of the solvated anticodon domain of yeast tRNA(Phe) in the presence and absence of the codon trinucleotides UUC and UUU, we find that, although at a gross level the structures are quite similar for the free and the bound domains, there are small but distinct differences in certain parts of the molecule, notably near the Y37 base. Comparison of the dynamics in terms of interatomic or inter-residual distance fluctuation for the free and the bound domains showed regions of enhanced rigidity in the loop region in the presence of codons. Because fluorescence experiments suggested the existence of multiple conformers of the anticodon domain, which interconvert on a much larger time scale than our simulations, we probed the conformational space using five independent trajectories of 500 ps duration. A generalized ergodic measure analysis of the trajectories revealed that at least for this time scale, all the trajectories populated separate parts of the conformational space, indicating a need for even longer simulations or enhanced sampling of the conformational space to give an unequivocal answer to this question.

A molecular dynamics simulation of the flavin mononucleotide-RNA aptamer complex

We report on an unrestrained molecular dynamics simulation of the flavin mononucleotide (FMN)-RNA aptamer. The simulated average structure maintains both cross-strand and intermolecular FMN-RNA nuclear Overhauser effects from the nmr experiments and has all qualitative features of the nmr structure including the G10-U12-A25 base triple and the A13-G24, A8-G28, and G9-G27 mismatches. However, the relative orientation of the hairpin loop to the remaining part of the molecule differs from the nmr structure. The simulation predicts that the flexible phosphoglycerol part of FMN moves toward G27 and forms hydrogen bonds. There are structurally long-lived water molecules in the FMN binding pocket forming hydrogen bonds within FMN and between FMN and RNA. In addition, long-lived water is found bridging primarily RNA backbone atoms. A general feature of the environment of long-lived "structural" water is at least two and in most cases three or four potential acceptor atoms. The 2'-OH group of RNA usually acts as an acceptor in interactions with the solvent. There are almost no intrastrand O2'H(n) vertical ...O4'(n + 1) hydrogen bonds within the RNA backbone. In the standard case the preferred orientation of the 2'-OH hydrogen atoms is approximately toward O3' of the same nucleotide. However, a relatively large number of conformations with the backbone torsional angle gamma in the trans orientation is found. A survey of all experimental RNA x-ray structures shows that this backbone conformation occurs but is less frequent than found in the simulation. Experimental nmr RNA aptamer structures have a higher fraction of this conformation as compared to the x-ray structures. The backbone conformation of nucleotide n + 1 with the torsional angle gamma in the trans orientation leads to a relatively short distance between 2'-OH(n) and O5'(n + 1), enabling hydrogen-bond formation. In this case the preferred orientation of the 2'-OH hydrogen atom is approximately toward O5'(n + 1). We find two relatively short and dynamically stable types of backbone-backbone next-neighbor contacts, namely C2'(H)(n) vertical ...O4'(n + 1) and C5'(H)(n + 1) vertical ...O2'(n). These interactions may affect both backbone rigidity and thermodynamic stability of RNA helical structures.

The effect of queuosine on tRNA structure and function

Computational modeling was performed to determine the potential function of the queuosine modification of tRNA found in wobble position 34 of tRNAasp, tRNAasn, tRNAhis, and tRNAtyr. Using the crystal structure of tRNAasp and a tRNA-tRNA-mRNA complex model, we show that the queuosine modification serves as a structurally restrictive base for tRNA anticodon loop flexibility. An extended intraresidue and intramolecular hydrogen bonding network is established by queuosine. The quaternary amine of the 7-aminomethyl side chain hydrogen bonds with the base's carbonyl oxygen. This positions the dihydroxycyclopentenediol ring of queuosine in proper orientation for hydrogen bonding with the backbone of the neighboring uridine 33 residue. The interresidue association stabilizes the formation of a cross-loop hydrogen bond between the uridine 33 base and the phosphoribosyl backbone of the cytosine at position 36. Additional interactions between RNAs in the translation complex were studied with regard to potential codon context and codon bias effects. Neither steric nor electrostatic interaction occurs between aminoacyl- and peptidyl-site tRNA anticodon loops that are modified with queuosine. However, there is a difference in the strength of anticodon/codon associations (codon bias) based on the presence or lack of queuosine in the wobble position of the tRNA. Unmodified (guanosine-containing) tRNAasp forms a very stable association with cytosine (GAC), but is much less stable in complex with a uridine-containing codon (GAU). Queuosine-modified tRNAasp exhibits no bias for either of cognate codons GAC or GAU and demonstrates a lower binding energy similar to the wobble pairing of guanosine-containing tRNA with a GAU codon. This is proposed to be due to the inflexibility of the queuosine-modified anticodon loop to accommodate proper positioning for optimal Watson-Crick type associations. A preliminary survey of codon usage patterns in oncodevelopmental versus housekeeping gene transcripts suggests a significant difference in bias for the queuosine-associated codons. Therefore, the queuosine modification may have the potential to influence cellular growth and differentiation by codon bias-based regulation of protein synthesis for discrete mRNA transcripts.

Simulations of reversible protein aggregate and crystal structure

We simulated the structure of reversible protein aggregates as a function of protein surface characteristics, protein-protein interaction energies, and the entropic penalty accompanying the immobilization of protein in a solid phase. These simulations represent an extension of our previous work on kinetically irreversible protein aggregate structure and are based on an explicit accounting of the specific protein-protein interactions that occur within reversible aggregates and crystals. We considered protein monomers with a mixture of hydrophobic and hydrophilic surface regions suspended in a polar solvent; the energetic driving force for aggregation is provided by the burial of solvent-exposed hydrophobic surface area. We analyzed the physical properties of the generated aggregates, including density, protein-protein contact distributions, solvent accessible surface area, porosity, and order, and compared our results with the protein crystallization literature as well as with the kinetically irreversible case. The physical properties of reversible aggregates were consonant with those observed for the irreversible aggregates, although in general, reversible aggregates were more stable energetically and were more crystal-like in their order content than their irreversible counterparts. The reversible aggregates were less dense than the irreversible aggregates, indicating that the increased energetic stability is derived primarily from the optimality rather than the density of the packing in the solid phase. The extent of hydrophobic protein-protein contacts and solvent-exposed surface area within the aggregate phase depended on the aggregation pathway: reversible aggregates tended to have a greater proportion of hydrophobic-hydrophobic contacts and a smaller fraction of hydrophobic solvent-exposed surface area. Furthermore, the arrangement of hydrophobic patches on the protein surface played a major role in the distribution of protein contacts and solvent content. This was readily reflected in the order of the aggregates: the greater the contiguity of the hydrophobic patches on the monomer surface, the less ordered the aggregates became, despite the opportunities for rearrangement offered by a reversible pathway. These simulations have enhanced our understanding of the impact of protein structural motifs on aggregate properties and on the demarcation between aggregation and crystallization.

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.

Internal dynamics of tRNA(Phe) studied by depolarized dynamic light scattering

The collective internal dynamics of transfer RNA(Phe) from brewer's yeast in solution was studied by depolarized dynamic light scattering (DDLS). Within the melting region of tRNA the depolarized spectra consist of two Lorentzian, where the narrow (slow) component describes the overall rotation of the macromolecule. The broad component is attributed to the collective reorientation of the bases within the biopolymer. At high temperature only this relaxation process is observed in the spectrum. The viscosity dependence of the collective internal relaxation process is described by the Stokes-Einstein-Debye equation for rotational diffusion. Estimates of the internal orientational pair correlation factor from the integral depolarized intensities of tRNA(Phe) solutions indicates that the observed dynamics correspond to the collective reorientation of approximately 5 bases. A comparison of the results presented with DDLS studies on the aggregation of the mononucleotide guanosine-5'-monophosphate confirms this result. For a further characterization of the relaxation process we studied the effect of hydrostatic pressure (1-1000 bar) on the depolarized spectra of tRNA. While other spectroscopic methods like nmr, fluorescence polarization anisotropy decay, or ESR give information about the very local motion of a single base within the DNA or RNA, this study shows that by DDLS one can characterize collective internal motions of macromolecules.

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.

Local mobility of nucleic acids as determined from crystallographic data. II. Z-form DNA

Directions and magnitudes of the local mobility of the Z-DNA hexamer duplex CpGpCpGpCpG have been determined by crystallographic refinement of anisotropic displacement parameters using the observed X-ray diffraction data. The cytidine and guanosine residues demonstrate different modes of mobility, implying that a dinucleotide is the smallest repeating unit in terms of flexibility as well as structure. Directions of librational and translational mobility of the cytidine and guanosine residues of Z-DNA are similar to those observed for the same nucleotides in B-DNA. This suggests that the local mobility of DNA is primarily determined by the individual nucleotide type and by the constraints of Watson-Crick base-pairing, rather than by helical form. Differences in the magnitudes of mobility may be responsible for some of the different physical properties of B-DNA and Z-DNA. The B to Z transition is discussed in terms of the observed flexibilities of these two helical forms.

Refinement of a molecular model for lamprey hemoglobin from Petromyzon marinus

A molecular model for the protein and ambient solvent of the complex of cyanide with methemoglobin V from the sea lamprey Petromyzon marinus yields an R-factor of 0.142 against X-ray diffraction data to 2.0 A resolution. The root-mean-square discrepancies from ideal bond length and angle are, respectively, 0.014 A and 1.5 degrees. Atoms that belong to planar groups deviate by 0.012 A from planes determined by a least-squares procedure. The average standard deviation for chiral volumes, peptide torsion angle and torsion angles of side-chains are 0.150 A3, 2.0 degrees and 19.4 degrees, respectively. The root-mean-square variation in the thermal parameters of bonded atoms of the polypeptide backbone is 1.21 A2; the variation in thermal parameters for side-chain atoms is 2.13 A2. The model includes multiple conformations for 11 side-chains of the 149 amino acid residues of the protein. We identify 231 locations as sites of water molecules in full or partial occupancy. The sum of occupancy factors for these sites is approximately 154, representing 28% of the 550 molecules of water within the crystallographic asymmetric unit. The environment of the heme in the cyanide complex of lamprey methemoglobin resembles the deoxy state of the mammalian tetramer. In particular, the bond between atom NE2 of the proximal histidine and the Fe lies 5.1 degrees from the normal of the heme plane. In deoxy- and carbonmonoxyhemoglobins, the deviations from the normal to the heme plane are 7 to 8 degrees and 1 degree, respectively. Furthermore, the inequality in the distance of atom CD2 of the proximal histidine from the pyrrole nitrogen of ring-C of the heme (distance = 3.29 A) and CE1 from the pyrrole nitrogen of ring-A (distance = 3.06 A) is characteristic of deoxyhemoglobin, not carbonmonoxyhemoglobin, where these distances are equal. Finally, a hydrogen bond exists between carbonyl 111 and the hydroxyl of tyrosine 149. The corresponding hydrogen link in the mammalian tetramer is central to the T to R state transition and is present in deoxyhemoglobin but absent in carbonmonoxyhemoglobin. We suggest that the low affinity of oxygen for lamprey hemoglobin may be a consequence of these T-state geometries.

Molecular dynamics simulations of d(C-G-C-G-A) X d(T-C-G-C-G) with and without "hydrated" counterions

We present the results of molecular dynamics simulations on d(C-G-C-G-A) X d(T-C-G-C-G) with fully charged phosphates with and without inclusion of counterions. The average structures found in the two simulations are similar, but the simulation with counterions does give an average helix repeat, tilt, and twist in better agreement with those found in the x-ray structure of d(C-G-C-G-A-A-T-T-C-G-C-G)2. The average sugar pucker phases and amplitudes are in qualitative agreement with those found in NMR studies of double-helical DNA, and a number of examples of sugar repuckering from C2' endo to C3' endo carbon conformations in the sugar ring are found. The hydrogen bond correlations as well as torsion correlations are analyzed, and some interesting long-range correlations between dihedral angles are found.

Influence of solvent accessibility and intermolecular contacts on atomic mobilities in hemerythrins

Thermal factor parameters (B values) have been compared from the refined crystal structures of the myohemerythrin from Themiste zostericola and of the octameric hemerythrin from Themiste dyscrita. These B values, which are directly related to atomic mobilities, were found to correlate rather closely with the solvent accessible areas within the respective crystals. Although protomeric units of the two molecules have exceptionally similar three-dimensional structures, there are marked differences between the patterns of relative atomic mobilities along the polypeptide chains. The differences correspond to lattice and oligomer contacts. An adjustment of the B values based on the fraction of accessible area occluded by contacts yields values that correlate well between the independent subunits and that should pertain more closely to those for the protomer free in solution.