Molecular Dynamics Simulations and Steered Molecular Dynamics Simulations of Glabridin Bound to Wild Type and V30A Mutant Transthyretin: Ligand-linked Perturbation of Tertiary Conformation

Transthyretin(TTR), as a tetrameric protein, functions as a neuroprotector. The native TTR homotetramer dissociates into dimers and monomers. Dimers and monomers self-assemble into amyloid fibrils, and this process can lead to some diseases. Native TTR homotetramer is a widely accepted model for TTR amyloid formation. In this study, simulations using molecular dynamics(MD) and steered MD(SMD) were performed to explore the mechanisms for glabridin(Glab), a specific inhibitor for TTR binding, for V30A mutant and wild-type(WT) TTR. MD simulation results indicate that, compared with Glab binding to WT and V30A mutant, the WT TTR could lead to the collapse of β-strands from Ser52 to His56 at chain A. This phenomenon facilitated the easy dissociation of chains A and C. Calculations of the binding free energy between the two chains showed that the V30A-Glab TTR complex displayed a lower binding energy than other systems(WT TTR and WT-Glab TTR). Then, SMD simulation was performed to explore the unbinding pathway for Glab through the WT and V30A mutant TTR. The results show that Lys15(chain A) produced a hydrogen bond with Glab at the force peak via the WT TTR tunnel. Meanwhile, in the V30A TTR mutant, the hydrogen bond between Lys15(chain A) and Glab was broken at the force peak. This condition was beneficial for Glab to be taken off from the protein. Our theoretical results will be useful in designing a new specific inhibitor of TTR protein to control the TTR homotetramer dissociation.     

Conjugated β-Cyclodextrin Enhances the Affinity of Folic Acid towards FRα: Molecular Dynamics Study

Drug targeting is a progressive area of research with folate receptor alpha (FRα) receiving significant attention as a biological marker in cancer drug delivery. The binding affinity of folic acid (FA) to the FRα active site provides a basis for recognition of FRα. In this study, FA was conjugated to beta-cyclodextrin (βCD) and subjected to in silico analysis (molecular docking and molecular dynamics (MD) simulation (100 ns)) to investigate the affinity and stability for the conjugated system compared to unconjugated and apo systems (ligand free). Docking studies revealed that the conjugated FA bound into the active site of FRα with a docking score (free binding energy < −15 kcal/mol), with a similar binding pose to that of unconjugated FA. Subsequent analyses from molecular dynamics (MD) simulations, root mean square deviation (RMSD), root mean square fluctuation (RMSF), and radius of gyration (Rg) demonstrated that FA and FA–βCDs created more dynamically stable systems with FRα than the apo-FRα system. All systems reached equilibrium with stable RMSD values ranging from 1.9–2.4 Å and the average residual fluctuation values of the FRα backbone atoms for all residues (except for terminal residues ARG8, THR9, THR214, and LEU215) were less than 2.1 Å with a consistent Rg value of around 16.8 Å throughout the MD simulation time (0–100 ns). The conjugation with βCD improved the stability and decreased the mobility of all the residues (except residues 149–151) compared to FA–FRα and apo-FRα systems. Further analysis of H-bonds, binding free energy (MM-PBSA), and per residue decomposition energy revealed that besides APS81, residues HIS20, TRP102, HIS135, TRP138, TRP140, and TRP171 were shown to have more favourable energy contributions in the holo systems than in the apo-FRα system, and these residues might have a direct role in increasing the stability of holo systems.     

分子动力学模拟研究质子化态在HIV-1 Protease-Indinavir复合物中的作用

I 型人体免疫缺陷病毒(HIV-1)蛋白酶中Asp25/Asp25'的质子化对于理论研究HIV-1 蛋白酶和抑制剂的作用机制
以及氨基酸变异对抗药性的影响有重要意义. 分别对Protease-Indinavir (PR-IDV)复合物的六种可能的质子化态进行了
5 ns 的分子动力学模拟, 分析了不同状态对动力学特征和结构的影响, 用molecular mechanics/Possion-Boltzman surface
area (MM-PBSA)方法计算了PR 和IDV 在各种状态下的结合自由能. 计算结果说明A 链Asp25 的OD2 的质子化是最
为可能的状态. 对PR-IDV 复合物中起到媒介作用的水分子与PR-IDV 复合物形成的氢键进行了分析, 分析结果说明不
同的质子化态对水分子在PR-IDV 复合物中所起的媒介作用没有影响, 这一结果与我们先前对PR-BEA369 复合物的研
究不同. 我们的研究结果为更高效的PR 抑制剂的设计以及PR 氨基酸变异对药物抗药性的研究提供了理论上的指导.     

乙酰胆碱酯酶AChE与1,7-二氮杂咔唑抑制剂的作用机理的分子动力学模拟

通过分子对接、分子动力学(MD)模拟以及成键自由能分析方法,从原子水平上模拟研究了3种1,7-二氮杂咔唑衍生物(分别记为M1、M2和M3)与ACh E的结合模式及相互作用机理,分析和讨论了研究体系的静电相互作用和范德华相互作用(vd W)。用MM-PBSA方法计算的3种抑制剂与ACh E之间的结合自由能与抑制剂的实验生物活性数据(IC50值)相对应。分析结果表明,残基S286与抑制剂之间形成的氢键作用有利于抑制剂与ACh E之间的结合。范德华相互作用,尤其是抑制剂与关键残基W279和Y334的作用,对抑制剂与ACh E之间的结合自由能有较大的贡献,在区分抑制剂M1(或M2)和M3的生物活性上发挥着重要的作用。     

p53-MDM2相互作用的分子力学和动力学研究

p53-MDM2相互作用的抑制已经成为治疗癌症的新方法. 本文将分子动力学模拟和MM-PBSA(molecular mechanics/passion-Boltzman surface area)方法结合起来研究MDM2-p53相互作用机制. 结果证明范德华相互作用驱动了MDM2与p53的结合. 基于残基-残基相互作用的计算不仅证明p53的三个残基Phe19′, Trp23′和Leu26′与MDM2有较强的相互作用,而且还发现另外两个残基Leu22′和Pro27′也与MDM2有较强的相互作用,这为抗癌药物的设计提供了新靶标. 同时也证明CH-CH,CH-π和π-π相互作用驱动了p53在MDM2疏水性裂缝中的结合.     

抑制剂8CA与脂肪细胞脂肪酸结合蛋白(A-FABP)结合模式的分子动力学研究

摘要：脂肪细胞脂肪酸结合蛋白A-FABP(Adipocyte fatty-acid binding protein)是治疗脂质调节生物过程相关疾病的重要靶标. 分子动力学模拟和MM-PBSA方法被采用研究抑制剂8CA与A-FABP结合模式. 研究结果表明静电相互作用和范德华作用驱动了抑制剂8CA与A-FABP的结合。基于残基的能量分解表明抑制剂8CA与R126间的极性相互作用为抑制剂与A-FABP的结合提供了重要贡献. 该残基与8CA的相互作用较好地稳定了抑制剂与A-FABP复合物的稳定性. 我们期望这个研究能为治疗炎症、动脉硬化和代谢病药物设计提供一定的理论指导。     

5,10-Methylene-5,6,7,8-tetrahydrofolate conformational transitions upon binding to thymidylate synthase: molecular mechanics and continuum solvent studies

Summary We applied the molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) approach to evaluate relative stability of the extended (flat) and C-shaped (bent) solution conformational forms of the 5,10-methylene-5,6,7,8-tetrahydrofolate (mTHF) molecule in aqueous solution. Calculations indicated that both forms have similar free energies in aqueous solution but detailed energy components are different. The bent solution form has lower intramolecular electrostatic and van der Waals interaction energies. The flat form has more favorable solvation free energy and lower contribution from the bond, angle and torsion angle molecular mechanical internal energies. We exploit these results and combine them with known crystallographic data to provide a model for the progressive binding of the mTHF molecule, a natural cofactor of thymidylate synthase (TS), to the complex forming in the TS-catalyzed reaction. We propose that at the time of initial weak binding in the open enzyme the cofactor molecule remains in a close balance between the flat and bent solution conformations, with neither form clearly favored. Later, thymidylate synthase undergoes conformational change leading to the closure of the active site and the mTHF molecule is withdrawn from the solvent. That effect shifts the thermodynamic equilibrium of the mTHF molecule toward the bent solution form. At the same time, burying the cofactor molecule in the closed active site produces numerous contacts between mTHF and protein that render change in the shape of the mTHF molecule. As a result, the bent solution conformer is converted to more strained L-shaped bent enzyme conformer of the mTHF molecule. The strain in the bent enzyme conformation allows for the tight binding of the cofactor molecule to the productive ternary complex that forms in the closed active site, and facilitates the protonation of the imidazolidine N10 atom, which promotes further reaction. *To whom correspondence should be addressed. Fax: +4822-822-5342; E-mail: a.jarmula@nencki.gov.pl     

分子动力学研究抑制剂pDI6W与MDM2的相互作用机制

p53-MDM2相互作用已经成为治疗癌症药物设计的重要靶标。本文采用分子动力学模拟和MM-PBSA(molecular mechanics-Possion-Boltzmann surface area)方法计算了肽类抑制剂pDI6W与MDM2的结合自由能,结果证明范德华作用是pDI6W与MDM2结合的主要力量。相关矩阵的计算结果表明pDI6W的结合主要诱导了MDM2内部的反相关运动。基于成对残基的自由能分解计算不仅证明pDI6W的5个残基Phe19′, Trp22′, Trp23′, Leu26′和Thr27′能够与MDM2产生较强的相互作用,而且确认了CH-π, CH-CH和π-π相互作用驱动了pDI6W在MDM2疏水裂缝中的结合。这为抗癌药物的设计提供了理论上的指导。     

Dimerization energetics of the G-protein coupled bile acid receptor TGR5 from all-atom simulations

We describe the first extensive energetic evaluation of GPCR dimerization on the atomistic level by means of potential of mean force (PMF) computations and implicit solvent/implicit membrane end-point free energy calculations (MM-PBSA approach). Free energies of association computed from the PMFs show that the formation of both the 1/8 and 4/5 interface is energetically favorable for TGR5, the first GPCR known to be activated by hydrophobic bile acids and neurosteroids. Furthermore, formation of the 1/8 interface is favored over that of the 4/5 interface. Both results are in line with our previous FRET experiments in live cells. Differences in lipid-protein interactions are identified to contribute to the observed differences in free energies of association. A per-residue decomposition of the MM-PBSA effective binding energy reveals hot spot residues specific for both interfaces that form clusters. This knowledge may be used to guide the design of dimerization inhibitors or perform mutational studies to explore physiological consequences of distorted TGR5 association. Finally, we characterized the role of Y111, located in the conserved (D/E)RY motif, as a facilitator of TGR5 interactions. The types of computations performed here should be transferable to other transmembrane proteins that form dimers or higher oligomers as long as good structural models of the dimeric or oligomeric states are available. Such computations may help to overcome current restrictions due to an imperfect energetic representation of protein association at the coarse-grained level. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.     

Computational Investigation on the Ethylene-induced Esterase from Citrus Sinensis

In order to understand the interaction between ethylene-induced esterase(EIE, EC 3.1.1.) from Citrus sinensis and α-naphthyl acetate, a 3D model of EIE was generated based on the crystal structure of the tobacco salicylic acid binding protein 2(SABP2). With the aid of the molecular mechanics and molecular dynamics methods, the final refined model was obtained and its reliability was further assessed. In this study, the docking results show that the main-chain amide of residue His85 and residue Val18 can form hydrogen bonds to the carbonyl oxygen group of α-naphthyl acetate. MM-PBSA method was applied to calculating the binding free energy between EIE mutants and α-naphthyl acetate. Our calculated binding free energy of each of the two mutant complexes is increased compared with that of the one of the wild type, which is unfavorable to the reaction. It is well consistent with the experimental data. The above results clearly indicate that His85 and Val18 in EIE function as the oxyanion role and take part in the catalytic reaction. The new structural insights obtained from this computational study are expected to stimulate further biochemical studies on the structures and mechanisms of EIE and other members of the plant α/β hydrolases.