Backbone conformational flexibility of the lipid modified membrane anchor of the human N-Ras protein investigated by solid-state NMR and molecular dynamics simulation

The lipid modified human N-Ras protein, implicated in human cancer development, is of particular interest due to its membrane anchor that determines the activity and subcellular location of the protein. Previous solid-state NMR investigations indicated that this membrane anchor is highly dynamic, which may be indicative of backbone conformational flexibility. This article aims to address if a dynamic exchange between three structural models exist that had been determined previously. We applied a combination of solid-state nuclear magnetic resonance (NMR) methods and replica exchange molecular dynamics (MD) simulations using a Ras peptide that represents the terminal seven amino acids of the human N-Ras protein. Analysis of correlations between the conformations of individual amino acids revealed that Cys 181 and Met 182 undergo collective conformational exchange. Two major structures constituting about 60% of all conformations could be identified. The two conformations found in the simulation are in rapid exchange, which gives rise to low backbone order parameters and nuclear spin relaxation as measured by experimental NMR methods. These parameters were also determined from two 300 ns conventional MD simulations, providing very good agreement with the experimental data.     

Structural and dynamic characterization of an unfolded state of poplar apo-plastocyanin formed under nondenaturing conditions

Plastocyanin is a predominantly beta-sheet protein containing a type I copper center. The conformational ensemble of a denatured state of apo-plastocyanin formed in solution under conditions of low salt and neutral pH has been investigated by multidimensional heteronuclear NMR spectroscopy. Chemical shift assignments were obtained by using three-dimensional triple-resonance NMR experiments to trace through-bond heteronuclear connectivities along the backbone and side chains. The (3)J(HN,Halpha) coupling constants, (15)N-edited proton-proton nuclear Overhauser effects (NOEs), and (15)N relaxation parameters were also measured for the purpose of structural and dynamic characterization. Most of the residues corresponding to beta-strands in the folded protein exhibit small upfield shifts of the (13)C(alpha) and (13)CO resonances relative to random coil values, suggesting a slight preference for backbone dihedral angles in the beta region of (phi,psi) space. This is further supported by the presence of strong sequential d(alphaN)(i, i + 1) NOEs throughout the sequence. The few d(NN)(i, i + 1) proton NOEs that are observed are mostly in regions that form loops in the native plastocyanin structure. No medium or long-range NOEs were observed. A short sequence, between residues 59 and 63, was found to populate a nonnative helical conformation in the unfolded state, as indicated by the shift of the (13)C(alpha), (13)CO, and (1)H(alpha) resonances relative to random coil values and by the decreased values of the (3)J(HN,Halpha) coupling constants. The (15)N relaxation parameters indicate restriction of motions on a nanosecond timescale in this region. Intriguingly, this helical conformation is present in a sequence that is close to but not in the same location as the single short helix in the native folded protein. The results are consistent with earlier NMR studies of peptide fragments of plastocyanin and confirm that the regions of the sequence that form beta-strands in the native protein spontaneously populate the beta-region of (phi,psi) space under folding conditions, even in the absence of stabilizing tertiary interactions. We conclude that the state of apo-plastocyanin present under nondenaturing conditions is a noncompact unfolded state with some evidence of nativelike and nonnative local structuring that may be initiation sites for folding of the protein.     

Conformational properties of alpha-synuclein in its free and lipid-associated states

alpha-Synuclein (alphaS) is a presynaptic terminal protein that is believed to play an important role in the pathogenesis of Parkinson's disease (PD). We have used NMR spectroscopy to characterize the conformational properties of alphaS in solution as a free monomer and when bound to lipid vesicles and lipid-mimetic detergent micelles. Free wild-type alphaS is largely unfolded in solution, but exhibits a region with a preference for helical conformations that may be important in the aggregation of alphaS into fibrils. The N-terminal region of alphaS binds to synthetic lipid vesicles and detergent micelles in vitro and adopts a highly helical conformation, consistent with predictions based on sequence analysis. The C-terminal part of the protein does not associate with either vesicles or micelles, remaining free and unfolded. These results suggest that one function of alphaS may be to tether as of yet unidentified partners to lipid surfaces via interactions with its C-terminal tail.     

Preferred proline puckerings in cis andtrans peptide groups: Implications for collagen stability

The interplay between side-chain and main-chain conformations is a distinctive characteristic of proline residues. Here we report the results of a statistical analysis of proline conformations using a large protein database. In particular, we found that proline residues with the preceding peptide bond in the cis state preferentially adopt a down puckering. Indeed, out of 178 cis proline residues, as many as 145 (81%) are down. By analyzing the 1-4 and 1-5 nonbonding distances between backbone atoms, we provide a structural explanation for the observed trend. The observed correlation between proline puckering and peptide bond conformation suggests a new mechanism to explain the reported shift of the cis-trans equilibrium in proline derivatives. The implications of these results for the current models of collagen stability are also discussed.     

Characterisation of the conformational properties of urea-unfolded Im7: implications for the early stages of protein folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.     

Stereospecific synthesis and bio-activity of novel beta(3)-adrenoceptor agonists and inverse agonists

Since it is widely distributed into the body, beta(3)-adrenoceptor is becoming an attractive target for the treatment of several pathologies such as obesity, type 2 diabetes, metabolic syndrome, cachexia, overactive bladder, ulcero-inflammatory disorder of the gut, preterm labour, anxiety and depressive disorders, and heart failure. New compounds belonging to the class of arylethanolamines bearing one or two stereogenic centres were prepared in good yields as racemates and optically active forms. They were, then, evaluated for their intrinsic activity towards beta(3)-adrenoceptor and their affinity for beta(1)- and beta(2)-adrenergic receptors. Stereochemical features were found to play a crucial role in determining the behaviour of such compounds. In particular, alpha-racemic, (alphaR)- and (alphaS)-2-{4-[2-(2-hydroxy-2-phenylethylamino)ethyl]phenoxy}-2- methylpropanoic acid, (alpha-rac, beta-rac)-, (alphaR, betaS)- and (alphaR, betaR)- 2-{4-[2-(2-hydroxy-2-phenylethylamino)ethyl]phenoxy}propanoic acid were found to be endowed with beta(3)-adrenoceptor agonistic activity. Whereas, (alphaS, betaS)- and (alphaS, betaR)-2-{4-[2-(2-hydroxy-2-phenylethylamino)ethyl]phenoxy}propanoic acid behaved as beta(3)-adrenoceptor inverse agonists. Such compounds showed no affinity for beta(1)- and beta(2)-adrenergic receptor, respectively. Thus, resulting highly selective beta(3)-adrenoceptor ligands.     

The beta-turn scaffold of tripeptide containing an azaphenylalanine residue

The conformational preferences of azaphenylalanine-containing peptide were investigated using a model compound, Ac-azaPhe-NHMe with ab initio method at the HF/3-21G and HF/6-31G(*) levels, and the seven minimum energy conformations with trans orientation of acetyl group and the 4 minimum energy conformations with cis orientation of acetyl group were found at the HF/6-31G(*) level if their mirror images were not considered. An average backbone dihedral angle of the 11 minimum energy conformations is phi=+/-91 degrees +/-24 degrees , psi =+/-18 degrees +/-10 degrees (or +/-169 degrees +/-8 degrees ), corresponding to the i+2 position of beta-turn (delta(R)) or polyproline II (beta(P)) structure, respectively. The chi(1) angle in the aromatic side chain of azaPhe residue adopts preferentially between +/-60 degrees and +/-130 degrees, which reflect a steric hindrance between the N-terminal carbonyl group or the C-terminal amide group and the aromatic side chain with respect to the configuration of the acetyl group. These conformational preferences of Ac-azaPhe-NHMe predicted theoretically were compared with those of For-Phe-NHMe to characterize the structural role of azaPhe residue. Four tripeptides containing azaPhe residue, Boc-Xaa-azaPhe-Ala-OMe [Xaa=Gly(1), Ala(2), Phe(3), Asn(4)] were designed and synthesized to verify whether the backbone torsion angles of azaPhe reside are still the same as compared with theoretical conformations and how the preceding amino acids of azaPhe residue perturb the beta-turn skeleton in solution. The solution conformations of these tripeptide models containing azaPhe residue were determined in CDCl(3) and DMSO solvents using NMR and molecular modeling techniques. The characteristic NOE patterns, the temperature coefficients of amide protons and small solvent accessibility for the azapeptides 1-4 reveal to adopt the beta-turn structure. The structures of azapeptides containing azaPhe residue from a restrained molecular dynamics simulation indicated that average dihedral angles [(phi(1), psi(1)), (phi(2), psi(2))] of Xaa-azaPhe fragment in azapeptide, Boc-Xaa-azaPhe-Ala-OMe were [(-68 degrees, 135 degrees ), (116 degrees, -1 degrees )], and this implies that the intercalation of an azaPhe residue in tripeptide induces the betaII-turn conformation, and the volume change of a preceding amino acid of azaPhe residue in tripeptides would not perturb seriously the backbone dihedral angle of beta-turn conformation. We believe such information could be critical in designing useful molecules containing azaPhe residue for drug discovery and peptide engineering.     

NMR analysis of a kinetically trapped intermediate of a disulfide-deficient mutant of the starch-binding domain of glucoamylase

A thermally unfolded disulfide-deficient mutant of the starch-binding domain of glucoamylase refolds into a kinetically trapped metastable intermediate when subjected to a rapid lowering of temperature. We attempted to characterise this intermediate using multidimensional NMR spectroscopy. The ¹H-¹⁵N heteronuclear single quantum coherence spectrum after a rapid temperature decrease (the spectrum of the intermediate) showed good chemical shift dispersion but was significantly different from that of the native state, suggesting that the intermediate adopts a nonnative but well-structured conformation. Large chemical shift changes for the backbone amide protons between the native and the intermediate states were observed for residues in the β-sheet consisting of strands 2, 3, 5, 6, and 7 as well as in the C-terminal region. These residues were found to be in close proximity to aromatic residues, suggesting that the chemical shift changes are mainly due to ring current shifts caused by the aromatic residues. The two-dimensional nuclear Overhauser enhancement (NOE) spectroscopy experiments showed that the intermediate contained substantial, native-like NOE connectivities, although there were fewer cross peaks in the spectrum of the intermediate compared with that of the native state. It was also shown that there were native-like interresidue NOEs for residues buried in the protein, whereas many of the NOE cross peaks were lost for the residues involved in a surface-exposed aromatic cluster. These results suggest that, in the intermediate, the aromatic cluster at the surface is structurally less organised, whereas the interior of the protein has relatively rigid, native-like side-chain packing.     

A structural and functional role for 11-mer repeats in alpha-synuclein and other exchangeable lipid binding proteins

We have used NMR spectroscopy and limited proteolysis to characterize the structural properties of the Parkinson's disease-related protein alpha-synuclein in lipid and detergent micelle environments. We show that the lipid or micelle surface-bound portion of the molecule adopts a continuously helical structure with a single break. Modeling alphaS as an ideal alpha-helix reveals a hydrophobic surface that winds around the helix axis in a right-handed fashion. This feature is typical of 11-mer repeat containing sequences that adopt right-handed coiled coil conformations. In order to bind a flat or convex lipid surface, however, an unbroken helical alphaS structure would need to adopt an unusual, slightly unwound, alpha11/3 helix conformation (three complete turns per 11 residues). The break we observe in the alphaS helix may allow the protein to avoid this unusual conformation by adopting two shorter stretches of typical alpha-helical structure. However, a quantitative analysis suggests the possibility that the alpha11/3 conformation may in fact exist in lipid-bound alphaS. We discuss how structural features of helical 11-mer repeats could play a role in the reversible lipid binding function of alpha-synuclein and generalize this argument to include the 11-mer repeat-containing apolipoproteins, which also require the ability to release readily from lipid surfaces. A search of protein sequence databases confirms that synuclein-like 11-mer repeats are present in other proteins that bind lipids reversibly and predicts such a role for a number of hypothetical proteins of unknown function.     

Solution structure of the state 1 conformer of GTP-bound H-Ras protein and distinct dynamic properties between the state 1 and state 2 conformers

Ras small GTPases undergo dynamic equilibrium of two interconverting conformations, state 1 and state 2, in the GTP-bound forms, where state 2 is recognized by effectors, whereas physiological functions of state 1 have been unknown. Limited information, such as static crystal structures and (31)P NMR spectra, was available for the study of the conformational dynamics. Here we determine the solution structure and dynamics of state 1 by multidimensional heteronuclear NMR analysis of an H-RasT35S mutant in complex with guanosine 5'-(β, γ-imido)triphosphate (GppNHp). The state 1 structure shows that the switch I loop fluctuates extensively compared with that in state 2 or H-Ras-GDP. Also, backbone (1)H,(15)N signals for state 2 are identified, and their dynamics are studied by utilizing a complex with c-Raf-1. Furthermore, the signals for almost all the residues of H-Ras·GppNHp are identified by measurement at low temperature, and the signals for multiple residues are found split into two peaks corresponding to the signals for state 1 and state 2. Intriguingly, these residues are located not only in the switch regions and their neighbors but also in the rigidly structured regions, suggesting that global structural rearrangements occur during the state interconversion. The backbone dynamics of each state show that the switch loops in state 1 are dynamically mobile on the picosecond to nanosecond time scale, and these mobilities are significantly reduced in state 2. These results suggest that multiconformations existing in state 1 are mostly deselected upon the transition toward state 2 induced by the effector binding.     

Structural effects of O-glycosylation on a 15-residue peptide from the mucin (MUC1) core protein

To study the effect of O-glycosylation on the conformational propensities of a peptide backbone, the 15-residue peptide PPAHGVTSAPDTRPA (PPA15) from the MUC1 protein core and its analogue PPA15(T7), glycosylated with alpha-N-acetylgalactosamine on Thr7, were prepared and investigated by NMR spectroscopy. The peptide contains both the GVTSAP sequence, which is an effective substrate for GalNAc-T1 and -T3 transferases, and the PDTRP fragment, which is a well-known immunodominant epitope recognized by several anti-MUC1 monoclonal antibodies. Useful structural results were obtained in water upon decreasing the temperature to 5-10 degrees C. The sugar attachment slightly affected the conformational equilibrium of the peptide backbone near the glycosylated Thr7 residue. The clustering of low-energy conformations for both PPA15 and PPA15(T7) within the GVTSAP and APDTRP fragments revealed structural similarities between glycosylated and nonglycosylated peptides. For the GVTSAP region, minor but distinct clusters formed by either PPA15 or PPA15(T7) conformers showed distinct structural propensities of the peptide backbone specific for either the nonglycosylated or the glycosylated peptide. The peptide backbone of the APDTRP fragment, which is a well-known immunodominant region, resembled an S-shaped bend. A similar structural motif was found in the GVTSAP fragment. The S-shaped structure of the peptide backbone is formed by consecutive inverse gamma-turn conformations partially stabilized by hydrogen bonding. A comparison of the solution structure of the APDTRP fragment with a crystal structure of the MUC1 peptide antigen bound to the breast tumor-specific antibody SM3 demonstrated significant structural similarities in the general shape.     

The Role of Aromatic–Aromatic Interactions in Strand–Strand Stabilization of β-Sheets

Aromatic–aromatic interactions have long been believed to play key roles in protein structure, folding, and binding functions. However, we still lack full understanding of the contributions of aromatic–aromatic interactions to protein stability and the timing of their formation during folding. Here, using an aromatic ladder in the β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), as a case study, we find that aromatic π stacking plays a greater role in the Phe65–Phe71 cross-strand pair, while in another pair, Phe50–Phe65, hydrophobic interactions are dominant. The Phe65–Phe71 pair spans β-strands 4 and 5 in the β-barrel, which lack interstrand hydrogen bonding, and we speculate that it compensates energetically for the absence of strand–strand backbone interactions. Using perturbation analysis, we find that both aromatic–aromatic pairs form after the transition state for folding of CRABP1, thus playing a role in the final stabilization of the β-sheet rather than in its nucleation as had been earlier proposed. The aromatic interaction between strands 4 and 5 in CRABP1 is highly conserved in the intracellular lipid-binding protein (iLBP) family, and several lines of evidence combine to support a model wherein it acts to maintain barrel structure while allowing the dynamic opening that is necessary for ligand entry. Lastly, we carried out a bioinformatics analysis and found 51 examples of aromatic–aromatic interactions across non-hydrogen-bonded β-strands outside the iLBPs, arguing for the generality of the role played by this structural motif.     

Determination of the Individual Roles of the Linker Residues in the Interdomain Motions of Calmodulin Using NMR Chemical Shifts

Many protein molecules are formed by two or more domains whose structures and dynamics are closely related to their biological functions. It is thus important to develop methods to determine the structural properties of these multidomain proteins. Here, we characterize the interdomain motions in the calcium-bound state of calmodulin (Ca^2 +-CaM) using NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations. We find that the conformational fluctuations of the interdomain linker, which are largely responsible for the overall interdomain motions of CaM, can be well described by exploiting the information provided by chemical shifts. We thus identify 10 residues in the interdomain linker region that change their conformations upon substrate binding. Five of these residues (Met76, Lys77, Thr79, Asp80 and Ser81) are highly flexible and cover the range of conformations observed in the substrate-bound state, while the remaining five (Arg74, Lys75, Asp78, Glu82 and Glu83) are much more rigid and do not populate conformations typical of the substrate-bound form. The ensemble of conformations representing the Ca^2 +-CaM state obtained in this study is in good agreement with residual dipolar coupling, paramagnetic resonance enhancement, small-angle X-ray scattering and fluorescence resonance energy transfer measurements, which were not used as restraints in the calculations. These results provide initial evidence that chemical shifts can be used to characterize the conformational fluctuations of multidomain proteins.     

N-terminal region of alpha-synuclein is essential for the fatty acid-induced oligomerization of the molecules

Exposure of alpha-synuclein (alphaS), a major component of Lewy bodies in Parkinson's disease, to polyunsaturated fatty acids (PUFAs) triggers the formation of soluble alphaS oligomers. Here, we demonstrate that PUFA binds recombinant alphaS protein through its N-terminal region (residues 2-60). In HEK293 cells, alphaS mutants lacking the N-terminal region failed to form oligomers in the presence of PUFA. The PUFA-induced alphaS oligomerization was accelerated by C-terminal truncation or Ser129 phosphorylation of alphaS; however, this effect was abolished by deletion of the N-terminus. The results indicate that the N-terminus of alphaS is essential for the PUFA-induced alphaS oligomerization.     

Amide proton exchange rates of oxidized and reduced Saccharomyces cerevisiae iso-1-cytochrome c.

Proton NMR spectroscopy was used to determine the rate constant, kobs, for exchange of labile protons in both oxidized (Fe(III)) and reduced (Fe(II)) iso-1-cytochrome c. We find that slowly exchanging backbone amide protons tend to lack solvent-accessible surface area, possess backbone hydrogen bonds, and are present in regions of regular secondary structure as well as in omega-loops. Furthermore, there is no correlation between kobs and the distance from a backbone amide nitrogen to the nearest solvent-accessible atom. These observations are consistent with the local unfolding model. Comparisons of the free energy change for denaturation, delta Gd, at 298 K to the free energy change for local unfolding, delta Gop, at 298 K for the oxidized protein suggest that certain conformations possessing higher free energy than the denatured state are detected at equilibrium. Reduction of the protein results in a general increase in delta Gop. Comparisons of delta Gd to delta Gop for the reduced protein show that the most open states of the reduced protein possess more structure than its chemically denatured form. This persistent structure in high-energy conformations of the reduced form appears to involve the axially coordinated heme.     

(1)H/(15)N heteronuclear NMR spectroscopy shows four dynamic domains for phospholamban reconstituted in dodecylphosphocholine micelles

We report the backbone dynamics of monomeric phospholamban in dodecylphosphocholine micelles using (1)H/(15)N heteronuclear NMR spectroscopy. Phospholamban is a 52-amino acid membrane protein that regulates Ca-ATPase in cardiac muscle. Phospholamban comprises three structural domains: a transmembrane domain from residues 22 to 52, a connecting loop from 17 to 21, and a cytoplasmic domain from 1 to 16 that is organized in an "L"-shaped structure where the transmembrane and the cytoplasmic domain form an angle of approximately 80 degrees (Zamoon et al., 2003; Mascioni et al., 2002). T(1), T(2), and (1)H/(15)N nuclear Overhauser effect values measured for the amide backbone resonances were interpreted using the model-free approach of Lipari and Szabo. The results point to the existence of four dynamic domains, revealing the overall plasticity of the cytoplasmic helix, the flexible loop, and part of the transmembrane domain (residues 22-30). In addition, using Carr-Purcell-Meiboom-Gill-based experiments, we have characterized phospholamban dynamics in the micros-ms timescale. We found that the majority of the residues in the cytoplasmic domain, the flexible loop, and the first ten residues of the transmembrane domain undergo dynamics in the micros-ms range, whereas minimal dynamics were detected for the transmembrane domain. Hydrogen/deuterium exchange factors measured at different temperatures support the existence of slow motion in both the loop and the cytoplasmic helix. We propose that these dynamic properties are critical factors in the biomolecular recognition of phospholamban by Ca-ATPase and other interacting proteins such as protein kinase A and protein phosphatase 1.     

Secondary structure and dynamics of micelle bound beta- and gamma-synuclein

We have used solution state NMR spectroscopy to characterize the secondary structure and backbone dynamics of the proteins beta- and gamma-synuclein in their detergent micelle-bound conformations. Comparison of the results with those previously obtained for the Parkinson's disease-linked protein alpha-synuclein shows that structural differences between the three homologous synuclein family members are directly related to variations in their primary amino acid sequences. An 11-residue deletion in the lipid-binding domain of beta-synuclein leads to the destabilization of an entire segment of the micelle-bound helical structure containing the deletion site. The acidic C-terminal tail region of gamma-synuclein, which displays extensive sequence divergence, is more highly disordered than the corresponding regions in the other two family members. The observed structural differences are likely to mediate functional variations between the three proteins, with differences between alpha- and beta-synuclein expected to revolve around their lipid interactions, while differences in gamma-synuclein function are expected to result from different protein-protein interactions mediated by its unique C-terminal tail.     

pK(a) values for side-chain carboxyl groups of a PGB1 variant explain salt and pH-dependent stability

Determination of pK(a) values of titrating residues in proteins provides a direct means of studying electrostatic coupling as well as pH-dependent stability. The B1 domain of protein G provides an excellent model system for such investigations. In this work, we analyze the observed pK(a) values of all carboxyl groups in a variant of PGB1 (T2Q, N8D, N37D) at low and high ionic strength as determined using (1)H-(13)C heteronuclear NMR in a structural context. The pK(a) values are used to calculate the pH-dependent stability in low and high salt and to investigate electrostatic coupling in the system. The observed pK(a) values can explain the pH dependence of protein stability but require pK(a) shifts relative to model values in the unfolded state, consistent with persistent residual structure in the denatured state. In particular, we find that most of the deviations from the expected random coil values can be explained by a significantly upshifted pK(a) value. We show also that (13)C backbone carbonyl data can be used to study electrostatic coupling in proteins and provide specific information on hydrogen bonding and electrostatic potential at nontitrating sites.     

Solution structure of the B form of oxidized rat microsomal cytochrome b5 and backbone dynamics via 15N rotating-frame NMR-relaxation measurements. Biological implications.

Cytochrome b5 in solution has two isomers (A and B) differing by a 180 degrees rotation of the protoporphyrin IX plane around the axis defined by the alpha and gamma meso protons. Homonuclear and heteronuclear NMR spectroscopy has been employed in order to solve the solution structure of the minor (B) form of the oxidized state of the protein and to probe its backbone dynamics in the microsecond--ms timescale in both oxidation states. A family of 40 conformers has been obtained using 1302 meaningful NOEs and 220 pseudocontact shifts and is characterized by high quality and good resolution (rmsd to the mean structure of 0.055 +/- 0.009 nm and 0.103 +/- 0.011 nm for backbone and heavy atoms, respectively). Extensive comparisons of the structural and dynamics changes associated with the A-to-B form interconversion for both oxidation states were subsequently performed. Propionate 6 experiences a redox-state-dependent reorientation as does propionate 7 in the A form. Significant insights are obtained into the role of the protein frame for efficient biological function and backbone mobility is proposed to be one of the factors that could control the reduction potential of the heme.