Change In Substrate Specificity Research Articles

The bioenergetic "CK Clamp" technique detects substrate-specific changes in mitochondrial respiration and membrane potential during early VML injury pathology.

Volumetric muscle loss (VML) injuries are characterized by non-recoverable loss of tissue resulting in contractile and metabolic dysfunction. The characterization of metabolic dysfunction in volumetric muscle loss-injured muscle has been interpreted from permeabilized myofiber respiration experiments involving saturating ADP levels and non-physiologic ATP:ADP concentration ratios. The extent to which this testing condition obscures the analysis of mitochondrial (dys) function after volumetric muscle loss injury is unclear. An alternative approach is described that leverages the enzymatic reaction of creatine kinase and phosphocreatine to assess mitochondrial respiration and membrane potential at clamped physiologic ATP:ADP ratios, "CK Clamp." The objective of this study was to validate the CK Clamp in volumetric muscle loss-injured muscle and to detect differences that may exist between volumetric muscle loss-injured and uninjured muscles at 1, 3, 5, 7, 10, and 14days post-injury. Volumetric muscle loss-injured muscle maintains bioenergetic features of the CK Clamp approach, i.e., mitochondrial respiration rate (JO2) titters down and mitochondrial membrane potential is more polarized with increasing ATP:ADP ratios. Pyruvate/malate/succinate-supported JO2 was significantly less in volumetric muscle loss-injured muscle at all timepoints compared to uninjured controls (-26% to -84%, p < 0.001) and electron conductance was less at day 1 (-60%), 5 (-52%), 7 (-35%), 10 (-59%), and 14 (-41%) (p < 0.001). Palmitoyl-carnitine/malate-supported JO2 and electron conductance were less affected following volumetric muscle loss injury. volumetric muscle loss-injury also corresponded with a more polarized mitochondrial membrane potential across the clamped ATP:ADP ratios at day 1 and 10 (pyruvate and palmitoyl-carnitine, respectively) (+5%, p < 0.001). This study supports previous characterizations of metabolic dysfunction and validates the CK Clamp as a tool to investigate bioenergetics in traumatically-injured muscle.

Experimental–theoretical study of laccase as a detoxifier of aflatoxins

We investigate laccase-mediated detoxification of aflatoxins, fungal carcinogenic food contaminants. Our experimental comparison between two aflatoxins with similar structures (AFB1 and AFG2) shows significant differences in laccase-mediated detoxification. A multi-scale modeling approach (Docking, Molecular Dynamics, and Density Functional Theory) identifies the highly substrate-specific changes required to improve laccase detoxifying performance. We employ a large-scale density functional theory-based approach, involving more than 7000 atoms, to identify the amino acid residues that determine the affinity of laccase for aflatoxins. From this study we conclude: (1) AFB1 is more challenging to degrade, to the point of complete degradation stalling; (2) AFG2 is easier to degrade by laccase due to its lack of side products and favorable binding dynamics; and (3) ample opportunities to optimize laccase for aflatoxin degradation exist, especially via mutations leading to π–π stacking. This study identifies a way to optimize laccase for aflatoxin bioremediation and, more generally, contributes to the research efforts aimed at rational enzyme optimization.

Enhanced Substrate Specificity of Thermostable Cytochrome P450 CYP175A1 by Site Saturation Mutation on Tyrosine 68.

Thermostable cytochrome P450 (CYP175A1) cloned from Thermus thermophilus shows mid-point unfolding temperature (Tm) of 88°C (361K) along with high thermodynamic stability making it a potential industrially viable biocatalyst. Molecular docking analyses, and structural superposition with steroidogenic and fatty acid metabolizing cytochrome P450 s suggested that the tyrosine 68 may have important role in binding as well as metabolism of substrates by the enzyme. Site-saturation mutation of the tyrosine 68 residue was carried out and several unique mutations were obtained that were properly folded and showed high thermostability. We investigated the effects of variation of the single residue, Tyr68 at the substrate binding pocket of the enzyme on the substrate specificity of CYP175A1. Screening of the mutant colonies of CYP175A1 obtained after saturation mutagenesis of Tyr68 using saturated fatty acid, myristic acid and poly unsaturated fatty acids showed that the Y68K had notable binding and catalytic activity for mono-oxygenation of the saturated fatty acid (myristic acid), which had no major detectable binding affinity towards the WT enzyme. The Y68R mutant of CYP175A1, on the other hand was found to selectively bind and catalyse reaction of cholesterol. The wild type as well as both the mutants of the enzyme however bind poly unsaturated fatty acids. The results thus show that saturation mutation of a single amino acid at the substrate binding pocket of the thermostable cytochrome P450 could induce sufficient changes in the substrate binding pocket of the enzyme that can efficiently change substrate specificity of the enzyme.

Developmental programming of mitochondrial substrate metabolism in skeletal muscle of adult sheep by cortisol exposure before birth.

Prenatal glucocorticoid overexposure causes adult metabolic dysfunction in several species but its effects on adult mitochondrial function remain largely unknown. Using respirometry, this study examined mitochondrial substrate metabolism of fetal and adult ovine biceps femoris (BF) and semitendinosus (ST) muscles after cortisol infusion before birth. Physiological increases in fetal cortisol concentrations pre-term induced muscle- and substrate-specific changes in mitochondrial oxidative phosphorylation capacity in adulthood. These changes were accompanied by muscle-specific alterations in protein content, fibre composition and abundance of the mitochondrial electron transfer system (ETS) complexes. In adult ST, respiration using palmitoyl-carnitine and malate was increased after fetal cortisol treatment but not with other substrate combinations. There were also significant increases in protein content and reductions in the abundance of all four ETS complexes, but not ATP synthase, in the ST of adults receiving cortisol prenatally. In adult BF, intrauterine cortisol treatment had no effect on protein content, respiratory rates, ETS complex abundances or ATP synthase. Activity of citrate synthase, a marker of mitochondrial content, was unaffected by intrauterine treatment in both adult muscles. In the ST but not BF, respiratory rates using all substrate combinations were significantly lower in the adults than fetuses, predominantly in the saline-infused controls. The ontogenic and cortisol-induced changes in mitochondrial function were, therefore, more pronounced in the ST than BF muscle. Collectively, the results show that fetal cortisol overexposure programmes mitochondrial substrate metabolism in specific adult muscles with potential consequences for adult metabolism and energetics.

Atomistic insight on structure and dynamics of spinach acyl carrier protein with substrate length

The plant acyl-acyl carrier protein (ACP) desaturases are a family of soluble enzymes that convert saturated fatty acyl-ACPs into their cis-monounsaturated equivalents in an oxygen-dependent reaction. These enzymes play a key role in biosynthesis of monounsaturated fatty acids in plants. ACPs are central proteins in fatty acid biosynthesis that deliver acyl chains to desaturases. They have been reported to show a varying degree of local dynamics and structural variability depending on the acyl chain size. It has been suggested that substrate-specific changes in ACP structure and dynamics have a crucial impact on the desaturase enzymatic activity. Using molecular dynamics simulations, we investigated the intrinsic solution structure and dynamics of ACP from spinach with four different acyl chains: capric (C10), myristic (C14), palmitic (C16), and stearic (C18) acids. We found that the fatty acids can adopt two distinct structural binding motifs, which feature different binding free energies and influence the ACP dynamics in a different manner. Docking simulations of ACP to castor Δ9-desaturase and ivy Δ4-desaturase suggest that ACP desaturase interactions could lead to a preferential selection between the motifs.

Molecular Design and Regulation of Metalloenzyme Activities through Two Novel Approaches: Ferritin and P450s

Abstract We have developed two novel approaches for the construction of artificial metalloenzymes showing either unique catalytic activities or substrate specificity. The first example is the use of a hollow cage of apo-ferritin as a reaction vessel for hydrogenation of olefins, Suzuki-Miyaura C-C coupling and phenylacetylene polymerization by employing Pd0 nano-clusters, Pd2+(η3-C3H5) complexes and Rh1+(nbd) (nbd = norbornadiene) complexes introduced in the hollow cage, respectively. The second approach is the use of “decoy molecules” to change substrate specificity of P450s, allowing epoxidation and hydroxylation activities toward nonnative organic substrates in P450SPα, P450BSβ and P450BM3 without the mutation of any amino acid. Finally, the decoy strategy has been applied to an in vivo system of P450, i.e., the use of P450BM3 expressed in the whole cell of E. coli to oxidize benzene to phenol.

APOBEC3A Loop 1 Is a Determinant for Single-Stranded DNA Binding and Deamination.

The apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3 (APOBEC3 or A3) family of proteins functions in the innate immune system. The A3 proteins are interferon inducible and hypermutate deoxycytidine to deoxyuridine in foreign single-stranded DNA (ssDNA). However, this deaminase activity cannot discriminate between foreign and host ssDNA at the biochemical level, which presents a significant danger when A3 proteins gain access to the nucleus. Interestingly, this A3 capability can be harnessed when coupled with novel CRISPR-Cas9 proteins to create a targeted base editor. Specifically, A3A has been used in vitro to revert mutations associated with disease states. Recent structural studies have shown the importance of loop regions of A3A and A3G in ssDNA recognition and positioning for deamination. In this work, we further examined loop 1 of A3A to determine how it affects substrate selection, as well as the efficiency of deamination, in the hopes of advancing the potential of A3A in base editing technology. We found that mutating residue H29 enhanced deamination activity without changing substrate specificity. Also interestingly, we found that increasing the length of loop 1 decreases substrate specificity. Overall, these results lead to a better understanding of substrate recognition and deamination by A3A and the A3 family of proteins.

Changing substrate specificity and iteration of amino acid chain elongation in glucosinolate biosynthesis through targeted mutagenesis of Arabidopsis methylthioalkylmalate synthase 1.

Methylthioalkylmalate synthases catalyse the committing step of amino acid chain elongation in glucosinolate biosynthesis. As such, this group of enzymes plays an important role in determining the glucosinolate composition of Brassicaceae species, including Arabidopsis thaliana. Based on protein structure modelling of MAM1 from A. thaliana and analysis of 57 MAM sequences from Brassicaceae species, we identified four polymorphic residues likely to interact with the 2-oxo acid substrate. Through site-directed mutagenesis, the natural variation in these residues and the effect on product composition were investigated. Fifteen MAM1 variants as well as the native MAM1 and MAM3 from A. thaliana were characterised by heterologous expression of the glucosinolate chain elongation pathway in Escherichia coli. Detected products derived from leucine, methionine or phenylalanine were elongated with up to six methylene groups. Product profile and accumulation were changed in 14 of the variants, demonstrating the relevance of the identified residues. The majority of the single amino acid substitutions decreased the length of methionine-derived products, while approximately half of the substitutions increased the phenylalanine-derived products. Combining two substitutions enabled the MAM1 variant to increase the number of elongation rounds of methionine from three to four. Notably, characterisation of the native MAMs indicated that MAM1 and not MAM3 is responsible for homophenylalanine production. This hypothesis was confirmed by glucosinolate analysis in mam1 and mam3 mutants of A. thaliana.

Conversion of the sensor kinase DcuS of Escherichia coli of the DcuB/DcuS sensor complex to the C4 -dicarboxylate responsive form by the transporter DcuB.

The sensor kinase DcuS of Escherichia coli co-operates under aerobic conditions with the C4 -dicarboxylate transporter DctA to form the DctA/DcuS sensor complex. Under anaerobic conditions C4 -dicarboxylate transport in fumarate respiration is catalyzed by C4 -dicarboxylate/fumarate antiporter DcuB. (i) DcuB interacted with DcuS as demonstrated by a bacterial two-hybrid system (BACTH) and by co-chromatography of the solubilized membrane-proteins (mHPINE assay). (ii) In the DcuB/DcuS complex only DcuS served as the sensor since mutations in the substrate site of DcuS changed substrate specificity of sensing, and substrates maleate or 3-nitropropionate induced DcuS response without affecting the fumarate site of DcuB. (iii) The half-maximal concentration for induction of DcuS by fumarate (1 to 2 mM) and the corresponding Km for transport (50 µM) differ by a factor of 20 to 40. Therefore, the fumarate sites are different in transport and sensing. (iv) Increasing levels of DcuB converted DcuS from the permanent ON (DcuB deficient) state to the fumarate responsive form. Overall, the data show that DcuS and DcuB form a DcuB/DcuS complex representing the C4 -dicarboxylate responsive form, and that the sensory site of the complex is located in DcuS whereas DcuB is required for converting DcuS to the sensory competent state.

Crystal structure and identification of a key amino acid for glucose tolerance, substrate specificity, and transglycosylation activity of metagenomic β-glucosidase Td2F2.

β-Glucosidase Td2F2 isolated from a compost metagenome has high glucose tolerance and transglycosylation activity. In this study, we determined the high-resolution crystal structure of Td2F2. It has a unique structure at the -1 subsite that is important for substrate specificity but not for glucose tolerance. To elucidate the mechanism(s) of glucose tolerance, we isolated a glucose-sensitive Td2F2 mutant using random mutagenesis. In this mutant, Asn223 residue located between subsites +1 and +2 was mutated. The Asn223 mutation resulted in reduced glucose tolerance and transglycosylation activity, and drastically changed substrate specificity. These results indicate that the structure between subsites +1 and +2 is critical for the glucose tolerance and substrate specificity of Td2F2. Our findings shed light on the glucose tolerance and transglycosylation activity mechanisms of glycoside hydrolase family 1 β-glucosidases. The atomic coordinates and structure factors (codes 3WH5, 3WH6, 3WH8, 3WH7, 5AYB, and 5AYI) have been deposited in the Protein Data Bank (http://wwpdb.org/).

Insights into electron leakage in the reaction cycle of cytochrome P450 BM3 revealed by kinetic modeling and mutagenesis.

As a single polypeptide, cytochrome P450 BM3 fuses oxidase and reductase domains and couples each domain's function to perform catalysis with exceptional activity upon binding of substrate for hydroxylation. Mutations introduced into the enzyme to change its substrate specificity often decrease coupling efficiency between the two domains, resulting in unproductive consumption of cofactors and formation of water and/or reactive species. This phenomenon can correlate with leakage, in which P450 BM3 uses electrons from NADPH to reduce oxygen to water and/or reactive species even without bound substrate. The physical basis for leakage is not yet well understood in this particular member of the cytochrome P450 family. To clarify the relationship between leakage and coupling, we used simulations to illustrate how different combinations of kinetic parameters related to substrate-free consumption of NADPH and substrate hydroxylation can lead to either minimal effects on coupling or a dramatic decrease in coupling as a result of leakage. We explored leakage in P450 BM3 by introducing leakage-enhancing mutations and combining these mutations to assess whether doing so increases leakage further. The variants in this study provide evidence that while a transition to high spin may be vital for coupled hydroxylation, it is not required for enhanced leakage; substrate binding and the consequent shift in spin state are not necessary as a redox switch for catalytic oxidation of NADPH. Additionally, the variants in this study suggest a tradeoff between leakage and stability and thus evolvability, as the mutations we investigated were far more deleterious than other mutations that have been used to change substrate specificity.

Function of MsiR on canavanine-mediated repression in Mesorhizobium tianshanense

Mesorhizobium tianshanense employs MsiA as canavanine exporter, which is upregulated by MsiR, to successfully form a symbiosis with the legume Glycyrrhiza uralensis. In this research, through gel-shift and bacterial two-hybrid examination, MsiR was found to spontaneously form dimers and bind to msiA promoter without additional canavanine. Six truncated forms of MsiR were constructed, and the conserved helix-turn-helix (HTH), substrate-binding, and surface-loop domains were found essential for MsiR functions. Random mutagenesis was used to study the functional sites of MsiR. Seven point mutants were selected, in which three mutants constitutively induced msiA expression without additional canavanine, two mutants partially changed substrate specificity, and the other two were almost null mutants. Results from the site mutation show that the functional subunits (HTH domain, dimerization interface domains, and C-terminal) are important in the conformation and induction ability of MsiR.

Genetic analysis of the Warburg effect in yeast.

We recently discovered that the Warburg effect, defined by the dramatically enhanced metabolism of glucose to pyruvate, even in well-oxygenated cancer cells, can occur as a consequence of mutations that enhance lipid biosynthesis at the expense of respiratory capacity. Specifically, mutations in the E1 subunit of either of two respiratory enzymes, pyruvate dehydrogenase (PDC) or α-ketoglutarate dehydrogenase (KGDC), change substrate specificity from the 3-carbon α-ketoacid pyruvate, or the 5-carbon α-ketoacid α-ketoglutarate, to the 4-carbon α-ketoacid oxaloacetate (OADC). These mutations result in OADC-catalyzed synthesis of malonyl-CoA (MaCoA), the essential precursor of all fatty acids. These mutants arose as spontaneous suppressors of a yeast acc1(cs) cold-sensitive mutation encoding an altered form of AcCoA carboxylase (Acc1) that fails to produce MaCoA at the restrictive temperature (16 °C). Notably, these suppressors are respiratory defective as a result of the same nuclear mutations that suppress acc1(cs). These mutants also suppress sensitivity to Soraphen A, a potent inhibitor of Acc1 activity, at normal temperature (30 °C). To our knowledge, OADC activity has never been identified in eukaryotic cells. Our results offer a novel perspective on the Warburg effect: the reprogramming of energy metabolism in cancer cells as a consequence of mutational impairment of respiration to meet the fatty acid requirements of rapidly proliferating cells. We suggest OADC activity is a common feature of cancer cells and represents a novel target for the development of chemotherapeutics.

The role of substrate strain in the mechanism of the carbon–carbon lyases

The carbon–carbon lyases, tryptophan indole lyase (TIL) and tyrosine phenol-lyase (TPL) are bacterial enzymes which catalyze the reversible elimination of indole and phenol from l-tryptophan and l-tyrosine, respectively. These PLP-dependent enzymes show high sequence homology (∼40% identity) and both form homotetrameric structures. Steady state kinetic studies with both enzymes show that an active site base is essential for activity, and α-deuterated substrates exhibit modest primary isotope effects on kcat and kcat/Km, suggesting that substrate deprotonation is partially rate-limiting. Pre-steady state kinetics with TPL and TIL show rapid formation of external aldimine intermediates, followed by deprotonation to give quinonoid intermediates absorbing at about 500nm. In the presence of phenol and indole analogues, 4-hydroxypyridine and benzimidazole, the quinonoid intermediates of TPL and TIL decay to aminoacrylate intermediates, with λmax at about 340nm. Surprisingly, there are significant kinetic isotope effects on both formation and subsequent decay of the quinonoid intermediates when α-deuterated substrates are used. The crystal structure of TPL with a bound competitive inhibitor, 4-hydroxyphenylpropionate, identified several essential catalytic residues: Tyr-71, Thr-124, Arg-381, and Phe-448. The active sites of TIL and TPL are highly conserved with the exceptions of these residues: Arg-381(TPL)/Ile-396 (TIL); Thr-124 (TPL)/Asp-137 (TIL), and Phe-448 (TPL)/His-463 (TIL). Mutagenesis of these residues results in dramatic decreases in catalytic activity without changing substrate specificity. The conserved tyrosine, Tyr-71 (TPL)/Tyr-74 (TIL) is essential for elimination activity with both enzymes, and likely plays a role as a proton donor to the leaving group. Mutation of Arg-381 and Thr-124 of TPL to alanine results in very low but measurable catalytic activity. Crystallography of Y71F and F448H TPL with 3-fluoro-l-tyrosine bound demonstrated that there are two quinonoid structures, relaxed and tense. In the relaxed structure, the substrate aromatic ring is in plane with the Cβ–Cγ bond, but in the tense structure, the substrate aromatic ring is about 20° out of plane with the Cβ–Cγ bond. In the tense structure, hydrogen bonds are formed between the substrate OH and the guanidinium of Arg-381 and the OH of Thr-124, and the phenyl rings of Phe-448 and 449 provide steric strain. Based on the effects of mutagenesis, the substrate strain is estimated to contribute about 108 to TPL catalysis. Thus, the mechanisms of TPL and TIL require both substrate strain and acid/base catalysis, and substrate strain is probably responsible for the very high substrate specificity of TPL and TIL.

Tailoring D‐Amino Acid Oxidase from the Pig Kidney to R‐Stereoselective Amine Oxidase and its Use in the Deracemization of α‐Methylbenzylamine

The deracemization of racemic amines to yield enantioenriched amines using S-stereoselective amine oxidases (AOx) has recently been attracting attention. However, R-stereoselective AOx that are suitable for deracemization have not yet been identified. An R-stereoselective AOx was now evolved from porcine kidney D-amino acid oxidase (pkDAO) and subsequently use for the deracemization of racemic amines. The engineered pkDAO, which was obtained by directed evolution, displayed a markedly changed substrate specificity towards R amines. The mutant enzyme exhibited a high preference towards the substrate α-methylbenzylamine and was used to synthesize the S amine through deracemization. The findings of this study indicate that further investigations on the structure-activity relationship of AOx are warranted and also provide a new method for biotransformations in organic synthesis.

Tailoring D‐Amino Acid Oxidase from the Pig Kidney to R‐Stereoselective Amine Oxidase and its Use in the Deracemization of α‐Methylbenzylamine

AbstractThe deracemization of racemic amines to yield enantioenriched amines using S‐stereoselective amine oxidases (AOx) has recently been attracting attention. However, R‐stereoselective AOx that are suitable for deracemization have not yet been identified. An R‐stereoselective AOx was now evolved from porcine kidney D‐amino acid oxidase (pkDAO) and subsequently use for the deracemization of racemic amines. The engineered pkDAO, which was obtained by directed evolution, displayed a markedly changed substrate specificity towards R amines. The mutant enzyme exhibited a high preference towards the substrate α‐methylbenzylamine and was used to synthesize the S amine through deracemization. The findings of this study indicate that further investigations on the structure–activity relationship of AOx are warranted and also provide a new method for biotransformations in organic synthesis.

Substrate-specific changes in mitochondrial respiration in skeletal and cardiac muscle of hibernating thirteen-lined ground squirrels

During torpor, the metabolic rate (MR) of thirteen-lined ground squirrels (Ictidomys tridecemlineatus) is considerably lower relative to euthermia, resulting in part from temperature-independent mitochondrial metabolic suppression in liver and skeletal muscle, which together account for ~40% of basal MR. Although heart accounts for very little (<0.5%) of basal MR, in the present study, we showed that respiration rates were decreased up to 60% during torpor in both subsarcolemmal (SS) and intermyofibrillar (IM) mitochondria from cardiac muscle. We further demonstrated pronounced seasonal (summer vs. winter [i.e., interbout] euthermia) changes in respiration rates in both mitochondrial subpopulations in this tissue, consistent with a shift in fuel use away from carbohydrates and proteins and towards fatty acids and ketones. By contrast, these seasonal changes in respiration rates were not observed in either SS or IM mitochondria isolated from hind limb skeletal muscle. Both populations of skeletal muscle mitochondria, however, did exhibit metabolic suppression during torpor, and this suppression was 2- to 3-fold greater in IM mitochondria, which provide ATP for Ca(2+)- and myosin ATPases, the activities of which are likely quite low in skeletal muscle during torpor because animals are immobile. Finally, these changes in mitochondrial respiration rates were still evident when standardized to citrate synthase activity rather than to total mitochondrial protein.

Evidence of Ternary Complex Formation in Trypanosoma cruzi trans-Sialidase Catalysis

Trypanosoma cruzi trans-sialidase (TcTS) is a key target protein for Chagas disease chemotherapy. In this study, we investigated the implications of active site flexibility on the biochemical mechanism of TcTS. Molecular dynamics studies revealed remarkable plasticity in the TcTS catalytic site, demonstrating, for the first time, how donor substrate engagement with the enzyme induces an acceptor binding site in the catalytic pocket that was not previously captured in crystal structures. Furthermore, NMR data showed cooperative binding between donor and acceptor substrates, supporting theoretical results. In summary, our data put forward a coherent dynamic framework to understand how a glycosidase evolved its highly efficient trans-glycosidase activity.

Multisite phosphorylation networks as signal processors for Cdk1

The order and timing of cell cycle events is controlled by changing substrate specificity and different activity thresholds of cyclin-dependent kinases (CDK). However, it is not understood how a single protein kinase can trigger hundreds of switches in a sufficiently time-resolved fashion. We show that the cyclin-Cdk1-Cks1-dependent phosphorylation of multisite targets in Saccharomyces cerevisiae is controlled by key substrate parameters including distances between phosphorylation sites, the distribution of serines and threonines as phospho-acceptors, and the positioning of cyclin-docking motifs. The component mediating the key interactions in this process is Cks1, the phospho-adaptor subunit of the cyclin-Cdk1-Cks1 complex. We propose that variation of these parameters within the networks of phosphorylation sites in different targets provides a wide range of possibilities for the differential amplification of Cdk1 signals, providing a mechanism to generate a wide range of thresholds in the cell cycle.

Rapid birth-and-death evolution of the xenobiotic metabolizing NAT gene family in vertebrates with evidence of adaptive selection

BackgroundThe arylamine N-acetyltransferases (NATs) are a unique family of enzymes widely distributed in nature that play a crucial role in the detoxification of aromatic amine xenobiotics. Considering the temporal changes in the levels and toxicity of environmentally available chemicals, the metabolic function of NATs is likely to be under adaptive evolution to broaden or change substrate specificity over time, making NATs a promising subject for evolutionary analyses. In this study, we trace the molecular evolutionary history of the NAT gene family during the last ~450 million years of vertebrate evolution and define the likely role of gene duplication, gene conversion and positive selection in the evolutionary dynamics of this family.ResultsA phylogenetic analysis of 77 NAT sequences from 38 vertebrate species retrieved from public genomic databases shows that NATs are phylogenetically unstable genes, characterized by frequent gene duplications and losses even among closely related species, and that concerted evolution only played a minor role in the patterns of sequence divergence. Local signals of positive selection are detected in several lineages, probably reflecting response to changes in xenobiotic exposure. We then put a special emphasis on the study of the last ~85 million years of primate NAT evolution by determining the NAT homologous sequences in 13 additional primate species. Our phylogenetic analysis supports the view that the three human NAT genes emerged from a first duplication event in the common ancestor of Simiiformes, yielding NAT1 and an ancestral NAT gene which in turn, duplicated in the common ancestor of Catarrhini, giving rise to NAT2 and the NATP pseudogene. Our analysis suggests a main role of purifying selection in NAT1 protein evolution, whereas NAT2 was predicted to mostly evolve under positive selection to change its amino acid sequence over time. These findings are consistent with a differential role of the two human isoenzymes and support the involvement of NAT1 in endogenous metabolic pathways.ConclusionsThis study provides unequivocal evidence that the NAT gene family has evolved under a dynamic process of birth-and-death evolution in vertebrates, consistent with previous observations made in fungi.