Collective variable modelling of nucleic acids

THz-Spectroscopy of Biological Molecules

The terahertz frequency absorption spectraof DNA molecules reflect low-frequencyinternal helical vibrations involvingrigidly bound subgroups that are connectedby the weakest bonds, including thehydrogen bonds of the DNA base pairs,and/or non-bonded interactions. Althoughnumerous difficulties make the directidentification of terahertz phonon modes inbiological materials very challenging, ourresearch has shown that such measurementsare both possible and fruitful. Spectra ofdifferent DNA samples reveal a large numberof modes and a reasonable level ofsequence-specific uniqueness. In an attemptto show that the long wavelength absorptionfeatures are intrinsic properties ofbiological materials determined by phononmodes, a normal mode analysis has been usedto predict the absorption spectra ofpolynucleotide RNA Poly[G]-Poly[C]. Directcomparison demonstrated a correlationbetween calculated and experimentallyobserved spectra of the RNA polymers, thusconfirming that the fundamental physicalnature of the observed resonance structureis caused by the internal vibration modesin the macromolecules.In this work we demonstrate results fromFourier-Transform Infrared (FTIR)spectroscopy of DNA macromolecules andrelated biological materials in theterahertz frequency range. Carefulattention was paid to the possibility ofinterference or etalon effects in thesamples, and phenomena were clearlydifferentiated from the actual phononmodes. In addition, we studied thedependence of transmission spectra ofaligned DNA and polynucleotide film sampleson molecule orientation relative to theelectromagnetic field, showing the expectedchange in mode strength as a function ofsample orientation. Further, the absorptioncharacteristics were extracted from thetransmission data using the interferencespectroscopy technique, and a stronganisotropy of terahertz characteristics wasdemonstrated.

THZ-FREQUENCY SPECTROSCOPIC SENSING OF DNA AND RELATED BIOLOGICAL MATERIALS

The terahertz frequency absorption spectra of DNA molecules reflect low-frequency internal helical vibrations involving rigidly bound subgroups that are connected by the weakest bonds, including the hydrogen bonds of the DNA base pairs, and/or non-bonded interactions. Although numerous difficulties make the direct identification of terahertz phonon modes in biological materials very challenging, recent studies have shown that such measurements are both possible and useful. Spectra of different DNA samples reveal a large number of modes and a reasonable level of sequence-specific uniqueness. This chapter utilizes computational methods for normal mode analysis and theoretical spectroscopy to predict the low-frequency vibrational absorption spectra of short artificial DNA and RNA. Here the experimental technique is described in detail, including the procedure for sample preparation. Careful attention was paid to the possibility of interference or etalon effects in the samples, and phenomena were clearly differentiated from the actual phonon modes. The results from Fourier-transform infrared spectroscopy of DNA macromolecules and related biological materials in the terahertz frequency range are presented. In addition, a strong anisotropy of terahertz characteristics is demonstrated. Detailed tests of the ability of normal mode analysis to reproduce RNA vibrational spectra are also conducted. A direct comparison demonstrates a correlation between calculated and experimentally observed spectra of the RNA polymers, thus confirming that the fundamental physical nature of the observed resonance structure is caused by the internal vibration modes in the macromolecules. Application of artificial neural network analysis for recognition and discrimination between different DNA molecules is discussed.

Protein-DNA binding specificity: a grid-enabled computational approach applied to single and multiple protein assemblies

We use a physics-based approach termed ADAPT to analyse the sequence-specific interactions of three proteins which bind to DNA on the side of the minor groove. The analysis is able to estimate the binding energy for all potential sequences, overcoming the combinatorial problem via a divide-and-conquer approach which breaks the protein-DNA interface down into a series of overlapping oligomeric fragments. All possible base sequences are studied for each fragment. Energy minimisation with an all-atom representation and a conventional force field allows for conformational adaptation of the DNA and of the protein side chains for each new sequence. As a result, the analysis depends linearly on the length of the binding site and complexes as large as the nucleosome can be treated, although this requires access to grid computing facilities. The results on the three complexes studied are in good agreement with experiment. Although they all involve significant DNA deformation, it is found that this does not necessarily imply that the recognition will be dominated by the sequence-dependent mechanical properties of DNA.

Computer simulation study of probe-target hybridization in model DNA microarrays: effect of probe surface density and target concentration

We use lattice Monte Carlo simulations to study the thermodynamics of hybridization of single-stranded "target" genes in solution with complementary "probe" DNA molecules immobilized on a microarray surface. The target molecules in our system contain 48 segments and the probes tethered on a hard surface contain 8-24 segments. The segments on the probe and target are distinct, with each segment representing a sequence of nucleotides that interacts exclusively with its unique complementary target segment with a single hybridization energy; all other interactions are zero. We examine how surface density (number of probes per unit surface area) and concentration of target molecules affect the extent of hybridization. For short probe lengths, as the surface density increases, the probability of binding long stretches of target segments increases at low surface density, reaches a maximum at an intermediate surface density, and then decreases at high surface density. Furthermore, as the surface density increases, the target is less likely to bind completely to one probe; instead, it binds simultaneously to multiple probes. At short probe lengths, as the target concentration increases, the fraction of targets binding completely to the probes (specificity) decreases. At long probe lengths, varying the target concentration does not affect the specificity. At all target concentrations as the probe length increases, the fraction of target molecules bound to the probes by at least one segment (sensitivity) increases while the fraction of target molecules completely bound to the probes (specificity) decreases. This work provides general guidelines to maximizing microarray sensitivity and specificity. Our results suggest that the sensitivity and specificity can be maximized by using probes 130-180 nucleotides long at a surface density in the range of 7 x 10(-5)- 3 x 10(-4) probe molecules per nm(2).

An analysis of the electrostatic interaction between nucleic acid bases

Results from several commonly used approximate methods of evaluating electrostatic interactions have been compared to the rigorous, nonexpanded electrostatic energies at both uncorrelated and correlated levels of theory. We examined a number of energy profiles for both hydrogen bonded and stacked configurations of the nucleic acid base pairs. We found that the penetration effects play an extremely important role and the expanded electrostatic energies are significantly underestimated with respect to the ab initio values. Apart from the inability to reproduce the magnitudes of the ab initio electrostatic energy, there are other problems with the available approximate electrostatic models. For example, the distributed multipole analysis, one of the most advanced methods, is extremely sensitive to the basis set and level of theory used to evaluate the multipole moments. Detailed ab initio results are provided that other researchers could use to test their approximate models.

Computer simulation study of molecular recognition in model DNA microarrays

DNA microarrays have been widely adopted by the scientific community for a variety of applications. To improve the performance of microarrays there is a need for a fundamental understanding of the interplay between the various factors that affect microarray sensitivity and specificity. We use lattice Monte Carlo simulations to study the thermodynamics and kinetics of hybridization of single-stranded target genes in solution with complementary probe DNA molecules immobilized on a microarray surface. The target molecules in our system contain 48 segments and the probes tethered on a hard surface contain 8-24 segments. The segments on the probe and target are distinct and each segment represents a sequence of nucleotides ( approximately 11 nucleotides). Each probe segment interacts exclusively with its unique complementary target segment with a single hybridization energy; all other interactions are zero. We examine how the probe length, temperature, or hybridization energy, and the stretch along the target that the probe segments complement, affect the extent of hybridization. For systems containing single probe and single target molecules, we observe that as the probe length increases, the probability of binding all probe segments to the target decreases, implying that the specificity decreases. We observe that probes 12-16 segments ( approximately 132-176 nucleotides) long gave the highest specificity and sensitivity. This agrees with the experimental results obtained by another research group, who found an optimal probe length of 150 nucleotides. As the hybridization energy increases, the longer probes are able to bind all their segments to the target, thus improving their specificity. The hybridization kinetics reveals that the segments at the ends of the probe are most likely to start the hybridization. The segments toward the center of the probe remain bound to the target for a longer time than the segments at the ends of the probe.

Large-amplitude fast motions in double-stranded DNA driven by solvent thermal fluctuations

The nature of the internal dynamics of double-stranded DNA in aqueous environment remains to be established. We consider the motions to stem from thermal fluctuations/dissipations of the harmonic modes of beads (bases and sugars) in a cylindrical geometry that are tracked through the stochastic Langevin trajectories; these are characterized by parameters obtained from published data. The present approach has allowed a comparative study of the dynamics for DNA lengths in the range of 20-600 base pairs. For this range, we find that rotational motions about directions parallel to the helix axis (opening, twist) and perpendicular to it (propeller-twist, roll) contribute significantly to the dynamics. For a 20-mer at a solvent viscosity of 1 cP, the calculated fluorescence anisotropy profile exhibits a fast decay in the subnanosecond range due to large-amplitude fluctuations at the mesoscopic level. This feature reproduces the experimental behavior well, and suggests a possible way for the initiation of biological processes: they may be suddenly triggered on this scale through the occurrence of favorable thermal fluctuations. This analysis also reveals that, as is the case for a 20-mer, the dynamics of longer N-mers are dominated by internal motions, and are modulated by the viscosity of the solvent, in agreement with our previous experimental observations. Moreover, the model indicates that occurrence of partially concerted rotations of the bases due to thermal fluctuations can possibly be sustained over a DNA length of the order of 100 A at 1 ns, suggesting a possible mechanism for action-at-a-distance in transcription.

DNA sequence-dependent deformability?insights from computer simulations

The article reviews some recent developments in studying DNA sequence-dependent deformability, with emphasis on computer modeling. After a brief outline of available experimental techniques, we proceed to computational methods and focus on atomic-resolution molecular dynamics (MD) simulations. A sequence-dependent local (base-pair step) force field inferred from MD is compared with force fields obtained by other techniques. Various methods for establishing global (flexible-rod) DNA elastic constants are reviewed, including an approach based on atomic resolution MD. The problem of defining the global deformation variables, as well as the question of anisotropy and nonlocal effects, are discussed. As an example, both local and global deformability calculations from atomic-resolution MD of EcoRI dodecamer are presented.

Computation of DNA backbone conformations

We present an algorithm for the computation of 2'-deoxyribose-phosphodiester backbone conformations that are stereochemically compatible with a given arrangement of nucleic acid bases in a DNA structure. The algorithm involves the sequential computation of 2'-deoxyribose and phosphodiester conformers (collectively referred to as a backbone 'segment'), beginning at the 5'-end of a DNA strand. Computation of the possible segment conformations is achieved by the initial creation of a fragment library, with each fragment representing a set of bond lengths, bond angles and torsion angles. Following exhaustive searching of sugar conformations, each segment conformation is reduced to a single vector, defined by a specific distance, angle and torsion angle, that allows calculation of the O(1)' position. A given 'allowed' conformation of a backbone segment is determined based on its compatibility with the base positions and with the position of the preceding backbone segment. Initial computation of allowable segment conformations of a strand is followed by the determination of continuous backbone solutions for the strand, beginning at the 3'-end. The algorithm is also able to detect repeating segment conformations that arise in structures containing geometrically repeating dinucleotide steps. To illustrate the utility and properties of the algorithm, we have applied it to a series of experimental DNA structures. Regardless of the conformational complexity of these structures, we are able to compute backbone conformations for each structure. Hence, the algorithm, which is currently implemented within a new computer program NASDAC (Nucleic Acids: Structure, Dynamics and Conformation), should have generally applicability to the computation of DNA structures.

Theoretical methods for the simulation of nucleic acids

Different theoretical methods for the description of nucleic acid structures are reviewed. Firstly, we introduce the concept of classical force-field in the context of nucleic acid structures, discussing their accuracy. We then examine theoretical approaches to the description of nucleic acids based on: i) a rigid or quasi-rigid description of the molecule, ii) molecular mechanics optimization, and iii) molecular dynamics. Special emphasis is made ion current state of the art molecular dynamics simulations of nucleic acids structures.

ADAPT: a molecular mechanics approach for studying the structural properties of long DNA sequences

We describe an original approach to determining sequence-structure relationships for DNA. This approach, termed ADAPT, combines all-atom molecular mechanics with a multicopy algorithm to build nucleotides that contain all four standard bases in variable proportions. These nucleotides enable us to search very rapidly for base sequences that energetically favor chosen types of DNA deformation or chosen DNA-protein or DNA-ligand interactions. Sequences satisfying the chosen criteria can be found by energy minimization, combinatorial sequence searching, or genome scanning, in a manner similar to the threading approaches developed for protein structure prediction. In the latter case, we are able to analyze roughly 2000 base pairs per second. Applications of the method to DNA allomorphic transitions, DNA deformation, and specific DNA interactions are presented.

Sequence-structure relationships in DNA oligomers: a computational approach

A collective-variable model for DNA structure is used to predict the conformation of a set of 30 octamer, decamer, and dodecamer oligomers for which high-resolution crystal structures are available. The model combines an all-atom base pair representation with an empirical backbone, emphasizing the role of base stacking in fixing sequence-dependent structure. We are able to reproduce trends in roll and twist to within 5 degrees across a large database of both A- and B-DNA oligomers. A genetic algorithm approach is used to search for global minimum structures and this is augmented by a grid search to identify local minimums. We find that the number of local minimums is highly sequence dependent, with certain sequences having a set of minimums that span the entire range between canonical A- and B-DNA conformations. Although the global minimum does not always agree with the crystal structure, for 24 of the 30 oligomers, we find low-energy local minimums that match the experimental step parameters. Discrepancies throw some light on the role of crystal packing in determining the solid-state conformation of double-helical DNA.

Modeling DNA deformations

Recent developments have been made in modeling double-helical DNA at four levels of three-dimensional structure: the all-atom level, whereby an oligonucleotide duplex is surrounded by a shroud of solvent molecules; the base-pair level, with explicit backbone atoms; the mesoscopic level, that is, a few hundred base pairs, with the local duplex conformation described by knowledge-based harmonic energy functions; and the scale of several thousand nucleotides, with the duplex described as an ideal elastic rod. Predictions of the sequence-dependent bending and twisting of the double helix, as well as solvent- and force-induced B-->A and over-stretching conformational transitions, are compared with experimental data. These subtle conformational changes are critical to the functioning of the double helix, including its packaging in the close confines of the cell, the mutual fit of DNA and protein in nucleoprotein complexes, and the effective recognition of base pairs in recombination and transcription.

Nucleic acids: theory and computer simulation, Y2K

Molecular dynamics simulations on DNA and RNA that include solvent are now being performed under realistic environmental conditions of water activity and salt. Improvements to force-fields and treatments of long-range interactions have significantly increased the reliability of simulations. New studies of sequence effects, axis bending, solvation and conformational transitions have appeared.