Molecular dynamics simulations of the influenza hemagglutinin fusion peptide in micelles and bilayers: conformational analysis of peptide and lipids

Molecular dynamics simulations of the influenza hemagglutinin fusion peptide in two differently sized dodecylphosphocholine micelles and a palmitoyl oleoyl phosphatidylcholine bilayer were generated to analyze the influence of the environment. Four independent trajectories (5 ns each for the bilayer, and 2 ns each for the micelles) were generated for each system. The peptide lies at the surface of the micelles, while its N-terminal region inserts deeply in the bilayer. This leads to a substantial increase of the solvation and rigidity of the peptide in micelles as compared to the bilayer. The average structures, nevertheless, are similar in all three systems and agree reasonably with micelle-based NMR structures. When in the bilayer, the peptide increases the chain gauche population and area of adjacent lipids in the same binding leaflet, while it has the opposite effect for the nearby lipids of the other leaflet. These changes, which occur spontaneously to fill voids and defects, cause a decrease in the thickness of the membrane in the neighborhood of the peptide. They would be expected to promote positive curvature, as consistent with the formation of the convex bulge, or "nipple", in the initial stage of membrane fusion. An extension of the classical surfactant theory of Israelachvili based on shapes is proposed to introduce the concept of a "dynamically induced shape" of the membrane lipids by the peptide.     

一种水蛭素类融合蛋白与凝血酶作用的动力学模拟

通过将凝血酶抑制剂水蛭素的C端20肽片段嫁接到血小板结合蛋白AnnexinⅤ上,可以期望获得既具有抗凝血活性,同时又具有导向性的新型工程蛋白质分子.利用计算机辅助分子设计手段模拟了该融合蛋白的分子结构,并对该融合蛋白对凝血酶的抑制作用进行了分子动力学模拟,得到了支持上述想法的结果.     

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Interaction between cystatin-derived peptides and papain

The interaction between papain and synthetic peptides which tentatively mimic cystatin surfaces was investigated both enzymatically and structurally. Measurements of dissociation equilibrium constants for the interaction of papain with these peptides modified by successive deletions or substitutions demonstrated that the QVVAG segment, which is highly conserved throughout members of the cystatin superfamily, is essential for the interaction. The glycylcontaining (N-terminal) fragments and PW-containing (C-terminal) fragments were found to be of lesser importance, since each could be deleted without significantly modifying the interaction. These fragments improved the stability of the interacting QVVAG region, which appeared to be substrate-like in all peptides tested, as it was cleaved at the A-G bond upon peptide-papain interaction. Replacement of the A residue at the scissile bond of the QVVAG by a blocked cysteinyl residue reduced the rate of cleavage of the susceptible bond and therefore shifted the resulting peptide from a substrate to an inhibitor. Derivatization of this substituted peptide at its N- and C-terminal ends by fluoresceinyl groups resulted in a dramatic decrease in theK_i to 0.5 µM. This improvement in the inhibitory properties of the substituted and derivatized peptides was correlated with structural changes as analyzed by molecular dynamic calculations. The results were compared to those proposed for the mechanism of inhibition by natural inhibitors of the cystatin superfamily.     

Interaction between dengue virus fusion peptide and lipid bilayers depends on peptide clustering

Dengue fever is one of the most widespread tropical diseases in the world. The disease is caused by a virus member of the Flaviviridae family, a group of enveloped positive sense single-stranded RNA viruses. Dengue virus infection is mediated by virus glycoprotein E, which binds to the cell surface. After uptake by endocytosis, this protein induces the fusion between viral envelope and endosomal membrane at the acidic environment of the endosomal compartment. In this work, we evaluated by steady-state and time-resolved fluorescence spectroscopy the interaction between the peptide believed to be the dengue virus fusion peptide and large unilamellar vesicles, studying the extent of partition, fusion capacity and depth of insertion in membranes. The roles of the bilayer composition (neutral and anionic phospholipids), ionic strength and pH of the medium were also studied. Our results indicate that dengue virus fusion peptide has a high affinity to vesicles composed of anionic lipids and that the interaction is mainly electrostatic. Both partition coefficient and fusion index are enhanced by negatively charged phospholipids. The location determined by differential fluorescence quenching using lipophilic probes demonstrated that the peptide is in an intermediate depth in the hemilayers, in-between the bilayer core and its surface. Ultimately, these data provide novel insights on the interaction between dengue virus fusion peptide and its target membranes, namely, the role of oligomerization and specific types of membranes.     

Computational screening of fatty acid synthase inhibitors against thioesterase domain

Thioesterase (TE) domain of fatty acid synthase (FAS) is an attractive therapeutic target for design and development of anticancer drugs. In this present work, we search for the potential FAS inhibitors of TE domain from the ZINC database based on similarity search using three natural compounds as templates, including flavonoids, terpenoids, and phenylpropanoids. Molecular docking was used to predict the interaction energy of each screened ligand compared to the reference compound, which is methyl γ-linolenylfluorophosphonate (MGLFP). Based on this computational technique, rosmarinic acid and its eight analogs were observed as a new series of potential FAS inhibitors, which showed a stronger binding affinity than MGLFP. Afterward, nine docked complexes were studied by molecular dynamics simulations for investigating protein–ligand interactions and binding free energies using MM-PB(GB)SA, MM-3DRISM-KH, and QM/MM-GBSA methods. The binding free energy calculation indicated that the ZINC85948835 (R34) displayed the strongest binding efficiency against the TE domain of FAS. There are eight residues (S2308, I2250, E2251, Y2347, Y2351, F2370, L2427, and E2431) mainly contributed for the R34 binding. Moreover, R34 could directly form hydrogen bonds with S2308, which is one of the catalytic triad of TE domain. Therefore, our finding suggested that R34 could be a potential candidate as a novel FAS-TE inhibitor for further drug design.     

The structure of the major transition state for folding of an FF domain from experiment and simulation

We have analysed the transition state of folding of the four-helix FF domain from HYPA/FBP11 by high-resolution experiment and simulation as part of a continuing effort to understand the principles of folding and the refinement of predictive methods. The major transition state for folding was subjected to a Phi-value analysis utilising 50 mutants. The transition state contained a nucleus for folding centred around the end of helix 1 (H1) and the beginning of helix 2 (H2). Secondary structure in this region was fully formed (PhiF=0.9-1) and tertiary interactions were well developed. Interactions in the distal part of the native structure were weak (PhiF=0-0.2). The hydrophobic core and other parts of the protein displayed intermediate Phi-values, suggesting that interactions coalesce as the end of H1 and beginning of H2 are in the process of being formed. The distribution of Phi-values resembled that of barnase, which folds via an intermediate, rather than that of CI2 which folds by a concerted nucleation-condensation mechanism. The overall picture of the transition state structure identified in molecular dynamics simulations is in qualitative agreement, with the turn connecting H1 and H2 being formed, a loosened core, and H4 partially unfolded and detached from the core. There are some differences in the details and interpretation of specific Phi-values.     

Docking and molecular dynamics studies of peptide inhibitors of ornithine decarboxylase: a rate-limiting enzyme for the metabolism of Fusarium solani

Fusarium solani causes a wide variety of diseases in plants. Polyamine biosynthesis is responsible for the growth and pathogenicity of the fungus. The initial step of this pathway involves the decarboxylation of ornithine to putrescine, and is catalyzed by the enzyme ornithine decarboxylase (ODC). Inhibiting this process may be a promising approach for the management of fungal disease in various crops. Therefore, there is a need to develop inhibitors of ODC that have higher binding capacity than ornithine. Fifteen peptides were designed and modeled based on physicochemical properties of residues in the active site of ODC. The peptide GLIWGNGPF showed the highest dock score. It is assumed that the de novo design of peptides could be a potential approach to inhibit polyamine biosynthesis. Molecular dynamics studies make an important contribution to understanding the effect of the binding of peptides and the stability of an ODC-peptide complex system.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:8.     

HIV-1跨膜蛋白gp41与N-取代吡咯衍生物的结合模式研究

采用分子对接,分子动力学(MD)模拟和分子力学/泊松-波尔兹曼溶剂可有面积方法与分子力学/广义伯恩溶剂可及面积方法(MM-PBSA/MM-GBSA),预测两种N-取代吡咯衍生物与HIV-1 跨膜蛋白gp41疏水口袋的结合模式与作用机理．分子对接采用多种受体构象,并从结果中选取几种可能的结合模式进行MD 模拟,然后通过MM-PBSA计算结合能的方法识别最优的结合模式. MM-PBSA计算结果表明,范德华相互作用是结合的主要驱动力,而极性相互作用决定了配体在结合过程中的取向．进一步的结合能分解显示,配体的羧基与gp41残基Arg579的静电相互作用对结合有重要贡献．上述工作为进一步优化N-取代吡咯衍生物类的HIV-1融合抑制剂建立了良好的理论基础．