Effects of Macromolecular Crowding on the Conformational Ensembles of Disordered Proteins

Impact of Guanidinium Hydrochloride on the Shapes of Prothymosin-α and α-Synuclein Is Dramatically Different

The effects of guanidinium hydrochloride (GdmCl) on two intrinsically disordered proteins (IDPs) are investigated using simulations of the self-organized polymer-IDP (SOP-IDP) model. The impact of GdmCl is taken into account using the molecular transfer model (MTM). We show that due to the dramatic reduction in the stiffness of the highly charged Prothymosin-α (ProTα) with increasing concentration of GdmCl ([GdmCl]), the radius of gyration (Rg) decreases sharply until about 1.0 M. Above 1.0 M, ProTα expands, caused by the swelling effect of GdmCl. In contrast, Rg of α-Synuclein (αSyn) swells as continuously as [GdmCl] increases, with most of the expansion occurring at concentrations less than 0.2 M. Strikingly, the amplitude of the small-angle X-ray scattering (SAXS) profiles for ProTα increases until [GdmCl] ≈ 1.0 M and decreases beyond 1.0 M. The [GdmCl]-dependent SAXS profiles for αSyn, which has a pronounced bump at small wave vector (q ∼ 0.5 nm-1) at low [GdmCl] (≤0.2 M), monotonically decrease at all values of [GdmCl]. The contrasting behavior predicted by the combination of MTM and SOP-IDP simulations may be qualitatively understood by modeling ProTα as a strongly charged polyelectrolyte with nearly uniform density of charges along the chain contour and αSyn as a nearly neutral polymer, except near the C-terminus, where the uncompensated negatively charged residues are located. The precise predictions for the SAXS profiles as a function of [GdmCl] can be readily tested.

Identifying protein conformational states in the Protein Data Bank: Toward unlocking the potential of integrative dynamics studies

Studying protein dynamics and conformational heterogeneity is crucial for understanding biomolecular systems and treating disease. Despite the deposition of over 215 000 macromolecular structures in the Protein Data Bank and the advent of AI-based structure prediction tools such as AlphaFold2, RoseTTAFold, and ESMFold, static representations are typically produced, which fail to fully capture macromolecular motion. Here, we discuss the importance of integrating experimental structures with computational clustering to explore the conformational landscapes that manifest protein function. We describe the method developed by the Protein Data Bank in Europe - Knowledge Base to identify distinct conformational states, demonstrate the resource's primary use cases, through examples, and discuss the need for further efforts to annotate protein conformations with functional information. Such initiatives will be crucial in unlocking the potential of protein dynamics data, expediting drug discovery research, and deepening our understanding of macromolecular mechanisms.

Design of functional intrinsically disordered proteins

Many proteins do not fold into a fixed three-dimensional structure, but rather function in a highly disordered state. These intrinsically disordered proteins pose a unique challenge to protein engineering and design: How can proteins be designed de novo if not by tailoring their structure? Here, we will review the nascent field of design of intrinsically disordered proteins with focus on applications in biotechnology and medicine. The design goals should not necessarily be the same as for de novo design of folded proteins as disordered proteins have unique functional strengths and limitations. We focus on functions where intrinsically disordered proteins are uniquely suited including disordered linkers, desiccation chaperones, sensors of the chemical environment, delivery of pharmaceuticals, and constituents of biomolecular condensates. Design of functional intrinsically disordered proteins relies on a combination of computational tools and heuristics gleaned from sequence-function studies. There are few cases where intrinsically disordered proteins have made it into industrial applications. However, we argue that disordered proteins can perform many roles currently performed by organic polymers, and that these proteins might be more designable due to their modularity.

Bridging Soft Interaction and Excluded Volume in Crowded Milieu through Subtle Protein Dynamics

The impact of macromolecular crowding on biological macromolecules has been elucidated through the excluded volume phenomenon and soft interactions. However, it has often been difficult to provide a clear demarcation between the two regions. Here, using temperature-dependent dynamics (local and global) of the multidomain protein human serum albumin (HSA) in the presence of commonly used synthetic crowders (Dextran 40, PEG 8, Ficoll 70, and Dextran 70), we have shown the presence of a transition that serves as a bridge between the soft and hard regimes. The bridging region is independent of the crowder identity and displays no apparent correlation with the critical overlap concentration of the polymeric crowding agents. Moreover, the dynamics of domains I and II and the protein gating motion respond differently, thereby bringing to the fore the asymmetry underlying the crowder influence on HSA. In addition, solvent-coupled and decoupled protein motions indicate the heterogeneity of the dynamic landscape in the crowded milieu. We also propose an intriguing correlation between protein stability and dynamics, with increased global stability being accompanied by eased local domain motion.

The conformational phase diagram of charged polymers in the presence of attractive bridging crowders

Using extensive molecular dynamics simulations, we obtain the conformational phase diagram of a charged polymer in the presence of oppositely charged counterions and neutral attractive crowders for monovalent, divalent, and trivalent counterion valencies. We demonstrate that the charged polymer can exist in three phases: (1) an extended phase for low charge densities and weak polymer-crowder attractive interactions [Charged Extended (CE)]; (2) a collapsed phase for high charge densities and weak polymer-crowder attractive interactions, primarily driven by counterion condensation [Charged Collapsed due to Intra-polymer interactions [(CCI)]; and (3) a collapsed phase for strong polymer-crowder attractive interactions, irrespective of the charge density, driven by crowders acting as bridges or cross-links [Charged Collapsed due to Bridging interactions [(CCB)]. Importantly, simulations reveal that the interaction with crowders can induce collapse, despite the presence of strong repulsive electrostatic interactions, and can replace condensed counterions to facilitate a direct transition from the CCI and CE phases to the CCB phase.

Excluded Volume and Weak Interactions in Crowded Solutions Modulate Conformations and RNA Binding of an Intrinsically Disordered Tail

The cellular milieu is a solution crowded with a significant concentration of different components (proteins, nucleic acids, metabolites, etc.). Such a crowded environment affects protein conformations, dynamics, and interactions. Intrinsically disordered proteins and regions are particularly sensitive to these effects. Here, we investigate the impact on an intrinsically disordered tail that flanks a folded domain, the N-terminal domain, and the RNA-binding domain of the SARS-CoV-2 nucleocapsid protein. We mimic the crowded environment of the cell using polyethylene glycol (PEG) and study its impact on protein conformations using single-molecule Förster resonance energy transfer. We found that high-molecular-weight PEG induces a collapse of the disordered N-terminal tail, whereas low-molecular-weight PEG induces a chain expansion. Our data can be explained by accounting for two opposing contributions: favorable interactions between the protein and crowder molecules and screening of excluded volume interactions. We further characterized the interaction between protein and RNA in the presence of crowding agents. While for all PEG molecules tested, we observed an increase in the binding affinity, the trend is not monotonic as a function of the degree of PEG polymerization. This points to the role of nonspecific protein-PEG interactions on binding in addition to the entropic effects due to crowding. To separate the enthalpic and entropic components of the effects, we investigated the temperature dependence of the association constants in the absence and presence of crowders. Finally, we compared the effects of crowding across mutations in the disordered region and found that the threefold difference in association constants for two naturally occurring variants of the SARS-CoV-2 nucleocapsid protein is reduced to almost identical affinities in the presence of crowders. Overall, our data provide new insights into understanding and modeling the contribution of crowding effects on disordered regions, including the impact of interactions between proteins and crowders and their interplay when binding a ligand.

Effects of Homogeneous and Heterogeneous Crowding on Translational Diffusion of Rigid Bovine Serum Albumin and Disordered Alfa-Casein

Intracellular environment includes proteins, sugars, and nucleic acids interacting in restricted media. In the cytoplasm, the excluded volume effect takes up to 40% of the volume available for occupation by macromolecules. In this work, we tested several approaches modeling crowded solutions for protein diffusion. We experimentally showed how the protein diffusion deviates from conventional Brownian motion in artificial conditions modeling the alteration of medium viscosity and rigid spatial obstacles. The studied tracer proteins were globular bovine serum albumin and intrinsically disordered α-casein. Using the pulsed field gradient NMR, we investigated the translational diffusion of protein probes of different structures in homogeneous (glycerol) and heterogeneous (PEG 300/PEG 6000/PEG 40,000) solutions as a function of crowder concentration. Our results showed fundamentally different effects of homogeneous and heterogeneous crowded environments on protein self-diffusion. In addition, the applied "tracer on lattice" model showed that smaller crowding obstacles (PEG 300 and PEG 6000) create a dense net of restrictions noticeably hindering diffusing protein probes, whereas the large-sized PEG 40,000 creates a "less restricted" environment for the diffusive motion of protein molecules.

Conformational Plasticity in α-Synuclein and How Crowded Environment Modulates It

A 140-residue intrinsically disordered protein (IDP), α-synuclein (αS), is known to adopt conformations that are vastly plastic and susceptible to environmental cues and crowders. However, the inherently heterogeneous nature of αS has precluded a clear demarcation of its monomeric precursor between aggregation-prone and functionally relevant aggregation-resistant states and how a crowded environment could modulate their mutual dynamic equilibrium. Here, we identify an optimal set of distinct metastable states of αS in aqueous media by dissecting a 73 μs-long molecular dynamics ensemble via building a comprehensive Markov state model (MSM). Notably, the most populated metastable state corroborates with the dimension obtained from PRE-NMR studies of αS monomer, and it undergoes kinetic transition at diverse time scales with a weakly populated random-coil-like ensemble and a globular protein-like state. However, subjecting αS to a crowded environment results in a nonmonotonic compaction of these metastable conformations, thereby skewing the ensemble by either introducing new tertiary contacts or by reinforcing the innate contacts. The early stage of dimerization process is found to be considerably expedited in the presence of crowders, albeit promoting nonspecific interactions. Together with this, using an extensively sampled ensemble of αS, this exposition demonstrates that crowded environments can potentially modulate the conformational preferences of IDP that can either promote or inhibit aggregation events.

Crowder-directed interactions and conformational dynamics in multistimuli-responsive intrinsically disordered protein

The consequences of crowding on the dynamic conformational ensembles of intrinsically disordered proteins (IDPs) remain unresolved because of their ultrafast motion. Here, we report crowder-induced interactions and conformational dynamics of a prototypical multistimuli-responsive IDP, Rec1-resilin. The effects of a range of crowders of varying sizes, forms, topologies, and concentrations were examined using spectroscopic, spectrofluorimetric, and contrast-matching small- and ultrasmall-angle neutron scattering investigation. To achieve sufficient neutron contrast against the crowders, deuterium-labeled Rec1-resilin was biosynthesized successfully. Moreover, the ab initio "shape reconstruction" approach was used to obtain three-dimensional models of the conformational assemblies. The IDP revealed crowder-specific systematic extension and compaction with the level of macromolecular crowding. Last, a robust extension-contraction model has been postulated to capture the fundamental phenomena governing the observed behavior of IDPs. The study provides insights and fresh perspectives for understanding the interactions and structural dynamics of IDPs in crowded states.

Disordered-Ordered Protein Binary Classification by Circular Dichroism Spectroscopy

Intrinsically disordered proteins lack a stable tertiary structure and form dynamic conformational ensembles due to their characteristic physicochemical properties and amino acid composition. They are abundant in nature and responsible for a large variety of cellular functions. While numerous bioinformatics tools have been developed for in silico disorder prediction in the last decades, there is a need for experimental methods to verify the disordered state. CD spectroscopy is widely used for protein secondary structure analysis. It is usable in a wide concentration range under various buffer conditions. Even without providing high-resolution information, it is especially useful when NMR, X-ray, or other techniques are problematic or one simply needs a fast technique to verify the structure of proteins. Here, we propose an automatized binary disorder-order classification method by analyzing far-UV CD spectroscopy data. The method needs CD data at only three wavelength points, making high-throughput data collection possible. The mathematical analysis applies the k-nearest neighbor algorithm with cosine distance function, which is independent of the spectral amplitude and thus free of concentration determination errors. Moreover, the method can be used even for strong absorbing samples, such as the case of crowded environmental conditions, if the spectrum can be recorded down to the wavelength of 212 nm. We believe the classification method will be useful in identifying disorder and will also facilitate the growth of experimental data in IDP databases. The method is implemented on a webserver and freely available for academic users.

Crowding affects structural dynamics and contributes to membrane association of the NS3/4A complex

Using molecular dynamics simulations, we describe how crowded environments affect the internal dynamics and diffusion of the hepatitis C virus proteases NS3/4A. This protease plays a key role in viral replication and is successfully used as a target for antiviral treatment. The NS3 enzyme requires a peptide cofactor, called NS4A, with its central part interacting with the NS3 β-sheet, and flexible, protruding terminal tails that are unstructured in water solution. The simulations describe the enzyme and water molecules at atomistic resolution, whereas crowders are modeled via either all-atom or coarse-grained models to emphasize different aspects of crowding. Crowders reflect the polyethylene glycol (PEG) molecules used in the experiments to mimic the crowded surrounding. A bead-shell model of folded coarse-grained PEG molecules considers mainly the excluded volume effect, whereas all-atom PEG models afford more protein-like crowder interactions. Circular dichroism spectroscopy experiments of the NS4A N-terminal tail show that a helical structure is formed in the presence of PEG crowders. The simulations suggest that crowding may assist in the formation of an NS4A helical fragment, positioned exactly where a transmembrane helix would fold upon the NS4A contact with the membrane. In addition, partially interactive PEGs help the NS4A N-tail to detach from the protease surface, thus enabling the process of helix insertion and potentially helping the virus establish a replication machinery needed to produce new viruses. Results point to an active role of crowding in assisting structural changes in disordered protein fragments that are necessary for their biological function.

Homo-oligomerization of the human adenosine A2A receptor is driven by the intrinsically disordered C-terminus

G protein-coupled receptors (GPCRs) have long been shown to exist as oligomers with functional properties distinct from those of the monomeric counterparts, but the driving factors of oligomerization remain relatively unexplored. Herein, we focus on the human adenosine A2A receptor (A2AR), a model GPCR that forms oligomers both in vitro and in vivo. Combining experimental and computational approaches, we discover that the intrinsically disordered C-terminus of A2AR drives receptor homo-oligomerization. The formation of A2AR oligomers declines progressively with the shortening of the C-terminus. Multiple interaction types are responsible for A2AR oligomerization, including disulfide linkages, hydrogen bonds, electrostatic interactions, and hydrophobic interactions. These interactions are enhanced by depletion interactions, giving rise to a tunable network of bonds that allow A2AR oligomers to adopt multiple interfaces. This study uncovers the disordered C-terminus as a prominent driving factor for the oligomerization of a GPCR, offering important insight into the effect of C-terminus modification on receptor oligomerization of A2AR and other GPCRs reconstituted in vitro for biophysical studies.

Crowding Effects on the Structure and Dynamics of the Intrinsically Disordered Nuclear Chromatin Protein NUPR1

The intracellular environment is crowded with macromolecules, including sugars, proteins and nucleic acids. In the cytoplasm, crowding effects are capable of excluding up to 40% of the volume available to any macromolecule when compared to dilute conditions. NUPR1 is an intrinsically disordered protein (IDP) involved in cell-cycle regulation, stress-cell response, apoptosis processes, DNA binding and repair, chromatin remodeling and transcription. Simulations of molecular crowding predict that IDPs can adopt compact states, as well as more extended conformations under crowding conditions. In this work, we analyzed the conformation and dynamics of NUPR1 in the presence of two synthetic polymers, Ficoll-70 and Dextran-40, which mimic crowding effects in the cells, at two different concentrations (50 and 150 mg/ml). The study was carried out by using a multi-spectroscopic approach, including: site-directed spin labelling electron paramagnetic resonance spectroscopy (SDSL-EPR), nuclear magnetic resonance spectroscopy (NMR), circular dichroism (CD), small angle X-ray scattering (SAXS) and dynamic light scattering (DLS). SDSL-EPR spectra of two spin-labelled mutants indicate that there was binding with the crowders and that the local dynamics of the C and N termini of NUPR1 were partially affected by the crowders. However, the overall disordered nature of NUPR1 did not change substantially in the presence of the crowders, as shown by circular dichroism CD and NMR, and further confirmed by EPR. The changes in the dynamics of the paramagnetic probes appear to be related to preferred local conformations and thus crowding agents partially affect some specific regions, further pinpointing that NUPR1 flexibility has a key physiological role in its activity.

Protein folding and assembly in confined environments: Implications for protein aggregation in hydrogels and tissues

In the biological milieu of a cell, soluble crowding molecules and rigid confined environments strongly influence whether the protein is properly folded, intrinsically disordered proteins assemble into distinct phases, or a denatured or aggregated protein species is favored. Such crowding and confinement factors act to exclude solvent volume from the protein molecules, resulting in an increased local protein concentration and decreased protein entropy. A protein's structure is inherently tied to its function. Examples of processes where crowding and confinement may strongly influence protein function include transmembrane protein dimerization, enzymatic activity, assembly of supramolecular structures (e.g., microtubules), nuclear condensates containing transcriptional machinery, protein aggregation in the contexts of disease and protein therapeutics. Historically, most protein structures have been determined from pure, dilute protein solutions or pure crystals. However, these are not the environments in which these proteins function. Thus, there has been an increased emphasis on analyzing protein structure and dynamics in more "in vivo-like" environments. Complex in vitro models using hydrogel scaffolds to study proteins may better mimic features of the in vivo environment. Therefore, analytical techniques need to be optimized for real-time analysis of proteins within hydrogel scaffolds.

Depletion interactions modulate the binding between disordered proteins in crowded environments

The molecular environment in a biological cell is much more crowded than the conditions commonly used in biochemical and biophysical experiments in vitro. It is therefore important to understand how the conformations and interactions of biological macromolecules are affected by such crowding. Addressing these questions quantitatively, however, has been challenging owing to a lack of sufficiently detailed experimental information and theoretical concepts suitable for describing crowding, especially when polymeric crowding agents and biomolecules are involved. Here, we use the combination of extensive single-molecule experiments with established and recent theoretical concepts to investigate the interaction between two intrinsically disordered proteins. We observe pronounced effects of crowding on their interactions and provide a quantitative framework for rationalizing these effects.

Intrinsically disordered proteins (IDPs) abound in cellular regulation. Their interactions are often transitory and highly sensitive to salt concentration and posttranslational modifications. However, little is known about the effect of macromolecular crowding on the interactions of IDPs with their cellular targets. Here, we investigate the influence of crowding on the interaction between two IDPs that fold upon binding, with polyethylene glycol as a crowding agent. Single-molecule spectroscopy allows us to quantify the effects of crowding on a comprehensive set of observables simultaneously: the equilibrium stability of the complex, the association and dissociation kinetics, and the microviscosity, which governs translational diffusion. We show that a quantitative and coherent explanation of all observables is possible within the framework of depletion interactions if the polymeric nature of IDPs and crowders is incorporated based on recent theoretical developments. The resulting integrated framework can also rationalize important functional consequences, for example, that the interaction between the two IDPs is less enhanced by crowding than expected for folded proteins of the same size.

Nanoscale Surface-Induced Unfolding of Single Fibronectin Is Restricted by Serum Albumin Crowding

Understanding nanoscale protein conformational changes at solid-liquid interfaces is critical for predicting how proteins will impact the performance of biomaterials in vivo. Crowding is an important contributor to conformational stability. Here we apply single-molecule high resolution imaging with photobleaching to directly measure dye-conjugated fibronectin's unfolding in varying conditions of crowding with human serum albumin on aminosilanized glass. Using this approach, we identify serum albumin's crowding mechanism. We find that fibronectin achieves larger degrees of unfolding when not crowded by coadsorbed serum albumin. Serum albumin does not as effectively constrict fibronectin's conformation if it is sequentially, rather than simultaneously, introduced, suggesting that serum albumin's crowding mechanism is dependent on its ability to sterically block fibronectin's unfolding during the process of adsorption. Because fibronectin's conformation is dependent on interfacial macromolecular crowding under in vitro conditions, it is important to consider the role of in vivo crowding on protein activity.

Modeling Crowded Environment in Molecular Simulations

Biomolecules perform their various functions in living cells, namely in an environment that is crowded by many macromolecules. Thus, simulating the dynamics and interactions of biomolecules should take into account not only water and ions but also other binding partners, metabolites, lipids and macromolecules found in cells. In the last decade, research on how to model macromolecular crowders around proteins in order to simulate their dynamics in models of cellular environments has gained a lot of attention. In this mini-review we focus on the models of crowding agents that have been used in computer modeling studies of proteins and peptides, especially via molecular dynamics simulations.

Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation

Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.

Intrinsically Disordered Protein Exhibits Both Compaction and Expansion under Macromolecular Crowding

Conformational malleability allows intrinsically disordered proteins (IDPs) to respond agilely to their environments, such as nonspecifically interacting with in vivo bystander macromolecules (or crowders). Previous studies have emphasized conformational compaction of IDPs due to steric repulsion by macromolecular crowders, but effects of soft attraction are largely unexplored. Here we studied the conformational ensembles of the IDP FlgM in both polymer and protein crowders by small-angle neutron scattering. As crowder concentrations increased, the mean radius of gyration of FlgM first decreased but then exhibited an uptick. Ensemble optimization modeling indicated that FlgM conformations under protein crowding segregated into two distinct populations, one compacted and one extended. Coarse-grained simulations showed that compacted conformers fit into an interstitial void and occasionally bind to a surrounding crowder, whereas extended conformers snake through interstitial crevices and bind multiple crowders simultaneously. Crowder-induced conformational segregation may facilitate various cellular functions of IDPs.

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.

Atomistic Modeling of Intrinsically Disordered Proteins Under Polyethylene Glycol Crowding: Quantitative Comparison with Experimental Data and Implication of Protein-Crowder Attraction

The malleability of intrinsically disordered proteins (IDPs) has generated great interest in understanding how their conformations respond to crowded cellular environments. Experiments can report gross properties such as fluorescence resonance energy transfer (FRET) efficiency but cannot resolve the conformational ensembles of IDPs and their interactions with macromolecular crowders. Computation can in principle provide the latter information but in practice has been hampered by the enormous expense for realistic modeling of IDPs and crowders and for sufficient conformational sampling. Here, taking advantage of a powerful method called FMAP (fast Fourier transform-based modeling of atomistic protein-crowder interactions), we computed how the conformational ensembles of three IDPs are modified in concentrated polyethylene glycol (PEG) 6000 solutions. We represented the IDPs at the all-atom level and the PEG molecules at a coarse-grained level and calculated the experimental observable, i.e., FRET efficiency. Whereas accounting for only steric repulsion of PEG led to overestimation of crowding effects, quantitative agreement with experimental data was obtained upon including mild IDP-PEG attraction. The present work demonstrates that realistic modeling of IDPs under crowded conditions for direct comparison with experiments is now achievable.

Molecular simulations of cellular processes

It is, nowadays, possible to simulate biological processes in conditions that mimic the different cellular compartments. Several groups have performed these calculations using molecular models that vary in performance and accuracy. In many cases, the atomistic degrees of freedom have been eliminated, sacrificing both structural complexity and chemical specificity to be able to explore slow processes. In this review, we will discuss the insights gained from computer simulations on macromolecule diffusion, nuclear body formation, and processes involving the genetic material inside cell-mimicking spaces. We will also discuss the challenges to generate new models suitable for the simulations of biological processes on a cell scale and for cell-cycle-long times, including non-equilibrium events such as the co-translational folding, misfolding, and aggregation of proteins. A prominent role will be played by the wise choice of the structural simplifications and, simultaneously, of a relatively complex energetic description. These challenging tasks will rely on the integration of experimental and computational methods, achieved through the application of efficient algorithms. Graphical abstract.

Predictions of Backbone Dynamics in Intrinsically Disordered Proteins Using De Novo Fragment-Based Protein Structure Predictions

Intrinsically disordaered proteins (IDPs) are a prevalent phenomenon with over 30% of human proteins estimated to have long disordered regions. Computational methods are widely used to study IDPs, however, nearly all treat disorder in a binary fashion, not accounting for the structural heterogeneity present in disordered regions. Here, we present a new de novo method, FRAGFOLD-IDP, which addresses this problem. Using 200 protein structural ensembles derived from NMR, we show that FRAGFOLD-IDP achieves superior results compared to methods which can predict related data (NMR order parameter, or crystallographic B-factor). FRAGFOLD-IDP produces very good predictions for 33.5% of cases and helps to get a better insight into the dynamics of the disordered ensembles. The results also show it is not necessary to predict the correct fold of the protein to reliably predict per-residue fluctuations. It implies that disorder is a local property and it does not depend on the fold. Our results are orthogonal to DynaMine, the only other method significantly better than the naïve prediction. We therefore combine these two using a neural network. FRAGFOLD-IDP enables better insight into backbone dynamics in IDPs and opens exciting possibilities for the design of disordered ensembles, disorder-to-order transitions, or design for protein dynamics.

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.

Protein Composition Determines the Effect of Crowding on the Properties of Disordered Proteins

Unlike dilute experimental conditions under which biological molecules are typically characterized, the cell interior is crowded by macromolecules, which affects both the thermodynamics and kinetics of in vivo processes. Although the excluded-volume effects of macromolecular crowding are expected to cause compaction of unfolded and disordered proteins, the extent of this effect is uncertain. We use a coarse-grained model to represent proteins with varying sequence content and directly observe changes in chain dimensions in the presence of purely repulsive spherical crowders. We find that the extent of crowding-induced compaction is dependent not only on crowder size and concentration, but also on the properties of the protein itself. In fact, we observe a nonmonotonic trend between the dimensions of the polypeptide chain in bulk and the degree of compaction: the most extended chains experience up to 24% compaction, the most compact chains show virtually no change, and intermediate chains compress by up to 40% in size at a 40% crowder volume fraction. Free-volume theory combined with an impenetrable ellipsoidal representation of the chains predicts the crowding effects only for collapsed protein chains. An additional scaling factor, which can be easily computed from protein-crowder potential of mean force, corrects for the penetrability of extended chains and is sufficient to capture the observed nonmonotonic trend in compaction.

Conditional Chaperone-Client Interactions Revealed by Genetically Encoded Photo-cross-linkers

The cell envelope is an integral and essential component of Gram-negative bacteria. As the front line during host-pathogen interactions, it is directly challenged by host immune responses as well as other harsh extracellular stimuli. The high permeability of the outer-membrane and the lack of ATP energy system render it difficult to maintain important biological activities within the periplasmic space under stress conditions. The HdeA/B chaperone machinery is the only known acid resistant system found in bacterial periplasm, enabling enteric pathogens to survive through the highly acidic human stomach and establish infections in the intestine. These two homologous chaperones belong to a fast growing family of conditionally disordered chaperones that conditionally lose their well-defined three-dimensional structures to exert biological activities. Upon losing ordered structures, these proteins commit promiscuous binding of diverse clients in response to environmental stimulation. For example, HdeA and HdeB are well-folded inactive dimers at neutral pH but become partially unfolded to protect a wide array of acid-denatured proteins upon acid stress. Whether these conditionally disordered chaperones possess client specificities remains unclear. This is in part due to the lack of efficient tools to investigate such versatile and heterogeneous protein-protein interactions under living conditions. Genetically encoded protein photo-cross-linkers have offered a powerful strategy to capture protein-protein interactions, showing great potential in profiling protein interaction networks, mapping binding interfaces, and probing dynamic changes in both physiological and pathological settings. Despite great success, photo-cross-linkers that can simultaneously capture the promiscuous binding partners and directly identify the interaction interfaces remain technically challenging. Furthermore, methods for side-by-side profiling and comparing the condition-dependent client pools from two homologous chaperones are lacking. Herein, we introduce our recent efforts in developing a panel of versatile genetically encoded photo-cross-linkers to study the disorder-mediated chaperone-client interactions in living cells. In particular, we have developed a series of proteomic-based strategies relying on these new photo-cross-linkers to systematically compare the client profiles of HdeA and HdeB, as well as to map their interaction interfaces. These studies revealed the mode-of-action, particularly the client specificity, of these two conditionally disordered chaperones. In the end, some recent elegant work from other groups that applied the genetically encoded photo-cross-linking strategy to illuminate important protein-protein interactions within bacterial cell envelope is also discussed.

Combining Diffusion NMR and Small-Angle Neutron Scattering Enables Precise Measurements of Polymer Chain Compression in a Crowded Environment

The effect of particles on the behavior of polymers in solution is important in a number of important phenomena such as the effect of "crowding" proteins in cells, colloid-polymer mixtures, and nanoparticle "fillers" in polymer solutions and melts. In this Letter, we study the effect of spherical inert nanoparticles (which we refer to as "crowders") on the diffusion coefficient and radius of gyration of polymers in solution using pulsed-field-gradient NMR and small-angle neutron scattering (SANS), respectively. The diffusion coefficients exhibit a plateau below a characteristic polymer concentration, which we identify as the overlap threshold concentration c^{⋆}. Above c^{⋆}, in a crossover region between the dilute and semidilute regimes, the (long-time) self-diffusion coefficients are found, universally, to decrease exponentially with polymer concentration at all crowder packing fractions, consistent with a structural basis for the long-time dynamics. The radius of gyration obtained from SANS in the crossover regime changes linearly with an increase in polymer concentration, and must be extrapolated to c^{⋆} in order to obtain the radius of gyration of an individual polymer chain. When the polymer radius of gyration and crowder size are comparable, the polymer size is very weakly affected by the presence of crowders, consistent with recent computer simulations. There is significant chain compression, however, when the crowder size is much smaller than the polymer radius gyration.

Simulations of disordered proteins and systems with conformational heterogeneity

Intrinsically disordered proteins (IDPs) and protein regions can facilitate a wide variety of complex physiological processes such as binding, signaling, and formation of membraneless organelles. They can however also play pathological roles by aggregating into cytotoxic oligomers and fibrils. Characterizing the structure and function of disordered proteins is an onerous task, primarily because these proteins adopt transient structures, which are difficult to capture in experiments. Simulations have emerged as a powerful tool for interpreting and augmenting experimental measurements of IDPs. In this review we focus on computer simulations of disordered protein structures, functions, assemblies, and emerging questions that, taken together, give an overview of the field as it exists today.

Macromolecular and Small Molecular Crowding Have Similar Effects on α-Synuclein Structure

The intracellular milieu contains upwards of 400 g of macromolecules per liter. This crowding is thought to have a larger influence on intrinsically disordered proteins, whose chains are expanded, than on compact globular proteins. Classic theories of macromolecular crowding predict that increasing excluded volume effects will lead disordered proteins to compaction, and a great deal of data, from both simulation and experiments support this idea. We used nuclear magnetic resonance, circular dichroism, and fluorescence spectroscopies to characterize the structure and fibrillation of α-synuclein, an intrinsically disordered protein implicated in Parkinson's disease, using Ficoll70, its monomer sucrose and bovine serum albumin as crowding agents. Surprisingly, volume exclusion induced by high concentrations of macromolecules may not be the main reason for the compaction of α-synuclein. Our results indicate that all aspects crowding must be considered to understand protein conformation under crowded conditions.

Folding propensity of intrinsically disordered proteins by osmotic stress

Proteins imparted with intrinsic disorder conduct a range of essential cellular functions. To better understand the folding and hydration properties of intrinsically disordered proteins (IDPs), we used osmotic stress to induce conformational changes in nuclear co-activator binding domain (NCBD) and activator for thyroid hormone and retinoid receptor (ACTR) separate from their mutual binding. Osmotic stress was applied by the addition of small and polymeric osmolytes, where we discovered that water contributions to NCBD folding always exceeded those for ACTR. Both NCBD and ACTR were found to gain α-helical structure with increasing osmotic stress, consistent with their folding upon NCBD/ACTR complex formation. Using small-angle neutron scattering (SANS), we further characterized NCBD structural changes with the osmolyte ethylene glycol. Here a large reduction in overall size initially occurred before substantial secondary structural change. By focusing on folding propensity, and linked hydration changes, we uncover new insights that may be important for how IDP folding contributes to binding.

Protein folding, binding, and droplet formation in cell-like conditions

The many bystander macromolecules in the crowded cellular environments present both steric repulsion and weak attraction to proteins undergoing folding or binding and hence impact the thermodynamic and kinetic properties of these processes. The weak but nonrandom binding with bystander macromolecules may facilitate subcellular localization and biological function. Weak binding also leads to the emergence of a protein-rich droplet phase, which has been implicated in regulating a variety of cellular functions. All these important problems can now be addressed by realistic modeling of intermolecular interactions. Configurational sampling of concentrated protein solutions is an ongoing challenge.

The Differential Response of Proteins to Macromolecular Crowding

The habitat in which proteins exert their function contains up to 400 g/L of macromolecules, most of which are proteins. The repercussions of this dense environment on protein behavior are often overlooked or addressed using synthetic agents such as poly(ethylene glycol), whose ability to mimic protein crowders has not been demonstrated. Here we performed a comprehensive atomistic molecular dynamic analysis of the effect of protein crowders on the structure and dynamics of three proteins, namely an intrinsically disordered protein (ACTR), a molten globule conformation (NCBD), and a one-fold structure (IRF-3) protein. We found that crowding does not stabilize the native compact structure, and, in fact, often prevents structural collapse. Poly(ethylene glycol) PEG500 failed to reproduce many aspects of the physiologically-relevant protein crowders, thus indicating its unsuitability to mimic the cell interior. Instead, the impact of protein crowding on the structure and dynamics of a protein depends on its degree of disorder and results from two competing effects: the excluded volume, which favors compact states, and quinary interactions, which favor extended conformers. Such a viscous environment slows down protein flexibility and restricts the conformational landscape, often biasing it towards bioactive conformations but hindering biologically relevant protein-protein contacts. Overall, the protein crowders used here act as unspecific chaperons that modulate the protein conformational space, thus having relevant consequences for disordered proteins.

Most in vitro and in silico biophysical experiments generally study proteins in an isolated environment, overlooking that their natural environment—the cell cytoplasm—is a solution that is highly populated by proteins. To address this knowledge gap, here we explored how a crowded environment alters the conformational sampling of three proteins, each with a different degree of disorder and flexibility. We simulated a crowded system composed by the three proteins and reaching a cell-like concentration and compared the protein behavior observed with that induced by PEG500, a synthetic crowding agent. Despite some similarities between the environments, protein crowders showed a number of characteristics that raise concerns about the use of diluted solutions or synthetic agents when studying protein behavior.

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.

Tick receptor for outer surface protein A from Ixodes ricinus - the first intrinsically disordered protein involved in vector-microbe recognition

The tick receptor for outer surface protein A (TROSPA) is the only identified factor involved in tick gut colonization by various Borrelia species. TROSPA is localized in the gut epithelium and can recognize and bind the outer surface bacterial protein OspA via an unknown mechanism. Based on earlier reports and our latest observations, we considered that TROSPA would be the first identified intrinsically disordered protein (IDP) involved in the interaction between a vector and a pathogenic microbe. To verify this hypothesis, we performed structural studies of a TROSPA mutant from Ixodes ricinus using both computational and experimental approaches. Irrespective of the method used, we observed that the secondary structure content of the TROSPA polypeptide chain is low. In addition, the collected SAXS data indicated that this protein is highly extended and exists in solution as a set of numerous conformers. These features are all commonly considered hallmarks of IDPs. Taking advantage of our SAXS data, we created structural models of TROSPA and proposed a putative mechanism for the TROSPA-OspA interaction. The disordered nature of TROSPA may explain the ability of a wide spectrum of Borrelia species to colonize the tick gut.

Concentrated Solutions of Single-Chain Nanoparticles: A Simple Model for Intrinsically Disordered Proteins under Crowding Conditions

By means of large-scale computer simulations and small-angle neutron scattering (SANS), we investigate solutions of single-chain nanoparticles (SCNPs), covering the whole concentration range from infinite dilution to melt density. The analysis of the conformational properties of the SCNPs reveals that these synthetic nano-objects share basic ingredients with intrinsically disordered proteins (IDPs), as topological polydispersity, generally sparse conformations, and locally compact domains. We investigate the role of the architecture of the SCNPs in their collapse behavior under macromolecular crowding. Unlike in the case of linear macromolecules, which experience the usual transition from self-avoiding to Gaussian random-walk conformations, crowding leads to collapsed conformations of SCNPs resembling those of crumpled globules. This behavior is already found at volume fractions (about 30%) that are characteristic of crowding in cellular environments. The simulation results are confirmed by the SANS experiments. Our results for SCNPs--a model system free of specific interactions--propose a general scenario for the effect of steric crowding on IDPs: collapse from sparse conformations at high dilution to crumpled globular conformations in cell environments.

Dynamic footprint of sequestration in the molecular fluctuations of osteopontin

The sequestration of calcium phosphate by unfolded proteins is fundamental to the stabilization of biofluids supersaturated with respect to hydroxyapatite, such as milk, blood or urine. The unfolded state of osteopontin (OPN) is thought to be a prerequisite for this activity, which leads to the formation of core-shell calcium phosphate nanoclusters. We report on the structures and dynamics of a native OPN peptide from bovine milk, studied by neutron spectroscopy and small-angle X-ray and neutron scattering. The effects of sequestration are quantified on the nanosecond- ångström resolution by elastic incoherent neutron scattering. The molecular fluctuations of the free phosphopeptide are in agreement with a highly flexible protein. An increased resilience to diffusive motions of OPN is corroborated by molecular fluctuations similar to those observed for globular proteins, yet retaining conformational flexibilities. The results bring insight into the modulation of the activity of OPN and phosphopeptides with a role in the control of biomineralization. The quantification of such effects provides an important handle for the future design of new peptides based on the dynamics-activity relationship.

Digested disorder, Quarterly intrinsic disorder digest (October-November-December, 2013)

This is the 4th issue of the Digested Disorder series that represents reader's digest of the scientific literature on intrinsically disordered proteins. The only 2 criteria for inclusion in this digest are the publication date (a paper should be published within the covered time frame) and topic (a paper should be dedicated to any aspect of protein intrinsic disorder). The current digest issue covers papers published during the fourth quarter of 2013; i.e. during the period of October, November, and December of 2013. Similar to previous issues, the papers are grouped hierarchically by topics they cover, and for each of the included paper a short description is given on its major findings.

Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding

The diverse biological functions of intrinsically disordered proteins (IDPs) have markedly raised our appreciation of protein conformational versatility, whereas the existence of energetically favorable yet functional detrimental nonnative interactions underscores the physical limitations of evolutionary optimization. Here we survey recent advances in using biophysical modeling to gain insight into experimentally observed nonnative behaviors and IDP properties. Simulations of IDP interactions to date focus mostly on coupled folding-binding, which follows essentially the same organizing principle as the local-nonlocal coupling mechanism in cooperative folding of monomeric globular proteins. By contrast, more innovative theories of electrostatic and aromatic interactions are needed for the conceptually novel but less-explored 'fuzzy' complexes in which the functionally bound IDPs remain largely disordered.

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.

Further Development of the FFT-based Method for Atomistic Modeling of Protein Folding and Binding under Crowding: Optimization of Accuracy and Speed

Recently, we (Qin, S.; Zhou, H. X. J. Chem. Theory Comput.2013, 9, 4633–4643) developed the FFT-based method for Modeling Atomistic Proteins–crowder interactions, henceforth FMAP. Given its potential wide use for calculating effects of crowding on protein folding and binding free energies, here we aimed to optimize the accuracy and speed of FMAP. FMAP is based on expressing protein–crowder interactions as correlation functions and evaluating the latter via fast Fourier transform (FFT). The numerical accuracy of FFT improves as the grid spacing for discretizing space is reduced, but at increasing computational cost. We sought to speed up FMAP calculations by using a relatively coarse grid spacing of 0.6 Å and then correcting for discretization errors. This strategy was tested for different types of interactions (hard-core repulsion, nonpolar attraction, and electrostatic interaction) and over a wide range of protein–crowder systems. We were able to correct for the numerical errors on hard-core repulsion and nonpolar attraction by an 8% inflation of atomic hard-core radii and on electrostatic interaction by a 5% inflation of the magnitudes of protein atomic charges. The corrected results have higher accuracy and enjoy a speedup of more than 100-fold over those obtained using a fine grid spacing of 0.15 Å. With this optimization of accuracy and speed, FMAP may become a practical tool for realistic modeling of protein folding and binding in cell-like environments.