Synergistic effects in the melting of DNA hydration shell: melting of the minor groove hydration spine in poly(dA).poly(dT) and its effect on base pair stability

Specific and highly efficient condensation of GC and IC DNA by polyaza pyridinophane derivatives

Two bis-polyaza pyridinophane derivatives and their monomeric reference compounds revealed strong interactions with ds-DNA and RNA. The bis-derivatives show a specific condensation of GC- and IC-DNA, which is almost two orders of magnitude more efficient than the well-known condensation agent spermine. The type of condensed DNA was identified as ψ-DNA, characterized by the exceptionally strong CD signals. At variance to the almost silent AT(U) polynucleotides, these strong CD signals allow the determination of GC-condensates at nanomolar nucleobase concentrations. Detailed thermodynamic characterisation by ITC reveals significant differences between the DNA binding of the bis-derivative compounds (enthalpy driven) and that of spermine and of their monomeric counterparts (entropy driven). Atomic force microscopy confirmed GC-DNA compaction by the bis-derivatives and the formation of toroid- and rod-like structures responsible for the ψ-type pattern in the CD spectra.

Investigation of changes in the arrangement of water molecules and salt ions surrounding different atoms of the DNA molecule during the melting process: a molecular dynamics simulation study

A molecular dynamics simulation was performed on a B-DNA duplex (CGCGAATTGCGC) at different temperatures. The DNA was immerged in a saltwater medium with 1 mol/L NaCl concentration. The arrangements of water molecules and cations around the different atoms of DNA on the melting pathway were investigated. Almost for all atoms of the DNA by double helix → single-stranded transition, the water molecules released from the DNA duplex and cations were close to single-stranded DNA, but this behavior was not clearly seen at melting temperatures. Therefore, release of water molecules and cations approaching the DNA by the increase of temperature does not have any effect on the sharpness of the transition curve. Most of the water molecules and cations were found to be around the negatively charged phosphate oxygen atoms. The number of water molecules released from the first shell hydration upon melting in the minor groove was higher than in the major groove, and intrusion of cations into the minor groove after melting was higher than into the major groove. The hydrations of imino protons were different from each other and were dependent on DNA bases.

Studies of base pair sequence effects on DNA solvation based on all-atom molecular dynamics simulations

Detailed analyses of the sequence-dependent solvation and ion atmosphere of DNA are presented based on molecular dynamics (MD) simulations on all the 136 unique tetranucleotide steps obtained by the ABC consortium using the AMBER suite of programs. Significant sequence effects on solvation and ion localization were observed in these simulations. The results were compared to essentially all known experimental data on the subject. Proximity analysis was employed to highlight the sequence dependent differences in solvation and ion localization properties in the grooves of DNA. Comparison of the MD-calculated DNA structure with canonical A- and B-forms supports the idea that the G/C-rich sequences are closer to canonical A- than B-form structures, while the reverse is true for the poly A sequences, with the exception of the alternating ATAT sequence. Analysis of hydration density maps reveals that the flexibility of solute molecule has a significant effect on the nature of observed hydration. Energetic analysis of solute-solvent interactions based on proximity analysis of solvent reveals that the GC or CG base pairs interact more strongly with water molecules in the minor groove of DNA that the AT or TA base pairs, while the interactions of the AT or TA pairs in the major groove are stronger than those of the GC or CG pairs. Computation of solvent-accessible surface area of the nucleotide units in the simulated trajectories reveals that the similarity with results derived from analysis of a database of crystallographic structures is excellent. The MD trajectories tend to follow Manning's counterion condensation theory, presenting a region of condensed counterions within a radius of about 17 A from the DNA surface independent of sequence. The GC and CG pairs tend to associate with cations in the major groove of the DNA structure to a greater extent than the AT and TA pairs. Cation association is more frequent in the minor groove of AT than the GC pairs. In general, the observed water and ion atmosphere around the DNA sequences is the MD simulation is in good agreement with experimental observations.

Collagen-DNA complex

Previously presented models of collagen-DNA (7) and collagen-siRNA (8) complexes point to a general description of delivery systems and indicate to what specific topology that system should be equipped with to effectively deliver the gene into the living body via in vivo and in vitro injection. We focused our interest on the nature of collagen-DNA complex structure and the molecular and environmental determinants of the self-association processes of collagen with the presence of DNA. In this aspect, the self-association of collagen-DNA complex offers an opportunity to characterize a unique system, which may be related to the general mechanisms of self-association of fiber macromolecules by water bridges. For characterizing the collagen-DNA interaction, we used FTIR-ATR, NMR, and AFM experiments done on a separate collagen film, DNA film, and on the peptide-DNA aqueous solution. We demonstrate that collagen-DNA spontaneously forms self-assembling complex systems in aqueous solution. Such self-association of the complex could be induced by electrostatic interactions between neutral collagen cylinders, having strong dipole moment, and negatively charged DNA cylinders. A final complex could be formed by hydrogen bonds between specified donor groups of collagen and phosphate acceptor groups of DNA. According to FTIR measurements, a collagen triple helix should not change global conformation during collagen-DNA complex formation.

On the origin of the decrease in stability of the DNA hairpin d(GCGAAGC) on complexation with aromatic drugs

Molecular dynamics simulations of drug-DNA complexes have been carried out in order to explain the experimentally observed decrease in thermal stability of the DNA hairpin d(GCGAAGC) on binding the aromatic drug molecules, daunomycin, ethidium bromide, novantrone and proflavine. This complexation behavior is in contrast to the stabilizing effect of the same aromatic drug molecules on DNA duplexes. Analysis of the energy parameters and the hydration properties of the complexes shows that the main factor correlating with the decrease in melting temperatures of the drug-hairpin complexes is the number of water bridges, with a reduction of at least 40% on ligand binding.

Effects of high hydrostatic pressures on living cells: a consequence of the properties of macromolecules and macromolecule-associated water

Sixty percent of the Earth's biomass is found in the sea, at depths greater than 1000 m, i.e., at hydrostatic pressures higher than 100 atm. Still more surprising is the fact that living cells can reversibly withstand pressure shifts of 1000 atm. One explanation lies in the properties of cellular water. Water forms a very thin film around macromolecules, with a heterogeneous structure that is an image of the heterogeneity of the macromolecular surface. The density of water in contact with macromolecules reflects the physical properties of their different domains. Therefore, any macromolecular shape variations involving the reorganization of water and concomitant density changes are sensitive to pressure (Le Chatelier's principle). Most of the pressure-induced changes to macromolecules are reversible up to 2000 atm. Both the effects of pressure shifts on living cells and the characteristics of pressure-adapted species are opening new perspectives on fundamental problems such as regulation and adaptation.

Molecular dynamics simulation study of DNA dodecamer d(CGCGAATTCGCG) in solution: conformation and hydration

A molecular dynamics simulation of the dodecamer duplex d(CGCGAATTCGCG) using the particle mesh Ewald sum assumed a B-conformation remarkably close to the observed X-ray structure. The Ewald summation method effectively eliminates the usual "cut-off" of long-range interactions and allowed us to evaluate the full effect of the electrostatic forces. This simulation showed remarkable agreement with the Dickerson X-ray structure in both average structure and B-factors; within the EcoRI site itself, the rms deviation between the average theoretical and observed structures was 1.1 A. The width of the minor groove fluctuated between a wide and narrow configuration with the latter corresponding closely to the X-ray structure. The simulation also suggested a strong sequence-dependent signature on the minor groove width in both wide and narrow conformers. Hydration shells in both the major and minor grooves were observed. The "spine of hydration" in the minor groove was clear. In the major groove the first hydration shell appears to be a ribbon-like structure that reproduces the principal features of observed X-ray structures; subtle variations of this hydration pattern suggest sequence dependencies. Sequence-dependent features were also examined for helical and other geometric parameters. The successful reproduction of many experimentally observed fine structural features shows that the Ewald summation significantly improves the fidelity of the calculations.

Temperature-dependent ultraviolet resonance Raman spectroscopy of the premelting state of dA.dT DNA

Poly(dA).poly(dT) and DNA duplex with four or more adenine bases in a row exhibits a broad, solid-state structural premelting transition at about 35 degrees C. The low-temperature structure is correlated with the phenomena of "bent DNA." We have conducted temperature-dependent ultraviolet resonance Raman measurements of the structural transition using poly(dA).poly(dT) at physiological salt conditions, and are able to identify, between the high and low temperature limits, changes in the vibrational frequencies associated with the C4 carbonyl stretching mode in the thymine ring and the N6 scissors mode of the amine in the adenine ring of poly(dA).poly(dT). This work supports the model that the oligo-dA tracts' solid-state structural premelting transition is due to a set of cross-stand bifurcated hydrogen bonds between consecutive dA. dT pairs.

Hydration of the DNA bases is local

Ordered hydration sites were determined for the nucleotide bases in B-type conformations using the crystal structure data on 14 B-DNA decamer structures. A method of density representation was extended so that positions, occupancies, and distributions of the hydration sites were predicted around a B-DNA double helix by a method analogous to crystallographic refinement. The predicted hydration sites correctly reproduce the main features of hydration around the B-DNA dodecamer. In contrast to the previous observations, the newly available crystal data show the same extent of hydration of guanine and adenine, and of cytosine and thymine.

Hydrostatic and osmotic pressure as tools to study macromolecular recognition

Clearly, hydrostatic and osmotic pressure techniques offer unique potential in the study of fundamental problems of molecular recognition in biological systems. With the recent advances in technology such investigations are rapidly becoming commonplace. We look forward to further advances and their report in succeeding compendiums such as this volume.

Premelting base pair opening probability and drug binding constant of a daunomycin-poly d(GCAT).poly d(ATGC) complex

We calculate room temperature thermal fluctuational base pair opening probability of a daunomycin-poly d(GCAT).poly d(ATGC) complex. This system is constructed at an atomic level of detail based on x-ray analysis of a crystal structure. The base pair opening probabilities are calculated from a modified self-consistent phonon approach of anharmonic lattice dynamics theory. We find that daunomycin binding substantially enhances the thermal stability of one of the base pairs adjacent the drug because of strong hydrogen bonding between the drug and the base. The possible effect of this enhanced stability on the drug inhibition of DNA transcription and replication is discussed. We also calculate the probability of drug dissociation from the helix based on the selfconsistent calculation of the probability of the disruption of drug-base H-bonds and the unstacking probability of the drug. The calculations can be used to determine the equilibrium drug binding constant which is found to be in good agreement with observations on similar daunomycin-DNA systems.