Phylogeographic analyses of an endemic Neotropical fox (Lycalopex vetulus) reveal evidence of hybridization with a different canid species (L. gymnocercus)

The hoary fox (Lycalopex vetulus) is the only species of the Canidae (Mammalia: Carnivora) endemic to Brazil, and so far has been the target of few genetic studies. Using microsatellites and mtDNA markers, we investigated its present genetic diversity and population structure. We also tested the hypothesis that this species currently hybridizes with the pampas fox (L. gymnocercus), as suggested by previous mtDNA data from two individuals. We collected tissue and blood samples from animals representing most of the two species' distributions in Brazil (n = 87), including their recently discovered geographic contact zone in São Paulo state. We observed that the hoary fox exhibits high levels of genetic diversity and low levels of population structure. We identified six individuals from São Paulo state with clear evidence of hybridization based on introgressed pampas fox mtDNA and/or admixed microsatellite genotypes (three individuals bore both types of evidence). These results demonstrate the existence of admixed individuals between hoary and pampas foxes in southeastern Brazil, representing the first identified case of inter-species admixture between native South American canids. We discuss our findings in the context of the evolutionary history of these foxes and address potential conservation implications of this interspecies hybridization process.

Brain volumetry from CT-scan endocasts of three neotropical carnivores

Studies on brain anatomy can clarify specific evolutionary and behavioural aspects of wild animals. The rich diversity in a broad range of habitats makes carnivorans especially eligible for studying the relations between the brain form and behaviour, cognitive, sensorial and motor traits. This study compared the brain's contour and total and segmented brain volumetry in three species of neotropical carnivorans. CT images of 17 skulls of three species were acquired: Conepatus chinga (n = 6), Galictis cuja (n = 6) and Lontra longicaudis (n = 5). Three-dimensional endocasts allowed for estimating the brain's total and segmented volumes (olfactory bulb, rostral cerebrum, caudal cerebrum and cerebellum/brain stem). The average volume percentage of the segments was compared interspecifically and intraspecifically between the sexes. The otter has a notably more complex gyrification, typical for semiaquatic carnivorans. Proportionally, the olfactory bulb was significantly larger in hog-nosed skunks, possibly due to a better sense of smell for capturing insects. The proportional volumes of the rostral cerebrum, caudal cerebrum and cerebellum/brain stem segments did not differ between these species. Social behaviour traits and tactile, motor and balance skills were probably not sufficiently distinct to reflect differences in the brain segments analysed in these three species.

Coronary anatomy in neotropical carnivores: A comparative analysis

This study aims to investigate the anatomy of coronary arteries and their primary branches in neotropical carnivores, given the functional significance of coronary circulation and the limited existing descriptions in this order. For this purpose, coronary arteries were dissected, and their branches were counted in 74 hearts from 20 species belonging to five families: Canidae, Mustelidae, Procyonidae, Mephitidae, and Felidae. In all specimens of the suborder Caniformia, the subsinuosal interventricular branch originated at the end of the circumflex branch of the left coronary artery (LCA). In contrast, in all hearts of Felidae specimens, the subsinuosal branch was derived from the right coronary artery. Among the 50 neotropical specimens of the suborder Caniformia, 13.1 ± 5.2 ventricular branches originated from the LCA and only 4.6 ± 1.4 from the right (p < 0.0001), characterizing a consistent left dominance. However, in the 24 specimens of the Felidae family, 10.5 ± 4.0 ventricular branches were found originating from the LCA and 10.4 ± 3.4 from the right, with no difference between the means (p = 0.82), defining a balanced pattern. It is suggested that the type of coronary circulation may correlate with the phylogeny of carnivorans.

NEOTROPICAL FRESHWATER FISHES: A dataset of occurrence and abundance of freshwater fishes in the Neotropics

The Neotropical region hosts 4225 freshwater fish species, ranking first among the world's most diverse regions for freshwater fishes. Our NEOTROPICAL FRESHWATER FISHES data set is the first to produce a large-scale Neotropical freshwater fish inventory, covering the entire Neotropical region from Mexico and the Caribbean in the north to the southern limits in Argentina, Paraguay, Chile, and Uruguay. We compiled 185,787 distribution records, with unique georeferenced coordinates, for the 4225 species, represented by occurrence and abundance data. The number of species for the most numerous orders are as follows: Characiformes (1289), Siluriformes (1384), Cichliformes (354), Cyprinodontiformes (245), and Gymnotiformes (135). The most recorded species was the characid Astyanax fasciatus (4696 records). We registered 116,802 distribution records for native species, compared to 1802 distribution records for nonnative species. The main aim of the NEOTROPICAL FRESHWATER FISHES data set was to make these occurrence and abundance data accessible for international researchers to develop ecological and macroecological studies, from local to regional scales, with focal fish species, families, or orders. We anticipate that the NEOTROPICAL FRESHWATER FISHES data set will be valuable for studies on a wide range of ecological processes, such as trophic cascades, fishery pressure, the effects of habitat loss and fragmentation, and the impacts of species invasion and climate change. There are no copyright restrictions on the data, and please cite this data paper when using the data in publications.

Ticks, mites, fleas, and vector-borne pathogens in free-ranging neotropical wild felids from southern Brazil

The study of vectors and vector-borne diseases associated with ectoparasites from free-living Neotropical small wild felids is scarce, and the few existing studies on this theme focused on the genus Panthera, Paleartic species or in captive animals. For this reason, the aim of this study was to identify the diversity of ectoparasites and potentially associated Rickettsia and Bartonella species in free-ranging neotropical wild cats collected in the State of Rio Grande do Sul, southern Brazil. A total of 82 ticks, 10 fleas and 22 mites were collected from 18 Geoffroy's cats (Leopardus geoffroyi), 11 margays (Leopardus wiedii), two southern tiger cats (Leopardus guttulus), two jaguarundis (Herpailurus yagouaroundi), three ocelots (Leopardus pardalis) and two pumas (Puma concolor). We identified four tick species: Rhipicephalus microplus and three species from genus Amblyomma, the most frequent being Amblyomma aureolatum; three flea species corresponding to Ctenocephalides felis, Xenopsylla cheopis, and Pulex irritans; and one mite of genus Eutrombicula. In ectoparasites we found DNA of Rickettsia parkeri and Rickettsia asembonensis in ticks and DNA of Bartonella clarridgeiae, Bartonella sp. and R. asembonensis in fleas. Our results highlight the evidences of vectors and vector-borne agents in wild animals, and their potential wide distribution in the Pampa biome and the southernmost portion of the Atlantic Rainforest in Brazil.

Geographic variation in the skull morphology of the lesser grison (Galictis cuja: Carnivora, Mustelidae) from two Brazilian ecoregions

Background

The lesser grison (Galictis cuja) is one of the least known carnivores in the Neotropical region. Its wide geographical occurrence and range of habitats could lead to morphological variations along its distribution. So, this study aimed to investigate the variation in skull shape and size of this species, by testing the existence of ecotypes adapted to their respective environments (Uruguayan savanna and Atlantic Forest), as well as its relationship with selected abiotic variables.

Methods

The skulls of 52 museum specimens were photographed in the ventral, dorsal, and lateral views, and were analyzed using geometric morphometric techniques.

Results

We found sexual size dimorphism, with males being larger than females. The shape variation between sexes, as well as between ecoregions, is mostly explained by the effect of allometry. The specimens from Uruguayan savanna are larger than the ones from the Atlantic Forest. Size variation was also significantly correlated to latitude, temperature and precipitation patterns. No correlation between skull shape with geographical distance was detected.

Discussion

Morphometric measurements and diet data of lesser grison in regions from higher latitudes than our sampling show a tendency to heavier individuals, and the consumption of bigger prey compared to Uruguayan savanna. The results indicated the smaller specimens associated to low variability in annual temperature, congruent to Atlantic Forest region. An explanation for observed variation may be related to the “resource rule” but, due the minimal natural history information regards this species, we can just speculate about this.

NEOTROPICAL CARNIVORES: a data set on carnivore distribution in the Neotropics

Mammalian carnivores are considered a key group in maintaining ecological health and can indicate potential ecological integrity in landscapes where they occur. Carnivores also hold high conservation value and their habitat requirements can guide management and conservation plans. The order Carnivora has 84 species from 8 families in the Neotropical region: Canidae; Felidae; Mephitidae; Mustelidae; Otariidae; Phocidae; Procyonidae; and Ursidae. Herein, we include published and unpublished data on native terrestrial Neotropical carnivores (Canidae; Felidae; Mephitidae; Mustelidae; Procyonidae; and Ursidae). NEOTROPICAL CARNIVORES is a publicly available data set that includes 99,605 data entries from 35,511 unique georeferenced coordinates. Detection/non-detection and quantitative data were obtained from 1818 to 2018 by researchers, governmental agencies, non-governmental organizations, and private consultants. Data were collected using several methods including camera trapping, museum collections, roadkill, line transect, and opportunistic records. Literature (peer-reviewed and grey literature) from Portuguese, Spanish and English were incorporated in this compilation. Most of the data set consists of detection data entries (n = 79,343; 79.7%) but also includes non-detection data (n = 20,262; 20.3%). Of those, 43.3% also include count data (n = 43,151). The information available in NEOTROPICAL CARNIVORES will contribute to macroecological, ecological, and conservation questions in multiple spatio-temporal perspectives. As carnivores play key roles in trophic interactions, a better understanding of their distribution and habitat requirements are essential to establish conservation management plans and safeguard the future ecological health of Neotropical ecosystems. Our data paper, combined with other large-scale data sets, has great potential to clarify species distribution and related ecological processes within the Neotropics. There are no copyright restrictions and no restriction for using data from this data paper, as long as the data paper is cited as the source of the information used. We also request that users inform us of how they intend to use the data.

NEOTROPICAL ALIEN MAMMALS: a data set of occurrence and abundance of alien mammals in the Neotropics

Biological invasion is one of the main threats to native biodiversity. For a species to become invasive, it must be voluntarily or involuntarily introduced by humans into a nonnative habitat. Mammals were among first taxa to be introduced worldwide for game, meat, and labor, yet the number of species introduced in the Neotropics remains unknown. In this data set, we make available occurrence and abundance data on mammal species that (1) transposed a geographical barrier and (2) were voluntarily or involuntarily introduced by humans into the Neotropics. Our data set is composed of 73,738 historical and current georeferenced records on alien mammal species of which around 96% correspond to occurrence data on 77 species belonging to eight orders and 26 families. Data cover 26 continental countries in the Neotropics, ranging from Mexico and its frontier regions (southern Florida and coastal-central Florida in the southeast United States) to Argentina, Paraguay, Chile, and Uruguay, and the 13 countries of Caribbean islands. Our data set also includes neotropical species (e.g., Callithrix sp., Myocastor coypus, Nasua nasua) considered alien in particular areas of Neotropics. The most numerous species in terms of records are from Bos sp. (n = 37,782), Sus scrofa (n = 6,730), and Canis familiaris (n = 10,084); 17 species were represented by only one record (e.g., Syncerus caffer, Cervus timorensis, Cervus unicolor, Canis latrans). Primates have the highest number of species in the data set (n = 20 species), partly because of uncertainties regarding taxonomic identification of the genera Callithrix, which includes the species Callithrix aurita, Callithrix flaviceps, Callithrix geoffroyi, Callithrix jacchus, Callithrix kuhlii, Callithrix penicillata, and their hybrids. This unique data set will be a valuable source of information on invasion risk assessments, biodiversity redistribution and conservation-related research. There are no copyright restrictions. Please cite this data paper when using the data in publications. We also request that researchers and teachers inform us on how they are using the data.

Anatomical variations in cervical vertebrae in two species of neotropical canids: What is the meaning?

Anatomical variation in wild canids is rarely described. Altered shapes, uni- or bilaterally situated, were observed in ventral laminae of C6 and/or in C7 vertebral body of two fox species: Lycalopex gymnocercus and Cerdocyon thous. The specimens were young adults, collected dead on highways in southern Brazil. Deformities were considered anatomical variations because apparently there was no interference in function. The occurrence of such variations in two related species from the same area can suggest an environmental common cause rather than genetic issues. Future studies may deepen the understanding between cause and effect of these vertebral deformities in wild canids.

Who's afraid of the big bad boar? Assessing the effect of wild boar presence on the occurrence and activity patterns of other mammals

Wild boar are considered one the world’s worst invasive species and linked to biodiversity loss, competition for resources, predation of native species, and habitat modifications. In this study, we use camera traps to evaluate whether the invasive wild boar had an effect on the medium-sized mammal community of a protected area in southern Brazil. Based on photographic records, we evaluated whether the presence and relative abundance of wild boar was associated with a decrease in diversity or change in activity of medium-sized mammals. All comparisons were made between samples where wild boar were present or absent. The records of each camera during a season were considered a sample. The wild boar was the fourth most common species in the study area being present in 7.8% of the photographic records. The species richness of mammals was not negatively affected by the occurrence of wild boar and most common species did not exhibit changes in the daily activity pattern. However, we recorded an increase in the time elapsed between an observation of wild boar and the record of the next species relative to the average latency period observed among other mammalian species. This average latency period was similar to that observed in the case of large predators such as Puma, and its increase could be reflective partly of the avoidance of native species to wild boar. Nevertheless, our results show that the effect of invasive wild boar on the mammal community is not large as expected.

NEOTROPICAL XENARTHRANS: a data set of occurrence of xenarthran species in the Neotropics

Xenarthrans-anteaters, sloths, and armadillos-have essential functions for ecosystem maintenance, such as insect control and nutrient cycling, playing key roles as ecosystem engineers. Because of habitat loss and fragmentation, hunting pressure, and conflicts with domestic dogs, these species have been threatened locally, regionally, or even across their full distribution ranges. The Neotropics harbor 21 species of armadillos, 10 anteaters, and 6 sloths. Our data set includes the families Chlamyphoridae (13), Dasypodidae (7), Myrmecophagidae (3), Bradypodidae (4), and Megalonychidae (2). We have no occurrence data on Dasypus pilosus (Dasypodidae). Regarding Cyclopedidae, until recently, only one species was recognized, but new genetic studies have revealed that the group is represented by seven species. In this data paper, we compiled a total of 42,528 records of 31 species, represented by occurrence and quantitative data, totaling 24,847 unique georeferenced records. The geographic range is from the southern United States, Mexico, and Caribbean countries at the northern portion of the Neotropics, to the austral distribution in Argentina, Paraguay, Chile, and Uruguay. Regarding anteaters, Myrmecophaga tridactyla has the most records (n = 5,941), and Cyclopes sp. have the fewest (n = 240). The armadillo species with the most data is Dasypus novemcinctus (n = 11,588), and the fewest data are recorded for Calyptophractus retusus (n = 33). With regard to sloth species, Bradypus variegatus has the most records (n = 962), and Bradypus pygmaeus has the fewest (n = 12). Our main objective with Neotropical Xenarthrans is to make occurrence and quantitative data available to facilitate more ecological research, particularly if we integrate the xenarthran data with other data sets of Neotropical Series that will become available very soon (i.e., Neotropical Carnivores, Neotropical Invasive Mammals, and Neotropical Hunters and Dogs). Therefore, studies on trophic cascades, hunting pressure, habitat loss, fragmentation effects, species invasion, and climate change effects will be possible with the Neotropical Xenarthrans data set. Please cite this data paper when using its data in publications. We also request that researchers and teachers inform us of how they are using these data.

Functional myology of the thoracic limb in Pampas fox (Lycalopex gymnocercus): a descriptive and comparative analysis

The characteristics of the muscles of the thoracic limb were evaluated in 22 specimens of Lycalopex gymnocercus. Descriptive and comparative analyses showed similarity with other canids in terms of topography and tendon insertions. Differences with the domestic dog were observed in the pectoralis profundus, triceps brachii and interflexorii muscles. Intraspecific variations were observed in the rhomboideus capitis, serratus ventralis cervicis, extensor carpi radialis, extensor digiti I and II, lumbricales, flexor digiti I brevis, abductor digiti I brevis, and flexor digiti V muscles. The analyses of muscle architecture carried out in nine specimens showed that there was no difference in muscle percentage mass in the thoracic limb of males and females, but a young specimen showed significant lower percentage mass. The triceps brachii caput longus muscle showed the greatest mass, the subscapularis muscle showed the greatest physiological cross-sectional area value, and the extrinsic muscles, in general, presented the longest fascicles and higher architectural indexes. Muscle architecture data were compatible with those of a thoracic limb adapted to fast cursorial locomotion that prioritizes movements in a sagittal plane instead of rotation or adduction/abduction. There was a high association between functional percentage mass of the muscles in the thoracic limb and phylogeny in the Carnivora order. It may be inferred that carnivoran muscle mass is largely determined by phylogeny.

BRAZIL ROAD-KILL: a data set of wildlife terrestrial vertebrate road-kills

Mortality from collision with vehicles is the most visible impact of road traffic on wildlife. Mortality due to roads (hereafter road-kill) can affect the dynamic of populations of many species and can, therefore, increase the risk of local decline or extinction. This is especially true in Brazil, where plans for road network upgrading and expansion overlaps biodiversity hotspot areas, which are of high importance for global conservation. Researchers, conservationists and road planners face the challenge to define a national strategy for road mitigation and wildlife conservation. The main goal of this dataset is a compilation of geo-referenced road-kill data from published and unpublished road surveys. This is the first Data Paper in the BRAZIL series (see ATLANTIC, NEOTROPICAL, and BRAZIL collections of Data Papers published in Ecology), which aims make public road-kill data for species in the Brazilian Regions. The dataset encompasses road-kill records from 45 personal communications and 26 studies published in peer-reviewed journals, theses and reports. The road-kill dataset comprises 21,512 records, 83% of which are identified to the species level (n = 450 species). The dataset includes records of 31 amphibian species, 90 reptile species, 229 bird species, and 99 mammal species. One species is classified as Endangered, eight as Vulnerable and twelve as Near Threatened. The species with the highest number of records are: Didelphis albiventris (n = 1,549), Volatinia jacarina (n = 1,238), Cerdocyon thous (n = 1,135), Helicops infrataeniatus (n = 802), and Rhinella icterica (n = 692). Most of the records came from southern Brazil. However, observations of the road-kill incidence for non-Least Concern species are more spread across the country. This dataset can be used to identify which taxa seems to be vulnerable to traffic, analyze temporal and spatial patterns of road-kill at local, regional and national scales and also used to understand the effects of road-kill on population persistence. It may also contribute to studies that aims to understand the influence of landscape and environmental influences on road-kills, improve our knowledge on road-related strategies on biodiversity conservation and be used as complementary information on large-scale and macroecological studies. No copyright or proprietary restrictions are associated with the use of this data set other than citation of this Data Paper.

Trophic ecology of sympatric small cats in the Brazilian Pampa

Information about resource partitioning among small cat species that live in sympatry in South America is fairly incomplete. Knowledge about feeding habits is essential for understanding the role of these predators in the environment, the impact on prey populations, and potential competition among themselves and with other carnivores. This study aimed to describe and compare the diet of four sympatric small cats in the grasslands of southern Brazil. We analysed the stomach contents of 37 Geoffroy’s cats (Leopardus geoffroyi), 27 margays (Leopardus wiedii), 14 pampas cats (Leopardus colocola), and 20 jaguarundis (Herpailurus yagouaroundi) obtained as road kill in the Brazilian Pampa in southern Brazil. Small mammals were the most representative class consumed by all cats, followed by Aves, Reptilia, and Amphibia. Some items, such as rodents Cavia aperea, Akodon sp., Oligoryzomys sp. and Passeriformes were consumed by all cat species. Niche overlap varied widely, from 10% (margay x jaguarundi) to 92% (jaguarundi x pampas cat). Niche breadth indicated that jaguarundi were the most specialized of the cats (B_sta = 0.24) in this region, with a diet closely associated to C. aperea. Margay consumed more items associated with arboreal behaviour than other cat species, but consumed more terrestrial items than arboreal ones. The pampas cat consumed mostly terrestrial species associated with open fields. Geoffroy’s cat consumed mammals found in a diversity of habitats, indicating high ecological flexibility. Species with more similarity in diet such as jaguarundi and pampas cat probably present temporal segregation in activity. In conclusion, despite their habitat and diet similarities, these four species explore distinct microhabitats by foraging different prey groups, what favor them to live in sympatry.

Activity pattern of medium and large sized mammals and density estimates of Cuniculus paca (Rodentia: Cuniculidae) in the Brazilian Pampa

Between July 2014 and April 2015, we conducted weekly inventories of the circadian activity patterns of mammals in Passo Novo locality, municipality of Alegrete, southern Brazil. The vegetation is comprised by a grassy-woody steppe (grassland). We used two camera traps alternately located on one of four 1 km transects, each separated by 1 km. We classified the activity pattern of species by the percentage of photographic records taken in each daily period. We identify Cuniculus paca individuals by differences in the patterns of flank spots. We then estimate the density 1) considering the area of riparian forest present in the sampling area, and 2) through capture/recapture analysis. Cuniculus paca, Conepatus chinga and Hydrochoerus hydrochaeris were nocturnal, Cerdocyon thous had a crepuscular/nocturnal pattern, while Mazama gouazoubira was cathemeral. The patterns of circadian activity observed for medium and large mammals in this Pampa region (southern grasslands) may reflect not only evolutionary, biological and ecological affects, but also human impacts not assessed in this study. We identified ten individuals of C. paca through skin spot patterns during the study period, which were recorded in different transects and months. The minimum population density of C. paca was 3.5 individuals per km2 (resident animals only) and the total density estimates varied from 7.1 to 11.8 individuals per km2, when considering all individuals recorded or the result of the capture/recapture analysis, respectively.

Home range and density of three sympatric felids in the Southern Atlantic Forest, Brazil

Home range and minimal population densities of Southern tiger cat (Leopardus guttulus), margay (Lepardus wiedii) and jaguarundi (Puma yagouaroundi) were estimated between 2005 and 2006 in Taquari Valley, near the southern edge of the Atlantic Rainforest in Brazil. Home range data were collected by conventional radio telemetry (VHF) locations in a highly fragmented landscape. The average home range size, calculated using 95% kernel density estimates, was 16.01 km2 for Southern tiger cat, 21.85 km2 for margay and 51.45 km2 for jaguarundi. Telemetry data were used to obtain minimal density estimates of 0.08 Southern tiger cats / km2, and 0.04 jaguarundi / km2. The density estimates arise from areas where ocelot (Leopardus pardalis) and other larger-bodied carnivores were locally extinct, and they suggest a specific type of mesopredator release known as the ocelot effect, which is likely enabling the increase in smaller felid populations in this area.

Trophic relationships of sympatric small carnivores in fragmented landscapes of southern Brazil: niche overlap and potential for competition

Between 2000 and 2010, digestive tracts collected from carnivore carcasses found in southern Brazil were analyzed to determine the frequency and proportion of items constituting the diets of each species. Material was collected and analyzed from 194 animals of 10 species:

Trophic ecology and the use of shelters and latrines by the Neotropical otter (Lontra longicaudis) in the Taquari Valley, Southern Brazil

This manuscript presents information about the ecology of Lontra longicaudis (Olfers, 1818) in the Taquari Valley, State of Rio Grande do Sul, southern Brazil. The study was carried out in two areas located in the Forquetinha Creek and in the Forqueta River from January to December 2003. The otters are specialist feeders (Bsta = 0.24), with a diet based mostly on fish, especially those of the families Loricariidae and Cichlidae. Most shelters used by the species were excavated burrows underneath tree roots, while shelters within rocks were used less frequently. The burrows showed great variation in size, being found on average 3.5 m (sd = 3.6 m) away from the margin and 2.5 m (sd = 1.2 m) above the water level. Scent marks were made preferentially on rocks and fallen tree trunks at the edge of the water. There was a tendency to increase the reutilization of latrines in detriment of using new sites throughout the sample period.

NADPH-cytochrome P450 oxidoreductase. Structural basis for hydride and electron transfer

NADPH-cytochrome P450 oxidoreductase catalyzes transfer of electrons from NADPH, via two flavin cofactors, to various cytochrome P450s. The crystal structure of the rat reductase complexed with NADP(+) has revealed that nicotinamide access to FAD is blocked by an aromatic residue (Trp-677), which stacks against the re-face of the isoalloxazine ring of the flavin. To investigate the nature of interactions between the nicotinamide, FAD, and Trp-677 during the catalytic cycle, three mutant proteins were studied by crystallography. The first mutant, W677X, has the last two C-terminal residues, Trp-677 and Ser-678, removed; the second mutant, W677G, retains the C-terminal serine residue. The third mutant has the following three catalytic residues substituted: S457A, C630A, and D675N. In the W677X and W677G structures, the nicotinamide moiety of NADP(+) lies against the FAD isoalloxazine ring with a tilt of approximately 30 degrees between the planes of the two rings. These results, together with the S457A/C630A/D675N structure, allow us to propose a mechanism for hydride transfer regulated by changes in hydrogen bonding and pi-pi interactions between the isoalloxazine ring and either the nicotinamide ring or Trp-677 indole ring. Superimposition of the mutant and wild-type structures shows significant mobility between the two flavin domains of the enzyme. This, together with the high degree of disorder observed in the FMN domain of all three mutant structures, suggests that conformational changes occur during catalysis.

Differential contributions of NADPH-cytochrome P450 oxidoreductase FAD binding site residues to flavin binding and catalysis

Transfer of reducing equivalents from NADPH to the cytochromes P450 is mediated by NADPH-cytochrome P450 oxidoreductase, which contains stoichiometric amounts of tightly bound FMN and FAD. Hydrogen bonding and van der Waals interactions between FAD and amino acid residues in the FAD binding site of the reductase serve to regulate both flavin binding and reactivity. The precise orientation of key residues (Arg(454), Tyr(456), Cys(472), Gly(488), Thr(491), and Trp(677)) has been defined by x-ray crystallography (Wang, M., Roberts, D. L., Paschke, R., Shea, T. M., Masters, B. S., Kim, J.-J. P. (1997) Proc. Natl. Acad. Sci. U. S. A. 94, 8411-8416). The current study examines the relative contributions of these residues to FAD binding and catalysis by site-directed mutagenesis and kinetic analysis. Mutation of either Tyr(456), which makes van der Waals contact with the FAD isoalloxazine ring and also hydrogen-bonds to the ribityl 4'-hydroxyl, or Arg(454), which bonds to the FAD pyrophosphate, decreases the affinity for FAD 8000- and 25,000-fold, respectively, with corresponding decreases in cytochrome c reductase activity. In contrast, substitution of Thr(491), which also interacts with the pyrophosphate grouping, had a relatively modest effect on both FAD binding (100-fold decrease) and catalytic activity (2-fold decrease), while the G488L mutant exhibited, respectively, 800- and 50-fold decreases in FAD binding and catalytic activity. Enzymic activity of each of these mutants could be restored by addition of FAD. Kinetic properties and the FMN content of these mutants were not affected by these substitutions, with the exception of a 3-fold increase in Y456S K(m)(cyt )(c) and a 70% decrease in R454E FMN content, suggesting that the FMN- and FAD-binding domains are largely, but not completely, independent. Even though Trp(677) is stacked against the re-face of FAD, suggesting an important role in FAD binding, deletion of both Trp(677) and the carboxyl-terminal Ser(678) decreased catalytic activity 50-fold without affecting FAD content.

Molecular basis for cell-specific regulation of the NADPH-cytochrome P450 oxidoreductase gene

NADPH-cytochrome P450 oxidoreductase (CYPOR), a flavoprotein localized in the nuclear envelope and endoplasmic reticulum of most cell types, is responsible for transferring electrons from NADPH to the cytochromes P450 as well as heme oxygenase, squalene epoxidase, and cytochrome b(5). CYPOR is encoded by a single gene and, similar to many housekeeping genes, has a TATA-less, GC-rich promoter with multiple Sp1 consensus sites. The current work has delineated the importance of multiple cis-acting elements contained within the proximal promoter for basal expression of the CYPOR gene. Transcription factor binding sites within this region included two upstream Sp1 motifs, a SEC element containing overlapping Sp1/Egr-1/CACCC box motifs, and a novel site designated the OxidoReductase Upstream element (ORU). Mutational modification of the ORU element, leading to a loss of protein binding, resulted in an approximately 90% decrease in transcriptional activity in H4IIE cells. Similarly, inactivation of the Egr-1/CACCC segment of the SEC element dramatically reduced promoter activity to less than 10% of wild-type, while mutagenesis of the contiguous Sp1 site did not affect basal transcription. Although both Sp1 sites contained within the minimal promoter were required for optimal expression in H4IIE cells, loss of these sites was compensated for by those Sp1 motifs located upstream of position 206, suggesting that Sp1 was acting as a position-independent enhancer. Hence, the CYPOR promoter was distinguished from the majority of TATA-less promoters in that Sp1 was not a primary transcriptional regulator and by the fact that the Sp1 binding site closest to the transcription start site was nonfunctional. Furthermore, both the SEC and ORU elements were essential for basal expression; however, the ORU element exhibited cell-specific differences in regulatory activity. Thus, several mechanisms appear to be in place to selectively alter the expression of the CYPOR gene.

Two-stage glucocorticoid induction of CYP3A23 through both the glucocorticoid and pregnane X receptors

Glucocorticoid inducibility of the CYP3A23 gene is conferred by a multisite unit comprising binding sites for several members of the nuclear receptor superfamily of transcription factors, including the chicken ovalbumin upstream promoter-transcription factor COUP-TF, pregnane X receptor (PXR), and hepatocyte nuclear factor 4 (HNF-4). The presence of three binding sites, each of which interacts with more than one factor, contributes to the complexity of the CYP3A23 glucocorticoid-responsive region. Despite the glucocorticoid sensitivity of this gene, direct binding of ligand-activated glucocorticoid receptor (GR) to the CYP3A23 dexamethasone-responsive region (DexRE) is not required for induction. This study demonstrates that DexRE-2 is the key element within the CYP3A23 proximal promoter, conferring ligand sensitivity via its interaction with the PXR/RXRalpha heterodimer. The DexRE-1 and HNF-4 sites are not ligand-responsive, but are essential accessory elements required for full promoter inducibility. In addition to ligand-mediated activation of PXR, the overall induction response involves a GR-mediated stimulation of PXR and RXRalpha expression. Hence, the induction pathway can be divided into two stages. In stage one, maximal induction requires a GR-dependent increase in PXR and RXRalpha expression, and stage two is characterized by direct transcriptional activation of CYP3A23, which is dependent on ligand-activated PXR as well as accessory factors bound at the DexRE-1 and HNF-4 sites. Because multiple proteins bind at each element within the glucocorticoid-responsive region, factors not contributing to ligand responsiveness, such as chicken ovalbumin upstream promoter-transcription factor, may modulate the response through competitive interactions.

Differential glucocorticoid responses of CYP3A23 and CYP3A2 are mediated by selective binding of orphan nuclear receptors

CYP3A2 and CYP3A23 are two cytochrome P450 genes in rat that are differentially regulated in both their constitutive activities and their responsiveness to glucocorticoids, the prototypic CYP3A inducers. CYP3A2 displays 20-25% of the response to glucocorticoids as CYP3A23 despite extensive sequence homology in their 5'-regulatory regions. Promoter deletion analyses revealed that the CYP3A2 -57 to -168 region, homologous to the CYP3A23 dexamethasone-responsive region, mediated its low level activation. When this region was analyzed by DNase I footprinting, three binding sites were shown to correspond to the functional elements described for CYP3A23: DexRE-1, DexRE-2, and Site A (J. M. Huss and C. B. Kasper (1998) J. Biol. Chem. 273: 16155-16162). The CYP3A2 DexRE-2 and Site A elements bear two mismatches each from the CYP3A23 elements but displayed similar binding patterns in footprinting and gel-shift analyses as their CYP3A23 counterparts. The region containing 3A2DexRE-1 has six mismatches and displayed unique footprinting and gel-shift patterns compared to 3A23DexRE-1. Functional assays revealed that four mismatches within the DexRE-1 and DexRE-2 elements accounted for the differential inducibility of the two isoforms. We propose that the reduced responsiveness of CYP3A2 is the result of preferential binding of COUP-TF at the CYP3A2 DexRE-1 site. In contrast, CYP3A23 DexRE-1 associates with an accessory factor(s) that acts in concert with downstream sites to mediate the strong glucocorticoid induction response observed for CYP3A23. Site A mismatches did not influence induction magnitude but were responsible for basal activity differences. Higher CYP3A23 basal activity appears to be due to an E-box in 3A23SiteA that interacts with USF1, a ubiquitous bHLH/leucine zipper transcription factor. This site is disrupted in the corresponding 3A2SiteA. Hence, 4 nucleotide mismatches within two elements account for the difference in glucocorticoid induction, and a single mismatch is responsible for the fivefold difference in the basal activities of CYP3A2 and CYP3A23.

Mechanistic studies on the reductive half-reaction of NADPH-cytochrome P450 oxidoreductase

Site-directed mutagenesis has been employed to study the mechanism of hydride transfer from NADPH to NADPH-cytochrome P450 oxidoreductase. Specifically, Ser457, Asp675, and Cys630 have been selected because of their proximity to the isoalloxazine ring of FAD. Substitution of Asp675 with asparagine or valine decreased cytochrome c reductase activities 17- and 677-fold, respectively, while the C630A substitution decreased enzymatic activity 49-fold. Earlier studies had shown that the S457A mutation decreased cytochrome c reductase activity 90-fold and also lowered the redox potential of the FAD semiquinone (Shen, A., and Kasper, C. B. (1996) Biochemistry 35, 9451-9459). The S457A/D675N and S457A/D675N/C630A mutants produced roughly multiplicative decreases in cytochrome c reductase activity (774- and 22000-fold, respectively) with corresponding decreases in the rates of flavin reduction. For each mutation, increases were observed in the magnitudes of the primary deuterium isotope effects with NADPD, consistent with decreased rates of hydride transfer from NADPH to FAD and an increase in the relative rate limitation of hydride transfer. Asp675 substitutions lowered the redox potential of the FAD semiquinone. In addition, the C630A substitution shifted the pKa of an ionizable group previously identified as necessary for catalysis (Sem, D. S., and Kasper, C. B. (1993) Biochemistry 32, 11539-11547) from 6.9 to 7.8. These results are consistent with a model in which Ser457, Asp675, and Cys630 stabilize the transition state for hydride transfer. Ser457 and Asp675 interact to stabilize both the transition state and the FAD semiquinone, while Cys630 interacts with the nicotinamide ring and the fully reduced FAD, functioning as a proton donor/acceptor to FAD.

Nuclear receptor involvement in the regulation of rat cytochrome P450 3A23 expression

Many genes of the cytochrome P450 3A (CYP3A) subfamily, including several human and rat isoforms, are inducible by glucocorticoids. In the rat CYP3A23 gene, a 110-base pair segment of the proximal 5'-flanking region mediates dexamethasone activation. Three binding sites (DexRE-1, DexRE-2, and Site A), identified by DNase I footprinting analysis, were characterized for their relative contribution to both basal activity and dexamethasone inducibility. Site-directed mutagenesis of DexRE-1 (-144 to -169) and DexRE-2 (-118 to -136) demonstrated that each contained a core imperfect AGGTCA direct repeat, which comprised a consensus nuclear receptor binding site, and was essential for dexamethasone responsiveness but was not required for basal activity. Competition gel shift and supershift analyses revealed that both sites can bind the orphan nuclear receptor chicken ovalbumin upstream promoter-transcription factor. Site A (-85 to -110) was shown to be important for both basal activity and dexamethasone responsiveness. Point mutants displayed a reduced (2-3-fold) induction response, compared with 15-fold for wild-type, which was accompanied by a 40-60% drop in basal activity. Site A was shown to bind the liver-enriched nuclear receptor hepatocyte nuclear factor 4. Our studies demonstrate that the mechanism mediating glucocorticoid-inducible transcriptional activity of CYP3A23 involves multiple binding sites for members of the nuclear receptor superfamily.

Thyroid regulation of NADPH:cytochrome P450 oxidoreductase: identification of a thyroid-responsive element in the 5'-flank of the oxidoreductase gene

The current study demonstrates that T3-activated transcription of the NADPH:cytochrome P450 oxidoreductase (P450R) gene is dependent on the thyroid hormonal status of the animal, with both transcriptional and post-transcriptional pathways being important in regulating the cellular P450R mRNA level. The region required for transcriptional activation of the P450R gene by T3 has been identified. Nuclear run-on experiments demonstrated that the effects of T3 on P450R transcription are dependent on thyroid status, with a transcriptional enhancement obtained in T3-treated hypothyroid rat liver (1.8-fold increase) but not in T3-treated euthyroid animals. Transient cotransfection of P450R promoter/chloramphenicol acetyl transferase (CAT) constructs and the thyroid hormone receptor beta1 (TR beta1) expression plasmid into rat hepatoma H4IIE cells resulted in a 2.4-fold induction of promoter activity that was both T3 and TR beta1 dependent. Analysis of promoter deletion constructs identified a P450R-thyroid response region (P450R-TRE; bases, -564 to -536) containing three imperfect direct repeats of the thyroid response motif, AGGTCA. Mutational analysis further established that T3 induction was dependent only on the upstream direct repeat, having the sequence AGGTGAgctgAGGCCA. Footprint analysis showed that all three motifs were protected by proteins present in rat liver nuclear extracts, and a direct interaction between P450R-TRE and T3 receptors TR alpha1 and TR beta1 was demonstrated by gel-shift analysis. In vitro binding studies with P450R-TRE revealed the formation of heterodimeric complexes when TR alpha1 was coincubated with either the retinoic X receptor alpha or nuclear extract from rat liver, COS, or H4IIE cells. In addition, placement of the P450R-TRE upstream of the T3-nonresponsive heterologous thymidine kinase promoter resulted in a 2.7-fold transcriptional enhancement that was both T3 and TR beta1 dependent. Previous studies have demonstrated that T3 augments P450R mRNA levels approximately 20-30-fold and approximately 12-fold, respectively, in hypothyroid and euthyroid rats. Hence, for the hypothyroid state, transcriptional and post-transcriptional events contribute to the T3-induced mRNA increases; however, the marked increase in message level in T3-treated euthyroid animals depends primarily on post-transcriptional pathways.

Role of Ser457 of NADPH-cytochrome P450 oxidoreductase in catalysis and control of FAD oxidation-reduction potential

Site-directed mutagenesis of Ser457 of NADPH-cytochrome P450 oxidoreductase demonstrates that this residue plays a major role in both hydride transfer from NADPH to FAD and modulation of FAD redox potential. Substitution of Ser457 with alanine or cysteine decreases the rates of reduction of the substrates cytochrome c and potassium ferricyanide approximately 100-fold, while substitution with threonine produces a 20-fold decrease in activity. No changes are observed in k(m)NADPH, KiNADP+, or flavin content, indicating that these substitutions have no effect on cofactor binding but affect catalysis only. k(m)cyt c values are decreased in parallel with the observed decreases in the rates of the reductive half-reaction. Stopped-flow studies with the S457A mutant show a 100-fold decrease in the rate of flavin reduction. The primary deuterium isotope effect on Kcat for cytochrome c reduction increases from 2.7 for the wild-type enzyme to 9.0 for the S457A mutant, consistent with a change in the rate-determining step from NADP+ release in the wild-type enzyme to hydride transfer in the S457A mutant. The primary deuterium isotope effect on K1 for flavin reduction at high ionic strength (I = 535 mM) increases from 12.2 for the wild-type enzyme to > 20 for the S457A mutant, consistent again with an increase in the relative rate limitation of hydride transfer. Furthermore, anaerobic titration of S457A indicates that the redox potential of the FAD semiquinone has been decreased. Data presented in this study support the hypothesis that Ser457 is involved in hydrogen bonding interactions which stabilize both the transition state for hydride transfer and the reduced FAD.

Transcriptional regulation of the TATA-less NADPH cytochrome P-450 oxidoreductase gene

Multiple cis-acting DNA sequences regulating expression of the rat liver NADPH-cytochrome P-450 oxidoreductase gene have been identified in transient transfection assays using promoter deletion constructs linked to the chloramphenicol acetyl transferase gene. The TATA-less promoter possesses nine GC-boxes which contain the consensus sequence for the transcription factor Sp1. While loss of the seven distal GC-boxes had minimal effect on transcriptional activity, deletion of the next 35 bp, from -206 to -172, resulted in approximately 90% loss of promoter activity. Contained within this region is an Sp1 binding site indicating that either (1) this particular consensus sequence was essential for transcription, (2) the two proximal GC boxes act in concert, or (3) a yet unidentified regulatory element resides within this 35-bp stretch. In addition, transfection experiments demonstrated that two separate distal regions (-622 to -1167 and -1500 to -2300) contain negative regulatory elements which down-regulate gene transcription in a position-independent manner. Mobility-shift analyses and DNase footprinting identified sequences in the proximal region of the promoter that bound proteins present in nuclear extracts. Four protected segments were observed within the first 100 bp upstream of the transcription start site; these include (1) the region encompassing the transcription start site (-7 to +4), (2) the region normally occupied by a TATA-box (-38 to -18), (3) the bases from -78 to -60 which contain the regulatory element CACC, and (4) bases -105 to -92 which include an Sp1 binding site. Hence, regulation of the NADPH-cytochrome P-450 oxidoreductase gene is controlled by both positive and negative regulatory elements, and, of the nine Sp1 consensus sites, the two proximal sites are sufficient to support basal transcription.

Dexamethasone responsiveness of a major glucocorticoid-inducible CYP3A gene is mediated by elements unrelated to a glucocorticoid receptor binding motif

Elements responsible for dexamethasone responsiveness of CYP3A23, a major glucocorticoid-inducible member of the CYP3A gene family, have been identified. DNase I footprint analysis of the proximal promoter region revealed three protected sites (sites A, B, and C) within the sequence defined by -167 to -60. Mutational analysis demonstrated that both sites B and C were necessary for maximum glucocorticoid responsiveness and functioned in a cooperative manner. Interestingly, neither site contained a glucocorticoid responsive element. Embedded in site C was an imperfect direct repeat (5'-AACTCAAAGGAGGTCA-3'), showing homology to an AGGTCA steroid receptor motif, typically recognized by the estrogen receptor family, while site B contained an ATGAACT direct repeat; these core sequences were designated dexamethasone response elements 1 and 2 (DexRE-1 and -2), respectively. Neither element has previously been associated with a glucocorticoid-activated transcriptional response. Conversion of the DexRE-1 to either a perfect thyroid hormone or vitamin D3 responsive element further enhanced induction by dexamethasone. Gel-shift analysis demonstrated that glucocorticoid receptor did not associate with either DexRE-1 or -2; hence, glucocorticoid receptor does not directly mediate glucocorticoid induction of CYP3A23. These unusual features suggest an alternate pathway through which glucocorticoids exert their effects.

Role of acidic residues in the interaction of NADPH-cytochrome P450 oxidoreductase with cytochrome P450 and cytochrome c

Site-directed mutagenesis of the acidic clusters 207Asp-Asp-Asp209 and 213Glu-Glu-Asp215 of NADPH-cytochrome P450 oxidoreductase demonstrates that both cytochrome c and cytochrome P450 interact with this region; however, the sites and mechanisms of interaction of the two substrates are clearly distinct. Substitutions in the first acidic cluster did not affect cytochrome c or ferricyanide reductase activity, but substitution of asparagine for aspartate at position 208 reduced cytochrome P450-dependent benzphetamine N-demethylase activity by 63% with no effect on KP450m or KNADPHm. Substitutions in the second acidic cluster affected cytochrome c reduction but not benzphetamine N-demethylase or ferricyanide reductase activity. The E213Q enzyme exhibited a 59% reduction in cytochrome c reductase activity and a 47% reduction in KCyt cm under standard conditions (x0.27 M potassium phosphate, pH 7.7), as well as a decreased KCyt cm at every ionic strength and a shift of the salt dependence of cytochrome c reductase activity toward lower ionic strengths. The E214Q substitution did not affect cytochrome c reductase activity under standard conditions, but shifted the salt dependence of cytochrome c reductase activity toward higher ionic strengths. Measurements of the effect of ionic strength on steady-state kinetic properties indicated that increasing ionic strength destabilized the reductase-cytochrome c3+ ground state and reductase-cytochrome c transition state complexes for the wild-type, E213Q, and E214Q enzymes, suggesting the presence of electrostatic interactions involving Glu213 and Glu214 as well as additional residues outside this region. The ionic strength dependence of kcat/KCyt cm for the wild-type and E214Q enzymes is consistent with the presence of charge-pairing interactions in the transition state and removal of a weak ionic interaction in the reductase-cytochrome c transition-state complex by the E214Q substitution. The ionic strength dependence of the E213Q enzyme, however, is not consistent with a simple electrostatic model. Effects of ionic strength on kinetic properties of E213Q suggest that substitution of glutamine stabilizes the reductase-cytochrome c3+ ground-state complex, leading to a net increase in activation energy and decrease in kcat. Glu213 is also involved in a repulsive interaction with cytochrome c3+. Cytochrome c2+ Ki for the wild-type enzyme was 82.4 microM at 118 mM ionic strength and 10.8 microM at 749 mM ionic strength; similar values were observed for the E214Q enzyme. Cytochrome c Ki for the E213Q enzyme was 17.6 microM at 118 mM and 15.7 microM at 749 mM ionic strength, consistent with removal of an electrostatic repulsion between the reductase and cytochrome c2+.

Effect of ionic strength on the kinetic mechanism and relative rate limitation of steps in the model NADPH-cytochrome P450 oxidoreductase reaction with cytochrome c

Although the kinetic mechanism of the NADPH-cytochrome P450 oxidoreductase (P450R) reaction with cytochrome c3+ has been determined at 850 mM ionic strength [Sem, D.S., & Kasper, C. B. (1994) Biochemistry 33, 12012-12021], this mechanism is no longer valid at lower ionic strength. At 850 mM ionic strength, the mechanism is two-site ping-pong, and reaction at the electron acceptor site is itself ping-pong. As the ionic strength is decreased below 850 mM, the initial velocity profiles begin to show curvature when cytochrome c3+ is the varied substrate. These data are consistent with a mechanism that is still two-site ping-pong, but now with random sequential binding of two molecules of cytochrome c3+ at the electron acceptor site. Decreasing ionic strength also causes a change in rate-limiting steps, with (V/K)cytc and (V/K)NADPH increasing while Vmax decreases (below 500 mM ionic strength). These results are consistent with favorable ionic interactions being important for binding NADPH and cytochrome c3+ and with product (NADP+) release becoming the rate-limiting step in Vmax at low ionic strength. Vmax decreases significantly at higher ionic strength (> 500 mM), while (V/K)NADPH decreases only slightly. The DV isotope effect is largest (2.4) at 500 mM ionic strength but decreases at both low and high ionic strength as steps other than hydride transfer become more rate-limiting. D(V/K)NADPH also decreases at both low and high ionic strength, but to a lesser extent than DV.(ABSTRACT TRUNCATED AT 250 WORDS)

Kinetic mechanism for the model reaction of NADPH-cytochrome P450 oxidoreductase with cytochrome c

The kinetic mechanism of NADPH-cytochrome P450 oxidoreductase (P450R) has been determined for the model reaction with cytochrome c3+. Although initial velocity studies show parallel patterns, consistent with a classical (one-site) ping-pong mechanism that precludes the formation of a ternary NADPH-P450R-cytochrome c3+ complex, product and dead-end inhibition results suggest a nonclassical (two-site) ping-pong mechanism [Northrop, D. B. (1969) J. Biol. Chem. 244, 5808-5819]. This mechanism is a hybrid of the random sequential (ternary complex) and ping-pong mechanisms, since ternary complexes can form as well as intermediate, modified forms of the enzyme that can be present in the absence of any bound substrate. The complete rate equation is derived for this mechanism, and values for Vmax, (V/K)NADPH, (V/K)cytc, and the corresponding Michaelis constants are presented in terms of microscopic rate constants along with the expected product inhibition patterns (Appendix). Inhibition by NADP+ is competitive versus NADPH and uncompetitive versus cytochrome c3+, while inhibition by cytochrome c2+ is competitive versus cytochrome c3+ and noncompetitive versus NADPH. These inhibition patterns are consistent with the proposed two-site mechanism. This mechanism would give the same initial velocity patterns as the classical one-site ping-pong mechanism, but it allows for the formation of a ternary complex, with NADPH and cytochrome c3+ reacting independently at two separate sites on P450R. The D(V/K)NADPH isotope effect is not affected by cytochrome c3+ concentration, consistent with our assumption (in deriving the rate equation) that binding at the two sites is independent. At the high ionic strength used in this study (850 mM), the mechanism is two-site ping-pong, with the electron acceptor site itself reacting with cytochrome c3+ in a tetra uni ping-pong manner.

NADPH cytochrome P-450 oxidoreductase gene: identification and characterization of the promoter region

The untranslated first exon and the 5'-flanking region for the rat liver NADPH-cytochrome P450 oxidoreductase gene has been isolated from a Wistar-Furth genomic library. The remainder of the gene is composed of 15 exons which code for the mature protein and a 3'-nontranslated segment (T. D. Porter et al. Biochemistry, 1990, 29, 9814-9818). The 56-bp first exon resides 30.5 kb upstream from exon two, making the total gene length approximately 50 kb. While the region surrounding the start site (TCAGAGAC) was found to be homologous to a eukaryotic cap signal, the 5' flanking region possesses neither a TATA nor a CCAAT box. Instead it contains five GC-rich hexanucleotide consensus sequences for the transcription factor Sp1. These features clearly distinguish it from genes encoding other members of the mixed-function oxidase system, the cytochromes P450. Primer extension analysis and S1 nuclease mapping identified multiple transcriptional start sites. In many respects, the TATA-less oxidoreductase promoter resembles the promoter regions of dihydrofolate reductase and other housekeeping genes. Northern blot analysis demonstrates that this promoter is modulated by phenobarbital and trans-stilbene oxide, known inducers of oxidoreductase.

Enzyme-substrate binding interactions of NADPH-cytochrome P-450 oxidoreductase characterized with pH and alternate substrate/inhibitor studies

The pH dependence of the kinetic parameters for the reaction catalyzed by NADPH-cytochrome P-450 oxidoreductase (P-450R) has been determined, using various substrates and inhibitors. All Vmax and (V/K) profiles show pKas of 6.2-7.3, for an acidic group that is preferentially unprotonated for catalysis, and of 8.1-9.6, for a basic group that is preferentially protonated for catalysis. The presence of the wrong ionization state for both of these groups is tolerated more at lower ionic strength (300 mM) than at higher ionic strength (850 mM). Ionization of the basic group has a more pronounced effect on binding of substrate (cytochrome c or dichloroindophenol) than on catalysis, since ionization has only a 2-fold effect on Vmax with cytochrome c, and only a 5-fold effect on Vmax with dichloroindophenol, while (V/K) for both substrates continues to drop at high pH with no sign of reaching a plateau. Therefore, this basic group affects predominantly substrate binding and, to a lesser extent, catalysis. It is most likely located on the surface of the protein at the cytochrome c/dichloroindophenol binding site, near the FMN prosthetic group. The NADP+ pKi profile shows a pKa of 5.95 for the 2'-phosphate of NADP+, which is bound to P-450R as the dianion, and a pKa of 9.53 for an enzyme group that must be protonated in order to bind NADP+.(ABSTRACT TRUNCATED AT 250 WORDS)

Interaction with arginine 597 of NADPH-cytochrome P-450 oxidoreductase is a primary source of the uniform binding energy used to discriminate between NADPH and NADH

Site-directed mutagenesis has been used in conjunction with pH and alternate substrate/inhibitor studies to characterize the interactions between NADPH-cytochrome P-450 oxidoreductase (P-450R) and the 2'-phosphate of NADP(H) that provide P-450R with its strong nicotinamide nucleotide specificity. It is known that the 2'-phosphate of NADP(H) is bound to P-450R as the dianion and that interactions between it and residues on P-450R provide 5 kcal/mol of essentially uniform binding energy (preceding paper in this issue). In order to probe these interactions further, Arg597 of P-450R, which is homologous to Arg235 of ferredoxin-NADP+ reductase that forms a salt bridge with the 2'-phosphate of 2'-phospho-AMP in the crystal structure of that complex [Karplus, P. A., Daniels, M. J., & Herriott, J. R. (1991) Science 251, 60], was mutated to methionine. The mutant protein, P-450R (R597M), does not appear to have a grossly perturbed tertiary structure on the basis of the observation of similar 31P-NMR chemical shifts for FAD (pyrophosphate) bound to it and wild-type (WT) P-450R, although it is more unstable to urea denaturation. P-450R (R597M) has a Km for NADPH that is 150 times that of P-450R (WT) and a Ki for NADP+ that is 240 times that of P-450R (WT). In contrast, the R597M mutation has only a modest effect on the Km for NADH (0.8 WT) and the Ki for NAD+ (2.9 WT), indicating that Arg597 must have been interacting specifically with the 2'-phosphate of NADP(H). The R597M mutation has relatively little effect on kcat for NADPH (1.2 WT) or NADH (0.6 WT), indicating that the mutation is affecting ground and transition states to essentially the same degree, by removing 3 kcal/mol of uniform binding energy. The NADP+ pKi profile for P-450R (R597M) shows a pKa of 5.78 for the 2'-phosphate of NADP+, which is bound to P-450R (R597M) as the dianion, but the pKa of 9.5 for the preferentially protonated enzymic group observed in the P-450R (WT) profile is no longer present. It is argued then that the 2'-phosphate binding pocket of P-450R (WT) has a high positive charge density (> + 2) and that Arg597, which is in this binding pocket, has a highly perturbed pKa of 9.5. Finally, a general theoretical treatment of the thermodynamic consequences of individual and combined perturbations to complementary interacting groups on enzyme and substrate is presented (see Appendix).(ABSTRACT TRUNCATED AT 400 WORDS)

Expression of rat microsomal epoxide hydrolase in Escherichia coli. Identification of a histidyl residue essential for catalysis

The cDNA containing the complete coding region for rat microsomal epoxide hydrolase (EC 3.3.2.3) was cloned into the expression/secretion vector pIN-III-OmpA3 and expressed in Escherichia coli strain TG1. Recombinant epoxide hydrolase was found to represent 4-9% of total bacterial protein and catalyzed the hydrolysis of styrene oxide and benzo[a]pyrene 4,5-oxide with specific activities of 421 and 734 nmol min-1 mg of epoxide hydrolase-1, respectively. Previous work implicated a histidyl residue at or near the active site of the enzyme (DuBois, G. C., Appella, E., Levin, W., Lu, A. Y. H., and Jerina, D. M. (1978) J. Biol. Chem. 253, 2932-2939). Comparison of the amino acid sequences of rat, human, and rabbit epoxide hydrolases revealed the presence of 14 conserved histidyl residues. To investigate the role of these residues in epoxide hydrolysis, site-specific mutants were generated and expressed in E. coli. Mutants H64L, H82L, H115N, H126N, H129L, H148N, H170L, H176L, H242L, H247L, H301L, H385L, K386M-H387L, delta 385-391, and H407L catalyzed the hydrolysis of benzo[a]pyrene 4,5-oxide with specific activities between 115 and 830 nmol min-1 mg-1. Mutants H431L, H431N, and H431R were all found to have activities of < 5 nmol min-1 mg-1, which is at least 150-fold less than the activity of the wild type enzyme. A Vm versus pH profile for the recombinant wild type epoxide hydrolase revealed a broad pH optimum of 6.5 to 8.5 and the presence of three ionizable groups with pKa values of 5.8 +/- 0.2, 9.2 +/- 0.1, and 9.7 +/- 0.4. The group with a pKa of 5.8 is preferentially unprotonated, while the other two groups are preferentially protonated for catalysis. We propose that histidine 431 corresponds to the group with a pKa of 5.8, while the others, with pKa values of 9.2 and 9.7 likely represent lysyl, cysteinyl, or tyrosyl residues. Thus, the data are consistent with a model where His-431 acts as a general base, abstracting a proton from water, while another residue(s), perhaps lysine, act as a general acid protonating the alkoxide anion that forms upon cleavage of the carbon-oxygen bond.

31P NMR spectroscopic studies on purified, native and cloned, expressed forms of NADPH-cytochrome P450 reductase

31P NMR spectroscopy has been utilized in conjunction with site-directed mutagenesis and phospholipid analysis to determine structural aspects of the prosthetic flavins, FAD and FMN, of NADPH-cytochrome P450 reductase. Comparisons are made among detergent-solubilized and protease (steapsin)-solubilized preparations of porcine liver reductases, showing unequivocally that the 31P NMR signals at approximately 0.0 ppm in the detergent-solubilized, hydrophobic form are attributable to phospholipids. By extraction and TLC analysis, the phospholipid contents of detergent-solubilized rat liver reductase, both tissue-purified and Escherichia coli-expressed, have been determined to reflect the membranes from which the enzyme was extracted. In addition, the cloned, wild-type NADPH-cytochrome P450 reductase exhibits an additional pair of signals downfield of the normal FAD pyrophosphate resonances reported by Otvos et al. [(1986) Biochemistry 25, 7220-7228], but these signals are not observed with tissue-purified or mutant enzyme preparations. The Tyr140----Asp140 mutant, which exhibits only 20% of wild-type activity, displays no gross changes in 31P NMR spectra. However, the Tyr178----Asp178 mutant, which has no catalytic activity and does not bind FMN, exhibits no FMN 31P NMR signal and a normal, but low intensity, pair of signals for FAD. The latter experiments, taking advantage of mutations in residues putatively on either side of the FMN isoalloxazine ring, suggest subtle to severe changes in the binding of the flavin prosthetic groups and, perhaps, cooperative interactions of flavin binding to NADPH-cytochrome P450 reductase.

Geometric relationship between the nicotinamide and isoalloxazine rings in NADPH-cytochrome P-450 oxidoreductase: implications for the classification of evolutionarily and functionally related flavoproteins

The stereospecificity of hydride abstraction from NADPH and the conformation of the nicotinamide ring around the glycosidic bond have been determined for the flavoprotein NADPH-cytochrome P-450 oxidoreductase (P-450R). The A-side (pro-R) hydrogen is abstracted from NADPH, and the nicotinamide ring is in the anti conformation. These results are consistent with the apparently strong correlation between A-side stereospecificity and anti conformation and between B-side stereospecificity and syn conformation [You, K. (1985) CRC Crit. Rev. Biochem. 17, 313]. This correlation reveals how the flavin and nicotinamide rings are oriented relative to each other. In P-450R, the flavin is then "on top of" (on the exo side of) the nicotinamide ring. In another flavoprotein dehydrogenase, glutathione reductase, which is a B-side/anti enzyme [Pai, E. F., & Schulz, G. E. (1983) J. Biol. Chem. 258, 1752], the flavin is "underneath" (on the endo side of) the nicotinamide ring. We argue that all enzymes that are evolutionarily related to these two flavoproteins should have their respective overall configurations. The overall configuration is defined by the following five properties: (1) relative orientation of the isoalloxazine and nicotinamide rings, (2) stereospecificity of hydride transfer to/from the nicotinamide ring, (3) conformation of the nicotinamide ring around the glycosidic bond, (4) stereospecificity of hydride transfer to/from the flavin, and (5) conformation of the flavin around its N5-N10 axis. There are only eight possible overall configurations, and a knowledge of only three of the five properties is needed to determine which one is present (as long as the combination of properties is not 1, 2, 3 or 1, 4, 5).(ABSTRACT TRUNCATED AT 250 WORDS)

NADPH-cytochrome P-450 oxidoreductase. The role of cysteine 566 in catalysis and cofactor binding

Site-directed mutagenesis was employed to investigate the role of Cys566 in the catalytic mechanism of rat liver NADPH-cytochrome P-450 oxidoreductase. Rat NADPH-cytochrome P-450 oxidoreductase and mutants containing either alanine or serine at position 566 were expressed in Escherichia coli and purified to homogeneity. Substitution of alanine at position 566 had no effect on enzymatic activity with the acceptors cytochrome c and ferricyanide but did increase trans-hydrogenase activity with 3-acetylpyridine adenine dinucleotide phosphate by 79%. The Km for NADPH was increased 2.5-fold, and the NADP+ KI was increased 4.8-fold compared with that found for the wild-type enzyme. The conservative substitution, Ser566, produced a 50% decrease in cytochrome c reductase activity whereas activity with ferricyanide was decreased 57%, and 3-acetylpyridine adenine dinucleotide phosphate activity was unaffected. The NADPH Km was increased 4.6-fold, and the NADP+ KI increased 7.6-fold. The dependence of cytochrome c reductase activity on the KCl concentration was markedly altered by the Cys566 substitutions. Maximum activity for the wild-type enzyme was observed at approximately 0.18 M KCl whereas maximum activity for the mutant enzymes was observed between 0.04 and 0.09 M KCl. The pH dependence of cytochrome c reductase activity, cytochrome c Km, and flavin content were unaffected by these substitutions. These results demonstrate that Cys566 is not essential for activity of rat liver NADPH-cytochrome P-450 oxidoreductase although the cysteine side chain does affect the interaction of NADPH with the enzyme.

Transcriptional regulation and localization of the tissue-specific induction of epoxide hydrolase by lead acetate in rat kidney

Intraperitoneal administration of lead acetate to male Sprague-Dawley rats resulted in the tissue-specific transcriptional activation of the microsomal epoxide hydrolase gene in kidney. This response was followed by a 15-fold increase in the level of kidney epoxide hydrolase mRNA, while no change in mRNA level was noted in liver. Treated animals also showed no increase in mRNA levels for NADPH-cytochrome P-450 oxidoreductase, cytochrome P-450b/e, cytochrome P-450PCN (where PCN is pregnenolone 16 alpha-carbonitrile), or serum albumin in either kidney or liver. Immunoquantitation of the enzyme revealed a 16-fold increase in kidney upon treatment with lead acetate, but significant changes in epoxide hydrolase levels were not noted in liver, heart, spleen, lung, small intestine, or testis. The enzymatic activity for liver and kidney paralleled the immunochemical results; however, the activity increase in kidney was only one-third of the increase noted for total enzyme protein. Immunohistochemical analysis of epoxide hydrolase protein in sections of rat kidney demonstrated that in lead acetate-treated animals there was a marked increase in staining of the cytoplasm of the proximal tubular cells in the outer cortex as compared with kidneys from control animals. In contrast, considerable protein was also localized to collecting ducts, but no change was evident in the content of the epoxide hydrolase gene product in these structures in control and lead acetate-treated animals. Immunohistochemical differences were not noted between livers from control and lead-treated animals. Furthermore, the staining patterns for NADPH-cytochrome P-450 oxidoreductase were the same for control and treated animals in both kidney and liver. Quantitative measurements of lead uptake by various rat tissues showed liver, spleen, and small intestine reaching a maximum of approximately 12,000 ng of lead/g of dry tissue at 8, 8, and 16 h, respectively, while kidney, lung, and testis peaked (approximately 3,000 ng of lead/g of dry tissue) at 16, 16, and 12 h, respectively.

Activation of mitomycin C by NADPH:cytochrome P-450 reductase

Mitomycin C is an alkylating agent used in cancer chemotherapy that shows some specificity towards hypoxic cells. The therapeutic effects of this compound are thought to result from its metabolic activation by enzymes such as NADPH:cytochrome P-450 reductase. In a previous report we described a Chinese hamster ovary cell line resistant to mitomycin C, which had a decreased NADPH:cytochrome P-450 reductase activity coupled with a lower rate of mitomycin C metabolism. In order to provide further evidence that the lower reductase activity is a factor in the resistance mechanism, we incorporated NADPH:cytochrome P-450 reductase into cytotoxicity assays and showed that it significantly sensitizes cells to mitomycin C. Also, the difference in drug sensitivity between the wild-type and drug-resistant Chinese hamster ovary cells was no longer observed. In addition to these studies, we expressed a rat liver NADPH:cytochrome P-450 reductase cDNA in a Salmonella typhimurium strain, LR5000. The bacteria expressing the rat NADPH: cytochrome P-450 reductase showed increased sensitivity to mitomycin C when incubated with this compound under aerobic conditions. However, under hypoxic conditions increased sensitivity was not observed. This parallels the previous finding with mitomycin C-resistant Chinese hamster ovary cells. These data provide direct evidence for the role of NADPH:cytochrome P-450 reductase in the cytotoxic action of this mitomycin C under aerobic but not hypoxic conditions and suggest that reduced levels of this enzyme can lead to drug resistance. P-450 reductase expressed in S. typhimurium may provide a valuable tool for evaluating the role of this enzyme in the toxicity of drugs activated through a one electron reduction pathway.

NADPH-cytochrome P-450 oxidoreductase gene organization correlates with structural domains of the protein

cDNA clones to rat liver NADPH-cytochrome P-450 oxidoreductase were used to isolate genomic clones from a Wistar-Furth inbred rat genomic DNA library. Fifteen exons containing the coding region and 3'-nontranslated segment of the P-450 reductase gene were identified, spanning 20 kilobases of DNA contained in 3 lambda-Charon 35 clones. The organization of this single copy gene reveals a general correspondence between exons and structural domains of the protein, with the segment responsible for anchoring the reductase to the microsomal membrane and several segments involved in FMN, FAD, and NADPH binding encoded by discrete exons.

Glucocorticoid repression and basal regulation of the epoxide hydrolase promoter

Through a series of promoter deletions and gene transfer experiments we have examined the basal regulation and glucocorticoid-mediated repression of the rat epoxide hydrolase gene. Three regions of the 5' flanking sequence were found to influence the basal level of promoter function in H4IIE hepatoma cells. Region A (-891 to -355 bp) contains an apparent repressor of epoxide hydrolase expression, while regions B (-271 to -171 bp) and C (-141 to -85) were found to contain important sequences required for optimal promoter activity. Previous work has demonstrated that dexamethasone represses epoxide hydrolase transcription by approximately 50% in isolated rat liver nuclei, and, in this study, we have demonstrated that the ability of the epoxide hydrolase promoter to drive CAT expression is similarly repressed in H4IIE cells treated with 1 microM dexamethasone. Furthermore, the level of endogenous epoxide hydrolase mRNA is decreased by 70-88% in nontransfected H4IIE cells treated with dexamethasone. Interestingly, promoter activity was not decreased by dexamethasone in COS cells, which lack glucocorticoid receptors. The current data show that sequences from -42 to +110 bp are sufficient to support the dexamethasone response, and, furthermore, they suggest that repression may not require direct interaction of the ligand-receptor complex with the promoter region.

Quantitation of mRNAs specific for the mixed-function oxidase system in rat liver and extrahepatic tissues during development

Evaluation of ontogenetic expression of the cytochrome P450PCN and cytochrome P450b gene families as well as the NADPH-cytochrome P450 oxidoreductase and epoxide hydrolase genes in Holtzmann rats showed that basal levels of mRNAs encoding these enzymes could be detected in most tissues. Distinct developmental patterns of mRNA expression are evident for these four proteins in liver and extrahepatic tissues. Levels of cytochrome P450b-like mRNA were comparable in adult lung and liver, while cytochrome P450PCN-homologous mRNA exhibited low levels in lung and approximately 100-fold higher levels in liver. Cytochrome P450PCN-homologous mRNA also reached substantial levels in adult intestine, and was also present in placenta, where it increased approximately 4-fold 24 h before birth. Epoxide hydrolase mRNA was demonstrated to be highest in liver followed by kidney, lung, and intestine but was extremely low in brain. NADPH-cytochrome P450 oxidoreductase mRNA in kidney, lung, prostate, adrenal, and intestine exhibited levels comparable to that found in liver; however, the pattern of expression for oxidoreductase mRNA was unique in that levels declined at maturity in liver, kidney, and intestine but not in lung and brain. Development of mixed-function oxidase and epoxide hydrolase activities in liver was distinct from that in other tissues in that mRNAs for all four proteins rose dramatically after parturition. Testis from immature males demonstrated low levels of all the mRNAs assayed, which ranged from 20% (oxidoreductase) to less than 1% (cytochrome P450PCN and epoxide hydrolase) of the levels found in liver.

Structural analysis of the FMN binding domain of NADPH-cytochrome P-450 oxidoreductase by site-directed mutagenesis

Comparison of the amino acid sequence of rat liver NADPH-cytochrome P-450 oxidoreductase with that of flavoproteins of known three-dimensional structure suggested that residues Tyr-140 and Tyr-178 are involved in binding of FMN to the protein. To test this hypothesis, NADPH-cytochrome P-450 oxidoreductase was expressed in Escherichia coli using the expression-secretion vector pIN-III-ompA3, and site-directed mutagenesis was employed to selectively alter these residues and demonstrate that they are major determinants of the FMN-binding site. Bacterial expression produced a membrane-bound 80-kDa protein containing 1 mol each of FMN and FAD per mol of enzyme, which reduced cytochrome c at a rate of 51.5 mumol/min/mg of protein and had absorption spectra and kinetic properties very similar to those of the rat liver enzyme. Replacement of Tyr-178 with aspartate abolished FMN binding and cytochrome c reductase activity. Incubation with FMN increased catalytic activity to a maximum of 8.6 mumol/min/mg of protein. Replacement of Tyr-140 with aspartate did not eliminate FMN binding, but reduced cytochrome c reductase activity about 5-fold, suggesting that FMN may be bound in a conformation which does not permit efficient electron transfer. Substitution of phenylalanine at either position 140 or 178 had no effect on FMN content or catalytic activity. The FAD level in the Asp-178 mutant was also decreased, suggesting that FAD binding is dependent upon FMN; FAD incorporation may occur co-translationally and require prior formation of an intact FMN domain.

Structure and organization of the microsomal xenobiotic epoxide hydrolase gene

The gene for the microsomal xenobiotic rat liver epoxide hydrolase has been isolated and characterized. Clones were obtained from a Wistar Furth Charon 35 genomic library by hybridization with a full-length epoxide hydrolase cDNA. The gene for the xenobiotic epoxide hydrolase is approximately 16 kilobases in length and consists of 9 exons ranging in size from 109 to 420 base pairs and 8 intervening sequences, the largest of which is 3.2 kilobases. S1-nuclease mapping, primer extension studies, and sequence analysis were used to determine the 5' cap site and the size of the first exon (170 base pairs). Regulatory sequences analogous to TATA, CCAAT, and core enhancer sequences were noted in the 5'-flanking region of the gene. The cDNA and gene for epoxide hydrolase displayed nucleotide sequence identity although they were isolated from different rat strains. Also, Southern blot analysis of restricted liver DNA from inbred Fischer 344 and Wistar Furth rat strains, and outbred Sprague-Dawley rats indicated a high degree of structural similarity for the epoxide hydrolase gene within these three strains. Only a single functional epoxide hydrolase gene was identified and no evidence of hybridization to the genes for the microsomal cholesterol epoxide hydrolase or the cytosolic epoxide hydrolase was observed. However, a pseudogene for the microsomal xenobiotic epoxide hydrolase was isolated and characterized from the genomic library.