An empirical correlation between 1JC.alpha.H.alpha. and protein backbone conformation

Quantitative evaluation of positive ϕ angle propensity in flexible regions of proteins from three-bond J couplings

(3)JHNHα and (3)JC'C' couplings can be readily measured in isotopically enriched proteins and were shown to contain precise information on the backbone torsion angles, ϕ, sampled in disordered regions of proteins. However, quantitative interpretation of these couplings required the population of conformers with positive ϕ angles to be very small. Here, we demonstrate that this restriction can be removed by measurement of (3)JC'Hα values. Even though the functional forms of the (3)JC'Hα and (3)JHNHα Karplus equations are the same, large differences in their coefficients enable accurate determination of the fraction of time that positive ϕ angles are sampled. A four-dimensional triple resonance HACANH[C'] E.COSY experiment is introduced to simultaneously measure (3)JC'Hα and (3)JHNC' in the typically very congested spectra of disordered proteins. High resolution in these spectra is obtained by non-uniform sampling (in the 0.1-0.5% range). Application to the intrinsically disordered protein α-synuclein shows that while most residues have close-to-zero positive ϕ angle populations, up to 16% positive ϕ population is observed for Asn residues. Positive ϕ angle populations determined with the new approach agree closely with consensus values from protein coil libraries and prior analysis of a large set of other NMR parameters. The combination of (3)JHNC' and (3)JC'C' provides information about the amplitude of ϕ angle dynamics.

A maximum entropy approach to the study of residue-specific backbone angle distributions in α-synuclein, an intrinsically disordered protein

α-Synuclein is an intrinsically disordered protein of 140 residues that switches to an α-helical conformation upon binding phospholipid membranes. We characterize its residue-specific backbone structure in free solution with a novel maximum entropy procedure that integrates an extensive set of NMR data. These data include intraresidue and sequential H(N) − H(α) and H(N) − H(N) NOEs, values for (3) JHNHα, (1) JHαCα, (2) JCαN, and (1) JCαN, as well as chemical shifts of (15)N, (13)C(α), and (13)C' nuclei, which are sensitive to backbone torsion angles. Distributions of these torsion angles were identified that yield best agreement to the experimental data, while using an entropy term to minimize the deviation from statistical distributions seen in a large protein coil library. Results indicate that although at the individual residue level considerable deviations from the coil library distribution are seen, on average the fitted distributions agree fairly well with this library, yielding a moderate population (20-30%) of the PPII region and a somewhat higher population of the potentially aggregation-prone β region (20-40%) than seen in the database. A generally lower population of the αR region (10-20%) is found. Analysis of (1)H − (1)H NOE data required consideration of the considerable backbone diffusion anisotropy of a disordered protein.

An introduction to biological NMR spectroscopy

NMR spectroscopy is a powerful tool for biologists interested in the structure, dynamics, and interactions of biological macromolecules. This review aims at presenting in an accessible manner the requirements and limitations of this technique. As an introduction, the history of NMR will highlight how the method evolved from physics to chemistry and finally to biology over several decades. We then introduce the NMR spectral parameters used in structural biology, namely the chemical shift, the J-coupling, nuclear Overhauser effects, and residual dipolar couplings. Resonance assignment, the required step for any further NMR study, bears a resemblance to jigsaw puzzle strategy. The NMR spectral parameters are then converted into angle and distances and used as input using restrained molecular dynamics to compute a bundle of structures. When interpreting a NMR-derived structure, the biologist has to judge its quality on the basis of the statistics provided. When the 3D structure is a priori known by other means, the molecular interaction with a partner can be mapped by NMR: information on the binding interface as well as on kinetic and thermodynamic constants can be gathered. NMR is suitable to monitor, over a wide range of frequencies, protein fluctuations that play a crucial role in their biological function. In the last section of this review, intrinsically disordered proteins, which have escaped the attention of classical structural biology, are discussed in the perspective of NMR, one of the rare available techniques able to describe structural ensembles. This Tutorial is part of the International Proteomics Tutorial Programme (IPTP 16 MCP).

Left-handed helical preference in an achiral peptide chain is induced by an L-amino acid in an N-terminal type II β-turn

Oligomers of the achiral amino acid Aib adopt helical conformations in which the screw-sense may be controlled by a single N-terminal residue. Using crystallographic and NMR techniques, we show that the left- or right-handed sense of helical induction arises from the nature of the β-turn at the N terminus: the tertiary amino acid L-Val induces a left-handed type II β-turn in both the solid state and in solution, while the corresponding quaternary amino acid L-α-methylvaline induces a right-handed type III β-turn.

One-bond and two-bond J couplings help annotate protein secondary-structure motifs: J-coupling indexing applied to human endoplasmic reticulum protein ERp18

NMR coupling constants, both direct one-bond ((1)J) and geminal two-bond ((2)J), are employed to analyze the protein secondary structure of human oxidized ERp18. Coupling constants collected and evaluated for the 18 kDa protein comprise 1268 values of (1)J(CαHα), (1)J(CαCβ), (1)J(CαC'), (1)J(C'N'), (1)J(N'Cα), (1)J(N') (HN), (2)J(CαN'), (2)J(HNCα), (2)J(C'HN), and (2)J(HαC'). Comparison with (1)J and (2)J data from reference proteins and pattern analysis on a per-residue basis permitted main-chain ϕ,ψ torsion-angle combinations of many of the 149 amino-acid residues in ERp18 to be narrowed to particular secondary-structure motifs. J-coupling indexing is here being developed on statistical criteria and used to devise a ternary grid for interpreting patterns of relative values of J. To account for the influence of the varying substituent pattern in different amino-acid sidechains, a table of residue-type specific threshold values was compiled for discriminating small, medium, and large categories of J. For the 15-residue insertion that distinguishes the ERp18 fold from that of thioredoxin, the J-coupling data hint at a succession of five isolated Type-I β turns at progressively shorter sequence intervals, in agreement with the crystal structure.

Correlation of (2)J couplings with protein secondary structure

Geminal two-bond couplings ((2)J) in proteins were analyzed in terms of correlation with protein secondary structure. NMR coupling constants measured and evaluated for a total six proteins comprise 3999 values of (2)J(CalphaN'), (2)J(C'HN), (2)J(HNCalpha), (2)J(C'Calpha), (2)J(HalphaC'), (2)J(HalphaCalpha), (2)J(CbetaC'), (2)J(N'Halpha), (2)J(N'Cbeta), and (2)J(N'C'), encompassing an aggregate 969 amino-acid residues. A seamless chain of pattern comparisons across the spectrum datasets recorded allowed the absolute signs of all (2)J coupling constants studied to be retrieved. Grouped by their mediating nucleus, C', N' or C(alpha), (2)J couplings related to C' and N' depend significantly on phi,psi torsion-angle combinations. beta turn types I, I', II and II', especially, can be distinguished on the basis of relative-value patterns of (2)J(CalphaN'), (2)J(HNCalpha), (2)J(C'HN), and (2)J(HalphaC'). These coupling types also depend on planar or tetrahedral bond angles, whereas such dependences seem insignificant for other types. (2)J(HalphaCbeta) appears to depend on amino-acid type only, showing negligible correlation with torsion-angle geometry. Owing to its unusual properties, (2)J(CalphaN') can be considered a "one-bond" rather than two-bond interaction, the allylic analog of (1)J(N'Calpha), as it were. Of all protein J coupling types, (2)J(CalphaN') exhibits the strongest dependence on molecular conformation, and among the (2)J types, (2)J(HNCalpha) comes second in terms of significance, yet was hitherto barely attended to in protein structure work.

A strained DNA binding helix is conserved for site recognition, folding nucleation, and conformational modulation

Nucleic acid recognition is often mediated by alpha-helices or disordered regions that fold into alpha-helix on binding. A peptide bearing the DNA recognition helix of HPV16 E2 displays type II polyproline (PII) structure as judged by pH, temperature, and solvent effects on the CD spectra. NMR experiments indicate that the canonical alpha-helix is stabilized at the N-terminus, while the PII forms at the C-terminus half of the peptide. Re-examination of the dihedral angles of the DNA binding helix in the crystal structure and analysis of the NMR chemical shift indexes confirm that the N-terminus half is a canonical alpha-helix, while the C-terminal half adopts a 3(10) helix structure. These regions precisely match two locally driven folding nucleii, which partake in the native hydrophobic core and modulate a conformational switch in the DNA binding helix. The peptide shows only weak and unspecific residual DNA binding, 10(4)-fold lower affinity, and 500-fold lower discrimination capacity compared with the domain. Thus, the precise side chain conformation required for modulated and tight physiological binding by HPV E2 is largely determined by the noncanonical strained alpha-helix conformation, "presented" by this unique architecture. (c) 2009 Wiley Periodicals, Inc. Biopolymers 91: 432-443, 2009.

Deuterium isotope effects on 15N backbone chemical shifts in proteins

Quantum mechanical calculations are presented that predict that one-bond deuterium isotope effects on the (15)N chemical shift of backbone amides of proteins, (1)Delta(15)N(D), are sensitive to backbone conformation and hydrogen bonding. A quantitative empirical model for (1)Delta(15)N(D) including the backbone dihedral angles, Phi and Psi, and the hydrogen bonding geometry is presented for glycine and amino acid residues with aliphatic side chains. The effect of hydrogen bonding is rationalized in part as an electric-field effect on the first derivative of the nuclear shielding with respect to N-H bond length. Another contributing factor is the effect of increased anharmonicity of the N-H stretching vibrational state upon hydrogen bonding, which results in an altered N-H/N-D equilibrium bond length ratio. The N-H stretching anharmonicity contribution falls off with the cosine of the N-H...O bond angle. For residues with uncharged side chains a very good prediction of isotope effects can be made. Thus, for proteins with known secondary structures, (1)Delta(15)N(D) can provide insights into hydrogen bonding geometries.

Variation in protein C(alpha)-related one-bond J couplings

Four types of polypeptide (1)J(C alpha X) couplings are examined, involving the main-chain carbon C(alpha) and either of four possible substituents. A total 3105 values of (1)J(C alpha H alpha), (1)J(C alpha C beta), (1)J(C alpha C'), and (1)J(C alpha N') were collected from six proteins, averaging 143.4 +/- 3.3, 34.9 +/- 2.5, 52.6 +/- 0.9, and 10.7 +/- 1.2 Hz, respectively. Analysis of variances (ANOVA) reveals a variety of factors impacting on (1)J and ranks their relative statistical significance and importance to biomolecular NMR structure refinement. Accordingly, the spread in the (1)J values is attributed, in equal proportions, to amino-acid specific substituent patterns and to polypeptide-chain geometry, specifically torsions phi, psi, and chi(1) circumjacent to C(alpha). The (1)J coupling constants correlate with protein secondary structure. For alpha-helical phi, psi combinations, (1)J(C alpha H alpha) is elevated by more than one standard deviation (147.8 Hz), while both (1)J(C alpha N') and (1)J(C alpha C beta) fall short of their grand means (9.5 and 33.7 Hz). Rare positive phi torsion angles in proteins exhibit concomitant small (1)J(C alpha H alpha) and (1)J(C alpha N') (138.4 and 9.6 Hz) and large (1)J(C alpha C beta) (39.9 Hz) values. The (1)J(C alpha N') coupling varies monotonously over the phi torsion range typical of beta-sheet secondary structure and is largest (13.3 Hz) for phi around -160 degrees. All four coupling types depend on psi and thus help determine a torsion that is notoriously difficult to assess by traditional approaches using (3)J. Influences on (1)J stemming from protein secondary structure and other factors, such as amino-acid composition, are largely independent.

Electric field effects on one-bond indirect spin-spin coupling constants and possible biomolecular perspectives

Electric field (EF) induced changes of one-bond indirect spin-spin coupling constants are investigated on a wide range of molecules including peptide models. EFs were both externally applied and internally calculated without external EF application by the hybrid density functional theory method. Reliable agreement with experimental data has been obtained for calculated one-bond J-couplings. The role of the EF sign and direction, internal and induced components, hydrogen bonding, internuclear distance and hyperconjugative interactions on the one-bond J-coupling vs EF interconnection is analyzed. A linear dependence of 1J on EF projection along the bond is obtained, if the bound atoms possess different enough electron densities and an EF determined by the electronic polarization exists along the bond. Accentuating the 1JNH couplings as possible EF sensitive parameters, a systematic study is done in two sets of molecules with a large variation of the native internal EF value. The most EF affected component of the 1JNH coupling constant is the spin-dipole term of Ramsey's formulation; however, in the total J-coupling formation, the EF influence on the Fermi contact term is the most significant. The induced EF projection along the bond is 6.7 times weaker in magnitude than the simulated external uniform field. The absolute EF dependence of the one-bond J-coupling involves only the internal field, which is the sum of the induced field (if the external field exists) and the internuclear field determined by the native polarization. That linear and universal dependence joins the corresponding couplings in a diverse set of molecules under various electrostatic conditions. Many types of the one-bond J-couplings can be potentially measured in biomolecules, and the study of their relation with the electrostatic properties at the corresponding sites opens a new avenue to the full exploitation of the NMR measurable parameters with novel and exciting applications.

J-GFT NMR for precise measurement of mutually correlated nuclear spin-spin couplings

G-matrix Fourier transform (GFT) NMR spectroscopy is presented for accurate and precise measurement of chemical shifts and nuclear spin-spin couplings correlated according to spin system. The new approach, named "J-GFT NMR", is based on a largely extended GFT NMR formalism and promises to have a broad impact on projection NMR spectroscopy. Specifically, constant-time J-GFT (6,2)D (HA-CA-CO)-N-HN was implemented for simultaneous measurement of five mutually correlated NMR parameters, that is, 15N backbone chemical shifts and the four one-bond spin-spin couplings 13Calpha-1Halpha, 13Calpha-13C', 15N-13C', and 15N-1HNu. The experiment was applied for measuring residual dipolar couplings (RDCs) in an 8 kDa protein Z-domain aligned with Pf1 phages. Comparison with RDC values extracted from conventional NMR experiments reveals that RDCs are measured with high precision and accuracy, which is attributable to the facts that (i) the use of constant time evolution ensures that signals do not broaden whenever multiple RDCs are jointly measured in a single dimension and (ii) RDCs are multiply encoded in the multiplets arising from the joint sampling. This corresponds to measuring the couplings multiple times in a statistically independent manner. A key feature of J-GFT NMR, i.e., the correlation of couplings according to spin systems without reference to sequential resonance assignments, promises to be particularly valuable for rapid identification of backbone conformation and classification of protein fold families on the basis of statistical analysis of dipolar couplings.

Metal-binding and nuclease activity of an antimicrobial peptide analogue of the salivary histatin 5

The salivary antimicrobial peptide histatin 5 is characterized by its cationic nature, structural flexibility, and the presence of two metal-binding sites (the ATCUN motif and a Zn-binding motif). These properties make this peptide a good model for the design of new drugs of low molecular weight. In this work, we have synthesized and studied a new peptide, an analogue of the histatin 5 named ATCUN-C16, which contains both metal-binding centers. The results show that our 20-residue-derived peptide preserves anticandidal activity and exhibits a higher propensity to assume a stable conformation in a hydrophobic environment than do histatin 5 and the C16 peptide that contains the 16 residues of the C-terminal part of histatin 5, although overall our peptide remains a flexible molecule. ACTUN-C16 was found to bind DNA in a gel retardation assay and to have a nuclease activity in the presence of copper and zinc ions and ascorbate. Its nuclease activity can be attributed to the synergistic action of oxidative and hydrolytic activities due to the Cu-ATCUN complex and to the zinc ion coordination, respectively. The results show a new property of this family of salivary peptides and suggest a novel use of this peptide as a small nuclease and biotechnological tool.

Hyperphosphorylation of tau induces local polyproline II helix

Alzheimer's disease is characterized by two protein precipitates, extracellular amyloid plaques and intracellular neurofibrillary tangles (NFTs). The primary constituent of NFTs is a hyperphosphorylated form of the microtubule-binding protein tau. Hyperphosphorylation of tau on over 30 residues, primarily within proline-rich sequences, is associated with conformational changes whose nature is poorly defined. Peptides derived from the proline-rich region of tau (residues 174-242) were synthesized, and the conformations were analyzed for the nonphosphorylated and phosphorylated peptides. CD and NMR data indicate that phosphorylation of serine and threonine residues in proline-rich sequences induces a conformational change to a type II polyproline helix. The largest phosphorylation-dependent conformational changes observed by CD were for tau peptides incorporating residues 174-183 or residues 229-238. Phosphoserine and phosphothreonine residues exhibited ordered values of (3)J(alphaN) (3.1-6.2 Hz; mean = 4.7 Hz) compared to nonphosphorylated serine and threonine. Phosphorylation of a tau peptide consisting of tau residues 196-209 resulted in the disruption of a nascent alpha-helix. These results suggest that global reorganization of tau may occur upon hyperphosphorylation of proline-rich sequences in tau.

Electronic control of amide cis-trans isomerism via the aromatic-prolyl interaction

The cis-trans isomerization of prolyl amide bonds results in large structural and functional changes in proteins and is a rate-determining step in protein folding. We describe a novel electronic strategy to control cis-trans isomerization, based on the demonstration that interactions between aromatic residues and proline are tunable by aromatic electronics. A series of peptides of sequence TXPN, X = Trp, pyridylalanine, pentafluorophenylalanine, or 4-Z-phenylalanine derivatives (Z = electron-donating, electron-withdrawing, or electron-neutral substituents), was synthesized and Ktrans/cis analyzed by NMR. Electron-rich aromatic residues stabilized cis amide bond formation, while electron-poor aromatics relatively favored trans amide bond formation. A Hammett correlation between aromatic electronics and cis-trans isomerization was observed. These results indicate that the interaction between aromatic residues and proline, which is observed to stabilize cis amide bonds and is also a general stabilizing interaction ubiquitous in proteins and protein-protein complexes, is not stabilized exclusively by a classical hydrophobic effect. To a large extent, the aromatic-prolyl interaction is driven and controllable by an electronic effect between the aromatic ring pi-electrons and the proline ring, consistent with a C-H-pi interaction as the key stabilizing force. The aromatic-prolyl interaction is electronically tunable by 0.9 kcal/mol and is enthalpic in nature. In addition, by combining aromatic ring electronics and stereoelectronic effects using 4-fluoroprolines, we demonstrate broad tuning (2.0 kcal/mol) of cis-trans isomerism in tetrapeptides. We demonstrate a simple tetrapeptide, TWflpN, that exhibits 60% cis amide bond and adopts a type VIa1 beta-turn conformation.

Molecular basis for phosphorylation-dependent, PEST-mediated protein turnover

Proteasomal-mediated rapid turnover of proteins is often modulated by phosphorylation of PEST sequences. The E2 protein from papillomavirus participates in gene transcription, DNA replication, and episomal genome maintenance. Phosphorylation of a PEST sequence located in a flexible region accelerates its degradation. NMR analysis of a 29 amino acid peptide fragment derived from this region shows pH-dependent polyproline II and alpha helix structures, connected by a turn. Phosphorylation, in particular that at serine 301, disrupts the overall structure, and point mutations have either stabilizing or destabilizing effects. There is an excellent correlation between the thermodynamic stability of different peptides and the half-life of E2 proteins containing the same mutations in vivo. The structure around the PEST region appears to have evolved a marginal stability that is finely tunable by phosphorylation. Thus, conformational stability, rather than recognition of a phosphate modification, modulates the degradation of this PEST sequence by the proteasome machinery.

Fast and accurate computation of the 13C chemical shifts for an alanine-rich peptide

The purpose of this work is, first, to present a fast and accurate technique to compute Boltzmann-averaged values of the quantum-chemical 13C chemical shifts for each amino acid in oligopeptides, demonstrated here by an application to the peptide Ac-XXAAAAAAAOO-NH2 (where X denotes diaminobutyric acid, A is alanine, and O is ornithine) [XAO] and, second, to discuss the capability of the 13Calpha and 13Cbeta chemical shifts to distinguish the PP(II) conformation from the alpha-helix and statistical-coil conformations. Use is made of a combination of approaches, summarized as follows: (1) derivation of an ensemble of conformations by using a molecular mechanics technique; (2) use of a clustering procedure to form families and build a reduced set of conformations consisting of the lowest-energy conformations of each family, and (3) computation of the 13C chemical shifts for the lowest-energy conformations of each family, using a quantum-chemical approach that treats a selected residue, or group of residues, with a 6-311+G(2d,p) locally-dense basis set, while the remaining residues in the sequence are treated with a 3-21G basis set. The whole procedure is quite accurate and speeds up the computation of the Boltzmann-averaged values of the 13C-chemical shifts by several orders of magnitude. The present application sheds some light on the conformational preference for alanine and non-alanine residues to occupy the PP(II) helical region of the Ramachandran map.

The NMR structure of the sensory domain of the membranous two-component fumarate sensor (histidine protein kinase) DcuS of Escherichia coli

The structure of the water-soluble, periplasmic domain of the fumarate sensor DcuS (DcuS-pd) has been determined by NMR spectroscopy in solution. DcuS is a prototype for a sensory histidine kinase with transmembrane signal transfer. DcuS belongs to the CitA family of sensors that are specific for sensing di- and tricarboxylates. The periplasmic domain is folded autonomously and shows helices at the N and the C terminus, suggesting direct linking or connection to helices in the two transmembrane regions. The structure constitutes a novel fold. The nearest structural neighbor is the Per-Arnt-Sim domain of the photoactive yellow protein that binds small molecules covalently. Residues Arg107, His110, and Arg147 are essential for fumarate sensing and are found clustered together. The structure constitutes the first periplasmic domain of a two component sensory system and is distinctly different from the aspartate sensory domain of the Tar chemotaxis sensor.

NMR identification of left-handed polyproline type II helices

NMR characteristics of a model left-handed 3(1)-helical peptide are reported in this study. With temperature and sequence corrections on the predicted random coil (15)N chemical shifts, a significant (15)N chemical shift deviation is observed for the model 3(1) peptide. The (15)N chemical shift differences also correlate well with the molar ellipticities (at 220 nm) of the CD spectra at different temperatures, indicating that the (15)N chemical shift is a sensitive probe for 3(1)-helices. The average (3)J(HNalpha) and (1)J(CalphaHalpha) values of the model peptide are determined to be 6.5 and 142.6 Hz, respectively, which are consistent with the values calculated from the geometry of 3(1)-helices. With careful measurements of amide (15)N chemical shifts and incorporating temperature and sequence effect corrections, the (15)N chemical shifts can be used together with (3)J(HNalpha) and (1)J(CalphaHalpha) to differentiate 3(1)-helices from random coils with high confidence. Based on the observed NMR characteristics, a strategy is developed for probing left-handed 3(1)-helical structures from other secondary structures.

Measurement of small one-bond proton-carbon residual dipolar coupling constants in partially oriented (13)C natural abundance oligosaccharide samples: analysis of heteronuclear (1)J(CH)-modulated spectra with the BIRD inversion pulse

Two 2D J-modulated HSQC-based experiments were designed for precise determination of small residual dipolar one-bond carbon-proton coupling constants in (13)C natural abundance carbohydrates. Crucial to the precision of a few hundredths of Hz achieved by these methods was the use of long modulation intervals and BIRD pulses, which acted as semiselective inversion pulses. The BIRD pulses eliminated effective evolution of all but (1)J(CH) couplings, resulting in signal modulation that can be described by simple modulation functions. A thorough analysis of such modulation functions for a typical four-spin carbohydrate spin system was performed for both experiments. The results showed that the evolution of the (1)H-(1)H and long-range (1)H-(13)C couplings during the BIRD pulses did not necessitate the introduction of more complicated modulation functions. The effects of pulse imperfections were also inspected. While weakly coupled spin systems can be analyzed by simple fitting of cross peak intensities, in strongly coupled spin systems the evolution of the density matrix needs to be considered in order to analyse data accurately. However, if strong coupling effects are modest the errors in coupling constants determined by the "weak coupling" analysis are of similar magnitudes in oriented and isotropic samples and are partially cancelled during dipolar coupling calculation. Simple criteria have been established as to when the strong coupling treatment needs to be invoked.

Solution structure and backbone dynamics of an engineered arginine-rich subdomain 2 of the hepatitis C virus NS3 RNA helicase

The NS3 protein of the hepatitis C virus (HCV) is a 631 amino acid residue bifunctional enzyme with a serine protease localized to the N-terminal 181 residues and an RNA helicase located in the C-terminal 450 residues. The HCV NS3 RNA helicase consists of three well-defined subdomains which all contribute to its helicase activity. The second subdomain of the HCV helicase is flexibly linked to the remainder of the NS3 protein and could undergo rigid-body movements during the unwinding of double-stranded RNA. It also contains several motifs that are implicated in RNA binding and in coupling NTP hydrolysis to nucleic acid unwinding and translocation. As part of our efforts to use NMR techniques to assist in deciphering the enzyme's structure-function relationships and developing specific small molecule inhibitors, we have determined the solution structure of an engineered subdomain 2 of the NS3 RNA helicase of HCV, d(2Delta)-HCVh, and studied the backbone dynamics of this protein by (15)N-relaxation experiments using a model-free approach. The NMR studies on this 142-residue construct reveal that overall subdomain 2 of the HCV helicase is globular and well structured in solution even in the absence of the remaining parts of the NS3 protein. Its solution structure is very similar to the corresponding parts in the X-ray structures of the HCV NS3 helicase domain and intact bifunctional HCV NS3 protein. Slow hydrogen-deuterium exchange rates map to a well-structured, stable hydrophobic core region away from the subdomain interfaces. In contrast, the regions facing the subdomain interfaces in the HCV NS3 helicase domain are less well structured in d(2Delta)-HCVh, show fast hydrogen-deuterium exchange rates, and the analysis of the dynamic properties of d(2Delta)-HCVh reveals that these regions of the protein show distinct dynamical features. In particular, residues in motif V, which may be involved in transducing allosteric effects of nucleotide binding and hydrolysis on RNA binding, exhibit slow conformational exchange on the milli- to microsecond time-scale. The intrinsic conformational flexibility of this loop region may facilitate conformational changes required for helicase function.

Solution structure of DinI provides insight into its mode of RecA inactivation

The Escherichia coli RecA protein triggers both DNA repair and mutagenesis in a process known as the SOS response. The 81-residue E. coli protein DinI inhibits activity of RecA in vivo. The solution structure of DinI has been determined by multidimensional triple resonance NMR spectroscopy, using restraints derived from two sets of residual dipolar couplings, obtained in bicelle and phage media, supplemented with J couplings and a moderate number of NOE restraints. DinI has an alpha/beta fold comprised of a three-stranded beta-sheet and two alpha-helices. The beta-sheet topology is unusual: the central strand is flanked by a parallel and an antiparallel strand and the sheet is remarkably flat. The structure of DinI shows that six negatively charged Glu and Asp residues on DinI's kinked C-terminal alpha-helix form an extended, negatively charged ridge. We propose that this ridge mimics the electrostatic character of the DNA phospodiester backbone, thereby enabling DinI to compete with single-stranded DNA for RecA binding. Biochemical data confirm that DinI is able to displace ssDNA from RecA.

High-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a hydroxamic acid inhibitor

The high-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a sulfonamide derivative of a hydroxamic acid compound (WAY-151693) has been determined by multidimensional heteronuclear NMR. A total of 30 structures were calculated for residues 7-164 by means of hybrid distance geometry-simulated annealing using a total of 3280 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for the 30 structures is 0.43(+/-0.05) A for the backbone atoms, 0.80(+/-0.09) A for all atoms, and 0.47(+/-0.04) A for all atoms excluding disordered side-chains. The overall structure of MMP-13 is composed of a beta-sheet consisting of five beta-strands in a mixed parallel and anti-parallel arrangement and three alpha-helices where its overall fold is consistent with previously solved MMP structures. A comparison of the NMR structure of MMP-13 with the published 1.6 A resolution X-ray structure indicates that the major differences between the structures is associated with loop dynamics and crystal-packing interactions. The side-chains of some active-site residues for the NMR and X-ray structures of MMP-13 adopt distinct conformations. This is attributed to the presence of unique inhibitors in the two structures that encounter distinct interactions with MMP-13. The major structural difference observed between the MMP-13 and MMP-1 NMR structures is the relative size and shape of the S1' pocket where this pocket is significantly longer for MMP-13, nearly reaching the surface of the protein. Additionally, MMP-1 and MMP-13 exhibit different dynamic properties for the active-site loop and the structural Zn-binding region. The inhibitor WAY-151693 is well defined in the MMP-13 active-site based on a total of 52 distance restraints. The binding motif of WAY-151693 in the MMP-13 complex is consistent with our previously reported MMP-1:CGS-27023A NMR structure and is similar to the MMP-13: RS-130830 X-ray structure.

NMR structure of free RGS4 reveals an induced conformational change upon binding Galpha

Heterotrimeric guanine nucleotide-binding proteins (G-proteins) are transducers in many cellular transmembrane signaling systems where regulators of G-protein signaling (RGS) act as attenuators of the G-protein signal cascade by binding to the Galpha subunit of G-proteins (G(i)(alpha)(1)) and increasing the rate of GTP hydrolysis. The high-resolution solution structure of free RGS4 has been determined using two-dimensional and three-dimensional heteronuclear NMR spectroscopy. A total of 30 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 2871 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for residues 5-134 for the 30 structures is 0.47 +/- 0.05 A for the backbone atoms, 0. 86 +/- 0.05 A for all atoms, and 0.56 +/- 0.04 A for all atoms excluding disordered side chains. The NMR solution structure of free RGS4 suggests a significant conformational change upon binding G(i)(alpha)(1) as evident by the backbone atomic rms difference of 1. 94 A between the free and bound forms of RGS4. The underlying cause of this structural change is a perturbation in the secondary structure elements in the vicinity of the G(i)(alpha)(1) binding site. A kink in the helix between residues K116-Y119 is more pronounced in the RGS4-G(i)(alpha)(1) X-ray structure relative to the free RGS4 NMR structure, resulting in a reorganization of the packing of the N-terminal and C-terminal helices. The presence of the helical disruption in the RGS4-G(i)(alpha)(1) X-ray structure allows for the formation of a hydrogen-bonding network within the binding pocket for G(i)(alpha)(1) on RGS4, where RGS4 residues D117, S118, and R121 interact with residue T182 from G(i)(alpha)(1). The binding pocket for G(i)(alpha)(1) on RGS4 is larger and more accessible in the free RGS4 NMR structure and does not present the preformed binding site observed in the RGS4-G(i)(alpha)(1) X-ray structure. This observation implies that the successful complex formation between RGS4 and G(i)(alpha)(1) is dependent on both the formation of the bound RGS4 conformation and the proper orientation of T182 from G(i)(alpha)(1). The observed changes for the free RGS4 NMR structure suggest a mechanism for its selectivity for the Galpha-GTP-Mg(2+) complex and a means to facilitate the GTPase cycle.

NMR structure of tissue inhibitor of metalloproteinases-1 implicates localized induced fit in recognition of matrix metalloproteinases

A high quality solution structure of the matrix metalloproteinase inhibitory N-terminal domain of recombinant human tissue inhibitor of metalloproteinases-1 (N-TIMP-1) has been determined. For the rigidly packed residues, the average RMSD to the mean structure is 0. 57 A for the backbone atoms and 1.00 A for all heavy atoms. Comparison of the solution structure of free N-TIMP-1 with the crystal structure of TIMP-1 bound to the catalytic domain of MMP-3 ( Gomis-R]uth et al., 1997 ) shows that the structural core of the beta barrel flanked by helices is nearly unchanged by the association with MMP-3, evident from a backbone RMSD of 1.15 A. However, clear differences in the conformation of the MMP-binding ridge of free and MMP-bound TIMP-1 suggest induced fit throughout the ridge. The MMP-dependent conformational changes in the ridge include a dramatic bending of AB loop residues Glu28 through Leu34, moderate hinge bending of the CD-loop about residues Ala65 and Cys70, and modest bending of the Cys1 through Pro6 segment. A large number of interresidue Nuclear Overhauser enhancements (NOEs) augmented by stereospecific assignments, torsion restraints, and dipolar couplings (an average of 18 non-trivial restraints per residue) engender confidence in these structural inferences. A tight cluster of three lysine residues and one arginine residue atop beta-strands A and B, and identical among TIMP sequences, form the heart of a highly conserved electropositive patch that may interact with anionic components of the extracellular matrix.

High-resolution solution structure of the inhibitor-free catalytic fragment of human fibroblast collagenase determined by multidimensional NMR

The high-resolution solution structure of the inhibitor-free catalytic fragment of human fibroblast collagenase (MMP-1), a protein of 18.7 kDa, which is a member of the matrix metalloproteinase family, has been determined using three-dimensional heteronuclear NMR spectroscopy. A total of 30 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 3333 experimental NMR restraints, consisting of 2409 approximate interproton distance restraints, 84 distance restraints for 42 backbone hydrogen bonds, 426 torsion angle restraints, 125 3JNH alpha restraints, 153 C alpha restraints, and 136 C beta restraints. The atomic rms distribution about the mean coordinate positions for the 30 structures for residues 7-137 and 145-163 is 0.42 +/- 0.04 A for the backbone atoms, 0.80 +/- 0.04 A for all atoms, and 0.50 +/- 0.03 A for all atoms excluding disordered side chains. The overall structure of MMP-1 is composed of a beta-sheet consisting of five beta-strands in a mixed parallel and anti-parallel arrangement and three alpha-helices. A best-fit superposition of the NMR structure of inhibitor-free MMP-1 with the 1.56 A resolution X-ray structure by Spurlino et al. [Spurlino, J. C., Smallwood, A. M., Carlton, D. D., Banks, T. M., Vavra, K. J., Johnson, J. S., Cook, E. R., Falvo, J., and Wahl, R. C., et al. (1994) Proteins: Struct., Funct., Genet. 19, 98-109] complexed with a hydroxamate inhibitor yields a backbone atomic rms difference of 1.22 A. The majority of differences between the NMR and X-ray structure occur in the vicinity of the active site for MMP-1. This includes an increase in mobility for residues 138-144 and a displacement for the Ca(2+)-loop (residues 74-80). Distinct differences were observed for side-chain torsion angles, in particular, the chi 1 for N80 is -60 degrees in the NMR structure compared to 180 degrees in the X-ray. This results in the side chain of N80 occupying and partially blocking access to the active site of MMP-1.

High-resolution solution structure of basic fibroblast growth factor determined by multidimensional heteronuclear magnetic resonance spectroscopy

The high-resolution solution structure of recombinant human basic fibroblast growth factor (FGF-2), a protein of 17.2 kDa that exhibits a variety of functions related to cell growth and differentiation, has been determined using three-dimensional heteronuclear NMR spectroscopy. A total of 30 structures were calculated by means of hybrid distance geometry--simulated annealing using a total of 2865 experimental NMR restraints, consisting of 2486 approximate inteproton distance restraints, 50 distance restraints for 25 backbone hydrogen bonds, and 329 torsion angle restraints. The atomic rms distribution about the mean coordinate positions for the 30 structures for residues 29-152 is 0.43 +/- 0.03 A for the backbone atoms, 0.83 +/- 0.05 A for all atoms, and 0.51 +/- 0.04 A for all atoms excluding disordered side chains. The overall structure of FGF-2 consists of 11 extended antiparallel beta-strands arranged in three groups of three or four strands connected by tight turns and loop regions creating a pseudo-3-fold symmetry. Two strands from each group come together to form a beta-sheet barrel of six antiparallel beta-strands. A helix-like structure was observed for residues 131-136, which is part of the heparin binding site (residues 128-138). The discovery of the helix-like region in the primary heparin binding site instead of the beta-strand conformation described in the X-ray structures may have important implications in understanding the nature of heparin--FGF-2 interactions. A total of seven tightly bound water molecules were found in the FGF-2 structure, two of which are located in the heparin binding site. The first 28 N-terminal residues appear to be disordered, which is consistent with previous X-ray structures. A best fit superposition of the NMR structure of FGF-2 with the 1.9 A resolution X-ray structure by Zhu et al. (1991) yields a backbone atomic rms difference of 0.94 A, indicative of a close similarity between the NMR and X-ray structures.

Recent progress in understanding chemical shifts

In the past three or four years computer hardware and software developments have reached the stage where the nuclear magnetic resonance (NMR) spectra of many molecular systems can now be accurately evaluated. Detailed analysis of chemical shifts may soon become a routine part of solid (and liquid) state NMR structure prediction in chemistry and biology, and this Article covers the development of the topic from its earliest beginnings.

Exhaustive enumeration of protein conformations using experimental restraints

We present an efficient new algorithm that enumerates all possible conformations of a protein that satisfy a given set of distance restraints. Rapid growth of all possible self-avoiding conformations on the diamond lattice provides construction of alpha-carbon representations of a protein fold. We investigated the dependence of the number of conformations on pairwise distance restraints for the proteins crambin, pancreatic trypsin inhibitor, and ubiquitin. Knowledge of between one and two contacts per monomer is shown to be sufficient to restrict the number of candidate structures to approximately 1,000 conformations. Pairwise RMS deviations of atomic position comparisons between pairs of these 1,000 structures revealed that these conformations can be grouped into about 25 families of structures. These results suggest a new approach to assessing alternative protein folds given a very limited number of distance restraints. Such restraints are available from several experimental techniques such as NMR, NOESY, energy transfer fluorescence spectroscopy, and crosslinking experiments. This work focuses on exhaustive enumeration of protein structures with emphasis on the possible use of NOESY-determined distance restraints.

Determination of the three-dimensional structure of margatoxin by 1H, 13C, 15N triple-resonance nuclear magnetic resonance spectroscopy

The solution structure of the 39-residue peptide margatoxin, a scorpion toxin that selectively blocks the voltage-gated potassium-channel Kv1.3, has been determined by NMR spectroscopy. The toxin was isotopically labeled with 13C and 15N and studied using two-dimensional homonuclear and three- and four-dimensional heteronuclear NMR spectroscopy. The final structure was determined using 501 constraints, comprising 422 NOE constraints, 60 dihedral angle constraints, 9 disulfide constraints, and 10 hydrogen bond constraints. Structures were initially determined with the program PEGASUS and subsequently refined with X-PLOR. The average rms deviation from a calculated average structure for the backbone atoms of residues 3-38 is 0.40 A. A helix is present from residues 11 to 20 and includes two proline residues at positions 15 and 16. A loop at residues 21-24 leads into a two-strand antiparallel sheet from residues 25 to 38 with a turn at residues 30-33. Residues 3-6 run adjacent to the 33-38 strand but do not form a canonical beta-strand. The two additional residues of margatoxin, relative to the related toxins charybdotoxin and iberiotoxin, insert in a manner that extends the beta-sheet by one residue. Otherwise, the global structure is very similar to that of these two other toxins. The longer sheet may have implications for channel selectivity.

Differential deuterium isotope shifts and one-bond 1H-13C scalar couplings in the conformational analysis of protein glycine residues

The one-bond deuterium isotope shift effect for glycine C alpha resonances exhibits a conformational dependence comparable to that of the corresponding 1JHC scalar coupling in both magnitude (approximately 11 Hz at 14.1 T) and dihedral angle dependence. The similarity in the conformational dependence of the 1JHC and deuterium isotope shift values suggests a common physical basis. Given the known distribution of (phi, psi) main-chain dihedral angles for glycine residues, the deuterium isotope shifts and the 1JHC scalar couplings can determine conformations in the left- and right-handed helical-to-bridge regions of the (phi, psi) plane to an accuracy of approximately 13 degrees. In the absence of stereochemical assignments, the differential deuterium isotope shifts and the 1JHC scalar couplings can be combined with limited independent structural information (e.g., the sign of phi) to determine the chirality of the deuterium substitution.

Calculations of one-, two- and three-bond nuclear spin-spin couplings in a model peptide and correlations with experimental data

We present ab initio calculations of the Fermi contact term and experimental correlations of six coupling constants, 3JHNH alpha, 1JC alpha H alpha, 2JC'H alpha, 1JC alpha N, 2JC alpha N and 1JC'N, in a peptide as functions of the backbone dihedral angles, phi and psi. Given estimates of experimental uncertainties, we find semiquantitative experimental correlations for 3JHNH alpha, 1JC alpha N and 2JC alpha N, qualitative correlations for 1JC alpha H alpha and 2JC'H alpha, but no experimental correlations of practical utility for 1JC'N, owing to its complex dependence on at least four dihedral angles. Errors in the estimation of dihedral angles from X-ray crystallographic data for proteins, which result from uncertainties in atom-to-atom distances, place substantial limitations on the quantitative reliability of coupling constant calculations fitted to such data. In the accompanying paper [Edison, A.S. et al., J. Biomol. NMR, 4, 543-551] we apply the results of the coupling constant calculations presented here to the estimation of phi and psi angles in staphylococcal nuclease from experimental coupling constants.

The NMR structure of the inhibited catalytic domain of human stromelysin-1

The three-dimensional structure of the catalytic domain of stromelysin-1 complexed with an N-carboxyl alkyl inhibitor has been determined by NMR methods. The global fold consists of three helices, a five stranded beta-sheet and a methionine located in a turn near the catalytic histidines, classifying stromelysin-1 as a metzincin. Stromelysin-1 is unique in having two independent zinc binding sites: a catalytic site and a structural site. The inhibitor binds in an extended conformation. The S1' subsite is a deep hydrophobic pocket, whereas S2' appears shallow and S3' open.

Coupling constants again: experimental restraints in structure refinement

Utilization of coupling constants as restraints in computational structure refinement is reviewed. In addition, we address the effect of conformational averaging and examine different approaches to apply the restraints when the experimental observable is obviously a result of averaging. Here, two different computational methods are compared. The simulation of a single structure with time-dependent restraints produces results very similar to those obtained with the calculation of numerous copies of the molecule (an ensemble of structures) and ensemble averaging. The advantages and disadvantages of the two methods are illustrated with simulations of cyclosporin A, for which 117 NOEs and 62 homo- and heteronuclear coupling constants have been measured.

The high-resolution, three-dimensional solution structure of human interleukin-4 determined by multidimensional heteronuclear magnetic resonance spectroscopy

The high-resolution three-dimensional solution structure of recombinant human interleukin-4 (IL-4), a protein of approximately 15 kDa which plays a key role in the regulation of B and T lymphocytes, has been determined using three- and four-dimensional heteronuclear NMR spectroscopy. The structure is based on a total of 2973 experimental NMR restraints, comprising 2515 approximate interproton distance restraints, 102 distance restraints for 51 backbone hydrogen bonds, and 356 torsion angle restraints. A total of 30 structures was calculated by means of hybrid distance geometry-simulated annealing, and the atomic rms distribution about the mean coordinate positions for residues 8-129 is 0.44 +/- 0.03 A for the backbone atoms, 0.83 +/- 0.03 A for all atoms, and 0.51 +/- 0.04 A for all atoms excluding disordered side chains. The N- and C-terminal residues (1-7 and 130-133, respectively) appear to be disordered. The structure of IL-4 is dominated by a left-handed four-helix bundle with an unusual topology comprising two overhand connections. The linker elements between the helices are formed by either long loops, small helical turns, or short strands. The latter include a mini anti-parallel beta-sheet. A best fit superposition of the NMR structure of IL-4 with the 2.25 A resolution crystal structure [Wlodawer, A., Pavlovsky, A., & Gutschina, A. (1992) FEBS Lett. 309, 59-64] yields a backbone atomic rms difference of 1.37 A which can be mainly attributed to tighter packing of the helices in the crystal structure. This is indicated by an approximately 20% reduction in the axial separation of three pairs of helices (alpha A-alpha C, alpha A-alpha D, and alpha C-alpha D) in the crystal structure relative to the NMR structure and may reflect the greater flexibility of the molecule in solution which is reduced in the crystal due to intermolecular contacts. Comparison of the NMR structure of IL-4 with the X-ray structures of two other related proteins, granulocyte-macrophage colony stimulating factor [Diedrichs, K., Boone, T., & Karplus, P. A. (1992) Science 254, 1779-1782] and human growth hormone [de Vos, A. M., Ultsch, M., & Kossiakoff, A. A. (1992) Science 255, 306-312], that bind to the same hematopoietic superfamily of cell surface receptors reveals a remarkably similar topological fold, despite the absence of any significant overall sequence identity, and substantial differences in the relative lengths of the helices, the lengths and the nature of the various connecting elements, and the pattern and number of disulfide bridges.(ABSTRACT TRUNCATED AT 400 WORDS)