Folding propensity of intrinsically disordered proteins by osmotic stress

Structural characterization of an intrinsically disordered protein complex using integrated small-angle neutron scattering and computing

Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure-function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Lastly, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.

Mechanosensitive ion channels MSL8, MSL9, and MSL10 have environmentally sensitive intrinsically disordered regions with distinct biophysical characteristics in vitro

Intrinsically disordered protein regions (IDRs) are highly dynamic sequences that rapidly sample a collection of conformations over time. In the past several decades, IDRs have emerged as a major component of many proteomes, comprising ~30% of all eukaryotic protein sequences. Proteins with IDRs function in a wide range of biological pathways and are notably enriched in signaling cascades that respond to environmental stresses. Here, we identify and characterize intrinsic disorder in the soluble cytoplasmic N-terminal domains of MSL8, MSL9, and MSL10, three members of the MscS-like (MSL) family of mechanosensitive ion channels. In plants, MSL channels are proposed to mediate cell and organelle osmotic homeostasis. Bioinformatic tools unanimously predicted that the cytosolic N-termini of MSL channels are intrinsically disordered. We examined the N-terminus of MSL10 (MSL10N) as an exemplar of these IDRs and circular dichroism spectroscopy confirms its disorder. MSL10N adopted a predominately helical structure when exposed to the helix-inducing compound trifluoroethanol (TFE). Furthermore, in the presence of molecular crowding agents, MSL10N underwent structural changes and exhibited alterations to its homotypic interaction favorability. Lastly, interrogations of collective behavior via in vitro imaging of condensates indicated that MSL8N, MSL9N, and MSL10N have sharply differing propensities for self-assembly into condensates, both inherently and in response to salt, temperature, and molecular crowding. Taken together, these data establish the N-termini of MSL channels as intrinsically disordered regions with distinct biophysical properties and the potential to respond uniquely to changes in their physiochemical environment.

Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations

Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.

Molecular simulations of IDPs: From ensemble generation to IDP interactions leading to disorder-to-order transitions

Intrinsically disordered proteins (IDPs) lack a well-defined three-dimensional structure but do exhibit some dynamical and structural ordering. The structural plasticity of IDPs indicates that entropy-driven motions are crucial for their function. Many IDPs undergo function-related disorder-to-order transitions upon by their interaction with specific binding partners. Approaches that are based on both experimental and theoretical tools enable the biophysical characterization of IDPs. Molecular simulations provide insights into IDP structural ensembles and disorder-to-order transition mechanisms. However, such studies depend strongly on the chosen force field parameters and simulation techniques. In this chapter, we provide an overview of IDP characteristics, review all-atom force fields recently developed for IDPs, and present molecular dynamics-based simulation methods that allow IDP ensemble generation as well as the characterization of disorder-to-order transitions. In particular, we introduce metadynamics, replica exchange molecular dynamics simulations, and also kinetic models resulting from Markov State modeling, and provide various examples for the successful application of these simulation methods to IDPs.

Revealing the Hidden Sensitivity of Intrinsically Disordered Proteins to their Chemical Environment

Intrinsically disordered protein-regions (IDRs) make up roughly 30% of the human proteome and are central to a wide range of biological processes. Given a lack of persistent tertiary structure, all residues in IDRs are, to some extent, solvent exposed. This extensive surface area, coupled with the absence of strong intramolecular contacts, makes IDRs inherently sensitive to their chemical environment. We report a combined experimental, computational, and analytical framework for high-throughput characterization of IDR sensitivity. Our framework reveals that IDRs can expand or compact in response to changes in their solution environment. Importantly, the direction and magnitude of conformational change depend on both protein sequence and cosolute identity. For example, some solutes such as short polyethylene glycol chains exert an expanding effect on some IDRs and a compacting effect on others. Despite this complex behavior, we can rationally interpret IDR responsiveness to solution composition changes using relatively simple polymer models. Our results imply that solution-responsive IDRs are ubiquitous and can provide an additional layer of regulation to biological systems.

NFAT5, which protects against hypertonicity, is activated by that stress via structuring of its intrinsically disordered domain

Nuclear Factor of Activated T cells 5 (NFAT5) is a transcription factor (TF) that mediates protection from adverse effects of hypertonicity by increasing transcription of genes, including those that lead to cellular accumulation of protective organic osmolytes. NFAT5 has three intrinsically ordered (ID) activation domains (ADs). Using the NFAT5 N-terminal domain (NTD), which contains AD1, as a model, we demonstrate by biophysical methods that the NTD senses osmolytes and hypertonicity, resulting in stabilization of its ID regions. In the presence of sufficient NaCl or osmolytes, trehalose and sorbitol, the NFAT5 NTD undergoes a disorder-to-order shift, adopting higher average secondary and tertiary structure. Thus, NFAT5 is activated by the stress that it protects against. In its salt and/or osmolyte-induced more ordered conformation, the NTD interacts with several proteins, including HMGI-C, which is known to protect against apoptosis. These findings raise the possibility that the increased intracellular ionic strength and elevated osmolytes caused by hypertonicity activate and stabilize NFAT5.

Characterization of an Extensive Interface on Vitronectin for Binding to Plasminogen Activator Inhibitor-1: Adoption of Structure in an Intrinsically Disordered Region

Small-angle neutron scattering (SANS) measurements were pursued to study human vitronectin, a protein found in tissues and the circulation that regulates cell adhesion/migration and proteolytic cascades that govern hemostasis and pericellular proteolysis. Many of these functions occur via interactions with its binding partner, plasminogen activator inhibitor-1 (PAI-1), the chief inhibitor of proteases that lyse and activate plasminogen. We focused on a region of vitronectin that remains uncharacterized from previous X-ray scattering, nuclear magnetic resonance, and computational modeling approaches and which we propose is involved in binding to PAI-1. This region, which bridges the N-terminal somatomedin B (SMB) domain with a large central β-propeller domain of vitronectin, appears unstructured and has characteristics of an intrinsically disordered domain (IDD). The effect of osmolytes was evaluated using circular dichroism and SANS to explore the potential of the IDD to undergo a disorder-to-order transition. The results suggest that the IDD favors a more ordered structure under osmotic pressure; SANS shows a smaller radius of gyration (R_g) and a more compact fold of the IDD upon addition of osmolytes. To test whether PAI-1 binding is also coupled to folding within the IDD structure, a set of SANS experiments with contrast variation were performed on the complex of PAI-1 with a vitronectin fragment corresponding to the N-terminal 130 amino acids (denoted the SMB-IDD because it contains the SMB domain and IDD in linear sequence). Analysis of the SANS data using the Ensemble Optimization Method confirms that the SMB-IDD adopts a more compact configuration when bound to PAI-1. Calculated structures for the PAI-1:SMB-IDD complex suggest that the IDD provides an interaction surface outside of the primary PAI-1-binding site located within the SMB domain; this binding is proposed to lead to the assembly of higher-order structures of vitronectin and PAI-1 commonly found in tissues.

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

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

Does water stress promote the proteome-wide adjustment of intrinsically disordered proteins in plants?

Plant response to water stress involves the activation of mechanisms expected to help them cope with water scarcity. Among these mechanisms, proteome-wide adjustment is well known. This includes actions to save energy, protect cellular and molecular components, and maintain vital functions of the cell. Intrinsically disordered proteins, which are proteins without a rigid three-dimensional structure, are seen as emerging multifunctional cellular components of proteomes. They are highly abundant in eukaryotic proteomes, and numerous functions for these proteins have been proposed. Here, we discuss several reasons why the collection of intrinsically disordered proteins in a proteome (disordome) could be subjected to an active regulation during conditions of water scarcity in plants. We also discuss the potential misinterpretations of disordome content estimations made so far due to bias-prone data and the need for reliable analysis based on experimental data in order to acknowledge the plasticity nature of the disordome.

Metabolic Profile of Oral Squamous Carcinoma Cell Lines Relies on a Higher Demand of Lipid Metabolism in Metastatic Cells

Tumor cells are subjected to a broad range of selective pressures. As a result of the imposed stress, subpopulations of surviving cells exhibit individual biochemical phenotypes that reflect metabolic reprograming. The present work aimed at investigating metabolic parameters of cells displaying increasing degrees of metastatic potential. The metabolites present in cell extracts fraction of tongue fibroblasts and of cell lines derived from human tongue squamous cell carcinoma lineages displaying increasing metastatic potential (SCC9 ZsG, LN1 and LN2) were analyzed by ¹H NMR (nuclear magnetic resonance) spectroscopy. Living, intact cells were also examined by the non-invasive method of fluorescence lifetime imaging microscopy (FLIM) based on the auto fluorescence of endogenous NADH. The cell lines reproducibly exhibited distinct metabolic profiles confirmed by Partial Least-Square Discriminant Analysis (PLS-DA) of the spectra. Measurement of endogenous free and bound NAD(P)H relative concentrations in the intact cell lines showed that ZsG and LN1 cells displayed high heterogeneity in the energy metabolism, indicating that the cells would oscillate between glycolysis and oxidative metabolism depending on the microenvironment’s composition. However, LN2 cells appeared to have more contributions to the oxidative status, displaying a lower NAD(P)H free/bound ratio. Functional experiments of energy metabolism, mitochondrial physiology, and proliferation assays revealed that all lineages exhibited similar energy features, although resorting to different bioenergetics strategies to face metabolic demands. These differentiated functions may also promote metastasis. We propose that lipid metabolism is related to the increased invasiveness as a result of the accumulation of malonate, methyl malonic acid, n-acetyl and unsaturated fatty acids (CH₂)_n in parallel with the metastatic potential progression, thus suggesting that the NAD(P)H reflected the lipid catabolic/anabolic pathways.