NMR analysis of cleaved Escherichia coli thioredoxin (1-73/74-108) and its P76A variant: cis/trans peptide isomerization

Enhanced sensitivity by nonuniform sampling enables multidimensional MAS NMR spectroscopy of protein assemblies

We report dramatic sensitivity enhancements in multidimensional MAS NMR spectra by the use of nonuniform sampling (NUS) and introduce maximum entropy interpolation (MINT) processing that assures the linearity between the time and frequency domains of the NUS acquired data sets. A systematic analysis of sensitivity and resolution in 2D and 3D NUS spectra reveals that with NUS, at least 1.5- to 2-fold sensitivity enhancement can be attained in each indirect dimension without compromising the spectral resolution. These enhancements are similar to or higher than those attained by the newest-generation commercial cryogenic probes. We explore the benefits of this NUS/MaxEnt approach in proteins and protein assemblies using 1-73-(U-(13)C,(15)N)/74-108-(U-(15)N) Escherichia coli thioredoxin reassembly. We demonstrate that in thioredoxin reassembly, NUS permits acquisition of high-quality 3D-NCACX spectra, which are inaccessible with conventional sampling due to prohibitively long experiment times. Of critical importance, issues that hinder NUS-based SNR enhancement in 3D-NMR of liquids are mitigated in the study of solid samples in which theoretical enhancements on the order of 3-4 fold are accessible by compounding the NUS-based SNR enhancement of each indirect dimension. NUS/MINT is anticipated to be widely applicable and advantageous for multidimensional heteronuclear MAS NMR spectroscopy of proteins, protein assemblies, and other biological systems.

Estimating the "steric clash" at cis peptide bonds

To account for the scarcity of cis peptide bonds in proteins, especially in nonproline (or secondary amide) cases, a steric-clash argument is often put forward, in a scheme where the R lateral chains are facing parallel one another, and the backbone is kept in an "all- trans"-like arrangement. Although such a steric conflict can be partly relieved through proper adjustment of the backbone dihedral angles, one can try to estimate its associated energy cost. To this end, quantum-chemistry approaches using a differential-torsion protocol and bond-separation-energy analyses are applied to N-ethyl propionamide CH3-CH2-CO-NH-CH2-CH3, regarded as a model capable of exhibiting C beta...C beta interaction as in alanine succession. The calculations provide an increment of 9 kcal/mol, quite close to that obtained in the nearly isostere (gsg) rotamer of n-hexane (10 kcal/mol), suggesting the local effects induced by methyl-methyl contact are similar in both cases. Analogous treatments on larger radicals as encountered in leucine or phenylalanine dimers do not change this increment much, which therefore defines the basic reference per-plaque quota to be overcome along all- cis chains. Explicit modeling indicated it can be reduced by up to a factor of 4.

Magic angle spinning NMR spectroscopy of thioredoxin reassemblies

Differentially isotopically enriched 1-73((13)C,(15)N)/74-108((15)N) and 1-73((15)N)/74-108((13)C,(15)N) Escherichia coli thioredoxin reassemblies prepared by fragment complementation were investigated by high-resolution magic angle spinning solid-state NMR spectroscopy. Nearly complete resonance assignments, secondary and tertiary structure analysis are reported for 1-73((13)C,(15)N)/74-108((15)N) reassembled thioredoxin. Temperature dependence of the dipolar-assisted rotational resonance (DARR) spectra reveals the residues undergoing intermediate timescale motions at temperatures below - 15 degrees C. Analysis of the DARR intensity buildups as a function of mixing time in these reassemblies indicates that at long mixing times medium- and long-range cross-peaks do not experience dipolar truncation, suggesting that isotopic dilution is not required for gaining nontrivial distance restraints for structure calculations.

On the NMR analysis of pKa values in the unfolded state of proteins by extrapolation to zero denaturant

Detailed knowledge of the pH-dependence in both folded and unfolded states of proteins is essential to understand the role of electrostatics in protein stability. The increasing number of natively disordered proteins constitutes an excellent source for the NMR analysis of pKa values in the unfolded state of proteins. However, the tendency of many natively disordered proteins to aggregate via intermolecular hydrophobic clusters limits their NMR analysis over a wide pH range. To assess whether the pKa values in natively disordered polypeptides can be extrapolated from NMR measurements in the presence of denaturants, the natively disordered backbone of the C-terminal fragment 75 to 105 of Human Thioredoxin was studied. First, assignments using triple resonance experiments were performed to confirm lack of secondary structure. Then the pH-dependence of the amides and carboxylate side chains of Glu residues (Glu88, Glu95, Glu98, and Glu103) in the pH range from 2.0 to 7.0 was monitored using 2D 1H15N HSQC and 3D C(CO)NH experiments, and the behavior of their amides and corresponding carboxyl groups was compared to confirm the absence of nonlocal interactions. Lastly, the effect of increasing dimethyl urea concentration on the pKa values of these Glu residues was monitored. The results indicate that: (i) the dispersion in the pKa of carboxyl groups and the pH midpoints of amides in Glu residues is about 0.5 pH units and 0.6 pH units, respectively; (ii) the backbone amides of the Glu residues exhibit pH midpoints which are within 0.2 pH units from those of their carboxylates; (iii) the addition of denaturant produces upshifts in the pKa values of Glu residues that are nearly independent of their position in the sequence; and (iv) these upshifts show a nonlinear behavior in denaturant concentration, complicating the extrapolation to zero denaturant. Nevertheless, the relative ordering of the pKa values of Glu residues is preserved over the whole range of denaturant concentrations indicating that measurements at high denaturant concentration (e.g. 4 M dimethyl urea) can yield a qualitatively correct ranking of the pKa of these residues in natively disordered proteins whose pH-dependence cannot be monitored directly by NMR.

Elucidating quantitative stability/flexibility relationships within thioredoxin and its fragments using a distance constraint model

Numerous quantitative stability/flexibility relationships, within Escherichia coli thioredoxin (Trx) and its fragments are determined using a minimal distance constraint model (DCM). A one-dimensional free energy landscape as a function of global flexibility reveals Trx to fold in a low-barrier two-state process, with a voluminous transition state. Near the folding transition temperature, the native free energy basin is markedly skewed to allow partial unfolded forms. Under native conditions the skewed shape is lost, and the protein forms a compact structure with some flexibility. Predictions on ten Trx fragments are generally consistent with experimental observations that they are disordered, and that complementary fragments reconstitute. A hierarchical unfolding pathway is uncovered using an exhaustive computational procedure of breaking interfacial cross-linking hydrogen bonds that span over a series of fragment dissociations. The unfolding pathway leads to a stable core structure (residues 22-90), predicted to act as a kinetic trap. Direct connection between degree of rigidity within molecular structure and non-additivity of free energy is demonstrated using a thermodynamic cycle involving fragments and their hierarchical unfolding pathway. Additionally, the model provides insight about molecular cooperativity within Trx in its native state, and about intermediate states populating the folding/unfolding pathways. Native state cooperativity correlation plots highlight several flexibly correlated regions, giving insight into the catalytic mechanism that facilitates access to the active site disulfide bond. Residual native cooperativity correlations are present in the core substructure, suggesting that Trx can function when it is partly unfolded. This natively disordered kinetic trap, interpreted as a molten globule, has a wide temperature range of metastability, and it is identified as the "slow intermediate state" observed in kinetic experiments. These computational results are found to be in overall agreement with a large array of experimental data.

The crystal structure of the cis-proline to glycine variant (P114G) of ribonuclease A

Replacement of a cis-proline by glycine at position 114 in ribonuclease A leads to a large decrease in thermal stability and simplifies the refolding kinetics. A crystallographic approach was used to determine whether the decrease in thermal stability results from the presence of a cis glycine peptide bond, or from a localized structural rearrangement caused by the isomerization of the mutated cis 114 peptide bond. The structure was solved at 2.0 A resolution and refined to an R-factor of 19.5% and an R(free) of 21.9%. The overall conformation of the protein was similar to that of wild-type ribonuclease A; however, there was a large localized rearrangement of the mutated loop (residues 110-117-a 9.3 A shift of the Calpha atom of residue 114). The peptide bond before Gly114 is in the trans configuration. Interestingly, a large anomalous difference density was found near residue 114, and was attributed to a bound cesium ion present in the crystallization experiment. The trans isomeric configuration of the peptide bond in the folded state of this mutant is consistent with the refolding kinetics previously reported, and the associated protein conformational change provides an explanation for the decreased thermal stability.

High-affinity fragment complementation of a fibronectin type III domain and its application to stability enhancement

The tenth fibronectin type III (FN3) domain of human fibronectin (FNfn10), a prototype of the ubiquitous FN3 domain, is a small, monomeric beta-sandwich protein. In this study, we have bisected FNfn10 in each loop to generate a total of six fragment pairs. We found that fragment pairs bisected at multiple loops of FNfn10 show complementation in vivo as tested with a yeast two-hybrid system. The dissociation constant of these fragment pairs determined in vitro were as low as 3 nM, resulting in one of the tightest fragment complementation systems reported so far. Furthermore, we show that the affinity of fragment complementation is correlated with the stability of the uncut parent protein. Exploring this correlation, we screened a yeast two-hybrid library of one fragment and identified mutations that suppress the effect of a destabilizing mutation in the other fragment. One of the identified mutations significantly increased the stability of the uncut wild-type protein, proving that fragment complementation can be used as a novel strategy for the selection of proteins with enhanced stability.

pH dependence of amide chemical shifts in natively disordered polypeptides detects medium-range interactions with ionizable residues

A growing number of natively disordered proteins undergo a folding/binding process that is essential for their biological function. An interesting question is whether these proteins have incompletely solvated regions that drive the folding/binding process. Although the presence of predominantly hydrophobic buried regions can be easily ascertained by high-sensitivity differential scanning calorimetry analysis, the identification of those residues implicated in the burial requires NMR analysis. We have selected a partially solvated natively disordered fragment of Escherichia coli, thioredoxin, C37 (38-108), for full NMR spectral assignment. The secondary chemical shifts, temperature coefficients, and relaxation rates (R(1) and R(2)) of this fragment indicate the presence of a flexible backbone without a stable hydrogen bond network near neutral pH. (1)H-(15)N heteronuclear single quantum coherence analysis of the pH dependence of amide chemical shifts in fragment C37 within pH 2.0 and 7.0 suggests the presence of interactions between nonionizable residues and the carboxylate groups of four Asp and four Glu residues. The pH midpoints (pH(m)) of the amides in the ionizable residues (Asp or Glu) and, consequently, the shifts in the pH(m) (DeltapH(m)) of these residues with respect to model tetrapeptides, are sequence-dependent; and the nonionizable residues that show pH dependence cluster around the ionizable ones. The same pH dependence has been observed in two fragments: M37 (38-73) and C73 (74-108), ruling out the participation of long-range interactions. Our studies indicate the presence of a 15-residue pH-dependent segment with the highest density of ionizable sites in the disordered ensembles of fragments C37 and M37. The observed correlations between ionizable and nonionizable residues in this segment suggest the organization of the backbone and side chains through local and medium-range interactions up to nine residues apart, in contrast to only a few interactions in fragment C73. These results agree qualitatively with the predominantly hydrophobic buried surface detected only in fragments C37 and M37 by highly sensitive differential scanning calorimetry analysis. This work offers a sensitive and rapid new tool to obtain clues about local and nonlocal interactions between ionizable and nonionizable residues in the growing family of natively disordered small proteins with full NMR assignments.

Reduced spectral density mapping of a partially folded fragment of E. coli thioredoxin

The backbone dynamics of a partially folded, N-terminal fragment of E. coli thioredoxin were investigated using nuclear magnetic resonance spectroscopy (NMR). Relaxation data were collected at three temperatures and analyzed using reduced spectral density mapping. As temperature was increased, the values for the viscosity normalized J(0) and for J(omegaH) increased, while J(omegaN) decreased. The global trend observed for the viscosity normalized J(0) was consistent with an increase in the hydrodynamic volume of the fragment and suggested the presence of correlated rotational motion in the absence of long range interactions. In addition, the residue specific variation observed for the viscosity normalized J(0) suggested contributions to J(omega) from a range of correlation times that are close to the global correlation time.

Regio- and stereoselective synthesis of (E)-alkene trans-Xaa-Pro dipeptide mimetics utilizing organocopper-mediated anti-S(N)2' reactions

Proline dipeptides (Xaa-Pro) exist as an equilibrium mixture of cis- and trans-rotamers, which depends on the energy barriers for imide isomerization. This conformation mixture contributes to both structure and function of proline-containing peptides and proteins. Structural motifs resembling these cis- or trans-conformers have served as useful tools for elucidating contributions of proline residues in the physicochemical and biological profiles of structures which contain them. Among such motifs are alkene dipeptide isosteres which mimic cis- or trans-imide using (Z)- or (E)-alkene, respectively. In this report, the first regio- and stereoselective syntheses of (E)-alkene dipeptide isosteres (20, 31, and 35) corresponding to trans-proline dipeptides are described. Key to the synthesis of these mimetics is the anti-S(N)2' reaction of vinyl aziridines such as 15 or vinyl oxazolidinones such as 28 and 32 with organocopper reagents "RCu" (R = CH(2)SiMe(2)(Oi-Pr)). Reaction of cis-vinylaziridine 15 derived from L-serine with organocopper reagent gave a precursor of the trans-L-Ser-D-Pro type alkene isosteres 20, accompanied by an S(N)2 side product. One limitation with the use of such aziridine-mediated methodology is formation of the corresponding trans-aziridine 22, which leads to L-L type isosteres, that is unstable and obtainable only in low yield. On the other hand, both isomers of oxazolidinone derivatives can be easily obtained from N-Boc-protected amino alcohols. The reaction of trans- 28 or cis-oxazolidinone derivative 32 with organocopper reagents proceeds quantitatively with high regio- and diastereoselectivities in anti-S(N)2' fashion. Subsequent oxidative treatment of the newly introduced isopropoxydimethylsilylmethyl group yields trans-L-Ser-L-Pro 31 or trans-L-Ser-D-Pro type isosteres 35, respectively. Of note, synthesized isostere 31 can also be converted to trans-phosphoSer-Pro 42 and trans-Cys-Pro mimetics 44. The present synthetic methodology affords trans-Xaa-Pro alkene-type dipeptide isosteres in high yield with relatively simple manipulation.

Local structural changes caused by peptidyl-prolyl cis/trans isomerization in the native state of proteins

Peptidyl-prolyl cis/trans isomerization, observed in the native state of an increasing number of proteins, is of considerable biological significance. The first evidence for an asymmetric transmission along the polypeptide chain of the structural effects of prolyl isomerization is now derived from the statistics of the C(alpha)/C(alpha)-atom distance distributions in the crystal structures of 848 non-homologous proteins. More detailed information on how isomerization affects segments adjacent to proline is obtained from crystal structures of proteins, that are more than 95% homologous, and that exhibit two different states of isomerization at a particular prolyl bond. The resulting 64 cases, which represent 3.8% of the database used, form pairs of coordinates which were analyzed for the existence of isomer-specific intramolecular nonbonded C(alpha)/C(alpha)-atom distances around the critical proline, and for the positional preferences for particular amino acids in the isomeric sequence segment. The probability that a native protein exhibits both prolyl isomers in the crystalline state increases in particular with a Pro at the third position N-terminal to the isomeric bond (-3 position), and with Ser, Gly and Asp at the position preceding the isomeric bond (-1 position). Structural alignment of matched pairs of isomeric proteins generates three classes with respect to position-specific distribution of C(alpha)-atom displacements around an isomeric proline imide bond. In the majority of cases the distribution of these intermolecular isomer-specific C(alpha)-atom distances shows a symmetric behavior for the N-terminal and C-terminal segment flanking the proline residue, and the magnitude did not exceed 1.3+/-0.6 A including the C(alpha) atoms in proximity to the prolyl bond. However, in the remaining 12 protein pairs the structural changes are unidirectional relative to the isomerizing bond whereby the magnitude of the isomer-specific effect exceeds 3.0+/-2.0 A even at positions remote to proline. Interestingly, the magnitude of the intramolecular isomer-specific C(alpha) atom displacements reveals a lever-arm amplification of the isomerization-mediated structural changes in a protein backbone. The observed backbone effects provide a structural basis for isomer-specific reactions of proline-containing polypeptides, and thus may play a role in biological recognition and regulation.

Reconstitution of a native-like SH2 domain from disordered peptide fragments examined by multidimensional heteronuclear NMR

The N-terminal SH2 domain from the p85alpha subunit of phosphatidylinositol 3' kinase is cleaved specifically into 9- and 5-kD fragments by limited proteolytic digestion with trypsin. The noncovalent SH2 domain complex and its constituent tryptic peptides have been investigated using high-resolution heteronuclear magnetic resonance (NMR). These studies have established the viability of the SH2 domain as a fragment complementation system. The individual peptide fragments are predominantly unstructured in solution. In contrast, the noncovalent 9-kD + 5-kD complex shows a native-like (1)H-(15)N HSQC spectrum, demonstrating that the two fragments fold into a native-like structure on binding. Chemical shift analysis of the noncovalent complex compared to the native SH2 domain reveals that the highest degree of perturbation in the structure occurs at the cleavage site within a flexible loop and along the hydrophobic interface between the two peptide fragments. Mapping of these chemical shift changes on the structure of the domain reveals changes consistent with the reduction in affinity for the target peptide ligand observed in the noncovalent complex relative to the intact protein. The 5-kD fragment of the homologous Src protein is incapable of structurally complementing the p85 9-kD fragment, either in complex formation or in the context of the full-length protein. These high-resolution structural studies of the SH2 domain fragment complementation features establish the suitability of the system for further protein-folding and design studies.

Heat capacity analysis of oxidized Escherichia coli thioredoxin fragments (1--73, 74--108) and their noncovalent complex. Evidence for the burial of apolar surface in protein unfolded states

We have calculated the absolute heat capacities of fragments 1--73 (N fragment) and 74--108 (C fragment) from thioredoxin, their complex and the uncleaved protein, from the concentration dependence of the apparent heat capacities of the solutions determined by differential scanning calorimetry. We find that, while the absolute heat capacities of uncleaved, unfolded thioredoxin and the C fragment are in good agreement with the theoretical values expected for fully solvated chains (calculated as the sum of the contributions of the constituent amino acids), the absolute heat capacities of the N fragment and the unfolded complex are about 2 kJ x K(-1) x mol(-1) lower than the fully solvated-chain values. We attribute this discrepancy to burial of the apolar surface in the N fragment (as burial of the polar area is expected to lead to an increase in heat capacity). Illustrative calculations suggest that burial of about 1000--1600 A(2) of apolar surface takes place in the N fragment (probably accompanied by the burial of a smaller amount of polar surface). In general, this work is supportive of heat capacity measurements on protein fragments being useful as probes of surface burial in studies to characterize protein unfolded states and the high regions of protein folding landscapes.