Structural Characterization of Unfolded States of Apomyoglobin using Residual Dipolar Couplings

The effects of biliverdin on pressure-induced unfolding of apomyoglobin: The specific role of Zn2+ ions

Apomyoglobin (apoMb), a model protein in biochemistry, exhibits a strong propensity to bind various ligands, which makes it a good candidate as a carrier of bioactive hydrophobic drugs. The stability of its hydrophobic pocket determines its potential as a carrier of bioactive compounds. High pressure (HP) is a potent tool for studying protein stability, revealing the specific role of hydrophobic cavities in unfolding. We probed the effects of biliverdin (BV) binding and its complex with Zn²⁺ ions on the structure and HP stability of apoMb. CD spectroscopy and SAXS measurements revealed that BV and BV-Zn²⁺ complexes make the apoMb structure more compact with higher α-helical content. We performed in-situ HP measurements of apoMb intrinsic fluorescence to demonstrate the ability of BV to stabilise apoMb structure at HP conditions. Furthermore, the presence of Zn²⁺ within the apoMb-BV complex significantly enhances the BV stabilisation effect. In-situ visible absorption study of BV chromophore confirmed the ability of Zn²⁺ to increase the stability of apoMb-BV complex under HP: the onset of complex dissociation is shifted by ~100 MPa in the presence of Zn²⁺. By combining HP-fluorescence and HP-visible absorption spectroscopy, our strategy highlights the crucial role of tetrapyrrole-metal complexes in stabilising apoMb hydrophobic pocket.

Simple and Effective Conformational Sampling Strategy for Intrinsically Disordered Proteins Using the UNRES Web Server

Intrinsically disordered proteins (IDPs) contain more charged amino acids than folded proteins, resulting in a lack of hydrophobic core(s) and a tendency to adopt rapidly interconverting structures rather than well-defined structures. The structural heterogeneity of IDPs, encoded by the amino acid sequence, is closely related to their unique roles in biological pathways, which require them to interact with different binding partners. Recently, Robustelli and co-workers have demonstrated that a balanced all-atom force field can be used to sample heterogeneous structures of disordered proteins ( Proc. Natl. Acad. Sci. U.S.A. 2018, 115, E4758-E4766). However, such a solution requires extensive computational resources, such as Anton supercomputers. Here, we propose a simple and effective solution to sample the conformational space of IDPs using a publicly available web server, namely, the UNited-RESidue (UNRES) web server. Our proposed solution requires no investment in computational resources and no prior knowledge of UNRES. UNRES Replica Exchange Molecular Dynamics (REMD) simulations were carried out on a set of eight disordered proteins at temperatures spanning from 270 to 430 K. Utilizing the latest UNRES force field designed for structured proteins, with proper selections of temperatures, we were able to produce comparable results to all-atom force fields as reported in work done by Robustelli and co-workers. In addition, NMR observables and the radius of gyration calculated from UNRES ensembles were directly compared with the experimental data to further evaluate the accuracy of the UNRES model at all temperatures. Our results suggest that carrying out the UNRES simulations at optimal temperatures using the UNRES web server can be a good alternative to sample heterogeneous structures of IDPs.

Transient On- and Off-Pathway Protein Folding Intermediate States Characterized with NMR Relaxation Dispersion

The earliest events in the folding of a protein are in general poorly understood. We used NMR R2 relaxation dispersion experiments to study transient local collapse events in the unfolded-state (U) conformational ensemble of apomyoglobin (apoMb). Local residual secondary structure (seen in regions corresponding to the A, D, E, and H helices of the folded protein) is largely unchanged over the pH range of 2.3-2.75, yet a significant pH-dependent increase in the conformational exchange contribution to the R2 relaxation rate (Rex) indicates that transient intramolecular contacts occur on a microsecond to millisecond time scale at pH 2.75. A comparison of 15N and 13CO relaxation dispersion data at pH 2.75 for residues in the A, B, G, and H regions, which participate in the earliest folding intermediates, indicates that chain collapse and secondary structure formation are rapid and concomitant. Increasingly stabilizing conditions (lower temperature, higher pH) result in the observation of a relaxation dispersion in the C, CD, and E regions of the protein, which are known to fold at later stages. Mutation of Trp14 in the A-helix region to Ala eliminates conformational exchange throughout the protein, and the mutation of hydrophobic residues in other regions results in the selective inhibition of conformational exchange in the B, G, or H regions. The R2 dispersion data for WT apoMb at pH 2.75 and 10 °C are best fit to a four-state model ABGH ⇆ AGH ⇆ U ⇆ ABCD that includes on-pathway (AGH and ABGH) and off-pathway (ABCD) transiently folded states, both of which are required to explain the behavior of the mutant proteins. The off-pathway intermediate is destabilized at higher temperatures. Our analysis provides insights into the earliest stages of apoMb folding where the collapsing polypeptide chain samples both productive and nonproductive states with stabilized secondary structure.

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔGu°'(H2O), of UBA(1) is 2.4 kcal mol-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, kf, in the absence of a denaturant of 13,000 s-1 and a Tanford β-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone 15NH, 13CO, and 13Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

New insights into the folding–unfolding mechanism and conformations of cytochrome C

Optical trapping experiments offer new insights into the folding and unfolding of cytochrome C.

Structural and dynamic characterization of the C-terminal tail of ErbB2: Disordered but not random

ErbB2 (or HER2) is a receptor tyrosine kinase overexpressed in some breast cancers and associated with poor prognosis. Treatments targeting the receptor extracellular and kinase domains have greatly improved disease outcome in the last 20 years. In parallel, the structures of these domains have been described, enabling better mechanistic understanding of the receptor function and targeted inhibition. However, the ErbB2 disordered C-terminal cytoplasmic tail (CtErbB2) remains very poorly characterized in terms of structure, dynamics, and detailed functional mechanism. Yet, it is where signal transduction is triggered via phosphorylation of tyrosine residues and carried out via interaction with adaptor proteins. Here, we report the first description, to our knowledge, of the ErbB2 disordered tail at atomic resolution using NMR, complemented by small-angle x-ray scattering. We show that although no part of CtErbB2 has any fully populated secondary or tertiary structure, it contains several transient α-helices and numerous transient polyproline II helices, populated up to 20 and 40%, respectively, and low but significant compaction. The presence of some structural elements suggests, along the lines of the results obtained for EGFR (ErbB1), that they may have a functional role in ErbB2's autoregulation processes. In addition, the transient formation of polyproline II helices is compliant with previously suggested interactions with SH3 domains. All in all, our in-depth structural study opens perspectives in the mechanistic understanding of ErbB2.

Intrinsically Disordered Bacterial Polar Organizing Protein Z, PopZ, Interacts with Protein Binding Partners Through an N-terminal Molecular Recognition Feature

The polar organizing protein Z (PopZ) is necessary for the formation of three-dimensional microdomains at the cell poles in Caulobacter crescentus, where it functions as a hub protein that recruits multiple regulatory proteins from the cytoplasm. Although a large portion of the protein is predicted to be natively unstructured, in reconstituted systems PopZ can self-assemble into a macromolecular scaffold that directly binds to at least ten different proteins. Here we report the solution NMR structure of PopZ^Δ134-177, a truncated form of PopZ that does not self-assemble but retains the ability to interact with heterologous proteins. We show that the unbound form of PopZ^Δ134-177 is unstructured in solution, with the exception of a small amphipathic α-helix in residues M10-I17, which is included within a highly conserved region near the N-terminal. In applying NMR techniques to map the interactions between PopZ^Δ134-177 and one of its binding partners, RcdA, we find evidence that the α-helix and adjoining amino acids extending to position E23 serve as the core of the binding motif. Consistent with this, a point mutation at position I17 severely compromises binding. Our results show that a partially structured Molecular Recognition Feature (MoRF) within an intrinsically disordered domain of PopZ contributes to the assembly of polar microdomains, revealing a structural basis for complex network assembly in Alphaproteobacteria that is analogous to those formed by intrinsically disordered hub proteins in other kingdoms.

Lessons Learned from 50 Years of Hemoglobin Research: Unstirred and Cell-Free Layers, Electrostatics, Baseball Gloves, and Molten Globules

Significance: Over the past 50 years, the mechanisms for O₂ storage and transport have been determined quantitatively on distance scales from millimeters to tenths of nanometers and timescales from seconds to picoseconds. Recent Advances: In this review, I have described four key conclusions from work done by my group and our close colleagues. (i) O₂ uptake by mammalian red cells is limited by diffusion through unstirred water layers adjacent to the cell surface and across cell-free layers adjacent to vessel walls. (ii) In most vertebrates, hemoglobins (Hbs) and myoglobins (Mbs), the distal histidine at the E7 helical position donates a strong hydrogen bond to bound O₂, which selectively enhances O₂ affinity, prevents carbon monoxide poisoning, and markedly slows autoxidation. (iii) O₂ binding to mammalian Hbs and Mbs occurs by migration of the ligand through a channel created by upward rotation of the His(E7) side chain, capture in the empty space of the distal pocket, and then coordination with the ferroprotoporphyrin IX (heme) iron atom. (iv) The assembly of Mbs and Hbs occurs by formation of molten globule intermediates, in which the N- and C-terminal helices have almost fully formed secondary structures, but the heme pockets are disordered and followed by high-affinity binding of heme. Critical Issues: These conclusions indicate that there are often compromises between O₂ transport function, holoprotein stability, and the efficiency of assembly. Future Directions: However, the biochemical mechanisms underlying these conclusions provide the framework for understanding globin evolution in greater detail and for engineering more efficient and stable globins.

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

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

The dormancy-specific regulator, SutA, is intrinsically disordered and modulates transcription initiation in Pseudomonas aeruginosa

Though most bacteria in nature are nutritionally limited and grow slowly, our understanding of core processes like transcription comes largely from studies in model organisms doubling rapidly. We previously identified a small protein of unknown function, SutA, in a screen of proteins synthesized in Pseudomonas aeruginosa during dormancy. SutA binds RNA polymerase (RNAP), causing widespread changes in gene expression, including upregulation of the ribosomal RNA genes. Here, using biochemical and structural methods, we examine how SutA interacts with RNAP and the functional consequences of these interactions. We show that SutA comprises a central α-helix with unstructured N- and C-terminal tails, and binds to the β1 domain of RNAP. It activates transcription from the rrn promoter by both the housekeeping sigma factor holoenzyme (Eσ⁷⁰ ) and the stress sigma factor holoenzyme (Eσ^S ) in vitro, but has a greater impact on Eσ^S . In both cases, SutA appears to affect intermediates in the open complex formation and its N-terminal tail is required for activation. The small magnitudes of in vitro effects are consistent with a role in maintaining activity required for homeostasis during dormancy. Our results add SutA to a growing list of transcription regulators that use their intrinsically disordered regions to remodel transcription complexes.

Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.

Mutations in the dimer interfaces of the dengue virus capsid protein affect structural stability and impair RNA-capsid interaction

The dengue virus 2 capsid protein (DENV2C) plays a primary structural role in the protection of the viral genome and is crucial for nucleocapsid assembly. In this study, we generated single mutants of DENV2C at L50 and L54 residues of the α2 helix, which was shown to interfere with the integration of the capsid into lipid droplets, and at residues L81 and I88 located in the α4 helix, which was shown to affect viral assembly. We demonstrated that the oligomeric states of DENV2C and its mutants exist primarily in the dimeric state in solution. All single-point mutations introduced in DENV2C promoted reduction in protein stability, an effect that was more pronounced for the L81N and I88N mutants, but not protein unfolding. All the single-point mutations affected the ability of DEN2C to interact with RNA. We concluded that mutations in the α2-α2' and α4-α4' dimer interfaces of DENV2C affect the structural stability of the protein and impair RNA-capsid interaction. These effects were more pronounced for mutations at the L81 and I88 residues in the α4 helix. These results indicate the importance of the α4-α4' dimer interface, which could be studied as a potential target for drug design in the future.

Characterization of intrinsically disordered proteins and their dynamic complexes: From in vitro to cell-like environments

Over the last two decades, it has become increasingly clear that a large fraction of the human proteome is intrinsically disordered or contains disordered segments of significant length. These intrinsically disordered proteins (IDPs) play important regulatory roles throughout biology, underlining the importance of understanding their conformational behavior and interaction mechanisms at the molecular level. Here we review recent progress in the NMR characterization of the structure and dynamics of IDPs in various functional states and environments. We describe the complementarity of different NMR parameters for quantifying the conformational propensities of IDPs in their isolated and phosphorylated states, and we discuss the challenges associated with obtaining structural models of dynamic protein-protein complexes involving IDPs. In addition, we review recent progress in understanding the conformational behavior of IDPs in cell-like environments such as in the presence of crowding agents, in membrane-less organelles and in the complex environment of the human cell.

High-Resolution 2D NMR of Disordered Proteins Enhanced by Hyperpolarized Water

This study demonstrates the usefulness derived from relying on hyperpolarized water obtained by dissolution DNP, for site-resolved biophysical NMR studies of intrinsically disordered proteins. Thanks to the facile amide-solvent exchange experienced by protons in these proteins, 2D NMR experiments that like HMQC rely on the polarization of the amide protons, can be enhanced using hyperpolarized water by several orders of magnitude over their conventional counterparts. Optimizations of the DNP procedure and of the subsequent injection into the protein sample are necessary to achieve these gains while preserving state-of-the-art resolution; procedures enabling this transfer of the hyperpolarized water and the achievement of foamless hyperpolarized protein solutions are demonstrated. These protocols are employed to collect 2D ¹⁵N-¹H HMQC NMR spectra of α-synuclein, showing residue-specific enhancements ≥100× over their thermal counterparts. These enhancements, however, vary considerably throughout the residues. The biophysics underlying this residue-specific behavior upon injection of hyperpolarized water is theoretically examined, the information that it carries is compared with results arising from alternative methods, and its overall potential is discussed.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign ¹H and ¹⁵N backbone amide and ¹³C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

Monitoring protein folding through high pressure NMR spectroscopy

High-pressure is a well-known perturbation method used to destabilize globular proteins. It is perfectly reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. In contrast to other perturbation methods such as heat or chemical denaturant that destabilize protein structures uniformly, pressure exerts local effects on regions or domains of a protein containing internal cavities. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility to monitor at a residue level the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. High-pressure NMR experiments can now be routinely performed, owing to the recent development of commercially available high-pressure sample cells. This review summarizes recent advances and some future directions of high-pressure NMR techniques for the characterization at atomic resolution of the energy landscape of protein folding.

How Does Your Protein Fold? Elucidating the Apomyoglobin Folding Pathway

Although each type of protein fold and in some cases individual proteins within a fold classification can have very different mechanisms of folding, the underlying biophysical and biochemical principles that operate to cause a linear polypeptide chain to fold into a globular structure must be the same. In an aqueous solution, the protein takes up the thermodynamically most stable structure, but the pathway along which the polypeptide proceeds in order to reach that structure is a function of the amino acid sequence, which must be the final determining factor, not only in shaping the final folded structure, but in dictating the folding pathway. A number of groups have focused on a single protein or group of proteins, to determine in detail the factors that influence the rate and mechanism of folding in a defined system, with the hope that hypothesis-driven experiments can elucidate the underlying principles governing the folding process. Our research group has focused on the folding of the globin family of proteins, and in particular on the monomeric protein apomyoglobin. Apomyoglobin (apoMb) folds relatively slowly (∼2 s) via an ensemble of obligatory intermediates that form rapidly after the initiation of folding. The folding pathway can be dissected using rapid-mixing techniques, which can probe processes in the millisecond time range. Stopped-flow measurements detected by circular dichroism (CD) or fluorescence spectroscopy give information on the rates of folding events. Quench-flow experiments utilize the differential rates of hydrogen-deuterium exchange of amide protons protected in parts of the structure that are folded early; protection of amides can be detected by mass spectrometry or proton nuclear magnetic resonance spectroscopy (NMR). In addition, apoMb forms an intermediate at equilibrium at pH ∼ 4, which is sufficiently stable for it to be structurally characterized by solution methods such as CD, fluorescence and NMR spectroscopies, and the conformational ensembles formed in the presence of denaturing agents and low pH can be characterized as models for the unfolded states of the protein. Newer NMR techniques such as measurement of residual dipolar couplings in the various partly folded states, and relaxation dispersion measurements to probe invisible states present at low concentrations, have contributed to providing a detailed picture of the apomyoglobin folding pathway. The research summarized in this Account was aimed at characterizing and comparing the equilibrium and kinetic intermediates both structurally and dynamically, as well as delineating the complete folding pathway at a residue-specific level, in order to answer the question: "What is it about the amino acid sequence that causes each molecule in the unfolded protein ensemble to start folding, and, once started, to proceed towards the formation of the correctly folded three-dimensional structure?"

Folding of apomyoglobin: Analysis of transient intermediate structure during refolding using quick hydrogen deuterium exchange and NMR

The structures of apomyoglobin folding intermediates have been widely analyzed using physical chemistry methods including fluorescence, circular dichroism, small angle X-ray scattering, NMR, mass spectrometry, and rapid mixing. So far, at least two intermediates (on sub-millisecond- and millisecond-scales) have been demonstrated for apomyoglobin folding. The combination of pH-pulse labeling and NMR is a useful tool for analyzing the kinetic intermediates at the atomic level. Its use has revealed that the latter-phase kinetic intermediate of apomyoglobin (6 ms) was composed of helices A, B, G and H, whereas the equilibrium intermediate, called the pH 4 molten-globule intermediate, was composed mainly of helices A, G and H. The improved strategy for the analysis of the kinetic intermediate was developed to include (1) the dimethyl sulfoxide method, (2) data processing with the various labeling times, and (3) a new in-house mixer. Particularly, the rapid mixing revealed that helices A and G were significantly more protected at the earlier stage (400 µs) of the intermediate (former-phase intermediate) than the other helices. Mutation studies, where each hydrophobic residue was replaced with an alanine in helices A, B, E, F, G and H, indicated that both non-native and native-like structures exist in the latter-phase folding intermediate. The N-terminal part of helix B is a weak point in the intermediate, and the docking of helix E residues to the core of the A, B, G and H helices was interrupted by a premature helix B, resulting in the accumulation of the intermediate composed of helices A, B, G and H. The prediction-based protein engineering produced important mutants: Helix F in a P88K/A90L/S92K/A94L mutant folded in the latter-phase intermediate, although helix F in the wild type does not fold even at the native state. Furthermore, in the L11G/W14G/A70L/G73W mutant, helix A did not fold but helix E did, which is similar to what was observed in the kinetic intermediate of apoleghemoglobin. Thus, this protein engineering resulted in a changed structure for the apomyoglobin folding intermediate.

Statistical Mechanical Model for pH-Induced Protein Folding: Application to Apomyoglobin

Despite the major role of pH in protein folding and stability, a quantitative understanding of the pH-induced protein folding mechanism remains elusive. Two conventional models, the Monod-Wyman-Changeux and Linderstrøm-Lang smeared charge models, respectively, have been used to analyze the formation/disruption of specific native structures and fluctuating non-native states. However, there are only a few models that can represent the overall kinetic events of folding/unfolding independent of the properties of relevant molecular species, which has hampered the efforts to systematically analyze pH-induced folding. Here, we constructed a statistical mechanical model that incorporates the protonation mechanism of conventional models along with a combined manual search and least-squares fitting procedure, which was used to investigate the folding of horse apomyoglobin over a wide pH range (2.2-6.7), with a time window ranging from ∼40 μs to ∼100 s, using continuous-/stopped-flow fluorescence at 8 °C. Quantitative analysis assuming a five-state sequential scheme indicated that (1) pH-induced folding/unfolding is represented by both specific binding and Coulombic interactions; (2) kinetic folding/unfolding intermediates share kinetic mechanisms with the equilibrium intermediate, indicating their equivalence; and (3) native-like properties are acquired successively during folding by intermediates and in transition states. This model could also be applied to a variety of association/dissociation processes.

Residual structures in the unfolded state of starch-binding domain of glucoamylase revealed by near-UV circular dichroism and protein engineering techniques

Protein folding is a thermodynamic process driven by energy gaps between the native and unfolded states. Although a wealth of information is available on the structure of folded species, there is a paucity of data on unfolded species. Here, we analyzed the structural properties of the unfolded state of the starch-binding domain of glucoamylase from Aspergillus niger (SBD) formed in the presence of guanidinium hydrochloride (GuHCl). Although far-UV CD and intrinsic tryptophan fluorescence spectra as well as small angle X-ray scattering suggested that SBD assumes highly unfolded structures in the presence of GuHCl, near-UV circular dichroism of wild-type SBD suggested the presence of residual structures in the unfolded state. Analyses of the unfolded states of tryptophan mutants (W543L, W563A, W590A and W615L) using Similarity Parameter, a modified version of root mean square deviation as a measure of similarity between two spectra, suggested that W543 and W563 have preferences to form native-like residual structures in the GuHCl-unfolded state. In contrast, W615 was entirely unstructured, while W590 tended to form non-native ordered structures in the unfolded state. These data and the amino acid sequence of SBD suggest that local structural propensities in the unfolded state can be determined by the probability of the presence of hydrophobic or charged residues nearby tryptophan residues.

Ensemble models of proteins and protein domains based on distance distribution restraints

Conformational ensembles of intrinsically disordered peptide chains are not fully determined by experimental observations. Uncertainty due to lack of experimental restraints and due to intrinsic disorder can be distinguished if distance distributions restraints are available. Such restraints can be obtained from pulsed dipolar electron paramagnetic resonance (EPR) spectroscopy applied to pairs of spin labels. Here, we introduce a Monte Carlo approach for generating conformational ensembles that are consistent with a set of distance distribution restraints, backbone dihedral angle statistics in known protein structures, and optionally, secondary structure propensities or membrane immersion depths. The approach is tested with simulated restraints for a terminal and an internal loop and for a protein with 69 residues by using sets of sparse restraints for underlying well-defined conformations and for published ensembles of a premolten globule-like and a coil-like intrinsically disordered protein.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

Mechanisms of amyloid formation revealed by solution NMR

Amyloid fibrils are proteinaceous elongated aggregates involved in more than fifty human diseases. Recent advances in electron microscopy and solid state NMR have allowed the characterization of fibril structures to different extents of refinement. However, structural details about the mechanism of fibril formation remain relatively poorly defined. This is mainly due to the complex, heterogeneous and transient nature of the species responsible for assembly; properties that make them difficult to detect and characterize in structural detail using biophysical techniques. The ability of solution NMR spectroscopy to investigate exchange between multiple protein states, to characterize transient and low-population species, and to study high molecular weight assemblies, render NMR an invaluable technique for studies of amyloid assembly. In this article we review state-of-the-art solution NMR methods for investigations of: (a) protein dynamics that lead to the formation of aggregation-prone species; (b) amyloidogenic intrinsically disordered proteins; and (c) protein-protein interactions on pathway to fibril formation. Together, these topics highlight the power and potential of NMR to provide atomic level information about the molecular mechanisms of one of the most fascinating problems in structural biology.

Apoglobin Stability Is the Major Factor Governing both Cell-free and in Vivo Expression of Holomyoglobin

Expression levels in animal muscle tissues and in Escherichia coli vary widely for naturally occurring mammalian myoglobins (Mb). To explore this variation, we developed an in vitro transcription and wheat germ extract-based translation assay to examine quantitatively the factors that govern expression of holoMb. We constructed a library of naturally occurring Mbs from two terrestrial and four deep-diving aquatic mammals and three distal histidine mutants designed to enhance apoglobin stability but decrease hemin affinity. A strong linear correlation is observed between cell-free expression levels of holo-metMb variants and their corresponding apoglobin stabilities, which were measured independently by guanidine HCl-induced unfolding titrations using purified proteins. In contrast, there is little dependence of expression on hemin affinity. Our results confirm quantitatively that deep diving mammals have highly stable Mbs that express to higher levels in animal myocytes, E. coli, and the wheat germ cell-free system than Mbs from terrestrial mammals. Our theoretical analyses show that the rate of aggregation of unfolded apoMb is very large, and as a result, the key factor for high level expression of holoMb, and presumably other heme proteins, is an ultra high fraction of folded, native apoglobin that is capable of rapidly binding hemin. This fraction is determined by the overall equilibrium folding constant and not hemin affinity. These results also demonstrate that the cell-free transcription/translation system can be used as a high throughput platform to screen for apoglobin stability without the need to generate large amounts of protein for in vitro unfolding measurements.

A hypothesis to reconcile the physical and chemical unfolding of proteins

High pressure (HP) or urea is commonly used to disturb folding species. Pressure favors the reversible unfolding of proteins by causing changes in the volumetric properties of the protein-solvent system. However, no mechanistic model has fully elucidated the effects of urea on structure unfolding, even though protein-urea interactions are considered to be crucial. Here, we provide NMR spectroscopy and 3D reconstructions from X-ray scattering to develop the "push-and-pull" hypothesis, which helps to explain the initial mechanism of chemical unfolding in light of the physical events triggered by HP. In studying MpNep2 from Moniliophthora perniciosa, we tracked two cooperative units using HP-NMR as MpNep2 moved uphill in the energy landscape; this process contrasts with the overall structural unfolding that occurs upon reaching a threshold concentration of urea. At subdenaturing concentrations of urea, we were able to trap a state in which urea is preferentially bound to the protein (as determined by NMR intensities and chemical shifts); this state is still folded and not additionally exposed to solvent [fluorescence and small-angle X-ray scattering (SAXS)]. This state has a higher susceptibility to pressure denaturation (lower p1/2 and larger ΔVu); thus, urea and HP share concomitant effects of urea binding and pulling and water-inducing pushing, respectively. These observations explain the differences between the molecular mechanisms that control the physical and chemical unfolding of proteins, thus opening up new possibilities for the study of protein folding and providing an interpretation of the nature of cooperativity in the folding and unfolding processes.

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering, Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.

Structure and dynamics of the MKK7-JNK signaling complex

Signaling specificity in the mitogen-activated protein kinase (MAPK) pathways is controlled by disordered domains of the MAPK kinases (MKKs) that specifically bind to their cognate MAPKs via linear docking motifs. MKK7 activates the c-Jun N-terminal kinase (JNK) pathway and is the only MKK containing three motifs within its regulatory domain. Here, we characterize the conformational behavior and interaction mechanism of the MKK7 regulatory domain. Using NMR spectroscopy, we develop an atomic resolution ensemble description of MKK7, revealing highly diverse intrinsic conformational propensities of the three docking sites, suggesting that prerecognition sampling of the bound-state conformation is not prerequisite for binding. Although the different sites exhibit similar affinities for JNK1, interaction kinetics differ considerably. Importantly, we determine the crystal structure of JNK1 in complex with the second docking site of MKK7, revealing two different binding modes of the docking motif correlating with observations from NMR exchange spectroscopy. Our results provide unique insight into how signaling specificity is regulated by linear motifs and, in general, into the role of conformational disorder in MAPK signaling.

NMR Spectroscopic Studies of the Conformational Ensembles of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are characterized by substantial conformational flexibility and thus not amenable to conventional structural biology techniques. Given their inherent structural flexibility NMR spectroscopy offers unique opportunities for structural and dynamic studies of IDPs. The past two decades have witnessed significant development of NMR spectroscopy that couples advances in spin physics and chemistry with a broad range of applications. This chapter will summarize key advances in NMR methodology. Despite the availability of efficient (multi-dimensional) NMR experiments for signal assignment of IDPs it is discussed that NMR of larger and more complex IDPs demands spectral simplification strategies capitalizing on specific isotope-labeling strategies. Prototypical applications of isotope labeling-strategies are described. Since IDP-ligand association and dissociation processes frequently occur on time scales that are amenable to NMR spectroscopy we describe in detail the application of CPMG relaxation dispersion techniques to studies of IDP protein binding. Finally, we demonstrate that the complementary usage of NMR and EPR data provide a more comprehensive picture about the conformational states of IDPs and can be employed to analyze the conformational ensembles of IDPs.

A relationship between the transient structure in the monomeric state and the aggregation propensities of α-synuclein and β-synuclein

α-Synuclein is an intrinsically disordered protein whose aggregation is implicated in Parkinson's disease. A second member of the synuclein family, β-synuclein, shares significant sequence similarity with α-synuclein but is much more resistant to aggregation. β-Synuclein is missing an 11-residue stretch in the central non-β-amyloid component region that forms the core of α-synuclein amyloid fibrils, yet insertion of these residues into β-synuclein to produce the βSHC construct does not markedly increase the aggregation propensity. To investigate the structural basis of these different behaviors, quantitative nuclear magnetic resonance data, in the form of paramagnetic relaxation enhancement-derived interatomic distances, are combined with molecular dynamics simulations to generate ensembles of structures representative of the solution states of α-synuclein, β-synuclein, and βSHC. Comparison of these ensembles reveals that the differing aggregation propensities of α-synuclein and β-synuclein are associated with differences in the degree of residual structure in the C-terminus coupled to the shorter separation between the N- and C-termini in β-synuclein and βSHC, making protective intramolecular contacts more likely.

p15PAF is an intrinsically disordered protein with nonrandom structural preferences at sites of interaction with other proteins

We present to our knowledge the first structural characterization of the proliferating-cell-nuclear-antigen-associated factor p15(PAF), showing that it is monomeric and intrinsically disordered in solution but has nonrandom conformational preferences at sites of protein-protein interactions. p15(PAF) is a 12 kDa nuclear protein that acts as a regulator of DNA repair during DNA replication. The p15(PAF) gene is overexpressed in several types of human cancer. The nearly complete NMR backbone assignment of p15(PAF) allowed us to measure 86 N-H(N) residual dipolar couplings. Our residual dipolar coupling analysis reveals nonrandom conformational preferences in distinct regions, including the proliferating-cell-nuclear-antigen-interacting protein motif (PIP-box) and the KEN-box (recognized by the ubiquitin ligase that targets p15(PAF) for degradation). In accordance with these findings, analysis of the (15)N R2 relaxation rates shows a relatively reduced mobility for the residues in these regions. The agreement between the experimental small angle x-ray scattering curve of p15(PAF) and that computed from a statistical coil ensemble corrected for the presence of local secondary structural elements further validates our structural model for p15(PAF). The coincidence of these transiently structured regions with protein-protein interaction and posttranslational modification sites suggests a possible role for these structures as molecular recognition elements for p15(PAF).

NMR contributions to structural dynamics studies of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are characterized by substantial conformational plasticity. Given their inherent structural flexibility X-ray crystallography is not applicable to study these proteins. In contrast, NMR spectroscopy offers unique opportunities for structural and dynamic studies of IDPs. The past two decades have witnessed significant development of NMR spectroscopy that couples advances in spin physics and chemistry with a broad range of applications. This article will summarize key advances in basic physical-chemistry and NMR methodology, outline their limitations and envision future R&D directions.

NMR characterization of the near native and unfolded states of the PTB domain of Dok1: alternate conformations and residual clusters

BACKGROUND

Phosphotyrosine binding (PTB) domains are critically involved in cellular signaling and diseases. PTB domains are categorized into three distinct structural classes namely IRS-like, Shc-like and Dab-like. All PTB domains consist of a core pleckstrin homology (PH) domain with additional structural elements in Shc and Dab groups. The core PH fold of the PTB domain contains a seven stranded β-sheet and a long C-terminal helix.

PRINCIPAL FINDINGS

In this work, the PTB domain of Dok1 protein has been characterized, by use of NMR spectroscopy, in solutions containing sub-denaturing and denaturing concentrations of urea. We find that the Dok1 PTB domain displays, at sub-denaturing concentrations of urea, alternate conformational states for residues located in the C-terminal helix and in the β5 strand of the β-sheet region. The β5 strand of PTB domain has been found to be experiencing significant chemical shift perturbations in the presence of urea. Notably, many of these residues in the helix and in the β5 strand are also involved in ligand binding. Structural and dynamical analyses at 7 M urea showed that the PTB domain is unfolded with islands of motionally restricted regions in the polypeptide chain. Further, the C-terminal helix appears to be persisted in the unfolded state of the PTB domain. By contrast, residues encompassing β-sheets, loops, and the short N-terminal helix lack any preferred secondary structures. Moreover, these residues demonstrated an intimate contact with the denaturant.

SIGNIFICANCE

This study implicates existence of alternate conformational states around the ligand binding pocket of the PTB domain either in the native or in the near native conditions. Further, the current study demonstrates that the C-terminal helical region of PTB domain may be considered as a potential site for the initiation of folding.

Structural characterization of intrinsically disordered proteins by NMR spectroscopy

Recent advances in NMR methodology and techniques allow the structural investigation of biomolecules of increasing size with atomic resolution. NMR spectroscopy is especially well-suited for the study of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) which are in general highly flexible and do not have a well-defined secondary or tertiary structure under functional conditions. In the last decade, the important role of IDPs in many essential cellular processes has become more evident as the lack of a stable tertiary structure of many protagonists in signal transduction, transcription regulation and cell-cycle regulation has been discovered. The growing demand for structural data of IDPs required the development and adaption of methods such as ¹³C-direct detected experiments, paramagnetic relaxation enhancements (PREs) or residual dipolar couplings (RDCs) for the study of ‘unstructured’ molecules in vitro and in-cell. The information obtained by NMR can be processed with novel computational tools to generate conformational ensembles that visualize the conformations IDPs sample under functional conditions. Here, we address NMR experiments and strategies that enable the generation of detailed structural models of IDPs.

Residual dipolar couplings measured in unfolded proteins are sensitive to amino-acid-specific geometries as well as local conformational sampling

Many functional proteins do not have well defined folded structures. In recent years, both experimental and computational approaches have been developed to study the conformational behaviour of this type of protein. It has been shown previously that experimental RDCs (residual dipolar couplings) can be used to study the backbone sampling of disordered proteins in some detail. In these studies, the backbone structure was modelled using a common geometry for all amino acids. In the present paper, we demonstrate that experimental RDCs are also sensitive to the specific geometry of each amino acid as defined by energy-minimized internal co-ordinates. We have modified the FM (flexible-Meccano) algorithm that constructs conformational ensembles on the basis of a statistical coil model, to account for these differences. The modified algorithm inherits the advantages of the FM algorithm to efficiently sample the potential energy landscape for coil conformations. The specific geometries incorporated in the new algorithm result in a better reproduction of experimental RDCs and are generally applicable for further studies to characterize the conformational properties of intrinsically disordered proteins. In addition, the internal-co-ordinate-based algorithm is an order of magnitude more efficient, and facilitates side-chain construction, surface osmolyte simulation, spin-label distribution sampling and proline cis/trans isomer simulation.

Structural characterization of intrinsically disordered proteins by the combined use of NMR and SAXS

In recent years, IDPs (intrinsically disordered proteins) have emerged as pivotal actors in biology. Despite IDPs being present in all kingdoms of life, they are more abundant in eukaryotes where they are involved in the vast majority of regulation and signalling processes. The realization that, in some cases, functional states of proteins were partly or fully disordered was in contradiction to the traditional view where a well defined three-dimensional structure was required for activity. Several experimental evidences indicate, however, that structural features in IDPs such as transient secondary-structural elements and overall dimensions are crucial to their function. NMR has been the main tool to study IDP structure by probing conformational preferences at residue level. Additionally, SAXS (small-angle X-ray scattering) has the capacity to report on the three-dimensional space sampled by disordered states and therefore complements the local information provided by NMR. The present review describes how the synergy between NMR and SAXS can be exploited to obtain more detailed structural and dynamic models of IDPs in solution. These combined strategies, embedded into computational approaches, promise the elucidation of the structure-function properties of this important, but elusive, family of biomolecules.

Flexible-meccano: a tool for the generation of explicit ensemble descriptions of intrinsically disordered proteins and their associated experimental observables

MOTIVATION

Intrinsically disordered proteins (IDPs) represent a significant fraction of the human proteome. The classical structure function paradigm that has successfully underpinned our understanding of molecular biology breaks down when considering proteins that have no stable tertiary structure in their functional form. One convenient approach is to describe the protein in terms of an equilibrium of rapidly inter-converting conformers. Currently, tools to generate such ensemble descriptions are extremely rare, and poorly adapted to the prediction of experimental data.

RESULTS

We present flexible-meccano-a highly efficient algorithm that generates ensembles of molecules, on the basis of amino acid-specific conformational potentials and volume exclusion. Conformational sampling depends uniquely on the primary sequence, with the possibility of introducing additional local or long-range conformational propensities at an amino acid-specific resolution. The algorithm can also be used to calculate expected values of experimental parameters measured at atomic or molecular resolution, such as nuclear magnetic resonance (NMR) and small angle scattering, respectively. We envisage that flexible-meccano will be useful for researchers who wish to compare experimental data with those expected from a fully disordered protein, researchers who see experimental evidence of deviation from 'random coil' behaviour in their protein, or researchers who are interested in working with a broad ensemble of conformers representing the flexibility of the IDP of interest.

AVAILABILITY

A fully documented multi-platform executable is provided, with examples, at http://www.ibs.fr/science-213/scientific-output/software/flexible-meccano/

CONTACT

martin.blackledge@ibs.fr.

Experiments and simulations show how long-range contacts can form in expanded unfolded proteins with negligible secondary structure

The sizes of unfolded proteins under highly denaturing conditions scale as N(0.59) with chain length. This suggests that denaturing conditions mimic good solvents, whereby the preference for favorable chain-solvent interactions causes intrachain interactions to be repulsive, on average. Beyond this generic inference, the broader implications of N(0.59) scaling for quantitative descriptions of denatured state ensembles (DSEs) remain unresolved. Of particular interest is the degree to which N(0.59) scaling can simultaneously accommodate intrachain attractions and detectable long-range contacts. Here we present data showing that the DSE of the N-terminal domain of the L9 (NTL9) ribosomal protein in 8.3 M urea lacks detectable secondary structure and forms expanded conformations in accord with the expected N(0.59) scaling behavior. Paramagnetic relaxation enhancements, however, indicate the presence of detectable long-range contacts in the denatured-state ensemble of NTL9. To explain these observations we used atomistic thermal unfolding simulations to identify ensembles whose properties are consistent with all of the experimental observations, thus serving as useful proxies for the DSE of NTL9 in 8.3 M urea. Analysis of these ensembles shows that residual attractions are present under mimics of good solvent conditions, and for NTL9 they result from low-likelihood, medium/long-range contacts between hydrophobic residues. Our analysis provides a quantitative framework for the simultaneous observation of N(0.59) scaling and low-likelihood long-range contacts for the DSE of NTL9. We propose that such low-likelihood intramolecular hydrophobic clusters might be a generic feature of DSEs that play a gatekeeping role to protect against aggregation during protein folding.

Investigation of the Polymeric Properties of α-Synuclein and Comparison with NMR Experiments: A Replica Exchange Molecular Dynamics Study

Intrinsically disordered proteins (IDPs) have been shown to be involved in a number of cellular functions, in addition to their predominance in diseased states. α-synuclein may be described as one such IDP implicated in the pathology of Parkinson's disease. Understanding the conformational characteristics of the monomeric state of α-synuclein is necessary for understanding the role of the monomer conformation in aggregation. Polymer theories have been applied to investigate the statistical properties of homopolymeric IDPs. Here we use Replica Exchange Molecular Dynamics (REMD) simulations using temperature as a proxy for solvent quality to examine how well these theories developed for homopolymeric chains describe heteropolymeric α-synuclein. Our results indicate that α-synuclein behaves like a homopolymer at the extremes of solvent quality, while in the intermediate solvent regime, the uneven distribution of charged residues along the sequence strongly influences the conformations adopted by the chain. We refine the ensemble extracted from the REMD simulations of α-synuclein, which shows the best qualitative agreement with experiment, by fitting to the experimental NMR Residual Dipolar Couplings (RDCs) and Paramagnetic Relaxation Enhancements (PREs). Our results demonstrate that the detailed shape of the RDC patterns are sensitive to the angular correlations that are local in sequence while longer range anti-correlations which arise from packing constraints affect the RDC magnitudes.

The interplay between transient α-helix formation and side chain rotamer distributions in disordered proteins probed by methyl chemical shifts

The peptide backbones of disordered proteins are routinely characterized by NMR with respect to transient structure and dynamics. Little experimental information is, however, available about the side chain conformations and how structure in the backbone affects the side chains. Methyl chemical shifts can in principle report the conformations of aliphatic side chains in disordered proteins and in order to examine this two model systems were chosen: the acid denatured state of acyl-CoA binding protein (ACBP) and the intrinsically disordered activation domain of the activator for thyroid hormone and retinoid receptors (ACTR). We find that small differences in the methyl carbon chemical shifts due to the γ-gauche effect may provide information about the side chain rotamer distributions. However, the effects of neighboring residues on the methyl group chemical shifts obscure the direct observation of γ-gauche effect. To overcome this, we reference the chemical shifts to those in a more disordered state resulting in residue specific random coil chemical shifts. The (13)C secondary chemical shifts of the methyl groups of valine, leucine, and isoleucine show sequence specific effects, which allow a quantitative analysis of the ensemble of χ(2)-angles of especially leucine residues in disordered proteins. The changes in the rotamer distributions upon denaturation correlate to the changes upon helix induction by the co-solvent trifluoroethanol, suggesting that the side chain conformers are directly or indirectly related to formation of transient α-helices.