Evidence for [2Fe-2S]2+ and Linear [3Fe-4S]1+ Clusters in a Unique Family of Glycine/Cysteine-Rich Fe-S Proteins from Megavirinae Giant Viruses.

We have discovered a protein with an amino acid composition exceptionally rich in glycine and cysteine residues in the giant virus mimivirus. This small 6 kDa protein is among the most abundant proteins in the icosahedral 0.75 µm viral particles; it has no predicted function but is probably essential for infection. The aerobically purified red-brownish protein overproduced inEscherichia coli contained both iron and inorganic sulfide. UV/vis, EPR, and Mössbauer studies revealed that the viral protein, coined GciS, accommodated two distinct Fe-S clusters: a diamagnetic S = 0 [2Fe-2S]2+ cluster and a paramagnetic S = 5/2 linear [3Fe-4S]1+ cluster, a geometry rarely stabilized in native proteins. Orthologs of mimivirus GciS were identified within all clades of Megavirinae, a Mimiviridae subfamily infecting Acanthamoeba, including the distantly related tupanviruses, and displayed the same spectroscopic features. Thus, these glycine/cysteine-rich proteins form a new family of viral Fe-S proteins sharing unique Fe-S cluster binding properties.

Iron Insertion at the Assembly Site of the ISCU Scaffold Protein Is a Conserved Process Initiating Fe-S Cluster Biosynthesis.

Iron-sulfur (Fe-S) clusters are prosthetic groups of proteins biosynthesized on scaffold proteins by highly conserved multi-protein machineries. Biosynthesis of Fe-S clusters into the ISCU scaffold protein is initiated by ferrous iron insertion, followed by sulfur acquisition, via a still elusive mechanism. Notably, whether iron initially binds to the ISCU cysteine-rich assembly site or to a cysteine-less auxiliary site via N/O ligands remains unclear. We show here by SEC, circular dichroism (CD), and Mössbauer spectroscopies that iron binds to the assembly site of the monomeric form of prokaryotic and eukaryotic ISCU proteins via either one or two cysteines, referred to the 1-Cys and 2-Cys forms, respectively. The latter predominated at pH 8.0 and correlated with the Fe-S cluster assembly activity, whereas the former increased at a more acidic pH, together with free iron, suggesting that it constitutes an intermediate of the iron insertion process. Iron not binding to the assembly site was non-specifically bound to the aggregated ISCU, ruling out the existence of a structurally defined auxiliary site in ISCU. Characterization of the 2-Cys form by site-directed mutagenesis, CD, NMR, X-ray absorption, Mössbauer, and electron paramagnetic resonance spectroscopies showed that the iron center is coordinated by four strictly conserved amino acids of the assembly site, Cys35, Asp37, Cys61, and His103, in a tetrahedral geometry. The sulfur receptor Cys104 was at a very close distance and apparently bound to the iron center when His103 was missing, which may enable iron-dependent sulfur acquisition. Altogether, these data provide the structural basis to elucidate the Fe-S cluster assembly process and establish that the initiation of Fe-S cluster biosynthesis by insertion of a ferrous iron in the assembly site of ISCU is a conserved mechanism.

OrpR is a σ54 -dependent activator using an iron-sulfur cluster for redox sensing in Desulfovibrio vulgaris Hildenborough.

Enhancer binding proteins (EBPs) are key players of σ54 -regulation that control transcription in response to environmental signals. In the anaerobic microorganism Desulfovibrio vulgaris Hildenborough (DvH), orp operons have been previously shown to be coregulated by σ54 -RNA polymerase, the integration host factor IHF and a cognate EBP, OrpR. In this study, ChIP-seq experiments indicated that the OrpR regulon consists of only the two divergent orp operons. In vivo data revealed that (i) OrpR is absolutely required for orp operons transcription, (ii) under anaerobic conditions, OrpR binds on the two dedicated DNA binding sites and leads to high expression levels of the orp operons, (iii) increasing the redox potential of the medium leads to a drastic down-regulation of the orp operons expression. Moreover, combining functional and biophysical studies on the anaerobically purified OrpR leads us to propose that OrpR senses redox potential variations via a redox-sensitive [4Fe-4S]2+ cluster in the sensory PAS domain. Overall, the study herein presents the first characterization of a new Fe-S redox regulator belonging to the σ54 -dependent transcriptional regulator family probably advantageously selected by cells adapted to the anaerobic lifestyle to monitor redox stress conditions.

Mechanism of inhibition of NiFe hydrogenase by nitric oxide.

Hydrogenases reversibly catalyze the oxidation of molecular hydrogen and are inhibited by several small molecules including O2, CO and NO. In the present work, we investigate the mechanism of inhibition by NO of the oxygen-sensitive NiFe hydrogenase from Desulfovibrio fructosovorans by coupling site-directed mutagenesis, protein film voltammetry (PFV) and EPR spectroscopy. We show that micromolar NO strongly inhibits NiFe hydrogenase and that the mechanism of inhibition is complex, with NO targeting several metallic sites in the protein. NO reacts readily at the NiFe active site according to a two-step mechanism. The first and faster step is the reversible binding of NO to the active site followed by a slower and irreversible transformation at the active site. NO also induces irreversible damage of the iron-sulfur centers chain. We give direct evidence of preferential nitrosylation of the medial [3Fe-4S] to form dinitrosyl-iron complexes.

A threonine stabilizes the NiC and NiR catalytic intermediates of [NiFe]-hydrogenase.

The heterodimeric [NiFe] hydrogenase from Desulfovibrio fructosovorans catalyzes the reversible oxidation of H2 into protons and electrons. The catalytic intermediates have been attributed to forms of the active site (NiSI, NiR, and NiC) detected using spectroscopic methods under potentiometric but non-catalytic conditions. Here, we produced variants by replacing the conserved Thr-18 residue in the small subunit with Ser, Val, Gln, Gly, or Asp, and we analyzed the effects of these mutations on the kinetic (H2 oxidation, H2 production, and H/D exchange), spectroscopic (IR, EPR), and structural properties of the enzyme. The mutations disrupt the H-bond network in the crystals and have a strong effect on H2 oxidation and H2 production turnover rates. However, the absence of correlation between activity and rate of H/D exchange in the series of variants suggests that the alcoholic group of Thr-18 is not necessarily a proton relay. Instead, the correlation between H2 oxidation and production activity and the detection of the NiC species in reduced samples confirms that NiC is a catalytic intermediate and suggests that Thr-18 is important to stabilize the local protein structure of the active site ensuring fast NiSI-NiC-NiR interconversions during H2 oxidation/production.

The megavirus chilensis Cu,Zn-superoxide dismutase: the first viral structure of a typical cellular copper chaperone-independent hyperstable dimeric enzyme.

UNLABELLED: Giant viruses able to replicate in Acanthamoeba castellanii penetrate their host through phagocytosis. After capsid opening, a fusion between the internal membranes of the virion and the phagocytic vacuole triggers the transfer in the cytoplasm of the viral DNA together with the DNA repair enzymes and the transcription machinery present in the particles. In addition, the proteome analysis of purified mimivirus virions revealed the presence of many enzymes meant to resist oxidative stress and conserved in the Mimiviridae. Megavirus chilensis encodes a predicted copper, zinc superoxide dismutase (Cu,Zn-SOD), an enzyme known to detoxify reactive oxygen species released in the course of host defense reactions. While it was thought that the metal ions are required for the formation of the active-site lid and dimer stabilization, megavirus chilensis SOD forms a very stable metal-free dimer. We used electron paramagnetic resonance (EPR) analysis and activity measurements to show that the supplementation of the bacterial culture with copper and zinc during the recombinant expression of Mg277 is sufficient to restore a fully active holoenzyme. These results demonstrate that the viral enzyme's activation is independent of a chaperone both for disulfide bridge formation and for copper incorporation and suggest that its assembly may not be as regulated as that of its cellular counterparts. A SOD protein is encoded by a variety of DNA viruses but is absent from mimivirus. As in poxviruses, the enzyme might be dispensable when the virus infects Acanthamoeba cells but may allow megavirus chilensis to infect a broad range of eukaryotic hosts. IMPORTANCE: Mimiviridae are giant viruses encoding more than 1,000 proteins. The virion particles are loaded with proteins used by the virus to resist the vacuole's oxidative stress. The megavirus chilensis virion contains a predicted copper, zinc superoxide dismutase (Cu,Zn-SOD). The corresponding gene is present in some megavirus chilensis relatives but is absent from mimivirus. This first crystallographic structure of a viral Cu,Zn-SOD highlights the features that it has in common with and its differences from cellular SODs. It corresponds to a very stable dimer of the apo form of the enzyme. We demonstrate that upon supplementation of the growth medium with Cu and Zn, the recombinant protein is fully active, suggesting that the virus's SOD activation is independent of a copper chaperone for SOD generally used by eukaryotic SODs.

Kinetics of substrate inhibition of periplasmic nitrate reductase.

Periplasmic nitrate reductase catalyzes the reduction of nitrate into nitrite using a mononuclear molybdenum cofactor that has nearly the same structure in all enzymes of the DMSO reductase family. In previous electrochemical investigations, we found that the enzyme exists in several inactive states, some of which may have been previously isolated and mistaken for catalytic intermediates. In particular, the enzyme slowly and reversibly inactivates when exposed to high concentrations of nitrate. Here, we study the kinetics of substrate inhibition and its dependence on electrode potential and substrate concentration to learn about the properties of the active and inactive forms of the enzyme. We conclude that the substrate-inhibited enzyme never significantly accumulates in the EPR-active Mo(+V) state. This conclusion is relevant to spectroscopic investigations where attempts are made to trap a Mo(+V) catalytic intermediate using high concentrations of nitrate.

Reductive activation in periplasmic nitrate reductase involves chemical modifications of the Mo-cofactor beyond the first coordination sphere of the metal ion.

In Rhodobacter sphaeroides periplasmic nitrate reductase NapAB, the major Mo(V) form (the "high g" species) in air-purified samples is inactive and requires reduction to irreversibly convert into a catalytically competent form (Fourmond et al., J. Phys. Chem., 2008). In the present work, we study the kinetics of the activation process by combining EPR spectroscopy and direct electrochemistry. Upon reduction, the Mo (V) "high g" resting EPR signal slowly decays while the other redox centers of the protein are rapidly reduced, which we interpret as a slow and gated (or coupled) intramolecular electron transfer between the [4Fe-4S] center and the Mo cofactor in the inactive enzyme. Besides, we detect spin-spin interactions between the Mo(V) ion and the [4Fe-4S](1+) cluster which are modified upon activation of the enzyme, while the EPR signatures associated to the Mo cofactor remain almost unchanged. This shows that the activation process, which modifies the exchange coupling pathway between the Mo and the [4Fe-4S](1+) centers, occurs further away than in the first coordination sphere of the Mo ion. Relying on structural data and studies on Mo-pyranopterin and models, we propose a molecular mechanism of activation which involves the pyranopterin moiety of the molybdenum cofactor that is proximal to the [4Fe-4S] cluster. The mechanism implies both the cyclization of the pyran ring and the reduction of the oxidized pterin to give the competent tricyclic tetrahydropyranopterin form.

Resolution of key roles for the distal pocket histidine in cytochrome C nitrite reductases.

Cytochrome c nitrite reductases perform a key step in the biogeochemical N-cycle by catalyzing the six-electron reduction of nitrite to ammonium. These multiheme cytochromes contain a number of His/His ligated c-hemes for electron transfer and a structurally differentiated heme that provides the catalytic center. The catalytic heme has proximal ligation from lysine, or histidine, and an exchangeable distal ligand bound within a pocket that includes a conserved histidine. Here we describe properties of a penta-heme cytochrome c nitrite reductase in which the distal His has been substituted by Asn. The variant is unable to catalyze nitrite reduction despite retaining the ability to reduce a proposed intermediate in that process, namely, hydroxylamine. A combination of electrochemical, structural and spectroscopic studies reveals that the variant enzyme simultaneously binds nitrite and electrons at the catalytic heme. As a consequence the distal His is proposed to play a key role in orienting the nitrite for N-O bond cleavage. The electrochemical experiments also reveal that the distal His facilitates rapid nitrite binding to the catalytic heme of the native enzyme. Finally it is noted that the thermodynamic descriptions of nitrite- and electron-binding to the active site of the variant enzyme are modulated by the prevailing oxidation states of the His/His ligated hemes. This behavior is likely to be displayed by other multicentered redox enzymes such that there are wide implications for considering the determinants of catalytic activity in this important and varied group of oxidoreductases.

O2-independent formation of the inactive states of NiFe hydrogenase.

We studied the mechanism of aerobic inactivation of Desulfovibrio fructosovorans nickel-iron (NiFe) hydrogenase by quantitatively examining the results of electrochemistry, EPR and FTIR experiments. They suggest that, contrary to the commonly accepted mechanism, the attacking O(2) is not incorporated as an active site ligand but, rather, acts as an electron acceptor. Our findings offer new ways toward the understanding of O(2) inactivation and O(2) tolerance in NiFe hydrogenases.

Relating diffusion along the substrate tunnel and oxygen sensitivity in hydrogenase.

In hydrogenases and many other redox enzymes, the buried active site is connected to the solvent by a molecular channel whose structure may determine the enzyme's selectivity with respect to substrate and inhibitors. The role of these channels has been addressed using crystallography and molecular dynamics, but kinetic data are scarce. Using protein film voltammetry, we determined and then compared the rates of inhibition by CO and O2 in ten NiFe hydrogenase mutants and two FeFe hydrogenases. We found that the rate of inhibition by CO is a good proxy of the rate of diffusion of O2 toward the active site. Modifying amino acids whose side chains point inside the tunnel can slow this rate by orders of magnitude. We quantitatively define the relations between diffusion, the Michaelis constant for H2 and rates of inhibition, and we demonstrate that certain enzymes are slowly inactivated by O2 because access to the active site is slow.

DFT investigation of the molybdenum cofactor in periplasmic nitrate reductases: structure of the Mo(V) EPR-active species.

The periplasmic nitrate reductase NAP belongs to the DMSO reductase family that regroups molybdoenzymes housing a bis-molybdopterin cofactor as the active site. Several forms of the Mo(V) state, an intermediate redox state in the catalytic cycle of the enzyme, have been evidenced by EPR spectroscopy under various conditions, but their structure and catalytic relevance are not fully understood. On the basis of structural data available from the literature, we built several models that reproduce the first coordination sphere of the molybdenum cofactor and used DFT methods to make magneto-structural correlations on EPR-detected species. "High-g" states, which are the most abundant Mo(V) species, are characterized by a low-anisotropy g tensor and a high g(min) value. We assign this signature to a six-sulfur coordination sphere in a pseudotrigonal prismatic geometry with a partial disulfide bond. The "very high-g" species is well described with a sulfido ion as the sixth ligand. The "low-g" signal can be successfully associated to a Mo(V) sulfite-oxidase-type active site with only one pterin moiety coordinated to the molybdenum ion with an oxo or sulfido axial ligand. For all these species we investigate their catalytic activity using a thermodynamic point of view on the molybdenum coordination sphere. Beyond the periplasmic nitrate reductase case, this work provides useful magneto-structural correlations to characterize EPR-detected species in mononuclear molybdoenzymes.

Rates of intra- and intermolecular electron transfers in hydrogenase deduced from steady-state activity measurements.

Electrons are transferred over long distances along chains of FeS clusters in hydrogenases, mitochondrial complexes, and many other respiratory enzymes. It is usually presumed that electron transfer is fast in these systems, despite the fact that there has been no direct measurement of rates of FeS-to-FeS electron transfer in any respiratory enzyme. In this context, we propose and apply to NiFe hydrogenase an original strategy that consists of quantitatively interpreting the variations of steady-state activity that result from changing the nature of the FeS clusters which connect the active site to the redox partner, and/or the nature of the redox partner. Rates of intra- and intermolecular electron transfer are deduced from such large data sets. The mutation-induced variations of electron transfer rates cannot be explained by changes in intercenter distances and reduction potentials. This establishes that FeS-to-FeS rate constants are extremely sensitive to the nature and coordination of the centers.

Original design of an oxygen-tolerant [NiFe] hydrogenase: major effect of a valine-to-cysteine mutation near the active site.

Hydrogenases are efficient biological catalysts of H(2) oxidation and production. Most of them are inhibited by O(2), and a prerequisite for their use in biotechnological applications under air is to improve their oxygen tolerance. We have previously shown that exchanging the residue at position 74 in the large subunit of the oxygen-sensitive [NiFe] hydrogenase from Desulfovibrio fructosovorans could impact the reaction of the enzyme with O(2) (Dementin, S.; J. Am. Chem. Soc. 2009, 131, 10156-10164; Liebgott, P. P.; Nat. Chem. Biol. 2010, 6, 63-70). This residue, a valine in the wild-type enzyme, located at the bottleneck of the gas channel near the active site, has here been exchanged with a cysteine. A thorough characterization using a combination of kinetic, spectroscopic (EPR, FTIR), and electrochemical studies demonstrates that the V74C mutant has features of the naturally occurring oxygen-tolerant membrane-bound hydrogenases (MBH). The mutant is functional during several minutes under O(2), has impaired H(2)-production activity, and has a weaker affinity for CO than the WT. Upon exposure to O(2), it is converted into the more easily reactivatable inactive form, Ni-B, and this inactive state reactivates about 20 times faster than in the WT enzyme. Control experiments carried out with the V74S and V74N mutants indicate that protonation of the position 74 residue is not the reason the mutants reactivate faster than the WT enzyme. The electrochemical behavior of the V74C mutant toward O(2) is intermediate between that of the WT enzyme from D. fructosovorans and the oxygen-tolerant MBH from Aquifex aeolicus.

Experimental approaches to kinetics of gas diffusion in hydrogenase.

Hydrogenases, which catalyze H(2) to H(+) conversion as part of the bioenergetic metabolism of many microorganisms, are among the metalloenzymes for which a gas-substrate tunnel has been described by using crystallography and molecular dynamics. However, the correlation between protein structure and gas-diffusion kinetics is unexplored. Here, we introduce two quantitative methods for probing the rates of diffusion within hydrogenases. One uses protein film voltammetry to resolve the kinetics of binding and release of the competitive inhibitor CO; the other is based on interpreting the yield in the isotope exchange assay. We study structurally characterized mutants of a NiFe hydrogenase, and we show that two mutations, which significantly narrow the tunnel near the entrance of the catalytic center, decrease the rates of diffusion of CO and H(2) toward and from the active site by up to 2 orders of magnitude. This proves the existence of a functional channel, which matches the hydrophobic cavity found in the crystal. However, the changes in diffusion rates do not fully correlate with the obstruction induced by the mutation and deduced from the x-ray structures. Our results demonstrate the necessity of measuring diffusion rates and emphasize the role of side-chain dynamics in determining these.

Dependence of catalytic activity on driving force in solution assays and protein film voltammetry: insights from the comparison of nitrate reductase mutants.

Rhodobacter sphaeroides periplasmic nitrate reductase (Rs NapAB) is one of the enzymes whose assays give odd results: in spectrophotometric assays with methyl viologen as the electron donor, the activity increases as the reaction progresses, whereas the driving force provided by the soluble redox partner decreases; in protein film voltammetry (PFV), whereby the enzyme directly exchanges electrons with an electrode, the activity of NapAB decreases at large overpotential, whereas a monotonic increase is expected [Elliott, S. J., et al. (2002) Biochim. Biophys. Acta 1555, 54-59]. The relations between these phenomena and the catalytic mechanism are still debated. By studying NapAB mutants, we found that the peculiar dependences of electrochemical and solution activities on driving force are greatly affected by substituting certain amino acids that are located in the vicinity of the active site (M153, Q384, R392); this led us to establish and discuss the relation between the experimental parameters of the electrochemical and spectrophotometric assays: we show that the rate of reduction of the enzyme (which depends on the electrode potential or on the concentration of reduced MV) modulates the activity of the enzyme, but the "solution potential" does not. Our results also support the view that the complex profiles of activity versus potential are fingerprints of the active site chemistry, rather than direct consequences of changes in the redox states of relays that are remote from the active site.

Introduction of methionines in the gas channel makes [NiFe] hydrogenase aero-tolerant.

Hydrogenases catalyze the conversion between 2H(+) + 2e(-) and H(2)(1). Most of these enzymes are inhibited by O(2), which represents a major drawback for their use in biotechnological applications. Improving hydrogenase O(2) tolerance is therefore a major contemporary challenge to allow the implementation of a sustainable hydrogen economy. We succeeded in improving O(2) tolerance, which we define here as the ability of the enzyme to resist for several minutes to O(2) exposure, by substituting with methionines small hydrophobic residues strongly conserved in the gas channel. Remarkably, the mutated enzymes remained active in the presence of an O(2) concentration close to that found in aerobic solutions in equilibrium with air, while the wild type enzyme is inhibited in a few seconds. Crystallographic and spectroscopic studies showed that the structure and the chemistry at the active site are not affected by the mutations. Kinetic studies demonstrated that the inactivation is slower and reactivation faster in these mutants. We propose that in addition to restricting O(2) diffusion to the active site of the enzyme, methionine may also interact with bound peroxide and provide an assisted escape route for H(2)O(2) toward the gas channel. These results show for the first time that it is possible to improve O(2)-tolerance of [NiFe] hydrogenases, making possible the development of biohydrogen production systems.

Tuning the redox properties of a [4Fe-4S] center to modulate the activity of Mo-bisPGD periplasmic nitrate reductase.

Molybdoenzymes are ubiquitous in living organisms and catalyze, for most of them, oxidation-reduction reactions using a large range of substrates. Periplasmic nitrate reductase (NapAB) from Rhodobacter sphaeroides catalyzes the 2-electron reduction of nitrate into nitrite. Its active site is a Mo bis-(pyranopterin guanine dinucleotide), or Mo-bisPGD, found in most prokaryotic molybdoenzymes. A [4Fe-4S] cluster and two c-type hemes form an intramolecular electron transfer chain that deliver electrons to the active site. Lysine 56 is a highly conserved amino acid which connects, through hydrogen-bonds, the [4Fe-4S] center to one of the pyranopterin ligands of the Mo-cofactor. This residue was proposed to be involved in the intramolecular electron transfer, either defining an electron transfer pathway between the two redox cofactors, and/or modulating their redox properties. In this work, we investigated the role of this lysine by combining site-directed mutagenesis, activity assays, redox titrations, EPR and HYSCORE spectroscopies. Removal of a positively-charged residue at position 56 strongly decreased the redox potential of the [4Fe-4S] cluster at pH 8 by 230 mV to 400 mV in the K56H and K56M mutants, respectively, thus affecting the kinetics of electron transfer from the hemes to the [4Fe-4S] center up to 5 orders of magnitude. This effect was partly reversed at acidic pH in the K56H mutant likely due to protonation of the imidazole ring of the histidine. Overall, our study demonstrates the critical role of a charged residue from the second coordination sphere in tuning the reduction potential of the [4Fe-4S] cluster in RsNapAB and related molybdoenzymes.

Major Mo(V) EPR signature of Rhodobacter sphaeroides periplasmic nitrate reductase arising from a dead-end species that activates upon reduction. Relation to other molybdoenzymes from the DMSO reductase family.

Enzymes of the DMSO reductase family use a mononuclear Mo-bis(molybdopterin) cofactor (MoCo) to catalyze a variety of oxo-transfer reactions. Much functional information on nitrate reductase, one of the most studied members of this family, has been gained from EPR spectroscopy, but this technique is not always conclusive because the signature of the MoCo is heterogeneous, and which signals correspond to active species is still unsure. We used site-directed mutagenesis, EPR and protein film voltammetry to demonstrate that the MoCo in R. sphaeroides periplasmic nitrate reductase (NapAB) is subject to an irreversible reductive activation process whose kinetics we precisely define. This activation quantitatively correlates with the disappearance of the so-called "Mo(V) high-g" EPR signal, but this reductive process is too slow to be part of the normal catalytic cycle. Therefore, in NapAB, this most intense and most commonly observed signature of the MoCo arises from a dead-end, inactive state that gives a catalytically competent species only after reduction. This activation proceeds, even without substrate, according to a reduction followed by an irreversible nonredox step, both of which are pH independent. An apparently similar process occurs in other nitrate reductases (both assimilatory and membrane bound), and this also recalls the redox cycling procedure, which activates periplasmic DMSO reductases and simplifies their spectroscopic signatures. Hence we propose that heterogeneity at the active site and reductive activation are common properties of enzymes from the DMSO reductase family. Regarding NapAB, the fact that we could detect no Mo EPR signal upon reoxidizing the fully reduced enzyme suggests that the catalytically active form of the Mo(V) is thermodynamically unstable, as is the case for other enzymes of the DMSO reductase family. Our original approach, which combines spectroscopy and protein film voltammetry, proves useful for discriminating the forms of the active site on the basis of their catalytic properties. This could be applied to other enzymes for which the question arises as to the catalytic relevance of certain long-lived, spectroscopically characterized species.