Structure and dynamics of the membrane-bound cytochrome P450 2C9

Insight into the effect of ibuprofen on the permeability of the membrane: a molecular dynamic simulation study

Studying interactions between drugs and cell membranes is of great interest to designing novel drugs, optimizing drug delivery, and discerning drug mechanism action. In this study, we investigated the physical properties of the bilayer membrane model of POPC upon interaction with ibuprofen (IBU) using molecular dynamics simulations. The area per lipid (APL) was calculated to describe the effect of ibuprofen on the packing properties of the lipid bilayer. The APL was 0.58 nm2 and 0.63 nm2 for the membrane in low and high IBU respectively, and 0.57 nm2 for the membrane without IBU. Our finding showed that the mean square deviation (MSD) increased with increased ibuprofen content. In addition, the order parameter for the hydrocarbon chain of lipids increased with increased ibuprofen content. There was an increment in the transfer free energy after the head group region while it was maximum in the hydrophobic core for hydrogen peroxide (H2O2) (∼6.2 kcal.mol-1) and H2O (∼3.4 kcal.mol-1) which then decreased to respective values of (∼4.6 kcal.mol-1), and (∼2.3 kcal.mol-1) at the center of the bilayer in the presence of IBU. It seems that in the presence of ibuprofen, the free energy profile of the permeability of water and H2O2 significantly decreased. These findings show that ibuprofen significantly influences the physical properties of the bilayer by decreasing the packing and intermolecular interaction in the hydrocarbon chain region and increasing the water permeability of the bilayer. These results may provide insights into the local cytotoxic side effects of ibuprofen and its underlying molecular mechanisms.Communicated by Ramaswamy H. Sarma.

Multiresolution molecular dynamics simulations reveal the interplay between conformational variability and functional interactions in membrane-bound cytochrome P450 2B4

Cytochrome P450 2B4 (CYP 2B4) is one of the best-characterized CYPs and serves as a key model system for understanding the mechanisms of microsomal class II CYPs, which metabolize most known drugs. The highly flexible nature of CYP 2B4 is apparent from crystal structures that show the active site with either a wide open or a closed heme binding cavity. Here, we investigated the conformational ensemble of the full-length CYP 2B4 in a phospholipid bilayer, using multiresolution molecular dynamics (MD) simulations. Coarse-grained MD simulations revealed two predominant orientations of CYP 2B4's globular domain with respect to the bilayer. Their refinement by atomistic resolution MD showed adaptation of the enzyme's interaction with the lipid bilayer, leading to open configurations that facilitate ligand access to the heme binding cavity. CAVER analysis of enzyme tunnels, AquaDuct analysis of water routes, and Random Acceleration Molecular Dynamics simulations of ligand dissociation support the conformation-dependent passage of molecules between the active site and the protein surroundings. Furthermore, simulation of the re-entry of the inhibitor bifonazole into the open conformation of CYP 2B4 resulted in binding at a transient hydrophobic pocket within the active site cavity that may play a role in substrate binding or allosteric regulation. Together, these results show how the open conformation of CYP 2B4 facilitates the binding of substrates from and release of products to the membrane, whereas the closed conformation prolongs the residence time of substrates or inhibitors and selectively allows the passage of smaller reactants via the solvent and water channels.

Effect of Lipid Bilayer Anchoring on the Conformational Properties of the Cytochrome P450 2D6 Binding Site

Human cytochrome P450 (CYP) proteins metabolize 75% of small-molecule pharmaceuticals, which makes structure-based modeling of CYP metabolism and inhibition, bolstered by the current availability of X-ray crystal structures of CYP globular catalytic domains, an attractive prospect. Accounting for this broad metabolic capacity is a combination of the existence of multiple different CYP proteins and the capacity of a single CYP protein to metabolize multiple different small molecules. It is thought that structural plasticity and flexibility contribute to this latter property; therefore, incorporating diverse conformational states of a particular CYP is likely an important consideration in structure-based CYP metabolism and inhibition modeling. While all-atom explicit-solvent molecular dynamics simulations can be used to generate conformational ensembles under biologically relevant conditions, existing CYP crystal structures are of the globular domain only, whereas human CYPs contain N-terminal transmembrane and linker peptides that anchor the globular catalytic domain to the endoplasmic reticulum. To determine whether this can cause significant differences in the sampled binding site conformations, microsecond scale all-atom explicit-solvent molecular dynamics simulations of the CYP2D6 globular domain in an aqueous environment were compared with those of the full-length protein anchored in a POPC lipid bilayer. While bilayer-anchoring damped some structural fluctuations in the globular domain relative to the aqueous simulations, none of the affected residues included binding site pocket residues. Furthermore, clustering of molecular dynamics snapshots based on either pairwise binding site pocket RMSD or volume differences demonstrated a lack of separation of snapshots from the two simulation conditions into different clusters. These results suggest the substantially simpler and computationally cheaper aqueous simulation approach can be used to generate a relevant conformational ensemble of the CYP2D6 binding site for structure-based metabolism and inhibition modeling.

Multilevel superposition for deciphering the conformational variability of protein ensembles

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

Equilibrium landscape of ingress/egress channels and gating residues of the Cytochrome P450 3A4

The Cytochrome P450 (CYP) enzymes metabolize a variety of drugs, which may potentially lead to toxicity or reduced efficacy when drugs are co-administered. These drug-drug interactions are often manifested by CYP3A4, the most prevalent of all CYP isozymes. We carried out multiple MD simulations employing CAVER to quantify the channels, and Hidden Markov Models (HMM) to characterize the behavior of the gating residues. We discuss channel properties, bottleneck residues with respect to their likelihood to deem the respective channel ingress or egress, gating residues regarding their open or closed states, and channel location relative to the membrane. Channels do not display coordinated motion and randomly transition between different conformations. Gateway residues also behave in a random fashion. Our findings shed light on the equilibrium behavior of the gating residues and channels in the apo state.

Direct interaction between the transmembrane helices stabilize cytochrome P450 2B4 and cytochrome b5 redox complex

The catalytic activity of cytochrome P450 2B4 (CYP2B4) is moderated by its cognate redox partner cytochrome b5 (Cyt-b5). The endoplasmic reticulum (ER) membrane and intermolecular transmembrane (TM) interaction between CYP2B4 and Cyt-b5 regulate the substrate catalysis and the reaction rate. This emphasizes the significance of elucidating the molecular basis of CYP2B4 and Cyt-b5 complexation in a membrane environment to better understand the enzymatic activity of CYP2B4. Our previous solid-state NMR studies revealed the membrane topology of the transmembrane domains of these proteins in the free and complex forms. Here, we show the cross-angle complex formation by the single-pass TM domains of CYP2B4 and Cyt-b5, which is mainly driven by several salt-bridges (E2-R128, R21-D104 and K25-D104), using a multi-microsecond molecular dynamic simulation. Additionally, the leucine-zipper residues (L8, L12, L15, L18 and L19 from CYP2B4) and π-stacking between H23 and F20 residues of CYP2B4 and W110 of Cyt-b5 are identified to stabilize the TM-TM complex in the ER membrane. The simulated tilts of the helices in the free and in the complex are in excellent agreement with solid-state NMR results. The TM-TM packing influences a higher order structural stability when compared to the complex formed by the truncated soluble domains of these two proteins. MM/PBSA based binding free energy estimates nearly 100-fold higher binding affinity (ΔG = -2810.68 ± 696.44 kJ/mol) between the soluble domains of the full-length CYP2B4 and Cyt-b5 when embedded in lipid membrane as compared to the TM-domain-truncated soluble domains (ΔG = -27.406 ± 10.32 kJ/mol). The high-resolution full-length CYP2B4-Cyt-b5 complex structure and its dynamics in a native ER membrane environment reported here could aid in the development of approaches to effectively modulate the drug-metabolism activity of CYP2B4.

Phenylalanine Residues in the Active Site of CYP2E1 Participate in Determining the Binding Orientation and Metabolism-Dependent Genotoxicity of Aromatic Compounds

The composition of amino acids forming the active site of a CYP enzyme is impactful in its substrate selectivity. For CYP2E1, the role of PHE residues in the formation of effective binding orientations for its aromatic substrates remains unclear. In this study, molecular docking and molecular dynamics analysis were performed to reflect the interactions between PHEs in the active site of human CYP2E1 and various aromatic compounds known as its substrates. The results indicated that the orientation of 1-methylpyrene (1-MP) in the active site was highly determined by the presence of PHEs, PHE478 contributing to the binding free energy most significantly. Moreover, by building a random forest model the relationship between each of 19 molecular descriptors of polychlorinated biphenyl (PCB) compounds (from molecular docking, quantum mechanics, and physicochemical properties) and their human CYP2E1-dependent mutagenicityas established mostly in our lab, was investigated. The presence of PHEs did not appear to significantly modify the electronic or structural feature of each bound ligand (PCB), instead, the flexibility of the conformation of PHEs contributed substantially to the effective binding energy and orientation. It is supposed that PHE residues adjust their own conformation to permit a suitablly shaped cavity for holding the ligand and forming its orientation as favorable for a biochemical reaction. This study has provided some insights into the role of PHEs in guiding the interactive adaptation of the active site of human CYP2E1 for the binding and metabolism of aromatic substrates.

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

Human Orphan Cytochrome P450 2U1 Catalyzes the ω-Hydroxylation of Leukotriene B4

Cytochrome P450 2U1 (CYP2U1) identified from the human genome remains poorly known since few data are presently available on its physiological function(s) and substrate(s) specificity. CYP2U1 mutations are associated with complicated forms of hereditary spastic paraplegia, alterations of mitochondrial architecture and bioenergetics. In order to better know the biological roles of CYP2U1, we used a bioinformatics approach. The analysis of the data invited us to focus on leukotriene B₄ (LTB₄), an important inflammatory mediator. Here, we show that CYP2U1 efficiently catalyzes the hydroxylation of LTB₄ predominantly on its ω-position. We also report docking experiments of LTB₄ in a 3D model of truncated CYP2U1 that are in agreement with this hydroxylation regioselectivity. The involvement of CYP2U1 in the metabolism of LTB₄ could have strong physiological consequences in cerebral pathologies including ischemic stroke because CYP2U1 is predominantly expressed in the brain.

Cholesterol's inhibition effect on entering of chlorzoxazone into phosphatidylethanolamine bilayer: Relevance to cytochrome P450 drug metabolism at endoplasmic reticulum membranes

Many drugs are metabolized by cytochrome P450 (CYP) in the endoplasmic reticulum (ER) membrane. Recent studies have shown that CYP-substrate drugs reach the CYP active site after entering the lipid hydrophobic part of the ER membrane. To clarify the role of cholesterol (Chol) in the CYP-related drug metabolic process, we investigated the lipid bilayer entry of CYP-substrate drugs using a model membrane system as follows. The model membrane system comprised palmitoyl-oleoyl-phosphatidylethanolamine (POPE) and Chol. Phosphatidylethanolamine is the second major phospholipid component of ER membranes. Chlorzoxazone (CZX) was used as the CYP-substrate drug. Calorimetric measurements showed that the addition of CZX to POPE bilayers decreased the gel-liquid crystal phase transition temperature; X-ray diffraction indicated that CZX distributes into the liquid crystal phase bilayers but not practically the gel phase POPE bilayers. In the presence of Chol, dialysis and X-ray structural analyses showed that Chol inhibited CZX entry into the bilayer with an increase in Chol concentration. The Chol concentration in the ER membrane (5-10 mol%) is much lower than that in the plasma membrane (approximately 30 mol%). This fact may allow CYP-substrate drugs to enter the hydrophobic portion of the ER membrane more easily than other organelle membranes, yielding efficient drug metabolism.

Exploring the Interactome of Cytochrome P450 2E1 in Human Liver Microsomes with Chemical Crosslinking Mass Spectrometry

Aiming to elucidate the system-wide effects of the alcohol-induced increase in the content of cytochrome P450 2E1 (CYP2E1) on drug metabolism, we explored the array of its protein-protein interactions (interactome) in human liver microsomes (HLM) with chemical crosslinking mass spectrometry (CXMS). Our strategy employs membrane incorporation of purified CYP2E1 modified with photoreactive crosslinkers benzophenone-4-maleimide and 4-(N-succinimidylcarboxy)benzophenone. Exposure of bait-incorporated HLM samples to light was followed by isolating the His-tagged bait protein and its crosslinked aggregates on Ni-NTA agarose. Analyzing the individual bands of SDS-PAGE slabs of thereby isolated protein with the toolset of untargeted proteomics, we detected the crosslinked dimeric and trimeric complexes of CYP2E1 with other drug-metabolizing enzymes. Among the most extensively crosslinked partners of CYP2E1 are the cytochromes P450 2A6, 2C8, 3A4, 4A11, and 4F2, UDP-glucuronosyltransferases (UGTs) 1A and 2B, fatty aldehyde dehydrogenase (ALDH3A2), epoxide hydrolase 1 (EPHX1), disulfide oxidase 1α (ERO1L), and ribophorin II (RPN2). These results demonstrate the exploratory power of the proposed CXMS strategy and corroborate the concept of tight functional integration in the human drug-metabolizing ensemble through protein-protein interactions of the constituting enzymes.

Role of P34S, G169R, R296C, and S486T Substitutions in Ligand Access and Catalysis for Cytochrome P450 2D6 Allelic Variants CYP2D6*14A and CYP2D6*14B

BACKGROUND

Genetic polymorphism of cytochrome P450 (CYP) contributes to variability in drug metabolism, clearance, and response. This study aimed to investigate the functional and molecular basis for altered ligand binding and catalysis in CYP2D6*14A and CYP2D6*14B, two unique alleles common in the Asian population.

METHODS

CYP proteins expressed in Escherichia coli were studied using the substrate 3-cyano-7- ethoxycoumarin (CEC) and inhibitor probes (quinidine, fluoxetine, paroxetine, terbinafine) in the enzyme assay. Computer modelling was additionally used to create three-dimensional structures of the CYP2D6*14 variants.

RESULTS

Kinetics data indicated significantly reduced intrinsic clearance in CYP2D6*14 variants, suggesting that P34S, G169R, R296C, and S486T substitutions worked cooperatively to alter the conformation of the active site that negatively impacted the deethylase activity of CYP2D6. For the inhibition studies, IC50 values decreased in quinidine, paroxetine, and terbinafine but increased in fluoxetine, suggesting a varied ligand-specific susceptibility to inhibition. Molecular docking further demonstrated the role of P34S and R296C in altering access channel dimensions, thereby affecting ligand access and binding and subsequently resulting in varied inhibition potencies.

CONCLUSION

In summary, the differential selectivity of CYP2D6*14 variants for the ligands (substrate and inhibitor) was governed by the alteration of the active site and access channel architecture induced by the natural mutations found in the alleles.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Differential Interactions of Selected Phytocannabinoids with Human CYP2D6 Polymorphisms

Cytochrome P450 2D6 (CYP2D6) is primarily expressed in the liver and in the central nervous system. It is known to be highly polymorphic in nature. It metabolizes several endogenous substrates such as anandamide (AEA). Concomitantly, it is involved in phase 1 metabolism of several antidepressants, antipsychotics, and other drugs. Research in the field of phytocannabinoids (pCBs) has recently accelerated owing to their legalization and increasing medicinal use for pain and inflammation. The primary component of cannabis is THC, which is well-known for its psychotropic effects. Since CYP2D6 is an important brain and liver P450 and is known to be inhibited by CBD, we investigated the interactions of four important highly prevalent CYP2D6 polymorphisms with selected phytocannabinoids (CBD, THC, CBDV, THCV, CBN, CBG, CBC, β-carophyllene) that are rapidly gaining popularity. We show that there is differential binding of CYP2D6*17 to pCBs as compared to WT CYP2D6. We also perform a more detailed comparison of WT and *17 CYP2D6, which reveals the possible regulation of AEA metabolism by CBD. Furthermore, we use molecular dynamics to delineate the mechanism of this binding, inhibition, and regulation. Taken together, we have found that the interactions of CYP2D6 with pCBs vary by polymorphism and by specific pCB class.

Scalable production and application of Pichia pastoris whole cell catalysts expressing human cytochrome P450 2C9

BACKGROUND

Currently, the numerous and versatile applications in pharmaceutical and chemical industry make the recombinant production of cytochrome P450 enzymes (CYPs) of great biotechnological interest. Accelerating the drug development process by simple, quick and scalable access of human drug metabolites is key for efficient and targeted drug development in response to new and sometimes unexpected medical challenges and needs. However, due its biochemical complexity, scalable human CYP (hCYP) production and their application in preparative biotransformations was still in its infancy.

RESULTS

A scalable bioprocess for fine-tuned co-expression of hCYP2C9 and its essential complementary human cytochrome P450 reductase (hCPR) in the yeast Pichia pastoris (Komagataella phaffii) is presented. High-throughput screening (HTS) of a transformant library employing a set of diverse bidirectional expression systems with different regulation patterns and a fluorimetric assay was used in order to fine-tune hCYP2C9 and hCPR co-expression, and to identify best expressing clonal variants. The bioprocess development for scalable and reliable whole cell biocatalyst production in bioreactors was carried out based on rational optimization criteria. Among the different alternatives studied, a glycerol carbon-limiting strategy at high µ showed highest production rates, while methanol co-addition together with a decrease of µ provided the best results in terms of product to biomass yield and whole cell activity. By implementing the mentioned strategies, up to threefold increases in terms of production rates and/or yield could be achieved in comparison with initial tests. Finally, the performance of the whole cell catalysts was demonstrated successfully in biotransformation using ibuprofen as substrate, demonstrating the expected high selectivity of the human enzyme catalyst for 3'hydroxyibuprofen.

CONCLUSIONS

For the first time a scalable bioprocess for the production of hCYP2C9 whole cell catalysts was successfully designed and implemented in bioreactor cultures, and as well, further tested in a preparative-scale biotransformation of interest. The catalyst engineering procedure demonstrated the efficiency of the employment of a set of differently regulated bidirectional promoters to identify transformants with most effective membrane-bound hCYP/hCPR co-expression ratios and implies to become a model case for the generation of other P. pastoris based catalysts relying on co-expressed enzymes such as other P450 catalysts or enzymes relying on co-expressed enzymes for co-factor regeneration.

An electron transfer competent structural ensemble of membrane-bound cytochrome P450 1A1 and cytochrome P450 oxidoreductase

Cytochrome P450 (CYP) heme monooxygenases require two electrons for their catalytic cycle. For mammalian microsomal CYPs, key enzymes for xenobiotic metabolism and steroidogenesis and important drug targets and biocatalysts, the electrons are transferred by NADPH-cytochrome P450 oxidoreductase (CPR). No structure of a mammalian CYP-CPR complex has been solved experimentally, hindering understanding of the determinants of electron transfer (ET), which is often rate-limiting for CYP reactions. Here, we investigated the interactions between membrane-bound CYP 1A1, an antitumor drug target, and CPR by a multiresolution computational approach. We find that upon binding to CPR, the CYP 1A1 catalytic domain becomes less embedded in the membrane and reorients, indicating that CPR may affect ligand passage to the CYP active site. Despite the constraints imposed by membrane binding, we identify several arrangements of CPR around CYP 1A1 that are compatible with ET. In the complexes, the interactions of the CPR FMN domain with the proximal side of CYP 1A1 are supplemented by more transient interactions of the CPR NADP domain with the distal side of CYP 1A1. Computed ET rates and pathways agree well with available experimental data and suggest why the CYP-CPR ET rates are low compared to those of soluble bacterial CYPs.

Crystal Structures of Drug-Metabolizing CYPs

The complex enzyme kinetics displayed by drug-metabolizing cytochrome P450 enzymes (CYPs) (see Chapter 9 ) can, in part, be explained by an examination of their crystallographic protein structures. Fortunately, despite low sequence similarity between different families of drug-metabolizing CYPs, there exists a high degree of structural homology within the superfamily. This similarity in the protein fold allows for a direct comparison of the structural features of CYPs that contribute toward differences in substrate binding, heterotropic and homotropic cooperativity, and genetic variability in drug metabolism. In this chapter, we first provide an overview of the nomenclature and the role of structural features that are common in all CYPs. We then apply these definitions to understand the different substrate specificities and functions in the CYP3A, CYP2C, and CYP2D families of enzymes.

Preparing Membrane Proteins for Simulation Using CHARMM-GUI

Molecular dynamics simulations of membrane proteins have grown dramatically in the last 20 years. Running these simulations first requires embedding the protein's three-dimensional structure in a lipid bilayer of a suitable composition, one that resembles its native environment. This step is far from trivial, especially for modeling heterogeneous mixtures of lipids. CHARMM-GUI, a webserver for simulation system preparation greatly simplifies this step, allowing for the construction of complex heterogeneous and/or asymmetric membranes. Here, we demonstrate how to use CHARMM-GUI to build the membrane for the outer-membrane protein BamA.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Advanced Understanding of the Electron Transfer Pathway of Cytochrome P450s

Cytochrome P450s are heme-thiolate enzymes that participate in carbon source assimilation, natural compound biosynthesis and xenobiotic metabolism in all kingdoms of life. P450s can catalyze various reactions by using a wide range of organic compounds, thus exhibiting great potential in biotechnological applications. The catalytic reactions of P450s are driven by electron equivalents that are sourced from pyridine nucleotides and delivered by cognate or matching redox partners (RPs). The electron transfer (ET) route from RPs to P450s involves one or more redox center-containing domains. As the rate of ET is one of the main determinants of P450 efficacy, an in-depth understanding of the P450 ET pathway should increase our knowledge of these important enzymes and benefit their further applications. Here, the various P450 RP systems along with current understanding of their ET routes will be reviewed. Notably, state-of-the-art structural studies of the two main types of self-sufficient P450 will also be summarized.

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.

Dissecting the Structural Plasticity and Dynamics of Cytochrome P450 2B4 by Molecular Dynamics Simulations

The plasticity of cytochromes P450 (P450s) is known to contribute significantly to their catalytic capacity of metabolizing various substrates. Although numerous studies have been performed, factors governing the plasticity and dynamics of P450s are still not fully understood. In this study, taking CYP2B4 as an example, we dissect the protein plasticity and dynamics in different environments. CYP2B4 is featured by a high degree of plasticity, which exhibits open, closed, and intermediate states. By analyzing the CYP2B4 crystal structures, we identified the structural features for the closed, open, and intermediate states. Interestingly, formation of the dimer structure was found in the open and intermediate states. The subsequent molecular dynamics (MD) simulations of the open structure in water confirmed the importance of the dimer form in stabilizing the open conformations. MD simulations of the closed and open structures in the membrane environment and the free energies for opening the F-G cassette obtained from the umbrella sampling calculations indicate that the membrane environment is important for stabilizing the F-G cassette. The dynamical network analysis indicates that Asp105 on the B-C loop plays an important role in transiting the structure from the open to the intermediate state. Our results thus unveil the mechanisms of dimer formation and open-to-intermediate transition for CYP2B4 in the water and membrane environments.

Rapid and convenient biotransformation procedure for human drug metabolizing enzymes using permeabilized fission yeast cells

The development of convenient assays for the in vitro study of drug metabolizing enzymes (DMEs) such as cytochromes P450 (CYPs) and UDP-glucuronosyltransferases (UGTs) greatly facilitates metabolism studies of candidate drug compounds and other xenobiotics. We have developed and optimized an experimental approach that combines the advantages of recombinant expression in yeast with a microsomal-like biotransformation and thus allows for rapid and convenient enzymatic assays. Recombinant strains of the fission yeast Schizosaccharomyces pombe have previously been demonstrated to functionally express human CYPs and UGTs. Permeabilization of such cells with Triton X-100 results in the formation of enzyme bags, which can be used as biocatalysts. This protocol describes the preparation of such enzyme bags (3 h) and their application in enzyme activity assays (4 h) utilizing either pro-luminescent substrates and luminescence measurements or non-luminescent substrates and liquid chromatography coupled to mass spectrometry (LC-MS). Both applications provide practical tools for investigating CYP and UGT reactions in vitro without the need for additional sophisticated instrumentation or expertise.

Human CYP2E1-dependent mutagenicity of benzene and its hydroxylated metabolites in V79-derived cells: Suppression and enhancement by ethanol pretreatment

Benzene is a human carcinogen that requires metabolic activation. We previously observed that benzene and its hydroxylated metabolites induce micronuclei in mammalian cells expressing human CYP2E1. This study was initially aimed to study another endpoint, the induction of gene mutations by those compounds in the same cell models. A V79-derived cell line expressing human CYP2E1 and sulfotransferase (SULT) 1A1 (V79-hCYP2E1-hSULT1A1) pretreated with ethanol (a CYP2E1 stabilizer) was used in the hprt gene mutagenicity assay. Phenol, hydroquinone, catechol, and 1,2,4-trihydroxybenzene all induced gene mutations, while they were inactive, or only weakly positive (hydroquinone), in parental V79-Mz cells. Unexpectedly, benzene was non-mutagenic in both cell lines, but it became positive in V79-hCYP2E1-hSULT1A1 cells using regimes of short exposure/long recovery without ethanol pretreatment, for both gene mutations and micronuclei formation. In silico molecular simulation showed binding energies and positions favorable for each compound to be oxidized by human CYP2E1, benzene demonstrating the highest affinity. By tunnel analysis, ethanol binding did not limit benzene to pass tunnel S, which was specifically active for benzene. However, its end product, acetic acid, decreased the occurrence of tunnel S from 5.4 to 2.2% and extended the length of its bottleneck from 5.5 to 9.0 Å. With residual ethanol molecules still being present in CYP2E1 for a period of time after benzene exposure, the acetic acid formed could limit the entrance of benzene, thus inhibit its metabolic activation. In summary, ethanol may interfere with the activation of benzene to mutagenic metabolites, at least in cultured cells.

The Effect of Force-Field Parameters on Cytochrome P450-Membrane Interactions: Structure and Dynamics

The simulation of membrane proteins requires compatible protein and lipid force fields that reproduce the properties of both the protein and the lipid bilayer. Cytochrome P450 enzymes are bitopic membrane proteins with a transmembrane helical anchor and a large cytosolic globular domain that dips into the membrane. As such, they are representative and challenging examples of membrane proteins for simulations, displaying features of both peripheral and integral membrane proteins. We performed molecular dynamics simulations of three cytochrome P450 isoforms (2C9, 2C19 and 1A1) in a 2-oleoyl-1-palmitoyl-sn-glycerol-3-phosphocholine bilayer using two AMBER force field combinations: GAFF-LIPID with ff99SB for the protein, and LIPID14 with ff14SB for the protein. Comparison of the structural and dynamic properties of the proteins, the lipids and the protein-membrane interactions shows differing sensitivity of the cytochrome P450 isoforms to the choice of force field, with generally better agreement with experiment for the LIPID14 + ff14SB combination.

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

Communicated by Ramaswamy H. Sarma.

Differing Membrane Interactions of Two Highly Similar Drug-Metabolizing Cytochrome P450 Isoforms: CYP 2C9 and CYP 2C19

The human cytochrome P450 (CYP) 2C9 and 2C19 enzymes are two highly similar isoforms with key roles in drug metabolism. They are anchored to the endoplasmic reticulum membrane by their N-terminal transmembrane helix and interactions of their cytoplasmic globular domain with the membrane. However, their crystal structures were determined after N-terminal truncation and mutating residues in the globular domain that contact the membrane. Therefore, the CYP-membrane interactions are not structurally well-characterized and their dynamics and the influence of membrane interactions on CYP function are not well understood. We describe herein the modeling and simulation of CYP 2C9 and CYP 2C19 in a phospholipid bilayer. The simulations revealed that, despite high sequence conservation, the small sequence and structural differences between the two isoforms altered the interactions and orientations of the CYPs in the membrane bilayer. We identified residues (including K72, P73, and I99 in CYP 2C9 and E72, R73, and H99 in CYP 2C19) at the protein-membrane interface that contribute not only to the differing orientations adopted by the two isoforms in the membrane, but also to their differing substrate specificities by affecting the substrate access tunnels. Our findings provide a mechanistic interpretation of experimentally observed effects of mutagenesis on substrate selectivity.

Crystal structure of bacterial CYP116B5 heme domain: New insights on class VII P450s structural flexibility and peroxygenase activity

Class VII cytochromes P450 are self-sufficient enzymes carrying a phthalate family oxygenase-like reductase domain and a P450 domain fused in a single polypeptide chain. The biocatalytic applications of CYP116B members are limited by the need of the NADPH cofactor and the lack of crystal structures as a starting point for protein engineering. Nevertheless, we demonstrated that the heme domain of CYP116B5 can use hydrogen peroxide as electron donor bypassing the need of NADPH. Here, we report the crystal structure of CYP116B5 heme domain in complex with histidine at 2.6 Å of resolution. The structure reveals the typical P450 fold and a closed conformation with an active site cavity of 284 ų in volume, accommodating a histidine molecule forming a hydrogen bond with the water molecule present as 6th heme iron ligand. MD simulations in the absence of any ligand revealed the opening of a tunnel connecting the active site to the protein surface through the movement of F-, G- and H-helices. A structural alignment with bacterial cytochromes P450 allowed the identification of amino acids in the proximal heme site potentially involved in peroxygenase activity. The availability of the crystal structure provides the bases for the structure-guided design of new biocatalysts.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Can multiscale simulations unravel the function of metallo-enzymes to improve knowledge-based drug discovery?

Metallo-enzymes are a large class of biomolecules promoting specialized chemical reactions. Quantum-classical quantum mechanics/molecular mechanics molecular dynamics, describing the metal site at quantum mechanics level, while accounting for the rest of system at molecular mechanics level, has an accessible time-scale limited by its computational cost. Hence, it must be integrated with classical molecular dynamics and enhanced sampling simulations to disentangle the functions of metallo-enzymes. In this review, we provide an overview of these computational methods and their capabilities. In particular, we will focus on some systems such as CYP19A1 a Fe-dependent enzyme involved in estrogen biosynthesis, and on Mg2+-dependent DNA/RNA processing enzymes/ribozymes and the spliceosome, a protein-directed ribozyme. This information may guide the discovery of drug-like molecules and genetic manipulation tools.

Four Major Channels Detected in the Cytochrome P450 3A4: A Step toward Understanding Its Multispecificity

We computed the network of channels of the 3A4 isoform of the cytochrome P450 (CYP) on the basis of 16 crystal structures extracted from the Protein Data Bank (PDB). The calculations were performed with version 2 of the CCCPP software that we developed for this research project. We identified the minimal cost paths (MCPs) output by CCCPP as probable ways to access to the buried active site. The algorithm of calculation of the MCPs is presented in this paper, with its original method of visualization of the channels. We found that these MCPs constitute four major channels in CYP3A4. Among the many channels proposed by Cojocaru et al. in 2007, we found that only four of them open in 3A4. We provide a refined description of these channels together with associated quantitative data.

Influence of Transmembrane Helix Mutations on Cytochrome P450-Membrane Interactions and Function

Human cytochrome P450 (CYP) enzymes play an important role in the metabolism of drugs, steroids, fatty acids, and xenobiotics. Microsomal CYPs are anchored in the endoplasmic reticulum membrane by an N-terminal transmembrane (TM) helix that is connected to the globular catalytic domain by a flexible linker sequence. However, the structural and functional importance of the TM-helix is unclear because it has been shown that CYPs can still associate with the membrane and have enzymatic activity in reconstituted systems after truncation or modification of the N-terminal sequence. Here, we investigated the effect of mutations in the N-terminal TM-helix residues of two human steroidogenic enzymes, CYP 17A1 and CYP 19A1, that are major drug targets for cancer therapy. These mutations were originally introduced to increase the expression of the proteins in Escherichia coli. To investigate the effect of the mutations on protein-membrane interactions and function, we carried out coarse-grained and all-atom molecular dynamics simulations of the CYPs in a phospholipid bilayer. We confirmed the orientations of the globular domain in the membrane observed in the simulations by linear dichroism measurements in a Nanodisc. Whereas the behavior of CYP 19A1 was rather insensitive to truncation of the TM-helix, mutations in the TM-helix of CYP 17A1, especially W2A and E3L, led to a gradual drifting of the TM-helix out of the hydrophobic core of the membrane. This instability of the TM-helix could affect interactions with the allosteric redox partner, cytochrome b5, required for CYP 17A1's lyase activity. Furthermore, the simulations showed that the mutant TM-helix influenced the membrane interactions of the CYP 17A1 globular domain. In some simulations, the mutated TM-helix obstructed the substrate access tunnel from the membrane to the CYP active site, indicating a possible effect on enzyme function.

The Catalytic Mechanism of Steroidogenic Cytochromes P450 from All-Atom Simulations: Entwinement with Membrane Environment, Redox Partners, and Post-Transcriptional Regulation

Cytochromes P450 (CYP450s) promote the biosynthesis of steroid hormones with major impact on the onset of diseases such as breast and prostate cancers. By merging distinct functions into the same catalytic scaffold, steroidogenic CYP450s enhance complex chemical transformations with extreme efficiency and selectivity. Mammalian CYP450s and their redox partners are membrane-anchored proteins, dynamically associating to form functional machineries. Mounting evidence signifies that environmental factors are strictly intertwined with CYP450s catalysis. Atomic-level simulations have the potential to provide insights into the catalytic mechanism of steroidogenic CYP450s and on its regulation by environmental factors, furnishing information often inaccessible to experimental means. In this review, after an introduction of computational methods commonly employed to tackle these systems, we report the current knowledge on three steroidogenic CYP450s—CYP11A1, CYP17A1, and CYP19A1—endowed with multiple catalytic functions and critically involved in cancer onset. In particular, besides discussing their catalytic mechanisms, we highlight how the membrane environment contributes to (i) regulate ligand channeling through these enzymes, (ii) modulate their interactions with specific protein partners, (iii) mediate post-transcriptional regulation induced by phosphorylation. The results presented set the basis for developing novel therapeutic strategies aimed at fighting diseases originating from steroid metabolism dysfunction.

Monotopic Membrane Proteins Join the Fold

Monotopic membrane proteins, classified by topology, are proteins that embed into a single face of the membrane. These proteins are generally underrepresented in the Protein Data Bank (PDB), but the past decade of research has revealed new examples that allow the description of generalizable features. This Opinion article summarizes shared characteristics including oligomerization states, modes of membrane association, mechanisms of interaction with hydrophobic or amphiphilic substrates, and homology to soluble folds. We also discuss how associations of monotopic enzymes in pathways can be used to promote substrate specificity and product composition. These examples highlight the challenges in structure determination specific to this class of proteins, but also the promise of new understanding from future study of these proteins that reside at the interface.

Molecular Dynamics Simulations Reveal Structural Differences among Allelic Variants of Membrane-Anchored Cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) is an enzyme that is involved in the metabolism of roughly 25% of all marketed drugs and therefore belongs to the most important enzymes in drug metabolism. CYP2D6 features a high degree of genetic polymorphism that can significantly affect the metabolic activity of an individual. In extreme cases, structural changes at the level of single amino acids can either increase its enzymatic activity abolishing the drug therapeutic effect or completely disable the enzyme and elevate drug plasma level potentially leading to adverse effects. In this study, starting from the crystal structure, we built a full-length membrane-anchored all-atom model of the wild-type CYP2D6 as well as five of its variants differing in the enzymatic activity. We validated our models with available experimental data and compared their structural properties with molecular dynamics simulations. The main focus of this study was to identify differences that could mechanistically explain the altered activity of the variants and improve our understanding of their functioning. We observed differences in the opening frequencies and minimal diameters of tunnels that connect the buried active site to the surrounding solvent environment. The variants CYP2D6*4 and CYP2D6*10 associated with missing or decreased activity showed less frequent opening of the tunnels compared to the wild-type. Both CYP2D6*10 and CYP2D6*17 showed a deprivation of an important ligand tunnel suggesting a feasible reason for their altered substrate specificity. Next, the altered fold at the N-terminal anchor region and the decreased active site volume caused by the amino acid mutations of the CYP2D6*4 variant offer an explanation for the absence of its metabolic activity. The mutations in CYP2D6*53 contributed to a significant enlargement of an important ligand tunnel and an extension of the active site cavity. This could explain the altered metabolic profile as well as the enhanced metabolic rates of this particular variant supporting its designation as a possible cause for the ultrarapid metabolizer phenotype. We believe these novel structural insights could advance the fields of personalized medicine and enzyme engineering. Furthermore, they could aid in guiding laboratory as well as computational experiments in the future.

Ligand Access Channels in Cytochrome P450 Enzymes: A Review

Quantitative structure-activity relationships may bring invaluable information on structural elements of both enzymes and substrates that, together, govern substrate specificity. Buried active sites in cytochrome P450 enzymes are connected to the solvent by a network of channels exiting at the distal surface of the protein. This review presents different in silico tools that were developed to uncover such channels in P450 crystal structures. It also lists some of the experimental evidence that actually suggest that these predicted channels might indeed play a critical role in modulating P450 functions. Amino acid residues at the entrance of the channels may participate to a first global ligand recognition of ligands by P450 enzymes before they reach the buried active site. Moreover, different P450 enzymes show different networks of predicted channels. The plasticity of P450 structures is also important to take into account when looking at how channels might play their role.

Insights into molecular mechanisms of drug metabolism dysfunction of human CYP2C9*30

Cytochrome P450 2C9 (CYP2C9) metabolizes about 15% of clinically administrated drugs. The allelic variant CYP2C9*30 (A477T) is associated to diminished response to the antihypertensive effects of the prodrug losartan and affected metabolism of other drugs. Here, we investigated molecular mechanisms involved in the functional consequences of this amino-acid substitution. Molecular dynamics (MD) simulations performed for the active species of the enzyme (heme in the Compound I state), in the apo or substrate-bound state, and binding energy analyses gave insights into altered protein structure and dynamics involved in the defective drug metabolism of human CYP2C9.30. Our data revealed an increased rigidity of the key Substrate Recognition Sites SRS1 and SRS5 and shifting of the β turn 4 of SRS6 toward the helix F in CYP2C9.30. Channel and binding substrate dynamics analyses showed altered substrate channel access and active site accommodation. These conformational and dynamic changes are believed to be involved in the governing mechanism of the reduced catalytic activity. An ensemble of representative conformations of the WT and A477T mutant properly accommodating drug substrates were identified, those structures can be used for prediction of new CYP2C9 and CYP2C9.30 substrates and drug-drug interactions.

Lipid molecules can induce an opening of membrane-facing tunnels in cytochrome P450 1A2

Cytochrome P450 1A2 (P450 1A2, CYP1A2) is a membrane-bound enzyme that oxidizes a broad range of hydrophobic substrates. The structure and dynamics of both the catalytic and trans-membrane (TM) domains of this enzyme in the membrane/water environment were investigated using a multiscale computational approach, including coarse-grained and all-atom molecular dynamics. Starting from the spontaneous self-assembly of the system containing the TM or soluble domain immersed in randomized dilauroyl phosphatidylcholine (DLPC)/water mixture into their respective membrane-bound forms, we reconstituted the membrane-bound structure of the full-length P450 1A2. This structure includes a TM helix that spans the membrane, while being connected to the catalytic domain by a short flexible loop. Furthermore, in this model, the upper part of the TM helix interacts directly with a conserved and highly hydrophobic N-terminal proline-rich segment of the catalytic domain; this segment and the FG loop are immersed in the membrane, whereas the remaining portion of the catalytic domain remains exposed to aqueous solution. The shallow membrane immersion of the catalytic domain induces a depression in the opposite intact layer of the phospholipids. This structural effect may help in stabilizing the position of the TM helix directly beneath the catalytic domain. The partial immersion of the catalytic domain also allows for the enzyme substrates to enter the active site from either aqueous solution or phospholipid environment via several solvent- and membrane-facing tunnels in the full-length P450 1A2. The calculated tunnel dynamics indicated that the opening probability of the membrane-facing tunnels is significantly enhanced when a DLPC molecule spontaneously penetrates into the membrane-facing tunnel 2d. The energetics of the lipid penetration process were assessed by the linear interaction energy (LIE) approximation, and found to be thermodynamically feasible.

Toward a systems approach to the human cytochrome P450 ensemble: interactions between CYP2D6 and CYP2E1 and their functional consequences

Functional cross-talk among human drug-metabolizing cytochrome P450 through their association is a topic of emerging importance. Here, we studied the interactions of human CYP2D6, a major metabolizer of psychoactive drugs, with one of the most prevalent human P450 enzymes, ethanol-inducible CYP2E1. Detection of P450-P450 interactions was accomplished through luminescence resonance energy transfer between labeled proteins incorporated into human liver microsomes and the microsomes of insect cells containing NADPH-cytochrome P450 reductase. The potential of CYP2D6 to form oligomers in the microsomal membrane is among the highest observed with human cytochrome P450 studied up to date. We also observed the formation of heteromeric complexes of CYP2D6 with CYP2E1 and CYP3A4, and found a significant modulation of these interactions by 3,4-methylenedioxymethylamphetamine, a widespread drug of abuse metabolized by CYP2D6. Our results demonstrate an ample alteration of the catalytic properties of CYP2D6 and CYP2E1 caused by their association. In particular, we demonstrated that preincubation of microsomes containing co-incorporated CYP2D6 and CYP2E1 with CYP2D6-specific substrates resulted in considerable time-dependent activation of CYP2D6, which presumably occurs via a slow substrate-induced reorganization of CYP2E1-CYP2D6 hetero-oligomers. Furthermore, we demonstrated that the formation of heteromeric complexes between CYP2E1 and CYP2D6 affects the stoichiometry of futile cycling and substrate oxidation by CYP2D6 by means of decreasing the electron leakage through the peroxide-generating pathways. Our results further emphasize the role of P450-P450 interactions in regulatory cross-talk in human drug-metabolizing ensemble and suggest a role of interactions of CYP2E1 with CYP2D6 in pharmacologically important instances of alcohol-drug interactions.

Thermodynamics of camphor migration in cytochrome P450cam by atomistic simulations

Understanding the mechanisms of ligand binding to enzymes is of paramount importance for the design of new drugs. Here, we report on the use of a novel biased molecular dynamics (MD) methodology to study the mechanism of camphor binding to cytochrome P450cam. Microsecond-long MD simulations allowed us to observe reaction coordinates characterizing ligand diffusion from the active site of cytochrome P450cam to solvent via three egress routes. These atomistic simulations were used to estimate thermodynamic quantities along the reaction coordinates and indicate diverse binding configurations. The results suggest that the diffusion of camphor along the pathway near the substrate recognition site (SRS) is thermodynamically preferred. In addition, we show that the diffusion near the SRS is triggered by a transition from a heterogeneous collection of closed ligand-bound conformers to the basin comprising the open conformations of cytochrome P450cam. The conformational change accompanying this switch is characterized by the retraction of the F and G helices and the disorder of the B' helix. These results are corroborated by experimental studies and provide detailed insight into ligand binding and conformational behavior of the cytochrome family. The presented methodology is general and can be applied to other ligand-protein systems.

The catalytic function of cytochrome P450 is entwined with its membrane-bound nature

Cytochrome P450, a family of monooxygenase enzymes, is organized as a catalytic metabolon, which requires enzymatic partners as well as environmental factors that tune its complex dynamic. P450 and its reducing counterparts—cytochrome P450-reductase and cytochrome b ₅—are membrane-bound proteins located in the cytosolic side of the endoplasmic reticulum. They are believed to dynamically associate to form functional complexes. Increasing experimental evidence signifies the role(s) played by both protein-protein and protein-lipid interactions in P450 catalytic function and efficiency. However, the biophysical challenges posed by their membrane-bound nature have severely limited high-resolution understanding of the molecular interfaces of these interactions. In this article, we provide an overview of the current knowledge on cytochrome P450, highlighting the environmental factors that are entwined with its metabolic function. Recent advances in structural biophysics are also discussed, setting up the bases for a new paradigm in the study of this important class of membrane-bound enzymes.

Clinical, genetic, and structural basis of congenital adrenal hyperplasia due to 11β-hydroxylase deficiency

Congenital adrenal hyperplasia (CAH), resulting from mutations in CYP11B1, a gene encoding 11β-hydroxylase, represents a rare autosomal recessive Mendelian disorder of aberrant sex steroid production. Unlike CAH caused by 21-hydroxylase deficiency, the disease is far more common in the Middle East and North Africa, where consanguinity is common often resulting in identical mutations. Clinically, affected female newborns are profoundly virilized (Prader score of 4/5), and both genders display significantly advanced bone ages and are oftentimes hypertensive. We find that 11-deoxycortisol, not frequently measured, is the most robust biochemical marker for diagnosing 11β-hydroxylase deficiency. Finally, computational modeling of 25 missense mutations of CYP11B1 revealed that specific modifications in the heme-binding (R374W and R448C) or substrate-binding (W116C) site of 11β-hydroxylase, or alterations in its stability (L299P and G267S), may predict severe disease. Thus, we report clinical, genetic, hormonal, and structural effects of CYP11B1 gene mutations in the largest international cohort of 108 patients with steroid 11β-hydroxylase deficiency CAH.

Physical Studies of P450-P450 Interactions: Predicting Quaternary Structures of P450 Complexes in Membranes from Their X-ray Crystal Structures

Cytochrome P450 enzymes, which catalyze oxygenation reactions of both exogenous and endogenous chemicals, are membrane bound proteins that require interaction with their redox partners in order to function. Those responsible for drug and foreign compound metabolism are localized primarily in the endoplasmic reticulum of liver, lung, intestine, and other tissues. More recently, the potential for P450 enzymes to exist as supramolecular complexes has been shown by the demonstration of both homomeric and heteromeric complexes. The P450 units in these complexes are heterogeneous with respect to their distribution and function, and the interaction of different P450s can influence P450-specific metabolism. The goal of this review is to examine the evidence supporting the existence of physical complexes among P450 enzymes. Additionally, the review examines the crystal lattices of different P450 enzymes derived from X-ray diffraction data to make assumptions regarding possible quaternary structures in membranes and in turn, to predict how the quaternary structures could influence metabolism and explain the functional effects of specific P450-P450 interactions.