The nature of stacking interations in polynucleotides. Molecular states in Oligo- and polyribocytidylic acids by relaxation analysis

Kinetic Mechanism of RNA Helix-Terminal Basepairing-A Kinetic Minima Network Analysis

RNA functions are often kinetically controlled. The folding kinetics of RNAs involves global structural changes and local nucleotide movement, such as base flipping. The most elementary step in RNA folding is the closing and opening of a basepair. By integrating molecular dynamics simulation, master equation, and kinetic Monte Carlo simulation, we investigate the kinetics mechanism of RNA helix-terminal basepairing. The study reveals a six-state folding scheme with three dominant folding pathways of tens, hundreds, and thousands of nanoseconds of folding timescales, respectively. The overall kinetics is rate limited by the detrapping of a misfolded state with the overall folding time of 10^-5 s. Moreover, the analysis examines the different roles of the various driving forces, such as the basepairing and stacking interactions and the ion binding/dissociation effects on structural changes. The results may provide useful insights for developing a basepair opening/closing rate model and further kinetics models of large RNAs.

Physics-based all-atom modeling of RNA energetics and structure

The database of RNA sequences is exploding, but knowledge of energetics, structures, and dynamics lags behind. All-atom computational methods, such as molecular dynamics, hold promise for closing this gap. New algorithms and faster computers have accelerated progress in improving the reliability and accuracy of predictions. Currently, the methods can facilitate refinement of experimentally determined nuclear magnetic resonance and x-ray structures, but are 'unreliable' for predictions based only on sequence. Much remains to be discovered, however, about the many molecular interactions driving RNA folding and the best way to approximate them quantitatively. The large number of parameters required means that a wide variety of experimental results will be required to benchmark force fields and different approaches. As computational methods become more reliable and accessible, they will be used by an increasing number of biologists, much as x-ray crystallography has expanded. Thus, many fundamental physical principles underlying the computational methods are described. This review presents a summary of the current state of molecular dynamics as applied to RNA. It is designed to be helpful to students, postdoctoral fellows, and faculty who are considering or starting computational studies of RNA. WIREs RNA 2017, 8:e1422. doi: 10.1002/wrna.1422.

Kinetics of the B-A transition of DNA: analysis of potential contributions to a reaction barrier

Because of open problems in the relation between results obtained by relaxation experiments and molecular dynamics simulations on the B-A transition of DNA, relaxation measurements of the B-A dynamics have been extended to a wider range of conditions. Field-induced reaction effects are measured selectively by the magic angle technique using a novel cell construction preventing perturbations from cell window anisotropy. The kinetics was recorded for the case of poly[d(AT)] up to the salt concentration limit of 4.4 mM, where aggregation does not yet interfere. Now experimental data on the B-A dynamics are available for poly[d(AT)] at salt concentrations of 0.18, 0.73, 2.44 and 4.4 mM. In all cases, a spectrum of time constants is found, ranging from ~ 10 μs up to components approaching ~ 1 ms. The relatively small dependence of these data on the salt concentration indicates that electrostatic effects on the kinetics are not as strong as may be expected. The ethanol content at the transition center is a linear function of the logarithm of the salt concentration, and the slope is close to that expected from polyelectrolyte theory. The B-A transition dynamics was also measured in D₂O at a salt concentration of 2.4 mM: the center of the transition is found at 20.0 mol/l H₂O and at 20.1 mol/l D₂O with an estimated accuracy of ± 0.1 mol/l; the spectrum of time constants at the respective transition centers is very similar. The experimental results are discussed regarding the data obtained by molecular dynamics simulations.

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.

The nuclear magnetic resonance of CCCC RNA reveals a right-handed helix, and revised parameters for AMBER force field torsions improve structural predictions from molecular dynamics

The sequence dependence of RNA energetics is important for predicting RNA structure. Hairpins with C(n) loops are consistently less stable than hairpins with other loops, which suggests the structure of C(n) regions could be unusual in the "unfolded" state. For example, previous nuclear magnetic resonance (NMR) evidence suggested that polycytidylic acid forms a left-handed helix. In this study, UV melting experiments show that the hairpin formed by r(5'GGACCCCCGUCC) is less stable than r(5'GGACUUUUGUCC). NMR spectra for single-stranded C(4) oligonucleotide, mimicking the unfolded hairpin loop, are consistent with a right-handed A-form-like helix. Comparisons between NMR spectra and molecular dynamics (MD) simulations suggest that recent reparametrizations, parm99χ_YIL and parm99TOR, of the AMBER parm99 force field improve the agreement between structural features for C(4) determined by NMR and predicted by MD. Evidently, the force field revisions to parm99 improve the modeling of RNA energetics and therefore structure.

Structural, mechanical, and thermodynamic properties of a coarse-grained DNA model

We explore in detail the structural, mechanical, and thermodynamic properties of a coarse-grained model of DNA similar to that recently introduced in a study of DNA nanotweezers [T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, Phys. Rev. Lett. 134, 178101 (2010)]. Effective interactions are used to represent chain connectivity, excluded volume, base stacking, and hydrogen bonding, naturally reproducing a range of DNA behavior. The model incorporates the specificity of Watson-Crick base pairing, but otherwise neglects sequence dependence of interaction strengths, resulting in an "average base" description of DNA. We quantify the relation to experiment of the thermodynamics of single-stranded stacking, duplex hybridization, and hairpin formation, as well as structural properties such as the persistence length of single strands and duplexes, and the elastic torsional and stretching moduli of double helices. We also explore the model's representation of more complex motifs involving dangling ends, bulged bases and internal loops, and the effect of stacking and fraying on the thermodynamics of the duplex formation transition.

Unique Physical Signature of DNA Curvature and Its Implications for Structure and Dynamics

A particularly sensitive birefringence technique is used to analyze a curved DNA fragment with 118 bp and a standard DNA with 119 bp. At salt concentrations from 0.5 to 10 mM, both fragments show the usual negative stationary birefringence and monotonic transients - differences are relatively small. At 100 mM salt the curved DNA shows a positive stationary birefringence and non-monotonic transients with processes having amplitudes of opposite sign, whereas signals of the standard DNA remain as usual. Transients induced by reversal of the field vector indicate the existence of a permanent dipole for the curved DNA. 2-MHz-ac pulses induce a negative stationary birefringence in both DNAs. These results are consistent with calculations on models for curved DNA predicting a quasi-permanent dipole and a positive dichroism/birefringence. The quasi-permanent dipole results from the loss of symmetry in the charge distribution of the curved polyelectrolyte. The appearance of the unique signature of curvature at high salt is mainly due to a strong decrease of the polarizability by about 2 orders of magnitude. The special mode of orientation resulting from the quasi-permanent dipole is expected to contribute to the gel migration anomaly. The time constants of birefringence decay for the curved fragment are shorter than those of the 119 bp fragment by a factor of approximately 1.10 at 0.6 mM salt, whereas this factor is approximately 1.20 at 100 mM Na+. If both fragments were normal DNA with 3.4 A rise per base pair, the factor would be approximately 1.02. At high salt and high electric field strengths the factor increases up to 1.37. The implications for the bending dynamics and the potential to distinguish static from dynamic persistence by field reversal experiments are discussed. The dependence of the curvature on the salt concentration indicated by the time constants is consistent with a clear decrease of the electrophoretic anomaly at decreasing salt concentration.

The nature of "unusual" electro-optical transients observed for DNA

Unusual electro-optical transients have been observed for many different polymers and colloidal systems. These effects provoked serious confusion, because a simple-minded interpretation can be completely misleading. The case of double helical DNA is of particular interest, because DNA has been studied in more detail than other systems and because of its biological function. DNA is subject to bending, which implies a loss of symmetry. Due to its high charge density, non-symmetric conformations must have a non-symmetric distribution of charges leading to a torque of considerable magnitude in the presence of external electric fields. The dipole moment describing this torque must be calculated in a coordinate system with its origin at the center of diffusion. The resulting dipole values are in the range of thousands of Debye units. Because the new dipole type is analogous to but not identical with permanent dipoles, the notation "quasi-permanent" dipole is suggested. Application of this concept, using commonly accepted parameters for DNA and established procedures for calculation of electro-optical transients, leads to "unusual" transients. Thus, these transients must be expected from well-known parameters of DNA double helices. The influence of the quasi-permanent dipole moment may be amplified considerably by hydrodynamic coupling. This effect has been demonstrated for the case of smoothly bent rods. Both model calculations and experiments illustrate the danger of getting data that may be completely misleading. For example, depending on pulse amplitudes and/or pulse lengths, electro-optical decays may be accelerated artificially due to superposition of decay components with opposite amplitudes. Experiments show that unusual transients and apparent acceleration effects disappear, when high frequency sine pulses are used for the electro-optical analysis of DNA. Electro-optical effects depend upon the internal dynamics of the object under investigation. In general, the dynamics of DNA bending was assumed to be fast compared to rotational diffusion. Because stacking rearrangements in single stranded nucleic acids are relatively slow and recently the dynamics of the B-A transition was observed in the time range >1 micros, it is likely that there are also relatively slow rearrangements between bending conformers. Bending transitions are expected to be relatively fast, when there are no activation barriers in the bending pathway, and may be slow, when activation barriers must be passed between bending conformers.

DNA folding and melting observed in real time redefine the energy landscape

We report real-time observations of the folding and melting of DNA by probing two active sites of a hairpin structure, the bases and the stem end, and using an ultrafast T-jump. Studies at different initial temperatures (before, during, and after melting) provide the time scale of water heating (<20 ps), single-strand destacking (700 ps to 2 ns), and hairpin destacking (microseconds and longer) in solutions of various ionic strengths and pH values. The behavior of transient changes gives direct evidence to the existence of intermediate collapsed structures, labile in destacking but compact in nature, and indicates that melting is not a two-state process. We propose a landscape that is defined by these nucleation structures and destacking for efficient folding and melting.

Dynamics of gapped DNA recognition by human polymerase beta

Kinetics of human polymerase beta binding to gapped DNA substrates having single stranded (ss) DNA gaps with five or two nucleotide residues in the ssDNA gap has been examined, using the fluorescence stopped-flow technique. The mechanism of the recognition does not depend on the length of the ssDNA gap. Formation of the enzyme complex with both DNA substrates occurs by a minimum three-step reaction, with the bimolecular step followed by two isomerization steps. The results indicate that the polymerase initiates the association with gapped DNA substrates through the DNA-binding subsite located on the 8-kDa domain of the enzyme. This first association step is independent of the length of the ssDNA gap and is characterized by similar rate constants for both examined DNA substrates. The subsequent, first-order transition occurs at the rate of approximately 600-1200 s(-1). This is the major docking step accompanied by favorable free energy changes in which the 31-kDa domain engages in interactions with the DNA. The 5'-terminal PO(4)(-) group downstream from the primer is not a specific recognition element of the gap. However, the phosphate group affects the enzyme orientation in the complex with the DNA, particularly, for the substrate with a longer gap.

Kinetic mechanisms of rat polymerase beta-ssDNA interactions. Quantitative fluorescence stopped-flow analysis of the formation of the (Pol beta)(16) and (Pol beta)(5) ssDNA binding mode

Kinetics of rat polymerase beta (pol beta) binding to the single-stranded DNA (ssDNA) in the (pol beta)(16) and (pol beta)(5) binding modes has been examined, using the fluorescence stopped-flow technique. Binding of the enzyme to the ssDNA containing fluorescein is characterized by a strong increase of the DNA fluorescence, which provides an excellent signal to quantitatively study the complex mechanism of the ssDNA recognition process. The experiments were performed with a 20-mer ssDNA, which can engage the enzyme in the (pol beta)(16) binding mode, i.e. it encompasses the entire, total DNA-binding site of rat pol beta, and with a 10-mer which binds the enzyme exclusively in the (pol beta)(5) binding mode where only the 8 kDa domain of the enzyme is engaged in interactions with the DNA. The data indicate that the formation of the (pol beta)(16) binding mode occurs by a minimum three-step mechanism with the bimolecular binding step followed by two isomerizations: [formula-see text] A similar mechanism is observed in the formation of the (pol beta)(5) binding mode, although at low salt concentrations there is an additional, slow step in the reaction. The data analysis was performed using the matrix projection operator technique, a powerful method to address stopped-flow kinetics, particularly, amplitudes. The binding modes differ in the free energy changes of the partial reactions and ion effects on transitions between intermediates that reflect different participation of the two structural domains. The formation of both binding modes is initiated by the fast association with the ssDNA through the 8 kDa domain, followed by transitions induced by interactions at the interface of the 8 kDa domain and the DNA. In the (pol beta)(16) binding mode, the subsequent intermediates are stabilized by the DNA binding to the DNA-binding subsite on the 31 kDa domain. The data indicate that interactions of the ssDNA-binding subsite of the 8 kDa domain with the ssDNA, controlled by the ion binding, induce conformational transitions of the formed complexes in both binding modes. The sequential nature of the determined mechanisms indicates a lack of kinetically significant conformational equilibrium of rat pol beta, prior to ssDNA binding.

Kinetic mechanism of the single-stranded DNA recognition by Escherichia coli replicative helicase DnaB protein. Application of the matrix projection operator technique to analyze stopped-flow kinetics

Kinetics of the Escherichia coli primary replicative helicase DnaB protein binding to a single-stranded DNA, in the presence of the ATP non-hydrolyzable analog AMP-PNP, have been performed, using the fluorescence stopped-flow technique. This is the first direct determination of the mechanism of the ssDNA recognition by a hexameric helicase. Binding of the fluorescent etheno-derivative of a ssDNA to the enzyme is characterized by a strong increase of the nucleic acid fluorescence, which provides an excellent signal to quantitatively study the mechanism of ssDNA recognition by the helicase. The kinetic experiments have been performed with a ssDNA 20-mer, depsilonA(pepsilonA)(19), that encompasses the entire, total ssDNA-binding site of the helicase and with the 10-mer depsilonA(pepsilonA)(9), which binds exclusively to the ssDNA strong subsite within the total ssDNA-binding site. Association of the DnaB helicase with the 20-mer is characterized by three relaxation times, which indicates that the binding occurs by the minimum three-step mechanism where the bimolecular binding step is followed by two isomerization steps. This mechanism is described by the equation: Helicase+ssDNAk1/(k1)<-->(k-1)(H-ssDNA)1(k2)<-->(k-2)(H-ssDNA)2 (k3)<-->(k-3)(H-ssDNA)3. The value of the bimolecular rate constant, k(1), is four to six orders of magnitude lower than the value expected for the diffusion-controlled reaction. Moreover, quantitative amplitude analysis suggests that the major conformational change of the ssDNA takes place in the formation of the (H-ssDNA)(1). These results indicate that the determined first step includes formation of the collision and an additional transition of the protein-ssDNA complex, most probably the local opening of the protein hexamer. The data indicate that the binding mechanism reflects the interactions of the ssDNA predominantly through the strong ssDNA-binding subsite. The analysis of the stopped-flow kinetics has been performed using the matrix-projection operator technique, which provides a powerful method to address stopped-flow kinetics, particularly, the amplitudes. The method allowed us to determine the specific fluorescence changes accompanying the formation of all the intermediates. The sequential nature of the determined mechanism indicates the lack of the kinetically significant conformational equilibrium of the DnaB hexamer as well as a transient dissociation of the hexamer prior to the ssDNA binding. The significance of these results for the functioning of the DnaB helicase is discussed.

On the stability of nucleic acid structures in solution: enthalpy-entropy compensations, internal rotations and reversibility

The thermodynamics of self-association (stacking) of free bases and nucleotides, intramolecular stacking in dinucleotides, nearest-neighbour base pair stacking interactions in duplex DNA and RNA, and the formation of hairpin loops illustrate enthalpy/entropy compensations. Large stacking exothermicities are associated with large negative entropy changes that ensure that delta G is small, permitting readily reversible associations in solution. We rationalise enthalpy/entropy compensations with reference to residual motions and torsional vibrations which make a larger entropic contribution to binding when - delta H approximately kT (thermal energy at room temperature), than when - delta H >> kT. We present a factorisation of experimental free energies for helix formation in terms of approximate contributions from the restriction of rotations, hydrophobic interactions, electrostatic interactions due to base stacking, and contributions from hydrogen bonding, and estimate the adverse free energy cost per rotor (mainly entropy) of ordering the phosphate backbone as between 1.9 and 5.4 kJ mol-1 [averaged over 12 rotors per base pair for A-U on A-U stacking (lower limit), and G-C on C-G stacking (upper limit)]. The largest cost is associated with the most exothermic stacking interactions, while the range of values is consistent with earlier conclusions from data on the fusion of hydrocarbon chains (lower value), and with entropy changes in covalent isomerisations of small molecules involving severe restrictions (upper value).

Kinetic and thermodynamic investigation of the (dien)Pd(II)-polycytidylic acid system

The reaction between (dien)PdCl+ and polycytidylic acid was studied using spectroscopic and stopped-flow methods. In neutral solution, the palladium complex binds at the N3 site of the cytosine base and causes a noncooperative disruption of the ordered helical structure of poly(C). Interaction at the phosphate group of the polynucleotide was also demonstrated by using the dye acridine orange as an indicator. The results of this study show that the mechanism previously proposed for cytidine and CMP can be applied to poly(C), taking into account particular features of the polymer (polyelectrolytic nature, structure, etc.). In particular, electrostatic effects seem to play a major role in the interaction with metal ion complexes like (dien)Pd(II).

Anticodon loop of tRNAPhe: structure, dynamics, and Mg2+ binding

The structure, dynamics, and Mg2+ binding reactions of the isolated anticodon hairpin loop from tRNAPhe (yeast) have been analyzed by fluorescence-detected temperature-jump relaxation, melting experiments, and equilibrium sedimentation. Most of the measurements were performed at an ionic strength of 0.15 M and at temperatures below 25 degrees C, where the hairpin loop proved to be stable. A relaxation effect with a time constant of approximately 100 microseconds, indicated by the Wye base fluorescence, is attributed to a conformational change of the anticodon loop and is very similar to a corresponding transition observed previously for the whole tRNAPhe molecule. A Mg2+ binding site reflected by an inner-sphere relaxation process and associated with a strong increase of the Wye base fluorescence closely resembles a corresponding site observed in the complete tRNAPhe and is attributed to a site in the anticodon loop identified by X-ray analysis. In addition to the Mg2+ site in the loop, which is associated with a binding constant of 2 X 10(3) M-1, the existence of sites with a higher affinity is demonstrated by an unusual relaxation effect, showing a minimum in the reciprocal time constant with increasing Mg2+ concentration. The experimental data can be described by a transition between two states and Mg2+ binding to both states resulting in a reaction cycle, which is extended by an additional Mg2+ binding reaction to one of the states. The unusual effect has not been observed for the complete tRNAPhe and is also not observed when Ca2+ is added instead of Mg2+. This result indicates the existence of a conformational change involving Mg2+ inner-sphere complexation.(ABSTRACT TRUNCATED AT 250 WORDS)

Anticodon-anticodon interactions in solution. Studies of the self-association of yeast or Escherichia coli tRNAAsp and of their interactions with Escherichia coli tRNAVal

The temperature-jump method was used to measure the thermodynamic and kinetic parameters of the yeast tRNAAsp (anticodon GUC) duplex, which involves a U/U mismatch in the middle position of the quasi self-complementary anticodon, and of the yeast tRNAAsp (GUC)-Escherichia coli tRNAVal (GAC) complex, in which the tRNAs have complementary anticodons. The existence of the tRNAAsp duplex involving GUC-GUC interactions as evidenced in the crystal structure has now been demonstrated in solution. However, the value of its association constant (Kass = 10(4)M-1 at 0 degrees C) is characteristic of a rather weak complex, when compared with that between tRNAAsp and tRNAVal (Kass = 4 X 10(6) M-1 at 0 degrees C), the effect being essentially linked to differences in the rate constant for dissociation. tRNAAsp split in the anticodon by T1 ribonuclease gives no relaxation signal, indicating that the effects observed with intact tRNA were entirely due to anticodon interactions. No duplex formation was observed with other tRNAs having quasi self-complementary GNC anticodons (where N is C, A or G), such as E. coli tRNAGly (GCC), E. coli tRNAVal (GAC) or E. coli tRNAAla (GGC). This is compatible with the idea that, probably as in the crystal structure, a short double helix is formed in solution between the two GUC anticodons. Because of steric effects, such a complex formation would be hindered if a cytosine, adenine or guanine residue were located in the middle position of the anticodon. Escherichia coli tRNAAsp possessing a modified G residue, the Q base, at the first position of the anticodon, showed a weaker self-association than yeast tRNAAsp but its complex with E. coli tRNAVal was found to be only 1.5 times less stable than that between yeast tRNAAsp and E. coli tRNAVal. Temperature-jump experiments conducted under conditions mimicking those used for the crystallization of yeast tRNAAsp (in the presence of 1.6 M-ammonium sulphate and 3mM-spermine) revealed an important stabilization of the yeast and E. coli tRNAAsp duplexes or of their complexes with E. coli tRNAVal. The effect is due exclusively to ammonium sulphate; it is entropy driven and its influence is reflected on the association rate constant; no influence on the dissociation rate constant was observed. For all tRNA-tRNA complexes, the melting temperature upon addition of ammonium sulphate was considerably increased. This study permits the definition of solution conditions in which tRNAs with appropriate anticodons exist mainly as anticodon-anticodon dimers.

The effect of point mutations affecting Escherichia coli tryptophan tRNA on anticodon-anticodon interactions and on UGA suppression

tRNA species in Escherichia coli that translate codons starting with U contain 2-methyl-thio-N6-isopentenyl-adenosine in position 37, 3' adjacent to the anticodon. The role of this hypermodification in protein synthesis and trp operon attenuation has been investigated. Temperature-jump relaxation methods have been applied to study the interaction between E. coli tRNAPro, with anticodon VGG (V is uridine-5-oxyacetic acid) complementary to that of tRNATrp, and three species of E. coli tRNATrp: wild type tRNATrp (with ms2i6A37 and G24), UGA suppressor tRNATrp (with ms2i6A37 and A24 in the dihydrouridine stem but the same anticodon CCA), and the same suppressor molecule but ms2i6A-deficient as a result of the mutation miaA. Complex formation between tRNAPro and ms2i6A-containing tRNATrp shows thermodynamic parameters close to those found for several other pairs of tRNA with complementary anticodons. However, ms2i6A-deficient tRNATrp makes less stable complexes with tRNAPro, which dissociate eightfold faster. No effect on the complementary anticodon interaction of the mutation in the dihydrouridine stem can be detected. When the tRNA analogous to the opal codon, E. coli tRNASerIV (anticodon VGA) replaces tRNAPro in similar experiments, very weak complexes are observed with both normally hypermodified species of tRNATrp, the wild type and UGA suppressor; these show a lifetime about 50-fold shorter than with tRNAPro, but are again similar. No complex formation is detectable with the ms2i6A-deficient species. This may explain why the hypermodification is necessary for the efficient suppression of the UGA terminator of Q beta coat protein in vitro. The data on complexes with tRNAPro suggest that deficiency in ms2i6A may also reduce the efficiency of UGG reading. Thus, miaA may affect trp operon attenuation by slowing translation of the tandem UGG codons in the leader sequence. Temperature-jump differential spectra suggest that ms2i6 stabilizes the anticodon interaction by improved stacking of base 37.

Thermodynamic studies of RNA stability

Enthalpies and entropies of helix stabilization due to addition of 3' terminal unpaired nucleotides to a CCGG or GGCC core double helix are derived from UV melting studies. The results suggest stacking provides a significant fraction of the free energy of a terminal base pair. The effects of temperature, aggregation, and ionic strength on the determination of thermodynamic parameters are considered. Helix propagation parameters are revised and extended based on recent additions to the data set.

Effect of threonylcarbamoyl modification (t6A) in yeast tRNA Arg III on codon-anticodon and anticodon-anticodon interactions. A thermodynamic and kinetic evaluation

The effect of N-[9-(beta-D-ribofuranosyl) purin-6-ylcarbamoyl]threonine (t6A) adjacent to anticodon U-C-U of yeast tRNA Arg III (where U is a modified U), compared to its unmodified adenosine counterpart, has been evaluated by three independent methods: (a) the polynucleotide-directed binding of tRNA on ribosomes, (b) the ribosome-free trinucleotide binding to the anticodon, (c) the anticodon-anticodon binding test. The results obtained by these three methods indicate a small but significant stabilization effect of t6A on the binding of yeast tRNA Arg III with (a) poly(A,G) in the presence of Escherichia coli ribosomes, (b) free A-G-A triplet, and (c) E. coli tRNA Ser V (anticodon G-G-A). We therefore conclude that the stabilization effect of t6A occurs on U x A and U x G base pairs adjacent to the 5' side of the modified nucleoside, most probably by stacking.

Solvent effects on thermodynamics of double-helix formation in (dG-dC)3

The thermodynamics of double-helix formation by (dG-dC)3 have been measured in aqueous solvent mixtures containing 10 mol % methanol, ethanol, 1-propanol, formamide, N,N-dimethylformamide, or urea and 20 mol % ethanol. Optical activity measurements indicate the conformation of the double helix at 3 degree C is the same in all the solvent mixtures except 20 mol % ethanol. All the cosolvents destabilize the helix relative to water. With 10 mol % alcohol cosolvents, this destabilization is associated with a more unfavorable entropy change averaging approximately 8% and a more favorable enthalpy change averaging approximately 5%. This is consistent with a small contribution of hydrophobic bonding to stability. In contrast, the destabilization by formamide, N,N-dimethylformamide, and urea is associated with a more unfavorable enthalpy change averaging approximately 23% and a more favorable entropy change averaging approximately 21%. Since all three of these cosolvents have dipole moments larger than water, this is consistent with increased competition for dipolar interactions between the nucleic acid bases. None of the results correlate with any one bulk solvent parameter such as surface tension, viscosity, or dielectric constant. With 20 mol % ethanol, optical activity measurements are consistent with a partial B to C form transition. This is associated with a 27% less favorable enthalpy and 25% more favorable entropy for helix formation relative to water. Since the B to C transition is associated with helix dehydration, this may imply a significant contribution of bound water to stability.

Solvent effects on the kinetics and thermodynamics of stacking in poly(cytidylic acid)

The Raman laser temperature-jump technique has been used to measure the kinetics of the coil to helix reaction of poly(cytidylic acid) [poly(C)] in aqueous cosolvent mixtures. The rate of helix formation has a low activation energy and is proportional to reciprocal solvent viscosity. The observations suggest helix formation is rotationally diffusion controlled. The rate of coil formation in poly(C) has an activation energy of approximately 11 kcal/mol, presumably reflecting the electronic stacking interactions which stabilize the helix. Viscous cosolvents, glycerol or sucrose, slow down the rate of coil formation; acetonitrile and formamide at 5 mol % increase the rate relative to that in water. The polar cosolvents may specifically attack a cytosine stack. The absorbance vs. temperature profiles for poly(C) are analyzed with the one-dimensional Ising model. When only optical data are used, the cooperativity parameter, sigma, and the enthalpy, delta H, cannot be uniquely determined. A method is proposed that allows determination of sigma by combining spectroscopic and calorimetric data. The values of sigma derived for poly(C) are between 0.8 and 1.0, and delta H is about -9 kcal/mol of stack. An alternative method using integration of the excess heat capacity curve and extrapolation to fully stacked and random coil species yields a delta H of -7 kcal/mol of stack.

Thermal perturbation differential spectra of ribonucleic acids. III. Chain length effect

Thermal perturbation differential spectra of several adenylic acid oligomers, two single-strand polyribonucleotides, poly(A) and poly(C), and five trinucleoside diphosphates, U-A-A, U-A-G, U-G-A, C-U-C and C-U-A, were obtained and analysed. It is shown that these differential spectra cannot be entirely described by the nearest-neighbour approximation devised from the appropriate mononucleotides and dinucleoside monophosphates. In an attempt to determine the origin of this discrepancy, we have examined the possible optical changes arising from increased electronic interactions, changes in conformation, solvent accessibility or thermodynamic properties. This study indicates that dinucleoside monophosphates on one the hand and trinucleoside diphosphates on the other hand, are separate classes from the conformational point of view. The capacity to assume stacks, absent or negligible in higher oligomers or polymers, makes them poor models for the stacking interaction in longer nucleic acids. It is also shown that in trinucleoside diphosphates, interaction between the two terminal bases arising from bulging out of the middle base is very likely to occur. This type of interaction has to be taken into account in the description of the temperature perturbation differential spectra of trimers.

Stacking properties of a highly hydrophobic dinucleotide sequence, N6, N6-dimethyladenylyl(3' leads to 5')N6, N6-dimethyladenosine, occurring in 16--18-S ribosomal RNA

The thermal denaturation ultraviolet absorption spectra of N6,N6-dimethyladenylyl(3' leads to 5')-N6,N6-dimethyladenosine (m2(6)Apm2(6)A), which is a common sequence in 16--18-S ribosomal RNA, in aqueous buffer at pH 7 have been measured over the temperature range 3-90 degrees C. These data have been used to determine the thermodynamic quantitites associated with the intramolecular stacking equilibria. At 25 degrees C in neutral aqueous solution m2(6)Apmw(6)A exists mainly (about 81%) as a stacked form, so that the stacking interactions are stronger than those in the parent unmethylated adenylyl-(3'-5')adenosine (ApA), where about 52% is stacked. From the parameters of delta H and delta S, it is concluded that 'hidden' hydrophobic inteactions are of prime importance in the enhanced stability of m2(6)Apm2(6)A. Transphosphorylation reaction of ApA and m2(6)Apm2(6)A to form the corresponding cyclic 2',3'-phosphates has been studied. First-order rate constants at 25 degrees C for the reactions, which are base-catalyzed, have been obtained. Insertion of two methyl groups at N-6 of ApA reduces the rate of transphosphorylation. Effects of stacking on rates are discussed in the light of reaction mechanisms.

RNA structure

This review is concerned primarily with the physical structure and changes in the structure of RNA molecules. It will be evident that we have not attempted comprehensive coverage of what amounts to a vast literature. We have tried to stay away from particular areas that have been recently reviewed elsewhere. Citations to and information from them are included, however, so that access to the literature is available. Much of what we treat in depth deals with the crystal structures and solution behaviour of model RNA compounds, including synthetic polymers and molecular fragments such as dinucleoside phosphates. Sequence data on natural RNA are cited, but not in detail. Similarly, apart from tRNA, natural RNAs the structural determinations of which are presently not so far advanced, are not dwelt upon. We have tried to present in detail the available structural data with scaled drawings that permit facile comparisons of molecular geometries.

Conformation of viroids

Viroids are uncoated infectious RNA molecules (MW 107 000-127 000) known as pathogens of certain higher plants. Thermodynamic and kinetic studies were carried out on highly purified viroid preparations by applying UV-absorption melting analysis and temperature jump methods. The thermal denaturation of viroids is characterized by high thermal stability, high cooperativity and a high degree of base pairing. Two relaxation processes could be resolved; a process in the sec range could be evaluated as an independent all-or-none-transition with the following properties: reaction enthalpy= 550 kcal/mol, activation enthalpy of the dissociation = 470 kcal/mol; G : C content = 72 %. These data indicate the existence of an uninterrupted double helix of 52 base pairs. A process in the msec range involves 15 - 25 base pairs which are most probably distributed over several short double helical stretches. A tentative model for the secondary structure of viroids isproposed and the possible functional implications of their physicochemical properties are discussed.

The kinetics of bisulphite modification of reactive residues in E. coli tRNA2Phe

E coli tRNA2Phe was modified at 25 degrees C with 3M sodium bisulphite, pH6.0, for periods of up to 48 hours, Three cytadinine residues, at position 17, 74 and 75 from the 5' end were each deaminated to uridine. The 2-methylthio-N6-isopentenyl adenosine at position 37 formed a 1:1 bi-sulphite addition product which was stable to alkaii. No other residues were permanently modified. The rate of modification of each residue was first order with respect to remaining unmodified nucleotide, the time of half reaction, t1/2, being different for each residue. C17 reaction reacted at twice the rate of cytidine in PolyC, indicating that it occupied a very exposed position in the tRNA.

The kinetics of binding of U-U-C-A to a dodecanucleotide anticodon fragment from yeast tRNA-Phe

The kinetics of U-U-C-A binding to the dodecanucleotide (A-Cm-U-Gm-A-A-Y-A-psi-m5C-U-Gp) isolated from the anticodon region of yeast tRNA-Phe are similar to the kinetics of binding of U-U-C-A to intact tRNA-Phe. A large enhancement in binding constant over that predicted for U-U-C-A-U-G-A-A is observed for both the complexes of dodecanucleotide and tRNA-Phe with U-U-C-A. This strongly suggests that both the anticodon loop in tRNA-Phe and the dodecanucleotide can form four base pairs with U-U-C-A. Furthermore, the enhanced stability cannot be attributed to a special conformation of the anticodon loop, but instead the anticodon loop is probably flexible. A likely explanation for the increased binding is the effect of non-base-paired ends. This increased thermodynamic stability comes from a larger entropy gain rather than a larger enthalpy decrease.

Thermodynamic and kinetic parameters of ion condensation to polynucleotides. Outer sphere complex formed by Mg++ions

The coupling of ion binding to the single strand helix-coil transition in poly (A) and poly (C) is used to obtain information about both processes by ion titration and field-jump relaxation methods. Characterisation of the field-jump relaxation in poly (C) at various concentrations of monovalent ions leads to the evaluation of a stability constant K = 71 M(-1) for the ion binding to the polymer. the rate constant of helix formation is found to be 1.3 X 10(7)s(-1), whereas the dissociation rate is 1.0 X 10(6)s(-1). Similar data are preseented for poly (A) and poly (dA). The interaction of Mg(++) and Ca(++) with poly (A) and poly (C) is measured by a titration method using the polymer absorbance for the indication of binding. The data can be represented by a model with independent binding "sites". The stability constants increase with decreasing salt concentration from 2.7 X 10(4) M(-1) at medium ionic strengths up to 2.7 X 10(7) M(-1) at low ionic strength. The number of ions bound per nucleotide residue is in the range 0.2 to 0.3. Relaxation time constants associated with Mg(++) binding are characterised over a broad range of Mg(++) concentrations form 5 muM to 500 muM. The observed concentration dependence supports the conclusion on the number of binding places inferred from equilibrium titrations. The rate of Mg(++) and Ca(++) association to the polymer is close to the limit of diffusion control kR =1 X 10(10) to 2 X 10(10) M(-1) s(-1)). This high rate demonstrates that Mg(++) and Ca(++) ions do not form inner-sphere complexes with the polynucleotides. Apparently the distance between two adjacent phosphates is too large for a simultaneous site binding of Mg(++) or Ca(++), and inner sphere complexation at a single phosphate seems to be too weak. The data support the view that the ions like Mg(++) and Ca(++) surround the polynucleotides in the form of a mobile ion cloud without site binding.