Conformation of valine side chains in ribonuclease T1 determined by NMR studies of homonuclear and heteronuclear 3J coupling constants

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.

Using 1H and 13C NMR chemical shifts to determine cyclic peptide conformations: a combined molecular dynamics and quantum mechanics approach

Conformational analysis from NMR and density-functional prediction of low-energy ensembles (CANDLE), a new approach for determining solution structures.

A J-modulated protonless NMR experiment characterizes the conformational ensemble of the intrinsically disordered protein WIP

Intrinsically disordered proteins (IDPs) are multi-conformational polypeptides that lack a single stable three-dimensional structure. It has become increasingly clear that the versatile IDPs play key roles in a multitude of biological processes, and, given their flexible nature, NMR is a leading method to investigate IDP behavior on the molecular level. Here we present an IDP-tailored J-modulated experiment designed to monitor changes in the conformational ensemble characteristic of IDPs by accurately measuring backbone one- and two-bond J(¹⁵N,¹³Cα) couplings. This concept was realized using a unidirectional (H)NCO ¹³C-detected experiment suitable for poor spectral dispersion and optimized for maximum coverage of amino acid types. To demonstrate the utility of this approach we applied it to the disordered actin-binding N-terminal domain of WASp interacting protein (WIP), a ubiquitous key modulator of cytoskeletal changes in a range of biological systems. One- and two-bond J(¹⁵N,¹³Cα) couplings were acquired for WIP residues 2-65 at various temperatures, and in denaturing and crowding environments. Under native conditions fitted J-couplings identified in the WIP conformational ensemble a propensity for extended conformation at residues 16-23 and 45-60, and a helical tendency at residues 28-42. These findings are consistent with a previous study of the based upon chemical shift and RDC data and confirm that the WIP^2-65 conformational ensemble is biased towards the structure assumed by this fragment in its actin-bound form. The effects of environmental changes upon this ensemble were readily apparent in the J-coupling data, which reflected a significant decrease in structural propensity at higher temperatures, in the presence of 8 M urea, and under the influence of a bacterial cell lysate. The latter suggests that crowding can cause protein unfolding through protein-protein interactions that stabilize the unfolded state. We conclude that J-couplings are a useful measureable in characterizing structural ensembles in IDPs, and that the proposed experiment provides a practical method for accurately performing such measurements, once again emphasizing the power of NMR in studying IDP behavior.

Side Chain Conformational Distributions of a Small Protein Derived from Model-Free Analysis of a Large Set of Residual Dipolar Couplings

Accurate quantitative measurement of structural dispersion in proteins remains a prime challenge to both X-ray crystallography and NMR spectroscopy. Here we use a model-free approach based on measurement of many residual dipolar couplings (RDCs) in differentially orienting aqueous liquid crystalline solutions to obtain the side chain χ₁ distribution sampled by each residue in solution. Applied to the small well-ordered model protein GB3, our approach reveals that the RDC data are compatible with a single narrow distribution of side chain χ₁ angles for only about 40% of the residues. For more than half of the residues, populations greater than 10% for a second rotamer are observed, and four residues require sampling of three rotameric states to fit the RDC data. In virtually all cases, sampled χ₁ values are found to center closely around ideal g^–, g⁺ and t rotameric angles, even though no rotamer restraint is used when deriving the sampled angles. The root-mean-square difference between experimental ³J_HαHβ couplings and those predicted by the Haasnoot-parametrized, motion-adjusted Karplus equation reduces from 2.05 to 0.75 Hz when using the new rotamer analysis instead of the 1.1-Å X-ray structure as input for the dihedral angles. A comparison between observed and predicted ³J_HαHβ values suggests that the root-mean-square amplitude of χ₁ angle fluctuations within a given rotamer well is ca. 20°. The quantitatively defined side chain rotamer equilibria obtained from our study set new benchmarks for evaluating improved molecular dynamics force fields, and also will enable further development of quantitative relations between side chain chemical shift and structure.

Transforming between discrete and continuous angle distribution models: application to protein χ₁ torsions

Two commonly employed angular-mobility models for describing amino-acid side-chain χ(1) torsion conformation, the staggered-rotamer jump and the normal probability density, are discussed and performance differences in applications to scalar-coupling data interpretation highlighted. Both models differ in their distinct statistical concepts, representing discrete and continuous angle distributions, respectively. Circular statistics, introduced for describing torsion-angle distributions by using a universal circular order parameter central to all models, suggest another distribution of the continuous class, here referred to as the elliptic model. Characteristic of the elliptic model is that order parameter and circular variance form complementary moduli. Transformations between the parameter sets that describe the probability density functions underlying the different models are provided. Numerical aspects of parameter optimization are considered. The issues are typified by using a set of χ(1) related (3) J coupling constants available for FK506-binding protein. The discrete staggered-rotamer model is found generally to produce lower order parameters, implying elevated rotatory variability in the amino-acid side chains, whereas continuous models tend to give higher order parameters that suggest comparatively less variation in angle conformations. The differences perceived regarding angular mobility are attributed to conceptually different features inherent to the models.

A method to explore protein side chain conformational variability using experimental data

Experimentally measured values of molecular properties or observables of biomolecules such as proteins are generally averages over time and space, which do not contain sufficient information to determine the underlying conformational distribution of the molecules in solution. The relationship between experimentally measured NMR (3)J-coupling values and the corresponding dihedral angle values is a particularly complicated case due to its nonlinear, multiple-valued nature. Molecular dynamics (MD) simulations at constant temperature can generate Boltzmann ensembles of molecular structures that are free from a priori assumptions about the nature of the underlying conformational distribution. They suffer, however, from limited sampling with respect to time and conformational space. Moreover, the quality of the obtained structures is dependent on the choice of force field and solvation model. A recently proposed method that uses time-averaging with local-elevation (LE) biasing of the conformational search provides an elegant means of overcoming these three problems. Using a set of side chain (3)J-coupling values for the FK506 binding protein (FKBP), we first investigate the uncertainty in the angle values predicted theoretically. We then propose a simple MD-based technique to detect inconsistencies within an experimental data set and identify degrees of freedom for which conformational averaging takes place or for which force field parameters may be deficient. Finally, we show that LE MD is the best method for producing ensembles of structures that, on average, fit the experimental data.

Accurate measurement of methyl 13C chemical shifts by solid-state NMR for the determination of protein side chain conformation: the influenza a M2 transmembrane peptide as an example

The use of side chain methyl (13)C chemical shifts for the determination of the rotameric conformation of Val and Leu residues in proteins by solid-state NMR spectroscopy is described. Examination of the solution NMR stereospecifically assigned methyl groups shows significant correlation between the difference in the two methyl carbons' chemical shifts and the side chain conformation. It is found that alpha-helical and beta-sheet backbones cause different side chain methyl chemical shift trends. In alpha-helical Leu's, a relatively large absolute methyl (13)C shift difference of 2.89 ppm is found for the most populated mt rotamer (chi(1) = -60 degrees, chi(2) = 180 degrees), while a much smaller value of 0.73 ppm is found for the next populated tp rotamer (chi(1) = 180 degrees, chi(2) = 60 degrees). For alpha-helical Val residues, the dominant t rotamer (chi(1) = 180 degrees) has more downfield Cgamma2 chemical shifts than Cgamma1 by 1.71 ppm, while the next populated m rotamer (chi(1) = -60 degrees) shows the opposite trend of more downfield Cgamma1 chemical shift by 1.23 ppm. These significantly different methyl (13)C chemical shifts exist despite the likelihood of partial rotameric averaging at ambient temperature. We show that these conformation-dependent methyl (13)C chemical shifts can be utilized for side chain structure determination once the methyl (13)C resonances are accurately measured by double-quantum (DQ) filtered 2D correlation experiments, most notably the dipolar DQ to single-quantum (SQ) correlation technique. The advantage of the DQ-SQ correlation experiment over simple 2D SQ-SQ correlation experiments is demonstrated on the transmembrane peptide of the influenza A M2 proton channel. The methyl chemical shifts led to predictions of the side chain rotameric states for several Val and Leu residues in this tetrameric helical bundle. The predicted Val rotamers were further verified by dipolar correlation experiments that directly measure the chi(1) torsion angles. It was found that the chemical-shift-predicted side chain conformations are fully consistent with the direct torsion angle results; moreover, the methyl (13)C chemical shifts are sensitive to approximately 5 degrees changes in the chi(1) torsion angle due to drug binding.

Mechanistic analysis of a suicide inactivator of histone demethylase LSD1

Lysine-specific demethylase 1 (LSD1) is a transcriptional repressor and a flavin-dependent amine oxidase that is responsible for the removal of methyl from lysine 4 of histone H3. In this study, we characterize the mechanism and scope of LSD1 inhibition by a propargylamine-derivatized histone H3 substrate (1). Unlike aziridinyl and cyclopropylamine-derivatized histone H3 peptide substrate analogues, compound 1 appears to covalently modify and irreversibly inactivate LSD1 with high potency. Accompanying this inactivation is a spectroscopic change, which shifts the absorbance maximum to 392 nm. Spectral changes associated with the 1-LSD1 complex and reactivity to decreased pH and sodium borohydride treatment were suggestive of a structure involving a flavin-linked inhibitor conjugate between N5 of the flavin and the terminal carbon of the inhibitor. Using a 13C-labeled inhibitor, NMR analysis of the 1-flavin conjugate was consistent with this structural assignment. Kinetic analysis of the spectroscopic shift induced by 1 showed that the flavin adduct formed in a reaction with kinetic constants similar to those of the LSD1 inactivation process. Taken together, these data support a mechanism of LSD1 inactivation by 1 involving amine oxidation followed by Michael addition to the propargylic imine. We further examined the potential for a biotinylated analogue of 1 (1-Btn) to be used as a tool in affinity pulldown experiments. Using 1-Btn, it was feasible to selectively pull down spiked and endogenous LSD1 from HeLa cell nuclear extracts, setting the stage for activity-based demethylase proteomics.

A versatile component-coupling model to account for substituent effects: application to polypeptide phi and chi(1) torsion related (3)J data

A model is proposed for collating fundamental and incremental component couplings to account for substituent effects on (3)J arising from, for example, amino-acid type variation. The unique topology patterns encountered in each of the common amino acids were modeled by assigning substituents on a (3)J coupling path to four simple categories comprising only relative positions: central (inner) vs. terminal (outer) and first-sphere vs. second-sphere. Associated increment values then reflect the influences on each (3)J coupling accessible for torsion-angle determination. Facility of use of this model, in comparison with previous ones, owes to its strict limitation to no more than three Karplus coefficients for each specific torsion-angle dependency derived. The model was integrated in the concept of self-consistent (3)J analysis and applied to polypeptide fragments X-N-C(alpha)-Y and X-C(alpha)-C(beta)-Y related to torsions phi and chi(1), respectively, yielding quantitative effects of both first- and second-sphere substituents. Regarding the polypeptide backbone, the model predicts first-sphere substituent effects on phi-related (3)J couplings to be within experimental uncertainty because main-chain topologies are identical in most amino-acid types, except for marginal effects in glycine and proline. However, effects in excess of standard errors in (3)J(phi) measurements are anticipated from second-sphere substituent variation. Regarding amino-acid side chains, first-sphere substituent effects on chi(1)-related (3)J couplings were previously found pivotal to accurate torsion-angle interpretation. Taking additional second-sphere effects on (3)J(chi(1)) into account is here demonstrated further to improve biomolecular structure analysis.

Asymmetric Karplus curves for the protein side-chain 3J couplings

The standard Karplus equation for calculating 3J coupling constants from any given dihedral angle requires three empirical coefficients be determined that relate to the magnitudes of three modes of the angle dependency of 3J. Considering cosine modes only (bimodal, unimodal and baseline component), Karplus curves are generally symmetric with respect to the sign of the angle argument. Typically, their primary and secondary maxima differ in amplitude, whereas the two minima are of equal depth. However, chiral molecular topologies, such as those surrounding the main-chain and side-chain torsions in amino-acid residues, preclude, as regards substituent positioning, exact mirror-image conformations from being formed--for any given torsion-angle value. It is therefore unlikely that 3J couplings assume identical values for the corresponding positive and negative dihedral angles. This suggests that a better empirical fit of the torsion-angle dependency of 3J could be obtained when removing the constraint of symmetrically identical coupling constants. A sine term added to the Karplus equation allows independent modelling of both curve minima typically located near dihedral-angle values of +90 degrees and -90 degrees. Revisiting an extensive 3J coupling dataset previously recorded to determine the side-chain torsions chi1 in the protein flavodoxin, the asymmetric Karplus model accomplishes a more accurate fit to the experimental data. Asymmetries revealed in the angle dependencies exceed the experimental precision in determining 3J. Accounting for these effects helps improve molecular models.

Stereospecific assignments of protein NMR resonances based on the tertiary structure and 2D/3D NOE data

In many cases of protein structure determination by NMR a high-quality structure is required. An important contribution to structural precision is stereospecific assignment of magnetically nonequivalent prochiral methylene and methyl groups, eliminating the need for introducing pseudoatoms and pseudoatom corrections in distance restraint lists. Here, we introduce the stereospecific assignment program that uses the resonance assignment, a preliminary 3D structure and 2D and/or 3D nuclear Overhauser effect spectroscopy peak lists for stereospecific assignment. For each prochiral group the algorithm automatically calculates a score for the two different stereospecific assignment possibilities, taking into account the presence and intensity of the nuclear Overhauser effect (NOE) peaks that are expected from the local environment of each prochiral group (i.e., the close neighbors). The performance of the algorithm has been tested and used on NMR data of alpha-helical and beta-sheet proteins using homology models and/or X-ray structures. The program produced no erroneous stereospecific assignments provided the NOEs were carefully picked and the 3D model was sufficiently accurate. The set of NOE distance restraints produced by nmr2st using the results of the SSA module was superior in generating good-quality ensembles of NMR structures (low deviations from upper limits in conjunction with low root-mean-square-deviation values) in the first round of structure calculations. The program uses a novel approach that employs the entire 3D structure of the protein to obtain stereospecific assignment; it can be used to speed up the NMR structure refinement and to increase the quality of the final NMR ensemble even when no scalar or residual dipolar coupling information is available.

Interpreting dynamically-averaged scalar couplings in proteins

The experimental determination of scalar three-bond coupling constants represents a powerful method to probe both the structure and dynamics of proteins. The detailed structural interpretation of such coupling constants is usually based on Karplus relationships, which allow the measured couplings to be related to the torsion angles of the molecules. As the measured couplings are sensitive to thermal fluctuations, the parameters in the Karplus relationships are better derived from ensembles representing the distributions of dihedral angles present in solution, rather than from single conformations. We present a method to derive such parameters that uses ensembles of conformations determined through dynamic-ensemble refinement--a method that provides structural ensembles that simultaneously represent both the structure and the associated dynamics of a protein.

Side chain dynamics in unfolded protein states: an NMR based 2H spin relaxation study of delta131delta

NMR relaxation data on disordered proteins can provide insight into both structural and dynamic properties of these molecules. Because of chemical shift degeneracy in correlation spectra, detailed site-specific analyses of side chain dynamics have not been possible. Here, we present new experiments for the measurement of side chain dynamics in methyl-containing residues in unfolded protein states. The pulse schemes are similar to recently proposed methods for measuring deuterium spin relaxation rates in (13)CH(2)D methyl groups in folded proteins.(1) However, because resolution in (1)H-(13)C correlation maps of unfolded proteins is limiting, relaxation data are recorded as a series of (1)H-(15)N spectra. The methodology is illustrated with an application to the study of side chain dynamics in delta131delta, a large disordered fragment of staphylococcal nuclease containing residues 1-3 and 13-140 of the wide-type protein. A good correlation between the order parameters of the symmetry axes of the methyl groups and the backbone (1)H-(15)N bond vectors of the same residue is observed. Simulations establish that such a correlation is only possible if the unfolded state is comprised of an ensemble of structures which are not equiprobable. A motional model, which combines wobbling-in-a-cone and Gaussian axial fluctuations, is proposed to estimate chi(1) torsion angle fluctuations, sigma(chi)()1, of Val and Thr residues on the basis of the backbone and side chain order parameters. Values of sigma(chi)()1 are approximately 10 degrees larger than what has previously been observed in folded proteins. Of interest, the value of sigma(chi)()1 for Val 104 is considerably smaller than for other Val or Thr residues, suggesting that it may be part of a hydrophobic cluster. Notably large (15)N transverse relaxation rates are observed in this region. To our knowledge, this is the first time that side chain dynamics in an unfolded state have been studied in detail by NMR.

Self-consistent Karplus parametrization of 3J couplings depending on the polypeptide side-chain torsion chi1

Recently proposed self-consistent 3J coupling analysis (Schmidt, J. M.; Blümel, M.; Löhr, F.; Rüterjans, H. J. Biomol. NMR 1999, 14, 1-12) has been carried out to calibrate Karplus parameters constituting the empirical dependence of 3J coupling constants on the chi1 dihedral angle in amino acid side chains. The procedure involves simultaneous least-squares optimization of six sets of three Karplus coefficients related to all six 3J coupling types accessible in 15N,13C-labeled proteins. A simple concept of fundamental and incremental component couplings is proposed to account for substituent effects, eventually yielding amino acid topology-specific Karplus parameters. The method is exemplified with recombinant Desulfovibrio vulgaris flavodoxin (147 amino acids, 16 kDa) with reference to a total of 749 experimental 3JHalpha,Hbeta, 3JN',Hbeta, 3JC',Hbeta, 3JHalpha,Cgamma, 3JN',Cgamma, and 3JC',Cgamma coupling constants. Unlike other parametrizations, the present method does not make reference to X-ray coordinates, so that the Karplus coefficients obtained are not influenced by differences between solution and crystal states. Cross validation using X-ray torsion angles demonstrates the improvement relative to previous parametrizations. The Karplus coefficients derived are applicable to other proteins, too. Parameter refinement also yields a series of chi1 torsion angles, providing valuable constraints for protein structure determination, as well as optional parameters of local angular mobility in the contexts of Gaussian random fluctuation or a three-site jump model. The procedure permits automatic stereospecific assignments of Hbeta and Cgamma chemical shifts. The majority of the flavodoxin side-chain conformations agrees with high-resolution X-ray structures of the protein. Marked deviations between NMR and X-ray datasets are attributed to different rotameric states due to crystal-packing effects and to conformational equilibria between multiple chi1 rotamers.

A refined solution structure of hen lysozyme determined using residual dipolar coupling data

A high resolution NMR structure of hen lysozyme has been determined using 209 residual 1H-15N dipolar coupling restraints from measurements made in two different dilute liquid crystalline phases (bicelles) in conjunction with a data set of 1632 NOE distance restraints, 110 torsion angle restraints, and 60 hydrogen bond restraints. The ensemble of 50 low-energy calculated structures has an average backbone RMSD of 0.50+/-0.13A to the mean structure and of 1.49+/-0.10A to the crystal structure of hen lysozyme. To assess the importance of the dipolar coupling data in the structure determination, the final structures are compared with an ensemble calculated using an identical protocol but excluding the dipolar coupling restraints. The comparison shows that structures calculated with the dipolar coupling data are more similar to the crystal structure than those calculated without, and have better stereochemical quality. The structures also show improved quality factors when compared with additional dipolar coupling data that were not included in the structure calculations, with orientation-dependent 15N chemical shift changes measured in the bicelle solutions, and with T1/T2 values obtained from 15N relaxation measurements. Analysis of the ensemble of NMR structures and comparisons with crystal structures, 15N relaxation data, and molecular dynamics simulations of hen lysozyme provides a detailed description of the solution structure of this protein and insights into its dynamical behavior.

NMR and molecular mechanics study of pyrethrins I and II

Bioassay-directed fractionation of the organic extract of the Kenyan pyrethrum flowers (Chrysanthemum cinerariaefolium Vissiani) resulted in the isolation of two natural pyrethrin esters, pyrethrin I (PI) and pyrethrin II (PII) as the major constituents. These esters elicited inhibition of the multiple drug resistant (MDR) Mycobacterium tuberculosis. The high-field (1)H and (13)C nuclear magnetic resonance (NMR) chemical shifts of PI and PII were unequivocally assigned using modern two-dimensional (2D) proton-detected heteronuclear multiple-quantum coherence (HMQC) and heteronuclear multiple-bond correlation (HMBC) experiments. The conformations of both esters were deduced from (1)H-(1)H vicinal coupling constants and confirmed by 2D nuclear Overhauser effect spectroscopy (NOESY). Computer molecular modeling (MM) studies revealed that PI and PII molecules adopt a "love-seat" conformation in chloroform (CDCl(3)) solution.

Side-chain conformations in an unfolded protein: chi1 distributions in denatured hen lysozyme determined by heteronuclear 13C, 15N NMR spectroscopy

Using a 13C and 15N-labelled sample, multi-dimensional heteronuclear NMR techniques have been carried out to characterise hen lysozyme denatured in 8 M urea at pH 2.0. The measurement of 3J(C',Cgamma) and 3J(N,Cgamma) coupling constants has enabled side-chain chi1 torsion angle populations to be probed in the denatured polypeptide chain. Analysis of the coupling constant data has allowed the relative populations of the three staggered rotamers about chi1 to be defined for 51 residues. The amino acids can broadly be divided into five classes that show differing side-chain conformational preferences in the denatured state. These range from a strong preference for the -60 degrees chi1 rotamer for methionine and leucine (74-79 % population) to a favouring of the +60 degrees chi1 rotamer for threonine (67 % population). The differences in behaviour reflect the steric and electrostatic characteristics of the side-chains concerned. A close agreement is seen between the chi1 populations calculated from the experimental coupling constant data and predictions from the statistical model for a random coil that uses the chi1 torsion angle distributions in a data base of native protein structures. Short-range interactions therefore dominate in determining the local conformational properties of side-chains in a denatured protein. Deviations are, however, observed for many of the aromatic residues involved in hydrophobic clusters within the denatured protein. For these residues the effects of additional non-local interactions in the clusters presumably play a major role in determining the chi1 preferences.

Conformation in solution and dynamics of a structurally constrained linear insect kinin pentapeptide analogue

The preferred conformations of the active diuretic insect kinin pentapeptide analogue Phe-Phe-Aib-Trp-Gly-NH2 were studied using nmr spectroscopy and molecular modeling. Structure sets consistent with rotating frame nuclear Overhauser effect spectroscopy distance constraints obtained by restrained simulated annealing in vacuo indicate a predominant population of a type II beta-turn involving the Phe1-Trp4 region. An equilibrium between this type II and a type I beta-turn formed by residues Phe2 and Gly5 was observed in a 5 ns restrained molecular dynamics simulation using the implicit generalized Born solvent accessible surface area (GB/SA) solvation model. When subjected to 500 ps dynamics with explicit water both beta-turn folds were conserved throughout the simulations. The results obtained with implicit and explicit solvation models are compared, and their consistency with the nmr observations is discussed. The behavior of the linear pentapeptide in this study is in agreement with an earlier report on the consensus conformation of the insect kinin active core derived from analysis of cyclic active analogues.

New methylene specific experiments for the measurement of scalar spin-spin coupling constants between protons attached to 13C

New two- and three-dimensional NMR methods are proposed for the measurement of 3J(H, H) coupling constants between two adjacent methylene moieties. The new experiment, which is based on a combination of the E.COSY principle and double/zero quantum heteronuclear spectroscopy, has been applied to diaceton-glucose and to the protein rhodniin. The coupling constants of CH-CH2 groups have been compared with those obtained from a HCCH-E.COSY experiment to check the reliability of the results. An analysis of the coupling constants derived by comparison between experimental and simulated spectra is presented. Simulations were done with the program wtest considering fully correlated dipolar relaxation. Side-chain conformations in amino acids with adjacent methylene groups can be determined by the new experiment.

Side-chains in native and random coil protein conformations. Analysis of NMR coupling constants and chi1 torsion angle preferences

The behaviour of amino acid side-chains in proteins in solution has been characterised by analysing NMR 3JHalphaH beta coupling constants and crystallographic chi1 torsion angles. Side-chains both in the core of native folded proteins and in situations where there is an absence of close packing including the random coil state have been considered. An analysis of experimental 3JHalphaH beta coupling constant data for ten proteins shows that in the core of native proteins a very close similarity is observed between the chi1 conformations adopted in solution and in crystals. There is clear evidence, however, for significant motional averaging about the chi1 torsion angles in solution. Using a model of a Gaussian distribution about the average torsion angles the extent of these fluctuations has been quantified; the standard deviation for the motion is 26 degrees, the fluctuations about chi1 in the protein core being similar in size to those found for main-chain phi torsion angles in solution. From the distribution of chi1 torsion angles in a data base of protein crystal structures, torsion angle populations and coupling constants have been predicted for a random coil polypeptide. Significant variations in the chi1 distributions for different amino acids give differences in the predicted coupling constants; for 3JHalphaH beta, for example, values of 5.1 and 5.7 Hz are predicted for serine compared with 4.9 and 9.9 Hz for leucine. Experimental data for short unstructured peptides show an excellent agreement with the predictions, indicating that the overall chi1 distributions in protein crystals reflect the local preferences of the amino acids. Predictions from the protein data base therefore provide an important framework for interpreting experimental data for non-native protein conformations and for residues on the surface of folded proteins.

Conformational studies of paclitaxel analogs modified at the C-2' position in hydrophobic and hydrophilic solvent systems

The conformations of two paclitaxel analogs modified at the C-2' position, 2'-deoxypaclitaxel and 2'-methoxypaclitaxel, were studied in hydrophobic and hydrophilic solvent systems by a combination of NMR spectroscopy, CD measurements, and molecular modeling. Both analogs have hydrophobic and hydrophilic conformations that resemble those of paclitaxel itself in the same media. Since the two have diminished biological activities in a number of bioactivity assays and the hydrogen-bonding capability of the 2'-hydroxyl group has been eliminated, we postulate that this group is involved in hydrogen bonding with tubulin and plays an important role in molecular recognition. The results of this study are in agreement with our earlier report on paclitaxel 2'-acetate, an analog in which the 2'-hydroxyl group hydrogen-bonding capacity has also been eliminated.

Limits of NMR structure determination using variable target function calculations: ribonuclease T1, a case study

Limits of NMR structure determination using multidimensional NMR spectroscopy, variable target function calculations and relaxation matrix analysis were explored using the model protein ribonuclease T1 (RNase T1). The enzyme consists of 104 amino acid residues and has a molecular mass of approximately 11 kDa. Primary experimental data comprise 1856 assigned NOE intensities, 493 3J coupling constants and 62 values of amid proton exchange rates. From these data, 2580 distance bounds, 168 allowed ranges for torsional angles and stereospecific assignments for 75% of beta-methylene protons as well as for 80% of diastereotopic methyl groups were derived. Whenever possible, the distance restraints were refined in a relaxation matrix analysis including amid proton exchange data for improvement of lower distance limits. Description of side-chain conformations were based on various models of motional averaging of 3J coupling constants. The final structure ensemble was selected from the starting ensemble comparing the global precision of structures with order parameters derived from 15N relaxation time measurements. Significant differences between the structure of RNase T1 in solution and in the crystal became apparent from a comparison of the two highly resolved structures.

The solution structure refinement of the paramagnetic reduced high-potential iron-sulfur protein I from Ectothiorhodospira halophila by using stable isotope labeling and nuclear relaxation

The reduced high-potential iron sulfur protein I from Ectothiorhodospira halophila which contains the [4Fe-4S]2+ polymetallic center has been fully labeled with 15N and 13C. The protein is paramagnetic, the nuclear relaxation times of nuclei close to the paramagnetic ion are drastically shortened and some strategic dipolar connectivities are lost. Notwithstanding, the solution structure has been reported [Banci, L., Bertini, I., Eltis, L. D., Felli, I. C., Kastrau, D. H. W., Luchinat, C., Piccioli, M., Pierattelli, R. & Smith, M. (1994) Eur. J. Biochem. 225, 715-725]. We have performed classical HNHA, HNCA soft-COSY, soft-HCCH E. COSY and 15N-1H correlated NOESY experiments in order to obtain a set of 3J scalar coupling constants. Some experiments have been optimized to counterbalance the effect of paramagnetism. From heteronuclear single-quantum experiments preceded by a 180 degrees pulse and variable delay times, the non-selective magnetization recovery has been followed from which the contribution to dipolar relaxation of nuclei due to the interaction with the paramagnetic metal ions (rho para) has been estimated. Finally, the intensities of NOEs have been corrected for the presence of paramagnetic metal ions and these corrected values together with 3J values and rho para data have been used to obtain a well defined solution structure. The aim is that of obtaining a structure with enough constraints to be well resolved all over the protein, including the vicinity of the paramagnetic metal cluster, which is anchored to the protein through the rho para constraints. In total, 1226 corrected NOESY crosspeaks (of which 945 were found to be meaningful), 37 one-dimensional NOEs, 39 3JHNH alpha and 37 3JHNC' (providing 45 phi dihedral angle constraints) 54 3JH alpha H beta and 31 3JNH beta (providing 26 chi 1 dihedral angle constraints), 4 chi 2 dihedral angle constraints of the coordinated cysteines, obtained from the hyperfine shifts of the beta CH protons, and 58 rho para constraints, have been used for structure calculation. Restrained molecular dynamics simulations have also been performed to provide the final family of structures. This research demonstrates that stable isotope labeling provides specific advantages for the NMR investigation of paramagnetic molecules, as the small magnetic moment of heteronuclei minimizes the paramagnetic influence of unpaired electrons.

Structure calculation from NMR data

NMR calculation methods have kept pace with the rapid extension of NMR experiments to larger molecules. By including additional data and effects of local dynamics in the refinement, we can obtain a more complete picture of the molecule in solution. The structure determination process is being aided by new methods to solve some aspects of spectral assignment during the structure calculation.

Determination of heteronuclear three-bond J-coupling constants in peptides by a simple heteronuclear relayed E.COSY experiment

A simple heteronuclear relayed E.COSY pulse sequence with a minimum number of pulses is proposed for the quantitative determination of heteronuclear three-bond J-coupling constants in uniformly 13C-enriched polypeptide samples. Numerous heteronuclear three-bond coupling constants, including 3JHNC, 3JHNC beta, 3JH beta C, and 3JH alpha C gamma, can be determined for each residue from a single heteronuclear relayed E.COSY spectrum. Couplings relevant for stereospecific assignments as well as for the determination of dihedral angles in the amino acid backbone and in side chains are obtained. The method is demonstrated on the uniformly 13C-enriched decapeptide antamanide (-Val1-Pro2-Pro3-Ala4-Phe5-Phe6-Pro7-Pro8-Phe9-Phe1 0-).