Molecular dynamics simulations of interaction between sub-bituminous coal and water

The high moisture content of sub-bituminous coal is associated with the interactions between coal and water. Because of complex composition and structure, the graphite surface modified by hydroxyl, carboxyl and carbonyl groups was used to represent the surface model of sub-bituminous coal according to XPS results. Density profiles for oxygen atoms and hydrogen atoms indicate that the coal surface properties affect the structural and dynamic characteristics of the interfacial water molecules. The interfacial water exhibits much more ordering than bulk water. The results of radial distribution functions, mean square displacement and local self-diffusion coefficient for water molecule related to three oxygen moieties confirmed that the water molecules prefer to absorb with carboxylic groups, and adsorption of water molecules at the hydroxy and carbonyl is similar.     

Molecular dynamics simulations of Leu-enkephalin in water and DMSO.

The structure of Leu-enkephalin (L-Enk) and Met-enkephalin (M-Enk) have frequently been studied, in particular by nuclear magnetic resonance spectroscopy. After more than 20 years of research, it was concluded that enkephalins have no preferred structure in aqueous solution, but that they may have in other solvents. We have performed molecular dynamics simulations of zwitterionic L-Enk in water, and zwitterionic as well as neutral L-Enk dimethyl sulfoxide (DMSO). In water the peptide is very flexible, although there seems to be a preference for compact conformations. In DMSO, the peptide forms a clear salt bridge in the zwitterionic form, but has no preferred conformation in the neutral form. This difference in conformation may provide an explanation for measurements in DMSO in which multiple conformations were found to exist. In this paper we introduce a new formulation for a dihedral angle autocorrelation function, and apply it to study side-chain dynamics in L-Enk. We find that the side-chain dynamics of the large Tyr and Phe residues cannot be adequately sampled in 2.0-ns simulations, while this does seem to be possible for the smaller Leu side chain.     

Molecular dynamic simulations of electric microfield distributions in a nonideal electron-positron plasma

A symmetric model of a two-component plasma is considered and the distributions of electric microfields acting on charged and neutral particles are calculated using the method of molecular dynamics at a fixed temperature of T = 30000 K and different values of the coupling parameter 0.2 ≤ Γ ≤ 1.2. Changes in these distributions with varying Γ are discussed. Special attention is paid to the behavior of the distribution tails. The behavior of these tails at a neutral point is shown to agree with the tails of the Holtsmark distribution, whereas the tails of the distribution at a charge are considerably heavier and are characterized by the exponent that varies within the range from −2.2 up to −1.8 as Γ increases.     

Molecular dynamics simulations of water within models of ion channels.

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.     

Molecular dynamics simulations

Molecular dynamics simulations have become a standard tool for the investigation of biomolecules. Simulations are performed of ever bigger systems using more realistic boundary conditions and better sampling due to longer sampling times. Recently, realistic simulations of systems as complex as transmembrane channels have become feasible. Simulations aid our understanding of biochemical processes and give a dynamic dimension to structural data; for example, the transformation of harmless prion protein into the disease-causing agent has been modeled.     

Molecular dynamics simulations of trihalomethanes removal from water using boron nitride nanosheets

Molecular dynamics simulations were performed to investigate the separation of trihalomethanes (THMs) from water using boron nitride nanosheets (BNNSs). The studied systems included THM molecules and a functionalized BNNS membrane immersed in an aqueous solution. An external pressure was applied to the z axis of the systems. Two functionalized BNNSs with large fluorinated-hydrogenated pore (F-H-pores) and small hydrogen-hydroxyl pore (H-OH-pores) were used. The pores of the BNNS membrane were obtained by passivating each nitrogen and boron atoms at the pore edges with fluorine and hydrogen atoms in the large pore or with hydroxyl and hydrogen atoms in the small pore. The results show that the BNNS with a small functionalized pore was impermeable to THM molecules, in contrast to the BNNS with a large functionalized pore. Using these membranes, water contaminants can be removed at lower cost.

Open image in new window

Graphical Abstract A snapshot of the simulation system. The BNNS membrane with the large functionalized pore is located in the middle of the box. The size of the box is 3 × 3 × 5 nm³. Green chlorine, cyan carbon, red oxygen, white hydrogen

     

Molecular dynamics simulations of a DMSO/water mixture using the AMBER force field

Due to its protective properties of biological samples at low temperatures and under desiccation, dimethyl sulfoxide (DMSO) in aqueous solutions has been studied widely by many experimental approaches and molecular dynamics (MD) simulations. In the case of the latter, AMBER is among the most commonly used force fields for simulations of biomolecular systems; however, the parameters for DMSO published by Fox and Kollman in 1998 have only been tested for pure liquid DMSO. We have conducted an MD simulation study of DMSO in a water mixture and computed several structural and dynamical properties such as of the mean density, self-diffusion coefficient, hydrogen bonding and DMSO and water ordering. The AMBER force field of DMSO is seen to reproduce well most of the experimental properties of DMSO in water, with the mixture displaying strong and specific water ordering, as observed in experiments and multiple other MD simulations with other non-polarizable force fields.

Open image in new window

Graphical abstract Hydration structure within hydrogen-bonding distance around a DMSOmolecule

     

Molecular dynamic simulations of the metallo-beta-lactamase from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor

The beta-lactam-based antibiotics are among the most prescribed and effective antibacterial agents. Widespread use of these antibiotics, however, has created tremendous pressure for the emergence of resistance mechanisms in bacteria. The most common cause of antibiotic resistance is bacterial production of actamases that efficiently degrade antibiotics. The metallo-beta-lactamases are of particular clinical concern due to their transference between bacterial strains. We used molecular dynamics (MD) simulations to further study the conformational changes that occur due to binding of an inhibitor to the dicanzinc metallo-beta-lactamase from Bacteroides fragilis. Our studies confirm previous findings that the major flap is a major source of plasticity within the active site, therefore its dynamic response should be considered in drug development. However, our results also suggest the need for care in using MD simulations in evaluating loop mobility, both due to relaxation times and to the need to accurately model the zinc active site. Our study also reveals two new robust responses to ligand binding. First, there are specific localized changes in the zinc active site—a local loop flip—due to ligand intercalation that may be critical to the function of this enzyme. Second, inhibitor binding perturbs the dynamics throughout the protein, without otherwise perturbing the enzyme structure. These dynamic perturbations radiate outward from the active site and their existence suggests that long-range communication and dynamics may be important in the activity of this enzyme. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Molecular dynamic simulations of the N-terminal receiver domain of NtrC reveal intrinsic conformational flexibility in the inactive state

The N-terminal receiver domain of NtrC is the molecular switch in the two-component signal transduction. It is the first protein where structures of both the active (phosphyroylated) and inactive (unphosphyroylated) states are determined experimentally. Phosphorylation of the NtrC at the active site induces large structural change. NMR experiments suggested that the wild type unphosphorylated NtrC adopts both the active and the inactive conformations and the phosphorylation stabilizes the active conformations. We applied free (unconstrained) molecular dynamic (MD) simulation to examine the intrinsic flexibilities and stabilities of the NtrC receiver domain in both the active and inactive conformations. Molecular dynamic simulations showed that the inactive state of NtrC receiver domain is more flexible than the active state. There were large movements in helix 4 and loop beta3-alpha3 which coincide with major structural differences between the inactive and active states. We observed large root-mean-square deviations from the initial starting structure and the large root-mean-square fluctuations during MD simulation for the inactive state. We then investigated the activation pathway with Targeted MD simulation. We show that the intrinsic flexibility in the loop beta3-alpha3 plays an important role in triggering the conformational change. Phosphorylation at the active site may serve to stabilize the conformational change. These results together suggest that the unphosphorylated NtrC receiver domain could be involved in a conformational equilibrium between two different states.     

Molecular dynamics simulations of the unfolding of an alpha-helical analogue of ribonuclease A S-peptide in water

Molecular dynamics simulations of the S-peptide analogue AETAAAKFLREHMDS have been conducted in aqueous solution for 300 ps at 278 K and for 500 ps in two different runs at 358 K. The results show agreement with experimental observations in that at low temperature, 5 degrees C, the helix is stable, while unfolding is observed at 85 degrees C. In the low-temperature simulation a solvent-separated ion pair was formed between Glu-2 and Arg-10, and the side chain of His-12 reoriented toward the C-terminal end of the alpha-helix. Detailed analyses of the unfolding pathways at high temperature have also revealed that the formation or disappearance of main-chain helical hydrogen bonds occurs frequently through an alpha in equilibrium with 3(10) in equilibrium with no hydrogen bond sequence.     

Molecular simulations of rhodopsin tetrameter

Di/oligomerization of G-protein coupled receptors (GPCRs) is well established, however very little is known regarding the interaction details. Current paper presents results of molecular dynamics simulations of theoretical model of rhodopsin tetramer with transducine (Gt) in lipid bilayer. Ligand-protein and receptor-receptor interactions have been proposed.     

Molecular mechanics and dynamics simulations of various dispersant models on the water surface (001)

Five dispersant-molecule models of succinimide, acrylate, imide, phenylsulfonic and salicyl were used to study their interactions with the water surface (001). The interaction energy, molecular configuration, charge distribution and radial distribution function (RDF) curve for each of the dispersant molecules were analyzed from the molecular mechanics (MM) and molecular dynamics (MD) simulation results. It can be seen that the system energies, mostly electrostatic and hydrogen bond energies, were reduced significantly when the dispersant molecules interacted with the water surface. The hydrophilic group of a dispersant molecule can attach itself to the water surface firmly and reach for a stable energy-minimized configuration, which is helpful to the dispersants' dispersancy. The influence exerted by the hydrophobic group of the dispersant molecule, which was the substituted hydrocarbon chain of n-octadecanyl in this paper, is discussed in comparison with the naked polar headgroup. Steric configuration, charge distribution and substitute hydrocarbon chain of the dispersant molecule influenced the interaction between dispersants and polar water surface.     

Investigation into the diffusion of water into HEMA-co-MOEP hydrogels

Cross-linked homopolymers and copolymers of 2-hydroxyethyl methacrylate, HEMA, and ethylene glycol methacrylate phosphate, MOEP, have been synthesized, and the diffusion of water into these systems has been investigated. Only polymers with 0-20 mol % MOEP exhibited ideal swelling behavior as extensive fracturing occurred in the systems with greater than 20 mol % MOEP as the polymers began to swell during water sorption. Gravimetric studies were used in conjunction with magnetic resonance imaging of the diffusion front to elucidate the diffusion mechanism for these systems. In the case of the cross-linked HEMA homopolymer gels, the water transport mechanism was determined to be concentration-independent Fickian diffusion. However, as the fraction of MOEP in the network increased, the transport mechanism became increasingly exponentially concentration-dependent but remained Fickian until the polymer consisted of 30 mol % MOEP where the water transport could no longer been described by Fickian diffusion.     

Molecular dynamic simulations analysis of ritonavir and lopinavir as SARS-CoV 3CL(pro) inhibitors

Since the emergence of the severe acute respiratory syndrome (SARS) to date, neither an effective antiviral drug nor a vaccine against SARS is available. However, it was found that a mixture of two HIV-1 proteinase inhibitors, lopinavir and ritonavir, exhibited some signs of effectiveness against the SARS virus. To understand the fine details of the molecular interactions between these proteinase inhibitors and the SARS virus via complexation, molecular dynamics simulations were carried out for the SARS-CoV 3CL^pro free enzyme (free SARS) and its complexes with lopinavir (SARS-LPV) and ritonavir (SARS-RTV). The results show that flap closing was clearly observed when the inhibitors bind to the active site of SARS-CoV 3CL^pro. The binding affinities of LPV and RTV to SARS-CoV 3CL^pro do not show any significant difference. In addition, six hydrogen bonds were detected in the SARS-LPV system, while seven hydrogen bonds were found in SARS-RTV complex.     

Molecular dynamics simulations in photosynthesis

Photosynthesis Research - Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular...