Backbone Hydrogen Bond Strengths Can Vary Widely in Transmembrane Helices

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.

Classifying tetraspanins: A universal system for numbering residues and a proposal for naming structural motifs and subfamilies

All tetraspanins have four transmembrane domains (TMs). The large extracellular loop (LEL) that connects the third and fourth TMs contains multiple secondary structures together with the family's signature Cys-Cys-Gly motif. These intriguing membrane proteins are involved in diverse and incompletely understood cellular processes including cell adhesion, tissue differentiation, immune cell maturation and host-parasite interactions. Here we present a classification system that accurately describes the position of each amino acid within its primary sequence based on both sequence and topological conservation of the TMs and LEL. This builds on the numbering systems that have been used in the G protein-coupled receptor (GPCR) field for nearly three decades and which have aided the understanding of GPCR structure/activity relationships and ligand interactions. The high-resolution structures of the tetraspanins CD81, CD9, CD53 and Tspan15 were used to validate the structural relevance of our new tetraspanin classification system. Modelling of all tetraspanin LELs highlighted flexibility in LEL disulfide bonding across the family and suggests that the structural arrangement of tetraspanin LELs is more complex than previously thought. We therefore propose a new subfamily naming system that addresses this added complexity and facilitates the systematic classification of human tetraspanins, shedding light on all structural motifs within the family. We anticipate that our universal tetraspanin classification system will enable progress in defining how sequence and structure inform function.

Implicit model to capture electrostatic features of membrane environment

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

Implicit model to capture electrostatic features of membrane environment

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

Lipid Bilayer Strengthens the Cooperative Network of a Membrane-Integral Enzyme

Lipid bilayer provides a two-dimensional hydrophobic solvent milieu for membrane proteins in cells. Although the native bilayer is widely recognized as an optimal environment for folding and function of membrane proteins, the underlying physical basis remains elusive. Here, employing the intramembrane protease GlpG of Escherichia coli as a model, we elucidate how the bilayer stabilizes a membrane protein and engages the protein's residue interaction network compared to the nonnative hydrophobic medium, micelles. We find that the bilayer enhances GlpG stability by promoting residue burial in the protein interior compared to micelles. Strikingly, while the cooperative residue interactions cluster into multiple distinct regions in micelles, the whole packed regions of the protein act as a single cooperative unit in the bilayer. Molecular dynamics (MD) simulation indicates that lipids less efficiently solvate GlpG than detergents. Thus, the bilayerinduced enhancement of stability and cooperativity likely stems from the dominant intraprotein interactions outcompeting the weak lipid solvation. Our findings reveal a foundational mechanism in the folding, function, and quality control of membrane proteins. The enhanced cooperativity benefits function facilitating propagation of local structural perturbation across the membrane. However, the same phenomenon can render the proteins' conformational integrity vulnerable to missense mutations causing conformational diseases1,2.

A thermodynamic investigation of amyloid precursor protein processing by human γ-secretase

Human γ-secretase cleaves the transmembrane domains (TMDs) of amyloid precursor protein (APP) into pathologically relevant amyloid-β peptides (Aβs). The detailed mechanisms of the unique endoproteolytic cleavage by the presenilin 1 domain (PS1) of γ-secretase are still poorly understood. Herein, we provide thermodynamic insights into how the α-helical APP TMD is processed by γ-secretase and elucidate the specificity of Aβ48/Aβ49 cleavage using unbiased molecular dynamics and bias-exchange metadynamics simulations. The thermodynamic data show that the unwinding of APP TMD is driven by water hydration in the intracellular pocket of PS1, and the scissile bond T32-L33 or L33-V34 of the APP TMD can slide down and up to interact with D257/D385 to achieve endoproteolysis. In the wild-type system, the L33-V34 scissile bond is more easily hijacked by D257/D385 than T32-L33, resulting in higher Aβ49 cleavage, while the T32N mutation on the APP TMD decreases the energy barrier of the sliding of the scissile bonds and increases the hydrogen bond occupancy for Aβ48 cleavage. In summary, the thermodynamic analysis elucidates possible mechanisms of APP TMD processing by PS1, which might facilitate rational drug design targeting γ-secretase.

Thermodynamic analysis from unbiased molecular dynamics and bias-exchange metadynamics simulations reveals possible mechanisms on how γ-secretase cleaves the transmembrane domains of amyloid precursor protein into amyloid-β peptides.

Free-energy landscapes and insertion pathways for peptides in membrane environment

Free-energy landscapes for short peptides-specifically for variants of the pH low insertion peptide (pHLIP)-in the heterogeneous environment of a lipid bilayer or cell membrane are constructed, taking into account a set of dominant interactions and the conformational preferences of the peptide backbone. Our methodology interprets broken internal H-bonds along the backbone of a polypeptide as statistically interacting quasiparticles, activated from the helix reference state. The favored conformation depends on the local environment (ranging from polar to nonpolar), specifically on the availability of external H-bonds (with H_{2}O molecules or lipid headgroups) to replace internal H-bonds. The dominant side-chain contribution is accounted for by residue-specific transfer free energies between polar and nonpolar environments. The free-energy landscape is sensitive to the level of pH in the aqueous environment surrounding the membrane. For high pH, we identify pathways of descending free energy that suggest a coexistence of membrane-adsorbed peptides with peptides in solution. A drop in pH raises the degree of protonation of negatively charged residues and thus increases the hydrophobicity of peptide segments near the C terminus. For low pH, we identify insertion pathways between the membrane-adsorbed state and a stable trans-membrane state with the C terminus having crossed the membrane.

How physical forces drive the process of helical membrane protein folding

Protein folding is a fundamental process of life with important implications throughout biology. Indeed, tens of thousands of mutations have been associated with diseases, and most of these mutations are believed to affect protein folding rather than function. Correct folding is also a key element of design. These factors have motivated decades of research on protein folding. Unfortunately, knowledge of membrane protein folding lags that of soluble proteins. This gap is partly caused by the greater technical challenges associated with membrane protein studies, but also because of additional complexities. While soluble proteins fold in a homogenous water environment, membrane proteins fold in a setting that ranges from bulk water to highly charged to apolar. Thus, the forces that drive folding vary in different regions of the protein, and this complexity needs to be incorporated into our understanding of the folding process. Here, we review our understanding of membrane protein folding biophysics. Despite the greater challenge, better model systems and new experimental techniques are starting to unravel the forces and pathways in membrane protein folding.

Lipid bilayer induces contraction of the denatured state ensemble of a helical-bundle membrane protein

Defining the denatured state ensemble (DSE) and disordered proteins is essential to understanding folding, chaperone action, degradation, and translocation. As compared with water-soluble proteins, the DSE of membrane proteins is much less characterized. Here, we measure the DSE of the helical membrane protein GlpG of Escherichia coli (E. coli) in native-like lipid bilayers. The DSE was obtained using our steric trapping method, which couples denaturation of doubly biotinylated GlpG to binding of two streptavidin molecules. The helices and loops are probed using limited proteolysis and mass spectrometry, while the dimensions are determined using our paramagnetic biotin derivative and double electron-electron resonance spectroscopy. These data, along with our Upside simulations, identify the DSE as being highly dynamic, involving the topology changes and unfolding of some of the transmembrane (TM) helices. The DSE is expanded relative to the native state but only to 15 to 75% of the fully expanded condition. The degree of expansion depends on the local protein packing and the lipid composition. E. coli's lipid bilayer promotes the association of TM helices in the DSE and, probably in general, facilitates interhelical interactions. This tendency may be the outcome of a general lipophobic effect of proteins within the cell membranes.

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.

Diverse Scientific Benchmarks for Implicit Membrane Energy Functions

Energy functions are fundamental to biomolecular modeling. Their success depends on robust physical formalisms, efficient optimization, and high-resolution data for training and validation. Over the past 20 years, progress in each area has advanced soluble protein energy functions. Yet, energy functions for membrane proteins lag behind due to sparse and low-quality data, leading to overfit tools. To overcome this challenge, we assembled a suite of 12 tests on independent data sets varying in size, diversity, and resolution. The tests probe an energy function's ability to capture membrane protein orientation, stability, sequence, and structure. Here, we present the tests and use the franklin2019 energy function to demonstrate them. We then identify areas for energy function improvement and discuss potential future integration with machine-learning-based optimization methods. The tests are available through the Rosetta Benchmark Server (https://benchmark.graylab.jhu.edu/) and GitHub (https://github.com/rfalford12/Implicit-Membrane-Energy-Function-Benchmark).

Interhelical H-Bonds Modulate the Activity of a Polytopic Transmembrane Kinase

DesK is a Histidine Kinase that allows Bacillus subtilis to maintain lipid homeostasis in response to changes in the environment. It is located in the membrane, and has five transmembrane helices and a cytoplasmic catalytic domain. The transmembrane region triggers the phosphorylation of the catalytic domain as soon as the membrane lipids rigidify. In this research, we study how transmembrane inter-helical interactions contribute to signal transmission; we designed a co-expression system that allows studying in vivo interactions between transmembrane helices. By Alanine-replacements, we identified a group of polar uncharged residues, whose side chains contain hydrogen-bond donors or acceptors, which are required for the interaction with other DesK transmembrane helices; a particular array of H-bond- residues plays a key role in signaling, transmitting information detected at the membrane level into the cell to finally trigger an adaptive response.

Membrane proteins enter the fold

Membrane proteins have historically been recalcitrant to biophysical folding studies. However, recent adaptations of methods from the soluble protein folding field have found success in their applications to transmembrane proteins composed of both α-helical and β-barrel conformations. Avoiding aggregation is critical for the success of these experiments. Altogether these studies are leading to discoveries of folding trajectories, foundational stabilizing forces and better-defined endpoints that enable more accurate interpretation of thermodynamic data. Increased information on membrane protein folding in the cell shows that the emerging biophysical principles are largely recapitulated even in the complex biological environment.

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.

Backbone Hydrogen Bond Energies in Membrane Proteins Are Insensitive to Large Changes in Local Water Concentration

A hallmark feature of biological lipid bilayer structure is a depth-dependent polarity gradient largely resulting from the change in water concentration over the angstrom length scale. This gradient is particularly steep as it crosses the membrane interfacial regions where the water concentration drops at least a million-fold along the direction of the bilayer normal. Although local water content is often assumed to be a major determinant of membrane protein stability, the effect of the water-induced polarity gradient upon backbone hydrogen bond strength has not been systematically investigated. We addressed this question by measuring the free energy change for a number of backbone hydrogen bonds in the transmembrane protein OmpW. These values were obtained at 33 backbone amides from hydrogen/deuterium fractionation factors by nuclear magnetic resonance spectroscopy. We surprisingly found that OmpW backbone hydrogen bond energies do not vary over a wide range of water concentrations that are characteristic of the solvation environment in the bilayer interfacial region. We validated the interpretation of our results by determining the hydrodynamic and solvation properties of our OmpW-micelle complex using analytical ultracentrifugation and molecular dynamics simulations. The magnitudes of the backbone hydrogen bond free energy changes in our study are comparable to those observed in water-soluble proteins, the H-segment of the leader peptidase helix used in the von Heijne and White biological scale experiments, and several interfacial peptides. Our results agree with those reported for the transmembrane α-helical portion of the amyloid precursor protein after the latter values were adjusted for kinetic isotope effects. Overall, our work suggests that backbone hydrogen bonds provide modest thermodynamic stability to membrane protein structures and that many amides are unaffected by dehydration within the bilayer.

Semi-conductive helical homochiral metal–organic frameworks based on enantiomeric proline derivatives

Chiral ligands (S)-H₂PCA and (R)-H₂PCA were synthesized by attaching (S)- and (R)-proline to 1,4-dicarboxybenzene. They further reacted with Cu(ii) and Co(ii) ions to obtain two pairs of helical homochiral MOFs with photocatalytic behaviours.

Superfine bacterial nanocellulose produced by reverse mutations in the bcsC gene during adaptive breeding of Komagataeibacter oboediens

Nonsense mutation in the bcsC gene occurred in the ethanol-adapted strain of Komagataeibacter oboediens MSKU 3, E3 strain, resulting in the loss of the function to produce BNC. In this study, we tried to restore the BNC-producing ability of E3 strain by the following adaptive mutation through repetitive static culture, and obtained four BNC-producing revertant strains, of which the bcsC gene had InDel mutations near the frameshift mutation region in E3 strain, resulting in several amino acid alterations compared with the BcsC of MSKU 3. Each revertant produced BNCs with different productivity on the static culture. Interestingly, one of the revertants, R37-9, produced BNC with a finer structure and narrower range of fibrils width, compared to others. The genome of R37-9 strain revealed only one amino acid substitution in the bcsC gene. Thus, we concluded that N713D mutation occurred in the bcsC gene is responsible for the finer fibrils structure.

Underestimated Noncovalent Interactions in Protein Data Bank

Noncovalent interactions (NCIs) play essential roles in the structure and function of biomacromolecules. There are various NCIs, e.g., hydrogen bonds (HBs), cation-π and π-π interactions, and ionic bonds, among which HBs are the most widespread and well-studied. By utilizing the ratio of the observed HBs over pseudo HBs (1.0 Å longer than the HB distance criteria without angle constraints), we demonstrated that HBs in both protein-ligand and protein-protein interfaces are overlooked in structures deposited in PDB. After the QM/MM optimization of 12 protein-ligand complexes, we showed that the overlooked HBs could be recovered. With a systematic search in the PDB, we found that the HB number per residue (N_HB/R) in proteins decreases as structural resolution becomes lower, implying that HBs are overlooked even today, regardless of the type of refinement approach used. Similarly, cation-π, π-π, and ionic interactions were found to be significantly lost, manifesting the universal underestimation of various NCIs. Considering the vital role of NCIs, it is important to recover the NCIs to facilitate drug design, to explore protein-protein interaction, and to study protein structure and function.

Harmonizing Experimental Data with Modeling to Predict Membrane Protein Insertion in Yeast

Membrane proteins must adopt their proper topologies within biological membranes, but achieving the correct topology is compromised by the presence of marginally hydrophobic transmembrane helices (TMHs). In this study, we report on a new model membrane protein in yeast that harbors two TMHs fused to an unstable nucleotide-binding domain. Because the second helix (TMH2) in this reporter has an unfavorable predicted free energy of insertion, we employed established methods to generate variants that alter TMH2 insertion free energy. We first found that altering TMH2 did not significantly affect the extent of protein degradation by the cellular quality control machinery. Next, we correlated predicted insertion free energies from a knowledge-based energy scale with the measured apparent free energies of TMH2 insertion. Although the predicted and apparent insertion energies showed a similar trend, the predicted free-energy changes spanned an unanticipated narrow range. By instead using a physics-based model, we obtained a broader range of free energies that agreed considerably better with the magnitude of the experimentally derived values. Nevertheless, some variants still inserted better in yeast than predicted from energy-based scales. Therefore, molecular dynamics simulations were performed and indicated that the corresponding mutations induced conformational changes within TMH2, which altered the number of stabilizing hydrogen bonds. Together, our results offer insight into the ability of the cellular quality control machinery to recognize conformationally distinct misfolded topomers, provide a model to assess TMH insertion in vivo, and indicate that TMH insertion energy scales may be limited depending on the specific protein and the mutation present.

Thermodynamics of Helix formation in small peptides of varying lengthin vacuo, implicit solvent and explicit solvent: Comparison between AMBER force fields

Helix formation is of great significance in protein folding. The helix-forming tendencies of amino acids are accumulated along the sequence to determine the helix-forming tendency of peptides. Computer simulation can be used to model this process in atomic details and give structural insights. In the current work, we employ equilibrate-state free energy simulation to systematically study the folding/unfolding thermodynamics of a series of mutated peptides. Two AMBER force fields including AMBER99SB and AMBER14SB are compared. The new 14SB force field uses refitted torsion parameters compared with 99SB and they share the same atomic charge scheme. We find that in vacuo the helix formation is mutation dependent, which reflects the different helix propensities of different amino acids. In general, there are helix formers, helix indifferent groups and helix breakers. The helical structure becomes more favored when the length of the sequence becomes longer, which arises from the formation of additional backbone hydrogen bonds in the lengthened sequence. Therefore, the helix indifferent groups and helix breakers will become helix formers in long sequences. Also, protonation-dependent helix formation is observed for ionizable groups. In 14SB, the helical structures are more stable than in 99SB and differences can be observed in their grouping schemes, especially in the helix indifferent group. In solvents, all mutations are helix indifferent due to protein–solvent interactions. The decrease in the number of backbone hydrogen bonds is the same with the increase in the number of protein–water hydrogen bonds. The 14SB in explicit solvent is able to capture the free energy minima in the helical state while 14SB in implicit solvent, 99SB in explicit solvent and 99SB in implicit solvent cannot. The helix propensities calculated under 14SB agree with the corresponding experimental values, while the 99SB results obviously deviate from the references. Hence, implicit solvent models are unable to correctly describe the thermodynamics even for the simple helix formation in isolated peptides. Well-developed force fields and explicit solvents are needed to correctly describe the protein dynamics. Aside from the free energy, differences in conformational ensemble under different force fields in different solvent models are observed. The numbers of hydrogen bonds formed under different force fields agree and they are mostly determined by the solvent model.

Modulating Hinge Flexibility in the APP Transmembrane Domain Alters γ-Secretase Cleavage

Intramembrane cleavage of the β-amyloid precursor protein C99 substrate by γ-secretase is implicated in Alzheimer's disease pathogenesis. Biophysical data have suggested that the N-terminal part of the C99 transmembrane domain (TMD) is separated from the C-terminal cleavage domain by a di-glycine hinge. Because the flexibility of this hinge might be critical for γ-secretase cleavage, we mutated one of the glycine residues, G38, to a helix-stabilizing leucine and to a helix-distorting proline. Both mutants impaired γ-secretase cleavage and also altered its cleavage specificity. Circular dichroism, NMR, and backbone amide hydrogen/deuterium exchange measurements as well as molecular dynamics simulations showed that the mutations distinctly altered the intrinsic structural and dynamical properties of the substrate TMD. Although helix destabilization and/or unfolding was not observed at the initial ε-cleavage sites of C99, subtle changes in hinge flexibility were identified that substantially affected helix bending and twisting motions in the entire TMD. These resulted in altered orientation of the distal cleavage domain relative to the N-terminal TMD part. Our data suggest that both enhancing and reducing local helix flexibility of the di-glycine hinge may decrease the occurrence of enzyme-substrate complex conformations required for normal catalysis and that hinge mobility can thus be conducive for productive substrate-enzyme interactions.

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.

Increased H-Bond Stability Relates to Altered ε-Cleavage Efficiency and Aβ Levels in the I45T Familial Alzheimer's Disease Mutant of APP

Cleavage of the amyloid precursor protein's (APP) transmembrane domain (TMD) by γ-secretase is a crucial step in the aetiology of Alzheimer's Disease (AD). Mutations in the APP TMD alter cleavage and lead to familial forms of AD (FAD). The majority of FAD mutations shift the preference of initial cleavage from ε49 to ε48, thus raising the AD-related Aβ42/Aβ40 ratio. The I45T mutation is among the few FAD mutations that do not alter ε-site preference, while it dramatically reduces the efficiency of ε-cleavage. Here, we investigate the impact of the I45T mutation on the backbone dynamics of the substrate TMD. Amide exchange experiments and molecular dynamics simulations in solvent and a lipid bilayer reveal an increased stability of amide hydrogen bonds at the ζ- and γ-cleavage sites. Stiffening of the H-bond network is caused by an additional H-bond between the T45 side chain and the TMD backbone, which alters dynamics within the cleavage domain. In particular, the increased H-bond stability inhibits an upward movement of the ε-sites in the I45T mutant. Thus, an altered presentation of ε-sites to the active site of γ-secretase as a consequence of restricted local flexibility provides a rationale for reduced ε-cleavage efficiency of the I45T mutant.

Improved free-energy landscape reconstruction of bacteriorhodopsin highlights local variations in unfolding energy

Precisely quantifying the energetics that drive the folding of membrane proteins into a lipid bilayer remains challenging. More than 15 years ago, atomic force microscopy (AFM) emerged as a powerful tool to mechanically extract individual membrane proteins from a lipid bilayer. Concurrently, fluctuation theorems, such as the Jarzynski equality, were applied to deduce equilibrium free energies (ΔG₀) from non-equilibrium single-molecule force spectroscopy records. The combination of these two advances in single-molecule studies deduced the free-energy of the model membrane protein bacteriorhodopsin in its native lipid bilayer. To elucidate this free-energy landscape at a higher resolution, we applied two recent developments. First, as an input to the reconstruction, we used force-extension curves acquired with a 100-fold higher time resolution and 10-fold higher force precision than traditional AFM studies of membrane proteins. Next, by using an inverse Weierstrass transform and the Jarzynski equality, we removed the free energy associated with the force probe and determined the molecular free-energy landscape of the molecule under study, bacteriorhodopsin. The resulting landscape yielded an average unfolding free energy per amino acid (aa) of 1.0 ± 0.1 kcal/mol, in agreement with past single-molecule studies. Moreover, on a smaller spatial scale, this high-resolution landscape also agreed with an equilibrium measurement of a particular three-aa transition in bacteriorhodopsin that yielded 2.7 kcal/mol/aa, an unexpectedly high value. Hence, while average unfolding ΔG₀ per aa is a useful metric, the derived high-resolution landscape details significant local variation from the mean. More generally, we demonstrated that, as anticipated, the inverse Weierstrass transform is an efficient means to reconstruct free-energy landscapes from AFM data.

Structure of APP-C991-99 and implications for role of extra-membrane domains in function and oligomerization

The 99 amino acid C-terminal fragment of Amyloid Precursor Protein APP-C99 (C99) is cleaved by γ-secretase to form Aβ peptide, which plays a critical role in the etiology of Alzheimer's Disease (AD). The structure of C99 consists of a single transmembrane domain flanked by intra and intercellular domains. While the structure of the transmembrane domain has been well characterized, little is known about the structure of the flanking domains and their role in C99 processing by γ-secretase. To gain insight into the structure of full-length C99, REMD simulations were performed for monomeric C99 in model membranes of varying thickness. We find equilibrium ensembles of C99 from simulation agree with experimentally-inferred residue insertion depths and protein backbone chemical shifts. In thin membranes, the transmembrane domain structure is correlated with extra-membrane structural states and the extra-membrane domain structural states become less correlated to each other. Mean and variance of the transmembrane and G₃₇G₃₈ hinge angles are found to increase with thinning membrane. The N-terminus of C99 forms β-strands that may seed aggregation of Aβ on the membrane surface, promoting amyloid formation. In thicker membranes the N-terminus forms α-helices that interact with the nicastrin domain of γ-secretase. The C-terminus of C99 becomes more α-helical as the membrane thickens, forming structures that may be suitable for binding by cytoplasmic proteins, while C-terminal residues essential to cytotoxic function become α-helical as the membrane thins. The heterogeneous but discrete extra-membrane domain states analyzed here open the path to new investigations of the role of C99 structure and membrane in amyloidogenesis. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.