Interaction of an Amphiphilic Peptide with a Phospholipid Bilayer Surface by Molecular Dynamics Simulation Study

Abstract

Corticotropin-releasing factor (CRF) is the principal neuroregulator of adrenocorticotropic hormone (ACTH) secretion. Previous experiments have demonstrated that CRF binds avidly to the surface of single egg phosphatidylcholine vesicles and its amphiphilic secondary structure might play an important role in the function. In this study, the interaction of the residues 13–41 in human CRF with the surface of a DOPC bilayer was investigated by molecular dynamics (MD) simulation in order to understand the role of the membrane surface in the formation of the amphiphilic α helix as well as to determine the effects of the peptide on the lipid bilayer. The model used included 60 DOPC molecules, 1 helical peptide (CRF_13–41) on the bilayer surface, and explicit waters of solvation in the lipid polar head group regions, together with constant-volume periodic boundary conditions in three dimensions. The MD simulation was carried out for 510 ps. In addition, CRF_13–41, initially in a helical form, was simulated in vacuo as a control. The results indicate that while it was completely unstable in vacuo, the peptide helical form was generally maintained on the bilayer surface, but with distortions near the terminal ends. The peptide was confined to the bilayer headgroup/water region, similar to that reported from neutron diffraction measurement of tripeptides bound to the phosphatidylcholine bilayer surface (Ref 1). The amphiphilicity of the peptide matched that of the bilayer headgroup environment, with the hydrophilic side oriented toward water and the hydrophobic side making contact with the bilayer hydrocarbon core. These results support the hypothesis that the amphiphilic environment of a membrane surface is important in the induction of peptide amphiphilic α-helical secondary structure. Two major effects of the peptide on the lipids were found: the first CH₂ segment in the lipid chains was significantly disordered and the lipid headgroup distribution was broadened towards the water region.     

Early Relaxation Dynamics in the LC 13 T Cell Receptor in Reaction to 172 Altered Peptide Ligands: A Molecular Dynamics Simulation Study

The interaction between the T cell receptor and the major histocompatibility complex is one of the most important events in adaptive immunology. Although several different models for the activation process of the T cell via the T cell receptor have been proposed, it could not be shown that a structural mechanism, which discriminates between peptides of different immunogenicity levels, exists within the T cell receptor. In this study, we performed systematic molecular dynamics simulations of 172 closely related altered peptide ligands in the same T cell receptor/major histocompatibility complex system. Statistical evaluations yielded significant differences in the initial relaxation process between sets of peptides at four different immunogenicity levels.     

Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K

Abstract

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.     

Molecular Dynamics Simulation of Structural Phase Transition of AlPO4 Induced by Pressure

Abstract

The mechanism of pressure-induced phase transition of AlPO₄ has been investigated by means of a molecular dynamics method of constant temperature and pressure. A new crystalline phase with space group C2, which has not yet been experimentally found, appears by an instantaneous compression of 60, 70 and 80 GPa at 300 K. At high temperature (2500 K) and pressure (58 GPa), another new phase of AlPo₄ (y-phase), which is composed of PO₆ and AlO₆ octahedra, has been observed.     

High-Glucose-Induced Endothelial Cell Injury Is Inhibited by a Peptide Derived from Human Apolipoprotein E

Although the importance of human apolipoprotein E (apoE) in vascular diseases has clearly been established, most of the research on apoE has focused on its role in cholesterol metabolism. In view of the observation that apoE and its functional domains impact extracellular matrix (ECM) remodeling, we hypothesized that apoE could also confer protection against ECM degradation by mechanisms independent of its role in cholesterol and lipoprotein transport. The ECM degrading enzyme, heparanase, is secreted by cells as pro-heparanase that is internalized through low-density lipoprotein (LDL) receptor-related protein-1 (LRP-1) to become enzymatically active. Both apoE and pro-heparanase bind the LRP-1. We further hypothesized that an apoE mimetic peptide (apoEdp) would inhibit the production of active heparanase by blocking LRP-1-mediated uptake of pro-heparanase and thereby decrease degradation of the ECM. To test this hypothesis, we induced the expression of heparanase by incubating human retinal endothelial cells (hRECs) with high glucose (30 mM) for 72 hours. We found that elevated expression of heparanase by high glucose was associated with increased shedding of heparan sulfate (ΔHS) and the tight junction protein occludin. Treatment of hRECs with 100 µM apoEdp in the presence of high glucose significantly reduced the expression of heparanase, shedding of ΔHS, and loss of occludin as detected by Western blot analysis. Either eye drop treatment of 1% apoEdp topically 4 times a day for 14 consecutive days or intraperitoneal injection (40 mg/kg) of apoEdp daily for 14 consecutive days in an in vivo mouse model of streptozotocin-induced diabetes inhibited the loss of tight junction proteins occludin and zona occludin- 1 (ZO-1). These findings imply a functional relationship between apoE and endothelial cell matrix because the deregulation of these molecules can be inhibited by a short peptide derived from the receptor-binding region of apoE. Thus, strategies targeting ECM-degrading enzymes could be therapeutically beneficial for treating diabetic retinopathy.     

Study of the Affinity between the Protein Kinase PKA and Peptide Substrates Derived from Kemptide Using Molecular Dynamics Simulations and MM/GBSA

We have carried out a protocol in computational biochemistry including molecular dynamics (MD) simulations and MM/GBSA free energy calculations on the complex between the protein kinase A (PKA) and the specific peptide substrate Kemptide (LRRASLG). We made the same calculations on other PKA complexes that contain Kemptide derivatives (with mutations of the arginines, and with deletions of N and C-terminal amino acids). We predicted shifts in the free energy changes from the free PKA to PKA-substrate complex (ΔΔG_E→ES) when Kemptide structure is modified (we consider that the calculated shifts correlate with the experimental shifts of the free energy changes from the free PKA to the transition states (ΔΔG_E→TS) determined by the catalytic efficiency (k_cat/K_M) changes). Our results demonstrate that it is possible to predict the kinetic properties of protein kinases using simple computational biochemistry methods. As an additional benefit, these methods give detailed molecular information that permit the analysis of the atomic forces that contribute to the affinity between protein kinases and their substrates.     

Mechano-Transduction Signals Derived from Self-Assembling Peptide Nanofibers Containing Long Motif of Laminin Influence Neurogenesis in In-Vitro and In-Vivo

Astroglial scaring and limited neurogenesis are two problematic issues in recovery of spinal cord injury (SCI). In the meantime, it seems that mechanical manipulations of scaffold to inhibit astroglial scarring and improve neurogenesis is worthy of value. In the present investigation, the effect of nanofiber (gel) concentration as a mechanical-stimuli in neurogenesis was investigated. Cell viability, membrane damage, and neural differentiation derived from endometrial stem cells encapsulated into self-assembling peptide nanofiber containing long motif of laminin were assessed. Then, two of their concentrations that had no significant difference of neural differentiation potential were selected for motor neuron investigation in SCI model of rat. MTT assay data showed that nanofibers at the concentrations of 0.125 and 0.25 % w/v induced higher and less cell viability than others, respectively, while cell viability derived from higher concentrations of 0.25 % w/v had ascending trend. Gene expression results showed that noggin along with laminin motif over-expressed TH gene and the absence of noggin or laminin motif did not in all concentrations. Bcl₂ over-expression is concomitant with the decrease of nanofiber stiffness, NF⁺ cells increment, and astrogenesis inhibition and dark neuron decrement in SCI model. It seems that stiffness affects on Bcl₂ gene expression and may through β-Catenin/Wnt signaling pathway and BMP-4 inhibition decreases astrogenesis and improves neurogenesis. However, stiffness had a significant effect on upregulation of GFAP⁺ cells and motor neuron recovery in in vivo. It might be concluded that eventually there is a critical definitive point concentration that at less or higher than of it changes cell behavior and neural differentiation through different molecular pathways.     

促细胞黏附生长的自组装肽水凝胶

随着细胞与组织工程的迅猛发展,能够促进细胞黏附、生长和分化的生物材料基质支架的研究日益重要。具有生物相容性且含水量超过99％的自组装肽水凝胶因其很好地符合理想的生物材料基质支架标准而备受重视。这类自我互补的两亲寡肽含50％的带电残基,并且以交替的离子亲水性和不带电的氨基酸残基周期性重复为特征;在其寡肽的氨基末端可用直接固相合成法修饰几个短序列生物活性模体进行功能化,用以促进不同细胞的黏附生长和靶向定位。现对自组装肽水凝胶的结构特征、自组装机制、对细胞黏附生长的影响以及未来自组装肽生物材料设计的目标进行综述．     

Molecular Dynamics Simulation of Model Lipid Membranes: Structural Effects of Impurities

Abstract

The gel to fluid phase transition or ordered to disordered phase transition observed in biological membranes are simulated by using constant energy Molecular Dynamics. The surface part of the membrane is modelled as a two-dimensional matrix formed by the head groups of the phospholipid molecules. Head molecules which are modelled as three spheres fused with three force centers, interact with each other via van der Waals and Coulomb type interactions. The -so called- impurity or foreign molecule embedded in the surface represents the protein type molecule which is present in biological membranes and control its activity. It is modelled as a pentagon having one force centers in each corner. It also interacts with the surface molecules again via van der Waals and Coulomb type interactions. The surface density is kept constant in the simulations of the systems with or without impurity. Structural and orientational changes due to impurity were observed and proved by monitoring two-dimensional order parameter. It has been shown that melting of the surface or breakage of the ordering of the surface molecules becomes easier and ordered to disordered phase transition temperature was lowered by 100 K if the impurity is present.     

Conformational Analysis of a Synthetic Antimicrobial Peptide in Water and Membrane-Mimicking Solvents: A Molecular Dynamics Simulation Study

We have investigated structural and dynamic properties of the synthetic peptide hlF1-11 (GRRRSVQWCA, i.e., the first 11 N-terminal amino acids of the human lactoferrin protein) in water, 250 mM NaCl solution, 50% (V/V) water–trifluoroethanol mixture, and in the membrane mimetic 4:4:1 methanol–chloroform–water mixture. For comparison, we have also performed analogous simulations for the biologically inactive control peptide featuring Ala substitutions in the 2, 3, 6 and 9 positions of the hlF1-11 sequence. Statistical analyses of the trajectories indicate that only in the membrane-mimicking medium hlF1-11 adopts preferentially a conformation suitable to interact effectively with the membrane. In this conformation the peptide cationic region is rather flexible and elongated, while the C-terminal hydrophobic moiety appears as a more rigid hairpin-shaped loop approximately perpendicular to the cationic region. No such conformation is statistically relevant for the control peptide.     

Large Scale Characterization of the LC13 TCR and HLA-B8 Structural Landscape in Reaction to 172 Altered Peptide Ligands: A Molecular Dynamics Simulation Study

The interplay between T cell receptors (TCRs) and peptides bound by major histocompatibility complexes (MHCs) is one of the most important interactions in the adaptive immune system. Several previous studies have computationally investigated their structural dynamics. On the basis of these simulations several structural and dynamical properties have been proposed as effectors of the immunogenicity. Here we present the results of a large scale Molecular Dynamics simulation study consisting of 100 ns simulations of 172 different complexes. These complexes consisted of all possible point mutations of the Epstein Barr Virus peptide FLRGRAYGL bound by HLA-B*08:01 and presented to the LC13 TCR. We compare the results of these 172 structural simulations with experimental immunogenicity data. We found that simulations with more immunogenic peptides and those with less immunogenic peptides are in fact highly similar and on average only minor differences in the hydrogen binding footprints, interface distances, and the relative orientation between the TCR chains are present. Thus our large scale data analysis shows that many previously suggested dynamical and structural properties of the TCR/peptide/MHC interface are unlikely to be conserved causal factors for peptide immunogenicity.     

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.     

分子动力学模拟与分子生物力学

生物大分子的微观结构动力学决定其生物学功能,其力学-化学耦合规律是分子生物力学的重点关注方向。分子动力学模拟是耦合生物大分子力学-化学性质微观结构动力学基础的有效手段,其结果可用于预测结构-功能关系、指导实验设计和诠释实验结果。本文简要介绍了分子动力学模拟的方法学特点、基本工作原理及其在分子生物力学中的应用,并展望了未来可能的发展方向和应用前景。     

Molecular Dynamics Simulation and Molecular Docking Studies of Angiotensin Converting Enzyme with Inhibitor Lisinopril and Amyloid Beta Peptide

Angiotensin converting enzyme (ACE) cleaves amyloid beta peptide. So far this cleavage mechanism has not been studied in detail at atomic level. Keeping this view in mind, we performed molecular dynamics simulation of crystal structure complex of testis truncated version of ACE (tACE) and its inhibitor lisinopril along with Zn²⁺ to understand the dynamic behavior of active site residues of tACE. Root mean square deviation results revealed the stability of tACE throughout simulation. The residues Ala 354, Glu 376, Asp 377, Glu 384, His 513, Tyr 520 and Tyr 523 of tACE stabilized lisinopril by hydrogen bonding interactions. Using this information in subsequent part of study, molecular docking of tACE crystal structure with Aβ-peptide has been made to investigate the interactions of Aβ-peptide with enzyme tACE. The residues Asp 7 and Ser 8 of Aβ-peptide were found in close contact with Glu 384 of tACE along with Zn²⁺. This study has demonstrated that the residue Glu 384 of tACE might play key role in the degradation of Aβ-peptide by cleaving peptide bond between Asp 7 and Ser 8 residues. Molecular basis generated by this attempt could provide valuable information towards designing of new therapies to control Aβ concentration in Alzheimer’s patient.     

Validating Solution Ensembles from Molecular Dynamics Simulation by Wide-Angle X-ray Scattering Data

Wide-angle x-ray scattering (WAXS) experiments of biomolecules in solution have become increasingly popular because of technical advances in light sources and detectors. However, the structural interpretation of WAXS profiles is problematic, partly because accurate calculations of WAXS profiles from structural models have remained challenging. In this work, we present the calculation of WAXS profiles from explicit-solvent molecular dynamics (MD) simulations of five different proteins. Using only a single fitting parameter that accounts for experimental uncertainties because of the buffer subtraction and dark currents, we find excellent agreement to experimental profiles both at small and wide angles. Because explicit solvation eliminates free parameters associated with the solvation layer or the excluded solvent, which would require fitting to experimental data, we minimize the risk of overfitting. We further find that the influence from water models and protein force fields on calculated profiles are insignificant up to q ≈ 15 nm^?1. Using a series of simulations that allow increasing flexibility of the proteins, we show that incorporating thermal fluctuations into the calculations significantly improves agreement with experimental data, demonstrating the importance of protein dynamics in the interpretation of WAXS profiles. In addition, free MD simulations up to one microsecond suggest that the calculated profiles are highly sensitive with respect to minor conformational rearrangements of proteins, such as an increased flexibility of a loop or an increase of the radius of gyration by  <  1%. The present study suggests that quantitative comparison between MD simulations and experimental WAXS profiles emerges as an accurate tool to validate solution ensembles of biomolecules.     

Molecular Dynamics Simulation Study of DNA Triplex Formed by Mixed Sequences in Solution

Abstract

The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC)?d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usual “cut-off” of the long—range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence—dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other legends.     

A Molecular Dynamics Simulation Study of Coaxial Stacking in RNA

Abstract

We report on unrestrained molecular dynamics simulations of an RNA tetramer binding to a tetra-nucleotide overhang at the 5′-end of an RNA hairpin (nicked structure) and of the corresponding continuous hairpin with Na⁺ as counterions. The simulations lead to stable structures and in this way a structural model for the coaxially stacked RNA hairpin is generated. The stacking interface in the coaxially stacked nicked hairpin structure is characterized by a reduced twist and shift and a slightly increased propeller twist as compared to the continuous system. This leads to an increased overlap between C22 and G23 in the stacking interface of the nicked structure. In the simulations the continuous RNA hairpin has an almost straight helical axis. On the other hand, the corresponding axis for the nicked structure exhibits a marked kink of 39°. The stacking interface exhibits no increased flexibility as compared to the corresponding base pair step in the continuous structure.     

Disruption Mechanism in the Helix of SPF Peptide by Interchanging E5 and K10 Residues: Inference from Molecular Dynamics Study

Abstract

Seminalplasmin (SPLN) is a 47-residue peptide (SDEKASPDKHHRFSLSRYAKLANR LANPKLLETFLSKWIGDRGNRSV) from bovine seminal plasma. It has broad spectrum antimicrobial activity, without any hemolytic activity. The 28–40 segment of SPLN with the sequence PKLLETFLSKWIG, designated as SPF, is the most hydrophobic stretch of SPLN and primarily responsible for the membrane-perturbing activity of SPLN. It was reported that SPF has a helical structure and the interchange of E5 and K10 residues disrupted the helical structure. The present paper reports a possible mechanism of disruption of the helical structure of SPF peptide during the interchange of E5 and K10 residues. The result is based on simulated annealing and molecular dynamics simulation studies on SPF and its four analogues with K10E, K10D, E5K, and E5K & K10E substitutions. It showed that K10 residue has a critical role in maintaining the highest helical content and the positions of charged residues are also very important for maintaining the helical structure of the SPF peptide. Formation of some new long-range hydrogen bonds and the rupture of some shortrange hydrogen bonds involving the tenth residue led to the disruption of helical structure of SPF peptide when E5 and K10 residues are interchanged.