Applying the stochastic difference equation to DNA conformational transitions: a study of B-Z and B-A DNA transitions

Despite the existence of numerous models to account for the B-Z DNA transition, experimenters have not yet arrived at a conclusive answer to the structural and dynamical features of the B-Z transition. By applying the stochastic difference equation to simulate the B-Z DNA transition, we have shown that the stretched intermediate model of the B-Z transition is more probable than other B-Z transition models such as the Harvey model. This is accomplished by comparing potential energy profiles of various B-Z DNA transition models and calculating relative probabilities based on the stochastic difference equation with respect to length (SDEL) formalism. The results garnered in this article allow for new approaches in determining the structural transition of B-DNA to Z-DNA experimentally. We have also simulated the B-A DNA transition using the stochastic difference equation. Unlike the B-Z DNA transition, the mechanism for the B-A DNA transition is well established. The variation in the pseudorotation angle during the transition is in good agreement with experimental results. Qualitative features of the simulated B-A transition also agree well with experimental data. The SDEL approach is thus a suitable numerical technique to compute long-time molecular dynamics trajectory for DNA molecules.     

Overview of simulation studies on the enzymatic activity and conformational dynamics of the GTPase Ras

Over the last 40 years, we have learnt a great deal about the Ras onco-proteins. These intracellular molecular switches are essential for the function of a variety of physiological processes, including signal transduction cascades responsible for cell growth and proliferation. Molecular simulations and free energy calculations have played an essential role in elucidating the conformational dynamics and energetics underlying the GTP hydrolysis reaction catalysed by Ras. Here we present an overview of the main lessons from molecular simulations on the GTPase reaction and conformational dynamics of this important anti-cancer drug target. In the first part, we summarise insights from quantum mechanical and combined quantum mechanical/molecular mechanical simulations as well as other free energy methods and highlight consensus viewpoints as well as remaining controversies. The second part provides a very brief overview of new insights emerging from large-scale molecular dynamics simulations. We conclude with a perspective regarding future studies of Ras where computational approaches will likely play an active role.     

持续性催化酶类机理研究及其分子动力学模拟

酶分子催化机理研究是生命科学研究领域一个重要的问题.近80年来,过渡态理论在解释酶催化机理问题上占据了主导地位,结合热力学循环、锁钥学说、诱导契合学说以及酶活性中心柔性学说等理论,可以很好地解释多种酶分子的催化过程.近年来,随着蛋白质结构解析方法、单分子分析检测技术及计算机模拟技术的发展,人们对酶分子催化机理的认识愈加深刻.但持续性催化酶类的催化动力研究表明,过渡态理论的解释并不充分.本文对酶催化机理研究的相关进展进行了综述,并针对持续性酶类催化动态过程的特点提出了可能的研究方向及可行的研究方法.     

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.     

On the Approximation of Solvent Effects on the Conformation and Dynamics of Cyclosporin A by Stochastic Dynamics Simulation Techniques

Abstract

The molecular simulation technique of stochastic dynamics (SD) is tested by application to the immunosuppressive drug cyclosporin A (CPA). Two stochastic dynamics simulations are performed, one (SD_CCl4 ) with atomic friction coefficients proportional to the viscosity of the nonpolar solvent CCl₄, and one (SD_H2O) with atomic friction coefficients corresponding to an aqueous solution. The atomic friction coefficients are also taken proportional to an approximate expression for the atomic accessible surface area. The properties of both stochastic dynamics simulations are compared to those of two full molecular dynamics (MD) simulations of cyclosporin A, one in a box with 591 CCl₄ molecules, and one in a box with 632 H₂O molecules.

The properties of cyclosporin A as found in the molecular dynamics simulation in CCl₄ are well reproduced by the SD_CCl4 simulation. This indicates that the neglect of a mean force reresenting the average solvent effects on the solute is justified in the case of nonpolar solvents. For polar solvents, like water, this mean force may not be neglected. The SD_H2O simulation of cyclosporin A clearly fails to reproduce the amount of hydrogen bonding found in the molecular dynamics stimulation of cyclosporin A in water.

A comparison with a molecular dynamics simulation of cyclosporin A in vacuo shows that both the SD_CCl4 and the SD_H2O simulation come closer to the properties of the molecular dynamics simulations in CCl₄ and in H₂O than a molecular dynamics simulation in vacuo.     

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.     

Potential Therapeutic Agents on Alzheimer's Disease through Molecular Docking and Molecular Dynamics Simulation Study of Plant-Based Compounds

Globally Alzheimer's disease (AD) is a highly complex, heterogeneous, and multifactorial neurological disease. AD is categorized clinically through a steady loss in memory and progressive decline of cognitive function. So far, there is no effective cure is available for the treatment of AD. Here, we identified Plant-based compounds (PBCs) from seven therapeutic plants through pharmacophore and pharmacokinetics approaches. Subsequently, we retrieved 65 AD associated proteins by Text Mining approach .We observed the interactions between 39 PBCs with 65 AD-associated targets by using molecular docking. Further, we carried out Molecular dynamics simulation analysis to predict the steady binding of top drug-target complexes. The entire MD simulation results analysis was evidence that seven drug-target complexes consistently interacted during the in silico experiment. The top complexes were the target CHLE interacted with 2 PBCs (Pseudojujubogenin and Anahygrine), target VDAC1 interacted with Withanolide R, target THOP1 interacted with Withaolide R, target AOFB interacted with 2 PBCs (Nardostachysin and Viscosalactone B), and target ACHE interacted with the drug (12-Deoxywithastramonolide). These PBCs have stably and flexibly interacted at the protein‘s active site region. Our results suggest that these PBCs and targets are potential therapeutic candidates for molecular development in AD.     

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.     

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.     

Role of electrostatic interactions for the stability and folding behavior of cold shock protein

The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Study on the binding characteristics of hydroxylated polybrominated diphenyl ethers and thyroid transporters using the multispectral technique and computational simulation

Hydroxylated polybrominated diphenyl ethers (OH-PBDEs) are a class of toxic environmental pollutants that are persistent, bioaccumulative, and difficult to degrade. Their structure is very similar to the thyroid hormone (T4) and uses the body’s thyroid transporter (TTR) binding to interfere with the endocrine balance, disrupting the body’s normal physiological activity. According to Fourier transform infrared spectroscopy and dynamics simulation of do_dssp module analysis, there are three kinds of OH-PBDEs that can induce TTR secondary structural changes. Fluorescence spectra and UV–Vis spectra show that for the three kinds of OH-PBDEs for TTR, the main methods of quenching are static quenching and non-radiative energy transfer. According to thermodynamic analysis, ΔG < 0, ΔH > 0, and ΔS > 0 combine to show that the hydrophobic interaction is the main driving force of the combination. From the molecular docking analysis, it was found that 4′-hydroxy-2,2′,4,5′- tetrabromodiphenyl ether (4′-OH-BDE49) and 4 hydroxy-2,2′,3,4′,5,6,6′- heptabromodiphenyl ether (4-OH-BDE188) had a cationic–π interaction with TTR, whereas 4 hydroxy-2,2′,3,4,5,5′,6- heptabromodiphenyl ether (4-OH-BDE187) was bonded to TTR by hydrogen bonds to form stable complexes. In this paper, we highlight the consistency of spectroscopic experiments and computer simulations so as to provide a reliable analytical method for the toxicological properties of small molecule contaminants.     

Effect of intrinsic and extrinsic factors on the simulated D‐band length of type I collagen

A signature feature of collagen is its axial periodicity visible in TEM as alternating dark and light bands. In mature, type I collagen, this repeating unit, D, is 67 nm long. This periodicity reflects an underlying packing of constituent triple‐helix polypeptide monomers wherein the dark bands represent gaps between axially adjacent monomers. This organization is visible distinctly in the microfibrillar model of collagen obtained from fiber diffraction. However, to date, no atomistic simulations of this diffraction model under zero‐stress conditions have reported a preservation of this structural feature. Such a demonstration is important as it provides the baseline to infer response functions of physiological stimuli. In contrast, simulations predict a considerable shrinkage of the D‐band (11–19%). Here we evaluate systemically the effect of several factors on D‐band shrinkage. Using force fields employed in previous studies we find that irrespective of the temperature/pressure coupling algorithms, assumed salt concentration or hydration level, and whether or not the monomers are cross‐linked, the D‐band shrinks considerably. This shrinkage is associated with the bending and widening of individual monomers, but employing a force field whose backbone dihedral energy landscape matches more closely with our computed CCSD(T) values produces a small D‐band shrinkage of < 3%. Since this force field also performs better against other experimental data, it appears that the large shrinkage observed in earlier simulations is a force‐field artifact. The residual shrinkage could be due to the absence of certain atomic‐level details, such as glycosylation sites, for which we do not yet have suitable data. Proteins 2015; 83:1800–1812. © 2015 Wiley Periodicals, Inc.     

Proton Countertransport and Coupled Gating in the Sarcoplasmic Reticulum Calcium Pump

The calcium pump of the sarcoplasmic reticulum (SERCA) is an ATP-driven active transporter of Ca²⁺ ions that functions via an “alternating-access” cycle mechanism. In each cycle, SERCA transports two Ca²⁺ ions toward the lumen of the sarcoplasmic reticulum and two to three protons to the cytoplasm. How the latter conformational transition is coupled to cytoplasmic release of protons remains poorly understood. The present computational study shows how the mechanism of proton countertransport is coupled to the alternating access gating process in SERCA. Molecular dynamics simulation trajectories are generated starting from a series of configurations taken along the E2 to E1 transition pathway determined by the string method with swarms-of-trajectories. Simulations of different protonation configurations at the binding sites reveal how deprotonation events affect the opening of the cytoplasmic gate. The results show that there is a strong coupling between the chronological order of deprotonation, the entry of water molecules into the TM region, and the opening of the cytoplasmic gate. Deprotonation of E309 and E771 is sequential with E309 being the first to lose the proton. The deprotonation promotes the opening of the cytoplasmic gate but leads to a productive gating transition only if it occurs after the transmembrane domain has reached an intermediate conformation. Deprotonation of E309 and E771 is unproductive when it occurs too early because it causes the re-opening of the luminal gate.     

Molecular dynamics simulation of film size and vacancy defect effects on the thermal conductivities of argon thin films

It is well known that there is a size effect for the thermal conductivity of thin films and that vacancy defects in film reduce the film's thermal conduction. In this paper, the film size and vacancy defect effects on the thermal conductivities of argon thin films were studied by molecular dynamics simulations. The results show the existence of phonon boundary scattering. The results also confirm that the theoretical model based on the Boltzmann equation can accurately model the thermal conduction of thin argon films. Both the theoretical and MD results illustrate that, although, both the defect and the thickness of the thin film deduce the thermal conductivity, their physical mechanisms differ.     

Nucleotide modifications and tRNA anticodon-mRNA codon interactions on the ribosome

We have carried out molecular dynamics simulations of the tRNA anticodon and mRNA codon, inside the ribosome, to study the effect of the common tRNA modifications cmo(5)U34 and m(6)A37. In tRNA(Val), these modifications allow all four nucleotides to be successfully read at the wobble position in a codon. Previous data suggest that entropic effects are mainly responsible for the extended reading capabilities, but detailed mechanisms have remained unknown. We have performed a wide range of simulations to elucidate the details of these mechanisms at the atomic level and quantify their effects: extensive free energy perturbation coupled with umbrella sampling, entropy calculations of tRNA (free and bound to the ribosome), and thorough structural analysis of the ribosomal decoding center. No prestructuring effect on the tRNA anticodon stem-loop from the two modifications could be observed, but we identified two mechanisms that may contribute to the expanded decoding capability by the modifications: The further reach of the cmo(5)U34 allows an alternative outer conformation to be formed for the noncognate base pairs, and the modification results in increased contacts between tRNA, mRNA, and the ribosome.     

山奈酚和槲皮素对大鼠体内Neu5Gc合成的影响及分子机制探究

非人源性唾液酸N-羟乙酰神经氨酸(N-glycolylneuraminic acid, Neu5Gc)是红肉中潜在的致癌性因子。唾液酸转移酶是涉及转运其前体物质的关键酶之一。大鼠用不同浓度的山奈酚(kaempferol,KA)和槲皮素(quercetin,Qu)灌胃,并模拟宰前对Neu5Gc合成影响的结果表明,不同浓度的KA、Qu对大鼠体内肌肉组织、肝和肾中Neu5Gc的含量均有一定的影响,最大抑制率分别为36.44%±0.11%和33.37%±0.08%。为探究其抑制机制,采用唾液酸转移酶(sialyltransferase,ST)与KA、Qu进行分子对接和分子动力学模拟分析。分子对接复合物二维相互作用图表明:KA和Qu能有效占据ST特异性抑制剂5′-胞苷单磷酸的活性位点氨基酸残基;对最佳复合物进行200 ns的动力学模拟结果显示,ST和KA产生稳定氢键作用的氨基酸残基主要是His301和Gly298,与Qu产生氢键作用的氨基酸残基主要是Gly293、Glu324、Thr272和Ser276。结合自由能分析表明,范德华力、静电吸引作用对抑制过程具有重要作用。本研究为合成和筛选高效ST抑制剂提供了一定实验依据。     

Study on the inhibitory mechanism and binding mode of the hydroxycoumarin compound NSC158393 to HIV‐1 integrase by molecular modeling

Human immunodeficiency virus type 1 integrase (IN) is an essential enzyme in the life cycle of this virus and also an important target for the study of anti‐HIV drugs. In this work, the binding modes of the wild type IN core domain and the two mutants, that is, W132G and C130S, with the 4‐hydroxycoumarin compound NSC158393 were evaluated by using the “relaxed complex” molecular docking approach combined with molecular dynamics (MD) simulations. Based on the monomer MD simulations, both of the two substitutions affect not only the stability of the 128–136 peptides, but also the flexibility of the functional 140s loop. In principle, NSC158393 binds the 128–136 peptides of IN; however, the specific binding modes for the three systems are various. According to the binding mode of NSC158393 with WT, NSC158393 can effectively interfere with the stability of the IN dimer by causing a steric hindrance around the monomer interface. Additionally, through the comparative analysis of the MD trajectories of the wild type IN and the IN‐NSC158393 complex, we found that NSC15893 may also exert its inhibitory function by diminishing the mobility of the function loop of IN. Three key binding residues, that is, W131, K136, and G134, were discovered by energy decomposition calculated with the Molecular Mechanics Generalized Born Surface Area method. Characterized by the largest binding affinity, W131 is likely to be indispensable for the ligand binding. All the above results are consistent with experiment data, providing us some helpful information for understanding the mechanism of the coumarin‐based inhibitors. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 700–709, 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     

Simulation of the different biological activities of diethylstilbestrol (DES) on estrogen receptor alpha and estrogen-related receptor gamma

We describe the application of the molecular dynamics (MD) and molecular mechanics-generalized Born/surface area (MM-GB/SA) approaches to the simulation of the different biological activity of diethylstilbestrol (DES) on two highly homologous nuclear receptors-estrogen receptor alpha (ER-alpha) and estrogen-related receptor gamma (ERR-gamma). DES exerts an agonistic effect against ER-alpha and an antagonistic effect against ERR-gamma. Using the x-ray crystal structures of ER-alpha in the canonical agonist bound form (PDB code: 3ERD) and antagonist bound form (PDB code: 3ERT), ERR-gamma homology models have been constructed for the receptor in two different conformations. MM-GB/SA binding free energy calculations of DES in the ER-alpha and ERR-gamma structures suggest that DES exhibits a greater free energy of binding in the agonist bound conformation of ER-alpha, while the antagonist bound conformation is preferred for ERR-gamma. Further dissection of the free energy contributions coupled with calculation of the ligand binding pocket volume suggests that the van der Waals interactions for DES within the smaller binding pocket volume of ERR-gamma are less favorable and this is the main factor for DES antagonism in ERR-gamma. This approach has potential general applicability to the prediction of the biological activity of nuclear receptor ligands.     

Effect of chirality,length and diameter of carbon nanotubes on the adsorption of 20 amino acids: a molecular dynamics simulation study

We carried out molecular dynamics simulations to study the adsorption of all the 20 amino acids (AAs; aromatic, polar and non-polar) on the surface of chiral, zigzag and armchair single-walled carbon nanotubes. The adsorption was occurring in all systems. In the aromatic AAs, the π–π stacking and the semi-hydrogen bond formation cause a strong interaction with the carbon nanotubes (CNTs). We also investigated the chirality, length and diameter dependencies on adsorption energies. We found that all AAs have more tendency to adsorption on the chiral and zigzag CNTs over the armchair. The results show that increasing both the diameter and the length causes the enhancement of the adsorption energy. But, the effect of the length is more than of the diameter. For example, the adsorption energy of Trp on the surface of CNT (4,1), with 2 nm length, is 20.4 kcal/mol. When the length of CNT becomes twice, the adsorption energy increases by 24 ± 0.3%. But by doubling the diameter, the adsorption energy increased only by 9.8 ± 0.25%.