Refinement of Peptide Conformational Ensembles by 2D IR Spectroscopy: Application to Ala‒Ala‒Ala

Characterizing ensembles of intrinsically disordered proteins is experimentally challenging because of the ill-conditioned nature of ensemble determination with limited data and the intrinsic fast dynamics of the conformational ensemble. Amide I two-dimensional infrared (2D IR) spectroscopy has picosecond time resolution to freeze structural ensembles as needed for probing disordered-protein ensembles and conformational dynamics. Also, developments in amide I computational spectroscopy now allow a quantitative and direct prediction of amide I spectra based on conformational distributions drawn from molecular dynamics simulations, providing a route to ensemble refinement against experimental spectra. We performed a Bayesian ensemble refinement method on Ala–Ala–Ala against isotope-edited Fourier-transform infrared spectroscopy and 2D IR spectroscopy and tested potential factors affecting the quality of ensemble refinements. We found that isotope-edited 2D IR spectroscopy provides a stringent constraint on Ala–Ala–Ala conformations and returns consistent conformational ensembles with the dominant ppII conformer across varying prior distributions from many molecular dynamics force fields and water models. The dominant factor influencing ensemble refinements is the systematic frequency uncertainty from spectroscopic maps. However, the uncertainty of conformer populations can be significantly reduced by incorporating 2D IR spectra in addition to traditional Fourier-transform infrared spectra. Bayesian ensemble refinement against isotope-edited 2D IR spectroscopy thus provides a route to probe equilibrium-complex protein ensembles and potentially nonequilibrium conformational dynamics.     

Solution structure and internal dynamics of NSCP, a compact calcium-binding protein

The solution structure of Nereis diversicolor sarcoplasmic calcium-binding protein (NSCP) in the calcium-bound form was determined by NMR spectroscopy, distance geometry and simulated annealing. Based on 1859 NOE restraints and 262 angular restraints, 17 structures were generated with a rmsd of 0.87 A from the mean structure. The solution structure, which is highly similar to the structure obtained by X-ray crystallography, includes two open EF-hand domains, which are in close contact through their hydrophobic surfaces. The internal dynamics of the protein backbone were determined by studying amide hydrogen/deuterium exchange rates and 15N nuclear relaxation. The two methods revealed a highly compact and rigid structure, with greatly restricted mobility at the two termini. For most of the amide protons, the free energy of exchange-compatible structural opening is similar to the free energy of structural stability, suggesting that isotope exchange of these protons takes place through global unfolding of the protein. Enhanced conformational flexibility was noted in the unoccupied Ca2+-binding site II, as well as the neighbouring helices. Analysis of the experimental nuclear relaxation and the molecular dynamics simulations give very similar profiles for the backbone generalized order parameter (S2), a parameter related to the amplitude of fast (picosecond to nanosecond) movements of N(H)-H vectors. We also noted a significant correlation between this parameter, the exchange rate, and the crystallographic B factor along the sequence.     

Structural interpretation of hydrogen exchange protection factors in proteins: characterization of the native state fluctuations of CI2

Protection factors obtained from equilibrium hydrogen exchange experiments are an important source of structural information on both native and nonnative states of proteins. We present a method for determining ensembles of protein structures by using hydrogen exchange data as restraints in molecular dynamics simulations in conjunction with an empirical force-field. The method is applied to determine the ensemble of structures representing the native state of chymotrypsin inhibitor 2 (CI2), including the rare, large fluctuations responsible for hydrogen exchange.     

A Correspondence Between Solution-State Dynamics of an Individual Protein and the Sequence and Conformational Diversity of its Family

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using “Backrub” motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.     

Measuring 1HN temperature coefficients in invisible protein states by relaxation dispersion NMR spectroscopy

A method based on the Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for measuring the temperature coefficients of amide proton chemical shifts of low populated ‘invisible’ protein states that exchange with a ‘visible’ ground state on the millisecond time-scale. The utility of the approach is demonstrated with an application to an I58D mutant of the Pfl6 Cro protein that undergoes exchange between the native, folded state and a cold denatured, unfolded conformational ensemble that is populated at a level of 6% at 2.5°C. A wide distribution of amide temperature coefficients is measured for the unfolded state. The distribution is centered about –5.6 ppb/K, consistent with an absence of intra-molecular hydrogen bonds, on average. However, the large range of values (standard deviation of 2.1 ppb/K) strongly supports the notion that the unfolded state of the protein is not a true random coil polypeptide chain.     

Paramagnetic relaxation enhancement‐assisted structural characterization of a partially disordered conformation of ubiquitin

Nuclear magnetic resonance (NMR) is a powerful tool to study three‐dimensional structures as well as protein conformational fluctuations in solution, but it is compromised by increases in peak widths and missing signals. We previously reported that ubiquitin has two folded conformations, N₁ and N₂ and plus another folded conformation, I, in which some amide group signals of residues 33–41 almost disappeared above 3 kbar at pH 4.5 and 273 K. Thus, well‐converged structural models could not be obtained for this region owing to the absence of distance restraints. Here, we reexamine the problem using the ubiquitin Q41N variant as a model for this locally disordered conformation, I. We demonstrate that the variant shows pressure‐induced loss of backbone amide group signals at residues 28, 33, 36, and 39–41 like the wild‐type, with a similar but smaller effect on CαH and CβH signals. In order to characterize this I structure, we measured paramagnetic relaxation enhancement (PRE) under high pressure to obtain distance restraints, and calculated the structure assisted by Bayesian inference. We conclude that the more disordered I conformation observed at pH 4.0, 278 K, and 2.5 kbar largely retained the N₂ conformation, although the amide groups at residues 33–41 have more heterogeneous conformations and more contact with water, which differ from the N₁ and N₂ states. The PRE‐assisted strategy has the potential to improve structural characterization of proteins that lack NMR signals, especially for relatively more open and hydrated protein conformations.     

The conformational characteristics in solution of the synthetic cyclic hexapeptide 5L-ala-D-ala

An attempt to elucidate the solution conformation(s) of the synthetic cyclic hexapeptide 5L -ala·D-ala is described. Nuclear magnetic resonance (nmr) spectra are recorded for the purpose of measuring the vicinal coupling constant between the amide and α-protons in each residue and to observe the deuterium exchange rate and temperature dependence of the chemical shift of each amide proton. Low-energy cyclic conformations, whose individual residues are in conformations consistent with the observed amide to α-proton coupling constant, are searched for in an approximate theoretical treatment. The two lowest energy, all trans peptide bond conformations generated are distinguishable by the presence or absence of a single intramolecular hydrogen bond. The observed temperature independence of the chemical shift of one of the amide protons is consistent with the presence of a single intramolecular hydrogen bond, while the observation of similar deuterium exchange rates for each of the amide protons indicates their comparable availability to solvent. Consequently, it is concluded that 5L -ala·D-ala is in rapid equilibrium between conformations with and without a single internal hydrogen bond and possesses considerable conformational flexibility in solution.     

Improved structural characterizations of the drkN SH3 domain unfolded state suggest a compact ensemble with native-like and non-native structure

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

Role of native-state structure in rubredoxin native-state hydrogen exchange

Amide exchange rates were measured for Pyrococcus furiosus (Pf) rubredoxin substituted with either Zn(II), Ga(III), or Ge(IV). Base-catalyzed exchange rate constants increase up to 3000-fold per unit charge for the highly protected amides surrounding the active site metal, yielding apparent residue-specific conformational energy decreases of more than 8 kcal/mol in a comparison of the Zn(II)- and Ge(IV)-substituted proteins. However, the exchange kinetics for many of the other amides of the protein are insensitive to these metal substitutions. These differential rates are inversely correlated with the distance between the amide nitrogen and the metal in the X-ray structure, out to a distance of at least 12 A, consistent with an electrostatic potential-dependent shifting of the amide nitrogen pK. This strongly correlated distance dependence is consistent with a nativelike structure for the exchange-competent conformations. The electric field potential within the interior of the rubredoxin structure gives rise to a change of as much as a million-fold in the rate for the exchange-competent state of the individual amide hydrogens. Nevertheless, the strength of these electrostatic interactions in Pf rubredoxin appears to be comparable to those previously reported within other proteins. As a result, contrary to the conventional analysis of hydrogen exchange data, for exchange processes that occur via nonglobal transitions, the residual conformational structure will often modulate the observed rates. Although this necessarily complicates the estimation of the conformational equilibria of these exchange-competent states, this dependence on residual structure can provide insight into the conformation of these transient states.     

Temperature-dependence of protein hydrogen bond properties as studied by high-resolution NMR

The temperature-dependence of a large number of NMR parameters describing hydrogen bond properties in the protein ubiquitin was followed over a range from 5 to 65 degrees C. The parameters comprise hydrogen bond (H-bond) scalar couplings, h3JNC', chemical shifts, amide proton exchange rates, 15N relaxation parameters as well as covalent 1JNC' and 1JNH couplings. A global weakening of the h3JNC' coupling with increasing temperature is accompanied by a global upfield shift of the amide protons and a decrease of the sequential 1JNC' couplings. If interpreted as a linear increase of the N...O distance, the change in h3JNC' corresponds to an average linear thermal expansion coefficient for the NH-->O hydrogen bonds of 1.7 x 10(-4)/K, which is in good agreement with overall volume expansion coefficients observed for proteins. A residue-specific analysis reveals that not all hydrogen bonds are affected to the same extent by the thermal expansion. The end of beta-sheet beta1/beta5 at hydrogen bond E64-->Q2 appears as the most thermolabile, whereas the adjacent hydrogen bond I3-->L15 connecting beta-strands beta1 and beta2 is even stabilized slightly at higher temperatures. Additional evidence for the stabilization of the beta1/beta2 beta-hairpin at higher temperatures is found in reduced hydrogen exchange rates for strand end residue V17. This reduction corresponds to a stabilizing change in free energy of 9.7 kJ/mol for the beta1/beta2 hairpin. The result can be linked to the finding that the beta1/beta2 hairpin behaves as an autonomously folding unit in the A-state of ubiquitin under changed solvent conditions. For several amide groups the temperature-dependencies of the amide exchange rates and H-bond scalar couplings are uncorrelated. Therefore, amide exchange rates are not a sole function of the hydrogen bond "strength" as given by the electronic overlap of donors and acceptors, but are clearly dependent on other blocking mechanisms.     

Representing an ensemble of NMR-derived protein structures by a single structure.

The usefulness of representing an ensemble of NMR-derived protein structures by a single structure has been investigated. Two stereochemical properties have been used to assess how a single structure relates to the ensemble from which it was derived, namely the distribution of phi psi torsion angles and the distribution of chi 1 torsion angles. The results show that the minimized average structure derived from the ensemble (a total of 11 ensembles from the Brookhaven Protein Data Bank were analyzed) does not always correspond well with this ensemble, particularly for those ensembles generated with a smaller number of experimentally derived restraints per residue. An alternative method that selects the member of the ensemble which is closest to the "average" of the ensemble has been investigated (a total of 23 ensembles from the Brookhaven Protein Data Bank were analyzed). Although this method selected a structure that on the whole corresponded more closely to the ensemble than did the minimized average structure, this is still not a totally reliable means of selecting a single structure to represent the ensemble. This suggests that it is advisable to study the ensemble as a whole. A study has also been made of the practice of selecting the "best" rather than the most representative member of the ensemble. This too suggests that the ensemble should be studied as a whole. A study of the conformational space occupied by the ensemble also suggests the need to consider the ensemble as a whole, particularly for those ensembles generated with a smaller number of experimentally derived restraints per residue.     

Iterative assembly of helical proteins by optimal hydrophobic packing

We present a method for the computer-based iterative assembly of native-like tertiary structures of helical proteins from alpha-helical fragments. For any pair of helices, our method, called MATCHSTIX, first generates an ensemble of possible relative orientations of the helices with various ways to form hydrophobic contacts between them. Those conformations having steric clashes, or a large radius of gyration of hydrophobic residues, or with helices too far separated to be connected by the intervening linking region, are discarded. Then, we attempt to connect the two helical fragments by using a robotics-based loop-closure algorithm. When loop closure is feasible, the algorithm generates an ensemble of viable interconnecting loops. After energy minimization and clustering, we use a representative set of conformations for further assembly with the remaining helices, adding one helix at a time. To efficiently sample the conformational space, the order of assembly generally proceeds from the pair of helices connected by the shortest loop, followed by joining one of its adjacent helices, always proceeding with the shorter connecting loop. We tested MATCHSTIX on 28 helical proteins each containing up to 5 helices and found it to heavily sample native-like conformations. The average rmsd of the best conformations for the 17 helix-bundle proteins that have 2 or 3 helices is less than 2 A; errors increase somewhat for proteins containing more helices. Native-like states are even more densely sampled when disulfide bonds are known and imposed as restraints. We conclude that, at least for helical proteins, if the secondary structures are known, this rapid rigid-body maximization of hydrophobic interactions can lead to small ensembles of highly native-like structures. It may be useful for protein structure prediction.     

Conformational studies of the N-terminal lipid-associating domain of human apolipoprotein C-I by CD and 1H NMR spectroscopy.

A peptide comprising the N-terminal 38 residues of human apolipoprotein C-I (apoC-I(1-38)) was synthesized using solid-phase methods and its solution conformation studied by CD and 1H NMR spectroscopy. The CD data indicate that apoC-I(1-38) has a similar helical content (55%) in the presence of saturating amounts of SDS or egg yolk lysophosphatidylcholine. A structural ensemble of SDS-bound apoC-I(1-38) was calculated from 464 NOE-based distance restraints using distance geometry methods. ApoC-I(1-38) adopts a helical structure between residues V4 and K30 and an extended C-terminus from Q31 when associated with SDS. The region K12-G15 undergoes slow conformational exchange as indicated by above-average amide resonance linewidths, large temperature coefficients, and fast exchange (< 2 h) of backbone amide protons with deuterium. The mobility of K12-G15 is reflected in the poorly defined dihedral angles of K12 and E13 in the calculated ensemble of structures. The average structure of apoC-I(1-38) is curved toward its hydrophobic face with bends of 125 degrees, centered at K12/E13, and 150 degrees, centered at K21. This curvature appears to be driven by the interaction of two hydrophobic clusters, one formed by residues L8, L11, F14, and L18, and the other by L25, I26, and I29, with the amphiphile SDS. Based on our present structural definition of apoC-I(1-38) and the previously obtained structure of the fragment apoC-I(35-53), we propose the secondary structure of intact apolipoprotein C-I.     

Determination of the transition state ensemble for the folding of ubiquitin from a combination of Phi and Psi analyses

Protein engineering techniques have emerged as powerful tools for characterizing transition states (TSs) for protein folding. Recently, the Ψ analysis, in which double-histidine mutations create the possibility of reversible crosslinking in the native state, has been proposed as an additional approach to the well-established Φ analysis. We present here a combination of these two procedures for defining the structure of the TS of ubiquitin, a small α/β protein that has been used extensively as a model system for both experimental and computational studies of the protein-folding process. We performed a series of molecular dynamics simulations in which Φ and Ψ values were used as ensemble-averaged structural restraints to determine an ensemble of structures representing the TS of ubiquitin. Although the available Ψ values for ubiquitin did not, by themselves, generate well-defined TS ensembles, the inclusion of the restricted set of zero or unity values, but not fractional ones, provided useful complementary information to the Φ analysis. Our results show that the TS of ubiquitin is formed by a relatively narrow ensemble of structures exhibiting an overall native-like topology in which the N-terminal and C-terminal regions are in close proximity.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

The Effect of a ΔK280 Mutation on the Unfolded State of a Microtubule-Binding Repeat in Tau

Tau is a natively unfolded protein that forms intracellular aggregates in the brains of patients with Alzheimer''s disease. To decipher the mechanism underlying the formation of tau aggregates, we developed a novel approach for constructing models of natively unfolded proteins. The method, energy-minima mapping and weighting (EMW), samples local energy minima of subsequences within a natively unfolded protein and then constructs ensembles from these energetically favorable conformations that are consistent with a given set of experimental data. A unique feature of the method is that it does not strive to generate a single ensemble that represents the unfolded state. Instead we construct a number of candidate ensembles, each of which agrees with a given set of experimental constraints, and focus our analysis on local structural features that are present in all of the independently generated ensembles. Using EMW we generated ensembles that are consistent with chemical shift measurements obtained on tau constructs. Thirty models were constructed for the second microtubule binding repeat (MTBR2) in wild-type (WT) tau and a ΔK280 mutant, which is found in some forms of frontotemporal dementia. By focusing on structural features that are preserved across all ensembles, we find that the aggregation-initiating sequence, PHF6*, prefers an extended conformation in both the WT and ΔK280 sequences. In addition, we find that residue K280 can adopt a loop/turn conformation in WT MTBR2 and that deletion of this residue, which can adopt nonextended states, leads to an increase in locally extended conformations near the C-terminus of PHF6*. As an increased preference for extended states near the C-terminus of PHF6* may facilitate the propagation of β-structure downstream from PHF6*, these results explain how a deletion at position 280 can promote the formation of tau aggregates.     

Cyclic AMP analog blocks kinase activation by stabilizing inactive conformation: conformational selection highlights a new concept in allosteric inhibitor design

The regulatory (R) subunit of protein kinase A serves to modulate the activity of protein kinase A in a cAMP-dependent manner and exists in two distinct and structurally dissimilar, end point cAMP-bound "B" and C-subunit-bound "H"-conformations. Here we report mechanistic details of cAMP action as yet unknown through a unique approach combining x-ray crystallography with structural proteomics approaches, amide hydrogen/deuterium exchange and ion mobility mass spectrometry, applied to the study of a stereospecific cAMP phosphorothioate analog and antagonist((Rp)-cAMPS). X-ray crystallography shows cAMP-bound R-subunit in the B form but surprisingly the antagonist Rp-cAMPS-bound R-subunit crystallized in the H conformation, which was previously assumed to be induced only by C-subunit-binding. Apo R-subunit crystallized in the B form as well but amide exchange mass spectrometry showed large differences between apo, agonist and antagonist-bound states of the R-subunit. Further ion mobility reveals the apo R-subunit as an ensemble of multiple conformations with collisional cross-sectional areas spanning both the agonist and antagonist-bound states. Thus contrary to earlier studies that explained the basis for cAMP action through "induced fit" alone, we report evidence for conformational selection, where the ligand-free apo form of the R-subunit exists as an ensemble of both B and H conformations. Although cAMP preferentially binds the B conformation, Rp-cAMPS interestingly binds the H conformation. This reveals the unique importance of the equatorial oxygen of the cyclic phosphate in mediating conformational transitions from H to B forms highlighting a novel approach for rational structure-based drug design. Ideal inhibitors such as Rp-cAMPS are those that preferentially "select" inactive conformations of target proteins by satisfying all "binding" constraints alone without inducing conformational changes necessary for activation.     

The structure and dynamics of rat apo-cellular retinol-binding protein II in solution: comparison with the X-ray structure

The structure and dynamics of rat apo-cellular retinol binding protein II (apo-CRBP II) in solution has been determined by multidimensional NMR analysis of uniformly enriched recombinant rat 13C, 15N-apo-CRBP II and 15N-apo-CRBP II. The final ensemble of 24 NMR structures has been calculated from 3274 conformational restraints or 24.4 restraints/residue. The average root-mean-square deviation of the backbone atoms for the final 24 structures relative to their mean structure is 1.06 A. Although the average solution structure is very similar to the crystal structure, it differs at the putative entrance to the binding cavity, which is formed by the helix-turn-helix motif, the betaC-betaD turn and the betaE-betaF turn. The mean coordinates of the main-chain atoms of amino acid residues 28-38 are displaced in the solution structure relative to the crystal structure. The side-chain of F58, located on the betaC-betaD turn, is reoriented such that it interacts with L37 and no longer blocks entry into the ligand-binding pocket. Residues 28-35, which form the second helix of the helix-turn-helix motif in the crystal structure, do not exhibit a helical conformation in the solution structure. The solution structure of apo-CRBP II exhibits discrete regions of backbone disorder which are most pronounced at residues 28-32, 37-38 and 73-76 in the betaE-betaF turn as evaluated by the consensus chemical shift index, the root-mean-square deviation, amide 1H exchange rates and 15N relaxation studies. These studies indicate that fluctuations in protein conformation occur on the microseconds to ms time-scale in these regions of the protein. Some of these exchange processes can be directly observed in the three-dimensional 15N-resolved NOESY spectrum. These results suggest that in solution, apo-CRBP II undergoes conformational changes on the microseconds to ms time-scale which result in increased access to the binding cavity.     

The effect of a DeltaK280 mutation on the unfolded state of a microtubule-binding repeat in Tau

Tau is a natively unfolded protein that forms intracellular aggregates in the brains of patients with Alzheimer's disease. To decipher the mechanism underlying the formation of tau aggregates, we developed a novel approach for constructing models of natively unfolded proteins. The method, energy-minima mapping and weighting (EMW), samples local energy minima of subsequences within a natively unfolded protein and then constructs ensembles from these energetically favorable conformations that are consistent with a given set of experimental data. A unique feature of the method is that it does not strive to generate a single ensemble that represents the unfolded state. Instead we construct a number of candidate ensembles, each of which agrees with a given set of experimental constraints, and focus our analysis on local structural features that are present in all of the independently generated ensembles. Using EMW we generated ensembles that are consistent with chemical shift measurements obtained on tau constructs. Thirty models were constructed for the second microtubule binding repeat (MTBR2) in wild-type (WT) tau and a DeltaK280 mutant, which is found in some forms of frontotemporal dementia. By focusing on structural features that are preserved across all ensembles, we find that the aggregation-initiating sequence, PHF6*, prefers an extended conformation in both the WT and DeltaK280 sequences. In addition, we find that residue K280 can adopt a loop/turn conformation in WT MTBR2 and that deletion of this residue, which can adopt nonextended states, leads to an increase in locally extended conformations near the C-terminus of PHF6*. As an increased preference for extended states near the C-terminus of PHF6* may facilitate the propagation of beta-structure downstream from PHF6*, these results explain how a deletion at position 280 can promote the formation of tau aggregates.