Modulation of calmodulin plasticity by the effect of macromolecular crowding

Contrasting effect of ficoll on apo and holo forms of bacterial chemotaxis protein Y: Selective destabilization of the conformationally altered holo form

Chemotaxis Y (CheY), upon metal binding, displays a drastic alteration in its structure and stability. This premise prompted us to study the effect of crowding on the two conformationally distinct states of the same test protein. A comparative analysis on the structure and thermal stability in the presence and absence of the macromolecular crowder, ficoll, and its monomeric unit, sucrose, revealed a contrasting effect of ficoll on the apo and holo forms. In the presence of ficoll while the thermal stability (T_m) of the apo form is enhanced, the thermal stability of the holo form is reduced. The selective lowering of T_m for the holo form in the combined presence of ficoll and sucrose and not in sucrose alone suggests that the contrasting effect is due to the macromolecular nature of ficoll. Since metal-protein interaction remains unperturbed in the presence of ficoll and Mg²⁺ sequestration is ruled out in a systematic manner the alternative possibility for the exclusive reduction in the thermal stability of the holo form is the ficoll-induced modulation of the relative population of apo and holo forms of CheY.

Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca2+-Calmodulin

Calmodulin (CaM) is a calcium-binding protein that transduces signals to downstream proteins through target binding upon calcium binding in a time-dependent manner. Understanding the target binding process that tunes CaM’s affinity for the calcium ions (Ca²⁺), or vice versa, may provide insight into how Ca²⁺-CaM selects its target binding proteins. However, modeling of Ca²⁺-CaM in molecular simulations is challenging because of the gross structural changes in its central linker regions while the two lobes are relatively rigid due to tight binding of the Ca²⁺ to the calcium-binding loops where the loop forms a pentagonal bipyramidal coordination geometry with Ca²⁺. This feature that underlies the reciprocal relation between Ca²⁺ binding and target binding of CaM, however, has yet to be considered in the structural modeling. Here, we presented a coarse-grained model based on the Associative memory, Water mediated, Structure, and Energy Model (AWSEM) protein force field, to investigate the salient features of CaM. Particularly, we optimized the force field of CaM and that of Ca²⁺ ions by using its coordination chemistry in the calcium-binding loops to match with experimental observations. We presented a “community model” of CaM that is capable of sampling various conformations of CaM, incorporating various calcium-binding states, and carrying the memory of binding with various targets, which sets the foundation of the reciprocal relation of target binding and Ca²⁺ binding in future studies.

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.

NMR-Based Methods for Protein Analysis

Nuclear magnetic resonance (NMR) spectroscopy is a well-established method for analyzing protein structure, interaction, and dynamics at atomic resolution and in various sample states including solution state, solid state, and membranous environment. Thanks to rapid NMR methodology development, the past decade has witnessed a growing number of protein NMR studies in complex systems ranging from membrane mimetics to living cells, which pushes the research frontier further toward physiological environments and offers unique insights in elucidating protein functional mechanisms. In particular, in-cell NMR has become a method of choice for bridging the huge gap between structural biology and cell biology. Herein, we review the recent developments and applications of NMR methods for protein analysis in close-to-physiological environments, with special emphasis on in-cell protein structural determination and the analysis of protein dynamics, both difficult to be accessed by traditional methods.

Folding of poly-amino acids and intrinsically disordered proteins in overcrowded milieu induced by pH change

pH-induced structural changes of the synthetic homopolypeptides poly-E, poly-K, poly-R, and intrinsically disordered proteins (IDPs) prothymosin α (ProTα) and linker histone H1, in concentrated PEG solutions simulating macromolecular crowding conditions within the membrane-less organelles, were characterized. The conformational transitions of the studied poly-amino acids in the concentrated PEG solutions depend on the polymerization degree of these homopolypeptides, the size of their side chains, the charge distribution of the side chains, and the crowding agent concentration. The results obtained for poly-amino acids are valid for IDPs having a significant total charge. The overcrowded conditions promote a significant increase in the cooperativity of the pH-induced coil-α-helix transition of ProTα and provoke histone H1 aggregation. The most favorable conditions for the pH-induced structural transitions in concentrated PEG solutions are realized when the charged residues are grouped in blocks, and when the distance between the end of the side group carrying charge and the backbone is small. Therefore, the block-wise distribution of charged residues within the IDPs not only plays an important role in the liquid-liquid phase transitions, but may also define the expressivity of structural transitions of these proteins in the overcrowded conditions of the membrane-less organelles.

A calmodulin-like protein (LCALA) is a new Leishmania amazonensis candidate for telomere end-binding protein

BACKGROUND

Leishmania spp. telomeres are composed of 5'-TTAGGG-3' repeats associated with proteins. We have previously identified LaRbp38 and LaRPA-1 as proteins that bind the G-rich telomeric strand. At that time, we had also partially characterized a protein: DNA complex, named LaGT1, but we could not identify its protein component.

METHODS AND RESULTS

Using protein-DNA interaction and competition assays, we confirmed that LaGT1 is highly specific to the G-rich telomeric single-stranded DNA. Three protein bands, with LaGT1 activity, were isolated from affinity-purified protein extracts in-gel digested, and sequenced de novo using mass spectrometry analysis. In silico analysis of the digested peptide identified them as a putative calmodulin with sequences identical to the T. cruzi calmodulin. In the Leishmania genome, the calmodulin ortholog is present in three identical copies. We cloned and sequenced one of the gene copies, named it LCalA, and obtained the recombinant protein. Multiple sequence alignment and molecular modeling showed that LCalA shares homology to most eukaryotes calmodulin. In addition, we demonstrated that LCalA is nuclear, partially co-localizes with telomeres and binds in vivo the G-rich telomeric strand. Recombinant LCalA can bind specifically and with relative affinity to the G-rich telomeric single-strand and to a 3'G-overhang, and DNA binding is calcium dependent.

CONCLUSIONS

We have described a novel candidate component of Leishmania telomeres, LCalA, a nuclear calmodulin that binds the G-rich telomeric strand with high specificity and relative affinity, in a calcium-dependent manner.

GENERAL SIGNIFICANCE

LCalA is the first reported calmodulin that binds in vivo telomeric DNA.

Calcium triggers reversal of calmodulin on nested anti-parallel sites in the IQ motif of the neuronal voltage-dependent sodium channel NaV1.2

Several members of the voltage-gated sodium channel family are regulated by calmodulin (CaM) and ionic calcium. The neuronal voltage-gated sodium channel NaV1.2 contains binding sites for both apo (calcium-depleted) and calcium-saturated CaM. We have determined equilibrium dissociation constants for rat NaV1.2 IQ motif [IQRAYRRYLLK] binding to apo CaM (~3nM) and (Ca2+)4-CaM (~85nM), showing that apo CaM binding is favored by 30-fold. For both apo and (Ca2+)4-CaM, NMR demonstrated that NaV1.2 IQ motif peptide (NaV1.2IQp) exclusively made contacts with C-domain residues of CaM (CaMC). To understand how calcium triggers conformational change at the CaM-IQ interface, we determined a solution structure (2M5E.pdb) of (Ca2+)2-CaMC bound to NaV1.2IQp. The polarity of (Ca2+)2-CaMC relative to the IQ motif was opposite to that seen in apo CaMC-Nav1.2IQp (2KXW), revealing that CaMC recognizes nested, anti-parallel sites in Nav1.2IQp. Reversal of CaM may require transient release from the IQ motif during calcium binding, and facilitate a re-orientation of CaMN allowing interactions with non-IQ NaV1.2 residues or auxiliary regulatory proteins interacting in the vicinity of the IQ motif.

Tracking the Energy Flow on Nanoscale via Sample-Transmitted Excitation Photoluminescence Spectroscopy

Tracking the energy flow in nanoscale materials is an important yet challenging goal. Experimental methods for probing the intermolecular energy transfer (ET) are often burdened by the spectral crosstalk between donor and acceptor species, which complicates unraveling their individual contributions. This issue is particularly prominent in inorganic nanoparticles and biological macromolecules featuring broad absorbing profiles. Here, we demonstrate a general spectroscopic strategy for measuring the ET efficiency between nanostructured or molecular dyes exhibiting a significant donor-acceptor spectral overlap. The reported approach is enabled through spectral shaping of the broadband excitation light with solutions of donor molecules, which inhibits the excitation of respective donor species in the sample. The resulting changes in the acceptor emission induced by the spectral modulation of the excitation beam are then used to determine the quantum efficiency and the rate of ET processes between arbitrary fluorophores (molecules, nanoparticles, polymers) with high accuracy. The feasibility of the reported method was demonstrated using a control donor-acceptor system utilizing a protein-bridged Cy3-Cy5 dye pair and subsequently applied for studying the energy flow in a CdSe₅₆₀-CdSe₆₀₀ binary nanocrystal film.

Structure and Conformational Properties of d-Glucose/d-Galactose-Binding Protein in Crowded Milieu

Conformational changes of d-glucose/d-galactose-binding protein (GGBP) were studied under molecular crowding conditions modeled by concentrated solutions of polyethylene glycols (PEG-12000, PEG-4000, and PEG-600), Ficoll-70, and Dextran-70, addition of which induced noticeable structural changes in the GGBP molecule. All PEGs promoted compaction of GGBP and lead to the increase in ordering of its structure. Concentrated solutions of PEG-12000 and PEG-4000 caused GGBP aggregation. Although Ficoll-70 and Dextran-70 also promoted increase in the GGBP ordering, the structural outputs were different for different crowders. For example, in comparison with the GGBP in buffer, the intrinsic fluorescence spectrum of this protein was shifted to short-wave region in the presence of PEGs but was red-shifted in the presence of Ficoll-70 and Dextran-70. It was hypothesized that this difference could be due to the specific interaction of GGBP with the sugar-based polymers (Ficoll-70 and Dextran-70), indicating that protein can adopt different conformations in solutions containing molecular crowders of different chemical nature. It was also shown that all tested crowding agents were able to stabilize GGBP structure shifting the GGBP guanidine hydrochloride (GdnHCl)-induced unfolding curves to higher denaturant concentrations, but their stabilization capabilities did not depend on the hydrodynamic dimensions of the polymers molecules. Refolding of GGBP was complicated by protein aggregation in all tested solutions of crowding agents. The lowest yield of refolded protein was achieved in the highly concentrated solutions of PEG-12000. These data support the previous notion that the influence of macromolecular crowders on proteins is rather complex phenomenon that extends beyond the excluded volume effects.

Relative Cosolute Size Influences the Kinetics of Protein-Protein Interactions

Protein signaling occurs in crowded intracellular environments, and while high concentrations of macromolecules are postulated to modulate protein-protein interactions, analysis of their impact at each step of the reaction pathway has not been systematically addressed. Potential cosolute-induced alterations in target association are particularly important for a signaling molecule like calmodulin (CaM), where competition among >300 targets governs which pathways are selectively activated. To explore how high concentrations of cosolutes influence CaM-target affinity and kinetics, we methodically investigated each step of the CaM-target binding mechanism under crowded or osmolyte-rich environments mimicked by ficoll-70, dextran-10, and sucrose. All cosolutes stabilized compact conformers of CaM and modulated association kinetics by affecting diffusion and rates of conformational change; however, the results showed that differently sized molecules had variable effects to enhance or impede unique steps of the association pathway. On- and off-rates were modulated by all cosolutes in a compensatory fashion, producing little change in steady-state affinity. From this work insights were gained on how high concentrations of inert crowding agents and osmolytes fit into a kinetic framework to describe protein-protein interactions relevant for cellular signaling.

Conformational frustration in calmodulin-target recognition

Calmodulin (CaM) is a primary calcium (Ca(2+) )-signaling protein that specifically recognizes and activates highly diverse target proteins. We explored the molecular basis of target recognition of CaM with peptides representing the CaM-binding domains from two Ca(2+) -CaM-dependent kinases, CaMKI and CaMKII, by employing experimentally constrained molecular simulations. Detailed binding route analysis revealed that the two CaM target peptides, although similar in length and net charge, follow distinct routes that lead to a higher binding frustration in the CaM-CaMKII complex than in the CaM-CaMKI complex. We discovered that the molecular origin of the binding frustration is caused by intermolecular contacts formed with the C-domain of CaM that need to be broken before the formation of intermolecular contacts with the N-domain of CaM. We argue that the binding frustration is important for determining the kinetics of the recognition process of proteins involving large structural fluctuations.

Beyond the excluded volume effects: mechanistic complexity of the crowded milieu

Macromolecular crowding is known to affect protein folding, binding of small molecules, interaction with nucleic acids, enzymatic activity, protein-protein interactions, and protein aggregation. Although for a long time it was believed that the major mechanism of the action of crowded environments on structure, folding, thermodynamics, and function of a protein can be described in terms of the excluded volume effects, it is getting clear now that other factors originating from the presence of high concentrations of “inert” macromolecules in crowded solution should definitely be taken into account to draw a more complete picture of a protein in a crowded milieu. This review shows that in addition to the excluded volume effects important players of the crowded environments are viscosity, perturbed diffusion, direct physical interactions between the crowding agents and proteins, soft interactions, and, most importantly, the effects of crowders on solvent properties.

Unusual domain movement in a multidomain protein in the presence of macromolecular crowders

FRET-based distance changes between domains I and II of HSA reveal significant differences in the manner in which dextran based macromolecular crowding agents influence the serum protein.

What macromolecular crowding can do to a protein

The intracellular environment represents an extremely crowded milieu, with a limited amount of free water and an almost complete lack of unoccupied space. Obviously, slightly salted aqueous solutions containing low concentrations of a biomolecule of interest are too simplistic to mimic the “real life” situation, where the biomolecule of interest scrambles and wades through the tightly packed crowd. In laboratory practice, such macromolecular crowding is typically mimicked by concentrated solutions of various polymers that serve as model “crowding agents”. Studies under these conditions revealed that macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, protein-protein interactions, protein-nucleic acid interactions, and pathological aggregation. The goal of this review is to systematically analyze currently available experimental data on the variety of effects of macromolecular crowding on a protein molecule. The review covers more than 320 papers and therefore represents one of the most comprehensive compendia of the current knowledge in this exciting area.

Minimal effects of macromolecular crowding on an intrinsically disordered protein: a small-angle neutron scattering study

Small-angle neutron scattering was used to study the effects of macromolecular crowding by two globular proteins, i.e., bovine pancreatic trypsin inhibitor and equine metmyoglobin, on the conformational ensemble of an intrinsically disordered protein, the N protein of bacteriophage λ. The λ N protein was uniformly labeled with (2)H, and the concentrations of D2O in the samples were adjusted to match the neutron scattering contrast of the unlabeled crowding proteins, thereby masking their contribution to the scattering profiles. Scattering from the deuterated λ N was recorded for samples containing up to 0.12 g/mL bovine pancreatic trypsin inhibitor or 0.2 g/mL metmyoglobin. The radius of gyration of the uncrowded protein was estimated to be 30 Å and was found to be remarkably insensitive to the presence of crowders, varying by <2 Å for the highest crowder concentrations. The scattering profiles were also used to estimate the fractal dimension of λ N, which was found to be ∼1.8 in the absence or presence of crowders, indicative of a well-solvated and expanded random coil under all of the conditions examined. These results are contrary to the predictions of theoretical treatments and previous experimental studies demonstrating compaction of unfolded proteins by crowding with polymers such as dextran and Ficoll. A computational simulation suggests that some previous treatments may have overestimated the effective volumes of disordered proteins and the variation of these volumes within an ensemble. The apparent insensitivity of λ N to crowding may also be due in part to weak attractive interactions with the crowding proteins, which may compensate for the effects of steric exclusion.

The effect of macromolecular crowding on the electrostatic component of barnase-barstar binding: a computational, implicit solvent-based study

Macromolecular crowding within the cell can impact both protein folding and binding. Earlier models of cellular crowding focused on the excluded volume, entropic effect of crowding agents, which generally favors compact protein states. Recently, other effects of crowding have been explored, including enthalpically-related crowder–protein interactions and changes in solvation properties. In this work, we explore the effects of macromolecular crowding on the electrostatic desolvation and solvent-screened interaction components of protein–protein binding. Our simple model enables us to focus exclusively on the electrostatic effects of water depletion on protein binding due to crowding, providing us with the ability to systematically analyze and quantify these potentially intuitive effects. We use the barnase–barstar complex as a model system and randomly placed, uncharged spheres within implicit solvent to model crowding in an aqueous environment. On average, we find that the desolvation free energy penalties incurred by partners upon binding are lowered in a crowded environment and solvent-screened interactions are amplified. At a constant crowder density (fraction of total available volume occupied by crowders), this effect generally increases as the radius of model crowders decreases, but the strength and nature of this trend can depend on the water probe radius used to generate the molecular surface in the continuum model. In general, there is huge variation in desolvation penalties as a function of the random crowder positions. Results with explicit model crowders can be qualitatively similar to those using a lowered “effective” solvent dielectric to account for crowding, although the “best” effective dielectric constant will likely depend on multiple system properties. Taken together, this work systematically demonstrates, quantifies, and analyzes qualitative intuition-based insights into the effects of water depletion due to crowding on the electrostatic component of protein binding, and it provides an initial framework for future analyses.

Protein conformational populations and functionally relevant substates

Functioning proteins do not remain fixed in a unique structure, but instead they sample a range of conformations facilitated by motions within the protein. Even in the native state, a protein exists as a collection of interconverting conformations driven by thermodynamic fluctuations. Motions on the fast time scale allow a protein to sample conformations in the nearby area of its conformational landscape, while motions on slower time scales give it access to conformations in distal areas of the landscape. Emerging evidence indicates that protein landscapes contain conformational substates with dynamic and structural features that support the designated function of the protein. Nuclear magnetic resonance (NMR) experiments provide information about conformational ensembles of proteins. X-ray crystallography allows researchers to identify the most populated states along the landscape, and computational simulations give atom-level information about the conformational substates of different proteins. This ability to characterize and obtain quantitative information about the conformational substates and the populations of proteins within them is allowing researchers to better understand the relationship between protein structure and dynamics and the mechanisms of protein function. In this Account, we discuss recent developments and challenges in the characterization of functionally relevant conformational populations and substates of proteins. In some enzymes, the sampling of functionally relevant conformational substates is connected to promoting the overall mechanism of catalysis. For example, the conformational landscape of the enzyme dihydrofolate reductase has multiple substates, which facilitate the binding and the release of the cofactor and substrate and catalyze the hydride transfer. For the enzyme cyclophilin A, computational simulations reveal that the long time scale conformational fluctuations enable the enzyme to access conformational substates that allow it to attain the transition state, therefore promoting the reaction mechanism. In the long term, this emerging view of proteins with conformational substates has broad implications for improving our understanding of enzymes, enzyme engineering, and better drug design. Researchers have already used photoactivation to modulate protein conformations as a strategy to develop a hypercatalytic enzyme. In addition, the alteration of the conformational substates through binding of ligands at locations other than the active site provides the basis for the design of new medicines through allosteric modulation.

Protein recognition and selection through conformational and mutually induced fit

Protein-protein interactions drive most every biological process, but in many instances the domains mediating recognition are disordered. How specificity in binding is attained in the absence of defined structure contrasts with well-established experimental and theoretical work describing ligand binding to protein. The signaling protein calmodulin presents a unique opportunity to investigate mechanisms for target recognition given that it interacts with several hundred different targets. By advancing coarse-grained computer simulations and experimental techniques, mechanistic insights were gained in defining the pathways leading to recognition and in how target selectivity can be achieved at the molecular level. A model requiring mutually induced conformational changes in both calmodulin and target proteins was necessary and broadly informs how proteins can achieve both high affinity and high specificity.

Reaching new levels of realism in modeling biological macromolecules in cellular environments

An increasing number of studies are aimed at modeling cellular environments in a comprehensive and realistic fashion. A major challenge in these efforts is how to bridge spatial and temporal scales over many orders of magnitude. Furthermore, there are additional challenges in integrating different aspects ranging from questions about biomolecular stability in crowded environments to the description of reactive processes on cellular scales. In this review, recent studies with models of biomolecules in cellular environments at different levels of detail are discussed in terms of their strengths and weaknesses. In particular, atomistic models, implicit representations of cellular environments, coarse-grained and spheroidal models of biomolecules, as well as the inclusion of reactive processes via reaction-diffusion models are described. Furthermore, strategies for integrating the different models into a comprehensive description of cellular environments are discussed.

Direct observation of protein unfolded state compaction in the presence of macromolecular crowding

Proteins fold and function in cellular environments that are crowded with other macromolecules. As a consequence of excluded volume effects, compact folded states of proteins should be indirectly stabilized due to destabilization of extended unfolded conformations. Here, we assess the role of excluded volume in terms of protein stability, structural dimensions and folding dynamics using a sugar-based crowding agent, dextran 20, and the small ribosomal protein S16 as a model system. To specifically address dimensions, we labeled the protein with BODIPY at two positions and measured Trp-BODIPY distances under different conditions. As expected, we found that dextran 20 (200 mg/ml) stabilized the variants against urea-induced unfolding. At conditions where the protein is unfolded, Förster resonance energy transfer measurements reveal that in the presence of dextran, the unfolded ensemble is more compact and there is residual structure left as probed by far-ultraviolet circular dichroism. In the presence of a crowding agent, folding rates are faster in the two-state regime, and at low denaturant concentrations, a kinetic intermediate is favored. Our study provides direct evidence for protein unfolded-state compaction in the presence of macromolecular crowding along with its energetic and kinetic consequences.

Where soft matter meets living matter--protein structure, stability, and folding in the cell

A protein is a biopolymer that self-assembles through the process of protein folding. A cell is a crowded space where the surrounding macromolecules of a protein can limit the number of ways of folding. These crowding macromolecules can also affect the shape and the size of a physically malleable, or 'soft, squishy', protein with regulatory purposes. In this review, we focus on the in silico approaches of coarse-grained molecular simulations that enable the investigation of protein folding in a cell-like environment. When these simulation results were compared with experimentally measured properties of a protein, such joint effort has yielded new ideas on the specific function of a protein in cells. We also highlighted the recent developments of computer modeling and simulations that encompass the importance of the shape of a macromolecule, the interactions between macromolecules, and the hydrodynamic interactions on the kinetics and thermodynamics of a protein in a high concentration of protein solution and in cytoplasmic environments.

Effects of macromolecular crowding agents on protein folding in vitro and in silico

Proteins fold and function inside cells which are environments very different from that of dilute buffer solutions most often used in traditional experiments. The crowded milieu results in excluded-volume effects, increased bulk viscosity and amplified chances for inter-molecular interactions. These environmental factors have not been accounted for in most mechanistic studies of protein folding executed during the last decades. The question thus arises as to how these effects-present when polypeptides normally fold in vivo-modulate protein biophysics. To address excluded volume effects, we use synthetic macromolecular crowding agents, which take up significant volume but do not interact with proteins, in combination with strategically selected proteins and a range of equilibrium and time-resolved biophysical (spectroscopic and computational) methods. In this review, we describe key observations on macromolecular crowding effects on protein stability, folding and structure drawn from combined in vitro and in silico studies. As expected based on Minton's early predictions, many proteins (apoflavodoxin, VlsE, cytochrome c, and S16) became more thermodynamically stable (magnitude depends inversely on protein stability in buffer) and, unexpectedly, for apoflavodoxin and VlsE, the folded states changed both secondary structure content and, for VlsE, overall shape in the presence of macromolecular crowding. For apoflavodoxin and cytochrome c, which have complex kinetic folding mechanisms, excluded volume effects made the folding energy landscapes smoother (i.e., less misfolding and/or kinetic heterogeneity) than in buffer.

Is buffer a good proxy for a crowded cell-like environment? A comparative NMR study of calmodulin side-chain dynamics in buffer and E. coli lysate

Biophysical studies of protein structure and dynamics are typically performed in a highly controlled manner involving only the protein(s) of interest. Comparatively fewer such studies have been carried out in the context of a cellular environment that typically involves many biomolecules, ions and metabolites. Recently, solution NMR spectroscopy, focusing primarily on backbone amide groups as reporters, has emerged as a powerful technique for investigating protein structure and dynamics in vivo and in crowded “cell-like” environments. Here we extend these studies through a comparative analysis of Ile, Leu, Val and Met methyl side-chain motions in apo, Ca²⁺-bound and Ca²⁺, peptide-bound calmodulin dissolved in aqueous buffer or in E. coli lysate. Deuterium spin relaxation experiments, sensitive to pico- to nano-second time-scale processes and Carr-Purcell-Meiboom-Gill relaxation dispersion experiments, reporting on millisecond dynamics, have been recorded. Both similarities and differences in motional properties are noted for calmodulin dissolved in buffer or in lysate. These results emphasize that while significant insights can be obtained through detailed “test-tube” studies, experiments performed under conditions that are “cell-like” are critical for obtaining a comprehensive understanding of protein motion in vivo and therefore for elucidating the relation between motion and function.

A physics-based approach of coarse-graining the cytoplasm of Escherichia coli (CGCYTO)

We have investigated protein stability in an environment of Escherichia coli cytoplasm using coarse-grained computer simulations. To coarse-grain a small slide of E. coli's cytoplasm consisting of over 16 million atoms, we have developed a self-assembled clustering algorithm (CGCYTO). CGCYTO uses the shape parameter and asphericity as well as a parameter λ (ranging from 0 to 1) that measures the covolume of a test protein and a macromolecule against the covolume of a test protein and a sphere of equal volume as that of a macromolecule for the criteria of coarse-graining a cytoplasmic model. A cutoff λ(c) = 0.8 was chosen based on the size of a test protein and computational resources and it determined the resolution of a coarse-grained cytoplasm. We compared the results from a polydisperse cytoplasmic model (PD model) produced by CGCYTO with two other coarse-grained hard-sphere cytoplasmic models: 1), F70 model, macromolecules in the cytoplasm were modeled by homogeneous hard spheres with a radius of 55 Å, the size of Ficoll70 and 2), HS model, each macromolecule in the cytoplasm was modeled by a hard sphere of equal volume. It was found that the folding temperature T(f) of a test protein (apoazurin) in a PD model is ~5° greater than that in a F70 model. In addition, the deviation of T(f) in a PD model is twice as much as that in a HS model when an apoazurin is randomly placed at different voids formed by particle fluctuations in PD models.

Folding dynamics of Trp-cage in the presence of chemical interference and macromolecular crowding. I

Proteins fold and function in the crowded environment of the cell's interior. In the recent years it has been well established that the so-called "macromolecular crowding" effect enhances the folding stability of proteins by destabilizing their unfolded states for selected proteins. On the other hand, chemical and thermal denaturation is often used in experiments as a tool to destabilize a protein by populating the unfolded states when probing its folding landscape and thermodynamic properties. However, little is known about the complicated effects of these synergistic perturbations acting on the kinetic properties of proteins, particularly when large structural fluctuations, such as protein folding, have been involved. In this study, we have first investigated the folding mechanism of Trp-cage dependent on urea concentration by coarse-grained molecular simulations where the impact of urea is implemented into an energy function of the side chain and/or backbone interactions derived from the all-atomistic molecular dynamics simulations with urea through a Boltzmann inversion method. In urea solution, the folding rates of a model miniprotein Trp-cage decrease and the folded state slightly swells due to a lack of contact formation between side chains at the terminal regions. In addition, the equilibrium m-values of Trp-cage from the computer simulations are in agreement with experimental measurements. We have further investigated the combined effects of urea denaturation and macromolecular crowding on Trp-cage's folding mechanism where crowding agents are modeled as hard-spheres. The enhancement of folding rates of Trp-cage is most pronounced by macromolecular crowding effect when the extended conformations of Trp-cast dominate at high urea concentration. Our study makes quantitatively testable predictions on protein folding dynamics in a complex environment involving both chemical denaturation and macromolecular crowding effects.

The effect of macromolecular crowding, ionic strength and calcium binding on calmodulin dynamics

The flexibility in the structure of calmodulin (CaM) allows its binding to over 300 target proteins in the cell. To investigate the structure-function relationship of CaM, we combined methods of computer simulation and experiments based on circular dichroism (CD) to investigate the structural characteristics of CaM that influence its target recognition in crowded cell-like conditions. We developed a unique multiscale solution of charges computed from quantum chemistry, together with protein reconstruction, coarse-grained molecular simulations, and statistical physics, to represent the charge distribution in the transition from apoCaM to holoCaM upon calcium binding. Computationally, we found that increased levels of macromolecular crowding, in addition to calcium binding and ionic strength typical of that found inside cells, can impact the conformation, helicity and the EF hand orientation of CaM. Because EF hand orientation impacts the affinity of calcium binding and the specificity of CaM's target selection, our results may provide unique insight into understanding the promiscuous behavior of calmodulin in target selection inside cells.

Proteins are workhorses for driving biological functions inside cells. Calmodulin (CaM) is a protein that can carry cellular signals by triggered conformational changes due to calcium binding that alters target binding. Interestingly, CaM is able to bind over 300 targets. One of the challenges in characterizing CaM's ability to bind multiple targets lies in that CaM is a flexible protein and its structure is easily modulated by the physicochemical changes in its surroundings, particularly inside a complex cellular milieu. In order to determine structure-function relationships of CaM, we employed a combined approach of experiments, computer simulations and statistical physics in the investigation of the effect of calcium-binding, salt concentration, and macromolecular crowding on CaM. The results revealed unique folding energy landscapes of CaM in the absence and presence of calcium ions and the structural implications of CaM are interpreted under cell-like conditions. Further, a large conformational change in CaM in response to environmental impacts, dictates the packing of local helices that may be critical to its function of target binding and recognition among vast target selections.

Structure, function, and folding of phosphoglycerate kinase are strongly perturbed by macromolecular crowding

We combine experiment and computer simulation to show how macromolecular crowding dramatically affects the structure, function, and folding landscape of phosphoglycerate kinase (PGK). Fluorescence labeling shows that compact states of yeast PGK are populated as the amount of crowding agents (Ficoll 70) increases. Coarse-grained molecular simulations reveal three compact ensembles: C (crystal structure), CC (collapsed crystal), and Sph (spherical compact). With an adjustment for viscosity, crowded wild-type PGK and fluorescent PGK are about 15 times or more active in 200 mg/ml Ficoll than in aqueous solution. Our results suggest a previously undescribed solution to the classic problem of how the ADP and diphosphoglycerate binding sites of PGK come together to make ATP: Rather than undergoing a hinge motion, the ADP and substrate sites are already located in proximity under crowded conditions that mimic the in vivo conditions under which the enzyme actually operates. We also examine T-jump unfolding of PGK as a function of crowding experimentally. We uncover a nonmonotonic folding relaxation time vs. Ficoll concentration. Theory and modeling explain why an optimum concentration exists for fastest folding. Below the optimum, folding slows down because the unfolded state is stabilized relative to the transition state. Above the optimum, folding slows down because of increased viscosity.

A didactic model of macromolecular crowding effects on protein folding

A didactic model is presented to illustrate how the effect of macromolecular crowding on protein folding and association is modeled using current analytical theory and discrete molecular dynamics. While analytical treatments of crowding may consider the effect as a potential of average force acting to compress a polypeptide chain into a compact state, the use of simulations enables the presence of crowding reagents to be treated explicitly. Using an analytically solvable toy model for protein folding, an approximate statistical thermodynamic method is directly compared to simulation in order to gauge the effectiveness of current analytical crowding descriptions. Both methodologies are in quantitative agreement under most conditions, indication that both current theory and simulation methods are capable of recapitulating aspects of protein folding even by utilizing a simplistic protein model.

Diffusion, crowding & protein stability in a dynamic molecular model of the bacterial cytoplasm

A longstanding question in molecular biology is the extent to which the behavior of macromolecules observed in vitro accurately reflects their behavior in vivo. A number of sophisticated experimental techniques now allow the behavior of individual types of macromolecule to be studied directly in vivo; none, however, allow a wide range of molecule types to be observed simultaneously. In order to tackle this issue we have adopted a computational perspective, and, having selected the model prokaryote Escherichia coli as a test system, have assembled an atomically detailed model of its cytoplasmic environment that includes 50 of the most abundant types of macromolecules at experimentally measured concentrations. Brownian dynamics (BD) simulations of the cytoplasm model have been calibrated to reproduce the translational diffusion coefficients of Green Fluorescent Protein (GFP) observed in vivo, and “snapshots” of the simulation trajectories have been used to compute the cytoplasm's effects on the thermodynamics of protein folding, association and aggregation events. The simulation model successfully describes the relative thermodynamic stabilities of proteins measured in E. coli, and shows that effects additional to the commonly cited “crowding” effect must be included in attempts to understand macromolecular behavior in vivo.

The interior of a typical bacterial cell is a highly crowded place in which molecules must jostle and compete with each other in order to carry out their biological functions. The conditions under which such molecules are typically studied in vitro, however, are usually quite different: one or a few different types of molecules are studied as they freely diffuse in a dilute, aqueous solution. There is therefore a significant disconnect between the conditions under which molecules can be most usefully studied and the conditions under which such molecules usually “live”, and developing ways to bridge this gap is likely to be important for properly understanding molecular behavior in vivo. Toward this end, we show in this work that computer simulations can be used to model the interior of bacterial cells at a near atomic level of detail: the rates of diffusion of proteins are matched to known experimental values, and their thermodynamic stabilities are found to be in good agreement with the few measurements that have so far been performed in vivo. While the simulation approach is certainly not free of assumptions, it offers a potentially important complement to experimental techniques and provides a vivid illustration of molecular behavior inside a biological cell that is likely to be of significant educational value.

Models of macromolecular crowding effects and the need for quantitative comparisons with experiment

In recent years significant effort has been devoted to exploring the potential effects of macromolecular crowding on protein folding and association phenomena. Theoretical calculations and molecular simulations have, in particular, been exploited to describe aspects of protein behavior in crowded and confined conditions and many aspects of the simulated behavior have reflected, at least at a qualitative level, the behavior observed in experiments. One major and immediate challenge for the theorists is to now produce models capable of making quantitatively accurate predictions of in vitro behavior. A second challenge is to derive models that explain results obtained from experiments performed in vivo, the results of which appear to call into question the assumed dominance of excluded-volume effects in vivo.