The P450 superfamily: updated listing of all genes and recommended nomenclature for the chromosomal loci

Gene-Environment Interactions: My Unique Journey

I am deeply honored to be invited to write this scientific autobiography. As a physician-scientist, pediatrician, molecular biologist, and geneticist, I have authored/coauthored more than 600 publications in the fields of clinical medicine, biochemistry, biophysics, pharmacology, drug metabolism, toxicology, molecular biology, cancer, standardized gene nomenclature, developmental toxicology and teratogenesis, mouse genetics, human genetics, and evolutionary genomics. Looking back, I think my career can be divided into four distinct research areas, which I summarize mostly chronologically in this article: (a) discovery and characterization of the AHR/CYP1 axis, (b) pharmacogenomics and genetic prediction of response to drugs and other environmental toxicants, (c) standardized drug-metabolizing gene nomenclature based on evolutionary divergence, and (d) discovery and characterization of the SLC39A8 gene encoding the ZIP8 metal cation influx transporter. Collectively, all four topics embrace gene-environment interactions, hence the title of my autobiography.

Human Cytochrome P450 1, 2, 3 Families as Pharmacogenes with Emphases on Their Antimalarial and Antituberculosis Drugs and Prevalent African Alleles

Precision medicine gives individuals tailored medical treatment, with the genotype determining the therapeutic strategy, the appropriate dosage, and the likelihood of benefit or toxicity. Cytochrome P450 (CYP) enzyme families 1, 2, and 3 play a pivotal role in eliminating most drugs. Factors that affect CYP function and expression have a major impact on treatment outcomes. Therefore, polymorphisms of these enzymes result in alleles with diverse enzymatic activity and drug metabolism phenotypes. Africa has the highest CYP genetic diversity and also the highest burden of malaria and tuberculosis, and this review presents current general information on CYP enzymes together with variation data concerning antimalarial and antituberculosis drugs, while focusing on the first three CYP families. Afrocentric alleles such as CYP2A6*17, CYP2A6*23, CYP2A6*25, CYP2A6*28, CYP2B6*6, CYP2B6*18, CYP2C8*2, CYP2C9*5, CYP2C9*8, CYP2C9*9, CYP2C19*9, CYP2C19*13, CYP2C19*15, CYP2D6*2, CYP2D6*17, CYP2D6*29, and CYP3A4*15 are implicated in diverse metabolic phenotypes of different antimalarials such as artesunate, mefloquine, quinine, primaquine, and chloroquine. Moreover, CYP3A4, CYP1A1, CYP2C8, CYP2C18, CYP2C19, CYP2J2, and CYP1B1 are implicated in the metabolism of some second-line antituberculosis drugs such as bedaquiline and linezolid. Drug-drug interactions, induction/inhibition, and enzyme polymorphisms that influence the metabolism of antituberculosis, antimalarial, and other drugs, are explored. Moreover, a mapping of Afrocentric missense mutations to CYP structures and a documentation of their known effects provided structural insights, as understanding the mechanism of action of these enzymes and how the different alleles influence enzyme function is invaluable to the advancement of precision medicine.

Placing human gene families into their evolutionary context

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Insights into the Cytochrome P450 Monooxygenase Superfamily in Osmanthus fragrans and the Role of OfCYP142 in Linalool Synthesis

Osmanthus fragrans flowers have long been used as raw materials in food, tea, beverage, and perfume industries due to their attractive and strong fragrance. The P450 superfamily proteins have been reported to widely participate in the synthesis of plant floral volatile organic compounds (VOCs). To investigate the potential functions of P450 superfamily proteins in the fragrance synthesis of O. fragrans, we investigated the P450 superfamily genome wide. A total of 276 P450 genes were identified belonging to 40 families. The RNA-seq data suggested that many OfCYP genes were preferentially expressed in the flower or other organs, and some were also induced by multiple abiotic stresses. The expression patterns of seven flower-preferentially expressed OfCYPs during the five different flower aroma content stages were further explored using quantitative real-time PCR, showing that the CYP94C subfamily member OfCYP142 had the highest positive correlation with linalool synthesis gene OfTPS2. The transient expression of OfCYP142 in O. fragrans petals suggested that OfCYP142 can increase the content of linalool, an important VOC of the O. fragrans floral aroma, and a similar result was also obtained in flowers of OfCYP142 transgenic tobacco. Combined with RNA-seq data of the transiently transformed O. fragrans petals, we found that the biosynthesis pathway of secondary metabolites was significantly enriched, and many 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway genes were also upregulated. This evidence indicated that the OfCYP proteins may play critical roles in the flower development and abiotic response of O. fragrans, and that OfCYP142 can participate in linalool synthesis. This study provides valuable information about the functions of P450 genes and a valuable guide for studying further functions of OfCYPs in promoting fragrance biosynthesis of ornamental plants.

Feedback Inhibition Might Dominate the Accumulation Pattern of BR in the New Shoots of Tea Plants (Camellia sinensis)

Brassinosteroid (BR), a kind of polyhydroxylated steroid hormone, plays an important role in physiological and biochemical processes in plants. Studies were mainly focused on BR signaling and its exogenous spraying to help enhance crop yields. Few research studies are centered on the accumulation pattern of BR and its mechanism. Yet, it is crucial to unlock the mystery of the function of BR and its cross action with other hormones. Tea (Camellia sinensis (L.) O. Kuntze) is one of the important economic crops in some countries, and new shoots are the raw materials for the preparation of various tea products. Different concentrations of exogenous BR were reported to have different effects on growth and development. New shoots of tea plants can thus be considered a valuable research object to study the accumulation pattern of BR. In this study, the quantity of five BR components (brassinolide, 28-norbrassinolide, 28-homobrassinolide, castasterone, and 28-norcastasterone) in different tissues of tea plants, including buds (Bud), different maturity of leaves (L1, L2), and stems (S1, S2) were determined by UPLC-MS/MS. A total of 15 cDNA libraries of the same tissue with three repetitions for each were constructed and sequenced. The BR-accumulation pattern and gene expression pattern were combined together for weighted gene co-expression network analysis (WGCNA). BR-accumulation-relative genes were then screened using two methods, based on the K.in value and BR biosynthetic pathway (ko00905), respectively. The result showed that photosynthesis-related genes and CYP450 family genes were actively involved and might play important roles in BR accumulation and/or its accumulation pattern. First and foremost, feedback inhibition was more likely to dominate the accumulation pattern of BR in the new shoots of tea plants. Moreover, three conserved miRNAs with their target transcriptional factors and target mRNAs had been figured out from negative correlation modules that might be strongly linked to the BR-accumulation pattern. Our study provided an experimental basis for the role of BR in tea plants. The excavation of genes related to the accumulation pattern of BR provided the possibility of cross-action studies on the regulation of BR biosynthesis and the study between BR and other hormones.

Mutation Types of CYP71P1 Cause Different Phenotypes of Mosaic Spot Lesion and Premature Leaf Senescence in Rice

Lesion mimic mutants (LMMs) are ideal materials for studying programmed cell death and defense response in plants. Here we report investigations on two LMMs (msl-1 and msl-2) from the indica rice cultivar JG30 treated by ethyl methyl sulfone. Both of the mutants showed similar mosaic spot lesions at seedling stage, but they displayed different phenotypes along with development of the plants. At tillering stage, larger orange spots appeared on leaves of msl-2, while only small reddish-brown spots exhibit on leaves of msl-1. At heading stage, the msl-2 plants were completely dead, while the msl-1 plants were still alive even if showed apparent premature senility. For both the mutants, the mosaic spot lesion formation was induced by light; DAB and trypan blue staining showed a large amount of hydrogen peroxide accumulated at the lesion sites, accompanied by a large number of cell death. Consequently, reactive oxygen species were enriched in leaves of the mutants; SOD and CAT activities in the scavenging enzyme system were decreased compared with the wild type. In addition, degraded chloroplasts, decreased photosynthetic pigment content, down-regulated expression of genes associated with chloroplast synthesis/photosynthesis and up-regulated expression of genes related to senescence were detected in the mutants, but the abnormality of msl-2 was more serious than that of msl-1 in general. Genetic analysis and map-based cloning revealed that the lesion mimic and premature senescence traits of both the mutants were controlled by recessive mutated alleles of the SL (Sekiguchi lesion) gene, which encodes the CYP71P1 protein belonging to cytochrome P450 monooxygenase family. The difference of mutation sites and mutation types (SNP-caused single amino acid change and SNP-caused early termination of translation) led to the different phenotypes in severity between msl-1 and msl-2. Taken together, this work revealed that the CYP71P1 is involved in regulation of both premature senescence and cell death in rice, and its different mutation sites and mutation types could cause different phenotypes in terms of severity.

Cytochromes P450 (P450s): A review of the class system with a focus on prokaryotic P450s

Cytochromes P450 (P450s) are a large superfamily of heme-containing monooxygenases. P450s are found in all Kingdoms of life and exhibit incredible diversity, both at sequence level and also on a biochemical basis. In the majority of cases, P450s can be assigned into one of ten classes based on their associated redox partners, domain architecture and cellular localization. Prokaryotic P450s now represent a large diverse collection of annotated/known enzymes, of which many have great potential biocatalytic potential. The self-sufficient P450 classes (Class VII/VIII) have been explored significantly over the past decade, with many annotated and biochemically characterized members. It is clear that the prokaryotic P450 world is expanding rapidly, as the number of published genomes and metagenome studies increases, and more P450 families are identified and annotated (CYP families).

The Use of Liver Cell Cultures Derived from Different Mammalian Species in In Vitro Toxicological Studies: Implementation in Extrapolation Models?

In vitro techniques can be applied to obtain insight into the organotoxic potential of compounds and in the underlying mechanisms of action. During the last few decades the applicability of these techniques in toxicology has been demonstrated for an increasing number of cell types. Cultured hepatocytes have been shown to provide convenient in vitro systems for studying the role of biotransformation in the toxicity of a compound. The use of in vitro systems derived from a range of species may lead to a better understanding of species differences in the toxicity of xenobiotics, thus improving interspecies extrapolations. The applicability of this approach will be illustrated with results from experiments on the biotransformation and cytotoxicity of bromobenzene, the in vitro inducibility of biotransformation enzyme systems and the effects of hypolipidemics in primary hepatocyte cultures derived from a numberof animals.

Sense of Smell: Structural, Functional, Mechanistic Advancements and Challenges in Human Olfactory Research

Olfaction, the sense of smell detects and discriminate odors as well as social cues which influence our innate responses. The olfactory system in human beings is found to be weak as compared to other animals; however, it seems to be very precise. It can detect and discriminate millions of chemical moieties (odorants) even in minuscule quantities. The process initiates with the binding of odorants to specialized olfactory receptors, encoded by a large family of Olfactory Receptor (OR) genes belonging to the G-protein-coupled receptor superfamily. Stimulation of ORs converts the chemical information encoded in the odorants, into respective neuronal action-potentials which causes depolarization of olfactory sensory neurons. The olfactory bulb relays this signal to different parts of the brain for processing. Odors are encrypted using a combinatorial approach to detect a variety of chemicals and encode their unique identity. The discovery of functional OR genes and proteins provided an important information to decipher the genomic, structural and functional basis of olfaction. ORs constitute 17 gene families, out of which 4 families were reported to contain more than hundred members each. The olfactory machinery is not limited to GPCRs; a number of non- GPCRs is also employed to detect chemosensory stimuli. The article provides detailed information about such olfaction machinery, structures, transduction mechanism, theories of odor perception, and challenges in the olfaction research. It covers the structural, functional and computational studies carried out in the olfaction research in the recent past.

Development and Uses of Offline and Web-Searchable Metabolism Databases - The Case of Benzo[a]pyrene

BACKGROUND

The present work describes development of offline and web-searchable metabolism databases for drugs, other chemicals, and physiological compounds using human and model species, prompted by the large amount of data published after year 1990. The intent was to provide a rapid and accurate approach to published data to be applied both in science and to assist therapy.

METHODS

Searches for the data were done using the Pub Med database, accessing the Medline database of references and abstracts. In addition, data presented at scientific conferences (e.g., ISSX conferences) are included covering the publishing period beginning with the year 1976.

RESULTS

Application of the data is illustrated by the properties of benzo[a]pyrene (B[a]P) and its metabolites. Analysis show higher activity of P450 1A1 for activation of the (-)- isomer of trans-B[a]P-7,8-diol, while P4501B1 exerts higher activity for the (+)- isomer. P450 1A2 showed equally low activity in the metabolic activation of both isomers.

CONCLUSION

The information collected in the databases is applicable in prediction of metabolic drug-drug and/or drug-chemical interactions in clinical and environmental studies. The data on the metabolism of searched compound (exemplified by benzo[a]pyrene and its metabolites) also indicate toxicological properties of the products of specific reactions. The offline and web-searchable databases had wide range of applications (e.g. computer assisted drug design and development, optimization of clinical therapy, toxicological applications) and adjustment in everyday life styles.

Structure and function of the cytochrome P450 peroxygenase enzymes

The cytochromes P450 (P450s or CYPs) constitute a large heme enzyme superfamily, members of which catalyze the oxidative transformation of a wide range of organic substrates, and whose functions are crucial to xenobiotic metabolism and steroid transformation in humans and other organisms. The P450 peroxygenases are a subgroup of the P450s that have evolved in microbes to catalyze the oxidative metabolism of fatty acids, using hydrogen peroxide as an oxidant rather than NAD(P)H-driven redox partner systems typical of the vast majority of other characterized P450 enzymes. Early members of the peroxygenase (CYP152) family were shown to catalyze hydroxylation at the α and β carbons of medium-to-long-chain fatty acids. However, more recent studies on other CYP152 family P450s revealed the ability to oxidatively decarboxylate fatty acids, generating terminal alkenes with potential applications as drop-in biofuels. Other research has revealed their capacity to decarboxylate and to desaturate hydroxylated fatty acids to form novel products. Structural data have revealed a common active site motif for the binding of the substrate carboxylate group in the peroxygenases, and mechanistic and transient kinetic analyses have demonstrated the formation of reactive iron-oxo species (compounds I and II) that are ultimately responsible for hydroxylation and decarboxylation of fatty acids, respectively. This short review will focus on the biochemical properties of the P450 peroxygenases and on their biotechnological applications with respect to production of volatile alkenes as biofuels, as well as other fine chemicals.

Roles of Cytochrome P450 in Metabolism of Ethanol and Carcinogens

Cytochrome P450 (P450) enzymes are involved in the metabolism of carcinogens, as well as drugs, steroids, vitamins, and other classes of chemicals. P450s also oxidize ethanol, in particular P450 2E1. P450 2E1 oxidizes ethanol to acetaldehyde and then to acetic acid, roles also played by alcohol and aldehyde dehydrogenases. The role of P450 2E1 in cancer is complex in that P450 2E1 is also induced by ethanol, P450 2E1 is involved in the bioactivation and detoxication of a number of chemical carcinogens, and ethanol is an inhibitor of P450 2E1. Contrary to some literature, P450 2E1 expression and induction itself does not cause global oxidative stress in vivo, as demonstrated in studies using isoniazid treatment and gene deletion studies with rats and mice. However, a major fraction of P450 2E1 is localized in liver mitochondria instead of the endoplasmic reticulum, and studies with site-directed rat P450 2E1 mutants and natural human P450 2E1 N-terminal variants have shown that P450 2E1 localized in mitochondria is catalytically active and more proficient in producing reactive oxygen species and damage. The role of the mitochondrial oxidative stress in ethanol toxicity is still under investigation, as is the mechanism of altered electron transport to P450s that localize inside mitochondria instead of their typical endoplasmic reticulum environment.

The Role of Metabolism in Chemical-Induced Pulmonary Toxicity

The lung is a target organ for the toxic effects of several chemical agents, including natural products, industrial chemicals, pesticides, environmental agents, and occasionally, drugs. Factors that establish the lung as a target organ include selective tissue exposure, high tissue oxygenation, and the presence of bioactivating systems that can generate toxic products from initially innocuous substances. Selective pulmonary exposure most often results from the fact that the lung serves as the major portal of entry for most airborne substances, but in some cases, selective exposure is the consequence of accumulation of agents, such as certain basic amines, from the circulation. Lung tumor development following long-term exposure to cigarette smoke or diesel engine exhaust is an example of pulmonary toxicity resulting from selective external exposure. Selective internal exposure, on the other hand, is exemplified by the pulmonary uptake of the herbicide paraquat from the circulation which is in part responsible for its lung-toxic effects. Although the lung contains drug metabolizing enzymes similar to those found in the liver, the enzymatic systems in the lung are sometimes highly concentrated in specific cell types. In the rabbit, for example, the lung-selective toxicity of the natural product ipomeanol is the consequence of relatively large amounts of cytochromes P450 2B1 and 4B1 in nonciliated bronchiolar epithelial cells (Clara cells) of the terminal airways. These P450 enzymes are highly proficient in vitro at converting ipomeanol to reactive products. Lung tissue contains other enzymic systems which are capable of catalyzing phase I biostransformation pathways (e.g., flavin-containing amine monooxygenase, amine oxidase, and prostaglandin synthase). Examples, however, where pulmonary metabolism by these pathways results in lung toxicity are less numerous than with P450 mediated reactions. Pulmonary prostaglandin H-synthase mediated cooxygenation has been shown to activate procarcinogens such as benzo(a)pyrene 7,8-dihydrodiol, aflatoxin B1, and monosubstituted hydrazines. The activities of pulmonary phase II (conjugation) pathways may also contribute to lung toxicity. Low glutathione transferase activity (relative to P450 mediated aryl hydrocarbon hydroxylase activity) in lung tissue has been suggested to correlate with elevated risk of lung cancer in smokers. Other examples of lung-specific toxic agents and possible causative roles of biotransformation are also discussed.

Sequence analysis of novel CYP4 transcripts from Mytilus galloprovincialis

Cytochrome P450 enzymes (CYPs) are essential components of cellular detoxification system. We identified and characterized seven new cytochrome P450 gene transcript clusters in the populations of bivalve mollusc Mytilus galloprovincialis from three different locations. The phylogenetic analysis identified all transcripts as clusters within the CYP4 branch. Identified clusters, each comprising a number of transcript variants, were designated CYP4Y1, Y2, Y3, Y4, Y5, Y6 and Y7. Transcript clusters CYP4Y2 and Y7, and CYP4Y5 and Y6 showed site specificity, while the transcript clusters CYP4Y1, Y3 and Y4 were present at all investigated locations. The comparison of transcripts deduced amino acid sequences with CYP4s from vertebrate and invertebrate species showed high conservation of the residues and domains essential to the putative function of the enzyme, as terminal ω-hydroxylation and prostaglandin hydroxylation. Our results suggest the great expansion of the CYP4Y cDNAs indicative of CYP4 proteins in the mussel M. galloprovincialis presumably as a response to different environmental conditions.

Metabolism of agrochemicals and related environmental chemicals based on cytochrome P450s in mammals and plants

A yeast gene expression system originally established for mammalian cytochrome P450 monooxygenase cDNAs was applied to functional analysis of a number of mammalian and plant P450 species, including 11 human P450 species (CYP1A1, CYP1A2, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C18, CYP2C19, CYP2D6, CYP2E1 and CYP3A4). The human P450 species CYP1A1, CYP1A2, CYP2B6, CYP2C18 and CYP2C19 were identified as P450 species metabolising various agrochemicals and environmental chemicals. CYP2C9 and CYP2E1 specifically metabolised sulfonylurea herbicides and halogenated hydrocarbons respectively. Plant P450 species metabolising phenylurea and sulfonylurea herbicides were also identified mainly as the CYP71 family, although CYP76B1, CYP81B1 and CYP81B2 metabolised phenylurea herbicides. The transgenic plants expressing these mammalian and plant P450 species were applied to herbicide tolerance as well as phytoremediation of agrochemical and environmental chemical residues. The combined use of CYP1A1, CYP2B6 and CYP2C19 belonging to two families and three subfamilies covered a wide variety of herbicide tolerance and phytoremediation of these residues. The use of 2,4-D-and bromoxynil-induced CYP71AH11 in tobacco seemed to enhance herbicide tolerance and selectivity.

Cytochrome P450-2D6*4 polymorphism seminal relationship in infertile men

This study aimed to assess cytochrome (CY) P450-2D6*4 polymorphism relationship with semen variables in infertile men. In all, 308 men were included; fertile normozoospermia (N) (n = 77), asthenozoospermia (A) (n = 70), asthenoteratozoospermia (AT) (n = 75) and oligoasthenoteratozoospermia (OAT) (n = 86). They were subjected to history taking, clinical examination, semen analysis, sperm acrosin activity, seminal malondialdehyde (MDA) and CYP450-2D6*4 genotyping. CYP450-2D6*4 wild-type allele was represented in 76.5% of N, 70% of A, 66.7% of AT and 57.7% of OAT men where homozygous gene mutation was present in 5.9% of N, 20% of A, 26.6% of AT and 26.9% of OAT men, respectively. Sperm acrosin activity, sperm concentration, sperm motility, linear sperm velocity and sperm normal forms were significantly higher, and seminal MDA level was significantly lower in men with CYP450-2D6*4 wild-type allele compared with men with homozygous mutation. It is concluded that CYP450-2D6*4 wild-type allele has higher frequency where homozygous-type allele has lower frequency in N men compared with A, AT and OAT men. Sperm acrosin activity index, sperm concentration, sperm motility, linear sperm velocity and sperm normal forms were significantly higher, and seminal MDA level was significantly lower in men with CYP450-2D6*4 wild-type allele compared with men with homozygous mutation.

Similarity in recombination rate and linkage disequilibrium at CYP2C and CYP2D cytochrome P450 gene regions among Europeans indicates signs of selection and no advantage of using tagSNPs in population isolates

OBJECTIVE

Linkage disequilibrium (LD) and recombination rate variations are known to vary considerably between human genome regions and populations mostly because of the combined effects of mutation, recombination, and demographic history. Thus, the pattern of LD is a key issue to disentangle variants associated with complex traits. Here, we aim to describe the haplotype structure and LD variation at the pharmacogenetically relevant cytochrome P450 CYP2C and CYP2D gene regions among European populations.

METHODS

To assess the haplotype structure, LD pattern, and recombination rate variations in the clinically significant CYP2C and CYP2D regions, we genotyped 143 single-nucleotide polymorphisms (SNPs) across these two genome regions in a diverse set of 11 European population samples and one sub-Saharan African sample.

RESULTS

Our results showed extended patterns of LD and in general a low rate of recombination at these loci, and a low degree of allele differentiation for the two cytochrome P450 regions among Europeans, with the exception of the Sami and the Finns as European outliers. The Sami sample showed reduced haplotype diversity and higher LD for the two cytochrome P450 regions than the other Europeans, a feature that is proposed to enhance the LD mapping of underlying common complex traits. However, recombination hotspots and LD blocks at these two regions showed highly consistent structures across Europeans including Finns and Sami. Moreover, we showed that the CEPH sample has significantly higher tag transferability among Europeans and a more efficient tagging of both the rare CYP2C9 and the common CYP2C19 functional variants than the Sami. Our data set included CYP2C9*3 (rs1057910) and CYP2C19*2 (rs4244285) enzyme activity-altering variants associated in a recent genome-wide study with acenocoumarol-induced and warfarin-induced anticoagulation or to the antiplatelet effect of clopidogrel, respectively. Including these known activity-altering variants, we showed the haplotype variation and high derived allele frequencies of novel recently identified acenocoumarol genome-wide associated SNPs at CYP2C9 (rs4086116) and CYP2C18 (rs12772169, rs1998591, rs2104543, rs1042194) loci in a comprehensive set of 11 European populations. Furthermore, a significant frequency difference of a CYP2C19*2 gene mutation causing variable drug reactions was observed among Europeans.

CONCLUSION

The CEPH sample representing the general European population as such in the HapMap project seems to be the optimal population sample for the LD mapping of common complex traits among Europeans. Nevertheless, it is still argued that the unique pattern of LD in the Sami may offer an advantage for further association mapping, especially if multiple rare variants play a role in disease etiology. However, besides the activity-altering CYP2C9*3 (rs1057910) and CYP2C19*2 (rs4244285) variants, the high derived allele frequencies of novel recently identified acenocoumarol genome-wide associated SNPs at CYP2C9 (rs4086116) and CYP2C18 (rs12772169, rs1998591, rs2104543, rs1042194) loci variants indicated that the CYP2C region may have been influenced by selection. Thus, this fine-scale haplotype map of the CYP2C and CYP2D regions may help to choose markers for further association mapping of complex pharmacogenetic traits at these loci.

Phytomonitoring and phytoremediation of agrochemicals and related compounds based on recombinant cytochrome P450s and aryl hydrocarbon receptors (AhRs)

Molecular mechanisms of metabolism and modes of actions of agrochemicals and related compounds are important for understanding selective toxicity, biodegradability, and monitoring of biological effects on nontarget organisms. It is well-known that in mammals, cytochrome P450 (P450 or CYP) monooxygenases metabolize lipophilic foreign compounds. These P450 species are inducible, and both CYP1A1 and CYP1A2 are induced by aryl hydrocarbon receptor (AhR) combined with a ligand. Gene engineering of P450 and NADPH cytochrome P450 oxidoreductase (P450 reductase) was established for bioconversion. Also, gene modification of AhRs was developed for recombinant AhR-mediated β-glucronidase (GUS) reporter assay of AhR ligands. Recombinant P450 genes were transformed into plants for phytoremediation, and recombinant AhR-mediated GUS reporter gene expression systems were each transformed into plants for phytomonitoring. Transgenic rice plants carrying CYP2B6 metabolized the herbicide metolachlor and remarkably reduced the residues in the plants and soils under paddy field conditions. Transgenic Arabidopsis plants carrying recombinant guinea pig (g) AhR-mediated GUS reporter genes detected PCB126 at the level of 10 ng/g soils in the presence of biosurfactants MEL-B. Both phytomonitoring and phytoremediation plants were each evaluated from the standpoint of practical uses.

Cytochrome P450 diversity and induction by gorgonian allelochemicals in the marine gastropod Cyphoma gibbosum

BACKGROUND

Intense consumer pressure strongly affects the structural organization and function of marine ecosystems, while also having a profound effect on the phenotype of both predator and prey. Allelochemicals produced by prey often render their tissues unpalatable or toxic to a majority of potential consumers, yet some marine consumers have evolved resistance to host chemical defenses. A key challenge facing marine ecologists seeking to explain the vast differences in consumer tolerance of dietary allelochemicals is understanding the biochemical and molecular mechanisms underlying diet choice. The ability of marine consumers to tolerate toxin-laden prey may involve the cooperative action of biotransformation enzymes, including the inducible cytochrome P450s (CYPs), which have received little attention in marine invertebrates despite the importance of allelochemicals in their evolution.

RESULTS

Here, we investigated the diversity, transcriptional response, and enzymatic activity of CYPs possibly involved in allelochemical detoxification in the generalist gastropod Cyphoma gibbosum, which feeds exclusively on chemically defended gorgonians. Twelve new genes in CYP family 4 were identified from the digestive gland of C. gibbosum. Laboratory-based feeding studies demonstrated a 2.7- to 5.1-fold induction of Cyphoma CYP4BK and CYP4BL transcripts following dietary exposure to the gorgonian Plexaura homomalla, which contains high concentrations of anti-predatory prostaglandins. Phylogenetic analysis revealed that C. gibbosum CYP4BK and CYP4BL were most closely related to vertebrate CYP4A and CYP4F, which metabolize pathophysiologically important fatty acids, including prostaglandins. Experiments involving heterologous expression of selected allelochemically-responsive C. gibbosum CYP4s indicated a possible role of one or more CYP4BL forms in eicosanoid metabolism. Sequence analysis further demonstrated that Cyphoma CYP4BK/4BL and vertebrate CYP4A/4F forms share identical amino acid residues at key positions within fatty acid substrate recognition sites.

CONCLUSIONS

These results demonstrate differential regulation of CYP transcripts in a marine consumer feeding on an allelochemical-rich diet, and significantly advance our understanding of both the adaptive molecular mechanisms that marine consumers use to cope with environmental chemical pressures and the evolutionary history of allelochemical-metabolizing enzymes in the CYP superfamily.

The cytochrome p450 homepage

The Cytochrome P450 Homepage is a universal resource for nomenclature and sequence information on cytochrome P450 (CYP) genes. The site has been in continuous operation since February 1995. Currently, naming information for 11,512 CYPs are available on the web pages. The P450 sequences are manually curated by David Nelson, and the nomenclature system conforms to an evolutionary scheme such that members of CYP families and subfamilies share common ancestors. The organisation and content of the Homepage are described.

Structural and biochemical evidence that a TEM-1 beta-lactamase N170G active site mutant acts via substrate-assisted catalysis

TEM-1 beta-lactamase is the most common plasmid-encoded beta-lactamase in Gram-negative bacteria and is a model class A enzyme. The active site of class A beta-lactamases share several conserved residues including Ser(70), Glu(166), and Asn(170) that coordinate a hydrolytic water involved in deacylation. Unlike Ser(70) and Glu(166), the functional significance of residue Asn(170) is not well understood even though it forms hydrogen bonds with both Glu(166) and the hydrolytic water. The goal of this study was to examine the importance of Asn(170) for catalysis and substrate specificity of beta-lactam antibiotic hydrolysis. The codon for position 170 was randomized to create a library containing all 20 possible amino acids. The random library was introduced into Escherichia coli, and functional clones were selected on agar plates containing ampicillin. DNA sequencing of the functional clones revealed that only asparagine (wild type) and glycine at this position are consistent with wild-type function. The determination of kinetic parameters for several substrates revealed that the N170G mutant is very efficient at hydrolyzing substrates that contain a primary amine in the antibiotic R-group that would be close to the Asn(170) side chain in the acyl-intermediate. In addition, the x-ray structure of the N170G enzyme indicated that the position of an active site water important for deacylation is altered compared with the wild-type enzyme. Taken together, the results suggest the N170G TEM-1 enzyme hydrolyzes ampicillin efficiently because of substrate-assisted catalysis where the primary amine of the ampicillin R-group positions the hydrolytic water and allows for efficient deacylation.

Human ATP-binding cassette (ABC) transporter family

There exist four fundamentally different classes of membrane-bound transport proteins: ion channels; transporters; aquaporins; and ATP-powered pumps. ATP-binding cassette (ABC) transporters are an example of ATP-dependent pumps. ABC transporters are ubiquitous membrane-bound proteins, present in all prokaryotes, as well as plants, fungi, yeast and animals. These pumps can move substrates in (influx) or out (efflux) of cells. In mammals, ABC transporters are expressed predominantly in the liver, intestine, blood-brain barrier, blood-testis barrier, placenta and kidney. ABC proteins transport a number of endogenous substrates, including inorganic anions, metal ions, peptides, amino acids, sugars and a large number of hydrophobic compounds and metabolites across the plasma membrane, and also across intracellular membranes. The human genome contains 49 ABC genes, arranged in eight subfamilies and named via divergent evolution. That ABC genes are important is underscored by the fact that mutations in at least I I of these genes are already known to cause severe inherited diseases (eg cystic fibrosis and X-linked adrenoleukodystrophy [X-ALD]). ABC transporters also participate in the movement of most drugs and their metabolites across cell surface and cellular organelle membranes; thus, defects in these genes can be important in terms of cancer therapy, pharmacokinetics and innumerable pharmacogenetic disorders.

Update on the olfactory receptor (OR) gene superfamily

The olfactory receptor gene (OR) superfamily is the largest in the human genome. The superfamily contains 390 putatively functional genes and 465 pseudogenes arranged into 18 gene families and 300 subfamilies. Even members within the same subfamily are often located on different chromosomes. OR genes are located on all autosomes except chromosome 20, plus the X chromosome but not the Y chromosome. The gene:pseudogene ratio is lowest in human, higher in chimpanzee and highest in rat and mouse -- most likely reflecting the greater need of olfaction for survival in the rodent than in the human. The OR genes undergo allelic exclusion, each sensory neurone expressing usually only one odourant receptor allele; the mechanism by which this phenomenon is regulated is not yet understood. The nomenclature system (based on evolutionary divergence of genes into families and subfamilies of the OR gene superfamily) has been designed similarly to that originally used for the CYP gene superfamily.

Versatile capacity of shuffled cytochrome P450s for dye production

DNA family shuffling is a relatively new method of directed evolution used to create novel enzymes in order to improve their existing properties or to develop new features. This method of evolution in vitro has one basic requirement: a high similarity of initial parental sequences. Cytochrome P450 enzymes are relatively well conserved in their amino acid sequences. Members of the same family can have more than 40% of sequence identity at the protein level and are therefore good candidates for DNA family shuffling. These xenobiotic-metabolising enzymes have an ability to metabolise a wide range of chemicals and produce a variety of products including blue pigments such as indigo. By applying the specifically designed DNA family shuffling approach, catalytic properties of cytochrome P450 enzymes were further extended in the chimeric progeny to include a new range of blue colour formations. This mini-review evokes the possibility of exploiting directed evolution of cytochrome P450s and the novel enzymes created by DNA family shuffling for the production of new dyes.

Methylglyoxal and glucose metabolism: a historical perspective and future avenues for research

Methylglyoxal, an alpha-oxoaldehyde discovered in the 1880s, has had a hectic scientific career, at times being considered of fundamental importance and at other times viewed as playing a very subordinate role. Much has been learned about methylglyoxal, but the function of its production in the metabolic machinery is still unknown. This paper gives an overview of the changing role of methylglyoxal from a historical aspect and arrives at the conclusion that methylglyoxal is tightly bound to glycolysis from an evolutionary perspective, its production therefore being inevitable. It is not situated in the main stream of the glycolytic sequence, but a role can be assigned to its production in the phosphate supply of operating glycolysis in some prokaryotes and yeast under conditions of phosphate deficiency. This function is presumed to be performed by the enzyme methylglyoxal synthase, which is specialized for the conversion of dihydroxyacetone-phosphate to methylglyoxal. However, it is still unknown whether this enzyme and this kind of regulation also exist in animals.

Hormone receptor-related gene polymorphisms and prostate cancer risk in North Indian population

The purpose of this study was to analyse the frequency and type of mutations in the coding region of androgen receptor (AR) and to determine the role of polymorphisms in the intron 1 of ERalpha, exon 5 of ERbeta, intron 7 of progesterone, exon 7 of the aromatase (CYP19) and exon 9 of VDR genes in the risk of prostate cancer. PCR-RFLP analysis of all above the genes was on 100 prostate cancer patients and an equal number of matching controls. The study also included PCR-SSCP analyses of exons 2-8 of AR gene. The genotype containing -/- allele of ERalpha gene was statistically significant for the risk of prostate cancer pose (OR, 2.70; 95% CI, 1.08-6.70, P = 0.032) Rr genotype of ERbeta gene also have a higher risk (OR, 1.65; 95% CI, 0.52-5.23) for prostate cancer. The Cys allele of CYP19 gene was also associated with statistically significant increased risk of prostate cancer (OR; 2.28, 95% CI, 1.20-4.35, P = 0.012). tt genotype of codon 352 of VDR gene showed an OR of 0.43 for (95% CI, 0.13-1.39) and an OR for Tt genotype was 0.65 (95% CI, 0.36-1.16). Taken together, the results showed that in North Indian population, ERalpha and CYP19 genes may be playing a role in the risk of prostate cancer.

Ontogenesis of brain aromatase P450 expression in the bovine hypothalamus

Aromatase P450 (P450(AROM)), converting testosterone (T) into estradiol (E), plays an important role in sexual differentiation of neural structures in the developing mammalian brain. The aim of the present study was to characterize the qualitative and quantitative profile of P450(AROM) mRNA expression in the bovine hypothalamus (the region of the central nervous system in which the enzyme is mainly localized) using RT-PCR and quantitative real-time RT-PCR analysis, respectively. P450(AROM) expression was examined in the developing hypothalamus in a series of experimental groups investigated at 10 weeks interval one from the other. Our data indicate that in the bovine fetal hypothalamus P450(AROM) expression peaks at the second quarter of gestation. The presence of neural cells containing P450(AROM) in the bovine fetal hypothalamus was confirmed by immunohistochemistry, and localized in the medial preoptic area. We conclude that second quarter of the gestation is the developmental stage which represents a critical period for hypothalamic differentiation in bovine ontogenesis, an important difference with the rat and mouse, short gestation species in which P450(AROM) activity peaks around delivery.

Inter-individual variation in the metabolic activation of the antimalarial biguanides

The aryl-biguanides proguanil and chlorproguanil were developed as part of a collaborative programme between ICI and the Liverpool School of Tropical Medicine during the 1940s. The compounds were characterized by their absence of host toxicity. However, the rapid development of parasite resistance to the actions of these drugs and the development of the 4-aminoquinoline, chloroquine, severely limited their use. The subsequent widespread development of parasite resistance to chloroquine, together with the observations that the magnitude of dihydrofolate reductase inhibitor resistance (the site of action of the biguanides) developed to pyrimethamine is not directly correlated with biguanide resistance(1,2). has resulted in renewed interest in these drugs. In particular, proguanil is now the drug of choice for malaria prophylaxis, in combination with chloroquine; used in combination with a suitable sulphonamide, it may be of value in malaria therapy.

Lack of Interaction of Milnacipran with the Cytochrome P450 Isoenzymes Frequently Involved in the Metabolism of Antidepressants

OBJECTIVE

To compare the pharmacokinetics of milnacipran in extensive metabolisers (EMs) and poor metabolisers (PMs) of sparteine and mephenytoin, and to assess the influence of multiple administrations of milnacipran on the activity of cytochrome P450 (CYP) isoenzymes through its own metabolism and through various probes, namely CYP2D6 (sparteine/dextromethorphan), CYP2C19 (mephenytoin), CYP1A2 (caffeine) and CYP3A4 (endogenous 6-beta-hydroxy-cortisol excretion).

METHODS

Twenty-five healthy subjects, 12 EMs for both sparteine/dextromethorphan and mephenytoin, nine EMs for mephenytoin and PMs for sparteine/dextromethorphan (PM(2D6)) and four PMs for mephenytoin and EMs for sparteine/dextromethorphan (PM(2C19)) were administered milnacipran as a single 50 mg capsule on day 1 followed by a 50 mg capsule twice daily for 7 days. The pharmacokinetics of milnacipran and its oxidative metabolites were assessed after the first dose (day 1) and after multiple administration (day 8), and were compared for differences between CYP2D6 and CYP2C19 PMs and EMs. Metabolic tests were performed before (day -2), during (days 1 and 8) and after (day 20) milnacipran administration.

RESULTS

Milnacipran steady state was rapidly achieved. Metabolism was limited: approximately 50% unchanged drug, 30% as glucuronide and 20% as oxidative metabolite (mainly F2800 the N-dealkyl metabolite). Milnacipran administration to PM2D6 and PM2C19 subjects did not increase parent drug exposure or decrease metabolite exposure. Milnacipran oxidative metabolism is not mediated through CYP2D6 or CYP2C19 polymorphic pathways nor does it significantly interact with CYP1A2, CYP2C19, CYP2D6 or CYP3A4 activities.

CONCLUSION

Limited reciprocal pharmacokinetic interaction between milnacipran and CYP isoenzymes would confer flexibility in the therapeutic use of the drug when combined with antidepressants. Drug-drug interaction risk would be low, even if the combined treatments were likely to inhibit CYP2D6 and CYP2C19 isoenzyme activities.

Mechanism of impairment of cytochrome P450-dependent metabolism in hamster liver during leishmaniasis

The mechanism of impairment of cytochrome P450 (P450)-dependent metabolism in hamster liver during leishmaniasis is reported. A significant decrease in the level of P450 was observed on the 20th day of infection when the parasite load in the liver was maximum. The decrease in P450 level was accompanied by a significant increase in the level of marker enzymes of liver and degeneration of liver tissue. The impairment was isozyme-specific and concomitant with the induction of nitric oxide synthase. The results of in vitro experiments with generated nitric oxide and with scavengers demonstrated that the impairment is mediated by NO. Treatment of the infected animals with a combination therapy showed reduction in parasite load, reversal of P450 impairment, and recovery of liver enzymes and tissue close to normal.

Effect of cyclophosphamide on serum cyclosporine levels at the conditioning of hematopoietic stem cell transplantation

We retrospectively analyzed the factors that affect serum cyclosporine (CsA) concentrations up to day 14 after allogeneic hematopoietic stem cell transplantation (HSCT). In all, 103 transplant recipients who received MTX and CsA for acute GVHD prophylaxis were analyzed. No significant relationships between serum CsA concentrations and gender, age, serum creatinine levels, AST/ALT levels, or antibiotic/fluconazole administration were found by comparing median CsA concentrations or by using longitudinal or regression multivariate analyses. However, the mean of the median serum CsA concentration in patients (n=54) receiving the regimen containing cyclophosphamide (CY) (149.7 ng/ml; 95% confidence interval (CI): 132.1-167.4) was significantly (P<0.0001) lower than that in patients (n=49) receiving the non-CY regimen (217.3 ng/ml; 95% CI: 198.9-235.6). Longitudinal analysis and regression multivariate analysis showed that only administration of CY had a significant effect on the serum CsA concentration. Our results suggest that administration of CY during conditioning can reduce the effects on serum CsA concentrations during the 2 weeks following HSCT. The mechanism of this effect is not clear, but it may be due to the autoinduction of CY.

Differentiation-dependent induction of CYP1A1 in cultured rat small intestinal epithelial cells, colonocytes, and human colon carcinoma cells: basement membrane-mediated apoptosis

Rat small intestinal epithelial cells and human colon adenocarcinoma cells cultured on Matrigel expressed the differentiation specific enzyme, sucrase-isomaltase, as determined by indirect immunofluorescence. Rat small intestinal epithelial cells, rat colonocytes, and human colon adenocarcinoma cells developed an altered morphology when cultured on Matrigel and became apoptotic within 24-48 h. Benzo[a]pyrene and 2,3,7,8-tetrachlorodibenzo-p-dioxin caused a 2- and 5-fold induction, respectively, of ethoxyresorufin-o-deethylase activity in rat small intestinal epithelial cells cultured on Matrigel. Benzo[a]pyrene- or 2,3,7,8-tetrachlorodibenzo-p-dioxin-induced ethoxyresorufin-o-deethylase activity in rat small intestinal epithelial cells cultured on plastic was not detected. 2,3,7,8-tetrachlorodibenzo-p-dioxin treatment caused a 14-fold induction of transfected, rat CYP1A1-promoter-luciferase activity in rat small intestinal epithelial cells cultured on Matrigel. Benzo[a]pyrene and 2,3,7,8-tetrachlorodibenzo-p-dioxin treatment induced ethoxyresorufin-o-deethylase activity by 6- and 1.6-fold, respectively in rat colonocytes cultured on Matrigel. Induction of ethoxyresorufin-o-deethylase activity was not observed in rat colonocytes cultured on plastic. CYP1A1-promoter-luciferase activity was induced 3-fold by 2,3,7,8-tetrachlorodibenzo-p-dioxin in rat colonocytes cultured on Matrigel. Induction of CYP1A1-promoter-luciferase activity in rat small intestinal epithelial cells or rat colonocytes cultured on plastic was not observed. Ethoxyresorufin-o-deethylase activity in human colon adenocarcinoma cells, cultured on either plastic or Matrigel, was induced 7-fold by benzo[a]pyrene. 2,3,7,8-Tetrachlorodibenzo-p-dioxin-induced ethoxyresorufin-o-deethylase activity was 2-fold greater in human colon adenocarcinoma cells cultured on Matrigel compared to cells cultured on plastic. Extracellular matrix-mediated differentiation and apoptosis of intestinal cells provide in vitro systems for study of the regulation of CYP1A1 expression, carcinogen activation in the gut and mechanism(s) of apoptosis of colon cancer cells.

Molecular Cloning of Monkey CYP2C43 cDNA and Expression in Yeast

A cDNA clone designated as CYP2C43 was isolated from the rhesus monkey liver cDNA library. The first 16 amino acid residues at the N-terminal region of this cDNA product were identical with those of P450 CMLd which have been purified and characterized as S-mephenytoin 4'-hydroxylase in monkey liver. The respective nucleotide and deduced amino acid sequences of CYP2C43 were 83% and 77%, identical to those of monkey CYP2C20. Antibody against CYP2C9 detected a protein in the microsomes of yeast transformed CYP2C43 expression plasmid. The specific content of recombinant CYP2C43 was 78.0 pmol/mg protein and the yield was 4.23 nmol/l of the culture. CYP2C43 was able to metabolize S-mephenytoin stereo-selectively. The activity for S-mephenytoin in the microsomes reconstituted with or without cytochrome b(5) was found to be 96.2 or 23.7 pmol/min/nmol P450, respectively. CYP2C43, however, did not show any oxidative activity for tolbutamide. These results indicate that CYP2C43 is the second identified member of the monkey CYP2C subfamily and a cDNA clone encoding P450 CMLd in monkey.

Genetic Polymorphism of Human Cytochrome P450 Involved in Drug Metabolism

Recent advances in human gene analysis promoted by the human genome project have brought us a massive amount of information. These data can be seen and analyzed by personal computer through individual Web sites. As a result, the best use of bioinformatic is essential for recent molecular biology research. Genetic polymorphism of drug-metabolizing enzymes influences individual drug efficacy and safety through the alteration of pharmacokinetics and disposition of drugs. Considerable amounts of data have now accumulated as allelic differences of various drug metabolizing enzymes. Current understanding of genotype information on cytochrome P450 is hereby summarized, based on the Web site for their use in individual optimization of drug therapy.

Role of LYS271 and LYS279 residues in the interaction of cytochrome P4501A1 with NADPH-cytochrome P450 reductase

It has been proposed that negatively charged amino acids on the surface of reductase and positively charged amino acids on the surface of P450 mediate the binding of both proteins through electrostatic interactions. In this study, we used a site-directed mutagenesis approach to determine a role for two lysine residues (Lys271 and Lys279) of cytochrome P4501A1 in the interaction of P4501A1 with reductase. We prepared two mutants P4501A1Ile271 and P4501A1Ile279 with a mutation of the lysine at positions 271 and 279, respectively. We observed a strong inhibition (>80%) of the 7-ethoxycoumarin and ethoxyresorufin deethylation activity in the reductase-supported system for both mutants. In the cumene hydroperoxide-supported system, P4501A1Ile279 exhibited wild-type activity, but the P4501A1Ile271 mutant activity remained low. The CD spectrum and substrate-binding assay indicated that the secondary structure of P4501A1Ile271 is perturbed. To evaluate further the involvement of these P4501A1 lysine residues in reductase binding, we measured the KM of reductase for wild type and mutants. Both wild type and P4501A1Ile271 reached saturation in the range of reductase concentrations tested with KM values 5.1 and 11.2 pM, respectively. The calculated KM value for P4501A1Ile279 increased 9-fold, 44.4 pM, suggesting that the mutation affected binding of reductase to P4501A1. Stopped-flow spectroscopy was employed to evaluate the effect of mutations on electron transfer from reductase to heme iron. Both wild type and P450Ile279 showed biphasic kinetics with a approximately 40% participation of the fast step in the total activity. On the other hand, only single-phase kinetics for iron reduction was observed for P450Ile271, suggesting that the low activity of this mutant can be attributed not only to major structural changes but also to a disturbance in the electron transport.

The complete human olfactory subgenome

Olfactory receptors likely constitute the largest gene superfamily in the vertebrate genome. Here we present the nearly complete human olfactory subgenome elucidated by mining the genome draft with gene discovery algorithms. Over 900 olfactory receptor genes and pseudogenes (ORs) were identified, two-thirds of which were not annotated previously. The number of extrapolated ORs is in good agreement with previous theoretical predictions. The sequence of at least 63% of the ORs is disrupted by what appears to be a random process of pseudogene formation. ORs constitute 17 gene families, 4 of which contain more than 100 members each. "Fish-like" Class I ORs, previously considered a relic in higher tetrapods, constitute as much as 10% of the human repertoire, all in one large cluster on chromosome 11. Their lower pseudogene fraction suggests a functional significance. ORs are disposed on all human chromosomes except 20 and Y, and nearly 80% are found in clusters of 6-138 genes. A novel comparative cluster analysis was used to trace the evolutionary path that may have led to OR proliferation and diversification throughout the genome. The results of this analysis suggest the following genome expansion history: first, the generation of a "tetrapod-specific" Class II OR cluster on chromosome 11 by local duplication, then a single-step duplication of this cluster to chromosome 1, and finally an avalanche of duplication events out of chromosome 1 to most other chromosomes. The results of the data mining and characterization of ORs can be accessed at the Human Olfactory Receptor Data Exploratorium Web site (http://bioinfo.weizmann.ac.il/HORDE).