Open-Boundary Molecular Dynamics of a DNA Molecule in a Hybrid Explicit/Implicit Salt Solution

The composition and electrolyte concentration of the aqueous bathing environment have important consequences for many biological processes and can profoundly affect the behavior of biomolecules. Nevertheless, because of computational limitations, many molecular simulations of biophysical systems can be performed only at specific ionic conditions: either at nominally zero salt concentration, i.e., including only counterions enforcing the system’s electroneutrality, or at excessive salt concentrations. Here, we introduce an efficient molecular dynamics simulation approach for an atomistic DNA molecule at realistic physiological ionic conditions. The simulations are performed by employing the open-boundary molecular dynamics method that allows for simulation of open systems that can exchange mass and linear momentum with the environment. In our open-boundary molecular dynamics approach, the computational burden is drastically alleviated by embedding the DNA molecule in a mixed explicit/implicit salt-bathing solution. In the explicit domain, the water molecules and ions are both overtly present in the system, whereas in the implicit water domain, only the ions are explicitly present and the water is described as a continuous dielectric medium. Water molecules are inserted and deleted into/from the system in the intermediate buffer domain that acts as a water reservoir to the explicit domain, with both water molecules and ions free to enter or leave the explicit domain. Our approach is general and allows for efficient molecular simulations of biomolecules solvated in bathing salt solutions at any ionic strength condition.     

Targeting Telomeres: Molecular Dynamics and Free Energy Simulation of Gold-Carbene Binding to DNA

DNA sequences in regulatory regions and in telomers at the ends of chromosomes frequently contain tandem repeats of guanine nucleotides that can form stacked structures stabilized by Hoogsten pairing and centrally bound monovalent cations. The replication and elongation of telomeres requires the disruption of these G-quadruplex structures. Hence, drug molecules such as gold (Au)-carbene that stabilize G-quadruplexes may also interfere with the elongation of telomeres and, in turn, could be used to control cell replication and growth. To better understand the molecular mechanism of Au-carbene binding to G-quadruplexes, we employed molecular dynamics simulations and free energy simulations. Whereas very restricted mobility of two Au-carbene ligands was found upon binding as a doublet to one side of the G-quadruplex, much larger translational and orientational mobility was observed for a single Au-carbene binding at the second G-quadruplex surface. Comparative simulations on duplex DNA in the presence of Au-carbene ligands indicates a preference for the minor groove and weaker unspecific and more salt-dependent binding than to the G-quadruplex surface. Analysis of energetic contributions reveals a dominance of nonpolar and van der Waals interactions to drive binding. The simulations can also be helpful for proposing possible modifications that could improve Au-carbene affinity and specificity for G-quadruplex binding.     

Molecular Mechanics and Dynamics of an Abasic Frameshift in DNA and Comparison to NMR Data

Abstract

In a previous publication (Ph. Cuniasse, L.C. Sowers, R. Eritja, B. Kaplan, M.F. Goodman, J.A.H. Cognet, M. Le Bret, W. Guschlbauer and G.V. Fazakerley, Biochemistry 28, 2018 (1989), we determined by two dimensional NMR studies and molecular mechanics calculations the three-dimensional structure of a non-selfcomplementary oligonucleotide:

5′d(C1 P1 G2 P2 G3 P3 dr4 P4 G5 P5 G6 P6 C7)3′

3′d(G13P12C12PllCll P10 C10 P9 C9 P8 G8)5′

where dr, at the center of the first strand, is a model abasic site. In order to explain all the results arising from NMR measurements, we found that an equilibrium between two conformations was necessary. These conformations differ mainly by the sugar pucker of G5 which is C2′ endo or C3′ endo. The latter is stabilized by addition of counterions between phosphate residues P3 and P4.

In this paper, we have constructed systematically, all possible structures as a function of torsion angles delta of dr4 and of G5 by molecular mechanics in the presence or absence of counterions. Since these conformations were not forced with NMR distance measurements, this method allows detailed comparisons between all possible conformations and NMR data. Maps of contour lines of the potential energy, of fits to NMR distance measurements, and of helical twist as a function of torsion angles delta of dr4 and of G5 unravel the difficulties associated with the study of the G5 sugar pucker conformation equilibrium.

Sugar puckers and proton distances are very sensitive criteria to monitor molecular dynamics. Relying on these experimental criteria, we have tested many molecular dynamics preparation phases and we propose a new warm-up and equilibration procedure for molecular dynamics. Thus we show with a 290 ps molecular dynamic run that G5 is in conformational equilibrium and that all NMR data are well reproduced.     

Molecular Dynamics Studies of Electrolyte Solution/Metal Interfaces

Abstract

The results of Molecular Dynamics simulations of pure water near a Pt(100) and a mercury surface as well as an aqueous LiI solution in contact with the Pt(100) surface are reported. The flexible BJH water model is employed in the simulations and the metal-water, ion-water and platinum-ion potentials are derived from molecular orbital calculations. It is shown that the structural and dynamical properties of water and the ions in the adsorbed water layer are significantly different from those in the bulk region.     

Structure and Dynamics of Double Helical DNA in Torsion Angle Hyperspace: A Molecular Mechanics Approach

Abstract

Analysis of the conformational space populated by the torsion angles and the correlation between the conformational energy and the sequence of DNA are important for fully understanding DNA structure and function. Presence of seven variable torsion angles about single covalent bonds in DNA main chain puts a big challenge for such analysis. We have carried out restrained energy minimization studies for four representative dinucleosides, namely d(ApA):d(TpT), d(CpG):d(CpG), d(GpC):d(GpC) and d(CpA):d(TpG) to determine the energy hyperspace of DNA in context to the values of the torsion angles and the structural properties of the DNA conformations populating the favorable regions of this energy hyperspace. The torsion angles were manipulated by constraining their values at the reference points and then performing energy minimization. The energy minima obtained on the potential energy contour plots mostly correspond to the conformations populated in crystal structures of DNA. Some novel favorable conformations that are not present in crystal structure data are also found. The plots also suggest few low energy routes for conformational transitions or the associated energy barrier heights. Analyses of base pairing and stacking possibility reveal structural changes accompanying these transitions as well as the flexibility of different base steps towards variations in different torsion angles.     

Molecular Dynamics Simulation Studies of dTTP Binding and Catalysis Mediated by YhdE Dimerization

YhdE is a Maf-like (multicopy associated filamentation) protein that primarily acts as dTTPase to hydrolyze dTTP into dTMP and two phosphate molecules in cell metabolism pathway. Two crystal structures of YhdE have been previously determined, representing the open and closed active site conformations, respectively. Based on the structures, we have carried out molecular dynamics simulations and free energy calculations to investigate dTTP binding to and hydrolysis by YhdE. Our results suggest that YhdE closed state is structurally more compact than its open state at room temperature. YhdE open state is a favorable conformation for dTTP binding and closed state is a structurally favorable conformation for catalytic reaction. This observation is supported by the structure of YhdE homolog in complex with a nucleotide analog. Free energy calculations reveal that YhdE dimerization occurs preferentially in dTTP binding and is favorable for successive cooperative reaction. The key residues R11, R12 and K80, are found to contribute to the substrate stabilization. Further, YhdE dimerization and binding of dTTP induce the cooperative effect through a direct allosteric communication network in YhdE from the dTTP binding sites in the catalytic center to the intermolecular β-strand in YhdE dimer.     

Energetics of Glutathione Binding to Human Eukaryotic Elongation Factor 1 Gamma: Isothermal Titration Calorimetry and Molecular Dynamics Studies

The energetics of ligand binding to human eukaryotic elongation factor 1 gamma (heEF1γ) was investigated using reduced glutathione (GSH), oxidised glutathione (GSSG), glutathione sulfonate and S-hexylglutathione as ligands. The experiments were conducted using isothermal titration calorimetry, and the findings were supported using computational studies. The data show that the binding of these ligands to heEF1γ is enthalpically favourable and entropically driven (except for the binding of GSSG). The full length heEF1γ binds GSSG with lower affinity (K _d = 115 μM), with more hydrogen-bond contacts (ΔH = ?73.8 kJ/mol) and unfavourable entropy (?TΔS = 51.7 kJ/mol) compared to the glutathione transferase-like N-terminus domain of heEF1γ, which did not show preference to any specific ligand. Computational free binding energy calculations from the 10 ligand poses show that GSSG and GSH consistently bind heEF1γ, and that both ligands bind at the same site with a folded bioactive conformation. This study reveals the possibility that heEF1γ is a glutathione-binding protein.     

Dynamics of Small Molecules in a Dense Polymer Matrix: Molecular Dynamics Studies

Abstract

The anomalous diffusion regime appearing in the self-diffusion of small molecules in bulk amorphous polymers has been extensively studied by molecular dynamics simulations. A rather long simulation of duration ～ 10 ^?8 s is performed on a polyethylene-like simple polymer model containing either oxygen molecules or helium atoms as a diffusant. Dynamic properties evaluated for these diffusants are the mean-square displacement, the van Hove self correlation function, and the self part of the density autocorrelation. It is first confirmed that the anomalous diffusion regime appears in a few hundred picoseconds for oxygen molecule, while the Einstein relation adopted beyond this regime results in the self-diffusion coefficient of the order of ～ 10^?5 cm²/s. This anomaly is still observed for helium that diffuses much faster than oxygen. In the anomalous diffusion regime, it is found that the correlation functions for the two diffusants show characteristic features and become essentially the same as time is scaled appropriately. These features allow the estimation of the two characteristic spatial scales which are probably dominated by the microstructure of the polymer matrix, namely, the cage size and the distance between adjacent cages. The time dependence of the mean-square displacements of the two diffusants can be well interpreted by these characteristic spatial scales as time is scaled with the self-diffusion coefficients. It is shown that the anomalous diffusion regime arises from the inhomogeneous microstructure of the polymer matrix.     

Molecular Dynamics Simulations of a Calmodulin-peptide Complex in Solution

Abstract

We have performed an 4-ns MD simulation of calmodulin complexed with a target peptide in explicit water, under realistic conditions of constant temperature and pressure, in the presence of a physiological concentration of counterions and using Ewald summation to avoid truncation of long-range electrostatic forces. During the simulation the system tended to perform small fluctuations around a structure similar to, but somewhat looser than the starting crystal structure. The calmodulin-peptide complex was quite rigid and did not exhibit any large amplitude domain motions such as previously seen in apo- and calcium-bound calmodulin. We analyzed the calmodulin-peptide interactions by calculating buried surface areas, CHARMM interaction energies and continuum model interaction free energies. In the trajectory, the protein surface area buried by contact with the peptide is 1373 Ų, approximately evenly divided between the calmodulin N-terminal, C-terminal and central linker regions. A majority of this buried surface, 803 ·A², comes from nonpolar residues, in contrast to the protein as a whole, for which the surface is made up of mostly polar and charged groups. Our continuum calculations indicate that the largest favorable contribution to pep- tide binding comes from burial of molecular surface upon complex formation. Electrostatic contributions are favorable but smaller in the trajectory structures, and actually unfavorable for binding in the crystal structure. Since nonpolar groups make up most of buried surface of the protein, our calculations suggest that the hydrophobic effect is the main driving force for binding the helical peptide to calmodulin, consistent with thermodynamic analysis of experimental data. Besides the burial of nonpolar surface area, secondary contributions to peptide binding come from burial of polar surface and electrostatic interactions. In the nonpolar interactions a crucial role is played by the nine methionines of calmodulin. In the electrostatic interactions the negatively charged protein residues and positively charged peptide residues play a dominant role.     

Structural and Dynamical Studies of Humanin in Water and TFE/Water Mixture: A Molecular Dynamics Simulation

Summary

Proline-rich peptides are known to adopt preferentially the extended polyproline II (PPII) helical conformation, which is involved in several protein-protein recognition events. By resorting to molecular modelling techniques, we wished to investigate the extent to which PPII helices could be used for the formation of isohelical peptide-DNA complexes leading to the selective recognition of the major groove of B-DNA. For that purpose, we have grafted to a cationic intercalator, 9-amino-acridine, an oligopeptide having the sequence: Pro-Arg-Pro-Pro-Arg-Pro-Pro-Arg-Pro-Pro-Asp-Pro-Pro. Each residue in the sequence was set in the D configuration, to prevent enzymatic hydrolysis, and each Arg residue was designed to target O6/N7 of a guanine base following the intercalation site. The Asp residue was designed to target a cytosine base, whilst simultaneously forming a bidentate complex with the Arg three residues upstream. Energy-minimization, using the JUMNA procedure, led to the following conclusions: 1) major groove binding is favoured over minor groove or exclusive binding to the phosphates by large energy differences, of over 50 and 90 kcal/mole, respectively; 2) the two best bound sequences are those having three successive guanine bases on the same DNA strand, immediately adjacent to the intercalation site. Sequence d(CGGGC G), encountered in the Primer Binding Site of the HIV retrovirus, thus ranks amongst the best-bound sequences; 3) replacement of an individual guanine amongst the three ones upstream of the intercalation site, by an adenine base, weakens by > 6 kcal/mole the binding energetics; 4) the conformational rigidity of the DNA-bound PPII helix should enable for a modulation of the base sequence selectivity, by appropriate replacements of the Arg and Asp residues. Thus sequence CGGCAAG, also encountered in the HIV genome, could be targeted by an oligopeptide having the sequence Pro-Arg-Pro-Pro-Asp-Pro-Pro- Asn-Pro-Pro-Asn-Pro-Pro-Arg-Ala.     

Molecular Dynamics Simulations of NaCl-type Solid Solution Crystals: The First Application of Molecular Dynamics to Solid Solutions

Abstract

Solid solution crystals appear widely in the fields of earth sciences and inorganic material sciences. The physical properties of solid solutions may vary continuously with chemical composition. Sometimes, linear relationships of the properties with composition are assumed. However, this approximation is not always applicable (e.g., [1], [2], [3]). In order to elucidate the properties of solid solutions, studies on the relation between the macroscopic properties and the atomic configurations (microscopic property) in the crystal are desirable. One of the most effective approaches to the subject is molecular dynamics (MD). However, as far as the authors are aware there have been no molecular dynamics studies on solid solution crystals.     

Molecular Dynamics Simulation Study of DNA Triplex Formed by Mixed Sequences in Solution

Abstract

The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC)?d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usual “cut-off” of the long—range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence—dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other legends.     

Structural Basis for the Binding Affinity of a Homologous Series of Synthetic Phenoxazone Drugs with DNA: NMR and Molecular Mechanics Analysis

Abstract

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5′-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear ¹H-¹H 2D-NOESY and heteronuclear ¹H-³¹P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH₂ groups in the dimethylami- noalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH₃)₂ groups of the side chains.     

The Molecular Mechanics of P- and L-Selectin Lectin Domains Binding to PSGL-1

A laser trap was used to compare the load-dependent binding kinetics between truncated P- and L-selectin to their natural ligand, P-selectin glycoprotein ligand-1 (PSGL-1) over the predicted physiological range of loading rates. Human PSGL-1 was covalently coupled to polystyrene beads. Chimeric selectins were adsorbed to nitrocellulose-coated glass beads on a coverslip. A PSGL-1 bead was held in a laser trap and touched to a vertical surface of the glass bead, allowing a bond to form between selectin and ligand. The surface was moved away from the microsphere, applying load at a constant rate until bond rupture. Rupture force for both selectins increased with loading rate, but significant differences in rupture force between P- and L-selectin were observed only above 460 pN/s. These data are best represented as two energy barriers to unbinding, with the transition from the low to high loading rate regime at 260–290 pN/s. The data also allow the first estimate of a two-dimensional specific on-rate for binding of these two selectins to their natural ligand (1.7 μm²/s). These data suggest that P- and L-selectin lectin domains have very similar kinetics under physiological conditions.     

Solution Conformation of [D-Pen2,D-Pen5]enkephalin in Water: A NMR and Molecular Dynamics Study

The solution conformation of [D -Pen²,D -Pen⁵] enkephalin (DPDPE), a highly potent δ-selective opioid agonist, was examined by means of NMR, molecular mechanics and molecular dynamics methods. The structural information in the solvent water was obtained employing one- and two-dimensional methods of ¹H and ¹³C-NMR spectroscopy. Based on the distance geometry technique using the ROE data as input, 400 conformers were obtained and considered in the structure analysis. Alternatively, about 2000 conformers were stochastically generated and related to the NMR data after energy minimization. The structure analysis provides one conformer in agreement with all NMR data, which belongs to the lowest energy conformation group. This structure may serve as a reference conformer for DPDPE analogues synthesized with the aim of activity increase.     

Molecular Determinants of Epidermal Growth Factor Binding: A Molecular Dynamics Study

The epidermal growth factor receptor (EGFR) is a member of the receptor tyrosine kinase family that plays a role in multiple cellular processes. Activation of EGFR requires binding of a ligand on the extracellular domain to promote conformational changes leading to dimerization and transphosphorylation of intracellular kinase domains. Seven ligands are known to bind EGFR with affinities ranging from sub-nanomolar to near micromolar dissociation constants. In the case of EGFR, distinct conformational states assumed upon binding a ligand is thought to be a determining factor in activation of a downstream signaling network. Previous biochemical studies suggest the existence of both low affinity and high affinity EGFR ligands. While these studies have identified functional effects of ligand binding, high-resolution structural data are lacking. To gain a better understanding of the molecular basis of EGFR binding affinities, we docked each EGFR ligand to the putative active state extracellular domain dimer and 25.0 ns molecular dynamics simulations were performed. MM-PBSA/GBSA are efficient computational approaches to approximate free energies of protein-protein interactions and decompose the free energy at the amino acid level. We applied these methods to the last 6.0 ns of each ligand-receptor simulation. MM-PBSA calculations were able to successfully rank all seven of the EGFR ligands based on the two affinity classes: EGF>HB-EGF>TGF-α>BTC>EPR>EPG>AR. Results from energy decomposition identified several interactions that are common among binding ligands. These findings reveal that while several residues are conserved among the EGFR ligand family, no single set of residues determines the affinity class. Instead we found heterogeneous sets of interactions that were driven primarily by electrostatic and Van der Waals forces. These results not only illustrate the complexity of EGFR dynamics but also pave the way for structure-based design of therapeutics targeting EGF ligands or the receptor itself.