Sensitivity of coccolithophores to carbonate chemistry and ocean acidification.

About one-third of the carbon dioxide (CO(2)) released into the atmosphere as a result of human activity has been absorbed by the oceans, where it partitions into the constituent ions of carbonic acid. This leads to ocean acidification, one of the major threats to marine ecosystems and particularly to calcifying organisms such as corals, foraminifera and coccolithophores. Coccolithophores are abundant phytoplankton that are responsible for a large part of modern oceanic carbonate production. Culture experiments investigating the physiological response of coccolithophore calcification to increased CO(2) have yielded contradictory results between and even within species. Here we quantified the calcite mass of dominant coccolithophores in the present ocean and over the past forty thousand years, and found a marked pattern of decreasing calcification with increasing partial pressure of CO(2) and concomitant decreasing concentrations of CO(3)(2-). Our analyses revealed that differentially calcified species and morphotypes are distributed in the ocean according to carbonate chemistry. A substantial impact on the marine carbon cycle might be expected upon extrapolation of this correlation to predicted ocean acidification in the future. However, our discovery of a heavily calcified Emiliania huxleyi morphotype in modern waters with low pH highlights the complexity of assemblage-level responses to environmental forcing factors.

Controls of primary production in two phytoplankton blooms in the Antarctic Circumpolar Current.

The Antarctic Circumpolar Current has a high potential for primary production and carbon sequestration through the biological pump. In the current study, two large-scale blooms observed in 2012 during a cruise with R.V. Polarstern were investigated with respect to phytoplankton standing stocks, primary productivity and nutrient budgets. While net primary productivity was similar in both blooms, chlorophyll a -specific photosynthesis was more efficient in the bloom closer to the island of South Georgia (39 °W, 50 °S) compared to the open ocean bloom further east (12 °W, 51 °S). We did not find evidence for light being the driver of bloom dynamics as chlorophyll standing stocks up to 165 mg m-2 developed despite mixed layers as deep as 90 m. Since the two bloom regions differ in their distance to shelf areas, potential sources of iron vary. Nutrient (nitrate, phosphate, silicate) deficits were similar in both areas despite different bloom ages, but their ratios indicated more pronounced iron limitation at 12 °W compared to 39 °W. While primarily the supply of iron and not the availability of light seemed to control onset and duration of the blooms, higher grazing pressure could have exerted a stronger control toward the declining phase of the blooms.

Molecular modelling of the Norrie disease protein predicts a cystine knot growth factor tertiary structure.

The X-lined gene for Norrie disease, which is characterized by blindness, deafness and mental retardation has been cloned recently. This gene has been thought to code for a putative extracellular factor; its predicted amino acid sequence is homologous to the C-terminal domain of diverse extracellular proteins. Sequence pattern searches and three-dimensional modelling now suggest that the Norrie disease protein (NDP) has a tertiary structure similar to that of transforming growth factor beta (TGF beta). Our model identifies NDP as a member of an emerging family of growth factors containing a cystine knot motif, with direct implications for the physiological role of NDP. The model also sheds light on sequence related domains such as the C-terminal domain of mucins and of von Willebrand factor.

Progress in protein structure prediction?

Prediction of protein secondary structure is an old problem and progress has been slow. Recently, spectacular success has been claimed in the blind prediction of the catalytic subunit of the cAMP-dependent protein kinase. When predictions in this and other test cases are assessed critically, some claims of prediction success turn out to be exaggerated, but a kernel of real progress remains: protein structure prediction can be improved substantially when a family of related sequences is available. Enough so that molecular biologists equipped with a new amino acid sequence and a multiple sequence alignment in hand may be tempted to test the new prediction methods.

Loss of classical transient receptor potential 6 channel reduces allergic airway response.

BACKGROUND: Non-selective cation influx through canonical transient receptor potential channels (TRPCs) is thought to be an important event leading to airway inflammation. TRPC6 is highly expressed in the lung, but its role in allergic processes is still poorly understood. OBJECTIVE: The purpose of this study was to evaluate the role of TRPC6 in airway hyperresponsiveness (AHR) and allergic inflammation of the lung. METHODS: Methacholine-induced AHR was assessed by head-out body plethysmography of wild type (WT) and TRPC6(-/-) mice. Experimental airway inflammation was induced by intraperitoneal ovalbumin (OVA) sensitization, followed by OVA aerosol challenges. Allergic inflammation and mucus production were analysed 24 h after the last allergen challenge. RESULTS: Methacholine-induced AHR and agonist-induced contractility of tracheal rings were increased in TRPC6(-/-) mice compared with WT mice, most probably due to compensatory up-regulation of TRPC3 in airway smooth muscle cells. Most interestingly, when compared with WT mice, TRPC6(-/-) mice exhibited reduced allergic responses after allergen challenge as evidenced by a decrease in airway eosinophilia and blood IgE levels, as well as decreased levels of T-helper type 2 (Th2) cytokines (IL-5, IL-13) in the bronchoalveolar lavage. However, lung mucus production after allergen challenge was not altered by TRPC6 deficiency. CONCLUSIONS: TRPC6 deficiency inhibits specific allergic immune responses, pointing to an important immunological function of this cation channel in Th2 cells, eosinophils, mast cells and B cells.

Resistance of Arctic phytoplankton to ocean acidification and enhanced irradiance.

The Arctic Ocean is a region particularly prone to ongoing ocean acidification (OA) and climate-driven changes. The influence of these changes on Arctic phytoplankton assemblages, however, remains poorly understood. In order to understand how OA and enhanced irradiances (e.g., resulting from sea-ice retreat) will alter the species composition, primary production, and eco-physiology of Arctic phytoplankton, we conducted an incubation experiment with an assemblage from Baffin Bay (71°N, 68°W) under different carbonate chemistry and irradiance regimes. Seawater was collected from just below the deep Chl a maximum, and the resident phytoplankton were exposed to 380 and 1000 µatm pCO2 at both 15 and 35% incident irradiance. On-deck incubations, in which temperatures were 6 °C above in situ conditions, were monitored for phytoplankton growth, biomass stoichiometry, net primary production, photo-physiology, and taxonomic composition. During the 8-day experiment, taxonomic diversity decreased and the diatom Chaetoceros socialis became increasingly dominant irrespective of light or CO2 levels. We found no statistically significant effects from either higher CO2 or light on physiological properties of phytoplankton during the experiment. We did, however, observe an initial 2-day stress response in all treatments, and slight photo-physiological responses to higher CO2 and light during the first five days of the incubation. Our results thus indicate high resistance of Arctic phytoplankton to OA and enhanced irradiance levels, challenging the commonly predicted stimulatory effects of enhanced CO2 and light availability for primary production.

Prediction of protein secondary structure at better than 70% accuracy.

We have trained a two-layered feed-forward neural network on a non-redundant data base of 130 protein chains to predict the secondary structure of water-soluble proteins. A new key aspect is the use of evolutionary information in the form of multiple sequence alignments that are used as input in place of single sequences. The inclusion of protein family information in this form increases the prediction accuracy by six to eight percentage points. A combination of three levels of networks results in an overall three-state accuracy of 70.8% for globular proteins (sustained performance). If four membrane protein chains are included in the evaluation, the overall accuracy drops to 70.2%. The prediction is well balanced between alpha-helix, beta-strand and loop: 65% of the observed strand residues are predicted correctly. The accuracy in predicting the content of three secondary structure types is comparable to that of circular dichroism spectroscopy. The performance accuracy is verified by a sevenfold cross-validation test, and an additional test on 26 recently solved proteins. Of particular practical importance is the definition of a position-specific reliability index. For half of the residues predicted with a high level of reliability the overall accuracy increases to better than 82%. A further strength of the method is the more realistic prediction of segment length. The protein family prediction method is available for testing by academic researchers via an electronic mail server.

Redefining the goals of protein secondary structure prediction.

Secondary structure prediction recently has surpassed the 70% level of average accuracy, evaluated on the single residue states helix, strand and loop (Q3). But the ultimate goal is reliable prediction of tertiary (three-dimensional, 3D) structure, not 100% single residue accuracy for secondary structure. A comparison of pairs of structurally homologous proteins with divergent sequences reveals that considerable variation in the position and length of secondary structure segments can be accommodated within the same 3D fold. It is therefore sufficient to predict the approximate location of helix, strand, turn and loop segments, provided they are compatible with the formation of 3D structure. Accordingly, we define here a measure of segment overlap (Sov) that is somewhat insensitive to small variations in secondary structure assignments. The new segment overlap measure ranges from an ignorance level of 37% (random protein pairs) via a current level of 72% for a prediction method based on sequence profile input to neural networks (PHD) to an average 90% level for homologous protein pairs. We conclude that the highest scores one can reasonably expect for secondary structure prediction are a single residue accuracy of Q3 > 85% and a fractional segment overlap of Sov > 90%.

Protein fold recognition by prediction-based threading.

In fold recognition by threading one takes the amino acid sequence of a protein and evaluates how well it fits into one of the known three-dimensional (3D) protein structures. The quality of sequence-structure fit is typically evaluated using inter-residue potentials of mean force or other statistical parameters. Here, we present an alternative approach to evaluating sequence-structure fitness. Starting from the amino acid sequence we first predict secondary structure and solvent accessibility for each residue. We then thread the resulting one-dimensional (1D) profile of predicted structure assignments into each of the known 3D structures. The optimal threading for each sequence-structure pair is obtained using dynamic programming. The overall best sequence-structure pair constitutes the predicted 3D structure for the input sequence. The method is fine-tuned by adding information from direct sequence-sequence comparison and applying a series of empirical filters. Although the method relies on reduction of 3D information into 1D structure profiles, its accuracy is, surprisingly, not clearly inferior to methods based on evaluation of residue interactions in 3D. We therefore hypothesise that existing 1D-3D threading methods essentially do not capture more than the fitness of an amino acid sequence for a particular 1D succession of secondary structure segments and residue solvent accessibility. The prediction-based threading method on average finds any structurally homologous region at first rank in 29% of the cases (including sequence information). For the 22% first hits detected at highest scores, the expected accuracy rose to 75%. However, the task of detecting entire folds rather than homologous fragments was managed much better; 45 to 75% of the first hits correctly recognised the fold.

Effective use of sequence correlation and conservation in fold recognition.

Protein families are a rich source of information; sequence conservation and sequence correlation are two of the main properties that can be derived from the analysis of multiple sequence alignments. Sequence conservation is related to the direct evolutionary pressure to retain the chemical characteristics of some positions in order to maintain a given function. Sequence correlation is attributed to the small sequence adjustments needed to maintain protein stability against constant mutational drift. Here, we showed that sequence conservation and correlation were each frequently informative enough to detect incorrectly folded proteins. Furthermore, combining conservation, correlation, and polarity, we achieved an almost perfect discrimination between native and incorrectly folded proteins. Thus, we made use of this information for threading by evaluating the models suggested by a threading method according to the degree of proximity of the corresponding correlated, conserved, and apolar residues. The results showed that the fold recognition capacity of a given threading approach could be improved almost fourfold by selecting the alignments that score best under the three different sequence-based approaches.

Adaptation of protein surfaces to subcellular location.

In vivo, proteins occur in widely different physio-chemical environments, and, from in vitro studies, we know that protein structure can be very sensitive to environment. However, theoretical studies of protein structure have tended to ignore this complexity. In this paper, we have approached this problem by grouping proteins by their subcellular location and looking at structural properties that are characteristic to each location. We hypothesize that, throughout evolution, each subcellular location has maintained a characteristic physio-chemical environment, and that proteins in each location have adapted to these environments. If so, we would expect that protein structures from different locations will show characteristic differences, particularly at the surface, which is directly exposed to the environment. To test this hypothesis, we have examined all eukaryotic proteins with known three-dimensional structure and for which the subcellular location is known to be either nuclear, cytoplasmic, or extracellular. In agreement with previous studies, we find that the total amino acid composition carries a signal that identifies the subcellular location. This signal was due almost entirely to the surface residues. The surface residue signal was often strong enough to accurately predict subcellular location, given only a knowledge of which residues are at the protein surface. The results suggest how the accuracy of prediction of location from sequence can be improved. We concluded that protein surfaces show adaptation to their subcellular location. The nature of these adaptations suggests several principles that proteins may have used in adapting to particular physio-chemical environments; these principles may be useful for protein design.

Structure prediction of proteins--where are we now?

Although the 'structure from sequence' prediction problem remains fundamentally unsolved, new and promising methods in one, two and three dimensions have reopened the field. Significantly improved one-dimensional prediction of secondary structure from multiple sequence alignments is now in routine use. In the two-dimensional approach, inter-residue contacts can be detected by analysis of correlated mutations, albeit with low accuracy. Finally, three-dimensional methods, in which pseudopotentials or information values are derived from the databases, are proving their value for distinguishing between correct and incorrect models.

Comparing function and structure between entire proteomes.

More than 30 organisms have been sequenced entirely. Here, we applied a variety of simple bioinformatics tools to analyze 29 proteomes for representatives from all three kingdoms: eukaryotes, prokaryotes, and archaebacteria. We confirmed that eukaryotes have relatively more long proteins than prokaryotes and archaes, and that the overall amino acid composition is similar among the three. We predicted that approximately 15%-30% of all proteins contained transmembrane helices. We could not find a correlation between the content of membrane proteins and the complexity of the organism. In particular, we did not find significantly higher percentages of helical membrane proteins in eukaryotes than in prokaryotes or archae. However, we found more proteins with seven transmembrane helices in eukaryotes and more with six and 12 transmembrane helices in prokaryotes. We found twice as many coiled-coil proteins in eukaryotes (10%) as in prokaryotes and archaes (4%-5%), and we predicted approximately 15%-25% of all proteins to be secreted by most eukaryotes and prokaryotes. Every tenth protein had no known homolog in current databases, and 30%-40% of the proteins fell into structural families with >100 members. A classification by cellular function verified that eukaryotes have a higher proportion of proteins for communication with the environment. Finally, we found at least one homolog of experimentally known structure for approximately 20%-45% of all proteins; the regions with structural homology covered 20%-30% of all residues. These numbers may or may not suggest that there are 1200-2600 folds in the universe of protein structures. All predictions are available at http://cubic.bioc.columbia.edu/genomes.

Topology prediction for helical transmembrane proteins at 86% accuracy.

Previously, we introduced a neural network system predicting locations of transmembrane helices (HTMs) based on evolutionary profiles (PHDhtm, Rost B, Casadio R, Fariselli P, Sander C, 1995, Protein Sci 4:521-533). Here, we describe an improvement and an extension of that system. The improvement is achieved by a dynamic programming-like algorithm that optimizes helices compatible with the neural network output. The extension is the prediction of topology (orientation of first loop region with respect to membrane) by applying to the refined prediction the observation that positively charged residues are more abundant in extra-cytoplasmic regions. Furthermore, we introduce a method to reduce the number of false positives, i.e., proteins falsely predicted with membrane helices. The evaluation of prediction accuracy is based on a cross-validation and a double-blind test set (in total 131 proteins). The final method appears to be more accurate than other methods published: (1) For almost 89% (+/-3%) of the test proteins, all HTMs are predicted correctly. (2) For more than 86% (+/-3%) of the proteins, topology is predicted correctly. (3) We define reliability indices that correlate with prediction accuracy: for one half of the proteins, segment accuracy raises to 98%; and for two-thirds, accuracy of topology prediction is 95%. (4) The rate of proteins for which HTMs are predicted falsely is below 2% (+/-1%). Finally, the method is applied to 1,616 sequences of Haemophilus influenzae. We predict 19% of the genome sequences to contain one or more HTMs. This appears to be lower than what we predicted previously for the yeast VIII chromosome (about 25%).

Transmembrane helices predicted at 95% accuracy.

We describe a neural network system that predicts the locations of transmembrane helices in integral membrane proteins. By using evolutionary information as input to the network system, the method significantly improved on a previously published neural network prediction method that had been based on single sequence information. The input data were derived from multiple alignments for each position in a window of 13 adjacent residues: amino acid frequency, conservation weights, number of insertions and deletions, and position of the window with respect to the ends of the protein chain. Additional input was the amino acid composition and length of the whole protein. A rigorous cross-validation test on 69 proteins with experimentally determined locations of transmembrane segments yielded an overall two-state per-residue accuracy of 95%. About 94% of all segments were predicted correctly. When applied to known globular proteins as a negative control, the network system incorrectly predicted fewer than 5% of globular proteins as having transmembrane helices. The method was applied to all 269 open reading frames from the complete yeast VIII chromosome. For 59 of these, at least two transmembrane helices were predicted. Thus, the prediction is that about one-fourth of all proteins from yeast VIII contain one transmembrane helix, and some 20%, more than one.

Carboxyl group protonation upon reduction of the Paracoccus denitrificans cytochrome c oxidase: direct evidence by FTIR spectroscopy.

The redox reactions of the cytochrome c oxidase from Paracoccus denitrificans were investigated in a thin-layer cell designed for the combination of electrochemistry under anaerobic conditions with UV/VIS and IR spectroscopy. Quantitative and reversible electrochemical reactions were obtained at a surface-modified electrode for all cofactors as indicated by the optical signals in the 400-700 nm range. Fourier transform infrared (FTIR) difference spectra of reduction and oxidation (reduced-minus-oxidized and oxidized-minus-reduced, respectively) obtained in the 1800-1000 cm(-1) range reveal highly structured band features with major contributions in the amide I (1620-1680 cm(-1)) and amide II (1580-1520 cm(-1)) range which indicate structural rearrangements in the cofactor vicinity. However, the small amplitude of the IR difference signals indicates that these conformational changes are small and affect only individual peptide groups. In the spectral region above 1700 cm(-1), a positive peak in the reduced state (1733 cm(-1)) and negative peak in the oxidized st ate (1745 cm(-1)) are characteristic for the formation and decay of a COOH mode upon reduction. The most obvious interpretation of this difference signal is proton uptake by one Asp or Glu side chain carboxyl group in the reduced state and deprotonation of another Asp or Glu residue. Moreover, both residues could well be coupled as a donor-acceptor pair in the proton transfer chain. An alternative interpretation is in terms of a protonated carboxyl group which shifts to a different environment in the reduced state. The relevance of this first direct observation of protein protonation changes in the cytochrome c oxidase for vectorial proton transfer and the catalytic reaction is discussed.

Role of transmembrane domains in the functions of B- and T-cell receptors.

The antigen receptors on the surface of B- and T-lymphocytes are complexes of several integral membrane proteins, essential for their proper expression and function. Recent studies demonstrated that transmembrane (TM) domains of the components of these receptors play a critical role in their association and function. It was specifically demonstrated that in many cases point mutations in the TM domains can partially or completely disrupt the receptor surface expression and function. Here we review studies of the TM domains of B- and T-cell receptors. Furthermore, we use a novel method, PHDtopology, to provide estimates of the exact locations and lengths of the TM domains of the subunit components of these receptors. Most previous studies used single residue hydrophobicity as a criterion for determining the position and length of the TM domains. In contrast, PHDtopology utilizes a system of neural networks and the evolutionary information contained in multiple alignments of related sequences to predict the location, length, and orientation of transmembrane helices. Present results significantly differ from most published estimates of the TM domains of the B- and T-cell receptor components, primarily in the length of the TM domains. These results may lead to modification of putative TM motifs and re-interpretation of the results of studies using mutated TM domains. The availability of PHDtopology on the Internet would make it a valuable tool in the future studies of the TM domains of integral membrane proteins.

Exercise training for rehabilitation and secondary prevention of falls in geriatric patients with a history of injurious falls.

OBJECTIVE: To determine the safety and efficacy of an exercise protocol designed to improve strength, mobility, and balance and to reduce subsequent falls in geriatric patients with a history of injurious falls. DESIGN: A randomized controlled 3-month intervention trial, with an additional 3-month follow-up. SETTING: Out-patient geriatric rehabilitation unit. PARTICIPANTS: Fifty-seven female geriatric patients (mean age 82 +/- 4.8 years; range 75-90) admitted to acute care or inpatient rehabilitation with a history of recurrent or injurious falls including patients with acute fall-related fracture. INTERVENTION: Ambulatory training of strength, functional performance, and balance 3 times per week for 3 months. Patients of the control group attended a placebo group 3 times a week for 3 months. Both groups received an identical physiotherapeutic treatment 2 times a week, in which strengthening and balance training were excluded. MEASUREMENTS: Strength, functional ability, motor function, psychological parameters, and fall rates were assessed by standardized protocols at the beginning (T1) and the end (T2) of intervention. Patients were followed up for 3 months after the intervention (T3). RESULTS: No training-related medical problems occurred in the study group. Forty-five patients (79%) completed all assessments after the intervention and follow-up period. Adherence was excellent in both groups (intervention 85.4 +/- 27.8% vs control 84.2 +/- 29.3%). The patients in the intervention group increased strength, functional motor performance, and balance significantly. Fall-related behavioral and emotional restrictions were reduced significantly. Improvements persisted during the 3-month follow-up with only moderate losses. For patients of the control group, no change in strength, functional performance, or emotional status could be documented during intervention and follow-up. Fall incidence was reduced nonsignificantly by 25% in the intervention group compared with the control group (RR:0.753 CI:0.455-1.245). CONCLUSIONS: Progressive resistance training and progressive functional training are safe and effective methods of increasing strength and functional performance and reducing fall-related behavioral and emotional restrictions during ambulant rehabilitation in frail, high-risk geriatric patients with a history of injurious falls.

[Value of contrast-enhanced ultrasound vs. CT and MRI in palpable enlarged lymph nodes of the head and neck]. / Wertigkeit der kontrastmittelverstärkten Sonographie von Lymphknotenvergrösserungen im Kopf-Hals-Bereich versus Computertomographie und Magnetresonanztomographie.

PURPOSE: To compare the value of B-mode-, plain and contrast enhanced color Doppler ultrasound, CT and MRI with respect to their diagnostic accuracy in palpable enlarged cervical lymph nodes. MATERIAL AND METHODS: Thirty patients (18 - 90 years old) with palpable enlarged lymph nodes of the head and neck underwent B-mode-ultrasound, plain and contrast enhanced color Doppler, CT and MRI (gold standard: histologic analysis in 22 and clinical follow up for at least six months in eight patients). The criteria of malignancy were maximal and minimal lymph node diameter, M/Q-ratio, various morphologic criteria (necrosis, hilus line, internal structure, contour, contrast enhancement), spectral Doppler indices, and vascular architecture in color Doppler. RESULTS: The highest sensitivity (= 1.00, specificity = 0.07 - 0.15) was obtained measuring the lymph node diameter independent on the used imaging modality (ultrasound, CT, MRI), the highest specificity (= 1.00, sensitivity = 0.71) analyzing the vascularity of the lymph node by plain color Doppler. The highest diagnostic (= 0.93) accuracy was delivered by contrast enhanced color Doppler analysis of the vascularity. Sensitivity (= 0.94) and specificity (= 0.92) of this imaging modality were only slightly inferior to the top values. Fisher's exact test revealed significant values in differentiating malignant from benign lymph nodes for B-mode- and MR-analysis of the M/Q-ratio (p < 0001/p < 0.05), B-mode morphology (p < 0.00005), plain and contrast enhanced color Doppler analysis of the vascularity (p < 0.0001/p < 0.000005), MR-morphology (p < 0.0001), and CT-morphology (p < 0.005). CONCLUSION: CT is inferior to MRI, B-mode-ultrasound and contrast-enhanced color Doppler in the differential diagnosis of selectively analyzed, palpably enlarged cervical lymph nodes using the criteria of our study. The analysis of the MR-morphology revealed a slightly inferior diagnostic accuracy to B-mode morphology and color Doppler analysis of the vascularity.

Redox dependent conformational changes in the mixed valence form of the cytochrome c oxidase from p. The reorganization of glutamic acid 278 is coupled to the electron transfer from/to heme a and the binuclear center. denitrificans.

In this work we present the separation of FTIR difference signals induced by electron transfer to/from the redox centers of the cytochrome c oxidase from P. denitrificans and compare electrochemically induced FTIR difference spectra with those induced by CO photolysis. FTIR difference spectra of rebinding of CO to the half reduced (mixed valence) form of the cytochrome c oxidase after photolysis reflect the conformational changes induced by the rebinding of CO and by electron transfer reactions from heme a3 to heme a and further on to CUA. During this process, heme a3 (and CUB) are oxidized, whereas heme a and CuA are reduced. By subtracting these difference spectra from an electrochemically induced FTIR difference spectrum, where all four cofactors are reduced, the contributions for heme a3 (and CuB) could be separated. Correspondingly, the spectral contributions of heme a and CuA have been separated. The comparison of these spectra with the spectra calculated for the hemes on the basis of their redox dependent changes previously published in Hellwig et al., (Biochemistry 38, (1999) 1685-1694) show a high degree of similarity, except for additional signals coupled to the reorganization of the binuclear center upon CO rebinding. The separated spectra clearly show that the signals attributed to Glu278, an amino acid discussed to be crucial for proton pumping, is coupled to electron transfer to/from heme a and the binuclear heme a3-CuB center.