| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Medicine (K.G., P.E., S.B., A.H., R.M.F.), Atherosclerosis Research Unit, King Gustaf V Research Institute, and Department of Medical Biochemistry and Biophysics (M.S.), Karolinska Institutet, S-171 76 Stockholm, Sweden; and Department of Organic Chemistry (M.B.), University of Padova, 35131 Padova, Italy
Address all correspondence and requests for reprints to: Rachel Fisher, King Gustaf V Research Institute, Karolinska Hospital, S-171 76 Stockholm, Sweden. E-mail: rachel.fisher{at}medks.ki.se.
| Abstract |
|---|
|
|
|---|
Heterozygotes for a Gly209Ser polymorphism (Ser allele frequency 0.05) had significantly lower body mass index and, correcting for body mass index, significantly lower triglyceride concentrations, systolic blood pressure, insulin concentrations, and homeostasis model assessment index compared with common homozygotes. A three-dimensional model of the FATP4 protein based on structural and functional similarity with adenylate-forming enzymes revealed that the variable residue 209 is exposed in a region potentially involved in protein-protein interactions. Furthermore, the model indicated functional regions with respect to NEFA transport and acyl-coenzyme A synthase activity and membrane association.
These findings propose FATP4 as a candidate gene for the insulin resistance syndrome and provide a structural basis for understanding FATP function in NEFA transport and metabolism.
| Introduction |
|---|
|
|
|---|
FATP4 belongs to a highly homologous and evolutionary conserved gene family with six members (FATP16) in man, each with a distinct expression pattern (7). FATP4 expression is most prominent in central nervous tissue, intermediate in intestine, heart, liver, and pancreas, is relatively low in skeletal muscle, and has been observed in adipose tissue (8). The human FATP4 gene consists of 12 coding exons spanning more than 17 kb of genomic DNA on chromosome 9q34, giving rise to a 71-kDa protein, containing an AMP-binding motif (also called phosphate-binding loop or p-loop) conserved in FATP16, acyl-coenzyme A (acyl-CoA) synthases and other members of the adenylate-forming enzyme superfamily (9, 10). The AMP-binding motif has been suggested to convey acyl-CoA synthase activity to FATPs and homologous proteins (9, 11, 12) as part of their proposed role in NEFA metabolism. The importance of FATP4 function in long-chain fatty acid uptake by enterocytes, in which FATP4 is the most abundant FATP, has been demonstrated in vitro (13).
We have previously shown that the rare allele of an intronic polymorphism in the FATP1 gene was associated with increased postprandial triglyceridemia and a shift toward a smaller, more dense LDL phenotype (14). The aims of this study were to investigate variation within the FATP4 gene with respect to associations with fasting and postprandial lipid and lipoprotein variables and markers of insulin resistance in healthy, middle-aged Swedish men and to evaluate possible mechanisms behind any associations observed.
| Subjects and Methods |
|---|
|
|
|---|
A total of 608 50-yr-old men living in the county of Stockholm were selected at random from a registry of permanent residents (15). Exclusion criteria were non-Caucasian descent, chronic disease, history of cardiovascular disease, familial hypercholesterolemia, alcohol abuse, psychiatric disorders, and participation in other ongoing studies. Samples were collected after an overnight fast. Studies were approved by the Karolinska Hospital local ethics committee, and all subjects gave their informed consent.
Oral fat tolerance test
A subset of 105 consecutive individuals from the cohort with the apolipoprotein (apo) E3/E3 genotype underwent an oral fat tolerance test (OFTT) of the mixed-meal type (total energy content, 1000 kcal; 60% energy from fat) (15). Blood samples were drawn before and up to 6 h after food intake.
Biochemical analyses
Triglyceride (TG) and cholesterol concentrations in plasma and major lipoproteins were determined enzymatically (15). Very low-density lipoprotein (VLDL), LDL, and HDL were isolated from fasting plasma by a combination of preparative ultracentrifugation and precipitation of apoB-containing lipoproteins (15). Fasting and postprandial TG-rich lipoproteins; Svedberg floatation rate [Sf > 400; Sf 60400; and Sf 2060] were subfractionated by cumulative density ultracentrifugation (15). ApoB48 and apoB100 concentrations in isolated lipoprotein subfractions were determined by SDS-PAGE, staining with Coomassie G-250, and subsequent scanning (15). Blood glucose, plasma NEFA, insulin, and proinsulin concentrations were measured with standard methods (15).
Sequencing
PCR was used to amplify regions surrounding reported singlenucleotide polymorphisms (SNPs) in the FATP4 gene, namely a C/T polymorphism at position -856 of the putative proximal promoter region, a G/A polymorphism at position 625 in exon 3 (Gly209Ser), a C/T polymorphism at position 1542 in exon 10, and a G/A polymorphism at position 149 of intron 11. For the proximal promoter region, the primers 5'-CCTGGATGTCAGAAGAGTGAACTG-3' and 5'-TGCTGGTAGAGAACATGAGGTCTG-3' were used in 1.5 mM MgCl2 with 2 U AmpliTaq Gold (Applied Biosystems, Piscataway, NJ). Hot-start touch-down thermocycling was performed with 25 final cycles of 94 C for 45 sec, 53 C for 30 sec, and 72 C for 90 sec. For the exon 3 fragment, the primers 5'-GTGAGGTCCATGCCAGCCTG-3' and 5'-CACCTGTGAAGCCCTTGTCAG-3' were used in 2.0 mM MgCl2 with 1 U of AmpliTaq Gold. Hot-start touch-down thermocycling was performed with 25 final cycles of 30-sec steps at 94 C, 56 C, and 72 C. For the exon 10 C/T and intron 11 G/A polymorphisms, the primers 5'-GACGCAGGGTACACCTGGTGACAG-3' and 5'-CCTTGTCCTAGGCCAGCCTCTCTC-3' were used in 1.5 mM MgCl2 with 2 U Taq DNA polymerase (Promega, Madison, WI) using touch-down thermocycling with 30 final cycles of 94 C for 30 sec, 56 C for 30 sec, and 72 C for 60 sec. Purified PCR products were sequenced with the primers 5'-CCAGACTTCGCTCGTTCTCGCATG-3', 5'-TTGTCAGGGCAACTGGGAAG-3', 5'-CCTTGTGGTCAAACAAATCCTTGG-3', and 5'-GCTCAGGTCTTGGAGAAGGAACTG-3' for the proximal promoter, exon 3, exon 10, and intron 11 polymorphisms, respectively, using the DYEnamic ET terminator cycle sequencing premix kit (Amersham Pharmacia Biotech, Piscataway, NJ). Sequencing reactions were run on an ABI 377 automated DNA sequencer and the results analyzed with the ABI Prism software (Amersham Pharmacia Biotech).
Genetic analysis
Fragments surrounding the reported FATP4 exon 3 G/A (Gly209Ser) and intron 11 G/A polymorphisms and a fragment surrounding a novel intron 11 4T/3T polymorphism that was detected by sequencing were PCR amplified for subsequent analysis by pyrosequencing (16). For the exon 3 polymorphism, fragments were amplified as above using 6 pmol of each primer 5'-biotin-GTGAGGTCCATGCCAGCCTG-3' and 5'-CACCTGTGAAGCCCTTGTCAG-3'. Fragments surrounding the intron 11 G/A and 4T/3T polymorphisms were amplified with 5 pmol of each primer 5'-biotin-GGCTGTTACTTGACCTCTGCACAC-3' and 5'-GCGCTGTGCCTGGCACCTAATATG-3' in 1.5 mM MgCl2 and 5 pmol of each primer 5'-biotin-GGCTGTTACTTGACCTCTGCACAC-3' and 5'-GCCTCTCTCCTGCTCTAAAACTCC-3' in 1.75 mM MgCl2, respectively, with 1 U Taq DNA polymerase (Promega) using touch-down thermocycling with 20 final cycles of 30-sec steps at 94 C, 56 C, and 72 C. Samples were prepared and analyzed on a PSQ96 instrument using the SNP reagent kit in accordance with the manufacturers recommendations (Pyrosequencing AB, Uppsala, Sweden), with 15 pmol of pyrosequencing primer (5'-GTGCTTGGAGGCACCGCAC-3', 5'-CAATGTGTGTCTTCTCTGCCAGCC-3' and 5'-CTGCTCTAAAACTCCAAGAAGGCAAA-3' for the exon 3 G/A, intron 11 G/A, and intron 11 4T/3T polymorphisms, respectively). ApoE genotyping was performed using PCR and restriction enzyme digestion (17).
Protein structure modeling
Modeling template protein structures were identified by the threading algorithm 3D-PSSM (18) using the FATP4 cDNA (GenBank accession no. NP_005085) as search query. The final model was derived using Modeller4 (19), in which the 3D-PSSM-derived alignment was adjusted to relieve local energetic stress and to optimize the alignment in regions of low similarity. Results were visualized with the RasMol (20) or MolMol (21) software. A consensus of several methods was used to predict FATP4 membrane topology (PSIPRED, PRED-TMR, HMMTOP, TMHMM, TMBASE, www.hgmp.mrc.ac.uk/GenomeWeb/prot-transmembrane.html) and secondary structure (22), respectively. Sequence conservation within the FATP family was visualized in the structural model by ConSurf (23). Sequence alignments were also performed with the BLAST and ClustalW algorithms.
Statistical analysis
The Statview (SAS Institute Inc., Cary, NC) software was used to perform ANOVA and repeated-measures ANOVA. Fishers protected least significant difference tests were used for post hoc analysis. Skewed data were logarithmically transformed before analysis. Homeostasis model assessment (HOMA) index was calculated as the product of fasting insulin and glucose concentrations divided by 22.5 (24). Linkage disequilibrium calculations were performed with the Associate software 2.35 (J. Ott, Rockefeller University, New York, NY). Statistical significance was assigned to a value of P < 0.05.
| Results |
|---|
|
|
|---|
In a cohort of 608 healthy 50-yr-old Swedish men, frequencies of the rare alleles were 0.05 for both the FATP4 exon 3 G/A (Gly209Ser) and intron 11 G/A polymorphisms reported in the National Center for Biotechnology Information SNP database (Table 1
). Notably, Ser209 is reported as wild-type in murine FATP4 (7, 25). A T/C polymorphism at position -856 in the FATP4 putative proximal promoter region and a C/T polymorphism at position 1542 in exon 10 reported in the NCBI SNP database were not found after direct sequencing of 48 alleles. A novel 4T/3T polymorphism at position 249 in intron 11, identified when screening for the intron 11 G/A polymorphism, exhibited a frequency of 0.10 for the rare 3T allele (Table 1
). No obvious impact on any consensus transcription factor binding site was judged to be conveyed by either the intron 11 G/A polymorphism or the intron 11 4T/3T polymorphism (Ref.26 ; http://www.cbrc.jp/research/db/TFSEARCH.html). A novel G/A polymorphism at position -908 identified when screening for the -856 T/C polymorphism showed a frequency of only 0.01 for the rare allele and was not observed to significantly alter any consensus transcription factor binding site (Ref.26 ; http://www.cbrc.jp/research/db/TFSEARCH.html) and was therefore not investigated further. The Gly209Ser and intron 11 4T/3T polymorphic loci were in strong linkage disequilibrium, whereas the intron 11 G/A locus was not in strong linkage disequilibrium with either of the other two loci (Table 1
). All genotyping assays gave identical results, compared with those from direct sequencing of a limited number of individuals. Sequencing FATP4 exon 3, codon 194 was found to read CTG (Leu) in all of a total of 88 alleles, in accordance with the GenBank genomic clone AC006312, and not CCG (Pro) as in the GenBank FATP4 mRNA AF055899.
|
For each of the three polymorphisms, the number of rare homozygotes was low, especially in the OFTT subset, and certain parameters exhibited considerable variability within the rare homozygote groups. Therefore, comparisons of common homozygotes and heterozygotes were focused on (using post hoc tests), and little weight was attached to observations regarding rare homozygotes [for example, the high mean insulin concentration in Ser209 homozygotes, which is hard to interpret due to the very small sample size (n = 3) and large variability of this parameter in these individuals; Table 2
].
|
The phenotypic modulation associated with the intron 11 4T/3T polymorphism was similar to that of the Gly209Ser polymorphism regarding some variables (BMI, SBP, and TG concentrations; Table 2
), as expected, due to the almost complete linkage disequilibrium between these two polymorphic loci. However, when analyzing the intron 11 4T/3T polymorphism in Gly209 homozygotes only, BMI was the only variable that remained significantly different between genotypes (4T/4T, 26.2 ± 0.15 kg/m2; n = 475; 4T/3T, 25.0 ± 0.39; n = 62; P < 0.01). Heterozygotes for the intron 11 G/A polymorphism had significantly lower LDL cholesterol concentrations than common homozygotes (Table 2
).
In a subset of 107 subjects with the apoE3/E3 genotype who underwent an OFTT allele frequencies were similar to those in the cohort as a whole (Table 1
). Plasma TG concentrations at each time point (measured every hour 06 h after food intake) tended to be lower in carriers of the rare allele of the Gly209Ser polymorphism (Fig. 1
) as well as the intronic polymorphisms (data not shown). Similarly, for all three polymorphisms, heterozygotes exhibited a tendency toward lower NEFA concentrations, chylomicron/chylomicron remnant (Sf > 400) apoB48, and VLDL-I (Sf 60400) and VLDL-II (Sf 2060) apoB100 concentrations measured at 0, 3, and 6 h after food intake, compared with common homozygotes (data not shown). No significant differences between heterozygotes and common homozygotes in the postprandial lipid or lipoprotein response pattern, as judged by repeated-measures ANOVA, were observed for any of the three polymorphisms.
|
To investigate structural-functional aspects of FATP4, including potential effects of the Gly209Ser substitution on FATP4 function, a structural model of the FATP4 protein structure was generated. The FATP4 protein was judged to adopt the fold of the adenylate-forming enzyme superfamily because several members of this family were picked up by both BLAST searches of the Protein Data Bank (PDB) and by the threading algorithm 3D-PSSM (18) using the FATP4 sequence as search query. The resultant candidate structural templates identified for FATP4 were the PDB entries 1PG4 (Salmonella enterica acetyl-CoA synthase) (seACS) (27), 1MD9 (Bacillus subtilis 2,3-dihydroxybenzoate-AMP ligase) (10), and 1BA3 [Photinus pyralis luciferase (28), represented by 1LCI (29) in 3D-PSSM]. Despite low overall sequence identities (1520%), the AMP-binding motif was nearly identical, as were a number of strictly invariant residues within the superfamily of adenylate-forming enzymes (29) (Fig. 2
). seACS was chosen as template for modeling of the FATP4 structure because it scored highest, both with respect to sequence similarity and structural conservation (>95% certainty in 3D-PSSM).
|
|
The variable residue 209 was located in an exposed loop connecting two ß-strands, in a region of the protein predicted to face the cytosol. A Gly-to-Ser mutation could well influence the loop structure due to structural limitations of adoptable backbone angles. Although flanked by sequence stretches similar with seACS, the exposed residues surrounding FATP4 residue 209 exhibited relatively high variability within the FATP family as well as within the structurally conserved adenylate-forming enzyme family (Figs. 2
and 3
) (10), suggesting that these residues may convey protein-specific properties. Furthermore, the variable residue 209 was flanked by proline residues (SWEPG/SAVPP), which is typical of short interactive protein segments (31), and other hydrophobic residues. Taken together, this could well indicate that these residues are involved in protein-specific interactions.
| Discussion |
|---|
|
|
|---|
An intron 11 4T/3T polymorphism exhibited phenotypic associations resembling those of the Gly209Ser polymorphism, and an intron 11 G/A polymorphism was associated with LDL cholesterol concentrations. When analyzing the intron 11 4T/3T polymorphism in Gly209 homozygotes only, most of the significant associations were lost, indicating that the phenotype associated with the intron 11 4T/3T polymorphism was to a large extent dependent on the strong linkage disequilibrium between the intron 11 3T and Ser209 alleles. This is supported by the potential functionality of the Gly209Ser substitution, as indicated by the structural evaluation of this amino acid change and by the fact that there was no clear implication of either of the intronic polymorphisms in differential transcriptional regulation. However, influence of the intronic polymorphisms on co- or posttranscriptional processes cannot be ruled out.
The striking phenotypic differences associated with the Gly209Ser polymorphism may readily be explained by variations in NEFA uptake and metabolism and more specifically by variations in FATP4 activity. NEFA concentrations are strongly related to insulin sensitivity, in that high NEFA concentrations are associated with reduced peripheral insulin-stimulated glucose uptake and increased hepatic glucose output (1). This might be explained by high NEFA concentrations metabolically influencing glucose homeostasis (1, 34), NEFA-mediated transcriptional regulation (35), and physical-chemical effects of NEFAs. The observed variations in both fasting and postprandial TG concentrations could be explained by differences in peripheral TG clearance and/or hepatic VLDL-TG secretion. The relatively more insulinresistant phenotypes associated with the common Gly209 allele may include impaired insulin-mediated inhibition of VLDL secretion and adipose tissue hormone-sensitive lipase activity, and impaired insulin-stimulated activation of peripheral lipoprotein lipase. Furthermore, high local NEFA concentrations per se may impair lipoprotein lipase activity, and high circulating NEFA concentrations may increase hepatic VLDL-TG synthesis and secretion. An insulin resistance-related hypertensive phenotype may be caused by high concentrations of TGs or NEFAs inducing endothelial dysfunction by interfering with endothelial nitric oxide synthase activity (36, 37). FATP4 expression levels in isolated enterocytes are strongly correlated with uptake of long-chain FAs (13), suggesting that the observations in the postprandial test could be explained by differences in gastrointestinal absorption of FAs. However, because the majority of the observations in this study concern fasting variables, we assume that the primary effects associated with the polymorphisms studied also involve tissues other than gastrointestinal, such as adipose tissue, skeletal muscle, or liver.
In the proposed FATP4 structural model, the variable residue 209 was located in an exposed hydrophobic loop structure flanked by proline residues, suggested to convey protein-specific interactive properties. The candidate structural template B. subtilis 2,3-dihydroxybenzoate-AMP ligase and several other adenylate-forming enzymes are part of multienzyme complexes (10), indicating that the FATP4 structure also may comprise regions for potential protein-protein interaction. Such interactions could involve cytosolic tissue-specific FA binding proteins (FABPs) (38), acyl-CoA synthases, or CD36, the latter two having been shown to colocalize with FATPs at the plasma membrane (39, 40). Interestingly, residues 191475 of murine FATP1 have been shown to direct formation of functional FATP1 homodimers (41). Thus, the phenotypic modulation associated with the Gly209Ser substitution might be due to allelic differences in protein-protein interactions required either for the formation of functional FATP4 dimers or for interaction of FATP4 with other proteins involved in cellular FA handling.
The current study suggests a structural and functional classification of FATP4, and thereby other FATPs, as membrane-associated members of the adenylate-forming enzyme superfamily (including acyl-CoA synthases). Despite low overall sequence similarities, functionally critical regions in this superfamily were very similar in FATP4, both at the levels of sequence and structure (Figs. 2
and 3
). A substantial amount of experimental data regarding intrinsic FATP acyl-CoA synthase activity, and the importance of the AMP-binding motif for FATP function (9, 11, 12, 25), is in excellent agreement with this classification.
The functionality of the proposed FATP4 catalytic cleft is supported by the fact that conserved residues, implicated in AMP binding in FATP4 (Figs. 2
and 3
), were shown important for azido-ATP binding, acyl-CoA synthase activity, and FA transport activity in FATP1 and the yeast FATP homolog Fat1p by directed mutagenesis (9, 11, 12), as described above. Furthermore, other mutations that severely affected Fat1p function (12) (corresponding to FATP4 Tyr499, Phe508, and Ser516) map to the cleft region (Fig. 3
), as does a conserved fatty acyl-CoA synthase (FACS) signature motif proposed to contribute to the FA binding site in the Escherichia coli FACS (42, 43) (Fig. 3
). Interestingly, certain of the above Fat1p mutations selectively affected either acyl-CoA synthase activity or FA import (12). The proximity of the proposed catalytic region and membrane-binding surface suggests that they might be functionally linked in the context of FA binding. Of note, the interdomain orientation of adenylate-forming enzymes has been proposed to be significantly altered during their two-step catalytic reaction (adenylation and thioesterification by acyl-CoA synthases) (27). This suggests that the interface between the two domains could be important for such a reorientation as a functional feature in FATP-mediated acyl-CoA formation and FA transport.
In summary, we conclude that alterations in NEFA uptake and/or metabolism may underlie the phenotypes observed in this study and that these alterations may be caused by Gly209Ser genotype-associated variations in FATP4 activity in metabolically relevant tissues. Furthermore, we propose a model of FATP4 structure, including identification of regions with specific functionality, that will prove useful in further studies to evaluate the physiological and pathological role(s) of the FATPs.
| Footnotes |
|---|
Present address for M.S.: Molecular Biotechnology, IFM, Linköping University, S-581 83 Linköping, Sweden.
Abbreviations: acyl-CoA, Acyl-coenzyme A; apo, apolipoprotein; BMI, body mass index; FA, fatty acid; FACS, fatty acyl-CoA synthase; FATP, fatty acid transport protein; HDL, high-density lipoprotein; HOMA, homeostasis model assessment; LDL, low-density lipoprotein; NEFA, nonesterified fatty acid; OFTT, oral fat tolerance test; SBP, systolic blood pressure; seACS, Salmonella enterica acetyl-CoA synthase; Sf, Svedberg floatation rate; SNP, single-nucleotide polymorphism; TG, triglyceride; VLDL, very low-density lipoprotein.
Received April 22, 2003.
Accepted September 29, 2003.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. Cunningham and L. McDermott Long Chain PUFA Transport in Human Term Placenta J. Nutr., April 1, 2009; 139(4): 636 - 639. [Full Text] [PDF] |
||||
![]() |
P. A. Watkins Very-long-chain Acyl-CoA Synthetases J. Biol. Chem., January 25, 2008; 283(4): 1773 - 1777. [Full Text] [PDF] |
||||
![]() |
S. Lobo, B. M. Wiczer, A. J. Smith, A. M. Hall, and D. A. Bernlohr Fatty acid metabolism in adipocytes: functional analysis of fatty acid transport proteins 1 and 4 J. Lipid Res., March 1, 2007; 48(3): 609 - 620. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Doege and A. Stahl Protein-mediated Fatty Acid uptake: novel insights from in vivo models. Physiology, August 1, 2006; 21: 259 - 268. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. M. Hall, B. M. Wiczer, T. Herrmann, W. Stremmel, and D. A. Bernlohr Enzymatic Properties of Purified Murine Fatty Acid Transport Protein 4 and Analysis of Acyl-CoA Synthetase Activities in Tissues from FATP4 Null Mice J. Biol. Chem., March 25, 2005; 280(12): 11948 - 11954. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Doron-Faigenboim, A. Stern, I. Mayrose, E. Bacharach, and T. Pupko Selecton: a server for detecting evolutionary forces at a single amino-acid site Bioinformatics, January 1, 2005; 21(9): 2101 - 2103. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |