Structural motif-based homology modeling of CYP27A1 and site-directed mutational analyses affecting vitamin D hydroxylation

Human CYP27A1 is a mitochondrial cytochrome P450, which is principally found in the liver and plays important roles in the biological activation of vitamin D(3) and in the biosynthesis of bile acids. We have applied a systematic analysis of hydrogen bonding patterns in 11 prokaryotic and mammalian CYP crystal structures to construct a homology-based model of CYP27A1. Docking of vitamin D(3) structures into the active site of this model identified potential substrate contact residues in the F-helix, the beta-3 sheet, and the beta-5 sheet. Site-directed mutagenesis and expression in COS-1 cells confirmed that these positions affect enzymatic activity, in some cases shifting metabolism of 1alpha-hydroxyvitamin D(3) to favor 25- or 27-hydroxylation. The results suggest that conserved hydrophobic residues in the beta-5 hairpin help define the shape of the substrate binding cavity and that this structure interacts with Phe-248 in the F-helix. Mutations directed toward the beta-3a strand suggested a possible heme-binding interaction centered on Asn-403 and a structural role for substrate contact residues Thr-402 and Ser-404.     

Generation of a homology model for the human cytochrome P450, CYP24A1, and the testing of putative substrate binding residues by site-directed mutagenesis and enzyme activity studies

A systematic analysis of conserved H-bonding patterns and tertiary structural motifs from 13 crystal structures was used to create a homology model for the human multicatalytic cytochrome P450, CYP24A1, involved in catabolism of 1alpha,25-dihydroxyvitamin D3. The substrate was docked in the active site and used to identify potential substrate contact residues in the B' helix, B'/C loop, F-helix and the beta-5 hairpin. Seven CYP24A1 mutants were created and studied by mammalian cell transfection and CYP24A1 activity assay. Mutants showed reduced metabolic rates and altered metabolite patterns compared to wild-type. We conclude that: Ile-131 positions substrate via A-ring and cis-triene contacts; Trp-134 and Gly-499 are determinants of substrate access; Leu-148 contacts the substrate side-chain; Met-246 is important in mediating regioselectivity. Our findings validate the new model of CYP24A1, which can now be used to predict structural modifications for rational vitamin D drug design.     

The dynamics of camphor in the cytochrome P450 CYP101D2

The recent crystal structures of CYP101D2, a cytochrome P450 protein from the oligotrophic bacterium Novosphingobium aromaticivorans DSM12444 revealed that both the native (substrate‐free) and camphor‐soaked forms have open conformations. Furthermore, two other potential camphor‐binding sites were also identified from electron densities in the camphor‐soaked structure, one being located in the access channel and the other in a cavity on the surface near the F‐helix side of the F‐G loop termed the substrate recognition site. These latter sites may be key intermediate positions on the pathway for substrate access to or product egress from the active site. Here, we show via the use of unbiased atomistic molecular dynamics simulations that despite the open conformation of the native and camphor‐bound crystal structures, the underlying dynamics of CYP101D2 appear to be very similar to other CYP proteins. Simulations of the native structure demonstrated that the protein is capable of sampling many different conformational substates. At the same time, simulations with the camphor positioned at various locations within the access channel or recognition site show that movement towards the active site or towards bulk solvent can readily occur on a short timescale, thus confirming many previously reported in silico studies using steered molecular dynamics. The simulations also demonstrate how the fluctuations of an aromatic gate appear to control access to the active site. Finally, comparison of camphor‐bound simulations with the native simulations suggests that the fluctuations can be of similar level and thus are more representative of the conformational selection model rather than induced fit.     

Structure of the human lung cytochrome P450 2A13

The human lung cytochrome P450 2A13 (CYP2A13) activates the nicotine-derived procarcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) into DNA-altering compounds that cause lung cancer. Another cytochrome P450, CYP2A6, is also present in human lung, but at much lower levels. Although these two enzymes are 93.5% identical, CYP2A13 metabolizes NNK with much lower K(m) values than does CYP2A6. To investigate the structural differences between these two enzymes the structure of CYP2A13 was determined to 2.35A by x-ray crystallography and compared with structures of CYP2A6. As expected, the overall CYP2A13 and CYP2A6 structures are very similar with an average root mean square deviation of 0.5A for the Calpha atoms. Like CYP2A6, the CYP2A13 active site cavity is small and highly hydrophobic with a cluster of Phe residues composing the active site roof. Active site residue Asn(297) is positioned to hydrogen bond with an adventitious ligand, identified as indole. Amino acid differences between CYP2A6 and CYP2A13 at positions 117, 300, 301, and 208 relate to different orientations of the ligand plane in the two protein structures and may underlie the significant variations observed in binding and catalysis of many CYP2A ligands. In addition, docking studies suggest that residues 365 and 366 may also contribute to differences in NNK metabolism.     

Covalent linkage of prosthetic heme to CYP4 family P450 enzymes.

An extensive body of research on the structural properties of cytochrome P450 enzymes has established that these proteins possess a b-type heme prosthetic group which is noncovalently bound at the active site. Coordinate, electrostatic, and hydrogen bond interactions between the protein backbone and heme functional groups are readily overcome upon mild acid treatment of the enzyme, which releases free heme from the protein. In the present study, we have used a combination of HPLC, LC/ESI-MS, and SDS-PAGE techniques to demonstrate that members of the mammalian CYP4B, CYP4F, and CYP4A subfamilies bind their heme in an unusually tight manner. HPLC chromatography of CYP4B1 on a POROS R2 column under mild acidic conditions caused dissociation of less than one-third of the heme from the protein. Moreover, heme was not substantially removed from CYP4B1 under electrospray or electrophoresis conditions that readily release the prosthetic group from other non-CYP4 P450 isoforms. This was evidenced by an intact protein mass value of 59,217 +/- 3 amu for CYP4B1 (i.e., apoprotein plus heme) and extensive staining of this approximately 60 kDa protein with tetramethylbenzidine/H(2)O(2) following SDS-PAGE. In addition, treatment of CYP4B1, CYP4F3, and CYP4A5/7 with strong base generated a new, chromatographically distinct, polar heme species with a mass of 632.3 amu rather than 616.2 amu. This mass shift is indicative of the incorporation of an oxygen atom into the heme nucleus and is consistent with the presence of a novel covalent ester linkage between the protein backbone of the CYP4 family of mammalian P450s and their heme catalytic center.     

Gaussian neighborhood: A new measure of accessibility for residues of protein molecules

We introduce a new method for assessing the extent of residue exposure in proteins. For each atom of every residue a Gaussian-weighted atomic surroundings value (the G-neighborhood) is calculated. A normalized sum of G-neighborhood values over all the atoms of a residue is complementary to conventional surface accessibility characteristics. The G-0neighborhood value of a residue is a sensitive indicator of its location, strongly dependent on the 3D structure of a the protein. Correlations between secondary structures and patterns of G-neighborhood values for six different protein molecules are discussed. Comparison of the distribution of hydrophobic and charged residues in the 3D structure for the alcohol-soluble protein crambin and that of five water-soluble proteins (cytochrome c, flavodoxin, myoglobin, rhodanese, and Bence–Jones protein) shows striking differences in their G-neighborhood patterns. Contacts between the prosthetic group and the peptide portion of a protein as well as protein interdomain contacts and monomer–monomer contacts are characterized. © 1993 Wiley-Liss, Inc.     

Homology model of 1alpha,25-dihydroxyvitamin D3 24-hydroxylase cytochrome P450 24A1 (CYP24A1): active site architecture and ligand binding

Homology models of cytochrome P450 24A1 (CYP24A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using molecular operating environment (MOE) software the lowest energy CYP24A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP24A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP24A1 model. The natural substrate, 1,25-dihydroxyvitamin D(3) (calcitriol) and the CYP24 inhibitor (R)-N-(2-(1H-imidazol-1-yl)-2-phenylethyl)-4'-chlorobiphenyl-4-carboxamide ((R)-VID-400) were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure.     

Comparative modeling of 25-hydroxycholesterol-7α-hydroxylase (CYP7B1): ligand binding and analysis of hereditary spastic paraplegia type 5 CYP7B1 mutations

CYP7B1 mutations have been linked directly with the neurodegenerative disease hereditary spastic paraplegia (HSP), with mutations in the CYP7B1 gene identified as being directly responsible for autosomal recessive HSP type 5A (SPG5). To evaluate the potential impact of CYP7B1 mutations identified in SPG5 on binding and protein function, a comparative model of cytochrome P450 7B1 (CYP7B1) was constructed using human CYP7A1 as a template during model construction. The secondary structure was predicted using the PSIPRED and GOR4 prediction methods, the lowest energy CYP7B1 model was generated using MOE, and then this model was assessed in terms of stereochemical quality and the side chain environment using RAMPAGE, Verify3D and ProSA. Evaluation of the active site residues of the CYP7B1 model and validation of the active site architecture were performed via molecular docking experiments: the docking of the substrates 25-hydroxycholesterol and 27-hydroxycholesterol and the inhibitor 3α-Adiol identified structurally and functionally important residues. Mutational analysis of CYP7B1 amino acid mutations related to hereditary spastic paraplegia type 5 considered phosphorylation, ligand/substrate binding and the structural roles of mutated amino acid residues, with R112, T297 and S363 mutations expected to have a direct impact on ligand binding, while mutations involving R417 would indirectly affect ligand binding as a result of impairment in catalytic function.     

Unusual regioselectivity and active site topology of human cytochrome P450 2J2

The oxidation of six derivatives of terfenadone by recombinant human CYP2J2 (CYP = cytochrome P450) was studied by high-performance liquid chromatography coupled to mass spectrometry (MS) using tandem MS techniques and by 1H NMR spectroscopy. CYP2J2 exhibited a surprising regioselectivity in favor of the hydroxylation of the substrate terminal chain at the weakly reactive homobenzylic position. In contrast, hydroxylation of the same substrates by CYP3A4 mainly occurred on the most chemically reactive sites of the substrates (N-oxidation and benzylic hydroxylation). A 3D homology model of CYP2J2 was constructed using recently published structures of CYP2A6, CYP2B4, CYP2C8, CYP2C9, and CYP2D6 as templates. In contrast with other CYP2 structures, it revealed an active site cavity with a severely restricted access of substrates to the heme through a narrow hydrophobic channel. Dynamic docking of terfenadone derivatives in the CYP2J2 active site allowed one to interpret the unexpected regioselectivity of the hydroxylation of these substrates by CYP2J2, which is mainly based on this restricted access to the iron. The structural features that have been found to be important for recognition of substrates or inhibitors by CYP2J2 were also interpreted on the basis of CYP2J2-substrate interactions in this model.     

Decreasing the CYP2D6 contribution to metabolism of a CK1ε inhibitor

Our internal casein kinase 1ε lead inhibitor, compound 1 was partially cleared by the polymorphic cytochrome P450 2D6. CYP2D6 involvement in metabolism implies more extensive clinical trials. We therefore wanted to reduce the contribution to clearance by this enzyme. We utilized metabolism reports for compound 1 performed in recombinant CYP2D6 together with structure-metabolism variation in structures of closely related analogs in order to see if we could incorporate similar substitution patterns in our lead compound. In addition, we utilized a previously established docking method using a modified CYP2D6 crystal structure to see if the metabolism patterns in CYP2D6 could be reproduced to afford the metabolites in the metabolism reports as well as those for the compounds used in the structure-metabolism relationship. All three of these steps, the metabolism report, the establishment of structure-metabolism relationships and the docking, lead to compound 10 where CYP2D6 was not involved in the clearance pathways.     

Building an automated classification of DNA-binding protein domains

Intensive growth in 3D structure data on DNA-protein complexes as reflected in the Protein Data Bank (PDB) demands new approaches to the annotation and characterization of these data and will lead to a new understanding of critical biological processes involving these data. These data and those from other protein structure classifications will become increasingly important for the modeling of complete proteomes. We propose a fully automated classification of DNA-binding protein domains based on existing 3D-structures from the PDB. The classification, by domain, relies on the Protein Domain Parser (PDP) and the Combinatorial Extension (CE) algorithm for structural alignment. The approach involves the analysis of 3D-interaction patterns in DNA-protein interfaces, assignment of structural domains interacting with DNA, clustering of domains based on structural similarity and DNA-interacting patterns. Comparison with existing resources on describing structural and functional classifications of DNA-binding proteins was used to validate and improve the approach proposed here. In the course of our study we defined a set of criteria and heuristics allowing us to automatically build a biologically meaningful classification and define classes of functionally related protein domains. It was shown that taking into consideration interactions between protein domains and DNA considerably improves the classification accuracy. Our approach provides a high-throughput and up-to-date annotation of DNA-binding protein families which can be found at http://spdc.sdsc.edu.     

An automated method for modeling proteins on known templates using distance geometry.

We present an automated method incorporated into a software package, FOLDER, to fold a protein sequence on a given three-dimensional (3D) template. Starting with the sequence alignment of a family of homologous proteins, tertiary structures are modeled using the known 3D structure of one member of the family as a template. Homologous interatomic distances from the template are used as constraints. For nonhomologous regions in the model protein, the lower and the upper bounds for the interatomic distances are imposed by steric constraints and the globular dimensions of the template, respectively. Distance geometry is used to embed an ensemble of structures consistent with these distance bounds. Structures are selected from this ensemble based on minimal distance error criteria, after a penalty function optimization step. These structures are then refined using energy optimization methods. The method is tested by simulating the alpha-chain of horse hemoglobin using the alpha-chain of human hemoglobin as the template and by comparing the generated models with the crystal structure of the alpha-chain of horse hemoglobin. We also test the packing efficiency of this method by reconstructing the atomic positions of the interior side chains beyond C beta atoms of a protein domain from a known 3D structure. In both test cases, models retain the template constraints and any additionally imposed constraints while the packing of the interior residues is optimized with no short contacts or bond deformations. To demonstrate the use of this method in simulating structures of proteins with nonhomologous disulfides, we construct a model of murine interleukin (IL)-4 using the NMR structure of human IL-4 as the template. The resulting geometry of the nonhomologous disulfide in the model structure for murine IL-4 is consistent with standard disulfide geometry.     

Biochemical Characterization of the Cytochrome P450 CYP107CB2 from <Emphasis Type="Italic">Bacillus lehensis</Emphasis> G1

The bioconversion of vitamin D3 catalyzed by cytochrome P450 (CYP) requires 25-hydroxylation and subsequent 1α-hydroxylation to produce the hormonal activated 1α,25-dihydroxyvitamin D3. Vitamin D3 25-hydroxylase catalyses the first step in the vitamin D3 biosynthetic pathway, essential in the de novo activation of vitamin D3. A CYP known as CYP107CB2 has been identified as a novel vitamin D hydroxylase in Bacillus lehensis G1. In order to deepen the understanding of this bacterial origin CYP107CB2, its detailed biological functions as well as biochemical characteristics were defined. CYP107CB2 was characterized through the absorption spectral analysis and accordingly, the enzyme was assayed for vitamin D3 hydroxylation activity. CYP-ligand characterization and catalysis optimization were conducted to increase the turnover of hydroxylated products in an NADPH-regenerating system. Results revealed that the over-expressed CYP107CB2 protein was dominantly cytosolic and the purified fraction showed a protein band at approximately 62 kDa on SDS–PAGE, indicative of CYP107CB2. Spectral analysis indicated that CYP107CB2 protein was properly folded and it was in the active form to catalyze vitamin D3 reaction at C25. HPLC and MS analysis from a reconstituted enzymatic reaction confirmed the hydroxylated products were 25-hydroxyitamin D3 and 1α,25-dihydroxyvitamin D3 when the substrates vitamin D3 and 1α-hydroxyvitamin D3 were used. Biochemical characterization shows that CYP107CB2 performed hydroxylation activity at 25 °C in pH 8 and successfully increased the production of 1α,25-dihydroxyvitamin D3 up to four fold. These findings show that CYP107CB2 has a biologically relevant vitamin D3 25-hydroxylase activity and further suggest the contribution of CYP family to the metabolism of vitamin D3.     

Direct simulation of plastocyanin and cytochrome f interactions in solution

Most biological functions, including photosynthetic activity, are mediated by protein interactions. The proteins plastocyanin and cytochrome f are reaction partners in a photosynthetic electron transport chain. We designed a 3D computer simulation model of diffusion and interaction of spinach plastocyanin and turnip cytochrome f in solution. It is the first step in simulating the electron transfer from cytochrome f to photosystem 1 in the lumen of thylakoid. The model is multiparticle and it can describe the interaction of several hundreds of proteins. In our model the interacting proteins are represented as rigid bodies with spatial fixed charges. Translational and rotational motion of proteins is the result of the effect of stochastic Brownian force and electrostatic force. The Poisson-Boltzmann formalism is used to determine the electrostatic potential field generated around the proteins. Using this model we studied the kinetic characteristics of plastocyanin-cytochrome f complex formation for plastocyanin mutants at pH 7 and a variety of ionic strength values.     

Expression,structure-function,and molecular modeling of vitamin D P450s

Although vitamin D(3) is a natural product of a sunlight-mediated process in the skin, the secosteroid's biological function is dependent upon specific cytochrome P450 enzymes that mediate the parent vitamin's bioactivation and inactivation. Cytochrome P450C1 (CYP27B1) is the regulatory rate-limiting enzyme that directs the bioactivation process through introduction of a C-1alpha hydroxyl group. The resultant 1,25-dihydroxyvitamin D(3) (1,25D) is the biologically active secosteroid hormone that directs the multitude of vitamin D-dependent actions involved with calcium homeostasis, cellular differentiation and growth, and the immune response. The circulating and cellular level of 1,25D is regulated through a coordinated process involving the hormone's synthesis and degradation. Central to the degradation and turnover of 1,25D is the regulatory multi-catalytic cytochrome P450C24 (CYP24) enzyme that directs the introduction of C-24R groups onto targeted 25-hydroxy substrates. Discussed in this article is the action of the rat CYP24 to catalyze the side-chain oxidation and cleavage of 25-hydroxylated vitamin D metabolites. Expression and characterization of purified recombinant rat CYP24 is discussed in light of mutations directed at the enzyme's active site.     

Covalently linked heme in cytochrome p4504a fatty acid hydroxylases

Three independent experimental methods, liquid chromatography, denaturing gel electrophoresis with heme staining, and mass spectrometry, establish that the CYP4A class of enzymes has a covalently bound heme group even though the heme is not cross-linked to the protein in other P450 enzymes. Covalent binding has been demonstrated for CYP4A1, -4A2, -4A3, -4A8, and -4A11 heterologously expressed in Escherichia coli. However, the covalent link is also present in CYP4A1 isolated from rat liver and is not an artifact of heterologous expression. The extent of heme covalent binding in the proteins as isolated varies and is substoichiometric. In CYP4A3, the heme is attached to the protein via an ester link to glutamic acid residue 318, which is on the I-helix, and is predicted to be within the active site. This is the first demonstration that a class of cytochrome P450 enzymes covalently binds their prosthetic heme group.     

Prediction of amino acid residues participated in substrate recognition by cytochrome P450 subfamilies with broad substrate specificity

Cytochromes P450 comprise a large superfamily and several of their isoforms play a crucial role in metabolism of xenobiotics, including drugs. Although these enzymes demonstrate broad and cross‐substrate specificity, different cytochrome P450 subfamilies exhibit certain selectivity for some types of substrates. Analysis of amino acid residues of the active sites of six cytochrome subfamilies (CYP1А, CYP2А, CYP2С, CYP2D, CYP2E and CYP3А) enables to define subfamily‐specific patterns that consist of four residues. These residues are located on the periphery of the active sites of these cytochromes. We suggest that they can form a primary binding site at the entrance to the active site, defining cytochrome substrate recognition. Copyright © 2013 John Wiley & Sons, Ltd.     

In silico identification and screening of CYP24A1 inhibitors: 3D QSAR pharmacophore mapping and molecular dynamics analysis

Vitamin D is a key signalling molecule that plays a vital role in the regulation of calcium phosphate homeostasis and bone remodelling. The circulating biologically active form of vitamin D is regulated by the catabolic mechanism of cytochrome P450 24-hydroxylase (CYP24A1) enzyme. The over-expression of CYP24A1 negatively regulates the vitamin D level, which is the causative agent of chronic kidney disease, osteoporosis and several types of cancers. In this study, we found three potential lead molecules adverse to CYP24A1 through structure-based, atom-based pharmacophore and e-pharmacophore-based screening methods. Analysis was done by bioinformatics methods and tools like binding affinity (binding free energy), chemical reactivity (DFT studies) and molecular dynamics simulation (protein–ligand stability). Combined computational investigation showed that the compounds NCI_95001, NCI_382818 and UNPD_141613 may have inhibitory effects against the CYP24A1 protein.     

Fluorescent photoaffinity labeling of cytochrome P450 3A4 by lapachenole: identification of modification sites by mass spectrometry

While photoaffinity ligands (PALs) have been widely used to probe the structures of many receptors and transporters, their effective use in the study of membrane-bound cytochrome P450s is less established. Here, lapachenole has been used as an effective photoaffinity ligand of human P450 3A4, and mass spectrometry data demonstrating the efficient and specific photoaffinity labeling of CYP3A4 by this naturally occurring benzochromene compound is presented. Without photolysis, lapachenole is a substrate of CYP3A4 and can be metabolized to hydroxylated products by this enzyme. A high-performance liquid chromatography/electrospray ionization mass spectrometry (HPLC/ESI-MS) procedure was developed to analyze small amounts of intact purified CYP3A4, and analysis of the labeled protein showed the presence of one molecule of lapachenole bound per monomer of protein. Photolabeled CYP3A4 peptide adducts were further characterized by mass spectrometric analysis after proteolytic digestion and isolation of fluorescent photolabeled peptides. Two peptide adducts accounting for >95% of the labeled peptides were isolated by HPLC, and both peptides, ECYSVFTNR (positions 97-105) and VLQNFSFKPCK (positions 459-469), were identified by nano-LC/ESI quadrupole time-of-flight (QTOF) and matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry. The sites of modification were further localized to positions Cys-98 and Cys-468 for each peptide by nano-LC/ESI QTOF tandem mass spectrometry (MS/MS). The results provided the first direct evidence for interaction between the PAL and the putative B-B' loop region, which may serve as a substrate access channel or as a part of the CYP3A4 active site. In conclusion, benzochromene analogues are effective PALs, which may be used in the study of other cytochrome P450 structures.     

Metabolism of cholesterol, vitamin D3 and 20-hydroxyvitamin D3 incorporated into phospholipid vesicles by human CYP27A1

CYP27A1 is a mitochondrial cytochrome P450 which can hydroxylate vitamin D3 and cholesterol at carbons 25 and 26, respectively. The product of vitamin D3 metabolism, 25-hydroxyvitamin D3, is the precursor to the biologically active hormone, 1α,25-dihydroxyvitamin D3. CYP27A1 is attached to the inner mitochondrial membrane and substrates appear to reach the active site through the membrane phase. We have therefore examined the ability of bacterially expressed and purified CYP27A1 to metabolize substrates incorporated into phospholipid vesicles which resemble the inner mitochondrial membrane. We also examined the ability of CYP27A1 to metabolize 20-hydroxyvitamin D3 (20(OH)D3), a novel non-calcemic form of vitamin D derived from CYP11A1 action on vitamin D3 which has anti-proliferative activity on keratinocytes, leukemic and myeloid cells. CYP27A1 displayed high catalytic activity towards cholesterol with a turnover number (k(cat)) of 9.8 min(-1) and K(m) of 0.49 mol/mol phospholipid (510 μM phospholipid). The K(m) value of vitamin D3 was similar for that of cholesterol, but the k(cat) was 4.5-fold lower. 20(OH)D3 was metabolized by CYP27A1 to two major products with a k(cat)/K(m) that was 2.5-fold higher than that for vitamin D3, suggesting that 20(OH)D3 could effectively compete with vitamin D3 for catalysis. NMR and mass spectrometric analyses revealed that the two major products were 20,25-dihydroxyvitamin D3 and 20,26-dihydroxyvitamin D3, in almost equal proportions. Thus, the presence of the 20-hydroxyl group on the vitamin D3 side chain enables it to be metabolized more efficiently than vitamin D3, with carbon 26 in addition to carbon 25 becoming a major site of hydroxylation. Our study reports the highest k(cat) for the 25-hydroxylation of vitamin D3 by any human cytochrome P450 suggesting that CYP27A1 might be an important contributor to the synthesis of 25-hydroxyvitamin D3, particularly in tissues where it is highly expressed.