Search PubMed⌕ Search

Biomedical subjects

D L Beveridge

Publications and source records attributed to D L Beveridge.

52 records · Page 3Linked to original sources

Theoretical account of the 'spine of hydration' in the minor groove of duplex d(CGCGAATTCGCG).

The hydration of B-form DNA involves a novel sequence of ordered water molecules in the minor groove known as the spine of hydration. The positions of the ordered waters in the AATT region have been determined and indicate water bridges between the A-N3 and T-O2 acceptor sites. Theoretical calculation on hydration based on Monte Carlo simulation show that the water bridges are more sensitive to the prescription for intermolecular potentials than the choice of canonical vs. crystallographic geometry. The GROMOS force field and the SPC model for water are shown to give an accurate theoretical account of the experimental data.

Chemical Phenomena↗

A theoretical study of the aqueous hydration of canonical B d(CGCGAATTCGCG): Monte Carlo simulation and comparison with crystallographic ordered water sites.

Monte Carlo computer simulation is described for the dodecamer d(CGCGAATTCGCG) together with 1777 water molecules at an environmental density of 1 gm/cc in a cubic cell under periodic boundary conditions. Water-water interactions were treated using the TIP4P potential and the solute water interactions by TIP4P spliced with the non-bonded interactions from the AMBER 3.0 force field. The stimulation was subjected to proximity analysis to obtain solute coordination numbers and pair interaction energies for each solute atom. Hydration density distributions partitioned into contributions from the major groove side, the minor groove side and the sugar-phosphate backbone were examined, and the probabilities of occurence for one- and two-water bridges in the simulation were enumerated. The results were compared with observations of crystallographic ordered water sites from x-ray diffraction studies on the native dodecamer by Dickerson and coworkers.

Base Sequence↗

Conformational and helicoidal analysis of 30 PS of molecular dynamics on the d(CGCGAATTCGCG) double helix: "curves", dials and windows.

A new procedure for the analysis of the structure and molecular dynamics of duplex DNA is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the values of the conformational and helicoidal parameters for each structure entering the analysis using the method "Curves" developed by Lavery and Sklenar, J. Biomol. Str. Dyn. 6, 63 (1988), followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels, or "dials", for each IUPAC torsional parameter and "windows" on the range of values assumed by the linear and angular helicoidal parameters, and is presented in a form isomorphous with the structure per se. The complete time evolution of the conformational and helicoidal parameters of a DNA double helix can then be depicted in a set of six composite figures. Dynamical aspects of helix bending are also subsumed in this analysis. The procedure is illustrated with an analysis of the structures of canonical A and B forms of DNA and the 300 degrees K native dodecamer duplex d(CGCGAATTCGCG). The "dials and windows" are then used for a comprehensive analysis of 30 psec of molecular dynamics on the dodecamer in the vicinity of a canonical B-DNA energy minimum. This involves presentation of the time evolution of 206 conformational and 230 helicoidal parameters for the dodecamer. A number of interesting structural features can be recognized in the analysis, including crankshaft motions, BI - BII transitions, sugar repuckerings, and a description of spontaneous helix bending at what corresponds to the 1 degrees and 2 degrees "hinge points" indicated in the crystal structure. Our approach is expected to be directly useful for critical analysis of the effects of various assumptions about force field parameters, hydration and electrostatic effects and thus contribute to the development of reliable simulation protocols for nucleic acid systems. Extension of the method to present differential changes in conformational and helicoidal parameters is expected to be valuable for the analysis of structural and molecular dynamics studies of the reorganization and adaptation of DNA on complexation with various drugs and regulatory proteins.

Algorithms↗

Theoretical considerations on the "spine of hydration" in the minor groove of d(CGCGAATTCGCG).d(GCGCTTAAGCGC): Monte Carlo computer simulation.

A theoretical description of aqueous hydration in the minor groove of a B-form DNA is presented on the basis of a liquid-state Monte Carlo computer simulation on a system consisting of the oligonucleotide duplex d(CGCGAATTCGCG).d(GCGCTTAAGCGC) in a canonical B-form together with 1777 water molecules contained in a hexagonal prism cell and treated under periodic boundary conditions. The results are analyzed in terms of solvent density distributions. The calculated minor-groove solvent density shows considerable localization, indicative of discrete solvation sites and providing theoretical evidence for a well-defined ordered water structure. In the AATT sequence, this corresponds to the "spine of hydration" described by H. R. Drew and R. E. Dickerson [(1981) J. Mol. Biol. 151, 535-556] based on the x-ray crystal structure of the dodecamer hydrate. We find, however, that the calculated ordered water structure also extends into the CGCG flanking sequences, supported by the N2 hydrogen bond donors of the guanine residues and indicating that the spine of hydration could thus extend throughout the minor groove of a B-form DNA. This provides a possible explanation of the positive binding entropies observed by L. A. Marky and K. J. Breslauer [(1984) Proc. Natl. Acad. Sci. USA 84, 4359-4363] for both A.T and C.G sequences on the complexation of netropsin to the minor groove of DNAs. Implications of these results with regard to the thermodynamic stability of DNA in water and the sequence specificity of the minor groove hydration are discussed.

Computer Simulation↗

Monte Carlo computer simulation of the aqueous hydration of the glycine zwitterion at 25 degree C.

Monte Carlo computer simulation on a dilute aqueous solution of the glycine zwitterion are reported. The results are presented in terms of the Quasi-Component Distribution Functions (QCDF) of Ben Naim and partitioned into atomic and functional group contributions using the Proximity Criterion. The Proximity Criterion analysis has been extended to orientational properties and a new normalization procedure has been introduced for the radial distribution functions obtained by the Proximity Criterion. The solvation environment of the glycine zwitterion is found to contain, on the average, 14.4 water molecules out of which 3.2 belong to the ammonium group, 6.1 to the methylene group and 5.1 to the carboxyl group. The importance of the many-body statistical mechanical approach to hydration is emphasized by our finding that the configuration corresponding to the absolute minimum of the glycine zwitterion-water potential surface was found to have negligible statistical weight in the aqueous simulation.

Computer Simulation↗

Aqueous hydration of nucleic acid constituents: Monte Carlo computer simulation studies.

Monte Carlo computer simulations were performed on dilute aqueous solutions of thymine, cytosine, uracil, adenine, guanine, the dimethyl phosphate anion in the gauche-gauche conformation and a ribose and deoxyribose derivative. The aqueous hydration of each molecule was analysed in terms of quasi-component distribution functions based on the Proximity Criterion, and partitioned into hydrophobic, hydrophilic and ionic contributions. Color stereo views of selected hydration complexes are also presented. A preliminary discussion of the transferability of functional group coordination numbers is given. The results enable to comment on two current problems related to the hydration of nucleic acids: a) the theory of Dickerson and coworkers on the role of water in the relative stability of the A and B form of DNA and b) the idea of water bridges and filaments emerging from the computer simulation results on the hydration of DNA fragments by Clementi.

Carbohydrates↗

Monte Carlo studies on water in the dCpG/proflavin crystal hydrate.

The extensive water network identified in the crystallographic studies of the dCpG/Proflavin hydrate by Neidle, Berman and Shieh (Nature 288, 129, 1980) forms an ideal test case for a) assessing the accuracy of theoretical calculations on nucleic acid--water systems based on statistical thermodynamic computer simulation, and b) the possible use of computer simulation in predicting the water positions in crystal hydrates for use in the further refinement and interpretation of diffraction data. Monte Carlo studies have been carried out on water molecules in the unit cell of dCpG/proflavin, with the nucleic acid complex fixed and the condensed phase environment of the system treated by means of periodic boundary conditions. Intermolecular interactions are described by potential functions representative of quantum mechanical calculations developed by Clementi and coworkers, and widely used in recent studies of the aqueous hydration of various forms of DNA fragments. The results are analyzed in terms of hydrogen bond topology, hydrogen bond distances and energies, mean water positions, and water crystal probability density maps. Detailed comparison of calculated and experimentally observed results are given, and the sensitivity of results to choice of potential is determined by comparison with simulation results based on a set of empirical potentials.

Acridines↗

A Monte Carlo simulation study of the aqueous hydration of r(GpC)2: comparison with crystallographic ordered water sites.

Monte Carlo computer simulation is described for the dinucleotide duplex rGpC together with 562 water molecules at an environmental density of 1 g/cc in a cubic cell under periodic boundary conditions. Water-water interactions were treated using the TIP4P potential and the solute water interactions by TIP4P spliced with the nonbonded interactions from the AMBER 3.0 force field. The simulation was subjected to proximity analysis to obtain solute coordinate numbers and pair interaction energies for each solute atom. Hydration density distributions partitioned into contributions from the major groove side, the minor groove side and the sugar-phosphate backbone were examined, and the probabilities of occurrence for one- and two-water bridges in the simulation were enumerated. The results were compared with observations of crystallographic ordered water sites from x-ray diffraction studies on G and C containing small molecules, and in crystal structure determinations of the sodium, calcium, and ammonium salts of rGpC. The calculated results are generally consistent with the observed sites, except for cytosine N4, where a hydration site is predicted yet none observed in rGpC salts, and for guanine N3, which appears in this calculation to compete unfavorably with the adjacent donor site at guanine N2. There is, however, a significant probability of finding a one-water G-N3-W-G-N2 bridge indicated in the simulation. An explanation for the guanine N3 discrepancy in terms of electrostatic potentials is also offered. The calculated one- and two-water bridges in the rGpC hydration complex coincide in a number of cases to those observed in the ordered water structure of the sodium rGpC crystal hydrate.

Computer Simulation↗