Analysis of conformational parameters in nucleic acid fragments. II. Co-crystal complexes of nucleic acid bases.

Studies have been made of conformational parameters in co-crystal complexes and compounds of nucleic acid bases in which there is the possibility of formation of hetero-base-pairs. Using published data extracted from the Cambridge structural database, a total of 37 base-pairs were found, of which 25 were hetero-pairs and 12 homo-pairs. These base-pairs were subject to analysis to reveal hydrogen bond parameters, propeller twist, buckle and C1'-C1' separation (or a similar parameter if C1' atoms were not present). Hetero-pairs were found to show larger twists than homo-pairs, the magnitude of twist being unrelated to hydrogen bond parameters or buckle value. The propeller twisting is less pronounced in these nucleic acid bases than in nucleosides, but still has a significant magnitude. Propeller twisting in hetero-pairs is found to be larger than in homo-pairs. Hetero-pairs appear to be formed preferentially in competitive situations.     

Hydration of oligonucleotides in crystals

The solvent molecules found around crystallized oligonucleotides after X-ray refinement are analysed in terms of interaction sites to bases, phosphates and sugars in the three main forms of nucleic acid structures, the A-form, the B-form and the Z-form. The average numbers of contacts to nucleic acid atoms made by solvent molecules are identical in the three forms, but it appears that the average number of contacts solvent molecules make with each other depends on the resolution of the structure. The phosphate anionic oxygen atoms are the most hydrated, while the O(3′) and O(5′) backbone atoms and the ring oxygen atom O(4′) are the least hydrated. Among the hydrophilic atoms of the bases, there is a modulation of the relative water affinities with the nucleic acid form. Numerous hydration sites are such that water molecules can bridge hydrophilic atoms of the same residue, of adjacent residues on the same strand, of distant residues on the two strands, or belonging to symmetry-related residues. Through the helical periodicity of the nucleic acid structure, those bridges can lead to regular and striking hydration networks involving several water molecules and characteristic of the nucleic acid form. Solvent dynamics, as seen by temperature factor versus occupancy plots, seems intimately related to nucleic acid structure and dynamics, since they depend on hydration sites around the nucleic acids.     

Nucleic Acid Structure Analysis: A Users Guide to a Collection of New Analysis Programs

Abstract

Common nomenclature describing the geometry of nucleic acid structures was established at a 1988 EMBO Workshop on DNA Curvature and Bending (Diekmann, S. (1988) J. Mol. Biol. 208, 787–791; Diekmann, S. (1989) The EMBO Journal 8, 1–4; Sarma, RH. (1988) J. Biomol. Structure & Dynamics 6, 391–395; Dickerson, R.E. (1989) J. Biomol. Structured Dynamics 6, 627–634; Dickerson, RE. et al. (1989)Nuc. Acids Res. 17, 1979–1803). We have subsequently developed and incorporated sophisticated mathematics in a computer program to calculate the parameters described by the guidelines. The program calculates all the local parameters relating complementary bases and neighboring base and base pairs in both Cartesian and helical coordinate frames. In addition, the main mathematical property requested by the EMBO guidelines—that the magnitude of the parameters be independent of strand or direction of measurement—is accomplished without the use of a midway coordinate frame for the rotational parameters. The mathematics preserve the physical intuition used in defining the parameters; in particular, the rotational parameters are true rotations based on a simple physical model (rotation at constant angular velocity for a unit amount of time), not Euler angles or angles between vectors and planes as is the case with other approaches. As a result, the mathematical equations are symmetric with the property that a 5° tilt is the same as a 5° roll or a 5° twist, except that the rotations take place about different axes. In other approaches, a 5° tilt can mean a different amount of net rotation from a 5° roll or a 5° twist. In addition, a great deal of flexibility is built into the program so that the user has control over the analysis, including the input format, the coordinate frame used for the base pairing relationship, the point about which the rotations are performed, and which geometric relationships are analyzed. While there is a great deal of flexibility, the program is easy to use. Interactive queries and user accessible files make the options in the program very convenient to tailor to individual needs. In addition, there is also a program that calculates bond lengths, valence angles, and torsion angles along the nucleic acid backbone, and within the sugar and base rings. Another program ‘learns’ the identities of the bond lengths, valence angles, and torsion angles that the user would like to determine. This last program is especially useful for calculating the hydrogen bonds between atoms in complementary strands as well as the unusual hydrogen bonds found in recently determined nucleic acid NMR structures or within protein/nucleic acid complexes.     

Model studies of interactions between nucleic acids and proteins: hydrogen bonding of amides with nucleic acid bases.

The formation of hydrogen bonded complexes between nucleic acid bases and acetamide has been studied by nuclear magnetic resonance in CDC13 at different temperatures. Pairs of hydrogen bonds are formed when acetamide binds to nucleic acid bases. Thermodynamic parameters have been computed and compared to those obtained for the association of carboxylic acids with nucleic acid bases. The role of hydrogen bonded complexes in the association of proteins with nucleic acids is discussed.     

NMR studies of dynamics in RNA and DNA by 13C relaxation

RNA and DNA molecules experience motions on a wide range of time scales, ranging from rapid localized motions to much slower collective motions of entire helical domains. The many functions of RNA in biology very often require this molecule to change its conformation in response to biological signals in the form of small molecules, proteins or other nucleic acids, whereas local motions in DNA may facilitate protein recognition and allow enzymes acting on DNA to access functional groups on the bases that would otherwise be buried in Watson-Crick base pairs. Although these statements make a compelling case to study the sequence dependent dynamics in nucleic acids, there are few residue-specific studies of nucleic acid dynamics. Fortunately, NMR studies of dynamics of nucleic acids and nucleic acids-protein complexes are gaining increased attention. The aim of this review is to provide an update of the recent progress in studies of nucleic acid dynamics by NMR based on the application of solution relaxation techniques.     

DIGICALC: a restriction fragment analysis program.

DIGICALC is a program designed to aid in the acquisition, storage, and analysis of nucleic acid restriction fragment data. The chief considerations during program design were (i) ease of use for people with varying degrees of computer experience, (ii) minimal hardware requirements (e.g. an IBM PC), (iii) portability and ease of modification, and (iv) improved functionality in sizing and comparing restriction fragments over manual methods. The program accepts manual or digitizer input of nucleic acid fragment mobility, calculates the fragments' sizes, and provides the means to search the fragment database and to produce charts of fragment sizes.     

Fluorescence energy transfer monitored competitive equilibria of nucleic acids: applications in thermodynamics and screening

Precise thermodynamic characterization of nucleic acid complex stability is required to understand a variety of biologically significant events as well as to exploit the specific recognition capabilities of nucleic acids in biotechnology, diagnostics, and therapeutics. The development of a database of nucleic acid thermodynamics with sufficient precision to foster further developments in these areas requires new and improved measurement techniques. The combination of a competitive equilibrium titration with fluorescence energy transfer based detection provides a method for precise measurement of differences in free energy values for nucleic acid duplexes that far exceeds in precision those accessible via conventional methods. The method can be applied to detect and to characterize any deviation in a nucleic acid that alters duplex stability. Such deviations include, but are not limited to, mismatches; single nucleotide polymorphisms (SNP); chemically modified nucleotide bases, sugars or phosphates; and conformational anomalies or folding motifs, such as, loops or hairpins.     

Sequence errors described in GenBank: a means to determine the accuracy of DNA sequence interpretation.

The accuracy of nucleic acid sequence data interpretation was determined by assessing and quantifying the discrepancies reported in the GenBank database. This permitted the calculation of an Error Rate (ER) for nucleic acid sequence determination. If one assumes that most entries (TB, Total Bases) were independently verified or those without reported discrepancies were correct, the ER is 0.368 errors per 1000 bases. However, if one assumes that only those sequences with reported discrepancies (TBIQ, Total Bases from entries In Question) were verified and are thus correct, the ER is 2.887 errors per 1000 bases. This establishes the first set of limit boundaries of the ER for sequence interpretation and sequence errors within the GenBank database and provides the foundation for future assessments and the monitoring of sequence data accumulation. In addition, the ER measure provides a basis to evaluate the efficiency and merit of present and future automated nucleic acid sequencing technologies which will have a direct impact upon the final outcome of the "Human Genome Initiative".     

Temperature-dependent molecular dynamics and restrained X-ray refinement simulations of a Z-DNA hexamer

Molecular dynamics simulations of the Z-DNA hexamer 5BrdC-dG-5BrdC-dG-5BrdC-dG were performed at several temperatures between 100 K and 300 K. Above 250 K, a strong sequence-dependent flexibility in the nucleic acid is observed, with the guanine sugar and the phosphate of GpC sequences much more mobile than the cytosine sugar and phosphate of CpG sequences. At 300 K, the hexamer is in dynamic equilibrium between several Z forms, including the crystallographically determined ZI and ZII forms. The local base-pair geometry, however, is not very variable, except for the roll of the base-pairs. The hexamer molecular dynamics trajectories have been used to test the restrained parameter crystallographic refinement model for nucleic acids. X-ray diffraction intensities corresponding to observed diffraction data were computed. The average structures obtained from the simulations were then refined against the calculated intensities, using a restrained least-squares program developed for nucleic acids in order to analyse the effects of the refinement model on the derived quantities. In general, the temperature dependence of the atomic fluctuations determined directly from the refined Debye-Waller factors is in reasonably good agreement with the results obtained by calculating the atomic fluctuations directly from the Z-DNA molecular dynamics trajectories. The agreement is best for refinement of temperature factors without restraints. At the highest temperature studied (300 K), the effect of the refinement on the most mobile atoms (phosphates) is to significantly reduce the mean-square atomic fluctuations estimated from the refined Debye-Waller factors below the actual values (less than (delta r)2 greater than congruent to 0.5 A2). Analysis of the temperature-dependence of the mean-square atomic fluctuations provides information concerning the conformational potential within which the atoms move. The calculated temperature-dependence and anharmonicity of the Z-DNA helix are compared with the results observed for proteins. The average structures from the simulations were refined against the experimental X-ray intensities. It is found that low-temperature molecular dynamics simulations provide a useful tool for optimizing the refinement of X-ray structures.     

Anthracycline-DNA interactions studied with linear dichroism and fluorescence spectroscopy

DNA-binding geometry and dynamics of a number of anthracyclines, including adriamycin and 4-demethoxydaunorubicin, interacting with DNA have been studied by means of linear dichroism and fluorescence techniques. The anthracycline chromophore is found to be approximately parallel to the plane of the DNA bases and to have a restricted mobility, as would be expected for an intercalative binding mode, but there are variations between different directions in the chromophore as well as between the drugs. From dichroic spectra of adriamycin in an anisotropic host of poly(vinyl alcohol), absorption components corresponding to transitions with mutually orthogonal polarizations have been resolved. These can be exploited to determine the orientations of the two chromophore axes in the DNA complex relative to the DNA helix axis. In a certain binding regime the long axis of the bound anthracycline chromophores (with the exception of 4-demethoxydaunorubicin) is found to be approximately 10 degrees closer to perpendicular to the helix axis than are the DNA bases. This demonstrates that the average base tilt is at least 10 degrees. By contrast, the short axis of the aglycon moiety is found to be tilted some 20-30 degrees from perpendicular. This may be because it is probing a base direction with a more pronounced, static or dynamic, inclination than the average in DNA. The drug orientation and the DNA orientation (reflecting flexibility) are observed to vary differently and nonmonotonically with binding ratio, suggesting specific binding and varying site geometries.(ABSTRACT TRUNCATED AT 250 WORDS)     

Fast and accurate fitting of relaxation dispersion data using the flexible software package GLOVE

Relaxation dispersion spectroscopy is one of the most widely used techniques for the analysis of protein dynamics. To obtain a detailed understanding of the protein function from the view point of dynamics, it is essential to fit relaxation dispersion data accurately. The grid search method is commonly used for relaxation dispersion curve fits, but it does not always find the global minimum that provides the best-fit parameter set. Also, the fitting quality does not always improve with increase of the grid size although the computational time becomes longer. This is because relaxation dispersion curve fitting suffers from a local minimum problem, which is a general problem in non-linear least squares curve fitting. Therefore, in order to fit relaxation dispersion data rapidly and accurately, we developed a new fitting program called GLOVE that minimizes global and local parameters alternately, and incorporates a Monte-Carlo minimization method that enables fitting parameters to pass through local minima with low computational cost. GLOVE also implements a random search method, which sets up initial parameter values randomly within user-defined ranges. We demonstrate here that the combined use of the three methods can find the global minimum more rapidly and more accurately than grid search alone.     

FRETmatrix: a general methodology for the simulation and analysis of FRET in nucleic acids

Förster resonance energy transfer (FRET) is a technique commonly used to unravel the structure and conformational changes of biomolecules being vital for all living organisms. Typically, FRET is performed using dyes attached externally to nucleic acids through a linker that complicates quantitative interpretation of experiments because of dye diffusion and reorientation. Here, we report a versatile, general methodology for the simulation and analysis of FRET in nucleic acids, and demonstrate its particular power for modelling FRET between probes possessing limited diffusional and rotational freedom, such as our recently developed nucleobase analogue FRET pairs (base–base FRET). These probes are positioned inside the DNA/RNA structures as a replacement for one of the natural bases, thus, providing unique control of their position and orientation and the advantage of reporting from inside sites of interest. In demonstration studies, not requiring molecular dynamics modelling, we obtain previously inaccessible insight into the orientation and nanosecond dynamics of the bases inside double-stranded DNA, and we reconstruct high resolution 3D structures of kinked DNA. The reported methodology is accompanied by a freely available software package, FRETmatrix, for the design and analysis of FRET in nucleic acid containing systems.     

Fluorescence energy transfer monitored competitive equilibria of nucleic acids: Applications in thermodynamics and screening

Precise thermodynamic characterization of nucleic acid complex stability is required to understand a variety of biologically significant events as well as to exploit the specific recognition capabilities of nucleic acids in biotechnology, diagnostics, and therapeutics. The development of a database of nucleic acid thermodynamics with sufficient precision to foster further developments in these areas requires new and improved measurement techniques. The combination of a competitive equilibrium titration with fluorescence energy transfer based detection provides a method for precise measurement of differences in free energy values for nucleic acid duplexes that far exceeds in precision those accessible via conventional methods. The method can be applied to detect and to characterize any deviation in a nucleic acid that alters duplex stability. Such deviations include, but are not limited to, mismatches; single nucleotide polymorphisms (SNP); chemically modified nucleotide bases, sugars or phosphates; and conformational anomalies or folding motifs, such as, loops or hairpins. © 2002 Wiley Periodicals, Inc. Biopoly (Nucleic Acid Sci) 61: 214–223, 2002     

Conformational and helicoidal analysis of the molecular dynamics of proteins: "curves," dials and windows for a 50 psec dynamic trajectory of BPTI

A new procedure for the graphic analysis of molecular dynamics (MD) simulations on proteins is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the conformational and helicoidal parameters for each structure entering the analysis via the method "Curves," developed for proteins by Sklenar, Etchebest, and Lavery (Proteins: Structure, Function Genet. 6:46-60, 1989) followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels ("dials") for each torsional parameter and "windows" on the range values assumed by the linear and angular helicoidal parameters, and is present in a form isomorphous with the primary structure per se. The complete time evolution of dynamic structure can then be depicted in a set of four composite figures. Dynamic aspects of secondary and tertiary structure are also provided. The procedure is illustrated with an analysis of a 50 psec in vacuo simulation on the 58 residue protein, bovine pancreatic trypsin inhibitor (BPTI), in the vicinity of the local minimum on the energy surface corresponding to a high resolution crystal structure. The time evolution of 272 conformational and 788 helicoidal parameters for BPTI is analyzed. A number of interesting features can be discerned in the analysis, including the dynamic range of conformational and helicoidal motions, the dynamic extent of 2 degrees structure motifs, and the calculated fluctuations in the helix axis. This approach is expected to be useful for a critical analysis of the effects of various assumptions about force field parameters, truncation of potentials, solvation, and electrostatic effects, and can thus contribute to the development of more reliable simulation protocols for proteins. Extensions of the analysis to present differential changes in conformational and helicoidal parameters is expected to be valuable in MD studies of protein complexes with substrates, inhibitors, and effectors and in determining the nature of structural changes in protein-protein interactions.     

Photochemical stability of aromatic amino acids under picosecond laser UV-irradiation

Two-step photochemical decomposition of aromatic amino acids under picosecond laser UV-irradiation was investigated. These results were compared with the photochemical stability of nucleic acid bases. Using the known ratio between the nucleic acid bases and aromatic amino acids in native bacteriophages lambda and phi X174 it was shown that picosecond laser UV-inactivation of viruses occurred due to the photodegradation of nucleic acid.     

Adsorption mechanism of single guanine and thymine on single-walled carbon nanotubes

Bio-nano hybrids introduce magnificent applications of nanomaterials to various fields. The choice of carbon nanotube as well as sequence selection of the nucleic acid bases play a crucial role in shaping DNA–carbon nanotube hybrids. To come up with a clear vision for the choice of carbon nanotube and nucleic acid bases to create bio-nano hybrids, we studied the adsorption mechanism of the nucleic acid bases guanine and thymine on four different types of nanotubes based on density functional theory. Nucleic acid bases exhibit differential binding strengths according to their structural geometry, inter-molecular distances, the carbon nanotube diameter, and charge transfer. The π–π interaction mechanism between the adsorbent and adsorbate is discussed in terms of charge density profile and electronic band structure analysis.     

Fluorescent probe and permeability to cells of isopoly (S-carboxymethyl-L-cysteine) derivative of nucleic acid bases.

Isopoly(S-carboxymethyl-L-cysteine) derivatives of nucleic acid bases were found to form stable complex with oligo-DNA in vitro. Fluorescent probed isopoly(S-carboxymethyl cysteine) derivatives of nucleic acid bases were prepared as antisense oligomers. The transfection of the oligomer into cells was carried out by HVJ-liposome method. Fluorescence was observed from the cells treated with HVJ-liposome including fluorescent probed oligomers.