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.

Crystal structure of mevalonate 3,5-bisphosphate decarboxylase reveals insight into the evolution of decarboxylases in the mevalonate metabolic pathways.

Mevalonate 3,5-bisphosphate decarboxylase is involved in the recently discovered Thermoplasma-type mevalonate pathway. The enzyme catalyzes the elimination of the 3-phosphate group from mevalonate 3,5-bisphosphate as well as concomitant decarboxylation of the substrate. This entire reaction of the enzyme resembles the latter half-reactions of its homologs, diphosphomevalonate decarboxylase and phosphomevalonate decarboxylase, which also catalyze ATP-dependent phosphorylation of the 3-hydroxyl group of their substrates. However, the crystal structure of mevalonate 3,5-bisphosphate decarboxylase and the structural reasons of the difference between reactions catalyzed by the enzyme and its homologs are unknown. In this study, we determined the X-ray crystal structure of mevalonate 3,5-bisphosphate decarboxylase from Picrophilus torridus, a thermoacidophilic archaeon of the order Thermoplasmatales. Structural and mutational analysis demonstrated the importance of a conserved aspartate residue for enzyme activity. In addition, although crystallization was performed in the absence of substrate or ligands, residual electron density having the shape of a fatty acid was observed at a position overlapping the ATP-binding site of the homologous enzyme, diphosphomevalonate decarboxylase. This finding is in agreement with the expected evolutionary route from phosphomevalonate decarboxylase (ATP-dependent) to mevalonate 3,5-bisphosphate decarboxylase (ATP-independent) through the loss of kinase activity. We found that the binding of geranylgeranyl diphosphate, an intermediate of the archeal isoprenoid biosynthesis pathway, evoked significant activation of mevalonate 3,5-bisphosphate decarboxylase, and several mutations at the putative geranylgeranyl diphosphate-binding site impaired this activation, suggesting the physiological importance of ligand binding as well as a possible novel regulatory system employed by the Thermoplasma-type mevalonate pathway.

A bio-inspired cell-free system for cannabinoid production from inexpensive inputs.

Moving cannabinoid production away from the vagaries of plant extraction and into engineered microbes could provide a consistent, purer, cheaper and environmentally benign source of these important therapeutic molecules, but microbial production faces notable challenges. An alternative to microbes and plants is to remove the complexity of cellular systems by employing enzymatic biosynthesis. Here we design and implement a new cell-free system for cannabinoid production with the following features: (1) only low-cost inputs are needed; (2) only 12 enzymes are employed; (3) the system does not require oxygen and (4) we use a nonnatural enzyme system to reduce ATP requirements that is generally applicable to malonyl-CoA-dependent pathways such as polyketide biosynthesis. The system produces ~0.5 g l-1 cannabigerolic acid (CBGA) or cannabigerovarinic acid (CBGVA) from low-cost inputs, nearly two orders of magnitude higher than yeast-based production. Cell-free systems such as this may provide a new route to reliable cannabinoid production.

Unfolding of a ClC chloride transporter retains memory of its evolutionary history.

ClC chloride channels and transporters are important for chloride homeostasis in species from bacteria to human. Mutations in ClC proteins cause genetically inherited diseases, some of which are likely to involve folding defects. The ClC proteins present a challenging and unusual biological folding problem because they are large membrane proteins possessing a complex architecture, with many reentrant helices that go only partway through membrane and loop back out. Here we were able to examine the unfolding of the Escherichia coli ClC transporter, ClC-ec1, using single-molecule forced unfolding methods. We found that the protein could be separated into two stable halves that unfolded independently. The independence of the two domains is consistent with an evolutionary model in which the two halves arose from independently folding subunits that later fused together. Maintaining smaller folding domains of lesser complexity within large membrane proteins may be an advantageous strategy to avoid misfolding traps.

Exome sequencing and CRISPR/Cas genome editing identify mutations of ZAK as a cause of limb defects in humans and mice.

The CRISPR/Cas technology enables targeted genome editing and the rapid generation of transgenic animal models for the study of human genetic disorders. Here we describe an autosomal recessive human disease in two unrelated families characterized by a split-foot defect, nail abnormalities of the hands, and hearing loss, due to mutations disrupting the SAM domain of the protein kinase ZAK. ZAK is a member of the MAPKKK family with no known role in limb development. We show that Zak is expressed in the developing limbs and that a CRISPR/Cas-mediated knockout of the two Zak isoforms is embryonically lethal in mice. In contrast, a deletion of the SAM domain induces a complex hindlimb defect associated with down-regulation of Trp63, a known split-hand/split-foot malformation disease gene. Our results identify ZAK as a key player in mammalian limb patterning and demonstrate the rapid utility of CRISPR/Cas genome editing to assign causality to human mutations in the mouse in <10 wk.

A molecular rheostat maintains ATP levels to drive a synthetic biochemistry system.

Synthetic biochemistry seeks to engineer complex metabolic pathways for chemical conversions outside the constraints of the cell. Establishment of effective and flexible cell-free systems requires the development of simple systems to replace the intricate regulatory mechanisms that exist in cells for maintaining high-energy cofactor balance. Here we describe a simple rheostat that regulates ATP levels by controlling the flow down either an ATP-generating or non-ATP-generating pathway according to the free-phosphate concentration. We implemented this concept for the production of isobutanol from glucose. The rheostat maintains adequate ATP concentrations even in the presence of ATPase contamination. The final system including the rheostat produced 24.1 ± 1.8 g/L of isobutanol from glucose in 91% theoretical yield with an initial productivity of 1.3 g/L/h. The molecular rheostat concept can be used in the design of continuously operating, self-sustaining synthetic biochemistry systems.

A synthetic biochemistry module for production of bio-based chemicals from glucose.

Synthetic biochemistry, the cell-free production of biologically based chemicals, is a potentially high-yield, flexible alternative to in vivo metabolic engineering. To limit costs, cell-free systems must be designed to operate continuously with minimal addition of feedstock chemicals. We describe a robust, efficient synthetic glucose breakdown pathway and implement it for the production of bioplastic. The system's performance suggests that synthetic biochemistry has the potential to become a viable industrial alternative.

Molecular dynamics simulation strategies for protein-micelle complexes.

The structure and stability of membrane proteins can vary widely in different detergents and this variability has great practical consequences for working with membrane proteins. Nevertheless, the mechanisms that operate to alter the behavior of proteins in micelles are poorly understood and not predictable. Atomic simulations could provide considerable insight into these mechanisms. Building protein-micelle complexes for simulation is fraught with uncertainty, however, in part because it is often unknown how many detergent molecules are present in the complex. Here, we describe several convenient ways to employ Micelle Builder in CHARMM-GUI to rapidly construct protein-micelle complexes and performed simulations of the isolated voltage-sensor domain of voltage-dependent potassium-selective channel and an antimicrobial peptide papiliocin with varying numbers of detergents. We found that once the detergent number exceeds a threshold, protein-detergent interactions change very little and remain very consistent with experimental observations. Our results provide a platform for future studies of the interplays between protein structure and detergent properties at the atomic level. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Backbone Hydrogen Bond Strengths Can Vary Widely in Transmembrane Helices.

Although backbone hydrogen bonds in transmembrane (TM) helices have the potential to be very strong due to the low dielectric and low water environment of the membrane, their strength has never been assessed experimentally. Moreover, variations in hydrogen bond strength might be necessary to facilitate the TM helix breaking and bending that is often needed to satisfy functional imperatives. Here we employed equilibrium hydrogen/deuterium fractionation factors to measure backbone hydrogen bond strengths in the TM helix of the amyloid precursor protein (APP). We find an enormous range of hydrogen bond free energies, with some weaker than water-water hydrogen bonds and some over 6 kcal/mol stronger than water-water hydrogen bonds. We find that weak hydrogen bonds are at or near preferred γ-secretase cleavage sites, suggesting that the sequence of APP and possibly other cleaved TM helices may be designed, in part, to make their backbones accessible for cleavage. The finding that hydrogen bond strengths in a TM helix can vary widely has implications for membrane protein function, dynamics, evolution, and design.

Mapping the energy landscape for second-stage folding of a single membrane protein.

Membrane proteins are designed to fold and function in a lipid membrane, yet folding experiments within a native membrane environment are challenging to design. Here we show that single-molecule forced unfolding experiments can be adapted to study helical membrane protein folding under native-like bicelle conditions. Applying force using magnetic tweezers, we find that a transmembrane helix protein, Escherichia coli rhomboid protease GlpG, unfolds in a highly cooperative manner, largely unraveling as one physical unit in response to mechanical tension above 25 pN. Considerable hysteresis is observed, with refolding occurring only at forces below 5 pN. Characterizing the energy landscape reveals only modest thermodynamic stability (ΔG = 6.5 kBT) but a large unfolding barrier (21.3 kBT) that can maintain the protein in a folded state for long periods of time (t1/2 ∼3.5 h). The observed energy landscape may have evolved to limit the existence of troublesome partially unfolded states and impart rigidity to the structure.

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.

Shifting hydrogen bonds may produce flexible transmembrane helices.

The intricate functions of membrane proteins would not be possible without bends or breaks that are remarkably common in transmembrane helices. The frequent helix distortions are nevertheless surprising because backbone hydrogen bonds should be strong in an apolar membrane, potentially rigidifying helices. It is therefore mysterious how distortions can be generated by the evolutionary currency of random point mutations. Here we show that we can engineer a transition between distinct distorted helix conformations in bacteriorhodopsin with a single-point mutation. Moreover, we estimate the energetic cost of the conformational transitions to be smaller than 1 kcal/mol. We propose that the low energy of distortion is explained in part by the shifting of backbone hydrogen bonding partners. Consistent with this view, extensive backbone hydrogen bond shifts occur during helix conformational changes that accompany functional cycles. Our results explain how evolution has been able to liberally exploit transmembrane helix bending for the optimization of membrane protein structure, function, and dynamics.

Evidence of a novel mevalonate pathway in archaea.

Isoprenoids make up a remarkably diverse class of more than 25000 biomolecules that include familiar compounds such as cholesterol, chlorophyll, vitamin A, ubiquinone, and natural rubber. The two essential building blocks of all isoprenoids, isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP), are ubiquitous in the three domains of life. In most eukaryotes and archaea, IPP and DMAPP are generated through the mevalonate pathway. We have identified two novel enzymes, mevalonate-3-kinase and mevalonate-3-phosphate-5-kinase from Thermoplasma acidophilum, which act sequentially in a putative alternate mevalonate pathway. We propose that a yet unidentified ATP-independent decarboxylase acts upon mevalonate 3,5-bisphosphate, yielding isopentenyl phosphate, which is subsequently phosphorylated by the known isopentenyl phosphate kinase from T. acidophilum to generate the universal isoprenoid precursor, IPP.

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.

Regulation of clathrin adaptor function in endocytosis: novel role for the SAM domain.

During clathrin-mediated endocytosis, adaptor proteins play central roles in coordinating the assembly of clathrin coats and cargo selection. Here we characterize the binding of the yeast endocytic adaptor Sla1p to clathrin through a variant clathrin-binding motif that is negatively regulated by the Sla1p SHD2 domain. The crystal structure of SHD2 identifies the domain as a sterile alpha-motif (SAM) domain and shows a propensity to oligomerize. By co-immunoprecipitation, Sla1p binds to clathrin and self-associates in vivo. Mutations in the clathrin-binding motif that abolish clathrin binding and structure-based mutations in SHD2 that impede self-association result in endocytosis defects and altered dynamics of Sla1p assembly at the sites of endocytosis. These results define a novel mechanism for negative regulation of clathrin binding by an adaptor and suggest a role for SAM domains in clathrin-mediated endocytosis.

Characterization of the SAM domain of the PKD-related protein ANKS6 and its interaction with ANKS3.

BACKGROUND: Autosomal dominant polycystic kidney disease (ADPKD) is the most common genetic disorder leading to end-stage renal failure in humans. In the PKD/Mhm(cy/+) rat model of ADPKD, the point mutation R823W in the sterile alpha motif (SAM) domain of the protein ANKS6 is responsible for disease. SAM domains are known protein-protein interaction domains, capable of binding each other to form polymers and heterodimers. Despite its physiological importance, little is known about the function of ANKS6 and how the R823W point mutation leads to PKD. Recent work has revealed that ANKS6 interacts with a related protein called ANKS3. Both ANKS6 and ANKS3 have a similar domain structure, with ankyrin repeats at the N-terminus and a SAM domain at the C-terminus. RESULTS: The SAM domain of ANKS3 is identified as a direct binding partner of the ANKS6 SAM domain. We find that ANKS3-SAM polymerizes and ANKS6-SAM can bind to one end of the polymer. We present crystal structures of both the ANKS3-SAM polymer and the ANKS3-SAM/ANKS6-SAM complex, revealing the molecular details of their association. We also learn how the R823W mutation disrupts ANKS6 function by dramatically destabilizing the SAM domain such that the interaction with ANKS3-SAM is lost. CONCLUSIONS: ANKS3 is a direct interacting partner of ANKS6. By structurally and biochemically characterizing the interaction between the ANKS3 and ANKS6 SAM domains, our work provides a basis for future investigation of how the interaction between these proteins mediates kidney function.

Improving the tolerance of Escherichia coli to medium-chain fatty acid production.

Microbial fatty acids are an attractive source of precursors for a variety of renewable commodity chemicals such as alkanes, alcohols, and biofuels. Rerouting lipid biosynthesis into free fatty acid production can be toxic, however, due to alterations of membrane lipid composition. Here we find that membrane lipid composition can be altered by the direct incorporation of medium-chain fatty acids into lipids via the Aas pathway in cells expressing the medium-chain thioesterase from Umbellularia californica (BTE). We find that deletion of the aas gene and sequestering exported fatty acids reduces medium-chain fatty acid toxicity, partially restores normal lipid composition, and improves medium-chain fatty acid yields.

Rampant exchange of the structure and function of extramembrane domains between membrane and water soluble proteins.

Of the membrane proteins of known structure, we found that a remarkable 67% of the water soluble domains are structurally similar to water soluble proteins of known structure. Moreover, 41% of known water soluble protein structures share a domain with an already known membrane protein structure. We also found that functional residues are frequently conserved between extramembrane domains of membrane and soluble proteins that share structural similarity. These results suggest membrane and soluble proteins readily exchange domains and their attendant functionalities. The exchanges between membrane and soluble proteins are particularly frequent in eukaryotes, indicating that this is an important mechanism for increasing functional complexity. The high level of structural overlap between the two classes of proteins provides an opportunity to employ the extensive information on soluble proteins to illuminate membrane protein structure and function, for which much less is known. To this end, we employed structure guided sequence alignment to elucidate the functions of membrane proteins in the human genome. Our results bridge the gap of fold space between membrane and water soluble proteins and provide a resource for the prediction of membrane protein function. A database of predicted structural and functional relationships for proteins in the human genome is provided at sbi.postech.ac.kr/emdmp.

Modest stabilization by most hydrogen-bonded side-chain interactions in membrane proteins.

Understanding the energetics of molecular interactions is fundamental to all of the central quests of structural biology including structure prediction and design, mapping evolutionary pathways, learning how mutations cause disease, drug design, and relating structure to function. Hydrogen-bonding is widely regarded as an important force in a membrane environment because of the low dielectric constant of membranes and a lack of competition from water. Indeed, polar residue substitutions are the most common disease-causing mutations in membrane proteins. Because of limited structural information and technical challenges, however, there have been few quantitative tests of hydrogen-bond strength in the context of large membrane proteins. Here we show, by using a double-mutant cycle analysis, that the average contribution of eight interhelical side-chain hydrogen-bonding interactions throughout bacteriorhodopsin is only 0.6 kcal mol(-1). In agreement with these experiments, we find that 4% of polar atoms in the non-polar core regions of membrane proteins have no hydrogen-bond partner and the lengths of buried hydrogen bonds in soluble proteins and membrane protein transmembrane regions are statistically identical. Our results indicate that most hydrogen-bond interactions in membrane proteins are only modestly stabilizing. Weak hydrogen-bonding should be reflected in considerations of membrane protein folding, dynamics, design, evolution and function.

Thermodynamic stability of bacteriorhodopsin mutants measured relative to the bacterioopsin unfolded state.

The stability of bacteriorhodopsin (bR) has often been assessed using SDS unfolding assays that monitor the transition of folded bR (bR(f)) to unfolded (bR(u)). While many criteria suggest that the unfolding curves reflect thermodynamic stability, slow retinal (RET) hydrolysis during refolding makes it impossible to perform the most rigorous test for equilibrium, i.e., superimposable unfolding and refolding curves. Here we made a new equilibrium test by asking whether the refolding rate in the transition zone is faster than RET hydrolysis. We find that under conditions we have used previously, refolding is in fact slower than hydrolysis, strongly suggesting that equilibrium is not achieved. Instead, the apparent free energy values reported previously are dominated by unfolding rates. To assess how different the true equilibrium values are, we employed an alternative method by measuring the transition of bR(f) to unfolded bacterioopsin (bO(u)), the RET-free form of unfolded protein. The bR(f)-to-bO(u) transition is fully reversible, particular when we add excess RET. We compared the difference in unfolding free energies for 13 bR mutants measured by both assays. For 12 of the 13 mutants with a wide range of stabilities, the results are essentially the same within experimental error. The congruence of the results is fortuitous and suggests the energetic effects of most mutations may be focused on the folded state. The bR(f)-to-bO(u) reaction is inconvenient because many days are required to reach equilibrium, but it is the preferable measure of thermodynamic stability. This article is part of a Special Issue entitled: Protein Folding in Membranes.