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.

Engineered solubility tag for solution NMR of proteins

The low solubility of many proteins hinders large scale expression and purification as well as biophysical measurements. Here, we devised a general strategy to solubilize a protein by conjugating it at a solvent-exposed position to a 6 kDa protein that was re-engineered to be highly soluble. We applied this method to the CARD domain of Apoptosis-associated speck-like protein containing a CARD (ASC), which represents one member of a class of proteins that are notoriously prone to aggregation. Attachment of the tag to a cysteine residue, introduced by site-directed mutagenesis at its self-association interface, improved the solubility of the ASC CARD over 50-fold under physiological conditions. Although it is not possible to use nuclear magnetic resonance (NMR) to obtain a high quality 2D correlation spectrum of the wild type domain under physiological conditions, we demonstrate that NMR relaxation parameters of the solubilized variant are sufficiently improved to facilitate virtually any demanding measurement. The method shown here represents a straightforward approach for dramatically increasing protein solubility, enabled by ease of labeling as well as flexibility in tag placement with minimal perturbation to the target.

Experimental determination of microsecond reorientation correlation times in protein solutions

Reorientation correlation times in protein solutions are key determinants for feasibility and quality of NMR experiments. Yet, their accurate estimate is not easy, especially in the case of very large proteins. We show that nuclear magnetic relaxation dispersion (NMRD) can accurately determine reorientation times up to the microsecond range. A theoretical description for the analysis of the NMRD profiles is provided, and the protein reorientation time is shown to be provided by the longest correlation time among those needed to reproduce the experimental profile. Measurements are performed using samples of the archaeal proteasome double ring α7α7 and of αB-Crystallin in glycerol solutions.

Cross-validation of the structure of a transiently formed and low populated FF domain folding intermediate determined by relaxation dispersion NMR and CS-Rosetta

We have recently reported the atomic resolution structure of a low populated and transiently formed on-pathway folding intermediate of the FF domain from human HYPA/FBP11 [Korzhnev, D. M.; Religa, T. L.; Banachewicz, W.; Fersht, A. R.; Kay, L.E. Science 2011, 329, 1312-1316]. The structure was determined on the basis of backbone chemical shift and bond vector orientation restraints of the invisible intermediate state measured using relaxation dispersion nuclear magnetic resonance (NMR) spectroscopy that were subsequently input into the database structure determination program, CS-Rosetta. As a cross-validation of the structure so produced, we present here the solution structure of a mimic of the folding intermediate that is highly populated in solution, obtained from the wild-type domain by mutagenesis that destabilizes the native state. The relaxation dispersion/CS-Rosetta structures of the intermediate are within 2 Å of those of the mimic, with the nonnative interactions in the intermediate also observed in the mimic. This strongly confirms the structure of the FF domain folding intermediate, in particular, and validates the use of relaxation dispersion derived restraints in structural studies of invisible excited states, in general.

Nonnative interactions in the FF domain folding pathway from an atomic resolution structure of a sparsely populated intermediate: an NMR relaxation dispersion study

Several all-helical single-domain proteins have been shown to fold rapidly (microsecond time scale) to a compact intermediate state and subsequently rearrange more slowly to the native conformation. An understanding of this process has been hindered by difficulties in experimental studies of intermediates in cases where they are both low-populated and only transiently formed. One such example is provided by the on-pathway folding intermediate of the small four-helix bundle FF domain from HYPA/FBP11 that is populated at several percent with a millisecond lifetime at room temperature. Here we have studied the L24A mutant that has been shown previously to form nonnative interactions in the folding transition state. A suite of Carr-Purcell-Meiboom-Gill relaxation dispersion NMR experiments have been used to measure backbone chemical shifts and amide bond vector orientations of the invisible folding intermediate that form the input restraints in calculations of atomic resolution models of its structure. Despite the fact that the intermediate structure has many features that are similar to that of the native state, a set of nonnative contacts is observed that is even more extensive than noted previously for the wild-type (WT) folding intermediate. Such nonnative interactions, which must be broken prior to adoption of the native conformation, explain why the transition from the intermediate state to the native conformer (millisecond time scale) is significantly slower than from the unfolded ensemble to the intermediate and why the L24A mutant folds more slowly than the WT.

13CHD2 methyl group probes of millisecond time scale exchange in proteins by 1H relaxation dispersion: an application to proteasome gating residue dynamics

A pulse scheme is presented for quantifying millisecond time scale chemical exchange processes in proteins by measuring (1)H CPMG relaxation dispersion profiles of (13)CHD(2) methyl groups. The use of (13)CHD(2) isotopomers for (1)H methyl dispersion experiments eliminates problems with interconversion between differentially relaxing proton transitions that complicate the extraction of accurate exchange parameters when (13)CH(3) probes are used. Good agreement is demonstrated between extracted chemical shift differences from fits of dispersion profiles and the corresponding differences measured independently on a model exchanging system, validating the experiment. The methodology is applied to the gating residues of the T. acidiphilium proteasome that are shown to undergo extensive motion on the millisecond time scale.

A Transient and Low-Populated Protein-Folding Intermediate at Atomic Resolution

Proteins can sample conformational states that are critical for function but are seldom detected directly because of their low occupancies and short lifetimes. In this work, we used chemical shifts and bond-vector orientation constraints obtained from nuclear magnetic resonance relaxation dispersion spectroscopy, in concert with a chemical shift–based method for structure elucidation, to determine an atomic-resolution structure of an “invisible” folding intermediate of a small protein module: the FF domain. The structure reveals non-native elements preventing formation of the native conformation in the carboxyl-terminal part of the protein. This is consistent with the kinetics of folding in which a well-structured intermediate forms rapidly and then rearranges slowly to the native state. The approach introduces a general strategy for structure determination of low-populated and transiently formed protein states.

Optimal methyl labeling for studies of supra-molecular systems

Selective methyl labeling combined with HMQC spectroscopy that exploits a TROSY effect in (13)CH(3) spin systems has significantly extended the utility of solution NMR spectroscopy in studies of high molecular weight particles. Herein we compare the utility of (13)CH(3)- versus (13)CHD(2)-labeling of Ile, Leu, Val probes in supra-molecular systems through quantification of relative signal-to-noise ratios in optimized spectra of highly deuterated, (13)CH(3)- and (13)CHD(2)-labeled samples of the half proteasome (alpha(7)alpha(7), 360 kDa). It is shown that the sensitivity of spectra recorded on Ile, Leu, Val (13)CH(3)-labeled samples is between 1.5 and 2 fold higher than the corresponding data sets obtained on alpha(7)alpha(7) with (13)CHD(2) probes. Thus, labeling of supra-molecules with (13)CH(3) isotopomers remains the method of choice, but in applications where (13)CHD(2) moieties are required, sensitivity will in general not be limiting.

Effect of noncompetitive proteasome inhibition on bortezomib resistance

BACKGROUND

Bortezomib and the other proteasome inhibitors that are currently under clinical investigation bind to the catalytic sites of proteasomes and are competitive inhibitors. We hypothesized that proteasome inhibitors that act through a noncompetitive mechanism might overcome some forms of bortezomib resistance.

METHODS

5-amino-8-hydroxyquinoline (5AHQ) was identified through a screen of a 27-compound chemical library based on the quinoline pharmacophore to identify proteasome inhibitors. Inhibition of proteasome activity by 5AHQ was tested by measuring 7-amino-4-methylcoumarin (AMC) release from the proteasome substrate Suc-LLVY-AMC in intact human and mouse leukemia and myeloma cells and in tumor cell protein extracts. Cytotoxicity was assessed in 5AHQ-treated cell lines and primary cells from myeloma and leukemia patients using AlamarBlue fluorescence and MTS assays, trypan blue staining, and annexin V staining. 5AHQ-proteasome interaction was assessed by nuclear magnetic resonance. 5AHQ efficacy was evaluated in three leukemia xenograft mouse models (9-10 mice per group per model). All statistical tests were two-sided.

RESULTS

5AHQ inhibited the proteasome when added to cell extracts and intact cells (the mean concentration inhibiting 50% [IC(50)] of AMC release in intact cells ranged from 0.57 to 5.03 microM), induced cell death in intact cells from leukemia and myeloma cell lines (mean IC(50) values for cell growth ranged from 0.94 to 3.85 microM), and preferentially induced cell death in primary myeloma and leukemia cells compared with normal hematopoietic cells. 5AHQ was equally cytotoxic to human myelomonocytic THP1 cells and to THP1/BTZ500 cells, which are 237-fold more resistant to bortezomib than wild-type THP1 cells because of their overexpression and mutation of the bortezomib-binding beta5 proteasome subunit (mean IC(50) for cell death in the absence of bortezomib, wild-type THP1: 3.7 microM, 95% confidence interval = 3.4 to 4.0 microM; THP1/BTZ500: 6.6 microM, 95% confidence interval = 5.9 to 7.5 microM). 5AHQ interacted with the alpha subunits of the 20S proteasome at noncatalytic sites. Orally administered 5AHQ inhibited tumor growth in all three mouse models of leukemia without overt toxicity (eg, OCI-AML2 model, median tumor weight [interquartile range], 5AHQ vs control: 95.7 mg [61.4-163.5 mg] vs 247.2 mg [189.4-296.2 mg], P = .002).

CONCLUSIONS

5AHQ is a noncompetitive proteasome inhibitor that is cytotoxic to myeloma and leukemia cells in vitro and inhibits xenograft tumor growth in vivo. 5AHQ can overcome some forms of bortezomib resistance in vitro.

Dynamic regulation of archaeal proteasome gate opening as studied by TROSY NMR

The proteasome catalyzes the majority of protein degradation in the cell and plays an integral role in cellular homeostasis. Control over proteolysis by the 20S core-particle (CP) proteasome is achieved by gated access of substrate; thus, an understanding of the molecular mechanism by which these gates regulate substrate entry is critical. We used methyl-transverse relaxation optimized nuclear magnetic resonance spectroscopy to show that the amino-terminal residues that compose the gates of the alpha subunits of the Thermoplasma acidophilum proteasome are highly dynamic over a broad spectrum of time scales and that gating termini are in conformations that extend either well inside (closed gate) or outside (open gate) of the antechamber. Interconversion between these conformers on a time scale of seconds leads to a dynamic regulation of 20S CP proteolysis activity.

Comparison of multiple crystal structures with NMR data for engrailed homeodomain

Two methods are currently available to solve high resolution protein structures-X-ray crystallography and nuclear magnetic resonance (NMR). Both methods usually produce highly similar structures, but small differences between both solutions are always observed. Here the raw NMR data as well as the solved NMR structure were compared to the multiple crystal structures solved for the WT 60 residue three helix bundle engrailed homeodomain (EnHD) and single point mutants. There was excellent agreement between TALOS-predicted and crystal structure-observed dihedral angles and a good agreement for the (3) J(H ( N ) H ( alpha )) couplings for the multiple crystal structures. Around 1% of NOEs were violated for any crystal structure, but no NOE was inconsistent with all of the crystal structures. Violations usually occurred for surface residues or for residues for which multiple discreet conformations were observed between the crystal structures. Comparison of the disorder shown in the multiple crystal structures shows little correlation with dynamics under native conditions for this protein.

Fractional 13C enrichment of isolated carbons using [1-13C]- or [2- 13C]-glucose facilitates the accurate measurement of dynamics at backbone Calpha and side-chain methyl positions in proteins

A simple labeling approach is presented based on protein expression in [1-(13)C]- or [2-(13)C]-glucose containing media that produces molecules enriched at methyl carbon positions or backbone C(alpha) sites, respectively. All of the methyl groups, with the exception of Thr and Ile(delta1) are produced with isolated (13)C spins (i.e., no (13)C-(13)C one bond couplings), facilitating studies of dynamics through the use of spin-spin relaxation experiments without artifacts introduced by evolution due to large homonuclear scalar couplings. Carbon-alpha sites are labeled without concomitant labeling at C(beta) positions for 17 of the common 20 amino acids and there are no cases for which (13)C(alpha)-(13)CO spin pairs are observed. A large number of probes are thus available for the study of protein dynamics with the results obtained complimenting those from more traditional backbone (15)N studies. The utility of the labeling is established by recording (13)C R (1rho) and CPMG-based experiments on a number of different protein systems.

Solution structure of the BRK domains from CHD7

CHD7 is a member of the chromodomain helicase DNA binding domain (CHD) family of ATP-dependent chromatin remodelling enzymes. It is mutated in CHARGE syndrome, a multiple congenital anomaly condition. CHD7 is one of a subset of CHD proteins, unique to metazoans that contain the BRK domain, a protein module also found in the Brahma/BRG1 family of helicases. We describe here the NMR solution structure of the two BRK domains of CHD7. Each domain has a compact betabetaalphabeta fold. The second domain has a C-terminal extension consisting of two additional helices. The structure differs from those of other domains present in chromatin-associated proteins.

The folding pathway of an FF domain: characterization of an on-pathway intermediate state under folding conditions by (15)N, (13)C(alpha) and (13)C-methyl relaxation dispersion and (1)H/(2)H-exchange NMR spectroscopy

The FF domain from the human protein HYPA/FBP11 folds via a low-energy on-pathway intermediate (I). Elucidation of the structure of such folding intermediates and denatured states under conditions that favour folding are difficult tasks. Here, we investigated the millisecond time-scale equilibrium folding transition of the 71-residue four-helix bundle wild-type protein by (15)N, (13)C(alpha) and methyl(13)C Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiments and by (1)H/(2)H-exchange measurements. The relaxation data for the wild-type protein fitted a simple two-site exchange process between the folded state (F) and I. Destabilization of F in mutants A17G and Q19G allowed the detection of the unfolded state U by (15)N CPMG relaxation dispersion. The dispersion data for these mutants fitted a three-site exchange scheme, U<-->I<-->F, with I populated higher than U. The kinetics and thermodynamics of the folding reaction were obtained via temperature and urea-dependent relaxation dispersion experiments, along with structural information on I from backbone (15)N, (13)C(alpha) and side-chain methyl (13)C chemical shifts, with further information from protection factors for the backbone amide groups from (1)H/(2)H-exchange. Notably, helices H1-H3 are at least partially formed in I, while helix H4 is largely disordered. Chemical shift differences for the methyl (13)C nuclei suggest a paucity of stable, native-like hydrophobic interactions in I. These data are consistent with Phi-analysis of the rate-limiting transition state between I and F. The combination of relaxation dispersion and Phi data can elucidate whole experimental folding pathways.

The helix-turn-helix motif as an ultrafast independently folding domain: the pathway of folding of Engrailed homeodomain

Helices 2 and 3 of Engrailed homeodomain (EnHD) form a helix-turn-helix (HTH) motif. This common motif is believed not to fold independently, which is the characteristic feature of a motif rather than a domain. But we found that the EnHD HTH motif is monomeric and folded in solution, having essentially the same structure as in full-length protein. It had a sigmoidal thermal denaturation transition. Both native backbone and local tertiary interactions were formed concurrently at 4 x 10(5) s(-1) at 25 degrees C, monitored by IR and fluorescence T-jump kinetics, respectively, the same rate constant as for the fast phase in the folding of EnHD. The HTH motif, thus, is an ultrafast-folding, natural protein domain. Its independent stability and appropriate folding kinetics account for the stepwise folding of EnHD, satisfy fully the criteria for an on-pathway intermediate, and explain the changes in mechanism of folding across the homeodomain family. Experiments on mutated and engineered fragments of the parent protein with different probes allowed the assignment of the observed kinetic phases to specific events to show that EnHD is not an example of one-state downhill folding.

A single-quantum methyl 13C-relaxation dispersion experiment with improved sensitivity

A pulse sequence is described for recording single-quantum (13)C-methyl relaxation dispersion profiles of (13)C-selectively labeled methyl groups in proteins that offers significant improvements in sensitivity relative to existing approaches where initial magnetization derives from (13)C polarization. Sensitivity gains in the new experiment are achieved by making use of polarization from (1)H spins and (1)H --> (13)C --> (1)H type magnetization transfers. Its utility has been established by applications involving three different protein systems ranging in molecular weight from 8 to 28 kDa, produced using a number of different selective labeling approaches. In all cases exchange parameters from both (13)C-->(1)H and (1)H --> (13)C --> (1)H classes of experiment are in good agreement, with gains in sensitivity of between 1.7 and 4-fold realized using the new scheme.

Phi-analysis of the folding of the B domain of protein A using multiple optical probes

We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.

Solution structure of a protein denatured state and folding intermediate

The most controversial area in protein folding concerns its earliest stages. Questions such as whether there are genuine folding intermediates, and whether the events at the earliest stages are just rearrangements of the denatured state or progress from populated transition states, remain unresolved. The problem is that there is a lack of experimental high-resolution structural information about early folding intermediates and denatured states under conditions that favour folding because competent states spontaneously fold rapidly. Here we have solved directly the solution structure of a true denatured state by nuclear magnetic resonance under conditions that would normally favour folding, and directly studied its equilibrium and kinetic behaviour. We engineered a mutant of Drosophila melanogaster Engrailed homeodomain that folds and unfolds reversibly just by changing ionic strength. At high ionic strength, the mutant L16A is an ultra-fast folding native protein, just like the wild-type protein; however, at physiological ionic strength it is denatured. The denatured state is a well-ordered folding intermediate, poised to fold by docking helices and breaking some non-native interactions. It unfolds relatively progressively with increasingly denaturing conditions, and so superficially resembles a denatured state with properties that vary with conditions. Such ill-defined unfolding is a common feature of early folding intermediate states and accounts for why there are so many controversies about intermediates versus compact denatured states in protein folding.

Determinants of functionality in the ubiquitin conjugating enzyme family

The E2 enzymes are key enzymes in the ubiquitin and ubiquitin-like protein ligation pathways. To understand the functionality of the different E2 enzymes, we analyzed 190 protein sequences and 211 structures and electrostatic potentials. Key findings include: The ScUbc1 orthologs are defined by a C-terminal UBA domain. An N-terminal sequence motif that is highly conserved in all E2s except for Cdc34 orthologs is important for the stabilization of the L7 loop and is likely to be involved in E1 binding. ScUbc11p has a different electrostatic potential from E2-Cp and other proteins with which it has high sequence similarity but different functionality. All the E2s known to ubiquitinate histones have a negative potential. The members of the NCUBE family have a positive electrostatic potential, although its form is different from that of the SUMO conjugating E2s. The specificities of only the ScUbc4/Ubc5 and ScUbc1p orthologs are reflected in their L4 and L7 loops.

Testing protein-folding simulations by experiment: B domain of protein A

We have assessed the published predictions of the pathway of folding of the B domain of protein A, the pathway most studied by computer simulation. We analyzed the transition state for folding of the three-helix bundle protein, by using experimental Phi values on some 70 suitable mutants. Surprisingly, the third helix, which has the most stable alpha-helical structure as a peptide fragment, is poorly formed in the transition state, especially at its C terminus. The protein folds around a nearly fully formed central helix, which is stabilized by extensive hydrophobic side chain interactions. The turn connecting the poorly structured first helix to the central helix is unstructured, but the turn connecting the central helix to the third is in the process of being formed as the N-terminal region of the third helix begins to coalesce. The transition state is inconsistent with a classical framework mechanism and is closer to nucleation-condensation. None of the published atomistic simulations are fully consistent with the experimental picture although many capture important features. There is a continuing need for combining simulation with experiment to describe folding pathways, and of continued testing to improve predictive methods.