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

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

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.     

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.     

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.     

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%.     

The mechanism for the complexation and dissociation between siRNA and PMAL: a molecular dynamics simulation study based on a coarse-grained model

Abstract

The capacity of silencing genes makes small interfering RNA (siRNA) becomes potential candidates for curing many fatal diseases. Due to the low stability and delivery efficiency of siRNA, the design of amphiphilic carrier for siRNA delivery is vital for the practical gene therapy. In the present work, we explored how the complexation and dissociation of siRNA with poly (maleic anhydride-alt-1-decene) substituted with 3-(dimethylamino) propylamine (PMAL), which is a recent synthesised amphiphilic polymer and can be used in delivery of siRNAs and proteins, using traditional molecular dynamics simulations, together with steered molecular dynamics simulations. It was shown that the complexation of siRNA with PMALs can spontaneously occur, no matter what unit number of PMAL is. PMALs of different unit numbers form micelle-like structures and interact with siRNA surface. With the increase of unit number, PMAL becomes more flexible and interacts with siRNA from attachment to entanglement. The dissociation of PMAL from siRNA is an energy-consuming process. The free energy difference increases with the unit number of PMAL. The free energy for dissociation involves both the stretch of PMAL and the separation of PMAL from siRNA. Therefore, an optimal unit number of PMAL is critical for the delivery efficiency of siRNA when PMAL is used as carrier. In present work, when the radius of gyration of PMAL approaches to that of siRNA, PMAL gives a favoured both complexation and dissociation between siRNA and PMAL. Finally, we propose the mechanism of complexation and dissociation of PMAL with siRNA. The above simulation established a molecular insight of the interaction between siRNA and PMAL and was helpful for the design and applications of new PMAL-based polymers as siRNA delivery carriers.     

双酚A(BPA)与黑腹果蝇雌激素相关受体(dERR)的结合模式及对其基因表达的影响

【目的】雌激素相关受体(estrogen-related receptors, ERRs)属于核受体(nuclear receptor, NR)超家族,在昆虫新陈代谢和能量转换调节中发挥重要作用。双酚A(bisphenol A, BPA)与昆虫生殖和神经系统疾病有关。本研究旨在探究BPA影响黑腹果蝇Drosophila melanogaster ERR(dERR)的作用机理。【方法】通过AutoDock Vina模拟小分子BPA与Modeller 9.25构建的dERR蛋白的分子对接,使用Gromacs 5.1.9进行分子动力学模拟,并结合自由能的计算探究BPA与dERR的结合模式。利用qRT-PCR方法检测3种浓度(0.1, 1和10 μg/L) BPA处理6, 12和24 h对黑腹果蝇成虫和2龄幼虫中dERR基因转录水平的影响。【结果】通过分子对接和极性溶剂化能发现,Phe370及Leu334等侧链氨基酸是BPA与dERR结合的关键氨基酸。qRT-PCR分析显示,不同浓度的BPA处理6和12 h时黑腹果蝇成虫和2龄幼虫中dERR基因转录水平均发生显著变化,而0.1 μg/L BPA处理24 h时的几乎接近对照。【结论】BPA能够影响黑腹果蝇体内dERR的表达,作用机理可能跟BPA与dERR的潜在特异性结合有关。     

Insight into the key active sites of afChiA1 based on molecular dynamics simulations and free energy calculations

In this study, the binding of the enzyme chitinase A1 (afChiA1) from the plant-type Aspergillus fumigatus with four potent inhibitors, allosamidin (ASM), acetazolamide (AZM), 8-chloro-theophylline (CTP) and kinetin (KIT) is investigated by molecular docking, molecular dynamics simulation and binding free energy calculation. The results reveal that the electrostatic interactions play an important role in the stabilisation of the binding of afChiA1 with inhibitors. Based on the binding energy of afChiA1-ligands, the key residues (Gln37 and Trp312) in the active binding pocket of the complex systems are confirmed by molecular mechanics/Poisson–Boltzmann surface area method, and the active inhibitors, ASM and AZM, both could form strong interaction with Gln37 and Trp312, and the non-active ligands, CTP and KIT, could not interact with these two residues, which is consistent with the result of experimental report. Then, it is identified that Gln37 and Trp312 should be one of the important active site residues of afChiA1.     

A coarse grained molecular dynamics study on the structure and stability of small-sized liposomes

The dependence of geometric structure and thermal stability of liposomes on their component phospholipid molecules and distribution of molecules in the inner and the outer layers of the liposome is investigated by conducting molecular simulations in explicit water for the eight types of liposomes constructed from different phospholipids. Using molecular mechanics structure-relaxation based on the coarse grained (CG) model, stable structures of the solvated liposomes are obtained. In addition, the molecular dynamics (MD) simulations based on the CG model are carried out at 310 and 360 K for elucidating the change in structure of the solvated liposomes. The MD simulations reveal that liposomes having the same number of lipids (SNL) in both the inner and the outer layers keep their spherical structures even at 360 K. In particular, the SNLs composed of palmitoyloleoyl-phosphatidyl-ethanolamine1 or dimyristoylglycero-phosphatidyl-choline lipid exhibit a compact spherical shape. In contrast, liposomes having the same density of lipids in the inner and the outer layers cannot keep their spherical shapes at 360 K. The obtained results contribute toward developing novel liposomes with enhanced thermal stability.     

ATP Synthases in the Year 2000: Defining the Different Levels of Mechanism and Getting a Grip on Each

ATP synthases are unusually complex molecules, which fractionate most readily into two major units, one a water soluble unit called F₁ and the other a detergent soluble unit called F₀. In almost all known species the F₁ unit consists of 5 subunit types in the stoichiometric ratio ₃₃ while the F₀ unit contains 3 subunit types (a, b, and c) in E. coli, and at least 10 subunit types (a, b, c, and others) in higher animals. It is now believed by many investigators that during the synthesis of ATP, protons derived from an electrochemical gradient generated by an electron transport chain are directed through the F₀ unit in such a way as to drive the rotation of the single subunit, which extends from an oligomeric ring of at least 10 c subunits in F₀ through the center of F₁. It is further believed by many that the rotating subunit, by interacting sequentially with the 3 pairs of F₁ (360° cycle) in the presence of ADP, P_i, and Mg⁺⁺, brings about via power strokes conformational/binding changes in these subunits that promote the synthesis of ATP and its release on each pair. In support of these views, studies in several laboratories either suggest or demonstrate that F₀ consists in part of a proton gradient driven motor while F₁ consists of an ATP hydrolysis driven motor, and that the subunit does rotate during F₁ function. Therefore, current implications are that during ATP synthesis the former motor drives the latter in reverse via the subunit. This would suggest that the process of understanding the mechanism of ATP synthases can be subdivided into three major levels, which include elucidating those chemical and/or biophysical events involved in (1) inducing rotation of the subunit, (2) coupling rotation of this subunit to conformational/binding changes in each of the 3 pairs, and (3) forming ATP and water (from ADP, P_i, and Mg⁺⁺) and then releasing these products from each of the 3 catalytic sites. Significantly, it is at the final level of mechanism where the bond breaking/making events of ATP synthesis occur in the transition state, with the former two levels of mechanism setting the stage for this critical payoff event. Nevertheless, in order to get a better grip in this new century on how ATP synthases make ATP and then release it, we must take on the difficult challenge of elucidating each of the three levels of mechanism.     

矿区可持续发展系统动力学模拟与调控

矿区是由资源、环境、经济和社会等子系统构成的复杂系统。矿区的可持续发展依赖于矿各子系统的合理结构和发展模式以及人们对系统的有效调控,将可持续发展思想 模拟方法相结合,在系统地分析矿区ＲＥＥＳ系统结构的基础上,构建了矿区ＲＥＥＳ系统动力学模型,并以铁法矿为例论述了矿区ＲＥＥＳ系统的模拟和调控等有关问题。     

The effect of structural parameters and positive charge distance on the interaction free energy of antimicrobial peptides with membrane surface

Many attempts have been made to find hints explaining the relationship between physicochemical and structural properties of antimicrobial peptides (AMPs) which are relevant to their antimicrobial activities. We here found that there is a difference in the percentages of hydrophobic, hydrophilic, and charged residues between AMPs killing both bacteria and fungi (Group A) and AMPs that only kill bacteria (Group B). The percentage of charged residues in Group A AMPs is highly elevated, while in Group B the percentage of hydrophobic residues is increased. This result suggests a sequence-based mechanism of selectivity for AMPs. Moreover, we examined how the distance between basic residues affects the interaction free energy of AMPs with the membrane surface, since most of the known AMPs act by membrane perturbation. We measured the average distance between basic residues throughout the 3D structure of AMPs by defining D_pr parameter and calculated the interaction free energy for 10 AMPs that interacted with the DPPC membrane using molecular dynamics simulation. We found that the changes of the interaction free energy correlates with the change of D_pr by a linear regression coefficient of r^2?=?.47 and a cubic regression coefficient of r^2?=?.70.     

Introducing a novel method based on the imperialistic competitive algorithm to determine fluorine intermolecular potential from ab initio calculations and calculation of some properties via MD simulations

In this work, the possibility of obtaining an accurate site-site potential model suitable for use in molecular dynamics (MD) simulations of fluorine from ab initio calculations has been explored. The exploration was made on ab initio calculations. To reduce the ab initio pair potentials into a site-site potential, a higher significance was assigned to the configuration which is more stable. For this purpose, the imperialistic competitive algorithm (ICA) was implemented as a powerful optimisation tool. The calculated second virial coefficients were compared to the experimental values to test the quality of the presented intermolecular potential. The relative error for the calculated second virial coefficient ranged from 0.1 to 5.6%. MD simulations were used to evaluate the ability of the proposed intermolecular potential function. The relative error for the MD simulations ranged from 0.5 to 5.2%. The results are in good agreement with experimental data.     

Effects of substrate binding on the dynamics of RNase A: Molecular dynamics simulations of UpA bound and native RNase A

RNase A has been extensively used as a model protein in several biophysical and biochemical studies. Using the available structural and biochemical results, RNase A-UpA interaction has been computationally modeled at an atomic level. In this study, the molecular dynamics (MD) simulations of native and UpA bound RNase A have been carried out. The gross dynamical behavior and atomic fluctuations of the free and UpA bound RNase A have been characterized. Principal component analysis is carried out to identify the important modes of collective motion and to analyze the changes brought out in these modes of RNase A upon UpA binding. The hydrogen bonds are monitored to study the atomic details of RNase A-UpA interactions and RNase A-water interactions. Based on these analysis, the stability of the free and UpA bound RNase A are discussed. © 1997 John Wiley & Sons, Inc. Biopoly 42: 505–520, 1997     

Quantum simulation for muoniated and deuterated methyl radicals in implicit water solvent: combined ab initio path integral molecular dynamics and the polarizable continuum model simulation study

We have combined ab initio path integral molecular dynamics (PIMD) simulation and the polarizable continuum model (PCM) method to efficiently incorporate solvent effects into nuclear quantum fluctuation of molecular systems. Our combined ab initio PIMD–PCM simulation was applied to muoniated and deuterated methyl radical immersed in implicit water solvent to gain information on solvent and isotope effects from one simulation run. We found that solvent effects lead to the bond elongation and a decrease in the magnitude of isotropic hyperfine coupling constants. These are consistent with the trends in conventional static calculations and experiments. In addition, the performance of cavity models (universal force field, united atom specified for Kohn–Sham and these hybrid models) and the conservation of the PIMD–PCM Hamiltonian were accessed. We confirmed that solvent effects on nuclear quantum fluctuation are efficiently computed using our combined simulation of quantum solute in implicit solvent.     

运用隔室理论研究中药材吸水膨胀动力学总室线性乳突数学模型及参数分析

目的:运用隔室模型研究中药材吸水膨胀动力学总室线性乳突数学模型并对参数进行分析;方法:运用隔室理论及动力学方法建立中药材吸水膨胀动力学数学模型,用拉氏变换法求解,研究总室与单独隔室的数学模型的关系,并对总室模型求算的各参数进行分析.结果:中药材吸水膨胀总室,与各隔室一样遵循一级线性乳突模型,V_T为e的指数函数之和形式,总室膨胀为各分室膨胀之和.V_T~∞、V_i~∞可通过非线性曲线参数估计算得,各参数(转运常数)由其与α_i的关系方程组算得.结论:中药材吸水膨胀动力学表现为一线性乳突模型.总室与单独室一样可求得各动力学参数,该模型在一定条件下可直接转变为药物动力学数学模型.     

Mathematical simulation of the interactions among cyanobacteria, purple sulfur bacteria and chemotrophic sulfur bacteria in microbial mat communities

Abstract: A deterministic one-dimensional reaction diffusion model was constructed to simulate benthic stratification patterns and population dynamics of cyanobacteria, purple and colorless sulfur bacteria as found in marine microbial mats. The model involves the major biogeochemical processes of the sulfur cycle and includes growth metabolism and their kinetic parameters as described from laboratory experimentation. Hence, the metabolic production and consumption processes are coupled to population growth. The model is used to calculate benthic oxygen, sulfide and light profiles and to infer spatial relationships and interactions among the different populations. Furthermore, the model is used to explore the effect of different abiotic and biotic environmental parameters on the community structure. A strikingly clear pattern emerged of the interaction between purple and colorless sulfur bacteria: either colorless sulfur bacteria dominate or a coexistence is found of colorless and purple sulfur bacteria. The model predicts that purple sulfur bacteria only proliferate when the studied environmental parameters surpass well-defined threshold levels. However, once the appropriate conditions do occur, the purple sulfur bacteria are extremely successful as their biomass outweighs that of colorless sulfur bacteria by a factor of up to 17. The typical stratification pattern predicted closely resembles the often described bilayer communities which comprise a layer of purple sulfur bacteria below a cyanobacterial top-layer; colorless sulfur bacteria are predicted to sandwich in between both layers. The profiles of oxygen and sulfide shift on a diel basis similarly as observed in real systems.     

EPR spectra and molecular dynamics agree that the nucleotide pocket of myosin V is closed and that it opens on binding actin

We have used EPR spectroscopy and computational modeling of nucleotide-analog spin probes to investigate conformational changes at the nucleotide site of myosin V. We find that, in the absence of actin, the mobility of a spin-labeled diphosphate analog [spin-labeled ADP (SLADP)] bound at the active site is strongly hindered, suggesting a closed nucleotide pocket. The mobility of the analog increases when the MV·SLADP complex (MV = myosin V) binds to actin, implying an opening of the active site in the A·MV·SLADP complex (A = actin). The probe mobilities are similar to those seen with myosin II, despite the fact that myosin V has dramatically altered kinetics. Molecular dynamics (MD) simulation was used to understand the EPR spectra in terms of the X-ray database. The X-ray structure of MV·ADP·BeFx shows a closed nucleotide site and has been proposed to be the detached state. The MV·ADP structure shows an open nucleotide site and has been proposed to be the A·MV·ADP state at the end of the working powerstroke. MD simulation of SLADP docked in the closed conformation gave a probe mobility comparable to that seen in the EPR spectrum of the MV·SLADP complex. The simulation of the open conformation gave a probe mobility that was 35-40° greater than that observed experimentally for the A·MV·SLADP state. Thus, EPR, X-ray diffraction, and computational analysis support the closed conformation as a myosin V state that is detached from actin. The MD results indicate that the MV·ADP crystal structure, which may correspond to the strained actin-bound post-powerstroke conformation resulting from head-head interaction in the dimeric processive motor, is superopened.