Molecular Dynamics Simulation Study of the Interaction of Piscidin 1 with DPPC Bilayers: Structure-Activity Relationship

Abstract

To study structure-activity relationship of antimicrobial peptides and to design novel antimicrobial peptides with selectivity for bacterial cells, we have performed molecular dynamics simulations of the interaction of Piscidin (Pis1) and its two analogues (Pis1-AA and Pis1-PG) with dipalmitoylphosphatidylcholine (DPPC) bilayer through 45 ns. Our results inform us of the detailed location and orientation of the peptide with respect to the bilayer as well as provide about hydrogen-bond-formation patterns and electrostatics interactions. Simulations show that Pis1 and Pis-AA form the most hydrogen bonds and Pis-PG forms the fewest hydrogen bonds with lipid. Thus, Pis1 and Pis-AA should have stronger interactions with the lipid head group when compared to Pis-PG. Experimental studies have shown that Pis1 and Pis1-AA have a high antimicrobial and hemolytic activities, and Pis1-PG has low hemolytic activity while keeps potent antimicrobial activity. Our results complement the previous experimental studies. According to our MD results and previous experimental studies, Pis1 and Pis1-AA are more effective at the zwitterionic bilayer comparing Pis1-PG. These properties of Pis1-PG could be accordance with its low hemolytic activities.     

Consequences of introducing a disulfide bond into an antibacterial and hemolytic peptide.

The effect of introducing a disulfide bridge between the N- and C-terminal ends on the structure and biological activities of the 13-residue linear peptide PKLLKTFLSKWIG(SPFK), which has both antibacterial and hemolytic activity, have been investigated. The terminal amino acids P and G in SPFK were replaced by cysteines to form a disulfide bridge. The linear peptides C(Acm)KLLKTFLSKWIC(Acm) and C(Acm) KLLKTFLSKWIC(Acm)-amide, where Acm is acetamidomethyl group, showed antibacterial activity but did not possess hemolytic activity unlike SPFK. Introduction of an S-S bridge resulted in enhanced hemolytic activity compared with SPFK. The hemolytic activity was particularly pronounced in the cyclic peptide CKLLKTFLSKWIC-amide. Circular dichroism studies indicate that the cyclic peptides tend to adopt distorted helical structures. The cyclic peptides also have a greater affinity for lipid vesicles, which could be the reason for the effective perturbation of the erythrocyte membrane.     

Molecular dynamic simulation studies on the effect of one residue chain staggering on the structure and stability of heterotrimeric collagen-like peptides with interruption

A systematic molecular dynamics (MD) simulation has been carried out on collagen-like peptides with different combinations of interruptions in the Gly-X(AA) -Y(AA) repeats. Although experimental studies have been carried out to elucidate the structural consequences of homotrimeric collagen-like peptides, this is the first report on the structural effect on the heterotrimeric models with G4G and G1G breaks present simultaneously in the constituent chains with difference in one residue chain staggering. The results reveal that the axial registry of the interrupted region changes significantly from that of conventional triple helical peptide without interruption. Further, results from MD simulations show the formation of a kink in the interrupted region of the triple-helical peptides. The conformational analysis reveals that the interruption in the Gly-X(AA) -Y(AA) pattern in these peptides induces β-strand conformation in triple helical peptides. The conventional hydrogen bonds in the interrupted triad are affected and new nonconventional H-bonds are formed in the triple helical structure, and as a result interrupted region becomes locally fragile. MM-PBSA calculations on the different systems clearly suggest that the binding affinity varies marginally due to one residue staggering. However, it is found from the structural parameters that hydrogen-bonding pattern differs significantly due to the difference in the staggering of chains.     

Insights into unbinding mechanisms upon two mutations investigated by molecular dynamics study of GSK3beta-axin complex: role of packing hydrophobic residues

Glycogen synthase kinase 3beta (GSK 3beta) is a key component of several cellular processes including Wnt and insulin signalling pathways. The interaction of GSK3beta with scaffolding peptide axin is thought to be responsible for the effective phosphorylation of beta-catenin, the core effector of Wnt signaling, which has been linked with the occurrence of colon cancer and melanoma. It has been demonstrated that the binding of axin to GSK3beta is abolished by the single-point mutation of Val267 to Gly (V267G) in GSK3beta or Leu392 to Pro (L392P) in axin. Molecular dynamics (MD) simulations were performed on wild type (WT), V267G mutant and L392P one to elucidate the two unbinding mechanisms that occur through different pathways. Besides, rough energy and residue-based energy decomposition were calculated by MM_GBSA (molecular mechanical Generalized_Born surface area) approach to illuminate the instability of the two mutants. The MD simulations of the two mutants and WT reveal that the structure of GSK3beta remains unchanged, while axin moves away from the interfacial hydrophobic pockets in both two mutants. Axin exhibits positional shift in V267G mutant, whereas, losing the hydrogen bonds that are indispensable for stabilizing the helix structure of wild type axin, the helix of axin is distorted in L392P mutant. To conclude, both two mutants destroy the hydrophobic interaction that is essential to the stability of GSK3beta-axin complex.     

Terminal ion pairs stabilize the second beta-hairpin of the B1 domain of protein G

The effects of terminal ion pairs on the stability of a beta-hairpin peptide corresponding to the C-terminal residues of the B1 domain of protein G were determined using thermal unfolding as monitored by nuclear magnetic resonance and circular dichroism spectroscopy. Molecular dynamics (MD) simulations were also performed to examine the effect of ion pairs on the structures. Eight peptides were studied including the wild type (G41) and the N-terminal modified sequences that had the first residue deleted (E42), replaced with a Lys (K41), or extended by an additional Gly (G40). Acetylated variants were made to examine the effect of removing the positive N-terminal charge on beta-hairpin stability. The rank in stability determined experimentally is K41 > E42 approximately G41 approximately G40 > Ac-K41 > Ac-E42 approximately Ac-G41 > Ac-G40. The Tm of the K41 peptide is 12 degrees C higher than G41, while the Tm values for the acetylated peptides are less than their unacetylated forms by more than 15 degrees C. NOE cross-peaks between side-chain methylene groups at the N- and C-termini and larger CalphaH shifts compared to random values are seen for K41. The addition of 20% methanol increases the stability in K41 and G41. The MD studies complement these results by showing that the charged N-terminus is important to stability. The type of ion pair observed varies with peptide, and when formed the simulations show that the ion pair can prevent fraying of the beta-strands through electrostatic and hydrophobic contacts. Therefore, introducing favorable electrostatic interactions at the N- and C-termini can substantially enhance beta-hairpin stability and help define the structure.     

Design and activity of novel lactoferrampin analogues against O157:H7 enterohemorrhagic escherichia coli

Lactoferrampin 265–284 (LFampin 265–284) is a peptide consisting of residues 265–284 of N1‐domain of bovine Lactoferrin (LF). This peptide has several cationic groups in the C‐terminal lobe, exhibiting an antibacterial activity against a wide range of microorganisms. However, LFampin 265–284 exhibits low antimicrobial activity against the O157:H7 enterohaemorrhagic Escherichia coli (EHEC O157:H7) when compared with Lactoferrin chimera and Lactoferricin. Here, we have designed three analogues of LFampin 265–284 based on the distribution of cationic groups, hydrophobicity, size, and sequence. Analogues were synthesized by solid phase chemistry using Fmoc methodology obtaining peptides with 95% purity. All peptides maintain the ability to adopt helical conformations (checked by circular dichroism spectra and molecular simulations). Some of these analogues exhibited a significant increase in antimicrobial activity by counting colony forming units against EHEC O157:H7 compared to native LFampin 265–284, with MIC of 10 and 40 µM for 264G‐D265K and 264G‐D265K/S272R, respectively. The incorporation of a GKLI sequence in the N‐terminal lobe increased dramatically its antibacterial activity, an effect which has been attributed to the addition of cationic groups in the N‐terminal side that may stabilize the helical conformation of the new designed peptides. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 319–328, 2014.     

Folding thermodynamics of β‐hairpins studied by replica‐exchange molecular dynamics simulations

We study the differences in folding stability of β‐hairpin peptides, including GB1 hairpin and a point mutant GB1 K10G, as well as tryptophan zippers (TrpZips): TrpZip1, TrpZip2, TrpZip3‐1, and TrpZip4. By performing replica‐exchange molecular dynamics simulations with Amber03* force field (a modified version of Amber ff03) in explicit solvent, we observe ab initio folding of all the peptides except TrpZip3‐1, which is experimentally known to be the least stable among the peptides studied here. By calculating the free energies of unfolding of the peptides at room temperature and folding midpoint temperatures for thermal unfolding of peptides, we find that TrpZip4 and GB1 K10G peptides are the most stable β‐hairpins followed by TrpZip1, GB1, and TrpZip2 in the given order. Hence, the proposed K10G mutation of GB1 peptide results in enhanced stability compared to wild‐type GB1. An important goal of our study is to test whether simulations with Amber 03* model can reproduce experimentally predicted folding stability differences between these peptides. While the stabilities of GB1 and TrpZip1 yield close agreement with experiment, TrpZip2 is found to be less stable than predicted by experiment. However, as heterogenous folding of TrpZip2 may yield divergent thermodynamic parameters by different spectroscopic methods, mismatching of results with previous experimental values are not conclusive of model shortcomings. For most of the cases, molecular simulations with Amber03* can successfully reproduce experimentally known differences between the mutated peptides, further highlighting the predictive capabilities of current state‐of‐the‐art all‐atom protein force fields. Proteins 2015; 83:1307–1315. © 2015 Wiley Periodicals, Inc.     

Role of amino acid residues within the disulfide loop of thanatin,a potent antibiotic peptide

Thanatin, a 21-residue peptide, is an inducible insect peptide with a broad range of activity against bacteria and fungi. It has a C-terminal disulfide loop, like the frog skin secretion antimicrobial peptides of the brevinin family. In this study, we tried to find the effect of a number of amino acids between the disulfide bond. Thanatin showed stronger antibacterial activity to Gram negative bacteria than other mutants, except Th1; whereas, the mutant peptides with deletion had higher activity to Gram positive bacteria than thanatin. An increase in the number of amino acid(s) using the alanine residue decreased the antibacterial activity in all of the bacteria. Th1 with deletion of threonine at position 15 (Thr(1)(2)) showed similar antibacterial activity against Gram-negative bacteria, but had higher activity against the Gram positive bacteria. In order to study the structure-function relationship, we measured liposome disruption by the peptides and CD spectra of the peptides. Th1 also showed the highest liposome leaking activity and alpha-helical propensity in the sodium dodecyl sulfate solution, compared with other peptides. Liposome disruption activity was closely correlated with the anti-Gram positive bacterial activity. All of the peptides showed no hemolytic activity. Th1 was considered to be useful as an antimicrobial peptide with broad spectrum without toxicity     

Molecular dynamic studies on the impact of mutations on the structure,stability, and N‐terminal orientation of annexin A1: Implications for membrane aggregation

Multiple MD simulations were performed for the full‐length wild‐type A1, the full length A1 mutations S27E and S27A, as well as the N‐terminal peptide (AMVSEFLKQAWFIDNEEQEYIKTVKG S ²⁷ KGGPGSAVSPYPTFN) of wild‐type A1 and mutations S27E and S27A. The MD simulation trajectories of about 350 ns were generated and analyzed to examine the changes of core domain calcium binding affinity, core domain and N‐terminal domain structures, and N‐terminal domain orientation. Our results indicated that S27A and S27E mutations caused little changes on the calcium‐binding affinity of the core domain of A1. However, the S27A mutation made the N‐terminal domain of A1 less helical, and made the N‐terminal domain migrate faster toward the core domain; these impacts on A1 are beneficial to the membrane aggregation process. On the contrary, the S27E mutation made the N‐terminal domain of A1 more stable, and hindered the migration to the core domain; these changes on A1 are antagonistic for the membrane aggregation process. Our results using MD simulations provide an atomistic explanation for experimental observations that the S27E mutant showed a higher calcium concentration requirement and lower maximal extent of aggregation, while the wild‐type and two mutants S27E and S27A required identical calcium concentrations for liposome binding. Proteins 2014; 82:3327–3334. © 2014 Wiley Periodicals, Inc.     

Streptococcal surface proteins activate the contact system and control its antibacterial activity

Group G streptococci (GGS) are important bacterial pathogens in humans. Here, we investigated the interactions between GGS and the contact system, a procoagulant and proinflammatory proteolytic cascade that, upon activation, also generates antibacterial peptides. Two surface proteins of GGS, protein FOG and protein G (PG), were found to bind contact system proteins. Experiments utilizing contact protein-deficient human plasma and isogenic GGS mutant strains lacking FOG or PG showed that FOG and PG both activate the procoagulant branch of the contact system. In contrast, only FOG induced cleavage of high molecular weight kininogen, generating the proinflammatory bradykinin peptide and additional high molecular weight kininogen fragments containing the antimicrobial peptide NAT-26. On the other hand, PG protected the bacteria against the antibacterial effect of NAT-26. These findings underline the significance of the contact system in innate immunity and demonstrate that GGS have evolved surface proteins to exploit and modulate its effects.     

Solution structure and cell selectivity of piscidin 1 and its analogues

Piscidin 1 (Pis-1) is a novel cytotoxic peptide with a cationic alpha-helical structure that was isolated from the mast cells of hybrid striped bass [Silphaduang, U., and Noga, E. J. (2001) Nature 414, 268-269]. Pis-1 is not selective for bacterial versus mammalian cells. In the present study, to develop novel antibiotic peptides with selectivity for bacterial cells, we examined the effect of substituting two glycine residues, Gly8 and Gly13, with Ala or Pro on this peptide's structure and biological activities. The bacterial cell selectivity of the peptides decreased in the following order: Gly-->Pro analogues > Gly-->Pro/Ala analogues > Pis-1 > Gly-->Ala analogues. The antimicrobial and hemolytic activities and abilities to permeabilize the model phospholipid membranes were higher for Pis-1 with Gly or Pro at position 8 than for its counterparts with either Gly or Pro at position 13. We determined the tertiary structure of Pis-1 and its analogues in the presence of SDS micelles by NMR spectroscopy. We found that Pis-1 has an alpha-helical structure from Phe2 to Thr21. Also, Pis-1 AA (Gly8, Gly13-->Ala8, Ala13) with higher antibacterial and hemolytic activity than Pis-1 has a stable alpha-helical structure from Phe2 to Thr21. Pis-1 PG (Gly-->Pro8) with bacterial cell selectivity has a hinge structure at Pro8, which provides flexibility in piscidin, followed by a three-turn helix from Val10 to Gly22 in the C-terminal region. Taken together, our results demonstrate that the conformational flexibility provided by introduction of a Pro at position 8, coupled with the primary anchoring of phenylalanines and histidines in the N-terminus to the cell membrane and the optimal length of the C-terminal amphipathic alpha-helix, are the critical factors that confer antibacterial activity and bacterial cell selectivity to Pis-1 PG. Pis-1 PG may be a good candidate for the development of a new drug with potent antibacterial activity but without cytotoxicity.     

Lipid-induced conformational transition of the amyloid core fragment Abeta(28-35) and its A30G and A30I mutants

The interaction of the beta-amyloid peptide (Abeta) with neuronal membranes could play a key role in the pathogenesis of Alzheimer's disease. Recent studies have focused on the interactions of Abeta oligomers to explain the neuronal toxicity accompanying Alzheimer's disease. In our study, we have investigated the role of lipid interactions with soluble Abeta(28-35) (wild-type) and its mutants A30G and A30I in their aggregation and conformational preferences. CD and Trp fluorescence spectroscopic studies indicated that, immediately on dissolution, these peptides adopted a random coil structure. Upon addition of negatively charged 1,2-dipalmitoyl-syn-glycero-3-phospho-rac-(glycerol) sodium salt (PG) lipid, the wild-type and A30I mutant underwent reorganization into a predominant beta-sheet structure. However, no conformational changes were observed in the A30G mutant on interaction with PG. In contrast, the presence of zwitterionic 1,2-dipalmitoyl-syn-glycero-3-phosphatidylcholine (PC) lipid had no effect on the conformation of these three peptides. These observations were also confirmed with atomic force microscopy and the thioflavin-T assay. In the presence of PG vesicles, both the wild-type and A30I mutant formed fibrillar structures within 2 days of incubation in NaCl/P(i), but not in their absence. Again, no oligomerization was observed with PC vesicles. The Trp studies also revealed that both ends of the three peptides are not buried deep in the vesicle membrane. Furthermore, fluorescence spectroscopy using the environment-sensitive probe 1,6-diphenyl-1,3,5-hexatriene showed an increase in the membrane fluidity upon exposure of the vesicles to the peptides. The latter effect may result from the lipid head group interactions with the peptides. Fluorescence resonance energy transfer experiments revealed that these peptides undergo a random coil-to-sheet conversion in solution on aging and that this process is accelerated by negatively charged lipid vesicles. These results indicate that aggregation depends on hydrophobicity and propensity to form beta-sheets of the amyloid peptide, and thus offer new insights into the mechanism of amyloid neurodegenerative disease.     

Cecropin D‐like antibacterial peptides from the sphingid moth,Agrius convolvuli

Two major antibacterial peptides were isolated and purified from immunized larval hemolymph of Agrius convolvuli. Acid extraction, gel filtration, ultrafiltration, and reversed‐phase FPLC were used for purification of peptides. These peptides had similar molecular mass and amino acid composition. Moreover, 21 of the first 23 N terminal residues were identical. The peptides were highly homologous with cecropin D in size and primary sequence, and named Agrius cecropin D1 and D2. The molecular masses of Agrius cecropin D1 and D2 were 3,879.39 and 3,839.27, respectively. In antibacterial and hemolytic assays, Agrius cecropin D showed potent antibacterial activities against a panel of Gram positive and negative bacteria without hemolytic activity against human red blood cells. Notably, our antibacterial assay revealed Agrius cecropin D possessed stronger or at least equivalent activities against B. megaterium than cecropin A. It suggests that Agrius cecropin D, which has an alternative structure from cecropin D, could be the model for the development of peptide antibiotics. Arch. Insect Biochem. Physiol. 41:178–185, 1999. © 1999 Wiley‐Liss, Inc.     

Molecular characterization of the β‐amyloid(4‐10) epitope of plaque specific Aβ antibodies by affinity‐mass spectrometry using alanine site mutation

Alzheimer disease is a neurodegenerative disease affecting an increasing number of patients worldwide. Current therapeutic strategies are directed to molecules capable to block the aggregation of the β‐amyloid(1‐42) (Aβ) peptide and its shorter naturally occurring peptide fragments into toxic oligomers and amyloid fibrils. Aβ‐specific antibodies have been recently developed as powerful antiaggregation tools. The identification and functional characterization of the epitope structures of Aβ antibodies contributes to the elucidation of their mechanism of action in the human organism. In previous studies, the Aβ(4‐10) peptide has been identified as an epitope for the polyclonal anti‐Aβ(1‐42) antibody that has been shown capable to reduce amyloid deposition in a transgenic Alzheimer disease mouse model. To determine the functional significance of the amino acid residues involved in binding to the antibody, we report here the effects of alanine single‐site mutations within the Aβ‐epitope sequence on the antigen‐antibody interaction. Specific identification of the essential affinity preserving mutant peptides was obtained by exposing a Sepharose‐immobilized antibody column to an equimolar mixture of mutant peptides, followed by analysis of bound peptides using high‐resolution MALDI‐Fourier transform‐Ion Cyclotron Resonance mass spectrometry. For the polyclonal antibody, affinity was preserved in the H6A, D7A, S8A, and G9A mutants but was lost in the F4, R5, and Y10 mutants, indicating these residues as essential amino acids for binding. Enzyme‐linked immunosorbent assays confirmed the binding differences of the mutant peptides to the polyclonal antibody. In contrast, the mass spectrometric analysis of the mutant Aβ(4‐10) peptides upon affinity binding to a monoclonal anti‐Aβ(1‐17) antibody showed complete loss of binding by Ala‐site mutation of any residue of the Aβ(4‐10) epitope. Surface plasmon resonance affinity determination of wild‐type Aβ(1‐17) to the monoclonal Aβ antibody provided a binding constant K_D in the low nanomolar range. These results provide valuable information in the elucidation of the binding mechanism and the development of Aβ‐specific antibodies with improved therapeutic efficacy.     

Molecular dynamics investigation of the influence of anionic and zwitterionic interfaces on antimicrobial peptides' structure: implications for peptide toxicity and activity

Molecular dynamics simulations of three related helical antimicrobial peptides have been carried out in zwitterionic diphosphocholine (DPC) micelles and anionic sodiumdodecylsulfate (SDS) micelles. These systems can be considered as model mammalian and bacterial membrane interfaces, respectively. The goal of this study is to dissect the differences in peptide composition which make the mutant peptides (novispirin-G10 and novispirin-T7) less toxic than the parent peptide ovispirin (OVIS), although all three peptides have highly antibacterial properties. Compared to G10 and T7, OVIS inserts deepest into the DPC micelle. This correlates well with the lesser toxicity of G10 and T7. There is strong evidence which suggests that synergistic binding of hydrophobic residues drives binding of OVIS to the micelle. The helical content of G10 and T7 is reduced in the presence of DPC, and this leads to less amphipathic peptide structures, which bind weakly to the micelle. Simulations in SDS were carried out to compare the influence of membrane electrostatics on peptide structure. All three peptides bound strongly to SDS, and retained helical form. This corresponds well with their equally potent antibacterial properties. Based on the simulations, we argue that secondary structure stability often leads to toxic properties. We also propose that G10 and T7 operate by the carpet mechanism of cell lysis. Toxicity of peptides operating by the carpet mechanism can be attenuated by reducing the peptide helical content. The simulations successfully capture experimental binding states, and the different depths of binding of the three peptides to the two micelles correlate with their antibacterial and toxic properties.     

Conolysin-Mt: a conus peptide that disrupts cellular membranes

Conus venoms are estimated to comprise over 100,000 distinct pharmacologically active peptides, the majority probably targeting ion channels. Through the characterization of a cytolytic peptide from the venom of Conus mustelinus, conolysin-Mt, we expand the known conopeptide mechanisms to include association with and destruction of cellular membranes. A new 23AA conopeptide, conolysin-Mt has potent hemolytic activity when tested on human erythrocytes. At a concentration of 0.25 microM, the peptide permeabilized both negatively charged prokaryotic (PE:PG) and zwitterionic eukaryotic (PC:cholesterol) model membranes. The affinity constants (KA) of conolysin-Mt for PE:PG and PC:cholesterol model membranes were 0.9 +/- 0.3 x 10(7) and 3 +/- 1 x 10(7) M-1, respectively. In contrast, conolysin-Mt exhibited low antimicrobial activity (MIC > 50 microM) against two Escherichia coli strains, with an MIC for the Gram-positive S. aureus of 25-50 microM. The specificity of conolysin-Mt for native eukaryotic membranes is a novel feature of the peptide compared to other well-characterized cytolytic peptides such as melittin.     

Structure-activity relationship of an antibacterial peptide, maculatin 1.1, from the skin glands of the tree frog, Litoria genimaculata.

Maculatin 1.1 (Mac) is a cationic antibacterial peptide isolated from the dorsal glands of the tree frog, Litoria genimaculata, and has a sequence of GLFGVLAKVAAHVVPAIAEHF-NH2. A short peptide lacking the N-terminal two residues of Mac was reported to have no activity. To investigate the structure-activity relationship in detail, several analogs and related short peptides of Mac were synthesized. CD measurement showed that all the peptides took more or less an alpha-helical structure in the presence of anionic lipid vesicles. Analogs which are more basic than Mac had strong antibacterial and hemolytic activities, while short peptides lacking one or two terminal residues exhibited weak or no activity. Outer and inner membrane permeabilization activities of the peptides were also reduced with shortening of the peptide chain. These results indicate that the entire chain length of Mac is necessary for full activity, and the basicity of the peptides greatly affects the activity.     

Line-tension controlled mechanism for influenza fusion

Our molecular simulations reveal that wild-type influenza fusion peptides are able to stabilize a highly fusogenic pre-fusion structure, i.e. a peptide bundle formed by four or more trans-membrane arranged fusion peptides. We rationalize that the lipid rim around such bundle has a non-vanishing rim energy (line-tension), which is essential to (i) stabilize the initial contact point between the fusing bilayers, i.e. the stalk, and (ii) drive its subsequent evolution. Such line-tension controlled fusion event does not proceed along the hypothesized standard stalk-hemifusion pathway. In modeled influenza fusion, single point mutations in the influenza fusion peptide either completely inhibit fusion (mutants G1V and W14A) or, intriguingly, specifically arrest fusion at a hemifusion state (mutant G1S). Our simulations demonstrate that, within a line-tension controlled fusion mechanism, these known point mutations either completely inhibit fusion by impairing the peptide's ability to stabilize the required peptide bundle (G1V and W14A) or stabilize a persistent bundle that leads to a kinetically trapped hemifusion state (G1S). In addition, our results further suggest that the recently discovered leaky fusion mutant G13A, which is known to facilitate a pronounced leakage of the target membrane prior to lipid mixing, reduces the membrane integrity by forming a 'super' bundle. Our simulations offer a new interpretation for a number of experimentally observed features of the fusion reaction mediated by the prototypical fusion protein, influenza hemagglutinin, and might bring new insights into mechanisms of other viral fusion reactions.     

Design and synthesis of cyclic disulfide-bonded antibacterial peptides on the basis of the alpha helical domain of Tenecin 1, an insect defensin

We synthesized cyclic disulfide-bonded (i, i+4) peptides with various net positive charges (+2-+5) from linear peptides derived from the alpha helical domain of Tenecin 1, an insect defensin, and investigated the effect of the intradisulfide bridge (i, i+4) on hydrophobicity, secondary structure, leakage activity and binding activity for large unilamellar vesicles, antimicrobial activity, and hemolytic activity. Intradisulfide bridge formation of the peptides resulted in the increase of amphiphilicity and hydrophobicity. Cyclic forms of the peptides did not deeply penetrate into PG/PC (1:1, mole ratio) large unilamellar vesicles and had a decreased lipid membrane perturbation activity for PG/PC LUVs. When the peptides interacted with PG/CL (2:1, mole ratio) LUVs, cyclic peptides with a high net positive charge (+4-+5) showed similar binding affinities and leakage activities for vesicles to those of linear forms, whereas cyclic peptides with a low net positive charge (+2-+3) exhibited lower leakage activity than their linear forms. CD spectra indicate that the intradisulfide bridge (i, i+4) provided little conformational constraint to linear peptides in buffer solution but resulted in the decrease of alpha helicity of the peptides in lipid membrane mimic conditions. The cyclic peptide with the highest net positive charge had a similar antibacterial activity to that of the linear peptide, whereas the cyclic peptides with a low net positive charge (+3-+4) exhibited lower antibacterial activity than their linear forms. The cyclic peptides of an appropriate net charge showed more potent activities against some bacteria than those of linear forms under high salt conditions.     

Prediction of binding affinities between the human amphiphysin-1 SH3 domain and its peptide ligands using homology modeling, molecular dynamics and molecular field analysis

The SH3 domain of the human protein amphiphysin-1, which plays important roles in clathrin-mediated endocytosis, actin function and signaling transduction, can recognize peptide motif PXRPXR (X is any amino acid) with high affinity and specificity. We have constructed a complex structure of the amphiphysin-1 SH3 domain and a high-affinity peptide ligand PLPRRPPRA using homology modeling and molecular docking, which was optimized by molecular dynamics (MD). Three-dimensional quantitative structure-affinity relationship (3D-QSAR) analyses on the 200 peptides with known binding affinities to the amphiphysin-1 SH3 domain was then performed using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMSIA model showed promising predictive power, giving good predictions for about 95% of the peptides in the test set (absolute prediction errors less than 1.0). It was used to validate peptide-SH3 binding structure and provide insight into the structural requirements for binding of peptides to SH3 domains. Finally, MD simulations were performed to analyze the interaction between the SH3 domain and another peptide GFPRRPPPRG that contains with the PXRPXsR (s represents residues with small side chains) motif. MD simulations demonstrated that the binding conformation of GFPRRPPPRG is quite different from that of PLPRRPPRAA especially the four residues at the C terminal, which may explain why the CoMSIA model cannot give good predictions on the peptides of the PXRPXsR motif. Because of its efficiency and predictive power, the 3D-QSAR model can be used as a scoring filter for predicting peptide sequences bound to SH3 domains.