1
|
Mandal N, Surpeta B, Brezovsky J. Reinforcing Tunnel Network Exploration in Proteins Using Gaussian Accelerated Molecular Dynamics. J Chem Inf Model 2024; 64:6623-6635. [PMID: 39143923 DOI: 10.1021/acs.jcim.4c00966] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/16/2024]
Abstract
Tunnels are structural conduits in biomolecules responsible for transporting chemical compounds and solvent molecules from the active site. They have been shown to be present in a wide variety of enzymes across all functional and structural classes. However, the study of such pathways is experimentally challenging, because they are typically transient. Computational methods, such as molecular dynamics (MD) simulations, have been successfully proposed to explore tunnels. Conventional MD (cMD) provides structural details to characterize tunnels but suffers from sampling limitations to capture rare tunnel openings on longer time scales. Therefore, in this study, we explored the potential of Gaussian accelerated MD (GaMD) simulations to improve the exploration of complex tunnel networks in enzymes. We used the haloalkane dehalogenase LinB and its two variants with engineered transport pathways, which are not only well-known for their application potential but have also been extensively studied experimentally and computationally regarding their tunnel networks and their importance in multistep catalytic reactions. Our study demonstrates that GaMD efficiently improves tunnel sampling and allows the identification of all known tunnels for LinB and its two mutants. Furthermore, the improved sampling provided insight into a previously unknown transient side tunnel (ST). The extensive conformational landscape explored by GaMD simulations allowed us to investigate in detail the mechanism of ST opening. We determined variant-specific dynamic properties of ST opening, which were previously inaccessible due to limited sampling of cMD. Our comprehensive analysis supports multiple indicators of the functional relevance of the ST, emphasizing its potential significance beyond structural considerations. In conclusion, our research proves that the GaMD method can overcome the sampling limitations of cMD for the effective study of tunnels in enzymes, providing further means for identifying rare tunnels in enzymes with the potential for drug development, precision medicine, and rational protein engineering.
Collapse
Affiliation(s)
- Nishita Mandal
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, Poznan 61-614, Poland
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, Warsaw 02-109, Poland
| | - Bartlomiej Surpeta
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, Poznan 61-614, Poland
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, Warsaw 02-109, Poland
| | - Jan Brezovsky
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, Poznan 61-614, Poland
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, Warsaw 02-109, Poland
| |
Collapse
|
2
|
Feng Y, Gong C, Zhu J, Liu G, Tang Y, Li W. Unraveling the Ligand-Binding Sites of CYP3A4 by Molecular Dynamics Simulations with Solvent Probes. J Chem Inf Model 2024; 64:3451-3464. [PMID: 38593186 DOI: 10.1021/acs.jcim.4c00089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
Cytochrome P450 3A4 (CYP3A4) is one of the most important drug-metabolizing enzymes in the human body and is well known for its complicated, atypical kinetic characteristics. The existence of multiple ligand-binding sites in CYP3A4 has been widely recognized as being capable of interfering with the active pocket through allosteric effects. The identification of ligand-binding sites other than the canonical active site above the heme is especially important for understanding the atypical kinetic characteristics of CYP3A4 and the intriguing association between the ligand and the receptor. In this study, we first employed mixed-solvent molecular dynamics (MixMD) simulations coupled with the online computational predictive tools to explore potential ligand-binding sites in CYP3A4. The MixMD approach demonstrates better performance in dealing with the receptor flexibility compared with other computational tools. From the sites identified by MixMD, we then picked out multiple sites for further exploration using ensemble docking and conventional molecular dynamics (cMD) simulations. Our results indicate that three extra sites are suitable for ligand binding in CYP3A4, including one experimentally confirmed site and two novel sites.
Collapse
Affiliation(s)
- Yanjun Feng
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Changda Gong
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Jieyu Zhu
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Guixia Liu
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Yun Tang
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Weihua Li
- Shanghai Frontiers Science Center of Optogenetic Techniques for Cell Metabolism, Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| |
Collapse
|
3
|
Ackad EM, Biggers L, Meister M, Kontoyianni M. Equilibrium landscape of ingress/egress channels and gating residues of the Cytochrome P450 3A4. PLoS One 2024; 19:e0298424. [PMID: 38498575 PMCID: PMC10947690 DOI: 10.1371/journal.pone.0298424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 01/24/2024] [Indexed: 03/20/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Edward Michael Ackad
- Department of Physics, Southern Illinois University Edwardsville, Edwardsville, IL, United States of America
| | - Laurence Biggers
- Department of Internal Medicine at UT Southwestern Medical Center, Dallas, TX, United States of America
| | - Mary Meister
- University of Illinois College of Medicine, Chicago, IL, United States of America
| | - Maria Kontoyianni
- Department of Pharmaceutical Sciences, Southern Illinois University Edwardsville, Edwardsville, IL, United States of America
| |
Collapse
|
4
|
Kaluzhskiy L, Yablokov E, Gnedenko O, Burkatovskii D, Maslov I, Bogorodskiy A, Ershov P, Tsybruk T, Zelepuga E, Rutckova T, Kozlovskaya E, Dmitrenok P, Gilep A, Borshchevskiy V, Strushkevich N, Ivanov A. The effect of membrane composition on the interaction between human CYP51 and its flavonoid inhibitor - luteolin 7,3'-disulfate. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2024; 1866:184286. [PMID: 38272204 DOI: 10.1016/j.bbamem.2024.184286] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Revised: 12/22/2023] [Accepted: 01/16/2024] [Indexed: 01/27/2024]
Abstract
Cytochromes P450 (CYP) are a family of membrane proteins involved in the production of endogenous molecules and the metabolism of xenobiotics. It is well-known that the composition of the membrane can influence the activity and orientation of CYP proteins. However, little is known about how membrane composition affects the ligand binding properties of CYP. In this study, we utilized surface plasmon resonance and fluorescence lifetime analysis to examine the impact of membrane micro-environment composition on the interaction between human microsomal CYP51 (CYP51A1) and its inhibitor, luteolin 7,3'-disulphate (LDS). We observed that membranes containing cholesterol or sphingomyelin exhibited the lowest apparent equilibrium dissociation constant for the CYP51A1-LDS complex. Additionally, the tendency for relation between kinetic parameters of the CYP51A1-LDS complex and membrane viscosity and overall charge was observed. These findings suggest that the specific composition of the membrane, particularly the presence of cholesterol and sphingomyelin, plays a vital role in regulating the interaction between CYP enzymes and their ligands.
Collapse
Affiliation(s)
- Leonid Kaluzhskiy
- Institute of Biomedical Chemistry, 10 Building 8, Pogodinskaya Street, 119121 Moscow, Russia.
| | - Evgeniy Yablokov
- Institute of Biomedical Chemistry, 10 Building 8, Pogodinskaya Street, 119121 Moscow, Russia.
| | - Oksana Gnedenko
- Institute of Biomedical Chemistry, 10 Building 8, Pogodinskaya Street, 119121 Moscow, Russia.
| | - Dmitrii Burkatovskii
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia.
| | - Ivan Maslov
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia.
| | - Andrey Bogorodskiy
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia
| | - Pavel Ershov
- Institute of Biomedical Chemistry, 10 Building 8, Pogodinskaya Street, 119121 Moscow, Russia
| | - Tatsiana Tsybruk
- Institute of Bioorganic Chemistry NASB, 5 Building 2, V.F. Kuprevich Street, 220141 Minsk, Belarus
| | - Elena Zelepuga
- G.B. Elyakov Pacific Institute of Bioorganic Chemistry, Far Eastern Branch of the Russian Academy of Science, 159 Prospect 100-letiya Vladivostoka, 690022 Vladivostok, Russia.
| | - Tatyana Rutckova
- G.B. Elyakov Pacific Institute of Bioorganic Chemistry, Far Eastern Branch of the Russian Academy of Science, 159 Prospect 100-letiya Vladivostoka, 690022 Vladivostok, Russia
| | - Emma Kozlovskaya
- G.B. Elyakov Pacific Institute of Bioorganic Chemistry, Far Eastern Branch of the Russian Academy of Science, 159 Prospect 100-letiya Vladivostoka, 690022 Vladivostok, Russia
| | - Pavel Dmitrenok
- G.B. Elyakov Pacific Institute of Bioorganic Chemistry, Far Eastern Branch of the Russian Academy of Science, 159 Prospect 100-letiya Vladivostoka, 690022 Vladivostok, Russia.
| | - Andrei Gilep
- Institute of Biomedical Chemistry, 10 Building 8, Pogodinskaya Street, 119121 Moscow, Russia; Institute of Bioorganic Chemistry NASB, 5 Building 2, V.F. Kuprevich Street, 220141 Minsk, Belarus.
| | - Valentin Borshchevskiy
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia; Joint Institute for Nuclear Research, Dubna 141980, Russia
| | - Natallia Strushkevich
- Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, bld. 1, 121205 Moscow, Russia.
| | - Alexis Ivanov
- Institute of Biomedical Chemistry, 10 Building 8, Pogodinskaya Street, 119121 Moscow, Russia.
| |
Collapse
|
5
|
Paço L, Hackett JC, Atkins WM. Nanodisc-embedded cytochrome P450 P3A4 binds diverse ligands by distributing conformational dynamics to its flexible elements. J Inorg Biochem 2023; 244:112211. [PMID: 37080138 PMCID: PMC10175226 DOI: 10.1016/j.jinorgbio.2023.112211] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Revised: 03/12/2023] [Accepted: 04/03/2023] [Indexed: 04/22/2023]
Abstract
Cytochrome P450 3A4 (CYP3A4) metabolizes a wide range of drugs and toxins. Interactions of CYP3A4 with ligands are difficult to predict due to promiscuity and conformational flexibility. To better understand CYP3A4 conformational responses to ligands we use hydrogen deuterium exchange mass spectrometry (HDX-MS) to investigate the effect of ligands on nanodisc-embedded CYP3A4. For a subset of CYP3A4-ligand complexes, differences in the low-frequency modes derived by principal component analyses of molecular dynamics trajectories mirrored the HDX-MS results. The effects of ligands are distributed to flexible elements of CYP3A4 between stretches of secondary structure. The largest effects occur in the F- and G-helices, where most ligands increase the flexibility of the F-helix and connecting loops and decrease the flexibility of the C-term of the G-helix. Most ligands affect the E-F-G, CD and HI regions of the protein. Ligand-dependent differences are observed in the A"-A' loop, BC region, E-helix, K-β1 region, proximal loop, and C-term loop. Correlated HDX responses were observed in the CD region and the C-term of the G-helix that were most pronounced for Type II ligands. Collectively, the HDX and molecular dynamics results suggest that CYP3A4 accommodates diverse binding partners by propagating local backbone fluctuations from the binding site onto the flexible regions of the enzyme via long-range interactions that are differentially modulated by ligands. In contrast to the paradigm wherein ligands decrease protein dynamics at their binding site, a wide range of ligands modestly increase CYP3A4 dynamics throughout the protein including effects remote from the active site.
Collapse
Affiliation(s)
- Lorela Paço
- Department of Medicinal Chemistry, University of Washington, Seattle, WA 98195-7610, United States of America
| | - John C Hackett
- Department of Chemistry and Biochemistry, Florida International University, Miami, FL 33199, United States of America
| | - William M Atkins
- Department of Medicinal Chemistry, University of Washington, Seattle, WA 98195-7610, United States of America.
| |
Collapse
|
6
|
Fischer A, Bardakci F, Sellner M, Lill MA, Smieško M. Ligand pathways in estrogen-related receptors. J Biomol Struct Dyn 2023; 41:1639-1648. [PMID: 35068382 DOI: 10.1080/07391102.2022.2027818] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The three subtypes of estrogen-related receptors ERRα, ERRβ, and ERRγ are nuclear receptors mediating metabolic processes in various tissues such as the skeletal muscle, fat tissue, bone, and liver. Although the knowledge on their physiological ligands is limited, they have been implicated as drug targets for important indications including diabetes, cardiovascular diseases, and osteoporosis. As in other nuclear receptors, their ligand binding pocket is buried within the core of the receptor and connected to its surrounding by ligand pathways. Here, we investigated these pathways with conventional molecular dynamics as well as metadynamics simulations to reveal their distribution and their capability to facilitate ligand translocation. Dependent on the ERR subtype and the conformational state of the receptor, we could detect different pathways to be favored. Overall, the results suggested pathways IIIa and IIIb to be favored in the agonistic conformation, while antagonists preferred pathways I, II, and V. Along the pathways, the ligands passed different gating mechanisms of the receptor, including groups of protein residues as well as whole secondary structure elements, to leave the binding site. Even though these pathways are suggested to influence ligand specificity of the receptors and their elucidation might advance rational drug design, they have not yet been studied in ERRs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- André Fischer
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, Basel, Switzerland
| | - Ferhat Bardakci
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, Basel, Switzerland
| | - Manuel Sellner
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, Basel, Switzerland
| | - Markus A Lill
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, Basel, Switzerland
| | - Martin Smieško
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, Basel, Switzerland
| |
Collapse
|
7
|
Zhou Y, Li J, Baryshnikov G, Tu Y. Unraveling the Abnormal Molecular Mechanism of Suicide Inhibition of Cytochrome P450 3A4. J Chem Inf Model 2022; 62:6172-6181. [PMID: 36457253 PMCID: PMC9749025 DOI: 10.1021/acs.jcim.2c01035] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
Suicide inhibition of the CYP3A4 enzyme by a drug inactivates the enzyme in the drug biotransformation process and often shows safety concerns about the drug. Despite extensive experimental studies, the abnormal molecular mechanism of a suicide inhibitor that forms a covalent bond with the residue far away from the catalytically active center of CYP3A4 inactivating the enzyme remains elusive. Here, the authors used molecular simulation approaches to study in detail how diquinone methide (DQR), the metabolite product of raloxifene, unbinds from CYP3A4 and inactivates the enzyme at the atomistic level. The results clearly indicate that in one of the intermediate states formed in its unbinding process, DQR covalently binds to Cys239, a residue far away from the catalytically active center of CYP3A4, and hinders the substrate from entering or leaving the enzyme. This work therefore provides an unprecedented way of clarifying the abnormal mechanism of suicide inhibition of the CYP3A4 enzyme.
Collapse
Affiliation(s)
- Yang Zhou
- School of Pharmacy, Jinan University, 601 Huangpu Avenue West, Guangzhou510632, China.,Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, 114 28Stockholm, Sweden
| | - Junhao Li
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, 114 28Stockholm, Sweden
| | - Glib Baryshnikov
- Laboratory of Organic Electronics, Department of Science and Technology, Linköping University, 60174Norrköping, Sweden
| | - Yaoquan Tu
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, 114 28Stockholm, Sweden
| |
Collapse
|
8
|
Fu T, Zhang H, Zheng Q. Molecular Insights into the Heterotropic Allosteric Mechanism in Cytochrome P450 3A4-Mediated Midazolam Metabolism. J Chem Inf Model 2022; 62:5762-5770. [PMID: 36342224 DOI: 10.1021/acs.jcim.2c01264] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Cytochrome P450 3A4 (CYP3A4) is the main P450 enzyme for drug metabolism and drug-drug interactions (DDIs), as it is involved in the metabolic process of approximately 50% of drugs. A detailed mechanistic elucidation of DDIs mediated by CYP3A4 is commonly believed to be critical for drug optimization and rational use. Here, two typical probes, midazolam (MDZ, substrate) and testosterone (TST, allosteric effector), are used to investigate the molecular mechanism of CYP3A4-mediated heterotropic allosteric interactions, through conventional molecular dynamics (cMD) and well-tempered metadynamics (WT-MTD) simulations. Distance monitoring shows that TST can stably bind in two potential peripheral sites (Site 1 and Site 2) of CYP3A4. The binding of TST at these two sites can induce conformational changes in CYP3A4 flexible loops on the basis of conformational analysis, thereby promoting the transition of the MDZ binding mode and affecting the ratio of MDZ metabolites. According to the results of the residue interaction network, multiple allosteric communication pathways are identified that can provide vivid and applicable insights into the heterotropic allostery of TST on MDZ metabolism. Comparing the regulatory effects and the communication pathways, the allosteric effect caused by TST binding in Site 2 seems to be more pronounced than in Site 1. Our findings could provide a deeper understanding of CYP3A4-mediated heterotropic allostery at the atomic level and would be helpful for rational drug use as well as the design of new allosteric modulators.
Collapse
Affiliation(s)
- Tingting Fu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Hongxing Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Qingchuan Zheng
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China.,Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, Jilin University, Changchun 130023, China
| |
Collapse
|
9
|
Srdič M, Fessner ND, Yildiz D, Glieder A, Spiertz M, Schwaneberg U. Preparative Production of Functionalized (N- and O-Heterocyclic) Polycyclic Aromatic Hydrocarbons by Human Cytochrome P450 3A4 in a Bioreactor. Biomolecules 2022; 12:biom12020153. [PMID: 35204652 PMCID: PMC8961652 DOI: 10.3390/biom12020153] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 01/12/2022] [Accepted: 01/13/2022] [Indexed: 11/16/2022] Open
Abstract
Polycyclic aromatic hydrocarbons (PAHs) and their N- and O-containing derivatives (N-/O-PAHs) are environmental pollutants and synthetically attractive building blocks in pharmaceuticals. Functionalization of PAHs can be achieved via C-H activation by cytochrome P450 enzymes (e.g., P450 CYP3A4) in an environmentally friendly manner. Despite its broad substrate scope, the contribution of CYP3A4 to metabolize common PAHs in humans was found to be small. We recently showcased the potential of CYP3A4 in whole-cell biocatalysis with recombinant yeast Komagataella phaffii (Pichia pastoris) catalysts for the preparative-scale synthesis of naturally occurring metabolites in humans. In this study, we aimed at exploring the substrate scope of CYP3A4 towards (N-/O)-PAHs and conducted a bioconversion experiment at 10 L scale to validate the synthetic potential of CYP3A4 for the preparative-scale production of functionalized PAH metabolites. Hydroxylated products were purified and characterized using HPLC and NMR analysis. In total, 237 mg of fluorenol and 48 mg of fluorenone were produced from 498 mg of fluorene, with peak productivities of 27.7 μmol/L/h for fluorenol and 5.9 μmol/L/h for fluorenone; the latter confirmed that CYP3A4 is an excellent whole-cell biocatalyst for producing authentic human metabolites.
Collapse
Affiliation(s)
- Matic Srdič
- SeSaM-Biotech GmbH, 52074 Aachen, Germany;
- Institute of Biotechnology, RWTH Aachen University, 52074 Aachen, Germany
| | - Nico D. Fessner
- Institute of Molecular Biotechnology, Graz University of Technology, NAWI Graz, 8010 Graz, Austria;
| | - Deniz Yildiz
- DWI—Leibniz Institute for Interactive Materials, 52074 Aachen, Germany;
- Institute for Technical and Macromolecular Chemistry, RWTH Aachen University, 52074 Aachen, Germany
| | | | - Markus Spiertz
- SeSaM-Biotech GmbH, 52074 Aachen, Germany;
- Correspondence: (M.S.); (U.S.)
| | - Ulrich Schwaneberg
- Institute of Biotechnology, RWTH Aachen University, 52074 Aachen, Germany
- DWI—Leibniz Institute for Interactive Materials, 52074 Aachen, Germany;
- Correspondence: (M.S.); (U.S.)
| |
Collapse
|
10
|
Fischer A, Smieško M. A Conserved Allosteric Site on Drug-Metabolizing CYPs: A Systematic Computational Assessment. Int J Mol Sci 2021; 22:13215. [PMID: 34948012 PMCID: PMC8707821 DOI: 10.3390/ijms222413215] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Revised: 11/29/2021] [Accepted: 12/01/2021] [Indexed: 11/17/2022] Open
Abstract
Cytochrome P450 enzymes (CYPs) are the largest group of enzymes involved in human drug metabolism. Ligand tunnels connect their active site buried at the core of the membrane-anchored protein to the surrounding solvent environment. Recently, evidence of a superficial allosteric site, here denoted as hotspot 1 (H1), involved in the regulation of ligand access in a soluble prokaryotic CYP emerged. Here, we applied multi-scale computational modeling techniques to study the conservation and functionality of this allosteric site in the nine most relevant mammalian CYPs responsible for approximately 70% of drug metabolism. In total, we systematically analyzed over 44 μs of trajectories from conventional MD, cosolvent MD, and metadynamics simulations. Our bioinformatic analysis and simulations with organic probe molecules revealed the site to be well conserved in the CYP2 family with the exception of CYP2E1. In the presence of a ligand bound to the H1 site, we could observe an enlargement of a ligand tunnel in several members of the CYP2 family. Further, we could detect the facilitation of ligand translocation by H1 interactions with statistical significance in CYP2C8 and CYP2D6, even though all other enzymes except for CYP2C19, CYP2E1, and CYP3A4 presented a similar trend. As the detailed comprehension of ligand access and egress phenomena remains one of the most relevant challenges in the field, this work contributes to its elucidation and ultimately helps in estimating the selectivity of metabolic transformations using computational techniques.
Collapse
Affiliation(s)
| | - Martin Smieško
- Computational Pharmacy, Departement of Pharmaceutical Sciences, University of Basel, Klingelbergstrasse 61, 4056 Basel, Switzerland;
| |
Collapse
|
11
|
Fessner ND, Grimm C, Srdič M, Weber H, Kroutil W, Schwaneberg U, Glieder A. Natural Product Diversification by One‐Step Biocatalysis using Human P450 3A4. ChemCatChem 2021. [DOI: 10.1002/cctc.202101564] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Nico D. Fessner
- Institute of Molecular Biotechnology NAWI Graz Graz University of Technology Petersgasse 14 8010 Graz Austria
| | - Christopher Grimm
- Institute of Chemistry NAWI Graz University of Graz Heinrichstraße 28 8010 Graz Austria
| | - Matic Srdič
- SeSaM-Biotech GmbH Forckenbeckstraße 50 52074 Aachen Germany
- Bisy GmbH Wuenschendorf 292 Hofstätten an der Raab 8200 Hofstaetten Austria
| | - Hansjörg Weber
- Institute of Organic Chemistry NAWI Graz Graz University of Technology Stremayrgasse 9 8010 Graz Austria
| | - Wolfgang Kroutil
- Institute of Chemistry NAWI Graz University of Graz Heinrichstraße 28 8010 Graz Austria
| | - Ulrich Schwaneberg
- Institute of Biotechnology RWTH Aachen University Worringerweg 3 52074 Aachen Germany
| | - Anton Glieder
- Institute of Molecular Biotechnology NAWI Graz Graz University of Technology Petersgasse 14 8010 Graz Austria
| |
Collapse
|
12
|
Róg T, Girych M, Bunker A. Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design. Pharmaceuticals (Basel) 2021; 14:1062. [PMID: 34681286 PMCID: PMC8537670 DOI: 10.3390/ph14101062] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 10/14/2021] [Accepted: 10/15/2021] [Indexed: 11/17/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Tomasz Róg
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Mykhailo Girych
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Alex Bunker
- Drug Research Program, Division of Pharmaceutical Biosciences, Faculty of Pharmacy, University of Helsinki, 00014 Helsinki, Finland;
| |
Collapse
|
13
|
Molecular probes for human cytochrome P450 enzymes: Recent progress and future perspectives. Coord Chem Rev 2021. [DOI: 10.1016/j.ccr.2020.213600] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
14
|
Fessner ND, Srdič M, Weber H, Schmid C, Schönauer D, Schwaneberg U, Glieder A. Preparative‐Scale Production of Testosterone Metabolites by Human Liver Cytochrome P450 Enzyme 3A4. Adv Synth Catal 2020. [DOI: 10.1002/adsc.202000251] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Nico D. Fessner
- Institute of Molecular BiotechnologyGraz University of Technology, NAWI Graz Petersgasse 14/3 Austria
| | - Matic Srdič
- SeSaM-Biotech GmbH Aachen Germany
- Bisy GmbH Hofstaetten Austria
| | - Hansjörg Weber
- Institute of Organic ChemistryGraz University of Technology, NAWI Graz Austria
| | - Christian Schmid
- Institute of Molecular BiotechnologyGraz University of Technology, NAWI Graz Petersgasse 14/3 Austria
- Austrian Centre of Industrial Biotechnology (ACIB) Graz Austria
| | | | | | - Anton Glieder
- Institute of Molecular BiotechnologyGraz University of Technology, NAWI Graz Petersgasse 14/3 Austria
| |
Collapse
|
15
|
Ritacco I, Saltalamacchia A, Spinello A, Ippoliti E, Magistrato A. All-Atom Simulations Disclose How Cytochrome Reductase Reshapes the Substrate Access/Egress Routes of Its Partner CYP450s. J Phys Chem Lett 2020; 11:1189-1193. [PMID: 31986051 DOI: 10.1021/acs.jpclett.9b03798] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Cytochromes P450 enzymes (CYP450s) promote the oxidative metabolism of a variety of substrates via the electrons supplied by the cytochrome P450 reductase (CPR) and upon formation of a CPR/CYP450 adduct. In spite of the pivotal regulatory importance of this process, the impact of CPR binding on the functional properties of its partner CYP450 remains elusive. By performing multiple microsecond-long all-atom molecular dynamics simulations of a 520 000-atom model of a CPR/CYP450 adduct embedded in a membrane mimic, we disclose the molecular terms for their interactions, considering the aromatase (HA) enzyme as a proxy of the CYP450 family. Our study strikingly unveils that CPR binding alters HA's functional motions, bolstering a change in the shape and type of the channels traveled by substrates/products during their access/egress to/from the enzyme's active site. Our outcomes unprecedentedly contribute to extricate the many entangled facets of the CYP450 metabolon, redrafting its intricate panorama from an atomic-level perspective.
Collapse
Affiliation(s)
- Ida Ritacco
- CNR-IOM c/o SISSA , via Bonomea 265 , 34136 Trieste , Italy
| | | | | | - Emiliano Ippoliti
- IAS-5/INM-9 Computational Biomedicine Institute and JARA-HPC , Forschungszentrum Jülich , Wilhelm-Johnen-Straße , 52425 Jülich , Germany
| | | |
Collapse
|
16
|
Fischer A, Smieško M. Spontaneous Ligand Access Events to Membrane-Bound Cytochrome P450 2D6 Sampled at Atomic Resolution. Sci Rep 2019; 9:16411. [PMID: 31712722 PMCID: PMC6848145 DOI: 10.1038/s41598-019-52681-w] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Accepted: 10/18/2019] [Indexed: 12/12/2022] Open
Abstract
The membrane-anchored enzyme Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of around 25% of marketed drugs and its metabolic performance shows a high interindividual variation. While it was suggested that ligands access the buried active site of the enzyme from the membrane, no proof from unbiased simulations has been provided to support this hypothesis. Laboratory experiments fail to capture the access process which is suspected to influence binding kinetics. Here, we applied unbiased molecular dynamics (MD) simulations to investigate the access of ligands to wild-type CYP2D6, as well as the allelic variant CYP2D6*53. In multiple simulations, substrates accessed the active site of the enzyme from the protein-membrane interface to ultimately adopt a conformation that would allow a metabolic reaction. We propose the necessary steps for ligand access and the results suggest that the increased metabolic activity of CYP2D6*53 might be caused by a facilitated ligand uptake.
Collapse
Affiliation(s)
- André Fischer
- University of Basel, Department of Pharmaceutical Sciences, Basel, 4056, Switzerland
| | - Martin Smieško
- University of Basel, Department of Pharmaceutical Sciences, Basel, 4056, Switzerland.
| |
Collapse
|
17
|
Escalante DE, Aksan A. Role of Water Hydrogen Bonding on Transport of Small Molecules inside Hydrophobic Channels. J Phys Chem B 2019; 123:6673-6685. [PMID: 31310534 DOI: 10.1021/acs.jpcb.9b03060] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We conducted a systematic analysis of water networking inside smooth hyperboloid hydrophobic structures (cylindrical, barrel, and hourglass shapes) to elucidate the role of water hydrogen bonding on the transport of small hydrophobic molecules (ligands). Through a series of molecular dynamics simulations, we established that a hydrogen-bonded network forming along the centerline resulted in a water exclusion zone adjacent to the walls. The size of the exclusion zone is a function of the geometry and the nonbonded interaction strength, defining the effective hydrophobicity of the structure. Exclusion of water molecules from this zone results in lower apparent viscosity, leading to acceleration of ligand transport up to 7 times faster than that measured in the bulk. Transport of ligands into and out of the hydrophobic structures was shown to be controlled by a single water molecule that capped the narrow regions in the structure. This mechanism provides physical insights into the behavior and role of water in the bottleneck regions of hydrophobic enzyme channels. These findings were then used in a sister publication [ Escalante , D. E. , Comput. Struct. Biotechnol. J. 2019 17 757 760 ] to develop a model that can accurately predict the transport of ligands along nanochannels of broad-substrate specificity enzymes.
Collapse
Affiliation(s)
- Diego E Escalante
- Department of Mechanical Engineering , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Alptekin Aksan
- Department of Mechanical Engineering , University of Minnesota , Minneapolis , Minnesota 55455 , United States.,BioTechnology Institute , University of Minnesota , St. Paul , Minnesota 55108 , United States
| |
Collapse
|
18
|
Pravda L, Sehnal D, Toušek D, Navrátilová V, Bazgier V, Berka K, Svobodová Vareková R, Koca J, Otyepka M. MOLEonline: a web-based tool for analyzing channels, tunnels and pores (2018 update). Nucleic Acids Res 2019; 46:W368-W373. [PMID: 29718451 PMCID: PMC6030847 DOI: 10.1093/nar/gky309] [Citation(s) in RCA: 196] [Impact Index Per Article: 39.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2018] [Accepted: 04/12/2018] [Indexed: 12/27/2022] Open
Abstract
MOLEonline is an interactive, web-based application for the detection and characterization of channels (pores and tunnels) within biomacromolecular structures. The updated version of MOLEonline overcomes limitations of the previous version by incorporating the recently developed LiteMol Viewer visualization engine and providing a simple, fully interactive user experience. The application enables two modes of calculation: one is dedicated to the analysis of channels while the other was specifically designed for transmembrane pores. As the application can use both PDB and mmCIF formats, it can be leveraged to analyze a wide spectrum of biomacromolecular structures, e.g. stemming from NMR, X-ray and cryo-EM techniques. The tool is interconnected with other bioinformatics tools (e.g., PDBe, CSA, ChannelsDB, OPM, UniProt) to help both setup and the analysis of acquired results. MOLEonline provides unprecedented analytics for the detection and structural characterization of channels, as well as information about their numerous physicochemical features. Here we present the application of MOLEonline for structural analyses of α-hemolysin and transient receptor potential mucolipin 1 (TRMP1) pores. The MOLEonline application is freely available via the Internet at https://mole.upol.cz.
Collapse
Affiliation(s)
- Lukáš Pravda
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic
| | - David Sehnal
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic
| | - Dominik Toušek
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17. Listopadu 12, 771 46 Olomouc, Czech Republic
| | - Veronika Navrátilová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tr. 17. Listopadu 12, 771 46 Olomouc, Czech Republic
| | - Václav Bazgier
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tr. 17. Listopadu 12, 771 46 Olomouc, Czech Republic
| | - Karel Berka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tr. 17. Listopadu 12, 771 46 Olomouc, Czech Republic
| | - Radka Svobodová Vareková
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic
| | - Jaroslav Koca
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tr. 17. Listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
19
|
Escalante DE, Aksan A. Prediction of Ligand Transport along Hydrophobic Enzyme Nanochannels. Comput Struct Biotechnol J 2019; 17:757-760. [PMID: 31303980 PMCID: PMC6606821 DOI: 10.1016/j.csbj.2019.06.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2019] [Revised: 05/08/2019] [Accepted: 06/01/2019] [Indexed: 12/22/2022] Open
Abstract
Buried active sites of enzymes are connected to the bulk solvent through a network of hydrophobic channels. We developed a discretized model that can accurately predict ligand transport along hydrophobic channels up to six orders of magnitude faster than any other existing method. The non-dimensional nature of the model makes it applicable to any hydrophobic channel/ligand combination.
Collapse
Affiliation(s)
- Diego E. Escalante
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN 55455, United States
| | - Alptekin Aksan
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN 55455, United States
- BioTechnology Institute, University of Minnesota, St. Paul, MN 55108, United States
- Corresponding author at: Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN 55455, United States.
| |
Collapse
|
20
|
Hoias Teixeira M, Menegon Arantes G. Balanced internal hydration discriminates substrate binding to respiratory complex I. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2019; 1860:541-548. [DOI: 10.1016/j.bbabio.2019.05.004] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/06/2018] [Revised: 05/16/2019] [Accepted: 05/28/2019] [Indexed: 12/16/2022]
|
21
|
Paul S, Nair NN, Vashisth H. Phase space and collective variable based simulation methods for studies of rare events. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2019.1634268] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, USA
| | - Nisanth N. Nair
- Department of Chemistry, Indian Institute of Technology, Kanpur, India
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, USA
| |
Collapse
|
22
|
Four Major Channels Detected in the Cytochrome P450 3A4: A Step toward Understanding Its Multispecificity. Int J Mol Sci 2019; 20:ijms20040987. [PMID: 30823507 PMCID: PMC6412807 DOI: 10.3390/ijms20040987] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2019] [Revised: 02/12/2019] [Accepted: 02/20/2019] [Indexed: 12/27/2022] Open
Abstract
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.
Collapse
|
23
|
Nunes-Alves A, Zuckerman DM, Arantes GM. Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways. Biophys J 2019. [PMID: 29539393 DOI: 10.1016/j.bpj.2018.01.014] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.
Collapse
Affiliation(s)
- Ariane Nunes-Alves
- Department of Biochemistry, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
| | - Daniel M Zuckerman
- Department of Biomedical Engineering, School of Medicine, Oregon Health & Science University, Portland, Oregon.
| | | |
Collapse
|
24
|
The Catalytic Mechanism of Steroidogenic Cytochromes P450 from All-Atom Simulations: Entwinement with Membrane Environment, Redox Partners, and Post-Transcriptional Regulation. Catalysts 2019. [DOI: 10.3390/catal9010081] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
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.
Collapse
|
25
|
Kreß N, Halder JM, Rapp LR, Hauer B. Unlocked potential of dynamic elements in protein structures: channels and loops. Curr Opin Chem Biol 2018; 47:109-116. [DOI: 10.1016/j.cbpa.2018.09.010] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2018] [Revised: 08/22/2018] [Accepted: 09/11/2018] [Indexed: 10/28/2022]
|
26
|
Awasthi S, Nair NN. Exploring high‐dimensional free energy landscapes of chemical reactions. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1398] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Shalini Awasthi
- Department of Chemistry Indian Institute of Technology Kanpur Uttar Pradesh India
| | - Nisanth N. Nair
- Department of Chemistry Indian Institute of Technology Kanpur Uttar Pradesh India
| |
Collapse
|
27
|
Šrejber M, Navrátilová V, Paloncýová M, Bazgier V, Berka K, Anzenbacher P, Otyepka M. Membrane-attached mammalian cytochromes P450: An overview of the membrane's effects on structure, drug binding, and interactions with redox partners. J Inorg Biochem 2018; 183:117-136. [DOI: 10.1016/j.jinorgbio.2018.03.002] [Citation(s) in RCA: 81] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Revised: 02/16/2018] [Accepted: 03/01/2018] [Indexed: 01/08/2023]
|
28
|
Ligand Access Channels in Cytochrome P450 Enzymes: A Review. Int J Mol Sci 2018; 19:ijms19061617. [PMID: 29848998 PMCID: PMC6032366 DOI: 10.3390/ijms19061617] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 05/27/2018] [Accepted: 05/28/2018] [Indexed: 12/22/2022] Open
Abstract
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.
Collapse
|
29
|
Rydzewski J, Jakubowski R, Nowak W, Grubmüller H. Kinetics of Huperzine A Dissociation from Acetylcholinesterase via Multiple Unbinding Pathways. J Chem Theory Comput 2018; 14:2843-2851. [PMID: 29715428 DOI: 10.1021/acs.jctc.8b00173] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
The dissociation of huperzine A (hupA) from Torpedo californica acetylcholinesterase ( TcAChE) was investigated by 4 μs unbiased and biased all-atom molecular dynamics (MD) simulations in explicit solvent. We performed our study using memetic sampling (MS) for the determination of reaction pathways (RPs), metadynamics to calculate free energy, and maximum-likelihood estimation (MLE) to recover kinetic rates from unbiased MD simulations. Our simulations suggest that the dissociation of hupA occurs mainly via two RPs: a front door along the axis of the active-site gorge (pwf) and through a new transient side door (pws), i.e., formed by the Ω-loop (residues 67-94 of TcAChE). An analysis of the inhibitor unbinding along the RPs suggests that pws is opened transiently after hupA and the Ω-loop reach a low free-energy transition state characterized by the orientation of the pyridone group of the inhibitor directed toward the Ω-loop plane. Unlike pws, pwf does not require large structural changes in TcAChE to be accessible. The estimated free energies and rates agree well with available experimental data. The dissociation rates along the unbinding pathways are similar, suggesting that the dissociation of hupA along pws is likely to be relevant. This indicates that perturbations to hupA- TcAChE interactions could potentially induce pathway hopping. In summary, our results characterize the slow-onset inhibition of TcAChE by hupA, which may provide the structural and energetic bases for the rational design of the next-generation slow-onset inhibitors with optimized pharmacokinetic properties for the treatment of Alzheimer's disease.
Collapse
Affiliation(s)
- J Rydzewski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University , Grudziadzka 5 , 87-100 Toruń , Poland
| | - R Jakubowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University , Grudziadzka 5 , 87-100 Toruń , Poland.,Centre of New Technologies, University of Warsaw , Banacha 2c , 02-097 Warsaw , Poland
| | - W Nowak
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University , Grudziadzka 5 , 87-100 Toruń , Poland
| | - H Grubmüller
- Max Planck Institute for Biophysical Chemistry , Am Faßberg 11 , 37077 Göttingen , Germany
| |
Collapse
|
30
|
Don CG, Smieško M. Out‐compute drug side effects: Focus on cytochrome P450 2D6 modeling. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1366] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Charleen G. Don
- Department of Pharmaceutical SciencesUniversity of BaselBaselSwitzerland
| | - Martin Smieško
- Department of Pharmaceutical SciencesUniversity of BaselBaselSwitzerland
| |
Collapse
|
31
|
Fan JR, Li H, Zhang HX, Zheng QC. Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine. Biopolymers 2018; 109:e23108. [PMID: 29484634 DOI: 10.1002/bip.23108] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2017] [Revised: 02/06/2018] [Accepted: 02/07/2018] [Indexed: 12/25/2022]
Abstract
The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure-function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.
Collapse
Affiliation(s)
- Jing-Rong Fan
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| | - Heng Li
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, College of Life Science, Jilin University, Changchun, 130012, People's Republic of China
| | - Hong-Xing Zhang
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| | - Qing-Chuan Zheng
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, 130023, People's Republic of China.,Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, College of Life Science, Jilin University, Changchun, 130012, People's Republic of China
| |
Collapse
|
32
|
Pravda L, Sehnal D, Svobodová Vařeková R, Navrátilová V, Toušek D, Berka K, Otyepka M, Koča J. ChannelsDB: database of biomacromolecular tunnels and pores. Nucleic Acids Res 2018; 46:D399-D405. [PMID: 29036719 PMCID: PMC5753359 DOI: 10.1093/nar/gkx868] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2017] [Revised: 09/08/2017] [Accepted: 09/28/2017] [Indexed: 01/21/2023] Open
Abstract
ChannelsDB (http://ncbr.muni.cz/ChannelsDB) is a database providing information about the positions, geometry and physicochemical properties of channels (pores and tunnels) found within biomacromolecular structures deposited in the Protein Data Bank. Channels were deposited from two sources; from literature using manual deposition and from a software tool automatically detecting tunnels leading to the enzymatic active sites and selected cofactors, and transmembrane pores. The database stores information about geometrical features (e.g. length and radius profile along a channel) and physicochemical properties involving polarity, hydrophobicity, hydropathy, charge and mutability. The stored data are interlinked with available UniProt annotation data mapping known mutation effects to channel-lining residues. All structures with channels are displayed in a clear interactive manner, further facilitating data manipulation and interpretation. As such, ChannelsDB provides an invaluable resource for research related to deciphering the biological function of biomacromolecular channels.
Collapse
Affiliation(s)
- Lukáš Pravda
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
- National Centre for Biomolecular Research, Faculty of Science, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
| | - David Sehnal
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
- National Centre for Biomolecular Research, Faculty of Science, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
| | - Radka Svobodová Vařeková
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
- National Centre for Biomolecular Research, Faculty of Science, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
| | - Veronika Navrátilová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Dominik Toušek
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Karel Berka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jaroslav Koča
- CEITEC - Central European Institute of Technology, Masaryk University Brno, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
- National Centre for Biomolecular Research, Faculty of Science, Kamenice 5, 625 00 Brno-Bohunice, Czech Republic
| |
Collapse
|
33
|
Navrátilová V, Paloncýová M, Berka K, Mise S, Haga Y, Matsumura C, Sakaki T, Inui H, Otyepka M. Molecular insights into the role of a distal F240A mutation that alters CYP1A1 activity towards persistent organic pollutants. Biochim Biophys Acta Gen Subj 2017; 1861:2852-2860. [DOI: 10.1016/j.bbagen.2017.08.002] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2017] [Revised: 06/30/2017] [Accepted: 08/01/2017] [Indexed: 01/12/2023]
|
34
|
Thermodynamics of camphor migration in cytochrome P450cam by atomistic simulations. Sci Rep 2017; 7:7736. [PMID: 28798338 PMCID: PMC5552751 DOI: 10.1038/s41598-017-07993-0] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2017] [Accepted: 07/05/2017] [Indexed: 01/12/2023] Open
Abstract
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.
Collapse
|
35
|
Magistrato A. Direct in silico visualization of ligands channelling through proteins: The next-generation frontier of computational biology: Comment on 'Ligand diffusion via enhanced sampling molecular dynamics' by Jakub Rydzewski and Wieslaw Nowak. Phys Life Rev 2017; 22-23:82-84. [PMID: 28818495 DOI: 10.1016/j.plrev.2017.08.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2017] [Accepted: 08/07/2017] [Indexed: 10/19/2022]
Affiliation(s)
- Alessandra Magistrato
- CNR-IOM-Democritos c/o, International School for Advanced Studies (SISSA), Via Bonomea 265, 34135, Trieste, Italy.
| |
Collapse
|
36
|
Benkaidali L, André F, Moroy G, Tangour B, Maurel F, Petitjean M. The Cytochrome P450 3A4 has three Major Conformations: New Clues to Drug Recognition by this Promiscuous Enzyme. Mol Inform 2017; 36. [PMID: 28685969 DOI: 10.1002/minf.201700044] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2017] [Accepted: 06/28/2017] [Indexed: 01/02/2023]
Abstract
We computed the channels of the 3A4 isoform of the cytochrome P450 3A4 (CYP) on the basis of 24 crystal structures extracted from the Protein Data Bank (PDB). We identified three major conformations (denoted C, O1 and O2) using an enhanced version of the CCCPP software that we developed for the present work, while only two conformations (C and O2 ) are considered in the literature. We established the flowchart of definition of these three conformations in function of the structural and physicochemical parameters of the ligand. The channels are characterized with qualitative and quantitative parameters, and not only with their surrounding secondary structures as it is usually done in the literature.
Collapse
Affiliation(s)
- Lydia Benkaidali
- ITODYS, CNRS UMR 7086, Université Paris Diderot, France.,Université de Carthage, Faculté des Sciences, Département de Chimie, Bizerte, Tunisie
| | - François André
- CEA/I2BC, CNRS UMR 9198, Université Paris-Saclay, France
| | - Gautier Moroy
- MTi, INSERM UMR-S 973, Université Paris Diderot, France
| | - Bahoueddine Tangour
- Unité de Recherche de Modélisation en Sciences Fondamentales et Didactique, BP244, Université de Tunis El Manar, 2092, Tunis, Tunisie
| | | | - Michel Petitjean
- MTi, INSERM UMR-S 973, Université Paris Diderot, France.,Epôle de génoinformatique, CNRS UMR 7592, Institut Jacques Monod, Paris, France
| |
Collapse
|
37
|
Magistrato A, Sgrignani J, Krause R, Cavalli A. Single or Multiple Access Channels to the CYP450s Active Site? An Answer from Free Energy Simulations of the Human Aromatase Enzyme. J Phys Chem Lett 2017; 8:2036-2042. [PMID: 28423275 DOI: 10.1021/acs.jpclett.7b00697] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Cytochromes P450 (CYP450s), in particular, CYP19A1 and CYP17A1, are key clinical targets of breast and prostate anticancer therapies, critical players in drug metabolism, and their overexpression in tumors is associated with drug resistance. In these enzymes, ligand (substrates, drugs) metabolism occurs in deeply buried active sites accessible only via several grueling channels, whose exact biological role remains unclear. Gaining direct insights on the mechanism by which ligands travel in and out is becoming increasingly important given that channels are involved in the modulation of binding/dissociation kinetics and the specificity of ligands toward a CYP450. This has profound implications for enzymatic efficiency and drug efficacy/toxicity. Here, by applying free energy methods, for a cumulative simulation time of 20 μs, we provide detailed atomistic characterization and free energy profiles of the entry/exit routes preferentially followed by a substrate (androstenedione) and a last-generation inhibitor (letrozole) to/from the catalytic site of CYP19A1 (the human aromatase (HA) enzyme), a key clinical target against breast cancer, studied here as prototypical CYP450. Despite the remarkably different size/shape/hydrophobicity of the ligands, two channels appear accessible to their entrance, while only one exit route appears to be preferential. Our study shows that the preferential paths may be conserved among different CYP450s. Moreover, our results highlight that, at least in the case of HA, ligand channeling is associated with large enzyme structural rearrangements. A wise choice of the computational method and very long simulations are, thus, required to obtain fully converged quantitative free energy profiles, which might be used to design novel biocatalysts or next-generation cytochrome inhibitors with an in silico tuned Km.
Collapse
Affiliation(s)
| | - Jacopo Sgrignani
- Institute for Research in Biomedicine (IRB), Università della Svizzera Italiana (USI) , Via Vincenzo Vela 6, CH-6500 Bellinzona, Switzerland
| | - Rolf Krause
- Institute of Computational Science, Faculty of Informatics, Università della Svizzera Italiana (USI) , Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Andrea Cavalli
- Institute for Research in Biomedicine (IRB), Università della Svizzera Italiana (USI) , Via Vincenzo Vela 6, CH-6500 Bellinzona, Switzerland
| |
Collapse
|
38
|
Hermans SM, Pfleger C, Nutschel C, Hanke CA, Gohlke H. Rigidity theory for biomolecules: concepts, software, and applications. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1311] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Susanne M.A. Hermans
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| | - Christopher Pfleger
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| | - Christina Nutschel
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| | - Christian A. Hanke
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry; Heinrich Heine University Düsseldorf; Düsseldorf Germany
| |
Collapse
|
39
|
Rydzewski J, Nowak W. Ligand diffusion in proteins via enhanced sampling in molecular dynamics. Phys Life Rev 2017; 22-23:58-74. [PMID: 28410930 DOI: 10.1016/j.plrev.2017.03.003] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2016] [Revised: 10/28/2016] [Accepted: 03/28/2017] [Indexed: 01/17/2023]
Abstract
Computational simulations in biophysics describe the dynamics and functions of biological macromolecules at the atomic level. Among motions particularly important for life are the transport processes in heterogeneous media. The process of ligand diffusion inside proteins is an example of a complex rare event that can be modeled using molecular dynamics simulations. The study of physical interactions between a ligand and its biological target is of paramount importance for the design of novel drugs and enzymes. Unfortunately, the process of ligand diffusion is difficult to study experimentally. The need for identifying the ligand egress pathways and understanding how ligands migrate through protein tunnels has spurred the development of several methodological approaches to this problem. The complex topology of protein channels and the transient nature of the ligand passage pose difficulties in the modeling of the ligand entry/escape pathways by canonical molecular dynamics simulations. In this review, we report a methodology involving a reconstruction of the ligand diffusion reaction coordinates and the free-energy profiles along these reaction coordinates using enhanced sampling of conformational space. We illustrate the above methods on several ligand-protein systems, including cytochromes and G-protein-coupled receptors. The methods are general and may be adopted to other transport processes in living matter.
Collapse
Affiliation(s)
- J Rydzewski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruń, Poland.
| | - W Nowak
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruń, Poland.
| |
Collapse
|
40
|
Escalante DE, Aukema KG, Wackett LP, Aksan A. Simulation of the Bottleneck Controlling Access into a Rieske Active Site: Predicting Substrates of Naphthalene 1,2-Dioxygenase. J Chem Inf Model 2017; 57:550-561. [PMID: 28170277 DOI: 10.1021/acs.jcim.6b00469] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Naphthalene 1,2-dioxygenase (NDO) has been computationally understudied despite the extensive experimental knowledge obtained for this enzyme, including numerous crystal structures and over 100 demonstrated substrates. In this study, we have developed a substrate prediction model that moves away from the traditional active-site-centric approach to include the energetics of substrate entry into the active site. By comparison with experimental data, the accuracy of the model for predicting substrate oxidation is 92%, with a positive predictive value of 93% and a negative predictive value of 98%. Also, the present analysis has revealed that the amino acid residues that provided the largest energetic barrier for compounds entering the active site are residues F224, L227, P234, and L235. In addition, F224 is proposed to play a role in controlling ligand entrance via π-π stacking stabilization as well as providing stabilization via T-shaped π-π interactions once the ligand has reached the active-site cavity. Overall, we present a method capable of being scaled to computationally discover thousands of substrates of NDO, and we present parameters to be used for expanding the prediction method to other members of the Rieske non-heme iron oxygenase family.
Collapse
Affiliation(s)
- Diego E Escalante
- Department of Mechanical Engineering, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - Kelly G Aukema
- Department of Biochemistry, Molecular Biology and Biophysics, University of Minnesota , Minneapolis, Minnesota 55455, United States.,BioTechnology Institute, University of Minnesota , St. Paul, Minnesota 55108, United States
| | - Lawrence P Wackett
- Department of Biochemistry, Molecular Biology and Biophysics, University of Minnesota , Minneapolis, Minnesota 55455, United States.,BioTechnology Institute, University of Minnesota , St. Paul, Minnesota 55108, United States
| | - Alptekin Aksan
- Department of Mechanical Engineering, University of Minnesota , Minneapolis, Minnesota 55455, United States.,BioTechnology Institute, University of Minnesota , St. Paul, Minnesota 55108, United States
| |
Collapse
|
41
|
Riccardi L, Arencibia JM, Bono L, Armirotti A, Girotto S, De Vivo M. Lid domain plasticity and lipid flexibility modulate enzyme specificity in human monoacylglycerol lipase. Biochim Biophys Acta Mol Cell Biol Lipids 2017; 1862:441-451. [PMID: 28088576 DOI: 10.1016/j.bbalip.2017.01.002] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2016] [Revised: 12/20/2016] [Accepted: 01/09/2017] [Indexed: 12/17/2022]
Abstract
Human monoacylglycerol lipase (MAGL) is a membrane-interacting enzyme that generates pro-inflammatory signaling molecules. For this reason, MAGL inhibition is a promising strategy to treat pain, cancer, and neuroinflammatory diseases. MAGL can hydrolyze monoacylglycerols bearing an acyl chain of different lengths and degrees of unsaturation, cleaving primarily the endocannabinoid 2-arachidonoylglycerol. Importantly, the enzymatic binding site of MAGL is confined by a 75-amino-acid-long, flexible cap domain, named 'lid domain', which is structurally similar to that found in several other lipases. However, it is unclear how lid domain plasticity affects catalysis in MAGL. By integrating extensive molecular dynamics simulations and free-energy calculations with mutagenesis and kinetic experiments, we here define a lid-domain-mediated mechanism for substrate selection and binding in MAGL catalysis. In particular, we clarify the key role of Phe159 and Ile179, two conserved residues within the lid domain, in regulating substrate specificity in MAGL. We conclude by proposing that other structurally related lipases may share this lid-domain-mediated mechanism for substrate specificity.
Collapse
Affiliation(s)
- Laura Riccardi
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Jose M Arencibia
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Luca Bono
- D3-PharmaChemistry, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Andrea Armirotti
- D3-PharmaChemistry, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Stefania Girotto
- CompuNet, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy; IAS-5/INM-9 Computational Biomedicine Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52428 Jülich, Germany.
| |
Collapse
|
42
|
Navrátilová V, Paloncýová M, Berka K, Otyepka M. Effect of Lipid Charge on Membrane Immersion of Cytochrome P450 3A4. J Phys Chem B 2016; 120:11205-11213. [DOI: 10.1021/acs.jpcb.6b10108] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Veronika Navrátilová
- Regional
Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacký University Olomouc, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Markéta Paloncýová
- Regional
Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacký University Olomouc, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Karel Berka
- Regional
Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacký University Olomouc, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Regional
Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacký University Olomouc, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
43
|
Di Meo F, Fabre G, Berka K, Ossman T, Chantemargue B, Paloncýová M, Marquet P, Otyepka M, Trouillas P. In silico pharmacology: Drug membrane partitioning and crossing. Pharmacol Res 2016; 111:471-486. [PMID: 27378566 DOI: 10.1016/j.phrs.2016.06.030] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/20/2016] [Revised: 06/30/2016] [Accepted: 06/30/2016] [Indexed: 01/09/2023]
Abstract
Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.
Collapse
Affiliation(s)
- Florent Di Meo
- INSERM UMR 850, Univ. Limoges, Faculty of Pharmacy, 2 rue du Dr Marcland, F-87025, Limoges, France
| | - Gabin Fabre
- LCSN, Univ. Limoges, Faculty of Pharmacy, 2 rue du Dr Marcland, F-87025, Limoges, France
| | - Karel Berka
- Regional Centre for Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky̿ University, Olomouc, Czech Republic
| | - Tahani Ossman
- INSERM UMR 850, Univ. Limoges, Faculty of Pharmacy, 2 rue du Dr Marcland, F-87025, Limoges, France
| | - Benjamin Chantemargue
- INSERM UMR 850, Univ. Limoges, Faculty of Pharmacy, 2 rue du Dr Marcland, F-87025, Limoges, France; Regional Centre for Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky̿ University, Olomouc, Czech Republic
| | - Markéta Paloncýová
- Regional Centre for Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky̿ University, Olomouc, Czech Republic
| | - Pierre Marquet
- INSERM UMR 850, Univ. Limoges, Faculty of Pharmacy, 2 rue du Dr Marcland, F-87025, Limoges, France
| | - Michal Otyepka
- Regional Centre for Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky̿ University, Olomouc, Czech Republic
| | - Patrick Trouillas
- INSERM UMR 850, Univ. Limoges, Faculty of Pharmacy, 2 rue du Dr Marcland, F-87025, Limoges, France; Regional Centre for Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky̿ University, Olomouc, Czech Republic.
| |
Collapse
|