Investigation of moisture diffusion in cross-linked epoxy moulding compound by molecular dynamics simulation

Molecular dynamics (MD) simulation was carried out to study the moisture diffusion in cross-linked epoxy resin, with the influence of temperature, water concentration and polymer conversion taken into account. The simulation results showed that the moisture diffusion coefficients increased with the increase in temperature. And generally, with the increased moisture concentration or decreased polymer conversion, the moisture diffusion coefficients reduced. However, the moisture diffusion was strongly inhibited when the number of epoxy groups in completely reacted epoxy resins was equal to the number of water molecules.     

Hydrogen-bond dynamics of water in presence of an amphiphile,tetramethylurea: signature of confinement-induced effects

Various experimental and simulation studies have suggested that the presence of amphiphilic molecules in aqueous solutions substantially perturbs the tetrahedral hydrogen-bond (H-bond) network of neat liquid water. Such structural perturbation is expected to impact H-bond lifetime of liquid water. Tetramethylurea (TMU) is an example of an amphiphile because it possesses both hydrophobic and hydrophilic moieties. Molecular dynamics simulations of (water+TMU) binary mixtures at various compositions have been performed in order to investigate the microscopic mechanism through which the amphiphiles influence the H-bond dynamics of liquid water at room temperature. Present simulations indicate lengthening of both water–water H-bond lifetime and H-bond structural relaxation time upon addition of TMU in aqueous solution. At the highest TMU mole fraction studied, H-bond lifetime and structural relaxation time are, respectively, ～4 and ～8 times longer than those in neat water. This is comparable with the slowing down of H-bond dynamics for water molecules confined in cyclodextrin cavities. Simulated relaxation profiles are multi-exponential in character at all mixture compositions, and simulated radial distribution functions suggest enhanced water–water and water–TMU interactions upon addition of TMU. No evidence for complete encapsulation of TMU by water H-bond network has been found.     

Study on thermomechanical properties of cross-linked epoxy resin

In this paper, molecular dynamics simulation was carried out to investigate the thermomechanical properties of cross-linked epoxy resin. The glass transition temperature, coefficients of thermal and moisture expansion, mechanical property parameters and so on are studied with the influence of temperature, water concentration and polymer conversion taken into account. The simulation results were in good agreement with existing experimental data.     

Long-term molecular dynamics simulation of copper plastocyanin in water

A long molecular dynamics simulation (1.1 ns) of fully hydrated plastocyanin has been performed and analysed to relate protein dynamics to structural elements and functional properties. The solvated structure is described in detail by the analysis of H-bond network. During all the simulation, the crystal H-bond network is maintained in the beta-sheet regions, while several H-bonds are broken or formed on the external surface of the protein. To evaluate whether such changes could be due to conformational rearrangements or to solvent competition, we have examined the average number of H-bonds between protein atoms and water molecules, and the root mean square deviations from crystal structure as a function of protein residues. Protein mobility and flexibility have been examined by positional and dihedral angle rms fluctuations. Finally, cross-correlation maps have revealed the existence of correlated motions among residues connected by hydrogen bonds.     

Molecular dynamics and quantum chemical studies of solvent effects on cyclo glycylglycine and glycylalanine dipeptides

Hydrogen bond (H-bond) interactions between the two cyclo dipeptides, cyclo(glycyl-glycine) (CGG) and cyclo(glycyl-alanine) (CGA), and water have been studied using molecular dynamics (MD) and quantum chemical methods. The MD studies have been carried out on CGG and CGA in water using fixed charge force field (AMBER ff03) for over 10 ns with a MD time step of 2 fs. The results of this study show that the solvation pattern influences the conformations of the cyclo dipeptides. Following molecular simulations, post Hartree–Fock and density functional theory methods have been used to explore the molecular properties of the cyclo dipeptides in gaseous and aqueous phase environments. The self-consistent reaction field theory has been used to optimise the cyclopeptides in diethyl ether (? = 4.3) and water (? = 78.5), and the solvent effects have been analysed. A cluster of eight water molecules leads to the formation of first solvation shell of CGG and CGA and the strong H-bonding mainly contributes to the interaction energies. The H-bond interactions have been analysed by the calculation of electron density ρ(r) and its Laplacian ▽²ρ(r) at bond critical points using atoms in molecules theory. The natural bond orbital analysis was carried out to reveal the nature of H-bond interactions. In the solvated complexes, the keto carbons registered the maximum NMR chemical shifts.     

Water and polypeptide conformations in the gramicidin channel. A molecular dynamics study.

Theoretical studies of ion channels address several important questions. The mechanism of ion transport, the role of water structure, the fluctuations of the protein channel itself, and the influence of structural changes are accessible from these studies. In this paper, we have carried out a 70-ps molecular dynamics simulation on a model structure of gramicidin A with channel waters. The backbone of the protein has been analyzed with respect to the orientation of the carbonyl and the amide groups. The results are in conformity with the experimental NMR data. The structure of water and the hydrogen bonding network are also investigated. It is found that the water molecules inside the channel act as a collective chain; whereas the conformation in which all the waters are oriented with the dipoles pointing along the axis of the channel is a preferred one, others are also accessed during the dynamics simulation. A collective coordinate involving the channel waters and some of the hydrogen bonding peptide partners is required to describe the transition of waters from one configuration to the other.     

DFT study of water adsorption on lignite molecule surface

High moisture content is a main characteristic of low-rank coal, such as lignite. Numerous oxygen containing functional groups in lignite make it represent some special properties, and these functional groups affect the adsorption mechanisms of water molecules on lignite surface. This study reports some typical water?·?·?·?lignite conformations, along with a detailed analysis of the geometry, electrostatic potential distribution, reduced density gradient of interaction, and interaction energy decomposition. The results show that water molecules tend to aggregate around functional groups, and hydrogen bonds play a dominant role in the interaction. The adsorption energy of water cluster on lignite surface is larger than that of isolated water molecule, a good linear relationship between the interaction distance and adsorption energy of layers has been found. Since water is a polar molecule, the local minima and maxima of electrostatic potential in conformations increase along with more water adsorbing on lignite surface. Reduced density gradient analysis shows that H-bonds, van der Waals interaction, and a little steric make up the interaction between water cluster and lignite molecule. In these studied conformations which mainly are H-bond complexes, electrostatic and exchange repulsion play a dominant role, whereas polarization and dispersion make relatively small contribution to the interaction. Attractive and repulsive interaction both affect the stability of water?·?·?·?lignite conformations.     

Recognition dynamics of trinuclear copper cluster and associated histidine residues through conserved or semi-conserved water molecules in human Ceruloplasmin: The involvement of aspartic and glutamic acid gates

Human Ceruloplasmin belongs to the family of multi-copper oxidases and it is involved in different physiological processes, copper ion transport, iron metabolism, iron homeostasis, and biogenic amine metabolism. MD-simulation studies have indicated the higher hydrophilic susceptibility of the trinuclear copper cluster in native CP compared to its oxygen bound form. The copper (T2/T3) atom Cu3047 of the cluster, which is close to T1 copper center Cu3052 (~13 Å) has a higher affinity for water molecules compared to other copper centers. The water molecules of W3, W4, W5, W9, and W12 conserved water sites are coordinated to Cu3047, where W3, W9, and W12 centers are found to play some crucial role in the stabilization of native trinuclear copper cluster. The hydrogen bonding interaction of Asp169, Glu112, Asp995, and Glu1032 residues with the copper-bound conserved water molecules (W3, W4, W5, W10, and W11) in native CP is observed to be unique. The conformational flexibility of Asp169 and Glu112 and their association with the copper-bound water molecules, but the absence of such interaction in O₂-bound simulated structure of the enzyme is indicating some plausible rational on the role of these acidic residues in the gating of O₂ molecule in the native trinuclear Cu cluster of CP. The simulation results may shade some new light on the biochemistry/chemistry of CP, specially on the hydration dynamics of the trinuclear copper cluster.     

Microstructures,interactions and dynamics properties studies of N-methyldiethanolamine + guanidinium triflate ionic liquid + water tertiary system at the standard temperature

Molecular dynamics simulations with an all-atom force field have been carried out in order to understand the phase equilibrium behaviour of ternary aqueous mixtures containing guanidinium triflate ionic liquid [gua][OTf] and water mixed with N-methyldiethanolamine (MDEA) in different function composition at the standard temperature of 298.15 K. A very good numerical agreement has been obtained for the prediction of the mixture densities. The analysis of structural and dynamic properties showed that the molecular level of ternary mixtures is slightly affected by the presence of MDEA and [gua][OTf] molar fractions. For MDEA–water interactions in [gua][OTf] media, we found that MDEA prefers to be surrounded by water molecules rather than by MDEA molecules even at a high MDEA molar fraction. While for [gua][OTf]–water interaction in MDEA media, as [gua][OTf] molar fraction increases, water molecules replace counterions in the coordination shell of both ions, thus weakening their interaction. On the other hand, for MDEA–[gua][OTf] interactions in water media, we have found that as the molar fraction of [gua][OTf] increases, a sulfonate group from anion appears to have a stronger association by making hydrogen bonding with MDEA molecules. The chemical process using ionic liquids (ILs) as solvents is commonly limited by their high viscosity. Based on their physical properties such as viscosities, these ternary solvents can be applied in natural gas industry, such as removing carbon dioxide using aqueous MDEA and IL at high pressure.     

Understanding adsorption-desorption dynamics of BMP-2 on hydroxyapatite (001) surface

The interaction between protein molecules and the hydroxyapatite (HAP) crystal is an important research topic in many fields. However, the nature of their noncovalent bonding is still not clear at the atomic level. In this work, molecular dynamics simulation, steered molecular dynamics simulation, and quantum chemistry calculations were used to study the adsorption-desorption dynamics of BMP-2 on HAP (001) surface. The results suggest that there are three types of functional groups through which BMP-2 can interact with HAP crystallite, and they are -OH, -NH(2), and -COO(-). Based on the different orientations of protein, each might interact with HAP crystallite individually, or, two or three of them can work cooperatively. Concerning the mechanisms of interaction, it is found that the water-bridged H-bond plays an important role, which is the main force for groups without net charges. If there were more than one set of adsorption groups for a certain orientation of protein, the adsorption-desorption process would likely be stepwise. On the contrary, if there were only one set, there would be only the key-adsorption period. The results of density functional theory calculations confirm the actual existence of this type of water-bridged H-bond. Furthermore, it is also found that the CHARMM27 force field could provide correct structural information qualitatively, although the data are slightly different from those obtained by UB3LYP/6-31G* method.     

On the properties of aqueous amide solutions through classical molecular dynamics simulations

The structure and dynamics of infinitely diluted aqueous amide solutions is studied for 13 compounds in the NVT ensemble using classical molecular dynamics simulations. The aim of this work is to provide valuable insights into the effect of amides on liquid water properties in order to understand the amides role in the kinetic inhibition of clathrate hydrate formation in natural gas mixtures. The OPLS-AA forcefield is used to describe the amides, with parameters obtained through fitting of computed B3LYP/6-311++g^{* *} data when not available in the literature, and the SPC-E model is applied for water molecules. Structural properties of the solutions are analyzed via calculated radial distribution functions and dynamic properties are studied with the computed mean square displacements and velocity autocorrelation functions. Most of the studied compounds show a remarkable structuring effect on the surrounding water with strong interactions resulting from hydrogen bonding between solute and solvent molecules. Hydrophobic and hydrophilic synergistic effects influence the amide–water interaction and the properties of the water solvation shells around amides.     

A new gating site in human aquaporin-4: Insights from molecular dynamics simulations

Aquaporin-4 (AQP4) is the predominant water channel in different organs and tissues. An alteration of its physiological functioning is responsible for several disorders of water regulation and, thus, is considered an attractive target with a promising therapeutic and diagnostic potential. Molecular dynamics (MD) simulations performed on the AQP4 tetramer embedded in a bilayer of lipid molecules allowed us to analyze the role of spontaneous fluctuations occurring inside the pore. Following the approach by Hashido et al. [Hashido M, Kidera A, Ikeguchi M (2007) Biophys J 93: 373–385], our analysis on 200 ns trajectory discloses three domains inside the pore as key elements for water permeation. Herein, we describe the gating mechanism associated with the well-known selectivity filter on the extracellular side of the pore and the crucial regulation ensured by the NPA motifs (asparagine, proline, alanine). Notably, on the cytoplasmic side, we find a putative gate formed by two residues, namely, a cysteine belonging to the loop D (C178) and a histidine from loop B (H95). We observed that the spontaneous reorientation of the imidazole ring of H95 acts as a molecular switch enabling H-bond interaction with C178. The occurrence of such local interaction seems to be responsible for the narrowing of the pore and thus of a remarkable decrease in water flux rate. Our results are in agreement with recent experimental observations and may represent a promising starting point to pave the way for the discovery of chemical modulators of AQP4 water permeability.     

Binding dynamics of two water molecules constrained within the scytalone dehydratase binding pocket

Two water molecules reside between inhibitors and active site residues of scytalone dehydratase. A molecular dynamics study is consistent with one water molecule binding less tightly than the other. Inhibitor binding studies with site-directed mutants indicate that the hydrogen bonding network around the less mobile water molecule contributes much greater binding energy than that around the more mobile one.     

Modeling and prediction of the mixed-mode retention mechanisms for puerarin and its analogues on n-octylamine modified poly(glycidyl methacrylate-co-ethylene glycol dimethacrylate) monoliths

n-Octylamine modified poly(glycidyl methacrylate-co-ethylene glycol dimethacrylate) (poly(GMA-co-EGDMA)) monoliths were prepared for the rapid screening and determination of puerarin content of a crude extract Radix puerariae. The mixed-mode retention mechanisms for puerarin and its analogues on n-octylamine modified monoliths were investigated using a variety of solvent systems, chromatographic evaluation and molecular dynamics (MDs) modeling. The equilibrated conformations between cross-linked polymers and target molecules were obtained from MD modeling. Both the polymer skeleton and functional groups played important roles in the recognition process. The cross-linker formed a structural network skeleton, in which recognition cavities were formed surrounded by functional groups. The polymer network structures provided good interaction access for isoflavones. The active groups recognized isoflavones by both intermolecular hydrogen bonding and hydrophobic interaction. The interaction energies and retention factors between polymers and target molecules were also evaluated and compared. A higher value of interaction energy corresponded to a higher value of retention factor. The potential of using modeling technology for predicting the chromatographic performances of target molecules was explored.     

Anti-insulin antibody structure and conformation. II. Molecular dynamics with explicit solvent.

Molecular dynamics at 300 K was used as a conformation searching tool to analyze a knowledge-based structure prediction of an anti-insulin antibody. Solvation effects were modeled by packing water molecules around the antigen binding loops. Some loops underwent backbone and side-chain conformational changes during the 95-ps equilibration, and most of these new, lower potential energy conformations were stable during the subsequent 200-ps simulation. Alterations to the model include changes in the intraloop, main-chain hydrogen bonding network of loop H3, and adjustments of Tyr and Lys side chains of H3 induced by hydrogen bonding to water molecules. The structures observed during molecular dynamics support the conclusion of the previous paper that hydrogen bonding will play the dominant role in antibody-insulin recognition. Determination of the structure of the antibody by x-ray crystallography is currently being pursued to provide an experimental test of these results. The simulation appears to improve the model, but longer simulations at higher temperatures should be performed.     

Dynamic water networks in cytochrome C oxidase from Paracoccus denitrificans investigated by molecular dynamics simulations

We present a molecular dynamics study of cytochrome c oxidase from Paracoccus denitrificans in the fully oxidized state, embedded in a fully hydrated dimyristoylphosphatidylcholine lipid bilayer membrane. Parallel simulations with different levels of protein hydration, 1.125 ns each in length, were carried out under conditions of constant temperature and pressure using three-dimensional periodic boundary conditions and full electrostatics to investigate the distribution and dynamics of water molecules and their corresponding hydrogen-bonded networks inside cytochrome c oxidase. The majority of the water molecules had residence times shorter than 100 ps, but a few water molecules are fixed inside the protein for up to 1.125 ns. The hydrogen-bonded network in cytochrome c oxidase is not uniformly distributed, and the degree of water arrangement is variable. The average number of solvent sites in the proton-conducting K- and D-pathways was determined. In contrast to single water files in narrow geometries we observe significant diffusion of individual water molecules along these pathways. The highly fluctuating hydrogen-bonded networks, combined with the significant diffusion of individual water molecules, provide a basis for the transfer of protons in cytochrome c oxidase, therefore leading to a better understanding of the mechanism of proton pumping.     

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

Vibrational Properties of Cyclodextrin�CWater Solutions Investigated by Low-Frequency Raman Scattering: Temperature and Concentration Effects

Depolarized low-frequency Raman spectra from several cyclodextrin–water solutions have been investigated as a function of both temperature and macrocycle concentration. The differences between the vibrational spectra of solutions and pure water have been discussed, focusing the attention on the modifications of the vibrational bands assigned to the H-bond bending and stretching intermolecular modes of water. These features are in turn related to the structural changes occurring in the H-bonded water molecules allowing us to evince a destructuring effect on the tetrahedral hydrogen bonding arrangements induced in solution by increasing temperature and solute concentration.     

Structure and dynamics of primary hydration shell of phosphatidylcholine bilayers at subzero temperatures.

Deuterium NMR relaxation and intensity measurements of the 2H-labeled H2O/dimyristoyl phosphatidylcholine bilayer were performed to understand the molecular origin of the freezing event of phospholipid headgroup and the structure and dynamics of unfrozen water molecules in the interbilayer space at subzero temperatures. The results suggest that about one to two water molecules associated with the phosphate group freeze during the freezing event of phospholipid headgroups, whereas about five to six waters near the trimethylammonium group behave as a water cluster and remain unfrozen at temperatures as low as -70 degrees C. In addition, temperature-dependent T1 and T2 relaxation times suggest that dynamic coupling occurs not only between the phosphate group and its bound water, but also between the methyl group and the adjacent water molecules. Based on these observations, the primary hydration shell of phosphatidylcholine headgroup at subzero temperatures is suggested to consist of two distinct regions: a clathrate-like water cluster, most likely a water pentamer, near the hydrophobic methyl group, and hydration water molecules associated with the phosphate group.     

A molecular dynamics simulation study of the solvent isotope effect on copper plastocyanin

The effect of heavy water on the structure and dynamics of copper plastocyanin as well as on some aspects of the solvent dynamics at the protein-solvent interfacial region have been investigated by molecular dynamics simulation. The simulated system has been analyzed in terms of the atomic root mean square deviation and fluctuations, intraprotein H-bond pattern, dynamical cross-correlation map and the results have been compared with those previously obtained for plastocyanin in H2O (Ciocchetti et al. Biophys. Chem. 69 (1997), 185-198). The simulated plastocyanin structure in the two solvents, averaging 1 ns, is very similar along the beta-structure regions, while the most significant differences are registered, analogous to the turns and the regions likely involved in the electron transfer pathway. Moreover, plastocyanin in D2O shows an increase in the number of both the intraprotein H-bonds and the residues involved in correlated motions. An analysis of the protein-solvent coupling evidenced that D2O makes the H-bond formation more difficult with the solvent molecules for positively charged and polar residues, while an opposite trend is observed for negatively charged residues. On the other hand, the frequency of exchange of the solvent molecules involved in the protein-solvent H-bond formation is significantly depressed in D2O. The results are discussed also in connection with protein functionality and briefly with some experimental results connected with the thermostability of proteins in D2O.