Analysis of backbone dynamics in cytochrome b5 using 15N-NMR relaxation measurements

GRPY: An Accurate Bead Method for Calculation of Hydrodynamic Properties of Rigid Biomacromolecules

Two main problems that arise in the context of hydrodynamic bead modeling are an inaccurate treatment of bead overlaps and the necessity of using volume corrections when calculating intrinsic viscosity. We present a formalism based on the generalized Rotne-Prager-Yamakawa approximation that successfully addresses both of these issues. The generalized Rotne-Prager-Yamakawa method is shown to be highly effective for the calculation of transport properties of rigid biomolecules represented as assemblies of spherical beads of different sizes, both overlapping and nonoverlapping. We test the method on simple molecular shapes as well as real protein structures and compare its performance with other computational approaches.

Mitochondrial and Microsomal Ferric b5 Cytochromes Exhibit Divergent Conformational Plasticity in the Context of a Common Fold

Native-state hydrogen-deuterium exchange (HDX) monitored by NMR spectroscopy has been used to compare conformational plasticity in ferric rat liver outer mitochondrial membrane cytochrome b5 (rOM b5) and ferric bovine liver microsomal cytochrome b5 (bMc b5). Analysis of the data indicated that rOM b5 is the less conformationally flexible protein on the time scale probed by the HDX experiments. The data also suggest a likely contributor to the much higher kinetic barrier for the release of hemin from OM b5s in comparison to Mc b5s, a characteristic that may be to a large extent the source of their divergent functional properties. Specifically, the data indicate that conformational mobility within helices alpha4 and alpha5, which flank the loop harboring axial ligand His63, is considerably more restricted in rOM b5 than in bMc b5. The lower conformational flexibility of alpha4 and alpha5 in rOM b5 can reasonably be attributed to more extensive hydrophobic packing in that region of the protein, arising from two conserved side chain packing motifs in OM cytochrome b5s [Altuve, A., Wang, L., Benson, D. R., and Rivera, M. (2004) Biochem. Biophys. Res. Commun. 314, 602-609].

Myoglobin and Cytochromeb5: A Nuclear Magnetic Resonance Study of a Highly Dynamic Protein Complex†

The transient complex of bovine myoglobin and cytochrome b(5) has been investigated using a combination of NMR chemical shift mapping, (15)N relaxation data, and protein docking simulations. Chemical shift perturbations observed for cytochrome b(5) amide resonances upon complex formation with either metmyoglobin (Fe(III)) or carbon monoxide-bound myoglobin (Fe(II)) are more than 10-fold smaller than in other transient redox protein complexes. From (15)N relaxation experiments, an increase in the overall correlation time of cytochrome b(5) in the presence of myoglobin is observed, confirming that complex formation is occurring. The chemical shift perturbations of proton and nitrogen amide nuclei as well as heme protons of cytochrome b(5) titrate with increasing myoglobin concentrations, also demonstrating the formation of a weak complex with a K(a) in the inverse millimolar range. The perturbed residues map over a wide surface area of cytochrome b(5), with patches of residues located around the exposed heme 6-propionate as well as at the back of the protein. The nature of the affected residues is mostly negatively charged contrary to perturbed residues in other transient complexes, which are mainly hydrophobic or polar. Protein docking simulations using the NMR data as constraints show several docking geometries both close to and far away from the exposed heme propionates of myoglobin. Overall, the data support the emerging view that this complex consists of a dynamic ensemble of orientations in which each protein constantly diffuses over the surface of the other. The characteristic NMR features may serve as a structural tool for the identification of such dynamic complexes.

Redox-coupled conformational alternations in cytochrome c(3) from D. vulgaris Miyazaki F on the basis of its reduced solution structure

Heteronuclear NMR spectroscopy was performed to determine the solution structure of (15)N-labeled ferrocytochrome c(3) from Desulfovibrio vulgaris Miyazaki F (DvMF). Although the folding of the reduced cytochrome c(3) in solution was similar to that of the oxidized one in the crystal structure, the region involving hemes 1 and 2 was different. The redox-coupled conformational change is consistent with the reported solution structure of D. vulgaris Hildenborough ferrocytochrome c(3), but is different from those of other cytochromes c(3). The former is homologous with DvMF cytochrome c(3) in amino acid sequence. Small displacements of hemes 1 and 2 relative to hemes 3 and 4 were observed. This observation is consistent with the unusual behavior of the 2(1)CH(3) signal of heme 3 reported previously. As shown by the (15)N relaxation parameters of the backbone, a region between hemes 1 and 2 has more flexibility than the other regions. The results of this work strongly suggest that the cooperative reduction of hemes 1 and 2 is based on the conformational changes of the C-13 propionate of heme 1 and the aromatic ring of Tyr43, and the interaction between His34 and His 35 through covalent and coordination bonds.

Structural and dynamic perturbations induced by heme binding in cytochrome b5

The water-soluble domain of rat hepatic cytochrome b(5) undergoes marked structural changes upon heme removal. The solution structure of apocytochrome b(5) shows that the protein is partially folded in the absence of the heme group, exhibiting a stable module and a disordered heme-binding loop. The quality of the apoprotein structure in solution was improved with the use of heteronuclear NMR data. Backbone amide hydrogen exchange was studied to characterize cooperative units in the protein. It was found that this criterion distinguished the folded module from the heme-binding loop in the apoprotein, in contrast to the holoprotein. The osmolyte trimethylamine N-oxide (TMAO) did not affect the structure of the apoprotein in the disordered region. TMAO imparted a small stabilization consistent with an unfolded state effect correlating with the extent of buried surface area in the folded region of the native apoprotein. The failure of the osmolyte to cause large conformational shifts in the disordered loop supported the view that the specificity of the local sequence for the holoprotein fold was best developed with the stabilization of the native state through heme binding. To dissect the role of the heme prosthetic group in forcing the disordered region into the holoprotein conformation, the axial histidine belonging to the flexible loop (His63) was replaced with an alanine, and the structural properties of the protein with carbon-monoxide-ligated reduced iron were studied. The His63Ala substitution resulted in a protein with lower heme affinity but nevertheless capable of complete refolding. This indicated that the coordination bond was not necessary to establish the structural features of the holoprotein. In addition, the weak binding of the heme in this protein resulted in conformational shifts at a location distant from the binding site. The data suggested an uneven distribution of cooperative elements in the structure of the cytochrome.

Assignment of 15N chemical shifts and 15N relaxation measurements for oxidized and reduced iso-1-cytochrome c

A protocol for complete isotopic labeling of iso-1-cytochrome c from the eukaryote Saccharomyces cerevisiae is reported. Assignments are reported for the vast majority of the 15N amide resonances in both oxidized and reduced states. 15N heteronuclear relaxation experiments were collected to study the picosecond-nanosecond backbone dynamics of this protein. Relaxation rates were computed and fit to spectral density functions by a model-free analysis. Backbone amides in the overlapping loop B/C region are the most flexible on the picosecond-nanosecond time scale in both forms of the protein. The results show that, on average, the protein backbone is slightly more dynamic in the oxidized than the reduced state, though not significantly so. Exchange terms, which suggest significant motion on a time scale at least an order of magnitude slower than the overall correlation time of 5.2 ns, were required for only two residues in the reduced state and 27 residues in the oxidized state. When analyzed on a per-residue basis, the lower order parameters found in the oxidized state were scattered throughout the protein, with a few continuous segments found in loop C and the C-terminal helix, suggesting greater flexibility of these regions in the oxidized state. The results provide dynamic interpretations for previously presented structural and functional data, including redox-dependent changes that occur in the protein. The way is now paved for extensive dynamic analysis of variant cytochromes c.

Backbone dynamics of apocytochrome b5 in its native, partially folded state

The backbone dynamics in the native state of apocytochrome b5 were studied using 15N nuclear magnetic spin relaxation measurements. The field (11.7 and 14.1 T) and temperature (10-25 degrees C) dependence of the relaxation parameters (R1, R2, and R1rho) and the 1H-15N NOE established that the protein undergoes multiple time scale internal motions related to the secondary structure. The relaxation data were analyzed with the reduced spectral density mapping approach and within the extended model-free framework. The apoprotein was confirmed to contain a disordered heme-binding loop of approximately 30 residues with dynamics on the sub-nanosecond time scale (0.6 < S2 < 0.7, 100 ps < taue < 500 ps). This loop is attached to a structured hydrophobic core, rigid on the picosecond time scale (S2 > 0.75, taue < 50 ps). The inability to fit the data for several residues with the model-free protocol revealed the presence of correlated motion. An exchange contribution was detected in the transverse relaxation rate (R2) of all residues. The differential temperature response of R2 along the backbone supported slower exchange rates for residues in the loop (tauex > 300 micros) than for the folded polypeptide chain (tauex < 150 micros). The distribution of the reduced spectral densities at the 1H and 15N frequencies followed the dynamic trend and predicted the slowing of the internal motions at 10 degrees C. Comparison of the dynamics with those of the holoprotein [Dangi, B., Sarma, S., Yan, C., Banville, D. L., and Guiles, R. D. (1998) Biochemistry 37, 8289-8302] demonstrated that binding of the heme alters the time scale of motions both in the heme-binding loop and in the structured hydrophobic core.

Probing the backbone dynamics of oxidized and reduced rat microsomal cytochrome b5 via 15N rotating frame NMR relaxation measurements: biological implications

Rotating frame 15N relaxation NMR experiments have been performed to study the local mobility of the oxidized and reduced forms of rat microsomal cytochrome b5, in the microsecond to millisecond time range. Measurements of rotating frame relaxation rates (R1rho) were performed as a function of the effective magnetic field amplitude by using off-resonance radio frequency irradiation. Detailed analysis of the two data sets resulted in the identification of slow motions along the backbone nitrogens for both oxidation states of the protein. The local mobility of reduced and oxidized cytochrome b5 turned out to be significantly different; 28 backbone nitrogens of the oxidized form were shown to participate in a conformational exchange process, while this number dropped to 12 in the reduced form. The correlation time, tauex, for the exchange processes could be determined for 21 and 9 backbone nitrogens for oxidized and reduced cytochrome b5, respectively, with their values ranging between 70 and 280 microseconds. The direct experimental evidence provided in this study for the larger mobility of the oxidized form of the protein is consistent with the different backbone NH solvent exchangeability recently documented for the two oxidation states [Arnesano, F., et al. (1998) Biochemistry 37, 173-184]. Our experimental observations may have significant biological implications. The differential local mobility between the two oxidation states is proposed to be an important factor controlling the molecular recognition processes in which cytochrome b5 is involved.

The origin of differences in the physical properties of the equilibrium forms of cytochrome b5 revealed through high-resolution NMR structures and backbone dynamic analyses

On the basis of a comparison of high-resolution solution structures calculated for both equilibrium forms of rat ferrocytochrome b5, differences in reduction potential and thermodyanmic stability have been characterized in terms of significant structural and dynamic differences between the two forms. The dominant difference between A and B conformations has long been known to be due to a 180 degrees rotation of the heme in the binding pocket about an axis defined by the alpha- and gamma-meso carbons, however, the B form has not been structurally characterized until now. The most significant differences observed between the two forms were the presence of a hydrogen bond between the 7-propionate and the S64 amide in the A form but not the B form and surprisingly a displacement of the heme out of the binding pocket by 0.9 A in the B form relative to the A form. The magnitude of other factors which could contribute to the known difference in reduction potentials in the bovine protein [Walker, F. A., Emrick, D., Rivera, J. E., Hanquet, B. J., and Buttlaire, D. H. (1988) J. Am. Chem. Soc. 110, 6234-6240], such as differences in the orientation of the axial imidazoles and differences in hydrogen bond strength to the imidazoles, have been evaluated. The dominant effector of the reduction potential would appear to be the lack of the hydrogen bond to the S64 amide in the B form which frees up the propionate to charge stabilize the iron in the oxidized state and thus lower the reduction potential of the B form. The structure we report for the A form, based on heteronuclear NMR restraints, involving a total of 1288 restraints strongly resembles both the X-ray crystal structure of the bovine protein and a recently reported structure for the A form of the rat protein based on homonuclear data alone [Banci, L., Bertini, I., Ferroni, F., and Rosato, A. (1997) Eur. J. Biochem. 249, 270-279]. The rmsd for the backbone atoms of the A form is 0.54 A (0.92 A for all non-hydrogens). The rmsd for the backbone of the B form is 0.51 A (0. 90 A for all non-hydrogen atoms). An analysis of backbone dynamics based on a model-free analysis of 15N relaxation data, which incorporated axially symmetric diffusion tensor modeling of the cytochrome, indicates that the protein is more rigid in the reduced state relative to the oxidized state, based on a comparison with order parameters reported for the bovine protein in the oxidized state [Kelly, G. P., Muskett, F. W., and Whitford, D. (1997) Eur. J. Biochem. 245, 349-354].

Paramagnetic Inversion of the Sign of the Interference Contribution to the Transverse Relaxation of the Imido Protons of the Coordinated Imidazoles in the Uniformly 15N-Labeled Cytochrome c3

In the spectrum of uniformly 15N-labeled cytochrome c3, the relative linewidths of the doublet peaks of the 15N-coupled imido proton of the coordinated imidazole group were reversed on oxidation. This inversion was explained by the interference relaxation process between the electron-proton dipolar and 15N-1H dipolar interactions. The inversion can be used to assign the imido protons of the coordinated imidazole groups in heme proteins. Copyright 1998 Academic Press.