Molecular dynamics simulation study of the interaction of trehalose with lipid membranes

The interactions of the cryoprotective agent trehalose with a lipid membrane made of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine at 323 K were studied by means of molecular dynamics simulations. It was observed that trehalose binds to the phospholipid headgroups with its main axis parallel to the membrane normal. Trehalose establishes hydrogen bonds with the carbonyl and phosphate groups and replaces water molecules from the lipid headgroup. Notably, the number of hydrogen bonds (HBs) that the membrane made with its environment was conserved after trehalose binding. The HBs between lipid and trehalose have a longer lifetime than those established between lipid and water. The binding of the sugar does not produce changes either in the lipid area or in the lipid order parameter. The effect of trehalose on the dipole potential is in agreement with experimental results. The contribution of the different components to the membrane dipole potential was analyzed. It was observed that the binding of trehalose produces changes in the different components and the sugar itself contributes to the surface potential due to the polarization of its hydroxyl in the interface.     

Molecular interaction energies and optimal configuration of a cubane dimer

We have studied the dependence of the binding energy of a cubane dimer on the mutual orientation of and the distance between the composing monomers employing the second‐order Møller–Plesset perturbation scheme (MP2) with the cc‐pVDZ molecular basis set. We have found that the MP2 contribution from the molecular correlations is responsible for the bound state of the cubane dimer, whereas the Hartree–Fock contribution remains anti‐bonding at all intermolecular distances. Starting with two molecules in the standard orientation and centers of mass at (0,0,0) and (0,0,d), respectively, the maximal binding energy is found at d = 5.125 Å and one of the monomers rotated by 45° about the z‐axis. This configuration implies that the hydrogen atoms belonging to different monomers tend to repel each other. The results are in agreement with experimental data on the optimal packing of cubane molecules in the solid state. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Molecular dynamics study of perovskite structures with modified interatomic interaction potentials

The structure of compounds with the perovskite structure ABX₃ (A and B are cations, X are anions O^2—, F^—, Cl^—, Br^—, and I^—), which are widely used in engineering due to unique electrical, optical, and photovoltaic properties, has been considered. Hybrid organic—inorganic halide perovskites important for photovoltaics of a new generation are worth mentioning; they contain cations of organic nitrogen bases as monovalent cations. A molecular dynamics (MD) study of the CaTiO₃ base structure (Ca²⁺, Ti⁴⁺, and O^2—) has been performed in order to develop the methodology of computer simulation and optimization of the shape and parameters of atomic potentials for perovskite systems.     

Theoretical study on Pd dimer and trimer interaction with the hydrogen molecule

The H₂ interaction with the Pd dimer and trimer were studied using multiconfigurational self-consistent field (MC-SCF) calculations with the relativistic effective core potential (RECP); the correlation energy correction was included in the extended multireference configuration interaction (MRCI), variational and perturbative to second order. Here, we considered the Pd₂ first six states: ³Σ⁺_u, ¹Σ⁺_g, ³Π_g, ³Δ_xy, ¹Σ⁺_u, and ³Σ⁺_g. For them, the four geometrical approaches included were the side-on H₂ toward Pd₂, for the hydrogen molecule in and out the Pd dimer plane; the perpendicular end-on H₂ toward Pd₂; and the perpendicular end-on Pd₂ to H₂. The Pd₂ ground state is ³Σ⁺_u, which only captures H₂ in the C_2v end-on approach, softly relaxing the H(SINGLE BOND)H bond. The closed-shell ¹Σ⁺_g captures the H₂ molecule in all the approaches considered: The side-on approach of this state presents deep wells and relaxes the H(SINGLE BOND)H bond, and the end-on approach captures H₂ with a relatively longer H(SINGLE BOND)H distance and also a deep well. The ³Π_g state was the only one which did not capture H₂. For the triangular Pd₃ clusters, H₂ was approached in the C_2v symmetry in and out of the Pd₃ plane. In the triangular case, H₂ was absorbed in both spin states, with deep wells and relaxing the H(SINGLE BOND)H distance. The linear Pd₃ singlet and triplet states capture outside of the Pd₃ and break the H(SINGLE BOND)H bond. © 1997 John Wiley & Sons, Inc.     

Theoretical study on intermolecular interaction of epoxyethane dimer

Ab initio calculations at Hartree–Fock and fourth‐order Mø ller–Plesset (MP4) correlation correction levels with 6‐31G* basis set have been performed on the epoxyethane dimer. Dimer binding energies have been corrected for the basis set superposition error (BSSE) and the zero‐point energy. The greatest corrected dimer binding energy is −8.36 kJ/mol at the MP4/6‐31G*//HF/6‐31G* level. The natural bond orbital analysis has been performed to trace the origin of the weak interactions that stabilize dimer. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 94–98, 2000     

Molecular dynamics studies on the interaction of 4-acetylamino-5-hydroxynapthalene-2,7-disulfonic acid with catalytic domain of avian sarcoma virus integrase dimer

We report here 500 ps molecular dynamics (MD) simulation results on two molecules of 4-acetylamino-5-hydroxy naphthalene 2-7-disulfonic acid (Y-3) in the catalytic domain (residues 54–199) of avian sarcoma virus integrase (ASV-IN) dimer. Starting model was obtained on the basis of PDB coordinates (PDB code 1A5X) of ASV-IN-Y3 monomer by Lubkowski et al. [Proc. Natl Acad. Sci. USA, 95 (1998) 4831–4836]. Two molecules of Y3 were docked in the active cavity of the dimer using IMF [Phys. Edu., 5 (1988) 169–176]. Energy minimization (EM) and MD simulation were carried out using Sander's module of [ : Assisted Model building with Energy Refinement: a computer simulation software developed by the University of California, USA, 1997] with all atom force field for Y3 and united atom force field for IN. Analysis of the ligand protein interaction and perturbative changes in the complex is presented in the paper. Salient features related to higher pharmacological activity of Y-3 compared to other analogs from acetyl-amino-naphthalene series are enumerated.     

A density-functional study on pi-aromatic interaction: benzene dimer and naphthalene dimer

The long-range correction (LC) scheme of density-functional theory (DFT) was applied to the calculation of the pi-aromatic interaction of the benzene dimer and naphthalene dimer. In previous calculations, it was confirmed that the LC scheme [Iikura et al., J. Chem. Phys. 115, 3540 (2001)] gives very accurate potential- energy surfaces (PESs) of small van der Waals (vdW) complexes by combining with the Anderson-Langreth-Lundqvist (ALL) vdW correlation functional [Andersson et al., Phys. Rev. Lett. 76, 102 (1996)] (LC-DFT + ALL). In this study, LC-DFT+ALL method was examined by calculating a wide range of PES of the benzene dimer including parallel, T-shaped, and parallel-displaced configurations. As a result, we succeeded in reproducing very accurate PES within the energy deviance of less than 1 kcalmol in comparison with the results of high-level ab initio molecular-orbital methods at all reference points on the PES. It was also found that LC-DFT + ALL gave accurate results independent of exchange-correlation functional used, in contrast with the strong functional dependencies of conventional pure functionals. This indicates that both exchange repulsion and van der Waals attractive interactions should be correctly incorporated in conventional pure functionals in order to calculate accurate pi-aromatic interactions. We also found that LC-DFT + ALL method has a low basis-set dependency in the calculations of pi-aromatic interactions. The present scheme was also successfully applied to the pi,[ellipsis (horizontal)],pi stacking interactions of naphthalene dimer. This may suggest that LC-DFT + ALL method would be a powerful tool in the calculations of large molecules such as biomolecules.     

Diastereochemically controlled porphyrin dimer formation on a DNA duplex scaffold

DNA-porphyrin conjugates were designed and synthesized for the preparation of the conformationally controlled porphyrin dimer structures constructed on a d(GCGTATACGC)2. Porphyrin derivatives were introduced to the central TATpA sequence where p represents the phosphoramidate for the attachment of the free-base porphyrin (FbP) and zinc-coordinated porphyrin (ZnP), which allows contact of the two porphyrins in the minor groove. The porphyrin dimers were characterized using CD, UV-vis, steady-state, and time-resolved fluorescence spectroscopies, indicating that the porphyrins form face-to-face conformations. Also the co-facial conformation was confirmed by comparison with spectra of the non-self-complementary duplex containing one porphyrin moiety. Introduction of zinc into porphyrin moiety destabilized the duplex formation. Two diastereomers showed different thermal stabilities and affected the conformations of porphyrin dimers. The temperature-dependent assembly and the conformational change of the porphyrin dimer on the duplex DNA were observed in the UV-vis spectra, indicating that the dynamic movement of the porphyrin dimer occurs on the duplex. The results indicate that the porphyrin dimers of DNA-FbP conjugates are overlapped clockwise and are located in the minor groove of the usual B-form DNA backbone. The interaction and conformation of two porphyrin moieties are controlled by the following three factors: (1) temperature change during and after formation of the duplex porphyrins at lower temperature; (2) diastereochemistry of the phosphoramidates where porphyrins are connected via a linker; and (3) zinc ion coordination that destabilizes the interaction of porphyrins as well duplex formation.     

Force-induced insulin dimer dissociation: a molecular dynamics study

Understanding the forces and dynamics of insulin dissociation is critical for devising formulations for the treatment of insulin-dependent diabetes. In earlier work, we applied AFM-based force spectroscopy to covalently tethered and oriented insulin monomers to assess the effect of molecular orientation on insulin-insulin binding forces. We report here on steered molecular dynamics simulations of the insulin dissociation force spectroscopy experiment. Consistent with our experiments, our simulation results suggest that insulin dimer dissociation occurs near the limit of extensibility of the B-chain. We have also found that the forced dissociation of the insulin dimer is a rate-dependent process, involving significant conformational changes to the monomer(s). The insulin dimer dissociation pathway also depends on the relative strength of the inter-monomer interactions across the antiparallel beta-sheet interface and the intra-monomer interactions of residues A1 and A30 with the insulin B-chain. Our simulation results strongly support the design of bioactive insulin analogues that involves altering hydrogen bonding and hydrophobic interactions across the beta-sheet dimer interface.     

A theoretical study on the intermolecular interaction of energetic system-Nitromethane dimer

Three optimized geometries of nitromethane dimer have been obtained at the HF/6-31G level.Dimer binding energies have been corrected for the basis set superposition error (BSSE) and the zero point energy.Computed results indicate that the cyclic structure of (CH3NO2)2 is the most stable of three optimized geometries,whose corrected binding energyis 17.29 kJ mol-1 at the MP4SDTQ/6-31G//HF/6-31G level.In the optimized structures of nitromethane dimer,the inter-molecular hydrogen bond has not been found; and the charge-transfer interaction between CH3NO2 subsystems is weak; and the correlation interaction energy makes a little contribution to the intermolecular interaction energy of the dimer.     

Molecular dynamics simulation of sub-and supercritical water with a new interaction potential

A molecular dynamics experiment was performed for water under sub-and supercritical conditions with a new interaction potential including the nonelectrostatic H-bond component. The internal energy, self-diffusion coefficient, mean number of H-bonds per water molecule, and distributions of molecules according to the number of H-bonds were calculated over a wide temperature range at pressures of 50 and 100 MPa. The temperature dependences of these properties were analyzed.     

Molecular dynamics simulation study of superhydrated perdeuterated natrolite using a new interaction potential model

To test a new interaction potential, molecular dynamics simulations of zeolite natrolite were performed for the structures under ambient conditions hydrated by perdeuterated water and at high pressure (1.87 GPa) in the superhydrated phase, which were recently studied by neutron diffraction. The experimental structures were reproduced with reasonable accuracy, and the hydrogen bond features are discussed. As in ordinary natrolite, a flip motion of water molecules around the HOH bisector is found, which, together with translational oscillations, gives rise to transient hydrogen bonds between water molecules, which do not appear from experimental equilibrium coordinates. The dynamics of water molecules can explain some problems encountered in refining the experimental structure. Vibrational spectra of natrolite containing perdeuterated water, which are not yet measured, were simulated, and their qualitative trend is discussed.     

Molecular dynamics simulation study of interaction between model rough hydrophobic surfaces

We study some aspects of hydrophobic interaction between molecular rough and flexible model surfaces. The model we use in this work is based on a model we used previously (Eun, C.; Berkowitz, M. L. J. Phys. Chem. B 2009, 113, 13222-13228), when we studied the interaction between model patches of lipid membranes. Our original model consisted of two graphene plates with attached polar headgroups; the plates were immersed in a water bath. The interaction between such plates can be considered as an example of a hydrophilic interaction. In the present work, we modify our previous model by removing the charge from the zwitterionic headgroups. As a result of this procedure, the plate character changes: it becomes hydrophobic. By separating the total interaction (or potential of mean force, PMF) between plates into the direct and the water-mediated interactions, we observe that the latter changes from repulsive to attractive, clearly emphasizing the important role of water as a medium. We also investigate the effect of roughness and flexibility of the headgroups on the interaction between plates and observe that roughness enhances the character of the hydrophobic interaction. The presence of a dewetting transition in a confined space between charge-removed plates confirms that the interaction between plates is strongly hydrophobic. In addition, we notice that there is a shallow local minimum in the PMF in the case of the charge-removed plates. We find that this minimum is associated with the configurational changes that flexible headgroups undergo as the two plates are brought together.     

Molecular dynamics simulation study on the interaction of KRN 7000 and three analogues with human CD1d

Many analogues of KRN 7000, a synthetic glycolipid (α-galactosylceramide) exhibiting immuno-stimulatory activity and antitumor properties, were previously synthesized and tested in order to understand the reasons for the resulting biological activity and Th1/Th2 cytokine profile. Principles have been established for the interaction of such glycolipids with the human CD1d molecule but the exact mechanism by which different ligands with the same polar head elicit distinct biological responses remains unclear. Based on these experiments and on the available crystal structures, protein-ligand interactions are explored using molecular dynamics simulations. Hydrogen bond interactions are examined with regard to the polar group orientation. The influence of modulations on the dynamic behavior of the CD1d-glycolipid complex is addressed. Overall, our data support the mechanism by which the shortening of the α-GalCer sphingosine chain causes a significant twist of the CD1d α1 helix structure from residue Phe84 that affects the position of CD1d residues involved in the TCR recognition.     

Molecular dynamics simulation study of the water-mediated interaction between zwitterionic and charged surfaces

We calculated the potential of mean force (PMF) for the interaction between a model zwitterionic bilayer and a model charged bilayer. To understand the role of water, we separated the PMF into two components: one due to direct interaction and the other due to water-mediated interaction. In our calculations, we observed that water-mediated interaction is attractive at larger distances and repulsive at shorter. The calculation of the entropic and enthalpic contributions to the solvent-mediated components of the PMF showed that attraction is entropically dominant, while repulsion is dominated by the enthalpy.     

Molecular dynamics simulations of formamide interaction with hydrocyanic acid on a catalytic surface TiO2

The behavior of water—formamide and hydrocyanic acid—formamide solutions on an anatase surface have been studied using molecular dynamics (MD) simulation method. The interaction activation energies have been estimated for the temperature range from 250 up to 400 K. The diffusion coefficients and structural radial distribution functions have been calculated for the formamide, water and hydrocyanic acid on an anatase surface. The calculated activation energies of the water—formamide—anatase and hydrocyanic acid—formamide—anatase systems were analyzed and compared. A comparative analysis of the systems under investigation was performed and a possible correlation between the obtained MD results and the molecular mechanism involving the formamide’s interaction with dioxide titan adsorbing surface were discussed.     

Molecular interaction of N-formylmorpholine with alkanols,a spectroscopic study

The interactions of N-formylmorpholine (NFM) with different alkanols were monitored using proton nuclear magnetic resonance and IR spectroscopy. It was found that the NFM interaction with alkanols decreases with increasing length and branching of the alkanols. A good correlation was obtained between the proton chemical shift and the volume of mixing of alkanols with NFM.     

Molecular spectroscopic study on the interaction of tetracyclines with serum albumins

A molecular spectroscopic investigation of the interaction between tetracyclines antibiotics and human serum albumin or bovine serum albumin was reported. The influences of some metal ions on the interaction were also studied. When tetracyclines drugs were added into the solution containing serum albumins, the fluorescence intensity of serum albumins decreased with the increasing of the drugs concentrations, which is due to the formation of new non-fluorescence complexes of drug-serum albumin. The tetracyclines acted as quenchers and quenched the fluorescence of the serum albumins. The binding constants and the number of the binding sites of the reaction of tetracyclines and serum albumins were obtained. The main sorts of acting force between the drugs and serum albumins were found and the action distances and the energy transfer efficiencies between donor-acceptor were calculated based on the Foster energy transference.     

Base flipping of the thymine dimer in duplex DNA

Exposure of two adjacent thymines in DNA to UV light of 260-320 nm can result in the formation of the cis,syn-cyclobutane pyrimidine dimer (CPD). The structure of DNA containing an intrahelical CPD lesion has been previously studied experimentally and computationally. However, the structure of the extrahelical, flipped-out, CPD lesion, which has been shown to be the structure that binds to the CPD repair enzyme, DNA photolyase, has yet to be reported. In this work the structure of both the flipped-in and the flipped-out CPD lesions in duplex DNA is reported. These structures were calculated using 8 ns molecular dynamics (MD) simulations. These structures are then used to define the starting and ending points for the base-flipping process for the CPD lesion. Using a complex, two-dimensional pseudodihedral coordinate, the potential of mean force (PMF) for the base-flipping process was calculcated using novel methodology. The free energy of the flipped-out CPD is roughly 6.5 kcal/mol higher than that of the flipped-in state, indicating that the barrier to flipping out is much lower for CPD than for undamaged DNA. This may indicate that the flipped-out CPD lesion may be recognized by its repair enzyme, DNA photolyase, whereas previous studies of other damaged, as well as nondamaged, bases indicate that they are recognized by enzymes in the intrahelical, flipped-in state.