Effects of surface amino acid replacements in cytochrome c peroxidase on complex formation with cytochrome c.

Site-directed mutagenesis was employed to examine the role played by specific surface residues in the activity of cytochrome c peroxidase. The double charge, aspartic acid to lysine, point mutations were constructed at positions 37, 79, and 217 on the surface of cytochrome c peroxidase, sites purported to be within or proximal to the recognition site for cytochrome c in an electron-transfer productive complex formed by the two proteins. The resulting mutant peroxidases were examined for catalytic activity by steady-state measurements and binding affinity by two methods, fluorescence binding titration and cytochrome c affinity chromatography. The cloned peroxidases exhibit similar UV-visible spectra to the wild-type yeast protein, indicating that there are no major structural differences between the cloned peroxidases and the wild-type enzyme. The aspartic acid to lysine mutations at positions 79 and 217 exhibited similar turnover numbers and binding affinities to that seen for the "wild type-like" cloned peroxidase. The same change at position 37 caused more than a 10-fold decrease in both turnover of and binding affinity for cytochrome c. This empirical finding localizes a primary recognition region critical to the dynamic complex. Models from the literature proposing structures for the complex between peroxidase and cytochrome c are discussed in light of these findings.     

A pathogenic cytochrome b mutation reveals new interactions between subunits of the mitochondrial bc1 complex

Energy transduction in mitochondria involves five oligomeric complexes embedded within the inner membrane. They are composed of catalytic and noncatalytic subunits, the role of these latter proteins often being difficult to assign. One of these complexes, the bc1 complex, is composed of three catalytic subunits including cytochrome b and seven or eight noncatalytic subunits. Recently, several mutations in the human cytochrome b gene have been linked to various diseases. We have studied in detail the effects of a cardiomyopathy generating mutation G252D in yeast. This mutation disturbs the biogenesis of the bc1 complex at 36 degrees C and decreases the steady-state level of the noncatalytic subunit Qcr9p. In addition, the G252D mutation and the deletion of QCR9 show synergetic defects that can be partially bypassed by suppressor mutations at position 252 and by a new cytochrome b mutation, P174T. Altogether, our results suggest that the supernumerary subunit Qcr9p enhances or stabilizes the interactions between the catalytic subunits, this role being essential at high temperature.     

Structural complementarity of distance constraints obtained from chemical cross-linking and amino acid coevolution

The analysis of amino acid coevolution has emerged as a practical method for protein structural modeling by providing structural contact information from alignments of amino acid sequences. In parallel, chemical cross-linking/mass spectrometry (XLMS) has gained attention as a universally applicable method for obtaining low-resolution distance constraints to model the quaternary arrangements of proteins, and more recently even protein tertiary structures. Here, we show that the structural information obtained by XLMS and coevolutionary analysis are effectively complementary: the distance constraints obtained by each method are almost exclusively associated with non-coincident pairs of residues, and modeling results obtained by the combination of both sets are improved relative to considering the same total number of constraints of a single type. The structural rationale behind the complementarity of the distance constraints is discussed and illustrated for a representative set of proteins with different sizes and folds.     

Chemical characterization of protein-protein interactions between cytochrome P-450 and cytochrome b5

Native cytochrome b5 interacts with either RLM5 or LM2 to form tight equimolar complexes (Kd = 250 and 540 nM, respectively) in which the content of high spin cytochrome P-450 was substantially increased. Cytochrome b5 caused 3- and 7-fold increases in the binding affinities of RLM5 and LM2 for benzphetamine, respectively, and benzphetamine decreased the apparent Kd for cytochrome b5 binding. Upon formation of the ternary complex between cytochromes P-450, b5, and benzphetamine the percentage of cytochrome P-450 in the high spin state was increased from 28 to 74 (RLM5) and from 9 to 85 (LM2). Cytochrome b5 caused 13- and 7-fold increases in the rate of RLM5- and LM2-dependent p-nitroanisole demethylation, respectively. Amino-modified (ethyl acetimidate or acetic anhydride) cytochrome b5 produced results similar to those obtained above with native cytochrome b5. In contrast, modification of as few as 5 mol of carboxyl groups/mol of amidinated cytochrome b5 resulted in both a substantial loss of the spectrally observed interactions with either cytochrome P-450 LM2 or cytochrome P-450 RLM5, and in a loss of the cytochrome b5-mediated stimulation of p-nitroanisole demethylation catalyzed by either monooxygenase. In further studies, native and fully acetylated cytochromes b5 reoxidized carbonmonoxy ferrous LM2 at least 20 times faster than amidinated, carboxyl-modified cytochrome b5 derivatives. In contrast, amidination, or acetylation of amino groups, or amidination of amino groups plus methylamidination of the carboxyl groups did not appreciably slow the rate of reduction of the cytochrome b5 by NADPH-cytochrome P-450 reductase. Collectively, the results provide strong evidence for an essential role of cytochrome b5 carboxyl groups in functional interactions with RLM5 and LM2.     

Rapid redox equilibrium between the mitochondrial Q pool and cytochrome b during triphasic reduction of cytochrome b by succinate

The reliability of monitoring the redox reactions of cytochrome b using the different wavelengths employed by different authors has been reexamined. It was found that 562-575 nm is suitable in succinate: cytochrome c reductase but not in mitochondria, in which case 562-540 nm is a better pair. Direct optical measurements of the redox reaction kinetics of the mitochondrial Q pool using a commercial dual-wavelength spectrophotometer are possible when succinate is used as the electron donor. Using the correct wavelength pair, and with malonate to slow down the electron input, the reduction course of cytochrome b was still triphasic but a plateau or a turn replaced the oxidation phase previously reported by several authors. At the same time, the reduction course of the Q pool was also triphasic, and in perfect match with that of cytochrome b. Destruction of the Rieske iron-sulfur cluster by British anti-Lewisite (BAL) + O2 treatment or prereduction of the high-potential components made the reduction of both Q and b monophasic. The plot of log (Q/QH2) against log (b3+/b2+) gave a straight line with an n value of 1.7 for cytochrome b at pH 7.4. This n value rose to 2.0 at pH 6.5 and dropped to 1.4 at pH 8.5. On the other hand, the mid-point potential of cytochrome b relative to that of the Q pool remained essentially unchanged between pH 6.5 and 8.4. BAL treatment had a small effect on the midpoint potential of cytochrome b relative to that of the Q pool and had no effect on the n value. Addition of quinone homologues and analogues extended the plateau phase in the reduction of cytochrome b, but exogenous quinones did not equilibrate rapidly with cytochrome b. It was concluded that the appearance of the plateau between the two reduction phases of Q and b is caused by the rapid delivery of electrons to the high-potential components of the respiratory chain as envisaged in the Q cycle; the unexpected n value for cytochrome b suggests a concerted reduction by QH2 of two species of cytochromes b-562.     

Genetic studies of the lac repressor. XIII. Extensive amino acid replacements generated by the use of natural and synthetic nonsense suppressors

We have altered the amino acid sequence of the lac repressor one residue at a time by utilizing a collection of nonsense suppressors that permit the insertion of 13 different amino acids in response to the amber (UAG) codon, as well as an additional amino acid in response to the UGA codon. We used this collection to suppress nonsense mutations at 141 positions in the lacI gene, which encodes the 360 amino acid long lac repressor, including 53 new nonsense mutations which we constructed by oligonucleotide-directed mutagenesis. This method has generated over 1600 single amino acid substitutions in the lac repressor. We have cataloged the effects of these replacements and have interpreted the results with the objective of gaining a better understanding of lac repressor structure, and protein structure in general. The DNA binding domain of the repressor, involving the amino-terminal 59 amino acids, is extremely sensitive to substitution, with 70% of the replacements resulting in the I- phenotype. However, the remaining 301 amino acid core of the repressor is strikingly tolerant of substitutions, with only 30% of the amino acids introduced causing the I- phenotype. This analysis reveals the location of sites in the protein involved in inducer binding, tighter binding to operator and thermal stability, and permits a virtual genetic image reconstruction of the lac repressor protein.     

Effects of charged amino-acid mutation on the solution structure of cytochrome b5 and binding between cytochrome b5 and cytochrome c

The solution structure of oxidized bovine microsomal cytochrome b(5) mutant (E48, E56/A, D60/A) has been determined through 1524 meaningful nuclear Overhauser effect constraints together with 190 pseudocontact shift constraints. The final family of 35 conformers has rmsd values with respect to the mean structure of 0.045+/-0.009 nm and 0.088+/-0.011 nm for backbone and heavy atoms, respectively. A characteristic of this mutant is that of having no significant changes in the whole folding and secondary structure compared with the X-ray and solution structures of wild-type cytochrome b(5). The binding of different surface mutants of cytochrome b(5) with cytochrome c shows that electrostatic interactions play an important role in maintaining the stability and specificity of the protein complex formed. The differences in association constants demonstrate the electrostatic contributions of cytochrome b(5) surface negatively charged residues, which were suggested to be involved in complex formation in the Northrup and Salemme models, have cumulative effect on the stability of cyt c-cyt b(5) complex, and the contribution of Glu48 is a little higher than that of Glu44. Moreover, our result suggests that the docking geometry proposed by Northrup, which is involved in the participation of Glu48, Glu56, Asp60, and heme propionate of cytochrome b(5), do occur in the association between cytochrome b(5) and cytochrome c.     

Calculation of the three-dimensional structure of Saccharomyces cerevisiae cytochrome b inserted in a lipid matrix

Cytochrome b is an integral membrane protein, which forms the core of the ubiquinol-cytochrome c oxidoreductase (cytochrome bc1) complex. A computer-aided three-dimensional modeling procedure was carried out in four steps. First, the candidate hydrophobic helices were searched for throughout the protein primary sequence by a computer procedure based upon the method of Eisenberg; second, a secondary helical structure was imposed to the transmembrane peptides; third, the helical segments at a lipid-water interface were oriented, and finally the possible interactions between helices with similar properties were investigated. This procedure enabled the identification of nine hydrophobic segments, of which eight are membrane-spanning helices while one has amphipathic properties. Three hydrophilic receptor-binding domains were also identified. Based upon their hydrophobicity profiles, the transmembrane helices could be associated in pairs inside the lipid bilayer. In our folding model proposed for cytochrome b, all mutation sites are not only located on the same side of the membrane but are also in close proximity in the three-dimensional structure. Inhibitor resistance mutational sites which were recently characterized (di Rago, J.-P., and Colson, A.-M. (1988) J. Biol. Chem. 263, 12564-12570) have been located on this model. Moreover, the receptor-binding domains and the mutation sites are close neighbors in the three-dimensional spatial representation.     

Effects of amino acid replacements in yeast iso-1 cytochrome c on heme accessibility and intracomplex electron transfer in complexes with cytochrome c peroxidase

The kinetics of reduction of wild type and several site-specific mutants of yeast iso-1 cytochrome c (Arg-13----Ile, Gln-16----Ser, Gln-16----Lys, Lys-27----Gln, Lys-72----Asp), both free and in 1:1 complexes with yeast cytochrome c peroxidase, by free flavin semiquinones have been studied. Intramolecular one-electron transfer from the ferrous cytochromes c to the H2O2-oxidized peroxidase at both low (8 mM) and high (275 mM) ionic strengths was also studied. The accessibility of the cytochrome c heme within the electrostatically stabilized complex and the rate constants for intramolecular electron transfer at both low and high ionic strength are highly dependent on the specific amino acids present at the protein-protein interface. Importantly, replacement by uncharged amino acids of Arg or Lys residues thought to be important in orientation and/or stabilization of the electron-transfer complex resulted in increased rates of electron transfer. In all cases, an increase in ionic strengths from 8 to 275 mM also produced increased intramolecular electron-transfer rate constants. The results suggest that the electrostatically stabilized 1:1 complex is not optimized for electron transfer and that by neutralization of key positively charged residues, or by an increase in the ionic strength thereby masking the ionic interactions, the two proteins can orient themselves to allow the formation of a more efficient electron-transfer complex.     

Differences in C-terminal amino acid sequences between erythrocyte and liver cytochrome b5 isolated from pig and human. Evidence for two tissue-specific forms of cytochrome b5

Two forms of cytochrome b5, a soluble erythrocyte form and a membrane-bound liver form, were purified from pig and human, and structural differences between them were analyzed. Porcine and human erythrocyte cytochrome b5 consisted of 97 amino acid residues and contained the same catalytic domain structure (residues 1-96) as that of the corresponding liver cytochrome b5, but had one amino acid replacement at the C-terminus (residue 97). These results suggest that erythrocyte cytochrome b5 is not derived from the liver protein by proteolysis but a translational product from another distinct mRNA of cytochrome b5.     

Location of haem-binding sites in the mitochondrial cytochrome b

Amino acid sequences of 6 mitochondrial cytochromes b are analysed to locate the binding sites of the two protohaems. Four invariant histidines are found in two protein segments which span the membrane. In each, two histidines are separated by 13 residues. This would place them on the same sides of alpha-helices, and the protohaems could be sandwiched as bis-imidazole complexes between the two transmembrane segments. In this model the haems are located in different halves of the bilayer; the Fe-Fe distance is about 20 A.     

In vitro interactions between the two mitochondrial membrane proteins VDAC and cytochrome c oxidase

VDAC, a mitochondrial outer membrane channel, is involved in the control of aerobic metabolism and in apoptotic processes via numerous protein-protein interactions. To unveil those interactions, we screened a human liver cDNA library with the phage display methodology optimized to target VDAC reconstituted into a membrane environment. One positively selected clone yielded a sequence matching a part of the subunit I of human cytochrome c oxidase (COX), a mitochondrial inner membrane enzyme. Such putative interaction was never reported before. This interaction proved to be functional as evidenced by the effect of the human and yeast isoforms of VDAC on the oxidation of cytochrome c by the pure holoenzyme and by the effect of the COX epitope on VDAC permeability. Our results providing four independently obtained evidences of VDAC-COX interaction in vitro, would support a novel and potentially important level of mitochondrial regulation given the respective locations and functions of both proteins.     

Probing the differences between rat liver outer mitochondrial membrane cytochrome b5 and microsomal cytochromes b5

Two distinct forms of cytochrome b5 exist in the rat hepatocyte. One is associated with the membrane of the endoplasmic reticulum (microsomal, or Mc, cyt b5) while the other is associated with the outer membrane of liver mitochondria (OM cyt b5). Rat OM cyt b5, the only OM cyt b5 identified so far, has a significantly more negative reduction potential and is substantially more stable toward chemical and thermal denaturation than Mc cytochromes b5. In addition, hemin is kinetically trapped in rat OM cyt b5 but not in the Mc proteins. As a result, no transfer of hemin from rat OM cyt b5 to apomyoglobin is observed at pH values as low as 5.2, nor can the thermodyamically favored ratio of hemin orientational isomers be achieved under physiologically relevant conditions. These differences are striking given the similarity of the respective protein folds. A combined theoretical and experimental study has been conducted in order to probe the structural basis behind the remarkably different properties of rat OM and Mc cytochromes b5. Molecular dynamics (MD) simulations starting from the crystal structure of bovine Mc cyt b5 revealed a conformational change that exposes several internal residues to the aqueous environment. The new conformation is equivalent to the "cleft-opened" intermediate observed in a previously reported MD simulation of bovine Mc cyt b5 [Storch, E. M., and Daggett, V. (1995) Biochemistry 34, 9682-9693]. The rat OM protein does not adopt a comparable conformation in MD simulations, thus restricting access of water to the protein interior. Subsequent comparisons of the protein sequences and structures suggested that an extended hydrophobic network encompassing the side chains of Ala-18, Ile-32, Leu-36, and Leu-47 might contribute to the inability of rat OM cyt b5 to adopt the cleft-opened conformation and, hence, stabilize its fold relative to the Mc isoforms. A corresponding network is not present in bovine Mc cyt b5 because positions 18, 32, and 47, are occupied by Ser, Leu, and Arg, respectively. To probe the roles played by Ala-18, Ile-32, and Leu-47 in endowing rat OM cyt b5 with its unusual structural properties, we have replaced them with the corresponding residues in bovine Mc cyt b5. Hence, the I32L (single), A18S/L47R (double), and A18S/L47R/I32L (triple) mutants of rat OM cyt b5 were prepared. The stability of these proteins was found to decrease in the following order: WT rat OM > rat OM I32L > rat OM A18S/L47R > rat OM A18S/L47R/I32L > bovine Mc cyt b5. The decrease in stability of the rat OM protein correlates with the extent to which the hydrophobic cluster involving the side chains of residues 18, 32, 36, and 47 has been disrupted. Complete disruption of the hydrophobic network in the triple mutant is confirmed in a 2.0 A resolution crystal structure of the protein. Disruption of the hydrophobic network also facilitates hemin loss at pH 5.2 for the double and triple mutants, with the less stable triple mutant exhibiting the greater rate of hemin transfer to apomyoglobin. Finally, 1H NMR spectroscopy and side-by-side comparisons of the crystal structures of bovine Mc, rat OM, and rat OM A18S/L47R/I32L cyt b5 allowed us to conclude that the nature of residue 32 plays a key role in controlling the relative stability of hemin orientational isomers A and B in rat OM cyt b5. A similar analysis led to the conclusion that Leu-70 and Ser-71 play a pivotal role in stabilizing isomer A relative to isomer B in Mc cytochromes b5.     

Evolutionary couplings of amino acid residues reveal structure and function of bacterial signaling proteins

The genomic era along with major advances in high‐throughput sequencing technology has led to a rapid expansion of the genomic and consequently the protein sequence space. Bacterial extracytoplasmic function sigma factors have emerged as an important group of signaling proteins in bacteria involved in many regulatory decisions, most notably the adaptation to cell envelope stress. Their wide prevalence and amplification among bacterial genomes has led to sub‐group classification and the realization of diverse signaling mechanisms. Mathematical frameworks have been developed to utilize extensive protein sequence alignments to extract co‐evolutionary signals of interaction. This has proven useful in a number of different biological fields, including de novo structure prediction, protein–protein partner identification and the elucidation of alternative protein conformations for signal proteins, to name a few. The mathematical tools, commonly referred to under the name ‘Direct Coupling Analysis’ have now been applied to deduce molecular mechanisms of activation for sub‐groups of extracytoplasmic sigma factors adding to previous successes on bacterial two‐component signaling proteins. The amplification of signal transduction protein genes in bacterial genomes made them the first to be amenable to this approach but the sequences are available now to aid the molecular microbiologist, no matter their protein pathway of interest.     

Fitness effects of amino acid replacements in the beta-galactosidase of Escherichia coli

Two genetic procedures were used to obtain amino acid replacements in the lacZ-encoded beta-galactosidase in Escherichia coli. Amino acid replacements could be obtained without regard to their effects on lactase activity by selecting spontaneous mutations that relieved the strong polarity of six nonsense mutations. When streaked on MacConkey- lactose indicator plates, approximately 75% of these mutants gave strong red lactose-fermenting colonies, and 25% gave white nonfermenting colonies. Mutants from 11 other nonsense codons were isolated directly using MacConkey-lactose indicator plates, on which positive color indication requires only 0.5% of the wildtype lactase activity. Among the total of 17 codons, 25 variant beta-galactosidases were identified using electrophoresis and thermal denaturation studies. The fitness effects of these variant beta-galactosidases were determined using competition experiments conducted with lactose as the sole nutrient limiting the growth rate in chemostat cultures. Three of the replacements were deleterious, one was selectively advantageous, and the selective effects of the remaining 21 were undetectable under conditions in which the smallest detectable selection coefficient was approximately 0.4%/generation.      

Relation between the structure of oligopeptides and the amino acid sequence

The present paper deals with determination of the relationship between the order of the arrangement of amino acids in comparatively short-range oligopeptides (tetrapeptides) and their conformational potentialities. It is shown that the spatial and conformational possibilities of the tetrapeptides composed of the same amino acid residues exhibit high sensibility to their mutual arrangement, i. e. to the amino acid sequence. A detailed conformation analysis vividly demonstrated that the difference in conformational possibilities is manly determined by different conditions of realization of residual interactions. It is shown convincingly that energetic differences of the fragments are due to different interaction contributions for each of the considered fragments.     

Resistance mutations reveal the atovaquone-binding domain of cytochrome b in malaria parasites

Atovaquone represents a class of antimicrobial agents with a broad-spectrum activity against various parasitic infections, including malaria, toxoplasmosis and Pneumocystis pneumonia. In malaria parasites, atovaquone inhibits mitochondrial electron transport at the level of the cytochrome bc1 complex and collapses mitochondrial membrane potential. In addition, this drug is unique in being selectively toxic to parasite mitochondria without affecting the host mitochondrial functions. A better understanding of the structural basis for the selective toxicity of atovaquone could help in designing drugs against infections caused by mitochondria-containing parasites. To that end, we derived nine independent atovaquone-resistant malaria parasite lines by suboptimal treatment of mice infected with Plasmodium yoelii; these mutants exhibited resistance to atovaquone-mediated collapse of mitochondrial membrane potential as well as inhibition of electron transport. The mutants were also resistant to the synergistic effects of atovaquone/ proguanil combination. Sequencing of the mitochondrially encoded cytochrome b gene placed these mutants into four categories, three with single amino acid changes and one with two adjacent amino acid changes. Of the 12 nucleotide changes seen in the nine independently derived mutants 11 replaced A:T basepairs with G:C basepairs, possibly because of reactive oxygen species resulting from atovaquone treatment. Visualization of the resistance-conferring amino acid positions on the recently solved crystal structure of the vertebrate cytochrome bc1 complex revealed a discrete cavity in which subtle variations in hydrophobicity and volume of the amino acid side-chains may determine atovaquone-binding affinity, and thereby selective toxicity. These structural insights may prove useful in designing agents that selectively affect cytochrome bc1 functions in a wide range of eukaryotic pathogens.