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.     

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.     

A Robust Water Potential Parameterisation

Abstract

We compare molecular dynamics simulation results for the properties of liquid water predicted by four novel water potential models. These models are designed as a combination of parameters taken from the dedicated but brittle TIP3P water potential, and the more flexible but less accurate parameterisations such as the Dreiding and Universal force fields. We find that a hybrid of Dreiding and TIP3P delivers the best results, yielding a density, diffusion coefficient and radial distribution function in good agreement with experiment, performing in some respects even better than the dedicated reference TIP3P model. Another Dreiding based force field predicts semi-quantitative results for the water structure and dynamics while the Universal force field based models are incapable of simulating a condensed phase of water at all, continuing to expand indefinitely. These observations are useful for selecting and designing robust water force field parameterisations that can be used for general simulation purposes.     

The Application of the Reaction-Field Method to the Calculation of Dielectric Constants

Abstract

The behaviour of the popular TIP3P water model has been investigated using both molecular dynamics and Monte Carlo simulation procedures. Long-range electrostatic interactions were included through a reaction-field treatment, and the nonbonded interactions were either truncated at the cutoff distance, or smoothly scaled to zero using a switching function. The thermodynamic observables, and in particular the dipole-dipole correlation functions, are found to differ between the two simulation techniques if a rigid nonbonded cutoff is applied. However, use of a switching function gives exact agreement between the simulation methodologies. This difference is ascribed to the effect of energy pumping in the molecular dynamics simulations, and suggests that dielectric constants calculated using this simulation method with the fluctuation procedure in conjunction with a reaction field should be reappraised. Thus the Monte Carlo simulation procedure offers a number of intrinsic advantages over molecular dynamics for the calculation of dielectric constants with a reaction field. The most precise value for the dielectric constant of TIP3P is calculated to be 102 ± 3 at 298 K.     

Hydration at the TD Damaged Site of DNA and its Role in the Formation of Complex with T4 Endonuclease V

Abstract

An analysis of the distribution of water around DNA surface focusing on the role of the distribution of water molecules in the proper recognition of damaged site by repair enzyme T4 Endonuclease V was performed. The native DNA dodecamer, dodecamer with the thymine dimer (TD) and complex of DNA and part of repair enzyme T4 Endonuclease V were examined throughout the 500 ps of molecular dynamics simulation. During simulation the number of water molecules close to the DNA atoms and the residence time were calculated. There is an increase in number of water molecules lying in the close vicinity to TD if compared with those lying close to two native thymines (TT). Densely populated area with water molecules around TD is one of the factors detected by enzyme during scanning process. The residence time was found higher for molecule of the complex and the six water molecules were found occupying the stabile positions between the TD and catalytic center close to atoms P, C3′ and N3. These molecules originate water mediated hydrogen bond network that contribute to the stability of complex required for the onset of repair process.     

Effect of Warmup Protocol and Sampling Time on Convergence of Molecular Dynamics Simulations of a DNA Dodecamer Using AMBER 4.1 and Particle-Mesh Ewald Method

Abstract

This report describes one 3000 ps and two 1500 ps molecular dynamic simulations on a TATA box containing dodecamer DNA duplex in a periodic box of TIP3P water molecules, using the AMBER 4.1 implementation of the particle-mesh Ewald method. We compare the effect of warmup protocol and simulation time length on the root-mean square deviation (RMSD) parameter. For the longer simulation, the RMSD computed for the 500–1000 ps time interval is representative of longer time intervals, including 500–3000 ps. The various warmup protocols do not appear to have a significant effect on the simulation results. Based on the present results, DNA sequence-dependent differences in RMSD, or related properties, should exceed two standard deviations before being attributed to non-simulation factors, such as warmup protocol and sampling time effects; we recommend a minimum criterion of at least a three standard deviation difference with a sampling period of at least 500–1000 ps. In addition, while end effects appear negligible there is a consistent dependence of RMSD on DNA helix length.     

Monte-Carlo Simulation of DNA Duplex Hydration. B and B′ Conformations of Poly(dA) · Poly(dT) Have Different Hydration Shells

Abstract

Monte-Carlo simulation of poly(dA) · poly(dT) hydration by 30 water molecules per nucleotide pair has been performed. Two B-family conformations, both with a 36° helical twist but with different minor groove widths, were considered. One conformation is Arnott's standard B form, the other one is specific for poly(dA) · poly(dT) B′ form with a narrowed minor groove. The mean energies and the mean numbers of water-water and water-DNA hydrogen bonds are close for the two conformations. Nevertheless, the hydration shell of the B' form differs drastically from that of the standard B form. The water arrangement in the minor groove of the B′ form resembles the spine of hydration in the central part of Dickerson's dodecamer d(CGCGAATTCGCG). No such spine is formed in the hydration shell of the usual B form with a wider minor groove. In this conformation water bridges between adenine N3 or thymine O2 and oxygen of the sugar ring of the neighbouring nucleotide along the chain can be formed (“strings” in Dickerson's decamer d(CCAAGATTGG)).     

An Unexpected Major Groove Binding of Netropsin and Distamycin A to tRNAphe

Abstract

Crystalline complexes of yeast tRNA^phe and the oligopeptide antibiotics netropsin and distamycin A were prepared by diffusing drugs into crystals of tRNA. X-ray structure analyses of these complexes reveal a single common binding site for both drugs which is located in the major or deep groove of the tRNA T-stem. The netropsin-tRNA complex is stabilized by specific hydrogen bonds between the amide groups of the drug and the tRNA bases G51 0(6), U52 0(4) and G53 N(7) on one strand, and is further stabilized by electrostatic interactions between the positively charges guanidino side chain of the drug and the tRNA phosphate P53 on the same strand and the positively charged amidino propyl side chain and the phosphates P61, P62 and P63 on the opposite strand of the double helix. These results are in contrast to the implicated minor groove binding of these drugs to non-guanine sequences in DNA. The binding to the GUG sequence in tRNA implies that major groove binding to certain DNA sequences is possible.     

Monte Carlo Study of Structural and Thermodynamic Properties of Liquid Chloroform Using a Five Site Model

A five site potential model combining Lennard–Jones plus Coulomb potential functions has been developed for chloroform molecule. The partial charges needed for Coulombic interactions were derived using the chelpg procedure implemented in the gaussian 92 program. These calculations were performed at the MP2 level with MC-311G* basis set for Cl and 6-311G** for C and H atoms. The parameters for the Lennard–Jones potentials were optimized to reproduce experimental values for the density and enthalpy of vaporization of the pure liquid at 298 K and 1 atm. The statistical mechanics calculations were performed with the Monte Carlo method in the isothermic and isobaric (NpT) ensemble. Besides the values obtained for density, ρ, and molar enthalpy of vaporization at constant pressure, Δ H_V, for liquid chloroform, results for molar volume, V_m, molar heat capacity, C_p, isobaric thermal expansivity, α_p, and isothermal compressibility, κ_T, for this pure liquid are also in very good agreement with experimental observations. Size effects on the values of thermodynamic properties were investigated. The potential model was also tested by computing the free energy for solvating one chloroform molecule into its own liquid at 298 K using a statistical perturbation approach. The result obtained compares well with the experimental value. Site–site pair correlation functions were calculated and are in good accordance with theoretical results available in the literature. Dipole–dipole correlation functions for the present five site model were also calculated at different carbon–carbon distances. These correlations were compared to those obtained using the four site model reported in the literature. An investigation of the solvent dependence of the relative free energy for cis/trans conversion of a hypothetical solute in TIP4P water and chloroform was accomplished. The results show strong interaction of water and chloroform molecules with the gauche conformer. The value obtained for the free energy barrier for cis/trans rotation in TIP4P water is higher than that for chloroform. This result is in agreement with the continuous theory for solvation as the conformer with higher dipole moment is more favoured by the solvent with higher dieletric constant. The results also show an increase in entropy as the solute goes from the cis to the trans geometry and this result is more appreciable in the aqueous solution. This revised version was published online in June 2006 with corrections to the Cover Date.     

Computer Simulation of the Dielectric Properties of Liquid Water

Abstract

Monte Carlo simulations of water in the NVT ensemble using three models (SPC, TIP4P and TIPS2) are reported. The internal energy, dielectric constant, and the site-site radial distribution functions of liquid water (temperature 300 K and mass density 1 gm cc^?1) were calculated and compared with experiment. It was found that of the three intermolecular potential models, SPC gives the best dielectric constant. Since SPC also yields acceptable results for the energy and structure, it is judged to be the best among the three models studied.     

Simulation of protein diffusion: a sensitive probe of protein–solvent interactions

Aqueous solutions of Candida antarctica lipase B (CALB) were simulated considering three different water models (SPC/E, TIP3P, TIP4P) by a series of molecular dynamics (MD) simulations of three different box sizes (L = 9, 14, and 19 nm) to determine the diffusion coefficient, the water viscosity and the protein density. The protein–water systems were equilibrated for 500 ns, followed by 100 ns production runs which were analysed. The diffusional properties of CALB were characterized by the Stokes radius (R_S), which was derived from the diffusion coefficient and the viscosity. R_S was compared to the geometric radius (R_G) of CALB, which was derived from the protein density. R_S and R_G differed by 0.27 nm for SPC/E and by 0.40 and 0.39 nm for TIP3P and TIP4P, respectively, which characterizes the thickness of the diffusive hydration layer on the protein surface. The simulated hydration layer of CALB resulted in agreement with those experimentally determined for other seven different proteins of comparable size. By avoiding the most common pitfalls, protein diffusion can be reliably simulated: simulating different box sizes to account for the finite size effect, equilibrating the protein–water system sufficiently, and using the complete production run for the determination of the diffusion coefficient.     

Theoretical Calculation of the Liquid—Vapor Coexistence Curve of Water,Chloroform and Methanol with the Cavity-Biased Monte Carlo Method in the Gibbs Ensemble

Abstract

The Gibbs ensemble computer simulation method of Panagiotopoulos is combined by the cavity-biased sampling technique used previously in the grand-canonical ensemble. The combined technique is applied to the determination of the liquid—vapor coexistence curve of the Lennard—Jones fluid as a test case, two water models (SPC and TIP4P) as well as methanol and chloroform, both described with the OPLS model. The application of the virial-based sampling technique, used earlier in the isobaric ensemble is also discussed.     

Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG

The crystal structure of the DNA dodecamer duplex CATGGGCCCATG lies on a structural continuum along the transition between A- and B-DNA. The dodecamer possesses the normal vector plot and inclination values typical of B-DNA, but has the crystal packing, helical twist, groove width, sugar pucker, slide and x-displacement values typical of A-DNA. The structure shows highly ordered water structures, such as a double spine of water molecules against each side of the major groove, stabilizing the GC base pairs in an A-like conformation. The different hydration of GC and AT base pairs provides a physical basis for solvent-dependent facilitation of the A↔B helix transition by GC base pairs. Crystal structures of CATGGGCCCATG and other A/B-DNA intermediates support a ‘slide first, roll later’ mechanism for the B→A helix transition. In the distribution of helical parameters in protein–DNA crystal structures, GpG base steps show A-like properties, reflecting their innate predisposition for the A conformation.     

Molecular dynamics study of substance P peptides in a biphasic membrane mimic

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in a TIP3P water/CCl4 biphasic solvent system as a mimic for the water-membrane system. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) in the equilibration stage of the simulation. The starting orientation/position of the peptides for the MD simulation was either parallel to the water/CCl4 interface or in a perpendicular/insertion mode. In both cases the peptides equilibrated and adopted a near-parallel orientation within approximately 250 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained in the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbone of nearly all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions with the hydrophobic phase come from the hydrophobic interactions of the side chains of Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, and Phe-7, Leu-10, Met-11 and, to a lesser extent, Tyr-8 in SP-Y8. Concerted conformational transitions took place in the time frame of hundreds of picoseconds. The concertedness of the transition was due to the tendency of the peptide to maintain the necessary secondary structure to position the peptide properly with respect to the water/CCl4 interface.     

Exocyclic groups in the minor groove influence the backbone conformation of DNA

Exocyclic groups in the minor groove of DNA modulate the affinity and positioning of nucleic acids to the histone protein. The addition of exocyclic groups decreases the formation of this protein–DNA complex, while their removal increases nucleosome formation. On the other hand, recent theoretical results show a strong correlation between the B_I/B_II phosphate backbone conformation and the hydration of the grooves of the DNA. We performed a simulation of the d(CGCGAATTCGCG)₂ Drew Dickerson dodecamer and one simulation of the d(CGCIAATTCGCG)₂ dodecamer in order to investigate the influence of the exocyclic amino group of guanine. The removal of the amino group introduces a higher intrinsic flexibility to DNA supporting the suggestions that make the enhanced flexibility responsible for the enlarged histone complexation affinity. This effect is attributed to changes in the destacking interactions of both strands of the DNA. The differences in the hydration of the minor groove could be the explanation of this flexibility. The changed hydration of the minor groove also leads to a different B_I/B_II substate pattern. Due to the fact that the histone preferentially builds contacts with the backbone of the DNA, we propose an influence of these B_I/B_II changes on the nucleosome formation process. Thus, we provide an additional explanation for the enhanced affinity to the histone due to removal of exocyclic groups. In terms of B_I/B_II we are also able to explain how minor groove binding ligands could affect the nucleosome assembly without disrupting the structure of DNA.     

Energetics of water proton configurations in gas hydrates: comparison of various water models

The total interaction energies for a large number of water proton configurations in the unit cell of hydrate structure I consisting of 46 molecules are compared for qualitatively different water models, such as SPC/E, TIP4P, TIP5P, TIP 3f and AMOEBA. All calculations were carried out using the TINKER molecular modelling package. The Ewald summation method with metallic tin-foil boundary conditions is used to account for long-range electrostatic interactions. It was established that there is a high correlation between the energies calculated using the five water models (interaction potentials). The average correlation coefficient for all pairs of potentials is equal to 0.91. Analogous calculations were carried out to evaluate the consistency of the different water models with respect to a new property of the ice-like system: the hydrogen-bond-reversal asymmetry. It was established that, for all water models, there is relatively high correlation between the energy differences for proton configurations with opposite direction of all hydrogen bonds. In this case, the average correlation coefficient is 0.77. Data for the TIP4P potential differ noticeably from the others, especially owing to the variation in the total interaction energy. The validity and usefulness of simple discrete models of inter-molecular interactions are discussed.     

Alteration of the groove width of DNA induced by the multimodal hydrogen bonding of denaturants with DNA bases in its grooves affects their stability

BackgroundDenaturants, namely, urea and guanidinium chloride (GdmCl) affect the stability as well as structure of DNA. Critical assessment of the role of hydrogen bonding of these denaturants with the different regions of DNA is essential in terms of its stability and structural aspect. However, the understanding of the mechanistic aspects of structural change of DNA induced by the denaturants is not yet well understood.MethodsIn this study, various spectroscopic along with molecular dynamics (MD) simulation techniques were employed to understand the role of hydrogen bonding of these denaturants with DNA bases in their stability and structural change.Results and conclusionIt has been found that both, GdmCl and urea intrude into groove region of DNA by striping surrounding water. The hydrogen bonding pattern of Gdm⁺ and urea with DNA bases in its groove region is multimodal and distinctly different from each other. The interaction of GdmCl with DNA is stabilized by electrostatic interaction whereas electrostatic and Lennard-Jones interactions both contribute for urea. Gdm⁺ forms direct hydrogen bond with the bases in the minor groove of DNA whereas direct and water assisted hydrogen bond takes place with urea. The hydrogen bond formed between Gdm⁺ with bases in the groove region of DNA is stronger than urea due to strong electrostatic interaction along with less self-aggregation of Gdm⁺ than urea. The distinct hydrogen bonding capability of Gdm⁺ and urea with DNA bases in its groove region affects its width differently. The interaction of Gdm⁺ decreases the width of the minor and major groove which probably increases the strength of hydrogen bond between the Watson-Crick base pairs of DNA leading to its stability. In contrast, the interaction of urea does not affect much to the width of the grooves except the marginal increase in the minor groove width which probably decreases the strength of hydrogen bond between Watson Crick base pairs leading to the destabilization of DNA.General significanceOur study clearly depicts the role of hydrogen bonding between DNA bases and denaturants in their stability and structural change which can be used further for designing of the guanidinium based drug molecules.     

Structure of Poly(dA)·Poly(dT) is not Identical to the AT Rich Regions of the Single Crystal Structure of CGCGAATTBrCGCG. The Consequence of this to Netropsin Binding to Poly(dA)·Poly(dT)

Abstract

The basic assumption of Dickerson and Kopka (J. Biomole. Str. Dyns. 2, 423, 1985) that the conformation of poly(dA)·poly(dT) in solution is identical to the AT rich region of the single crystal structure of the Dickerson dodecamer is not supported by any experimental data. In poly(dA)·poly(dT), NOE and Raman studies indicate that the dA and dT units are conformationally equivalent and display the (anti-S-type sugar)-conformation; incorporation of this nucleotide geometry into a double helix leads to a conventional regular B-helix in which the width of the minor groove is 8A. The derived structure is consistent with all available experimental data on poly(dA)·poly(dT) obtained under solution conditions. In the crystal structure of the dodecamer, the dA and dT units have distinctly different conformations—dA residues adopt (anti, S-type sugar pucker), while dT residues belong to (low anti, N-type sugar pucker). These different conformations of the dA and dT units along with the large propeller twist can be accommodated in a double helix in which the minor groove is shrunk from 8A to less than 4A. In the conventional right handed B-form of poly(dA)·poly(dT) with the 8A wide minor groove, netropsin has to bind asymmetrically along the dA strand to account for the NOE and chemical shift data and to generate a stereochemically sound structure (Sarma et al, J. Biomole. Str. Dyns. 2, 1085, 1985).     

上调OPN和TIP30表达对口腔鳞癌cal-27细胞的影响

目的:探讨上调骨桥蛋白(OPN)和TIP30表达对口腔鳞状细胞癌cal-27细胞增殖和迁移能力的影响。方法:采用过表达OPN载体病毒和过表达TIP30质粒分别感染或转染cal-27细胞株后,通过实时荧光定量PCR法和蛋白印迹western blot法检测OPN和TIP30的表达,CCK-8法检测细胞增殖能力,细胞划痕实验检测细胞的迁移能力。结果:过表达OPN载体病毒感染的cal-27-sh OPN细胞株OPN m RNA和蛋白的表达均明显高于cal-27-sh N细胞株(P0.001,P=0.002),且细胞的增殖能力和迁移能力强(P0.001,P=0.002)。过表达TIP30质粒转染的cal-27-p TIP30细胞株TIP30 m RNA和蛋白的表达均明显高于cal-27-pc3.0细胞株(P=0.001,P=0.003),且细胞的增殖能力和迁移能力弱(P0.001,P=0.005)。结论:OPN能促进cal-27细胞的增殖和迁移,TIP30抑制cal-27细胞的增殖和迁移。