computational prediction of the sequence dependent structure , dynamics , and thermodynamics of rna is important for extracting information from genome sequences .
this is particularly important for rna because it has many biological functions,(1 ) including catalysis , control of gene expression , roles in enzymatic activities,(6 ) and perhaps affecting protein folding.(7 ) yet , the database of three - dimensional structures of rna is limited in comparison to that of proteins .
experimental methods such as x - ray crystallography and nmr are unable to keep pace with the rapid discovery of important rna sequences .
even when experimental data are available , computational methods employing force fields are used for refinement of structure . due to computational expense ,
these methods include residue - centered force fields ( coarse - grained(8 ) ) , atom - centered force fields ( amber,(9 ) charmm , and gromos ) , approximate quantum mechanics , and mixed quantum - mechanics / molecular mechanics ( qm / mm ) .
molecular mechanics force fields are inexpensive enough that molecular dynamics ( md ) simulations on the microsecond and beyond time scale are now possible , enabling the comparison of structural and thermodynamic properties of nucleic acids with experiment .
for some systems such as rna with tandem ga base pairs and unnatural ribonucleotides , and quadruplex dna with single stranded loops , the current , widely used amber99 force field does not reproduce experimental results .
it was originally parametrized in 1995 on the basis of qm calculations on small model systems because computer power was limited.(9 ) moreover , the biological importance of rna was only beginning to be discovered so that model systems for qm computations focused on proteins and dna .
, we test the effect of the  reparameterization as substituted into the amber99 force field to give amber99.(44 ) briefly , the parametrization for  used qm calculations on each of the complete nucleosides , adenosine , cytidine , guanosine , and uridine .
this was prompted by the observation that md calculations for single stranded tetramers produced an unlikely fraction of syn conformers .
thus , for each nucleoside , 288 data points were fit to two dihedrals with four cosine terms.(44 ) comparison between nmr spectra and md simulations for cytidine and uridine revealed that amber99 improved predictions of the fractions of nucleosides having c3-endo sugar puckers and anti orientations of bases.(44 ) recently , the effect of amber99 and other  revisions on structural predictions for uucg and gnra rna hairpins were analyzed.(45 ) the  revisions were crucial to maintain the structural characteristics of these hairpins.(45 ) experimental benchmarks can facilitate the improvement of force fields .
for example , interpretation of nmr spectra of short single - stranded oligoribonucleotides can be complicated due to multiple conformations and dynamics , but can reveal structural propensities . in particular , the nuclear overhauser effect ( noe ) differs by roughly 20-fold between syn and anti conformations of a base relative to its ribose , the ribose pucker can be determined from j12 coupling constants , and various noes report on base  base stacking .
base stacking , especially for the middle nearest neighbor of a - c.(50 ) similarly , single - stranded gaaac has adenine  adenine stacking both as dna and rna.(51 ) nmr spectra of caau , caau , and caau , where m and l stand for 2o - methyl and locked nucleic acid ( lna ) sugars ( i.e. , ribose with a methylene bridge between the 2 oxygen and 4 carbon ) , respectively , revealed helical preorganization , especially for oligonucleotides with lna residues.(52 ) short single strands provide a stringent test of rna force fields for torsions because the conformational space is not limited by hydrogen bonding between bases and is minimally limited by steric exclusions . moreover ,
structural transitions in single stranded rna typically occur in 100 ns so that comparisons can be made to md calculations on accessible time scales . in this paper ,
h nuclear overhauser effect spectroscopy ( noesy ) , h double - quantum - filtered correlation spectroscopy ( dqf - cosy ) , total correlation spectroscopy ( tocsy ) , and h  p correlation spectroscopy ( hetcor ) nmr spectra reveal that single - stranded gacc rna ( figure 1 ) has a preference for a - form - like conformations as defined by c3-endo sugar puckers , anti orientation of bases relative to sugars and base 
these results are compared to predictions from md simulations with amber99(9 ) and amber99(44 ) force fields .
thus , amber99 should improve modeling of the structures and dynamics of unpaired regions and the thermodynamics of folding an rna from an unpaired single strand .
( left ) single - stranded gacc rna and atom notations used ( right , top to bottom ) in guanine , adenine , cytosine , and ribose and phosphate . for g / a , and c ,  dihedral angle is defined by o4-c1-n9-c4 and o4-c1-n1-c2 , respectively .
the nmr sample had 3.2 mm gacc in 80 mm nacl , 10 mm sodium phosphate , and 0.5 mm disodium edta at ph 7.0 .
two deuterium exchanges were performed on the sample with 99.9% d2o ( cambridge isotopes laboratories ) .
p hetcor , dqf - cosy , and tocsy at 2 and 20 c ( table s1 ) .
noesy spectra were recorded with mixing times ( m ) of 200 and 800 ms .
spin coupling constants were used to deduce the sugar puckers for each residue ( table s3 ) .
proton distance restraints were calculated from magnitudes of noes in 200 ms mixing time noesy spectra at 2 c ( table s4 ) using the two - spin approximation.(57 ) restraint limits were calculated by allowing for a  4-fold range in each noe volume to account for relaxation effects , spin diffusion , and/or baseline distortions .
the magnitudes of noes are proportional to 1/(rij ) , where rij is the distance between protons i and j. h1 to h2 noes were used as reference noes in the distance calculations with h1 to h2 distances taken as 2.75 and 3.00  , respectively , for c3-endo and c2-endo conformations ( see supporting information for detail ) .
the upper limit of some of the distance restraints was set to 10.0  in order to search a broader conformational space even though the magnitudes of these noe peaks are medium in size indicating a distance < 5  .
similarly , the lower and upper limits of some distance restraints were set to 4.5 and 99.0  , respectively , when no noe was seen ( table s4 ) .
simulated annealing md simulations with nmr restraints were done with sander(58 ) using the amber99 force field .
the first md simulation was initiated from a starting structure of gacc in a - form .
another 199 simulations were performed with the starting structure for each taken from the final structure of the previous simulation .
the generalized born implicit solvent model with a 1 m salt concentration was used in the simulated annealing protocol . for each simulation ,
a total of 50 ps of md were run with a 1 fs time step .
berendsen temperature control(62 ) with a coupling time of 1 ps was used in all simulated annealing md simulations .
distance and sugar pucker restraints extracted from noe data were applied using square bottom wells with parabolic sides with force constants of 50.0 kcal / mol  and 50.0 kcal / mol rad , respectively ( table s5 ) .
chirality restraints with 50.0 kcal / mol rad were used to keep chiral atoms in the correct orientations so that the covalent bonds do not flip at high temperatures in the simulated annealing protocol ( table s5 ) .
temperature was increased to 3000 k in 5 ps and gradually cooled to 0 k in 45 ps . during the first 5 ps
, relative weights of the nmr restraints were set to 0.1 to make sure that the structure was independent of the starting conformation ( see supporting information for structural parameters as a function of time during simulation ) .
the nmr restraint weight was set to 1.0 during the rest of the simulation ( table s6 ) .
a - form single - stranded gacc rna was created with the nucgen module of amber9.(58 ) gacc was solvated with tip3p water molecules(63 ) in a truncated octahedral box with 1076 water molecules .
the system was minimized in two steps : ( 1 ) with gacc held fixed with a restraint force of 500 kcal / mol  , steepest descent minimization of 500 steps was followed by a conjugate gradient minimization of 500 steps in order to minimize the potential energy due to water molecules .
( 2 ) with all restraints removed , steepest descent minimization of 1000 steps was followed by a conjugate gradient minimization of 1500 steps .
after minimization , two steps of pressure equilibration were done with the sander module in amber9 : ( 1 ) gacc was held fixed with a restraint force of 10 kcal / mol -  .
constant volume dynamics with an atom - based cutoff of 10.0  was used for the direct space sum .
a total of 20 ps of md were run with a 2 fs time step .
( 2 ) the above conditions were chosen , except that constant pressure dynamics with isotropic position scaling was turned on .
a total of 100 ps of md were run with a 2 fs time step .
each production run was similar to the second step of the pressure equilibration described above .
a total of 50 ns of md was run at 300 k with a 1 fs time step . in production runs , simulations were carried out with the pmemd module in amber9.(58 ) trajectory files were written at each 500 fs time step .
different random seed numbers were used for each simulation . in order to examine the time evolution of gacc with the amber99 force field
, a single simulation of 1.9 s was also run at 300 k. this simulation took over six months on a dell precision t5400 workstation using the eight 2.33 mhz intel xeon dual quad cores , and was stopped after no new stable structures appeared during the last 800 ns .
root - mean - square deviations of whole structure and backbone from a - form gacc single strand as well as  and  dihedral angles of each residue were analyzed . measured rate constants for stacking and unstacking of nucleotides in ac , ca , cc , a7 , poly a , and poly c range from 5  10 to 2  10 s , indicating that structural transitions occur in roughly 100 ns .
coupling constants from 1d h nmr spectra of gacc revealed that from 2 to 55 c all the residues prefer c3-endo sugar pucker , especially the middle nucleotides , a2 and c3 ( tables 1 and s3 and figure s1 ) .
figure 2 shows the noesy walk for gacc in an 800 ms mixing time noesy spectrum . in the 200 ms noesy spectrum ,
the noes of g1h1-h8 , a2h1-h8 , c3h1-h6 , and c4h1-h6 are all weak , implying anti base orientation for the residues ( table 1 and figure s2 ) . noes of g1h2/h3-a2h8 , a2h2/h3-c3h6 , and c3h2/h3-c4h6 are consistent with base stacking in a - form structures ( figure 3 ) .
weak c3h6-c4h2 and c3h3-c4h2 cross - peaks are seen in 200 and 800 ms noesy spectra ( figures 3 and s3 ) , however , which suggests a population of conformations with a flipped , upside down c4 sugar
thus , no c3h6-c4h2 and c3h3-c4h2 distance restraints were used in the simulated annealing protocol .
the chemical shifts of all phosphorus resonances were within 0.2 ppm indicating no unusual backbone conformations ( table s1 ) .
noesy walk of gacc from 800 ms mixing time noesy experiment ( 2 c ) .
200 ms 2d noesy of gacc at 2 c showing cross - peaks of the sugar and h6/h8 protons .
note that intense peaks are seen between the inter - residue protons , g1h2/h3-a2h8 , a2h2/h3-c3h6 , and c3h2/h3-c4h6 , implying base 
cross - peak of c3h6-c4h2 is due to minor conformation ( see also figure s3 ) .
values for amber99 and amber99 not in parentheses are on the basis of fifteen individual 50 ns simulations with errors corresponding to sample standard deviations .
the percentages are listed for different time domains in the simulation in the following order : 1700 ns/7011100 ns/11011900 ns/11900 ns .
fractions of sugar pucker in c3endo conformation are calculated from plots such as figure s5 . for nmr values ,
the same thing is done for % anti using plots such as figure s6 .  angles of 200 , 300 , and 60 correspond to anti , high - anti , and syn conformations , respectively .
intranucleoside h1 to h8/h6 distances for anti and syn conformations are typically 3.5  and 2.1  , respectively , so that the noe volumes for anti and syn differ by 20-fold .
values labeled  are standard deviations , but are not due to random error . because gacc is dynamic , a perfect force field would still give different percentages for different 50 ns simulations .
no standard deviations are given for the 1.9 s simulation because it shows different structural behavior in the 1700 ns , 7011100 ns , and 11011900 ns time domains .
these values include both anti and high - anti conformations ( see figure s6 ) .
most of the syn conformations occur after 1.1 s when gacc is not in a stable conformation ( figure 4 ) . prior to 1.1 s , g1 and a2
the percentage of structures consistent with nmr spectra was calculated by integrating the number of structures with distances between 1.50  and the upper limit of the noe distance restraint and dividing by the total number of structures generated .
this reflects the fact that nmr noes depend on ( 1/rij ) as weighted by the fraction of population with distance , rij .
time evolution ( in ns ) of gacc single - strand with amber99 force field .
top two plots show rms deviations ( in  ) of whole structure and backbone from a - form gacc .
 dihedral angles around 180 , 300 , and 90 represent anti , high - anti , and syn base orientations , respectively .
 dihedral angles around 90 and 120 represent sugar puckering of c3-endo and c2-endo , respectively .
red and green colors in top plot show sampling of major and minor conformations , respectively , while blue color shows sampling of a non - a - form conformation ( see figure s4 ) .
for each of the 50 ns unrestrained md simulations , snapshots were extracted from the trajectories at intervals of 5 ns and analyzed in order to cluster them into structural groups .
no preferred conformation was found for the fifteen 50 ns amber99 md simulations , and g1 of gacc went to a syn base orientation within 0.5 ns on average .
coupled with the force field s preference for c2-endo sugar pucker,(44 ) this creates distortions in the backbone and the other residues . as a result , random - coil conformations are observed within less than 10 ns with the amber99 force field . while simulations with amber99 generated rapid and large structural fluctuations of gacc within 50 ns , simulations with amber99 did not . for the fifteen 50 ns amber99 md simulations , 107 and 11 out of 150 structures
cluster into a - form - like major and minor conformations , respectively . to further explore the time dependence of dynamics with amber99 ,
time evolutions of rms deviations of all atoms and backbone heavy atoms with respect to a pure a - form gacc conformation and of the  and  torsion angles of each residue are shown in figure 4 .
this decomposition provides insight into the individual conformational transformations of the structural dynamics of gacc .
cluster analysis was done on this single 1.9 s md simulation and revealed that 63 and 18 out of 380 structures represent major and minor conformations , respectively , while 77 structures represent a stable non - a - form conformation ( see figures 4 and s4 ) , which is not consistent with nmr spectra .
the rest of the structures in the 1.9 s md simulation do not cluster to any conformation and are in non - a - form conformations .
the results show that for 700 ns amber99 generates structures consistent with nmr but that additional revisions are necessary . in particular
, gacc did not return to an a - form - like structure during the 800 ns between 1.1 and 1.9 s when the predicted structures were random coil .
structures extracted every 0.5 ps from the unrestrained md trajectories were analyzed to quantify the sugar pucker and base orientation of gacc residues predicted by each force field .
population distribution plots of pseudorotation phase angles and  torsions for each residue were created ( figures s5 and s6 ) to obtain the percentages of c3-endo sugar pucker and anti type base orientation ( table 1 ) .
the sugar pucker of a five - membered ring can be described by the pseudorotation phase angle . in rna and dna , pseudorotation phase angles of 18 and 162 that correspond to c3-endo and c2-endo sugar pucker , respectively
these percentages can be quantitatively compared to nmr results ( table 1 ) . at 25 c ,
nmr indicates around 65% , 85% , 85% , and 75% c3-endo sugar pucker for g1 , a2 , c3 , and c4 , respectively . for a series of fifteen 50 ns simulations that start in a - form ,
the amber99 force field predicts 49% , 33% , 64% , and 41% c3-endo sugar pucker for g1 , a2 , c3 , and c4 , respectively ( table 1 and figure s5 ) .
the 1.9 s amber99 simulation predicts 76% , 81% , 83% , and 67% c3-endo sugar pucker for g1 , a2 , c3 , and c4 , respectively , which agrees reasonably well with the nmr data ( table 1 ) .
the predicted percentages of c3-endo sugar pucker , however , depend on the time domain within the 1.9 s simulation ( table 1 ) .
for example , they are 90% , 85% , 92% , and 74% for g1 , a2 , c3 , and c4 , respectively , in the first 700 ns , where gacc samples a - form - like and random coil conformations ( figure 4 ) . in contrast , they are 52% , 69% , 67% , and 65% in the 11011900 ns region where md predicts random coil .
quantitative comparisons can also be made between predicted and nmr determined orientation of bases relative to sugars .
for a series of fifteen 50 ns simulations that start in a - form , the amber99 force field predicts gacc to have 37% , 97% , 90% , and 92% anti / high - anti conformations for g1 , a2 , c3 , and c4 , respectively .
the anti range of the  torsions observed in amber99 simulations was 180360 and 180300 for g1/a2 and c3/c4 , respectively ( see figure s6 ) .
the broad range for amber99 is due to generation of both anti and high - anti conformations .
the 1.9 s amber99 simulation predicts 53% , 57% , 93% , and 90% , anti conformations for g1 , a2 , c3 , and c4 , respectively , with only 26% of high - anti conformations for these residues .
most of the syn conformations , however , occur after 1.1 s when gacc is not in a stable structure ( table 1 and figure 4 ) . prior to 1.1 s ,
g1 and a2 are 80% anti on average , with the rest either syn or high - anti .
for each force field , the structures used in part 3.3 were also used to calculate the percentages of structures consistent with the nmr distance restraints for some of the h  h noes ( table 1 ) .
these predictions were extracted from the population distribution plots shown in figures s7s9 and provide another measure of prediction accuracy . in 50 ns simulations starting with a - form and using the amber99 force field , the percentages of structures consistent with nmr distance restraints for g1h2-a2h8 , a2h2-c3h6 , a2h3-c3h6 , a2h2-c3h1 , c3h2-c4h6 , and c3h3-c4h6 are 19% , 24% , 49% , 9% , 41% , and 65% , respectively . with the 1.9 s amber99
simulation , the fractions are 31% , 54% , 60% , 63% , 35% , and 72% , respectively ( table 1 ) .
the noe cross - peak a2h2-c3h1 represents stacking of the middle bases a2 and c3 , and this cross - peak is underestimated much more by amber99 than by the amber99 force field , even when the amber99 simulation is run to 1.9 s .
as expected , the percentages of structures in the 1.9 s simulation that are consistent with nmr distance restraints depend on the time domain ( table 1 ) . for example ,
between 701 and 1100 ns , the g1h2-a2h8 and c3h2-c4h6 distances never fall in the range consistent with nmr data . as discussed below
, the stable structure generated between 701 and 1100 ns is inconsistent with nmr spectra .
comparison of nmr measured and md predicted noe  distances  ( r ) is a more direct test of the force fields . these  distances  are not true distances if gacc populates more than one conformation .
the structures used in part 3.3 were used to calculate the noes ( 1/rij ) for each snapshot , which then yielded average noe and , therefore , average  distance 
, r. that is , the  distance  is an ensemble average with each exact distance weighted by its inverse sixth power . for the combined fifteen 50 ns md simulations , rms deviations of six noe  distances  between nmr and , respectively , amber99 and amber99 are 0.84 and 0.39  . for the single 1.9 s md simulation of amber99 ,
r is the  distance  calculated from measured or predicted noe ( 1/r ) .
if more than one conformation is present , then r does not represent any physical distance .
the exp , a99 and a99 subscripts stand for experiment , amber99 , and amber99 .
predictions are from single 1.9 s md simulation of amber99 and for combined fifteen 50 ns md of both force fields ( 15  50 ns ) using structures extracted every 0.5 ps ( see table s7 for detail ) .
 distances  in parentheses are for the following time domains of the 1.9 s simulation : 1700 ns/7011100 ns/11011900 ns .
 distances  , r , for amber99 simulations of 50 ns starting from a - form are consistently shorter than experimental  distances  from nmr , presumably because gacc is not entirely a - form in solution , but remains a - form in 50 ns simulations ( see comparisons of noes in table s11 ) . the shorter  distances  , however , fall within the range consistent with noe distance restraints ( table 1 ) .
the 1.9 s simulation allows more time for sampling structures , so that the systematic effect of short simulation time is reduced ( see figure s16 ) .
again as expected , the distance comparison for the 1.9 s simulation depends on the time domain ( table 2 ) .
the distances predicted from the first 700 ns agree better with the nmr than those predicted from the 701 to 1100 ns and 1101 to 1900 ns domains .
figure 4 shows the time evolution for the 1.9 s md simulation . in figure 4 , regions highlighted with red and green represent sampling of major and minor a - form - like conformations , respectively , whereas regions highlighted with blue represent a stable non - a - form conformation ( figure s4 and table s8 ) .
the rest of the regions highlighted with black represent sampling of random - coil non - a - form conformations . during the first 700 ns md
, there are multiple major to minor and a - form - like to non - a - form conformational transitions ( figure 4 ) .
the switch from major to minor conformation is due to the conformational change of the backbone ( figure 4 ) . during these periods ,  dihedrals of g1 , a2 , c3 , and c4 are mostly around 180 , which is the anti region .
g1 goes to the syn conformation a few times during these periods , but returns to anti conformation almost instantly ( figure 4 ) .
between 200 and 375 ns , 550600 ns , and around 520 ns , however , the structure resembles neither major nor minor conformations . at the beginning of these times ,
base orientation of g1 goes to either high - anti or syn , which distorts the backbone .
this allows a2 to change its base orientation to either high - anti or syn ( figure 4 ) . during these times ,
c3 and c4 mostly stay in anti base orientation , but g1 , a2 , and c3 start sampling c2-endo sugar puckering ( figure 4 ) . during the whole simulation time
between 701 and 1100 ns , a stable and unique non - a - form conformation is sampled by amber99 ( figures 4 and s4 ) . at around 700 ns ,
the base orientation and sugar puckering of all four residues are anti and c3-endo , respectively , but the backbone changes dramatically and creates a conformation with g1 between c3 and c4 ( figure s4 ) .
this stable conformation has unique potential noe cross - peaks such as g1h8-c3/c4h2/h3 , g1h8-c3/c4h5/h6 , and c3h6-g1h2/h3 due to proton  proton distances < 5  .
these potential cross - peaks would be easy to detect because they would appear in unique and uncrowded noesy regions but no such cross - peaks have been detected . after 1100 ns , the structure transforms to random coil conformations .
again , first g1 goes to a syn base orientation that gives enough flexibility to the backbone to search other possible configurations .
c3 and c4 rarely change their base orientations but sample syn a couple of times ( figure 4 ) .
similar to the random - coil conformations seen within 700 ns , fluctuations of the  torsions are high after 1100 ns where all the residues start to sample c2-endo sugar puckering significantly ( figure 4 ) .
the terminal g1 and c4 nucleotides are flexible enough to change their conformations from anti to syn and/or c3-endo to c2-endo .
these changes are seen in figure 4 , but are not the main reason that the structure samples random - coil conformations ; this broader sampling is due primarily to conformational changes of a2 .
c3 and c4 , however , are found primarily in anti conformations as expected for cytidine,(44 ) even when the structure is in a random - coil conformation .
one aspect of figure 4 is that in order for the single - stranded gacc to have a stable conformation , it has to have all the bases in anti orientations ( figure 4 ) . otherwise , high fluctuations on the  and  torsions yield random - coil conformations .
gacc presumably has an ensemble of structures in rapid exchange that contribute to the nmr spectrum .
simulated annealing with nmr restraints was used to see if any particular structures are consistent with the nmr data . the protocol described in the methods section generated 182 out of 200 structures with constraint energies less than 0.1 kcal / mol . of these structures ,
cluster analyses showed that 116 and 32 represent a - form - like major and first minor conformations , respectively ( figures 5 and s4 and tables s9 and s10 ) .
another 18 structures represent a second a - form - like minor conformation with an average energy less favorable than both the major and first minor conformations ( figure s4 ) .
both minor conformations have a flipped , upside down c4 sugar , but with different backbone orientations ( figure s4 ) .
the flipped c4 sugar conformation is consistent with weak c3h6c4h2 and c3h3-c4h2 noes that were observed , but not used as restraints .
the major and first minor conformations have stacking between g1/a2 , a2/c3 , and c3/c4 , while the c3/c4 stacking is missing in the second minor conformation . with the exception of the flipped c4 ,
overlap of the nmr modeled ( black ) and amber99 force field predicted ( red ) conformations of major ( left ) and minor ( right ) species . predicted structures are average from fifteen 50 ns unrestrained simulations starting from a - form .
heavy atom rmsd of major and minor conformations with respect to nmr modeled structures are 0.6 and 0.9  , respectively ( see text and figure s4 for detail ) .
major and minor conformations generated during 1.9 s simulation are similar to these ( see figure s4 ) .
in molecular mechanics simulations of rna , the torsional potential for the  torsion angle has a significant effect on the preferred orientations of the bases . for the amber force field , this potential
was originally parametrized in 1995.(9 ) advances in computers since 1995 allow more thorough qm computations to parametrize the  torsions for rna .
initial comparisons were made between amber99 and amber99 simulations and nmr spectra for cytidine and uridine , but comparisons for adenosine and guanosine were prevented due to aggregation at nmr concentrations.(44 ) nmr results for gacc provide a test for parameters for g and a , and for inclusion of parameters in simulations of an unpaired oligonucleotide with a ribose - phosphate backbone .
amber99 and other revised  parameters have also been tested against structures for uucg and gnra hairpins where the nucleotides are more constrained.(45 ) accurate simulations of unpaired regions of rna are important for predicting structures of loops and thermodynamics of folding .
the nmr results for gacc can be qualitatively compared to predictions of the amber99 and amber99 force fields . in the amber99 force field
unrestrained simulations , there are large fluctuations within 50 ns in rms deviations of individual residues from anti / c3-endo conformation and in backbone from a - form ( table 1 ) . evidently , anti / c3-endo conformations and a - form are not preferred by the amber99 force field and transition barriers allow rapid conformational changes ( see figures s10 and s11 ) . in contrast , all residues of gacc in the fifteen 50 ns amber99 force field simulations show stable conformations similar to anti / c3-endo . because few transitions are seen within 50 ns with amber99 , the 1.9 s simulation provides a better comparison . in that simulation , a stable backbone similar to a - form
after 700 ns , various random - coil and a non - a - form stable structure inconsistent with nmr spectra are seen .
the non - a - form conformation is likely an artifact of the force field , and suggests that reparameterization beyond  is necessary .
the high - dimensional character of the conformational space , however , does not allow definitive identification of the torsions causing this non - a - form conformation to be stable enough to be sampled significantly by amber99. histogram analyses of rna backbones in crystal structures show preferred regions for  ,  ,  ,  ,  , and  , and most dihedrals in the 1.9 s md occur in these regions ( figure s12 ) .
all the residues of the stable non - a - form conformation have  between 7 and 71 ( table s8 ) , whereas an a - form rna has  around 70 and in crystal structures ,  is rarely between 40 and + 40.(68 ) it is possible that a revision of the  torsion for rna might improve predictions . alternatively ,
the stable non - a - form and the random - coil conformations may actually exist in solution , but are not prevalent enough to influence the nmr spectra . testing this hypothesis , however
quantitative comparisons show that the sugar pucker , base orientation , and proton  proton noe distances derived from nmr results are predicted better by amber99 than amber99 ( tables 1 and 2 ) .
even weak c3h6- c4h2 and c3h3- c4h2 noes are consistent with the flipped c4 generated by unrestrained amber99 simulations . in order to test the effects of reparameterizing only  , the parmbsc0 force field,(46 ) which has revised parameters for  and  , was not used in md simulations .
preliminary results on a different benchmark , however , indicate that combining parmbsc0 with amber99 further improves agreement with experiment ( i.y .
a different reparameterization of  produced the ode force field.(43 ) although the ode force field has more parameters than amber99 , it predicts similar potential energy surfaces for nucleosides.(44 ) thus , it should provide predictions similar to those of amber99 , although tests on uucg and gnra tetraloop hairpins indicate that the ode force field overpredicts high - anti conformations while amber99 does not.(45 ) additional tests are required to see if this difference is general .
the structural fluctuations observed with amber99 occur in less than 10 ns , which is much faster than the 100 ns expected on the basis of experimental measurements on single - stranded oligo and polynucleotides . in contrast , amber99 generated structures consistent with nmr spectra are typically stable for 50100 ns ( figure 4 ) .
evidently , amber99 reasonably approximates barriers to conformational change . because the dynamics of unpaired oligonucleotides occur in less than 1 s , nmr reports a population weighted average of spectral properties for the ensemble of structures
. somewhat surprisingly , the nmr spectra for gacc are consistent with expectations for a high population of a - form - like conformations .
simulated annealing with amber99 and nmr restraints generates two a - form - like structures with the minor structure having an upside down terminal ribose that is consistent with minor noesy peaks .
this suggests that the ensemble has largely a - form conformations , rather than a high population of random conformations .
the first 700 ns of amber99 unrestrained md is largely consistent with this interpretation ( figure 4 ) . after 700 ns , however , a structure inconsistent with nmr spectra is generated , suggesting that further refinement of the force field will be important .
experimental nmr restraints , however , appear able to compensate for limitations of the current amber force field . in order to sample physical regions of conformational space in md simulations , the force field
such sampling is important for predicting the structure and dynamics of rna and the thermodynamics of folding an rna from an unpaired single strand .
modification of  parameters in the amber99 force field improved agreement between simulations and nmr spectra of single - stranded gacc rna , as summarized in figure 6 .
the modified  torsion parameters are therefore expected to improve sampling of non - watson 
crick regions and terminal base pairs in rna simulations . in watson  crick regions of an rna duplex , the  torsions are restricted by cross - strand hydrogen bonding , so the amber99 torsion parameters should have less effect .
a stable non - a - form and random - coil conformations are observed significantly during the 1.9 s md simulation , suggesting refinement of other parameters will be important . in principle , such refinement can rely on qm calculations and the parmbsc0 revisions(46 ) will likely help .
many subtle interactions , however , define rna.(42 ) thus , some combination of knowledge based force fields(74 ) with qm may prove useful .
finally , nmr spectra for single - stranded regions of rna provide useful benchmarks for testing any revisions .
 accurately modeling unpaired regions of rna is important for predicting structure , dynamics , and thermodynamics of folded rna . comparisons between nmr data and molecular dynamics simulations provide a test of force fields used for modeling . here , 
 nmr spectroscopy , including noesy , 1h31p hetcor , dqf - cosy , and tocsy , was used to determine conformational preferences for single - stranded gacc rna . 
 the spectra are consistent with a conformational ensemble containing major and minor a - form - like structures . in a series of 50 ns molecular dynamics ( md ) simulations with the amber99 force field in explicit solvent , initial a - form - like structures 
 rapidly evolve to disordered conformations . 
 a set of 50 ns simulations with revised  torsions ( amber99 force field ) gives two primary conformations , consistent with the nmr spectra . a single 1.9 s md simulation with the amber99 force field showed that the major and minor conformations are retained for almost 68% of the time in the first 700 ns , with multiple transformations from a - form to non - a - form conformations . 
 for the rest of the simulation , random - coil structures and a stable non - a - form conformation inconsistent with nmr spectra were seen . 
 evidently , the amber99 force field improves structural predictions for single - stranded gacc rna compared to the amber99 force field , but further force field improvements are needed . 
Introduction
Methods
Results
Discussion