Computational Insight Into Vitamin K1 ω-Hydroxylation by Cytochrome P450 4F2

The fatty acid omega hydroxylase genes (CYP4 family) in the progression of metabolic dysfunction-associated steatotic liver disease (MASLD): An RNA sequence database analysis and review

Fatty acid omega hydroxylase P450s consist of enzymes that hydroxylate various chain-length saturated and unsaturated fatty acids (FAs) and bioactive eicosanoid lipids. The human cytochrome P450 gene 4 family (CYP4) consists of 12 members that are associated with several human diseases. However, their role in the progression of metabolic dysfunction-associated fatty liver disease (MASLD) remains largely unknown. It has long been thought that the induction of CYP4 family P450 during fasting and starvation prevents FA-related lipotoxicity through FA metabolism to dicarboxylic acids that are chain-shortened in peroxisomes and then transported to the mitochondria for complete oxidation. Several studies have revealed that peroxisome succinate transported to the mitochondria is used for gluconeogenesis during fasting and starvation, and recent evidence suggests that peroxisome acetate can be utilized for lipogenesis and lipid droplet formation as well as epigenetic modification of gene transcription. In addition, omega hydroxylation of the bioactive eicosanoid arachidonic acid to 20-Hydroxyeicosatetraenoic acid (20-HETE) is essential for activating the GPR75 receptor, leading to vasoconstriction and cell proliferation. Several mouse models of diet-induced MASLD have revealed the induction of selective CYP4A members and the suppression of CYP4F during steatosis and steatohepatitis, suggesting a critical metabolic role in the progression of fatty liver disease. Thus, to further investigate the functional roles of CYP4 genes, we analyzed the differential gene expression of 12 members of CYP4 gene family in datasets from the Gene Expression Omnibus (GEO) from patients with steatosis, steatohepatitis, fibrosis, cirrhosis, and hepatocellular carcinoma. We also observed the differential expression of various CYP4 genes in the progression of MASLD, indicating that different CYP4 members may have unique functional roles in the metabolism of specific FAs and eicosanoids at various stages of fatty liver disease. These results suggest that targeting selective members of the CYP4A family is a viable therapeutic approach for treating and managing MASLD.

Functional Characterization of 29 Cytochrome P450 4F2 Variants Identified in a Population of 8380 Japanese Subjects and Assessment of Arachidonic Acid ω-Hydroxylation

Cytochrome P450 4F2 (CYP4F2) is an enzyme that is involved in the metabolism of arachidonic acid (AA), vitamin E and K, and xenobiotics including drugs. CYP4F2*3 polymorphism (rs2108622; c.1297G>A; p.Val433Met) has been associated with hypertension, ischemic stroke, and variation in the effectiveness of the anticoagulant drug warfarin. In this study, we characterized wild-type CYP4F2 and 28 CYP4F2 variants, including a Val433Met substitution, detected in 8380 Japanese subjects. The CYP4F2 variants were heterologously expressed in 293FT cells to measure the concentrations of CYP4F2 variant holoenzymes using carbon monoxide-reduced difference spectroscopy, where the wild type and 18 holoenzyme variants showed a peak at 450 nm. Kinetic parameters [Vmax , substrate concentration producing half of Vmax (S50 ), and intrinsic clearance (CL int ) as Vmax /S50 ] of AA ω-hydroxylation were determined for the wild type and 21 variants with enzyme activity. Compared with the wild type, two variants showed significantly decreased CL int values for AA ω-hydroxylation. The values for seven variants could not be determined because no enzymatic activity was detected at the highest substrate concentration used. Three-dimensional structural modeling was performed to determine the reason for reduced enzymatic activity of the CYP4F2 variants. Our findings contribute to a better understanding of CYP4F2 variant-associated diseases and possible future therapeutic strategies. SIGNIFICANCE STATEMENT: CYP4F2 is involved in the metabolism of arachidonic acid and vitamin K, and CYP4F2*3 polymorphisms have been associated with hypertension and variation in the effectiveness of the anticoagulant drug warfarin. This study presents a functional analysis of 28 CYP4F2 variants identified in Japanese subjects, demonstrating that seven gene polymorphisms cause loss of CYP4F2 function, and proposes structural changes that lead to altered function.

Computational analysis of missense variant CYP4F2*3 (V433M) in association with human CYP4F2 dysfunction: a functional and structural impact

BACKGROUND

Cytochrome P450 4F2 (CYP4F2) enzyme is a member of the CYP4 family responsible for the metabolism of fatty acids, therapeutic drugs, and signaling molecules such as arachidonic acid, tocopherols, and vitamin K. Several reports have demonstrated that the missense variant CYP4F2*3 (V433M) causes decreased activity of CYP4F2 and inter-individual variations in warfarin dose in different ethnic groups. However, the molecular pathogenicity mechanism of missense V433M in CYP4F2 at the atomic level has not yet been completely elucidated.

METHODS AND RESULTS

In the current study, we evaluated the effect of the V433M substitution on CYP4F2 using 14 different bioinformatics tools. Further molecular dynamics (MD) simulations were performed to assess the impact of the V433M mutation on the CYP4F2 protein structure, stability, and dynamics. In addition, molecular docking was used to illustrate the effect of V433M on its interaction with vitamin K1. Based on our results, the CYP4F2*3 variant was a damaging amino acid substitution with a destabilizing nature. The simulation results showed that missense V433M affects the dynamics and stability of CYP4F2 by reducing its compactness and stability, which means that it tends to change the overall structural conformation and flexibility of CYP4F2. The docking results showed that the CYP4F2*3 variant decreased the binding affinity between vitamin K1 and CYP4F2, which reduced the activity of CYP4F2*3 compared to native CYP4F2.

CONCLUSIONS

This study determined the molecular pathogenicity mechanism of the CYP4F2*3 variant on the human CYP4F2 protein and provided new information for understanding the structure-function relationship of CYP4F2 and other CYP4 enzymes. These findings will aid in the development of effective drugs and treatment options.

Insights into the binding mechanism between α-TOH and CYP4F2: A homology modeling, molecular docking, and molecular dynamics simulation study

α-Tocopherol (α-TOH) is a potent antioxidant. The concentrations of α-TOH in plasma are closely related to human health. α-TOH can be regulated by the metabolism of cytochrome P450 4F2 (CYP4F2). However, the atomic-level basis for this regulation process remains elusive. Here, we successfully constructed the structure of CYP4F2 by homology modeling and obtained the α-TOH-CYP4F2 complex models using molecular docking. Three parallel 500 ns molecular dynamics simulations were performed on each complex model to investigate the details of the interaction between α-TOH and CYP4F2. MM-GBSA method combined with principal component analysis shows that 8 key residues establish a hydrophobic cavity stabilizing α-TOH in the pocket of CYP4F2 and S423 forms an important hydrogen bond with α-TOH anchoring α-TOH in the favorable position for ω-hydroxylation. Based on our simulation results and the experimental facts, we designed mutation simulation experiments to clarify the important role of two key residues (S423 and V433) in the binding of α-TOH with CYP4F2. The results show that the mutations directly or indirectly change the binding mode of α-TOH and decrease its binding affinity with CYP4F2, which is unfavorable for ω-hydroxylation. Our results could enrich the information on structure-function relationships of CYP4F2 and provide valuable insights into the regulatory mechanism of CYP4F2 on the metabolism of α-TOH.

Structural insights into understudied human cytochrome P450 enzymes

Human cytochrome P450 (CYP) enzymes are widely known for their pivotal role in the metabolism of drugs and other xenobiotics as well as of endogenous chemicals. In addition, CYPs are involved in numerous pathophysiological pathways and, hence, are therapeutically relevant. Remarkably, a portion of promising CYP targets is still understudied and, as a consequence, untargeted, despite their huge therapeutic potential. An increasing number of X-ray and cryo-electron microscopy (EM) structures for CYPs have recently provided new insights into the structural basis of CYP function and potential ligand binding. This structural knowledge of CYP functionality is essential for both understanding metabolism and exploiting understudied CYPs as drug targets. In this review, we summarize and highlight structural knowledge about this enzyme class, with a focus on understudied CYPs and resulting opportunities for structure-based drug design. Teaser: This review summarizes recent structural insights into understudied cytochrome P450 enzymes. We highlight the impact of molecular modeling for mechanistically explaining pathophysiological effects establishing understudied CYPs as promising drug targets.

Deciphering Reaction Determinants of Altered-Activity CYP2D6 Variants by Well-Tempered Metadynamics Simulation and QM/MM Calculations

The xenobiotic metabolizing enzyme CYP2D6 is the P450 cytochrome family member with the highest rate of polymorphism. This causes changes in the enzyme activity and specificity, which can ultimately lead to adverse reactions during drug treatment. To avoid or lower CYP-related toxicity risks, prediction of the most likely positions within a molecule where a metabolic reaction might occur is paramount. In order to obtain accurate predictions, it is crucial to understand all phenomena within the active site of the enzyme that contribute to an efficient substrate recognition and the subsequent catalytic reaction together with their relative weight within the overall thermodynamic context. This study aims to define the weight of the driving forces upon the C-H bond activation within CYP2D6 wild-type and a clinically relevant allelic variant with increased activity (CYP2D6*53) featuring two amino acid mutations in close vicinity of the heme. First, we investigated the steric and electrostatic complementarity of the substrate bufuralol using well-tempered metadynamics simulations with the aim to obtain the free energy profiles for each site of metabolism (SoM) within the different active sites. Second, the stereoelectronic complementarity was determined for each SoM within the two different active-site environments. Relying on the well-tempered metadynamics simulation energy profiles of each SoM, we identified the binding mode that was closest to the preferred transition-state geometry for efficient C-H bond activation. The binding modes were then used as starting structures for the quantum mechanics/molecular mechanics calculations performed to quantify the corresponding activation barriers. Our results show the relevance of the steric component in orienting the SoM in an energetically accessible position toward the heme. However, the corresponding intrinsic reactivity and electronic complementarity within the active site must be accurately evaluated in order to obtain a meaningful reaction prediction, from which the predominant SoM can be determined. The F120I mutation lowered the activation barrier for the major site and one of the minor SoMs. However, it had an impact neither on the CYP2D6 enantioselectivity preference of the oxidation reaction nor on the stereoselectivity from the substrate point of view.