Folding profiles of antimicrobial scorpion venom-derived peptides on hydrophobic surfaces: a molecular dynamics study

Abstract

Most helical antimicrobial peptides (AMPs) are usually unfolded in aqueous solution; however they acquire their secondary structure in the presence of a hydrophobic environment such as lipid membranes. Being the biological membranes the main target of many AMPs it is necessary to understand their way of action. Pandinin 2 (Pin2) is an alpha-helical AMP isolated from the venom of the African scorpion Pandinus imperator which shows high antimicrobial activity against Gram-positive bacteria and it is less active against Gram-negative bacteria, nevertheless, it has strong hemolytic activity. Its chemically synthesized Pin2GVG analog has low hemolytic activity while keeping its antimicrobial activity. With the aim of exploring the partition and subsequent folding of these peptides, in this work we report the results of extensive molecular dynamics simulations of Pin2 and Pin2GVG peptides in the presence of 2 hydrophobic environments such as dodecyl-phosphocholine (DPC) micelle and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocoline (POPC) membrane. Our results indicate that Pin2 folds in DPC with a 79% of alpha-helical content, which is in agreement with the experimental results, while in POPC it has 62.5% of alpha-helical content. On the other hand, Pin2GVG presents a higher percentage of alpha-helical structure in POPC and a smaller content in DPC when compared with Pin2. These results can help to better choose the starting structures in future molecular dynamics simulations of AMPs, because these peptides can adopt slightly different conformations depending on the hydrophobic environment.

Communicated by Ramaswamy H. Sarma     

Interaction of alginate single-chain polyguluronate segments with mono- and divalent metal cations: a comparative molecular dynamics study

The interaction of a set of monovalent (Na⁺, K⁺) and divalent (Mg²⁺, Ca²⁺) metal cations with single-chain polyguluronate (periodic chain based on a dodecameric repeat unit, 2₁-helical conformation) is investigated using explicit-solvent molecular dynamics simulations (at 300 K and 1 bar). A total of 14 (neutralising) combinations of the different ions are considered (single type of cation or simultaneous presence of two types of cation, either in the presence or absence of chloride anions). The main observations are: (1) the chain conformation and intramolecular hydrogen bonding is insensitive to the counter-ion environment; (2) the binding of the cations is essentially non-specific for all ions considered (counter-ion atmosphere confined within a cylinder of high ionic density, but no well-defined binding sites); (3) the density and tightness of the distributions of the different cations within the counter-ion atmosphere follow the approximate sequence Ca²⁺>Mg²⁺>K⁺～Na⁺; (4) the solvent-separated binding of the cations to the carboxylate groups of the chain is frequent, and its occurrence follows the approximate sequence K⁺>Na⁺>Ca²⁺>Mg²⁺ (contact-binding events as well as the binding of a cation to multiple carboxylate groups are very infrequent); and (5) the counter-ion atmosphere typically leads to a complete screening of the chain charge within 1.0–1.2 nm of the chain axis and, for most systems, to a charge reversal at about 1.5 nm (i.e. the effective chain charge becomes positive at this distance and as high in magnitude as one-quarter of the bare chain charge, before slowly decreasing to zero). These findings agree well (in a qualitative sense) with available experimental data and predictions from simple analytical models, and provide further insight concerning the nature of alginate–cation interactions in aqueous solution.     

Antimicrobial peptide RP-1 structure and interactions with anionic versus zwitterionic micelles

Topologically, platelet factor-4 kinocidins consist of distinct N-terminal extended, C-terminal helical, and interposing gamma-core structural domains. The C-terminal alpha-helices autonomously confer direct microbicidal activity, and the synthetic antimicrobial peptide RP-1 is modeled upon these domains. In this study, the structure of RP-1 was assessed using several complementary techniques. The high-resolution structure of RP-1 was determined by NMR in anionic sodium dodecyl sulfate (SDS) and zwitterionic dodecylphosphocholine (DPC) micelles, which approximate prokaryotic and eukaryotic membranes, respectively. NMR data indicate the peptide assumes an amphipathic alpha-helical backbone conformation in both micelle environments. However, small differences were observed in the side-chain orientations of lysine, tyrosine, and phenylalanine residues in SDS versus DPC environments. NMR experiments with a paramagnetic probe indicated differences in positioning of the peptide within the two micelle types. Molecular dynamics (MD) simulations of the peptide in both micelle types were also performed to add insight into the peptide/micelle interactions and to assess the validity of this technique to predict the structure of peptides in complex with micelles. MD independently predicted RP-1 to interact only peripherally with the DPC micelle, leaving its spherical shape intact. In contrast, RP-1 entered deeply into and significantly distorted the SDS micelle. Overall, the experimental and MD results support a preferential specificity of RP-1 for anionic membranes over zwitterionic membranes. This specificity likely derives from differences in RP-1 interaction with distinct lipid systems, including subtle differences in side chain orientations, rather than gross changes in RP-1 structure in the two lipid environments.     

Molecular dynamics studies of the antimicrobial peptides piscidin 1 and its mutants with a DOPC lipid bilayer

Piscidin 1 (Pis‐1) has a high broad‐spectrum activity against bacteria, fungi, and viruses but it also has a moderate hemolytic activities. To improve the antibacterial activity and to reduce toxicity, mutants Pis‐1AA (G8A/G13A double mutant) and Pis‐1PG (G8P mutant) have been designed based on the crystal structure of Pis‐1. Eighteen independent molecular dynamics (MD) simulations of Pis‐1 and its mutants with membranes are conducted in this article. Furthermore, 60 independent MD simulations of three peptides in water box have also been discussed for comparison. The results indicate that the unfolding process starts at the middle of the peptide. Pis‐1 disrupts easily in the region of Val10‐Lys14. Pis‐1PG has a flexible N‐terminal region, and the interaction between N‐terminal and C‐terminal is very weak. Pis‐1AA has the most stable helical structure. In addition, percentage of native contacts and hydrogen bonds analysis are also performed. Lipid‐peptide interaction analysis suggests that Pis‐1 and Pis‐1AA has a stronger interaction with the zwitterionic dioleoylphosphatidylcholine (DOPC) lipid bilayer than Pis‐1PG. When compared with the results of peptide with membrane, peptides are unstable and unfolding quickly in water solution. Our results are applicable in examining diversities on hemolytic, antibacterial, and selectivity of antimicrobial peptides. © 2012 Wiley Periodicals, Inc. Biopolymers 97:998–1009, 2012.     

Effects of cholesterol concentration on the interaction of cytarabine with lipid membranes: a molecular dynamics simulation study

Liposomal cytarabine, DepoCyt, is a chemotherapy agent which is used in cancer treatment. This form of cytarabine has more efficacy and fewer side effects relative to the other forms. Since DepoCyt contains the cytarabine encapsulated within phosphatidylcholine and the sterol molecules, we modeled dioleoylphosphatidylcholine (DOPC)/cholesterol bilayer membrane as a carrier for cytarabine to study drug–bilayer interactions. For this purpose, we performed a series of united-atom molecular dynamics (MD) simulations for 25?ns to investigate the interactions between cytarabine and cholesterol-containing DOPC lipid bilayers. Only the uncharged form of cytarabine molecule was investigated. In this study, different levels of the cholesterol content (0, 20, and 40%) were used. MD simulations allowed us to determine dynamical and structural properties of the bilayer membrane and to estimate the preferred location and orientation of the cytarabine molecule inside the bilayer membrane. Properties such as membrane thickness, area per lipid, diffusion coefficient, mass density, bilayer packing, order parameters, and intermolecular interactions were examined. The results show that by increasing the cholesterol concentration in the lipid bilayers, the bilayer thickness increases and area per lipid decreases. Moreover, in accordance with the experiments, our calculations show that cholesterol molecules have ordering effect on the hydrocarbon acyl chains. Furthermore, the cytarabine molecule preferentially occupies the polar region of the lipid head groups to form specific interactions (hydrogen bonds). Our results fully support the experimental data. Our finding about drug–bilayer interaction is crucial for the liposomal drug design.     

Interfacial water on crystalline silica: a comparative molecular dynamics simulation study

Understanding the properties of interfacial water at solid–liquid interfaces is important in a wide range of applications. Molecular dynamics is becoming a widespread tool for this purpose. Unfortunately, however, the results of such studies are known to strongly depend on the selection of force fields. It is, therefore, of interest to assess the extent by which the implemented force fields can affect the predicted properties of interfacial water. Two silica surfaces, with low and high surface hydroxyl density, respectively, were simulated implementing four force fields. These force fields yield different orientation and flexibility of surface hydrogen atoms, and also different interaction potentials with water molecules. The properties for interfacial water were quantified by calculating contact angles, atomic density profiles, surface density distributions, hydrogen bond density profiles and residence times for water near the solid substrates. We found that at low surface density of hydroxyl groups, the force field strongly affects the predicted contact angle, while at high density of hydroxyl groups, water wets all surfaces considered. From a molecular-level point of view, our results show that the position and intensity of peaks observed from oxygen and hydrogen atomic density profiles are quite different when different force fields are implemented, even when the simulated contact angles are similar. Particularly, the surfaces simulated by the CLAYFF force field appear to attract water more strongly than those simulated by the Bródka and Zerda force field. It was found that the surface density distributions for water strongly depend on the orientation of surface hydrogen atoms. In all cases, we found an elevated number of hydrogen bonds formed between interfacial water molecules. The hydrogen bond density profile does not depend strongly on the force field implemented to simulate the substrate, suggesting that interfacial water assumes the necessary orientation to maximise the number of water–water hydrogen bonds irrespectively of surface properties. Conversely, the residence time for water molecules near the interface strongly depends on the force field and on the flexibility of surface hydroxyl groups. Specifically, water molecules reside for longer times at contact with rigid substrates with high density of hydroxyl groups. These results should be considered when comparisons between simulated and experimental data are attempted.     

Principles of carbopeptoid folding: a molecular dynamics simulation study.

The conformational spaces of five oligomers of tetrahydrofuran-based carbopeptoids in chloroform and dimethyl sulfoxide were investigated through nine molecular dynamics simulations. Prompted by nuclear magnetic resonance experiments that indicated various stable folds for some but not all of these carbopeptoids, their folding behaviour was investigated as a function of stereochemistry, chain length and solvent. The conformational distributions of these molecules were analysed in terms of occurrence of hydrogen bonds, backbone torsional-angle distributions, conformational clustering and solute configurational entropy. While a cis-linkage across the tetrahydrofuran ring favours right-handed helical structures, a trans-linkage results in a larger conformational variability. Intra-solute hydrogen bonding is reduced with increasing chain length and with increasing solvent polarity. Solute configurational entropies confirm the picture obtained: they are smaller for cis- than for trans-linked peptides, for chloroform than for dimethyl sulfoxide as solvent and for shorter peptide chains. The simulations provide an atomic picture of molecular conformational variability that is consistent with the available experimental data.     

Membrane protein dynamics versus environment: simulations of OmpA in a micelle and in a bilayer

The bacterial outer membrane protein OmpA is one of the few membrane proteins whose structure has been solved both by X-ray crystallography and by NMR. Crystals were obtained in the presence of detergent, and the NMR structure is of the protein in a detergent micelle. We have used 10 ns duration molecular dynamics simulations to compare the behaviour of OmpA in a detergent micelle and in a phospholipid bilayer. The dynamic fluctuations of the protein structure seem to be ca 1.5 times greater in the micelle environment than in the lipid bilayer. There are subtle differences between the nature of OmpA-detergent and OmpA-lipid interactions. As a consequence of the enhanced flexibility of the OmpA protein in the micellar environment, side-chain torsion angle changes are such as to lead to formation of a continuous pore through the centre of the OmpA molecule. This may explain the experimentally observed channel formation by OmpA.     

Mechanical properties of armchair silicene nanoribbons with edge cracks: a molecular dynamics study

Silicene has been proven to be a promising material with attractive electronic properties. During the synthesis of silicene, structural defects such as edge crack are likely to be generated and such defects in silicene have impacts on its properties. Herein, molecular dynamics simulations were performed to investigate the mechanical properties of the armchair silicene nanoribbons (ASiNRs) with edge cracks. Our results showed that the mechanical properties of the ASiNRs decrease because of the existence of edge crack. Both the pristine ASiNRs and the ASiNRs with edge cracks show brittle fracture behaviours. The crack length plays an important role in determining the critical strain and fracture strength of the ASiNRs. Moreover, we investigated the effects of strain rate and temperature on the mechanical properties of the ASiNRs with edge cracks. We observed that the increasing strain rate increases the critical strain and fracture strength while decreasing the Young’s modulus. Low-strain rates also changes the expanded directions of cracks in the ASiNRs. We also found that the increasing temperature could significantly decrease the mechanical properties of the ASiNRs with edge cracks.     

The molecular dynamics of cholesterol in bilayer membranes: A deuterium NMR study

The ²H-NMR spectra of selectively deuterated cholesterol, intercalated in egg phosphatidyl-choline, were examined. The orientation of the axis of motional averaging was calculated using the observed quadrupole splittings and the atomic coordinates. With the known orientation of the rotation axis, quadrupole splittings observed for deuterium labels on cholesterol can be related to the molecular order parameter of the sterol. In addition, knowledge of the axis orientation allows prediction of the magnitudes of quadrupole splittings for deuterium at other positions, which is useful in the choice of labelling for particular applications. Finally, preliminary relaxation time measurements yield information on the rates of anisotropic motion of cholesterol in bilayer membranes.     

NMR structures and molecular dynamics simulation of hylin‐a1 peptide analogs interacting with micelles

Antimicrobial peptides are recognized candidates with pharmaceutical potential against epidemic emerging multi‐drug resistant bacteria. In this study, we use nuclear magnetic resonance spectroscopy and molecular dynamics simulations to determine the unknown structure and evaluate the interaction with dodecylphosphatidylcholine (DPC) and sodium dodecylsulphate (SDS) micelles with three W⁶‐Hylin‐a1 analogs antimicrobial peptides (HyAc, HyK, and HyD). The HyAc, HyK, and HyD bound to DPC micelles are all formed by a unique α‐helix structure. Moreover, all peptides reach the DPC micelles' core, which thus suggests that the N‐terminal modifications do not influence the interaction with zwiterionic surfaces. On the other hand, only HyAc and HyK peptides are able to penetrate the SDS micelle core while HyD remains always at its surface. The stability of the α‐helical structure, after peptide‐membrane interaction, can also be important to the second step of peptide insertion into the membrane hydrophobic core during permeabilization. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.     

A study of comparative modelling,simulation and molecular dynamics of CXCR3 receptor with lipid bilayer

The G-coupled receptors seen on the cell surface are composites with a lipid bilayer. The chemokines are kind of G-coupled receptor which majorly involved in the activation and downstream signalling of the cell. In general, many G-coupled receptors lack their 3D structures which become a hurdle in the drug designing process. In this study, comparative modelling of the CXCR3 receptor was carried out, structure evaluation was done using various tools and softwares. Additionally, molecular dynamics and docking were performed to prove the structural quality and architecture. Interestingly, the studies like toggle switch mechanism, lipid dynamics, virtual screening were carried out to find the potent antagonist for the CXCR3 receptor. During virtual screening 14,303 similar molecules were retrieved among them only four compounds have an ability to interact with a crucial amino acid residue of an antagonist. Hence, these screened compounds can serve as a drug candidate for a CXCR3 receptor, but further in vitro and in vivo studies are ought to do to prove its same efficacy.     

Conformational dynamics of the mitochondrial ADP/ATP carrier: a simulation study

The mitochondrial ADP/ATP carrier is a six helix bundle membrane transport protein, which couples the exit of ATP from the mitochondrial matrix to the entry of ADP. Extended (4×20 ns) molecular dynamics simulations of the carrier, in the presence and absence of bound inhibitor (carboxyatractyloside), have been used to explore the conformational dynamics of the protein in a lipid bilayer environment, in the presence and absence of the carboxyatractyloside inhibitor. The dynamic flexibility (measured as conformational drift and fluctuations) of the protein is reduced in the presence of bound inhibitor. Proline residues in transmembrane helices H1, H3 and H5 appear to form dynamic hinges. Fluctuations in inter-helix salt bridges are also observed over the time course of the simulations. Inhibitor-protein and lipid-protein interactions have been characterised in some detail. Overall, the simulations support a transport mechanism in which flexibility about the proline hinges enables a transition between a ‘closed’ and an ‘open’ pore-like state of the carrier protein.     

Epigallocatechin-3-gallate location and interaction with late endosomal and plasma membrane model membranes by molecular dynamics

Epigallocatechin-3-gallate (EGCG) is the most abundant polyphenol in green tea and it has been reported to have many beneficial properties against many different types of illnesses and infections. However, the exact mechanism/s underlying its biological effects are unknown. It has been previously shown that EGCG is capable of binding to and disrupting the membrane, so that some of its effects on biological systems could be ascribed to its capacity to incorporate into the biological membrane and modulate its structure. In this work, we have used atomistic molecular dynamics (MD) to discern the location and orientation of EGCG in model membranes and the possible existence of specific interactions with membrane lipids. For that goal, we have used in our simulation two complex model membranes, one resembling the plasma membrane (PM) and the other one the late endosome (LE) membrane. Our results support that EGCG tends to associate with the membrane and exists inside it in a relatively stable and steady location with a low propensity to be associated with other EGCG molecules. Interestingly, EGCG forms hydrogen bonds with POPC and POPE in the PM system but POPC and BMP and no POPE in the LE. These data suggest that the broad beneficial effects of EGCG could be mediated, at least in part, through its membranotropic effects and therefore membrane functioning.

Communicated by Ramaswamy H. Sarma     

One membrane protein,two structures and six environments: a comparative molecular dynamics simulation study of the bacterial outer membrane protein PagP

PagP is a bacterial outer membrane protein consisting of an 8 stranded transmembrane β-barrel and an N-terminal α-helix. It is an enzyme which catalyses transfer of a palmitoyl chain from a phospholipid to lipid A. Molecular dynamics simulations have been used to compare the dynamic behaviour in simulations starting from two different structures (X-ray vs. NMR) and in six different environments (detergent micelles formed by dodecyl phosphocholine and by octyl glucoside, vs. four species of phospholipid bilayer). Analysis of interactions between the protein and its environment reveals the role played by the N-terminal α-helix, which interacts with the lipid headgroups to lock the PagP molecule into the bilayer. The PagP β-barrel adopts a tilted orientation in lipid bilayers, facilitating access of lipid tails into the mouth of the central binding pocket. In simulations starting from the X-ray structure in lipid bilayer, the L1 and L2 loops move towards one another, leading to the formation of a putative active site by residues H33, D76 and S77 coming closer together.     

A comparative molecular dynamics study on BACE1 and BACE2 flap flexibility

Beta-amyloid precursor protein cleavage enzyme1 (BACE1) and beta-amyloid precursor protein cleavage enzyme2 (BACE2), members of aspartyl protease family, are close homologs and have high similarity in their protein crystal structures. However, their enzymatic properties are different, which leads to different clinical outcomes. In this study, we performed sequence analysis and all-atom molecular dynamic (MD) simulations for both enzymes in their ligand-free states in order to compare their dynamical flap behaviors. This is to enhance our understanding of the relationship between sequence, structure and the dynamics of this protein family. Sequence analysis shows that in BACE1 and BACE2, most of the ligand-binding sites are conserved, indicative of their enzymatic property as aspartyl protease members. The other conserved residues are more or less unsystematically localized throughout the structure. Herein, we proposed and applied different combined parameters to define the asymmetric flap motion; the distance, d₁, between the flap tip and the flexible region; the dihedral angle, φ, to account for the twisting motion and the TriCα angle, θ₂ and θ₁. All four combined parameters were found to appropriately define the observed “twisting” motion during the flaps different conformational states. Additional analysis of the parameters indicated that the flaps can exist in an ensemble of conformations, i.e. closed, semi-open and open conformations for both systems. However, the behavior of the flap tips during simulations is different between BACE1 and BACE2. The BACE1 active site cavity is more spacious as compared to that of BACE2. The analysis of 10S loop and 113S loop showed a similar trend to that of flaps, with the BACE1 loops being more flexible and less stable than those of BACE2. We believe that the results, methods and perspectives highlighted in this report would assist researchers in the discovery of BACE inhibitors as potential Alzheimer’s disease therapies.