Exploring free-energy landscapes of intrinsically disordered proteins at atomic resolution using NMR spectroscopy

Extensive tests and evaluation of the CHARMM36IDPSFF force field for intrinsically disordered proteins and folded proteins

Intrinsically disordered proteins (IDPs) have received increasing attention in recent studies due to their structural heterogeneity and critical biological functions. To fully understand the structural properties and determine accurate ensembles of IDPs, molecular dynamics (MD) simulation was widely used to sample diverse conformations and reveal the structural dynamics. However, the classical state-of-the-art force fields perform well for folded proteins while being unsatisfactory for the simulations of disordered proteins reported in many previous studies. Thus, improved force fields were developed to precisely describe both folded proteins and disordered proteins. Preliminary tests show that our newly developed CHARMM36IDPSFF (C36IDPSFF) force field can well reproduce the experimental observables of several disordered proteins, but more tests on different types of proteins are needed to further evaluate the performance of C36IDPSFF. Here, we extensively simulate short peptides, disordered proteins, and fast-folding proteins as well as folded proteins, and compare the simulated results with the experimental observables. The simulation results show that C36IDPSFF could substantially reproduce the experimental observables for most of the tested proteins but some limitations are also found in the radius of gyration of large disordered proteins and the stability of fast-folding proteins. This force field will facilitate large scale studies of protein structural dynamics and functions using MD simulations.

Perspective: the essential role of NMR in the discovery and characterization of intrinsically disordered proteins

The 2019 ISMAR Prize recognized NMR studies of disordered proteins. Here we provide a highly personal perspective on the discovery of intrinsically disordered proteins and the development and application of NMR methods to characterize their conformational ensembles, dynamics, and interactions.

A Unified Description of Intrinsically Disordered Protein Dynamics under Physiological Conditions Using NMR Spectroscopy

Intrinsically disordered proteins (IDPs) are flexible biomolecules whose essential functions are defined by their dynamic nature. Nuclear magnetic resonance (NMR) spectroscopy is ideally suited to the investigation of this behavior at atomic resolution. NMR relaxation is increasingly used to detect conformational dynamics in free and bound forms of IDPs under conditions approaching physiological, although a general framework providing a quantitative interpretation of these exquisitely sensitive probes as a function of experimental conditions is still lacking. Here, measuring an extensive set of relaxation rates sampling multiple-time-scale dynamics over a broad range of crowding conditions, we develop and test an integrated analytical description that accurately portrays the motion of IDPs as a function of the intrinsic properties of the crowded molecular environment. In particular we observe a strong dependence of both short-range and long-range motional time scales of the protein on the friction of the solvent. This tight coupling between the dynamic behavior of the IDP and its environment allows us to develop analytical expressions for protein motions and NMR relaxation properties that can be accurately applied over a vast range of experimental conditions. This unified dynamic description provides new insight into the physical behavior of IDPs, extending our ability to quantitatively investigate their conformational dynamics under complex environmental conditions, and accurately predicting relaxation rates reporting on motions on time scales up to tens of nanoseconds, both in vitro and in cellulo.

Quantitative Conformational Analysis of Functionally Important Electrostatic Interactions in the Intrinsically Disordered Region of Delta Subunit of Bacterial RNA Polymerase

Electrostatic interactions play important roles in the functional mechanisms exploited by intrinsically disordered proteins (IDPs). The atomic resolution description of long-range and local structural propensities that can both be crucial for the function of highly charged IDPs presents significant experimental challenges. Here, we investigate the conformational behavior of the δ subunit of RNA polymerase from Bacillus subtilis whose unfolded domain is highly charged, with 7 positively charged amino acids followed by 51 acidic amino acids. Using a specifically designed analytical strategy, we identify transient contacts between the two regions using a combination of NMR paramagnetic relaxation enhancements, residual dipolar couplings (RDCs), chemical shifts, and small-angle scattering. This strategy allows the resolution of long-range and local ensemble averaged structural contributions to the experimental RDCs, and reveals that the negatively charged segment folds back onto the positively charged strand, compacting the conformational sampling of the protein while remaining highly flexible in solution. Mutation of the positively charged region abrogates the long-range contact, leaving the disordered domain in an extended conformation, possibly due to local repulsion of like-charges along the chain. Remarkably, in vitro studies show that this mutation also has a significant effect on transcription activity, and results in diminished cell fitness of the mutated bacteria in vivo. This study highlights the importance of accurately describing electrostatic interactions for understanding the functional mechanisms of IDPs.

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

Generation of the configurational ensemble of an intrinsically disordered protein from unbiased molecular dynamics simulation

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.

Bayesian-Maximum-Entropy Reweighting of IDP Ensembles Based on NMR Chemical Shifts

Bayesian and Maximum Entropy approaches allow for a statistically sound and systematic fitting of experimental and computational data. Unfortunately, assessing the relative confidence in these two types of data remains difficult as several steps add unknown error. Here we propose the use of a validation-set method to determine the balance, and thus the amount of fitting. We apply the method to synthetic NMR chemical shift data of an intrinsically disordered protein. We show that the method gives consistent results even when other methods to assess the amount of fitting cannot be applied. Finally, we also describe how the errors in the chemical shift predictor can lead to an incorrect fitting and how using secondary chemical shifts could alleviate this problem.

Novel phosphorylation-dependent regulation in an unstructured protein

This work explores how phosphorylation of an unstructured protein region in inhibitor-2 (I2) regulates protein phosphatase-1 (PP1) enzyme activity using molecular dynamics (MD). Free I2 is largely unstructured; however, when bound to PP1, three segments adopt a stable structure. In particular, an I2 helix (i-helix) blocks the PP1 active site and inhibits phosphatase activity. I2 phosphorylation in the PP1-I2 complex activates phosphatase activity without I2 dissociation. The I2 Thr74 regulatory phosphorylation site is in an unstructured domain in PP1-I2. PP1-I2 MD demonstrated that I2 phosphorylation promotes early steps of PP1-I2 activation in explicit solvent models. Moreover, phosphorylation-dependent activation occurred in PP1-I2 complexes derived from I2 orthologs with diverse sequences from human, yeast, worm, and protozoa. This system allowed exploration of features of the 73-residue unstructured human I2 domain critical for phosphorylation-dependent activation. These studies revealed that components of I2 unstructured domain are strategically positioned for phosphorylation responsiveness including a transient α-helix. There was no evidence that electrostatic interactions of I2 phosphothreonine74 influenced PP1-I2 activation. Instead, phosphorylation altered the conformation of residues around Thr74. Phosphorylation uncurled the distance between I2 residues Glu71 to Tyr76 to promote PP1-I2 activation, whereas reduced distances reduced activation. This I2 residue Glu71 to Tyr76 distance distribution, independently from Thr74 phosphorylation, controls I2 i-helix displacement from the PP1 active site leading to PP1-I2 activation.

Polymer effects modulate binding affinities in disordered proteins

Structural disorder is widespread in regulatory protein networks. Weak and transient interactions render disordered proteins particularly sensitive to fluctuations in solution conditions such as ion and crowder concentrations. How this sensitivity alters folding coupled binding reactions, however, has not been fully understood. Here, we demonstrate that salt jointly modulates polymer properties and binding affinities of 5 disordered proteins from a transcription factor network. A combination of single-molecule Förster resonance energy transfer experiments, polymer theory, and molecular simulations shows that all 5 proteins expand with increasing ionic strengths due to Debye-Hückel charge screening. Simultaneously, pairwise affinities between the proteins increase by an order of magnitude within physiological salt limits. A quantitative analysis shows that 50% of the affinity increase can be explained by changes in the disordered state. Disordered state properties therefore have a functional relevance even if these states are not directly involved in biological functions. Numerical solutions of coupled binding equilibria with our results show that networks of homologous disordered proteins can function surprisingly robustly in fluctuating cellular environments, despite the sensitivity of its individual proteins.

Structure determination using solution NMR: Is it worth the effort?

It has been almost 40 years since solution NMR joined X-ray crystallography as a technique for determining high-resolution structures of proteins. Since then NMR derived structure has contributed in fundamental ways to our understanding of the function of biomolecules. With the already existing mature field of X-ray crystallography and the emergence of cryo-EM as techniques to tackle high-resolution structures of large protein complexes, the role of NMR in structure determination has been questioned. However, NMR has the unique ability to recapitulate the dynamic motion of proteins in their structures, while size limitations of the biomolecular systems that can be routinely studied still present challenges. The field has continually developed methodology and instrumentation since its introduction, pushing its frontiers and redefining its limits. Here we present a brief overview of NMR-based structure determination over the past 40 years. We outline the current state of the field and look ahead to the challenges that still need to be addressed to realize the future potential of NMR as a structural technique.

Structure Determination by Single-Particle Cryo-Electron Microscopy: Only the Sky (and Intrinsic Disorder) is the Limit

Traditionally, X-ray crystallography and NMR spectroscopy represent major workhorses of structural biologists, with the lion share of protein structures reported in protein data bank (PDB) being generated by these powerful techniques. Despite their wide utilization in protein structure determination, these two techniques have logical limitations, with X-ray crystallography being unsuitable for the analysis of highly dynamic structures and with NMR spectroscopy being restricted to the analysis of relatively small proteins. In recent years, we have witnessed an explosive development of the techniques based on Cryo-electron microscopy (Cryo-EM) for structural characterization of biological molecules. In fact, single-particle Cryo-EM is a special niche as it is a technique of choice for the structural analysis of large, structurally heterogeneous, and dynamic complexes. Here, sub-nanometer atomic resolution can be achieved (i.e., resolution below 10 Å) via single-particle imaging of non-crystalline specimens, with accurate 3D reconstruction being generated based on the computational averaging of multiple 2D projection images of the same particle that was frozen rapidly in solution. We provide here a brief overview of single-particle Cryo-EM and show how Cryo-EM has revolutionized structural investigations of membrane proteins. We also show that the presence of intrinsically disordered or flexible regions in a target protein represents one of the major limitations of this promising technique.

The Nucleoprotein and Phosphoprotein of Measles Virus

Measles virus is a negative strand virus and the genomic and antigenomic RNA binds to the nucleoprotein (N), assembling into a helical nucleocapsid. The polymerase complex comprises two proteins, the Large protein (L), that both polymerizes RNA and caps the mRNA, and the phosphoprotein (P) that co-localizes with L on the nucleocapsid. This review presents recent results about N and P, in particular concerning their intrinsically disordered domains. N is a protein of 525 residues with a 120 amino acid disordered C-terminal domain, N_tail. The first 50 residues of N_tail extricate the disordered chain from the nucleocapsid, thereby loosening the otherwise rigid structure, and the C-terminus contains a linear motif that binds P. Recent results show how the 5′ end of the viral RNA binds to N within the nucleocapsid and also show that the bases at the 3′ end of the RNA are rather accessible to the viral polymerase. P is a tetramer and most of the protein is disordered; comprising 507 residues of which around 380 are disordered. The first 37 residues of P bind N, chaperoning against non-specific interaction with cellular RNA, while a second interaction site, around residue 200 also binds N. In addition, there is another interaction between C-terminal domain of P (XD) and N_tail. These results allow us to propose a new model of how the polymerase binds to the nucleocapsid and suggests a mechanism for initiation of transcription.

Conjugation of NMR and SAXS for flexible and multidomain protein structure determination: From sample preparation to model refinement

Experimental information from small angle X-ray scattering (SAXS) is conjugated with nuclear magnetic resonance (NMR) spectroscopy data for the improvement of protein structure determination, particularly for flexible, multidomain or intrinsically disordered proteins. Individually, each of these techniques presents capabilities and limitations: NMR excels in local information, providing atomic resolution, but is limited by protein size, whereas SAXS yields a global envelope of the protein with lower resolution, but revealing domain positions. Different conjugation methodologies use the complementarity of both techniques' independent constraints to achieve comprehensive protein structure determination and resolve dynamics at a moderate computational expense.

Free Ligand 1D NMR Conformational Signatures To Enhance Structure Based Drug Design of a Mcl-1 Inhibitor (AZD5991) and Other Synthetic Macrocycles

The three-dimensional conformations adopted by a free ligand in solution impact bioactivity and physicochemical properties. Solution 1D NMR spectra inherently contain information on ligand conformational flexibility and three-dimensional shape, as well as the propensity of the free ligand to fully preorganize into the bioactive conformation. Herein we discuss some key learnings, distilled from our experience developing potent and selective synthetic macrocyclic inhibitors, including Mcl-1 clinical candidate AZD5991. Case studies have been selected from recent oncology research projects, demonstrating how 1D NMR conformational signatures can complement X-ray protein-ligand structural information to guide medicinal chemistry optimization. Learning to extract free ligand conformational information from routinely available 1D NMR signatures has proven to be fast enough to guide medicinal chemistry decisions within design cycles for compound optimization.

Using Small-Angle Scattering Data and Parametric Machine Learning to Optimize Force Field Parameters for Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) play important roles in many aspects of normal cell physiology, such as signal transduction and transcription, as well as pathological states, including Alzheimer's, Parkinson's, and Huntington's disease. Unlike their globular counterparts that are defined by a few structures and free energy minima, IDP/IDR comprise a large ensemble of rapidly interconverting structures and a corresponding free energy landscape characterized by multiple minima. This aspect has precluded the use of structural biological techniques, such as X-ray crystallography and nuclear magnetic resonance (NMR) for resolving their structures. Instead, low-resolution techniques, such as small-angle X-ray or neutron scattering (SAXS/SANS), have become a mainstay in characterizing coarse features of the ensemble of structures. These are typically complemented with NMR data if possible or computational techniques, such as atomistic molecular dynamics, to further resolve the underlying ensemble of structures. However, over the past 10–15 years, it has become evident that the classical, pairwise-additive force fields that have enjoyed a high degree of success for globular proteins have been somewhat limited in modeling IDP/IDR structures that agree with experiment. There has thus been a significant effort to rehabilitate these models to obtain better agreement with experiment, typically done by optimizing parameters in a piecewise fashion. In this work, we take a different approach by optimizing a set of force field parameters simultaneously, using machine learning to adapt force field parameters to experimental SAXS scattering profiles. We demonstrate our approach in modeling three biologically IDP ensembles based on experimental SAXS profiles and show that our optimization approach significantly improve force field parameters that generate ensembles in better agreement with experiment.

Proteasome Activation to Combat Proteotoxicity

Loss of proteome fidelity leads to the accumulation of non-native protein aggregates and oxidatively damaged species: hallmarks of an aged cell. These misfolded and aggregated species are often found, and suggested to be the culpable party, in numerous neurodegenerative diseases including Huntington's, Parkinson's, Amyotrophic Lateral Sclerosis (ALS), and Alzheimer's Diseases (AD). Many strategies for therapeutic intervention in proteotoxic pathologies have been put forth; one of the most promising is bolstering the efficacy of the proteasome to restore normal proteostasis. This strategy is ideal as monomeric precursors and oxidatively damaged proteins, so called "intrinsically disordered proteins" (IDPs), are targeted by the proteasome. This review will provide an overview of disorders in proteins, both intrinsic and acquired, with a focus on susceptibility to proteasomal degradation. We will then examine the proteasome with emphasis on newly published structural data and summarize current known small molecule proteasome activators.

bHLH-PAS Proteins: Their Structure and Intrinsic Disorder

The basic helix–loop–helix/Per-ARNT-SIM (bHLH–PAS) proteins are a class of transcriptional regulators, commonly occurring in living organisms and highly conserved among vertebrates and invertebrates. These proteins exhibit a relatively well-conserved domain structure: the bHLH domain located at the N-terminus, followed by PAS-A and PAS-B domains. In contrast, their C-terminal fragments present significant variability in their primary structure and are unique for individual proteins. C-termini were shown to be responsible for the specific modulation of protein action. In this review, we present the current state of knowledge, based on NMR and X-ray analysis, concerning the structural properties of bHLH–PAS proteins. It is worth noting that all determined structures comprise only selected domains (bHLH and/or PAS). At the same time, substantial parts of proteins, comprising their long C-termini, have not been structurally characterized to date. Interestingly, these regions appear to be intrinsically disordered (IDRs) and are still a challenge to research. We aim to emphasize the significance of IDRs for the flexibility and function of bHLH–PAS proteins. Finally, we propose modern NMR methods for the structural characterization of the IDRs of bHLH–PAS proteins.

Entropy and Information within Intrinsically Disordered Protein Regions

Bioinformatics and biophysical studies of intrinsically disordered proteins and regions (IDRs) note the high entropy at individual sequence positions and in conformations sampled in solution. This prevents application of the canonical sequence-structure-function paradigm to IDRs and motivates the development of new methods to extract information from IDR sequences. We argue that the information in IDR sequences cannot be fully revealed through positional conservation, which largely measures stable structural contacts and interaction motifs. Instead, considerations of evolutionary conservation of molecular features can reveal the full extent of information in IDRs. Experimental quantification of the large conformational entropy of IDRs is challenging but can be approximated through the extent of conformational sampling measured by a combination of NMR spectroscopy and lower-resolution structural biology techniques, which can be further interpreted with simulations. Conformational entropy and other biophysical features can be modulated by post-translational modifications that provide functional advantages to IDRs by tuning their energy landscapes and enabling a variety of functional interactions and modes of regulation. The diverse mosaic of functional states of IDRs and their conformational features within complexes demands novel metrics of information, which will reflect the complicated sequence-conformational ensemble-function relationship of IDRs.

NMR characterization of solvent accessibility and transient structure in intrinsically disordered proteins

In order to understand the conformational behavior of intrinsically disordered proteins (IDPs) and their biological interaction networks, the detection of residual structure and long-range interactions is required. However, the large number of degrees of conformational freedom of disordered proteins require the integration of extensive sets of experimental data, which are difficult to obtain. Here, we provide a straightforward approach for the detection of residual structure and long-range interactions in IDPs under near-native conditions using solvent paramagnetic relaxation enhancement (sPRE). Our data indicate that for the general case of an unfolded chain, with a local flexibility described by the overwhelming majority of available combinations, sPREs of non-exchangeable protons can be accurately predicted through an ensemble-based fragment approach. We show for the disordered protein α-synuclein and disordered regions of the proteins FOXO4 and p53 that deviation from random coil behavior can be interpreted in terms of intrinsic propensity to populate local structure in interaction sites of these proteins and to adopt transient long-range structure. The presented modification-free approach promises to be applicable to study conformational dynamics of IDPs and other dynamic biomolecules in an integrative approach.

Solvent-dependent segmental dynamics in intrinsically disordered proteins

Protein and water dynamics have a synergistic relationship, which is particularly important for intrinsically disordered proteins (IDPs), although the details of this coupling remain poorly understood. Here, we combine temperature-dependent molecular dynamics simulations using different water models with extensive nuclear magnetic resonance (NMR) relaxation to examine the importance of distinct modes of solvent and solute motion for the accurate reproduction of site-specific dynamics in IDPs. We find that water dynamics play a key role in motional processes internal to "segments" of IDPs, stretches of primary sequence that share dynamic properties and behave as discrete dynamic units. We identify a relationship between the time scales of intrasegment dynamics and the lifetime of hydrogen bonds in bulk water. Correct description of these motions is essential for accurate reproduction of protein relaxation. Our findings open important perspectives for understanding the role of hydration water on the behavior and function of IDPs in solution.

General Purpose Water Model Can Improve Atomistic Simulations of Intrinsically Disordered Proteins

Unconstrained atomistic simulations of intrinsically disordered proteins and peptides (IDP) remain a challenge: widely used, "general purpose" water models tend to favor overly compact structures relative to experiment. Here we have performed a total of 93 μs of unrestrained MD simulations to explore, in the context of IDPs, a recently developed "general-purpose" 4-point rigid water model OPC, which describes liquid state of water close to experiment. We demonstrate that OPC, together with a popular AMBER force field ff99SB, offers a noticeable improvement over TIP3P in producing more realistic structural ensembles of three common IDPs benchmarks: 55-residue apo N-terminal zinc-binding domain of HIV-1 integrase ("protein IN"), amyloid β-peptide (Aβ₄₂) (residues 1-42), and 26-reside H4 histone tail. As a negative control, computed folding profile of a regular globular miniprotein (CLN025) in OPC water is in appreciably better agreement with experiment than that obtained in TIP3P, which tends to overstabilize the compact native state relative to the extended conformations. We employed Aβ₄₂ peptide to investigate the possible influence of the solvent box size on simulation outcomes. We advocate a cautious approach for simulations of IDPs: we suggest that the solvent box size should be at least four times the radius of gyration of the random coil corresponding to the IDP. The computed free energy landscape of protein IN in OPC resembles a shallow "tub" - conformations with substantially different degrees of compactness that are within 2 k_B T of each other. Conformations with very different secondary structure content coexist within 1 k_B T of the global free energy minimum. States with higher free energy tend to have less secondary structure. Computed low helical content of the protein has virtually no correlation with its degree of compactness, which calls into question the possibility of using the helicity as a metric for assessing performance of water models for IDPs, when the helicity is low. Predicted radius of gyration ( R _g) of H4 histone tail in OPC water falls in-between that of a typical globular protein and a fully denatured protein of the same size; the predicted R _g is consistent with two independent predictions. In contrast, H4 tail in TIP3P water is as compact as the corresponding globular protein. The computed free energy landscape of H4 tail in OPC is relatively flat over a significant range of compactness, which, we argue, is consistent with its biological function as facilitator of internucleosome interactions.

Quo Vadis Biomolecular NMR Spectroscopy?

In-cell nuclear magnetic resonance (NMR) spectroscopy offers the possibility to study proteins and other biomolecules at atomic resolution directly in cells. As such, it provides compelling means to complement existing tools in cellular structural biology. Given the dominance of electron microscopy (EM)-based methods in current structure determination routines, I share my personal view about the role of biomolecular NMR spectroscopy in the aftermath of the revolution in resolution. Specifically, I focus on spin-off applications that in-cell NMR has helped to develop and how they may provide broader and more generally applicable routes for future NMR investigations. I discuss the use of ‘static’ and time-resolved solution NMR spectroscopy to detect post-translational protein modifications (PTMs) and to investigate structural consequences that occur in their response. I argue that available examples vindicate the need for collective and systematic efforts to determine post-translationally modified protein structures in the future. Furthermore, I explain my reasoning behind a Quinary Structure Assessment (QSA) initiative to interrogate cellular effects on protein dynamics and transient interactions present in physiological environments.

Concomitant disorder and high-affinity zinc binding in the human zinc- and iron-regulated transport protein 4 intracellular loop

The human zinc- and iron-regulated transport protein 4 (hZIP4) protein is the major plasma membrane protein responsible for the uptake of zinc in the body, and as such it plays a key role in cellular zinc homeostasis. hZIP4 plasma membrane levels are regulated through post-translational modification of its large, disordered, histidine-rich cytosolic loop (ICL2) in response to intracellular zinc concentrations. Here, structural characteristics of the isolated disordered loop region, both in the absence and presence of zinc, were investigated using nuclear magnetic resonance (NMR) spectroscopy. NMR chemical shifts, coupling constants and temperature coefficients of the apoprotein, are consistent with a random coil with minor propensities for transient polyproline Type II helices and β-strand in regions implicated in post-translational modifications. The ICL2 protein remains disordered upon zinc binding, which induces exchange broadening. Paramagnetic relaxation enhancement experiments reveal that the histidine-rich region in the apoprotein makes transient tertiary contacts with predicted post-translational modification sites. The residue-specific data presented here strengthen the relationship between hZIP4 post-translational modifications, which impact its role in cellular zinc homeostasis, and zinc sensing by the intracellular loop. Furthermore, the zinc sensing mechanism employed by the ICL2 protein demonstrates that high-affinity interactions can occur in the presence of conformational disorder.

Fast NMR method to probe solvent accessibility and disordered regions in proteins

Understanding protein structure and dynamics, which govern key cellular processes, is crucial for basic and applied research. Intrinsically disordered protein (IDP) regions display multifunctionality via alternative transient conformations, being key players in disease mechanisms. IDP regions are abundant, namely in small viruses, allowing a large number of functions out of a small proteome. The relation between protein function and structure is thus now seen from a different perspective: as IDP regions enable transient structural arrangements, each conformer can play different roles within the cell. However, as IDP regions are hard and time-consuming to study via classical techniques (optimized for globular proteins with unique conformations), new methods are required. Here, employing the dengue virus (DENV) capsid (C) protein and the immunoglobulin-binding domain of streptococcal protein G, we describe a straightforward NMR method to differentiate the solvent accessibility of single amino acid N-H groups in structured and IDP regions. We also gain insights into DENV C flexible fold region biological activity. The method, based on minimal pH changes, uses the well-established ¹H-¹⁵N HSQC pulse sequence and is easily implementable in current protein NMR routines. The data generated are simple to interpret, with this rapid approach being an useful first-choice IDPs characterization method.

Solution Ensemble of the C-Terminal Domain from the Transcription Factor Pdx1 Resembles an Excluded Volume Polymer

The pancreatic and duodenal homeobox 1 (Pdx1) is an essential pancreatic transcription factor. The C-terminal intrinsically disordered domain of Pdx1 (Pdx1-C) has a heavily biased amino acid composition; most notably, 18 of 83 residues are proline, including a hexaproline cluster near the middle of the chain. For these reasons, Pdx1-C is an attractive target for structure characterization, given the availability of suitable methods. To determine the solution ensembles of disordered proteins, we have developed a suite of ¹³C direct-detect NMR experiments that provide high spectral quality, even in the presence of strong proline enrichment. Here, we have extended our suite of NMR experiments to include four new pulse programs designed to record backbone residual dipolar couplings in a ¹³C,¹⁵N-CON detection format. Using our NMR strategy, in combination with small-angle X-ray scattering measurements and Monte Carlo simulations, we have determined that Pdx1-C is extended in solution, with a radius of gyration and internal scaling similar to that of an excluded volume polymer, and a subtle tendency toward a collapsed structure to the N-terminal side of the hexaproline sequence. This structure leaves Pdx1-C exposed for interactions with trans-regulatory co-factors that contribute with Pdx1 to transcription control in the cell.

NMR Characterization of Long-Range Contacts in Intrinsically Disordered Proteins from Paramagnetic Relaxation Enhancement in 13 C Direct-Detection Experiments

Intrinsically disordered proteins (IDPs) carry out many biological functions. They lack a stable 3D structure and are able to adopt many different conformations in dynamic equilibrium. The interplay between local dynamics and global rearrangements is key for their function. A widely used experimental NMR spectroscopy approach to study long-range contacts in IDPs exploits paramagnetic effects, and ¹ H detection experiments are generally used to determine paramagnetic relaxation enhancement (PRE) for amide protons. However, under physiological conditions, exchange broadening hampers the detection of solvent-exposed amide protons, which reduces the content of information available. Herein, we present an experimental approach based on direct carbon detection of PRE that provides improved resolution, reduced sensitivity to exchange broadening, and complementary information derived from the use of different starting polarization sources.

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.

AMBER and CHARMM Force Fields Inconsistently Portray the Microscopic Details of Phosphorylation

Phosphorylation of serine, threonine, and tyrosine is one of the most frequently occurring and crucial post-translational modifications of proteins often associated with important structural and functional changes. We investigated the direct effect of phosphorylation on the intrinsic conformational preferences of amino acids as a potential trigger of larger structural events. We conducted a comparative study of force fields on terminally capped amino acids (dipeptides) as the simplest model for phosphorylation. Our bias-exchange metadynamics simulations revealed that all model dipeptides sampled a great heterogeneity of ensembles affected by introduction of mono- and dianionic phosphate groups. However, the detected changes in populations of backbone conformers and side-chain rotamers did not reveal a strong discriminatory shift in preferences, as could be anticipated for the bulky, charged phosphate group. Furthermore, the AMBER and CHARMM force fields provided inconsistent populations of individual conformers as well as net structural trends upon phosphorylation. Detailed analysis of ensembles revealed competition between hydration and formation of internal hydrogen bonds involving amide hydrogens and the phosphate group. The observed difference in hydration free energy and potential for hydrogen bonding in individual force fields could be attributed to the different partial atomic charges used in each force field and, hence, the different parametrization strategies. Nevertheless, conformational propensities and net structural changes upon phosphorylation are difficult to extract from experimental measurements, and existing experimental data provide limited guidance for force field assessment and further development.

Towards a Stochastic Paradigm: From Fuzzy Ensembles to Cellular Functions

The deterministic sequence → structure → function relationship is not applicable to describe how proteins dynamically adapt to different cellular conditions. A stochastic model is required to capture functional promiscuity, redundant sequence motifs, dynamic interactions, or conformational heterogeneity, which facilitate the decision-making in regulatory processes, ranging from enzymes to membraneless cellular compartments. The fuzzy set theory offers a quantitative framework to address these problems. The fuzzy formalism allows the simultaneous involvement of proteins in multiple activities, the degree of which is given by the corresponding memberships. Adaptation is described via a fuzzy inference system, which relates heterogeneous conformational ensembles to different biological activities. Sequence redundancies (e.g., tandem motifs) can also be treated by fuzzy sets to characterize structural transitions affecting the heterogeneous interaction patterns (e.g., pathological fibrillization of stress granules). The proposed framework can provide quantitative protein models, under stochastic cellular conditions.

Structure and pro-toxic mechanism of the human Hsp90/PPIase/Tau complex

The molecular chaperone Hsp90 is critical for the maintenance of cellular homeostasis and represents a promising drug target. Despite increasing knowledge on the structure of Hsp90, the molecular basis of substrate recognition and pro-folding by Hsp90/co-chaperone complexes remains unknown. Here, we report the solution structures of human full-length Hsp90 in complex with the PPIase FKBP51, as well as the 280 kDa Hsp90/FKBP51 complex bound to the Alzheimer’s disease-related protein Tau. We reveal that the FKBP51/Hsp90 complex, which synergizes to promote toxic Tau oligomers in vivo, is highly dynamic and stabilizes the extended conformation of the Hsp90 dimer resulting in decreased Hsp90 ATPase activity. Within the ternary Hsp90/FKBP51/Tau complex, Hsp90 serves as a scaffold that traps the PPIase and nucleates multiple conformations of Tau’s proline-rich region next to the PPIase catalytic pocket in a phosphorylation-dependent manner. Our study defines a conceptual model for dynamic Hsp90/co-chaperone/client recognition.

The chaperone Hsp90 plays a key role in maintaining cellular homeostasis. Here the authors provide structural insights into substrate recognition and the pro-folding mechanism of Hsp90/co-chaperone complexes by studying the complex of Hsp90 with its co-chaperone FKBP51 and the substrate Tau bound Hsp90/FKBP51 ternary complex using a NMR based integrative approach.

Characterization of Dynamic IDP Complexes by NMR Spectroscopy

NMR spectroscopy has proven to be a key method for studying intrinsically disordered proteins (IDPs). Nonetheless, traditional NMR methods developed for solving structures of ordered protein complexes are insufficient for the full characterization of dynamic IDP complexes, where the energy landscape is broader and more rugged. Furthermore, due to their high sensitivity to environmental changes, NMR studies of IDP complexes must be conducted with extra care and the observed NMR parameters thoroughly evaluated to enable disentanglement of binding events from ensemble distribution changes. In this chapter, written for the non-NMR expert, we start out by outlining sample preparation for IDP complexes, guide through the recording and evaluation of diagnostic ¹H,¹⁵N-HSQC spectra, and delineate more sophisticated NMR strategies to follow for the particular type of complex. The most relevant experiments are then described in terms of aims, needs, pitfalls, analysis, and expected outcomes, with references to recent examples.

Fold or not to fold upon binding - does it really matter?

Protein interactions are usually determined by well-defined contact patterns. In this scenario, structuring of the interface is a prerequisite, which takes place prior or coupled to binding. Recent data, however, indicate plasticity of the templated folding pathway as well as considerable variations: polymorphism or dynamics in the bound-state. Conformational fluctuations in both cases are modulated by non-native, transient contacts, which complement suboptimal binding motifs to improve affinity. Here I discuss both templated folding and fuzzy binding mechanisms and propose a uniform scheme.

Elucidating binding mechanisms and dynamics of intrinsically disordered protein complexes using NMR spectroscopy

Advances in characterizing complexes of intrinsically disordered proteins (IDPs) have led to the discovery of a remarkably diverse interaction landscape that includes folding-upon-binding, highly dynamic complexes, multivalent interactions as well as regulatory switches controlled by post-translational modifications. Nuclear magnetic resonance (NMR) spectroscopy has in recent years made significant contributions to this field by describing the binding mechanisms and mapping conformational dynamics on multiple time scales. Importantly, this progress has been associated with specific methodological developments in NMR, for example in exchange techniques, allowing challenging biological systems to be studied at atomic resolution. In general, the level of dynamics observed in IDP complexes does not correlate with binding affinities, demonstrating the intricate relationship between conformational dynamics and IDP regulatory function.

A Method for Determining Structure Ensemble of Large Disordered Protein: Application to a Mechanosensing Protein

Structure characterization of intrinsically disordered proteins (IDPs) remains a key obstacle in understanding their functional mechanisms. Due to the highly dynamic feature of IDPs, structure ensembles instead of static unique structures are often derived from experimental data. Several state-of-the-art computational methods have been developed to select an optimal ensemble from a pregenerated structure pool, but they suffer from low efficiency for large IDPs. Here we present a matching pursuit genetic algorithm (MPGA) for structure ensemble determination, which takes advantages from both matching pursuit (MP) to reduce the search space and genetic algorithm (GA) to reduce the restriction on constraint types. The MPGA method is validated using a reference ensemble with predefined structures. In comparison with the conventional GA, MPGA takes much less computational time for large IDPs. The utility of the method is demonstrated by application to structure ensemble determination of a mechanosensing protein domain with 306 amino acids. The structure ensemble determined reveals that the N-terminal region 1-240 is more compact than the C-terminal region 240-306. The unique structural feature explains why only a small portion of YXXP tyrosine residues can be phosphorylated easily by kinases in the absence of extension force and why the phosphorylation is force-dependent.

An ultraweak interaction in the intrinsically disordered replication machinery is essential for measles virus function

Measles virus genome encapsidation is essential for viral replication and is controlled by the intrinsically disordered phosphoprotein (P) maintaining the nucleoprotein in a monomeric form (N) before nucleocapsid assembly. All paramyxoviruses harbor highly disordered amino-terminal domains (PNTD) that are hundreds of amino acids in length and whose function remains unknown. Using nuclear magnetic resonance (NMR) spectroscopy, we describe the structure and dynamics of the 90-kDa N0PNTD complex, comprising 450 disordered amino acids, at atomic resolution. NMR relaxation dispersion reveals the existence of an ultraweak N-interaction motif, hidden within the highly disordered PNTD, that allows PNTD to rapidly associate and dissociate from a specific site on N while tightly bound at the amino terminus, thereby hindering access to the surface of N. Mutation of this linear motif quenches the long-range dynamic coupling between the two interaction sites and completely abolishes viral transcription/replication in cell-based minigenome assays comprising integral viral replication machinery. This description transforms our understanding of intrinsic conformational disorder in paramyxoviral replication. The essential mechanism appears to be conserved across Paramyxoviridae, opening unique new perspectives for drug development against this family of pathogens.

Accurate Electron-Nucleus Distances from Paramagnetic Relaxation Enhancements

Measurements of paramagnetic relaxation enhancements (PREs) in ¹H NMR spectra are an important tool to obtain long-range distance information in proteins, but quantitative interpretation is easily compromised by nonspecific intermolecular PREs. Here we show that PREs generated by lanthanides with anisotropic magnetic susceptibilities offer a route to accurate calibration-free distance measurements. As these lanthanides change ¹H chemical shifts due to pseudocontact shifts, the relaxation rates in the paramagnetic and diamagnetic state can be measured with a single sample that simultaneously contains the protein labeled with a paramagnetic and a diamagnetic lanthanide ion. Nonspecific intermolecular PREs are thus automatically subtracted when calculating the PREs as the difference in nuclear relaxation rates between paramagnetic and diamagnetic protein. Although PREs from lanthanides with anisotropic magnetic susceptibilities are complicated by additional cross-correlation effects and residual dipolar couplings (RDCs) in the paramagnetic state, these effects can be controlled by the choice of lanthanide ion and experimental conditions. Using calbindin D_9k with erbium, we succeeded in measuring intramolecular PREs with unprecedented accuracy, resulting in distance predictions with a root-mean-square-deviation of <0.9 Å in the range 11-24 Å.

Crowding Stabilizes DMSO-Water Hydrogen-Bonding Interactions

Up to 40% of intracellular water is confined due to the dense packing of macromolecules, ions, and osmolytes. Despite the large body of work concerning the effect of additives on the biomolecular structure and stability, the role of crowding and heterogeneity in these interactions is not well understood. Here, infrared spectroscopy and molecular dynamics simulations are used to describe the mechanisms by which crowding modulates hydrogen bonding interactions between water and dimethyl sulfoxide (DMSO). Specifically, we use formamide and dimethylformamide (DMF) as molecular crowders and show that the S═O hydrogen bond populations in aqueous mixtures are increased by both amides. These additives increase the amount of water within the DMSO first solvation shell through two mechanisms: (a) directly stabilizing water-DMSO hydrogen bonds; (b) increasing water exposure by destabilizing DMSO-DMSO self-interactions. Further, we quantified the hydrogen bond enthalpies between the different components: DMSO-water (61 kJ/mol) > DMSO-formamide (32 kJ/mol) > water-water (23 kJ/mol) ≫ formamide-water (4.7 kJ/mol). Spectra of carbonyl stretching vibrations show that DMSO induces the dehydration of amides as a result of strong DMSO-water interactions, which has been suggested as the main mechanism of protein destabilization.

Quantum Mechanics Study on Hydrophilic and Hydrophobic Interactions in the Trivaline-Water System

With the aim to elucidate hydrophobic effects in the unfolded state of peptides, DFT-M062X computations on the Val₃H⁺· nH₂O ( n up to 22) clusters have been accomplished. As far as the main chain is concerned, four conformers with β-strand and/or polyproline type II conformations, PPII (indicated as β-β, β-PPII, PPII-β, and PPII-PPII), have been found by changing the ϕ and ψ angles. For bare peptide, the side chain (isopropyl) of each residue can independently take on three different orientations with negligible effects on energetics. The great isopropyl spatial separations in β-β and β-PPII conformers allow for the construction of synergic and extensive water-water and water-peptide H-bonding in the minimal hydration Val₃H⁺·22H₂O models without significant steric encumbrance. Conversely, due to the proximity of the isopropyl of the central residue with the other two, some restrictions in the water shell construction around the peptide become evident for the PPII-PPII conformer and the number of energetically accessible structures decreases. This is indicative of correlated motion involving isopropyls and backbone mediated by water molecules, the origin of the nearest neighbor effects. Comparing the thermodynamic data of Ala₃H⁺·22H₂O and Val₃H⁺·22H₂O, what emerges is that both hydration enthalpy and entropy drive the β-strand stability of the latter.

Accurately Predicting Disordered Regions of Proteins Using Rosetta ResidueDisorder Application

Although many proteins necessitate well-folded structures to properly instigate their biological functions, a large fraction of functioning proteins contain regions-known as intrinsically disordered protein regions-where stable structures are not likely to form. Notable functional roles of intrinsically disordered proteins are in transcriptional regulation, translation, and cellular signal transduction. Moreover, intrinsically disordered protein regions are highly abundant in many proteins associated with various human diseases, therefore these segments have become attractive drug targets for potential therapeutics. Over the past decades, numerous computational methods have been developed to accurately predict disordered regions of proteins. Here we introduce a user-friendly and reliable approach for the prediction of disordered protein regions using the structure prediction software Rosetta. Using 245 proteins from a benchmark data set (16 DisProt database proteins) and a test data set (229 proteins with NMR data), we use Rosetta to predict the global protein structures and then show that there is a statistically significant difference between Rosetta scores in disordered and ordered regions, with scores being less favorable in disordered regions. Furthermore, the difference in scores between ordered and disordered protein regions is sufficient to accurately identify disordered protein regions. As a result, our Rosetta ResidueDisorder method (benchmark data set prediction accuracy of 71.77% and independent test data set prediction accuracy of 65.37%) outperformed other established disorder prediction tools and did not exhibit a biased prediction toward either ordered or disordered regions. To facilitate usage, a Rosetta application has been developed for the Rosetta ResidueDisorder method.

Fuzziness in Protein Interactions-A Historical Perspective

The proposal that coupled folding to binding is not an obligatory mechanism for intrinsically disordered (ID) proteins was put forward 10 years ago. The notion of fuzziness implies that conformational heterogeneity can be maintained upon interactions of ID proteins, which has a functional impact either on regulated assembly or activity of the corresponding complexes. Here I review how the concept has evolved in the past decade, via increasing experimental data providing insights into the mechanisms, pathways and regulatory modes. The effects of structural diversity and transient contacts on protein assemblies have been collected and systematically analyzed (Fuzzy Complexes Database, http://protdyn-database.org). Fuzziness has also been exploited as a framework to decipher molecular organization of higher-order protein structures. Quantification of conformational heterogeneity opens exciting future perspectives for drug discovery from small molecule-ID protein interactions to supramolecular assemblies.

Hypothesis: structural heterogeneity of the unfolded proteins originating from the coupling of the local clusters and the long-range distance distribution

We propose a hypothesis that explains two apparently contradicting observations for the heterogeneity of the unfolded proteins. First, the line confocal method of the single-molecule Förster resonance energy transfer (sm-FRET) spectroscopy revealed that the unfolded proteins possess broad peaks in the FRET efficiency plot, implying the significant heterogeneity that lasts longer than milliseconds. Second, the fluorescence correlation method demonstrated that the unfolded proteins fluctuate in the time scale shorter than 100 ns. To formulate the hypothesis, we first summarize the recent consensus for the structure and dynamics of the unfolded proteins. We next discuss the conventional method of the sm-FRET spectroscopy and its limitations for the analysis of the unfolded proteins, followed by the advantages of the line confocal method that revealed the heterogeneity. Finally, we propose that the structural heterogeneity formed by the local clustering of hydrophobic residues modulates the distribution of the long-range distance between the labeled chromophores, resulting in the broadening of the peak in the FRET efficiency plot. A clustering of hydrophobic residues around the chromophore might further contribute to the broadening. The proposed clusters are important for the understanding of protein folding mechanism.

Conformational sampling of the intrinsically disordered dsRBD-1 domain from Arabidopsis thaliana DCL1

Partial folding and stability of DCL1-dsRBD1.

Selecting Conformational Ensembles Using Residual Electron and Anomalous Density (READ)

Heterogeneous and dynamic biomolecular complexes play a central role in many cellular processes but are poorly understood due to experimental challenges in characterizing their structural ensembles. To address these difficulties, we developed a hybrid methodology that combines X-ray crystallography with ensemble selections typically used in NMR studies to determine structural ensembles of heterogeneous biomolecular complexes. The method, termed READ, for residual electron and anomalous density, enables the visualization of heterogeneous conformational ensembles of complexes within crystals. Here we present a detailed protocol for performing the ensemble selections to construct READ ensembles. From a diverse pool of binding poses, a selection scheme is used to determine a subset of conformations that maximizes agreement with the X-ray data. Overall, READ is a general approach for obtaining a high-resolution view of dynamic protein-protein complexes.

Force field development and simulations of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) play important roles in many physiological processes such as signal transduction and transcriptional regulation. Computer simulations that are based on empirical force fields have been increasingly used to understand the biophysics of disordered proteins. In this review, we focus on recent improvement of protein force fields, including polarizable force fields, concerning their accuracy in modeling intrinsically disordered proteins. Some recent benchmarks and applications of these force fields are also overviewed.

Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental approaches for investigating the conformational behaviour of intrinsically disordered proteins (IDPs). IDPs represent a significant fraction of all proteomes, and, despite their importance for understanding fundamental biological processes, the molecular basis of their activity still remains largely unknown. The functional mechanisms exploited by IDPs in their interactions with other biomolecules are defined by their intrinsic dynamic modes and associated timescales, justifying the considerable interest over recent years in the development of technologies adapted to measure and describe this behaviour. NMR spin relaxation delivers information-rich, site-specific data reporting on conformational fluctuations occurring throughout the molecule. Here we review recent progress in the use of ¹⁵N relaxation to identify local backbone dynamics and long-range chain-like motions in unfolded proteins.

Prediction of nearest neighbor effects on backbone torsion angles and NMR scalar coupling constants in disordered proteins

Using fine-tuned hydrogen bonding criteria, a library of coiled peptide fragments has been generated from a large set of high-resolution protein X-ray structures. This library is shown to be an improved representation of ϕ/ψ torsion angles seen in intrinsically disordered proteins (IDPs). The ϕ/ψ torsion angle distribution of the library, on average, provides good agreement with experimentally observed chemical shifts and ³ J_HN-Hα coupling constants for a set of five disordered proteins. Inspection of the coil library confirms that nearest-neighbor effects significantly impact the ϕ/ψ distribution of residues in the coil state. Importantly, ³ J_HN-Hα coupling constants derived from the nearest-neighbor modulated backbone ϕ distribution in the coil library show improved agreement to experimental values, thereby providing a better way to predict ³ J_HN-Hα coupling constants for IDPs, and for identifying locations that deviate from fully random behavior.