Protein contents in biological membranes can explain abnormal solvation of charged and polar residues

The Pathway of a Transmembrane Helix Insertion into the Membrane Assisted by Sec61α Channel

The significant inconsistency between the experimental and simulation results of the free energy for the translocon-assisted insertion of the transmembrane helix (TMH) has not been reasonably explained. Understanding the mechanism of TMH insertion through the translocon is the key to solving this problem. In this study, we performed a series of coarse-grained molecular dynamics simulations and calculated the potential mean forces (PMFs) for three insertion processes of a hydrophobic TMH. The simulations reveal the pathway of the TMH insertion assisted by a translocon. The results indicate that the TMH contacts the top of the lateral gate first and then inserts down the lateral gate, which agrees with the sliding model. The TMH begins to transfer laterally to the bilayer when it is blocked by the plug and reaches the exit of the lateral gate, where there is a free energy minimum point. We also found that the connecting section between TM2 and TM3 of Sec61α prevented TMH from leaving the lateral gate and directly transitioning to the surface-bound state. These findings provide insight into the mechanism of the insertion of TMH through the translocon.

Antimicrobial and Cell-Penetrating Peptides: Understanding Penetration for the Design of Novel Conjugate Antibiotics

Antimicrobial peptides (AMPs) are short oligopeptides that can penetrate the bacterial inner and outer membranes. Together with cell-penetrating peptides (CPPs), they are called membrane active peptides; peptides which can translocate across biological membranes. Over the last fifty years, attempts have been made to understand the molecular features that drive the interactions of membranes with membrane active peptides. This review examines the features of a membrane these peptides exploit for translocation, as well as the physicochemical characteristics of membrane active peptides which are important for translocation. Moreover, it presents examples of how these features have been used in recent years to create conjugates consisting of a membrane active peptide, called a "vector", attached to either a current or novel antibiotic, called a "cargo" or "payload". In addition, the review discusses what properties may contribute to an ideal peptide vector able to deliver cargoes across the bacterial outer membrane as the rising issue of antimicrobial resistance demands new strategies to be employed to combat this global public health threat.

A protet-based, protonic charge transfer model of energy coupling in oxidative and photosynthetic phosphorylation

Textbooks of biochemistry will explain that the otherwise endergonic reactions of ATP synthesis can be driven by the exergonic reactions of respiratory electron transport, and that these two half-reactions are catalyzed by protein complexes embedded in the same, closed membrane. These views are correct. The textbooks also state that, according to the chemiosmotic coupling hypothesis, a (or the) kinetically and thermodynamically competent intermediate linking the two half-reactions is the electrochemical difference of protons that is in equilibrium with that between the two bulk phases that the coupling membrane serves to separate. This gradient consists of a membrane potential term Δψ and a pH gradient term ΔpH, and is known colloquially as the protonmotive force or pmf. Artificial imposition of a pmf can drive phosphorylation, but only if the pmf exceeds some 150-170mV; to achieve in vivo rates the imposed pmf must reach 200mV. The key question then is 'does the pmf generated by electron transport exceed 200mV, or even 170mV?' The possibly surprising answer, from a great many kinds of experiment and sources of evidence, including direct measurements with microelectrodes, indicates it that it does not. Observable pH changes driven by electron transport are real, and they control various processes; however, compensating ion movements restrict the Δψ component to low values. A protet-based model, that I outline here, can account for all the necessary observations, including all of those inconsistent with chemiosmotic coupling, and provides for a variety of testable hypotheses by which it might be refined.

Molecular Dynamics Simulations Based on Polarizable Models Show that Ion Permeation Interconverts between Different Mechanisms as a Function of Membrane Thickness

Different mechanisms have been proposed to explain the permeation of charged compounds through lipid membranes. Overall, it is expected that an ion-induced defect permeation mechanism, where substantial membrane deformations accompany ion movement, should be dominant in thin membranes but that a solubility-diffusion mechanism, where ions partition into the membrane core with large associated dehydration energy costs, becomes dominant in thicker membranes. However, while this physical picture is intuitively reasonable, capturing the interconversion between these two permeation mechanisms in molecular dynamics (MD) simulations based on atomic models is challenging. In particular, simulations relying on nonpolarizable force fields are artificially unfavorable to the solubility-diffusion mechanism, as induced polarization of the nonpolar hydrocarbon is ignored, causing overestimated free energy costs for charged molecules to enter into this region of the membrane. In this study, all-atom MD simulations based on nonpolarizable and polarizable force fields are used to quantitatively characterize the permeation process for the arginine side chain analog methyl-guanidinium through bilayer membranes of mono-unsaturated phosphatidylcholine lipids with and without cholesterol, resulting in thicknesses spanning from ∼24 to ∼42 Å. With simulations based on a nonpolarizable force field, ion translocation can take place solely through an ion-induced defect mechanism, with free energy barriers increasing linearly from 14 to 40 kcal/mol, depending on the thickness. However, with simulations based on a polarizable force field, ion translocation is predominantly dominated by an ion-induced defect mechanism in thin membranes, which progressively converts to a solubility-diffusion mechanism as the membranes get thicker. The transition between the two mechanisms occurs at a thickness of ∼29 Å, with lipid tails of 22 or more carbon atoms. This situation appears to represent the upper limit for ion-induced defect permeation within the current polarizable models. Beyond this thickness, it becomes energetically preferable for the ion to dehydrate and partition into the membrane core-a phenomenon that cannot be captured using the nonpolarizable models. Induced electronic polarizability therefore leads not just to a shift in permeation energetics but to an interconversion between two strikingly different physical mechanisms. The result highlights the importance of induced polarizability in modeling lipid membranes.

Thermal conductivity of the cell membrane in the presence of cholesterol and amyloid precursor protein

The cell membrane is responsible for the transportation of heat between inside and outside the cell. Whether the thermal properties of the cell membrane are affected by the cholesterol concentration or the membrane proteins has not been investigated so far. Although the experimental measurement of the membrane thermal conductivity was not available until very recently, computational methods have been widely used for this purpose. In this study, we carry out molecular dynamics simulations to investigate the relation between the concentration of cholesterol and the thermal conductivity of a model membrane. Our results suggest an increase in the membrane thermal conductivity upon increasing the concentration of cholesterol in the membrane. Moreover, we find that the asymmetric distribution of cholesterol in the two membrane leaflets decreases thermal conductivity. We also find a rectification effect when heat flows in opposite directions through a model membrane decorated with the amyloid precursor protein. The results of this study apply to the advancement of selective treatment methods, as well as the development of new materials such as biological rectifiers.

Insertion of Bacteriorhodopsin Helix C Variants into Biological Membranes

A peptide corresponding to bacteriorhodopsin (bR) helix C, later named pHLIP, inserts across lipid bilayers as a monomeric α-helix at acidic pH, but is an unstructured surface-bound monomer at neutral pH. As a result of such pH-responsiveness, pHLIP targets acidic tumors and has been used as a vehicle for imaging and drug-delivery cargoes. To gain insights about the insertion of bR helix C into biological membranes, we replaced two key aspartic residues that control the topological transition from the aqueous phase into a lipid bilayer. Here, we used an in vitro transcription-translation system to study the translocon-mediated insertion of helix C-derived segments into rough microsomes. Our data provide the first quantitative biological understanding of this effect. Interestingly, replacing the aspartic residues by glutamic residues does not significantly alters the insertion propensity, while replacement by alanines promotes a transmembrane orientation. These results are consistent with mutational data obtained in synthetic liposomes by manipulating pH conditions. Our findings support the notion that the translocon facilitates topogenesis under physiological pH conditions.

Site-Specific Peptide Probes Detect Buried Water in a Lipid Membrane

Transmembrane peptides contain polar residues in the interior of the membrane, which may alter the electrostatic environment and favor hydration in the otherwise nonpolar environment of the membrane core. Here, we demonstrate a general, nonperturbative strategy to probe hydration of the peptide backbone at specific depths within the bilayer using a combination of site-specific isotope labels, ultrafast two-dimensional infrared spectroscopy, and spectral modeling based on molecular dynamics simulations. Our results show that the amphiphilic pH-low insertion peptide supports a highly heterogeneous environment, with significant backbone hydration of nonpolar residues neighboring charged residues. For example, a leucine residue located as far as 1 nm into the hydrophobic bulk reports hydrogen-bonded populations as high as ∼20%. These findings indicate that the polar nature of these residues may facilitate the transport of water molecules into the hydrophobic core of the membrane.

Characterizing the Molecular Mechanisms for Flipping Charged Peptide Flanking Loops across a Lipid Bilayer

The cell membrane largely prevents the passive diffusion of charged molecules due to the large free energy barrier associated with translocating charged groups across the hydrophobic lipid bilayer core. Despite this barrier, some peptides can interconvert between transmembrane and surface-adsorbed states by "flipping" charged flanking loops across the bilayer on a surprisingly rapid second-minute time scale. The transmembrane helices of some multispanning membrane proteins undergo similar reorientation processes, suggesting that loop-flipping may be a mechanism for regulating membrane protein topology; however, the molecular mechanisms underlying this behavior remain unknown. In this work, we study the loop-flipping behavior exhibited by a peptide with a hydrophobic transmembrane helix, charged flanking loops, and a central, membrane-exposed aspartate residue of varying protonation state. We utilize all-atom temperature accelerated molecular dynamics simulations to predict the likelihood of loop-flipping without predefining specific loop-flipping pathways. We demonstrate that this approach can identify multiple possible flipping pathways, with the prevalence of each pathway depending on the protonation state of the central residue. In particular, we find that a charged central residue facilitates loop-flipping by stabilizing membrane water defects, enabling the "self-catalysis" of charge translocation. These findings provide detailed molecular-level insights into charged loop-flipping pathways that may generalize to other charge translocation processes, such as lipid flip-flop or the large-scale conformational rearrangements of multispanning membrane proteins.

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.

Free Energy Calculations of Membrane Permeation: Challenges Due to Strong Headgroup-Solute Interactions

Understanding how different classes of molecules move across biological membranes is a prerequisite to predicting a solute's permeation rate, which is a critical factor in the fields of drug design and pharmacology. We use biased molecular dynamics computer simulations to calculate and compare the free energy profiles of translocation of several small molecules across 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC) lipid bilayers as a first step toward determining the most efficient method for free energy calculations. We study the translocation of arginine, a sodium ion, alanine, and a single water molecule using the metadynamics, umbrella sampling, and replica exchange umbrella sampling techniques. Within the fixed lengths of our simulations, we find that all methods produce similar results for charge-neutral permeants, but not for polar or positively charged molecules. We identify the long relaxation time scale of electrostatic interactions between lipid headgroups and the solute to be the principal cause of this difference and show that this slow process can lead to an erroneous dependence of computed free energy profiles on the initial system configuration. We demonstrate the use of committor analysis to validate the proper sampling of the presumed transition state, which in our simulations is achieved only in replica exchange calculations. On the basis of these results we provide some useful guidance to perform and evaluate free energy calculations of membrane permeation.

Insight into the interactions, residue snorkeling, and membrane disordering potency of a single antimicrobial peptide into different lipid bilayers

Pardaxin, with a bend-helix-bend-helix structure, is a membrane-active antimicrobial peptide that its membrane activity depends on the lipid bilayer composition. Herein, all-atom molecular dynamics (MD) simulations were performed to provide further molecular insight into the interactions, structural dynamics, orientation behavior, and cationic residues snorkeling of pardaxin in the DMPC, DPPC, POPC, POPG, POPG/POPE (3:1), and POPG/POPE (1:3) lipid bilayers. The results showed that the C-terminal helix of the peptide was maintained in all six types of the model-bilayers and pardaxin was tilted into the DMPC, DPPC, and POPG/POPE mixed bilayers more than the POPC and POPG bilayers. As well as, the structure of zwitterionic membranes was more affected by the peptide than the anionic bilayers. Taken together, the study demonstrated that the cationic residues of pardaxin snorkeled toward the interface of lipid bilayers and all phenylalanine residues of the peptide played important roles in the peptide-membrane interactions. We hope that this work will provide a better understanding of the interactions of antimicrobial peptides with the membranes.

Grafting Charged Species to Membrane-Embedded Scaffolds Dramatically Increases the Rate of Bilayer Flipping

The cell membrane is a barrier to the passive diffusion of charged molecules due to the chemical properties of the lipid bilayer. Surprisingly, recent experiments have identified processes in which synthetic and biological charged species directly transfer across lipid bilayers on biologically relevant time scales. In particular, amphiphilic nanoparticles have been shown to insert into lipid bilayers, requiring the transport of charged species across the bilayer. The molecular factors facilitating this rapid insertion process remain unknown. In this work, we use atomistic molecular dynamics simulations to calculate the free energy barrier associated with "flipping" charged species across a lipid bilayer for species that are grafted to a membrane-embedded scaffold, such as a membrane-embedded nanoparticle. We find that the free energy barrier for flipping a grafted ligand can be over 7 kcal/mol lower than the barrier for translocating an isolated, equivalent ion, yielding a 5 order of magnitude decrease in the corresponding flipping time scale. Similar results are found for flipping charged species grafted to either nanoparticle or protein scaffolds. These results reveal new mechanistic insight into the flipping of charged macromolecular components that might play an important, yet overlooked, role in signaling and charge transport in biological settings. Furthermore, our results suggest guidelines for the design of synthetic materials capable of rapidly flipping charged moieties across the cell membrane.

Comprehensive analysis of the numbers, lengths and amino acid compositions of transmembrane helices in prokaryotic, eukaryotic and viral integral membrane proteins of high-resolution structure

We report a comprehensive analysis of the numbers, lengths and amino acid compositions of transmembrane helices in 235 high-resolution structures of integral membrane proteins. The properties of 1551 transmembrane helices in the structures were compared with those obtained by analysis of the same amino acid sequences using topology prediction tools. Explanations for the 81 (5.2%) missing or additional transmembrane helices in the prediction results were identified. Main reasons for missing transmembrane helices were mis-identification of N-terminal signal peptides, breaks in α-helix conformation or charged residues in the middle of transmembrane helices and transmembrane helices with unusual amino acid composition. The main reason for additional transmembrane helices was mis-identification of amphipathic helices, extramembrane helices or hairpin re-entrant loops. Transmembrane helix length had an overall median of 24 residues and an average of 24.9 ± 7.0 residues and the most common length was 23 residues. The overall content of residues in transmembrane helices as a percentage of the full proteins had a median of 56.8% and an average of 55.7 ± 16.0%. Amino acid composition was analysed for the full proteins, transmembrane helices and extramembrane regions. Individual proteins or types of proteins with transmembrane helices containing extremes in contents of individual amino acids or combinations of amino acids with similar physicochemical properties were identified and linked to structure and/or function. In addition to overall median and average values, all results were analysed for proteins originating from different types of organism (prokaryotic, eukaryotic, viral) and for subgroups of receptors, channels, transporters and others.

Interplay between hydrophobicity and the positive-inside rule in determining membrane-protein topology

The energetics of membrane-protein interactions determine protein topology and structure: hydrophobicity drives the insertion of helical segments into the membrane, and positive charges orient the protein with respect to the membrane plane according to the positive-inside rule. Until recently, however, quantifying these contributions met with difficulty, precluding systematic analysis of the energetic basis for membrane-protein topology. We recently developed the dsTβL method, which uses deep sequencing and in vitro selection of segments inserted into the bacterial plasma membrane to infer insertion-energy profiles for each amino acid residue across the membrane, and quantified the insertion contribution from hydrophobicity and the positive-inside rule. Here, we present a topology-prediction algorithm called TopGraph, which is based on a sequence search for minimum dsTβL insertion energy. Whereas the average insertion energy assigned by previous experimental scales was positive (unfavorable), the average assigned by TopGraph in a nonredundant set is -6.9 kcal/mol. By quantifying contributions from both hydrophobicity and the positive-inside rule we further find that in about half of large membrane proteins polar segments are inserted into the membrane to position more positive charges in the cytoplasm, suggesting an interplay between these two energy contributions. Because membrane-embedded polar residues are crucial for substrate binding and conformational change, the results implicate the positive-inside rule in determining the architectures of membrane-protein functional sites. This insight may aid structure prediction, engineering, and design of membrane proteins. TopGraph is available online (topgraph.weizmann.ac.il).

Amino acid composition analysis of human secondary transport proteins and implications for reliable membrane topology prediction

Secondary transporters in humans are a large group of proteins that transport a wide range of ions, metals, organic and inorganic solutes involved in energy transduction, control of membrane potential and osmotic balance, metabolic processes and in the absorption or efflux of drugs and xenobiotics. They are also emerging as important targets for development of new drugs and as target sites for drug delivery to specific organs or tissues. We have performed amino acid composition (AAC) and phylogenetic analyses and membrane topology predictions for 336 human secondary transport proteins and used the results to confirm protein classification and to look for trends and correlations with structural domains and specific substrates and/or function. Some proteins showed statistically high contents of individual amino acids or of groups of amino acids with similar physicochemical properties. One recurring trend was a correlation between high contents of charged and/or polar residues with misleading results in predictions of membrane topology, which was especially prevalent in Mitochondrial Carrier family proteins. We demonstrate how charged or polar residues located in the middle of transmembrane helices can interfere with their identification by membrane topology tools resulting in missed helices in the prediction. Comparison of AAC in the human proteins with that in 235 secondary transport proteins from Escherichia coli revealed similar overall trends along with differences in average contents for some individual amino acids and groups of similar amino acids that are presumed to result from a greater number of functions and complexity in the higher organism.

Decrypting protein insertion through the translocon with free-energy calculations

Protein insertion into a membrane is a complex process involving numerous players. The most prominent of these players is the Sec translocon complex, a conserved protein-conducting channel present in the cytoplasmic membrane of bacteria and the membrane of the endoplasmic reticulum in eukaryotes. The last decade has seen tremendous leaps forward in our understanding of how insertion is managed by the translocon and its partners, coming from atomic-detailed structures, innovative experiments, and well-designed simulations. In this review, we discuss how experiments and simulations, hand-in-hand, teased out the secrets of the translocon-facilitated membrane insertion process. In particular, we focus on the role of free-energy calculations in elucidating membrane insertion. Amazingly, despite all its apparent complexity, protein insertion into membranes is primarily driven by simple thermodynamic and kinetic principles. This article is part of a Special Issue entitled: Membrane proteins edited by J.C. Gumbart and Sergei Noskov.

Mutational scanning reveals the determinants of protein insertion and association energetics in the plasma membrane

Insertion of helix-forming segments into the membrane and their association determines the structure, function, and expression levels of all plasma membrane proteins. However, systematic and reliable quantification of membrane-protein energetics has been challenging. We developed a deep mutational scanning method to monitor the effects of hundreds of point mutations on helix insertion and self-association within the bacterial inner membrane. The assay quantifies insertion energetics for all natural amino acids at 27 positions across the membrane, revealing that the hydrophobicity of biological membranes is significantly higher than appreciated. We further quantitate the contributions to membrane-protein insertion from positively charged residues at the cytoplasm-membrane interface and reveal large and unanticipated differences among these residues. Finally, we derive comprehensive mutational landscapes in the membrane domains of Glycophorin A and the ErbB2 oncogene, and find that insertion and self-association are strongly coupled in receptor homodimers.

Protein-fluctuation-induced water-pore formation in ion channel voltage-sensor translocation across a lipid bilayer membrane

We have applied a combined fluorescence microscopy and single-ion-channel electric current recording approach, correlating with molecular dynamics (MD) simulations, to study the mechanism of voltage-sensor domain translocation across a lipid bilayer. We use the colicin Ia ion channel as a model system, and our experimental and simulation results show the following: (1) The open-close activity of an activated colicin Ia is not necessarily sensitive to the amplitude of the applied cross-membrane voltage when the cross-membrane voltage is around the resting potential of excitable membranes; and (2) there is a significant probability that the activation of colicin Ia occurs by forming a transient and fluctuating water pore of ∼15 Å diameter in the lipid bilayer membrane. The location of the water-pore formation is nonrandom and highly specific, right at the insertion site of colicin Ia charged residues in the lipid bilayer membrane, and the formation is intrinsically associated with the polypeptide conformational fluctuations and solvation dynamics. Our results suggest an interesting mechanistic pathway for voltage-sensitive ion channel activation, and specifically for translocation of charged polypeptide chains across the lipid membrane under a transmembrane electric field: the charged polypeptide domain facilitates the formation of hydrophilic water pore in the membrane and diffuses through the hydrophilic pathway across the membrane; i.e., the charged polypeptide chain can cross a lipid membrane without entering into the hydrophobic core of the lipid membrane but entirely through the aqueous and hydrophilic environment to achieve a cross-membrane translocation. This mechanism sheds light on the intensive and fundamental debate on how a hydrophilic and charged peptide domain diffuses across the biologically inaccessible high-energy barrier of the hydrophobic core of a lipid bilayer: The peptide domain does not need to cross the hydrophobic core to move across a lipid bilayer.

Effects of Lys to Glu mutations in GsMTx4 on membrane binding, peptide orientation, and self-association propensity, as analyzed by molecular dynamics simulations

GsMTx4, a gating modifier peptide acting on cationic mechanosensitive channels, has a positive charge (+5e) due to six Lys residues. The peptide does not have a stereospecific binding site on the channel but acts from the boundary lipids within a Debye length of the pore probably by changing local stress. To gain insight into how these Lys residues interact with membranes, we performed molecular dynamics simulations of Lys to Glu mutants in parallel with our experimental work. In silico, K15E had higher affinity for 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine bilayers than wild-type (WT) peptide or any other mutant tested, and showed deeper penetration than WT, a finding consistent with the experimental data. Experimentally, the inhibitory activities of K15E and K25E were most compromised, whereas K8E and K28E inhibitory activities remained similar to WT peptide. Binding of WT in an interfacial mode did not influence membrane thickness. With interfacial binding, the direction of the dipole moments of K15E and K25E was predicted to differ from WT, whereas those of K8E and K28E oriented similarly to that of WT. These results support a model in which binding of GsMTx4 to the membrane acts like an immersible wedge that serves as a membrane expansion buffer reducing local stress and thus inhibiting channel activity. In simulations, membrane-bound WT attracted other WT peptides to form aggregates. This may account for the positive cooperativity observed in the ion channel experiments. The Lys residues seem to fine-tune the depth of membrane binding, the tilt angle, and the dipole moments.

Evaluating Force Fields for the Computational Prediction of Ionized Arginine and Lysine Side-Chains Partitioning into Lipid Bilayers and Octanol

Abundant peptides and proteins containing arginine (Arg) and lysine (Lys) amino acids can apparently permeate cell membranes with ease. However, the mechanisms by which these peptides and proteins succeed in traversing the free energy barrier imposed by cell membranes remain largely unestablished. Precise thermodynamic studies (both theoretical and experimental) on the interactions of Arg and Lys residues with model lipid bilayers can provide valuable clues to the efficacy of these cationic peptides and proteins. We have carried out molecular dynamics simulations to calculate the interactions of ionized Arg and Lys side-chains with the zwitterionic 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayer for 10 widely used lipid/protein force fields: CHARMM36/CHARMM36, SLIPID/AMBER99SB-ILDN, OPLS-AA/OPLS-AA, Berger/OPLS-AA, Berger/GROMOS87, Berger/GROMOS53A6, GROMOS53A6/GROMOS53A6, nonpolarizable MARTINI, polarizable MARTINI, and BMW MARTINI. We performed umbrella sampling simulations to obtain the potential of mean force for Arg and Lys side-chains partitioning from water to the bilayer interior. We found significant differences between the force fields, both for the interactions between side-chains and bilayer surface, as well as the free energy cost for placing the side-chain at the center of the bilayer. These simulation results were compared with the Wimley-White interfacial scale. We also calculated the free energy cost for transferring ionized Arg and Lys side-chains from water to both dry and wet octanol. Our simulations reveal rapid diffusion of water molecules into octanol whereby the equilibrium mole fraction of water in the wet octanol phase was ∼25%. Surprisingly, our free energy calculations found that the high water content in wet octanol lowered the water-to-octanol partitioning free energies for cationic residues by only 0.6 to 0.7 kcal/mol.

Folding of Aquaporin 1: multiple evidence that helix 3 can shift out of the membrane core

The folding of most integral membrane proteins follows a two-step process: initially, individual transmembrane helices are inserted into the membrane by the Sec translocon. Thereafter, these helices fold to shape the final conformation of the protein. However, for some proteins, including Aquaporin 1 (AQP1), the folding appears to follow a more complicated path. AQP1 has been reported to first insert as a four-helical intermediate, where helix 2 and 4 are not inserted into the membrane. In a second step, this intermediate is folded into a six-helical topology. During this process, the orientation of the third helix is inverted. Here, we propose a mechanism for how this reorientation could be initiated: first, helix 3 slides out from the membrane core resulting in that the preceding loop enters the membrane. The final conformation could then be formed as helix 2, 3, and 4 are inserted into the membrane and the reentrant regions come together. We find support for the first step in this process by showing that the loop preceding helix 3 can insert into the membrane. Further, hydrophobicity curves, experimentally measured insertion efficiencies and MD-simulations suggest that the barrier between these two hydrophobic regions is relatively low, supporting the idea that helix 3 can slide out of the membrane core, initiating the rearrangement process.

Sidekick for Membrane Simulations: Automated Ensemble Molecular Dynamics Simulations of Transmembrane Helices

The interactions of transmembrane (TM) α-helices with the phospholipid membrane and with one another are central to understanding the structure and stability of integral membrane proteins. These interactions may be analyzed via coarse grained molecular dynamics (CGMD) simulations. To obtain statistically meaningful analysis of TM helix interactions, large (N ca. 100) ensembles of CGMD simulations are needed. To facilitate the running and analysis of such ensembles of simulations, we have developed Sidekick, an automated pipeline software for performing high throughput CGMD simulations of α-helical peptides in lipid bilayer membranes. Through an end-to-end approach, which takes as input a helix sequence and outputs analytical metrics derived from CGMD simulations, we are able to predict the orientation and likelihood of insertion into a lipid bilayer of a given helix of a family of helix sequences. We illustrate this software via analyses of insertion into a membrane of short hydrophobic TM helices containing a single cationic arginine residue positioned at different positions along the length of the helix. From analyses of these ensembles of simulations, we estimate apparent energy barriers to insertion which are comparable to experimentally determined values. In a second application, we use CGMD simulations to examine the self-assembly of dimers of TM helices from the ErbB1 receptor tyrosine kinase and analyze the numbers of simulation repeats necessary to obtain convergence of simple descriptors of the mode of packing of the two helices within a dimer. Our approach offers a proof-of-principle platform for the further employment of automation in large ensemble CGMD simulations of membrane proteins.

Why is the biological hydrophobicity scale more accurate than earlier experimental hydrophobicity scales?

The recognition of transmembrane helices by the translocon is primarily guided by the average hydrophobicity of the potential transmembrane helix. However, the exact hydrophobicity of each amino acid can be identified in several different ways. The free energy of transfer for amino acid analogues between a hydrophobic media, for example, octanol and water can be measured or obtained from simulations, the hydrophobicity can also be estimated by statistical properties from known transmembrane segments and finally the contribution of each amino acid type for the probability of translocon recognition has recently been measured directly. Although these scales correlate quite well, there are clear differences between them and it is not well understood which scale represents neither the biology best nor what the differences are. Here, we try to provide some answers to this by studying the ability of different scales to recognize transmembrane helices and predict the topology of transmembrane proteins. From this analysis it is clear that the biological hydrophobicity scale as well scales created from statistical analysis of membrane helices perform better than earlier experimental scales that are mainly based on measurements of amino acid analogs and not directly on transmembrane helix recognition. Using these results we identified the properties of the scales that perform better than other scales. We find, for instance, that the better performing scales consider proline more hydrophilic. This shows that transmembrane recognition is not only governed by pure hydrophobicity but also by the helix preferences for amino acids, as proline is a strong helix breaker.

The effect of structural parameters and positive charge distance on the interaction free energy of antimicrobial peptides with membrane surface

Many attempts have been made to find hints explaining the relationship between physicochemical and structural properties of antimicrobial peptides (AMPs) which are relevant to their antimicrobial activities. We here found that there is a difference in the percentages of hydrophobic, hydrophilic, and charged residues between AMPs killing both bacteria and fungi (Group A) and AMPs that only kill bacteria (Group B). The percentage of charged residues in Group A AMPs is highly elevated, while in Group B the percentage of hydrophobic residues is increased. This result suggests a sequence-based mechanism of selectivity for AMPs. Moreover, we examined how the distance between basic residues affects the interaction free energy of AMPs with the membrane surface, since most of the known AMPs act by membrane perturbation. We measured the average distance between basic residues throughout the 3D structure of AMPs by defining Dpr parameter and calculated the interaction free energy for 10 AMPs that interacted with the DPPC membrane using molecular dynamics simulation. We found that the changes of the interaction free energy correlates with the change of Dpr by a linear regression coefficient of r(2 )= .47 and a cubic regression coefficient of r(2 )= .70.

Implicit membrane treatment of buried charged groups: application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

On the nature of the apparent free energy of inserting amino acids into membrane through the translocon

The nature of the biological free energy scale (ΔGapp), obtained from translocon mediated insertion studies, has been a major puzzle and the subject of major controversies. Part of the problem has been the complexity of the insertion process that discouraged workers from considering the feasible kinetics schemes and left the possible impression that ΔGapp presents some simple partition. Here we extend and clarify our recent analysis of the insertion problem using well-defined kinetics schemes and a free energy profile. We point out that although the rate constants of some steps are far from being obvious, it is essential to consider explicitly such schemes in order to advance in analyzing the meaning of ΔGapp. It is then shown that under some equilibrium conditions the kinetics scheme leads to a simple formula that allows one to relate ΔGapp to the actual free energy of partitioning between the water, the membrane, and the translocon. Other options are also considered (including limits with irreversible transitions that can be described by linear free energy relationships (LFERs)). It is concluded that it is unlikely that a kinetics plus thermodynamic based analysis can lead to a result that identifies ΔGapp with the partition between the membrane and the translocon. Thus, we argue that unless such analysis is presented, it is unjustified to assume that ΔGapp corresponds to the membrane translocon equilibrium or to some other arbitrary definition. Furthermore, we point out that the presumption that it is sufficient to just calculate the PMF for going from the translocon (TR) to the membrane and then to assume irreversible diffusive motion to water and for further entrance to the membrane is not a valid analysis. Overall, we point out that it is important to try to relate ΔGapp to a well-defined kinetics scheme (regardless of the complication of the system) in order to determine whether the energies of inserting positively charged residues to the membrane are related to the corresponding ΔGapp. It is also suggested that deviations from our simple formula for equilibrium conditions can help in identifying and analyzing kinetics barriers.

The different interactions of lysine and arginine side chains with lipid membranes

The basic amino acids lysine (Lys) and arginine (Arg) play important roles in membrane protein activity, the sensing of membrane voltages, and the actions of antimicrobial, toxin, and cell-penetrating peptides. These roles are thought to stem from the strong interactions and disruptive influences of these amino acids on lipid membranes. In this study, we employ fully atomistic molecular dynamics simulations to observe, quantify, and compare the interactions of Lys and Arg with saturated phosphatidylcholine membranes of different thickness. We make use of both charged (methylammonium and methylguanidinium) and neutral (methylamine and methylguanidine) analogue molecules, as well as Lys and Arg side chains on transmembrane helix models. We find that the free energy barrier experienced by a charged Lys crossing the membrane is strikingly similar to that of a charged Arg (to within 2 kcal/mol), despite the two having different chemistries, H-bonding capability, and hydration free energies that differ by ∼10 kcal/mol. In comparison, the barrier for neutral Arg is higher than that for neutral Lys by around 5 kcal/mol, being more selective than that for the charged species. This can be explained by the different transport mechanisms for charged or neutral amino acid side chains in the membrane, involving membrane deformations or simple dehydration, respectively. As a consequence, we demonstrate that Lys would be deprotonated in the membrane, whereas Arg would maintain its charge. Our simulations also reveal that Arg attracts more phosphate and water in the membrane, and can form extensive H-bonding with its five H-bond donors to stabilize Arg-phosphate clusters. This leads to enhanced interfacial binding and membrane perturbations, including the appearance of a trans-membrane pore in a thinner membrane. These results highlight the special role played by Arg as an amino acid to bind to, disrupt, and permeabilize lipid membranes, as well as to sense voltages for a range of peptide and protein activities in nature and in engineered bionanodevices.

The simulation approach to lipid-protein interactions

The interactions between lipids and proteins are crucial for a range of biological processes, from the folding and stability of membrane proteins to signaling and metabolism facilitated by lipid-binding proteins. However, high-resolution structural details concerning functional lipid/protein interactions are scarce due to barriers in both experimental isolation of native lipid-bound complexes and subsequent biophysical characterization. The molecular dynamics (MD) simulation approach provides a means to complement available structural data, yielding dynamic, structural, and thermodynamic data for a protein embedded within a physiologically realistic, modelled lipid environment. In this chapter, we provide a guide to current methods for setting up and running simulations of membrane proteins and soluble, lipid-binding proteins, using standard atomistically detailed representations, as well as simplified, coarse-grained models. In addition, we outline recent studies that illustrate the power of the simulation approach in the context of biologically relevant lipid/protein interactions.

Assembly and Stability of α-Helical Membrane Proteins

Grease to grease - this is how one might begin to describe the tendency of hydrophobic stretches in protein amino acid sequences to form transmembrane domains. While this simple rule contains a lot of truth, the mechanisms of membrane protein folding, the insertion of hydrophobic protein domains into the lipid bilayer, and the apparent existence of highly polar residues in some proteins in the hydrophobic membrane core are subjects of lively debate - an indication that many details remain unresolved. Here, we present a historical survey of recent insights from experiments and computational studies into the rules and mechanisms of α-helical membrane protein assembly and stability.

Direct simulation of early-stage Sec-facilitated protein translocation

Direct simulations reveal key mechanistic features of early-stage protein translocation and membrane integration via the Sec-translocon channel. We present a novel computational protocol that combines non-equilibrium growth of the nascent protein with microsecond timescale molecular dynamics trajectories. Analysis of multiple, long timescale simulations elucidates molecular features of protein insertion into the translocon, including signal-peptide docking at the translocon lateral gate (LG), large lengthscale conformational rearrangement of the translocon LG helices, and partial membrane integration of hydrophobic nascent-protein sequences. Furthermore, the simulations demonstrate the role of specific molecular interactions in the regulation of protein secretion, membrane integration, and integral membrane protein topology. Salt-bridge contacts between the nascent-protein N-terminus, cytosolic translocon residues, and phospholipid head groups are shown to favor conformations of the nascent protein upon early-stage insertion that are consistent with the Type II (N(cyt)/C(exo)) integral membrane protein topology, and extended hydrophobic contacts between the nascent protein and the membrane lipid bilayer are shown to stabilize configurations that are consistent with the Type III (N(exo)/C(cyt)) topology. These results provide a detailed, mechanistic basis for understanding experimentally observed correlations between integral membrane protein topology, translocon mutagenesis, and nascent-protein sequence.

Protein Structure in Membrane Domains

Of great interest to the academic and pharmaceutical research communities, helical transmembrane proteins are characterized by their ability to dissolve and fold in lipid bilayers—properties conferred by polypeptide spans termed transmembrane domains (TMDs). The apolar nature of TMDs necessitates the use of membrane-mimetic solvents for many structure and folding studies. This review examines the relationship between TMD structure and solvent environment, focusing on principles elucidated largely in membrane-mimetic environments with single-TMD protein and peptide models. Following a brief description of TMD sequence and conformational characteristics gleaned from the structural database, we present an overview of the conceptual models used to study folding in vitro. The impact of sequence and solvent context on the incorporation of TMDs into membranes, and its role in measurements of TMD self-assembly strengths, is then described. We conclude with a discussion of the nonspecific effects of membrane components on TMD stability.

Determination of membrane-insertion free energies by molecular dynamics simulations

The accurate prediction of membrane-insertion probability for arbitrary protein sequences is a critical challenge to identifying membrane proteins and determining their folded structures. Although algorithms based on sequence statistics have had moderate success, a complete understanding of the energetic factors that drive the insertion of membrane proteins is essential to thoroughly meeting this challenge. In the last few years, numerous attempts to define a free-energy scale for amino-acid insertion have been made, yet disagreement between most experimental and theoretical scales persists. However, for a recently resolved water-to-bilayer scale, it is found that molecular dynamics simulations that carefully mimic the conditions of the experiment can reproduce experimental free energies, even when using the same force field as previous computational studies that were cited as evidence of this disagreement. Therefore, it is suggested that experimental and simulation-based scales can both be accurate and that discrepancies stem from disparities in the microscopic processes being considered rather than methodological errors. Furthermore, these disparities make the development of a single universally applicable membrane-insertion free energy scale difficult.

The role of membrane thickness in charged protein-lipid interactions

Charged amino acids are known to be important in controlling the actions of integral and peripheral membrane proteins and cell disrupting peptides. Atomistic molecular dynamics studies have shed much light on the mechanisms of membrane binding and translocation of charged protein groups, yet the impact of the full diversity of membrane physico-chemical properties and topologies has yet to be explored. Here we have performed a systematic study of an arginine (Arg) side chain analog moving across saturated phosphatidylcholine (PC) bilayers of variable hydrocarbon tail length from 10 to 18 carbons. For all bilayers we observe similar ion-induced defects, where Arg draws water molecules and lipid head groups into the bilayers to avoid large dehydration energy costs. The free energy profiles all exhibit sharp climbs with increasing penetration into the hydrocarbon core, with predictable shifts between bilayers of different thickness, leading to barrier reduction from 26 kcal/mol for 18 carbons to 6 kcal/mol for 10 carbons. For lipids of 10 and 12 carbons we observe narrow transmembrane pores and corresponding plateaus in the free energy profiles. Allowing for movements of the protein and side chain snorkeling, we argue that the energetic cost for burying Arg inside a thin bilayer will be small, consistent with recent experiments, also leading to a dramatic reduction in pK(a) shifts for Arg. We provide evidence that Arg translocation occurs via an ion-induced defect mechanism, except in thick bilayers (of at least 18 carbons) where solubility-diffusion becomes energetically favored. Our findings shed light on the mechanisms of ion movement through membranes of varying composition, with implications for a range of charged protein-lipid interactions and the actions of cell-perturbing peptides. This article is part of a Special Issue entitled: Membrane protein structure and function.

Transfer of arginine into lipid bilayers is nonadditive

Computer simulations suggest that the translocation of arginine through the hydrocarbon core of a lipid membrane proceeds by the formation of a water-filled defect that keeps the arginine molecule hydrated even at the center of the bilayer. We show here that adding additional arginine molecules into one of these water defects causes only a small change in free energy. The barrier for transferring multiple arginines through the membrane is approximately the same as for a single arginine and may even be lower depending on the exact geometry of the system. We discuss these results in the context of arginine-rich peptides such as antimicrobial and cell-penetrating peptides.

Computational studies of membrane proteins: models and predictions for biological understanding

We discuss recent progresses in computational studies of membrane proteins based on physical models with parameters derived from bioinformatics analysis. We describe computational identification of membrane proteins and prediction of their topology from sequence, discovery of sequence and spatial motifs, and implications of these discoveries. The detection of evolutionary signal for understanding the substitution pattern of residues in the TM segments and for sequence alignment is also discussed. We further discuss empirical potential functions for energetics of inserting residues in the TM domain, for interactions between TM helices or strands, and their applications in predicting lipid-facing surfaces of the TM domain. Recent progresses in structure predictions of membrane proteins are also reviewed, with further discussions on calculation of ensemble properties such as melting temperature based on simplified state space model. Additional topics include prediction of oligomerization state of membrane proteins, identification of the interfaces for protein-protein interactions, and design of membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Hydrophobicity scales: a thermodynamic looking glass into lipid-protein interactions

The partitioning of amino acid sidechains into the membrane is a key aspect of membrane protein folding. However, lipid bilayers exhibit rapidly changing physicochemical properties over their nanometer-scale thickness, which complicates understanding the thermodynamics and microscopic details of membrane partitioning. Recent data from diverse approaches, including protein insertion by the Sec translocon, folding of a small beta-barrel membrane protein and computer simulations of the exact distribution of a variety of small molecules and peptides, have joined older hydrophobicity scales for membrane protein prediction. We examine the correlations among the scales and find that they are remarkably correlated even though there are large differences in magnitude. We discuss the implications of these scales for understanding membrane protein structure and function.

Efficient molecular mechanics simulations of the folding, orientation, and assembly of peptides in lipid bilayers using an implicit atomic solvation model

Membrane proteins comprise a significant fraction of the proteomes of sequenced organisms and are the targets of approximately half of marketed drugs. However, in spite of their prevalence and biomedical importance, relatively few experimental structures are available due to technical challenges. Computational simulations can potentially address this deficit by providing structural models of membrane proteins. Solvation within the spatially heterogeneous membrane/solvent environment provides a major component of the energetics driving protein folding and association within the membrane. We have developed an implicit solvation model for membranes that is both computationally efficient and accurate enough to enable molecular mechanics predictions for the folding and association of peptides within the membrane. We derived the new atomic solvation model parameters using an unbiased fitting procedure to experimental data and have applied it to diverse problems in order to test its accuracy and to gain insight into membrane protein folding. First, we predicted the positions and orientations of peptides and complexes within the lipid bilayer and compared the simulation results with solid-state NMR structures. Additionally, we performed folding simulations for a series of host-guest peptides with varying propensities to form alpha helices in a hydrophobic environment and compared the structures with experimental measurements. We were also able to successfully predict the structures of amphipathic peptides as well as the structures for dimeric complexes of short hexapeptides that have experimentally characterized propensities to form beta sheets within the membrane. Finally, we compared calculated relative transfer energies with data from experiments measuring the effects of mutations on the free energies of translocon-mediated insertion of proteins into lipid bilayers and of combined folding and membrane insertion of a beta barrel protein.

Exploring peptide-membrane interactions with coarse-grained MD simulations

The interaction of α-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an α-helix within a membrane, or that change the probability that the α-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an α-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.

Membrane proteins: from bench to bits

Membrane proteins currently receive a lot of attention, in large part thanks to a steady stream of high-resolution X-ray structures. Although the first few structures showed proteins composed of tightly packed bundles of very hydrophobic more or less straight transmembrane α-helices, we now know that helix-bundle membrane proteins can be both highly flexible and contain transmembrane segments that are neither very hydrophobic nor necessarily helical throughout their lengths. This raises questions regarding how membrane proteins are inserted into the membrane and fold in vivo, and also complicates life for bioinformaticians trying to predict membrane protein topology and structure.

Multiscale modeling of biological functions

Recent years have witnessed a tremendous explosion in computational power, which in turn has resulted in great progress in the complexity of the biological and chemical problems that can be addressed by means of all-atom simulations. Despite this, however, our computational time is not infinite, and in fact many of the key problems of the field were resolved long before the existence of the current levels of computational power. This review will start by presenting a brief historical overview of the use of multiscale simulations in biology, and then present some key developments in the field, highlighting several cases where the use of a physically sound simplification is clearly superior to a brute-force approach. Finally, some potential future directions will be discussed.

Membrane insertion of a voltage sensor helix

Most membrane proteins contain a transmembrane (TM) domain made up of a bundle of lipid-bilayer-spanning α-helices. TM α-helices are generally composed of a core of largely hydrophobic amino acids, with basic and aromatic amino acids at each end of the helix forming interactions with the lipid headgroups and water. In contrast, the S4 helix of ion channel voltage sensor (VS) domains contains four or five basic (largely arginine) side chains along its length and yet adopts a TM orientation as part of an independently stable VS domain. Multiscale molecular dynamics simulations are used to explore how a charged TM S4 α-helix may be stabilized in a lipid bilayer, which is of relevance in the context of mechanisms of translocon-mediated insertion of S4. Free-energy profiles for insertion of the S4 helix into a phospholipid bilayer suggest that it is thermodynamically favorable for S4 to insert from water to the center of the membrane, where the helix adopts a TM orientation. This is consistent with crystal structures of Kv channels, biophysical studies of isolated VS domains in lipid bilayers, and studies of translocon-mediated S4 helix insertion. Decomposition of the free-energy profiles reveals the underlying physical basis for TM stability, whereby the preference of the hydrophobic residues of S4 to enter the bilayer dominates over the free-energy penalty for inserting charged residues, accompanied by local distortion of the bilayer and penetration of waters. We show that the unique combination of charged and hydrophobic residues in S4 allows it to insert stably into the membrane.

Interaction of diverse voltage sensor homologs with lipid bilayers revealed by self-assembly simulations

Voltage sensors (VS) domains couple the activation of ion channels/enzymes to changes in membrane voltage. We used molecular dynamics simulations to examine interactions with lipids of several VS homologs. VSs in intact channels in the activated state are exposed to phospholipids, leading to a characteristic local distortion of the lipid bilayer which decreases its thickness by ∼10 Å. This effect is mediated by a conserved hydrophilic stretch in the S4-S5 segment linking the VS and the pore domains, and may favor gating charges crossing the membrane. In cationic lipid bilayers lacking phosphate groups, VSs form fewer contacts with lipid headgroups. The S3-S4 paddle motifs show persistent interactions of individual lipid molecules, influenced by the hairpin loop. In conclusion, our results suggest common interactions with phospholipids for various VS homologs, providing insights into the molecular basis of their stabilization in the membrane and how they are altered by lipid modification.

Free-energy cost for translocon-assisted insertion of membrane proteins

Nascent membrane proteins typically insert in a sequential fashion into the membrane via a protein-conducting channel, the Sec translocon. How this process occurs is still unclear, although a thermodynamic partitioning between the channel and the membrane environment has been proposed. Experiment- and simulation-based scales for the insertion free energy of various amino acids are, however, at variance, the former appearing to lie in a narrower range than the latter. Membrane insertion of arginine, for instance, requires 14-17 kcal/mol according to molecular dynamics simulations, but only 2-3 kcal/mol according to experiment. We suggest that this disagreement is resolved by assuming a two-stage insertion process wherein the first step, the insertion into the translocon, is energized by protein synthesis and, therefore, has an effectively zero free-energy cost; the second step, the insertion into the membrane, invokes the translocon as an intermediary between the fully hydrated and the fully inserted locations. Using free-energy perturbation calculations, the effective transfer free energies from the translocon to the membrane have been determined for both arginine and leucine amino acids carried by a background polyleucine helix. Indeed, the insertion penalty for arginine as well as the insertion gain for leucine from the translocon to the membrane is found to be significantly reduced compared to direct insertion from water, resulting in the same compression as observed in the experiment-based scale.

The energetics of transmembrane helix insertion into a lipid bilayer

Free energy profiles for insertion of a hydrophobic transmembrane protein α-helix (M2 from CFTR) into a lipid bilayer have been calculated using coarse-grained molecular dynamics simulations and umbrella sampling to yield potentials of mean force along a reaction path corresponding to translation of a helix across a lipid bilayer. The calculated free energy of insertion is smaller when a bilayer with a thinner hydrophobic region is used. The free energies of insertion from the potentials of mean force are compared with those derived from a number of hydrophobicity scales and with those derived from translocon-mediated insertion. This comparison supports recent models of translocon-mediated insertion and in particular suggests that: 1), helices in an about-to-be-inserted state may be located in a hydrophobic region somewhat thinner than the core of a lipid bilayer; and/or 2), helices in a not-to-be-inserted state may experience an environment more akin (e.g., in polarity/hydrophobicity) to the bilayer/water interface than to bulk water.