Understanding the energetics of helical peptide orientation in membranes

A Critical Review of Short Antimicrobial Peptides from Scorpion Venoms, Their Physicochemical Attributes, and Potential for the Development of New Drugs

Scorpion venoms have proven to be excellent sources of antimicrobial agents. However, although many of them have been functionally characterized, they remain underutilized as pharmacological agents, despite their evident therapeutic potential. In this review, we discuss the physicochemical properties of short scorpion venom antimicrobial peptides (ssAMPs). Being generally short (13-25 aa) and amidated, their proven antimicrobial activity is generally explained by parameters such as their net charge, the hydrophobic moment, or the degree of helicity. However, for a complete understanding of their biological activities, also considering the properties of the target membranes is of great relevance. Here, with an extensive analysis of the physicochemical, structural, and thermodynamic parameters associated with these biomolecules, we propose a theoretical framework for the rational design of new antimicrobial drugs. Through a comparison of these physicochemical properties with the bioactivity of ssAMPs in pathogenic bacteria such as Staphylococcus aureus or Acinetobacter baumannii, it is evident that in addition to the net charge, the hydrophobic moment, electrostatic energy, or intrinsic flexibility are determining parameters to understand their performance. Although the correlation between these parameters is very complex, the consensus of our analysis suggests that there is a delicate balance between them and that modifying one affects the rest. Understanding the contribution of lipid composition to their bioactivities is also underestimated, which suggests that for each peptide, there is a physiological context to consider for the rational design of new drugs.

Implicit model to capture electrostatic features of membrane environment

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Implicit model to capture electrostatic features of membrane environment

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

Modeling of the Passive Permeation of Mercury and Methylmercury Complexes Through a Bacterial Cytoplasmic Membrane

Cellular uptake and export are important steps in the biotransformation of mercury (Hg) by microorganisms. However, the mechanisms of transport across biological membranes remain unclear. Membrane-bound transporters are known to be relevant, but passive permeation may also be involved. Inorganic Hg^II and methylmercury ([CH₃Hg^II]⁺) are commonly complexed with thiolate ligands. Here, we have performed extensive molecular dynamics simulations of the passive permeation of Hg^II and [CH₃Hg^II]⁺ complexes with thiolate ligands through a model bacterial cytoplasmic membrane. We find that the differences in free energy between the individual complexes in bulk water and at their most favorable position within the membrane are ∼2 kcal mol^-1. We provide a detailed description of the molecular interactions that drive the membrane crossing process. Favorable interactions with carbonyl and tail groups of phospholipids stabilize Hg-containing solutes in the tail-head interface region of the membrane. The calculated permeability coefficients for the neutral compounds CH₃S-Hg^II-SCH₃ and CH₃Hg^II-SCH₃ are on the order of 10^-5 cm s^-1. We conclude that small, nonionized Hg-containing species can permeate readily through cytoplasmic membranes.

Thermodynamic and kinetic characterization of transmembrane helix association

The transient dimerization of transmembrane proteins is an important event in several cellular processes and here we use coarse-grain and meso-scale modeling methods to quantify their underlying dynamics.

Mechanical properties of lipid bilayers and regulation of mechanosensitive function: from biological to biomimetic channels

Material properties of lipid bilayers, including thickness, intrinsic curvature and compressibility regulate the function of mechanosensitive (MS) channels. This regulation is dependent on phospholipid composition, lateral packing and organization within the membrane. Therefore, a more complete framework to understand the functioning of MS channels requires insights into bilayer structure, thermodynamics and phospholipid structure, as well as lipid-protein interactions. Phospholipids and MS channels interact with each other mainly through electrostatic forces and hydrophobic matching, which are also crucial for antimicrobial peptides. They are excellent models for studying the formation and stabilization of membrane pores. Importantly, they perform equivalent responses as MS channels: (1) tilting in response to tension and (2) dissipation of osmotic gradients. Lessons learned from pore forming peptides could enrich our knowledge of mechanisms of action and evolution of these channels. Here, the current state of the art is presented and general principles of membrane regulation of mechanosensitive function are discussed.

Monitoring membrane binding and insertion of peptides by two-color fluorescent label

Herein, we developed an approach for monitoring membrane binding and insertion of peptides using a fluorescent environment-sensitive label of the 3-hydroxyflavone family. For this purpose, we labeled the N-terminus of three synthetic peptides, melittin, magainin 2 and poly-l-lysine capable to interact with lipid membranes. Binding of these peptides to lipid vesicles induced a strong fluorescence increase, which enabled to quantify the peptide-membrane interaction. Moreover, the dual emission of the label in these peptides correlated well with the depth of its insertion measured by the parallax quenching method. Thus, in melittin and magainin 2, which show deep insertion of their N-terminus, the label presented a dual emission corresponding to a low polar environment, while the environment of the poly-l-lysine N-terminus was rather polar, consistent with its location close to the bilayer surface. Using spectral deconvolution to distinguish the non-hydrated label species from the hydrated ones and two photon fluorescence microscopy to determine the probe orientation in giant vesicles, we found that the non-hydrated species were vertically oriented in the bilayer and constituted the best indicators for evaluating the depth of the peptide N-terminus in membranes. Thus, this label constitutes an interesting new tool for monitoring membrane binding and insertion of peptides.

Coarse grained molecular dynamics simulations of transmembrane protein-lipid systems

Many biological cellular processes occur at the micro- or millisecond time scale. With traditional all-atom molecular modeling techniques it is difficult to investigate the dynamics of long time scales or large systems, such as protein aggregation or activation. Coarse graining (CG) can be used to reduce the number of degrees of freedom in such a system, and reduce the computational complexity. In this paper the first version of a coarse grained model for transmembrane proteins is presented. This model differs from other coarse grained protein models due to the introduction of a novel angle potential as well as a hydrogen bonding potential. These new potentials are used to stabilize the backbone. The model has been validated by investigating the adaptation of the hydrophobic mismatch induced by the insertion of WALP-peptides into a lipid membrane, showing that the first step in the adaptation is an increase in the membrane thickness, followed by a tilting of the peptide.

Peptide nanopores and lipid bilayers: interactions by coarse-grained molecular-dynamics simulations

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (alpha-helical and beta-barrel), and the seven different bilayer systems range in thickness from approximately 28 to approximately 43 A. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.

Determining the orientation of protegrin-1 in DLPC bilayers using an implicit solvent-membrane model

Continuum models that describe the effects of solvent and biological membrane molecules on the structure and behavior of antimicrobial peptides, holds a promise to improve our understanding of the mechanisms of antimicrobial action of these peptides. In such methods, a lipid bilayer model membrane is implicitly represented by multiple layers of relatively low dielectric constant embedded in a high dielectric aqueous solvent, while an antimicrobial peptide is accounted for by a dielectric cavity with fixed partial charge at the center of each one of its atoms. In the present work, we investigate the ability of continuum approaches to predict the most probable orientation of the β-hairpin antimicrobial peptide Protegrin-1 (PG-1) in DLPC lipid bilayers by calculating the difference in the transfer free energy from an aqueous environment to a membrane-water environment for multiple orientations. The transfer free energy is computed as a sum of two terms; polar/electrostatic and non-polar. They both include energetic and entropic contributions to the free energy. We numerically solve the Poisson-Boltzmann equation to calculate the electrostatic contribution to the transfer free energy, while the non-polar contribution to the free energy is approximated using a linear solvent accessible surface area relationships. The most probable orientation of PG-1 is that with the lowest relative transfer free energy. Our simulation results indicate that PG-1 assumes an oblique orientation in DLPC lipid bilayers. The predicted most favorable orientation was with a tilt angle of 19°, which is in qualitative agreement with the experimentally observed orientations derived from solid-state NMR data.

Packing density of the erythropoietin receptor transmembrane domain correlates with amplification of biological responses

The formation of signal-promoting dimeric or oligomeric receptor complexes at the cell surface is modulated by self-interaction of their transmembrane (TM) domains. To address the importance of TM domain packing density for receptor functionality, we examined a set of asparagine mutants in the TM domain of the erythropoietin receptor (EpoR). We identified EpoR-T242N as a receptor variant that is present at the cell surface similar to wild-type EpoR but lacks visible localization in vesicle-like structures and is impaired in efficient activation of specific signaling cascades. Analysis by a molecular modeling approach indicated an increased interhelical distance for the EpoR-T242N TM dimer. By employing the model, we designed additional mutants with increased or decreased packing volume and confirmed a correlation between packing volume and biological responsiveness. These results propose that the packing density of the TM domain provides a novel layer for fine-tuned regulation of signal transduction and cellular decisions.

Monitoring orientation and dynamics of membrane-bound melittin utilizing dansyl fluorescence

Melittin is a cationic hemolytic peptide isolated from the European honey bee, Apis mellifera. In spite of a number of studies, there is no consensus regarding the orientation of melittin in membranes. In this study, we used a melittin analogue that is covalently labeled at its amino terminal (Gly-1) with the environment-sensitive 1-dimethylamino-5-sulfonylnaphthalene (dansyl) group to obtain information regarding the orientation and dynamics of the amino terminal region of membrane-bound melittin. Our results show that the dansyl group in Dns-melittin exhibits red edge excitation shift in vesicles of 1,2-dioleoyl-sn-glycero-3-phosphocholine, implying its localization in a motionally restricted region of the membrane. This is further supported by wavelength-dependent anisotropy and lifetime changes and time-resolved emission spectra characterized by dynamic Stokes shift, which indicates relatively slow solvent relaxation in the excited state. Membrane penetration depth analysis using the parallax method shows that the dansyl group is localized at a depth of approximately 18 A from the center of the bilayer in membrane-bound Dns-melittin. Further analysis of dansyl and tryptophan depths in Dns-melittin shows that the tilt angle between the helix axis of membrane-bound melittin and the bilayer normal is approximately 70 degrees. Our results therefore suggest that melittin adopts a pseudoparallel orientation in DOPC membranes at low concentration.

Investigating the mechanisms of photosynthetic proteins using continuum electrostatics

Computational methods based on continuum electrostatics are widely used in theoretical biochemistry to analyze the function of proteins. Continuum electrostatic methods in combination with quantum chemical and molecular mechanical methods can help to analyze even very complex biochemical systems. In this article, applications of these methods to proteins involved in photosynthesis are reviewed. After giving a short introduction to the basic concepts of the continuum electrostatic model based on the Poisson-Boltzmann equation, we describe the application of this approach to the docking of electron transfer proteins, to the comparison of isofunctional proteins, to the tuning of absorption spectra, to the analysis of the coupling of electron and proton transfer, to the analysis of the effect of membrane potentials on the energetics of membrane proteins, and to the kinetics of charge transfer reactions. Simulations as those reviewed in this article help to analyze molecular mechanisms on the basis of the structure of the protein, guide new experiments, and provide a better and deeper understanding of protein functions.

Pore formation in a lipid bilayer under a tension ramp: modeling the distribution of rupture tensions

The rupture of fluid membrane vesicles with a steady ramp of micropipette suction has been shown to produce a distribution of breakage tensions, with a mean that rises rapidly with tension rate. Starting from a lattice model that incorporates the essential features of the lipid bilayers held together with hydrophobic forces, and developing it to handle varying tension rates, we reproduce the main features of the experimental results. In essence, we show that the rupture kinetics are driven by the nucleation and growth of pores, with two limiting kinetics-growth-limited and nucleation-limited. The model has been extended to address the role of peptides in solution that can adsorb and insert themselves into the bilayer. At concentrations below those required to spontaneously rupture the membrane, the effect of the peptides is to lower the rupture tensions systematically for all tension rates.

A generalized born implicit-membrane representation compared to experimental insertion free energies

An implicit-membrane representation based on the generalized Born theory of solvation has been developed. The method was parameterized against the water-to-cyclohexane insertion free energies of hydrophobic side-chain analogs. Subsequently, the membrane was compared with experimental data from translocon inserted polypeptides and validated by comparison with an independent dataset of six membrane-associated peptides and eight integral membrane proteins of known structure and orientation. Comparison of the insertion energy of alpha-helical model peptides with the experimental values from the biological hydrophobicity scale of Hessa et al. gave a correlation of 93% with a mean unsigned error of 0.64 kcal/mol, when charged residues were ignored. The membrane insertion energy was found to be dependent on residue position. This effect is particularly pronounced for charged and polar residues, which strongly prefer interfacial locations. All integral membrane proteins investigated orient and insert correctly into the implicit-membrane model. Remarkably, the membrane model correctly predicts a partially inserted configuration for the monotopic membrane protein cyclooxygenase, matching experimental and theoretical predictions. To test the applicability and usefulness of the implicit-membrane method, molecular simulations of influenza A M2 as well as the glycophorin A dimer were performed. Both systems remain structurally stable and integrated into the membrane.

Influence of the Membrane Potential on the Protonation of Bacteriorhodopsin: Insights from Electrostatic Calculations into the Regulation of Proton Pumping

Proton binding and release are elementary steps for the transfer of protons within proteins, which is a process that is crucial in biochemical catalysis and biological energy transduction. Local electric fields in proteins affect the proton binding energy compared to aqueous solution. In membrane proteins, also the membrane potential affects the local electrostatics and can thus be crucial for protein function. In this paper, we introduce a procedure to calculate the protonation probability of titratable sites of a membrane protein in the presence of a membrane potential. In the framework of continuum electrostatics, we use a modified Poisson-Boltzmann equation to include the influence of the membrane potential. Our method considers that in a transmembrane protein each titratable site is accessible for protons from only one side of the membrane depending on the hydrogen bond pattern of the protein. We show that the protonation of sites receiving their protons from different sides of the membrane is differently influenced by the membrane potential. In addition, the effect of the membrane potential is combined with the effect of the pH gradient resulting from proton pumping. Our method is applied to bacteriorhodopsin, a light-activated proton pump. We find that the proton pumping of this protein might be regulated by Asp115, a conserved residue for which no function has been identified yet. According to our calculations, the interaction of Asp115 with Asp85 leads to the protonation of the latter if the pH gradient or the membrane potential becomes too large. Since Asp85 is the primary proton acceptor in the photocycle, bacteriorhodopsin molecules in which Asp85 is protonated cannot pump protons. Furthermore, we estimate how the membrane potential affects the energetics of the individual proton-transfer reactions of the photocycle. Most reactions, except the initial proton transfer from the Schiff base to Asp85, are influenced. Our calculations give new insights into the mechanism with which bacteriorhodopsin senses the membrane potential and the pH gradient and how the proton pumping is regulated by these parameters.

Conformation and environment of channel-forming peptides: a simulation study

Ion channel-forming peptides enable us to study the conformational dynamics of a transmembrane helix as a function of sequence and environment. Molecular dynamics simulations are used to study the conformation and dynamics of three 22-residue peptides derived from the second transmembrane domain of the glycine receptor (NK4-M2GlyR-p22). Simulations are performed on the peptide in four different environments: trifluoroethanol/water; SDS micelles; DPC micelles; and a DMPC bilayer. A hierarchy of alpha-helix stabilization between the different environments is observed such that TFE/water < micelles < bilayers. Local clustering of trifluoroethanol molecules around the peptide appears to help stabilize an alpha-helical conformation. Single (S22W) and double (S22W,T19R) substitutions at the C-terminus of NK4-M2GlyR-p22 help to stabilize a helical conformation in the micelle and bilayer environments. This correlates with the ability of the W22 and R19 side chains to form H-bonds with the headgroups of lipid or detergent molecules. This study provides a first atomic resolution comparison of the structure and dynamics of NK4-M2GlyR-p22 peptides in membrane and membrane-mimetic environments, paralleling NMR and functional studies of these peptides.

Evaluating tilt angles of membrane-associated helices: comparison of computational and NMR techniques

A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.

Energetics and stability of transmembrane helix packing: A replica-exchange simulation with a knowledge-based membrane potential

The energetics and stability of the packing of transmembrane helices were investigated by Monte Carlo simulations with the replica-exchange method. The helices were modeled with a united atom representation, and the CHARMM19 force field was employed. Based on known experimental structures of membrane proteins, an implicit knowledge-based potential was developed to describe the helix-membrane interactions at the residue level, whose validity was tested through prediction of the orientations when single helices were inserted into a membrane. Two systems were studied in this article, namely the glycophorin A dimer, and helices A and B of Bacteriorhodopsin. For the glycophorin A dimer, the most stable structure (0.5 A away from the experimental structure) is mainly stabilized by the favorable helix-helix interactions, and has the most population regardless of the helix-membrane interaction. However, for helices A and B of Bacteriorhodopsin, it was found that the packing determined by helix-helix interactions is nonspecific, and a native-like structure (0.2 A from the experimental one) can be identified from several structural analogs as the most stable one only after applying the membrane potential. Our results suggest that the contribution from the helix-membrane interaction could be critical in the correct packing of transmembrane helices in the membrane.

The α Helix Dipole: Screened Out?

Aligned alpha helix peptide dipoles sum to a "macroscopic" dipole parallel to the helix axis that has been implicated in protein folding and function. However, in aqueous solution the dipole is counteracted by an electrostatic reaction field generated by the solvent, and the strength of the helix dipole may reduce drastically from its value in vacuum. Here, using atomic-detail helix models and Poisson-Boltzmann continuum electrostatics calculations, the net effective dipole moment, mu(eff), is calculated. Some initially surprising results are found. Whereas in vacuum mu(eff) increases with helix length, the opposite is found to be the case for transmembrane helices. In soluble proteins, mu(eff) is found to vary strongly with the orientation and position of the helix relative to the aqueous medium. A set of rules is established to estimate of the strength of mu(eff) from graphical inspection of protein structures.