Structural Basis for the Mutation-Induced Dysfunction of Human CYP2J2: A Computational Study

CYP2J2-mediated metabolism of arachidonic acid in heart: A review of its kinetics, inhibition and role in heart rhythm control

Cytochrome P450 2 J2 (CYP2J2) is primarily expressed extrahepatically and is the predominant epoxygenase in human cardiac tissues. This highlights its key role in the metabolism of endogenous substrates. Significant scientific interest lies in cardiac CYP2J2 metabolism of arachidonic acid (AA), an omega-6 polyunsaturated fatty acid, to regioisomeric bioactive epoxyeicosatrienoic acid (EET) metabolites that show cardioprotective effects including regulation of cardiac electrophysiology. From an in vitro perspective, the accurate characterization of the kinetics of CYP2J2 metabolism of AA including its inhibition and inactivation by drugs could be useful in facilitating in vitro-in vivo extrapolations to predict drug-AA interactions in drug discovery and development. In this review, background information on the structure, regulation and expression of CYP2J2 in human heart is presented alongside AA and EETs as its endogenous substrate and metabolites. The in vitro and in vivo implications of the kinetics of this endogenous metabolic pathway as well as its perturbation via inhibition and inactivation by drugs are elaborated. Additionally, the role of CYP2J2-mediated metabolism of AA to EETs in cardiac electrophysiology will be expounded.

A Case of Symptomatic Supratherapeutic International Normalized Ratio on Rivaroxaban

Rivaroxaban is a direct oral anticoagulant that works by inhibiting factor Xa. Direct anticoagulants have largely replaced direct vitamin K inhibitors (VKAs) due to the decreased risk of major hemorrhages and the lack of need for regular monitoring and dose adjustments. However, there have been multiple reports of elevated international normalized ratio (INR) and bleeding incidents in patients on rivaroxaban, which brings into question the potential need for monitoring. We report a case of an INR of 4.8 in a rivaroxaban-naïve patient who presented with gastrointestinal bleeding and a significant drop in hemoglobin four days after starting rivaroxaban. We present possible pharmacologic explanations. We propose the idea that specific subgroups of patients may be at risk for true INR elevations and may benefit from routine monitoring of their INR while on rivaroxaban.

Cardiomyocyte-specific CYP2J2 and its therapeutic implications

INTRODUCTION

Cytochrome P450s (CYPs) are a superfamily of monooxygenases with diverse biological roles. CYP2J2 is an isozyme highly expressed in the heart where it metabolizes endogenous substrates such as N-3/N-6 polyunsaturated fatty acids (PUFA) to produce lipid mediators involved in homeostasis and cardioprotective responses. Expanding our knowledge of the role CYP2J2 has within the heart is important for understanding its impact on cardiac health and disease.

AREAS COVERED

The objective of this review was to assess the state of knowledge regarding cardiac CYP2J2. A literature search was conducted using PubMed-MEDLINE (from 2022 and earlier) to evaluate relevant studies regarding CYP2J2 mediated cardioprotection, small molecule modulators, effects of CYP2J2 substrates toward biologically relevant effects and implications of CYP2J2 polymorphisms and sexual dimorphism in the heart.

EXPERT OPINION

Cardiac CYP2J2-mediated metabolism of endogenous and exogenous substrates have been shown to impact cardiac function. Identifying individual factors, like sex and age, that affect CYP2J2 require further elucidation to better understand CYP2J2's clinical relevance. Resolving the biological targets and activities of CYP2J2-derived PUFA metabolites will be necessary to safely target CYP2J2 and design novel analogues. Targeting CYP2J2 for therapeutic aims offers a potential novel approach to regulating cardiac homeostasis, drug metabolism and cardioprotection.

Atypical kinetics of cytochrome P450 2J2: Epoxidation of arachidonic acid and reversible inhibition by xenobiotic inhibitors

Extrahepatic CYP2J2 metabolism of arachidonic acid (AA) to bioactive regioisomeric epoxyeicosatrienoic acids (EETs) is implicated in both physiological and pathological conditions. Here, we aimed to characterize atypical substrate inhibition kinetics of this endogenous metabolic pathway and its reversible inhibition by xenobiotic inhibitors when AA is used as the physiologically-relevant substrate vis-à-vis conventional probe substrate astemizole (AST). As compared to typical Michaelis-Menten kinetics observed for AST, complete substrate inhibition was observed for CYP2J2 metabolism of AA to 14,15-EET whereby velocity of the reaction declined significantly at concentrations of AA above 20-30 µM with an estimated substrate inhibition constant (K_s) of 31 µM. In silico sequential docking of two AA substrates to orthosteric (OBS) and adjacent secondary binding sites (SBS) within a 3-dimensional homology model of CYP2J2 revealed favorable and comparable binding poses of glide-scores -3.1 and -3.8 respectively. Molecular dynamics (MD) simulations ascertained CYP2J2 conformational stability with dual AA substrate binding as time-dependent root mean squared deviation (RMSD) of protein Cα atoms and ligand heavy atoms stabilized to a plateau in all but one trajectory (n=6). The distance between heme-iron and ω6 (C14, C15) double bond of AA in OBS also increased from 7.5 ± 1.4 Å to 8.5 ± 1.8 Å when CYP2J2 was simulated with only AA in OBS versus the presence of AA in both OBS and SBS (p<0.001), supporting the observed in vitro substrate inhibition phenomenon. Poor correlation was observed between inhibitory constants (K_i) determined for a panel of nine competitive and mixed mode xenobiotic inhibitors against CYP2J2 metabolism of AA as compared to AST, whereby 4 out of 9 drugs had a greater than 5-fold difference between K_i values. Nonlinear Eadie-Hofstee plots illustrated that complete substrate inhibition of CYP2J2 by AA was not attenuated even at high concentrations of xenobiotic inhibitors which further corroborates that CYP2J2 may accommodate three or more ligands simultaneously. In light of the atypical kinetics, our results highlight the importance of using physiologically-relevant substrates in in vitro enzymatic inhibition assays for the characterization of xenobiotic-endobiotic interactions which is applicable to other complex endogenous metabolic pathways beyond CYP2J2 metabolism of AA to EETs. The accurate determination of K_i would further facilitate the association of xenobiotic-endobiotic interactions to observed therapeutic or toxic outcomes.

Molecular determinant of substrate binding and specificity of cytochrome P450 2J2

Cytochrome P450 2J2 (CYP2J2) is responsible for the epoxidation of endogenous arachidonic acid, and is involved in the metabolism of exogenous drugs. To date, no crystal structure of CYP2J2 is available, and the proposed structural basis for the substrate recognition and specificity in CYP2J2 varies with the structural models developed using different computational protocols. In this study, we developed a new structural model of CYP2J2, and explored its sensitivity to substrate binding by molecular dynamics simulations of the interactions with chemically similar fluorescent probes. Our results showed that the induced-fit binding of these probes led to the preferred active poses ready for the catalysis by CYP2J2. Divergent conformational dynamics of CYP2J2 due to the binding of each probe were observed. However, a stable hydrophobic clamp composed of residues I127, F310, A311, V380, and I487 was identified to restrict any substrate access to the active site of CYP2J2. Molecular docking of a series of compounds including amiodarone, astemizole, danazol, ebastine, ketoconazole, terfenadine, terfenadone, and arachidonic acid to CYP2J2 confirmed the role of those residues in determining substrate binding and specificity of CYP2J2. In addition to the flexibility of CYP2J2, the present work also identified other factors such as electrostatic potential in the vicinity of the active site, and substrate strain energy and property that have implications for the interpretation of CYP2J2 metabolism.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

Impact of peripheral mutations on the access channels of human cytochrome P450 1A2

Communicated by Ramaswamy H. Sarma.

Role of Arginine 117 in Substrate Recognition by Human Cytochrome P450 2J2

The influence of Arginine 117 of human cytochrome P450 2J2 in the recognition of ebastine and a series of terfenadone derivatives was studied by site-directed mutagenesis. R117K, R117E, and R117L mutants were produced, and the behavior of these mutants in the hydroxylation of ebastine and terfenadone derivatives was compared to that of wild-type CYP2J2. The data clearly showed the importance of the formation of a hydrogen bond between R117 and the keto group of these substrates. The data were interpreted on the basis of 3D homology models of the mutants and of dynamic docking of the substrates in their active site. These modeling studies also suggested the existence of a R117-E222 salt bridge between helices B’ and F that would be important for maintaining the overall folding of CYP2J2.

Discovery of a regioselectivity switch in nitrating P450s guided by molecular dynamics simulations and Markov models

The dynamic motions of protein structural elements, particularly flexible loops, are intimately linked with diverse aspects of enzyme catalysis. Engineering of these loop regions can alter protein stability, substrate binding, and even dramatically impact enzyme function. When these flexible regions are structurally unresolvable, computational reconstruction in combination with large-scale molecular dynamics simulations can be used to guide the engineering strategy. Here, we present a collaborative approach consisting of both experiment and computation that led to the discovery of a single mutation in the F/G loop of the nitrating cytochrome P450 TxtE that simultaneously controls loop dynamics and completely shifts the enzyme's regioselectivity from the C4 to the C5 position of L-tryptophan. Furthermore, we find that this loop mutation is naturally present in a subset of homologous nitrating P450s and confirm that these uncharacterized enzymes exclusively produce 5-nitro-L-tryptophan, a previously unknown biosynthetic intermediate.

Bridging scales through multiscale modeling: a case study on protein kinase A

The goal of multiscale modeling in biology is to use structurally based physico-chemical models to integrate across temporal and spatial scales of biology and thereby improve mechanistic understanding of, for example, how a single mutation can alter organism-scale phenotypes. This approach may also inform therapeutic strategies or identify candidate drug targets that might otherwise have been overlooked. However, in many cases, it remains unclear how best to synthesize information obtained from various scales and analysis approaches, such as atomistic molecular models, Markov state models (MSM), subcellular network models, and whole cell models. In this paper, we use protein kinase A (PKA) activation as a case study to explore how computational methods that model different physical scales can complement each other and integrate into an improved multiscale representation of the biological mechanisms. Using measured crystal structures, we show how molecular dynamics (MD) simulations coupled with atomic-scale MSMs can provide conformations for Brownian dynamics (BD) simulations to feed transitional states and kinetic parameters into protein-scale MSMs. We discuss how milestoning can give reaction probabilities and forward-rate constants of cAMP association events by seamlessly integrating MD and BD simulation scales. These rate constants coupled with MSMs provide a robust representation of the free energy landscape, enabling access to kinetic, and thermodynamic parameters unavailable from current experimental data. These approaches have helped to illuminate the cooperative nature of PKA activation in response to distinct cAMP binding events. Collectively, this approach exemplifies a general strategy for multiscale model development that is applicable to a wide range of biological problems.

Lipid-metabolizing CYPs in the regulation and dysregulation of metabolism

The cytochrome P450s (CYPs) represent a highly divergent class of enzymes involved in the oxidation of organic compounds. A subgroup of CYPs metabolize ω3-arachidonic and linoleic acids and ω6-docosahexaenoic and eicosapentaenoic polyunsaturated fatty acids (PUFAs) into a series of related biologically active mediators. Over the past 20 years, increasing evidence has emerged for a role of these PUFA-derived mediators in physiological and pathophysiological processes in the vasculature, during inflammation, and in the regulation of metabolism. With recent technological advances and increased availability of lipid mass spectroscopy, we are now starting to discern the patterns of these CYP-PUFA products in health and disease. These analyses not only are revealing the diverse spectrum of lipid nutrients regulated by CYPs, but also clearly indicate that the balance of these mediators changes with dietary intake of different PUFA classes. These findings suggest that we are only just beginning to understand all of the relevant lipid species produced by CYP pathways. Moreover, we are still a long way from understanding the nature and presence of their receptors, their tissue expression, and the pathophysiological processes they regulate. This review highlights these future issues in the context of lipid-metabolizing CYP enzymes, focusing particularly on the CYP450 family of epoxygenases and the lipid mediators they produce.

Advances in molecular modeling of human cytochrome P450 polymorphism

Cytochrome P450 (CYP) is a supergene family of metabolizing enzymes involved in the phase I metabolism of drugs and endogenous compounds. CYP oxidation often leads to inactive drug metabolites or to highly toxic or carcinogenic metabolites involved in adverse drug reactions (ADR). During the last decade, the impact of CYP polymorphism in various drug responses and ADR has been demonstrated. Of the drugs involved in ADR, 56% are metabolized by polymorphic phase I metabolizing enzymes, 86% among them being CYP. Here, we review the major CYP polymorphic forms, their impact for drug response and current advances in molecular modeling of CYP polymorphism. We focus on recent studies exploring CYP polymorphism performed by the use of sequence-based and/or protein-structure-based computational approaches. The importance of understanding the molecular mechanisms related to CYP polymorphism and drug response at the atomic level is outlined.