Simulation of voltage-dependent interactions of α-helical peptides with lipid bilayers

Pore formation in lipid bilayers by channel-forming peptides and toxins is thought to follow voltage-dependent insertion of amphipathic α-helices into lipid bilayers. We have developed an approximate potential for use within the CHARMm molecular mechanics program which enables one to simulate voltage-dependent interaction of such helices with a lipid bilayer. Two classes of helical peptides which interact with lipid bilayers have been studied: (a) δ-toxin, a 26 residue channel-forming peptide from Staphylococcus aureus; and (b) synthetic peptides corresponding to the α5 and α7 helices of the pore-forming domain of Bacillus thuringiensis CryIIIA δ-endotoxin. Analysis of δ-toxin molecular dynamics (MD) simulations suggested that the presence of a transbilayer voltage stabilized the inserted location of δ-toxin helices, but did not cause insertion per se. A series of simulations for the α5 and α7 peptides revealed dynamic switching of the α5 helix between a membrane-associated and a membrane-inserted state in response to a transbilayer voltage. In contrast the α7 helix did not exhibit such switching but instead retained a membrane associated state. These results are in agreement with recent experimental studies of the interactions of synthetic α5 and α7 peptides with lipid bilayers.     

Free energy landscapes of two model peptides: α-helical and β-hairpin peptides explored with Brownian dynamics simulation

We applied an atomistic Brownian dynamics (BD) simulation with multiple time step method for the folding simulation of a 13-mer α-helical peptide and a 12-mer β-hairpin peptide, giving successful folding simulations. In this model, the driving energy contribution towards folding came from both electrostatic and van der Waals interactions for the α-helical peptide and from van der Waals interactions for the β-hairpin peptide. Although, many non-native structures having the same or lower energy than that of native structure were observed, the folded states formed the most populated cluster when the structures obtained by the BD simulations were subjected to the cluster analysis based on distance-based root mean square deviation of side-chains between different structures. This result indicates that we can predict the native structures from conformations sampled by BD simulation.     

Effect of hydrocarbon stapling on the properties of α-helical antimicrobial peptides isolated from the venom of hymenoptera

The impact of inserting hydrocarbon staples into short α-helical antimicrobial peptides lasioglossin III and melectin (antimicrobial peptides of wild bee venom) on their biological and biophysical properties has been examined. The stapling was achieved by ring-closing olefin metathesis, either between two S-2-(4′-pentenyl) alanine residues (S ₅) incorporated at i and i + 4 positions or between R-2-(7′-octenyl) alanine (R ₈) and S ₅ incorporated at the i and i + 7 positions, respectively. We prepared several lasioglossin III and melectin analogs with a single staple inserted into different positions within the peptide chains as well as analogs with double staples. The stapled peptides exhibited a remarkable increase in hemolytic activity, while their antimicrobial activities decreased. Some single stapled peptides showed a higher resistance against proteolytic degradation than native ones, while the double stapled analogs were substantially more resistant. The CD spectra of the singly stapled peptides measured in water showed only a slightly better propensity to form α-helical structure when compared to native peptides, whereas the doubly stapled analogs exhibited dramatically enhanced α-helicity.     

Structure and dynamics of the lipid modifications of a transmembrane α-helical peptide determined by ²H solid-state NMR spectroscopy

The fusion of biological membranes is mediated by integral membrane proteins with α-helical transmembrane segments. Additionally, those proteins are often modified by the covalent attachment of hydrocarbon chains. Previously, a series of de novo designed α-helical peptides with mixed Leu/Val sequences was presented, mimicking fusiogenically active transmembrane segments in model membranes (Hofmann et al., Proc. Natl. Acad. Sci. USA 101 (2004) 14776-14781). From this series, we have investigated the peptide LV16 (KKKW LVLV LVLV LVLV LVLV KKK), which was synthesized featuring either a free N-terminus or a saturated N-acylation of 2, 8, 12, or 16 carbons. We used 2H and 31P NMR spectroscopy to investigate the structure and dynamics of those peptide lipid modifications in POPC and DLPC bilayers and compared them to the hydrocarbon chains of the surrounding membrane. Except for the C2 chain, all peptide acyl chains were found to insert well into the membrane. This can be explained by the high local lipid concentrations the N-terminal lipid chains experience. Further, the insertion of these peptides did not influence the membrane structure and dynamics as seen from the 2H and 31P NMR data. In spite of the fact that the longer acyl chains insert into the membrane, they do not adapt their lengths to the thickness of the bilayer. Even the C16 lipid chain on the peptide, which could match the length of the POPC palmitoyl chain, exhibited lower order parameters in the upper chain, which get closer and finally reach similar values in the lower chain region. 2H NMR square law plots reveal motions of slightly larger amplitudes for the peptide lipid chains compared to the surrounding phospholipids. In spite of the significantly different chain lengths of the acylations, the fraction of gauche defects in the inserted chains is constant.     

The effect of phosphorylation on the conformation of oligo-peptides with Ser–Pro motif: a molecular dynamics simulation

Model tetrapeptide system was designed to investigate the cis/trans isomerization of peptidyl-prolyl imide bond of Ser–Pro motif. To establish the side-chain O-phosphorylation effect in regulating the peptides conformations, molecular dynamics (MD) simulations where carried out on the designed tetrapeptides and their corresponding phosphorylated forms by MD Insight II Discovery3 approach. The most stable configurations and the statistic cis/trans concentration distribution demonstrated that the phosphorylation evidently influences the peptidyl-prolyl imide bond isomerization and works as a key effect in regulating the peptide conformations. The charge state and the site provided for the charge of the phosphate moiety might be an important key. The results also demonstrated that phosphorylation changes the cis conformation ratio of the peptide and the maximum cis value is obtained when the phosphate group has no negative charge.     

Effect of disulphide bond position on salt resistance and LPS-neutralizing activity of α-helical homo-dimeric model antimicrobial peptides

To investigate the effects of disulphide bond position on the salt resistance and lipopolysaccharide (LPS)-neutralizing activity of α-helical homo-dimeric antimicrobial peptides (AMPs), we synthesized an α-helical model peptide (K6L4W1) and its homo-dimeric peptides (di-K(6)L(4)W(1)-N, di-K(6)L(4)W(1)-M, and di-K(6)L(4)W(1)-C) with a disulphide bond at the N-terminus, the central position, and the C-terminus of the molecules, respectively. Unlike (6)L(4)W(1) and di-K(6)L(4)W(1)-M, the antimicrobial activity of di-K(6)L(4)W(1)-N and di-K(6)L(4)W(1)-C was unaffected by 150 mM NaCl. Both di-K(6)L(4)W(1)-N and di-K(6)L(4)W(1)-C caused much greater inhibitory effects on nitric oxide (NO) release in LPS-induced mouse macrophage RAW 264.7 cells, compared to di-K(6)L(4)W(1)-M. Taken together, our results indicate that the presence of a disulphide bond at the N- or C-terminus of the molecule, rather than at the central position, is more effective when designing salt-resistant α-helical homo-dimeric AMPs with potent antimicrobial and LPS-neutralizing activities. [BMB reports 2011; 44(11): 747-752].     

Structure of the antimicrobial β-hairpin peptide protegrin-1 in a DLPC lipid bilayer investigated by molecular dynamics simulation

All atom molecular dynamics simulations of the 18-residue β-hairpin antimicrobial peptide protegrin-1 (PG-1, RGGRLCYCRRRFCVCVGR-NH₂) in a fully hydrated dilauroylphosphatidylcholine (DLPC) lipid bilayer have been implemented. The goal of the reported work is to investigate the structure of the peptide in a membrane environment (previously solved only in solution [R.L. Fahrner, T. Dieckmann, S.S.L. Harwig, R.I. Lehrer, D. Eisenberg, J. Feigon, Solution structure of protegrin-1, a broad-spectrum antimicrobial peptide from porcine leukocytes. Chemistry and Biology, 3 (1996) 543-550]), and to delineate specific peptide-membrane interactions which are responsible for the peptide's membrane binding properties. A novel, previously unknown, “kick” shaped conformation of the peptide was detected, where a bend at the C-terminal β-strand of the peptide caused the peptide backbone at residues 16-18 to extend perpendicular to the β-hairpin plane. This bend was driven by a highly persistent hydrogen-bond between the polar peptide side-chain of TYR7 and the unshielded backbone carbonyl oxygen atom of GLY17. The H-bond formation relieves the unfavorable free energy of insertion of polar groups into the hydrophobic membrane core. PG-1 was anchored to the membrane by strong electrostatic binding of the protonated N-terminus of the peptide to the lipid head group phosphate anions. The orientation of the peptide in the membrane, and its influence on bilayer structural and dynamic properties are in excellent agreement with solid state NMR measurements [S. Yamaguchi, T. Hong, A. Waring, R.I. Lehrer, M. Hong, Solid-State NMR Investigations of Peptide-Lipid Interaction and Orientation of a b-Sheet Antimicrobial Peptide, Protegrin, Biochemistry, 41 (2002) 9852-9862]. Importantly, two simulations which started from different initial orientations of the peptide converged to the same final equilibrium orientation of the peptide relative to the bilayer. The kick-shaped conformation was observed only in one of the two simulations.     

Scrutiny of chain-length and N-terminal effects in α-helix folding: a molecular dynamics study on polyalanine peptides

Protein folding remains an unsolved problem as main-chain, side-chain, and solvent interactions remain entangled and have been hard to resolve. Polyalanines are promising models to analyze protein folding initiation and propagation structurally as well as energetically. In the present work, the effect of chain-length and N-terminal residue stereochemistry in polyalanine peptides are investigated for their role in the nucleation of α-helical conformation. The end-protected polyalanine peptides, tetra-alanine, Ac-^LAla₄-NHMe (Ia) and Ac-^DAla-^LAla₃-NHMe (Ib), hexa-alanine, Ac-^LAla₆-NHMe (IIa) and Ac-^DAla-^LAla₅-NHMe (IIb), and octa-alanine, Ac-^LAla₈-NHMe (IIIa) and Ac-^DAla-^LAla₇-NHMe (IIIb), are assessed as chain-length and stereochemical-structure perturbed models. The appreciable variations in the sampling of α-helical conformation, including a sampling of α-helix folds, due to the cooperative effect of chain-length and N-terminal residue stereochemistry have been noted. The electrostatics of α-helical conformation rather than the conformational entropy of the main-chain appear to be decisive in the initiation of α-helix folding. The results of the present work will enhance our understanding on the nucleation of α-helical conformation in short peptides and aid in the design of novel peptides with α-helical structure that can modulate disease-related protein–protein interactions.     

Effect of flavonoids on the destabilization of α-synuclein fibrils and their conversion to amorphous aggregate: A molecular dynamics simulation and experimental study

The second most prevalent neurodegenerative disease, Parkinson's disease (PD), is caused by the accumulation and deposition of fibrillar aggregates of the α-Syn into the Lewy bodies. To create a potent pharmacological candidate to destabilize the preformed α-Syn fibril, it is important to understand the precise molecular mechanism underlying the destabilization of the α-Syn fibril. Through molecular dynamics simulations and experiments, we have examined the molecular mechanisms causing the destabilization and suppression of a newly discovered α-Syn fibril with a Greek-key-like shape and an aggregation prone state (APS) of α-Syn in the presence and absence of various Flvs. According to MD simulation and experimental evidence, morin, quercetin, and myricetin are the Flvs, most capable of destabilizing the fibrils and converting them into amorphous aggregates. Compared to galangin and kaempferol, they have more hydroxyl groups and form more hydrogen bonds with fibrils.The processes by which morin and myricetin prevent new fibril production from APS and destabilize the fibrils are different. According to linear interaction energy analysis, van der Waals interaction predominates with morin, and electrostatic interaction dominates with myricetin. Our MD simulation and experimental findings provide mechanistic insights into how Flvs destabilize α-Syn fibrils and change their morphology, opening the door to developing structure-based drugs for treating Parkinson's disease.     

Inhibitory effect of α-ketoglutaric acid on α-glucosidase: integrating molecular dynamics simulation and inhibition kinetics

Abstract

The inhibition of α-glucosidase is used as a key clinical approach to treat type 2 diabetes mellitus and thus, we assessed the inhibitory effect of α-ketoglutaric acid (AKG) on α-glucosidase with both an enzyme kinetic assay and computational simulations. AKG bound to the active site and interacted with several key residues, including ASP68, PHE157, PHE177, PHE311, ARG312, TYR313, ASN412, ILE434 and ARG439, as detected by protein–ligand docking and molecular dynamics simulations. Subsequently, we confirmed the action of AKG on α-glucosidase as mixed-type inhibition with reversible and rapid binding. The relevant kinetic parameter IC₅₀ was measured (IC₅₀ = 1.738?±?0.041?mM), and the dissociation constant was determined (K_{i Slope} = 0.46?±?0.04?mM). Regarding the relationship between structure and activity, a high AKG concentration induced the slight modulation of the shape of the active site, as monitored by hydrophobic exposure. This tertiary conformational change was linked to AKG inhibition and mostly involved regional changes in the active site. Our study provides insight into the functional role of AKG due to its structural property of a hydroxyphenyl ring that interacts with the active site. We suggest that similar hydroxyphenyl ring-containing compounds targeting key residues in the active site might be potential α-glucosidase inhibitors. Abbreviations AKG alpha-ketoglutaric acid

pNPG 4-nitrophenyl-α-d-glucopyranoside

ANS 1-anilinonaphthalene-8-sulfonate

MD molecular dynamics.

Communicated by Ramaswamy H. Sarma     

Effect of the A30P mutation on the structural dynamics of micelle-bound αSynuclein released in water: a molecular dynamics study

Atomistic molecular dynamics simulation has been used to probe the effect of the A30P mutation on the structural dynamics of micelle-bound, helical αSynuclein when released in an aqueous environment. On the timescales simulated, the effect of the mutation on the secondary structure is restricted to local changes close to the mutation site in the N-terminal helical domain. The changes are transient, and all residues except Lys23 recover their initial structure. The local behavior due to the mutation gives rise to a global difference in the A30P mutant in the form of a permanent kink in the N-terminal helical domain.     

Differential activity of lytic α-helical peptides on lactobacilli and lactobacilli-derived liposomes

Eukaryotic antimicrobial peptides (AMPs) interact with plasma membrane of bacteria, fungi and eukaryotic parasites. Noteworthy, Lactobacillus delbrueckii subsp. lactis (CIDCA 133) and L. delbrueckii subsp. bulgaricus (CIDCA 331) show different susceptibility to human beta-defensins (β-sheet peptides). In the present work we extended the study to α-helical peptides from anuran amphibian (Aurein 1.2, Citropin 1.1 and Maculatin 1.1). We studied the effect on whole bacteria and liposomes formulated with bacterial lipids through growth kinetics, flow cytometry, leakage of liposome content and studies of peptide insertion in lipid monolayers.Growth of strain CIDCA 331 was dramatically inhibited in the presence of all three peptides and minimal inhibitory concentrations were lower than those for strain CIDCA 133. Flow cytometry revealed that AMPs lead to the permeabilization of bacteria.In addition, CIDCA 331-derived liposomes showed high susceptibility, leading to content leakage and structural disruption. Accordingly, peptide insertion in lipid monolayers demonstrated spontaneous interaction of AMPs with CIDCA 331 lipids. In contrast, lipids monolayers from strain CIDCA 133 were less susceptible.Summarizing we demonstrate that the high resistance of the probiotic strain CIDCA 133 to AMPs extends to α helix peptides Aurein, Citropin and Maculatin. This behavior could be ascribed in part to differences in membrane composition. These findings, along with the previously demonstrated resistance to β defensins from human origin, suggest that strain CIDCA 133 is well adapted to host innate immune effectors from both mammals and amphibians thus indicating conserved mechanisms of interaction with key components of the innate immune system.     

Binding of amphipathic α-helical antimicrobial peptides to lipid membranes: Lessons from temporins B and L

Temporins constitute a family of amphipathic α-helical antimicrobial peptides (AMP) and contain some of the shortest cytotoxic peptides, comprised of only 10-14 residues. General characteristics of temporins parallel those of other AMP, both in terms of structural features and biophysical properties relating to their interactions with membrane lipids, with selective lipid-binding properties believed to underlie the discrimination between target vs host cells. Lipid-binding properties also contribute to the cytotoxicity AMP, causing permeabilization of their target cell membranes. The latter functional property of AMP involves highly interdependent acidic phospholipid-induced conformational changes, aggregation, and formation of toxic oligomers in the membrane. These oligomers are subsequently converted to amyloid-type fibers, as demonstrated for e.g. temporins B and L in our laboratory, and more recently for dermaseptins by Auvynet et al. Amyloid state represents the generic minimum in the folding/aggregation free energy landscape, and for AMP its formation most likely serves to detoxify the peptides, in keeping with the current consensus on mature amyloid being inert and non-toxic. The above scenario is supported by sequence analyses of temporins as well as other amphipathic α-helical AMP belonging to diverse families. Accordingly, sequence comparison identifies ‘conformational switches’, domains with equal probabilities for adopting random coil, α-helical and β-sheet structures. These regions were further predicted also to aggregate and assemble into amyloid β-sheets. Taken together, the lipid-binding properties and structural characterization lend support to the notion that the mechanism of membrane permeabilization by temporins B and L and perhaps of most AMP could be very similar, if not identical, to that of the paradigm amyloid forming cytotoxic peptides, responsible for degenerative cell loss in e.g. prion, Alzheimer's and Parkinson's disease, and type 2 diabetes.     

Prediction of the binding mode between GSK3β and a peptide derived from GSKIP using molecular dynamics simulation

GSK3β plays an important role in many physiological functions; dysregulated GSK3β is involved in human diseases such as diabetes, cancer, and Alzheimer's disease. This study uses MD simulations to determine the interaction between GSK3β and a peptide derived from GSKIP, a novel GSK3β interacting protein. Results show that GSKIPtide is inlaid in a binding pocket consisting of an α-helix and an extended loop near the carboxy-terminal end. This binding pocket is hydrophobic, and is responsible for the protein-protein interaction of two other GSK3β interacting proteins: FRAT and Axin. The GSKIPtide binding mode is closer to that of AxinGID (in the Axin-GSK3-interacting domain). The single-point mutations of V267G and Y288F in GSK3β differentiate the binding modes between GSK3 and GSKIPtide, AxinGID, and FRATide. The V2677G mutation of GSK3β reduces the GSKIPtide binding affinity by 70% and abolishes the binding affinity with AxinGID, but has no effect on FRATide. However, GSK3β Y288F completely abolishes the FRATide binding without affecting GSKIPtide or AxinGID binding. An analysis of the GSK3β-GSKIPtide complex structure and the X-ray crystal structures of GSK3β-FRATide and GSK3β-AxinGID complexes suggests that the hydroxyl group of Y288 is crucial to maintaining a hydrogen bond network in GSK3β-FRATide. The hydrophobic side chain of V267 maintains the integrity of helix-helix ridge-groove hydrophobic interaction for GSK3β-GSKIPtide and GSK3β-AxinGID. This study simulates these two mutant systems to provide atomic-level evidence of the aforementioned experimental results and validate the wild-type complex structure prediction.     

A molecular dynamics simulation study of the effect of water–graphene interaction on the properties of confined water

The role of water in determining the structure and stability of biomacromolecules has been well studied. In this work, molecular dynamics simulations have been applied to investigate the effect of surface hydrophobicity on the structure and dynamics of water confined between graphene surfaces. In order to evaluate this effect, we apply various attractive/repulsive water–graphene interaction potentials (hydrophobicity). The properties of confined water are studied by applying a purely repulsive interaction potential between water–graphene (modelled as a repulsive r^?12 potential) and repulsive–attractive forces (modelled as an LJ(12-6) potential). Compared to the case of a purely repulsive graphene–water potential, the inclusion of repulsive–attractive forces leads to formation of sharp peaks for density and the number of hydrogen bonds. Also, it was found that repulsive–attractive graphene–water potential caused slower hydrogen bonds dynamics and restricted the diffusion coefficient of water. Consequently, it was found that hydrogen bond breakage and formation rate with the repulsive r^?12 potential model, will increase compared to the corresponding water confined with the LJ(12-6) potential.     

Effect of collagenase–gelatinase ratio on the mechanical properties of a collagen fibril: a combined Monte Carlo–molecular dynamics study

Loading in cartilage is supported primarily by fibrillar collagen, and damage will impair the function of the tissue, leading to pathologies such as osteoarthritis. Damage is initiated by two types of matrix metalloproteinases, collagenase and gelatinase, that cleave and denature the collagen fibrils in the tissue. Experimental and modeling studies have revealed insights into the individual contributions of these two types of MMPs, as well as the mechanical response of intact fibrils and fibrils that have experienced random surface degradation. However, no research has comprehensively examined the combined influences of collagenases and gelatinases on collagen degradation nor studied the mechanical consequences of biological degradation of collagen fibrils. Such preclinical examinations are required to gain insights into understanding, treating, and preventing degradation-related cartilage pathology. To develop these insights, we use sequential Monte Carlo and molecular dynamics simulations to probe the effect of enzymatic degradation on the structure and mechanics of a single collagen fibril. We find that the mechanical response depends on the ratio of collagenase to gelatinase—not just the amount of lost fibril mass—and we provide a possible mechanism underlying this phenomenon. Overall, by characterizing the combined influences of collagenases and gelatinases on fibril degradation and mechanics at the preclinical research stage, we gain insights that may facilitate the development of targeted interventions to prevent the damage and loss of mechanical integrity that can lead to cartilage pathology.

     

Glycines: role in α-helical membrane protein structures and a potential indicator of native conformation

Among the growing number of membrane protein structures in the Protein Data Bank, there are many transmembrane domains that appear to be native-like; at the same time, there are others that appear to have less than complete native-like character. Hence, there is an increasing need for validation tools that distinguish native-like from non-native-like structures. Membrane mimetics used in protein structural characterizations differ in numerous physicochemical properties from native membranes and provide many opportunities for introducing non-native-like features into membrane protein structures. One possible approach for validating membrane protein structures is based on the use of glycine residues in transmembrane domains. Here, we have reviewed the membrane protein structure database and identified a set of benchmark proteins that appear to be native-like. In these structures, conserved glycine residues rarely face the lipid interstices, and many of them participate in close helix-helix packing. Glycine-based validation allowed the identification of non-native-like features in several membrane proteins and also shows the potential for verifying the native-like character for numerous other membrane protein structures.     

Effect of aluminum on the binding properties of α-chymotrypsin

 The effect of aluminum ions on the binding properties of α-chymotrypsin has been studied. The results show that aluminum does not affect the catalytic rate constant k _cat, but it acts as an enzyme activator favoring the binding of the substrate to the catalytic site (i.e. decreasing K _m). Furthermore, aluminum binding to α-chymotrypsin displays about a threefold decrease in its affinity for the macromolecular inhibitor bovine pancreatic trypsin inhibitor (BPTI). Altogether, the different effect of aluminum on the binding of synthetic substrates (e.g. N-α-benzoyl-l-tyrosine ethyl ester, BTEE) and macromolecular inhibitors (e.g. BPTI) to α-chymotrypsin suggests the occurrence of an aluminum-linked conformational change in the enzyme molecule which brings about a marked structural change at the primary and secondary recognition sites for substrates and inhibitors. The modulative effect exerted by aluminum on the enzyme hydrolytic activity has been investigated also as a function of pH. The ion-linked effect appears to be dependent on the pH in a complex fashion, which suggests that aluminum binding is controlled by the protonation of at least two classes of residues on the enzyme molecule. Received: 5 December 1996 / Accepted: 11 March 1997     

Interaction models of substrate peptides and β-secretase studied by NMR spectroscopy and molecular dynamics simulation

The formation of β-amyloid peptide (Aβ) is initiated from cleavage of amyloid precursor protein (APP) by a family of protease, α-, β-, and γ-secretase. Sub W, a substrate peptide, consists of 10 amino acids, which are adjacent to the β-cleavage site of wild-type APP, and Sub M is Swedish mutant with double mutations on the left side of the β-cleavage site of APP. Sub W is a normal product of the metabolism of APP in the secretary pathway. Sub M is known to increase the efficiency of β-secretase activity, resulting in a more specific binding model compared to Sub W. Three-dimensional structures of Sub W and Sub M were studied by CD and NMR spectroscopy in water solution. On the basis of these structures, interaction models of β-secretase and substrate peptides were determined by molecular dynamics simulation. Four hydrogen bonds and one water-mediated interaction were formed in the docking models. In particular, the hydrogen bonding network of Sub M-BACE formed spread over the broad region of the active site of β-secretase (P5-P3′), and the side chain of P2-Asn formed a hydrogen bond specifically with the side chain of Arg235. These are more favorable to the cleavage of Sub M by β-secretase than Sub W. The two substrate peptides showed different tendency to bind to β-secretase and this information may useful for drug development to treat and prevent Alzheimer’s disease.