The solvation of bromide anion in acetonitrile: a structural study based on the combination of theoretical calculations and X-ray absorption spectroscopy

This work studies the solvation of bromide in acetonitrile by combining quantum mechanics, computer simulations and X-ray absorption near edge structure (XANES) spectroscopy. Three different sets of interaction potentials were tested, one of them derived from literature and the other two are simple modifications of the previous one to include specificities of the bromide–acetonitrile interactions. Results for microsolvation of bromide were obtained by quantum mechanical optimization and classical minimization of small clusters [Br(ACN) _n ]^?  (n = 9, 20). Analysis of molecular dynamics (MD) simulations has provided structural, dynamic and energetic aspects of the solvation phenomenon. The theoretical computation of Br K-edge XANES spectrum in solution using the structural information obtained from the different simulations has allowed the comparison among the three different potentials, as well as the examination of the main structural and dynamic factors determining the shape of the experimental spectrum.     

Simulated dynamics of Ne@C60 aggregates beyond dissociation

Molecular dynamics (MD) computer simulations are utilized to better understand the dynamics of small (N = 5) endohedral Ne@C₆₀ aggregates. Multiple runs at various temperatures are used to increase the reliability of our statistics. The aggregate holds together until somewhere between T = 1150 and 1200 K, where it dissociates, showing no intermediate sign of melting or fullerene disintegration. When the temperature is increased to around T = 4000 K, the encapsulated neon atoms begin to leave the aggregate, with the fullerene molecules still remaining intact. At temperatures near T = 4400 K, thermal disintegration of the fullerenes preempts the aggregate dissociation. Above this temperature neon atoms are more quickly released and the fullerenes form a larger connected structure, with bonding taking place in atom pairs from different original fullerene molecules. Escape constants and half lives are calculated for the temperature range 4000 K ≤ T ≤ 5000 K. The agreements and disagreements of results of this work with experiments suggest that classical MD simulations are useful in describing fullerene systems at low temperatures and near disintegration, but require development of new techniques before it is possible to accurately model windowing at temperatures below T = 3000 K.     

Multishelled Gold Nanowires

Abstract

The current miniaturization of electronic devices raises many questions about the properties of various materials at nanometre-scales. Recent molecular dynamics computer simulations have shown that small finite nanowires of gold exist as multishelled structures of lasting stability. These classical simulations are based on a well-tested embedded atom potential. Molecular dynamics simulation studies of metallic nanowires should help in developing methods for their fabrication, such as electron-beam litography and scanning tunneling microscopy.     

Helix stability of oligoglycine,oligoalanine, and oligo‐β‐alanine dodecamers reflected by hydrogen‐bond persistence

Helices are important structural/recognition elements in proteins and peptides. Stability and conformational differences between helices composed of α‐ and β‐amino acids as scaffolds for mimicry of helix recognition has become a theme in medicinal chemistry. Furthermore, helices formed by β‐amino acids are experimentally more stable than those formed by α‐amino acids. This is paradoxical because the larger sizes of the hydrogen‐bonding rings required by the extra methylene groups should lead to entropic destabilization. In this study, molecular dynamics simulations using the second‐generation force field, AMOEBA (Ponder, J.W., et al., Current status of the AMOEBA polarizable force field. J Phys Chem B, 2010. 114 (8): p. 2549–64.) explored the stability and hydrogen‐bonding patterns of capped oligo‐β‐alanine, oligoalanine, and oligoglycine dodecamers in water. The MD simulations showed that oligo‐β‐alanine has strong acceptor+2 hydrogen bonds, but surprisingly did not contain a large content of 3₁₂‐helical structures, possibly due to the sparse distribution of the 3₁₂‐helical structure and other structures with acceptor+2 hydrogen bonds. On the other hand, despite its backbone flexibility, the β‐alanine dodecamer had more stable and persistent <3.0 Å hydrogen bonds. Its structure was dominated more by multicentered hydrogen bonds than either oligoglycine or oligoalanine helices. The 3₁ (PII) helical structure, prevalent in oligoglycine and oligoalanine, does not appear to be stable in oligo‐β‐alanine indicating its competition with other structures (stacking structure as indicated by MD analyses). These differences are among the factors that shape helical structural preferences and the relative stabilities of these three oligopeptides. Proteins 2014; 82:3043–3061. © 2014 Wiley Periodicals, Inc.     

Molecular dynamics simulation of chitinase I from Thermomyces lanuginosus SSBP to ensure optimal activity

Abstract

The fungal chitinase I obtained from Thermomyces lanuginosus SSBP, a thermophilic deuteromycete, has an optimum growth temperature and pH of 323.15 K and 6.5, respectively. This enzyme plays an important task in the defence mechanism of organisms against chitin-containing parasites by hydrolysing β-1, 4-linkages in chitin. It acts as both anti-fungal and biofouling agents, with some being thermostable and suitable for the industrial applications. Three-dimensional model of chitinase I enzyme was predicted and analysed using various bioinformatics tools. The structure of chitinase I exhibited a well-defined TIM barrel topology with an eight-stranded α/β domain. Structural analysis and folding studies at temperatures ranging from 300 to 375 K using 10 ns molecular dynamics simulations clearly showed the stability of the protein was evenly distributed even at higher temperatures, in accordance with the experimental results. We also carried out a number of 20 ns constant pH molecular dynamics simulations of chitinase I at a pH range 2–6 in a solvent. This work was aimed at establishing the optimum activity and stability profiles of chitinase I. We observed a strong conformational pH dependence of chitinase I and the enzyme retained their characteristic TIM barrel topology at low pH.     

Ab initio and molecular dynamics studies of solid β-HMX: effects of hydrostatic pressure and high temperature

Using first-principles density functional theory and classical molecular dynamics (MD), the structural, electronic and mechanical properties of the energetic material β-HMX have been studied. The crystal structure optimised by the local density approximation calculations compares reasonably with the experimental data. Electronic band structure and density of states indicate that β-HMX is an insulator with a band gap of 3.059 eV. The pressure effect on the crystal structure and physical properties has been investigated in the range of 0–40 GPa. The crystal structure and electronic properties change slightly as the pressure increases from 0 to 2.5 GPa; when the pressure is above 2.5 GPa, further increment of the pressure results in significant changes in crystal structure. There is a larger compression along the b-axis than along the a- and c-axes. Isothermal–isobaric MD simulations on β-HMX were performed in the temperature range of 5–400 K. Phase transition at 360 K, corresponding to a volume interrupt, was found. The computed thermal expansion coefficients show anisotropic behaviour with a slightly larger expansion along the b- and c-axes than along the a-axis. In the temperature range of 5–360 K, β-HMX possesses good plasticity and its stiffness decreases with increasing the temperature.     

Exploring the structure and dynamics of nano-confined water molecules using molecular dynamics simulations

Confinement effects can lead to drastic changes in the structural and dynamical properties of water molecules. In this work, we have performed classical molecular dynamics simulations of endohedral fullerenes of type (H₂O)_n@C_m (n = 1, 12, 21, 62, 108 and m = 60, 180, 240, 500 and 720) to explore the effects of spherical confinement on water properties. It is shown that these confined water molecules can form distinct solvation pattern depending upon the available space inside the fullerene cavity. For the systems with smaller diameter, cage-like structure is predominant whereas bulk-like structure is observed for larger fullerenes. The orientational relaxation of these confined water molecules showed slower relaxation as the cavity diameter increases except for the (H₂O)₂₁@C₂₄₀. In this case, stable cage-like structure hinders the overall dynamics of the trapped water molecules. Finally, we have calculated the hydrogen bond lifetimes from the hydrogen bond time correlation functions and compared with that of bulk water.     

Structure prediction and functional analyses of a thermostable lipase obtained from Shewanella putrefaciens

Previous experimental studies on thermostable lipase from Shewanella putrefaciens suggested the maximum activity at higher temperatures, but with little information on its conformational profile. In this study, the three-dimensional structure of lipase was predicted and a 60 ns molecular dynamics (MD) simulation was carried out at temperatures ranging from 300 to 400 K to better understand its thermostable nature at the molecular level. MD simulations were performed in order to predict the optimal activity of thermostable lipase. The results suggested strong conformational temperature dependence. The thermostable lipase maintained its bio-active conformation at 350 K during the 60 ns MD simulations.     

Molecular dynamics simulation and binding free energy studies of novel leads belonging to the benzofuran class inhibitors of Mycobacterium tuberculosis Polyketide Synthase 13

In this work, the binding mechanism of new Polyketide Synthase 13 (Pks13) inhibitors has been studied through molecular dynamics simulation and free energy calculations. The drug Tam1 and its analogs, belonging to the benzofuran class, were submitted to 100 ns simulations, and according to the results obtained for root mean square deviation, all the simulations converged from approximately 30 ns. For the analysis of backbone flotation, the root mean square fluctuations were plotted for the Cα atoms; analysis revealed that the greatest fluctuation occurred in the residues that are part of the protein lid domain. The binding free energy value (ΔG_bind) obtained for the Tam16 lead molecule was of ?51.43 kcal/mol. When comparing this result with the ΔG_bind values for the remaining analogs, the drug Tam16 was found to be the highest ranked: this result is in agreement with the experimental results obtained by Aggarwal and collaborators, where it was verified that the IC₅₀ for Tam16 is the smallest necessary to inhibit the Pks13 (IC₅₀ = 0.19 μM). The energy decomposition analysis suggested that the residues which most interact with inhibitors are: Ser1636, Tyr1637, Asn1640, Ala1667, Phe1670, and Tyr1674, from which the greatest energy contribution to Phe1670 was particularly notable. For the lead molecule Tam16, a hydrogen bond with the hydroxyl of the phenol not observed in the other analogs induced a more stable molecular structure. Aggarwal and colleagues reported this hydrogen bonding as being responsible for the stability of the molecule, optimizing its physic-chemical, toxicological, and pharmacokinetic properties.     

A Combined 2D-NMR and Molecular Dynamics Analysis of the Structure of the Actinomycin D: d(ATGCAT)2 Complex

Abstract

We present a comparative analysis of an NMR experiment and molecular and harmonic dynamics simulations of an actinomycin D: d(ATGCAT)₂ complex. A comparison of NOE measurements and 1/R⁶ weighted proton-proton distances confirm the general correctness of the Actinomycin D-DNA model proposed by Sobell. There are, however, some substantial differences between the proton-proton distances inferred from the NOE results and the molecular and harmonic dynamics simulations. The remaining discrepancies could either come from contributions of other conformations to the average properties of the complex or from uncertainties in the NMR distance analysis. An analysis of the molecular dynamics helix properties, sugar puckers, hydrogen bonding, rms fluctuations and torsional properties are qualitatively consistent with those from previous simulations, but the presence of an intercalated drug leads to some new structural and dynamical features.     

Molecular Dynamics Simulations of Rhodopsin Point Mutants at the Cytoplasmic Side of Helices 3 and 6

Abstract

The present work reports on a structural analysis carried out through different computer simulations of a set of rhodopsin mutants with differential functional features in regard to the wild type. Most of these mutants, whose experimental features had previously been reported [Ramon et al. J Biol Chem 282, 14272–14282 (2007)], were designed to perturb a network of electrostatic interactions located at the cytoplasmic sides of transmembrane helices 3 and 6. Geometric and energetic features derived from the detailed analysis of a series of molecular dynamics simulations of the different rhodopsin mutants, involving positions 134(3.49), 247(6.30), and 251(6.34), suggest that the protein structure is sensitive to these mutations through the local changes induced that extend further to the secondary structure of neighboring helices and, ultimately, to the packing of the helical bundle. Overall, the results obtained highlight the complexity of the analyzed network of electrostatic interactions where the effect of each mutation on protein structure can produce rather specific features.     

Effects of hesperidin,a flavanone glycoside interaction on the conformation,stability, and aggregation of lysozyme: multispectroscopic and molecular dynamic simulation studies?

Hesperidin (HESP), a flavanone glycoside, shows high antioxidant properties and posses ability to go through the blood–brain barrier. Therefore, it could be a potential drug molecule against aggregation based diseases such as Alzheimer’s, Parkinson’s, and systemic amyloidoses. In this work, we investigated the potential of HESP to interact with hen egg-white lysozyme (HEWL) monomer and prevent its aggregation. The HESP–HEWL binding studies were performed using a fluorescence quenching technique, molecular docking and molecular dynamics simulations. We found a strong interaction of HESP with the lysozyme monomer (K_a, ~ 5 × 10⁴ M^?1) mainly through hydrogen bonding, water bridges, and hydrophobic interactions. We showed that HESP molecule spanned the highly aggregation prone region (amino acid residues 48-101) of HEWL and prevented its fibrillar aggregation. Further, we found that HESP binding completely inhibited amorphous aggregation of the protein induced by disulfide-reducing agent tries-(2-carboxyethyl) phosphine. Conformational and stability studies as followed by various tertiary and secondary structure probes revealed that HESP binding only marginally affected the lysozyme monomer conformation and increased both stability and reversibility of the protein against thermal denaturation. Future studies should investigate detail effects of HESP on solvent dynamics, structure, and toxicity of various aggregates. The answers to these questions will not only target the basic sciences, but also have application in biomedical and biotechnological sciences.     

Solvation free energies of nucleic acid bases in ionic liquids

The solvation free energies of five nucleic acid bases in [C_nbim]Br (where n = 2, 4, 6) ionic liquids (ILs) were computed using the Bennett acceptance ratio (BAR) method employing molecular dynamics simulations. The computed free energies using BAR were in agreement with other methods. The large and negative predicted free energies of the bases in ILs indicated that the bases were better solvated in the ILs rather than in water. Hydrogen bonding interactions between polar sites of the bases and ILs’ ions significantly contributed to the solvation mechanism.     

Conformation dependence of the CαDα stretch mode in peptides. II. Explicitly hydrated alanine peptide structures

Our previous studies of the potential utility of the C^αD^α stretch frequency, ν(CD), as a tool for determining conformation in peptide systems (Mirkin and Krimm, J Phys Chem A 2004, 108, 10923–10924; 2007, 111, 5300–5303) dealt with the spectroscopic characteristics of isolated alanine peptides with α_R, β, and polyproline II structures. We have now extended these ab initio calculations to include various explicit‐water environments interacting with such conformers. We find that the structure‐discriminating feature of this technique is in fact enhanced as a result of the conformation‐specific interactions of the bonding waters, in part due to our finding (Mirkin and Krimm, J Phys Chem B 2008, 112, 15268) that C^α? D^α…O(water) hydrogen bonds can be present in addition to those expected between water and the CO and NH of the peptide groups. In fact, ν(CD) is hardly affected by the latter bonding but can be shifted by up to 70 cm^?1 by the former hydrogen bonds. We also discuss the factors that will have to be considered in developing the molecular dynamics (MD) treatment needed to satisfactorily take account of the influence of outer water layers on the structure of the first‐layer water molecules that hydrogen bond to the peptide backbone. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 791–800, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Non-equilibrium molecular dynamics study of nanoscale thermal contact resistance

Interfaces play an important role in microscale and nanoscale heat transfer processes with molecular dynamics (MD) simulations often used to study these interfacial phenomena. In this study, two models were used to simulate thermal conduction across micro contact points and the thermal contact resistance using non-equilibrium molecular dynamics simulations with consideration of the near field radiation. When the ratio of the length of the micro contact to the length of the conduction region is less than 0.125, the influence of the near field radiation should be considered; but when the ratio is larger than 0.2, it can be neglected. When the computational domain sizes are 8.50 × 10.62 × 8.50 nm and 10.62 × 10.62 × 10.62 nm, the MD results show that the thermal contact resistance exponentially increases with decreasing area of the micro contact point and increases with increasing micro contact layer thickness. The MD thermal contact resistances in nanoscale are much larger than that of the classical thermal analysis since the material thermal conductivity reduction is ignored in the classical model. The results also show that material defects increase the thermal resistance.     

On the Occurrence of Three-Center Hydrogen Bonds in Cyclodextrins in Crystalline Form and in Aqueous Solution: Comparison of Neutron Diffraction and Molecular Dynamics Results

Abstract

Three-center (bifurcated) hydrogen bonds may play a role by serving as an intermediate state between different dynamically changing hydrogen bonding patterns. Hydrogen bonding configurations can be studied experimentally by neutron diffraction and theoretically by computer simulation techniques. Here, both methods are used to analyse the occurrence of three-center hydrogen bonds in crystals of cyclodextrins.

Almost all experimentally observed three-center hydrogen bonds in the crystal are reproduced in the molecular dynamics (MD) simulations, even as far as the detailed asymmetric geometry is concerned. On the basis of this result a MD simulation of cyclodextrin in aqueous solution is searched for the occurrence of three-center hydrogen bonds. Significant differences are found. In solution more different three-center hydrogen bonds per α-cyclodextrin molecule are observed than in the crystal but the population (existence as percent of the simulation period) of each three-center hydrogen bond is lower in solution than in crystal. These may indeed serve as intermediate states in the process of changing one hydrogen bonding pattern into another.     

3D-QSAR (CoMFA,CoMSIA), molecular docking and molecular dynamics simulations study of 6-aryl-5-cyano-pyrimidine derivatives to explore the structure requirements of LSD1 inhibitors

Recently, Histone Lysine Specific Demethylase 1 (LSD1) was regarded as a promising anticancer target for the novel drug discovery. And several small molecules as LSD1 inhibitors in different structures have been reported. In this work, we carried out a molecular modeling study on the 6-aryl-5-cyano-pyrimidine fragment LSD1 inhibitors using three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to generate 3D-QSAR models. The results show that the best CoMFA model has q² = 0.802, r²_ncv = 0.979, and the best CoMSIA model has q² = 0.799, r²_ncv = 0.982. The electrostatic, hydrophobic and H-bond donor fields play important roles in the models. Molecular docking studies predict the binding mode and the interactions between the ligand and the receptor protein. Molecular dynamics simulations results reveal that the complex of the ligand and the receptor protein are stable at 300 K. All the results can provide us more useful information for our further drug design.     

Study on properties of liquid ammonia via molecular dynamics simulation based on ABEEMσπ polarisable force field

Abstract

The molecular dynamics simulation has been performed to investigate the charge distribution, structural and dynamical properties of liquid ammonia at 273 K using a polarisable force field of the atom-bond electronegativity equalisation method (ABEEMσπ). One ammonia molecule in this model has eight charge sites, one N atomic site, three H atomic sites, three N–H bond sites and one lone-pair electron site. ABEEMσπ model can present the quantitative site charges of molecular ammonias in liquid and their changing in response to their surroundings. The radial distribution functions and dynamical properties are in fair agreement with the available experimental data. The first peak of g_NN(r) appears at N–N distance of ~3.50 ± 0.05 Å where most hydrogen bonds are formed. The average coordination number of the first shell is 13.0 ± 0.1 among which a central ammonia molecule intimately connects 3 ~ 4 ammonia molecules by hydrogen bonds. The power spectrum shows the vibrations of hydrogen bonds. For a reference, a simple estimation of the average hydrogen bonding energy in liquid ammonia is 6.5 ± 0.1 kcal/mol larger than 3.8 ± 0.3 kcal/mol in dimer ammonia. Our simulation results provide more detailed information about liquid ammonia.     

Vapor–liquid equilibria and thermophysical behavior of the SPC-HW model for heavy water

The vapor–liquid coexistence curve of the simple point charge heavy-water model (SPC-HW), [J. Chem. Phys., 114, 8064–8067 (2001)] is determined by Gibbs Ensemble Monte-Carlo (GEMC) simulation. The estimated critical conditions of the model based on the Wegner-type expansion for the order parameters and the rectilinear diameter are ρ_c = 0.300 g/cc, T _c = 661 K and P _c = 156 bars. The dielectric constant determined by isothermal–isochoric molecular dynamics is underpredicted along the coexistence curve by 29–44% in comparison with the experimental values. The analysis of the orthobaric temperature dependence of the system microstructure, in terms of the three site–site radial distribution functions, indicates that the first coordination numbers for the oxygen–oxygen and the oxygen–deuterium interactions are ～4.3 ± 0.1 and ～1.9 ± 0.1 at T = 300 K, and decrease by 15 and 55%, respectively, at criticality. The dipole–dipole correlation functions show that the orientational order in heavy water is quickly lost beyond the first oxygen–oxygen coordination shell. The model's second virial coefficient is determined by Monte-Carlo integration and used to aid the interpretation of the predicted phase equilibrium results.