Membrane-bound ARF1 peptide: interpretation of neutron diffraction data by molecular dynamics simulation methods

Adenosine diphosphate ribosylation factor-1 (ARF1) is activated by cell membrane binding of a self-folding N-terminal domain. We have previously presented four possible conformations of the membrane bound, human ARF1 N-terminal peptide in planar lipid bilayers of DOPC and DOPG (7:3 molar ratio), determined from lamellar neutron diffraction and circular dichroism data. In this paper we analyse the four possible conformations by molecular dynamics simulations. The aim of these simulations was to use MD to distinguish which of the four possible membrane bound structures was the most likely. The most likely conformation was determined according to the following criteria: (a) location of label positions on the peptide in relation to the bilayer, (b) lowest mean square displacement from the initial structure, (c) lowest system energy, (d) most peptide-lipid headgroup hydrogen bonding, (e) analysis of phi/psi angles of the peptide. These findings demonstrate the application of molecular dynamics simulations to explore neutron diffraction data.     

Exploring single-domain antibody thermostability by molecular dynamics simulation

Abstract

Single-domain antibodies also known as nanobodies are recombinant antigen-binding domains that correspond to the heavy-chain variable region of camelid antibodies. Previous experimental studies showed that the nanobodies have stable and active structures at high temperatures. In this study, the thermal stability and dynamics of nanobodies have been studied by employing molecular dynamics simulation at different temperatures. Variations in root mean square deviation, native contacts, and solvent-accessible surface area of the nanobodies during the simulation were calculated to analyze the effect of different temperatures on the overall conformation of the nanobody. Then, the thermostability mechanism of this protein was studied through calculation of dynamic cross-correlation matrix, principal component analyses, native contact analyses, and root mean square fluctuation. Our results manifest that the side chain conformation of some residues in the complementarity-determining region 3 (CDR3) and also the interaction between α-helix region of CDR3 and framework2 play a critical role to stabilize the protein at a high temperature.

Communicated by Ramaswamy H. Sarma     

Adsorption and desorption surface dynamics of gaseous adsorbate on silicate-1 by molecular dynamics simulation

The dynamics of adsorption and desorption of gaseous molecules on the external surface of a crystal and a membrane of zeolite silicate-1 is investigated by molecular dynamics simulation. The gases are argon and three hydrocarbons, n-heptane, n-butane and ethylene. The sticking coefficient and the desorption coefficient are calculated for different coverages. The results clearly show that the desorption coefficients increase with the coverage contrary to the sticking coefficients. To have a better insight in the process, the desorption and adsorption time are computed, they are very similar and they show an increase with the coverage except for n-heptane which exhibit a specific decreasing behaviour at high loading.     

Cation redistribution upon dehydration of Na58Y faujasite zeolite: a joint neutron diffraction and molecular simulation study

Using neutron scattering and Monte Carlo simulation, we investigate the distribution of cations in Na₅₈Y faujasite upon (de)hydration. We introduce a new method for the assignment of cations to specific sites in molecular simulations from their local environment. This allows us to bypass the need of the coordinates of crystallographic sites, which vary as water adsorption induces changes in the zeolite framework structure. Although the agreement between experiments and simulation is excellent at high temperature, some differences are observed below 150°C. We show that these differences are due to the presence of water and that temperature itself as well as adsorption-induced deformation of the framework play a less important role. We demonstrate the migration of sodium to sites III upon water adsorption, not observed for other Si:Al ratios.     

Exploring the effect of inhibitor AKB-9778 on VE-PTP by molecular docking and molecular dynamics simulation

Diabetic macular edema, also known as diabetic eye disease, is mainly caused by the overexpression of vascular endothelial protein tyrosine phosphatase (VE-PTP) at hypoxia/ischemic. AKB-9778 is a known VE-PTP inhibitor that can effectively interact with the active site of VE-PTP to inhibit the activity of VE-PTP. However, the binding pattern of VE-PTP with AKB-9778 and the dynamic implications of AKB-9778 on VE-PTP system at the molecular level are poorly understood. Through molecular docking, it was found that the AKB-9778 was docked well in the binding pocket of VE-PTP by the interactions of hydrogen bond and Van der Waals. Furthermore, after molecular dynamic simulations on VE-PTP system and VE-PTP ^AKB-9778 system, a series of postdynamic analyses found that the flexibility and conformation of the active site undergone an obvious transition after VE-PTP binding with AKB-9778. Moreover, by constructing the RIN, it was found that the different interactions in the active site were the detailed reasons for the conformational differences between these two systems. Thus, the finding here might provide a deeper understanding of AKB-9778 as VE-PTP Inhibitor.     

Effect of hydrogen coverage on the buckling of penta-graphene by molecular dynamics simulation

By employing a series of MD simulations, buckling behaviour of penta-graphene and functionalised penta-graphene having different hydrogen (H) coverage is presented in this study. To this end, the buckling onset strain is determined for different systems. The results reveal that the new allotrope is slightly stiffer than graphene. Moreover, the effect of H adatoms in the range 0–56% on buckling behaviour is investigated. Finally it is shown that the H coverage has an influence on stability, and ripple-type destortion of the penta-graphene under compression.     

Molecular dynamics simulation of hydration in myoglobin

This study was carried out to evaluate the stability of the 89 bound water molecules that were observed in the neutron diffraction study of CO myoglobin. The myoglobin structure derived from the neutron analysis was used as the starting point in the molecular dynamics simulation using the software package CHARMM. After solvation of the protein, energy minimization and equilibration of the system, 50 ps of Newtonian dynamics was performed. This data showed that only 4 water molecules are continously bound during the length of this simulation while the other solvent molecules exhibit considerable mobility and are breaking and reforming hydrogen bonds with the protein. At any instant during the simulation, 73 of the hydration sites observed in the neutron structure are occupied by water. © 1995 Wiley-Liss, Inc.     

Water activity coefficients in aqueous amino acid solutions by molecular dynamics simulation: 1. Force field development

New force fields for molecular dynamics (MD) simulation of aqueous zwitterionic amino acid simulations were developed. These were especially designed to calculate activity coefficient of water in amino acid solutions with high accuracy. For example, aqueous solutions of the following amino acids were considered: glycine, alanine, α-aminobutyric acid, α-aminovalerianic acid, valine and leucine. The force fields were obtained by quantum chemical calculations using B3LYP/6-31G and MP2/6-311(d,p) model theories in combination with the Merz–Kollmann–Singh scheme. To further increase the accuracy of the force field, a polarised continuum was considered in all quantum chemical calculations. Water activity coefficients obtained from MD using different all-purpose literature force fields, namely, OPLS, AMBER ff03 and GROMOS 53A6 as well as experimental data are compared with the results utilising the new force field. The new force field is shown to give better results compared with experimental data than existing force fields.     

An efficiently extendable and fine-grain parallelised molecular dynamics simulation program: Mid

A new implementation of molecular dynamics simulation is presented. We employed policy-based design to achieve static polymorphism within our simulation programs. This technique provides flexibility and extensibility without additional if-statement branching in the simulation program development. It is shown that policy-based implementation prevents computational performance degradation. We used a fine-grained domain decomposition scheme to parallelise the simulation program. The smaller size decomposition reduces the total amount of inter-processing-core communication and affords good scalability for parallel calculation of short-range forces. The calculation of long-range interactions limits the total scalability. For enhanced performance at high levels of parallelism, the calculation methods for long-range interactions should be improved.     

Analysis of dissociation process for gas hydrates by molecular dynamics simulation

The dissociation processes of methane and carbon dioxide hydrates were investigated by molecular dynamics simulation. The simulations were performed with 368 water molecules and 64 gas molecules using NPT ensembles. The TraPPE (single-site) and 5-site models were adopted for methane molecules. The EPM2 (3-site) and SPC/E models were used for carbon dioxide and water molecules, respectively. The simulations were carried out at 270 K and 5.0 MPa for hydrate stabilisation. Then, temperature was increased up to 370 K. The temperature increasing rates were 0.1–20 TK/s. The gas hydrates dissociated during increasing temperature or at 370 K. The potential models of methane molecule did not much influence the dissociation process of methane hydrate. The mechanisms of dissociation process were analysed with the coordination numbers and mean square displacements. It was found that the water cages break down first, then the gas molecules escape from the water cages. The methane hydrate was more stable than the carbon dioxide hydrate at the calculated conditions.     

Study on the interactions between diketo-acid inhibitors and prototype foamy virus integrase-DNA complex via molecular docking and comparative molecular dynamics simulation methods

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an important drug target for anti-acquired immune deficiency disease (AIDS) treatment and diketo-acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN. Due to lack of three-dimensional structures including detail interactions between HIV-1 IN and its substrate viral DNA, the drug design and screening platform remains incompleteness and deficient. In addition, the action mechanism of DKA inhibitors with HIV-1 IN is not well understood. In view of the high homology between the structure of prototype foamy virus (PFV) IN and that of HIV-1 IN, we used PFV IN as a surrogate model for HIV-1 IN to investigate the inhibitory mechanism of raltegravir (RLV) and the binding modes with a series of DKA inhibitors. Firstly, molecular dynamics simulations of PFV IN, IN-RLV, IN-DNA, and IN-DNA-RLV systems were performed for 10?ns each. The interactions and inhibitory mechanism of RLV to PFV IN were explored through overall dynamics behaviors, catalytic loop conformation distribution, and hydrogen bond network analysis. The results show that the coordinated interactions of RLV with IN and viral DNA slightly reduce the flexibility of catalytic loop region of IN, and remarkably restrict the mobility of the CA end of viral DNA, which may lead to the partial loss of the inhibitory activity of IN. Then, we docked a series of DKA inhibitors into PFV IN-DNA receptor and obtained the IN-DNA-inhibitor complexes. The docking results between PFV IN-DNA and DKA inhibitors agree well with the corresponding complex of HIV-1 IN, which proves the dependability of PFV IN-DNA used for the anti-AIDS drug screening. Our study may help to make clear some theoretical questions and to design anti-AIDS drug based on the structure of IN.     

Replica-exchange molecular dynamics simulation for understanding the initial process of amyloid peptide aggregation

We applied replica-exchange molecular dynamics simulation to five fragments of amyloid-β peptide in order to study the mechanism of amyloid fibril formation. In this study, we calculated the free energy by focusing on how to form the β-structures to obtain the dominant structures. We classify the obtained β-structures and elucidate the order of β-structure assembly.     

Extending the capabilities of targeted molecular dynamics: simulation of a large conformational transition in plasminogen activator inhibitor 1

Plasminogen activator inhibitor type 1 (PAI-1) is an inhibitor of plasminogen activators such as tissue-type plasminogen activator or urokinase-type plasminogen activator. For this molecule, different conformations are known. The inhibiting form that interacts with the proteinases is called the active form. The noninhibitory, noncleavable form is called the latent form. X-ray and modeling studies have revealed a large change in position of the reactive center loop (RCL), responsible for the interaction with the proteinases, that is inserted into a beta-sheet (s4A) in the latent form. The mechanism underlying this spontaneous conformational change (half-life = 2 h at 37 degrees C) is not known in detail. This investigation attempts to predict a transition path from the active to the latent structure at the atomic level, by using simulation techniques. Together with targeted molecular dynamics (TMD), a plausible assumption on a rigid body movement of the RCL was applied to define an initial guess for an intermediate. Different pathways were simulated, from the active to the intermediate, from the intermediate to the latent structure and vice versa under different conditions. Equilibrium simulations at different steps of the path also were performed. The results show that a continuous pathway from the active to the latent structure can be modeled. This study also shows that this approach may be applied in general to model large conformational changes in any kind of protein for which the initial and final three-dimensional structure is known.     

Understanding interaction and dynamics of water molecules in the epoxy via molecular dynamics simulation

The robust structural integrity of the epoxy plays an important role in ensuring the long-term service life of its applications, which is affected by the absorbed moisture. In order to understand the mechanism of the moisture effect, the knowledge of the interaction and dynamics of the water molecules inside the epoxy is of great interest. Molecular dynamics simulation is used in this work to investigate the structure and bonding behaviour of the water molecules in the highly cross-linked epoxy network. When the moisture concentration is low, the water molecules are well dispersed in the cross-linked structure and located in the vicinity of the epoxy functional groups, which predominantly form the hydrogen bond (H-bond) with the epoxy network, resulting in the low water mobility in the epoxy. At the high concentration, the water favourably forms the large cluster due to the predominant water–water H-bond interaction, and the water molecules diffuse primarily inside the cluster, which leads to the high water mobility and the accelerated H-bond dynamics. The variation of the bonding behaviour and dynamics of the water molecules reported here could be exploited to understand the material change and predict the long-term performance of the epoxy-based products during the intended service life.     

Investigation of the role of disulphide bond in modulating internal motions of BmK AGAP protein by molecular dynamics simulation

Elucidating structural determinants in the functional regions of toxins can provide useful knowledge for designing novel analgesic peptides. A series of 100 ns MD simulations were performed on the scorpion toxin BmK AGAP in native and disulphide bond broken states. The comparison of disulphide bond broken states with the native state showed the α-helix was found to be the key to the analgesic activity. Furthermore, our results revealed disulphide bonds have considerable influence on the functionally important essential modes of motions and the correlations between the motions of the Core domain and the C-terminal region which are involved in the analgesic activity. Therefore, we can conclude that disulphide bonds have a crucial role in modulating the function via adjusting the dynamics of scorpion toxin BmK AGAP molecule.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Protein under pressure: molecular dynamics simulation of the arc repressor

Experimental nuclear magnetic resonance results for the Arc Repressor have shown that this dimeric protein dissociates into a molten globule at high pressure. This structural change is accompanied by a modification of the hydrogen-bonding pattern of the intermolecular beta-sheet: it changes its character from intermolecular to intramolecular with respect to the two monomers. Molecular dynamics simulations of the Arc Repressor, as a monomer and a dimer, at elevated pressure have been performed with the aim to study this hypothesis and to identify the major structural and dynamical changes of the protein under such conditions. The monomer appears less stable than the dimer. However, the complete dissociation has not been seen because of the long timescale needed to observe this phenomenon. In fact, the protein structure altered very little when increasing the pressure. It became slightly compressed and the dynamics of the side-chains and the unfolding process slowed down. Increasing both, temperature and pressure, a tendency of conversion of intermolecular into intramolecular hydrogen bonds in the beta-sheet region has been detected, supporting the mentioned hypothesis. Also, the onset of denaturation of the separated chains was observed.