NMR Conformers Research Articles

TROSY experiment for refinement of backbone psi and phi by simultaneous measurements of cross-correlated relaxation rates and 3,4J(H alpha HN) coupling constants.

The TROSY principle has been introduced into a HNCA experiment, which is designed for measurements of the intraresidual and sequential H(alpha)-C(alpha)/H(N)-N dipole/dipole and H(alpha)-C(alpha)/N dipole/CSA cross-correlated relaxation rates. In addition, the new experiment provides values of the (3,4)J(H alpha HN) coupling constants measured in an E.COSY manner. The conformational restraints for the psi and phi angles are obtained through the use of the cross-correlated relaxation rates together with the Karplus-type dependencies of the coupling constants. Improved signal-to-noise is achieved through preservation of all coherence transfer pathways and application of the TROSY principle. The application of the [(15)N,(13)C]-DQ/ZQ-[(15)N,(1)H]-TROSY-E.COSY experiment to the 16 kDa apo-form of the E. coli Heme Chaperon protein CcmE is described. Overall good agreement is achieved between psi and phi angles measured with the new experiment and the average values determined from an ensemble of 20 NMR conformers.

Solution structure of the catalytic domain of γδ resolvase. Implications for the mechanism of catalysis

The site-specific DNA recombinase, γδ resolvase, from Escherichia coli catalyzes recombination of res site-containing plasmid DNA to two catenated circular DNA products. The catalytic domain (residues 1-105), lacking a C-terminal dimerization interface, has been constructed and the NMR solution structure of the monomer determined. The RMSD of the NMR conformers for residues 2–92 excluding residues 37–45 and 64–73 is 0.41 Å for backbone atoms and 0.88 Å for all heavy atoms. The NMR solution structure of the monomeric catalytic domain (residues 1–105) was found to be formed by a four-stranded parallel β-sheet surrounded by three helices. The catalytic domain (residues 1–105), deficient in the C-terminal dimerization domain, was monomeric at high salt concentration, but displayed unexpected dimerization at lower ionic strength. The unique solution dimerization interface at low ionic strength was mapped by NMR. With respect to previous crystal structures of the dimeric catalytic domain (residues 1-140), differences in the average conformation of active-site residues were found at loop 1 containing the catalytic S10 nucleophile, the β1 strand containing R8, and at loop 3 containing D67, R68 and R71, which are required for catalysis. The active-site loops display high-frequency and conformational backbone dynamics and are less well defined than the secondary structures. In the solution structure, the D67 side-chain is proximal to the S10 side-chain making the D67 carboxylate group a candidate for activation of S10 through general base catalysis. Four conserved Arg residues can function in the activation of the phosphodiester for nucleophilic attack by the S10 hydroxyl group. A mechanism for covalent catalysis by this class of recombinases is proposed that may be related to dimer interface dissociation.

Fluctuations in ion pairs and their stabilities in proteins.

This report investigates the effect of systemic protein conformational flexibility on the contribution of ion pairs to protein stability. Toward this goal, we use all NMR conformer ensembles in the Protein Data Bank (1) that contain at least 40 conformers, (2) whose functional form is monomeric, (3) that are nonredundant, and (4) that are large enough. We find 11 proteins adhering to these criteria. Within these proteins, we identify 22 ion pairs that are close enough to be classified as salt bridges. These are identified in the high-resolution crystal structures of the respective proteins or in the minimized average structures (if the crystal structures are unavailable) or, if both are unavailable, in the "most representative" conformer of each of the ensembles. We next calculate the electrostatic contribution of each such ion pair in each of the conformers in the ensembles. This results in a comprehensive study of 1,201 ion pairs, which allows us to look for consistent trends in their electrostatic contributions to protein stability in large sets of conformers. We find that the contributions of ion pairs vary considerably among the conformers of each protein. The vast majority of the ion pairs interconvert between being stabilizing and destabilizing to the structure at least once in the ensembles. These fluctuations reflect the variabilities in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other, and with respect to other charged groups in the remainder of the protein. The higher crystallographic B-factors for the respective side-chains are consistent with these fluctuations. The major conclusion from this study is that salt bridges observed in crystal structure may break, and new salt bridges may be formed. Hence, the overall stabilizing (or, destabilizing) contribution of an ion pair is conformer population dependent.

Backbone dynamics and refined solution structure of the N-terminal domain of DNA polymerase β. Correlation with DNA binding and dRP lyase activity

Mammalian DNA polymerase β functions in the base excision DNA repair pathway filling in short patches (1–5 nt) in damaged DNA and removing deoxyribose 5′-phosphate from the 5′-side of damaged DNA. The backbone dynamics and the refined solution structure of the N-terminal domain of β-Pol have been characterized in order to establish the potential contribution(s) of backbone motion to the DNA binding and deoxyribose 5′-phosphate lyase function of this domain. The N-terminal domain is formed from four helices packed as two antiparallel pairs with a 60° crossing between the pairs. The RMSD of the NMR conformers (residues 13–80) is 0.37 Å for the backbone heavy atoms and 0.78 Å for all heavy atoms. NMR characterization of the binding site(s) for a ssDNA-5mer, ssDNA-8mer, ssDNA-9mer, and dsDNA-12mer shows a consensus surface for the binding of these various DNA oligomers, that surrounds and includes the deoxyribose 5′-phosphate lyase active site region. Connection segments between helices 1 and 2 and between helices 3 and 4 each contribute to DNA binding. Helix-3-turn-helix-4 forms a helix-hairpin-helix motif. The highly conserved hairpin sequence (LPGVG) displays a significant degree of picosecond time-scale motion within the backbone, that is possibly important for DNA binding at the phosphodiester backbone. An Ω-loop connecting helices 1 and 2 and helix-2 itself display significant exchange contributions (Rex) at the backbone amides due to apparent conformational type motion on a millisecond time-scale. This motion is likely important in allowing the Ω-loop and helix-2 to shift toward, and productively interact with, gapped DNA. The deoxyribose 5′-phosphate lyase catalytic residues that include K72 which forms the Schiff’s base, Y39 which is postulated to promote proton transfer to the aldehyde, and K35 which assists in phosphate elimination, show highly restricted backbone motion. H34, which apparently participates in detection of the abasic site hole and assists in the opening of the hemiacetal, shows conformational exchange.

Fluctuations between stabilizing and destabilizing electrostatic contributions of ion pairs in conformers of the c-Myc-Max leucine zipper

In solution proteins often exhibit backbone and side-chain flexibility. Yet electrostatic interactions in proteins are sensitive to motions. Hence, here we study the contribution of ion pairs toward protein stability in a range of conformers which sample the conformational space in solution. Specifically, we focus on the electrostatic contributions of ion pairs to the stability of each of the conformers in the NMR ensemble of the c-Myc-Max leucine zipper and to their average energy minimized structure. We compute the electrostatic contributions of inter- and intra-helical ion pairs and of an ion pair network. We find that the electrostatic contributions vary considerably among the 40 NMR conformers. Each ion pair, and the network, fluctuates between being stabilizing and being destabilizing. This fluctation reflects the variability in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other and with respect to other charged groups in the rest of the protein. Ion pair interactions in the c-Myc-Max leucine zipper in solution depend on the protein conformer which is analyzed. Hence, the overall stabilizing (or destabilizing) contribution of an ion pair is conformer population-dependent. This study indicates that free energy calculations performed using the continuum electrostatics methodology are sensitive to protein conformational details.

Solution structure of glutathione bound to human glutathione transferase P1-1: comparison of NMR measurements with the crystal structure.

The conformation of the bound glutathione (GSH) in the active site of the human glutathione transferase P1-1 (EC 2.5.1.18) has been studied by transferred NOE measurements and compared with those obtained by X-ray diffraction data. Two-dimensional TRNOESY and TRROESY experiments have been performed under fast-exchange conditions. The family of GSH conformers, compatible with TRNOE distance constraints, shows a backbone structure very similar to the crystal model. Interesting differences have been found in the side chain regions. After restrained energy minimization of a representative NMR conformer in the active site, the sulfur atom is not found in hydrogen-bonding distance of the hydroxyl group of Tyr 7. This situation is similar to the one observed in an "atypical" crystal complex grown at low pH and low temperature. The NMR conformers display also a poorly defined structure of the glutamyl moiety, and the presence of an unexpected intermolecular NOE could indicate a different interaction of this substrate portion with the G-site. The NMR data seem to provide a snapshot of GSH in a precomplex where the GSH glutamyl end is bound in a different fashion. The existence of this precomplex is supported by pre-steady-state kinetic experiments [Caccuri, A. M., Lo Bello, M., Nuccetelli, M., Nicotra, M., Rossi, P., Antonini, G., Federici, G., and Ricci, G. (1998) Biochemistry 37, 3028-3034] and preliminary time-resolved fluorescence data.

High-resolution solution NMR structure of the Z domain of staphylococcal protein A

Staphylococcal protein A (SpA) is a cell-wall-bound pathogenicity factor from the bacterium Staphylococcus aureus. Because of their small size and immunoglobulin (IgG)-binding activities, domains of protein A are targets for protein engineering efforts and for the development of computational approaches for de novo protein folding. The NMR solution structure of an engineered IgG-binding domain of SpA, the Z domain (an analog of the B domain of SpA), has been determined by simulated annealing with restrained molecular dynamics on the basis of 671 conformational constraints. The Z domain contains three well-defined α-helices corresponding to polypeptide segments Lys7 to Leu17 (helix 1), Glu24 to Asp36 (helix 2), and Ser41 to Ala54 (helix 3). A family of ten conformers representing the solution structure of the Z domain was computed by simulated annealing of restrained molecular dynamics using the program CONGEN. The average of the root-mean-square deviations (r.m.s.d.) of the individual NMR conformers, relative to the mean coordinates, for the backbone atoms N, C α and C′ of residues Phe5 through Ala56 is 0.69 Å; the corresponding backbone r.m.s.d. for the three-helical core is 0.44 Å. Helices 1, 2 and 3 are antiparallel in orientation (Ω 12 = −170(±4)°, Ω 13 = +16(±3)°, Ω 23 = +173(±7)°). A comparison of backbone amide hydrogen/deuterium exchange rates in free and IgG-bound Z domains demonstrates that the amide protons of helices 1, 2 and 3 are protected from rapid exchange in both states, indicating that all three helices are also intact in the IgG-bound state. These solution NMR results differ from the previously determined X-ray structure of the similar SpA B domain in complex with the F c fragment of a human IgG antibody, where helix 3 is not observed in the electron density map and from the solution NMR structure of the B domain, where helix 3 is observed but helix 1 is tilted by ∼30° with respect to helices 2 and 3. Hydrogen-bonded N-cap and C-cap formation is observed for all three helices of the Z domain; these capping interactions appear to be highly conserved in the five homologous domains of SpA.

Determination of the solution structure of the SH3 domain of human p56 Lck tyrosine kinase.

The solution structure of the SH3 domain of human p56 Lck tyrosine kinase (Lck-SH3) has been determined by multidimensional heteronuclear NMR spectroscopy. The structure was calculated from a total of 935 experimental restraints comprising 785 distance restraints derived from 1017 assigned NOE cross peaks and 150 dihedral angel restraints derived from 160 vicinal coupling constants. A novel combination of the constant-time 1H-13C NMR correlation experiment recorded with various delays of the constant-time refocusing delays and a fractionally 13C-labelled sample was exploited for the stereospecific assignment of prochiral methyl groups. Additionally, 28 restraints of 14 identified hydrogen bonds were included. A family of 25 conformers was selected to characterize the solution structure. The average root-mean-square deviations of the backbone atoms (N, C alpha, C', O) among the 25 conformers is 0.42 A for residues 7 to 63. The N- and C-terminal residues, 1 to 6 and 64 to 81, are disordered, while the well-converged residues 7 to 63 correspond to the conserved sequences of other SH3 domains. The topology of the SH3 structure comprises five anti-parallel beta-strands arranged to form two perpendicular beta-sheets, which are concave and twisted in the middle part. The overall secondary structure and the backbone conformation of the core beta-strands are almost identical to the X-ray structure of the fragment containing the SH2-SH3 domains of p56 Lck [Eck et al. (1994) Nature, 368, 764-769]. The X-ray structure of the SH3 domain in the tandem SH2-SH3 fragment is spatially included within the ensemble of the 25 NMR conformers, except for the segment of residues 14 to 18, which makes intermolecular contacts with an adjacent SH2 molecule and the phosphopeptide ligand in the crystal lattice. Local structural differences from other known SH3 domains are also observed, the most prominent of which is the absence in Lck-SH3 of the two additional short beta-strands in the regions Ser15 to Glu17 and Gly25 to Glu27 flanking the so-called "RT-Src' loop. This loop (residues Glu17 to Leu24), together with the "n-Src' loop (residues Gln37 to Ser46) confines the ligand interaction site which is formed by a shallow patch of hydrophobic amino acids (His14, Tyr16, Trp41, Phe54 and Phe59). Both loops are flexible and belong to the most mobile regions of the protein, which is assessed by the heteronuclear 15N, 1H-NOE values characterizing the degree of internal backbone motions. The aromatic residues of the ligand binding site are arranged such that they form three pockets for interactions with the polyproline ligand.

Solution structure of a naturally-occurring zinc-peptide complex demonstrates that the N-terminal zinc-binding module of the Lasp-1 LIM domain is an independent folding unit.

The three-dimensional solution structure of the 1:1 complex between the synthetic peptide ZF-1 and zinc was determined by 1H NMR spectroscopy. The peptide, initially isolated from pig intestines, is identical in sequence to the 30 N-terminal amino acid residues of the human protein Lasp-1 belonging to the LIM domain protein family. The final set of 20 energy-refined NMR conformers has an average rmsd relative to the mean structure of 0.55 A for the backbone atoms of residues 3-30. Calculations without zinc atom constraints unambiguously identified Cys 5, Cys 8, His 26, and Cys 29 as the zinc-coordinating residues. LIM domains consist of two sequential zinc-binding modules and the NMR structure of the ZF-1-zinc complex is the first example of a structure of an isolated module. Comparison with the known structures of the N-terminal zinc-binding modules of both the second LIM domain of chicken CRP and rat CRIP with which ZF-1 shares 50% and 43% sequence identity, respectively, supports the notion that the zinc-binding modules of the LIM domain have a conserved structural motif and identifies local regions of structural diversity. The similarities include conserved zinc-coordinating residues, a rubredoxin knuckle involving Cys 5 and Cys 8, and the coordination of the zinc ion by histidine N delta in contrast to the more usual coordination by N epsilon observed for other zinc-finger domains. The present structure determination of the ZF-1-zinc complex establishes the N-terminal half of a LIM domain as an independent folding unit. The structural similarities of N- and C-terminal zinc-binding modules of the LIM domains, despite limited sequence identity, lead to the proposal of a single zinc-binding motif in LIM domains. The coordinates are available from the Brookhaven protein data bank, entry 1ZFO.

Simultaneous Refinement of the Structure of BPTI Against NMR Data Measured in Solution and X-ray Diffraction Data Measured in Single Crystals

The structure of the bovine pancreatic trypsin inhibitor (BPTI) has been determined to high resolution by both NMR spectroscopy in solution and X-ray diffraction in crystals. The root-mean-square difference calculated between the two structures for the polypeptide backbone is 0·9 Å. Several amino acid side-chains, of which all but one are charged or polar, have different conformations. We find that by refining one structure simultaneously against both the NMR and crystallographic data sets, it can accommodate both. Different starting configurations were used, including the X-ray structure 5pti, an NMR conformer, and the X-ray structure in the full unit cell with extra solvent placed in the bulk solvent region. The X-ray structures quickly converged to accommodate the NMR data in addition to the crystallographic data. Starting from an NMR conformer, however, the convergence to accommodate the more abundant X-ray data in addition to the NMR data is much slower.