Local mobility of nucleic acids as determined from crystallographic data. I. RNA and B form DNA

A semi-rigid isoindoline-derived nitroxide spin label for RNA

A new isoindoline-derived benzimidazole nitroxide spin label, ImUm, was synthesized and incorporated into RNA oligoribonucleotides. ImUm is the first example of a conformationally unambiguous spin label for RNA, in which the nitroxide N-O bond lies on the same axis as the single bond used to attach the rigid isoindoline-based spin label to a uridine base. This results in minimal displacement of the nitroxide upon rotation of this single bond, which is a useful property for a label to be used for distance measurements. Continuous-wave (CW) EPR measurements of RNA duplexes containing ImUm indicate a restricted rotation around this single bond, presumably due to an intramolecular hydrogen bond between the benzimidazole N-H and O4 of the uracil. Orientation-selective pulsed electron-electron double resonance (PELDOR, also called double electron-electron resonance, or DEER) distance measurements between two spin labels in two RNA duplexes showed in one case a strong orientation dependence, further confirming the restricted motion of the spin labels in RNA duplexes.

An enhanced set of displacement parameter restraints in CRYSTALS

Crystallographic restraints are widely used during refinement of small-molecule and macromolecular crystal structures. They can be especially useful for introducing additional observations and information into structure refinements against low-quality or low-resolution data (e.g. data obtained at high pressure) or to retain physically meaningful parameter values in disordered or unstable refinements. However, despite the fact that the anisotropic displacement parameters (ADPs) often constitute more than half of the total model parameters determined in a structure analysis, there are relatively few useful restraints for them, examples being Hirshfeld rigid-bond restraints, direct equivalence of parameters and SHELXL RIGU-type restraints. Conversely, geometric parameters can be subject to a multitude of restraints (e.g. absolute or relative distance, angle, planarity, chiral volume, and geometric similarity). This article presents a series of new ADP restraints implemented in CRYSTALS [Parois, Cooper & Thompson (2015), Chem. Cent. J. 9, 30] to give more control over ADPs by restraining, in a variety of ways, the directions and magnitudes of the principal axes of the ellipsoids in locally defined coordinate systems. The use of these new ADPs results in more realistic models, as well as a better user experience, through restraints that are more efficient and faster to set up. The use of these restraints is recommended to preserve physically meaningful relationships between displacement parameters in a structural model for rigid bodies, rotationally disordered groups and low-completeness data.

Fusion RNAs in crystallographic studies of double-stranded RNA from trypanosome RNA editing

Head-to-head fusions of two identical double-stranded fragments of RNA can be designed to self-assemble from a single RNA species and form a double-stranded helix with a twofold rotation axis relating the two strands. These symmetrical RNA molecules are more likely to crystallize without end-on-end statistical packing disorder because the two halves of the molecule are identical. This approach can be used to study many fragments of double-stranded RNA or many isolated helical domains from large single-stranded RNAs that may not yet be amenable to high-resolution studies by crystallography or NMR. We used fusion RNAs to study one (the U-helix) of three functional domains formed when guide RNA binds to its cognate pre-edited mRNA from the trypanosome RNA editing system. The U-helix forms when the 3' oligo(U) tail of the guide RNA (gRNA) binds to the purine-rich, pre-edited mRNA upstream from the current RNA editing site. Fusion RNAs 16-and 32-base pairs in length formed crystals that gave diffraction to 1.37 and 1.05 Å respectively. We provide the composition of a fusion RNA crystallization screen and describe the X-ray data collection, structure determination, and refinement of the crystal structures of fusion RNAs.

TLS from fundamentals to practice

The Translation-Libration-Screw-rotation (TLS) model of rigid-body harmonic displacements introduced in crystallography by Schomaker & Trueblood (1968) is now a routine tool in macromolecular studies and is a feature of most modern crystallographic structure refinement packages. In this review we consider a number of simple examples that illustrate important features of the TLS model. Based on these examples simplified formulae are given for several special cases that may occur in structure modeling and refinement. The derivation of general TLS formulae from basic principles is also provided. This manuscript describes the principles of TLS modeling, as well as some select algorithmic details for practical application. An extensive list of applications references as examples of TLS in macromolecular crystallography refinement is provided.

Models for energy and charge transport and storage in biomolecules

Two models for energy and charge transport and storage in biomolecules are considered. A model based on the discrete nonlinear Schrödinger equation with long-range dispersive interactions (LRI's) between base pairs of DNA is offered for the description of nonlinear dynamics of the DNA molecule. We show that LRI's are responsible for the existence of an interval of bistability where two stable stationary states, a narrow, pinned state and a broad, mobile state, coexist at each value of the total energy. The possibility of controlled switching between pinned and mobile states is demonstrated. The mechanism could be important for controlling energy storage and transport in DNA molecules. Another model is offered for the description of nonlinear excitations in proteins and other anharmonic biomolecules. We show that in the highly anharmonic systems a bound state of Davydov and Boussinesq solitons can exist.

Structural changes of an abasic site in duplex DNA affect noncovalent binding of the spin label ç

The influence of structural changes of an abasic site in duplex DNA on noncovalent and site-directed spin labeling (NC-SDSL) of the spin label ç were examined with electron paramagnetic resonance (EPR) spectroscopy. The binding affinities of ç to sixteen different DNA duplexes containing all possible sequences immediately flanking the abasic site were determined and the results showed that the binding of ç is highly flanking-sequence dependent. In general, a 5′-dG nucleotide favors the binding of the spin label. In particular, 5′-d(G__T) was the best binding sequence whereas 5′-d(C__T) showed the lowest affinity. Changing the structure of the abasic site linker from a tetrahydrofuran analog (F) to the anucleosidic C₃-spacer (C₃) does not appreciably affect the binding of ç to the abasic site. For efficient binding of ç, the abasic site needs to be located at least four base pairs away from the duplex end. Introducing a methyl substituent at N3 of ç did not change the binding affinity, but a decreased binding was observed for both N3-ethyl and -propyl groups. These results will guide the design of abasic site receptors and spin label ligands for NC-SDSL of nucleic acids.

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.

Nucleic-acid structural deformability deduced from anisotropic displacement parameters

The growing numbers of very well resolved nucleic-acid crystal structures with anisotropic displacement parameters provide an unprecedented opportunity to learn about the natural motions of DNA and RNA. Here we report a new Monte-Carlo approach that takes direct account of this information to extract the distortions of covalent structure, base pairing, and dinucleotide geometry intrinsic to regularly organized double-helical molecules. We present new methods to test the validity of the anisotropic parameters and examine the apparent deformability of a variety of structures, including several A, B, and Z DNA duplexes, an AB helical intermediate, an RNA, a ligand-DNA complex, and an enzyme-bound DNA. The rigid-body parameters characterizing the positions of the bases in the structures mirror the mean parameters found when atomic motion is taken into account. The base-pair fluctuations intrinsic to a single structure, however, differ from those extracted from collections of nucleic-acid structures, although selected base-pair steps undergo conformational excursions along routes suggested by the ensembles. The computations reveal surprising new molecular insights, such as the stiffening of DNA and concomitant separation of motions of contacted nucleotides on opposite strands by the binding of Escherichia coli endonuclease VIII, which suggest how the protein may direct enzymatic action.

A Coarse Grained Model for Atomic-Detailed DNA Simulations with Explicit Electrostatics

Coarse-grain (CG) techniques allow considerable extension of the accessible size and time scales in simulations of biological systems. Although many CG representations are available for the most common biomacromolecules, very few have been reported for nucleic acids. Here, we present a CG model for molecular dynamics simulations of DNA on the multi-microsecond time scale. Our model maps the complexity of each nucleotide onto six effective superatoms keeping the "chemical sense" of specific Watson-Crick recognition. Molecular interactions are evaluated using a classical Hamiltonian with explicit electrostatics calculated under the framework of the generalized Born approach. This CG representation is able to accurately reproduce experimental structures, breathing dynamics, and conformational transitions from the A to the B form in double helical fragments. The model achieves a good qualitative reproduction of temperature-driven melting and its dependence on size, ionic strength, and sequence specificity. Reconstruction of atomistic models from CG trajectories give remarkable agreement with structural, dynamic, and energetic features obtained from fully atomistic simulation, opening the possibility to acquire nearly atomic detail data from CG trajectories.

On the relationship between diffraction patterns and motions in macromolecular crystals

The quality of many macromolecular crystal structures published recently has been enhanced through the use of new methods for treating the effects of molecular motion and disorder on diffraction patterns, among them a technique called translation, libration, screw-axis (TLS) parameterization. TLS parameterization rationalizes those effects in terms of domain-scale, rigid-body motions and, interestingly, the models for molecular motion that emerge when macromolecular diffraction data are analyzed this way often make sense biochemically. Here it is pointed out that all such models should be treated with caution until it is shown that they are consistent with the diffuse scatter produced by the crystals that provided the diffraction data from which they derive.

Effects of oligoDNA template length and sequence on binary self-assembly of a nucleotide bolaamphiphile

Templated self-assembly of nucleotide bolaamphiphile 1 (in which a 3'-phosphorylated thymidine moiety is connected to each end of a long oligomethylene chain) with a 10-, 20-, 30-, or 40-meric single-stranded oligoadenylic acid (2, 3, 4, or 5) led to the formation of right-handed helical nanofibers in 0.1x Tris/EDTA (TE) buffer solutions. The helical pitch increased as the length of the oligoadenylic acid template increased. DNA composed of oligoadenylic and oligocytidylic acid sequences (6, 7, and 8) also acted as templates to induce the formation of helical nanofiber structures. The diameter of the nanofibers remained constant (6-6.6 nm) irrespective of the template used. The binary self-assembly of 1 with 4 also produced higher-order, double-stranded nanofibers.

Progress towards a DNA-based model system for the study of electron spin-spin interactions

A DNA-based model system is described for studying electron spin-spin interactions between a paramagnetic metal ion and a nitroxide spin label. The modified base deoxythymidine-EDTA (dT-EDTA) chelates the divalent or trivalent metal ion and produces a new feature in the circular dichroism (CD) spectra that makes it possible to monitor local DNA melting. Based on the results of optical and electron paramagnetic resonance (EPR) experiments, we find that the terminus of the DNA duplex that incorporates dT-EDTA and the spin-label melts at a higher temperature than the rest of the DNA duplex. EPR microwave progressive power saturation experiments performed at 77 K are consistent with the specific binding of Dy(III) at the EDTA site and an intramolecular dipole-dipole interaction between the nitroxide spin-label and the chelated Dy(III). This model system should be suitable for studying the relaxation properties of metal ions by saturation-recovery EPR.

Ultrafast dynamics in DNA: "fraying" at the end of the helix

The dynamics of the electric fields in the interior of DNA are measured by using oligonucleotides in which a native base pair is replaced by a dye molecule (coumarin 102) whose emission spectrum is sensitive to the local electric field. Time-resolved measurements of the emission spectrum have been extended to a six decade time range (40 fs to 40 ns) by combining results from time-correlated photon counting, fluorescence up-conversion, and transient absorption. Recent results showed that when the reporter is placed in the center of the oligonucleotide, the dynamics are very broadly distributed over this entire time range and do not show specific time constants associated with individual processes (Andreatta, D.; et al. J. Am. Chem. Soc. 2005, 127, 7270). This paper examines an oligonucleotide with the reporter near its end. The broadly distributed relaxation seen before remains with little attenuation. In addition, a new relaxation with a well-defined relaxation time of 5 ps appears. This process is assigned to the rapid component of "fraying" at the end of the helix.

Modeling furanose ring dynamics in DNA

Determination of the conformational flexibility of the furanose ring is of vital importance in understanding the structure of DNA. In this work we have applied a model of furanose ring motion to the analysis of deuterium line shape data obtained from sugar rings in solid hydrated DNA. The model describes the angular trajectories of the atoms in the furanose ring in terms of pseudorotation puckering amplitude (q) and the pseudorotation puckering phase phi. Fixing q, the motion is thus treated as Brownian diffusion through an angular-dependent potential U(phi). We have simulated numerous line shapes varying the adjustable parameters, including the diffusion coefficient D, pseudorotation puckering amplitude q, and the form of the potential U(phi). We have used several forms of the potential, including equal double-well potentials, unequal double-well potentials, and a potential truncated to "second order" in the Fourier series. To date, we have obtained best simulations for both equilibrium and nonequilibrium (partially relaxed) solid-state deuterium NMR line shapes for the sample [2' '-2H]-2'-deoxycytidine at the position C3 (underlined) in the DNA sequence [d(CGCGAATTCGCG)]2, using a double-well potential with an equal barrier height of U(0) = 5.5k(B)T ( approximately 3.3 kcal/mol), a puckering amplitude of q = 0.4 A, and a diffusion coefficient characterizing the underlying stochastic jump rate D = 9.9 x 10(8) Hz. Then the rate of flux for the C-D bond over the barrier, i.e., the escape velocity or the overall rate of puckering between modes, was found to be 0.7 x 10(7) Hz.

Searching for the magic bullet against cancer: the butyrate saga

n-Butyric acid and its "polymorphic" derivatives have been largely but somehow "blindly" studied in oncology and in red cell diseases with consistent results through decades indicating a strong maturative effect determined by enhancement of gene transcription. Although these effects have been observed mainly in vitro, the relative absence of systemic toxicity of butyrates render these compounds appealing as specific therapeutic agents. More interestingly, their specific mechanism of action, i.e. inhibition of histone deacetylase and de-repression of transcription represents at present an unique tool for diseases such as acute leukemias which are characterised by a disregulation of co-repressors and co-activators of gene transcription. More insight into specificity and modalities of action of different butyrate derivatives may be a guarantee for excellent tailored antileukemic therapy in the future.

RNA crystallography

The current state of three-dimensional structure analysis of RNA by x-ray crystallography is summarized. The methods of sample preparation, crystallization, data collection, and structure solution are discussed, followed by a review of the RNA structures that have been determined and of common structural features, and finally, an appraisal of future prospects for x-ray crystal structure analysis of RNA.

The dynamic impact of CpG methylation in DNA

Solid-state deuterium NMR is used to investigate perturbations of the local, internal dynamics in the EcoRI restriction binding site, -GAATTC- induced by cytidine methylation. Methylation of the cytidine base in this sequence is known to suppress hydrolysis by the EcoRI restriction enzyme. Previous solid-state deuterium NMR studies have detected large amplitude motions of the phosphate-sugar backbone at the AT-CG junction of the unmethylated DNA sequence. This study shows that methylation of the cytidine base in a CpG dinucleotide reduces the amplitudes of motions of the phosphate-sugar backbone. These observations suggest a direct link between suppression of the amplitudes of localized, internal motions of the sugar-phosphate backbone of the DNA and inhibition of restriction enzyme cleavage.

The opening of a single base without perturbations of neighboring nucleotides: a study on crystal B-DNA duplex d(CGCGAATTCGCG)2

In this work we explore the possibility of the opening of a single base without perturbation of its neighboring nucleotides. Low energy base opening into the grooves can be accomplished by rotation of the relevant backbone and glycosidic bond torsion angles. The pathway has been determined by identifying zeta torsion angle as the reaction coordinate together with the accompanying geometric requirement that guides the displacement of other torsion angles. Our study on Dickerson dodecamer duplex d(CGCGAATTCGCG)2 showed that all bases with normal equilibrium zeta can be rotated by approximately 30 degrees, corresponding to approximately 3.5A base displacement, towards the major groove. Such an opening extent is comparable with estimated amplitudes of local angular motions in DNA bases from NMR experiments, which might facilitate proton exchange. The computed base opening energy barrier is also comparable with measured base pair opening enthalpy. These results indicate possible relevance of the pathway studied in this work with experimentally observed base pair opening process. Our analysis also showed a preference for base opening along the major groove and an abnormal opening behavior for bases with unusual equilibrium zeta torsion angle.

Site-specific dynamics in DNA: experiments

This chapter reviews the dynamics information obtained from experimental magnetic resonance studies of site-specifically labeled duplex DNA. A previous review (43) discusses the dynamics of duplex DNA; it develops a theory that shows how magnetic resonance experiments are used to detect those dynamics. The methods for obtaining information about dynamics as well as a summary of what is now known about the site-specific dynamics of DNA are presented. This review contains two methods sections which present results using electron paramagnetic resonance and nuclear magnetic resonance active probes.

Overall and internal dynamics of DNA as monitored by five-atom-tethered spin labels

Electron paramagnetic resonance (EPR) spectra of the two-atom-tethered six-membered ring thymidylate spin label (DUMTA) incorporated into duplexes of different sizes were found to display a helix length dependence and a local-order parameter S = 0.32 +/- 0.01 for B-DNA based on the dynamic cylinder model (Keyes, R. S., and A. M. Bobst. 1995. Detection of internal and overall dynamics of a two-atom-tethered spin-labeled DNA. Biochemistry. 34:9265-9276). This sensitivity to size, which reflects global tumbling, is now reported for the more flexible five-atom-tethered five-membered ring thymidylate spin label (DUAP) that can be readily incorporated enzymatically and sequence specifically into nucleic acids of different sizes. The DUAPs containing B-DNA systems were simulated with the same dynamic cylinder model, giving S = 0.20 +/- 0.01 for the more flexibly tethered spin label. This shows that S is dependent on tether length but not on global motion. An analysis with the same motional model of the B-Z transition in a (dG-dC)n polymer containing the five-atom-tethered six-membered ring cytidylate spin label (DCAT) (Strobel, O. K., R. S. Keyes, and A. M. Bobst. 1990b. Base dynamics of local Z-DNA conformations as detected by electron paramagnetic resonance with spin-labeled deoxycytidine analogues. Biochemistry. 29:8522-8528) revealed an increase in S from 0.15 +/- 0.01 to 0.26 +/- 0.01 in response to the B- to Z-DNA transition. This indicates that S is not only sensitive to tether length, but also to conformational changes in DNA. Both the DUAP- and the DCAT-labeled systems were also simulated with a base disk model. From the DUAP spectral series, the perpendicular component of the correlation time tau perpendicular describing the spin-labeled base diffusion was found to be sensitive to global tumbling, confirming earlier results obtained with DUMTA. The DCAT polymer results demonstrated that tau perpendicular monitors a conformational change from B- to Z-DNA, indicating that tau perpendicular is also sensitive to local base dynamics. These results confirm that the dynamics of five-atom-tethered nitroxides are coupled to the nucleic acid dynamics and, as with two-atom-tethered spin labels, can be characterized by S and tau perpendicular. The analyses of both spin-labeled systems provide good evidence for spin-labeled base motions within double-stranded DNA occurring on the nanosecond time scale, and establish that both labels can be used to monitor changes in global tumbling and local order parameter due to variations in DNA conformation and protein-DNA interactions.

Distance measurements in nucleic acids using windowless dipolar recoupling solid state NMR

A windowless, homonuclear dipolar recoupling pulse sequence (DRAWS) is described and a theoretical basis for describing its recoupling performance is developed using numerical techniques. It is demonstrated that DRAWS recouples weak dipolar interactions over a broad range of experimental and molecular conditions. We discuss two spectroscopic control experiments, which help to take into account effects due to insufficient proton decoupling, relaxation, and static dipolar couplings to nearby 13C spins at natural abundance. Finally DRAWS is used in combination with selective 13C labeling to measure 13C-13C distances in five doubly labeled DNA dodecamers, [d(CGCGAAT*T*CGCG)]2, which contain the binding site for the restriction enzyme EcoRI. The longest distance reported is 4.8 A. In most cases the distances agree well with those derived from X-ray crystallographic data, although small changes in hydration level can result in relatively large changes in internuclear distances.

Large-amplitude picosecond anisotropy decay of the intrinsic fluorescence of double-stranded DNA

The conformational flexibility of the DNA double helix is of great interest because of its potential role in protein recognition, packaging into chromosomes, formation of photodefects, and interaction with drugs. Theory finds that DNA is very flexible; however, there is a scarcity of experimental results that examine intrinsic properties of the DNA bases for the inherent flexibility in solution. We have studied the dynamics of poly(dA).poly(dT) and (dA)20.(dT)20 in a 50 mM cacodylate, 0.1 M NaCl, pH 7 buffer by using the time-correlated picosecond fluorescence anisotropy of thymine selectively excited at 293 nm. For both nucleic acids, a large-amplitude biphasic decrease in the anisotropy is observed that has a very fast, large-amplitude component on the picosecond time scale and a slower, smaller-amplitude component on the nanosecond time scale. These modes are sensitive to sucrose concentration, and are greatly attenuated at 77% sucrose by volume. This observation suggests that motions of the bases make a significant contribution to the observed fluorescence depolarization (in the absence of sucrose). Measurements on the single-stranded systems poly(dT) and (dT)20 reveal a much smaller amplitude of the very fast depolarization mode. These observations are consistent with a mechanism that involves concerted motions in the interior of the double-stranded systems.

How long does DNA keep the memory of its conformation? A time-dependent canonical correlation analysis of molecular dynamics simulation

The time dependence of the correlation between motions of different parts of DNA is analyzed from a 200 ps molecular dynamics simulation of the double-stranded self-complementary d(CTGATCAG) in the B form. Each nucleotide is decomposed into three subunits corresponding to the furanose ring (SU), the base (BA), and the backbone (SK). The motion of each subunit is considered as the superimposition of rigid body translation, rigid body rotation, and internal deformation. Canonical time-dependent correlation functions calculated with coordinates describing the different components of the subunits motion are defined and computed. This allows us to probe how long a particular type of motion of one subunit influences the other types of motions of other subunits (cross correlation functions) or how long a particular subunit keeps the memory of its own conformation or location (autocorrelation functions). From autocorrelation analysis it is found that deformation decorrelates within a few tenths of picoseconds, rotational correlation times are on the order of 8 ps, while translational motions are long-time correlated. The deformation of a subunit is not correlated to the deformation of another one (at the 200 ps time scale of our simulation), but influences slightly their translation and orientation as time increases.

Rigid-body motions of sub-units in DNA: a correlation analysis of a 200 ps molecular dynamics simulation

A 200 ps molecular dynamics simulation of the B-form double stranded self-complementary octanucleotide d(CTGATCAG) is analyzed in terms of correlated motions using the canonical analysis approach. Each nucleotide is decomposed in three sub-units corresponding to the base, the sugar ring and the backbone respectively. The correlation between the full dynamics of two sub-units was found to decrease as their mutual distance increases. The interpretation of the full dynamics of sub-units as the superimposition of rigid-body motions (translation and orientation) and deformation shows that the main source of correlation is rigid-body motions. Correlation between sub-units deformation is weak and practically vanishes for sub-units belonging to non-adjacent nucleotides. It is also shown that the correlation is much more important for sub-units of the same strand than of opposite strands. We conclude that the internal dynamics of the octanucleotide may be well described by rigid-body motions, the sub-units deformation having only local influence whereas sub-units translation and rotation have repercussion to long distances. The results presented in this study suggest how the number of degrees of freedom may be reduced for simulating long-time dynamics of oligonucleotides.

Molecular dynamics simulations of dinucleoside and dinucleoside-drug crystal hydrates

Molecular dynamics simulations have been performed on the dinucleoside monophosphates rGpC and dCpG, the latter in its intercalation complex with the acridine drug proflavine. The simulations were performed on the crystal structures, with crystallographically-located solvent molecules. It was found that satisfactory results were best obtained with restraints placed on the movements of the water molecules. Motions of individual atoms have been examined in terms of rms fluctuations and anisotropy and correlation functions. Relative motions of groups (phosphates, sugars, bases and proflavine molecules) have been analysed.

Simulation analysis of the stability mutant R96H of T4 lysozyme

Free energy simulation methods are used to analyze the effects of the mutation Arg 96----His on the stability of T4 lysozyme. The calculated stability change and the lack of significant structural rearrangement in the folded state due to the mutation are in agreement with experimental studies [Kitamura, S., & Sturtevant, J. M. (1989) Biochemistry 28, 3788-3792; Weaver, L. H., et al. (1989) Biochemistry 28, 3793-3797]. By use of thermodynamic integration, the contributions of specific interactions to the free energy change are evaluated. It is shown that a number of contributions that stabilize the wild type or the mutant partially cancel in the overall free energy difference; some of these involve the unfolded state. Comparison of the results with conclusions based on structural and thermodynamic data leads to new insights into the origin of the stability difference between wild-type and mutant proteins. Of particular interest is the importance of the contributions of more distant residues, solvent water, and the covalent linkage of the mutated amino acid. Also, the analysis of the interactions of Arg/His 96 with the C-terminal end of a helix (residues 82-90) makes it clear that the nearby carbonyl groups (Tyr 88 and Asp 89) make the dominant contribution, that the amide groups do not contribute significantly, and that the helix-dipole model is inappropriate for this case.

Dynamics in synthetic oligonucleotides. A solid-state deuterium NMR study

Solid-state 2H NMR spectroscopy was used to investigate dynamics in the [methyl-2H]thymidine-labeled oligonucleotide [d(CGCGAAT*T*CGCG)]2. Quadrupole echo line shapes, spin-lattice relaxation, and quadrupolar echo decay rates were investigated as a function of hydration W (moles of water/moles of nucleotide) between 0 and approximately 30. The amplitude of the base motion, modeled as a fast four-site libration, or diffusion in a cone, increased slightly with higher levels of hydration. A slower component of motion about the helix axis appears at W approximately 10 and increases in rate and amplitude, leading to the intermediate rate line shape observed at W approximately 21.

Solution structure of the parallel-stranded duplex oligonucleotide alpha-d(TCTAAAC)-beta-d(AGATTTG) via complete relaxation matrix analysis of the NOE effects and molecular mechanics calculations

The solution structure of the duplex formed by the association of the unnatural oligonucleotide alpha-d(TCTAAAC) with its natural and parallel complementary sequence beta-d(AGATTTG) was investigated by nuclear magnetic resonance spectroscopy and constrained molecular mechanics calculations. The structure was refined on the basis of interproton distances determined by NOE measurements for a series of mixing times. The NOE values were converted to distances by using the complete 134 x 134 relaxation matrix including all proton dipole-dipole interactions and spin diffusion. The computation of the relaxation matrix requires the Cartesian coordinates of the oligonucleotide, which are not known, a priori. To avoid this ambiguity, we used an iterative procedure in which the new distance constraints are obtained by using the complete relaxation matrix calculated from the previous structure. After three iterations, the process converged. The unnatural duplex alpha-d(TCTAAAC)-beta-d(AGATTTG) adopts in solution a right-helical structure with Watson-Crick base pairing, an anti conformation on the glycosyl linkage on the beta-strand, a syn conformation on the alpha-strand, and a 3'-exo conformation of the deoxyriboses for both sugar anomers. The three-dimensional structure obtained allowed us to describe the local heterogeneity of the duplex.

Crystal lattice packing is important in determining the bend of a DNA dodecamer containing an adenine tract

The crystal structure of a DNA duplex dodecamer d(CGCAAAAATGCG) and its complementary strand has been determined at 2.6-A resolution. Although our goal was to deduce the structural features of the static bending of the helical axis exhibited by adenine-tract structures in solution, we conclude that the overall bend of 20 degrees in the direction of the major groove observed here arises from the forces associated with crystal packing. An isomorphous dodecamer brominated on one strand provides experimental evidence that this asymmetric sequence is positioned in two orientations in the crystal lattice that are related by a 180 degrees rotation around the pseudodyad axis of the sequence. The bend in these two differently positioned DNA molecules depends on their orientation in the crystal, not on their sequence. As with previously determined structures containing adenine tracts, the adenine and thymine base pairs are highly propeller twisted. The N-6 of the adenine comes within hydrogen bonding distance of the O-4 of thymine one step down the helix, facilitating the formation of a series of bifurcated hydrogen bonds within the adenine tract. The adenine tract is relatively straight and the bending is localized outside this region.

Molecular mechanics model of supercoiled DNA

We describe a pseudo-atomic model of supercoiled DNA. Each base-pair of the DNA is represented in the model by three particles placed in a plane. The particle triplets are stacked to model stacked base-pairs in double-helical DNA, and closed circular conformations are generated to investigate supercoiling. This model is less detailed than all-atom models, which are too computationally demanding to be used to study supercoiling. On the other hand, this model contains details at the base-pair level and is therefore more elaborate than elastomechanical models. A potential energy function is written in terms of a set of internal co-ordinates defined to resemble a limited number of helical parameters. The modeled helical parameters, helical twist, base-roll, tilt and rise, are the most important parameters of the global shape of DNA. Experimentally measured mechanical properties of DNA are used to define the forces holding the particles together. We then use a procedure incorporating energy minimization and molecular dynamics to locate low energy conformations of the model DNA. The model was found to behave very much like rubber-tubing and elastomechanical models. The conformations and the effects of supercoiling pressure (a number proportional to the degree to which the total twist of the DNA has been altered from its natural value) on these conformations are all very similar to those observed in the latter two models. We also used this model to examine the effects of supercoiling pressure, base-sequence and mechanical properties on the conformations and energies of five sequences. The sequences studied include models of naturally straight DNA and DNA with static or natural bends.

Static disorder and librational motions of the purine bases in films of oriented Li-DNA

Solid-state 2H nuclear magnetic resonance line shapes have been obtained from folded films of oriented Li-DNA molecules with the purine bases selectively labeled with deuterium at the 8-position. From line shape simulations, the static base tilts as well as the anisotropic motional amplitudes were determined as a function of hydration level and temperature. It was found that the average tilt angle of the bases is close to 0 degrees and at a hydration of ten water molecules per nucleotide the distribution width of tilt angles about this average cannot be larger than 9 degrees (standard deviation). A slightly increased distribution width is observed at low hydration levels. The motional amplitudes are hydration dependent, with the tilting motion ranging from 4 degrees for the driest, up to 15 degrees for the wettest sample, and slightly larger amplitudes are observed for the twisting motion. The amplitude of the twisting motion is unaffected by a temperature decrease down to -60 degrees C, in contrast to the tilting motion that is suppressed at low temperatures.

Local mobility of nucleic acids as determined from crystallographic data. III. A daunomycin-DNA complex

The local mobility of the complex between the anti-tumor drug daunomycin and a DNA hexanucleotide duplex of sequence d(CpGpTpApCpG)2 has been determined by anisotropic refinement of single crystal X-ray diffraction data of 1.2 A resolution (1 A = 0.1 nm). The directions and amplitudes of the local motion indicate that changes in mobility of DNA due to daunomycin binding are primarily limited to the residues forming the intercalation site and do not propagate to the neighboring residues. The intercalated daunomycin ring system (aglycone) is rigidly fixed in the base stack, apparently serving as an anchor for the amino sugar segment of the drug which is one of the most mobile regions of the entire complex. The high flexibility of this amino sugar may be important for inhibition of replication and transcription not only by sterically blocking the minor groove, but also by allowing nonproductive interactions to be formed with various polymerases or other DNA-binding proteins. The crystallographic model is improved sufficiently by the rigid group anisotropic refinement to allow additional bound water molecules to be located.

Base sequence effects in double helical DNA. I. Potential energy estimates of local base morphology

A series of potential energy calculations have been carried out to estimate base sequence dependent structural differences in B-DNA. Attention has been focused on the simplest dimeric fragments that can be used to build long chains, computing the energy as a function of the orientation and displacement of the 16 possible base pair combinations within the double helix. Calculations have been performed, for simplicity, on free base pairs rather than complete nucleotide units. Conformational preferences and relative flexibilities are reported for various combinations of the roll, tilt, twist, lateral displacement, and propeller twist of individual residues. The predictions are compared with relevant experimental measures of conformation and flexibility, where available. The energy surfaces are found to fit into two distinct categories, some dimer duplexes preferring to bend in a symmetric fashion and others in a skewed manner. The effects of common chemical substitutions (uracil for thymine, 5-methyl cytosine for cytosine, and hypoxanthine for guanine) on the preferred arrangements of neighboring residues are also examined, and the interactions of the sugar-phosphate backbone are included in selected cases. As a first approximation, long range interactions between more distant neighbors, which may affect the local chain configuration, are ignored. A rotational isomeric state scheme is developed to describe the average configurations of individual dimers and is used to develop a static picture of overall double helical structure. The ability of the energetic scheme to account for documented examples of intrinsic B-DNA curvature is presented, and some new predictions of sequence directed chain bending are offered.