
==== Front
Nucleic Acids ResNucleic Acids ResearchNucleic Acids Research0305-10481362-4962Oxford University Press 10.1093/nar/gkl58017012271RNAMolecular dynamics simulations of human tRNAUUULys,3: the role of modified bases in mRNA recognition McCrate Nina E. Varner Mychel E. Kim Kenneth I. Nagan Maria C. *Division of Science, Truman State University, 100 East Normal, KirksvilleMO 63501, USA*To whom correspondence should be addressed. Tel: +1 660 785 4084; Fax: +1 660 785 4045; Email: mnagan@truman.edu11 2006 11 2006 29 9 2006 34 19 5361 5368 18 5 2006 21 7 2006 24 7 2006 © 2006 The Author(s)2006This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License () which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Accuracy in translation of the genetic code into proteins depends upon correct tRNA–mRNA recognition in the context of the ribosome. In human tRNAUUULys,3 three modified bases are present in the anticodon stem–loop—2-methylthio-N6-threonylcarbamoyladenosine at position 37 (ms2t6A37), 5-methoxycarbonylmethyl-2-thiouridine at position 34 (mcm5s2U34) and pseudouridine (ψ) at position 39—two of which, ms2t6A37 and mcm5s2U34, are required to achieve wild-type binding activity of wild-type human tRNAUUULys,3 [C. Yarian, M. Marszalek, E. Sochacka, A. Malkiewicz, R. Guenther, A. Miskiewicz and P. F. Agris (2000) Biochemistry, 39, 13390–13395]. Molecular dynamics simulations of nine tRNA anticodon stem–loops with different combinations of nonstandard bases were performed. The wild-type simulation exhibited a canonical anticodon stair-stepped conformation. The ms2t6 modification at position 37 is required for maintenance of this structure and reduces solvent accessibility of U36. Ms2t6A37 generally hydrogen bonds across the loop and may prevent U36 from rotating into solution. A water molecule does coordinate to ψ39 most of the simulation time but weakly, as most of the residence lifetimes are <40 ps.
==== Body
INTRODUCTION
Nonstandard bases are most prevalently found in tRNA molecules and may involve a variety of modifications from simple methylation to addition of a peptide moiety (1,2). Some tRNA molecules require posttranscriptionally modified nucleic acid bases for proper translation (3,4), recognition by proteins (5,6) and initiation of human immunodeficiency virus reverse transcription (7). Specifically, ribosomal binding (8,9), proper translocation from the aminoacyl site (A) to peptidyl site (P site) (10,11) and translational speed (12) may be dependent upon nonstandard bases, while the lack of modified bases may result in frameshifting errors (13). Most tRNAs contain at least one modified base (14) and modifications often occur at positions 34 and 37 (4). Modified bases are generally thought to increase stability of tRNA structure by providing additional hydrogen-bonding contacts, better base–base stacking interactions (4,15–19) and metal ion binding sites (20–22). Specifically, base modifications at the 34th and 37th positions of the tRNA are often important for facilitating cognate codon recognition (4). Still, the role specific modifications play in determining RNA structure, stability and function are not well understood.

Functional groups in modified nucleosides can dramatically affect the ability of tRNA anticodons to recognize their cognate mRNA codon and therefore participate in accurate translation of the genetic code (15,16). The third human tRNA coding for lysine (tRNALys,3) contains three posttranscriptionally modified bases in the anticodon stem–loop (ASL, Figure 1). Biochemical assays have established that 5-methoxycarbonylmethyl-2-thiouridine at position 34 (mcm3s2U34) and 2-methylthio-N6-threonylcarbamoyl-adenosine at position 37 (ms2t6A37) are required to mimic wild-type human tRNALys,3 binding to programmed ribosomes (9). Pseudouridine (ψ) presence at position 39 is not critical to anticodon–codon recognition but has been shown to contribute to helical stability (17). Escherichia coli tRNALys,, analogous to its human counterpart, contains a 5-methylaminomethyluridine-2-thiouridine at position 34 (mnm5s2U34) and N6-threonylcarbamoyl-adenosine at position 37 (t6A37) that are both important for ribosomal A site binding (8,9,23), translocation from the A site (11) and discrimination between AAA and AAG codons (24,25).

Several structures of human tRNALys,3 and E.coli tRNALys, variants with variously modified bases in a variety of contexts have been determined by both NMR spectroscopy (18,26–29) and X-ray crystallography (19,30). The expensive and often complex synthesis of nonstandard bases further complicates incorporation of all fully modified bases. Structural studies indicate that preservation of a U-turn motif in E.coli anticodon loop is important for recognition of cognate mRNA (28). Systematic, detailed NMR studies of fully modified E.coli and partially modified human tRNALys,3 indicate flexibility in the U36 position (29). In the context of the 30S ribosomal subunit, an anticodon loop containing t6A37 has been found to recognize the cognate mRNA codon (19) in a canonical Fuller–Hodgson tRNA form (31).

The contributions of each modified base to free tRNA and ribosomally bound tRNA structures remains elusive (4,16,32). Molecular dynamics (MD) studies considering specific contributions of modified bases mcm5s2U34, ms2t6A37, ψ39 to the structure and function of the human ASLLys,3 are reported. Specifically, the role of modified bases in preserving canonical tRNA anticodon structure is examined.

METHODS
Initial structures
Fully refined NMR structures are currently not available for fully modified tRNALys,3; therefore all simulations were begun from X-ray crystal structures. MD simulations of the tRNAUUULys,3 anticodon stem–loop with bases 27–43 (ASLLys,3 wild type) derived from the fully modified tRNALys,3 X-ray crystal structure solved to 3.3 Å resolution (PDB code 1FIR) (30) were performed with different combinations of mcm5s2U34, ms2t6A37, ψ39 and unmodified bases (Table 1). All variations of ASLLys,3 were derived from the Bénas et al. (30) structure. While this study was being conducted, the crystal structure of E.coli ASLLys,, which contains t6A37, with the 30S ribosomal subunit and mRNA was solved to 3.0 Å resolution (PDB code 1XMQ) (19). In E.coli ASLLys,, the ureido ring of t6A37 stacks with A38 stabilizing anticodon–codon interactions. To ensure adequate sampling of relevant structures, an additional system, wild-type 2, was derived by grafting the E.coli ASLLys, anticodon loop containing an ureido ring in t6A onto the Bénas ASLLys,3 stem and applying the appropriate full-human ASLLys,3 modifications.

Molecular dynamics simulations
All MD simulations were performed according to the same protocol. Trajectories were propagated with the Cornell et al. force field (33), as implemented in the SANDER module of the program AMBER 6.0 (34). ASLs were solvated by a 10 Å rectangular box of TIP3P water (35) and Na+ ions were added to neutralize the system. Simulations were performed under periodic boundary conditions, employing the particle mesh Ewald (36) formalism to account for long range electrostatics.

MD trajectories were propagated with a 2.0 fs time step, employing the SHAKE algorithm (37) on all hydrogen atoms, with a nonbonded cutoff of 9.0 Å. The pairlist was updated every 25 steps. Constant pressure (1 atm) and temperature (300 K) were maintained according to the Berendsen coupling algorithm (38). The equilibration protocol of Cheatham was followed (39) similar to previous simulations of nucleic acids (40–45). Simulations were allowed to equilibrate under production-run conditions for at least 2.0 ns before collection of statistics over 4.0 ns.

Force-field parameters
Parameters for nonstandard bases were derived to be consistent with the Cornell et al. force field (33) and specific parameters are available in Supplementary Data. Atom types, force constants and equilibrium values were assigned by analogy to the Cornell et al. nucleic acid, threonine side chain (t6) or methionine (ms2) force-field parameters (33), or the General Amber Force Field (46) for sulfur-related parameters. Partial atomic charges were derived from Hartree–Fock (47)/6-31G* (48–50) electrostatic potential charges fitted with a two-step restrained electrostatic potential procedure (51). Electronic structure calculations were performed with Gaussian 03M, revision B.04 (52).

Ms2t6A rotational barrier
Quantum mechanical calculations employing the modified Perdew-Wang 1-parameter model for kinetics (MPW1K) (53,54) and the 6-31+G (d,p) basis set (55,56) were used to determine the rotational barrier around the C10-N11 amide bond involved in ureido ring formation.

Data analysis
Trajectories were analyzed using the CARNAL and PTRAJ modules of AMBER 6.0 (34). All statistics were collected on 1 ps snapshots over 4.0 ns of data. Helical parameters were measured using the 3DNA program (57). Trajectories and structures were visualized using VMD (58) and UCSF chimera (59). Criteria for hydrogen bonding was a 3.5 Å distance between heavy atoms and a heavy atom-H-heavy atom angle of 180° ± 60°. Solvent accessible surface area (SASA) was calculated as a Connolly surface (60), as implemented in the MM-PBSA module of AMBER 8.0 (61). When determining if bases were parallel, the angle between normal vectors to the plane of the respective bases was used. Planes were defined using three points; N1, N3 and C5 for pyrimidines and N3, C6 and C8 for purines and the normal vectors were defined to point in the 3′ direction.

RESULTS
Trajectories
Trajectory stability was based upon steady fluctuations around constant values for root mean square deviation (r.m.s.d.) values from the average (Figure 2) and starting structures (Supplementary Data) as well as monitoring of helical parameters of the double-stranded region such as local twist, x-displacement from the helical axis and base pair inclination (Supplementary Data).

Overall, the helical portion of all ASL structures exhibited classic A-form geometry (62) with x-displacement values ranging from −4.5 ± 0.6 to −5.3 ± 0.9 Å and base pair rise parameters of 2.5 ± 0.3 to 2.8 ± 0.3 Å. The majority of structural analysis occurred in the loop region of the ASLs.

The wild-type simulation was unusually stable with an average r.m.s.d. from the starting structure of only 1.8 ± 0.3 Å. To ensure that the wild-type simulation was indeed naturally stable and that periodic effects or lack of sampling time did not influence the results, data for two additional simulations were collected (N. E. McCrate, unpublished data). In the first simulation, the water box surrounding the RNA was increased to 15 Å beyond the solute. In the other simulation, with 10 Å of water added beyond the solute, data were collected for 40 ns. Both additional simulations retained the structure observed in the reported wild-type simulation.

Presence of mcm5s2U34 and/or ms2t6A37 widens the loop region
All simulations contained a hydrogen bond between C32 O2 or N3 and the A38 exocyclic amine to varying extents of the simulation but averaging 74% occupancy. Simulations containing no modified bases or only ψ39, exhibited an additional hydrogen bond between the exocyclic amine of A37 and the O2 of U33 56 and 62% of the simulations time, respectively. Addition of either mcm5s2U34 and/or ms2t6A37 decreased the A37–U33 interaction to <4% of the simulation time.

The side chain of mcm5s2U34 did not form hydrogen bonds often with other parts of the ASL. The only detectable hydrogen bond observed occurred between the mcm5 carbonyl oxygen and the hydroxyl of ms2t6A37 in the wild-type simulation and the simulation containing both mcm5s2U34 and ms2t6A37 for <20% of the simulation time.

In the models presented herein there seems to be no correlation between presence of a U33 NH3 hydrogen bond to U36 O2P and any one conformation. Simulations containing mcm5s2U34 seem to lose a U33 NH3 hydrogen bond with U36 O2P (<3%) while simulations lacking a modification at position 34 close to form this hydrogen bond >90% of the simulation time. The exceptions to this trend are both wild-type simulations, which contain this hydrogen bond 100 and 91% for wild-type and wild-type 2 simulations, respectively.

Characterization of the anticodon stair-step
In the fully modified ASL structure including mcm5s2U34, ms2t6A37 and ψ39 (wild-type ASLLys,3), the structure retained a stair-stepped conformation, in which the uridines (Figure 3) are offset and nearly parallel to each other, observed previously in the X-ray crystal structures of free human tRNALys,3 (30) and when E.coli t6A ASLLys, pairs with its cognate mRNA codon in the context of the ribosome (19).

To characterize the stair-stepped conformation, the relative positions of the anticodon bases to each other were determined. To analyze the positions of the anticodon bases relative to each other, the distances between the center of mass (COM) of anticodon bases' aromatic rings were measured. The distance vectors were transformed into a relative coordinate system utilizing a rotation matrix and decomposed into x (dx), y (dy) and z (dz) components, to describe the relative positions of the bases. By analogy, the dx, dy and dz components describe the height width and length of ‘the stair’. From visual inspection, the wild-type simulation was determined to exhibit the prototypical stair-stepped conformation. To determine whether the geometry of the anticodon for each model matched that of the wild-type, dx, dy and dz between each set of rings (U34–U35 and U35–U36) were classified as ‘stair stepped’ if the values for both pairs of rings fell within 4 SDs of the average values calculated for the wild-type model and the angle between the normal vectors of the planes of the bases was 0 ± 45° (63) (Table 1).

To aid in data analysis, snapshots of each simulation were assessed for presence of a stair-stepped conformation. Based upon the criteria outlined above, the percent occupancies of stair-stepped interactions for positions 34–36 were determined (Table 2). Simulations that contained an ms2t6 modification at position 37 exhibited the stair-stepped anticodon configuration 75–100% of the simulation time, while those ASLs lacking base modifications at position 37 were not stair-stepped for most, if not all, of the simulation. In the wild-type 2 model, in which the anticodon appears to be stair-stepped, the U36 base is positioned closer to the N9 of ms2t6A37 than is observed in the other models. As a result of this difference in position, the wild-type 2 simulation is classified as having a less stair-stepped codon (46%) than the other models because it is stair-stepped in a slightly different conformation.

When ms2t6A37 is substituted with A37, U36 rotates around the glycosidic bond and is exposed to solvent. This was observed in all simulations lacking ms2t6 modifications at position 37. Quantification of U36 solvent exposure through the SASA of the U36 base indicates that the presence of the ms2t6 modification at position 37 reduces the SASA of U36 (Table 2). Relatively smaller U36 SASA values correspond to the stair-stepped conformation, while larger SASA values correspond to U36 unstacking. In addition, mcm5s2U at position 34 appears to slightly lower the SASA of U36 as all models lacking mcm5s2 modifications have slightly higher U36 SASA than their modified counterparts (10–30 Å2, Table 1).

Base stacking analysis in the anticodon loop
From visual inspection of the wild-type simulations, it seemed that bases 36–38 were interacting, perhaps by base–base stacking interactions. All simulations were examined to determine the extent of base–base stacking in the anticodon loop. A stacking interaction was defined to occur when the distance measured from the COM between two aromatic rings was <4 Å and the angle between the normal vectors of planes formed by the bases was ±20°, as determined appropriate by previous studies (63,64).

Based on the stacking criteria outlined above, the percent occupancy of stacking interactions for positions 36–38 was determined (Table 2). In all simulations containing ms2t6A37, for a portion of the simulation (3–34%), U36 was in close proximity to the A37 5-membered ring. The stacking data indicate that the presence of ms2t6A37 may slightly enhance base–base stacking between U36 and position 37.

Ms2t6A side chains
Modification of A37 encourages the base to remain intercalated between A38 and U36. In general, A37, modified or unmodified, hydrogen bonds with functional groups in the anticodon loop. Simulations containing ms2t6A37 tended to interact across the loop with base and sugar atoms in residues 32–35 for relatively short lifetimes. For instance, in the wild-type trajectory, the hydroxyl group from t6 was hydrogen bonded for most of the simulation time to some group in the tRNA. It preferred to interact with sugar oxygens in residue 33 but also hydrogen bonded with base oxygens in U34 and U35 carbonyls. When ms2t6 modifications were removed, the A37 disrupted the U36–U35 interaction and formed new h-bonds through N6 with bases in closer proximity, such as residues U35, U36 and A38. In the absence of ms2t6, the ASL rearranges to provide new hydrogen-bonding partners for the truncated N6 position of A37.

In the context of the ribosome, E.coli t6A37 contains an ureido ring which stacks with the 6-membered ring of A38 (19). The crystal structure of t6A alone contains an ureido ring (65). In the fully refined tRNALys3 X-ray structure (PDB code 1FIR), ms2t6A37 seems to adopt a conformation without an ureido ring although initial attempts were to try to have a clearer electron density showing an ureido ring-like conformation [P. Bénas, personal communication and figures in Ref. (30)]. Quantum mechanical calculations on N9-methylated ms2t6A at the MPW1K/6-31+G (d,p) level of theory indicate that the barrier to rotation around the C10-N11 bond to form an ureido ring is ∼9.0 kcal/mol at room temperature.

Pseudouridine stabilization
When ψ39 is present, a water bridge forms between A38 O2P, ψ39 NH3 and ψ39 O1P. Although the full water bridge is only present for ∼30% of the production run (detailed tables in Supplementary Data), the component hydrogen bonds between ψ39 NH3-water-A38 O2P and ψ39 NH3-water-ψ39 O1P are present for larger percentages of the production run, typically 50–80%. The waters involved in the water bridge are short lived with the bulk of the bridging interaction being attributable to water molecules having lifetimes at or below 40 ps and nearly all residence lifetimes being attributable to waters with lifetimes ≤120 ps.

DISCUSSION
As has been observed in the E.coli tRNAPhe ASL (66) and various studies involving tRNALys, variants (18,26,28,29), modifications act sterically in the anticodon loop to widen the loop and prevent a closed structure. This is consistent with the previously observed closed-loop conformation observed experimentally by NMR spectroscopy (18,29). Addition of mcm5s2U at position 34 and/or ms2t6A at position 37 introduces enough steric hindrance that the loop widens and interactions between U33 and A37 are no longer observed.

Pseudouridine and water
MD simulations with the Cornell et al. force field (33) and TIP3P (35) water have been shown to identify long-lived, structural water molecules in other systems (41,44). The water molecule that coordinates with ψ39 is present most of the time but is not long lived. Pseudouridine has been shown to stabilize tRNA structure through a water-mediated hydrogen-bonding network (67). The crystal structure of yeast tRNAPhe, which contains a ψ55, has a water molecule bridging the N1-H position to the neighboring phosphate oxygens of T54 and ψ55 (68). Additionally, a long-lived water molecule coordinated to ψ has been found in tRNAAsp through MDs simulations (41). This water-mediated hydrogen bond with ψ is believed to restrict ψ motion and stabilize A–ψ base pairs (41,69,70), which are most prevalently found at positions 31–39 (14).

Canonical tRNA
Overall the wild-type simulation retained the stair-stepped conformation with dx, dy and dz values all similar to that found in the X-ray crystal structure of human tRNALys,3 (30) and E.coli t6A-ASL bound in the 30S ribosome (19). Modifications at position 37 stabilize the anticodon conformation. From this data it appears that the presence of ms2t6A37 encourages retention of a highly ordered stair-stepped conformation.

Base modifications at position 37 seem to stabilize U36 participation in the stair-stepped conformation. U36 is the first base of the anticodon and is required for proper reading of the cognate codon base. The ms2t6A37 base is required for proper positioning of U36 to recognize its cognate codon base. Uridine with its hydrophilic carbonyls and N3, in general, is predisposed to an extra-helical position resulting in an increase in solvent exposure as noted in a systematic study of RNAs containing a single-base bulge (71). The highly ordered stair-stepped conformation of the anticodon bases seems to be important to codon recognition as originally hypothesized (31) and observed in crystal structures of tRNAs alone (72), interacting with other tRNAs through crystal packing effects (30,73–75), and in the context of the ribosome (19,76).

Base–base stacking data indicate that modifications on A37 may position the base to better interact with U36. However, this stacking interaction is not long lived in simulations. In order to have favorable base–base stacking interactions, the interaction of two hydrophobic bases must be enthalpically favorable (77,78). In this case, the negative electron densities of the aromatic bases may repel each other and electrostatic interactions, in which positive portions of the base above are interacting with the negative electron density of the base below, may predominate. This favorable electrostatic interaction has been observed in simple systems, such as benzene and water. In these benzene–water pairs, the partially positive hydrogens of the water molecules are interacting with the negative π-cloud of the benzene ring (79). It may be more favorable for a canonical anticodon to be stair-stepped because the positive portions of bases are aligned to interact with the negative π-cloud of the neighboring, lower base. When A37 is modified, it can partake in hydrogen bonds across the anticodon loop, effectively holding it up above U36. When modifications are removed, the base no longer has enthalpically favorable hydrogen-bonding interactions unless U36 rotates into solution, leaving unmodified A37 to favorably interact with U35 through electrostatic interactions and hydrogen bond with closer bases through its truncated 6-amino group.

Conformation of ms2t6A
The bulky ms2t6 modification containing hydrogen bond donors and acceptors may act sterically to restrict A37 movement while discouraging the displacement of U36 through weak stacking interactions and stabilizing the position of A37 through hydrogen bond contacts across the loop. In the wild-type 2 simulation, where the bulk of the t6A substituent's hydrogen bond donors and acceptors are solvent exposed, U36 exhibits markedly increased conformational freedom. This behavior suggests a delicate balance of forces, of which cross-strand hydrogen bonding is a portion, preventing the displacement of U36 (Figure 4).

NMR studies indicate that the SCH3 of ms2t6A weakens the ureido hydrogen-bonding interaction as compared to t6A (29). Taken together with quantum calculations indicating facile interconversion at room temperature between an ureido ring conformation and a non-ureido ring structure, the conformation around C10-N11 may change based upon the context of the ribosome. In the E.coli 30S ribosomal subunit, primarily responsible for decoding (80), A1493, A1492 and G530 have been found to recognize the first two base pairs of the anticodon–codon complex in the minor groove (76), triggering a closure of the small subunit around the tRNA (81). In this system, the ureido ring allows greater stacking interactions with A38 (19), possibly giving U36 more flexibility to adjust to the codon. The prokaryotic and eukaryotic ribosomes contain conserved regions of rRNA around the decoding site (82) and mRNA crosslinks to A1493 and A1492 corresponding residues in human 18S subunits (83). However, human tRNALys,3 may rotate the t6 moiety around the C10-N11 bond allowing the t6 side chain to hydrogen bond across the loop, thereby stabilizing U36 for correct positioning in the ribosome in a canonical fashion. Accuracy in tRNA–mRNA recognition is both dependent upon a kinetic proofreading mechanism and the stability of the cognate codon–anticodon complex (84). The role of modified bases in stablizing codon–anticodon complexes as opposed to destabilizing the loop conformation for an induced-fit in the ribosome seems to be base and system dependent (4).

SUPPLEMENTARY DATA
Supplementary Data are available at NAR Online.

Supplementary Material
[Supplementary Data]
 The authors would like to thank Philippe Bénas for discussions on the structure of human tRNALys,3. This work was supported by the Petroleum Research Fund (PRF 41701-GB7), the National Science Foundation (DUE-0431664 and CHE-0521063) and Truman State University. Funding to pay the Open Access publication charges for this article was provided by the Petroleum Research Fund.

Conflict of interest statement. None declared.

Figures and Tables
Figure 1 Secondary structure of human tRNALys,3 ASL. Chemical structures of modified bases found in ASLLys,3.

Figure 2 Plots of r.m.s.d. from the average structures. For each simulation, the average structure was calculated solely for the production run. The r.m.s.d. values have an arbitrary zero.

Figure 3 Stair-stepped conformation. (A) Stair-stepped conformation of the anticodon bases from the side (left) and from above (right) with dx, dy and dz parameters shown. (B) Snapshot of the wild-type simulation showing ms2t6A37 hydrogen bonding across the loop (left) and a snapshot of the ASLLys,3 ψ39 simulation when U36 is rotated into solution (right).

Figure 4 Conformations of ms2t6A. (A) with an ureido ring (B) without an ureido ring as calculated at the MPW1K/6-31+G (d,p) level of theory.

Table 1 The role of ms2t6A37

ASLLys,3	% Stair-stepped	SASA (Å2)	
wild-type ASLLys,3—mcm5s2U34, ms2t6A37, ψ39	96	150 (10)	
wild type 2—mcm5s2U34, ms2t6A37, ψ39	46	140 (20)	
ASLLys,3—mcm5s2U34, ms2t6A37	76	160 (20)	
ASLLys,3—ms2t6A37, ψ39	99	180 (20)	
ASLLys,3—ms2t6A37	93	160 (10)	
ASLLys,3—mcm5s2U34, ψ39	0	230 (30)	
ASLLys,3—mcm5s2U34	7	230 (10)	
ASLLys,3—ψ39	0	250 (10)	
ASLLys,3—unmodified	0	250 (10)	
Percentage of the simulation anticodons exhibited the stair-stepped conformation (left). SASA for the U36 base (right). Modified bases included in each simulation are indicated. Absence of a modification indicates the unmodified base was substituted at that position.

Table 2 Base–base stacking interactions

ASLLys,3	U36– A37 (5)	A37(5)– A38(5)	A37(6)– A38(5)	A37(u)– A38(6)	
wild-type ASLLys,3(mcm5s2U34, ms2t6A37, ψ39)	34	16	85		
wild-type 2 (mcm5s2U34, ms2t6A37, ψ39)	21	25	33	67	
ASLLys,3—mcm5s2U34, ms2t6A37	3	85	31		
ASLLys,3—ms2t6A37, ψ39	18	33	24		
ASLLys,3—ms2t6A37	28	67	32		
ASLLys,3—mcm5s2U34, ψ39	1	21	83		
ASLLys,3—mcm5s2U34	0	0	0		
ASLLys,3—ψ39	0	16	96		
ASLLys,3—unmodified	0	25	95		
The percent of the simulation time, base stacking was present for interactions between bases U36, A37 and A38 are shown. COM distances between bases were determined for the aromatic rings indicated; 5 indicates the 5-membered ring of A, 6 indicates the 6-membered ring of A and u indicates the ureido ring of ms2t6A37. Modified bases included in each simulation are indicated. Absence of a modification indicates the unmodified base was substituted at that position.
==== Refs
REFERENCES
1 Limbach P.A. Crain P.F. Pomerantz S.C. McCloskey J.A.  Structures of posttranscriptionally modified nucleosides from RNA Biochimie 1995 77 135 138 7541251 
2 Motorin Y. Grosjean H.  Grosjean H. Benne R.  Appendix 1: chemical structures and classification of posttranscriptionally modified nucleosides in RNA Modification and Editing of RNA 1998 Washington, DC American Society for Microbiology Press 534 549 
3 Yokoyama S. Nishimura S.  Soll D. RajBhandary U.L.  Modified nucleosides and codon recognition tRNA: Structure, Biosynthesis and Function 1995 Washington, DC American Society for Microbiology 207 223 
4 Agris P.F.  Decoding the genome: a modified view Nucleic Acids Res. 2004 32 223 238 14715921 
5 Sakurai M. Ohtsuki T. Watanabe K.  Modification at position 9 with 1-methyladenosine is crucial for structure and function of nematode mitochondrial tRNAs lacking the entire T-arm Nucleic Acids Res. 2005 33 1653 1661 15781491 
6 Blaise M. Becker H.D. Keith G. Cambillau C. Lapointe J. Giege R. Kern D.  A minimalist glutamyl-tRNA synthetase dedicated to aminoacylation of the tRNAAsp QUC anticodon Nucleic Acids Res. 2004 32 2768 2775 15150343 
7 Isel C. Marquet R. Keith G. Ehresmann C. Ehresmann B.  Modified nucleotides of tRNA3lys modulate primer/template loop–loop interaction in the initiation complex of HIV-1 reverse transcription J. Biol. Chem. 1993 268 25269 25272 7503978 
8 Ashraf S.S. Sochacka E. Cain R. Guenther R. Malkiewicz A. Agris P.F.  Single atom modification (O→S) of tRNA confers ribosome binding RNA 1999 5 188 194 10024171 
9 Yarian C. Marszalek M. Sochacka E. Malkiewicz A. Guenther R. Miskiewicz A. Agris P.F.  Modified nucleoside dependent Watson–Crick and wobble codon binding by tRNAUUULys species Biochemistry 2000 39 13390 13395 11063576 
10 Phelps S.S. Jeinic O. Joseph S.  Universally conserved interactions between the ribosome and the anticodon stem–loop of A site tRNA important for translocation Mol. Cell 2002 10 799 807 12419224 
11 Phelps S.S. Malkiewicz A. Agris P.F. Joseph S.  Modified nucleotides in tRNALys and tRNAVal are important for translocation J. Mol. Biol. 2004 338 439 444 15081802 
12 Kruger M. Pedersen S. Hagervall T.G.  The modification of the wobble base tRNAGlu modulates the translation rate of glutamic acid codons in vivo J. Mol. Biol. 1998 284 621 631 9826503 
13 Brierley I. Meredith M.R. Bloys A.J. Hagervall T.G.  Expression of a coronavirus ribosomal frameshift signal in Escherichia coli: influence of tRNA anticodon modification on frameshifting J. Mol. Biol. 1997 270 360 373 9237903 
14 Sprinzl M. Vassilenko K.S.  Compilation of tRNA sequences and sequences of tRNA genes Nucleic Acids Res. 2005 33 D139 D149 15608164 
15 Bjork G.R.  Söll D. RajBhandary U.L.  Biosynthesis and function of modified nucleosides tRNA: Structure, Biosynthesis, and Function 1995 Washington, DC American Society for Microbiology Press 165 206 
16 Agris P.F.  The importance of being modified: roles of modified nucleosides and Mg2+ in RNA structure and function Prog. Nucleic Acid Res. Mol. Biol. 1996 53 79 129 8650309 
17 Davis D.R. Veltri C.A. Nielsen L.  An RNA model system for investigation of pseudouridine stabilization of the codon–anticodon interaction in tRNALys, tRNAHis and tRNATyr J. Biomol. Struct. Dyn. 1998 15 1121 1132 9669557 
18 Durant P.C. Davis D.R.  Stabilization of the anticodon stem–loop of tRNALys,3 by an A+–C base pair and by pseudouridine J. Mol. Biol. 1999 285 115 131 9878393 
19 Murphy F.V. Ramakrishnan V. Malkiewicz A. Agris P.F.  The role of modifications in codon discrimination by tRNAUUULys Nature Struct. Mol. Biol 2004 11 1186 1191 15558052 
20 Chen Y. Sierzputowska-Gracz H. Guenther R. Everett K. Agris P.F.  5-Methylcytidine is required for cooperative binding of magnesium(2+) and a conformational transition at the anticodon stem–loop of yeast phenylalanine tRNA Biochemistry 1993 32 10249 10253 8399153 
21 Nobles K.N. Yarian C.S. Liu G. Guenther R.H. Agris P.F.  Highly conserved modified nucleosides influence Mg2+-dependent tRNA folding Nucleic Acids Res. 2002 30 4751 1760 12409466 
22 Cabello-Villegas J. Tworowska I. Nikonowicz E.P.  Metal ion stabilization of the U-turn of the A37 N6-dimethylallyl-modified anticodon stem–loop of Escherichia coli tRNAPhe Biochemistry 2004 43 55 66 14705931 
23 Von Ahsen U. Green R. Schroeder R. Noller H.F.  Identification of 2-hydroxyl groups required for interaction of a tRNA anticodon stem–loop region with the ribosome RNA 1997 3 49 56 8990398 
24 Lustig F. Elias P. Axberg T. Samuelsson T. Tittawella I. Lagerkvist U.  Codon reading and translational error J. Biol. Chem. 1981 256 2635 2643 6782093 
25 Yarian C. Townsend H. Czestkowski W. Sochacka E. Malkiewicz A. Guenther R. Miskiewicz A. Agris P.F.  Accurate translation of the genetic code depends on tRNA modified nucleosides J. Biol. Chem. 2002 227 16391 16395 11861649 
26 Durant P.C. Davis D.R.  The effects of pseudouridine and pH on the structure and dynamics of the anticodon stem–loop of tRNALys,3 Nucleic Acids Symp. Ser. 1997 36 56 57 9478205 
27 Stuart J.W. Gdaniec Z. Guenther R. Marszalek M. Sochacka E. Malkiewicz A. Agris P.F.  Functional anticodon architecture of human tRNALys3 includes disruption of intraloop hydrogen bonding by the naturally occurring amino acid modification, t6A Biochemistry 2000 39 13396 13404 11063577 
28 Sundaram M. Durant P.C. Davis D.R.  Hypermodified nucleosides in the anticodon of tRNALys stabilize a canonical U-turn structure Biochemistry 2000 39 12575 12584 11027137 
29 Durant P.C. Baji A.C. Sundaram M. Kumar R.K. Davis D.R.  Structural effects of hypermodified nucleosides in the Escherichia coli and human tRNALys anticodon loop: the effect of nucleosides s2U, mcm5U, mcm5s2U, mnm5s2U, t6A, and ms2t6A Biochemistry 2005 44 8078 8089 15924427 
30 Bénas P. Bec G. Keith G. Marquet R. Ehresmann C. Ehresmann B. Dumas P.  The crystal structure of HIV reverse-transcription primer tRNALys,3 shows a canonical anticodon loop RNA 2000 6 1347 1355 11073212 
31 Fuller W. Hodgson A.  Conformation of the anticodon loop in tRNA Nature 1967 215 817 821 6049728 
32 Grosjean H. Benne R.  Modification and Editing of RNA 1998 Washington, DC American Society for Microbiology Press 
33 Cornell W.D. Cieplak P. Bayly C.I. Gould I.R. Merz K.M. JrFerguson D.M. Spellmeyer D.C. Fox T. Caldwell J.W. Kollman P.A.  A second generation force field for the simulation of proteins, nucleic acids, and organic molecules J. Am. Chem. Soc. 1995 117 5179 5197 
34 Case D.A. Pearlman D.A. Caldwell J.W. Cheatham T.E. IIIRoss W.S. Simmerling C.L. Darden T.A. Mertz K.M. Stanton R.V. Cheng A.L.  1999 6.0 edn San Francisco, CA University of California 
35 Jorgensen W.L. Chandrasekhar J. Madura J.D.  Comparison of simple potential functions for simulating liquid water J. Chem. Phys. 1983 79 926 935 
36 York D.M. Darden T.A. Pedersen L.G.  The effect of long-range electrostatic interactions in simulations of macromolecular crystals: a comparison of the Ewald and truncated list methods J. Chem. Phys. 1993 99 8345 8348 
37 Ryckaert J.P. Ciccotti G. Berendsen H.J.C.  Numerical integration of the Cartesian equations of motion of a systems with constraints: molecular dynamics of n-alkanes J. Comput. Phys. 1977 23 327 341 
38 Berendsen H.J.C. Postma J.P.M. van Gunsteren W.F. A. DiNola A. Haak J.R.  Molecular dynamics with coupling to an external bath J. Comput. Phys. 1984 81 3684 3690 
39 Cheatham T.E. IIIKollman P.A.  Molecular dynamics simulations highlight the structural differences among DNA:DNA, RNA:RNA, and DNA:RNA hybrid duplexes J. Am. Chem. Soc. 1997 119 4805 4825 
40 Auffinger P. Bielecki L. Westhof E.  Symmetric K+ and Mg2+ ion-binding sites in the 5S rRNA loop E inferred from molecular dynamics simulations J. Mol. Biol. 2004 335 555 571 14672663 
41 Auffinger P. Westhof E.  RNA hydration: three nanoseconds of multiple molecular dynamics simulations of the solvated tRANAsp anticodon hairpin J. Mol. Biol. 1997 269 326 341 9199403 
42 Dixit S.B. Beveridge D.L. Case D.A. Cheatham T.E. IIIGiudice E. Lankas F. Lavery R. Maddocks J.H. Osman R. Sklenar H.  Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps Biophys. J. 2005 89 3721 3740 16169978 
43 Mokdad A. Krasovska M.V. Sponer J. Leontis N.B.  Structural and evolutionary classification of G/U wobble basepairs in the ribosome Nucleic Acids Res. 2006 34 1326 1341 16522645 
44 Nagan M.C. Kerimo S.S. Musier-Forsyth K. Cramer C.J.  Wild-type RNA microhelixAla and 3:70 variants: molecular dynamics analysis of local helical structure and tightly bound water J. Am. Chem. Soc. 1999 121 7310 7317 
45 Spackova N. Sponer J.  Molecular dynamics simulations of sarcin–ricin rRNA motif Nucleic Acids Res. 2006 34 697 708 16456030 
46 Wang J. Wolf R.M. Caldwell J.W. Kollman P.A. Case D.A.  Development and testing of a general amber force field J. Comput. Chem. 2003 25 1157 1174 15116359 
47 Hehre W.J. Radom L. Schleyer P.v.R. Pople J.A.  Ab Initio Molecular Orbital Theory. 1986 Wiley, NY 
48 Ditchfield R. Hehre W.J. Pople J.A.  Self-consistent molecular-orbital methods. IX. An extended Gaussian-type basis for molecular-orbital studies of organic molecules J. Chem. Phys. 1971 54 724 728 
49 Hariharan P.C. Pople J.A.  The effects of d-functions on molecular orbital energies for hydrocarbons Chem. Phys. Lett. 1972 66 217 219 
50 Hehre W.J. Ditchfield R. Pople J.A.  Self-consistent molecular orbital methods. XII. Further extensions of Gaussian-type basis sets for use in molecular orbital studies of organic molecules J. Am. Chem. Soc. 1972 56 2257 2261 
51 Bayly C.I. Cieplak P. Cornell W.D. Kollman P.A.  A well-behaved electrostatic potential based method using charge restraints for determining atom-centered charges: the RESP model J. Phys. Chem. 1993 97 10269 10280 
52 Frisch M.J. Trucks G.W. Schlegel H.B. Scuseria G.E. Robb M.A. Cheeseman J.R. Montgomery J.J.A. Vreven T. Kudin K.N. Burant J.C.  2003 Revision B.04 ed Wallingsford, CT Gaussian, Inc 
53 Lynch B.J. Fast P.L. Harris M. Truhlar D.G.  Adiabatic connection for kinetics J. Phys. Chem. A 2000 104 4811 4815 
54 Lynch B.J. Zhao Y. Truhlar D.G.  Effectiveness of diffuse basis functions for calculating relative energies by density functional theory J. Phys. Chem. A 2003 107 1384 1388 
55 Hariharan P.C. Pople J.A.  Influence of polarization functions on MO hydrogenation energies Theor. Chim. Acta 1973 28 213 222 
56 Francl M.M. Pietro W.J. Hehre W.J. Binkley J.S. Gordon M.S. DeFrees D.J. Pople J.A.  Self-consistent molecular orbital methods. XXIII. A polarization-type basis set for second-row elements J. Chem. Phys. 1982 77 3654 3665 
57 Lu X. Olson W.K.  3DNA: a software package for the analysis, rebuilding and visualization of three-dimensional nucleic acid structures Nucleic Acids Res. 2003 31 5108 5121 12930962 
58 Humphrey W. Dalke A. Schulten K.  VMD: visual molecular dynamics J. Mol. Graph 1996 14 33 38 8744570 
59 Pettersen E.F. Goddard T.D. Huang C.C. Couch G.S. Greenblatt D.M. Meng E.C. Ferrin T.E.  UCSF Chimera—a visualization system for exploratory research and analysis J. Comput. Chem. 2004 25 1605 1612 15264254 
60 Connolly M.L.  Analytical molecular surface calculation J. Appl. Cryst. 1983 16 548 558 
61 Case D.A. Darden T.A. Cheatham T.E. IIISimmerling C.L. Wang J. Duke R.E. Luo R. Merz K.M. Wang B. Pearlman D.A.  2004 La Jolla, CA Scripps Research Institute 
62 Lavery R. Zakrzewska K.  Neidle S.  Base and base pair morphologies, helical parameters, and definitions Oxford Handbook of Nucleic Acid Structure 1999 NY Oxford University Press 39 76 
63 Olson W.K. Dasika R.D.  Spatial configuration of ordered polynucleotide chains. 3. Polycyclonucleotides J. Am. Chem. Soc. 1976 98 5371 5380 956562 
64 Nagan M.C. Beuning P. Musier-Forsyth K. Cramer C.J.  Importance of discriminator base stacking interactions: molecular dynamics analysis of A73 microhelixAla variants Nucleic Acids Res. 2000 28 2527 2534 10871402 
65 Parthasarathy R. Ohrt J.M. Chheda G.B.  Modified nucleosides and conformation of anticodon loops: crystal structure of t6A and g6A Biochemistry 1977 16 4999 5008 911810 
66 Cabello-Villegas J. Winkler M.E. Nikonowicz E.P.  Solution conformations of unmodified and A37N6-dimethylallyl modified anticodon stem–loops of Escherichia coli tRNAPhe J. Mol. Biol. 2002 319 1015 1034 12079344 
67 Auffinger P. Westhof E.  Grosjean H. Benne R.  Effects of pseudouridylation Modification and Editing of RNA 1998 Washington, DC ASM Press 103 112 
68 Westhof E. Dumas P. Moras D.  Hydration of transfer RNA molecules: a crystallographic study Biochimie 1988 70 145 165 3134940 
69 Griffey R.H. Davis D.R. Yamaizumi Z. Nishimura S. Bax A. Hawkins B. Poulter C.D.  15N-labeled Escherichia coli tRNAfMet, tRNAGlu, tRNATyr, and tRNAPhe J. Biol. Chem. 1985 260 9734 9741 3894360 
70 Davis D.R. Poulter C.D.  1H-15N NMR studies of Escherichia coli tRNAPhe from hisT mutants: a structural role for pseudouridine Biochemistry 1991 30 4223 4231 2021615 
71 Zacharias M. Skelnar H.  Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model J. Mol. Biol. 1999 289 261 275 10366504 
72 Ladner J.E. Jack A. Robertus J.D. Brown R.S. Rhodes D. Clark B.F.C. Klug A.  A structure of yeast phenylalanine transfer RNA at 2.5 Å resolution Proc. Natl Acad. Sci. USA 1975 72 4414 4418 1105583 
73 Westhof E. Dumas P. Moras D.  Crystallographic refinement of yeast aspartic acid transfer RNA J. Mol. Biol. 1985 194 119 145 3897553 
74 Moras D. Dock A.C. Dumas P. Westhof E. Romby P. Ebel J.P. Giege R.  Anticodon-anticodon interaction induces conformational changes in tRNA: yeast tRNAAsp, a model for tRNA–mRNA recognition Proc. Natl Acad. Sci. USA 1986 83 932 936 3513167 
75 Basavappa R. Sigler P.B.  The 3 Å crystal structure of yeast initiator tRNA: functional implications in initiator/elongator discrimination EMBO J. 1991 10 2105 3111 2065656 
76 Ogle J.M. Brodersen D.E. Clemons W.M. Tarry M.J. Carter A.P. Ramakrishnan V.  Recognition of cognate transfer RNA by the 30S ribosomal unit Science 2001 292 897 902 11340196 
77 Guckian K.M. Ren R.X.-F. Chaudhuri N.C. Tahmassebi D.T. Kool E.T.  Factors contributing to aromatic stacking in water: evaluation in the context of DNA J. Am. Chem. Soc. 2000 122 2213 2222 
78 Kool E.T.  Replacing the nucleobases in DNA with designer molecules Acc. Chem. Res. 2002 35 936 943 12437318 
79 Fredericks S.Y. Jorda K.D. Zwier T.S.  Theoretical characterization of the structures and vibrational spectra of benzene-(H2O)n (n = 1–3) clusters J. Phys. Chem. 1996 100 7810 7821 
80 Carter A.P. Clemons W.M. JrBrodersen D.E. Morgan-Warren R.J. Wimberly B.T. Ramakrishnan V.  Functional insights from the structure of the 30S ribosomal subunit and its interactions with antibiotics Nature 2000 407 340 348 11014183 
81 Ogle J.M. Murphy F.V. IV, Tarry M.J. Ramakrishnan V.  Selection of tRNA by the ribosome requires a transition from an open to a closed form Cell 2002 111 721 732 12464183 
82 Spahn C.M.T. Beckmann R. Eswar N. Penczek P.A. Sali A. Blobel G. Frank J.  Structure of the 80S ribosome from Saccharomyces cerevisiae–tRNA–ribosome and subunit–subunit interactions Cell 2001 107 373 386 11701127 
83 Graifer D. Molotkov M. Styazhkina V. Demeshkina N. Bulygin K. Eremina A. Ivanov A. Laletina E. Ven'Yaminova A. Karpova G.  Variable and conserved elements of human ribosomes surrounding the mRNA at the decoding and upstream sites Nucleic Acids Res. 2004 32 3282 3293 15208366 
84 Rodnina M.V. Wintermeyer W.  Ribosome fidelity: tRNA discrimination, proofreading and induced fit TIBS 2001 26 124 130 11166571

