There is an active interest in peptides that readily cross cell membranes without the assistance of cell membrane receptors¹. Many of these are referred to as cell-penetrating peptides, which are frequently noted for their potential as drug delivery vectors^1-3. Moreover, there is increasing interest in antimicrobial peptides that operate via non-membrane lytic mechanisms^4,5, particularly those that cross bacterial membranes without causing cell lysis and kill cells by interfering with intracellular processes^6,7. In fact, authors have increasingly pointed out the relationship between cell-penetrating and antimicrobial peptides^1,8. A firm understanding of the process of membrane translocation and the relationship between peptide structure and its ability to translocate requires effective, reproducible assays for translocation. Several groups have proposed methods to measure translocation into large unilamellar lipid vesicles (LUVs)^9-13. LUVs serve as useful models for bacterial and eukaryotic cell membranes and are frequently used in peptide fluorescent studies^14,15. Here, we describe our application of the method first developed by Matsuzaki and co-workers to consider antimicrobial peptides, such as magainin and buforin II^16,17. In addition to providing our protocol for this method, we also present a straightforward approach to data analysis that quantifies translocation ability using this assay. The advantages of this translocation assay compared to others are that it has the potential to provide information about the rate of membrane translocation and does not require the addition of a fluorescent label, which can alter peptide properties¹⁸, to tryptophan-containing peptides. Briefly, translocation ability into lipid vesicles is measured as a function of the Foster Resonance Energy Transfer (FRET) between native tryptophan residues and dansyl phosphatidylethanolamine when proteins are associated with the external LUV membrane (Figure 1). Cell-penetrating peptides are cleaved as they encounter uninhibited trypsin encapsulated with the LUVs, leading to disassociation from the LUV membrane and a drop in FRET signal. The drop in FRET signal observed for a translocating peptide is significantly greater than that observed for the same peptide when the LUVs contain both trypsin and trypsin inhibitor, or when a peptide that does not spontaneously cross lipid membranes is exposed to trypsin-containing LUVs. This change in fluorescence provides a direct quantification of peptide translocation over time.

Biophysical Journal

Previous experimental work has shown that the functional properties of the mechanosensitive channel of large conductance (MscL) are affected by variations in lipid composition. Here, we utilize molecular dynamics simulations of Mycobacterium tuberculosis MscL to investigate such lipid composition effects on a molecular level. In particular, two sets of simulations were performed. In the first, trajectories using lipids with different headgroups (phosphatidylcholine and phosphatidylethanolamine) were compared. Protein-lipid interactions were clearly altered by the headgroup changes, leading to conformational differences in the C-terminal region of M. tuberculosis MscL. In the second set of simulations, lipid tails were gradually shortened, thinning the membrane over a molecular dynamics trajectory. These simulations showed evidence of hydrophobic matching between MscL and the lipid membrane, as previously proposed. For all simulations, protein-lipid interaction energies in the second transmembrane region were correlated to mutagenic data, emphasizing the importance of lipid interactions for proper MscL function.

FEBS Letters

Although molecular dynamics simulations are an important tool for studying membrane systems, relatively few simulations have used anionic lipids. This paper reports the first simulation of a pure phosphatidylglycerol (PG) bilayer. The properties of this equilibrated palmitoyloleoylphosphatidylglycerol membrane agree with experimental observations of PG membranes and with previous simulations of monolayers and mixed bilayers containing PG lipids. These simulations also provide interesting insights into hydrogen bonding interactions in PG membranes. This equilibrated membrane will be a useful starting point for simulations of membrane proteins interacting with PG lipids.

Biophysical Journal

A detailed picture of water and ion properties in small pores is important for understanding the behavior of biological ion channels. Several recent modeling studies have shown that small, hydrophobic pores exclude water and ions even if they are physically large enough to accommodate them, a mechanism called hydrophobic gating. This mechanism has been implicated in the gating of several channels, including the mechanosensitive channel of small conductance (MscS). Although the pore in the crystal structure of MscS is wide and was initially hypothesized to be open, it is lined by hydrophobic residues and may represent a nonconducting state. Molecular dynamics simulations were performed on MscS to determine whether or not the structure can conduct ions. Unlike previous simulations of hydrophobic nanopores, electric fields were applied to this system to model the transmembrane potential, which proved to be important. Although simulations without a potential resulted in a dehydrated, occluded pore, the application of a potential increased the hydration of the pore and resulted in current flow through the channel. The calculated channel conductance was in good agreement with experiment. Therefore, it is likely that the MscS crystal structure is closer to a conducting than a nonconducting state.

Biophysical Journal

The structure of the C-terminal domain of the mechanosensitive channel of large conductance (MscL) has generated significant controversy. As a result, several structures have been proposed for this region: the original crystal structure (1MSL) of the Mycobacterium tuberculosis homolog (Tb), a model of the Escherichia coli homolog, and, most recently, a revised crystal structure of Tb-MscL (2OAR). To understand which of these structures represents a physiological conformation, we measured the impact of mutations to the C-terminal domain on the thermal stability of Tb-MscL using circular dichroism and performed molecular dynamics simulations of the original and the revised crystal structures of Tb-MscL. Our results imply that this region is helical and adopts an alpha-helical bundle conformation similar to that observed in the E. coli MscL model and the revised Tb-MscL crystal structure.

FEBS Letters

Buforin II (BF2) is an antimicrobial peptide that is hypothesized to kill bacteria by entering cells and binding nucleic acids. To further investigate this proposed mechanism, we used computer modeling and experimental measurements to consider the interactions between BF2 and DNA. Computational and experimental results imply that the peptide forms specific interactions with DNA. Moreover, we observe a general correlation between DNA affinity and antimicrobial activity for a series of BF2 variants. Thus, our results support the proposed mechanism for BF2 and provide a useful approach for evaluating the nucleic acid interactions of other antimicrobial peptides.

Proteins

Buforin II is a 21-amino acid polycationic antimicrobial peptide derived from a peptide originally isolated from the stomach tissue of the Asian toad Bufo bufo gargarizans. It is hypothesized to target a wide range of bacteria by translocating into cells without membrane permeabilization and binding to nucleic acids. Previous research found that the structure and membrane interactions of buforin II are related to lipid composition. In this study, we used molecular dynamics (MD) simulations along with lipid vesicle experiments to gain insight into how buforin II interacts differently with phosphatidylcholine (PC), phosphatidylglycerol (PG), and phosphatidylethanolamine (PE) lipids. Fluorescent spectroscopic measurements agreed with the previous assertion that buforin II does not interact with pure PC vesicles. Nonetheless, the reduced entry of the peptide into anionic PG membranes versus neutral PC membranes during simulations correlates with the experimentally observed reduction in BF2 translocation through pure PG membranes. Simulations showing membrane entry into PC also provide insight into how buforin II may initially penetrate cell membranes. Our MD simulations also allowed us to consider how neutral PE lipids affect the peptide differently than PC. In particular, the peptide had a more helical secondary structure in simulations with PE lipids. A change in structure was also apparent in circular dichroism measurements. PE also reduced membrane entry in simulations, which correlates with decreased translocation in the presence of PE observed in previous studies. Together, these results provide molecular-level insight into how lipid composition can affect buforin II structure and function and will be useful in efforts to design peptides with desired antimicrobial and cell-penetrating properties.

The FEBS Journal

Thimet oligopeptidase (EC 3.4.24.15) is a zinc(II) endopeptidase implicated in the processing of numerous physiological peptides. Although its role in selecting and processing peptides is not fully understood, it is believed that flexible loop regions lining the substrate-binding site allow the enzyme to conform to substrates of varying structure. This study describes mutant forms of thimet oligopeptidase in which Gly or Tyr residues in the 599-611 loop region were replaced, individually and in combination, to elucidate the mechanism of substrate selection by this enzyme. Decreases in k(cat) observed on mutation of Tyr605 and Tyr612 demonstrate that these residues contribute to the efficient cleavage of most substrates. Modeling studies showing that a hinge-bend movement brings both Tyr612 and Tyr605 within hydrogen bond distance of the cleaved peptide bond supports this role. Thus, molecular modeling studies support a key role in transition state stabilization of this enzyme by Tyr605. Interestingly, kinetic parameters show that a bradykinin derivative is processed distinctly from the other substrates tested, suggesting that an alternative catalytic mechanism may be employed for this particular substrate. The data demonstrate that neither Tyr605 nor Tyr612 is necessary for the hydrolysis of this substrate. Relative to other substrates, the bradykinin derivative is also unaffected by Gly mutations in the loop. This distinction suggests that the role of Gly residues in the loop is to properly orientate these Tyr residues in order to accommodate varying substrate structures. This also opens up the possibility that certain substrates may be cleaved by an open form of the enzyme.

Peptides

Previous studies have identified several naturally occurring antimicrobial peptides derived from histone proteins. This research aimed to design novel histone-derived antimicrobial peptides (HDAPs). To this end, three novel peptides (DesHDAP1, DesHDAP2, and DesHDAP3) were designed based on a histone-DNA crystal structure and structural properties of buforin II, the best characterized naturally occurring HDAP. Molecular dynamics simulations and circular dichroism spectroscopy were used to further support the predicted structure and potential nucleic acid interactions of these three designed peptides. The antibacterial activity of the three peptides was then verified experimentally against a series of bacterial strains using a radial diffusion assay. One of these peptides is the first known fragment of histone H3 with antibacterial properties. Optical density measurements of bacterial cells exposed to the designed peptides implied that at least two of the novel peptides can induce cell death without causing significant membrane permeabilization, as observed for buforin II. The antibacterial potency of these designed HDAPs does not appear to correlate with their overall alpha-helical content, unlike previous observations for analogs of buforin II. However, the most potent designed peptide, DesHDAP1, shares a markedly similar circular dichroism spectrum with buforin II. These results demonstrate the potential of using histone structures as a framework for designing novel antimicrobial peptides. As well, these studies represent an important starting point for a broader characterization of properties shared by HDAPs.

Biochimica Et Biophysica Acta

Studies of bacterial ion channels have provided significant insights into the structure-function relationships of mechanosensitive and voltage-gated ion channels. However, to date, very few bacterial channels that respond to small molecules have been identified, cloned, and characterized. Here, we use bioinformatics to identify a novel family of bacterial cyclic nucleotide-gated (bCNG) ion channels containing a channel domain related by sequence homology to the mechanosensitive channel of small conductance (MscS). In this initial report, we clone selected members of this channel family, use electrophysiological measurements to verify their ability to directly gate in response to cyclic nucleotides, and use osmotic downshock to demonstrate their lack of mechanosensitivity. In addition to providing insight into bacterial physiology, these channels will provide researchers with a useful model system to investigate the role of ligand-gated ion channels (LGICs) in the signaling processes of higher organisms. The identification of these channels provides a foundation for structural and functional studies of LGICs that would be difficult to perform on mammalian channels. Moreover, the discovery of bCNG channels implies that bacteria have cyclic nucleotide-gated and cyclic nucleotide-modulated ion channels, which are analogous to the ion channels involved in eukaryotic secondary messenger signaling pathways.

Peptides

Buforin II (BF2) is a histone-derived antimicrobial peptide that causes cell death by translocating across membranes and interacting with nucleic acids. It contains one proline residue critical for its function. Previous research found that mutations replacing proline lead to decreased membrane translocation and antimicrobial activity as well as increased membrane permeabilization. This study further investigates the role of proline in BF2's antimicrobial mechanism by considering the effect of changing proline position on membrane translocation, membrane permeabilization, and antimicrobial activity. For this purpose, four mutants were made with proline substitution (P11A) or relocation (P11A/G7P, P11A/V12P, P11A/V15P). These mutations altered the amount of helical content. Although antimicrobial activity correlated with the α-helical content for the peptides containing proline, membrane translocation did not. This observation suggests that factors in BF2's bactericidal mechanism other than translocation must be altered by these mutations. To better explain these trends we also measured the nucleic acid binding and membrane permeabilization of the mutant peptides. A comparison of mutant and wild type BF2 activity revealed that BF2 relies principally on membrane translocation and nucleic acid binding for antimicrobial activity, although membrane permeabilization may play a secondary role for some BF2 variants. A better understanding of the role of proline in the BF2 antimicrobial mechanism will contribute to the further design and development of BF2 analogs. Moreover, since proline residues are prevalent among other antimicrobial peptides, this systematic characterization of BF2 provides general insights that can promote our understanding of other systems.

Future Medicinal Chemistry

Biochemistry and Molecular Biology Education : a Bimonthly Publication of the International Union of Biochemistry and Molecular Biology

As computational modeling plays an increasingly central role in biochemical research, it is important to provide students with exposure to common modeling methods in their undergraduate curriculum. This article describes a series of computer labs designed to introduce undergraduate students to energy minimization, molecular dynamics simulations, and homology modeling. These labs were created as part of a one-semester course on the molecular modeling of biochemical systems. Students who completed these activities felt that they were an effective component of the course, reporting improved comfort with the conceptual background and practical implementation of the computational methods. Although created as a component of a larger course, these activities could be readily adapted for a variety of other educational contexts. As well, all of these labs utilize software that is freely available in an academic environment and can be run on fairly common computer hardware, making them accessible to teaching environments without extensive computational resources.

Biophysical Journal

Mutations that alter the phenotypic behavior of the Escherichia coli mechanosensitive channel of small conductance (MscS) have been identified; however, most of these residues play critical roles in the transition between the closed and open states of the channel and are not directly involved in lipid interactions that transduce the tension response. In this study, we use molecular dynamic simulations to predict critical lipid interacting residues in the closed state of MscS. The physiological role of these residues was then investigated by performing osmotic downshock assays on MscS mutants where the lipid interacting residues were mutated to alanine. These experiments identified seven residues in the first and second transmembrane helices as lipid-sensing residues. The majority of these residues are hydrophobic amino acids located near the extracellular interface of the membrane. All of these residues interact strongly with the lipid bilayer in the closed state of MscS, but do not face the bilayer directly in structures associated with the open and desensitized states of the channel. Thus, the position of these residues relative to the lipid membrane appears related to the ability of the channel to sense tension in its different physiological states.

Biochemical and Biophysical Research Communications

We have recently identified and characterized the bacterial cyclic nucleotide gated (bCNG) subfamily of the larger mechanosensitive channel of small conductance (MscS) superfamily of ion channels. The channel domain of bCNG channels exhibits significant sequence homology to the mechanosensitive subfamily of MscS in the regions that have previously been used as a hallmark for channels that gate in response to mechanical stress. However, we have previously demonstrated that three of these channels are unable to rescue Escherichiacoli from osmotic downshock. Here, we examine an additional nine bCNG homologues and further demonstrate that the full-length bCNG channels are unable to rescue E. coli from hypoosmotic stress. However, limited mechanosensation is restored upon removal of the cyclic nucleotide binding domain. This indicates that the C-terminal domain of the MscS superfamily can drive channel gating and further highlight the ability of a superfamily of ion channels to be gated by multiple stimuli.

Biochimica Et Biophysica Acta

The increase in multidrug resistant bacteria has sparked an interest in the development of novel antibiotics. Antimicrobial peptides that operate by crossing the cell membrane may also have the potential to deliver drugs to intracellular targets. Buforin 2 (BF2) is an antimicrobial peptide that shares sequence identity with a fragment of histone subunit H2A and whose bactericidal mechanism depends on membrane translocation and DNA binding. Previously, novel histone-derived antimicrobial peptides (HDAPs) were designed based on properties of BF2, and DesHDAP1 and DesHDAP3 showed significant antibacterial activity. In this study, their DNA binding, permeabilization, and translocation abilities were assessed independently and compared to antibacterial activity to determine whether they share a mechanism with BF2. To investigate the importance of proline in determining the peptides' mechanisms of action, proline to alanine mutants of the novel peptides were generated. DesHDAP1, which shows significant similarities to BF2 in terms of secondary structure, translocates effectively across lipid vesicle and bacterial membranes, while the DesHDAP1 proline mutant shows reduced translocation abilities and antimicrobial potency. In contrast, both DesHDAP3 and its proline mutant translocate poorly, though the DesHDAP3 proline mutant is more potent. Our findings suggest that a proline hinge can promote membrane translocation in some peptides, but that the extent of its effect on permeabilization depends on the peptide's amphipathic properties. Our results also highlight the different antimicrobial mechanisms exhibited by histone-derived peptides and suggest that histones may serve as a source of novel antimicrobial peptides with varied properties.