From the test tube to the cell: exploring the folding and aggregation of a beta-clam protein

From Rat Tails to Glycoproteostasis: Motivated by Biology, Enabled by Biophysics, and Lucky

In this article I tell the story of my career path and how I have come to focus my research on protein folding in the cell. My early fascination with protein folding began during my undergraduate research. My graduate work exploited reductionist approaches to explore structural features in proteins by using cyclic peptide models ofβ-turns. My career trajectory from these early days to present, described in the first section of this article, illustrates the importance of pursuing the scientific questions that one finds most exciting and seizing professional opportunities that enable these questions to be tackled productively. In addition, this trajectory shows how serendipity can shape a career path. The second section describes the extraordinary scientific discoveries I have witnessed in protein folding during my career. Here I explain how I was drawn into the world of protein folding in thecell. This turning point allowed me to participate in the explosion of research on molecular chaperones in the early 90's and to help elucidate the nature of chaperone-substrate recognition, a problem I continue to focus on. Examples of our research contributions are presented in the third section, with a perspective on major challenges for the future offered in the last section. Throughout my career I have engaged in many collaborations;each has opened new scientific doors. Importantly, I seek to instill in my trainees the same excitement about research that I feel and to foster their growth as scientists and their discovery of their own passions and talents.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Structure and molecular evolution of the barcode fragment of cytochrome oxidase I (COI) in Macrocheles (Acari: Mesostigmata: Macrochelidae)

Consisting of approximately 320 species, Macrocheles is the most widely distributed genus in the family Macrochelidae. Though some studies have focused on the description of Macrochelidae using molecular techniques (e.g., RAPD) and sequencing of some genes, the interspecies relationships within Macrocheles still remain uncertain. As such, in the present study, we examine all publicly available data in GenBank to explore the evolutionary relationships, divergence times, and amino acid variations within Macrocheles. Exploring the patterns of variation in the secondary protein structure shows high levels of conservation in the second and last helices, emphasizing their involvement in the energy metabolism function of the cytochrome oxidase subunit I enzyme. According to our phylogenetic analysis, all available Macrocheles species are clustered in a monophyletic group. However, in the reconstructed trees, we subdivided M. merdarius and M. willowae into two well-supported intraspecific clades that are driven by geographic separation and host specificity. We also estimate the divergence time of selected species using calibration evidence from available fossils and previous studies. Thus, we estimate that the age of the Parasitiformes is 320.4 (273.3-384.3) Mya (Permian), and the Mesostigmata is 285.1 (270.8-286.4) Mya (Carboniferous), both with likely origins in the Paleozoic era. We also estimate that Macrocheles diverged from other Mesostigmata mites during the Mesozoic, approximately 222.9 Mya.

Protein stability [determination] problems

Human health depends on the correct folding of proteins, for misfolding and aggregation lead to diseases. An unfolded (denatured) protein can refold to its original folded state. How does this occur is known as the protein folding problem. One of several related questions to this problem is that how much more stable is the folded state than the unfolded state. There are several measures of protein stability. In this article, protein stability is given a thermodynamic definition and is measured by Gibbs free energy change ( Δ G D 0 ) associated with the equilibrium, native (N) conformation ↔ denatured (D) conformation under the physiological condition usually taken as dilute buffer (or water) at 25 °C. We show that this thermodynamic quantity ( Δ G D 0 ), where subscript D represents transition between N and D states, and superscript 0 (zero) represents the fact that the transition occurs in the absence of denaturant, can be neither measured nor predicted under physiological conditions. However, Δ G D can be measured in the presence of strong chemical denaturants such as guanidinium chloride and urea which are shown to destroy all noncovalent interactions responsible for maintaining the folded structure. A problem with this measurement is that the estimate of Δ G D 0 comes from the analysis of the plot of Δ G D versus denaturant concentration, which requires a long extrapolation of values of Δ G D , and all the three methods of extrapolation give three different values of Δ G D 0 for a protein. Thus, our confidence in the authentic value of Δ G D 0 is eroded. Another problem with this in vitro measurement of Δ G D 0 is that it is done on the pure protein sample in dilute buffer which is a very large extrapolation of the in vivo conditions, for the crowding effect on protein stability is ignored.

On the Effects of Disordered Tails, Supertertiary Structure and Quinary Interactions on the Folding and Function of Protein Domains

The vast majority of our current knowledge about the biochemical and biophysical properties of proteins derives from in vitro studies conducted on isolated globular domains. However, a very large fraction of the proteins expressed in the eukaryotic cell are structurally more complex. In particular, the discovery that up to 40% of the eukaryotic proteins are intrinsically disordered, or possess intrinsically disordered regions, and are highly dynamic entities lacking a well-defined three-dimensional structure, revolutionized the structure–function paradigm and our understanding of proteins. Moreover, proteins are mostly characterized by the presence of multiple domains, influencing each other by intramolecular interactions. Furthermore, proteins exert their function in a crowded intracellular milieu, transiently interacting with a myriad of other macromolecules. In this review we summarize the literature tackling these themes from both the theoretical and experimental perspectives, highlighting the effects on protein folding and function that are played by (i) flanking disordered tails; (ii) contiguous protein domains; (iii) interactions with the cellular environment, defined as quinary structures. We show that, in many cases, both the folding and function of protein domains is remarkably perturbed by the presence of these interactions, pinpointing the importance to increase the level of complexity of the experimental work and to extend the efforts to characterize protein domains in more complex contexts.

Computational Models for the Study of Protein Aggregation

Protein aggregation has been studied by many groups around the world for many years because it can be the cause of a number of neurodegenerative diseases that have no effective treatment. Obtaining the structure of related fibrils and toxic oligomers, as well as describing the pathways and main factors that govern the self-organization process, is of paramount importance, but it is also very difficult. To solve this problem, experimental and computational methods are often combined to get the most out of each method. The effectiveness of the computational approach largely depends on the construction of a reasonable molecular model. Here we discussed different versions of the four most popular all-atom force fields AMBER, CHARMM, GROMOS, and OPLS, which have been developed for folded and intrinsically disordered proteins, or both. Continuous and discrete coarse-grained models, which were mainly used to study the kinetics of aggregation, are also summarized.

An "Onion-like" Model of Protein Unfolding: Collective versus Site Specific Approaches

Approximating protein unfolding by an all-or-none cooperative event is a convenient assumption that can provide precious global information on protein stability. It is however quickly emerging that the scenario is far more complex and that global denaturation curves often hide a rich heterogeneity of states that are largely probe dependent. In this review, we revisit the importance of gaining site-specific information on the unfolding process. We focus on nuclear magnetic resonance, as this is the main technique able to provide site-specific information. We review historical and most modern approaches that have allowed an appreciable advancement of the field of protein folding. We also demonstrate how unfolding is a reporter dependent event, suggesting the outmost importance of selecting the reporter carefully.

Specific Interactions and Environment Flexibility Tune Protein Stability under Extreme Crowding

Macromolecular crowding influences protein mobility and stability in vivo. A precise description of the crowding effect on protein thermal stability requires the estimate of the combined effects of excluded volume, specific protein-environment interactions, as well as the thermal response of the crowders. Here, we explore an ideal model system, the lysozyme protein in powder state, to dissect the factors controlling the melting of the protein under extreme crowding. By deploying state-of-the art molecular simulations, supported by calorimetric experiments, we assess the role of the environment flexibility and of intermolecular electrostatic interactions. In particular, we show that the temperature-dependent flexibility of the macromolecular crowders, along with specific interactions, significantly alleviates the stabilizing contributions of the static volume effect.

Cleavage-Resistant Protein Labeling With Hydrophilic Trityl Enables Distance Measurements In-Cell

Sensitive in-cell distance measurements in proteins using pulsed-electron spin resonance (ESR) require reduction-resistant and cleavage-resistant spin labels. Among the reduction-resistant moieties, the hydrophilic trityl core known as OX063 is promising due to its long phase-memory relaxation time (T_m). This property leads to a sufficiently intense ESR signal for reliable distance measurements. Furthermore, the T_m of OX063 remains sufficiently long at higher temperatures, opening the possibility for measurements at temperatures above 50 K. In this work, we synthesized deuterated OX063 with a maleimide linker (mOX063-d₂₄). We show that the combination of the hydrophilicity of the label and the maleimide linker enables high protein labeling that is cleavage-resistant in-cells. Distance measurements performed at 150 K using this label are more sensitive than the measurements at 80 K. The sensitivity gain is due to the significantly short longitudinal relaxation time (T₁) at higher temperatures, which enables more data collection per unit of time. In addition to in vitro experiments, we perform distance measurements in Xenopus laevis oocytes. Interestingly, the T_m of mOX063-d₂₄ is sufficiently long even in the crowded environment of the cell, leading to signals of appreciable intensity. Overall, mOX063-d₂₄ provides highly sensitive distance measurements both in vitro and in-cells.

Polyanions Cause Protein Destabilization Similar to That in Live Cells

The structural stability of proteins is found to markedly change upon their transfer to the crowded interior of live cells. For some proteins, the stability increases, while for others, it decreases, depending on both the sequence composition and the type of host cell. The mechanism seems to be linked to the strength and conformational bias of the diffusive in-cell interactions, where protein charge is found to play a decisive role. Because most proteins, nucleotides, and membranes carry a net-negative charge, the intracellular environment behaves like a polyanionic (Z:1) system with electrostatic interactions different from those of standard 1:1 ion solutes. To determine how such polyanion conditions influence protein stability, we use negatively charged polyacetate ions to mimic the net-negatively charged cellular environment. The results show that, per Na+ equivalent, polyacetate destabilizes the model protein SOD1barrel significantly more than monoacetate or NaCl. At an equivalent of 100 mM Na+, the polyacetate destabilization of SOD1barrel is similar to that observed in live cells. By the combined use of equilibrium thermal denaturation, folding kinetics, and high-resolution nuclear magnetic resonance, this destabilization is primarily assigned to preferential interaction between polyacetate and the globally unfolded protein. This interaction is relatively weak and involves mainly the outermost N-terminal region of unfolded SOD1barrel. Our findings point thus to a generic influence of polyanions on protein stability, which adds to the sequence-specific contributions and needs to be considered in the evaluation of in vivo data.

Polyethylene glycol perturbs the unfolding of CRABP I: A correlation between experimental and theoretical approach

The importance of macromolecules paves the way towards a detailed molecular level investigation as all most all cellular processes occurring at the interior of cells in the form of proteins, enzymes, and other biological molecules are significantly affected because of their crowding. Thus, exploring the role of crowding environment on the denaturation and renaturation kinetics of protein molecules is of great importance. Here, CRABP I (cellular retinoic acid binding protein I) is employed as a model protein along with different molecular weights of Polyethylene glycol (PEG) as molecular crowders. The experimental evaluations are done by accessing the protein secondary structure analysis using circular dichroism (CD) spectroscopy and unfolding kinetics using intrinsic fluorescence of CRABP I at 37 °C to mimic the in vivo crowding environment. The unfolding kinetics results indicated that both PEG 2000 and PEG 4000 act as stabilizers by retarding the unfolding kinetic rates. Both kinetic and stability outcomes presented the importance of crowding environment on stability and kinetics of CRABP I. The molecular dynamics (MD) studies revealed that thirteen PEG 2000 molecules assembled during the 500 ns simulation, which increases the stability and percentage of β-sheet. The experimental findings are well supported by the molecular dynamics simulation results.

Insights into Protein Stability in Cell Lysate by 19 F NMR Spectroscopy

In living organisms, protein folding and function take place in an inhomogeneous, highly crowded environment possessing a concentration of diverse macromolecules of up to 400 g/L. It has been shown that the intracellular environment has a pronounced effect on the stability, dynamics and function of the protein under study, and has for this reason to be considered. However, most protein studies neglect the presence of these macromolecules. Consequently, we probe here the overall thermodynamic stability of cold shock protein B from Bacillus subtilis (BsCspB) in cell lysate. We found that an increase in cell lysate concentration causes a monotonic increase in the thermodynamic stability of BsCspB. This result strongly underlines the importance of considering the biological environment when inherent protein parameters are quantitatively determined. Moreover, we demonstrate that targeted application of 19 F NMR spectroscopy operates as an ideal tool for protein studies performed in complex cellular surroundings.

Proteins: "Boil 'Em, Mash 'Em, Stick 'Em in a Stew"

Cells of the vast majority of organisms are subject to temperature, pressure, pH, ionic strength, and other stresses. We discuss these effects in the light of protein folding and protein interactions in vitro, in complex environments, in cells, and in vivo. Protein phase diagrams provide a way of organizing different structural ensembles that occur under stress and how one can move among ensembles. Experiments that perturb biomolecules in vitro or in cells by stressing them have revealed much about the underlying forces that are competing to control protein stability, folding, and function. Two phenomena that emerge and serve to broadly classify effects of the cellular environment are crowding (mainly due to repulsive forces) and sticking (mainly due to attractive forces). The interior of cells is closely balanced between these emergent effects, and stress can tip the balance one way or the other. The free energy scale involved is small but significant on the scale of the "on/off switches" that control signaling in cells or of protein-protein association with a favorable function such as increased enzyme processivity. Quantitative tools from biophysical chemistry will play an important role in elucidating the world of crowding and sticking under stress.

Whole-Cell Models and Simulations in Molecular Detail

Comprehensive data about the composition and structure of cellular components have enabled the construction of quantitative whole-cell models. While kinetic network-type models have been established, it is also becoming possible to build physical, molecular-level models of cellular environments. This review outlines challenges in constructing and simulating such models and discusses near- and long-term opportunities for developing physical whole-cell models that can connect molecular structure with biological function.

Protecting activity of desiccated enzymes

Protein-based biological drugs and many industrial enzymes are unstable, making them prohibitively expensive. Some can be stabilized by formulation with excipients, but most still require low temperature storage. In search of new, more robust excipients, we turned to the tardigrade, a microscopic animal that synthesizes cytosolic abundant heat soluble (CAHS) proteins to protect its cellular components during desiccation. We find that CAHS proteins protect the test enzymes lactate dehydrogenase and lipoprotein lipase against desiccation-, freezing-, and lyophilization-induced deactivation. Our data also show that a variety of globular and disordered protein controls, with no known link to desiccation tolerance, protect our test enzymes. Protection of lactate dehydrogenase correlates, albeit imperfectly, with the charge density of the protein additive, suggesting an approach to tune protection by modifying charge. Our results support the potential use of CAHS proteins as stabilizing excipients in formulations and suggest that other proteins may have similar potential.

Bacterial inclusion bodies are industrially exploitable amyloids

Understanding the structure, functionalities and biology of functional amyloids is an issue of emerging interest. Inclusion bodies, namely protein clusters formed in recombinant bacteria during protein production processes, have emerged as unanticipated, highly tunable models for the scrutiny of the physiology and architecture of functional amyloids. Based on an amyloidal skeleton combined with varying amounts of native or native-like protein forms, bacterial inclusion bodies exhibit an unusual arrangement that confers mechanical stability, biological activity and conditional protein release, being thus exploitable as versatile biomaterials. The applicability of inclusion bodies in biotechnology as enriched sources of protein and reusable catalysts, and in biomedicine as biocompatible topographies, nanopills or mimetics of endocrine secretory granules has been largely validated. Beyond these uses, the dissection of how recombinant bacteria manage the aggregation of functional protein species into structures of highly variable complexity offers insights about unsuspected connections between protein quality (conformational status compatible with functionality) and cell physiology.

Dynamics of proteins aggregation. II. Dynamic scaling in confined media

In this paper, the second in a series devoted to molecular modeling of protein aggregation, a mesoscale model of proteins together with extensive discontinuous molecular dynamics simulation is used to study the phenomenon in a confined medium. The medium, as a model of a crowded cellular environment, is represented by a spherical cavity, as well as cylindrical tubes with two aspect ratios. The aggregation process leads to the formation of β sheets and eventually fibrils, whose deposition on biological tissues is believed to be a major factor contributing to many neuro-degenerative diseases, such as Alzheimer's, Parkinson's, and amyotrophic lateral sclerosis diseases. Several important properties of the aggregation process, including dynamic evolution of the total number of the aggregates, the mean aggregate size, and the number of peptides that contribute to the formation of the β sheets, have been computed. We show, similar to the unconfined media studied in Paper I [S. Zheng et al., J. Chem. Phys. 145, 134306 (2016)], that the computed properties follow dynamic scaling, characterized by power laws. The existence of such dynamic scaling in unconfined media was recently confirmed by experiments. The exponents that characterize the power-law dependence on time of the properties of the aggregation process in spherical cavities are shown to agree with those in unbounded fluids at the same protein density, while the exponents for aggregation in the cylindrical tubes exhibit sensitivity to the geometry of the system. The effects of the number of amino acids in the protein, as well as the size of the confined media, have also been studied. Similarities and differences between aggregation in confined and unconfined media are described, including the possibility of no fibril formation, if confinement is severe.

Specific Ion Effects on Protein Thermal Aggregation from Dilute Solutions to Crowded Environments

We have investigated specific ion effects on protein thermal aggregation from dilute solutions to crowded environments. Ovalbumin and poly(ethylene glycol) have been employed as the model protein and crowding agent, respectively. Our studies demonstrate that the rate-limiting step of ovalbumin thermal aggregation is changed from the aggregation of unfolded protein molecules to the unfolding of the protein molecules, when the solution conditions are varied from a dilute solution to a crowded environment. The specific ion effects acting on the thermal aggregation of ovalbumin generated by kosmotropic and chaotropic ions are different. The thermal aggregation of ovalbumin molecules is promoted by kosmotropic anions in dilute solutions via an increase in protein hydrophobic interactions. In contrast, ovalbumin thermal aggregation is facilitated by chaotropic ions in crowded environments through accelerated unfolding of protein molecules. Therefore, there are distinct mechanisms causing the ion specificities of protein thermal aggregation between dilute solutions and crowded environments. The ion specificities are dominated by ion-specific hydrophobic interactions between protein molecules and ion-specific unfolding of protein molecules in dilute solutions and crowded environments, respectively.

Directed evolution to improve protein folding in vivo

Recently, several innovative approaches have been developed that allow one to directly screen or select for improved protein folding in the cellular context. These methods have the potential of not just leading to a better understanding of the in vivo folding process, they may also allow for improved production of proteins of biotechnological interest.

Enzyme Efficiency but Not Thermostability Drives Cefotaxime Resistance Evolution in TEM-1 β-Lactamase

A leading intellectual challenge in evolutionary genetics is to identify the specific phenotypes that drive adaptation. Enzymes offer a particularly promising opportunity to pursue this question, because many enzymes' contributions to organismal fitness depend on a comparatively small number of experimentally accessible properties. Moreover, on first principles the demands of enzyme thermostability stand in opposition to the demands of catalytic activity. This observation, coupled with the fact that enzymes are only marginally thermostable, motivates the widely held hypothesis that mutations conferring functional improvement require compensatory mutations to restore thermostability. Here, we explicitly test this hypothesis for the first time, using four missense mutations in TEM-1 β-lactamase that jointly increase cefotaxime Minimum Inhibitory Concentration (MIC) ∼1500-fold. First, we report enzymatic efficiency (kcat/KM) and thermostability (Tm, and thence ΔG of folding) for all combinations of these mutations. Next, we fit a quantitative model that predicts MIC as a function of kcat/KM and ΔG. While kcat/KM explains ∼54% of the variance in cefotaxime MIC (∼92% after log transformation), ΔG does not improve explanatory power of the model. We also find that cefotaxime MIC rises more slowly in kcat/KM than predicted. Several explanations for these discrepancies are suggested. Finally, we demonstrate substantial sign epistasis in MIC and kcat/KM, and antagonistic pleiotropy between phenotypes, in spite of near numerical additivity in the system. Thus constraints on selectively accessible trajectories, as well as limitations in our ability to explain such constraints in terms of underlying mechanisms are observed in a comparatively "well-behaved" system.

Detection of oligonucleotides by PNA-peptide conjugates recognizing the biarsenical fluorescein complex FlAsH-EDT2

We report the application of the arsenical complex FlAsH-EDT₂ for the identification of oligonucleotide sequences. We designed PNA sequences conjugated to either a tetracysteine motif and to split tetracysteine sequences, that are recognized by FlAsH. The effect of conjugation of the PNA to the tetracysteine peptide and RNA hybridization on the fluorescence of the arsenical complex has been investigated. The reconstitution of the tetracysteine motif, starting from 15-mer PNAs conjugated to split tetracysteine sequences and hybridized to a complementary oligonucleotide was also explored.

Crowding in Cellular Environments at an Atomistic Level from Computer Simulations

The effects of crowding in biological environments on biomolecular structure, dynamics, and function remain not well understood. Computer simulations of atomistic models of concentrated peptide and protein systems at different levels of complexity are beginning to provide new insights. Crowding, weak interactions with other macromolecules and metabolites, and altered solvent properties within cellular environments appear to remodel the energy landscape of peptides and proteins in significant ways including the possibility of native state destabilization. Crowding is also seen to affect dynamic properties, both conformational dynamics and diffusional properties of macromolecules. Recent simulations that address these questions are reviewed here and discussed in the context of relevant experiments.

Physicochemical code for quinary protein interactions in Escherichia coli

How proteins sense and navigate the cellular interior to find their functional partners remains poorly understood. An intriguing aspect of this search is that it relies on diffusive encounters with the crowded cellular background, made up of protein surfaces that are largely nonconserved. The question is then if/how this protein search is amenable to selection and biological control. To shed light on this issue, we examined the motions of three evolutionary divergent proteins in the Escherichia coli cytoplasm by in-cell NMR. The results show that the diffusive in-cell motions, after all, follow simplistic physical-chemical rules: The proteins reveal a common dependence on (i) net charge density, (ii) surface hydrophobicity, and (iii) the electric dipole moment. The bacterial protein is here biased to move relatively freely in the bacterial interior, whereas the human counterparts more easily stick. Even so, the in-cell motions respond predictably to surface mutation, allowing us to tune and intermix the protein's behavior at will. The findings show how evolution can swiftly optimize the diffuse background of protein encounter complexes by just single-point mutations, and provide a rational framework for adjusting the cytoplasmic motions of individual proteins, e.g., for rescuing poor in-cell NMR signals and for optimizing protein therapeutics.

Osmotic Shock Induced Protein Destabilization in Living Cells and Its Reversal by Glycine Betaine

Many organisms can adapt to changes in the solute content of their surroundings (i.e., the osmolarity). Hyperosmotic shock causes water efflux and a concomitant reduction in cell volume, which is countered by the accumulation of osmolytes. This volume reduction increases the crowded nature of the cytoplasm, which is expected to affect protein stability. In contrast to traditional theory, which predicts that more crowded conditions can only increase protein stability, recent work shows that crowding can destabilize proteins through transient attractive interactions. Here, we quantify protein stability in living Escherichia coli cells before and after hyperosmotic shock in the presence and absence of the osmolyte, glycine betaine. The 7-kDa N-terminal src-homology 3 domain of Drosophila signal transduction protein drk is used as the test protein. We find that hyperosmotic shock decreases SH3 stability in cells, consistent with the idea that transient attractive interactions are important under physiologically relevant crowded conditions. The subsequent uptake of glycine betaine returns SH3 to the stability observed without osmotic shock. These results highlight the effect of transient attractive interactions on protein stability in cells and provide a new explanation for why stressed cells accumulate osmolytes.

Towards understanding cellular structure biology: In-cell NMR

To watch biological macromolecules perform their functions inside the living cells is the dream of any biologists. In-cell nuclear magnetic resonance is a branch of biomolecular NMR spectroscopy that can be used to observe the structures, interactions and dynamics of these molecules in the living cells at atomic level. In principle, in-cell NMR can be applied to different cellular systems to achieve biologically relevant structural and functional information. In this review, we summarize the existing approaches in this field and discuss its applications in protein interactions, folding, stability and post-translational modifications. We hope this review will emphasize the effectiveness of in-cell NMR for studies of intricate biological processes and for structural analysis in cellular environments.

Incorporation of Hard and Soft Protein-Protein Interactions into Models for Crowding Effects in Binary and Ternary Protein Mixtures. Comparison of Approximate Analytical Solutions with Numerical Simulation

In order to better understand how nonspecific interactions between solutes can modulate specific biochemical reactions taking place in complex media, we introduce a simplified model aimed at elucidating general principles. In this model, solutions containing two or three species of interacting globular proteins are modeled as a fluid of spherical particles interacting through square well potentials that qualitatively capture both steric hard core repulsion and longer-ranged attraction or repulsion. The excess chemical potential, or free energy of solvation, of each particle species is calculated as a function of species concentrations, particle radii, and square well interaction range and depth. The results of analytical models incorporating two-body and three-body interactions are compared with the estimates of free energy obtained via Widom insertion into simulated equilibrium square-well fluids. The analytical models agree well with results of numeric simulations carried out for a variety of model parameters and fluid compositions up to a total particle volume fraction of ca. 0.2.

Macromolecular Crowding In Vitro, In Vivo, and In Between

Biochemical processes take place in heterogeneous and highly volume-occupied or crowded environments that can considerably influence the reactivity and distribution of participating macromolecules. We summarize here the thermodynamic consequences of excluded-volume and long-range nonspecific intermolecular interactions for macromolecular reactions in volume-occupied media. In addition, we summarize and compare the information content of studies of crowding in vitro and in vivo. We emphasize the importance of characterizing the behavior not only of labeled tracer macromolecules but also the composition and behavior of unlabeled macromolecules in the immediate vicinity of the tracer. Finally, we propose strategies for extending quantitative analyses of crowding in simple model systems to increasingly complex media up to and including intact cells.

The macromolecular crowding effect--from in vitro into the cell

The influence of the cellular milieu, a complex and crowded solvent, is often neglected when biomolecular structure and function are studied in vitro. To mimic the cellular environment, crowding effects are commonly induced in vitro using artificial crowding agents like Ficoll or dextran. However, it is unclear if such effects are also observed in cellulo. Diverging results on protein stability in living cells point out the need for new quantitative methods to investigate the contributions of excluded volume and nonspecific interactions to the cellular crowding effect. We show how new crowding sensitive probes may be utilized to directly investigate crowding effects in living cells. Moreover, we discuss processes where crowding effects could play a crucial role in molecular cell biology.

Subcellular modulation of protein VlsE stability and folding kinetics

The interior of a cell interacts differently with proteins than a dilute buffer because of a wide variety of macromolecules, chaperones, and osmolytes that crowd and interact with polypeptide chains. We compare folding of fluorescent constructs of protein VlsE among three environments inside cells. The nucleus increases the stability of VlsE relative to the cytoplasm, but slows down folding kinetics. VlsE is also more stable in the endoplasmic reticulum, but unlike PGK, tends to aggregate there. Although fluorescent-tagged VlsE and PGK show opposite stability trends from in vitro to the cytoplasm, their trends from cytoplasm to nucleus are similar.

Peptide folding in the presence of interacting protein crowders

Using Monte Carlo methods, we explore and compare the effects of two protein crowders, BPTI and GB1, on the folding thermodynamics of two peptides, the compact helical trp-cage and the β-hairpin-forming GB1m3. The thermally highly stable crowder proteins are modeled using a fixed backbone and rotatable side-chains, whereas the peptides are free to fold and unfold. In the simulations, the crowder proteins tend to distort the trp-cage fold, while having a stabilizing effect on GB1m3. The extent of the effects on a given peptide depends on the crowder type. Due to a sticky patch on its surface, BPTI causes larger changes than GB1 in the melting properties of the peptides. The observed effects on the peptides stem largely from attractive and specific interactions with the crowder surfaces, and differ from those seen in reference simulations with purely steric crowder particles.

The Effect of Fluorescent Protein Tags on Phosphoglycerate Kinase Stability Is Nonadditive

It is frequently assumed that fluorescent protein tags used in biological imaging experiments are minimally perturbing to their host protein. As in-cell experiments become more quantitative and measure rates and equilibrium constants, rather than just "on-off" activity or the presence of a protein, it becomes more important to understand such perturbations. One criterion for a protein modification to be a perturbation is additivity of two perturbations (a linear effect on the protein free energy). Here we show that adding fluorescent protein tags to a host protein in vitro has a large nonadditive effect on its folding free energy. We compare an unlabeled, three singly labeled, and a doubly labeled enzyme (phosphoglycerate kinase). We propose two mechanisms for nonadditivity. In the "quinary interaction" mechanism, two tags interact transiently with one another, relieving the host protein from unfavorable tag-protein interactions. In the "crowding" mechanism, adding two tags provides the minimal crowding necessary to overcome destabilizing interactions of individual tags with the host protein. Both of these mechanisms affect protein stability in cells; we show here that they must also be considered for tagged proteins used for reference in vitro.

In-cell thermodynamics and a new role for protein surfaces

There is abundant, physiologically relevant knowledge about protein cores; they are hydrophobic, exquisitely well packed, and nearly all hydrogen bonds are satisfied. An equivalent understanding of protein surfaces has remained elusive because proteins are almost exclusively studied in vitro in simple aqueous solutions. Here, we establish the essential physiological roles played by protein surfaces by measuring the equilibrium thermodynamics and kinetics of protein folding in the complex environment of living Escherichia coli cells, and under physiologically relevant in vitro conditions. Fluorine NMR data on the 7-kDa globular N-terminal SH3 domain of Drosophila signal transduction protein drk (SH3) show that charge-charge interactions are fundamental to protein stability and folding kinetics in cells. Our results contradict predictions from accepted theories of macromolecular crowding and show that cosolutes commonly used to mimic the cellular interior do not yield physiologically relevant information. As such, we provide the foundation for a complete picture of protein chemistry in cells.

Associative Interactions in Crowded Solutions of Biopolymers Counteract Depletion Effects

The cytosol of Escherichia coli is an extremely crowded environment, containing high concentrations of biopolymers which occupy 20-30% of the available volume. Such conditions are expected to yield depletion forces, which strongly promote macromolecular complexation. However, crowded macromolecule solutions, like the cytosol, are very prone to nonspecific associative interactions that can potentially counteract depletion. It remains unclear how the cytosol balances these opposing interactions. We used a FRET-based probe to systematically study depletion in vitro in different crowded environments, including a cytosolic mimic, E. coli lysate. We also studied bundle formation of FtsZ protofilaments under identical crowded conditions as a probe for depletion interactions at much larger overlap volumes of the probe molecule. The FRET probe showed a more compact conformation in synthetic crowding agents, suggesting strong depletion interactions. However, depletion was completely negated in cell lysate and other protein crowding agents, where the FRET probe even occupied slightly more volume. In contrast, bundle formation of FtsZ protofilaments proceeded as readily in E. coli lysate and other protein solutions as in synthetic crowding agents. Our experimental results and model suggest that, in crowded biopolymer solutions, associative interactions counterbalance depletion forces for small macromolecules. Furthermore, the net effects of macromolecular crowding will be dependent on both the size of the macromolecule and its associative interactions with the crowded background.

Thermodynamics of protein destabilization in live cells

Although protein folding and stability have been well explored under simplified conditions in vitro, it is yet unclear how these basic self-organization events are modulated by the crowded interior of live cells. To find out, we use here in-cell NMR to follow at atomic resolution the thermal unfolding of a β-barrel protein inside mammalian and bacterial cells. Challenging the view from in vitro crowding effects, we find that the cells destabilize the protein at 37 °C but with a conspicuous twist: While the melting temperature goes down the cold unfolding moves into the physiological regime, coupled to an augmented heat-capacity change. The effect seems induced by transient, sequence-specific, interactions with the cellular components, acting preferentially on the unfolded ensemble. This points to a model where the in vivo influence on protein behavior is case specific, determined by the individual protein's interplay with the functionally optimized "interaction landscape" of the cellular interior.

Quantitative Description of a Protein Fitness Landscape Based on Molecular Features

Understanding the driving forces behind protein evolution requires the ability to correlate the molecular impact of mutations with organismal fitness. To address this issue, we employ here metallo-β-lactamases as a model system, which are Zn(II) dependent enzymes that mediate antibiotic resistance. We present a study of all the possible evolutionary pathways leading to a metallo-β-lactamase variant optimized by directed evolution. By studying the activity, stability and Zn(II) binding capabilities of all mutants in the preferred evolutionary pathways, we show that this local fitness landscape is strongly conditioned by epistatic interactions arising from the pleiotropic effect of mutations in the different molecular features of the enzyme. Activity and stability assays in purified enzymes do not provide explanatory power. Instead, measurement of these molecular features in an environment resembling the native one provides an accurate description of the observed antibiotic resistance profile. We report that optimization of Zn(II) binding abilities of metallo-β-lactamases during evolution is more critical than stabilization of the protein to enhance fitness. A global analysis of these parameters allows us to connect genotype with fitness based on quantitative biochemical and biophysical parameters.

Quinary structure modulates protein stability in cells

Protein quinary interactions organize the cellular interior and its metabolism. Although the interactions stabilizing secondary, tertiary, and quaternary protein structure are well defined, details about the protein-matrix contacts that comprise quinary structure remain elusive. This gap exists because proteins function in the crowded cellular environment, but are traditionally studied in simple buffered solutions. We use NMR-detected H/D exchange to quantify quinary interactions between the B1 domain of protein G and the cytosol of Escherichia coli. We demonstrate that a surface mutation in this protein is 10-fold more destabilizing in cells than in buffer, a surprising result that firmly establishes the significance of quinary interactions. Remarkably, the energy involved in these interactions can be as large as the energies that stabilize specific protein complexes. These results will drive the critical task of implementing quinary structure into models for understanding the proteome.

Probing the kinetic stabilities of Friedreich's ataxia clinical variants using a solid phase GroEL chaperonin capture platform

Numerous human diseases are caused by protein folding defects where the protein may become more susceptible to degradation or aggregation. Aberrant protein folding can affect the kinetic stability of the proteins even if these proteins appear to be soluble in vivo. Experimental discrimination between functional properly folded and misfolded nonfunctional conformers is not always straightforward at near physiological conditions. The differences in the kinetic behavior of two initially folded frataxin clinical variants were examined using a high affinity chaperonin kinetic trap approach at 25 °C. The kinetically stable wild type frataxin (FXN) shows no visible partitioning onto the chaperonin. In contrast, the clinical variants FXN-p.Asp122Tyr and FXN-p.Ile154Phe kinetically populate partial folded forms that tightly bind the GroEL chaperonin platform. The initially soluble FXN-p.Ile154Phe variant partitions onto GroEL more rapidly and is more kinetically liable. These differences in kinetic stability were confirmed using differential scanning fluorimetry. The kinetic and aggregation stability differences of these variants may lead to the distinct functional impairments described in Friedreich's ataxia, the neurodegenerative disease associated to frataxin functional deficiency. This chaperonin platform approach may be useful for identifying small molecule stabilizers since stabilizing ligands to frataxin variants should lead to a concomitant decrease in chaperonin binding.

Cytosolic selection systems to study protein stability

Here we describe biosensors that provide readouts for protein stability in the cytosolic compartment of prokaryotes. These biosensors consist of tripartite sandwich fusions that link the in vitro stability or aggregation susceptibility of guest proteins to the in vivo resistance of host cells to the antibiotics kanamycin, spectinomycin, and nourseothricin. These selectable markers confer antibiotic resistance in a wide range of hosts and are easily quantifiable. We show that mutations within guest proteins that affect their stability alter the antibiotic resistances of the cells expressing the biosensors in a manner that is related to the in vitro stabilities of the mutant guest proteins. In addition, we find that polyglutamine tracts of increasing length are associated with an increased tendency to form amyloids in vivo and, in our sandwich fusion system, with decreased resistance to aminoglycoside antibiotics. We demonstrate that our approach allows the in vivo analysis of protein stability in the cytosolic compartment without the need for prior structural and functional knowledge.

Residue level quantification of protein stability in living cells

The intracellular milieu differs from the dilute conditions in which most biophysical and biochemical studies are performed. This difference has led both experimentalists and theoreticians to tackle the challenging task of understanding how the intracellular environment affects the properties of biopolymers. Despite a growing number of in-cell studies, there is a lack of quantitative, residue-level information about equilibrium thermodynamic protein stability under nonperturbing conditions. We report the use of NMR-detected hydrogen-deuterium exchange of quenched cell lysates to measure individual opening free energies of the 56-aa B1 domain of protein G (GB1) in living Escherichia coli cells without adding destabilizing cosolutes or heat. Comparisons to dilute solution data (pH 7.6 and 37 °C) show that opening free energies increase by as much as 1.14 ± 0.05 kcal/mol in cells. Importantly, we also show that homogeneous protein crowders destabilize GB1, highlighting the challenge of recreating the cellular interior. We discuss our findings in terms of hard-core excluded volume effects, charge-charge GB1-crowder interactions, and other factors. The quenched lysate method identifies the residues most important for folding GB1 in cells, and should prove useful for quantifying the stability of other globular proteins in cells to gain a more complete understanding of the effects of the intracellular environment on protein chemistry.

An osmolyte mitigates the destabilizing effect of protein crowding

Most theories predict that macromolecular crowding stabilizes globular proteins, but recent studies show that weak attractive interactions can result in crowding-induced destabilization. Osmolytes are ubiquitous in biology and help protect cells against stress. Given that dehydration stress adds to the crowded nature of the cytoplasm, we speculated that cells might use osmolytes to overcome the destabilization caused by the increased weak interactions that accompany desiccation. We used NMR-detected amide proton exchange experiments to measure the stability of the test protein chymotrypsin inhibitor 2 under physiologically relevant crowded conditions in the presence and absence of the osmolyte glycine betaine. The osmolyte overcame the destabilizing effect of the cytosol. This result provides a physiologically relevant explanation for the accumulation of osmolytes by dehydration-stressed cells.