1
|
Zerroug E, Belaidi S, Chtita S, Tuffaha G, AbulQais F, Kciuk M, Dubey A, Taha MO. Structure‐Based Approaches for the Prediction of Alzheimer's Disease Inhibitors: Comparative Interactions Analysis, Pharmacophore Modeling and Molecular Dynamics Simulations. ChemistrySelect 2024; 9. [DOI: 10.1002/slct.202303307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 12/18/2023] [Indexed: 07/10/2024]
Abstract
AbstractDue to its significant role in neurodegeneration, Cyclin‐dependent kinase 5 (CDK5) has emerged as a potential target for addressing neuropathological disorders, including Alzheimer's disease (AD). The application of CDK5 inhibitors has demonstrated promise in the treatment of AD. This prompted us to model this interesting target using a computational workflow named Docking‐based Comparative Intermolecular Contacts Analysis (dbCICA). Approaches including 3D‐QSAR, genetic algorithm, and pharmacophore modeling were employed to discover new CDK inhibitors.
Collapse
Affiliation(s)
- Enfale Zerroug
- Group of Computational and Pharmaceutical Chemistry LMCE Laboratory University of Biskra, BP 145 Biskra 07000 Algeria
| | - Salah Belaidi
- Group of Computational and Pharmaceutical Chemistry LMCE Laboratory University of Biskra, BP 145 Biskra 07000 Algeria
| | - Samir Chtita
- Faculty of Sciences Ben M'Sik Hassan II University of Casablanca Sidi Othman, Casablanca Morocco
| | | | - Faizan AbulQais
- Department of Agricultural Microbiology Aligarh Muslim University Aligarh UP 202002 India
| | - Mateusz Kciuk
- Doctoral School of Exact and Natural Sciences University of Lodz Banacha Street 12/16 90-237 Lodz Poland
- Department of Molecular Biotechnology and Genetics University of Lodz Banacha 12/16 90-237 Lodz Poland
| | - Amit Dubey
- Department of Pharmacology Saveetha Dental College and Hospital Saveetha Institute of Medical and Technical Sciences Chennai Tamil Nadu India
| | - Mutasem O. Taha
- Drug Discovery Unit Department of Pharmaceutical Sciences Faculty of Pharmacy University of Jordan 11942 Amman Jordan
| |
Collapse
|
2
|
Zsidó BZ, Bayarsaikhan B, Börzsei R, Szél V, Mohos V, Hetényi C. The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering. Int J Mol Sci 2023; 24:11784. [PMID: 37511543 PMCID: PMC10381018 DOI: 10.3390/ijms241411784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Revised: 07/18/2023] [Accepted: 07/20/2023] [Indexed: 07/30/2023] Open
Abstract
Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.
Collapse
Affiliation(s)
- Balázs Zoltán Zsidó
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Bayartsetseg Bayarsaikhan
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Rita Börzsei
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Viktor Szél
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Violetta Mohos
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Csaba Hetényi
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| |
Collapse
|
3
|
Ashworth MA, Bombino E, de Jong RM, Wijma HJ, Janssen DB, McLean KJ, Munro AW. Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin. ACS Catal 2022; 12:15028-15044. [PMID: 36570080 PMCID: PMC9764288 DOI: 10.1021/acscatal.2c03974] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2022] [Revised: 10/22/2022] [Indexed: 11/29/2022]
Abstract
CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.
Collapse
Affiliation(s)
- Mark A. Ashworth
- Manchester
Institute of Biotechnology, School of Chemistry, The University of Manchester, Manchester M1 7DN, United Kingdom
| | - Elvira Bombino
- Department
of Biochemistry, Groningen Biomolecular Sciences and Biotechnology
Institute, University of Groningen, Nijenborgh 4, Groningen 9747 AG, Netherlands
| | - René M. de Jong
- DSM
Food & Beverage, Alexander Fleminglaan 1, 2613 AX Delft, the Netherlands
| | - Hein J. Wijma
- Department
of Biochemistry, Groningen Biomolecular Sciences and Biotechnology
Institute, University of Groningen, Nijenborgh 4, Groningen 9747 AG, Netherlands
| | - Dick B. Janssen
- Department
of Biochemistry, Groningen Biomolecular Sciences and Biotechnology
Institute, University of Groningen, Nijenborgh 4, Groningen 9747 AG, Netherlands,
| | - Kirsty J. McLean
- Manchester
Institute of Biotechnology, School of Chemistry, The University of Manchester, Manchester M1 7DN, United Kingdom,Department
of Biological and Geographical Sciences, School of Applied Sciences, University of Huddersfield, Huddersfield HD1 3DH, United Kingdom
| | - Andrew W. Munro
- Manchester
Institute of Biotechnology, School of Chemistry, The University of Manchester, Manchester M1 7DN, United Kingdom
| |
Collapse
|
4
|
Gingrich PW, Siegel JB, Tantillo DJ. Assessing Alkene Reactivity toward Cytochrome P450-Mediated Epoxidation through Localized Descriptors and Regression Modeling. J Chem Inf Model 2022; 62:1979-1987. [PMID: 35421306 DOI: 10.1021/acs.jcim.1c01567] [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/29/2022]
Abstract
The prediction of sites of epoxidation by cytochrome P450s during metabolism is particularly important in drug design, as epoxides are capable of alkylating biological macromolecules. Reliable methods are needed to quantitatively predict P450-mediated epoxidation barriers for inclusion in high-throughput screening campaigns alongside protein-ligand docking. Utilizing the fractional occupation number weighted density (FOD) and orbital-weighted Fukui index (fw+) as descriptors of local reactivity and a data set of 36 alkene epoxidation barriers computed with density functional theory (DFT), we developed and validated a multiple linear regression model for the reliable estimation of epoxidation barriers using only substrate structures as input. Using our recommended level of theory (GFN2-xTB//GFN-FF), mean absolute errors in the training and test sets were found to be 0.66 and 0.70 kcal/mol, respectively, with coefficients of determination of ca. 0.80. We demonstrate the utility of this approach on three known substrates of CYP101A1 and further show that this approach is inappropriate for particularly electron-rich alkenes. By employing a modern semiempirical method on force-field-generated geometries, the required descriptors can be calculated on the millisecond timescale per structure, making the approach well suited for incorporation into high-throughput methodologies alongside docking.
Collapse
Affiliation(s)
- Phillip W Gingrich
- Department of Chemistry, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Justin B Siegel
- Department of Chemistry, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Dean J Tantillo
- Department of Chemistry, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| |
Collapse
|
5
|
Zhang X, Xu M, Wu Z, Liu G, Tang Y, Li W. Assessment of CYP2C9 Structural Models for Site of Metabolism Prediction. ChemMedChem 2021; 16:1754-1763. [PMID: 33600055 DOI: 10.1002/cmdc.202000964] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Revised: 02/07/2021] [Indexed: 11/07/2022]
Abstract
Structure-based prediction of a compound's potential sites of metabolism (SOMs) mediated by cytochromes P450 (CYPs) is highly advantageous in the early stage of drug discovery. However, the accuracy of the SOMs prediction can be influenced by several factors. CYP2C9 is one of the major drug-metabolizing enzymes in humans and is responsible for the metabolism of ∼13 % of clinically used drugs. In this study, we systematically evaluated the effects of protein crystal structure models, scoring functions, heme forms, conserved active-site water molecules, and protein flexibility on SOMs prediction of CYP2C9 substrates. Our results demonstrated that, on average, ChemScore and GlideScore outperformed four other scoring functions: Vina, GoldScore, ChemPLP, and ASP. The performance of the crystal structure models with pentacoordinated heme was generally superior to that of the hexacoordinated iron-oxo heme (referred to as Compound I) models. Inclusion of the conserved active-site water molecule improved the prediction accuracy of GlideScore, but reduced the accuracy of ChemScore. In addition, the effect of the conserved water on SOMs prediction was found to be dependent on the receptor model and the substrate. We further found that one of snapshots from molecular dynamics simulations on the apo form can improve the prediction accuracy when compared to the crystal structural model.
Collapse
Affiliation(s)
- Xiaoxiao Zhang
- Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 20023, P. R. China
| | - Minjie Xu
- Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 20023, P. R. China
| | - Zengrui Wu
- Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 20023, P. R. China
| | - Guixia Liu
- Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 20023, P. R. China
| | - Yun Tang
- Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 20023, P. R. China
| | - Weihua Li
- Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai, 20023, P. R. China
| |
Collapse
|
6
|
Roy H, Nandi S. In-Silico Modeling in Drug Metabolism and Interaction: Current Strategies of Lead Discovery. Curr Pharm Des 2020; 25:3292-3305. [PMID: 31481001 DOI: 10.2174/1381612825666190903155935] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2019] [Accepted: 09/01/2019] [Indexed: 12/21/2022]
Abstract
BACKGROUND Drug metabolism is a complex mechanism of human body systems to detoxify foreign particles, chemicals, and drugs through bio alterations. It involves many biochemical reactions carried out by invivo enzyme systems present in the liver, kidney, intestine, lungs, and plasma. After drug administration, it crosses several biological membranes to reach into the target site for binding and produces the therapeutic response. After that, it may undergo detoxification and excretion to get rid of the biological systems. Most of the drugs and its metabolites are excreted through kidney via urination. Some drugs and their metabolites enter into intestinal mucosa and excrete through feces. Few of the drugs enter into hepatic circulation where they go into the intestinal tract. The drug leaves the liver via the bile duct and is excreted through feces. Therefore, the study of total methodology of drug biotransformation and interactions with various targets is costly. METHODS To minimize time and cost, in-silico algorithms have been utilized for lead-like drug discovery. Insilico modeling is the process where a computer model with a suitable algorithm is developed to perform a controlled experiment. It involves the combination of both in-vivo and in-vitro experimentation with virtual trials, eliminating the non-significant variables from a large number of variable parameters. Whereas, the major challenge for the experimenter is the selection and validation of the preferred model, as well as precise simulation in real physiological status. RESULTS The present review discussed the application of in-silico models to predict absorption, distribution, metabolism, and excretion (ADME) properties of drug molecules and also access the net rate of metabolism of a compound. CONCLUSION It helps with the identification of enzyme isoforms; which are likely to metabolize a compound, as well as the concentration dependence of metabolism and the identification of expected metabolites. In terms of drug-drug interactions (DDIs), models have been described for the inhibition of metabolism of one compound by another, and for the compound-dependent induction of drug-metabolizing enzymes.
Collapse
Affiliation(s)
- Harekrishna Roy
- Nirmala College of Pharmacy, Mangalagiri, Guntur, Affiliated to Acharya Nagarjuna University, Andhra Pradesh-522503, India
| | - Sisir Nandi
- Department of Pharmaceutical Chemistry, Global Institute of Pharmaceutical Education and Research, Affiliated to Uttarakhand Technical University, Kashipur-244713, India
| |
Collapse
|
7
|
Hu X, Maffucci I, Contini A. Advances in the Treatment of Explicit Water Molecules in Docking and Binding Free Energy Calculations. Curr Med Chem 2020; 26:7598-7622. [DOI: 10.2174/0929867325666180514110824] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2018] [Revised: 02/26/2018] [Accepted: 04/18/2018] [Indexed: 12/30/2022]
Abstract
Background:
The inclusion of direct effects mediated by water during the ligandreceptor
recognition is a hot-topic of modern computational chemistry applied to drug discovery
and development. Docking or virtual screening with explicit hydration is still debatable,
despite the successful cases that have been presented in the last years. Indeed, how to select
the water molecules that will be included in the docking process or how the included waters
should be treated remain open questions.
Objective:
In this review, we will discuss some of the most recent methods that can be used in
computational drug discovery and drug development when the effect of a single water, or of a
small network of interacting waters, needs to be explicitly considered.
Results:
Here, we analyse the software to aid the selection, or to predict the position, of water
molecules that are going to be explicitly considered in later docking studies. We also present
software and protocols able to efficiently treat flexible water molecules during docking, including
examples of applications. Finally, we discuss methods based on molecular dynamics
simulations that can be used to integrate docking studies or to reliably and efficiently compute
binding energies of ligands in presence of interfacial or bridging water molecules.
Conclusions:
Software applications aiding the design of new drugs that exploit water molecules,
either as displaceable residues or as bridges to the receptor, are constantly being developed.
Although further validation is needed, workflows that explicitly consider water will
probably become a standard for computational drug discovery soon.
Collapse
Affiliation(s)
- Xiao Hu
- Università degli Studi di Milano, Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Generale e Organica “A. Marchesini”, Via Venezian, 21 20133 Milano, Italy
| | - Irene Maffucci
- Pasteur, Département de Chimie, École Normale Supérieure, PSL Research University, Sorbonne Universités, UPMC Univ. Paris 06, CNRS, 75005 Paris, France
| | - Alessandro Contini
- Università degli Studi di Milano, Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Generale e Organica “A. Marchesini”, Via Venezian, 21 20133 Milano, Italy
| |
Collapse
|
8
|
Shadrack DM, Swai HS, Hassanali A. A computational study on the role of water and conformational fluctuations in Hsp90 in response to inhibitors. J Mol Graph Model 2019; 96:107510. [PMID: 31877402 DOI: 10.1016/j.jmgm.2019.107510] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2019] [Revised: 12/04/2019] [Accepted: 12/11/2019] [Indexed: 10/25/2022]
Abstract
Molecular chaperone Heat Shock Protein 90 (Hsp90) represents an interesting chemotherapeutic target for cancer treatments as it plays a role in cancer proliferation. Thus, continued effort to identify novel inhibitors of this target is an important task. Drug design using computational approach has gained significant attention in recent years. This work aims to propose docking protocols to re-purpose FDA-approved drugs targeting Hsp90. Sensitivity of results to different docking protocols such apo, holo and receptor ensembles (relaxed complex) structures, the role of water and conformational changes of Hsp90, are described. We show that the protein conformation and water have effects on drug binding. Holo relaxed complex receptors ensembles improves the binding energy of ligands to the protein. We also compare and contrast structural stability of three drugs namely: ezetimibe, pitavastatin and vilazodon in the Hsp90 protein. The results obtained serves as a possible basis towards developing Hsp90 inhibitors.
Collapse
Affiliation(s)
- Daniel M Shadrack
- Department of Health and Biomedical Sciences, School of Life Science and Bioengieering, The Nelson Mandela African Institution of Science and Technology, P.O.Box 447, Arusha, Tanzania; International Centre for Theoretical Physics, Strada Costiera, 11, 34151, Trieste, Italy; Department of Chemistry, FaNAS, St John's University of Tanzania, P.O.Box 47, Dodoma, Tanzania.
| | - Hulda S Swai
- Department of Health and Biomedical Sciences, School of Life Science and Bioengieering, The Nelson Mandela African Institution of Science and Technology, P.O.Box 447, Arusha, Tanzania
| | - Ali Hassanali
- International Centre for Theoretical Physics, Strada Costiera, 11, 34151, Trieste, Italy.
| |
Collapse
|
9
|
The role of hydration effects in 5-fluorouridine binding to SOD1: insight from a new 3D-RISM-KH based protocol for including structural water in docking simulations. J Comput Aided Mol Des 2019; 33:913-926. [PMID: 31686367 DOI: 10.1007/s10822-019-00239-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2019] [Accepted: 10/17/2019] [Indexed: 12/13/2022]
Abstract
Misfolded Cu/Zn superoxide dismutase enzyme (SOD1) shows prion-like propagation in neuronal cells leading to neurotoxic aggregates that are implicated in amyotrophic lateral sclerosis (ALS). Tryptophan-32 (W32) in SOD1 is part of a potential site for templated conversion of wild type SOD1. This W32 binding site is located on a convex, solvent exposed surface of the SOD1 suggesting that hydration effects can play an important role in ligand recognition and binding. A recent X-ray crystal structure has revealed that 5-Fluorouridine (5-FUrd) binds at the W32 binding site and can act as a pharmacophore scaffold for the development of anti-ALS drugs. In this study, a new protocol is developed to account for structural (non-displaceable) water molecules in docking simulations and successfully applied to predict the correct docked conformation binding modes of 5-FUrd at the W32 binding site. The docked configuration is within 0.58 Å (RMSD) of the observed configuration. The docking protocol involved calculating a hydration structure around SOD1 using molecular theory of solvation (3D-RISM-KH, 3D-Reference Interaction Site Model-Kovalenko-Hirata) whereby, non-displaceable water molecules are identified for docking simulations. This protocol was also used to analyze the hydrated structure of the W32 binding site and to explain the role of solvation in ligand recognition and binding to SOD1. Structural water molecules mediate hydrogen bonds between 5-FUrd and the receptor, and create an environment favoring optimal placement of 5-FUrd in the W32 binding site.
Collapse
|
10
|
Bruno A, Costantino G, Sartori L, Radi M. The In Silico Drug Discovery Toolbox: Applications in Lead Discovery and Optimization. Curr Med Chem 2019; 26:3838-3873. [PMID: 29110597 DOI: 10.2174/0929867324666171107101035] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2017] [Revised: 09/27/2017] [Accepted: 09/28/2017] [Indexed: 01/04/2023]
Abstract
BACKGROUND Discovery and development of a new drug is a long lasting and expensive journey that takes around 20 years from starting idea to approval and marketing of new medication. Despite R&D expenditures have been constantly increasing in the last few years, the number of new drugs introduced into market has been steadily declining. This is mainly due to preclinical and clinical safety issues, which still represent about 40% of drug discontinuation. To cope with this issue, a number of in silico techniques are currently being used for an early stage evaluation/prediction of potential safety issues, allowing to increase the drug-discovery success rate and reduce costs associated with the development of a new drug. METHODS In the present review, we will analyse the early steps of the drug-discovery pipeline, describing the sequence of steps from disease selection to lead optimization and focusing on the most common in silico tools used to assess attrition risks and build a mitigation plan. RESULTS A comprehensive list of widely used in silico tools, databases, and public initiatives that can be effectively implemented and used in the drug discovery pipeline has been provided. A few examples of how these tools can be problem-solving and how they may increase the success rate of a drug discovery and development program have been also provided. Finally, selected examples where the application of in silico tools had effectively contributed to the development of marketed drugs or clinical candidates will be given. CONCLUSION The in silico toolbox finds great application in every step of early drug discovery: (i) target identification and validation; (ii) hit identification; (iii) hit-to-lead; and (iv) lead optimization. Each of these steps has been described in details, providing a useful overview on the role played by in silico tools in the decision-making process to speed-up the discovery of new drugs.
Collapse
Affiliation(s)
- Agostino Bruno
- Experimental Therapeutics Unit, IFOM - The FIRC Institute for Molecular Oncology Foundation, Via Adamello 16 - 20139 Milano, Italy
| | - Gabriele Costantino
- Dipartimento di Scienze degli Alimenti e del Farmaco, Universita degli Studi di Parma, Viale delle Scienze, 27/A, 43124 Parma, Italy
| | - Luca Sartori
- Experimental Therapeutics Unit, IFOM - The FIRC Institute for Molecular Oncology Foundation, Via Adamello 16 - 20139 Milano, Italy
| | - Marco Radi
- Dipartimento di Scienze degli Alimenti e del Farmaco, Universita degli Studi di Parma, Viale delle Scienze, 27/A, 43124 Parma, Italy
| |
Collapse
|
11
|
Lu Q, Qi LW, Liu J. Improving protein–ligand binding prediction by considering the bridging water molecules in Autodock. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2019. [DOI: 10.1142/s0219633619500275] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Water plays a significant role in determining the protein–ligand binding modes, especially when water molecules are involved in mediating protein–ligand interactions, and these important water molecules are receiving more and more attention in recent years. Considering the effects of water molecules has gradually become a routine process for accurate description of the protein–ligand interactions. As a free docking program, Autodock has been most widely used in predicting the protein–ligand binding modes. However, whether the inclusion of water molecules in Autodock would improve its docking performance has not been systematically investigated. Here, we incorporate important bridging water molecules into Autodock program, and systematically investigate the effectiveness of these water molecules in protein–ligand docking. This approach was evaluated using 18 structurally diverse protein–ligand complexes, in which several water molecules bridge the protein–ligand interactions. Different treatment of water molecules were tested by using the fixed and rotatable water molecules, and a considerable improvement in successful docking simulations was found when including these water molecules. This study illustrates the necessity of inclusion of water molecules in Autodock docking, and emphasizes the importance of a proper treatment of water molecules in protein–ligand binding predictions.
Collapse
Affiliation(s)
- Qiangna Lu
- State Key Laboratory of Natural Medicines, Department of Chinese Medicine, Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, P. R. China
| | - Lian-Wen Qi
- State Key Laboratory of Natural Medicines, Department of Chinese Medicine, Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, P. R. China
| | - Jinfeng Liu
- State Key Laboratory of Natural Medicines, Department of Chinese Medicine, Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, P. R. China
| |
Collapse
|
12
|
Dmitriev AV, Lagunin AA, Karasev DА, Rudik AV, Pogodin PV, Filimonov DA, Poroikov VV. Prediction of Drug-Drug Interactions Related to Inhibition or Induction of Drug-Metabolizing Enzymes. Curr Top Med Chem 2019; 19:319-336. [PMID: 30674264 DOI: 10.2174/1568026619666190123160406] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Revised: 01/02/2019] [Accepted: 01/07/2019] [Indexed: 02/07/2023]
Abstract
Drug-drug interaction (DDI) is the phenomenon of alteration of the pharmacological activity of a drug(s) when another drug(s) is co-administered in cases of so-called polypharmacy. There are three types of DDIs: pharmacokinetic (PK), pharmacodynamic, and pharmaceutical. PK is the most frequent type of DDI, which often appears as a result of the inhibition or induction of drug-metabolising enzymes (DME). In this review, we summarise in silico methods that may be applied for the prediction of the inhibition or induction of DMEs and describe appropriate computational methods for DDI prediction, showing the current situation and perspectives of these approaches in medicinal and pharmaceutical chemistry. We review sources of information on DDI, which can be used in pharmaceutical investigations and medicinal practice and/or for the creation of computational models. The problem of the inaccuracy and redundancy of these data are discussed. We provide information on the state-of-the-art physiologically- based pharmacokinetic modelling (PBPK) approaches and DME-based in silico methods. In the section on ligand-based methods, we describe pharmacophore models, molecular field analysis, quantitative structure-activity relationships (QSAR), and similarity analysis applied to the prediction of DDI related to the inhibition or induction of DME. In conclusion, we discuss the problems of DDI severity assessment, mention factors that influence severity, and highlight the issues, perspectives and practical using of in silico methods.
Collapse
Affiliation(s)
| | - Alexey A Lagunin
- Institute of Biomedical Chemistry, Moscow, Russian Federation.,Pirogov Russian National Research Medical University, Moscow, RussiaN Federation
| | | | | | - Pavel V Pogodin
- Institute of Biomedical Chemistry, Moscow, Russian Federation
| | | | | |
Collapse
|
13
|
Microsecond MD simulations of human CYP2D6 wild-type and five allelic variants reveal mechanistic insights on the function. PLoS One 2018; 13:e0202534. [PMID: 30133539 PMCID: PMC6104999 DOI: 10.1371/journal.pone.0202534] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2018] [Accepted: 08/03/2018] [Indexed: 11/19/2022] Open
Abstract
Characterization of cytochrome P450 2D6 (CYP2D6) and the impact of the major identified allelic variants on the activity of one of the most dominating drug-metabolising enzymes is essential to increase drug safety and avoid adverse reactions. Microsecond molecular dynamics simulations have been performed to capture the dynamic signatures of this complex enzyme and five allelic variants with diverse enzymatic activity. In addition to the apo simulations, three substrates (bufuralol, veliparib and tamoxifen) and two inhibitors (prinomastat and quinidine) were included to explore their influence on the structure and dynamical features of the enzyme. Our results indicate that the altered enzyme activity can be attributed to changes in the hydrogen bonding network within the active site, and local structural differences in flexibility, position and shape of the binding pocket. In particular, the increased (CYP2D6*53) or the decreased (CYP2D6*17) activity seems to be related to a change in dynamics of mainly the BC loop due to a modified hydrogen bonding network around this region. In addition, the smallest active site volume was found for CYP2D6*4 (no activity). CYP2D6*2 (normal activity) showed no major differences in dynamic behaviour compared to the wild-type.
Collapse
|
14
|
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
|
15
|
Rudling A, Orro A, Carlsson J. Prediction of Ordered Water Molecules in Protein Binding Sites from Molecular Dynamics Simulations: The Impact of Ligand Binding on Hydration Networks. J Chem Inf Model 2018; 58:350-361. [PMID: 29308882 PMCID: PMC6716772 DOI: 10.1021/acs.jcim.7b00520] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Water
plays a major role in ligand binding and is attracting increasing
attention in structure-based drug design. Water molecules can make
large contributions to binding affinity by bridging protein–ligand
interactions or by being displaced upon complex formation, but these
phenomena are challenging to model at the molecular level. Herein,
networks of ordered water molecules in protein binding sites were
analyzed by clustering of molecular dynamics (MD) simulation trajectories.
Locations of ordered waters (hydration sites) were first identified
from simulations of high resolution crystal structures of 13 protein–ligand
complexes. The MD-derived hydration sites reproduced 73% of the binding
site water molecules observed in the crystal structures. If the simulations
were repeated without the cocrystallized ligands, a majority (58%)
of the crystal waters in the binding sites were still predicted. In
addition, comparison of the hydration sites obtained from simulations
carried out in the absence of ligands to those identified for the
complexes revealed that the networks of ordered water molecules were
preserved to a large extent, suggesting that the locations of waters
in a protein–ligand interface are mainly dictated by the protein.
Analysis of >1000 crystal structures showed that hydration sites
bridged
protein–ligand interactions in complexes with different ligands,
and those with high MD-derived occupancies were more likely to correspond
to experimentally observed ordered water molecules. The results demonstrate
that ordered water molecules relevant for modeling of protein–ligand
complexes can be identified from MD simulations. Our findings could
contribute to development of improved methods for structure-based
virtual screening and lead optimization.
Collapse
Affiliation(s)
- Axel Rudling
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Adolfo Orro
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Jens Carlsson
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC , Box 596, SE-751 24 Uppsala, Sweden
| |
Collapse
|
16
|
Dixit VA, Deshpande S. Advances in Computational Prediction of Regioselective and Isoform-Specific Drug Metabolism Catalyzed by CYP450s. ChemistrySelect 2016. [DOI: 10.1002/slct.201601051] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Vaibhav A. Dixit
- Department of Pharmaceutical Chemistry; School of Pharmacy and Technology Management (SPTM), Shri Vile Parle Kelavani Mandal's (SVKM's) Narsee Monjee Institute of Management Studies (NMIMS), Mukesh Patel Technology Park, Babulde, Bank of Tapi River; Mumbai-Agra Road Shirpur, Dist. Dhule−425405 India
| | - Shirish Deshpande
- Department of Pharmaceutical Chemistry; School of Pharmacy and Technology Management (SPTM), Shri Vile Parle Kelavani Mandal's (SVKM's) Narsee Monjee Institute of Management Studies (NMIMS), Mukesh Patel Technology Park, Babulde, Bank of Tapi River; Mumbai-Agra Road Shirpur, Dist. Dhule−425405 India
| |
Collapse
|
17
|
Nair PC, McKinnon RA, Miners JO. Cytochrome P450 structure–function: insights from molecular dynamics simulations. Drug Metab Rev 2016; 48:434-52. [DOI: 10.1080/03602532.2016.1178771] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
18
|
Li J, Cai J, Su H, Du H, Zhang J, Ding S, Liu G, Tang Y, Li W. Effects of protein flexibility and active site water molecules on the prediction of sites of metabolism for cytochrome P450 2C19 substrates. MOLECULAR BIOSYSTEMS 2016; 12:868-78. [DOI: 10.1039/c5mb00784d] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Structure-based prediction of sites of metabolism (SOMs) mediated by cytochrome P450s (CYPs) is of great interest in drug discovery and development.
Collapse
Affiliation(s)
- Junhao Li
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Jinya Cai
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Haixia Su
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Hanwen Du
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Juan Zhang
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Shihui Ding
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Guixia Liu
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Yun Tang
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| | - Weihua Li
- Shanghai Key Laboratory of New Drug Design
- School of Pharmacy
- East China University of Science and Technology
- Shanghai 200237
- China
| |
Collapse
|
19
|
Li L, Xu W, Lü Q. Improving protein-ligand docking with flexible interfacial water molecules using SWRosettaLigand. J Mol Model 2015; 21:294. [DOI: 10.1007/s00894-015-2834-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2014] [Accepted: 10/09/2015] [Indexed: 01/07/2023]
|
20
|
Spyrakis F, Cavasotto CN. Open challenges in structure-based virtual screening: Receptor modeling, target flexibility consideration and active site water molecules description. Arch Biochem Biophys 2015; 583:105-19. [DOI: 10.1016/j.abb.2015.08.002] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2015] [Revised: 08/03/2015] [Accepted: 08/03/2015] [Indexed: 01/05/2023]
|
21
|
Cheng CY, Olijve LLC, Kausik R, Han S. Cholesterol enhances surface water diffusion of phospholipid bilayers. J Chem Phys 2015; 141:22D513. [PMID: 25494784 DOI: 10.1063/1.4897539] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
Elucidating the physical effect of cholesterol (Chol) on biological membranes is necessary towards rationalizing their structural and functional role in cell membranes. One of the debated questions is the role of hydration water in Chol-embedding lipid membranes, for which only little direct experimental data are available. Here, we study the hydration dynamics in a series of Chol-rich and depleted bilayer systems using an approach termed (1)H Overhauser dynamic nuclear polarization (ODNP) NMR relaxometry that enables the sensitive and selective determination of water diffusion within 5-10 Å of a nitroxide-based spin label, positioned off the surface of the polar headgroups or within the nonpolar core of lipid membranes. The Chol-rich membrane systems were prepared from mixtures of Chol, dipalmitoyl phosphatidylcholine and/or dioctadecyl phosphatidylcholine lipid that are known to form liquid-ordered, raft-like, domains. Our data reveal that the translational diffusion of local water on the surface and within the hydrocarbon volume of the bilayer is significantly altered, but in opposite directions: accelerated on the membrane surface and dramatically slowed in the bilayer interior with increasing Chol content. Electron paramagnetic resonance (EPR) lineshape analysis shows looser packing of lipid headgroups and concurrently tighter packing in the bilayer core with increasing Chol content, with the effects peaking at lipid compositions reported to form lipid rafts. The complementary capability of ODNP and EPR to site-specifically probe the hydration dynamics and lipid ordering in lipid membrane systems extends the current understanding of how Chol may regulate biological processes. One possible role of Chol is the facilitation of interactions between biological constituents and the lipid membrane through the weakening or disruption of strong hydrogen-bond networks of the surface hydration layers that otherwise exert stronger repulsive forces, as reflected in faster surface water diffusivity. Another is the concurrent tightening of lipid packing that reduces passive, possibly unwanted, diffusion of ions and water across the bilayer.
Collapse
Affiliation(s)
- Chi-Yuan Cheng
- Department of Chemistry and Biochemistry and Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
| | - Luuk L C Olijve
- Laboratory of Macromolecular and Organic Chemistry and Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600 MB, Eindhoven, The Netherlands
| | - Ravinath Kausik
- Department of Chemistry and Biochemistry and Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
| | - Songi Han
- Department of Chemistry and Biochemistry and Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
| |
Collapse
|
22
|
Olsen L, Oostenbrink C, Jørgensen FS. Prediction of cytochrome P450 mediated metabolism. Adv Drug Deliv Rev 2015; 86:61-71. [PMID: 25958010 DOI: 10.1016/j.addr.2015.04.020] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2014] [Revised: 03/30/2015] [Accepted: 04/27/2015] [Indexed: 10/23/2022]
Abstract
Cytochrome P450 enzymes (CYPs) form one of the most important enzyme families involved in the metabolism of xenobiotics. CYPs comprise many isoforms, which catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be formed. However, it is often hard to rationalize what metabolites these enzymes generate. In recent years, many different in silico approaches have been developed to predict binding or regioselective product formation for the different CYP isoforms. These comprise ligand-based methods that are trained on experimental CYP data and structure-based methods that consider how the substrate is oriented in the active site or/and how reactive the part of the substrate that is accessible to the heme group is. We will review key aspects for various approaches that are available to predict binding and site of metabolism (SOM), what outcome can be expected from the predictions, and how they could potentially be improved.
Collapse
|
23
|
Scharkoi O, Becker R, Esslinger S, Weber M, Nehls I. Predicting sites of cytochrome P450-mediated hydroxylation applied to CYP3A4 and hexabromocyclododecane. MOLECULAR SIMULATION 2015. [DOI: 10.1080/08927022.2014.898845] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
24
|
Danishuddin M, Khan AU. Structure based virtual screening to discover putative drug candidates: Necessary considerations and successful case studies. Methods 2015; 71:135-45. [DOI: 10.1016/j.ymeth.2014.10.019] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2014] [Revised: 09/25/2014] [Accepted: 10/17/2014] [Indexed: 12/19/2022] Open
|
25
|
Practical Considerations in Virtual Screening and Molecular Docking. EMERGING TRENDS IN COMPUTATIONAL BIOLOGY, BIOINFORMATICS, AND SYSTEMS BIOLOGY 2015. [PMCID: PMC7173576 DOI: 10.1016/b978-0-12-802508-6.00027-2] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
Molecular docking has become an important common component of the drug discovery toolbox, and its relative low-cost implications and perceived simplicity of use has stimulated an everincreasing popularity within academic communities. The inherent “garbage-in-garbage-out” defect of molecular docking, however, leads a lot of researchers to dedicate countless hours to the identification of hit compounds that later prove to be inactive. Several considerations that can greatly improve the success and enrichment of true bioactive hit compounds are commonly overlooked at the initial stages of a molecular docking study. This chapter will cover several of these considerations, including protonation states, active site waters, separating actives from decoys, consensus docking and molecular mechanics generalized-Born/surface area (MM-GBSA) rescoring, and incorporation of pharmacophoric constraints, in an attempt to clarify what is, in fact, very complicated and inherent difficulties of a structure-based drug design study.
Collapse
|
26
|
Sheng Y, Chen Y, Wang L, Liu G, Li W, Tang Y. Effects of protein flexibility on the site of metabolism prediction for CYP2A6 substrates. J Mol Graph Model 2014; 54:90-9. [DOI: 10.1016/j.jmgm.2014.09.005] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2014] [Revised: 09/25/2014] [Accepted: 09/30/2014] [Indexed: 12/01/2022]
|
27
|
Lenselink EB, Beuming T, Sherman W, van Vlijmen HWT, IJzerman AP. Selecting an optimal number of binding site waters to improve virtual screening enrichments against the adenosine A2A receptor. J Chem Inf Model 2014; 54:1737-46. [PMID: 24835542 DOI: 10.1021/ci5000455] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
A major challenge in structure-based virtual screening (VS) involves the treatment of explicit water molecules during docking in order to improve the enrichment of active compounds over decoys. Here we have investigated this in the context of the adenosine A2A receptor, where water molecules have previously been shown to be important for achieving high enrichment rates with docking, and where the positions of some binding site waters are known from a high-resolution crystal structure. The effect of these waters (both their presence and orientations) on VS enrichment was assessed using a carefully curated set of 299 high affinity A2A antagonists and 17,337 decoys. We show that including certain crystal waters greatly improves VS enrichment and that optimization of water hydrogen positions is needed in order to achieve the best results. We also show that waters derived from a molecular dynamics simulation - without any knowledge of crystallographic waters - can improve enrichments to a similar degree as the crystallographic waters, which makes this strategy applicable to structures without experimental knowledge of water positions. Finally, we used decision trees to select an ensemble of structures with different water molecule positions and orientations that outperforms any single structure with water molecules. The approach presented here is validated against independent test sets of A2A receptor antagonists and decoys from the literature. In general, this water optimization strategy could be applied to any target with waters-mediated protein-ligand interactions.
Collapse
Affiliation(s)
- Eelke B Lenselink
- Division of Medicinal Chemistry, Leiden Academic Centre for Drug Research, Leiden University , Einsteinweg 55, 2333 CC Leiden, The Netherlands
| | | | | | | | | |
Collapse
|
28
|
Hayes C, Ansbro D, Kontoyianni M. Elucidating Substrate Promiscuity in the Human Cytochrome 3A4. J Chem Inf Model 2014; 54:857-69. [DOI: 10.1021/ci4006782] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Christina Hayes
- Department of Pharmaceutical
Sciences, School of Pharmacy, Southern Illinois University Edwardsville, Edwardsville, Illinois 62034, United States
| | - Daniel Ansbro
- Department of Pharmaceutical
Sciences, School of Pharmacy, Southern Illinois University Edwardsville, Edwardsville, Illinois 62034, United States
| | - Maria Kontoyianni
- Department of Pharmaceutical
Sciences, School of Pharmacy, Southern Illinois University Edwardsville, Edwardsville, Illinois 62034, United States
| |
Collapse
|
29
|
Rudik AV, Dmitriev AV, Lagunin AA, Filimonov DA, Poroikov VV. Metabolism site prediction based on xenobiotic structural formulas and PASS prediction algorithm. J Chem Inf Model 2014; 54:498-507. [PMID: 24417355 DOI: 10.1021/ci400472j] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
A new ligand-based method for the prediction of sites of metabolism (SOMs) for xenobiotics has been developed on the basis of the LMNA (labeled multilevel neighborhoods of atom) descriptors and the PASS (prediction of activity spectra for substances) algorithm and applied to predict the SOMs of the 1A2, 2C9, 2C19, 2D6, and 3A4 isoforms of cytochrome P450. An average IAP (invariant accuracy of prediction) of SOMs calculated by the leave-one-out cross-validation procedure was 0.89 for the developed method. The external validation was made with evaluation sets containing data on biotransformations for 57 cardiovascular drugs. An average IAP of regioselectivity for evaluation sets was 0.83. It was shown that the proposed method exceeds accuracy of SOM prediction by RS-Predictor for CYP 1A2, 2D6, 2C9, 2C19, and 3A4 and is comparable to or better than SMARTCyp for CYP 2C9 and 2D6.
Collapse
Affiliation(s)
- Anastasia V Rudik
- Orekhovich Institute of Biomedical Chemistry of the Russian Academy of Medical Sciences , Building 10/8, Pogodinskaya Str., Moscow, 119121, Russia
| | | | | | | | | |
Collapse
|
30
|
Perić-Hassler L, Stjernschantz E, Oostenbrink C, Geerke DP. CYP 2D6 binding affinity predictions using multiple ligand and protein conformations. Int J Mol Sci 2013; 14:24514-30. [PMID: 24351831 PMCID: PMC3876125 DOI: 10.3390/ijms141224514] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2013] [Revised: 11/28/2013] [Accepted: 12/04/2013] [Indexed: 01/18/2023] Open
Abstract
Because of the large flexibility and malleability of Cytochrome P450 enzymes (CYPs), in silico prediction of CYP binding affinities to drugs and other xenobiotic compounds is a true challenge. In the current work, we use an iterative linear interaction energy (LIE) approach to compute CYP binding affinities from molecular dynamics (MD) simulation. In order to improve sampling of conformational space, we combine results from simulations starting with different relevant protein-ligand geometries. For calculated binding free energies of a set of thiourea compounds binding to the flexible CYP 2D6 isoform, improved correlation with experiment was obtained by combining results of MD runs starting from distinct protein conformations and ligand-binding orientations. This accuracy was obtained from relatively short MD simulations, which makes our approach computationally attractive for automated calculations of ligand-binding affinities to flexible proteins such as CYPs.
Collapse
Affiliation(s)
- Lovorka Perić-Hassler
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands; E-Mail:
| | - Eva Stjernschantz
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands; E-Mail:
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria; E-Mail:
| | - Daan P. Geerke
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands; E-Mail:
| |
Collapse
|
31
|
Huang TW, Zaretzki J, Bergeron C, Bennett KP, Breneman CM. DR-predictor: incorporating flexible docking with specialized electronic reactivity and machine learning techniques to predict CYP-mediated sites of metabolism. J Chem Inf Model 2013; 53:3352-66. [PMID: 24261543 DOI: 10.1021/ci4004688] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Computational methods that can identify CYP-mediated sites of metabolism (SOMs) of drug-like compounds have become required tools for early stage lead optimization. In recent years, methods that combine CYP binding site features with CYP/ligand binding information have been sought in order to increase the prediction accuracy of such hybrid models over those that use only one representation. Two challenges that any hybrid ligand/structure-based method must overcome are (1) identification of the best binding pose for a specific ligand with a given CYP and (2) appropriately incorporating the results of docking with ligand reactivity. To address these challenges we have created Docking-Regioselectivity-Predictor (DR-Predictor)--a method that incorporates flexible docking-derived information with specialized electronic reactivity and multiple-instance-learning methods to predict CYP-mediated SOMs. In this study, the hybrid ligand-structure-based DR-Predictor method was tested on substrate sets for CYP 1A2 and CYP 2A6. For these data, the DR-Predictor model was found to identify the experimentally observed SOM within the top two predicted rank-positions for 86% of the 261 1A2 substrates and 83% of the 100 2A6 substrates. Given the accuracy and extendibility of the DR-Predictor method, we anticipate that it will further facilitate the prediction of CYP metabolism liabilities and aid in in-silico ADMET assessment of novel structures.
Collapse
Affiliation(s)
- Tao-wei Huang
- Department of Chemistry and Chemical Biology and §Department of Mathematics, Rensselaer Polytechnic Institute , Troy, New York 12180, United States
| | | | | | | | | |
Collapse
|
32
|
Accounting for Target Flexibility and Water Molecules by Docking to Ensembles of Target Structures: The HCV NS5B Palm Site I Inhibitors Case Study. J Chem Inf Model 2013; 54:481-97. [DOI: 10.1021/ci400367m] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
33
|
Kumar A, Zhang KYJ. Investigation on the Effect of Key Water Molecules on Docking Performance in CSARdock Exercise. J Chem Inf Model 2013; 53:1880-92. [DOI: 10.1021/ci400052w] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ashutosh Kumar
- Zhang Initiative
Research Unit, RIKEN, 2-1 Hirosawa,
Wako, Saitama 351-0198, Japan
| | - Kam Y. J. Zhang
- Zhang Initiative
Research Unit, RIKEN, 2-1 Hirosawa,
Wako, Saitama 351-0198, Japan
| |
Collapse
|
34
|
Madhavi Sastry G, Adzhigirey M, Day T, Annabhimoju R, Sherman W. Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments. J Comput Aided Mol Des 2013; 27:221-34. [DOI: 10.1007/s10822-013-9644-8] [Citation(s) in RCA: 2913] [Impact Index Per Article: 264.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2012] [Accepted: 04/03/2013] [Indexed: 12/11/2022]
|
35
|
Yuriev E, Ramsland PA. Latest developments in molecular docking: 2010-2011 in review. J Mol Recognit 2013; 26:215-39. [PMID: 23526775 DOI: 10.1002/jmr.2266] [Citation(s) in RCA: 193] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2012] [Revised: 01/16/2013] [Accepted: 01/19/2013] [Indexed: 12/28/2022]
Affiliation(s)
- Elizabeth Yuriev
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences; Monash University; Parkville; VIC; 3052; Australia
| | | |
Collapse
|
36
|
Songtawee N, Gleeson MP, Choowongkomon K. Computational study of EGFR inhibition: molecular dynamics studies on the active and inactive protein conformations. J Mol Model 2012; 19:497-509. [PMID: 22955422 DOI: 10.1007/s00894-012-1559-0] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2011] [Accepted: 08/02/2012] [Indexed: 10/27/2022]
Abstract
The structural diversity observed across protein kinases, resulting in subtly different active site cavities, is highly desirable in the pursuit of selective inhibitors, yet it can also be a hindrance from a structure-based design perspective. An important challenge in structure-based design is to better understand the dynamic nature of protein kinases and the underlying reasons for specific conformational preferences in the presence of different inhibitors. To investigate this issue, we performed molecular dynamics simulation on both the active and inactive wild type epidermal growth factor receptor (EGFR) protein with both type-I and type-II inhibitors. Our goal is to better understand the origin of the two distinct EGFR protein conformations, their dynamic differences, and their relative preference for Type-I inhibitors such as gefitinib and Type-II inhibitors such as lapatinib. We discuss the implications of protein dynamics from a structure-based design perspective.
Collapse
Affiliation(s)
- Napat Songtawee
- Department of Biochemistry, Faculty of Science, Kasetsart University, 50 Phaholyothin Rd, Chatuchak, Bangkok, 10900, Thailand
| | | | | |
Collapse
|
37
|
Zaretzki J, Rydberg P, Bergeron C, Bennett KP, Olsen L, Breneman CM. RS-Predictor models augmented with SMARTCyp reactivities: robust metabolic regioselectivity predictions for nine CYP isozymes. J Chem Inf Model 2012; 52:1637-59. [PMID: 22524152 DOI: 10.1021/ci300009z] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
RS-Predictor is a tool for creating pathway-independent, isozyme-specific, site of metabolism (SOM) prediction models using any set of known cytochrome P450 (CYP) substrates and metabolites. Until now, the RS-Predictor method was only trained and validated on CYP 3A4 data, but in the present study, we report on the versatility the RS-Predictor modeling paradigm by creating and testing regioselectivity models for substrates of the nine most important CYP isozymes. Through curation of source literature, we have assembled 680 substrates distributed among CYPs 1A2, 2A6, 2B6, 2C19, 2C8, 2C9, 2D6, 2E1, and 3A4, the largest publicly accessible collection of P450 ligands and metabolites released to date. A comprehensive investigation into the importance of different descriptor classes for identifying the regioselectivity mediated by each isozyme is made through the generation of multiple independent RS-Predictor models for each set of isozyme substrates. Two of these models include a density functional theory (DFT) reactivity descriptor derived from SMARTCyp. Optimal combinations of RS-Predictor and SMARTCyp are shown to have stronger performance than either method alone, while also exceeding the accuracy of the commercial regioselectivity prediction methods distributed by Optibrium and Schrödinger, correctly identifying a large proportion of the metabolites in each substrate set within the top two rank-positions: 1A2 (83.0%), 2A6 (85.7%), 2B6 (82.1%), 2C19 (86.2%), 2C8 (83.8%), 2C9 (84.5%), 2D6 (85.9%), 2E1 (82.8%), 3A4 (82.3%), and merged (86.0%). Comprehensive datamining of each substrate set and careful statistical analyses of the predictions made by the different models revealed new insights into molecular features that control metabolic regioselectivity and enable accurate prospective prediction of likely SOMs.
Collapse
Affiliation(s)
- Jed Zaretzki
- Department of Chemistry and Chemical Biology, Rensselaer Polytechnic Institute, Troy, New York 12180, USA
| | | | | | | | | | | |
Collapse
|
38
|
Berka K, Anzenbacherová E, Hendrychová T, Lange R, Mašek V, Anzenbacher P, Otyepka M. Binding of quinidine radically increases the stability and decreases the flexibility of the cytochrome P450 2D6 active site. J Inorg Biochem 2012; 110:46-50. [DOI: 10.1016/j.jinorgbio.2012.02.010] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2011] [Revised: 01/13/2012] [Accepted: 02/15/2012] [Indexed: 11/25/2022]
|
39
|
Cheng T, Li Q, Zhou Z, Wang Y, Bryant SH. Structure-based virtual screening for drug discovery: a problem-centric review. AAPS J 2012; 14:133-41. [PMID: 22281989 PMCID: PMC3282008 DOI: 10.1208/s12248-012-9322-0] [Citation(s) in RCA: 352] [Impact Index Per Article: 29.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2011] [Accepted: 01/04/2012] [Indexed: 11/30/2022] Open
Abstract
Structure-based virtual screening (SBVS) has been widely applied in early-stage drug discovery. From a problem-centric perspective, we reviewed the recent advances and applications in SBVS with a special focus on docking-based virtual screening. We emphasized the researchers' practical efforts in real projects by understanding the ligand-target binding interactions as a premise. We also highlighted the recent progress in developing target-biased scoring functions by optimizing current generic scoring functions toward certain target classes, as well as in developing novel ones by means of machine learning techniques.
Collapse
Affiliation(s)
- Tiejun Cheng
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, 8600 Rockville Pike, Bethesda, Maryland 20894 USA
| | - Qingliang Li
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, 8600 Rockville Pike, Bethesda, Maryland 20894 USA
| | - Zhigang Zhou
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, 8600 Rockville Pike, Bethesda, Maryland 20894 USA
| | - Yanli Wang
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, 8600 Rockville Pike, Bethesda, Maryland 20894 USA
| | - Stephen H. Bryant
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, 8600 Rockville Pike, Bethesda, Maryland 20894 USA
| |
Collapse
|
40
|
Ross GA, Morris GM, Biggin PC. Rapid and accurate prediction and scoring of water molecules in protein binding sites. PLoS One 2012; 7:e32036. [PMID: 22396746 PMCID: PMC3291545 DOI: 10.1371/journal.pone.0032036] [Citation(s) in RCA: 130] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2011] [Accepted: 01/18/2012] [Indexed: 12/21/2022] Open
Abstract
Water plays a critical role in ligand-protein interactions. However, it is still challenging to predict accurately not only where water molecules prefer to bind, but also which of those water molecules might be displaceable. The latter is often seen as a route to optimizing affinity of potential drug candidates. Using a protocol we call WaterDock, we show that the freely available AutoDock Vina tool can be used to predict accurately the binding sites of water molecules. WaterDock was validated using data from X-ray crystallography, neutron diffraction and molecular dynamics simulations and correctly predicted 97% of the water molecules in the test set. In addition, we combined data-mining, heuristic and machine learning techniques to develop probabilistic water molecule classifiers. When applied to WaterDock predictions in the Astex Diverse Set of protein ligand complexes, we could identify whether a water molecule was conserved or displaced to an accuracy of 75%. A second model predicted whether water molecules were displaced by polar groups or by non-polar groups to an accuracy of 80%. These results should prove useful for anyone wishing to undertake rational design of new compounds where the displacement of water molecules is being considered as a route to improved affinity.
Collapse
Affiliation(s)
- Gregory A. Ross
- Structural Bioinformatics and Computational Biochemistry, University of Oxford, Oxford, United Kingdom
| | | | - Philip C. Biggin
- Structural Bioinformatics and Computational Biochemistry, University of Oxford, Oxford, United Kingdom
- * E-mail:
| |
Collapse
|
41
|
Kirchmair J, Williamson MJ, Tyzack JD, Tan L, Bond PJ, Bender A, Glen RC. Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms. J Chem Inf Model 2012; 52:617-48. [PMID: 22339582 PMCID: PMC3317594 DOI: 10.1021/ci200542m] [Citation(s) in RCA: 187] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Metabolism of xenobiotics remains a central challenge
for the discovery
and development of drugs, cosmetics, nutritional supplements, and
agrochemicals. Metabolic transformations are frequently related to
the incidence of toxic effects that may result from the emergence
of reactive species, the systemic accumulation of metabolites, or
by induction of metabolic pathways. Experimental investigation of
the metabolism of small organic molecules is particularly resource
demanding; hence, computational methods are of considerable interest
to complement experimental approaches. This review provides a broad
overview of structure- and ligand-based computational methods for
the prediction of xenobiotic metabolism. Current computational approaches
to address xenobiotic metabolism are discussed from three major perspectives:
(i) prediction of sites of metabolism (SOMs), (ii) elucidation of
potential metabolites and their chemical structures, and (iii) prediction
of direct and indirect effects of xenobiotics on metabolizing enzymes,
where the focus is on the cytochrome P450 (CYP) superfamily of enzymes,
the cardinal xenobiotics metabolizing enzymes. For each of these domains,
a variety of approaches and their applications are systematically
reviewed, including expert systems, data mining approaches, quantitative
structure–activity relationships (QSARs), and machine learning-based
methods, pharmacophore-based algorithms, shape-focused techniques,
molecular interaction fields (MIFs), reactivity-focused techniques,
protein–ligand docking, molecular dynamics (MD) simulations,
and combinations of methods. Predictive metabolism is a developing
area, and there is still enormous potential for improvement. However,
it is clear that the combination of rapidly increasing amounts of
available ligand- and structure-related experimental data (in particular,
quantitative data) with novel and diverse simulation and modeling
approaches is accelerating the development of effective tools for
prediction of in vivo metabolism, which is reflected by the diverse
and comprehensive data sources and methods for metabolism prediction
reviewed here. This review attempts to survey the range and scope
of computational methods applied to metabolism prediction and also
to compare and contrast their applicability and performance.
Collapse
Affiliation(s)
- Johannes Kirchmair
- Unilever Centre for Molecular Science Informatics, Department of Chemistry, University of Cambridge, Lensfield Road, CB2 1EW, Cambridge, United Kingdom
| | | | | | | | | | | | | |
Collapse
|
42
|
Affiliation(s)
- David J Huggins
- Department of Oncology, Hutchison/MRC Research Centre, University of Cambridge, Hills Road, Cambridge, CB2 0XZ, United Kingdom.
| | | | | |
Collapse
|
43
|
Moroy G, Martiny VY, Vayer P, Villoutreix BO, Miteva MA. Toward in silico structure-based ADMET prediction in drug discovery. Drug Discov Today 2011; 17:44-55. [PMID: 22056716 DOI: 10.1016/j.drudis.2011.10.023] [Citation(s) in RCA: 164] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2011] [Revised: 10/07/2011] [Accepted: 10/21/2011] [Indexed: 12/12/2022]
Abstract
Quantitative structure-activity relationship (QSAR) methods and related approaches have been used to investigate the molecular features that influence the absorption, distribution, metabolism, excretion and toxicity (ADMET) of drugs. As the three-dimensional structures of several major ADMET proteins become available, structure-based (docking-scoring) computations can be carried out to complement or to go beyond QSAR studies. Applying docking-scoring methods to ADMET proteins is a challenging process because they usually have a large and flexible binding cavity; however, promising results relating to metabolizing enzymes have been reported. After reviewing current trends in the field we applied structure-based methods in the context of receptor flexibility in a case study involving the phase II metabolizing sulfotransferases. Overall, the explored concepts and results suggested that structure-based ADMET profiling will probably join the mainstream during the coming years.
Collapse
Affiliation(s)
- Gautier Moroy
- Inserm UMR-S 973, Molécules Thérapeutiques In Silico, Université Paris Diderot, Sorbonne Paris Cité, 35 Rue Helene Brion, 75013 Paris, France
| | | | | | | | | |
Collapse
|
44
|
Braga RC, Alves VM, Fraga CAM, Barreiro EJ, de Oliveira V, Andrade CH. Combination of docking, molecular dynamics and quantum mechanical calculations for metabolism prediction of 3,4-methylenedioxybenzoyl-2-thienylhydrazone. J Mol Model 2011; 18:2065-78. [DOI: 10.1007/s00894-011-1219-9] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2011] [Accepted: 08/09/2011] [Indexed: 11/29/2022]
|
45
|
Moors SLC, Vos AM, Cummings MD, Van Vlijmen H, Ceulemans A. Structure-Based Site of Metabolism Prediction for Cytochrome P450 2D6. J Med Chem 2011; 54:6098-105. [DOI: 10.1021/jm2006468] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Samuel L. C. Moors
- Department of Chemistry, Katholieke Universiteit Leuven, Celestijnenlaan 200F, 3001 Heverlee, Belgium
| | - Ann M. Vos
- Tibotec BVBA, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Maxwell D. Cummings
- Johnson & Johnson Pharmaceutical Research & Development, Welsh and McKean Roads, Springhouse, Pennsylvania 19477, United States
| | | | - Arnout Ceulemans
- Department of Chemistry, Katholieke Universiteit Leuven, Celestijnenlaan 200F, 3001 Heverlee, Belgium
| |
Collapse
|
46
|
Schyman P, Lai W, Chen H, Wang Y, Shaik S. The directive of the protein: how does cytochrome P450 select the mechanism of dopamine formation? J Am Chem Soc 2011; 133:7977-84. [PMID: 21539368 DOI: 10.1021/ja201665x] [Citation(s) in RCA: 200] [Impact Index Per Article: 15.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Dopamine can be generated from tyramine via arene hydroxylation catalyzed by a cytochrome P450 enzyme (CYP2D6). Our quantum mechanical/molecular mechanical (QM/MM) results reveal the decisive impact of the protein in selecting the 'best' reaction mechanism. Instead of the traditional Meisenheimer-complex mechanism, the study reveals a mechanism involving an initial hydrogen atom transfer from the phenolic hydroxyl group of the tyramine to the iron-oxo of the compound I (Cpd I), followed by a ring-π radical rebound that eventually leads to dopamine by keto-enol rearrangement. This mechanism is not viable in the gas phase since the O-H bond activation by Cpd I is endothermic and the process does not form a stable intermediate. By contrast, the in-protein reaction has a low barrier and is exothermic. It is shown that the local electric field of the protein environment serves as a template that stabilizes the intermediate of the H-abstraction step and thereby mediates the catalysis of dopamine formation at a lower energy cost. Furthermore, it is shown that external electric fields can either catalyze or inhibit the process depending on their directionality.
Collapse
Affiliation(s)
- Patric Schyman
- Institute of Chemistry and the Lise Meitner-Minerva Center for Computational Quantum Chemistry, The Hebrew University of Jerusalem, 91904 Jerusalem, Israel
| | | | | | | | | |
Collapse
|
47
|
Waszkowycz B, Clark DE, Gancia E. Outstanding challenges in protein–ligand docking and structure‐based virtual screening. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.18] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
| | - David E. Clark
- Argenta, 8/9 Spire Green Centre, Flex Meadow, Harlow CM19 5TR, UK
| | - Emanuela Gancia
- Argenta, 8/9 Spire Green Centre, Flex Meadow, Harlow CM19 5TR, UK
| |
Collapse
|
48
|
Conner KP, Woods C, Atkins WM. Interactions of cytochrome P450s with their ligands. Arch Biochem Biophys 2011; 507:56-65. [PMID: 20939998 PMCID: PMC3041843 DOI: 10.1016/j.abb.2010.10.006] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2010] [Revised: 10/01/2010] [Accepted: 10/04/2010] [Indexed: 01/12/2023]
Abstract
Cytochrome P450s (CYPs) are heme-containing monooxygenases that contribute to an enormous range of enzymatic function including biosynthetic and detoxification roles. This review summarizes recent studies concerning interactions of CYPs with ligands including substrates, inhibitors, and diatomic heme-ligating molecules. These studies highlight the complexity in the relationship between the heme spin state and active site occupancy, the roles of water in directing protein-ligand and ligand-heme interactions, and the details of interactions between heme and gaseous diatomic CYP ligands. Both kinetic and thermodynamic aspects of ligand binding are considered.
Collapse
Affiliation(s)
- Kip P. Conner
- Department of Medicinal Chemistry, Box 357610, University of Washington, Seattle, WA 98195-7610
| | - Caleb Woods
- Department of Medicinal Chemistry, Box 357610, University of Washington, Seattle, WA 98195-7610
| | - William M. Atkins
- Department of Medicinal Chemistry, Box 357610, University of Washington, Seattle, WA 98195-7610
| |
Collapse
|
49
|
Biarnés X, Bongarzone S, Vargiu AV, Carloni P, Ruggerone P. Molecular motions in drug design: the coming age of the metadynamics method. J Comput Aided Mol Des 2011; 25:395-402. [PMID: 21327922 DOI: 10.1007/s10822-011-9415-3] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2010] [Accepted: 01/28/2011] [Indexed: 01/25/2023]
Abstract
Metadynamics is emerging as a useful free energy method in physics, chemistry and biology. Recently, it has been applied also to investigate ligand binding to biomolecules of pharmacological interest. Here, after introducing the basic idea of the method, we review applications to challenging targets for pharmaceutical intervention. We show that this methodology, especially when combined with a variety of other computational approaches such as molecular docking and/or molecular dynamics simulation, may be useful to predict structure and energetics of ligand/target complexes even when the targets lack a deep binding cavity, such as DNA and proteins undergoing fibrillation in neurodegenerative diseases. Furthermore, the method allows investigating the routes of molecular recognition and the associated binding energy profiles, providing a molecular interpretation to experimental data.
Collapse
Affiliation(s)
- Xevi Biarnés
- International School for Advanced Studies, Trieste, Italy
| | | | | | | | | |
Collapse
|
50
|
Tarcsay Á, Keserű GM. In silicosite of metabolism prediction of cytochrome P450-mediated biotransformations. Expert Opin Drug Metab Toxicol 2011; 7:299-312. [DOI: 10.1517/17425255.2011.553599] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|