Elsevier

Gene

Volume 387, Issues 1–2, 31 January 2007, Pages 58-66
Gene

Structural divergence and adaptive evolution in mammalian cytochromes P450 2C

https://doi.org/10.1016/j.gene.2006.08.017Get rights and content

Abstract

Cytochromes P450 (CYPs) comprise a superfamily of enzymes involved in various physiological functions, including the metabolism of drugs and carcinogenic compounds present in food, making them of great importance for human health. The possibility that CYPs could be broadening or changing substrate specificity in accordance to the high diversity of xenobiotics compounds environmentally available suggests that their metabolic function could be under adaptive evolution. We evaluated the existence of functional divergence and signatures of selection on mammalian genes from the drug-metabolizing CYP2 family. Thirteen of the sites found to be functionally divergent and the eight found to be under strong positive selection occurred in important functional domains, namely on the substrate entrance channel and within the active site. Our results provide insight into CYPs evolution and the role of molecular adaptation in enzyme substrate-specificity diversification.

Introduction

Cytochromes P450 (CYPs) can be found in almost all organisms and present a high intra and interspecies diversity (Werck-Reichhart and Feyereisen, 2000). One of their tasks is the oxidation of xenobiotic compounds to facilitate their excretion from the organism (Omiecinski et al., 1999, Guengerich, 2001). These include drugs and carcinogenic compounds present in food (Anzenbacher and Anzenbacherova, 2001), which explains their importance for pharmaceutical research.

The possibility that the CYPs performing such tasks could be broadening or changing substrate specificity in accordance to environmental changes suggests that their metabolic function could be under adaptive evolution. In fact, mosquito CYPs seem to be involved in insecticide resistance (Ranson et al., 2002). Different evolutionary pressures on particular sites of CYPs should be revealed by the analysis of the rate of substitutions occurring within coding regions.

Protein functional divergence after gene duplication can be tested by calculating the replacement rates that occur in the originated subfamilies (Gaucher et al., 2001). Likelihood-ratio tests (LRTs) can be used to assess whether the rates differ between the subfamilies and/or are accelerated in one/all of them (Knudsen and Miyamoto, 2001, Gu, 2003, Knudsen et al., 2003).

Nonsynonymous substitutions may influence the fitness of an individual or population. Thus, adaptive molecular evolution may cause the nonsynonymous substitution rate (dN) to be higher than the synonymous rate (dS), with the ratio ω (dN/dS) being higher than 1 (Yang et al., 2000). The methods that provide statistical measures of such mutation rates fail, however, in evaluating the physicochemical importance of the correspondent amino acid changes and its consequences for the protein function (McClellan et al., 2005). Further analyses are necessary to determine if there is a statistically significant change in the amino acid properties at particular sites. CYPs chemical and structural variations in areas that are in contact with the ligands (substrate recognition sites, SRSs (Gotoh, 1992)) will have consequences on the size, shape, and chemical characteristics of both substrates and products.

In this study, we assessed functional divergence among mammalian CYP2C and used two types of methods to determine if these genes are under adaptive evolution: (i) a gene level approach, testing for functional divergence and positive selection using statistical methods and (ii) a protein level approach, evaluating statistically significant physicochemical amino acid changes and verifying the impact of the previous analyses on the protein three-dimensional (3D) structure. So far, positive selection has been mostly associated with protein recognition domains, such as those involved in immune response and reproduction (for a list of examples see Table S1 in supporting information; henceforth all supporting information will be indicated by an S preceding the correspondent numbering). We retrieved unequivocal evidence that diversifying selection is acting on CYP's active site, thus providing insight to understand rapid enzyme function diversification.

Section snippets

Sequences and structures

The CYP superfamily is classified into families and subfamilies (> 40% or 55% amino acid sequence identity, respectively) (Nelson et al., 2004). The enzymes used in this study belong to the largest and most diverse of CYP families, CYP2 (Omiecinski et al., 1999). Its members are responsible for the metabolism of a variety of different pharmaceutical agents and represent more than 20% of the human liver total CYPs content (Omiecinski et al., 1999). Mutations occurring both within and outside the

Selection analyses

The first evidence of selection in CYPs was retrieved by the SWAAP sliding window approach (Fig. 1A). Very low values of dN/dS were observed for the heme binding areas. This is consistent with the fact that these areas are associated with the catalytic oxidative mechanism, which is common to all CYPs (Poulos, 1995). By contrast, the substrate binding areas of the active site presented some of the higher values of dN/dS. These include SRS-1, SRS-3, SRS-2, SRS-5 and SRS-6 (Fig. 1B).

We used

Acknowledgments

We thank the FCT (Fundação para a Ciência e Tecnologia) for a doctoral scholarship (SFRH/BD/7089/2001) for R.F. and the NFCR (National Foundation for Cancer Research) Centre for Drug Discovery, University of Oxford, U.K., for financial support.

References (46)

  • M.A. Fares

    SWAPSC: sliding window analysis procedure to detect selective constraints

    Bioinformatics

    (2004)
  • M.A. Fares et al.

    The evolution of the heat-shock protein GroEL from Buchnera, the primary endosymbiont of aphids, is governed by positive selection

    Mol. Biol. Evol.

    (2002)
  • M.A. Fares et al.

    A sliding window-based method to detect selective constraints in protein-coding genes and its application to RNA viruses

    J. Mol. Evol.

    (2002)
  • N. Galtier et al.

    SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny

    Comput. Appl. Biosci.

    (1996)
  • E.A. Gaucher et al.

    Function-structure analysis of proteins using covarion-based evolutionary approaches: elongation factors

    PNAS

    (2001)
  • X. Gu

    Functional divergence in protein (family) sequence evolution

    Genetica

    (2003)
  • F.P. Guengerich

    Common and uncommon cytochrome P450 reactions related to metabolism and chemical toxicity

    Chem. Res. Toxicol.

    (2001)
  • D.T. Jones et al.

    The rapid generation of mutation data matrices from protein sequences

    Comput. Appl. Biosci.

    (1992)
  • B. Knudsen et al.

    A likelihood ratio test for evolutionary rate shifts and functional divergence among proteins

    Proc. Natl. Acad. Sci. U. S. A.

    (2001)
  • R.M. Laethem et al.

    Identification of rabbit cytochromes P450 2C1 and 2C2 as arachidonic acid epoxygenases

    Molecular

    (1992)
  • W.H. Li

    Unbiased estimation of the rates of synonymous and nonsynonymous substitution

    J. Mol. Evol.

    (1993)
  • D.A. McClellan et al.

    Physicochemical evolution and molecular adaptation of the cetacean and artiodactyl cytochrome b proteins

    Mol. Biol. Evol.

    (2005)
  • J.O. Miners et al.

    Cytochrome P4502C9: an enzyme of major importance in human drug metabolism

    Br. J. Clin. Pharmacol.

    (1998)
  • Cited by (0)

    View full text