Molecular dynamics simulation of the graphene–water interface: comparing water models

In this work, different water models (i.e. SPC/E, TIP3P, TIP4P/2005, TIP5P, SPC/Fw, TIP4P/2005f and SWM4_DP) are implemented to simulate water on neutral, negatively charged and positively charged graphene. In all cases ambient conditions are considered. Structural and dynamical properties for water are calculated to quantify the differences among various water models. The results show that SPC/E, TIP4P/2005, SPC/Fw, TIP4P/2005f and SWM4_DP water models yield a similar structure for interfacial water on graphene, whether it is neutral, negatively charged or positively charged. TIP5P is the model whose predictions for the structure of the interface deviate the most from those of the other models. Although qualitatively the results are for the most part similar, a large quantitative variation is observed among the dynamical properties predicted when various water models are implemented. Although experimental data are not available to discriminate the most/least accurate of the model predictions, our results could be useful for comparing results for interfacial water obtained implementing different models. Such critical comparison will benefit practical applications such as the development of energy-storage and water-desalination devices (e.g. electric double-layer capacitors), among others.     

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.     

Studies of ice melting using molecular dynamics

In this paper, the melting of ice 1h is studied using molecular dynamics (MD). Common potential functions employed in the MD simulations include SPC/E, TIP4P, TIP5P, TIP4/ice and TIP5P/E. We first conducted melting of ice bulks and then studied the melting speed of the ice/water interface during ice melting. It is found that various potentials result in similar ice-melting phenomena. The result is compared with the analytical solution for phase change problem. We also studied size effects and temperature gradient effects on ice melting.     

Molecular dynamics of water in the neighborhood of aquaporins

Present knowledge obtained by molecular dynamics (MD) simulation studies regarding the dynamics of water, both in the vicinity of biological membranes and within the proteinaceous water channels, also known as aquaporins (AQPs), is reviewed. A brief general summary of the water models most extensively employed in MD simulations (SPC, SPC/E, TIP3P, TIP4P), indicating their most relevant pros and cons, is likewise provided. Structural considerations of water are also discussed, based on different order parameters, which can be extracted from MD simulations as well as from experiments. Secondly, the behaviour of water in the neighbourhood of membranes by means of molecular dynamics simulations is addressed. Consequently, the comparison with previous experimental evidence is pointed out. In living cells, water is transported across the plasma membrane through the lipid bilayer and the aforementioned AQPs, which motivates this review to focus mostly on MD simulation studies of water within AQPs. Relevant contributions explaining peculiar properties of these channels are discussed, such as selectivity and gating. Water models used in these studies are also summarised. Finally, based on the information presented here, further MD studies are encouraged.     

Molecular Dynamics Simulation of Limiting Conductance for Li + Ion in Supercritical Water using Polarizable Models

We report results of molecular dynamics simulations of the limiting conductance of Li + ion in ambient water and in supercritical water using polarizable models for water and Li + . The limiting conductances of Li + in ambient water calculated from mean square displacement (MSD) using four points transferable intermolecular potential model (TIP4P), extended simple point charge model (SPC/E), and revised polarizable model 1 (RPOL1) are larger than the experimental value. The behavior of the limiting conductance of Li + in supercritical water using the RPOL models results in good agreement with experimental results for the limiting conductance of LiCl. The agreement of the RPOL1 model with the experimental results is much better than the RPOL2 model in the higher-density regime, whereas that of the RPOL2 model is much better than the RPOL1 model in the lower-density regime. Using the RPOL models (in contrast to the SPC/E model), the number of hydration water molecules around Li + is the dominating contributor to the limiting conductance in the higher-density regime. In agreement with the SPC/E model, the interaction strength between Li + and the hydration water molecules is a non-factor in the lower-density region since the potential energy per hydration water molecule decreases with decreasing water density at the lowest water densities.     

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.     

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.     

Global Potential Energy Minima of SPC/E Water Clusters without and with Induced Polarization Using a Genetic Algorithm

The global minimum potential energy structures of water clusters, (H₂O) _n , n = 2-14, have been calculated for the SPC/E (Simple Point Charge/Extended) model and a recent fluctuating charge version of the latter using a simple genetic algorithm. The SPC/E cluster geometries are in good agreement with previous TIP3P (Transferable Intermolecular Potential-3 Point) and TIP4P (Transferable Intermolecular Potential-4 Point) calculations as well as the interpretation of experimental measurements. In contrast to this, the polarizable version of the SPC/E model, which is based on the fluctuating charge approach, deviates rather strongly for n=6 with few exceptions. However, comparing the polarizable model to ab initio results for identical cluster geometries we find reasonable agreement for the magnitude of the average molecular dipole moment, the corresponding energy per molecule, and the average oxygen-oxygen distance as functions of n.Electronic Supplementary Material available.     

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.     

Study of conformational mobility of insulin,proinsulin, and insulin-like growth factors

The goal of the work was to perform by the method of molecular dynamics a comparative analysis of conformational mobility of evolutionary related peptides—insulin, proinsulin, IGF1, and IGF2. The proinsulin molecule has been shown to have the highest mobility, whereas IGF1—the lowest. Rotation radius (Rg) of insulin, IGF1, and IGF2 changes insignificantly, Rg of proinsulin decreases by reaching plateau after 6000 ps. The graphs of the mean square deviations (RMSD) from initial positions for A- and B-domains are practically identical, which indicates integrity of the carcass formed by A- and B-domains. The proinsulin C-domain behaves sufficiently chaotically. The IGF1 and IGF2 C-domains form ordered structures resembling horseshoe and an elongated hairpin. The D-domain makes the greatest contribution to the IGF2 mobility, but remains virtually immobile in IGF1, which might be the cause of high IGF1 stability. The obtained data can be used in development of new effective drugs for treatment of diabetes mellitus, which are based on the principle of evolutionary relationship of peptides.     

Structural models of viral insulin-like peptides and their analogs

The insulin receptor (IR), the insulin-like growth factor-1 receptor (IGF1R), and the insulin/IGF1 hybrid receptors (hybR) are homologous transmembrane receptors. The peptide ligands, insulin and IGF1, exhibit significant structural homology and can bind to each receptor via site-1 and site-2 residues with distinct affinities. The variants of the Iridoviridae virus family show capability in expressing single-chain insulin/IGF1 like proteins, termed viral insulin-like peptides (VILPs), which can stimulate receptors from the insulin family. The sequences of VILPs lacking the central C-domain (dcVILPs) are known, but their structures in unbound and receptor-bound states have not been resolved to date. We report all-atom structural models of three dcVILPs (dcGIV, dcSGIV, and dcLCDV1) and their complexes with the receptors (μIR, μIGF1R, and μhybR), and probed the peptide/receptor interactions in each system using all-atom molecular dynamics (MD) simulations. Based on the nonbonded interaction energies computed between each residue of peptides (insulin and dcVILPs) and the receptors, we provide details on residues establishing significant interactions. The observed site-1 insulin/μIR interactions are consistent with previous experimental studies, and a residue-level comparison of interactions of peptides (insulin and dcVILPs) with the receptors revealed that, due to sequence differences, dcVILPs also establish some interactions distinct from those between insulin and IR. We also designed insulin analogs and report enhanced interactions between some analogs and the receptors.     

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5-C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.     

Molecular dynamics simulation of a decasaccharide fragment of heparin in aqueous solution

Molecular dynamics (MD) simulations on heparin-water-sodium systems were carried out in order to establish a simulation protocol able to represent heparin solution conformation under physiological conditions. Atomic charges suitable for heparin oligosaccharides were obtained from ab initio quantum-mechanical computations, at the 6-31G(**) level. The GROMACS forcefield, the SPC, and SPC/E water models were employed. Also heparin was simulated with IdoA residues in 1C(4) or 2S(0) conformational states. The results of the performed MD simulations are in agreement with the available experimental data, suggesting that this approach can be applied for the study of heparin interactions with its target proteins and thus play a role in the development of new antithrombotic agents.     

Long time dynamics of Met-enkephalin: comparison of explicit and implicit solvent models

Met-enkephalin is one of the smallest opiate peptides. Yet, its dynamical structure and receptor docking mechanism are still not well understood. The conformational dynamics of this neuron peptide in liquid water are studied here by using all-atom molecular dynamics (MD) and implicit water Langevin dynamics (LD) simulations with AMBER potential functions and the three-site transferable intermolecular potential (TIP3P) model for water. To achieve the same simulation length in physical time, the full MD simulations require 200 times as much CPU time as the implicit water LD simulations. The solvent hydrophobicity and dielectric behavior are treated in the implicit solvent LD simulations by using a macroscopic solvation potential, a single dielectric constant, and atomic friction coefficients computed using the accessible surface area method with the TIP3P model water viscosity as determined here from MD simulations for pure TIP3P water. Both the local and the global dynamics obtained from the implicit solvent LD simulations agree very well with those from the explicit solvent MD simulations. The simulations provide insights into the conformational restrictions that are associated with the bioactivity of the opiate peptide dermorphin for the delta-receptor.     

Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins

Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 μs. Four combinations of protein and water force fields were tested: ff19SB/OPC , ff19SB/TIP4P-D , ff03CMAP/TIP4P-D , and a99SB-disp/TIP4P-disp , with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 μs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and ³J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.     

Leptin directly controls secretory activity of human ovarian granulosa cells: possible inter-relationship with the IGF/IGFBP system

AIMS: The aim of our in vitro studies was to understand the role of leptin and the insulin-like growth factor I/insulin-like growth factor protein (IGF/IGFBP) system in controlling human ovarian function. METHODS: We studied the action of leptin (0, 1, 10, or 100 ng/ml) and immunoneutralization of IGF-I using specific antiserum (0.1%) on the release of progesterone (P), estradiol (E), oxytocin (OT), IGF-I, IGFBP-3, and prostaglandins F (PGF) by these cells using radioimmunoassay/immunoradiometric assay. RESULTS: It was found that leptin stimulated the secretion of OT, IGFBP-3, and PGF. It suppressed the secretion of E and IGF-I, but not P, into the medium. The addition of antiserum against IGF-I decreased IGF-I output, increased P, OT, IGFBP-3, and PGF secretion, and had no effect on E release. Immunoneutralization of IGF-I also prevented or reversed the effects of leptin on P, E, IGF-I, IGFBP-3, PGF, but not on OT. CONCLUSIONS: These observations (1) demonstrate that leptin directly controls the secretory activity of human ovarian cells, (2) confirm the involvement of IGF-I in the regulation of ovarian cells, and (3) suggest an inter-relationship between leptin and the IGF/IGFBP system in the control of these functions and the involvement of IGF/IGFBP system in mediating leptin action on the ovary.     

Oocyte-mediated suppression of follicle-stimulating hormone- and insulin-like growth factor-induced secretion of steroids and inhibin-related proteins by bovine granulosa cells in vitro: possible role of transforming growth factor alpha

The objective was to investigate the potential role of the oocyte in modulating proliferation and basal, FSH-induced and insulin-like growth factor (IGF)-induced secretion of inhibin A (inh A), activin A (act A), follistatin (FS), estradiol (E(2)), and progesterone (P(4)) by mural bovine granulosa cells. Cells from 4- to 6-mm follicles were cultured in serum-free medium containing insulin and androstenedione, and the effects of ovine FSH and IGF analogue (LR3-IGF-1) were tested alone and in the presence of denuded bovine oocytes (2, 8, or 20 per well). Medium was changed every 48 h, cultures were terminated after 144 h, and viable cell number was determined. Results are based on combined data from four independent cultures and are presented for the last time period only when responses were maximal. Both FSH and IGF increased (P < 0.001) secretion of inh A, act A, FS, E(2), and P(4) and raised cell number. In the absence of FSH or IGF, coculture with oocytes had no effect on any of the measured hormones, although cell number was increased up to 1.8-fold (P < 0.0001). Addition of oocytes to FSH-stimulated cells dose-dependently suppressed (P < 0.0001) inh A (6-fold maximum suppression), act A (5.5-fold), FS (3.6-fold), E(2) (4.6-fold), and P(4) (2.4-fold), with suppression increasing with FSH dose. Likewise, oocytes suppressed (P < 0.001) IGF-induced secretion of inh A, act A, FS, and E(2) (P < 0.05) but enhanced IGF-induced P(4) secretion (1.7-fold; P < 0.05). Given the similarity of these oocyte-mediated actions to those we observed previously following epidermal growth factor (EGF) treatment, we used immunocytochemistry to determine whether bovine oocytes express EGF or transforming growth factor (TGF) alpha. Intense staining with TGFalpha antibody (but not with EGF antibody) was detected in oocytes both before and after coculture. Experiments involving addition of TGFalpha to granulosa cells confirmed that the peptide mimicked the effects of oocytes on cell proliferation and on FSH- and IGF-induced hormone secretion. These experiments indicate that bovine oocytes secrete a factor(s) capable of modulating granulosa cell proliferation and responsiveness to FSH and IGF in terms of steroidogenesis and production of inhibin-related peptides, bovine oocytes express TGFalpha but not EGF, and TGFalpha is a prime candidate for mediating the actions of oocytes on bovine granulosa cells.