Structural Comparison of the Unstable drkN SH3 Domain and a Stable Mutant†,‡

Interactions of the N- and C-Terminal SH3 Domains of Drosophila Drk with the Proline-Rich Peptides from Sos and Dos

Drk, a homologue of human GRB2 in Drosophila, receives signals from outside the cells through the interaction of its SH2 domain with the phospho-tyrosine residues in the intracellular regions of receptor tyrosine kinases (RTKs) such as Sevenless, and transduces the signals downstream through the association of its N- and C-terminal SH3 domains (Drk-NSH3 and Drk-CSH3, respectively) with proline-rich motifs (PRMs) in Son of Sevenless (Sos) or Daughter of Sevenless (Dos). Isolated Drk-NSH3 exhibits a conformational equilibrium between the folded and unfolded states, while Drk-CSH3 adopts only a folded confirmation. Drk interacts with PRMs of the PxxPxR motif in Sos and the PxxxRxxKP motif in Dos. Our previous study has shown that Drk-CSH3 can bind to Sos, but the interaction between Drk-NSH3 and Dos has not been investigated. To assess the affinities of both SH3 domains towards Sos and Dos, we conducted NMR titration experiments using peptides derived from Sos and Dos. Sos-S1 binds to Drk-NSH3 with the highest affinity, strongly suggesting that the Drk-Sos multivalent interaction is initiated by the binding of Sos-S1 and NSH3. Our results also revealed that the two Sos-derived PRMs clearly favour NSH3 for binding, whereas the two Dos-derived PRMs show almost similar affinity for NSH3 and CSH3. We have also performed docking simulations based on the chemical shift perturbations caused by the addition of Sos- and Dos-derived peptides. Finally, we discussed the various modes in the interactions of Drk with Sos/Dos.

The Analytical Flory Random Coil Is a Simple-to-Use Reference Model for Unfolded and Disordered Proteins

Denatured, unfolded, and intrinsically disordered proteins (collectively referred to here as unfolded proteins) can be described using analytical polymer models. These models capture various polymeric properties and can be fit to simulation results or experimental data. However, the model parameters commonly require users' decisions, making them useful for data interpretation but less clearly applicable as stand-alone reference models. Here we use all-atom simulations of polypeptides in conjunction with polymer scaling theory to parameterize an analytical model of unfolded polypeptides that behave as ideal chains (ν = 0.50). The model, which we call the analytical Flory random coil (AFRC), requires only the amino acid sequence as input and provides direct access to probability distributions of global and local conformational order parameters. The model defines a specific reference state to which experimental and computational results can be compared and normalized. As a proof-of-concept, we use the AFRC to identify sequence-specific intramolecular interactions in simulations of disordered proteins. We also use the AFRC to contextualize a curated set of 145 different radii of gyration obtained from previously published small-angle X-ray scattering experiments of disordered proteins. The AFRC is implemented as a stand-alone software package and is also available via a Google Colab notebook. In summary, the AFRC provides a simple-to-use reference polymer model that can guide intuition and aid in interpreting experimental or simulation results.

The analytical Flory random coil is a simple-to-use reference model for unfolded and disordered proteins

Denatured, unfolded, and intrinsically disordered proteins (collectively referred to here as unfolded proteins) can be described using analytical polymer models. These models capture various polymeric properties and can be fit to simulation results or experimental data. However, the model parameters commonly require users' decisions, making them useful for data interpretation but less clearly applicable as stand-alone reference models. Here we use all-atom simulations of polypeptides in conjunction with polymer scaling theory to parameterize an analytical model of unfolded polypeptides that behave as ideal chains (ν = 0.50). The model, which we call the analytical Flory Random Coil (AFRC), requires only the amino acid sequence as input and provides direct access to probability distributions of global and local conformational order parameters. The model defines a specific reference state to which experimental and computational results can be compared and normalized. As a proof-of-concept, we use the AFRC to identify sequence-specific intramolecular interactions in simulations of disordered proteins. We also use the AFRC to contextualize a curated set of 145 different radii of gyration obtained from previously published small-angle X-ray scattering experiments of disordered proteins. The AFRC is implemented as a stand-alone software package and is also available via a Google colab notebook. In summary, the AFRC provides a simple-to-use reference polymer model that can guide intuition and aid in interpreting experimental or simulation results.

Disordered proteins mitigate the temperature dependence of site-specific binding free energies

Biophysical characterization of protein-protein interactions involving disordered proteins is challenging. A common simplification is to measure the thermodynamics and kinetics of disordered site binding using peptides containing only the minimum residues necessary. We should not assume, however, that these few residues tell the whole story. Son of sevenless, a multidomain signaling protein from Drosophila melanogaster, is critical to the mitogen-activated protein kinase pathway, passing an external signal to Ras, which leads to cellular responses. The disordered 55 kDa C-terminal domain of Son of sevenless is an autoinhibitor that blocks guanidine exchange factor activity. Activation requires another protein, Downstream of receptor kinase (Drk), which contains two Src homology 3 domains. Here, we utilized NMR spectroscopy and isothermal titration calorimetry to quantify the thermodynamics and kinetics of the N-terminal Src homology 3 domain binding to the strongest sites incorporated into the flanking disordered sequences. Comparing these results to those for isolated peptides provides information about how the larger domain affects binding. The affinities of sites on the disordered domain are like those of the peptides at low temperatures but less sensitive to temperature. Our results, combined with observations showing that intrinsically disordered proteins become more compact with increasing temperature, suggest a mechanism for this effect.

Elucidating the mechanisms underlying protein conformational switching using NMR spectroscopy

How proteins switch between various ligand-free and ligand-bound structures has been a key biophysical question ever since the postulation of the Monod-Wyman-Changeux and Koshland-Nemethy-Filmer models over six decades ago. The ability of NMR spectroscopy to provide structural and kinetic information on biomolecular conformational exchange places it in a unique position as an analytical tool to interrogate the mechanisms of biological processes such as protein folding and biomolecular complex formation. In addition, recent methodological developments in the areas of saturation transfer and relaxation dispersion have expanded the scope of NMR for probing the mechanics of transitions in systems where one or more states constituting the exchange process are sparsely populated and 'invisible' in NMR spectra. In this review, we highlight some of the strategies available from NMR spectroscopy for examining the nature of multi-site conformational exchange, using five case studies that have employed NMR, either in isolation, or in conjunction with other biophysical tools.

Structural basis of protein substrate processing by human mitochondrial high-temperature requirement A2 protease

Protein aggregates are often toxic, leading to impaired cellular activities and disease. The human HtrA2 trimeric enzyme cleaves such aggregates, and mutations in HtrA2 are causative for various neurodegenerative disorders, such as Parkinson’s disease and essential tremor. The mechanism by which cleavage occurs has been studied using small peptides, but little information is available as to how HtrA2 protects cells from the pathologic effects of aggregation involving protein molecules that can form well-folded structures. Using solution NMR spectroscopy, we investigated the structural dynamics of the interaction between HtrA2 and a model protein substrate, demonstrating that HtrA2 preferentially binds to an unfolded substrate ensemble and providing insights into how HtrA2 function is regulated.

The human high-temperature requirement A2 (HtrA2) protein is a trimeric protease that cleaves misfolded proteins to protect cells from stresses caused by toxic, proteinaceous aggregates, and the aberrant function of HtrA2 is closely related to the onset of neurodegenerative disorders. Our methyl-transverse relaxation optimized spectroscopy (TROSY)–based NMR studies using small-peptide ligands have previously revealed a stepwise activation mechanism involving multiple distinct conformational states. However, very little is known about how HtrA2 binds to protein substrates and if the distinct conformational states observed in previous peptide studies might be involved in the processing of protein clients. Herein, we use solution-based NMR spectroscopy to investigate the interaction between the N-terminal Src homology 3 domain from downstream of receptor kinase (drk) with an added C-terminal HtrA2-binding motif (drkN SH3-PDZbm) that exhibits marginal folding stability and serves as a mimic of a physiological protein substrate. We show that drkN SH3-PDZbm binds to HtrA2 via a two-pronged interaction, involving both its C-terminal PDZ-domain binding motif and a central hydrophobic region, with binding occurring preferentially via an unfolded ensemble of substrate molecules. Multivalent interactions between several clients and a single HtrA2 trimer significantly stimulate the catalytic activity of HtrA2, suggesting that binding avidity plays an important role in regulating substrate processing. Our results provide a thermodynamic, kinetic, and structural description of the interaction of HtrA2 with protein substrates and highlight the importance of a trimeric architecture for function as a stress-protective protease that mitigates aggregation.

Atomic view of cosolute-induced protein denaturation probed by NMR solvent paramagnetic relaxation enhancement

The cosolvent effect arises from the interaction of cosolute molecules with a protein and alters the equilibrium between native and unfolded states. Denaturants shift the equilibrium toward the latter, while osmolytes stabilize the former. The molecular mechanism whereby cosolutes perturb protein stability is still the subject of considerable debate. Probing the molecular details of the cosolvent effect is experimentally challenging as the interactions are very weak and transient, rendering them invisible to most conventional biophysical techniques. Here, we probe cosolute-protein interactions by means of NMR solvent paramagnetic relaxation enhancement together with a formalism we recently developed to quantitatively describe, at atomic resolution, the energetics and dynamics of cosolute-protein interactions in terms of a concentration normalized equilibrium average of the interspin distance, [Formula: see text], and an effective correlation time, τ_c The system studied is the metastable drkN SH3 domain, which exists in dynamic equilibrium between native and unfolded states, thereby permitting us to probe the interactions of cosolutes with both states simultaneously under the same conditions. Two paramagnetic cosolute denaturants were investigated, one neutral and the other negatively charged, differing in the presence of a carboxyamide group versus a carboxylate. Our results demonstrate that attractive cosolute-protein backbone interactions occur largely in the unfolded state and some loop regions in the native state, electrostatic interactions reduce the [Formula: see text] values, and temperature predominantly impacts interactions with the unfolded state. Thus, destabilization of the native state in this instance arises predominantly as a consequence of interactions of the cosolutes with the unfolded state.

Danio rerio Oocytes for Eukaryotic In-Cell NMR

Understanding how the crowded and complex cellular milieu affects protein stability and dynamics has only recently become possible by using techniques such as in-cell nuclear magnetic resonance. However, the combination of stabilizing and destabilizing interactions makes simple predictions difficult. Here we show the potential of Danio rerio oocytes as an in-cell nuclear magnetic resonance model that can be widely used to measure protein stability and dynamics. We demonstrate that in eukaryotic oocytes, which are 3-6-fold less crowded than other cell types, attractive chemical interactions still dominate effects on protein stability and slow tumbling times, compared to the effects of dilute buffer.

A Unified De Novo Approach for Predicting the Structures of Ordered and Disordered Proteins

As recognition of the abundance and relevance of intrinsically disordered proteins (IDPs) continues to grow, demand increases for methods that can rapidly predict the conformational ensembles populated by these proteins. To date, IDP simulations have largely been dominated by molecular dynamics (MD) simulations, which require significant compute times and/or complex hardware. Recent developments in MD have afforded methods capable of simulating both ordered and disordered proteins, yet to date, accurate fold prediction from a sequence has been dominated by Monte Carlo (MC)-based methods such as Rosetta. To overcome the limitations of current approaches in IDP simulation using Rosetta while maintaining its utility for modeling folded domains, we developed PyRosetta-based algorithms that allow for the accurate de novo prediction of proteins across all degrees of foldedness along with structural ensembles of disordered proteins. Our simulations have accuracy comparable to state-of-the-art MD with vastly reduced computational demands.

Extended Experimental Inferential Structure Determination Method in Determining the Structural Ensembles of Disordered Protein States

Proteins with intrinsic or unfolded state disorder comprise a new frontier in structural biology, requiring the characterization of diverse and dynamic structural ensembles. We introduce a comprehensive Bayesian framework, the Extended Experimental Inferential Structure Determination (X-EISD) method, that calculates the maximum log-likelihood of a disordered protein ensemble. X-EISD accounts for the uncertainties of a range of experimental data and back-calculation models from structures, including NMR chemical shifts, J-couplings, Nuclear Overhauser Effects (NOEs), paramagnetic relaxation enhancements (PREs), residual dipolar couplings (RDCs), hydrodynamic radii (R _h ), single molecule fluorescence Förster resonance energy transfer (smFRET) and small angle X-ray scattering (SAXS). We apply X-EISD to the joint optimization against experimental data for the unfolded drkN SH3 domain and find that combining a local data type, such as chemical shifts or J-couplings, paired with long-ranged restraints such as NOEs, PREs or smFRET, yields structural ensembles in good agreement with all other data types if combined with representative IDP conformers.

Assessing Site-Specific Enhancements Imparted by Hyperpolarized Water in Folded and Unfolded Proteins by 2D HMQC NMR

Hyperpolarized water can be a valuable aid in protein NMR, leading to amide group ¹H polarizations that are orders of magnitude larger than their thermal counterparts. Suitable procedures can exploit this to deliver 2D ¹H–¹⁵N correlations with good resolution and enhanced sensitivity. These enhancements depend on the exchange rates between the amides and the water, thereby yielding diagnostic information about solvent accessibility. This study applied this “HyperW” method to four proteins exhibiting a gamut of exchange behaviors: PhoA^(350–471), an unfolded 122-residue fragment; barstar, a fully folded ribonuclease inhibitor; R17, a 13.3 kDa system possessing folded and unfolded forms under slow interconversion; and drkN SH3, a protein domain whose folded and unfolded forms interchange rapidly and with temperature-dependent population ratios. For PhoA4^(350–471) HyperW sensitivity enhancements were ≥300×, as expected for an unfolded protein sequence. Though fully folded, barstar also exhibited substantial enhancements; these, however, were not uniform and, according to CLEANEX experiments, reflected the solvent-exposed residues. R17 showed the expected superposition of ≥100-fold enhancements for its unfolded form, coexisting with more modest enhancements for their folded counterparts. Unexpected, however, was the behavior of drkN SH3, for which HyperW enhanced the unfolded but, surprisingly, enhanced even more certain folded protein sites. These preferential enhancements were repeatedly and reproducibly observed. A number of explanations—including three-site exchange magnetization transfers between water and the unfolded and folded states; cross-correlated relaxation processes from hyperpolarized “structural” waters and labile side-chain protons; and the possibility that faster solvent exchange rates characterize certain folded sites over their unfolded counterparts—are considered to account for them.

Rapid Quantification of Protein-Ligand Binding via 19F NMR Lineshape Analysis

Fluorine incorporation is ideally suited to many NMR techniques, and incorporation of fluorine into proteins and fragment libraries for drug discovery has become increasingly common. Here, we use one-dimensional ¹⁹F NMR lineshape analysis to quantify the kinetics and equilibrium thermodynamics for the binding of a fluorine-labeled Src homology 3 (SH3) protein domain to four proline-rich peptides. SH3 domains are one of the largest and most well-characterized families of protein recognition domains and have a multitude of functions in eukaryotic cell signaling. First, we showe that fluorine incorporation into SH3 causes only minor structural changes to both the free and bound states using amide proton temperature coefficients. We then compare the results from lineshape analysis of one-dimensional ¹⁹F spectra to those from two-dimensional ¹H-¹⁵N heteronuclear single quantum coherence spectra. Their agreement demonstrates that one-dimensional ¹⁹F lineshape analysis is a robust, low-cost, and fast alternative to traditional heteronuclear single quantum coherence-based experiments. The data show that binding is diffusion limited and indicate that the transition state is highly similar to the free state. We also measured binding as a function of temperature. At equilibrium, binding is enthalpically driven and arises from a highly positive activation enthalpy for association with small entropic contributions. Our results agree with those from studies using different techniques, providing additional evidence for the utility of ¹⁹F NMR lineshape analysis, and we anticipate that this analysis will be an effective tool for rapidly characterizing the energetics of protein interactions.

Sequence Reversal Prevents Chain Collapse and Yields Heat-Sensitive Intrinsic Disorder

Sequence patterns of charge, hydrophobicity, hydrogen bonding, and other amino acid physicochemical properties contribute to mechanisms of protein folding, but how sequence composition and patterns influence the conformational dynamics of the denatured state ensemble is not fully understood. To investigate structure-sequence relationships in the denatured state, we reversed the sequence of staphylococcal nuclease and characterized its structure, thermodynamic character, and hydrodynamic radius using circular dichroism spectroscopy, dynamic light scattering, analytical ultracentrifugation, and size-exclusion chromatography as a function of temperature. The macromolecular size of "Retro-nuclease" is highly expanded in solution with characteristics similar to biological intrinsically disordered proteins. In contradistinction to a disordered state, Retro-nuclease exhibits a broad sigmoid transition of its hydrodynamic dimensions as temperature is increased, indicating a thermodynamically controlled compaction. Counterintuitively, the magnitude of these temperature-induced hydrodynamic changes exceed that observed from thermal denaturation of folded unaltered staphylococcal nuclease. Undetectable by calorimetry and intrinsic tryptophan fluorescence, the lack of heat capacity or fluorescence changes throughout the thermal transition indicate canonical hydrophobic collapse did not drive the Retro-nuclease structural transitions. Temperature-dependent circular dichroism spectroscopy performed on Retro-nuclease and computer simulations correlate to temperature sensitivity in the intrinsic sampling of backbone conformations for polyproline II and α-helix. The experimental results indicate a role for sequence direction in mediating the collapse of the polypeptide chain, whereas the simulation trends illustrate the generality of the observed heat effects on disordered protein structure.

Conserved conformational selection mechanism of Hsp70 chaperone-substrate interactions

Molecular recognition is integral to biological function and frequently involves preferred binding of a molecule to one of several exchanging ligand conformations in solution. In such a process the bound structure can be selected from the ensemble of interconverting ligands a priori (conformational selection, CS) or may form once the ligand is bound (induced fit, IF). Here we focus on the ubiquitous and conserved Hsp70 chaperone which oversees the integrity of the cellular proteome through its ATP-dependent interaction with client proteins. We directly quantify the flux along CS and IF pathways using solution NMR spectroscopy that exploits a methyl TROSY effect and selective isotope-labeling methodologies. Our measurements establish that both bacterial and human Hsp70 chaperones interact with clients by selecting the unfolded state from a pre-existing array of interconverting structures, suggesting a conserved mode of client recognition among Hsp70s and highlighting the importance of molecular dynamics in this recognition event.

Proteins are the workhorses of a cell and are involved in almost all biological processes. Newly made proteins need to ‘fold’ into precise three-dimensional shapes in order to carry out their roles. However, proteins sometimes fold incorrectly or unfold. These protein forms are not able to work effectively and in some cases may even cause diseases.

Chaperone proteins help other proteins to fold correctly and are found in living organisms ranging in complexity from bacteria to humans. There are many different types of chaperones that play different roles inside cells. One, called Hsp70, binds to proteins that are incorrectly folded to help them to mature into their correct structures. However, it was not clear whether Hsp70 can also associate with the mature, correctly folded form of the proteins.

A technique called Nuclear Magnetic Resonance (NMR) spectroscopy can distinguish between mature, unfolded and chaperone-bound forms of the same protein. Sekhar et al. therefore used NMR to investigate which forms of a protein Hsp70 binds to. This revealed that both the bacterial and human versions of the Hsp70 chaperone interact only with unfolded proteins.

The results presented by Sekhar et al. also explain why Hsp70 does not disrupt the routine workings of the cell: because it does not bind to mature forms of proteins. These observations extend our understanding of how chaperones assist in folding proteins, and fit into a broader research theme exploring how proteins recognize one another. It will now be interesting to see whether the same mechanism holds for more complex forms of proteins, such as aggregates, or larger protein structures with regions of both folded and unfolded elements.

Cosolutes, Crowding, and Protein Folding Kinetics

Long accepted as the most important interaction, recent work shows that steric repulsions alone cannot explain the effects of macromolecular cosolutes on the equilibrium thermodynamics of protein stability. Instead, chemical interactions have been shown to modulate, and even dominate, crowding-induced steric repulsions. Here, we use 19F NMR to examine the effects of small and large cosolutes on the kinetics of protein folding and unfolding using the metastable 7 kDa N-terminal SH3 domain of the Drosophila signaling protein drk (SH3), which folds by a two-state mechanism. The small cosolutes consist of trimethylamine N-oxide and sucrose, which increase equilibrium protein stability, and urea, which destabilizes proteins. The macromolecules comprise the stabilizing sucrose polymer, Ficoll, and the destabilizing globular protein, lysozyme. We assessed the effects of these cosolutes on the differences in free energy between the folded state and the transition state and between the unfolded ensemble and the transition state. We then examined the temperature dependence to assess changes in activation enthalpy and entropy. The enthalpically mediated effects are more complicated than suggested by equilibrium measurements. We also observed enthalpic effects with the supposedly inert sucrose polymer, Ficoll, that arise from its macromolecular nature. Assessment of activation entropies shows important contributions from solvent and cosolute, in addition to the configurational entropy of the protein that, again, cannot be gleaned from equilibrium data. Comparing the effects of Ficoll to those of the more physiologically relevant cosolute lysozyme reveals that synthetic polymers are not appropriate models for understanding the kinetics of protein folding in cells.

Conformations of a Metastable SH3 Domain Characterized by smFRET and an Excluded-Volume Polymer Model

Conformational states of the metastable drkN SH3 domain were characterized using single-molecule fluorescence techniques. Under nondenaturing conditions, two Förster resonance energy transfer (FRET) populations were observed that corresponded to a folded and an unfolded state. FRET-estimated radii of gyration and hydrodynamic radii estimated by fluorescence correlation spectroscopy of the two coexisting conformations are in agreement with previous ensemble x-ray scattering and NMR measurements. Surprisingly, when exposed to high concentrations of urea and GdmCl denaturants, the protein still exhibits two distinct FRET populations. The dominant conformation is expanded, showing a low FRET efficiency, consistent with the expected behavior of a random chain with excluded volume. However, approximately one-third of the drkN SH3 conformations showed high, nearly 100%, FRET efficiency, which is shown to correspond to denaturation-induced looped conformations that remain stable on a timescale of at least 100 μs. These loops may contain interconverting conformations that are more globally collapsed, hairpin-like, or circular, giving rise to the observed heterogeneous broadening of this population. Although the underlying mechanism of chain looping remains elusive, FRET experiments in formamide and dimethyl sulfoxide suggest that interactions between hydrophobic groups in the distal regions may play a significant role in the formation of the looped state.

Understanding Protein Structure Deformation on the Surface of Gold Nanoparticles of Varying Size

Gold nanoparticles (AuNPs) have been of recent interest due to their unique optical properties and their biocompatibility. Biomolecules spontaneously adsorb to their surface, a trait that could potentially be exploited for drug targeting. Currently, it is unclear whether protein-AuNP interactions at the nanoparticle surface are dependent on nanoparticle size. In this work, we investigate whether varying surface curvature can induce protein unfolding and multilayer binding in citrate-coated AuNPs of various sizes. A recently developed NMR-based approach was utilized to determine the adsorption capacity, and protein NMR spectra were compared to determine whether nanoparticle size influences protein interactions at the surface. In addition, transmission electron microscopy (TEM) and dynamic light scattering (DLS) were employed to corroborate the NMR studies. Over a broad range of AuNP sizes (14-86 nm), we show that adsorption capacity can be predicted by assuming that proteins are compact and globular on the nanoparticle surface. Additionally, roughly one layer of proteins is adsorbed regardless of AuNP size. Our results hold for two proteins of significantly different sizes, GB3 (6 kDa) and bovine carbonic anhydrase (BCA, 29 kDa). However, the unstable drkN SH3 domain (ΔḠ0 ≈ 0, 7 kDa) does not appear to follow the same trend seen for stable, globular proteins. This observation suggests that unstable proteins can deform significantly when bound to AuNP surfaces. Taken together, the results of this work can be used to improve our knowledge of the mechanism of protein-AuNP interactions to optimize their use in the biomedical field.

Engineering Proteins for Thermostability with iRDP Web Server

Engineering protein molecules with desired structure and biological functions has been an elusive goal. Development of industrially viable proteins with improved properties such as stability, catalytic activity and altered specificity by modifying the structure of an existing protein has widely been targeted through rational protein engineering. Although a range of factors contributing to thermal stability have been identified and widely researched, the in silico implementation of these as strategies directed towards enhancement of protein stability has not yet been explored extensively. A wide range of structural analysis tools is currently available for in silico protein engineering. However these tools concentrate on only a limited number of factors or individual protein structures, resulting in cumbersome and time-consuming analysis. The iRDP web server presented here provides a unified platform comprising of iCAPS, iStability and iMutants modules. Each module addresses different facets of effective rational engineering of proteins aiming towards enhanced stability. While iCAPS aids in selection of target protein based on factors contributing to structural stability, iStability uniquely offers in silico implementation of known thermostabilization strategies in proteins for identification and stability prediction of potential stabilizing mutation sites. iMutants aims to assess mutants based on changes in local interaction network and degree of residue conservation at the mutation sites. Each module was validated using an extensively diverse dataset. The server is freely accessible at http://irdp.ncl.res.in and has no login requirements.

Heterogeneous binding of the SH3 client protein to the DnaK molecular chaperone

The molecular chaperone heat shock protein 70 (Hsp70) plays a vital role in cellular processes, including protein folding and assembly, and helps prevent aggregation under physiological and stress-related conditions. Although the structural changes undergone by full-length client proteins upon interaction with DnaK (i.e., Escherichia coli Hsp70) are fundamental to understand chaperone-mediated protein folding, these changes are still largely unexplored. Here, we show that multiple conformations of the SRC homology 3 domain (SH3) client protein interact with the ADP-bound form of the DnaK chaperone. Chaperone-bound SH3 is largely unstructured yet distinct from the unfolded state in the absence of DnaK. The bound client protein shares a highly flexible N terminus and multiple slowly interconverting conformations in different parts of the sequence. In all, there is significant structural and dynamical heterogeneity in the DnaK-bound client protein, revealing that proteins may undergo some conformational sampling while chaperone-bound. This result is important because it shows that the surface of the Hsp70 chaperone provides an aggregation-free environment able to support part of the search for the native state.

Towards a true protein movie: a perspective on the potential impact of the ensemble-based structure determination using exact NOEs

Confined by the Boltzmann distribution of the energies of the states, a multitude of structural states are inherent to biomolecules. For a detailed understanding of a protein's function, its entire structural landscape at atomic resolution and insight into the interconversion between all the structural states (i.e. dynamics) are required. Whereas dedicated trickery with NMR relaxation provides aspects of local dynamics, and 3D structure determination by NMR is well established, only recently have several attempts been made to formulate a more comprehensive description of the dynamics and the structural landscape of a protein. Here, a perspective is given on the use of exact NOEs (eNOEs) for the elucidation of structural ensembles of a protein describing the covered conformational space.

Sub-diffusion decays in fluorescence correlation spectroscopy: dye photophysics or protein dynamics?

Transitions between bright and dark fluorescent states of several rhodamine dyes were investigated by fluorescence correlation spectroscopy. We resolved two sub-diffusion exponential decays for free rhodamines in aqueous solutions, of which the slower component scales linearly with the viscosity of the solution. Correlation data for proteins and DNA labeled with tetramethylrhodamine were fitted with three to four exponential decays describing flickering dynamics on a time scale between 0.5 and 100 μs. We investigated the nature of these processes by performing experiments under different experimental conditions and for different samples. On the basis of how their population and lifetime change with viscosity, the oxygen content of the solution, the laser irradiance, and the detection geometry, we assigned these states, in the order of increasing lifetimes, to a triplet state, a hybrid between twisted-intramolecular-charge-transfer state and a ground state lactonic state, a lactonic state, and a photoionized state, respectively. Our data suggests that none of the observed sub-diffusion correlation decays can be directly assigned to the intramolecular dynamics of the labeled biomolecules. However, we found evidence that the intrinsic conformational dynamics of the biomolecule appears in the correlation curves as a modulation of the photophysics of the dye label. This shows the importance of accurate control measurements and appropriate modeling of the dye photophysics in fluorescence correlation studies, and it cautions against direct assignments of dark-state relaxation times to folding kinetics in proteins and nucleic acids.

Ensemble modeling of protein disordered states: experimental restraint contributions and validation

Disordered states of proteins include the biologically functional intrinsically disordered proteins and the unfolded states of normally folded proteins. In recent years, ensemble-modeling strategies using various experimental measurements as restraints have emerged as powerful means for structurally characterizing disordered states. However, these methods are still in their infancy compared with the structural determination of folded proteins. Here, we have addressed several issues important to ensemble modeling using our ENSEMBLE methodology. First, we assessed how calculating ensembles containing different numbers of conformers affects their structural properties. We find that larger ensembles have very similar properties to smaller ensembles fit to the same experimental restraints, thus allowing a considerable speed improvement in our calculations. In addition, we analyzed the contributions of different experimental restraints to the structural properties of calculated ensembles, enabling us to make recommendations about the experimental measurements that should be made for optimal ensemble modeling. The effects of different restraints, most significantly from chemical shifts, paramagnetic relaxation enhancements and small-angle X-ray scattering, but also from other data, underscore the importance of utilizing multiple sources of experimental data. Finally, we validate our ENSEMBLE methodology using both cross-validation and synthetic experimental restraints calculated from simulated ensembles. Our results suggest that secondary structure and molecular size distribution can generally be modeled very accurately, whereas the accuracy of calculated tertiary structure is dependent on the number of distance restraints used.

15N H/D-SOLEXSY experiment for accurate measurement of amide solvent exchange rates: application to denatured drkN SH3

Amide solvent exchange rates are regarded as a valuable source of information on structure/dynamics of unfolded (disordered) proteins. Proton-based saturation transfer experiments, normally used to measure solvent exchange, are known to meet some serious difficulties. The problems mainly arise from the need to (1) manipulate water magnetization and (2) discriminate between multiple magnetization transfer pathways that occur within the proton pool. Some of these issues are specific to unfolded proteins. For example, the compensation scheme used to cancel the Overhauser effect in the popular CLEANEX experiment is not designed for use with unfolded proteins. In this report we describe an alternative experimental strategy, where amide (15)N is used as a probe of solvent exchange. The experiment is performed in 50% H(2)O-50% D(2)O solvent and is based on the (HACACO)NH pulse sequence. The resulting spectral map is fully equivalent to the conventional HSQC. To fulfill its purpose, the experiment monitors the conversion of deuterated species, (15)N(D), into protonated species, (15)N(H), as effected by the solvent exchange. Conceptually, this experiment is similar to EXSY which prompted the name of (15)N(H/D)-SOLEXSY (SOLvent EXchange SpectroscopY). Of note, our experimental scheme, which relies on nitrogen rather than proton to monitor solvent exchange, is free of the complications described above. The developed pulse sequence was used to measure solvent exchange rates in the chemically denatured state of the drkN SH3 domain. The results were found to correlate well with the CLEANEX-PM data, r = 0.97, thus providing a measure of validation for both techniques. When the experimentally measured exchange rates are converted into protection factors, most of the values fall in the range 0.5-2, consistent with random-coil behavior. However, elevated values, ca. 5, are obtained for residues R38 and A39, as well as the side-chain indole of W36. This is surprising, given that high protection factors imply hydrogen bonding or hydrophobic burial not expected to occur in a chemically denatured state of a protein. We, therefore, hypothesized that elevated protection factors are an artefact arising from the calculation of the reference (random-coil) exchange rates. To confirm this hypothesis, we prepared samples of several short peptides derived from the sequence of the drkN SH3 domain; these samples were used to directly measure the reference exchange rates. The revised protection factors obtained in this manner proved to be close to 1.0. These results also have implications for the more compact unfolded state of drkN SH3, which appears to be fully permeable to water as well, with no manifestations of hydrophobic burial.

Structure and disorder in an unfolded state under nondenaturing conditions from ensemble models consistent with a large number of experimental restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include (1)H-(15)N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O(2)-induced (13)C paramagnetic shifts, hydrogen-exchange protection factors, and (15)N R(2) data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu(2+)-Ni(2+) binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative alpha-helical structure and both native-like and nonnative tertiary structure.

Improved Structural Characterizations of the drkN SH3 Domain Unfolded State Suggest a Compact Ensemble with Native-like and Non-native Structure

Due to their dynamic ensemble nature and a deficiency of experimental restraints, disordered states of proteins are difficult to characterize structurally. Here, we have expanded upon our previous work on the unfolded state of the Drosophila drk N-terminal (drkN) SH3 domain with our program ENSEMBLE, which assigns population weights to pregenerated conformers in order to calculate ensembles of structures whose properties are collectively consistent with experimental measurements. The experimental restraint set has been enlarged with newly measured paramagnetic relaxation enhancements from Cu(2+) bound to an amino terminal Cu(2+)-Ni(2+) binding (ATCUN) motif as well as nuclear Overhauser effect (NOE) and hydrogen exchange data from recent studies. In addition, two new pseudo-energy minimization algorithms have been implemented that have dramatically improved the speed of ENSEMBLE population weight assignment. Finally, we have greatly improved our conformational sampling by utilizing a variety of techniques to generate both random structures and structures that are biased to contain elements of native-like or non-native structure. Although it is not possible to uniquely define a representative structural ensemble, we have been able to assess various properties of the drkN SH3 domain unfolded state by performing ENSEMBLE minimizations of different conformer pools. Specifically, we have found that the experimental restraint set enforces a compact structural distribution that is not consistent with an overall native-like topology but shows preference for local non-native structure in the regions corresponding to the diverging turn and the beta5 strand of the folded state and for local native-like structure in the region corresponding to the beta6 and beta7 strands. We suggest that this approach could be generally useful for the structural characterization of disordered states.

Oxygen as a Paramagnetic Probe of Clustering and Solvent Exposure in Folded and Unfolded States of an SH3 Domain

The N-terminal SH3 domain of the Drosophila modular protein Drk undergoes slow exchange between a folded (Fexch) and highly populated unfolded (Uexch) state under nondenaturing buffer conditions, enabling both Fexch and Uexch states to be simultaneously monitored. The addition of dissolved oxygen, equilibrated to a partial pressure of either 30 atm or 60 atm, provides the means to study solvent exposure with atomic resolution via 13C NMR paramagnetic shifts in 1H,13C HSQC (heteronuclear single quantum coherence) spectra. Absolute differences in these paramagnetic shifts between the Fexch and Uexch states allow the discrimination of regions of the protein which undergo change in solvent exposure upon unfolding. Contact with dissolved oxygen for both the Fexch and Uexch states could also be assessed through 13C paramagnetic shifts which were normalized based on the corresponding paramagnetic shifts seen in the free amino acids. In the Fexch state, the 13C nuclei belonging to the hydrophobic core of the protein exhibited very weak normalized paramagnetic shifts while those with greater solvent accessible surface area exhibited significantly larger normalized shifts. The Uexch state displayed less varied 13C paramagnetic shifts although distinct regions of protection from solvent exposure could be identified by a lack of such shifts. These regions, which included Phe9, Thr12, Ala13, Lys21, Thr22, Ile24, Ile27, and Arg38, overlapped with those found to have residual nativelike and non-native structures in previous studies and in some cases provided novel information. Thus, the paramagnetic shifts from dissolved oxygen are highly useful in the study of a transient structure or clustering in disordered systems, where conventional NMR measurements (couplings, chemical shift deviations from random coil values, and NOEs) may give little information.

Partial unfolding of diverse SH3 domains on a wide timescale

SH3 domains are small, modular domains that are found in many proteins, especially signal transduction proteins such as tyrosine kinases. While much is known about the sequences and tertiary structures of SH3 domains, far less is known about their solution dynamics. A slow, partial unfolding event that occurs under physiological conditions was previously identified in the Hck SH3 domain using hydrogen exchange (HX) mass spectrometry (MS). To determine if this unfolding was unique to Hck SH3, HX MS was used to analyze 11 other SH3 domains: seven SH3 domains from Src-family kinases and five SH3 domains from various proteins. A wide variety of unfolding rates were found, with unfolding half-lives ranging from 1s to 1h. The Lyn and alpha-spectrin SH3 domains exhibited slow, partial unfolding in beta strands D and E and part of the RT-loop. Hck SH3 also underwent partial unfolding in the same region, implying that a unique feature in this area of the domains is responsible for the partial unfolding. Partial unfolding was, however, not a function of sequence conservation. Although the Fyn and Yes SH3 domains are very similar to Hck SH3 in sequence, they exhibited no evidence of partial unfolding. Overall, the results suggest that while the tertiary structure of SH3 domains is highly conserved, the dynamics of SH3 domains are variable.