NMR study of the conformational transition of cytochrome c upon the displacement of Met80 by exogenous ligand: structural and magnetic characterization of azidoferricytochrome c

As the exogenous ligand-cytochrome c complexes were purported to represent models for the unfolding intermediate of cytochrome c, NMR spectroscopy has been utilized to study the azide adduct of horse heart cytochrome c. The structure of azidoferricytochrome c was modeled by restrained energy minimization using paramagenetic pseudocontact shifts as constraints. The bound azide moiety was found to be tilted approximately 15 degrees from the heme normal. The displacement of Met80 by the exogenous azide molecule causes large structural rearrangement in the distal cavity. Furthermore, the conformation transition associated with the swing out of the loop containing Met80 and the shift of the 50s-helix increases the solvent accessibility of the heme group. To elucidate the heme electronic structure of the complex, the paramagnetic 13C shifts were analyzed in terms of a model based on the pi molecular orbitals of the heme under perturbed D(4) symmetry. It turned out that the His-Fe bonding provides the protein constraint that orients the in-plane anisotropy in the complex. The electronic properties are in accordance with the calculated magnetic susceptibility anisotropy and the structural information.     

The use of pseudocontact shifts to refine solution structures of paramagnetic metalloproteins: Met80Ala cyano-cytochrome c as an example

 The availability of NOE constraints and of the relative solution structure of a paramagnetic protein permits the use of pseudocontact shifts as further structural constraints. We have developed a strategy based on: (1) determination of the χ tensor anisotropy parameters from the starting structure; (2) recalculation of a new structure by using NOE and pseudocontact shift constraints simultaneously; (3) redetermination of the χ tensor anisotropy parameters from the new structure, and so on until self-consistency. The system investigated is the cyanide derivative of a variant of the oxidized Saccharomyces cerevisiae iso-1-cytochrome c containing the Met80Ala mutation. The structure has been substantially refined. It is shown that the analysis of the deviation of the experimental pseudocontact shifts from those calculated using the starting structure may be unsound, as may the simple structure refinement based on the pseudocontact shift constraints only. Received: 11 July 1995 / Accepted: 30 October 1995     

Redox-related conformational changes in Rhodobacter capsulatus cytochrome c2

WEFT-NOESY and transfer WEFT-NOESY NMR spectra were used to determine the heme proton assignments for Rhodobacter capsulatus ferricytochrome c2. The Fermi contact and pseudo-contact contributions to the paramagnetic effect of the unpaired electron in the oxidized state were evaluated for the heme and ligand protons. The chemical shift assignments for the 1H and 15N NMR spectra were obtained by a combination of 1H-1H and 1H-15N two-dimensional NMR spectroscopy. The short-range nuclear Overhauser effect (NOE) data are consistent with the view that the secondary structure for the oxidized state of this protein closely approximates that of the reduced form, but with redox-related conformational changes between the two redox states. To understand the decrease in stability of the oxidized state of this cytochrome c2 compared to the reduced form, the structural difference between the two redox states were analyzed by the differences in the NOE intensities, pseudo-contact shifts and the hydrogen-deuterium exchange rates of the amide protons. We find that the major difference between redox states, although subtle, involve heme protein interactions, orientation of the heme ligands, differences in hydrogen bond networks and, possible alterations in the position of some internal water molecules. Thus, it appears that the general destabilization of cytochrome c2, which occurs on oxidation, is consistent with the alteration of hydrogen bonds that result in changes in the internal dynamics of the protein.     

Orientation of the axial ligands and magnetic properties of the hemes in the triheme ferricytochrome PpcA from G. sulfurreducens determined by paramagnetic NMR

The geometry of the axial ligands of the hemes in the triheme cytochrome PpcA from Geobacter sulfurreducens was determined in solution for the ferric form using the unambiguous assignment of the NMR signals of the α-substituents of the hemes. The paramagnetic ¹³C shifts of the hemes can be used to define the heme electronic structure, the geometry of the axial ligands, and the magnetic susceptibility tensor. The latter establishes the magnitude and geometrical dependence of the pseudocontact shifts, which are crucial to warrant reliable structural constraints for a detailed structural characterization of this paramagnetic protein in solution.     

An optimized g-tensor for Rhodobacter capsulatus cytochrome c2 in solution: a structural comparison of the reduced and oxidized states.

The optimized g-tensor parameters for the oxidized form of Rhodobacter capsulatus cytochrome c2 in solution were obtained using a set (50) of backbone amide protons. Dipolar shifts for more than 500 individual protons of R. capsulatus cytochrome c2 have been calculated by using the optimized g-tensor and the X-ray crystallographic coordinates of the reduced form of R. capsulatus cytochrome c2. The calculated results for dipolar shifts are compared with the observed paramagnetic shifts. The calculated and the observed data are in good agreement throughout the entire protein, but there are significant differences between calculated and experimental results localized to the regions in the immediate vicinity of the heme ligand and the region of the front crevice of the protein (residues 44-50, 53-57, and 61-68). The results not only indicate that the overall solution structures are very similar in both the reduced and oxidized states, but that these structures in solution are similar to the crystal structure. However, there are small structural changes near the heme and the rearrangement of certain residues that result in changes in their hydrogen bonding concomitant with the change in the oxidation states; this was also evident in the data for the NH exchange rate measurements for R. capsulatus cytochrome c2.     

Pseudocontact shifts as constraints for energy minimization and molecular dynamics calculations on solution structures of paramagnetic metalloproteins

The pseudocontact shifts of NMR signals, which arise from the magnetic susceptibility anisotropy of paramagnetic molecules, have been used as structural constraints under the form of a pseudopotential in the SANDER module of the AMBER 4.1 molecular dynamics software package. With this procedure, restrained energy minimization (REM) and restrained molecular dynamics (RMD) calculations can be performed on structural models by using pseudocontact shifts. The structure of the cyanide adduct of the Met80Ala mutant of the yeast iso-1-cytochrome c has been used for successfully testing the calculations. For this protein, a family of structures is available, which was obtained by using NOE and pseudocontact shifts as constraints in a distance geometry program. The structures obtained by REM and RMD calculations with the inclusion of pseudocontact shifts are analyzed. Proteins 29:68–76, 1997. © 1997 Wiley-Liss, Inc.     

Chemical shift-based constraints for solution structure determination of paramagnetic low-spin heme proteins with bis-His and His-CN axial ligands: the cases of oxidized cytochrome b 5 and Met80Ala cyano-cytochrome c

The chemical shifts of the methyl protons of protoporphyrin IX, which are readily assigned, are related to the structural features of the axial histidine ligands in heme proteins with bis-His or His-CN axial coordination (Bertini I, Luchinat C, Parigi G, Walker FA (1999) JBIC 4:515-519). In the present paper, a module is developed which transforms the chemical shifts into a pseudo-potential energy that is a function of the dihedral angles defining the orientation of the axial ligand planes. Minimization of this pseudo-potential energy, together with the energetic contributions provided by the other constraints, yields structures consistent with the heme methyl chemical shifts. Oxidized cytochrome b(5) from rat and the cyanide derivative of the M80A mutant of yeast cytochrome c are used for test calculations. In the case of scarcity of NOEs for the axial ligands, owing to the presence of the paramagnetic center, the above structural constraints are shown to be quite precious. The newly refined structures are deposited in the PDB.     

The solution structure of the oxidized bovine microsomal cytochrome b(5) mutant V61H

Using 1488 NOE constraints, 19 stereo-specific assignments, 13 pairs of H-bond constraints, and 140 pseudo-contact shift constraints, a family of 35 structures of bovine microsomal cytochrome b(5) mutant V61H has been obtained through the program PSEUDYANA. The family has been further refined by restrained energy minimization to give a family of final structures. The RMSD values of final structures with respect to the average structure are 0.45+/-0.11 and 0.96+/-0.10A for backbone and heavy atoms, respectively. The final Deltachi(ax) and Deltachi(rh) values are 2.34 x 10(-32) and -0.67 x 10(-32)m(3), respectively. The comparisons between the solution structures of mutant V61H and WT cytochrome b(5), and X-ray structure of the mutant V61H show that the global folding of the molecule in solution is unchanged and the side-chain of His61 deviates from the heme pocket and extends into the solvent like in its crystal structure. However, the helices around the heme pocket undergo outward global displacement while their local conformations are well maintained. Meanwhile, the heme ring shows a little off the heme pocket, which accounts for the lower stability of the mutant. Additionally, the axial ligand rings counterclockwise rotate around His39 N-Fe axis due to the mutation, which is confirmed by variation of the hyperfine shifts of the heme protons of V61H compared to those of WT cytochrome b(5).     

Redox-dependent structure change and hyperfine nuclear magnetic resonance shifts in cytochrome c

Proton nuclear magnetic resonance assignments for reduced and oxidized equine cytochrome c show that many individual protons exhibit different chemical shifts in the two protein forms, reflecting diamagnetic shift effects due to structure change, and in addition contact and pseudocontact shifts that occur only in the paramagnetic oxidized form. To evaluate the chemical shift differences (delta delta) for structure change, we removed the pseudocontact shift contribution by a calculation based on knowledge of the electron spin g tensor. The g-tensor parameters were determined from the delta delta values of a large set (64) of C alpha H protons at well-defined spatial positions in the oxidized horse protein. The g-tensor calculation, when repeated using only 12 available C alpha H proton resonances for cytochrome c from tuna, proved to be remarkably stable. The largest principal value of the g tensor (gz) falls precisely along the ligand bond between the heme iron and methionine-80 sulfur, while gx and gy closely match the natural heme axes defined by the pyrrole nitrogens. The derived g tensor was then used together with spatial coordinates for the oxidized form to calculate the pseudocontact shift contribution (delta pc) to proton resonances at 400 identifiable sites throughout the protein, so that the redox-dependent chemical shift discrepancy, delta delta-delta pc, could be evaluated. Large residual changes in chemical shift define the Fermi contact shifts, which are found as expected to be limited to the immediate covalent structure of the heme and its ligands and to be asymmetrically distributed over the heme. Smaller chemical shift discrepancies point to a concerted change, involving residues 39-43 and 50-60 (bottom of the protein), and to other changes in the immediate vicinity of the heme ligands. Also, the three internal water molecules are implicated in redox sensitivity. The residues found to change are in good but not perfect agreement with prior X-ray diffraction observations of subangstrom redox-related displacements in the tuna protein. The chemical shift discrepancies observed appear in the main to reflect structure-dependent diamagnetic shifts rather than hyperfine effects due to displacements in the pseudocontact shift field. Although 51 protons in 29 different residues exhibit significant chemical shift changes, the general impression is one of small structural adjustments to redox-dependent strain rather than sizeable structural displacements or rearrangements.     

Structure and stability effects of the mutation of glycine 34 to serine in Rhodobacter capsulatus cytochrome c(2)

Gly 34 and the adjacent Pro 35 of Rhodobacter capsulatus cytochrome c(2) (or Gly 29 and Pro 30 in vertebrate cytochrome c) are highly conserved side chains among the class I c-type cytochromes. The mutation of Gly 34 to Ser in Rb. capsulatus cytochrome c(2) has been characterized in terms of physicochemical properties and NMR in both redox states. A comparison of the wild-type cytochrome c(2), the G34S mutation, and the P35A mutation is presented in the context of differences in chemical shifts, the differences in NOE patterns, and structural changes resulting from oxidation of the reduced cytochrome. G34S is substantially destabilized relative to wild-type (2.2 kcal/mol in the oxidized state) but similarly destabilized relative to P35A. Nevertheless, differences in terms of the impact of the mutations on specific structural regions are found when comparing G34S and P35A. Although available data indicates that the overall secondary structure of G34S and wild-type cytochrome c(2) are similar, a number of both perturbations of hydrogen bond networks and interactions with internal waters are found. Thus, the impact of the mutation at position 35 is propagated throughout the cytochrome but with alterations at defined sites within the molecule. Interestingly, we find that the substitution of serine at position 34 results in a perturbation of the heme beta meso and the methyl-5 protons. This suggests that the hydroxyl and beta carbon are positioned away from the solvent and toward the heme. This has the consequence of preferentially stabilizing the oxidized state in G34S, thus, altering hydrogen bond networks which involve the heme propionate, internal waters, and key amino acid side chains. The results presented provide important new insights into the stability and solution structure of the cytochrome c(2).     

Nuclear magnetic resonance signal chemical shifts and molecular simulations: a multidisciplinary approach to modeling copper protein structures

Nuclear magnetic resonance (NMR) chemical shifts are experimental observables that are available during the first stage of the protein structure determination process. Recently, some methodologies for building structural models of proteins using only these experimental data have been implemented. To assess the potential of these methods for modeling metalloproteins (generally considered a challenging benchmark), we determined the structures of the yeast copper chaperone Atx1 and the CuA domain of Thermus thermophilus cytochrome c oxidase starting from the available chemical shift data. The metal centers were modeled using molecular dynamics simulations with molecular mechanics potentials. The results obtained are evaluated and discussed.     

Paramagnetism-Based Restraints for Xplor-NIH

Modules that use paramagnetism-based NMR restraints have been developed and integrated in the well known program for solution structure determination Xplor-NIH; the complete set of such modules is called PARArestraints for Xplor-NIH. Paramagnetism-based restraints are paramagnetic relaxation enhancements, pseudocontact shifts, residual dipolar couplings due to metal and overall magnetic anisotropy, and cross correlation between Curie relaxation and nuclear-nuclear dipolar relaxation. The complete program has been tested by back-calculating NOEs and paramagnetism-based restraints from the X-ray structure of cytochrome c (553) from B. pasteurii. Furthermore, the same experimental restraints previously used to determine the solution structure of cytochrome c (553) itself, of cytochrome b (5), and of calbindin D(9k) with the program PARAMAGNETIC DYANA, have been used for structure calculations by using PARArestraints for Xplor-NIH. The agreement between the two programs is quite satisfactory and validates both protocols.     

Determination of solution structures of paramagnetic proteins by NMR

Standard procedures for using nuclear Overhauser enhancements (NOE) between protons to generate structures for diamagnetic proteins in solution from NMR data may be supplemented by using dipolar shifts if the protein is paramagnetic. This is advantageous since the electron-nuclear dipolar coupling provides relatively long-range geometric information with respect to the paramagnetic centre which complements the short-range distance constraints from NOEs. Several different strategies have been developed to date, but none of these attempts to combine data from NOEs and dipolar shifts in the initial stages of structure calculation or to determine three dimensional protein structures together with their magnetic properties. This work shows that the magnetic and atomic structures are highly correlated and that it is important to have additional constraints both to provide starting parameters for the magnetic properties and to improve the definition of the best fit. Useful parameters can be obtained for haem proteins from Fermi contact shifts; this approach is compared with a new method based on the analysis of dipolar shifts in haem methyl groups with respect to data from horse and tuna ferricytochromes c. The methods developed for using data from NOEs and dipolar shifts have been incorporated in a new computer program, PARADYANA, which is demonstrated in application to a model data set for the sequence of the haem octapeptide known as microperoxidase-8. Received: 13 October 1997 / Accepted: 19 December 1997     

15N-1H Residual dipolar coupling analysis of native and alkaline-K79A Saccharomyces cerevisiae cytochrome c

Residual dipolar couplings (RDCs) and pseudocontact shifts are experimentally accessible properties in nuclear magnetic resonance that are related to structural parameters and to the magnetic susceptibility anisotropy. We have determined RDCs due to field-induced orientation of oxidized-K79A and reduced cytochrome c at pH 7.0 and oxidized-K79A cytochrome c at pH 11.1 through measurements of amide (15)N-(1)H (1)J couplings at 800 and 500 MHz. The pH 7.0 RDCs for Fe(III)- and Fe(II)-cytochrome c together with available nuclear Overhauser effects were used to recalculate solution structures that were consistent with both sets of constraints. Molecular magnetic susceptibility anisotropy values were calculated for both redox states of the protein. By subtracting the residual dipolar couplings (RDCs) of the reduced form from those of the oxidized form measured at the same magnetic field (800 MHz), we found the RDC contribution of the paramagnetic metal ion in the oxidized protein. The magnetic susceptibility anisotropy, which was calculated from the structure, was found to be the same as that of the paramagnetic metal ion obtained independently from pseudocontact shifts, thereby indicating that the elements of secondary structure either are rigid or display the same mobility in both oxidation states. The residual dipolar coupling values of the alkaline-K79A form are small with respect to those of oxidized native cytochrome, whereas the pseudocontact shifts are essentially of the same magnitude, indicating local mobility. Importantly, this is the first time that mobility has been found through comparison of RDCs with pseudocontact shifts.     

De novo determination of protein structure by NMR using orientational and long-range order restraints

Orientational and novel long-range order restraints available from paramagnetic systems have been used to determine the backbone solution structure of the cytochrome c' protein to atomic resolution in the complete absence of restraints derived from the nuclear Overhauser effect. By exploiting the complementary geometric dependence of paramagnetic pseudocontact shifts and the recently proposed Curie-dipolar cross correlated relaxation effect, in combination with orientational constraints derived from residual dipolar coupling, autorelaxation rate ratios and secondary structure constraints, it is possible to define uniquely the fold and refine the tertiary structure of the protein (0.73 A backbone rmsd for 82/129 amino acid residues) starting from random atomic Cartesian coordinates. The structure calculation protocol, developed using specific models to describe the novel constraint interactions, is robust, requiring no precise a priori estimation of the various interaction strengths, and provides unambiguous convergence based only on the value of the target function. Tensor eigenvalues and their component orientations are allowed to float freely, and are thus simultaneously determined, and found to converge, during the structure calculation.     

De novo protein structure generation from incomplete chemical shift assignments

NMR chemical shifts provide important local structural information for proteins. Consistent structure generation from NMR chemical shift data has recently become feasible for proteins with sizes of up to 130 residues, and such structures are of a quality comparable to those obtained with the standard NMR protocol. This study investigates the influence of the completeness of chemical shift assignments on structures generated from chemical shifts. The Chemical-Shift-Rosetta (CS-Rosetta) protocol was used for de novo protein structure generation with various degrees of completeness of the chemical shift assignment, simulated by omission of entries in the experimental chemical shift data previously used for the initial demonstration of the CS-Rosetta approach. In addition, a new CS-Rosetta protocol is described that improves robustness of the method for proteins with missing or erroneous NMR chemical shift input data. This strategy, which uses traditional Rosetta for pre-filtering of the fragment selection process, is demonstrated for two paramagnetic proteins and also for two proteins with solid-state NMR chemical shift assignments. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structure determination of a Galectin-3-carbohydrate complex using paramagnetism-based NMR constraints

The determination of the location and conformation of a natural ligand bound to a protein receptor is often a first step in the rational design of molecules that can modulate receptor function. NMR observables, including NOEs, often provide the basis for these determinations. However, when ligands are carbohydrates, interactions mediated by extensive hydrogen-bonding networks often reduce or eliminate NOEs between ligand and protein protons. In these cases, it is useful to look to other distance- and orientation-dependent observables that can constrain the geometry of ligand-protein complexes. Here we illustrate the use of paramagnetism-based NMR constraints, including pseudo-contact shifts (PCS) and field-induced residual dipolar couplings (RDCs). When a paramagnetic center can be attached to the protein, field-induced RDCs and PCS reflect only bound-state properties of the ligand, even when averages over small fractions of bound states and large fractions of free states are observed. The effects can also be observed over a long range, making it possible to attach a paramagnetic center to a remote part of the protein. The system studied here is a Galectin-3-lactose complex. A lanthanide-binding peptide showing minimal flexibility with respect to the protein was integrated into the C terminus of an expression construct for the Galectin-3-carbohydrate-binding domain. Dysprosium ion, which has a large magnetic susceptibility anisotropy, was complexed to the peptide, making it possible to observe both PCSs and field-induced RDCs for the protein and the ligand. The structure determined from these constraints shows agreement with a crystal structure of a Galectin-3-N-acetyllactosamine complex.     

Strategies for the uses of lanthanide NMR shift probes in the determination of protein structure in solutio. Application to the EF calcium binding site of carp parvalbumin.

The homologous sequences observed for many calcium binding proteins such as parvalbumin, troponin C, the myosin light chains, and calmodulin has lead to the hypothesis that these proteins have homologous structures at the level of their calcium binding sites. This paper discusses the development of a nuclear magnetic resonance (NMR) technique which will enable us to test this structural hypothesis in solution. The technique involves the substitution of a paramagnetic lanthanide ion for the calcium ion which results in lanthanide induced shifts and broadening in the 1H NMR spectrum of the protein. These shifts are sensitive monitors of the precise geometrical orientation of each proton nucleus relative to the metal. The values of several parameters in the equation relating the NMR shifts to the structure are however known as priori. We have attempted to determine these parameters, the orientation and principal elements of the magnetic susceptibility tensor of the protein bound metal, by studying the lanthanide induced shifts for the protein parvalbumin whose structure has been determined by x-ray crystallographic techniques. The interaction of the lanthanide ytterbium with parvalbumin results in high resolution NMR spectra exhibiting a series of resonances with shifts spread over the range 32 to -19 ppm. The orientation and principal elements of the ytterbium magnetic susceptibility tensor have been determined using three assigned NMR resonances, the His-26 C2 and C4 protons and the amino terminal acetyl protons, and seven methyl groups; all with known geometry relative to the EF calcium binding site. The elucidation of these parameters has allowed us to compare the observed spectrum of the nuclei surrounding the EF calcium binding site of parvalbumin with that calculated from the x-ray structure. A significant number of the calculated shifts are larger than any of the observed shifts. We feel that a refinement of the x-ray based proton coordinates will be possible utilizing the geometric information contained in the lanthanide shifted NMR spectrum.