Hydrophobic Interactions Are a Key to MDM2 Inhibition by Polyphenols as Revealed by Molecular Dynamics Simulations and MM/PBSA Free Energy Calculations

p53, a tumor suppressor protein, has been proven to regulate the cell cycle, apoptosis, and DNA repair to prevent malignant transformation. MDM2 regulates activity of p53 and inhibits its binding to DNA. In the present study, we elucidated the MDM2 inhibition potential of polyphenols (Apigenin, Fisetin, Galangin and Luteolin) by MD simulation and MM/PBSA free energy calculations. All polyphenols bind to hydrophobic groove of MDM2 and the binding was found to be stable throughout MD simulation. Luteolin showed the highest negative binding free energy value of -173.80 kJ/mol followed by Fisetin with value of -172.25 kJ/mol. It was found by free energy calculations, that hydrophobic interactions (vdW energy) have major contribution in binding free energy.     

The protein-protein interactions between SMPI and thermolysin studied by molecular dynamics and MM/PBSA calculations

Thermolysin is a zinc-metalloendopeptidase secreted by the gram-positive thermophilic bacterium Bacillus thermoproteolyticus. Thermolysin belongs to the gluzinicin family of enzymes, which is selectively inhibited by Steptomyces metalloproteinase inhibitor (SMPI). Very little is known about the interaction between SMPI and thermolysin. Knowledge about the protein-protein interactions is very important for designing new thermolysin inhibitors with possible industrial or pharmaceutical applications. In the present study, two binding modes between SMPI and thermolysin were studied by 2300 picoseconds (ps) of comparative molecular dynamics (MD) simulations and calculation of the free energy of binding using the molecular mechanics-Poisson-Boltmann surface area (MM/PBSA) method. One of the positions, the 'horizontal arrow head docking' (HAHD) was similar to the previously proposed binding mode by Tate et al. (Tate, S., Ohno, A., Seeram, S. S., Hiraga, K., Oda, K., and Kainosho, M. J. Mol. Biol. 282, 435-446 (1998)). The other position, the 'vertical arrow head docking' (VAHD) was obtained by a manual docking guided by the shape and charge distribution of SMPI and the binding pocket of thermolysin. The calculations showed that SMPI had stronger interactions with thermolysin in the VAHD than in the HAHD complex, and the VAHD complex was considered more realistic than the HAHD complex. SMPI interacted with thermolysin not only at the active site but had auxiliary binding sites contributing to proper interactions. The VAHD complex can be used for designing small molecule inhibitors mimicking the SMPI-thermolysin binding interfaces.     

Revisiting Free Energy Calculations: A Theoretical Connection to MM/PBSA and Direct Calculation of the Association Free Energy

The prediction of absolute ligand-receptor binding affinities is essential in a wide range of biophysical queries, from the study of protein-protein interactions to structure-based drug design. End-point free energy methods, such as the Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) model, have received much attention and widespread application in recent literature. These methods benefit from computational efficiency as only the initial and final states of the system are evaluated, yet there remains a need for strengthening their theoretical foundation. Here a clear connection between statistical thermodynamics and end-point free energy models is presented. The importance of the association free energy, arising from one molecule's loss of translational and rotational freedom from the standard state concentration, is addressed. A novel method for calculating this quantity directly from a molecular dynamics simulation is described. The challenges of accounting for changes in the protein conformation and its fluctuations from separate simulations are discussed. A simple first-order approximation of the configuration integral is presented to lay the groundwork for future efforts. This model has been applied to FKBP12, a small immunophilin that has been widely studied in the drug industry for its potential immunosuppressive and neuroregenerative effects.     

Predicting Binding Free Energy Change Caused by Point Mutations with Knowledge-Modified MM/PBSA Method

A new methodology termed Single Amino Acid Mutation based change in Binding free Energy (SAAMBE) was developed to predict the changes of the binding free energy caused by mutations. The method utilizes 3D structures of the corresponding protein-protein complexes and takes advantage of both approaches: sequence- and structure-based methods. The method has two components: a MM/PBSA-based component, and an additional set of statistical terms delivered from statistical investigation of physico-chemical properties of protein complexes. While the approach is rigid body approach and does not explicitly consider plausible conformational changes caused by the binding, the effect of conformational changes, including changes away from binding interface, on electrostatics are mimicked with amino acid specific dielectric constants. This provides significant improvement of SAAMBE predictions as indicated by better match against experimentally determined binding free energy changes over 1300 mutations in 43 proteins. The final benchmarking resulted in a very good agreement with experimental data (correlation coefficient 0.624) while the algorithm being fast enough to allow for large-scale calculations (the average time is less than a minute per mutation).     

Structural Dynamics of Human Telomeric G-Quadruplex Loops Studied by Molecular Dynamics Simulations

Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.     

Synaptobrevin Transmembrane Domain Dimerization Studied by Multiscale Molecular Dynamics Simulations

Synaptic vesicle fusion requires assembly of the SNARE complex composed of SNAP-25, syntaxin-1, and synaptobrevin-2 (sybII) proteins. The SNARE proteins found in vesicle membranes have previously been shown to dimerize via transmembrane (TM) domain interactions. While syntaxin homodimerization is supposed to promote the transition from hemifusion to complete fusion, the role of synaptobrevin’s TM domain association in the fusion process remains poorly understood. Here, we combined coarse-grained and atomistic simulations to model the homodimerization of the sybII transmembrane domain and of selected TM mutants. The wild-type helix is shown to form a stable, right-handed dimer with the most populated helix-helix interface, including key residues predicted in a previous mutagenesis study. In addition, two alternative binding interfaces were discovered, which are essential to explain the experimentally observed higher-order oligomerization of sybII. In contrast, only one dimerization interface was found for a fusion-inactive poly-Leu mutant. Moreover, the association kinetics found for this mutant is lower as compared to the wild-type. These differences in dimerization between the wild-type and the poly-Leu mutant are suggested to be responsible for the reported differences in fusogenic activity between these peptides. This study provides molecular insight into the role of TM sequence specificity for peptide aggregation in membranes.     

Lipid Membrane Structure and Dynamics Studied by All-Atom Molecular Dynamics Simulations of Hydrated Phospholipid Bilayers

Abstract

The structure and dynamics of phosphatidylcholine bilayers are examined by reviewing the results of several nanoseconds of molecular dynamics simulations on a number of bilayer and monolayer models. The lengths of these simulations, the longest single one of which was 2 nanoseconds, were sufficiently long to effectively sample many of the longer-scale motions governing the behaviour of biomembranes. These simulations reproduce many experimental observables well and provide a degree of resolution currently unavailable experimentally.     

Cryo-Cooling Effect on DHFR Crystal Studied by Replica-Exchange Molecular Dynamics Simulations

Cryo-cooling is routinely performed before x-ray diffraction image collection to reduce the damage to crystals due to ionizing radiation. It has been suggested that although backbone structures are usually very similar between room temperature and cryo-temperature, cryo-cooling may hamper biologically relevant dynamics. In this study, the crystal of Escherichia coli dihydrofolate reductase is studied with replica-exchange molecular dynamics simulation, and the results are compared with the crystal structure determined at cryo-temperature and room temperature with the time-averaged ensemble method. Although temperature dependence of unit cell compaction and root mean-square fluctuation of Cα is found in accord with experiment, it is found that the protein structure at low temperature can be more heterogeneous than the ensemble of structures reported by using the time-averaged ensemble method, encouraging further development of the time-averaged ensemble method and indicating that data should be examined carefully to avoid overinterpretation of one average structure.     

Intrinsic Dynamics of Restriction Endonuclease EcoO109I Studied by Molecular Dynamics Simulations and X-Ray Scattering Data Analysis

EcoO109I is a type II restriction endonuclease that functions as a dimer in solution. Upon DNA binding to the enzyme, the two subunits rotate counterclockwise relative to each other, as the two catalytic domains undergo structural changes to capture the cognate DNA. Using a 150-ns molecular dynamics simulation, we investigated the intrinsic dynamics of the DNA-free enzyme in solution to elucidate the relationship between enzyme dynamics and structural changes. The simulation revealed that the enzyme is considerably flexible, and thus exhibits large fluctuations in the radius of gyration. The small-angle x-ray scattering profile calculated from the simulation, including scattering from explicit hydration water, was in agreement with the experimentally observed profile. Principal component analysis revealed that the major dynamics were represented by the open-close and counterclockwise motions: the former is required for the enzyme to access DNA, whereas the latter corresponds to structural changes upon DNA binding. Furthermore, the intrinsic dynamics in the catalytic domains were consistent with motions capturing the cognate DNA. These results indicate that the structure of EcoO109I is intrinsically flexible in the direction of its functional movement, to facilitate effective structural changes for sequence-specific DNA recognition and processing.     

Exploring the RING-Catalyzed Ubiquitin Transfer Mechanism by MD and QM/MM Calculations

Ubiquitylation is a universal mechanism for controlling cellular functions. A large family of ubiquitin E3 ligases (E3) mediates Ubiquitin (Ub) modification. To facilitate Ub transfer, RING E3 ligases bind both the substrate and ubiquitin E2 conjugating enzyme (E2) linked to Ub via a thioester bond to form a catalytic complex. The mechanism of Ub transfer catalyzed by RING E3 remains elusive. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations, we characterized this catalytic mechanism in detail. The three-dimensional model of dimeric RING E3 ligase RNF4 RING, E2 ligase UbcH5A, Ub and the substrate SUMO2 shows close contact between the substrate and Ub transfer catalytic center. Deprotonation of the substrate lysine by D117 on UbcH5A occurs with almost no energy barrier as calculated by MD and QM/MM calculations. Then, the side chain of the activated lysine gets close to the thioester bond via a conformation change. The Ub transfer pathway begins with a nucleophilic addition that forms an oxyanion intermediate of a 4.23 kcal/mol energy barrier followed by nucleophilic elimination, resulting in a Ub modified substrate by a 5.65 kcal/mol energy barrier. These results provide insight into the mechanism of RING-catalyzed Ub transfer guiding the discovery of Ub system inhibitors.     

Molecular Dynamics of a Schiff Base Tetramacrocycle Studied by NMR and MD Simulations

Schiff base condensation of the pentadentate tetrapodal amine 1 with 2,6-diformyl-4-methylphenol 2 leads in a [3 + 6] condensation to a previously not described macrocycle 3 which shows C₃-symmetry. X-ray analysis shows a truncated cone shape for 3. At T S 100°C, the ¹H-ROESY/EXSY spectrum reveals pairwise exchange of corresponding sites, indicative of inversion of the whole molecule in an umbrella-like fashion. Molecular dynamics simulations support this hypothesis.     

Interactions between Fengycin and Model Bilayers Quantified by Coarse-Grained Molecular Dynamics

Bacteria, particularly of the genus Bacillus, produce a wide variety of antifungal compounds. They act by affecting the lipid bilayers of fungal membranes, causing curvature-induced strain and eventual permeabilization. One class of these, known as fengycins, has been commercialized for treating agricultural infections and shows some promise as a possible antifungal pharmaceutical. Understanding the mechanism by which fengycins damage lipid bilayers could prove useful to the future development of related antifungal treatments. In this work, we present multi-microsecond-long simulations of fengycin interacting with different lipid bilayer systems. We see fengycin aggregation and uncover a clear aggregation pattern that is partially influenced by bilayer composition. We also quantify some local bilayer perturbations caused by fengycin binding, including curvature of the lipid bilayer and local electrostatic-driven reorganization.     

Interactions between Fengycin and Model Bilayers Quantified by Coarse-Grained Molecular Dynamics

Bacteria, particularly of the genus Bacillus, produce a wide variety of antifungal compounds. They act by affecting the lipid bilayers of fungal membranes, causing curvature-induced strain and eventual permeabilization. One class of these, known as fengycins, has been commercialized for treating agricultural infections and shows some promise as a possible antifungal pharmaceutical. Understanding the mechanism by which fengycins damage lipid bilayers could prove useful to the future development of related antifungal treatments. In this work, we present multi-microsecond-long simulations of fengycin interacting with different lipid bilayer systems. We see fengycin aggregation and uncover a clear aggregation pattern that is partially influenced by bilayer composition. We also quantify some local bilayer perturbations caused by fengycin binding, including curvature of the lipid bilayer and local electrostatic-driven reorganization.     

Tautomerism in 8-Nitroguanosine Studied by NMR and Theoretical Calculations

The guanine base in DNA, due to its low oxidation potential, is particularly sensitive to chemical modifications. A large number of guanine lesions have been characterized and studied in some detail due to their relationship with tissue inflammations. Nevertheless, one example of these lesions is the formation of 8-nitro-guanosine, but the NMR data of this compound was only partially interpreted. A comprehensive study of the two possible tautomeric forms, through a detailed characterization of this compound, has implications for its base pairing properties. The target compound was obtained through a synthetic sequence of five steps, where all intermediates were fully characterized using spectral data. The analysis of the two tautomers was then evaluated through NMR spectroscopy and theoretical calculations of the chemical shifts and NH coupling constants, which were also compared with the data from guanosine.     

嗜温冷休克蛋白力致去折叠研究

嗜温冷休克蛋白拥有一个由五个β股形成的反平行β桶结构,目前已被用于蛋白质去折叠的研究。当使用机械力对嗜温冷休克蛋白进行拉伸研究时,发现嗜温冷休克蛋白的去折叠过程具有明显的中间态。在常速和常力两种情况下对嗜温冷休克蛋白进行拉伸分子动力学模拟,发现其在两种情况下具有相同的去折叠次序,即C端β片层首先去折叠,随后N端β片层去折叠;同时这两种模拟都表现出明确的中间态。研究结果表明,嗜温冷休克蛋白抵抗外力作用除了依赖链间氢键外,分子内的静电相互作用也发挥着重要的作用。     

Protein-protein Interactions Visualized by Bimolecular Fluorescence Complementation in Tobacco Protoplasts and Leaves

Many proteins interact transiently with other proteins or are integrated into multi-protein complexes to perform their biological function. Bimolecular fluorescence complementation (BiFC) is an in vivo method to monitor such interactions in plant cells. In the presented protocol the investigated candidate proteins are fused to complementary halves of fluorescent proteins and the respective constructs are introduced into plant cells via agrobacterium-mediated transformation. Subsequently, the proteins are transiently expressed in tobacco leaves and the restored fluorescent signals can be detected with a confocal laser scanning microscope in the intact cells. This allows not only visualization of the interaction itself, but also the subcellular localization of the protein complexes can be determined. For this purpose, marker genes containing a fluorescent tag can be coexpressed along with the BiFC constructs, thus visualizing cellular structures such as the endoplasmic reticulum, mitochondria, the Golgi apparatus or the plasma membrane. The fluorescent signal can be monitored either directly in epidermal leaf cells or in single protoplasts, which can be easily isolated from the transformed tobacco leaves. BiFC is ideally suited to study protein-protein interactions in their natural surroundings within the living cell. However, it has to be considered that the expression has to be driven by strong promoters and that the interaction partners are modified due to fusion of the relatively large fluorescence tags, which might interfere with the interaction mechanism. Nevertheless, BiFC is an excellent complementary approach to other commonly applied methods investigating protein-protein interactions, such as coimmunoprecipitation, in vitro pull-down assays or yeast-two-hybrid experiments.     

壳聚糖与细菌细胞膜相互作用的分子动力学研究

目的 壳聚糖（chitosan,CS）是一种天然的广谱抗菌活性物质。现有研究表明,壳聚糖与细菌细胞膜的相互作用是其发挥抗菌功能的关键。受限于传统实验技术的表征能力,壳聚糖与细菌细胞膜相互作用的具体机制仍有待研究。本文旨在研究壳聚糖与细菌细胞膜相互作用的分子机制。方法 本研究利用全原子分子动力学模拟技术主要探究了完全脱乙酰化的不同聚合度壳聚糖（八聚糖、十二聚糖和十六聚糖）与革兰氏阴性菌外膜（outer membrane,OM）和革兰氏阳性菌质膜（cytoplasmic membrane,CM）相互作用的动态过程。结果 壳聚糖主要依靠其氨基、碳6位羟基和碳3位羟基与OM和CM的头部极性区发生快速结合。随后壳聚糖末端糖基单元倾向于插入OM内部,深度约1 nm,并与脂质分子脂肪酸链上的羰基形成稳定的氢键相互作用。与之相比,壳聚糖分子难以稳定地插入CM内部。壳聚糖结合对膜结构性质产生影响,主要表现在降低OM和CM的单分子脂质面积,显著减少OM和CM极性区的Ca2+和Na+数目,破坏阳离子介导的脂质间相互作用。结论 本研究证明,壳聚糖带正电的氨基基团是介导其与膜相互作用的关键,并破环脂质间的相互作...     

Analysis of the Interactions between Arp2/3 Complex and an Inhibitor Arpin by Molecular Dynamics Simulation

The Arp2/3 complex is one of the main regulators of the actin cytoskeleton and a basic molecular machine that nucleates the branched actin filaments. In this work, we studied the interaction of the Arp2/3 complex with its inhibitor, arpin, and revealed the amino-acid residues that are responsible for complex formation. The free-energy calculation for arpin binding to the Arp2/3 complex was performed using umbrella sampling. It has been shown that the dissociation constant of the Arp2/3–arpin complex is higher on average than that of Arp2/3 complexes with other inhibitors. Two arpin binding sites with different affinities were identified on the surface of the Arp2/3 complex. The mechanism of the inhibition of the Arp2/3 complex by arpin is discussed.     

抗菌肽BLFcin6与不同磷脂膜之间相互作用的分子动力学模拟研究

抗菌肽具有广谱抗菌特性,有望成为抗生素较好的替代产品.研究抗菌肽的抗菌机制,可以为新型抗菌肽的设计提供指导.无论抗菌肽采用哪种抗菌机制,其首先要稳定地吸附到细胞膜之上.因此,本文利用分子动力学模拟方法比较了抗菌小肽BLFcin6与5种不同细胞膜之间的相互作用.对这5种细胞膜而言,小肽会很快结合在POPG膜和DPPC-CHOL膜的表面,倾向于进入DPPC膜的疏水内部,与POPC膜和POPC-CHOL膜的接触很少.考察相互作用能,小肽与POPG膜的相互作用最强,主要是小肽与细胞膜亲水头部存在静电相互作用;小肽与DPPC膜的疏水尾部的相互作用较强,但受胆固醇影响,小肽只结合在DPPC-CHOL膜表面.在结合过程中,小肽N端的Arg会先结合到细胞膜上,静电相互作用在小肽锚定细胞膜的过程中起关键作用.以上研究从原子水平上解释了为什么BLFcin6小肽具有抗菌作用,哪些残基起关键作用,也为进一步开展BLFcin6小肽及其衍生小肽的研究奠定基础.