Structural divergence and adaptive evolution in mammalian cytochromes P450 2C
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)
Substrate recognition sites in cytochrome P450 family 2 (CYP2) proteins inferred from comparative analyses of amino acid and coding nucleotide sequences
J. Biol. Chem.
(1992)- et al.
Characterization of A 2Nd gene-product related to rabbit cytochrome-P-450-1
J. Biol. Chem.
(1987) - et al.
Evolutionary divergence of thyrotropin receptor structure
Mol. Genet. Metab.
(2004) - et al.
Using evolutionary rates to investigate protein functional divergence and conservation: a case study of the carbonic anhydrases (vol 164, pg 1261, 2002)
Genetics
(2003) - et al.
Cloning and expression of murine CYP2Cs and their ability to metabolize arachidonic acid
Arch. Biochem. Biophys.
(1998) Cytochrome P450
Curr. Opin. Struct. Biol.
(1995)- et al.
Arachidonic-acid metabolism by human cytochrome P450S-2C8, cytochrome P450S-2C9, cytochrome P450S-2E1, and cytochrome P450S-1A2 — regioselective oxygenation and evidence for a role for Cyp2C enzymes in arachidonic-acid epoxygenation in human liver-microsomes
Arch. Biochem. Biophys.
(1995) - et al.
Microsomal cytochrome P450 2C5: comparison to microbial P450s and unique features
J. Inorg. Biochem.
(2000) - et al.
Cytochromes P450 and metabolism of xenobiotics
Cell. Mol. Life Sci.
(2001) - et al.
Cytochrome P450-mediated oxidation of glucuronide derivatives: example of estradiol-17 beta-glucuronide oxidation to 2-hydroxyestradiol-17 beta-glucuronide by CYP2C8
Drug Metab. Dispos.
(2005)