Reconciling the roles of kinetic and thermodynamic factors in membrane-protein insertion

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.

Folding and Insertion of Transmembrane Helices at the ER

In eukaryotic cells, the endoplasmic reticulum (ER) is the entry point for newly synthesized proteins that are subsequently distributed to organelles of the endomembrane system. Some of these proteins are completely translocated into the lumen of the ER while others integrate stretches of amino acids into the greasy 30 Å wide interior of the ER membrane bilayer. It is generally accepted that to exist in this non-aqueous environment the majority of membrane integrated amino acids are primarily non-polar/hydrophobic and adopt an α-helical conformation. These stretches are typically around 20 amino acids long and are known as transmembrane (TM) helices. In this review, we will consider how transmembrane helices achieve membrane integration. We will address questions such as: Where do the stretches of amino acids fold into a helical conformation? What is/are the route/routes that these stretches take from synthesis at the ribosome to integration through the ER translocon? How do these stretches ‘know’ to integrate and in which orientation? How do marginally hydrophobic stretches of amino acids integrate and survive as transmembrane helices?

Molecular Modeling of Signal Peptide Recognition by Eukaryotic Sec Complexes

Here, we review recent molecular modelling and simulation studies of the Sec translocon, the primary component/channel of protein translocation into the endoplasmic reticulum (ER) and bacterial periplasm, respectively. Our focus is placed on the eukaryotic Sec61, but we also mention modelling studies on prokaryotic SecY since both systems operate in related ways. Cryo-EM structures are now available for different conformational states of the Sec61 complex, ranging from the idle or closed state over an inhibited state with the inhibitor mycolactone bound near the lateral gate, up to a translocating state with bound substrate peptide in the translocation pore. For all these states, computational studies have addressed the conformational dynamics of the translocon with respect to the pore ring, the plug region, and the lateral gate. Also, molecular simulations are addressing mechanistic issues of insertion into the ER membrane vs. translocation into the ER, how signal-peptides are recognised at all in the translocation pore, and how accessory proteins affect the Sec61 conformation in the co- and post-translational pathways.

The structure of the TOM core complex in the mitochondrial outer membrane

In the past three decades, significant advances have been made in providing the biochemical background of TOM (translocase of the outer mitochondrial membrane)-mediated protein translocation into mitochondria. In the light of recent cryoelectron microscopy-derived structures of TOM isolated from Neurospora crassa and Saccharomyces cerevisiae, the interpretation of biochemical and biophysical studies of TOM-mediated protein transport into mitochondria now rests on a solid basis. In this review, we compare the subnanometer structure of N. crassa TOM core complex with that of yeast. Both structures reveal remarkably well-conserved symmetrical dimers of 10 membrane protein subunits. The structural data also validate predictions of weakly stable regions in the transmembrane β-barrel domains of the protein-conducting subunit Tom40, which signal the existence of β-strands located in interfaces of protein-protein interactions.

Extension of the Variational Free Energy Profile and Multistate Bennett Acceptance Ratio Methods for High-Dimensional Potential of Mean Force Profile Analysis

We redevelop the variational free energy profile (vFEP) method using a cardinal B-spline basis to extend the method for analyzing free energy surfaces (FESs) involving three or more reaction coordinates. We also implemented software for evaluating high-dimensional profiles based on the multistate Bennett acceptance ratio (MBAR) method which constructs an unbiased probability density from global reweighting of the observed samples. The MBAR method takes advantage of a fast algorithm for solving the unbinned weighted histogram (UWHAM)/MBAR equations which replaces the solution of simultaneous equations with a nonlinear optimization of a convex function. We make use of cardinal B-splines and multiquadric radial basis functions to obtain smooth, differentiable MBAR profiles in arbitrary high dimensions. The cardinal B-spline vFEP and MBAR methods are compared using three example systems that examine 1D, 2D, and 3D profiles. Both methods are found to be useful and produce nearly indistinguishable results. The vFEP method is found to be 150 times faster than MBAR when applied to periodic 2D profiles, but the MBAR method is 4.5 times faster than vFEP when evaluating unbounded 3D profiles. In agreement with previous comparisons, we find the vFEP method produces superior FESs when the overlap between umbrella window simulations decreases. Finally, the associative reaction mechanism of hammerhead ribozyme is characterized using 3D, 4D, and 6D profiles, and the higher-dimensional profiles are found to have smaller reaction barriers by as much as 1.5 kcal/mol. The methods presented here have been implemented into the FE-ToolKit software package along with new methods for network-wide free energy analysis in drug discovery.

Characterization of the Features of Water Inside the SecY Translocon

Despite extended experimental and computational studies, the mechanism regulating membrane protein folding and stability in cell membranes is not fully understood. In this review, I will provide a personal and partial account of the scientific efforts undertaken by Dr. Stephen White to shed light on this topic. After briefly describing the role of water and the hydrophobic effect on cellular processes, I will discuss the physical chemistry of water confined inside the SecY translocon pore. I conclude with a review of recent literature that attempts to answer fundamental questions on the pathway and energetics of translocon-guided membrane protein insertion.

Recent Insights from Molecular Dynamics Simulations for G Protein-Coupled Receptor Drug Discovery

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. Some ligands that bind to GPCRs activate different downstream signaling pathways. G protein activation, or -arrestin biased signaling, involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. -arrestin biased signaling has become a hot target for structure-based drug discovery. However, challenges include that there are few crystal structures available in the Protein Data Bank and that GPCRs are highly dynamic. Hence, molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions with receptors and ligands. Here, we highlight recent MD simulation studies and enhanced sampling methods used to study biased G protein-coupled receptor signaling and their conformational dynamics as well as applications to drug discovery.

Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis

Advances over the past 25 years have revealed much about how the structural properties of membranes and associated proteins are linked to the thermodynamics and kinetics of membrane protein (MP) folding. At the same time biochemical progress has outlined how cellular proteostasis networks mediate MP folding and manage misfolding in the cell. When combined with results from genomic sequencing, these studies have established paradigms for how MP folding and misfolding are linked to the molecular etiologies of a variety of diseases. This emerging framework has paved the way for the development of a new class of small molecule "pharmacological chaperones" that bind to and stabilize misfolded MP variants, some of which are now in clinical use. In this review, we comprehensively outline current perspectives on the folding and misfolding of integral MPs as well as the mechanisms of cellular MP quality control. Based on these perspectives, we highlight new opportunities for innovations that bridge our molecular understanding of the energetics of MP folding with the nuanced complexity of biological systems. Given the many linkages between MP misfolding and human disease, we also examine some of the exciting opportunities to leverage these advances to address emerging challenges in the development of therapeutics and precision medicine.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

A molecular dynamics study on the resilience of Sec61 channel from open to closed state

When the nascent chain is released from the ribosome, its packing into the apolar environment of the lipid bilayer in the endoplasmic reticulum is facilitated by the Sec61 translocon. In this process, coupling of the conformational change of the channel is essential to transport the nascent chain and meanwhile maintain the membrane permeability barrier. Two molecular dynamics simulations were performed in the current work to investigate the resilience of the lateral gate and the linkage mechanism of the lateral gate, pore ring and plug. The results affirmed that the lateral gate is able to recover its partially-closed state rapidly after the nascent chain segment enters the bilayer. This triggers subsequent motions of the pore ring and plug, which prevent the small molecules passing through the pore. The pore diameter in the partially-closed state is about 6–7 Å. The plug would move upward ∼2 Å if the lateral gate could not close. Two waters permeate through the channel when the lateral gate was open. Water molecules could go across the bilayer via the gap of the open lateral gate due to the occluding of the pore ring and plug.

The lateral gate of Sec61 is able to recover its partially-closed state rapidly after the nascent chain segment enters the bilayer, which triggers subsequent motions of the pore ring and plug.

Interaction mapping of the Sec61 translocon identifies two Sec61α regions interacting with hydrophobic segments in translocating chains

Many proteins in organelles of the secretory pathway, as well as secretory proteins, are translocated across and inserted into the endoplasmic reticulum membrane by the Sec61 translocon, a protein-conducting channel. The channel consists of 10 transmembrane (TM) segments of the Sec61α subunit and possesses an opening between TM2b and TM7, termed the lateral gate. Structural and biochemical analyses of complexes of Sec61 and its ortholog SecY have revealed that the lateral gate is the exit for signal sequences and TM segments of translocating polypeptides to the lipid bilayer and also involved in the recognition of such hydrophobic sequences. Moreover, even marginally hydrophobic (mH) segments insufficient for membrane integration can be transiently stalled in surrounding Sec61α regions and cross-linked to them, but how the Sec61 translocon accommodates these mH segments remains unclear. Here, we used Cys-scanned variants of human Sec61α expressed in cultured 293-H cells to examine which channel regions associate with mH segments. A TM segment in a ribosome-associated polypeptide was mainly cross-linked to positions at the lateral gate, whereas an mH segment in a nascent chain was cross-linked to the Sec61α pore-interior positions at TM5 and TM10, as well as the lateral gate. Of note, cross-linking at position 180 in TM5 of Sec61α was reduced by an I179A substitution. We therefore conclude that at least two Sec61α regions, the lateral gate and the pore-interior site around TM5, interact with mH segments and are involved in accommodating them.

Simulations of lipid bilayers using the CHARMM36 force field with the TIP3P-FB and TIP4P-FB water models

The CHARMM36 force field for lipids is widely used in simulations of lipid bilayers. The CHARMM family of force fields were developed for use with the mTIP3P water model. This water model has an anomalously high dielectric constant and low viscosity, which limits its accuracy in the calculation of quantities like permeability coefficients. The TIP3P-FB and TIP4P-FB water models are more accurate in terms of the dielectric constant and transport properties, which could allow more accurate simulations of systems containing water and lipids. To test whether the CHARMM36 lipid force field is compatible with the TIP3P-FB and TIP4P-FB water models, we have performed simulations of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers. The calculated headgroup area, compressibility, order parameters, and X-ray form factors are in good agreement with the experimental values, indicating that these improved water models can be used with the CHARMM36 lipid force field without modification when calculating membrane physical properties. The water permeability predicted by these models is significantly different; the mTIP3P-model diffusion in solution and at the lipid-water interface is anomalously fast due to the spuriously low viscosity of mTIP3P-model water, but the potential of mean force of permeation is higher for the TIP3P-FB and TIP4P-FB models due to their high excess chemical potentials. As a result, the rates of water permeation calculated the FB water models are slower than the experimental value by a factor of 15-17, while simulations with the mTIP3P model only underestimate the water permeability by a factor of 3.

Cotranslational protein targeting to the membrane: Nascent-chain transfer in a quaternary complex formed at the translocon

Membrane proteins in bacteria are cotranslationally inserted into the plasma membrane through the SecYEG translocon. Ribosomes exposing the signal-anchor sequence (SAS) of a membrane protein are targeted to the translocon by the signal recognition particle (SRP) pathway. SRP scans translating ribosomes and forms high-affinity targeting complexes with those exposing a SAS. Recognition of the SAS activates SRP for binding to its receptor, FtsY, which, in turn, is primed for SRP binding by complex formation with SecYEG, resulting in a quaternary targeting complex. Here we examine the effect of SecYEG docking to ribosome-nascent-chain complexes (RNCs) on SRP binding and SAS transfer, using SecYEG embedded in phospholipid-containing nanodiscs and monitoring FRET between fluorescence-labeled constituents of the targeting complex. SecYEG–FtsY binding to RNC–SRP complexes lowers the affinity of SRP to both ribosome and FtsY, indicating a general weakening of the complex due to partial binding competition near the ribosomal peptide exit. The rearrangement of the quaternary targeting complex to the pre-transfer complex requires an at least partially exposed SAS. The presence of SecYEG-bound FtsY and the length of the nascent chain strongly influence nascent-chain transfer from SRP to the translocon and repositioning of SRP in the post-transfer complex.

Computed Free Energies of Peptide Insertion into Bilayers are Independent of Computational Method

We show that the free energy of inserting hydrophobic peptides into lipid bilayer membranes from surface-aligned to transmembrane inserted states can be reliably calculated using atomistic models. We use two entirely different computational methods: high temperature spontaneous peptide insertion calculations as well as umbrella sampling potential-of-mean-force (PMF) calculations, both yielding the same energetic profiles. The insertion free energies were calculated using two different protein and lipid force fields (OPLS protein/united-atom lipids and CHARMM36 protein/all-atom lipids) and found to be independent of the simulation parameters. In addition, the free energy of insertion is found to be independent of temperature for both force fields. However, we find major difference in the partitioning kinetics between OPLS and CHARMM36, likely due to the difference in roughness of the underlying free energy surfaces. Our results demonstrate not only a reliable method to calculate insertion free energies for peptides, but also represent a rare case where equilibrium simulations and PMF calculations can be directly compared.

Mechanism of the Spontaneous and Directional Membrane Insertion of a 2-Transmembrane Ion Channel

Protein insertion into membranes is a process occurring in every cell and every cellular compartment. Yet, many thermodynamic aspects of this fundamental biophysical process are not well understood. We investigated physicochemical parameters that influence protein insertion using the model protein KcsA, a 2-transmembrane ion channel. To understand what drives insertion and to identify individual steps of protein integration into a highly apolar environment, we investigated the contribution of electrostatic interactions and lipid composition on protein insertion on a single molecule level. We show that insertion of KcsA is spontaneous and directional as the cytosolic part of the protein does not translocate across the membrane barrier. Surprisingly, not hydrophobic residues but charged amino acids are crucial for the insertion of the unfolded protein into the membrane. Our results demonstrate the importance of electrostatic interactions between membrane and protein during the insertion process of hydrophobic polypeptides into the apolar membrane. On the basis of the observation that negatively charged lipids increase insertion events while high ionic strength in the surrounding aqueous phase decreases insertion events, a two-step mechanism is proposed. Here, an initial electrostatic attraction between membrane and protein represents the first step prior to insertion of hydrophobic residues into the hydrocarbon core of the membrane.

Insights into Cotranslational Nascent Protein Behavior from Computer Simulations

Regulation of protein stability and function in vivo begins during protein synthesis, when the ribosome translates a messenger RNA into a nascent polypeptide. Cotranslational processes involving a nascent protein include folding, binding to other macromolecules, enzymatic modification, and secretion through membranes. Experiments have shown that the rate at which the ribosome adds amino acids to the elongating nascent chain influences the efficiency of these processes, with alterations to these rates possibly contributing to diseases, including some types of cancer. In this review, we discuss recent insights into cotranslational processes gained from molecular simulations, how different computational approaches have been combined to understand cotranslational processes at multiple scales, and the new scenarios illuminated by these simulations. We conclude by suggesting interesting questions that computational approaches in this research area can address over the next few years.

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.

Molecular dynamics simulations of biological membranes and membrane proteins using enhanced conformational sampling algorithms

This paper reviews various enhanced conformational sampling methods and explicit/implicit solvent/membrane models, as well as their recent applications to the exploration of the structure and dynamics of membranes and membrane proteins. Molecular dynamics simulations have become an essential tool to investigate biological problems, and their success relies on proper molecular models together with efficient conformational sampling methods. The implicit representation of solvent/membrane environments is reasonable approximation to the explicit all-atom models, considering the balance between computational cost and simulation accuracy. Implicit models can be easily combined with replica-exchange molecular dynamics methods to explore a wider conformational space of a protein. Other molecular models and enhanced conformational sampling methods are also briefly discussed. As application examples, we introduce recent simulation studies of glycophorin A, phospholamban, amyloid precursor protein, and mixed lipid bilayers and discuss the accuracy and efficiency of each simulation model and method. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Integration of transmembrane domains is regulated by their downstream sequences

The Sec61 translocon catalyzes translocation of proteins into the endoplasmic reticulum and the lateral integration of transmembrane segments into the lipid bilayer. Integration is mediated by the hydrophobicity of a polypeptide segment consistent with thermodynamic equilibration between the translocon and the lipid membrane. Integration efficiency of a generic series of increasingly hydrophobic sequences (H-segments) was found to diverge significantly in different reporter constructs as a function of the ∼100 residues carboxyterminal of the H-segments. The hydrophobicity threshold of integration was considerably lowered by insertion of generic ∼20-residue peptides either made of flexible glycine-serine repeats, containing multiple negative charges, or consisting of an oligo-proline stretch. A highly flexible, 100-residue glycine-serine stretch maximally enhanced this effect. The apparent free energy of integration was found to be changed by more than 3 kcal/mol with the downstream sequences tested. The C-terminal sequences could also be shown to affect integration of natural mildly hydrophobic sequences. The results suggest that the conformation of the nascent polypeptide in the protected cavity between ribosome and translocon significantly influences the release of the H-segment into the bilayer.

How transmembrane peptides insert and orientate in biomembranes: a combined experimental and simulation study

After the synthesis of transmembrane peptides/proteins (TMPs), their insertion into a lipid bilayer is a fundamental biophysical process.

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.

Structure of the native Sec61 protein-conducting channel

In mammalian cells, secretory and membrane proteins are translocated across or inserted into the endoplasmic reticulum (ER) membrane by the universally conserved protein-conducting channel Sec61, which has been structurally studied in isolated, detergent-solubilized states. Here we structurally and functionally characterize native, non-solubilized ribosome-Sec61 complexes on rough ER vesicles using cryo-electron tomography and ribosome profiling. Surprisingly, the 9-Å resolution subtomogram average reveals Sec61 in a laterally open conformation, even though the channel is not in the process of inserting membrane proteins into the lipid bilayer. In contrast to recent mechanistic models for polypeptide translocation and insertion, our results indicate that the laterally open conformation of Sec61 is the only conformation present in the ribosome-bound translocon complex, independent of its functional state. Consistent with earlier functional studies, our structure suggests that the ribosome alone, even without a nascent chain, is sufficient for lateral opening of Sec61 in a lipid environment.

The protein-conducting channel Sec61 is responsible for protein transport and membrane insertion at the endoplasmic reticulum. Here, the authors determine the structure of ribosome-bound Sec61 in a native context, in which it adopts a laterally open conformation, irrespective of its functional state.

Anomalous behavior of water inside the SecY translocon

The heterotrimeric SecY translocon complex is required for the cotranslational assembly of membrane proteins in bacteria and archaea. The insertion of transmembrane (TM) segments during nascent-chain passage through the translocon is generally viewed as a simple partitioning process between the water-filled translocon and membrane lipid bilayer, suggesting that partitioning is driven by the hydrophobic effect. Indeed, the apparent free energy of partitioning of unnatural aliphatic amino acids on TM segments is proportional to accessible surface area, which is a hallmark of the hydrophobic effect [Öjemalm K, et al. (2011) Proc Natl Acad Sci USA 108(31):E359-E364]. However, the apparent partitioning solvation parameter is less than one-half the value expected for simple bulk partitioning, suggesting that the water in the translocon departs from bulk behavior. To examine the state of water in a SecY translocon complex embedded in a lipid bilayer, we carried out all-atom molecular-dynamics simulations of the Pyrococcus furiosus SecYE, which was determined to be in a "primed" open state [Egea PF, Stroud RM (2010) Proc Natl Acad Sci USA 107(40):17182-17187]. Remarkably, SecYE remained in this state throughout our 450-ns simulation. Water molecules within SecY exhibited anomalous diffusion, had highly retarded rotational dynamics, and aligned their dipoles along the SecY transmembrane axis. The translocon is therefore not a simple water-filled pore, which raises the question of how anomalous water behavior affects the mechanism of translocon function and, more generally, the partitioning of hydrophobic molecules. Because large water-filled cavities are found in many membrane proteins, our findings may have broader implications.

Timing is everything: unifying codon translation rates and nascent proteome behavior

Experiments have demonstrated that changing the rate at which the ribosome translates a codon position in an mRNA molecule's open reading frame can alter the behavior of the newly synthesized protein. That is, codon translation rates can govern nascent proteome behavior. We emphasize that this phenomenon is a manifestation of the nonequilibrium nature of cotranslational processes, and as such, there exist theoretical tools that offer a potential means to quantitatively predict the influence of codon translation rates on the broad spectrum of nascent protein behaviors including cotranslational folding, aggregation, and translocation. We provide a review of the experimental evidence for the impact that codon translation rates can have, followed by a discussion of theoretical methods that can describe this phenomenon. The development and application of these tools are likely to provide fundamental insights into protein maturation and homeostasis, codon usage bias in organisms, the origins of translation related diseases, and new rational design methods for biotechnology and biopharmaceutical applications.

Lateral opening of the bacterial translocon on ribosome binding and signal peptide insertion

Proteins are co-translationally inserted into the bacterial plasma membrane via the SecYEG translocon by lateral release of hydrophobic transmembrane segments into the phospholipid bilayer. The trigger for lateral opening of the translocon is not known. Here we monitor lateral opening by photo-induced electron transfer (PET) between two fluorophores attached to the two SecY helices at the rim of the gate. In the resting translocon, the fluorescence is quenched, consistent with a closed conformation. Ribosome binding to the translocon diminishes PET quenching, indicating opening of the gate. The effect is larger with ribosomes exposing hydrophobic transmembrane segments and vanishes at low temperature. We propose a temperature-dependent dynamic equilibrium between closed and open conformations of the translocon that is shifted towards partially and fully open by ribosome binding and insertion of a hydrophobic peptide, respectively. The combined effects of ribosome and peptide binding allow for co-translational membrane insertion of successive transmembrane segments.

Integral membrane proteins laterally partition from the SecYEG translocon into the phospholipid bilayer. Here, the authors use photo-induced electron transfer to show that ribosome binding induces the opening of the lateral gate, and demonstrate that lateral opening does not happen at low temperature.

Lateral opening and exit pore formation are required for BamA function

The outer membrane of Gram-negative bacteria is replete with a host of β-barrel outer membrane proteins (OMPs). Despite serving a variety of essential functions, including immune response evasion, the exact mechanism of OMP folding and membrane insertion remains largely unclear. The β-barrel assembly machinery complex is required for OMP biogenesis. Crystal structures and molecular dynamics (MD) simulations of the central and essential component, BamA, suggest a mechanism involving lateral opening of its barrel domain. MD simulations reported here reveal an additional feature of BamA: a substrate exit pore positioned above the lateral opening site. Disulfide crosslinks that prevent lateral opening and exit pore formation result in a loss of BamA function, which can be fully rescued by the reductant tris(2-carboxyethyl)phosphine. These data provide strong evidence that lateral opening and exit pore formation are required for BamA function.

A computational kinetic model of diffusion for molecular systems

Regulation of biomolecular transport in cells involves intra-protein steps like gating and passage through channels, but these steps are preceded by extra-protein steps, namely, diffusive approach and admittance of solutes. The extra-protein steps develop over a 10-100 nm length scale typically in a highly particular environment, characterized through the protein's geometry, surrounding electrostatic field, and location. In order to account for solute energetics and mobility of solutes in this environment at a relevant resolution, we propose a particle-based kinetic model of diffusion based on a Markov State Model framework. Prerequisite input data consist of diffusion coefficient and potential of mean force maps generated from extensive molecular dynamics simulations of proteins and their environment that sample multi-nanosecond durations. The suggested diffusion model can describe transport processes beyond microsecond duration, relevant for biological function and beyond the realm of molecular dynamics simulation. For this purpose the systems are represented by a discrete set of states specified by the positions, volumes, and surface elements of Voronoi grid cells distributed according to a density function resolving the often intricate relevant diffusion space. Validation tests carried out for generic diffusion spaces show that the model and the associated Brownian motion algorithm are viable over a large range of parameter values such as time step, diffusion coefficient, and grid density. A concrete application of the method is demonstrated for ion diffusion around and through the Eschericia coli mechanosensitive channel of small conductance ecMscS.

Surface-Tension Replica-Exchange Molecular Dynamics Method for Enhanced Sampling of Biological Membrane Systems

Conformational sampling is fundamentally important for simulating complex biomolecular systems. The generalized-ensemble algorithm, especially the temperature replica-exchange molecular dynamics method (T-REMD), is one of the most powerful methods to explore structures of biomolecules such as proteins, nucleic acids, carbohydrates, and also of lipid membranes. T-REMD simulations have focused on soluble proteins rather than membrane proteins or lipid bilayers, because explicit membranes do not keep their structural integrity at high temperature. Here, we propose a new generalized-ensemble algorithm for membrane systems, which we call the surface-tension REMD method. Each replica is simulated in the NPγT ensemble, and surface tensions in a pair of replicas are exchanged at certain intervals to enhance conformational sampling of the target membrane system. We test the method on two biological membrane systems: a fully hydrated DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine) lipid bilayer and a WALP23-POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane system. During these simulations, a random walk in surface tension space is realized. Large-scale lateral deformation (shrinking and stretching) of the membranes takes place in all of the replicas without collapse of the lipid bilayer structure. There is accelerated lateral diffusion of DPPC lipid molecules compared with conventional MD simulation, and a much wider range of tilt angle of the WALP23 peptide is sampled due to large deformation of the POPC lipid bilayer and through peptide-lipid interactions. Our method could be applicable to a wide variety of biological membrane systems.

Functional asymmetry within the Sec61p translocon

The Sec61 translocon forms a pore to translocate polypeptide sequences across the membrane and offers a lateral gate for membrane integration of hydrophobic (H) segments. A central constriction of six apolar residues has been shown to form a seal, but also to determine the hydrophobicity threshold for membrane integration: Mutation of these residues in yeast Sec61p to glycines, serines, aspartates, or lysines lowered the hydrophobicity required for integration; mutation to alanines increased it. Whereas four leucines distributed in an oligo-alanine H segment were sufficient for 50% integration, we now find four leucines in the N-terminal half of the H segment to produce significantly more integration than in the C-terminal half, suggesting functional asymmetry within the translocon. Scanning a cluster of three leucines through an oligo-alanine H segment showed high integration levels, except around the position matching that of the hydrophobic constriction in the pore where integration was strongly reduced. Both asymmetry and the position effect of H-segment integration disappeared upon mutation of the constriction residues to glycines or serines, demonstrating that hydrophobicity at this position within the translocon is responsible for the phenomenon. Asymmetry was largely retained, however, when constriction residues were replaced by alanines. These results reflect on the integration mechanism of transmembrane domains and show that membrane insertion of H segments strongly depends not only on their intrinsic hydrophobicity but also on the local conditions in the translocon interior. Thus, the contribution of hydrophobic residues in the H segment is not simply additive and displays cooperativeness depending on their relative position.

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.

In vitro reconstitution of lipid-dependent dual topology and postassembly topological switching of a membrane protein

Phospholipids could exert their effect on membrane protein topology either directly by interacting with topogenic signals of newly inserted proteins or indirectly by influencing the protein assembly machinery. In vivo lactose permease (LacY) of Escherichia coli displays a mixture of topological conformations ranging from complete inversion of the N-terminal helical bundle to mixed topology and then to completely native topology as phosphatidylethanolamine (PE) is increased from 0% to 70% of membrane phospholipids. These topological conformers are interconvertible by postassembly synthesis or dilution of PE in vivo. To investigate whether coexistence of multiple topological conformers is dependent solely on the membrane lipid composition, we determined the topological organization of LacY in an in vitro proteoliposome system in which lipid composition can be systematically controlled before (liposomes) and after (fliposomes) reconstitution using a lipid exchange technique. Purified LacY reconstituted into preformed liposomes of increasing PE content displayed inverted topology at low PE and then a mixture of inverted and proper topologies with the latter increasing with increasing PE until all LacY adopted its native topology. Interconversion between topological conformers of LacY was observed in a PE dose-dependent manner by either increasing or decreasing PE levels in proteoliposomes postreconstitution of LacY, clearly demonstrating that membrane protein topology can be changed simply by changing membrane lipid composition independent of other cellular factors. The results provide a thermodynamic-based lipid-dependent model for shifting the equilibrium between different conformational states of a membrane protein.

The bacterial translocon SecYEG opens upon ribosome binding

In co-translational translocation, the ribosome funnel and the channel of the protein translocation complex SecYEG are aligned. For the nascent chain to enter the channel immediately after synthesis, a yet unidentified signal triggers displacement of the SecYEG sealing plug from the pore. Here, we show that ribosome binding to the resting SecYEG channel triggers this conformational transition. The purified and reconstituted SecYEG channel opens to form a large ion-conducting channel, which has the conductivity of the plug deletion mutant. The number of ion-conducting channels inserted into the planar bilayer per fusion event roughly equals the number of SecYEG channels counted by fluorescence correlation spectroscopy in a single proteoliposome. Thus, the open probability of the channel must be close to unity. To prevent the otherwise lethal proton leak, a closed post-translational conformation of the SecYEG complex bound to a ribosome must exist.