From the Cover: Charge interactions can dominate the dimensions of intrinsically disordered proteins

Kinetic consequences of native state optimization of surface-exposed electrostatic interactions in the Fyn SH3 domain

Optimization of surface exposed charge-charge interactions in the native state has emerged as an effective means to enhance protein stability; but the effect of electrostatic interactions on the kinetics of protein folding is not well understood. To investigate the kinetic consequences of surface charge optimization, we characterized the folding kinetics of a Fyn SH3 domain variant containing five amino acid substitutions that was computationally designed to optimize surface charge-charge interactions. Our results demonstrate that this optimized Fyn SH3 domain is stabilized primarily through an eight-fold acceleration in the folding rate. Analyses of the constituent single amino acid substitutions indicate that the effects of optimization of charge-charge interactions on folding rate are additive. This is in contrast to the trend seen in folded state stability, and suggests that electrostatic interactions are less specific in the transition state compared to the folded state. Simulations of the transition state using a coarse-grained chain model show that native electrostatic contacts are weakly formed, thereby making the transition state conducive to nonspecific, or even nonnative, electrostatic interactions. Because folding from the unfolded state to the folding transition state for small proteins is accompanied by an increase in charge density, nonspecific electrostatic interactions, that is, generic charge density effects can have a significant contribution to the kinetics of protein folding. Thus, the interpretation of the effects of amino acid substitutions at surface charged positions may be complicated and consideration of only native-state interactions may fail to provide an adequate picture.

Length-dependent compaction of intrinsically disordered proteins

This work investigates the effect of chain length on the degree of compaction of intrinsically disordered proteins (IDPs). The three main IDP types, native coil (NC), pre-molten globule (PMG) and molten globule (MG), are compared by means of a compaction index (CI) normalized for chain length. The results point out a strong variability of compactness as a function of chain length within each group, with larger proteins populating more compact states. While qualitative sequence features are responsible for the main differences among groups, chain length seems to have an unspecific effect modulating the extent of compaction within each group. The results are consistent with a cooperative character of the weak interactions responsible for chain collapse.

How, when and why proteins collapse: the relation to folding

Unfolded proteins under strongly denaturing conditions are highly expanded. However, when the conditions are more close to native, an unfolded protein may collapse to a compact globular structure distinct from the folded state. This transition is akin to the coil-globule transition of homopolymers. Single-molecule FRET experiments have been particularly conducive in revealing the collapsed state under conditions of coexistence with the folded state. The collapse can be even more readily observed in natively unfolded proteins. Time-resolved studies, using FRET and small-angle scattering, have shown that the collapse transition is a very fast event, probably occurring on the submicrosecond time scale. The forces driving collapse are likely to involve both hydrophobic and backbone interactions. The loss of configurational entropy during collapse makes the unfolded state less stable compared to the folded state, thus facilitating folding.

Fractal dimension of an intrinsically disordered protein: small-angle X-ray scattering and computational study of the bacteriophage λ N protein

Small-angle X-ray scattering (SAXS) was used to characterize the bacteriophage λ N protein, a 107 residue intrinsically disordered protein (IDP) that functions as a transcriptional antitermination factor. The SAXS data were used to estimate both the average radius of gyration and the fractal dimension, a measure of the protein's internal scaling properties, under a variety of solution conditions. In the absence of denaturants, the radius of gyration was 38 ± 3.5 Å and the fractal dimension was 1.76 ± 0.05, slightly larger than the value predicted for a well-solvated polymer with excluded volume (1.7). Neither the radius of gyration nor the fractal dimension changed significantly on the addition of urea, further indicating that the protein is extensively unfolded and well solvated in the absence of denaturant. The addition of NaCl or D(2) O was found to promote aggregation, but did not appear to affect the properties of the monomeric form. The experimental SAXS profiles were also compared with those predicted by a computational model for a random-coil polypeptide, with an adjustable solvation energy term. The experimental data were well fit to the model with the solvation energy close to zero. These results indicate that the λ N protein is among the more expanded members of the broad class of IDPs, most likely because of its high content of charged residues and a large net charge (+15 at neutral pH). The expanded nature of the conformational ensemble may play a role in facilitating the interactions of the protein with other components of the dynamic transcriptional complex.

Modulation of an IDP binding mechanism and rates by helix propensity and non-native interactions: association of HIF1α with CBP

Intrinsically disordered proteins that acquire their three dimensional structures only upon binding to their targets are very important in cellular signal regulation. While experimental studies have been made on the structures of both bound (structured) and unbound (disordered) states, less is known about the actual folding-binding transition. Coarse grained simulations using native-centric (i.e. Gō) potentials have been particularly useful in addressing this problem, given the large search space for IDP binding, but have well-known deficiencies in reproducing the unfolded state structure and dynamics. Here, we investigate the interaction of HIF1α with CBP using a hierarchy of coarse-grained models, in each case matching the binding affinity at 300 K to the experimental value. Starting from a pure Gō-like model based on the native structure of the complex we go on to consider a more realistic model of helix propensity in the HIF1α, and finally the effect of non-native interactions between binding partners. We find structural disorder (i.e."fuzziness") in the bound state of HIF1α in all models which is supported by the results of atomistic simulations. Correcting the over-stabilized helices in the unbound state gives rise to a more cooperative folding-binding transition (destabilizing partially bound intermediates). Adding non-native contacts lowers the free energy barrier for binding to an almost barrierless scenario, leading to higher binding/unbinding rates relative to the other models, in better agreement with the near diffusion-limited binding rates measured experimentally. Transition state structures for the three models are highly disordered, supporting a fly-casting mechanism for binding.

Compaction properties of an intrinsically disordered protein: Sic1 and its kinase-inhibitor domain

IDPs in their unbound state can transiently acquire secondary and tertiary structure. Describing such intrinsic structure is important to understand the transition between free and bound state, leading to supramolecular complexes with physiological interactors. IDP structure is highly dynamic and, therefore, difficult to study by conventional techniques. This work focuses on conformational analysis of the KID fragment of the Sic1 protein, an IDP with a key regulatory role in the cell-cycle of Saccharomyces cerevisiae. FT-IR spectroscopy, ESI-MS, and IM measurements are used to capture dynamic and short-lived conformational states, probing both secondary and tertiary protein structure. The results indicate that the isolated Sic1 KID retains dynamic helical structure and populates collapsed states of different compactness. A metastable, highly compact species is detected. Comparison between the fragment and the full-length protein suggests that chain length is crucial to the stabilization of compact states of this IDP. The two proteins are compared by a length-independent compaction index.

Beyond the random coil: stochastic conformational switching in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) participate in critical cellular functions that exploit the flexibility and rapid conformational fluctuations of their native state. Limited information about the native state of IDPs can be gained by the averaging over many heterogeneous molecules that is unavoidable in ensemble approaches. We used single molecule fluorescence to characterize native state conformational dynamics in five synaptic proteins confirmed to be disordered by other techniques. For three of the proteins, SNAP-25, synaptobrevin and complexin, their conformational dynamics could be described with a simple semiflexible polymer model. Surprisingly, two proteins, neuroligin and the NMDAR-2B glutamate receptor, were observed to stochastically switch among distinct conformational states despite the fact that they appeared intrinsically disordered by other measures. The hop-like intramolecular diffusion found in these proteins is suggested to define a class of functionality previously unrecognized for IDPs.

Expanding the proteome: disordered and alternatively folded proteins

Proteins provide much of the scaffolding for life, as well as undertaking a variety of essential catalytic reactions. These characteristic functions have led us to presuppose that proteins are in general functional only when well structured and correctly folded. As we begin to explore the repertoire of possible protein sequences inherent in the human and other genomes, two stark facts that belie this supposition become clear: firstly, the number of apparent open reading frames in the human genome is significantly smaller than appears to be necessary to code for all of the diverse proteins in higher organisms, and secondly that a significant proportion of the protein sequences that would be coded by the genome would not be expected to form stable three-dimensional (3D) structures. Clearly the genome must include coding for a multitude of alternative forms of proteins, some of which may be partly or fully disordered or incompletely structured in their functional states. At the same time as this likelihood was recognized, experimental studies also began to uncover examples of important protein molecules and domains that were incompletely structured or completely disordered in solution, yet remained perfectly functional. In the ensuing years, we have seen an explosion of experimental and genome-annotation studies that have mapped the extent of the intrinsic disorder phenomenon and explored the possible biological rationales for its widespread occurrence. Answers to the question 'why would a particular domain need to be unstructured?' are as varied as the systems where such domains are found. This review provides a survey of recent new directions in this field, and includes an evaluation of the role not only of intrinsically disordered proteins but also of partially structured and highly dynamic members of the disorder-order continuum.

Effect of Src kinase phosphorylation on disordered C-terminal domain of N-methyl-D-aspartic acid (NMDA) receptor subunit GluN2B protein

NMDA receptors are ligand-gated ion channels with a regulatory intracellular C-terminal domain (CTD). In GluN2B, the CTD is the largest domain in the protein but is intrinsically disordered. The GluN2B subunit is the major tyrosine-phosphorylated protein in synapses. Src kinase phosphorylates the GluN2B CTD, but it is unknown how this affects channel activity. In disordered proteins, phosphorylation can tip the balance between order and disorder. Transitions can occur in both directions, so it is not currently possible to predict the effects of phosphorylation. We used single molecule fluorescence to characterize the effects of Src phosphorylation on GluN2B. Scanning fluorescent labeling sites throughout the domain showed no positional dependence of the energy transfer. Instead, efficiency only scaled with the separation between labeling sites suggestive of a relatively featureless conformational energy landscape. Src phosphorylation led to a general expansion of the polypeptide, which would result in greater exposure of known protein-binding sites and increase the physical separation between contiguous sites. Phosphorylation makes the CTD more like a random coil leaving open the question of how Src exerts its effects on the NMDA receptor.

Organic matrix-related mineralization of sea urchin spicules, spines, test and teeth

The camarodont echinoderms have five distinct mineralized skeletal elements: embryonic spicules, mature test, spines, lantern stereom and teeth. The spicules are transient structural elements whereas the spines, and test plates are permanent. The teeth grow continuously. The mineral is a high magnesium calcite, but the magnesium content is different in each type of skeletal element, varying from 5 to 40 mole% Mg. The organic matrix creates the spaces and environments for crystal initiation and growth. The detailed mechanisms of crystal regulation are not known, but acidic and phosphorylated matrix proteins may be of special importance. Biochemical studies, sequencing of the complete genome, and high-throughput proteomic analysis have not yet provided insight into the mechanisms of crystallization, calcite composition, and orientation applicable to all skeletal elements. The embryonic spicules are not representative of the mature skeletal elements. The next phase of research will have to focus on the specific localization of the proteins and individual biochemistries of each system with regard to mineral content and placement.

Dual role of protein phosphorylation in DNA activator/coactivator binding

Binding free energies are calculated for the phosphorylated and unphosphorylated complexes between the kinase inducible domain (KID) of the DNA transcriptional activator cAMP response element binding (CREB) protein and the KIX domain of its coactivator, CREB-binding protein (CBP). To our knowledge, this is the first application of a method based on a potential of mean force (PMF) with restraining potentials to compute the binding free energy of protein-protein complexes. The KID:KIX complexes are chosen here because of their biological relevance to the DNA transcription process and their relatively small size (81 residues for the KIX domain of CBP, and 28 residues for KID). The results for pKID:KIX and KID:KIX are -9.55 and -4.96 kcal/mol, respectively, in good agreement with experimental estimates (-8.8 and -5.8 kcal/mol, respectively). A comparison between specific contributions to protein-protein binding for the phosphorylated and unphosphorylated complexes reveals a dual role for the phosphorylation of KID at Ser-133 in effecting a more favorable free energy of the bound system: 1), stabilization of the unbound conformation of phosphorylated KID due to favorable intramolecular interactions of the phosphate group of Ser-133 with the charged groups of an arginine-rich region spanning both α-helices, which lowers the configurational entropy; and 2), more favorable intermolecular electrostatic interactions between pSer-133 and Arg-131 of KID, and Lys-662, Tyr-658, and Glu-666 of KIX. Charge reduction through ligand phosphorylation emerges as a possible mechanism for controlling the unbound state conformation of KID and, ultimately, gene expression. This work also demonstrates that the PMF-based method with restraining potentials provides an added benefit in that important elements of the binding pathway are evidenced. Furthermore, the practicality of the PMF-based method for larger systems is validated by agreement with experiment. In addition, we provide a somewhat differently structured exposition of the PMF-based method with restraining potentials and outline its generalization to systems in which both protein and ligand may adopt unbound conformations that are different from those of the bound state.

Minireview: structural insights into early folding events using continuous-flow time-resolved small-angle X-ray scattering

Small-angle X-ray scattering (SAXS) is a powerful method for obtaining quantitative structural information on the size and shape of proteins, and it is increasingly used in kinetic studies of folding and association reactions. In this minireview, we discuss recent developments in using SAXS to obtain structural information on the unfolded ensemble and early folding intermediates of proteins using continuous-flow mixing devices. Interfacing of these micromachined devices to SAXS beamlines has allowed access to the microsecond time regime. The experimental constraints in implementation of turbulence and laminar flow-based mixers with SAXS detection and a comparison of the two approaches are presented. Current improvements and future prospects of microsecond time-resolved SAXS and the synergy with ab initio structure prediction and molecular dynamics simulations are discussed.

The N(0)-binding region of the vesicular stomatitis virus phosphoprotein is globally disordered but contains transient α-helices

The phosphoprotein (P) of vesicular stomatitis virus (VSV) interacts with nascent nucleoprotein (N), forming the N(0)-P complex that is indispensable for the correct encapsidation of newly synthesized viral RNA genome. In this complex, the N-terminal region (P(NTR)) of P prevents N from binding to cellular RNA and keeps it available for encapsidating viral RNA genomes. Here, using nuclear magnetic resonance (NMR) spectroscopy and small-angle X-ray scattering (SAXS), we show that an isolated peptide corresponding to the 60 first N-terminal residues of VSV P (P(60)) and encompassing P(NTR) has overall molecular dimensions and a dynamic behavior characteristic of a disordered protein but transiently populates conformers containing α-helices. The modeling of P(60) as a conformational ensemble by the ensemble optimization method using SAXS data correctly reproduces the α-helical content detected by NMR spectroscopy and suggests the coexistence of subensembles of different compactness. The populations and overall dimensions of these subensembles are affected by the addition of stabilizing (1M trimethylamine-N-oxide) or destabilizing (6M guanidinium chloride) cosolvents. Our results are interpreted in the context of a scenario whereby VSV P(NTR) constitutes a molecular recognition element undergoing a disorder-to-order transition upon binding to its partner when forming the N(0)-P complex.

Protein folding at single-molecule resolution

The protein folding reaction carries great significance for cellular function and hence continues to be the research focus of a large interdisciplinary protein science community. Single-molecule methods are providing new and powerful tools for dissecting the mechanisms of this complex process by virtue of their ability to provide views of protein structure and dynamics without associated ensemble averaging. This review briefly introduces common FRET and force methods, and then explores several areas of protein folding where single-molecule experiments have yielded insights. These include exciting new information about folding landscapes, dynamics, intermediates, unfolded ensembles, intrinsically disordered proteins, assisted folding and biomechanical unfolding. Emerging and future work is expected to include advances in single-molecule techniques aimed at such investigations, and increasing work on more complex systems from both the physics and biology standpoints, including folding and dynamics of systems of interacting proteins and of proteins in cells and organisms. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Determining serpin conformational distributions with single molecule fluorescence

Conformational plasticity is key to inhibitory serpin function, and this plasticity gives serpins relatively easy access to alternative, dysfunctional conformations. Thus, a given serpin population may contain both functional and dysfunctional proteins. Single molecule fluorescence (SMF), with its ability to interrogate one fluorescently labeled protein at a time, is a powerful method for elucidating conformational distributions and monitoring how these distributions change over time. SMF and related methods have been particularly valuable for characterizing serpin polymerization. Fluorescence correlation spectroscopy experiments have revealed a second lag phase during in vitro α(1)-antitrypsin polymerization associated with the formation of smaller oligomers that then condense to form longer polymers [Purkayastha, P., Klemke, J. W., Lavender, S., Oyola, R., Cooperman, B. S., and Gai, F. (2005). Alpha 1-antitrypsin polymerization: A fluorescence correlation spectroscopic study. Biochemistry44, 2642-2649.]. SMF studies of in vitro neuroserpin polymerization have confirmed that a monomeric intermediate is required for polymer formation while providing a test of proposed polymerization mechanisms [Chiou, A., Hägglöf, P., Orte, A., Chen, A. Y., Dunne, P. D., Belorgey, D., Karlsson-Li, S., Lomas, D., and Klenerman, D. (2009). Probing neuroserpin polymerization and interaction with amyloid-beta peptides using single molecule fluorescence. Biophys. J.97, 2306-2315.]. SMF has also been used to monitor protease-serpin interactions. Single pair Förster resonance energy transfer studies of covalent protease-serpin complexes suggest that the extent of protease structural disruption in the complex is protease dependent [Liu, L., Mushero, N., Hedstrom, L., and Gershenson, A. (2006). Conformational distributions of protease-serpin complexes: A partially translocated complex. Biochemistry45, 10865-10872.]. SMF techniques are still evolving and the combination of SMF with encapsulation methods has the potential to provide more detailed information on the conformational changes associated with serpin polymerization, protease-serpin complex formation, and serpin folding.

The folding of single domain proteins--have we reached a consensus?

Rather than stressing the most recent advances in the field, this review highlights the fundamental topics where disagreement remains and where adequate experimental data are lacking. These topics include properties of the denatured state and the role of residual structure, the nature of the fundamental steps and barriers, the extent of pathway heterogeneity and non-native interactions, recent comparisons between theory and experiment, and finally, dynamical properties of the folding reaction.

DNA search efficiency is modulated by charge composition and distribution in the intrinsically disordered tail

Intrinsically disordered tails are common in DNA-binding proteins and can affect their search efficiency on nonspecific DNA by promoting the brachiation dynamics of intersegment transfer. During brachiation, the protein jumps between distant DNA regions via an intermediate state in which the tail and globular moieties are bound to different DNA segments. While the disordered tail must be long and positively charged to facilitate DNA search, the effect of its residue sequence on brachiation is unknown. We explored this issue using the NK-2 and Antp homeodomain transcription factors. We designed 566 NK-2 tail-variants and 55 Antp tail-variants having different net charges and positive charge distributions and studied their dynamics and DNA search efficiencies using coarse-grained molecular dynamics simulations. More intersegment transfers occur when the tail is moderately positively charged and the positive charges are clustered together in the middle of the tail or towards its N terminus. The presence of a negatively charged residue does not significantly affect protein brachiation, although it is likely that the presence of many negatively charged residues will complicate the DNA search mechanism. A bioinformatic analysis of 1,384 wild-type homeodomains illustrates that the charge composition and distribution in their N-tail sequences are consistent with an optimal charge pattern to promote intersegment transfer. Our study thus indicates that the residue sequence of the disordered tails of DNA-binding proteins has unique characteristics that were evolutionarily selected to achieve optimized function and suggests that the sequence-structure-function paradigm known for structured proteins is valid for intrinsically disordered proteins as well.