Side-chain Contributions to Membrane Protein Structure and Stability

Effects of stand-alone polar residue on membrane protein stability and structure

Helical membrane proteins generally have a hydrophobic nature, with apolar side chains comprising the majority of the transmembrane (TM) helices. However, whenever polar side chains are present in the TM domain, they often exert a crucial role in structural interactions with other polar residues, such as TM helix associations and oligomerization. Moreover, polar residues in the TM region also often participate in protein functions, such as the Schiff base bonding between Lys residues and retinal in rhodopsin-like membrane proteins. Although many studies have focused on these functional polar residues, our understanding of stand-alone polar residues that are energetically unfavored in TM helixes is limited. Here, we adopted bacteriorhodopsin (bR) as a model system and systematically mutated 17 of its apolar Leu or Phe residues to polar Asn. Stability measurements of the resulting mutants revealed that all of these polar substitutions reduced bR stability to various extents, and the extent of destabilization of each mutant bR is also correlated to different structural factors, such as the relative accessible surface area and membrane depth of the mutation site. Structural analyses of these Asn residues revealed that they form sidechain-to-backbone hydrogen bonds that alleviate the unfavorable energetics in hydrophobic and apolar surroundings. Our results indicate that membrane proteins are able to accommodate certain stand-alone polar residues in the TM region without disrupting overall structures.

Steric trapping strategy for studying the folding of helical membrane proteins

Elucidating the folding energy landscape of membrane proteins is essential to the understanding of the proteins' stabilizing forces, folding mechanisms, biogenesis, and quality control. This is not a trivial task because the reversible control of folding is inherently difficult in a lipid bilayer environment. Recently, novel methods have been developed, each of which has a unique strength in investigating specific aspects of membrane protein folding. Among such methods, steric trapping is a versatile strategy allowing a reversible control of membrane protein folding with minimal perturbation of native protein-water and protein-lipid interactions. In a nutshell, steric trapping exploits the coupling of spontaneous denaturation of a doubly biotinylated protein to the simultaneous binding of bulky monovalent streptavidin molecules. This strategy has been evolved to investigate key elements of membrane protein folding such as thermodynamic stability, spontaneous denaturation rates, conformational features of the denatured states, and cooperativity of stabilizing interactions. In this review, we describe the critical methodological advancement, limitation, and outlook of the steric trapping strategy.

Gold nanoparticle-powered screening of membrane protein-specific lipids from complex lipid mixtures

Membrane proteins (MPs) are affected by binding of specific lipids. We previously developed a methodology for systematically analyzing MP-lipid interactions leveraging surface plasmon resonance (SPR). In this method, the gold sensor chip surface was modified with a self-assembled monolayer (SAM), which allowed for a larger amount of MP-immobilization. However, the laborious lipid purification step remained a bottleneck. To address this issue, a new strategy has been developed utilizing gold nanoparticles (AuNPs) instead of the gold sensor chip. AuNPs were coated with SAM, on which MP was covalently anchored. The MP-immobilized AuNPs were mixed with a lipid mixture, and the recovered lipids were quantified by LC-MS. Bacteriorhodopsin (bR) was used as an MP to demonstrate this concept. We optimized immobilization conditions and confirmed the efficient immobilization of bR by dynamic light scattering and electron micrographs. Washing conditions for pulldown experiments were optimized to efficiently remove non-specific lipids. A new binding index was introduced to qualitatively reproduce the known affinity of lipids for bR. Consequently, the low-abundant and least-studied lipid S-TeGD was identified as a candidate for bR-specific lipids. This technique can skip the laborious lipid purification process, accelerating the screening of MP-specific lipids from complex lipid mixtures.

Quantifying a light-induced energetic change in bacteriorhodopsin by force spectroscopy

Ligand-induced conformational changes are critical to the function of many membrane proteins and arise from numerous intramolecular interactions. In the photocycle of the model membrane protein bacteriorhodopsin (bR), absorption of a photon by retinal triggers a conformational cascade that results in pumping a proton across the cell membrane. While decades of spectroscopy and structural studies have probed this photocycle in intricate detail, changes in intramolecular energetics that underlie protein motions have remained elusive to experimental quantification. Here, we measured these energetics on the millisecond time scale using atomic-force-microscopy-based single-molecule force spectroscopy. Precisely, timed light pulses triggered the bR photocycle while we measured the equilibrium unfolding and refolding of the terminal 8-amino-acid region of bR's G-helix. These dynamics changed when the EF-helix pair moved ~9 Å away from this end of the G helix during the "open" portion of bR's photocycle. In ~60% of the data, we observed abrupt light-induced destabilization of 3.4 ± 0.3 kcal/mol, lasting 38 ± 3 ms. The kinetics and pH-dependence of this destabilization were consistent with prior measurements of bR's open phase. The frequency of light-induced destabilization increased with the duration of illumination and was dramatically reduced in the triple mutant (D96G/F171C/F219L) thought to trap bR in its open phase. In the other ~40% of the data, photoexcitation unexpectedly stabilized a longer-lived putative misfolded state. Through this work, we establish a general single-molecule force spectroscopy approach for measuring ligand-induced energetics and lifetimes in membrane proteins.

Sensitivity optimization of a rhodopsin-based fluorescent voltage indicator

The ability to optically image cellular transmembrane voltages at millisecond-timescale resolutions can offer unprecedented insight into the function of living brains in behaving animals. Here, we present a point mutation that increases the sensitivity of Ace2 opsin-based voltage indicators. We use the mutation to develop Voltron2, an improved chemigeneic voltage indicator that has a 65% higher sensitivity to single APs and 3-fold higher sensitivity to subthreshold potentials than Voltron. Voltron2 retained the sub-millisecond kinetics and photostability of its predecessor, although with lower baseline fluorescence. In multiple in vitro and in vivo comparisons with its predecessor across multiple species, we found Voltron2 to be more sensitive to APs and subthreshold fluctuations. Finally, we used Voltron2 to study and evaluate the possible mechanisms of interneuron synchronization in the mouse hippocampus. Overall, we have discovered a generalizable mutation that significantly increases the sensitivity of Ace2 rhodopsin-based sensors, improving their voltage reporting capability.

How do surfactants unfold and refold proteins?

Although the anionic surfactant sodium dodecyl sulfate, SDS, has been used for more than half a century as a versatile and efficient protein denaturant for protein separation and size estimation, there is still controversy about its mode of interaction with proteins. The term "rod-like" structures for the complexes that form between SDS and protein, originally introduced by Tanford, is not sufficiently descriptive and does not distinguish between the two current vying models, namely protein-decorated micelles a.k.a. the core-shell model (in which denatured protein covers the surface of micelles) versus beads-on-a-string model (where unfolded proteins are surrounded by surfactant micelles). Thanks to a combination of structural, kinetic and computational work particularly within the last 5-10 years, it is now possible to rule decisively in favor of the core-shell model. This is supported unambiguously by a combination of calorimetric and small-angle X-ray scattering (SAXS) techniques and confirmed by increasingly sophisticated molecular dynamics simulations. Depending on the SDS:protein ratio and the protein molecular mass, the formed structures can range from multiple partly unfolded protein molecules surrounding a single shared micelle to a single polypeptide chain decorating multiple micelles. We also have much new insight into how this species forms. It is preceded by the binding of small numbers of SDS molecules which subsequently grow by accretion. Time-resolved SAXS analysis reveals an asymmetric attack by SDS micelles followed by distribution of the increasingly unfolded protein around the micelle. The compactness of the protein chain continues to evolve at higher SDS concentrations according to single-molecule studies, though the protein remains completely denatured on the tertiary structural level. SDS denaturation can be reversed by addition of nonionic surfactants that absorb SDS forming mixed micelles, leaving the protein free to refold. Refolding can occur in parallel tracks if only a fraction of the protein is initially stripped of SDS. SDS unfolding is nearly always reversible unless carried out at low pH, where charge neutralization can lead to superclusters of protein-surfactant complexes. With the general mechanism of SDS denaturation now firmly established, it largely remains to explore how other ionic surfactants (including biosurfactants) may diverge from this path.

De novo-designed transmembrane domains tune engineered receptor functions

De novo-designed receptor transmembrane domains (TMDs) present opportunities for precise control of cellular receptor functions. We developed a de novo design strategy for generating programmed membrane proteins (proMPs): single-pass α-helical TMDs that self-assemble through computationally defined and crystallographically validated interfaces. We used these proMPs to program specific oligomeric interactions into a chimeric antigen receptor (CAR) that we expressed in mouse primary T cells and found that both in vitro CAR T cell cytokine release and in vivo antitumor activity scaled linearly with the oligomeric state encoded by the receptor TMD, from monomers up to tetramers. All programmed CARs stimulated substantially lower T cell cytokine release relative to the commonly used CD28 TMD, which we show elevated cytokine release through lateral recruitment of the endogenous T cell costimulatory receptor CD28. Precise design using orthogonal and modular TMDs thus provides a new way to program receptor structure and predictably tune activity for basic or applied synthetic biology.

Universal platform for the generation of thermostabilized GPCRs that crystallize in LCP

Structural studies of G-protein-coupled receptors (GPCRs) are often limited by difficulties in obtaining well-diffracting crystals suitable for high-resolution structure determination. During the past decade, crystallization in lipidic cubic phase (LCP) has become the most successful and widely used technique for obtaining such crystals. Despite often intense efforts, many GPCRs remain refractory to crystallization, even if receptors can be purified in sufficient amounts. To address this issue, we have developed a highly efficient screening and stabilization strategy for GPCRs, based on a fluorescence thermal stability assay readout, which seems to correlate particularly well with those GPCR constructs that remain native during incorporation into the LCP. Detailed protocols are provided for rapid and cost-efficient mutant and construct generation using sequence- and ligation-independent cloning, high-throughput magnetic bead-based protein purification from small-scale expressions in mammalian cells, the screening and optimal combination of mutations for increased receptor thermostability and the rapid identification of suitable chimeric fusion protein constructs for successful crystallization in LCP. We exemplify the method on three receptors from two different classes: the neurokinin 1 receptor, the oxytocin receptor and the parathyroid hormone 1 receptor.

Untangling the complexity of membrane protein folding

Delineating the folding steps of helical-bundle membrane proteins has been a challenging task. Many questions remain unanswered, including the conformation and stability of the states populated during folding, the shape of the energy barriers between the states, and the role of lipids as a solvent in mediating the folding. Recently, theoretical frames have matured to a point that permits detailed dissection of the folding steps, and advances in experimental techniques at both single-molecule and ensemble levels enable selective modulation of specific steps for quantitative determination of the folding energy landscapes. We also discuss how lipid molecules would play an active role in shaping the folding energy landscape of membrane proteins, and how folding of multi-domain membrane proteins can be understood based on our current knowledge. We conclude this review by offering an outlook for emerging questions in the study of membrane protein folding.

Free-energy changes of bacteriorhodopsin point mutants measured by single-molecule force spectroscopy

Single amino acid mutations provide quantitative insight into the energetics that underlie the dynamics and folding of membrane proteins. Chemical denaturation is the most widely used assay and yields the change in unfolding free energy (ΔΔG). It has been applied to >80 different residues of bacteriorhodopsin (bR), a model membrane protein. However, such experiments have several key limitations: 1) a nonnative lipid environment, 2) a denatured state with significant secondary structure, 3) error introduced by extrapolation to zero denaturant, and 4) the requirement of globally reversible refolding. We overcame these limitations by reversibly unfolding local regions of an individual protein with mechanical force using an atomic-force-microscope assay optimized for 2 μs time resolution and 1 pN force stability. In this assay, bR was unfolded from its native bilayer into a well-defined, stretched state. To measure ΔΔG, we introduced two alanine point mutations into an 8-amino-acid region at the C-terminal end of bR's G helix. For each, we reversibly unfolded and refolded this region hundreds of times while the rest of the protein remained folded. Our single-molecule-derived ΔΔG for mutant L223A (-2.3 ± 0.6 kcal/mol) quantitatively agreed with past chemical denaturation results while our ΔΔG for mutant V217A was 2.2-fold larger (-2.4 ± 0.6 kcal/mol). We attribute the latter result, in part, to contact between Val²¹⁷ and a natively bound squalene lipid, highlighting the contribution of membrane protein-lipid contacts not present in chemical denaturation assays. More generally, we established a platform for determining ΔΔG for a fully folded membrane protein embedded in its native bilayer.

Functional Oligomeric Forms of Uncoupling Protein 2: Strong Evidence for Asymmetry in Protein and Lipid Bilayer Systems

Stoichiometry of uncoupling proteins (UCPs) and their coexistence as functional monomeric and associated forms in lipid membranes remain intriguing open questions. In this study, tertiary and quaternary structures of UCP2 were analyzed experimentally and through molecular dynamics (MD) simulations. UCP2 was overexpressed in the inner membrane of Escherichia coli, then purified and reconstituted in lipid vesicles. Structure and proton transport function of UCP2 were characterized by circular dichroism (CD) spectroscopy and fluorescence methods. Findings suggest a tetrameric functional form for UCP2. MD simulations conclude that tetrameric UCP2 is a dimer of dimers, is more stable than its monomeric and dimeric forms, is asymmetrical and induces asymmetry in the membrane's lipid structure, and a biphasic on-off switch between the dimeric units is its possible mode of transport. MD simulations also show that the water density inside the UCP2 monomer is asymmetric, with the cytoplasmic side having a higher water density and a wider radius. In contrast, the structurally comparable adenosine 5'-diphosphate (ADP)/adenosine 5'-triphosphate (ATP) carrier (AAC1) did not form tetramers, implying that tetramerization cannot be generalized to all mitochondrial carriers.

Optimizing the linker length for fusing an alcohol dehydrogenase with a cyclohexanone monooxygenase

The use of enzymes in organic synthesis is highly appealing due their remarkably high chemo-, regio- and enantioselectivity. Nevertheless, for biosynthetic routes to be industrially useful, the enzymes must fulfill several requirements. Particularly, in case of cofactor-dependent enzymes self-sufficient systems are highly valuable. This can be achieved by fusing enzymes with complementary cofactor dependency. Such bifunctional enzymes are also relatively easy to handle, may enhance stability, and promote product intermediate channeling. However, usually the characteristics of the linker, fusing the target enzymes, are not thoroughly evaluated. A poor linker design can lead to detrimental effects on expression levels, enzyme stability and/or enzyme performance. In this chapter, the effect of the length of a glycine-rich linker was explored for the case study of ɛ-caprolactone synthesis through an alcohol dehydrogenase-cyclohexanone monooxygenase fusion system. The procedure includes cloning of linker variants, expression analysis, determination of thermostability and effect on activity and conversion levels of 15 variants of different linker sizes. The protocols can also be used for the creation of other protein-protein fusions.

Structural cavities are critical to balancing stability and activity of a membrane-integral enzyme

Packing interaction is a critical driving force in the folding of helical membrane proteins. Despite the importance, packing defects (i.e., cavities including voids, pockets, and pores) are prevalent in membrane-integral enzymes, channels, transporters, and receptors, playing essential roles in function. Then, a question arises regarding how the two competing requirements, packing for stability vs. cavities for function, are reconciled in membrane protein structures. Here, using the intramembrane protease GlpG of Escherichia coli as a model and cavity-filling mutation as a probe, we tested the impacts of native cavities on the thermodynamic stability and function of a membrane protein. We find several stabilizing mutations which induce substantial activity reduction without distorting the active site. Notably, these mutations are all mapped onto the regions of conformational flexibility and functional importance, indicating that the cavities facilitate functional movement of GlpG while compromising the stability. Experiment and molecular dynamics simulation suggest that the stabilization is induced by the coupling between enhanced protein packing and weakly unfavorable lipid desolvation, or solely by favorable lipid solvation on the cavities. Our result suggests that, stabilized by the relatively weak interactions with lipids, cavities are accommodated in membrane proteins without severe energetic cost, which, in turn, serve as a platform to fine-tune the balance between stability and flexibility for optimal activity.

Quantifying the Native Energetics Stabilizing Bacteriorhodopsin by Single-Molecule Force Spectroscopy

We quantified the equilibrium (un)folding free energy ΔG_{0} of an eight-amino-acid region starting from the fully folded state of the model membrane-protein bacteriorhodopsin using single-molecule force spectroscopy. Analysis of equilibrium and nonequilibrium data yielded consistent, high-precision determinations of ΔG_{0} via multiple techniques (force-dependent kinetics, Crooks fluctuation theorem, and inverse Boltzmann analysis). We also deduced the full 1D projection of the free-energy landscape in this region. Importantly, ΔG_{0} was determined in bacteriorhodopsin's native bilayer, an advance over traditional results obtained by chemical denaturation in nonphysiological detergent micelles.

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.

Protein structure predictions by enhanced conformational sampling methods

In this Special Festschrift Issue for the celebration of Professor Nobuhiro Gō's 80th birthday, we review enhanced conformational sampling methods for protein structure predictions. We present several generalized-ensemble algorithms such as multicanonical algorithm, replica-exchange method, etc. and parallel Monte Carlo or molecular dynamics method with genetic crossover. Examples of the results of these methods applied to the predictions of protein tertiary structures are also presented.

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.

A concise method for quantitative analysis of interactions between lipids and membrane proteins

Although interactions between lipids and membrane proteins (MPs) have been considered crucially important for understanding the functions of lipids, lack of useful and convincing experimental methods has hampered the analysis of the interactions. Here, we developed a surface plasmon resonance (SPR)-based concise method for quantitative analysis of lipid-MP interactions, coating the sensor chip surface with self-assembled monolayer (SAM) with C₆-chain. To develop this method, we used bacteriorhodopsin (bR) as an MP, and examined its interaction with various types of lipids. The merits of using C₆-SAM-modified sensor chip are as follows: (1) alkyl-chains of SAM confer a better immobilization of MPs because of the efficient preconcentration due to hydrophobic contacts; (2) SAM provides immobilized MPs with a partial membranous environment, which is important for the stabilization of MPs; and (3) a thinner C₆-SAM layer (1 nm) compared with MP size forces the MP to bulge outward from the SAM surface, allowing extraneously injected lipids to be accessible to the hydrophobic transmembrane regions. Actually, the amount of bR immobilized on C₆-SAM is 10 times higher than that on a hydrophilic CM5 sensor chip, and AFM observations confirmed that bR molecules are exposed on the SAM surface. Of the lipids tested, S-TGA-1, a halobacterium-derived glycolipid, had the highest specificity to bR with a nanomolar dissociation constant. This is consistent with the reported co-crystal structure that indicates the formation of several intermolecular hydrogen bonds. Therefore, we not only reproduced the specific lipid-bR recognition, but also succeeded in its quantitative evaluation, demonstrating the validity and utility of this method.

Consensus designs and thermal stability determinants of a human glutamate transporter

Human excitatory amino acid transporters (EAATs) take up the neurotransmitter glutamate in the brain and are essential to maintain excitatory neurotransmission. Our understanding of the EAATs’ molecular mechanisms has been hampered by the lack of stability of purified protein samples for biophysical analyses. Here, we present approaches based on consensus mutagenesis to obtain thermostable EAAT1 variants that share up to ~95% amino acid identity with the wild type transporters, and remain natively folded and functional. Structural analyses of EAAT1 and the consensus designs using hydrogen-deuterium exchange linked to mass spectrometry show that small and highly cooperative unfolding events at the inter-subunit interface rate-limit their thermal denaturation, while the transport domain unfolds at a later stage in the unfolding pathway. Our findings provide structural insights into the kinetic stability of human glutamate transporters, and introduce general approaches to extend the lifetime of human membrane proteins for biophysical analyses.

Characterization of Denatured States and Reversible Unfolding of Sensory Rhodopsin II

Our understanding on the folding of membrane proteins lags behind that of soluble proteins due to challenges posed by the exposure of hydrophobic regions during in vitro chemical denaturation and refolding experiments. While different folding models are accepted for soluble proteins, only the two-stage model and the long-range interactions model have been proposed so far for helical membrane proteins. To address our knowledge gap on how different membrane proteins traverse their folding pathways, we have systematically investigated the structural features of SDS-denatured states and the kinetics for reversible unfolding of sensory rhodopsin II (pSRII), a retinal-binding photophobic receptor from Natronomonas pharaonis. pSRII is difficult to denature, and only SDS can dislodge the retinal chromophore without rapid aggregation. Even in 30% SDS (0.998 Χ_SDS), pSRII retains the equivalent of six out of seven transmembrane helices, while the retinal-binding pocket is disrupted, with transmembrane residues becoming more solvent exposed. Folding of pSRII from an SDS-denatured state harboring a covalently bound retinal chromophore shows deviations from an apparent two-state behavior. SDS denaturation to form the sensory opsin apo-protein is reversible. We report pSRII as a new model protein which is suitable for membrane protein folding studies and has a unique folding mechanism that differs from those of bacteriorhodopsin and bovine rhodopsin.

Bicelle size modulates the rate of bacteriorhodopsin folding

The conformational equilibria of integral membrane proteins have proven extremely difficult to characterize within native lipid bilayers. To circumvent technical issues, investigations of the structure and stability of α-helical membrane proteins are often carried out in mixed micelle or bicelle solvents that mimic the membrane and facilitate measurements of reversible folding. Under these conditions, the energetics of membrane protein folding are typically proportional to the mole fraction of an anionic detergent in the micelle. However, investigations of the folding and unfolding of bacteriorhodopsin (bR) surprisingly revealed that the folding rate is also highly sensitive to the bulk molar concentration of lipids and detergents. We show here that this rate enhancement coincides with changes in bicelle size and suggest this effect arises through restriction of the conformational search space during folding. In conjunction with previous mutagenic studies, these results provide additional evidence that a topological search limits the rate of bR folding. Furthermore, this finding provides insights into the manner by which micellar and bicellar environments influence the conformational stability of polytopic membrane proteins.

Optimization and 13CH3 methionine labeling of a signaling competent neurotensin receptor 1 variant for NMR studies

Neurotensin is a 13-residue peptide that acts as a neuromodulator of classical neurotransmitters such as dopamine and glutamate in the mammalian central nervous system, mainly by activating the G protein-coupled receptor (GPCR), neurotensin receptor 1 (NTS₁). Agonist binding to GPCRs shifts the conformational equilibrium of the transmembrane helices towards distinct, thermodynamically favorable conformations that favor effector protein interactions and promotes cell signaling. The introduction of site specific labels for NMR spectroscopy has proven useful for investigating this dynamic process, but the low expression levels and poor stability of GPCRs is a hindrance to solution NMR experiments. Several thermostabilized mutants of NTS₁ have been engineered to circumvent this, with the crystal structures of four of these published. The conformational dynamics of NTS₁ however, has not been thoroughly investigated with NMR. It is generally accepted that stabilized GPCRs exhibit attenuated signaling, thus we thoroughly characterized the signaling characteristics of several thermostabilized NTS₁ variants to identify an optimal variant for protein NMR studies. A variant termed enNTS₁ exhibited the best combination of signaling capability and stability upon solubilization with detergents. enNTS₁ was subsequently labeled with ¹³CH₃-methionine in E. coli and purified to homogeneity in the absence of bound ligands. Using solution NMR spectroscopy we observed several well dispersed ¹³CH₃-methionine resonances, many of which exhibited chemical shift changes upon the addition of the high affinity agonist peptide, NT8-13. Thus, enNTS₁ represents a novel tool for investigating ligand induced conformational changes in NTS₁ to gain insights into the molecular mechanisms underlying neurotensin signaling.

Kinetic stability of membrane proteins

Although membrane proteins constitute an important class of biomolecules involved in key cellular processes, study of the thermodynamic and kinetic stability of their structures is far behind that of soluble proteins. It is known that many membrane proteins become unstable when removed by detergent extraction from the lipid environment. In addition, most of them undergo irreversible denaturation, even under mild experimental conditions. This process was found to be associated with partial unfolding of the polypeptide chain exposing hydrophobic regions to water, and it was proposed that the formation of kinetically trapped conformations could be involved. In this review, we will describe some of the efforts toward understanding the irreversible inactivation of membrane proteins. Furthermore, its modulation by phospholipids, ligands, and temperature will be herein discussed.

Role of cardiolipin in stability of integral membrane proteins

Cardiolipin (CL) is a unique phospholipid with a dimeric structure having four acyl chains and two phosphate groups found almost exclusively in certain membranes of bacteria and of mitochondria of eukaryotes. CL interacts with numerous proteins and has been implicated in function and stabilization of several integral membrane proteins (IMPs). While both functional and stabilization roles of CL in IMPs has been generally acknowledged, there are, in fact, only limited number of quantitative analysis that support this function of CL. This is likely caused by relatively complex determination of parameters characterizing stability of IMPs and particularly intricate assessment of role of specific phospholipids such as CL in IMPs stability. This review aims to summarize quantitative findings regarding stabilization role of CL in IMPs reported up to now.

Discrimination of Native-like States of Membrane Proteins with Implicit Membrane-based Scoring Functions

A scoring protocol based on implicit membrane-based scoring functions and a new protocol for optimizing the positioning of proteins inside the membrane was evaluated for its capacity to discriminate native-like states from misfolded decoys. A decoy set previously established by the Baker lab (Proteins: Struct., Funct., Genet. 2006, 62, 1010-1025) was used along with a second set that was generated to cover higher resolution models. The Implicit Membrane Model 1 (IMM1), IMM1 model with CHARMM 36 parameters (IMM1-p36), generalized Born with simple switching (GBSW), and heterogeneous dielectric generalized Born versions 2 (HDGBv2) and 3 (HDGBv3) were tested along with the new HDGB van der Waals (HDGBvdW) model that adds implicit van der Waals contributions to the solvation free energy. For comparison, scores were also calculated with the distance-scaled finite ideal-gas reference (DFIRE) scoring function. Z-scores for native state discrimination, energy vs root-mean-square deviation (RMSD) correlations, and the ability to select the most native-like structures as top-scoring decoys were evaluated to assess the performance of the scoring functions. Ranking of the decoys in the Baker set that were relatively far from the native state was challenging and dominated largely by packing interactions that were captured best by DFIRE with less benefit of the implicit membrane-based models. Accounting for the membrane environment was much more important in the second decoy set where especially the HDGB-based scoring functions performed very well in ranking decoys and providing significant correlations between scores and RMSD, which shows promise for improving membrane protein structure prediction and refinement applications. The new membrane structure scoring protocol was implemented in the MEMScore web server ( http://feiglab.org/memscore ).

Reversible Unfolding of Rhomboid Intramembrane Proteases

Denaturant-induced unfolding of helical membrane proteins provides insights into their mechanism of folding and domain organization, which take place in the chemically heterogeneous, anisotropic environment of a lipid membrane. Rhomboid proteases are intramembrane proteases that play key roles in various diseases. Crystal structures have revealed a compact helical bundle with a buried active site, which requires conformational changes for the cleavage of transmembrane substrates. A dimeric form of the rhomboid protease has been shown to be important for activity. In this study, we examine the mechanism of refolding for two distinct rhomboids to gain insight into their secondary structure-activity relationships. Although helicity is largely abolished in the unfolded states of both proteins, unfolding is completely reversible for HiGlpG but only partially reversible for PsAarA. Refolding of both proteins results in reassociation of the dimer, with a 90% regain of catalytic activity for HiGlpG but only a 70% regain for PsAarA. For both proteins, a broad, gradual transition from the native, folded state to the denatured, partly unfolded state was revealed with the aid of circular dichroism spectroscopy as a function of denaturant concentration, thus arguing against a classical two-state model as found for many globular soluble proteins. Thermal denaturation has irreversible destabilizing effects on both proteins, yet reveals important functional details regarding substrate accessibility to the buried active site. This concerted biophysical and functional analysis demonstrates that HiGlpG, with a simple six-transmembrane-segment organization, is more robust than PsAarA, which has seven predicted transmembrane segments, thus rendering HiGlpG amenable to in vitro studies of membrane-protein folding.

Convergence and Sampling in Determining Free Energy Landscapes for Membrane Protein Association

Potential of mean force (PMF) calculations are used to characterize the free energy landscape of protein–lipid and protein–protein association within membranes. Coarse-grained simulations allow binding free energies to be determined with reasonable statistical error. This accuracy relies on defining a good collective variable to describe the binding and unbinding transitions, and upon criteria for assessing the convergence of the simulation toward representative equilibrium sampling. As examples, we calculate protein–lipid binding PMFs for ANT/cardiolipin and Kir2.2/PIP₂, using umbrella sampling on a distance coordinate. These highlight the importance of replica exchange between windows for convergence. The use of two independent sets of simulations, initiated from bound and unbound states, provide strong evidence for simulation convergence. For a model protein–protein interaction within a membrane, center-of-mass distance is shown to be a poor collective variable for describing transmembrane helix–helix dimerization. Instead, we employ an alternative intermolecular distance matrix RMS (D_RMS) coordinate to obtain converged PMFs for the association of the glycophorin transmembrane domain. While the coarse-grained force field gives a reasonable K_d for dimerization, the majority of the bound population is revealed to be in a near-native conformation. Thus, the combination of a refined reaction coordinate with improved sampling reveals previously unnoticed complexities of the dimerization free energy landscape. We propose the use of replica-exchange umbrella sampling starting from different initial conditions as a robust approach for calculation of the binding energies in membrane simulations.

Observation of helix associations for insertion of a retinal molecule and distortions of helix structures in bacteriorhodopsin

We applied a newly proposed prediction method for membrane protein structures to bacteriorhodopsin that has distorted transmembrane helices in the native structure. This method uses an implicit membrane model, which restricts sampling space during folding in a membrane region, and includes helix bending. Replica-exchange simulations were performed with seven transmembrane helices only without a retinal molecule. Obtained structures were classified into clusters of similar structures, which correspond to local-minimum free energy states. The two lowest free energy states corresponded to a native-like structure with the correct empty space for retinal and a structure with this empty space filled with a helix. Previous experiments of bacteriorhodopsin suggested that association of transmembrane helices enables them to make a room for insertion of a retinal. Our results are consistent with these results. Moreover, distortions of helices in the native-like structures were successfully reproduced. In the distortions, whereas the locations of kinks for all helices were similar to those of Protein Data Bank's data, the amount of bends was more similar for helices away from the retinal than for those close to the retinal in the native structure. This suggests a hypothesis that the amino-acid sequence specifies the location of kinks in transmembrane helices and that the amount of distortions depends on the interactions with the surrounding molecules such as neighboring helices, lipids, and retinal.

Documentation of an Imperative To Improve Methods for Predicting Membrane Protein Stability

There is a compelling and growing need to accurately predict the impact of amino acid mutations on protein stability for problems in personalized medicine and other applications. Here the ability of 10 computational tools to accurately predict mutation-induced perturbation of folding stability (ΔΔG) for membrane proteins of known structure was assessed. All methods for predicting ΔΔG values performed significantly worse when applied to membrane proteins than when applied to soluble proteins, yielding estimated concordance, Pearson, and Spearman correlation coefficients of <0.4 for membrane proteins. Rosetta and PROVEAN showed a modest ability to classify mutations as destabilizing (ΔΔG < −0.5 kcal/mol), with a 7 in 10 chance of correctly discriminating a randomly chosen destabilizing variant from a randomly chosen stabilizing variant. However, even this performance is significantly worse than for soluble proteins. This study highlights the need for further development of reliable and reproducible methods for predicting thermodynamic folding stability in membrane proteins.

A many-body term improves the accuracy of effective potentials based on protein coevolutionary data

The study of correlated mutations in alignments of homologous proteins proved to be successful not only in the prediction of their native conformation but also in the development of a two-body effective potential between pairs of amino acids. In the present work, we extend the effective potential, introducing a many-body term based on the same theoretical framework, making use of a principle of maximum entropy. The extended potential performs better than the two-body one in predicting the energetic effect of 308 mutations in 14 proteins (including membrane proteins). The average value of the parameters of the many-body term correlates with the degree of hydrophobicity of the corresponding residues, suggesting that this term partly reflects the effect of the solvent.

Principal component analysis of binding energies for single-point mutants of hT2R16 bound to an agonist correlate with experimental mutant cell response

Directed evolution is a technique that enables the identification of mutants of a particular protein that carry a desired property by successive rounds of random mutagenesis, screening, and selection. This technique has many applications, including the development of G protein-coupled receptor-based biosensors and designer drugs for personalized medicine. Although effective, directed evolution is not without challenges and can greatly benefit from the development of computational techniques to predict the functional outcome of single-point amino acid substitutions. In this article, we describe a molecular dynamics-based approach to predict the effects of single amino acid substitutions on agonist binding (salicin) to a human bitter taste receptor (hT2R16). An experimentally determined functional map of single-point amino acid substitutions was used to validate the whole-protein molecular dynamics-based predictive functions. Molecular docking was used to construct a wild-type agonist-receptor complex, providing a starting structure for single-point substitution simulations. The effects of each single amino acid substitution in the functional response of the receptor to its agonist were estimated using three binding energy schemes with increasing inclusion of solvation effects. We show that molecular docking combined with molecular mechanics simulations of single-point mutants of the agonist-receptor complex accurately predicts the functional outcome of single amino acid substitutions in a human bitter taste receptor.

A delocalized proton-binding site within a membrane protein

The role of protein-bound water molecules in protein function and catalysis is an emerging topic. Here, we studied the solvation of an excess proton by protein-bound water molecules and the contribution of the surrounding amino acid residues at the proton release site of the membrane protein bacteriorhodopsin. It hosts an excess proton within a protein-bound water cluster, which is hydrogen bonded to several surrounding amino acids. Indicative of delocalization is a broad continuum absorbance experimentally observed by time-resolved Fourier transform infrared spectroscopy. In combination with site-directed mutagenesis, the involvement of several amino acids (especially Glu-194 and Glu-204) in the delocalization was elaborated. Details regarding the contributions of the glutamates and water molecules to the delocalization mode in biomolecular simulations are controversial. We carried out quantum mechanics/molecular mechanics (QM/MM) self-consistent charge density functional tight-binding simulations for all amino acids that have been experimentally shown to be involved in solvation of the excess proton, and systematically investigated the influence of the quantum box size. We compared calculated theoretical infrared spectra with experimental ones as a measure for the correct description of excess proton delocalization. A continuum absorbance can only be observed for small quantum boxes containing few amino acids and/or water molecules. Larger quantum boxes, including all experimentally shown involved amino acids, resulted in narrow absorbance bands, indicating protonation of a single binding site in contradiction to experimental results. We conclude that small quantum boxes seem to reproduce representative extreme cases of proton delocalization modes: proton delocalization only on water molecules or only between Glu-194 and Glu-204. Extending the experimental spectral region to lower wave numbers, a water-delocalized proton reproduces the observed continuum absorbance better than a glutamate-shared delocalized proton. However, a full agreement between QM simulations and experimental results on the delocalized excess proton will require a larger quantum box as well as more sophisticated QM/MM methods.

Transmembrane helix straightening and buckling underlies activation of mechanosensitive and thermosensitive K(2P) channels

Mechanical and thermal activation of ion channels is central to touch, thermosensation, and pain. The TRAAK/TREK K(2P) potassium channel subfamily produces background currents that alter neuronal excitability in response to pressure, temperature, signaling lipids, and anesthetics. How such diverse stimuli control channel function is unclear. Here we report structures of K(2P)4.1 (TRAAK) bearing C-type gate-activating mutations that reveal a tilting and straightening of the M4 inner transmembrane helix and a buckling of the M2 transmembrane helix. These conformational changes move M4 in a direction opposite to that in classical potassium channel activation mechanisms and open a passage lateral to the pore that faces the lipid bilayer inner leaflet. Together, our findings uncover a unique aspect of K(2P) modulation, indicate a means for how the K(2P) C-terminal cytoplasmic domain affects the C-type gate which lies ∼40Å away, and suggest how lipids and bilayer inner leaflet deformations may gate the channel.

Crystal structure of Escherichia coli-expressed Haloarcula marismortui bacteriorhodopsin I in the trimeric form

Bacteriorhodopsins are a large family of seven-helical transmembrane proteins that function as light-driven proton pumps. Here, we present the crystal structure of a new member of the family, Haloarcula marismortui bacteriorhodopsin I (HmBRI) D94N mutant, at the resolution of 2.5 Å. While the HmBRI retinal-binding pocket and proton donor site are similar to those of other archaeal proton pumps, its proton release region is extended and contains additional water molecules. The protein's fold is reinforced by three novel inter-helical hydrogen bonds, two of which result from double substitutions relative to Halobacterium salinarum bacteriorhodopsin and other similar proteins. Despite the expression in Escherichia coli and consequent absence of native lipids, the protein assembles as a trimer in crystals. The unique extended loop between the helices D and E of HmBRI makes contacts with the adjacent protomer and appears to stabilize the interface. Many lipidic hydrophobic tail groups are discernible in the membrane region, and their positions are similar to those of archaeal isoprenoid lipids in the crystals of other proton pumps, isolated from native or native-like sources. All these features might explain the HmBRI properties and establish the protein as a novel model for the microbial rhodopsin proton pumping studies.

Bacteriorhodopsin folds through a poorly organized transition state

The folding mechanisms of helical membrane proteins remain largely uncharted. Here we characterize the kinetics of bacteriorhodopsin folding and employ φ-value analysis to explore the folding transition state. First, we developed and confirmed a kinetic model that allowed us to assess the rate of folding from SDS-denatured bacteriorhodopsin (bRU) and provides accurate thermodynamic information even under influence of retinal hydrolysis. Next, we obtained reliable φ-values for 16 mutants of bacteriorhodopsin with good coverage across the protein. Every φ-value was less than 0.4, indicating the transition state is not uniquely structured. We suggest that the transition state is a loosely organized ensemble of conformations.

The safety dance: biophysics of membrane protein folding and misfolding in a cellular context

Most biological processes require the production and degradation of proteins, a task that weighs heavily on the cell. Mutations that compromise the conformational stability of proteins place both specific and general burdens on cellular protein homeostasis (proteostasis) in ways that contribute to numerous diseases. Efforts to elucidate the chain of molecular events responsible for diseases of protein folding address one of the foremost challenges in biomedical science. However, relatively little is known about the processes by which mutations prompt the misfolding of α-helical membrane proteins, which rely on an intricate network of cellular machinery to acquire and maintain their functional structures within cellular membranes. In this review, we summarize the current understanding of the physical principles that guide membrane protein biogenesis and folding in the context of mammalian cells. Additionally, we explore how pathogenic mutations that influence biogenesis may differ from those that disrupt folding and assembly, as well as how this may relate to disease mechanisms and therapeutic intervention. These perspectives indicate an imperative for the use of information from structural, cellular, and biochemical studies of membrane proteins in the design of novel therapeutics and in personalized medicine.

Sequence-specific dimerization of a transmembrane helix in amphipol A8-35

As traditional detergents might destabilize or even denature membrane proteins, amphiphilic polymers have moved into the focus of membrane-protein research in recent years. Thus far, Amphipols are the best studied amphiphilic copolymers, having a hydrophilic backbone with short hydrophobic chains. However, since stabilizing as well as destabilizing effects of the Amphipol belt on the structure of membrane proteins have been described, we systematically analyze the impact of the most commonly used Amphipol A8-35 on the structure and stability of a well-defined transmembrane protein model, the glycophorin A transmembrane helix dimer. Amphipols are not able to directly extract proteins from their native membranes, and detergents are typically replaced by Amphipols only after protein extraction from membranes. As Amphipols form mixed micelles with detergents, a better understanding of Amphipol-detergent interactions is required. Therefore, we analyze the interaction of A8-35 with the anionic detergent sodium dodecyl sulfate and describe the impact of the mixed-micelle-like system on the stability of a transmembrane helix dimer. As A8-35 may highly stabilize and thereby rigidify a transmembrane protein structure, modest destabilization by controlled addition of detergents and formation of mixed micellar systems might be helpful to preserve the function of a membrane protein in Amphipol environments.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Bacteriorhodopsin: Would the real structural intermediates please stand up?

BACKGROUND

Bacteriorhodopsin (bR) is the simplest known light driven proton pump and has been heavily studied using structural methods: eighty four X-ray diffraction, six electron diffraction and three NMR structures of bR are deposited within the protein data bank. Twenty one X-ray structures report light induced structural changes and changes induced by mutation, changes in pH, thermal annealing or X-ray induced photo-reduction have also been examined.

SCOPE OF REVIEW

We argue that light-induced structural changes that are replicated across several studies by independent research groups are those most likely to represent what is happening in reality. We present both internal distance matrix analyses that sort deposited bR structures into hierarchal trees, and difference Fourier analysis of deposited X-ray diffraction data.

MAJOR CONCLUSIONS

An internal distance matrix analysis separates most wild-type bR structures according to their different crystal forms, indicating how the protein's structure is influenced by crystallization conditions. A similar analysis clusters eleven studies of illuminated bR crystals as one branch of a hierarchal tree with reproducible movements of the extracellular portion of helix C towards helix G, and of the cytoplasmic portion of helix F away from helices A, B and G. All crystallographic data deposited for illuminated crystals show negative difference density on a water molecule (Wat402) that forms H-bonds to the retinal Schiff Base and two aspartate residues (Asp85, Asp212) in the bR resting state. Other recurring difference density features indicated reproducible side-chain, backbone and water molecule displacements. X-ray induced radiation damage also disorders Wat402 but acts via cleaving the head-groups of Asp85 and Asp212.

GENERAL SIGNIFICANCE

A remarkable level of agreement exists when deposited structures and crystallographic observations are viewed as a whole. From this agreement a unified picture of the structural mechanism of light-induced proton pumping by bR emerges. This article is part of a Special Issue entitled Structural biochemistry and biophysics of membrane proteins.

Intrinsic conformational plasticity of native EmrE provides a pathway for multidrug resistance

EmrE is a multidrug resistance efflux pump with specificity to a wide range of antibiotics and antiseptics. To obtain atomic-scale insight into the attributes of the native state that encodes the broad specificity, we used a hybrid of solution and solid-state NMR methods in lipid bilayers and bicelles. Our results indicate that the native EmrE dimer oscillates between inward and outward facing structural conformations at an exchange rate (k(ex)) of ~300 s(-1) at 37 °C (millisecond motions), which is ~50-fold faster relative to the tetraphenylphosphonium (TPP(+)) substrate-bound form of the protein. These observables provide quantitative evidence that the rate-limiting step in the TPP(+) transport cycle is not the outward-inward conformational change in the absence of drug. In addition, using differential scanning calorimetry, we found that the width of the gel-to-liquid crystalline phase transition was 2 °C broader in the absence of the TPP(+) substrate versus its presence, which suggested that changes in transporter dynamics can impact the phase properties of the membrane. Interestingly, experiments with cross-linked EmrE showed that the millisecond inward-open to outward-open dynamics was not the culprit of the broadening. Instead, the calorimetry and NMR data supported the conclusion that faster time scale structural dynamics (nanosecond-microsecond) were the source and therefore impart the conformationally plastic character of native EmrE capable of binding structurally diverse substrates. These findings provide a clear example how differences in membrane protein transporter structural dynamics between drug-free and bound states can have a direct impact on the physical properties of the lipid bilayer in an allosteric fashion.

Unfolding study of a trimeric membrane protein AcrB

The folding of a multi-domain trimeric α-helical membrane protein, Escherichia coli inner membrane protein AcrB, was investigated. AcrB contains both a transmembrane domain and a large periplasmic domain. Protein unfolding in sodium dodecyl sulfate (SDS) and urea was monitored using the intrinsic fluorescence and circular dichroism spectroscopy. The SDS denaturation curve displayed a sigmoidal profile, which could be fitted with a two-state unfolding model. To investigate the unfolding of separate domains, a triple mutant was created, in which all three Trp residues in the transmembrane domain were replaced with Phe. The SDS unfolding profile of the mutant was comparable to that of the wild type AcrB, suggesting that the observed signal change was largely originated from the unfolding of the soluble domain. Strengthening of trimer association through the introduction of an inter-subunit disulfide bond had little effect on the unfolding profile, suggesting that trimer dissociation was not the rate-limiting step in unfolding monitored by fluorescence emission. Under our experimental condition, AcrB unfolding was not reversible. Furthermore, we experimented with the refolding of a monomeric mutant, AcrBΔloop , from the SDS unfolded state. The CD spectrum of the refolded AcrBΔloop superimposed well onto the spectra of the original folded protein, while the fluorescence spectrum was not fully recovered. In summary, our results suggested that the unfolding of the trimeric AcrB started with a local structural rearrangement. While the refolding of secondary structure in individual monomers could be achieved, the re-association of the trimer might be the limiting factor to obtain folded wild-type AcrB.

Kinetics and thermodynamics of membrane protein folding

Understanding protein folding has been one of the great challenges in biochemistry and molecular biophysics. Over the past 50 years, many thermodynamic and kinetic studies have been performed addressing the stability of globular proteins. In comparison, advances in the membrane protein folding field lag far behind. Although membrane proteins constitute about a third of the proteins encoded in known genomes, stability studies on membrane proteins have been impaired due to experimental limitations. Furthermore, no systematic experimental strategies are available for folding these biomolecules in vitro. Common denaturing agents such as chaotropes usually do not work on helical membrane proteins, and ionic detergents have been successful denaturants only in few cases. Refolding a membrane protein seems to be a craftsman work, which is relatively straightforward for transmembrane β-barrel proteins but challenging for α-helical membrane proteins. Additional complexities emerge in multidomain membrane proteins, data interpretation being one of the most critical. In this review, we will describe some recent efforts in understanding the folding mechanism of membrane proteins that have been reversibly refolded allowing both thermodynamic and kinetic analysis. This information will be discussed in the context of current paradigms in the protein folding field.

Small-angle scattering gives direct structural information about a membrane protein inside a lipid environment

Monomeric bacteriorhodopsin (bR) reconstituted into POPC/POPG-containing nanodiscs was investigated by combined small-angle neutron and X-ray scattering. A novel hybrid approach to small-angle scattering data analysis was developed. In combination, these provided direct structural insight into membrane-protein localization in the nanodisc and into the protein-lipid interactions. It was found that bR is laterally decentred in the plane of the disc and is slightly tilted in the phospholipid bilayer. The thickness of the bilayer is reduced in response to the incorporation of bR. The observed tilt of bR is in good accordance with previously performed theoretical predictions and computer simulations based on the bR crystal structure. The result is a significant and essential step on the way to developing a general small-angle scattering-based method for determining the low-resolution structures of membrane proteins in physiologically relevant environments.

The free energy landscape of dimerization of a membrane protein, NanC

Membrane proteins are frequently present in crowded environments, which favour lateral association and, on occasions, two-dimensional crystallization. To better understand the non-specific lateral association of a membrane protein we have characterized the free energy landscape for the dimerization of a bacterial outer membrane protein, NanC, in a phospholipid bilayer membrane. NanC is a member of the KdgM-family of bacterial outer membrane proteins and is responsible for sialic acid transport in E. coli. Umbrella sampling and coarse-grained molecular dynamics were employed to calculate the potentials of mean force (PMF) for a variety of restrained relative orientations of two NanC proteins as the separation of their centres of mass was varied. We found the free energy of dimerization for NanC to be in the range of to . Differences in the depths of the PMFs for the various orientations are related to the shape of the proteins. This was quantified by calculating the lipid-inaccessible buried surface area of the proteins in the region around the minimum of each PMF. The depth of the potential well of the PMF was shown to depend approximately linearly on the buried surface area. We were able to resolve local minima in the restrained PMFs that would not be revealed using conventional umbrella sampling. In particular, these features reflected the local organization of the intervening lipids between the two interacting proteins. Through a comparison with the distribution of lipids around a single freely-diffusing NanC, we were able to predict the location of these restrained local minima for the orientational configuration in which they were most pronounced. Our ability to make this prediction highlights the important role that lipid organization plays in the association of two NanCs in a bilayer.

Cells are surrounded by selectively-permeable bilayer membranes, enabling the cell to control its internal environment. Embedded within these membranes are a variety of membrane proteins, many of which facilitate this environmental control and are integral to numerous metabolic processes. Their location within the membrane and their mutual association are controlled by many factors. We use molecular dynamics simulations to investigate the free energy of association for a pair of relatively simple membrane proteins. By doing so, we are able to characterize the effect that the geometrical properties of the protein have on their mutual association in a bilayer environment, showing that there is a correlation between the buried surface area of two proteins when in contact and the strength of their interaction. We also observe the effect of protein-lipid-protein interactions in this free energy characterization. Such interactions are related to the preferential distribution of lipids around proteins in the membrane.

Measuring membrane protein stability under native conditions

The thermodynamic stability of proteins is typically measured at high denaturant concentrations and then extrapolated back to zero denaturant conditions to obtain unfolding free energies under native conditions. For membrane proteins, the extrapolations are fraught with considerable uncertainty as the denaturants may have complex effects on the membrane or micellar structure. We therefore sought to measure stability under native conditions, using a method that does not perturb the properties of the membrane or membrane mimetics. We use a technique called steric trapping to measure the thermodynamic stability of bacteriorhodopsin in bicelles and micelles. We find that bacteriorhodopsin has a high thermodynamic stability, with an unfolding free energy of ∼11 kcal/mol in dimyristoyl phosphatidylcholine bicelles. Nevertheless, the stability is much lower than predicted by extrapolation of measurements made at high denaturant concentrations. We investigated the discrepancy and found that unfolding free energy is not linear with denaturant concentration. Apparently, long extrapolations of helical membrane protein unfolding free energies must be treated with caution. Steric trapping, however, provides a method for making these measurements.

Retinal proteins as model systems for membrane protein folding

Experimental folding studies of membrane proteins are more challenging than water-soluble proteins because of the higher hydrophobicity content of membrane embedded sequences and the need to provide a hydrophobic milieu for the transmembrane regions. The first challenge is their denaturation: due to the thermodynamic instability of polar groups in the membrane, secondary structures in membrane proteins are more difficult to disrupt than in soluble proteins. The second challenge is to refold from the denatured states. Successful refolding of membrane proteins has almost always been from very subtly denatured states. Therefore, it can be useful to analyze membrane protein folding using computational methods, and we will provide results obtained with simulated unfolding of membrane protein structures using the Floppy Inclusions and Rigid Substructure Topography (FIRST) method. Computational methods have the advantage that they allow a direct comparison between diverse membrane proteins. We will review here both, experimental and FIRST studies of the retinal binding proteins bacteriorhodopsin and mammalian rhodopsin, and discuss the extension of the findings to deriving hypotheses on the mechanisms of folding of membrane proteins in general. This article is part of a Special Issue entitled: Retinal Proteins-You can teach an old dog new tricks.

Hepatitis C virus RNA replication and virus particle assembly require specific dimerization of the NS4A protein transmembrane domain

Hepatitis C virus (HCV) NS4A is a single-pass transmembrane (TM) protein essential for viral replication and particle assembly. The sequence of the NS4A TM domain is highly conserved, suggesting that it may be important for protein-protein interactions. To test this hypothesis, we measured the potential dimerization of the NS4A TM domain in a well-characterized two-hybrid TM protein interaction system. The NS4A TM domain exhibited a strong homotypic interaction that was comparable in affinity to glycophorin A, a well-studied human blood group antigen that forms TM homodimers. Several mutations predicted to cluster on a common surface of the NS4A TM helix caused significant reductions in dimerization, suggesting that these residues form an interface for NS4A dimerization. Mutations in the NS4A TM domain were further examined in the JFH-1 genotype 2a replicon system; importantly, all mutations that destabilized NS4A dimers also caused defects in RNA replication and/or virus assembly. Computational modeling of NS4A TM interactions suggests a right-handed dimeric interaction of helices with an interface that is consistent with the mutational effects. Furthermore, defects in NS4A oligomerization and virus particle assembly of two mutants were rescued by NS4A A15S, a TM mutation recently identified through forward genetics as a cell culture-adaptive mutation. Together, these data provide the first example of a functionally important TM dimer interface within an HCV nonstructural protein and reveal a fundamental role of the NS4A TM domain in coordinating HCV RNA replication and virus particle assembly.