The Conformational Landscape of an Intrinsically Disordered DNA-Binding Domain of a Transcription Regulator

Functional regulation of an intrinsically disordered protein via a conformationally excited state

A longstanding goal in the field of intrinsically disordered proteins (IDPs) is to characterize their structural heterogeneity and pinpoint the role of this heterogeneity in IDP function. Here, we use multinuclear chemical exchange saturation (CEST) nuclear magnetic resonance to determine the structure of a thermally accessible globally folded excited state in equilibrium with the intrinsically disordered native ensemble of a bacterial transcriptional regulator CytR. We further provide evidence from double resonance CEST experiments that the excited state, which structurally resembles the DNA-bound form of cytidine repressor (CytR), recognizes DNA by means of a "folding-before-binding" conformational selection pathway. The disorder-to-order regulatory switch in DNA recognition by natively disordered CytR therefore operates through a dynamical variant of the lock-and-key mechanism where the structurally complementary conformation is transiently accessed via thermal fluctuations.

Flexible Target Recognition of the Intrinsically Disordered DNA-Binding Domain of CytR Monitored by Single-Molecule Fluorescence Spectroscopy

The intrinsically disordered DNA-binding domain of cytidine repressor (CytR-DBD) folds in the presence of target DNA and regulates the expression of multiple genes in E. coli. To explore the conformational rearrangements in the unbound state and the target recognition mechanisms of CytR-DBD, we carried out single-molecule Förster resonance energy transfer (smFRET) measurements. The smFRET data of CytR-DBD in the absence of DNA show one major and one minor population assignable to an expanded unfolded state and a compact folded state, respectively. The population of the folded state increases and decreases upon titration with salt and denaturant, respectively, in an apparent two-state manner. The peak FRET efficiencies of both the unfolded and folded states change continuously with denaturant concentration, demonstrating the intrinsic flexibility of the DNA-binding domain and the deviation from a strict two-state transition. Remarkably, the CytR-DBD exhibits a compact structure when bound to both the specific and nonspecific DNA; however, the peak FRET efficiencies of the two structures are slightly but consistently different. The observed conformational heterogeneity highlights the potential structural changes required for CytR to bind variably spaced operator sequences.

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.

Quantification of Entropic Excluded Volume Effects Driving Crowding-Induced Collapse and Folding of a Disordered Protein

We investigate the conformational properties of the intrinsically disordered DNA-binding domain of CytR in the presence of the polymeric crowder polyethylene glycol (PEG). Integrating circular dichroism, nuclear magnetic resonance, and single-molecule Förster resonance energy transfer measurements, we demonstrate that disordered CytR populates a well-folded minor conformation in its native ensemble, while the unfolded ensemble collapses and folds with an increase in crowder density independent of the crowder size. Employing a statistical-mechanical model, the effective reduction in the accessible conformational space of a residue in the unfolded state is estimated to be 10% at 300 mg/mL PEG8000, relative to dilute conditions. The experimentally consistent PEG-temperature phase diagram thus constructed reveals that entropic effects can stabilize disordered CytR by 10 kJ mol^-1, driving the equilibrium toward folded conformations under physiological conditions. Our work highlights the malleable conformational landscape of CytR, the presence of a folded conformation in the disordered ensemble, and proposes a scaling relation for quantifying excluded volume effects on protein stability.

Structural-Energetic Basis for Coupling between Equilibrium Fluctuations and Phosphorylation in a Protein Native Ensemble

The functioning of proteins is intimately tied to their fluctuations in the native ensemble. The structural-energetic features that determine fluctuation amplitudes and hence the shape of the underlying landscape, which in turn determine the magnitude of the functional output, are often confounded by multiple variables. Here, we employ the FF1 domain from human p190A RhoGAP protein as a model system to uncover the molecular basis for phosphorylation of a buried tyrosine, which is crucial to the transcriptional activity associated with transcription factor TFII-I. Combining spectroscopy, calorimetry, statistical-mechanical modeling, molecular simulations, and in vitro phosphorylation assays, we show that the FF1 domain samples a diverse array of conformations in its native ensemble, some of which are phosphorylation-competent. Upon eliminating unfavorable charge-charge interactions through a single charge-reversal (K53E) or charge-neutralizing (K53Q) mutation, we observe proportionately lower phosphorylation extents due to the altered structural coupling, damped equilibrium fluctuations, and a more compact native ensemble. We thus establish a conformational selection mechanism for phosphorylation in the FF1 domain with K53 acting as a "gatekeeper", modulating the solvent exposure of the buried tyrosine. Our work demonstrates the role of unfavorable charge-charge interactions in governing functional events through the modulation of native ensemble characteristics, a feature that could be prevalent in ordered protein domains.

Non-specific DNA-driven quinary interactions promote structural transitions in proteins

The nature and distribution of charged residues on the surface of proteins play a vital role in determining the binding affinity, selectivity and kinetics of association to ligands. When it comes to DNA-binding domains (DBDs), these functional features manifest as anisotropic distribution of positively charged residues on the protein surface driven by the requirement to bind DNA, a highly negatively charged polymer. In this work, we compare the thermodynamic behavior of nine different proteins belonging to three families - LacR, engrailed and Brk - some of which are disordered in solution in the absence of DNA. Combining detailed electrostatic calculations and statistical mechanical modeling of folding landscapes at different distances and relative orientations with respect to DNA, we show that non-specific electrostatic interactions between the protein and DNA can promote structural transitions in DBDs. Such quinary interactions that are strictly agnostic to the DNA sequence induce varied behaviors including folding of disordered domains, partial unfolding of ordered proteins and (de-)population of intermediate states. Our work highlights that the folding landscape of proteins can be tuned as a function of distance from DNA and hints at possible reasons for DBDs exhibiting complex kinetic-thermodynamic behaviors in the absence of DNA.

Thermodynamics and folding landscapes of large proteins from a statistical mechanical model

Statistical mechanical models that afford an intermediate resolution between macroscopic chemical models and all-atom simulations have been successful in capturing folding behaviors of many small single-domain proteins. However, the applicability of one such successful approach, the Wako-Saitô-Muñoz-Eaton (WSME) model, is limited by the size of the protein as the number of conformations grows exponentially with protein length. In this work, we surmount this size limitation by introducing a novel approximation that treats stretches of 3 or 4 residues as blocks, thus reducing the phase space by nearly three orders of magnitude. The performance of the 'bWSME' model is validated by comparing the predictions for a globular enzyme (RNase H) and a repeat protein (IκBα), against experimental observables and the model without block approximation. Finally, as a proof of concept, we predict the free-energy surface of the 370-residue, multi-domain maltose binding protein and identify an intermediate in good agreement with single-molecule force-spectroscopy measurements. The bWSME model can thus be employed as a quantitative predictive tool to explore the conformational landscapes of large proteins, extract the structural features of putative intermediates, identify parallel folding paths, and thus aid in the interpretation of both ensemble and single-molecule experiments.

Controlling Structure and Dimensions of a Disordered Protein via Mutations

The dimensions of intrinsically disordered proteins (IDPs) are sensitive to small energetic-entropic differences between intramolecular and protein–solvent interactions. This is commonly observed on modulating solvent composition and temperature. However, the inherently heterogeneous conformational landscape of IDPs is also expected to be influenced by mutations that can (de)stabilize pockets of local and even global structure, native and non-native, and hence the average dimensions. Here, we show experimental evidence for the remarkably tunable landscape of IDPs by employing the DNA-binding domain of CytR, a high-sequence-complexity IDP, as a model system. CytR exhibits a range of structure and compactness upon introducing specific mutations that modulate microscopic terms, including main-chain entropy, hydrophobicity, and electrostatics. The degree of secondary structure, as monitored by far-UV circular dichroism (CD), is strongly correlated to average ensemble dimensions for 14 different mutants of CytR and is consistent with the Uversky–Fink relation. Our experiments highlight how average ensemble dimensions can be controlled via mutations even in the disordered regime, the prevalence of non-native interactions and provide testable controls for molecular simulations.

Protein plasticity driven by disorder and collapse governs the heterogeneous binding of CytR to DNA

The amplitude of thermodynamic fluctuations in biological macromolecules determines their conformational behavior, dimensions, nature of phase transitions and effectively their specificity and affinity, thus contributing to fine-tuned molecular recognition. Unique among large-scale conformational changes in proteins are temperature-induced collapse transitions in intrinsically disordered proteins (IDPs). Here, we show that CytR DNA-binding domain, an IDP that folds on binding DNA, undergoes a coil-to-globule transition with temperature in the absence of DNA while exhibiting energetically decoupled local and global structural rearrangements, and maximal thermodynamic fluctuations at the optimal bacterial growth temperature. The collapse is shown to be a continuous transition through a combination of statistical-mechanical modeling and all-atom implicit solvent simulations. Surprisingly, CytR binds single-site cognate DNA with negative cooperativity, described by Hill coefficients less than one, resulting in a graded binding response. We show that heterogeneity arising from varying binding-competent CytR conformations or orientations at the single-molecular level contributes to negative binding cooperativity at the level of bulk measurements due to the conflicting requirements of collapse transition, large fluctuations and folding-upon-binding. Our work reports strong evidence for functionally driven thermodynamic fluctuations in determining the extent of collapse and disorder with implications in protein search efficiency of target DNA sites and regulation.

Tunable order-disorder continuum in protein-DNA interactions

DNA-binding protein domains (DBDs) sample diverse conformations in equilibrium facilitating the search and recognition of specific sites on DNA over millions of energetically degenerate competing sites. We hypothesize that DBDs have co-evolved to sense and exploit the strong electric potential from the array of negatively charged phosphate groups on DNA. We test our hypothesis by employing the intrinsically disordered DBD of cytidine repressor (CytR) as a model system. CytR displays a graded increase in structure, stability and folding rate on increasing the osmolarity of the solution that mimics the non-specific screening by DNA phosphates. Electrostatic calculations and an Ising-like statistical mechanical model predict that CytR exhibits features of an electric potential sensor modulating its dimensions and landscape in a unique distance-dependent manner, while DNA plays the role of a non-specific macromolecular chaperone. Accordingly, CytR binds its natural half-site faster than the diffusion-controlled limit and even random DNA conforming to an electrostatic-steering binding mechanism. Our work unravels for the first time the synergistic features of a natural electrostatic potential sensor, a novel binding mechanism driven by electrostatic frustration and disorder, and the role of DNA in promoting distance-dependent protein structural transitions critical for switching between specific and non-specific DNA-binding modes.

Engineering Order and Cooperativity in a Disordered Protein

Structural disorder in proteins arises from a complex interplay between weak hydrophobicity and unfavorable electrostatic interactions. The extent to which the hydrophobic effect contributes to the unique and compact native state of proteins is, however, confounded by large compensation between multiple entropic and energetic terms. Here we show that protein structural order and cooperativity arise as emergent properties upon hydrophobic substitutions in a disordered system with non-intuitive effects on folding and function. Aided by sequence-structure analysis, equilibrium, and kinetic spectroscopic studies, we engineer two hydrophobic mutations in the disordered DNA-binding domain of CytR that act synergistically, but not in isolation, to promote structure, compactness, and stability. The double mutant, with properties of a fully ordered domain, exhibits weak cooperativity with a complex and rugged conformational landscape. The mutant, however, binds cognate DNA with an affinity only marginally higher than that of the wild type, though nontrivial differences are observed in the binding to noncognate DNA. Our work provides direct experimental evidence of the dominant role of non-additive hydrophobic effects in shaping the molecular evolution of order in disordered proteins and vice versa, which could be generalized to even folded proteins with implications for protein design and functional manipulation.

Molecular dynamics ensemble refinement of the heterogeneous native state of NCBD using chemical shifts and NOEs

Many proteins display complex dynamical properties that are often intimately linked to their biological functions. As the native state of a protein is best described as an ensemble of conformations, it is important to be able to generate models of native state ensembles with high accuracy. Due to limitations in sampling efficiency and force field accuracy it is, however, challenging to obtain accurate ensembles of protein conformations by the use of molecular simulations alone. Here we show that dynamic ensemble refinement, which combines an accurate atomistic force field with commonly available nuclear magnetic resonance (NMR) chemical shifts and NOEs, can provide a detailed and accurate description of the conformational ensemble of the native state of a highly dynamic protein. As both NOEs and chemical shifts are averaged on timescales up to milliseconds, the resulting ensembles reflect the structural heterogeneity that goes beyond that probed, e.g., by NMR relaxation order parameters. We selected the small protein domain NCBD as object of our study since this protein, which has been characterized experimentally in substantial detail, displays a rich and complex dynamical behaviour. In particular, the protein has been described as having a molten-globule like structure, but with a relatively rigid core. Our approach allowed us to describe the conformational dynamics of NCBD in solution, and to probe the structural heterogeneity resulting from both short- and long-timescale dynamics by the calculation of order parameters on different time scales. These results illustrate the usefulness of our approach since they show that NCBD is rather rigid on the nanosecond timescale, but interconverts within a broader ensemble on longer timescales, thus enabling the derivation of a coherent set of conclusions from various NMR experiments on this protein, which could otherwise appear in contradiction with each other.

Entropic Control of an Excited Folded-Like Conformation in a Disordered Protein Ensemble

Many intrinsically disordered proteins switch between unfolded and folded-like forms in the presence of their binding partner. The possibility of a pre-equilibrium between the two macrostates is challenging to discern given the complex conformational landscape. Here, we show that CytR, a disordered DNA-binding domain, samples a folded-like excited state in its native ensemble through equilibrium multi-probe spectroscopy, kinetics and an Ising-like statistical mechanical model. The population of the excited state increases upon stabilization of the native ensemble with an osmolyte, while decreasing with increasing temperatures. A conserved proline residue, the mutation of which weakens the binding affinity to the target promoter, is found to uniquely control the population of the minor excited state. Semi-quantitative statistical mechanical modeling reveals that the conformational diffusion coefficient of disordered CytR is an order of magnitude slower than the estimates from folded domains. The osmolyte and proline mutation smoothen and roughen up the landscape, respectively, apart from modulation of populations. Our work uncovers general strategies to probe for excited structured states in disordered ensembles, and to measure and modulate the roughness of the disordered landscapes, inter-conversion rates of species and their populations.

Energetic and topological determinants of a phosphorylation-induced disorder-to-order protein conformational switch

We show that the phosphorylation of 4E-BP2 acts as a triggering event to shape its folding-function landscape that is delicately balanced between conflicting favorable energetics and intrinsically unfavorable topological connectivity. We further provide first evidence that the fitness landscapes of proteins at the threshold of disorder can differ considerably from ordered domains.

Imprints of function on the folding landscape: functional role for an intermediate in a conserved eukaryotic binding protein

In the computational characterization of single domain protein folding, the effective free energies of numerous microstates are projected onto few collective degrees of freedom that in turn serve as well-defined reaction coordinates. In this regard, one-dimensional (1D) free energy profiles are widely used mainly for their simplicity. Since folding and functional landscapes are interlinked, how well can these reduced representations capture the structural and dynamic features of functional states while being simultaneously consistent with experimental observables? We investigate this issue by characterizing the folding of the four-helix bundle bovine acyl-CoA binding protein (bACBP), which exhibits complex equilibrium and kinetic behaviours, employing an Ising-like statistical mechanical model and molecular simulations. We show that the features of the 1D free energy profile are sufficient to quantitatively reproduce multiple experimental observations including millisecond chevron-like kinetics and temperature dependence, a microsecond fast phase, barrier heights, unfolded state movements, the intermediate structure and average ϕ-values. Importantly, we find that the structural features of the native-like intermediate (partial disorder in helix 1) are intricately linked to a unique interplay between packing and electrostatics in this domain. By comparison with available experimental data, we propose that this intermediate determines the promiscuous functional behaviour of bACBP that exhibits broad substrate specificity. Our results present evidence to the possibility of employing the statistical mechanical model and the resulting 1D free energy profile to not just understand folding mechanisms but to even extract features of functionally relevant states and their energetic origins.

A human transcription factor in search mode

Transcription factors (TF) can change shape to bind and recognize DNA, shifting the energy landscape from a weak binding, rapid search mode to a higher affinity recognition mode. However, the mechanism(s) driving this conformational change remains unresolved and in most cases high-resolution structures of the non-specific complexes are unavailable. Here, we investigate the conformational switch of the human mitochondrial transcription termination factor MTERF1, which has a modular, superhelical topology complementary to DNA. Our goal was to characterize the details of the non-specific search mode to complement the crystal structure of the specific binding complex, providing a basis for understanding the recognition mechanism. In the specific complex, MTERF1 binds a significantly distorted and unwound DNA structure, exhibiting a protein conformation incompatible with binding to B-form DNA. In contrast, our simulations of apo MTERF1 revealed significant flexibility, sampling structures with superhelical pitch and radius complementary to the major groove of B-DNA. Docking these structures to B-DNA followed by unrestrained MD simulations led to a stable complex in which MTERF1 was observed to undergo spontaneous diffusion on the DNA. Overall, the data support an MTERF1-DNA binding and recognition mechanism driven by intrinsic dynamics of the MTERF1 superhelical topology.

Demonstration of a folding after binding mechanism in the recognition between the measles virus NTAIL and X domains

In the past decade, a wealth of experimental data has demonstrated that a large fraction of proteins, while functional, are intrinsically disordered at physiological conditions. Many intrinsically disordered proteins (IDPs) undergo a disorder-to-order transition upon binding to their biological targets, a phenomenon known as induced folding. Induced folding may occur through two extreme mechanisms, namely conformational selection and folding after binding. Although the pre-existence of ordered structures in IDPs is a prerequisite for conformational selection, it does not necessarily commit to this latter mechanism, and kinetic studies are needed to discriminate between the two possible scenarios. So far, relatively few studies have addressed this issue from an experimental perspective. Here, we analyze the interaction kinetics between the intrinsically disordered C-terminal domain of the measles virus nucleoprotein (NTAIL) and the X domain (XD) of the viral phosphoprotein. Data reveal that NTAIL recognizes XD by first forming a weak encounter complex in a disordered conformation, which is subsequently locked-in by a folding step; i.e., binding precedes folding. The implications of our kinetic results, in the context of previously reported equilibrium data, are discussed. These results contribute to enhancing our understanding of the molecular mechanisms by which IDPs recognize their partners and represent a paradigmatic example of the need of kinetic methods to discriminate between reaction mechanisms.

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.

Evidence for the sequential folding mechanism in RNase H from an ensemble-based model

The number of distinct protein folding pathways starting from an unfolded ensemble, and hence, the folding mechanism is an intricate function of protein size, sequence complexity, and stability conditions. This has traditionally been a contentious issue particularly because of the ensemble nature of conventional experiments that can mask the complexity of the underlying folding landscape. Recent hydrogen-exchange experiments combined with mass spectrometry (HX-MS) reveal that the folding of RNase H proceeds in a hierarchical fashion with distinct intermediates and following a single discrete path. In our current work, we provide computational evidence for this unique folding mechanism by employing a structure-based statistical mechanical model. Upon calibrating the energetic terms of the model with equilibrium measurements, we predict multiple intermediate states in the folding of RNase H that closely resemble experimental observations. Remarkably, a simplified landscape representation adequately captures the folding complexity and predicts the possibility of a well-defined sequence of folding events. We supplement the statistical model study with both explicit solvent molecular simulations of the helical units and electrostatic calculations to provide structural and energetic insights into the early and late stages of RNase H folding that hint at the frustrated nature of its folding landscape.