The importance of the membrane interface as the reference state for membrane protein stability

The insertion of nascent polypeptide chains into lipid bilayer membranes and the stability of membrane proteins crucially depend on the equilibrium partitioning of polypeptides. For this, the transfer of full sequences of amino-acid residues into the bilayer, rather than individual amino acids, must be understood. Earlier studies have revealed that the most likely reference state for partitioning very hydrophobic sequences is the membrane interface. We have used μs-scale simulations to calculate the interface-to-transmembrane partitioning free energies ΔG_S→TM for two hydrophobic carrier sequences in order to estimate the insertion free energy for all 20 amino acid residues when bonded to the center of a partitioning hydrophobic peptide. Our results show that prior single-residue scales likely overestimate the partitioning free energies of polypeptides. The correlation of ΔG_S→TM with experimental full-peptide translocon insertion data is high, suggesting an important role for the membrane interface in translocon-based insertion. The choice of carrier sequence greatly modulates the contribution of each single-residue mutation to the overall partitioning free energy. Our results demonstrate the importance of quantifying the observed full-peptide partitioning equilibrium, which is between membrane interface and transmembrane inserted, rather than combining individual water-to-membrane amino acid transfer free energies.     

Determining Peptide Partitioning Properties via Computer Simulation

The transfer of polypeptide segments into lipid bilayers to form transmembrane helices represents the crucial first step in cellular membrane protein folding and assembly. This process is driven by complex and poorly understood atomic interactions of peptides with the lipid bilayer environment. The lack of suitable experimental techniques that can resolve these processes both at atomic resolution and nanosecond timescales has spurred the development of computational techniques. In this review, we summarize the significant progress achieved in the last few years in elucidating the partitioning of peptides into lipid bilayer membranes using atomic detail molecular dynamics simulations. Indeed, partitioning simulations can now provide a wealth of structural and dynamic information. Furthermore, we show that peptide-induced bilayer distortions, insertion pathways, transfer free energies, and kinetic insertion barriers are now accurate enough to complement experiments. Further advances in simulation methods and force field parameter accuracy promise to turn molecular dynamics simulations into a powerful tool for investigating a wide range of membrane active peptide phenomena.     

Coarse-grained molecular dynamics simulations of the energetics of helix insertion into a lipid bilayer

Experimental and computational studies have indicated that hydrophobicity plays a key role in driving the insertion of transmembrane alpha-helices into lipid bilayers. Molecular dynamics simulations allow exploration of the nature of the interactions of transmembrane alpha-helices with their lipid bilayer environment. In particular, coarse-grained simulations have considerable potential for studying many aspects of membrane proteins, ranging from their self-assembly to the relation between their structure and function. However, there is a need to evaluate the accuracy of coarse-grained estimates of the energetics of transmembrane helix insertion. Here, three levels of complexity of model system have been explored to enable such an evaluation. First, calculated free energies of partitioning of amino acid side chains between water and alkane yielded an excellent correlation with experiment. Second, free energy profiles for transfer of amino acid side chains along the normal to a phosphatidylcholine bilayer were in good agreement with experimental and atomistic simulation studies. Third, estimation of the free energy profile for transfer of an arginine residue, embedded within a hydrophobic alpha-helix, to the center of a lipid bilayer gave a barrier of approximately 15 kT. Hence, there is a substantial barrier to membrane insertion for charged amino acids, but the coarse-grained model still underestimates the corresponding free energy estimate (approximately 29 kT) from atomistic simulations (Dorairaj, S., and Allen, T. W. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 4943-4948). Coarse-grained simulations were then used to predict the free energy profile for transfer of a simple model transmembrane alpha-helix (WALP23) across a lipid bilayer. The results indicated that a transmembrane orientation was favored by about -70 kT.     

Interaction of cardiotoxins with membranes: a molecular modeling study

Incorporation of beta-sheet proteins into membrane is studied theoretically for the first time, and the results are validated by the direct experimental data. Using Monte Carlo simulations with implicit membrane, we explore spatial structure, energetics, polarity, and mode of insertion of two cardiotoxins with different membrane-destabilizing activity. Both proteins, classified as P- and S-type cardiotoxins, are found to retain the overall "three-finger" fold interacting with membrane core and lipid/water interface by the tips of the "fingers" (loops). The insertion critically depends upon the structure, hydrophobicity, and electrostatics of certain regions. The simulations reveal apparently distinct binding modes for S- and P-type cardiotoxins via the first loop or through all three loops, respectively. This rationalizes an earlier empirical classification of cardiotoxins into S- and P-type, and provides a basis for the analysis of experimental data on their membrane affinities. Accomplished with our previous simulations of membrane alpha-helices, the computational method may be used to study partitioning of proteins with diverse folds into lipid bilayers.     

Membrane adsorption, folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Membrane adsorption,folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Counterion-mediated membrane penetration: cationic cell-penetrating peptides overcome Born energy barrier by ion-pairing with phospholipids

Arginine-rich cell-penetrating peptides (CPPs) can enter cells non-endocytotically, despite that transport of charge across a membrane should be formally associated with an extremely high Born energy barrier. We studied partitioning of several derivatives of the CPP penetratin in a water-octanol two-phase system in presence of natural phospholipids to explore if solvation by ion-pairing to hydrophobic counter-ions may serve as a mechanism for cell internalisation. We demonstrate that anionic lipids can aid peptide partitioning into octanol. Particularly efficient partitioning into octanol is observed with an arginine-rich penetratin compared to a lysine-rich derivative. Substituting tryptophans for phenylalanines results in poor partitioning into octanol, due to decreased overall peptide hydrophobicity. Partitioning into octanol is dependent of phospholipid type and the peptides induced structural changes in the lipid assemblies found in octanol. Attachment of carboxyfluorescein as a model cargo was found to enhance peptide partitioning into octanol. We discuss our results with respect to theoretical electrostatic energies, empirical hydrophobicity scales and in terms of implications for CPP uptake mechanisms. An important improvement of the theoretical transfer energies is obtained when, instead of singular ions, the insertion of ion-paired dipolar species is considered.     

Counterion-mediated membrane penetration: Cationic cell-penetrating peptides overcome Born energy barrier by ion-pairing with phospholipids

Arginine-rich cell-penetrating peptides (CPPs) can enter cells non-endocytotically, despite that transport of charge across a membrane should be formally associated with an extremely high Born energy barrier. We studied partitioning of several derivatives of the CPP penetratin in a water-octanol two-phase system in presence of natural phospholipids to explore if solvation by ion-pairing to hydrophobic counter-ions may serve as a mechanism for cell internalisation. We demonstrate that anionic lipids can aid peptide partitioning into octanol. Particularly efficient partitioning into octanol is observed with an arginine-rich penetratin compared to a lysine-rich derivative. Substituting tryptophans for phenylalanines results in poor partitioning into octanol, due to decreased overall peptide hydrophobicity. Partitioning into octanol is dependent of phospholipid type and the peptides induced structural changes in the lipid assemblies found in octanol. Attachment of carboxyfluorescein as a model cargo was found to enhance peptide partitioning into octanol. We discuss our results with respect to theoretical electrostatic energies, empirical hydrophobicity scales and in terms of implications for CPP uptake mechanisms. An important improvement of the theoretical transfer energies is obtained when, instead of singular ions, the insertion of ion-paired dipolar species is considered.     

Interfacial folding and membrane insertion of a designed helical peptide

Nonconstitutive membrane-active proteins, such as diphtheria toxin, must refold on membrane interfaces in the course of membrane penetration. A useful step in deciphering this process is to understand quantitatively the energetics of interface-mediated insertion of model transmembrane helices. A difficulty is that peptides that are sufficiently hydrophobic to span a lipid bilayer have a strong tendency to aggregate in the aqueous phase. To learn how to control the aqueous and membrane behavior of model peptides, we designed a 31-residue peptide (TMX-3) whose properties are described here. TMX-3 has two important structural features: a proline residue in the hydrophobic core that discourages the formation of highly helical aggregates in solution and two histidine residues that allow control of membrane and solution interactions by means of pH changes. The partitioning of TMX-3 into membranes followed complex kinetics, induced helicity, and shifted the histidine pK(a) from 6.8 to approximately 6. Topology measurements disclosed two general modes of TMX-3 binding: interfacial (IF) at low peptide concentrations and partial transmembrane (TM) insertion at higher concentrations. Both modes were reversible and, consequently, suitable for thermodynamic analysis. The free energies of IF partitioning of TMX-3 with deprotonated (pH 7.6) and protonated histidines (pH 4.5) were estimated by fluorescence titration to be -6.7 and -5.0 kcal/mol, respectively. These results show that histidine titration is likely to be important in the pH-dependent refolding of toxins on membrane interfaces and that the most favored state of TMX-3 under any conditions is the IF folded state, which emphasizes the importance of such states in the spontaneous refolding and insertion of diphtheria and other membrane toxins.     

Membrane partitioning of the cleavage peptide in flock house virus

Membrane translocation of the ssRNA genome of nodaviruses has been proposed to be mediated by direct lipid-protein interactions between a postassembly autocatalytic cleavage product from the capsomere and the target membrane. We have recently shown that the 21-residue Met-->Nle variant of the N-terminal helical domain (denoted gamma(1)) of the cleavage peptide in flock house nodavirus increases membrane permeability to hydrophilic solutes and can alter both membrane structure and function, suggesting the possibility of peptide-triggered disruption of the endosomal membrane as a prelude to viral uncoating in the host cytoplasm. Elucidation of partitioning energetics would allow an assessment of the likelihood of this mechanism. We report herein complete thermodynamic characterization of the partitioning of gamma(1) to phospholipids by lipid-peptide titrations following changes in ellipticity, fluorescence signature, or calorimetric response. These experiments revealed a partitioning energy comparable to natural membrane-active peptide toxins, suggesting that the proposed mechanism may be possible. Additionally, a novel switch in the balance of partitioning forces was found: when the lipid headgroup was changed from zwitterionic to negatively charged, membrane association of the peptide became completely entropy-driven.     

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations

Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at L_o/L_d interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.     

Folding and membrane insertion of the pore-forming peptide gramicidin occur as a concerted process

Many antibiotic peptides function by binding and inserting into membranes. Understanding this process provides an insight into the fundamentals of both membrane protein folding and antibiotic peptide function. For the first time, in this work, flow-aligned linear dichroism (LD) is used to study the folding of the antibiotic peptide gramicidin. LD provides insight into the combined processes of peptide folding and insertion and has the advantage over other similar techniques of being insensitive to off-membrane aggregation events. By combining LD data with conventional measurements of protein fluorescence and circular dichroism, the mechanism of gramicidin insertion is elucidated. The mechanism consists of five separately assignable steps that include formation of a water-insoluble gramicidin aggregate, dissociation from the aggregate, partitioning of peptide to the membrane surface, oligomerisation on the surface and concerted insertion and folding of the peptide to the double-helical form of gramicidin. Measurement of the rates of each step shows that although changes in the fluorescence signal cease 10 s after the initiation of the process, the insertion of the peptide into the membrane is actually not complete for a further 60 min. This last membrane insertion phase is only apparent by measurement of LD and circular dichroism signal changes. In summary, this study demonstrates the importance of multi-technique approaches, including LD, in studies of membrane protein folding.     

Implicit solvent model estimates of the stability of model structures of the alamethicin channel

Alamethicin is a hydrophobic helical peptide of 20 residues, which oligomerizes to form ion-conducting channels in membranes. The behavior of an intact alamethicin channel in POPC bilayers was recently studied, using 2 ns molecular dynamics (MD) simulations of a model hexameric channel. These simulations produced numerous conformations of the channel. In the present study, we used 11 of these channel conformations and carried out continuum-solvent model calculations, similar to those used for the monomers in our previous studies, to investigate the energetics of the channel inside the lipid bilayer. Our results suggest that, out of the 11 channel conformations produced by the MD simulations, only four are stable inside the lipid bilayer, with water-to-membrane free energies of transfer ranging from ~–6 to ~–10 kcal/mol. Analysis of the results suggests two causes for the apparent instability of the remainder of the structures inside the lipid bilayer, both resulting from the desolvation of channel polar groups (i.e. their transfer from the aqueous phase into the bilayer). The first is specific, uncompensated backbone hydrogen bonds, which exist in the region of the channel exposed to the hydrocarbon of the lipid bilayer. The second is exposure of intra-pore water molecules to the surrounding lipid. Thus, the association of these structures with the membrane involves a large electrostatic desolvation free-energy penalty. The apparent conflict between continuum-solvent and MD calculations, and its significance for the interpretation of membrane proteins simulations, are discussed.     

Arginine in Membranes: The Connection Between Molecular Dynamics Simulations and Translocon-Mediated Insertion Experiments

Several laboratories have carried out molecular dynamics (MD) simulations of arginine interactions with lipid bilayers and found that the energetic cost of placing arginine in lipid bilayers is an order of magnitude greater than observed in molecular biology experiments in which Arg-containing transmembrane helices are inserted across the endoplasmic reticulum membrane by the Sec61 translocon. We attempt here to reconcile the results of the two approaches. We first present MD simulations of guanidinium groups alone in lipid bilayers, and then, to mimic the molecular biology experiments, we present simulations of hydrophobic helices containing single Arg residues at different positions along the helix. We discuss the simulation results in the context of molecular biology results and show that the energetic discrepancy is reduced, but not eliminated, by considering free energy differences between Arg at the interface and at the center of the model helices. The reduction occurs because Arg snorkeling to the interface prevents Arg from residing in the bilayer center where the energetic cost of desolvation is highest. We then show that the problem with MD simulations is that they measure water-to-bilayer free energies, whereas the molecular biology experiments measure the energetics of partitioning from translocon to bilayer, which raises the fundamental question of the relationship between water-to-bilayer and water-to-translocon partitioning. We present two thermodynamic scenarios as a foundation for reconciliation of the simulation and molecular biology results. The simplest scenario is that translocon-to-bilayer partitioning is independent of water-to-bilayer partitioning; there is no thermodynamic cycle connecting the two paths.     

Sequence analysis and membrane partitioning energies of alpha-helical antimicrobial peptides

Sequences of 221 alpha-helical antimicrobial peptides (alphaAMPs) were compared and 63-166 of them were selected and analyzed using Perl programs. The results showed that aliphatic amino acids Gly, Leu, Ala, Ile and two positively charged amino acids Lys and Arg were composed of more than 63% of the first 20 residues of alphaAMPs. The weighed mean membrane partitioning energies at positions from 1 to 25 of alphaAMPs were calculated. Profile of the partitioning energies suggests oblique membrane insertion and an amphipathic alpha-helical structure of the N-terminus of alphaAMP (residues from 1 to 13), a bend structure at positions 13 and 14, and a less structured C-terminus that parallels the surface of the membrane. These structural features are in good agreement with the experimentally determined membrane structure of hemagglutinin fusion peptide from influenza virus. We hypothesize that this (N-terminal oblique alpha-helix)-central bend-(C-terminus) could be a common structural motif of membrane-disruptive peptides.     

A Monte Carlo Study of Peptide Insertion into Lipid Bilayers: Equilibrium Conformations and Insertion Mechanisms

The membrane insertion behavior of two peptides, Magainin2 and M2δ, was investigated by applying the Monte Carlo simulation technique to a theoretical model. The model included many novel aspects, such as a new semi-empirical lipid bilayer model and a new set of semi-empirical transfer energies, which reproduced the experimental insertion behavior of Magainin2 and M2δ without parameter fitting. Additionally, we have taken into account diminished internal (intramolecular) hydrogen bonding at the N- and C-termini of helical peptides. All simulations were carried out at 305 K, above the membrane thermal phase transition temperature, and at pH 7.0. The peptide equilibrium conformations are discussed for a range of bilayers with tail polarities varying from octanol-like to alkane-like. Probability distributions of the individual amino-acid-residue positions show the dynamic nature of these equilibrium conformations. Two different insertion mechanisms for M2δ, and a translocation mechanism for Magainin2, are described. A study of the effect of bilayer thickness on M2δ insertion suggests a critical thickness above which insertion is unfavorable. Additionally, we did not need to use an orientational potential or array of hard cylinders to persuade M2δ to insert perpendicular to the membrane surface. Instead, we found that diminished internal hydrogen bonding in the helical conformation anchored the termini in the headgroups and resulted in a nearly perpendicular orientation.     

pH-dependent self-association of influenza hemagglutinin fusion peptides in lipid bilayers

We have recently designed a host-guest peptide system that allows us to quantitatively measure the energetics of interaction of viral fusion peptides with lipid bilayers. Here, we show that fusion peptides of influenza hemagglutinin reversibly associate with one another at membrane surfaces above critical surface concentrations, which range from one to five peptides per 1000 lipids in the systems that we investigated. It is further demonstrated by using circular dichroism and Fourier transform infrared spectroscopy that monomeric peptides insert into the bilayers in a predominantly alpha-helical conformation, whereas self-associated fusion peptides adopt predominantly antiparallel beta-sheet structures at the membrane surface. The two forms are readily interconvertible and the equilibrium between them is determined by the pH and ionic strength of the surrounding solution. Lowering the pH favors the monomeric alpha-helical conformation, whereas increasing the ionic strength shifts the equilibrium towards the membrane-associated beta-aggregates. The binding data are interpreted in terms of a cooperative binding model that yields free energies of insertion and free energies of self-association for each of the peptides studied at pH 7.4 and pH 5. At pH 5 and 35 mM ionic strength, the insertion energy of the 20 residue influenza hemagglutinin fusion peptide is -7.2 kcal/mol and the self-association energy is -1.9 kcal/mol. We propose that self-association of fusion peptides could be a major driving force for recruiting a small number of hemagglutinin trimers into a fusion site.     

Vectorial insertion of a β-helical peptide into membrane: a theoretical study on polytheonamide B

Spontaneous unidirectional, or vectorial, insertion of transmembrane peptides is a fundamental biophysical process for toxin and viral actions. Polytheonamide B (pTB) is a potent cytotoxic peptide with a β^6.3-helical structure. Previous experimental studies revealed that the pTB inserts into the membrane in a vectorial fashion and forms a channel with its single molecular length long enough to span the membrane. Also, molecular dynamics simulation studies demonstrated that the pTB is prefolded in aqueous solution. These are unique features of pTB because most of the peptide toxins form channels through oligomerization of transmembrane helices. Here, we performed all-atom molecular dynamics simulations to examine the dynamic mechanism of the vectorial insertion of pTB, providing underlying elementary processes of the membrane insertion of a prefolded single transmembrane peptide. We find that the insertion of pTB proceeds with only the local lateral compression of the membrane in three successive phases: “landing,” “penetration,” and “equilibration” phases. The free energy calculations using the replica-exchange umbrella sampling simulations present an energy cost of 4.3 kcal/mol at the membrane surface for the membrane insertion of pTB from bulk water. The trajectories of membrane insertion revealed that the insertion process can occur in two possible pathways, namely “trapped” and “untrapped” insertions; in some cases, pTB is trapped in the upper leaflet during the penetration phase. Our simulations demonstrated the importance of membrane anchoring by the hydrophobic N-terminal blocking group in the landing phase, leading to subsequent vectorial insertion.