1
|
Dulay ANG, de Guzman JCC, Marquez ZYD, Santana ESD, Arce J, Orosco FL. The potential of Chlorella spp. as antiviral source against African swine fever virus through a virtual screening pipeline. J Mol Graph Model 2024; 132:108846. [PMID: 39151375 DOI: 10.1016/j.jmgm.2024.108846] [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: 04/01/2024] [Revised: 06/26/2024] [Accepted: 08/02/2024] [Indexed: 08/19/2024]
Abstract
African swine fever (ASF) causes high mortality in pigs and threatens global swine production. There is still a lack of therapeutics available, with two vaccines under scrutiny and no approved small-molecule drugs. Eleven (11) viral proteins were used to identify potential antivirals in in silico screening of secondary metabolites (127) from Chlorella spp. The metabolites were screened for affinity and binding selectivity. High-scoring compounds were assessed through in silico ADMET (Absorption, Distribution, Metabolism, Excretion, Toxicity) predictions, compared to structurally similar drugs, and checked for off-target docking with prepared swine receptors. Molecular dynamics (MD) simulations determined binding stability while binding energy was measured in Molecular Mechanics - Generalized Born Surface Area (MMGBSA) or Poisson-Boltzmann Surface Area (MMPBSA). Only six (6) compounds passed until MD analyses, of which five (5) were stable after 100 ns of MD runs. Of these five compounds, only three had binding affinities that were comparable to or stronger than controls. Specifically, phytosterols 24,25-dihydrolanosterol and CID 4206521 that interact with the RNA capping enzyme (pNP868R), and ergosterol which bound to the Erv-like thioreductase (pB119L). The compounds identified in this study can be used as a theoretical basis for in vitro screening to develop potent antiviral drugs against ASFV.
Collapse
Affiliation(s)
- Albert Neil G Dulay
- Virology and Vaccine Research Program, Industrial Technology Development Institute, Department of Science and Technology, Taguig, 1632, Philippines
| | - John Christian C de Guzman
- Virology and Vaccine Research Program, Industrial Technology Development Institute, Department of Science and Technology, Taguig, 1632, Philippines
| | - Zyra Ysha D Marquez
- Department of Biology, College of Arts and Sciences, University of the Philippines - Manila, Manila, 1000, Philippines
| | - Elisha Sofia D Santana
- Department of Biology, College of Arts and Sciences, University of the Philippines - Manila, Manila, 1000, Philippines
| | - Jessamine Arce
- Department of Biology, College of Arts and Sciences, University of the Philippines - Manila, Manila, 1000, Philippines
| | - Fredmoore L Orosco
- Virology and Vaccine Research Program, Industrial Technology Development Institute, Department of Science and Technology, Taguig, 1632, Philippines; Department of Biology, College of Arts and Sciences, University of the Philippines - Manila, Manila, 1000, Philippines; S&T Fellows Program, Department of Science and Technology, Taguig, 1632, Philippines.
| |
Collapse
|
2
|
Zhao H. Structural Basis of Conformational Dynamics in the PROTAC-Induced Protein Degradation. ChemMedChem 2024; 19:e202400171. [PMID: 38655701 DOI: 10.1002/cmdc.202400171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2024] [Revised: 04/19/2024] [Accepted: 04/23/2024] [Indexed: 04/26/2024]
Abstract
Pronounced conformational dynamics is unveiled upon analyzing multiple crystal structures of the same proteins recruited to the same E3 ligases by PROTACs, and yet, is largely permissive for targeted protein degradation due to the intrinsic mobility of E3 assemblies creating a large ubiquitylation zone. Mathematical modelling of ternary dynamics on ubiquitylation probability confirms the experimental finding that ternary complex rigidification need not correlate with enhanced protein degradation. Salt bridges are found to prevail in the PROTAC-induced ternary complexes, and may contribute to a positive cooperativity and prolonged half-life. The analysis highlights the importance of presenting lysines close to the active site of the E2 enzyme while constraining ternary dynamics in PROTAC design to achieve high degradation efficiency.
Collapse
Affiliation(s)
- Hongtao Zhao
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| |
Collapse
|
3
|
Ugrani S. Inhibitor design for TMPRSS2: insights from computational analysis of its backbone hydrogen bonds using a simple descriptor. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2024; 53:27-46. [PMID: 38157015 PMCID: PMC10853362 DOI: 10.1007/s00249-023-01695-4] [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: 09/05/2023] [Revised: 12/04/2023] [Accepted: 12/07/2023] [Indexed: 01/03/2024]
Abstract
Transmembrane protease serine 2 (TMPRSS2) is an important drug target due to its role in the infection mechanism of coronaviruses including SARS-CoV-2. Current understanding regarding the molecular mechanisms of known inhibitors and insights required for inhibitor design are limited. This study investigates the effect of inhibitor binding on the intramolecular backbone hydrogen bonds (BHBs) of TMPRSS2 using the concept of hydrogen bond wrapping, which is the phenomenon of stabilization of a hydrogen bond in a solvent environment as a result of being surrounded by non-polar groups. A molecular descriptor which quantifies the extent of wrapping around BHBs is introduced for this. First, virtual screening for TMPRSS2 inhibitors is performed by molecular docking using the program DOCK 6 with a Generalized Born surface area (GBSA) scoring function. The docking results are then analyzed using this descriptor and its relationship to the solvent-accessible surface area term ΔGsa of the GBSA score is demonstrated with machine learning regression and principal component analysis. The effect of binding of the inhibitors camostat, nafamostat, and 4-guanidinobenzoic acid (GBA) on the wrapping of important BHBs in TMPRSS2 is also studied using molecular dynamics. For BHBs with a large increase in wrapping groups due to these inhibitors, the radial distribution function of water revealed that certain residues involved in these BHBs, like Gln438, Asp440, and Ser441, undergo preferential desolvation. The findings offer valuable insights into the mechanisms of these inhibitors and may prove useful in the design of new inhibitors.
Collapse
Affiliation(s)
- Suraj Ugrani
- Purdue University, West Lafayette, IN, 47907, USA.
| |
Collapse
|
4
|
Nilewski C, Labadie S, Wei B, Malhotra S, Do S, Gazzard L, Liu L, Shao C, Murray J, Izrayelit Y, Gustafson A, Endres NF, Ma F, Ye X, Zou J, Evangelista M. Structure-Based Design and Evaluation of Reversible KRAS G13D Inhibitors. ACS Med Chem Lett 2024; 15:21-28. [PMID: 38229748 PMCID: PMC10788945 DOI: 10.1021/acsmedchemlett.3c00478] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Revised: 11/18/2023] [Accepted: 11/28/2023] [Indexed: 01/18/2024] Open
Abstract
Oncogenic KRAS mutations were identified decades ago, yet the selective inhibition of specific KRAS mutant proteins represents an ongoing challenge. Recent progress has been made in targeting certain P-loop mutant proteins, in particular KRAS G12C, for which the covalent inhibition of the GDP state via the Switch II pocket is now a clinically validated strategy. Inhibition of other KRAS mutant proteins such as KRAS G13D, on the other hand, still requires clinical validation. The remoteness of the D13 residue relative to the Switch II pocket in combination with the solvent exposure and conformational flexibility of the D13 side chain, as well as the difficulties of targeting carboxylate residues covalently, renders this specific protein particularly challenging to target selectively. In this report, we describe the design and evaluation of potent and KRAS G13D-selective reversible inhibitors. Subnanomolar binding to the GDP state Switch II pocket and biochemical selectivity over WT KRAS are achieved by leveraging a salt bridge with D13.
Collapse
Affiliation(s)
- Christian Nilewski
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Sharada Labadie
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Binqing Wei
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Sushant Malhotra
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Steven Do
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Lewis Gazzard
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Li Liu
- Pharmaron-Beijing
Co. Ltd., 6 Taihe Road, BDA, Beijing 100176, P. R. China
| | - Cheng Shao
- Pharmaron-Beijing
Co. Ltd., 6 Taihe Road, BDA, Beijing 100176, P. R. China
| | - Jeremy Murray
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Yevgeniy Izrayelit
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Amy Gustafson
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Nicholas F. Endres
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Fang Ma
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Xin Ye
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Jun Zou
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Marie Evangelista
- Genentech
Inc., 1 DNA Way, South San Francisco, California 94080, United States
| |
Collapse
|
5
|
Zhao H. Kinetic modelling of the P-glycoprotein mediated efflux with a large-scale matched molecular pair analysis. Eur J Med Chem 2023; 261:115830. [PMID: 37774507 DOI: 10.1016/j.ejmech.2023.115830] [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: 05/30/2023] [Revised: 09/19/2023] [Accepted: 09/21/2023] [Indexed: 10/01/2023]
Abstract
P-glycoprotein (Pgp) mediated efflux impacts on the drug absorption, distribution, metabolism and excretion, and confers multidrug resistance to cancer cells. Kinetic modelling provides mechanistic insights into the relationship between the substrate-Pgp interactions and efflux, and bridges the gap between the preference of polar compounds as Pgp substrates and the hydrophobic nature of its drug-binding site. Matched molecular pair analysis supports the guidelines of controlling H-bond donors and polar surface area in the efflux mitigation, but also reveals insufficiency of this type of rule-based approach. Contrary to the rule-of-five compliant compounds, proteolysis-targeting chimeras (PROTACs) have shown the opposite preference of physicochemical properties to evade efflux. Our analysis reiterates the critical role of intrinsic passive permeability in the efflux ratio, and indeed, its mitigation is often driven by increased passive permeability. It is thus useful to separate the passive permeability from the structural context-specific substrate-Pgp interactions in the design cycle.
Collapse
Affiliation(s)
- Hongtao Zhao
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden.
| |
Collapse
|
6
|
Kariuki M, Rho JY, Hall SCL, Perrier S. Investigating the Impact of Hydrophobic Polymer Segments on the Self-Assembly Behavior of Supramolecular Cyclic Peptide Systems via Asymmetric-Flow Field Flow Fractionation. Macromolecules 2023; 56:6618-6632. [PMID: 37720562 PMCID: PMC10501196 DOI: 10.1021/acs.macromol.3c00442] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Revised: 08/08/2023] [Indexed: 09/19/2023]
Abstract
The present study examines the behavior of cyclic peptide polymer conjugates that have been designed to combine their self-assembling ability via H-bonding with the properties of amphiphilic diblock copolymers. Using a combination of asymmetric flow-field flow fractionation (AF4) and small-angle neutron scattering (SANS), we have uncovered unique insight based on the population of structures established at a 24 h equilibrium profile. Our results determine that by introducing a small quantity of hydrophobicity into the conjugated polymer corona, the resulting nanotube structures exhibit low unimer dissociation which signifies enhanced stability. Furthermore, as the hydrophobicity of the polymer corona is increased, the elongation of the nanotubes is observed due to an increase in the association of unimers. This encompasses not only the H-bonding of unimers into nanotubes but also the self-assembly of single nanotubes into segmented-nanotube structures with high aspect ratios. However, this influence relies on a subtle balance between the hydrophobicity and hydrophilicity of the polymer corona. This balance is proposed to determine the solvent entropic penalty of hydrating the system, whereby the cost scales with the hydrophobic quantity. Consequently, it has been suggested that at a critical hydrophobic quantity, the solvation penalty becomes high enough such that the self-assembly of the system deviates from ordered hydrogen bonding. The association behavior is instead dominated by the hydrophobic effect which results in the undesirable formation of disordered aggregates.
Collapse
Affiliation(s)
- Maria Kariuki
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, U.K.
| | - Julia Y. Rho
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, U.K.
| | - Stephen C. L. Hall
- ISIS
Neutron and Muon Source, Rutherford Appleton
Laboratory, Didcot OX11 0QX, U.K.
| | - Sébastien Perrier
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, U.K.
- Warwick
Medical School, University of Warwick, Coventry CV4 7AL, U.K.
- Faculty
of Pharmacy and Pharmaceutical Sciences, Monash University, Parkville, VIC 3052, Australia
| |
Collapse
|
7
|
Bhuia MS, Rokonuzzman M, Hossain MI, Ansari SA, Ansari IA, Islam T, Al Hasan MS, Mubarak MS, Islam MT. Anxiolytic- like Effects by trans-Ferulic Acid Possibly Occur through GABAergic Interaction Pathways. Pharmaceuticals (Basel) 2023; 16:1271. [PMID: 37765079 PMCID: PMC10535412 DOI: 10.3390/ph16091271] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Revised: 09/02/2023] [Accepted: 09/05/2023] [Indexed: 09/29/2023] Open
Abstract
Numerous previous studies reported that ferulic acid exerts anxiolytic activity. However, the mechanisms have yet to be elucidated. The current study aimed to investigate the anxiolytic effect of trans-ferulic acid (TFA), a stereoisomer of ferulic acid, and evaluated its underlying mechanism using in vivo and computational studies. For this, different experimental doses of TFA (25, 50, and 75 mg/kg) were administered orally to Swiss albino mice, and various behavioral methods of open field, hole board, swing box, and light-dark tests were carried out. Diazepam (DZP), a positive allosteric modulator of the GABAA receptor, was employed as a positive control at a dose of 2 mg/kg, and distilled water served as a vehicle. Additionally, molecular docking was performed to estimate the binding affinities of the TFA and DZP toward the GABAA receptor subunits of α2 and α3, which are associated with the anxiolytic effect; visualizations of the ligand-receptor interaction were carried out using various computational tools. Our findings indicate that TFA dose-dependently reduces the locomotor activity of the animals in comparison with the controls, calming their behaviors. In addition, TFA exerted the highest binding affinity (-5.8 kcal/mol) to the α2 subunit of the GABAA receptor by forming several hydrogen and hydrophobic bonds. Taken together, our findings suggest that TFA exerts a similar effect to DZP, and the compound exerts moderate anxiolytic activity through the GABAergic interaction pathway. We suggest further clinical studies to develop TFA as a reliable anxiolytic agent.
Collapse
Affiliation(s)
- Md. Shimul Bhuia
- Department of Pharmacy, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj 8100, Bangladesh; (M.S.B.); (M.R.); (M.I.H.); (T.I.); (M.S.A.H.)
| | - Md. Rokonuzzman
- Department of Pharmacy, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj 8100, Bangladesh; (M.S.B.); (M.R.); (M.I.H.); (T.I.); (M.S.A.H.)
| | - Md. Imran Hossain
- Department of Pharmacy, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj 8100, Bangladesh; (M.S.B.); (M.R.); (M.I.H.); (T.I.); (M.S.A.H.)
| | - Siddique Akber Ansari
- Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, P.O. Box 2457, Riyadh 11451, Saudi Arabia;
| | - Irfan Aamer Ansari
- Department of Drug Science and Technology, University of Turin, 10124 Turin, Italy;
| | - Tawhida Islam
- Department of Pharmacy, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj 8100, Bangladesh; (M.S.B.); (M.R.); (M.I.H.); (T.I.); (M.S.A.H.)
| | - Md. Sakib Al Hasan
- Department of Pharmacy, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj 8100, Bangladesh; (M.S.B.); (M.R.); (M.I.H.); (T.I.); (M.S.A.H.)
| | | | - Muhammad Torequl Islam
- Department of Pharmacy, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj 8100, Bangladesh; (M.S.B.); (M.R.); (M.I.H.); (T.I.); (M.S.A.H.)
| |
Collapse
|
8
|
Mishra V, Pathak AK, Bandyopadhyay T. Binding of human serum albumin with uranyl ion at various pH: an all atom molecular dynamics study. J Biomol Struct Dyn 2023; 41:7318-7328. [PMID: 36099177 DOI: 10.1080/07391102.2022.2120080] [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: 07/28/2022] [Accepted: 08/26/2022] [Indexed: 10/14/2022]
Abstract
Uranium is routinely handled in various stages of nuclear fuel cycle and its association with human serum albumin (HSA) has been reported in literature, however, their binding characteristics still remains obscure. The present study aims to understand interaction of uranium with HSA by employing all atom molecular dynamics simulation of the HSA-metal ion complex. His67, His247 and Asp249 residues constitute the major binding site of HSA, which capture the uranyl ion (UO22+). A total of six sets of initial coordinates are used for Zn2+-HSA and UO22+-HSA system at pH = 4, 7.4 and 9, respectively. Enhance sampling method, namely, well-tempered meta-dynamics (WT-MtD) is employed to study the binding and un-binding processes of UO22+ and Zn2+ ions. Potential of mean force (PMF) profiles are generated for all the six sets of complexes from the converged WT-MtD run. Various basins and barriers are observed along the (un)binding pathways. Hydrogen bond dynamics and short-range Coulomb interactions are evaluated from the equilibrium run at each basins and barriers for both the ions at all pH values. The binding of UO22+ ion with HSA is the result of the dynamical balance between UO22+-HSA and UO22+-water short range Coulomb interactions. Zn2+ ion interact more strongly than UO22+ at all pH through short range Coulomb interactions. PMF values further concludes that UO22+ cannot associate to the Zn2+ bound HSA protein but can be captured by free HSA at all pH values i.e. endosomal, alkaline and physiological pH.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Vijayakriti Mishra
- Radiation Safety Systems Division, Bhabha Atomic Research Centre, Mumbai, India
- Homi Bhabha National Institute, Mumbai, India
| | - Arup Kumar Pathak
- Homi Bhabha National Institute, Mumbai, India
- Chemistry Division, Bhabha Atomic Research Centre, Mumbai, India
| | - Tusar Bandyopadhyay
- Homi Bhabha National Institute, Mumbai, India
- Chemistry Division, Bhabha Atomic Research Centre, Mumbai, India
| |
Collapse
|
9
|
Aloke C, Iwuchukwu EA, Achilonu I. Exploiting Copaifera salikounda compounds as treatment against diabetes: An insight into their potential targets from a computational perspective. Comput Biol Chem 2023; 104:107851. [DOI: 10.1016/j.compbiolchem.2023.107851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2022] [Revised: 02/25/2023] [Accepted: 03/19/2023] [Indexed: 03/29/2023]
|
10
|
Vázquez J, Ginex T, Herrero A, Morisseau C, Hammock BD, Luque FJ. Screening and Biological Evaluation of Soluble Epoxide Hydrolase Inhibitors: Assessing the Role of Hydrophobicity in the Pharmacophore-Guided Search of Novel Hits. J Chem Inf Model 2023; 63:3209-3225. [PMID: 37141492 PMCID: PMC10207366 DOI: 10.1021/acs.jcim.3c00301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2023] [Indexed: 05/06/2023]
Abstract
The human soluble epoxide hydrolase (sEH) is a bifunctional enzyme that modulates the levels of regulatory epoxy lipids. The hydrolase activity is carried out by a catalytic triad located at the center of a wide L-shaped binding site, which contains two hydrophobic subpockets at both sides. On the basis of these structural features, it can be assumed that desolvation is a major factor in determining the maximal achievable affinity that can be attained for this pocket. Accordingly, hydrophobic descriptors may be better suited to the search of novel hits targeting this enzyme. This study examines the suitability of quantum mechanically derived hydrophobic descriptors in the discovery of novel sEH inhibitors. To this end, three-dimensional quantitative structure-activity relationship (3D-QSAR) pharmacophores were generated by combining electrostatic and steric or alternatively hydrophobic and hydrogen-bond parameters in conjunction with a tailored list of 76 known sEH inhibitors. The pharmacophore models were then validated by using two external sets chosen (i) to rank the potency of four distinct series of compounds and (ii) to discriminate actives from decoys, using in both cases datasets taken from the literature. Finally, a prospective study was performed including a virtual screening of two chemical libraries to identify new potential hits, which were subsequently experimentally tested for their inhibitory activity on human, rat, and mouse sEH. The use of hydrophobic-based descriptors led to the identification of six compounds as inhibitors of the human enzyme with IC50 < 20 nM, including two with IC50 values of 0.4 and 0.7 nM. The results support the use of hydrophobic descriptors as a valuable tool in the search of novel scaffolds that encode a proper hydrophilic/hydrophobic distribution complementary to the target's binding site.
Collapse
Affiliation(s)
- Javier Vázquez
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia i Ciències de l′Alimentació, Institut de Biomedicina (IBUB), Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
- Pharmacelera,
Parc Científic de Barcelona (PCB), Baldiri Reixac 4-8, 08028 Barcelona, Spain
| | - Tiziana Ginex
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia i Ciències de l′Alimentació, Institut de Biomedicina (IBUB), Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
| | - Albert Herrero
- Pharmacelera,
Parc Científic de Barcelona (PCB), Baldiri Reixac 4-8, 08028 Barcelona, Spain
| | - Christophe Morisseau
- Department
of Entomology and Nematology, and Comprehensive Cancer Center, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Bruce D. Hammock
- Department
of Entomology and Nematology, and Comprehensive Cancer Center, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - F. Javier Luque
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia i Ciències de l′Alimentació, Institut de Biomecidina (IBUB) and Institut de Química
Teòrica i Computacional (IQTCUB), Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
| |
Collapse
|
11
|
Miñarro-Lleonar M, Bertran-Mostazo A, Duro J, Barril X, Juárez-Jiménez J. Lenalidomide Stabilizes Protein-Protein Complexes by Turning Labile Intermolecular H-Bonds into Robust Interactions. J Med Chem 2023; 66:6037-6046. [PMID: 37083375 PMCID: PMC10184122 DOI: 10.1021/acs.jmedchem.2c01692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/22/2023]
Abstract
Targeted protein degradation is a promising therapeutic strategy, spearheaded by the anti-myeloma drugs lenalidomide and pomalidomide. These drugs stabilize very efficiently the complex between the E3 ligase Cereblon (CRBN) and several non-native client proteins (neo-substrates), including the transcription factors Ikaros and Aiolos and the enzyme Caseine Kinase 1α (CK1α,), resulting in their degradation. Although the structures for these complexes have been determined, there are no evident interactions that can account for the high efficiency of formation of the ternary complex. We show that lenalidomide's stabilization of the CRBN-CK1α complex is largely due to hydrophobic shielding of intermolecular hydrogen bonds. We also find a quantitative relationship between hydrogen bond robustness and binding affinities of the ternary complexes. These results pave the way to further understand cooperativity effects in drug-induced protein-protein complexes and could help in the design of improved molecular glues and more efficient protein degraders.
Collapse
Affiliation(s)
- Marina Miñarro-Lleonar
- Unitat de Fisicoquímica, Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l'Alimentació, Universitat de Barcelona (UB), Av. Joan XXIII, 27-31, 08028 Barcelona, Spain
- Institut de Química Teòrica i Computacional (IQTC), Facultat de Química i Física, Universitat de Barcelona (UB), C. Martí i Franquès, 1, 08028 Barcelona, Spain
- Institut de Biomedicina, Facultat de Biologia, Universitat de Barcelona (UB), Av. Diagonal, 643, 08028 Barcelona, Spain
| | - Andrea Bertran-Mostazo
- Unitat de Fisicoquímica, Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l'Alimentació, Universitat de Barcelona (UB), Av. Joan XXIII, 27-31, 08028 Barcelona, Spain
- Institut de Biomedicina, Facultat de Biologia, Universitat de Barcelona (UB), Av. Diagonal, 643, 08028 Barcelona, Spain
| | - Jorge Duro
- Unitat de Fisicoquímica, Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l'Alimentació, Universitat de Barcelona (UB), Av. Joan XXIII, 27-31, 08028 Barcelona, Spain
- Institut de Química Teòrica i Computacional (IQTC), Facultat de Química i Física, Universitat de Barcelona (UB), C. Martí i Franquès, 1, 08028 Barcelona, Spain
| | - Xavier Barril
- Unitat de Fisicoquímica, Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l'Alimentació, Universitat de Barcelona (UB), Av. Joan XXIII, 27-31, 08028 Barcelona, Spain
- Institut de Química Teòrica i Computacional (IQTC), Facultat de Química i Física, Universitat de Barcelona (UB), C. Martí i Franquès, 1, 08028 Barcelona, Spain
- Institut de Biomedicina, Facultat de Biologia, Universitat de Barcelona (UB), Av. Diagonal, 643, 08028 Barcelona, Spain
- Catalan Institution for Research and Advanced Studies (ICREA), Pg. Lluís Companys, 23 08010 Barcelona, Spain
| | - Jordi Juárez-Jiménez
- Unitat de Fisicoquímica, Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l'Alimentació, Universitat de Barcelona (UB), Av. Joan XXIII, 27-31, 08028 Barcelona, Spain
- Institut de Química Teòrica i Computacional (IQTC), Facultat de Química i Física, Universitat de Barcelona (UB), C. Martí i Franquès, 1, 08028 Barcelona, Spain
| |
Collapse
|
12
|
Aziz S, Waqas M, Iqbal A, Halim SA, Abdellattif MH, Khan A, Al-Harrasi A. Structure-based identification of potential substrate antagonists for isethionate sulfite-lyase enzyme of Bilophila Wadsworthia: Towards novel therapeutic intervention to curb gut-associated illness. Int J Biol Macromol 2023; 240:124428. [PMID: 37062383 DOI: 10.1016/j.ijbiomac.2023.124428] [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: 02/28/2023] [Revised: 04/03/2023] [Accepted: 04/09/2023] [Indexed: 04/18/2023]
Abstract
Bilophila wadsworthia is one of the prominent sources of hydrogen sulfide (H2S) production in appendices, excessive levels of which can result in a weaker colonic mucus barrier, inflammatory bowel disease, and colorectal cancer. Isethionate sulfite-lyase (IslA) enzyme catalyzes H2S production by cleaving CS bond in isethionate, producing acetaldehyde and sulfite. In this study, we aimed to identify potential substrate antagonists for IsIA using a structure-based drug design. Initially, pharmacophore-based computational screening of the ZINC20 database yielded 66 hits that were subjected to molecular docking targeting the isethionate binding site of IsIA. Based on striking docking scores, nine compounds showed strong interaction with critical IsIA residues (Arg189, Gln193, Glu470, Cys468, and Arg678), drug-like features, appropriate adsorption, metabolism, excretion, and excretion profile with non-toxicity. Molecular dynamics simulations uncovered the significant impact of binding the compounds on protein conformational dynamics. Finally, binding free energies revealed substantial binding affinity (ranging from -35.23 to -53.88 kcal/mol) of compounds (ZINC913876497, ZINC913856647, ZINC914263733, ZINC914137795, ZINC915757996, ZINC914357083, ZINC913934833, ZINC9143362047, and ZINC913854740) for IsIA. The compounds proposed herein through a multi-faceted computational strategy can be experimentally validated as potential substrate antagonists of B. wadsworthia's IsIA for developing new medications to curb gut-associated illness in the future.
Collapse
Affiliation(s)
- Shahkaar Aziz
- Institute of Biotechnology and Genetic Engineering, The University of Agriculture, Peshawar, Pakistan
| | - Muhammad Waqas
- Department of Biotechnology and Genetic Engineering, Hazara University, Mansehra, Pakistan; Natural and Medical Sciences Research Center, University of Nizwa, Birkat-ul-Mouz, Nizwa, Oman
| | - Aqib Iqbal
- Institute of Biotechnology and Genetic Engineering, The University of Agriculture, Peshawar, Pakistan; Department of Biotechnology, Abdul Wali Khan University Mardan, Mardan, Pakistan.
| | - Sobia Ahsan Halim
- Natural and Medical Sciences Research Center, University of Nizwa, Birkat-ul-Mouz, Nizwa, Oman
| | - Magda H Abdellattif
- Department of Chemistry, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
| | - Ajmal Khan
- Natural and Medical Sciences Research Center, University of Nizwa, Birkat-ul-Mouz, Nizwa, Oman.
| | - Ahmed Al-Harrasi
- Natural and Medical Sciences Research Center, University of Nizwa, Birkat-ul-Mouz, Nizwa, Oman.
| |
Collapse
|
13
|
Integrated approach for studying bioactive compounds from Cladosporium spp. against estrogen receptor alpha as breast cancer drug target. Sci Rep 2022; 12:22446. [PMID: 36575224 PMCID: PMC9794773 DOI: 10.1038/s41598-022-22038-x] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2022] [Accepted: 10/07/2022] [Indexed: 12/28/2022] Open
Abstract
Cladosporium spp. have been reported for their great diversity of secondary metabolites which represent as a prominent base material for verifying the biological activities. Several bioactive compounds which have antimicrobial, cytotoxic, quorum sensing inhibitory and phytotoxic activities have been isolated from Cladosporium species. Most of them are still needed to be explored for their anticancer properties. Therefore, the present study is focused on screening and identifying the bioactive compounds of Cladosporium spp. for their anticancer activity via the integrated approaches of Molecular Docking (MD), Molecular Dynamics Simulation (MDS) and Density Functional Theory (DFT) studies. A total of 123 bioactive compounds of Cladosporium spp. were explored for their binding affinity with the selected breast cancer drug target receptor such as estrogen receptor alpha (PDB:6CBZ). The Molecular Docking studies revealed that amongst the bioactive compounds screened, Altertoxin X and Cladosporol H showed a good binding affinity of - 10.5 kcal/mol and - 10.3 kcal/mol, respectively, with the estrogen receptor alpha when compared to the reference compound (17[Formula: see text]-Estradiol: - 10.2 kcal/mol). The MDS study indicated the stable binding patterns and conformation of the estrogen receptor alpha-Altertoxin X complex in a stimulating environment. In addition, in silico absorption, distribution, metabolism, excretion and toxicity (ADMET) study suggested that Altertoxin X has a good oral bioavailability with a high LD[Formula: see text] value of 2.375 mol/kg and did not cause any hepatotoxicity and skin sensitization. In summary, the integrated approaches revealed that Altertoxin X possesses a promising anticancer activity and could serve as a new therapeutic drug for breast cancer treatment.
Collapse
|
14
|
Bekker GJ, Kamiya N. Advancing the field of computational drug design using multicanonical molecular dynamics-based dynamic docking. Biophys Rev 2022; 14:1349-1358. [PMID: 36659995 PMCID: PMC9842809 DOI: 10.1007/s12551-022-01010-z] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Accepted: 10/14/2022] [Indexed: 11/20/2022] Open
Abstract
Multicanonical molecular dynamics (McMD)-based dynamic docking is a powerful tool to not only predict the native binding configuration between two flexible molecules, but it can also be used to accurately simulate the binding/unbinding pathway. Furthermore, it can also predict alternative binding sites, including allosteric ones, by employing an exhaustive sampling approach. Since McMD-based dynamic docking accurately samples binding/unbinding events, it can thus be used to determine the molecular mechanism of binding between two molecules. We developed the McMD-based dynamic docking methodology based on the powerful, but woefully underutilized McMD algorithm, combined with a toolset to perform the docking and to analyze the results. Here, we showcase three of our recent works, where we have applied McMD-based dynamic docking to advance the field of computational drug design. In the first case, we applied our method to perform an exhaustive search between Hsp90 and one of its inhibitors to successfully predict the native binding configuration in its binding site, as we refined our analysis methods. For our second case, we performed an exhaustive search of two medium-sized ligands and Bcl-xL, which has a cryptic binding site that differs greatly between the apo and holo structures. Finally, we performed a dynamic docking simulation between a membrane-embedded GPCR molecule and a high affinity ligand that binds deep within its receptor's pocket. These advanced simulations showcase the power that the McMD-based dynamic docking method has, and provide a glimpse of the potential our methodology has to unravel and solve the medical and biophysical issues in the modern world. Supplementary Information The online version contains supplementary material available at 10.1007/s12551-022-01010-z.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871 Japan
| | - Narutoshi Kamiya
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047 Japan
| |
Collapse
|
15
|
Ziada S, Diharce J, Raimbaud E, Aci-Sèche S, Ducrot P, Bonnet P. Estimation of Drug-Target Residence Time by Targeted Molecular Dynamics Simulations. J Chem Inf Model 2022; 62:5536-5549. [PMID: 36350238 DOI: 10.1021/acs.jcim.2c00852] [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/11/2022]
Abstract
Drug-target residence time has emerged as a key selection factor in drug discovery since the binding duration of a drug molecule to its protein target can significantly impact its in vivo efficacy. The challenge in studying the residence time, in early drug discovery stages, lies in how to cost-effectively determine the residence time for the systematic assessment of compounds. Currently, there is still a lack of computational protocols to quickly estimate such a measure, particularly for large and flexible protein targets and drugs. Here, we report an efficient computational protocol, based on targeted molecular dynamics, to rank drug candidates by their residence time and to obtain insights into ligand-target dissociation mechanisms. The method was assessed on a dataset of 10 arylpyrazole inhibitors of CDK8, a large, flexible, and clinically important target, for which the experimental residence time of the inhibitors ranges from minutes to hours. The compounds were correctly ranked according to their estimated residence time scores compared to their experimental values. The analysis of protein-ligand interactions along the dissociation trajectories highlighted the favorable contribution of hydrophobic contacts to residence time and revealed key residues that strongly affect compound residence time.
Collapse
Affiliation(s)
- Sonia Ziada
- Institut de Chimie Organique et Analytique (ICOA), UMR CNRS-Université d'Orléans 7311, Université d'Orléans BP 6759, Orléans Cedex 245067, France
| | - Julien Diharce
- Institut de Chimie Organique et Analytique (ICOA), UMR CNRS-Université d'Orléans 7311, Université d'Orléans BP 6759, Orléans Cedex 245067, France
| | - Eric Raimbaud
- Institut de Recherches Servier, 125 Chemin de Ronde, Croissy-sur-Seine78290, France
| | - Samia Aci-Sèche
- Institut de Chimie Organique et Analytique (ICOA), UMR CNRS-Université d'Orléans 7311, Université d'Orléans BP 6759, Orléans Cedex 245067, France
| | - Pierre Ducrot
- Institut de Recherches Servier, 125 Chemin de Ronde, Croissy-sur-Seine78290, France
| | - Pascal Bonnet
- Institut de Chimie Organique et Analytique (ICOA), UMR CNRS-Université d'Orléans 7311, Université d'Orléans BP 6759, Orléans Cedex 245067, France
| |
Collapse
|
16
|
Tateing S, Suree N. Decoding molecular recognition of inhibitors targeting HDAC2 via molecular dynamics simulations and configurational entropy estimation. PLoS One 2022; 17:e0273265. [PMID: 35981056 PMCID: PMC9387782 DOI: 10.1371/journal.pone.0273265] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2022] [Accepted: 08/05/2022] [Indexed: 11/18/2022] Open
Abstract
Molecular recognition by enzymes is a complicated process involving thermodynamic energies governing protein-ligand interactions. In order to aid the estimation of inhibitory activity of compounds targeting an enzyme, several computational methods can be employed to dissect this intermolecular contact. Herein, we report a structural dynamics investigation of an epigenetic enzyme HDAC2 in differentiating its binding to various inhibitors within the sub-sites of its active site. Molecular dynamics (MD) simulation was employed to elucidate the intermolecular interactions as well as the dynamics behavior of ligand binding. MD trajectories of five distinct HDAC2-inhibitor complexes reveal that compounds lacking adequate contacts with the opening rim of the active site possess high fluctuation along the cap portion, thus weakening the overall affinity. Key intermolecular interactions determining the effective binding of inhibitors include hydrogen bonds with Gly154, Asp181, and Tyr308; hydrophobic interactions between Phe155/Phe210 and the linker region; and a pi-stacking with Arg39 at the foot pocket. Decomposition of the binding free energy calculated per-residue by MM/PBSA also indicates that the interactions within the internal foot pocket, especially with residues Met35, Leu144, Gly305, and Gly306, can contribute significantly to the ligand binding. Additionally, configurational entropy of the binding was estimated and compared to the scale of the binding free energy in order to assess its contribution to the binding and to differentiate various ligand partners. It was found that the levels of entropic contribution are comparable among a set of structurally similar carbamide ligands, while it is greatly different for the set of unrelated ligands, ranging from 2.75 to 16.38 kcal/mol for the five inhibitors examined. These findings exemplify the importance of assessing molecular dynamics as well as estimating the entropic contribution in evaluating the ligand binding mechanism.
Collapse
Affiliation(s)
- Suriya Tateing
- Interdisciplinary Program in Biotechnology, Graduate School, Chiang Mai University, Chiang Mai, Thailand
- Division of Biochemistry and Biochemical Innovation, Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand
- Department of Plant and Soil Sciences, Faculty of Agriculture, Chiang Mai University, Chiang Mai, Thailand
| | - Nuttee Suree
- Division of Biochemistry and Biochemical Innovation, Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand
- Center of Excellence in Materials Science and Technology, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand
- * E-mail:
| |
Collapse
|
17
|
Zhuang Y, Thota N, Quirk S, Hernandez R. Implementation of Telescoping Boxes in Adaptive Steered Molecular Dynamics. J Chem Theory Comput 2022; 18:4649-4659. [PMID: 35830368 DOI: 10.1021/acs.jctc.2c00498] [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/28/2022]
Abstract
Long-time dynamical processes, such as those involving protein unfolding and ligand interactions, can be accelerated and realized through steered molecular dynamics (SMD). The challenge has been the extraction of information from such simulations that generalize for complex nonequilibrium processes. The use of Jarzynski's equality opened the possibility of determining the free energy along the steered coordinate, but sampling over the nonequilibrium trajectories is slow to converge. Adaptive steered molecular dynamics (ASMD) and other related techniques have been introduced to overcome this challenge through the use of stages. Here, we take advantage of these stages to address the numerical cost that arises from the required use of very large solvent boxes. We introduce telescoping box schemes within adaptive steered molecular dynamics (ASMD) in which we adjust the solvent box between stages and thereby vary (and optimize) the required number of solvent molecules. We have benchmarked the method on a relatively long α-helical peptide, Ala30, with respect to the potential of mean force and hydrogen bonds. We show that the use of telescoping boxes introduces little numerical error while significantly reducing the computational cost.
Collapse
Affiliation(s)
- Yi Zhuang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Nikhil Thota
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Stephen Quirk
- Kimberly-Clark Corporation, Atlanta, Georgia 30076-2199, United States
| | - Rigoberto Hernandez
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States.,Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States.,Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
18
|
Ahmad K, Rizzi A, Capelli R, Mandelli D, Lyu W, Carloni P. Enhanced-Sampling Simulations for the Estimation of Ligand Binding Kinetics: Current Status and Perspective. Front Mol Biosci 2022; 9:899805. [PMID: 35755817 PMCID: PMC9216551 DOI: 10.3389/fmolb.2022.899805] [Citation(s) in RCA: 30] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2022] [Accepted: 05/09/2022] [Indexed: 12/12/2022] Open
Abstract
The dissociation rate (k off) associated with ligand unbinding events from proteins is a parameter of fundamental importance in drug design. Here we review recent major advancements in molecular simulation methodologies for the prediction of k off. Next, we discuss the impact of the potential energy function models on the accuracy of calculated k off values. Finally, we provide a perspective from high-performance computing and machine learning which might help improve such predictions.
Collapse
Affiliation(s)
- Katya Ahmad
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
| | - Andrea Rizzi
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
- Atomistic Simulations, Istituto Italiano di Tecnologia, Genova, Italy
| | - Riccardo Capelli
- Department of Applied Science and Technology (DISAT), Politecnico di Torino, Torino, Italy
| | - Davide Mandelli
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
| | - Wenping Lyu
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen, China
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, China
| | - Paolo Carloni
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
- Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich, Jülich, Germany
| |
Collapse
|
19
|
Ferenczy GG, Kellermayer M. Contribution of Hydrophobic Interactions to Protein Mechanical Stability. Comput Struct Biotechnol J 2022; 20:1946-1956. [PMID: 35521554 PMCID: PMC9062142 DOI: 10.1016/j.csbj.2022.04.025] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2021] [Revised: 04/07/2022] [Accepted: 04/17/2022] [Indexed: 11/26/2022] Open
Abstract
The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.
Collapse
|
20
|
Debnath J, Parrinello M. Computing Rates and Understanding Unbinding Mechanisms in Host-Guest Systems. J Chem Theory Comput 2022; 18:1314-1319. [PMID: 35200023 DOI: 10.1021/acs.jctc.1c01075] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The long time scale associated with ligand residence times renders their computation challenging. Therefore, the influence of factors like solvation and steric hindrance on residence times is not fully understood. Here, we demonstrate in a set of model host-guest systems that the recently developed Gaussian mixture based enhanced sampling allows residence times to be computed and enables an understanding of their unbinding mechanism. We observe that guest unbinding often proceeds via a series of intermediate states that can be labeled by the number of water molecules present in the binding cavity. In several cases the residence time is correlated to the water trapping times in the cavity.
Collapse
Affiliation(s)
- Jayashrita Debnath
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland.,Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| | | |
Collapse
|
21
|
Trisciuzzi D, Siragusa L, Baroni M, Autiero I, Nicolotti O, Cruciani G. Getting Insights into Structural and Energetic Properties of Reciprocal Peptide-Protein Interactions. J Chem Inf Model 2022; 62:1113-1125. [PMID: 35148095 DOI: 10.1021/acs.jcim.1c01343] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Peptide-protein interactions play a key role for many cellular and metabolic processes involved in the onset of largely spread diseases such as cancer and neurodegenerative pathologies. Despite the progress in the structural characterization of peptide-protein interfaces, the in-depth knowledge of the molecular details behind their interactions is still a daunting task. Here, we present the first comprehensive in silico morphological and energetic study of peptide binding sites by focusing on both peptide and protein standpoints. Starting from the PixelDB database, a nonredundant benchmark collection of high-quality 3D crystallographic structures of peptide-protein complexes, a classification analysis of the most representative categories based on the nature of each cocrystallized peptide has been carried out. Several interpretable geometrical and energetic descriptors have been computed both from peptide and target protein sides in the attempt to unveil physicochemical and structural causative correlations. Finally, we investigated the most frequent peptide-protein residue pairs at the binding interface and made extensive energetic analyses, based on GRID MIFs, with the aim to study the peptide affinity-enhancing interactions to be further exploited in rational drug design strategies.
Collapse
Affiliation(s)
- Daniela Trisciuzzi
- Department of Pharmacy, Pharmaceutical Sciences, Università degli Studi di Bari "Aldo Moro", 70125 Bari, Italy.,Molecular Horizon s.r.l., Via Montelino, 30, 06084 Bettona (PG), Italy
| | - Lydia Siragusa
- Molecular Horizon s.r.l., Via Montelino, 30, 06084 Bettona (PG), Italy.,Molecular Discovery Ltd., Kinetic Business Centre, Theobald Street, Elstree, Borehamwood, Hertfordshire WD6 4PJ, United Kingdom
| | - Massimo Baroni
- Molecular Discovery Ltd., Kinetic Business Centre, Theobald Street, Elstree, Borehamwood, Hertfordshire WD6 4PJ, United Kingdom
| | - Ida Autiero
- Molecular Horizon s.r.l., Via Montelino, 30, 06084 Bettona (PG), Italy.,National Research Council, Institute of Biostructures and Bioimaging, 80138 Naples, Italy
| | - Orazio Nicolotti
- Department of Pharmacy, Pharmaceutical Sciences, Università degli Studi di Bari "Aldo Moro", 70125 Bari, Italy
| | - Gabriele Cruciani
- Department of Chemistry, Biology and Biotechnology, Università degli Studi di Perugia, via Elce di Sotto, 8, 06123 Perugia (PG), Italy
| |
Collapse
|
22
|
Weerakoon D, Carbajo RJ, De Maria L, Tyrchan C, Zhao H. Impact of PROTAC Linker Plasticity on the Solution Conformations and Dissociation of the Ternary Complex. J Chem Inf Model 2022; 62:340-349. [PMID: 35018781 DOI: 10.1021/acs.jcim.1c01036] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The conformational behavior of a small molecule free in solution is important to understand the free energy of binding to its target. This could be of special interest for proteolysis-targeting chimeras (PROTACs) due to their often flexible and lengthy linkers and the need to induce a ternary complex. Here, we report on the molecular dynamics (MD) simulations of two PROTACs, MZ1 and dBET6, revealing different linker conformational behaviors. The simulation of MZ1 in dimethyl sulfoxide (DMSO) agrees well with the nuclear magnetic resonance study, providing strong support for the relevance of our simulations. To further understand the role of linker plasticity in the formation of a ternary complex, the dissociation of the complex von Hippel-Lindau-MZ1-BRD4 is investigated in detail by steered simulations and is shown to follow a two-step pathway. Interestingly, both MZ1 and dBET6 display in water, a tendency toward an intramolecular lipophilic interaction between the two warheads. The hydrophobic contact of the two warheads would prevent them from binding to their respective proteins and might have an effect on the efficacy of induced cellular protein degradation. However, conformations featuring this hydrophobic contact of the two warheads are calculated to be marginally more favorable.
Collapse
Affiliation(s)
- Dhanushka Weerakoon
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg 43183, Sweden
| | - Rodrigo J Carbajo
- Chemistry, Oncology R&D, AstraZeneca, Cambridge CB4 0QA, United Kingdom
| | - Leonardo De Maria
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg 43183, Sweden
| | - Christian Tyrchan
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg 43183, Sweden
| | - Hongtao Zhao
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg 43183, Sweden
| |
Collapse
|
23
|
Potterton A, Heifetz A, Townsend-Nicholson A. Predicting Residence Time of GPCR Ligands with Machine Learning. Methods Mol Biol 2022; 2390:191-205. [PMID: 34731470 DOI: 10.1007/978-1-0716-1787-8_8] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Drug-target residence time, the duration of binding at a given protein target, has been shown in some protein families to be more significant for conferring efficacy than binding affinity. To carry out efficient optimization of residence time in drug discovery, machine learning models that can predict that value need to be developed. One of the main challenges with predicting residence time is the paucity of data. This chapter outlines all of the currently available ligand kinetic data, providing a repository that contains the largest publicly available source of GPCR-ligand kinetic data to date. To help decipher the features of kinetic data that might be beneficial to include in computational models for the prediction of residence time, the experimental evidence for properties that influence residence time are summarized. Finally, two different workflows for predicting residence time with machine learning are outlined. The first is a single-target model trained on ligand features; the second is a multi-target model trained on features generated from molecular dynamics simulations.
Collapse
Affiliation(s)
- Andrew Potterton
- Structural and Molecular Biology, University College London, London, UK
- Evotec (U.K.) Ltd., Abingdon, Oxfordshire, UK
| | | | | |
Collapse
|
24
|
Kleijwegt RJT, Winkenwerder W, Baan W, van der Schaaf J. Degradation kinetics and solvent effects of various long‐chain quaternary ammonium salts. INT J CHEM KINET 2021. [DOI: 10.1002/kin.21537] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Roel J. T. Kleijwegt
- Laboratory of Chemical Reactor Engineering Department of Chemical Engineering and Chemistry Eindhoven University of Technology Eindhoven The Netherlands
| | | | - Wim Baan
- Nouryon, Zutphenseweg 10 Deventer The Netherlands
| | - John van der Schaaf
- Laboratory of Chemical Reactor Engineering Department of Chemical Engineering and Chemistry Eindhoven University of Technology Eindhoven The Netherlands
| |
Collapse
|
25
|
Takedomi K, Ohta M, Ekimoto T, Ikeguchi M. Effect of Water Molecules on the Activating S810L Mutation of the Mineralocorticoid Receptor. J Chem Inf Model 2021; 61:3583-3592. [PMID: 34228431 DOI: 10.1021/acs.jcim.1c00389] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The mineralocorticoid receptor (MR) is a nuclear receptor whose endogenous ligands are mineralocorticoids, a type of steroid hormone. The activating S810L mutation is known to cause severe early-onset and pregnancy-related hypertension. Progesterone binds to the wild-type (WT) MR as a passive antagonist with fast dissociation; however, it binds to the S810L mutant as a full agonist with slow dissociation. The switch in the biological activity of progesterone is considered to be one of the causes of the disease. First, we used steered molecular dynamics simulations to analyze the dissociation process of progesterone for the WT and the S810L mutant. Progesterone in the WT dissociated from the ligand-binding pocket with a weak force in comparison with progesterone in the S810L mutant due to the large inflow of water molecules into the pocket. Therefore, we used conventional molecular dynamics simulations for the ligand-free structures of the WT and the S810L mutant to investigate the effect of the mutation on the inflow of water. In the WT, water molecules enter the ligand-binding pocket in two ways: in the vicinity of (i) Arg817 and (ii) Ser810. In contrast, few water molecules enter the pocket in the S810L mutant because of the large size and hydrophobic nature of the Leu810 side chain. Fast dissociation is a common feature among passive antagonists of MR; therefore, we inferred that the water inflow could be responsible for the dissociation kinetics of progesterone in the WT and the S810L mutant.
Collapse
Affiliation(s)
- Kei Takedomi
- Graduate School of Medicinal Life Science, Yokohama City University, Yokohama 230-0045, Japan.,Modality Laboratories, Innovative Research Division, Mitsubishi Tanabe Pharma Corporation, 1000 Kamoshida, Aoba-ku, Yokohama 227-0033, Japan
| | - Masateru Ohta
- HPC- and AI-driven Drug Development Platform Division, Center for Computational Science, RIKEN, Yokohama 230-0045, Japan
| | - Toru Ekimoto
- Graduate School of Medicinal Life Science, Yokohama City University, Yokohama 230-0045, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medicinal Life Science, Yokohama City University, Yokohama 230-0045, Japan.,HPC- and AI-driven Drug Development Platform Division, Center for Computational Science, RIKEN, Yokohama 230-0045, Japan
| |
Collapse
|
26
|
Wang Z, Bosma R, Kuhne S, van den Bor J, Garabitian W, Vischer HF, Wijtmans M, Leurs R, de Esch IJ. Exploring the Effect of Cyclization of Histamine H 1 Receptor Antagonists on Ligand Binding Kinetics. ACS OMEGA 2021; 6:12755-12768. [PMID: 34056427 PMCID: PMC8154229 DOI: 10.1021/acsomega.0c06358] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/28/2021] [Accepted: 04/16/2021] [Indexed: 06/12/2023]
Abstract
There is an increasing interest in guiding hit optimization by considering the target binding kinetics of ligands. However, compared to conventional structure-activity relationships, structure-kinetics relationships have not been as thoroughly explored, even for well-studied archetypical drug targets such as the histamine H1 receptor (H1R), a member of the family A G-protein coupled receptor. In this study, we show that the binding kinetics of H1R antagonists at the H1R is dependent on the cyclicity of both the aromatic head group and the amine moiety of H1R ligands, the chemotypes that are characteristic for the first-generation H1R antagonists. Fusing the two aromatic rings of H1R ligands into one tricyclic aromatic head group prolongs the H1R residence time for benchmark H1R ligands as well as for tailored synthetic analogues. The effect of constraining the aromatic rings and the basic amines is systematically explored, leading to a coherent series and detailed discussions of structure-kinetics relationships. This study shows that cyclicity has a pronounced effect on the binding kinetics.
Collapse
Affiliation(s)
| | | | | | - Jelle van den Bor
- Amsterdam Institute of Molecular
and Life Sciences (AIMMS), Division of Medicinal Chemistry, Faculty
of Science, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Wrej Garabitian
- Amsterdam Institute of Molecular
and Life Sciences (AIMMS), Division of Medicinal Chemistry, Faculty
of Science, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Henry F. Vischer
- Amsterdam Institute of Molecular
and Life Sciences (AIMMS), Division of Medicinal Chemistry, Faculty
of Science, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Maikel Wijtmans
- Amsterdam Institute of Molecular
and Life Sciences (AIMMS), Division of Medicinal Chemistry, Faculty
of Science, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Rob Leurs
- Amsterdam Institute of Molecular
and Life Sciences (AIMMS), Division of Medicinal Chemistry, Faculty
of Science, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Iwan J.P. de Esch
- Amsterdam Institute of Molecular
and Life Sciences (AIMMS), Division of Medicinal Chemistry, Faculty
of Science, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| |
Collapse
|
27
|
Lim KS, Lee DY, Han S, Bull DA, Won YW. Targeted delivery of heat shock protein 90 inhibitors prevents growth of HER2-positive tumor. Biomaterials 2021; 273:120817. [PMID: 33894402 DOI: 10.1016/j.biomaterials.2021.120817] [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: 11/17/2020] [Revised: 03/29/2021] [Accepted: 04/06/2021] [Indexed: 12/30/2022]
Abstract
Heat shock protein 90 (HSP90) plays a crucial role in the survival of cancer cells. When an inhibitor blocks the signaling pathway of HSP90, its client proteins are degraded, destabilized, and inactivated. Although HSP90 inhibitors are in various clinical trials, there are no HSP90 inhibitor-immunoconjugates due to the difficulty in chemical modification of HSP90 inhibitors. Here we show that biological affinity binding enables the incorporation of HSP90 inhibitors to an antibody without the need for chemical conjugation. We constructed a recombinant fusion protein composed of an anti-HER2 scFv and an HSP90 inhibitor-binding domain (HER2 scFv-HBD). The HBD spontaneously captures a HSP90 inhibitor, resulting in the formation of an HER2 scFv-HBD/HSP90 inhibitor complex. In an HER2-positive cancer mouse model, targeted delivery of HSP90 inhibitors was confirmed and improved anti-cancer efficacy was observed. We have proven the promise of tumor-directed HSP90 inhibition as a new form of targeted therapy.
Collapse
Affiliation(s)
- Kwang Suk Lim
- Department of Biotechnology and Bioengineering, Interdisciplinary Program in Biohealth-machinery Convergence Engineering, College of Art, Culture and Engineering, Kangwon National University, South Korea
| | - Daniel Y Lee
- Division of Cardiothoracic Surgery, Department of Surgery, University of Arizona College of Medicine, Tucson, AZ, USA
| | - Seungmin Han
- Division of Cardiothoracic Surgery, Department of Surgery, University of Arizona College of Medicine, Tucson, AZ, USA
| | - David A Bull
- Division of Cardiothoracic Surgery, Department of Surgery, University of Arizona College of Medicine, Tucson, AZ, USA
| | - Young-Wook Won
- Division of Cardiothoracic Surgery, Department of Surgery, University of Arizona College of Medicine, Tucson, AZ, USA.
| |
Collapse
|
28
|
Liu W, Jin H, Chen T, Zhang G, Lai S, Liu G. Investigating the Role of the N-Terminal Loop of PD-1 in Binding Process Between PD-1 and Nivolumab via Molecular Dynamics Simulation. Front Mol Biosci 2020; 7:574759. [PMID: 33102523 PMCID: PMC7522605 DOI: 10.3389/fmolb.2020.574759] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2020] [Accepted: 08/14/2020] [Indexed: 12/16/2022] Open
Abstract
The blockade of immune checkpoints, such as programmed death receptor 1 (PD-1) and programmed death ligand 1 protein (PD-L1), is a promising therapeutic approach in cancer immunotherapy. Nivolumab, a humanized IgG4 antibody targeting PD-1, was approved by the US Food and Drug Administration for several cancers in 2014. Crystal structures of the nivolumab/PD-1 complex show that the epitope of PD-1 locates at the IgV domain (including the FG and BC loops) and the N-terminal loop. Although the N-terminal loop of PD-1 has been shown to play a dominant role in the complex interface of the static structure, its role in the dynamic binding process has not been illustrated clearly. Here, eight molecular systems were established for nivolumab/PD-1 complex, and long-time molecular dynamics simulations were performed for each. Results showed that the N-terminal loop of PD-1 prefers to bind with nivolumab to stabilize the interface between IgV and nivolumab. Furthermore, the binding of the N-terminal loop with nivolumab induces the rebinding between the IgV domain and nivolumab. Thus, we proposed a two-step binding model for the nivolumab/PD-1 binding, where the interface switches to a high-affinity state with the help of the N-terminal loop. This finding suggests that the N-terminal loop of PD-1 might be a potential target for anti-PD-1 antibody design, which could serve as an important gatekeeper for the anti-PD-1 antibody binding.
Collapse
Affiliation(s)
- Wenping Liu
- Guangdong Food and Drug Vocational College, Guangzhou, China
| | - Haoyu Jin
- Guangdong Food and Drug Vocational College, Guangzhou, China
| | - Ting Chen
- Guangdong Food and Drug Vocational College, Guangzhou, China
| | - Gangping Zhang
- Guangdong Food and Drug Vocational College, Guangzhou, China
| | - Shengsheng Lai
- Guangdong Food and Drug Vocational College, Guangzhou, China
| | - Guangjian Liu
- Clinical Data Center, Guangzhou Women and Children's Medical Center, Guangzhou Medical University, Guangzhou, China
| |
Collapse
|
29
|
Mishra L, Sundararajan M, Bandyopadhyay T. MD simulation reveals differential binding of Cm(III) and Th(IV) with serum transferrin at acidic pH. Proteins 2020; 89:193-206. [PMID: 32892408 DOI: 10.1002/prot.26006] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Revised: 07/13/2020] [Accepted: 08/27/2020] [Indexed: 01/21/2023]
Abstract
The iron carrier human serum transferrin (sTf) is known to transport other metals, including some actinides (An). Radiotoxic An are routinely involved in the nuclear fuel cycle and the possibility of their accidental exposure cannot be ruled out. Understanding An interaction with sTf assumes a greater significance for the development of safe and efficacious chelators for their removal from the blood stream. Here we report several 100 ns equilibrium MD simulations of Cm(III)- and Th(IV)-loaded sTf at various protonation states of the protein to explore the possibility of the two An ions release and speciation. The results demonstrate variation in protonation state of dilysine pair (K206 and K296) and the tyrosine (Y188) residue is necessary for the opening of Cm(III)-bound protein and the release of the ion. For the tetravalent thorium, protonation of dilysine pair suffices to cause conformational changes of protein. However, in none of the protonation states, Th(IV) releases from sTf because of its strong electrostatic interaction with D63 in the first shell of the sTf binding cleft. Analysis of hydrogen bond, water bridge, and the evaluation of potential of mean forces of the An ions' release from sTf, substantiate the differential behavior of Cm(III) and Th(IV) at endosomal pH. The results provide insight in the regulation of Cm(III) and Th(IV) bioavailability that may prove useful for effective design of their decorporating agents and as well may help the future design of radiotherapy based on tetravalent ions.
Collapse
Affiliation(s)
- Lokpati Mishra
- Radiation Safety Systems Division, Bhabha Atomic Research Centre, Mumbai, India.,Department of Chemical Sciences, Homi Bhabha National Institute, Mumbai, India
| | | | - Tusar Bandyopadhyay
- Department of Chemical Sciences, Homi Bhabha National Institute, Mumbai, India.,Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, India
| |
Collapse
|
30
|
Bissaro M, Sturlese M, Moro S. The rise of molecular simulations in fragment-based drug design (FBDD): an overview. Drug Discov Today 2020; 25:1693-1701. [PMID: 32592867 PMCID: PMC7314695 DOI: 10.1016/j.drudis.2020.06.023] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2020] [Revised: 05/24/2020] [Accepted: 06/19/2020] [Indexed: 12/31/2022]
Abstract
Fragment-based drug discovery (FBDD) is an innovative approach, progressively more applied in the academic and industrial context, to enhance hit identification for previously considered undruggable biological targets. In particular, FBDD discovers low-molecular-weight (LMW) ligands (<300Da) able to bind to therapeutically relevant macromolecules in an affinity range from the micromolar (μM) to millimolar (mM). X-ray crystallography (XRC) and nuclear magnetic resonance (NMR) spectroscopy are commonly the methods of choice to obtain 3D information about the bound ligand-protein complex, but this can occasionally be problematic, mainly for early, low-affinity fragments. The recent development of computational fragment-based approaches provides a further strategy for improving the identification of fragment hits. In this review, we summarize the state of the art of molecular dynamics simulations approaches used in FBDD, and discuss limitations and future perspectives for these approaches.
Collapse
Affiliation(s)
- Maicol Bissaro
- Molecular Modeling Section, Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Padova, Italy
| | - Mattia Sturlese
- Molecular Modeling Section, Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Padova, Italy
| | - Stefano Moro
- Molecular Modeling Section, Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Padova, Italy.
| |
Collapse
|
31
|
Torrens-Fontanals M, Stepniewski TM, Aranda-García D, Morales-Pastor A, Medel-Lacruz B, Selent J. How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs. Int J Mol Sci 2020; 21:E5933. [PMID: 32824756 PMCID: PMC7460635 DOI: 10.3390/ijms21165933] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 08/11/2020] [Accepted: 08/15/2020] [Indexed: 01/08/2023] Open
Abstract
G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.
Collapse
Affiliation(s)
- Mariona Torrens-Fontanals
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Tomasz Maciej Stepniewski
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
- InterAx Biotech AG, PARK innovAARE, 5234 Villigen, Switzerland
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, 02-093 Warsaw, Poland
| | - David Aranda-García
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Adrián Morales-Pastor
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Brian Medel-Lacruz
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Jana Selent
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| |
Collapse
|
32
|
Capelli R, Lyu W, Bolnykh V, Meloni S, Olsen JMH, Rothlisberger U, Parrinello M, Carloni P. Accuracy of Molecular Simulation-Based Predictions of koff Values: A Metadynamics Study. J Phys Chem Lett 2020; 11:6373-6381. [PMID: 32672983 DOI: 10.1021/acs.jpclett.0c00999] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
The koff values of ligands unbinding to proteins are key parameters for drug discovery. Their predictions based on molecular simulation may under- or overestimate experiment in a system- and/or technique-dependent way. Here we use an established method-infrequent metadynamics, based on the AMBER force field-to compute the koff of the ligand iperoxo (in clinical use) targeting the muscarinic receptor M2. The ligand charges are calculated by either (i) the Amber standard procedure or (ii) B3LYP-DFT. The calculations using (i) turn out not to provide a reasonable estimation of the transition-state free energy. Those using (ii) differ from experiment by 2 orders of magnitude. On the basis of B3LYP DFT QM/MM simulations, we suggest that the observed discrepancy in (ii) arises, at least in part, from the lack of electronic polarization and/or charge transfer in biomolecular force fields. These issues might be present in other systems, such as DNA-protein complexes.
Collapse
Affiliation(s)
- Riccardo Capelli
- Computational Biomedicine Section, IAS-5/INM-9, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
- JARA-HPC, Forschungszentrum Jülich, D-54245 Jülich, Germany
| | - Wenping Lyu
- Computational Biomedicine Section, IAS-5/INM-9, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
| | - Viacheslav Bolnykh
- Computational Biomedicine Section, IAS-5/INM-9, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Simone Meloni
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Ferrara, Via Luigi Borsari 46, I-44121, Ferrara, Italy
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana (USI), Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Istituto Italiano di Tecnologia, Via Morego 30, I-16163 Genova, Italy
| | - Paolo Carloni
- Computational Biomedicine Section, IAS-5/INM-9, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
- JARA-Institute INM-11: Molecular Neuroscience and Neuroimaging, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
| |
Collapse
|
33
|
Bekker GJ, Araki M, Oshima K, Okuno Y, Kamiya N. Exhaustive search of the configurational space of heat-shock protein 90 with its inhibitor by multicanonical molecular dynamics based dynamic docking. J Comput Chem 2020; 41:1606-1615. [PMID: 32267975 DOI: 10.1002/jcc.26203] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2020] [Revised: 03/18/2020] [Accepted: 03/25/2020] [Indexed: 01/02/2023]
Abstract
Multicanonical molecular dynamics based dynamic docking was used to exhaustively search the configurational space of an inhibitor binding to the N-terminal domain of heat-shock protein 90 (Hsp90). The obtained structures at 300 K cover a wide structural ensemble, with the top two clusters ranked by their free energy coinciding with the native binding site. The representative structure of the most stable cluster reproduced the experimental binding configuration, but an interesting conformational change in Hsp90 could be observed. The combined effects of solvation and ligand binding shift the equilibrium from a preferred loop-in conformation in the unbound state to an α-helical one in the bound state for the flexible lid region of Hsp90. Thus, our dynamic docking method is effective at predicting the native binding site while exhaustively sampling a wide configurational space, modulating the protein structure upon binding.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute for Protein Research, Osaka University, Osaka, Japan
| | - Mitsugu Araki
- Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Kanji Oshima
- Biotechnology Research Laboratories, Kaneka Corporation, Takasago, Hyogo, Japan
| | - Yasushi Okuno
- Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Narutoshi Kamiya
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Hyogo, Japan
| |
Collapse
|
34
|
Distinct binding of cetirizine enantiomers to human serum albumin and the human histamine receptor H 1. J Comput Aided Mol Des 2020; 34:1045-1062. [PMID: 32572668 DOI: 10.1007/s10822-020-00328-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2020] [Accepted: 06/18/2020] [Indexed: 02/02/2023]
Abstract
Cetirizine, a major metabolite of hydroxyzine, became a marketed second-generation H1 antihistamine that is orally active and has a rapid onset of action, long duration of effects and a very good safety record at recommended doses. The approved drug is a racemic mixture of (S)-cetirizine and (R)-cetirizine, the latter being the levorotary enantiomer that also exists in the market as a third-generation, non-sedating and highly selective antihistamine. Both enantiomers bind tightly to the human histamine H1 receptor (hH1R) and behave as inverse agonists but the affinity and residence time of (R)-cetirizine are greater than those of (S)-cetirizine. In blood plasma, cetirizine exists in the zwitterionic form and more than 90% of the circulating drug is bound to human serum albumin (HSA), which acts as an inactive reservoir. Independent X-ray crystallographic work has solved the structure of the hH1R:doxepin complex and has identified two drug-binding sites for cetirizine on equine serum albumin (ESA). Given this background, we decided to model a membrane-embedded hH1R in complex with either (R)- or (S)-cetirizine and also the complexes of both ESA and HSA with these two enantiomeric drugs to analyze possible differences in binding modes between enantiomers and also among targets. The ensuing molecular dynamics simulations in explicit solvent and additional computational chemistry calculations provided structural and energetic information about all of these complexes that is normally beyond current experimental possibilities. Overall, we found very good agreement between our binding energy estimates and extant biochemical and pharmacological evidence. A much higher degree of solvent exposure in the cetirizine-binding site(s) of HSA and ESA relative to the more occluded orthosteric binding site in hH1R is translated into larger positional fluctuations and considerably lower affinities for these two nonspecific targets. Whereas it is demonstrated that the two known pockets in ESA provide enough stability for cetirizine binding, only one such site does so in HSA due to a number of amino acid replacements. At the histamine-binding site in hH1R, the distinct interactions established between the phenyl and chlorophenyl moieties of the two enantiomers with the amino acids lining up the pocket and between their free carboxylates and Lys179 in the second extracellular loop account for the improved pharmacological profile of (R)-cetirizine.
Collapse
|
35
|
Salmaso V, Jacobson KA. In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors. Biomolecules 2020; 10:E812. [PMID: 32466404 PMCID: PMC7356333 DOI: 10.3390/biom10060812] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 05/21/2020] [Accepted: 05/22/2020] [Indexed: 12/14/2022] Open
Abstract
Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.
Collapse
Affiliation(s)
| | - Kenneth A. Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA;
| |
Collapse
|
36
|
Moesgaard L, Reinholdt P, Wüstner D, Kongsted J. Modeling the Sterol-Binding Domain of Aster-A Provides Insight into Its Multiligand Specificity. J Chem Inf Model 2020; 60:2268-2281. [PMID: 32233488 DOI: 10.1021/acs.jcim.0c00086] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Intracellular transport of cholesterol and related sterols relies to a large degree on nonvesicular mechanisms, which are only partly understood at the molecular level. Aster proteins belonging to the Lam family of sterol transfer proteins have recently been identified as important catalysts of nonvesicular sterol exchange between the plasma membrane (PM) and endoplasmic reticulum (ER). Here, we used a range of computational tools to study the molecular mechanisms underlying sterol binding as well as multisterol ligand specificity of Aster-A. Our study focused primarily on gaining atomistic insight into the bound ligand-protein complex and was, on this basis, performed in the absence of any membrane. Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations provide a rationale for the experimentally found ranking of binding affinities of various sterols to Aster-A. In particular, the polarity of the sterols and the length of their alkyl chain could be identified as being critical determinants of ligand affinity. A Gibbs free energy decomposition identified a charged residue, Glu444, at the base of the binding pose as an important control point for sterol binding. Removing its net charge via protonation was found to cause significant changes to the environment surrounding this residue. In addition, the protonation of Glu444 was found to be paralleled by a large redistribution of molecular flexibility in the Aster domain. This finding was supplemented by multiple branched adaptive steered molecular dynamics (MB-ASMD) simulations by which we defined a possible molecular path for sterol release and demonstrated the importance of Glu444 in this process.
Collapse
Affiliation(s)
- Laust Moesgaard
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Daniel Wüstner
- Department of Biochemistry and Molecular Biology, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| |
Collapse
|
37
|
Inuki S. [Elucidation of Biological Mechanisms Using Synthetic Natural Products and Their Derivatives]. YAKUGAKU ZASSHI 2020; 140:455-470. [PMID: 32238625 DOI: 10.1248/yakushi.19-00206] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Natural products are useful sources in the search for biochemical probes and drug leads because of their unique biological activities. However, synthetic studies or functional analyses of polycyclic complex natural products or conjugated lipids (e.g., glycolipids) are often hampered because of their synthesis and handling are challenging. On the basis of rational designs, synthetic studies, and chemical modifications, natural products need to be optimized to more potent compounds with improved activities, selectivities and/or physical properties. We have been synthesizing natural products and their derivatives for the elucidation of their biological mechanisms and discovery of drug leads. This review describes three topics for developing functional compounds derived from natural products for prospective involvement in pharmaceutical research: 1) direct construction of the ergot alkaloid scaffold by palladium catalyzed domino cyclization of amino allenes; 2) identification of novel sphingosine kinase inhibitors through a structure-activity relationship study of jaspine B; and 3) design, synthesis and biological evaluation of novel CD1d glycolipid ligands containing modified lipid moieties.
Collapse
Affiliation(s)
- Shinsuke Inuki
- Graduate School of Pharmaceutical Sciences, Kyoto University
| |
Collapse
|
38
|
Glöckner S, Ngo K, Wagner B, Heine A, Klebe G. The Influence of Varying Fluorination Patterns on the Thermodynamics and Kinetics of Benzenesulfonamide Binding to Human Carbonic Anhydrase II. Biomolecules 2020; 10:E509. [PMID: 32230853 PMCID: PMC7226267 DOI: 10.3390/biom10040509] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2020] [Revised: 03/25/2020] [Accepted: 03/26/2020] [Indexed: 12/12/2022] Open
Abstract
The fluorination of lead-like compounds is a common tool in medicinal chemistry to alter molecular properties in various ways and with different goals. We herein present a detailed study of the binding of fluorinated benzenesulfonamides to human Carbonic Anhydrase II by complementing macromolecular X-ray crystallographic observations with thermodynamic and kinetic data collected with the novel method of kinITC. Our findings comprise so far unknown alternative binding modes in the crystalline state for some of the investigated compounds as well as complex thermodynamic and kinetic structure-activity relationships. They suggest that fluorination of the benzenesulfonamide core is especially advantageous in one position with respect to the kinetic signatures of binding and that a higher degree of fluorination does not necessarily provide for a higher affinity or more favorable kinetic binding profiles. Lastly, we propose a relationship between the kinetics of binding and ligand acidity based on a small set of compounds with similar substitution patterns.
Collapse
Affiliation(s)
- Steffen Glöckner
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35037 Marburg, Germany; (S.G.); (K.N.); (A.H.)
| | - Khang Ngo
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35037 Marburg, Germany; (S.G.); (K.N.); (A.H.)
| | - Björn Wagner
- F. Hoffmann-La Roche AG, Pharmaceutical Research & Early Development, Roche Innovation Center Basel, 4070 Basel, Switzerland;
| | - Andreas Heine
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35037 Marburg, Germany; (S.G.); (K.N.); (A.H.)
| | - Gerhard Klebe
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35037 Marburg, Germany; (S.G.); (K.N.); (A.H.)
| |
Collapse
|
39
|
Deganutti G, Moro S, Reynolds CA. A Supervised Molecular Dynamics Approach to Unbiased Ligand–Protein Unbinding. J Chem Inf Model 2020; 60:1804-1817. [DOI: 10.1021/acs.jcim.9b01094] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Giuseppe Deganutti
- School of Life Sciences, University of Essex, Wivenhoe Park, Colchester, CO4 3SQ, United Kingdom
| | - Stefano Moro
- Molecular Modeling Section, Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Via Marzolo 5, 35131, Padova, Italy
| | - Christopher A. Reynolds
- School of Life Sciences, University of Essex, Wivenhoe Park, Colchester, CO4 3SQ, United Kingdom
| |
Collapse
|
40
|
Hüfner-Wulsdorf T, Klebe G. Role of Water Molecules in Protein–Ligand Dissociation and Selectivity Discrimination: Analysis of the Mechanisms and Kinetics of Biomolecular Solvation Using Molecular Dynamics. J Chem Inf Model 2020; 60:1818-1832. [DOI: 10.1021/acs.jcim.0c00156] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Tobias Hüfner-Wulsdorf
- Institut für Pharmazeutische Chemie, Philipps Universität Marburg, Marbacher Weg 6, 35037 Marburg, Germany
| | - Gerhard Klebe
- Institut für Pharmazeutische Chemie, Philipps Universität Marburg, Marbacher Weg 6, 35037 Marburg, Germany
| |
Collapse
|
41
|
Mishra L, Sundararajan M, Bandyopadhyay T. Molecular dynamics simulations of plutonium binding and its decorporation from the binding-cleft of human serum transferrin. J Biol Inorg Chem 2020; 25:213-231. [PMID: 31980924 DOI: 10.1007/s00775-020-01753-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Accepted: 01/03/2020] [Indexed: 11/29/2022]
Abstract
The possibility of plutonium (Pu) intake by radiation workers can not be ruled out. Transportation of Pu(IV) to various organs/cells is mainly carried through iron-carrying protein, serum transferrin (sTf), by receptor-mediated endocytosis. Understanding the Pu-sTf interaction is a primary step toward future design of its decorporating agents. We report MD simulations of Pu(IV) binding with sTf and look out for its decorporation at extracellular pH using suitable ligands. MD simulations were carried out in polarizable water environment at different protonation states of the protein. Results unravel the binding motif of Pu(IV): (1) sTf binds the ion in closed conformation at extracellular serum pH with carbonate as synergistic anions, (2) change in protonation state of dilysine (K206 and K296)-trigger and that of the carbonate ion at acidic endosomal pH is found to cause conformational changes of protein, conducive for the heavy ion to be released, although; (3) strong electrostatic interaction between D63 in the binding-cleft and Pu(IV) is found not to ever set free the ion. In an endeavour to decorporate Pu(IV), fragmented molecular form of hydroxypyridinone (HOPO) and catechol (CAM)-based ligands are docked at the binding site (BS) of the protein and metadynamics simulations are conducted. Pu(IV) binding at BS is found to be so strong that it was not detached from BS with the docked HOPO. However, for the identical set of simulation parameters, CAM is found to facilitate dislodging the heavy ion from the protein's binding influence. Differential behaviour of the two chelators is further explored. Fragmented molecular form of hydroxy-pyridinone (HOPO) and catecholamide (CAM) ligands were docked at the binding-site (BS) of human serum transferrin (sTf) to explore their feasibility as plausible Pu(IV) decorporating agents by employing metadynamics method. CAM was found to dislodge Pu from the sTf BS, while HOPO could not.
Collapse
Affiliation(s)
- Lokpati Mishra
- Radiation Safety Systems Division, Bhabha Atomic Research Centre, Mumbai, 400 085, India.,Homi Bhabha National Institute, Mumbai, 400094, India
| | - Mahesh Sundararajan
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, 400 085, India
| | - Tusar Bandyopadhyay
- Homi Bhabha National Institute, Mumbai, 400094, India. .,Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, 400 085, India.
| |
Collapse
|
42
|
Borisov DV, Veselovsky AV. [Ligand-receptor binding kinetics in drug design]. BIOMEDITSINSKAIA KHIMIIA 2020; 66:42-53. [PMID: 32116225 DOI: 10.18097/pbmc20206601042] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Traditionally, the thermodynamic values of affinity are considered as the main criterion for the development of new drugs. Usually, these values for drugs are measured <i>in vitro</i> at steady concentrations of the receptor and ligand, which are differed from <i>in vivo</i> environment. Recent studies have shown that the kinetics of the process of drug binding to its receptor make significant contribution in the drug effectiveness. This has increased attention in characterizing and predicting the rate constants of association and dissociation of the receptor ligand at the stage of preclinical studies of drug candidates. A drug with a long residence time can determine ligand-receptor selectivity (kinetic selectivity), maintain pharmacological activity of the drug at its low concentration in vivo. The paper discusses the theoretical basis of protein-ligand binding, molecular determinants that control the kinetics of the drug-receptor binding. Understanding the molecular features underlying the kinetics of receptor-ligand binding will contribute to the rational design of drugs with desired properties.
Collapse
Affiliation(s)
- D V Borisov
- Institute of Biomedical Chemistry, Moscow, Russia
| | | |
Collapse
|
43
|
Poongavanam V, Kongsted J, Wüstner D. Computational Modeling Explains the Multi Sterol Ligand Specificity of the N-Terminal Domain of Niemann-Pick C1-Like 1 Protein. ACS OMEGA 2019; 4:20894-20904. [PMID: 31867479 PMCID: PMC6921270 DOI: 10.1021/acsomega.9b01668] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Accepted: 10/18/2019] [Indexed: 06/10/2023]
Abstract
Niemann-Pick C1 like 1 (NPC1L1) is a sterol transporter expressed in the apical membrane of enterocytes and hepatocytes. NPC1L1 resembles the lysosomal NPC1 protein including an N-terminal domain (NTD), which binds a variety of sterols. The molecular mechanisms underlying this multiligand specificity of the NTD of NPC1L1 (NPC1L1-NTD) are not known. On the basis of the crystal structure of NPC1L1-NTD, we have investigated the structural details of protein-sterol interactions using molecular mechanics Poisson Boltzmann surface area calculations here. We found a good agreement between experimental and calculated binding affinities with similar ranking of various sterol ligands. We defined hydrogen bonding of sterol ligands via the 3'-β-hydroxy group inside the binding pose as instrumental in stabilizing the interaction. A leucine residue (LEU213) at the mouth of the binding pocket transiently opens to allow for the access of sterol into the binding pose. Our calculations also predict that NPC1L1-NTD binds polyene sterols, such as dehydroergosterol or cholestatrienol with high affinity, which validates their use in future experiments as close intrinsically fluorescent cholesterol analogs. A free energy decomposition and computational mutation analysis revealed that the binding of various sterols to NPC1L1-NTD depends critically on specific amino acid residues within the binding pocket. Some of these residues were previously detected as being relevant for intestinal cholesterol absorption. We show that clinically known mutations in the NPC1L1-NTD associated with lowered risk of coronary heart disease result in strongly reduced binding energies, providing a molecular explanation for the clinical phenotype.
Collapse
Affiliation(s)
- Vasanthanathan Poongavanam
- Department of Physics, Chemistry and Pharmacy and Department of Biochemistry and
Molecular Biology, University of Southern
Denmark, DK-5230 Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy and Department of Biochemistry and
Molecular Biology, University of Southern
Denmark, DK-5230 Odense M, Denmark
| | - Daniel Wüstner
- Department of Physics, Chemistry and Pharmacy and Department of Biochemistry and
Molecular Biology, University of Southern
Denmark, DK-5230 Odense M, Denmark
| |
Collapse
|
44
|
Majewski M, Ruiz-Carmona S, Barril X. An investigation of structural stability in protein-ligand complexes reveals the balance between order and disorder. Commun Chem 2019. [DOI: 10.1038/s42004-019-0205-5] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Abstract
The predominant view in structure-based drug design is that small-molecule ligands, once bound to their target structures, display a well-defined binding mode. However, structural stability (robustness) is not necessary for thermodynamic stability (binding affinity). In fact, it entails an entropic penalty that counters complex formation. Surprisingly, little is known about the causes, consequences and real degree of robustness of protein-ligand complexes. Since hydrogen bonds have been described as essential for structural stability, here we investigate 469 such interactions across two diverse structure sets, comprising of 79 drug-like and 27 fragment ligands, respectively. Completely constricted protein-ligand complexes are rare and may fulfill a functional role. Most complexes balance order and disorder by combining a single anchoring point with looser regions. 25% do not contain any robust hydrogen bond and may form loose structures. Structural stability analysis reveals a hidden layer of complexity in protein-ligand complexes that should be considered in ligand design.
Collapse
|
45
|
Horváth I, Jeszenői N, Bálint M, Paragi G, Hetényi C. A Fragmenting Protocol with Explicit Hydration for Calculation of Binding Enthalpies of Target-Ligand Complexes at a Quantum Mechanical Level. Int J Mol Sci 2019; 20:ijms20184384. [PMID: 31489952 PMCID: PMC6770515 DOI: 10.3390/ijms20184384] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2019] [Revised: 09/03/2019] [Accepted: 09/04/2019] [Indexed: 12/22/2022] Open
Abstract
Optimization of the enthalpy component of binding thermodynamics of drug candidates is a successful pathway of rational molecular design. However, the large size and missing hydration structure of target-ligand complexes often hinder such optimizations with quantum mechanical (QM) methods. At the same time, QM calculations are often necessitated for proper handling of electronic effects. To overcome the above problems, and help the QM design of new drugs, a protocol is introduced for atomic level determination of hydration structure and extraction of structures of target-ligand complex interfaces. The protocol is a combination of a previously published program MobyWat, an engine for assigning explicit water positions, and Fragmenter, a new tool for optimal fragmentation of protein targets. The protocol fostered a series of fast calculations of ligand binding enthalpies at the semi-empirical QM level. Ligands of diverse chemistry ranging from small aromatic compounds up to a large peptide helix of a molecular weight of 3000 targeting a leukemia protein were selected for systematic investigations. Comparison of various combinations of implicit and explicit water models demonstrated that the presence of accurately predicted explicit water molecules in the complex interface considerably improved the agreement with experimental results. A single scaling factor was derived for conversion of QM reaction heats into binding enthalpy values. The factor links molecular structure with binding thermodynamics via QM calculations. The new protocol and scaling factor will help automated optimization of binding enthalpy in future molecular design projects.
Collapse
Affiliation(s)
- István Horváth
- Chemistry Doctoral School, University of Szeged, Dugonics tér 13, 6720 Szeged, Hungary.
| | - Norbert Jeszenői
- Institute of Physiology, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Mónika Bálint
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Gábor Paragi
- MTA-SZTE Biomimetic Systems Research Group, Dóm tér 8, 6720 Szeged, Hungary.
- Institute of Physics, University of Pécs, Ifjúság útja 6, 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
|
46
|
Bruno A, Barresi E, Simola N, Da Pozzo E, Costa B, Novellino E, Da Settimo F, Martini C, Taliani S, Cosconati S. Unbinding of Translocator Protein 18 kDa (TSPO) Ligands: From in Vitro Residence Time to in Vivo Efficacy via in Silico Simulations. ACS Chem Neurosci 2019; 10:3805-3814. [PMID: 31268683 DOI: 10.1021/acschemneuro.9b00300] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Translocator protein 18 kDa (TSPO) is a validated pharmacological target for the development of new treatments for neurological disorders. N,N-Dialkyl-2-phenylindol-3-ylglyoxylamides (PIGAs) are effective TSPO modulators and potentially useful therapeutics for the treatment of anxiety, central nervous system pathologies featuring astrocyte loss, and inflammatory-based neuropathologies. For this class of compounds, no correlation exists between the TSPO binding affinity and the corresponding functional efficacy. Rather, their biological effectiveness correlates with the kinetics of the unbinding events and more specifically with the residence time (RT). So far, the structural reasons for the different recorded RT of congeneric PIGAs remain elusive. Here, to understand the different kinetics of PIGAs, their unbinding paths were studied by employing enhanced-sampling molecular dynamics simulations. Results of these studies revealed how subtle structural differences between PIGAs have a substantial effect on the unbinding energetics. In particular, during the egress from the TSPO binding site, slow-dissociating PIGAs find tight interactions with the protein LP1 region thereby determining a long RT. Further support to these findings was achieved by in vivo studies, which demonstrated how the anxiolytic effect observed for the inspected PIGAs correlated with their RT to TSPO.
Collapse
Affiliation(s)
- Agostino Bruno
- Department of Pharmacy, University Federico II of Naples, Via D. Montesano 49, 80131 Naples, Italy
| | - Elisabetta Barresi
- Department of Pharmacy, Università di Pisa, Via Bonanno 6, 56126 Pisa, Italy
| | - Nicola Simola
- Department of Biomedical Sciences, University of Cagliari, Monserrato University Campus, 09042 Monserrato, Italy
| | - Eleonora Da Pozzo
- Department of Pharmacy, Università di Pisa, Via Bonanno 6, 56126 Pisa, Italy
| | - Barbara Costa
- Department of Pharmacy, Università di Pisa, Via Bonanno 6, 56126 Pisa, Italy
| | - Ettore Novellino
- Department of Pharmacy, University Federico II of Naples, Via D. Montesano 49, 80131 Naples, Italy
| | - Federico Da Settimo
- Department of Pharmacy, Università di Pisa, Via Bonanno 6, 56126 Pisa, Italy
| | - Claudia Martini
- Department of Pharmacy, Università di Pisa, Via Bonanno 6, 56126 Pisa, Italy
| | - Sabrina Taliani
- Department of Pharmacy, Università di Pisa, Via Bonanno 6, 56126 Pisa, Italy
| | - Sandro Cosconati
- DiSTABiF, University of Campania “Luigi Vanvitelli”, Via Vivaldi 43, 81100 Caserta, Italy
| |
Collapse
|
47
|
Sykes DA, Jain P, Charlton SJ. Investigating the Influence of Tracer Kinetics on Competition-Kinetic Association Binding Assays: Identifying the Optimal Conditions for Assessing the Kinetics of Low-Affinity Compounds. Mol Pharmacol 2019; 96:378-392. [PMID: 31436538 DOI: 10.1124/mol.119.116764] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2019] [Accepted: 06/22/2019] [Indexed: 12/16/2022] Open
Abstract
An increased appreciation of the importance of optimizing drug-binding kinetics has lead to the development of various techniques for measuring the kinetics of unlabeled compounds. One approach is the competition-association kinetic binding method first described in the 1980s. The kinetic characteristics of the tracer employed greatly affects the reliability of estimated kinetic parameters, a barrier to successfully introducing these kinetic assays earlier in the drug discovery process. Using a modeling and Monte Carlo simulation approach, we identify the optimal tracer characteristics for determining the kinetics of the range of unlabeled ligands typically encountered during the different stages of a drug discovery program (i.e., rapidly dissociating, e.g., k off = 10 minute-1 low-affinity "hits" through to slowly dissociating e.g., k off = 0.01 minute-1 high-affinity "candidates"). For more rapidly dissociating ligands (e.g., k off = 10 minute-1), the key to obtaining accurate kinetic parameters was to employ a tracer with a relatively fast off-rate (e.g., k off = 1 minute-1) or, alternatively, to increase the tracer concentration. Reductions in assay start-time ≤1second and read frequency ≤5 seconds significantly improved the reliability of curve fitting. Timing constraints are largely dictated by the method of detection, its inherent sensitivity (e.g., TR-FRET versus radiometric detection), and the ability to inject samples online. Furthermore, we include data from TR-FRET experiments that validate this simulation approach, confirming its practical utility. These insights into the optimal experimental parameters for development of competition-association assays provide a framework for identifying and testing novel tracers necessary for profiling unlabeled competitors, particularly rapidly dissociating low-affinity competitors.
Collapse
Affiliation(s)
- David A Sykes
- School of Life Sciences, Queen's Medical Centre, University of Nottingham, Nottingham, United Kingdom (D.A.S., P.J., S.J.C.); Centre of Membrane and Protein and Receptors (COMPARE), University of Birmingham and University of Nottingham, Midlands, United Kingdom (D.A.S., P.J., S.J.C.); and Excellerate Bioscience Ltd, Discovery Building, BioCity, Nottingham, United Kingdom (S.J.C.)
| | - Palash Jain
- School of Life Sciences, Queen's Medical Centre, University of Nottingham, Nottingham, United Kingdom (D.A.S., P.J., S.J.C.); Centre of Membrane and Protein and Receptors (COMPARE), University of Birmingham and University of Nottingham, Midlands, United Kingdom (D.A.S., P.J., S.J.C.); and Excellerate Bioscience Ltd, Discovery Building, BioCity, Nottingham, United Kingdom (S.J.C.)
| | - Steven J Charlton
- School of Life Sciences, Queen's Medical Centre, University of Nottingham, Nottingham, United Kingdom (D.A.S., P.J., S.J.C.); Centre of Membrane and Protein and Receptors (COMPARE), University of Birmingham and University of Nottingham, Midlands, United Kingdom (D.A.S., P.J., S.J.C.); and Excellerate Bioscience Ltd, Discovery Building, BioCity, Nottingham, United Kingdom (S.J.C.)
| |
Collapse
|
48
|
Bosma R, Wang Z, Kooistra AJ, Bushby N, Kuhne S, van den Bor J, Waring MJ, de Graaf C, de Esch IJ, Vischer HF, Sheppard RJ, Wijtmans M, Leurs R. Route to Prolonged Residence Time at the Histamine H 1 Receptor: Growing from Desloratadine to Rupatadine. J Med Chem 2019; 62:6630-6644. [PMID: 31274307 PMCID: PMC6750840 DOI: 10.1021/acs.jmedchem.9b00447] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
![]()
Drug–target
binding kinetics are an important predictor of in vivo drug efficacy,
yet the relationship
between ligand structures and their binding kinetics is often poorly
understood. We show that both rupatadine (1) and desloratadine
(2) have a long residence time at the histamine H1 receptor (H1R). Through development of a [3H]levocetirizine radiolabel, we find that the residence time
of 1 exceeds that of 2 more than 10-fold.
This was further explored with 22 synthesized rupatadine and desloratadine
analogues. Methylene-linked cycloaliphatic or β-branched substitutions
of desloratadine increase the residence time at the H1R,
conveying a longer duration of receptor antagonism. However, cycloaliphatic
substituents directly attached to the piperidine amine (i.e., lacking
the spacer) have decreased binding affinity and residence time compared
to their methylene-linked structural analogues. Guided by docking
studies, steric constraints within the binding pocket are hypothesized
to explain the observed differences in affinity and binding kinetics
between analogues.
Collapse
Affiliation(s)
- Reggie Bosma
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Zhiyong Wang
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Albert J Kooistra
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Nick Bushby
- Operations, BioPharmaceuticals R&D , AstraZeneca , Alderley Park , Macclesfield SK10 4TG , United Kingdom
| | - Sebastiaan Kuhne
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Jelle van den Bor
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Michael J Waring
- Medicinal Chemistry, Research and Early Development, Oncology R&D , AstraZeneca , Alderley Park , Macclesfield SK10 4TG , United Kingdom
| | - Chris de Graaf
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Iwan J de Esch
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Henry F Vischer
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Robert J Sheppard
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D , AstraZeneca , Gothenburg 431 50 , Sweden
| | - Maikel Wijtmans
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| | - Rob Leurs
- Amsterdam Institute for Molecules, Medicines and Systems, Division of Medicinal Chemistry, Faculty of Science , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands
| |
Collapse
|
49
|
Sánchez-Murcia PA, Mills A, Cortés-Cabrera Á, Gago F. Unravelling the covalent binding of zampanolide and taccalonolide AJ to a minimalist representation of a human microtubule. J Comput Aided Mol Des 2019; 33:627-644. [DOI: 10.1007/s10822-019-00208-w] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Accepted: 05/24/2019] [Indexed: 01/27/2023]
|
50
|
IJzerman AP, Guo D. Drug-Target Association Kinetics in Drug Discovery. Trends Biochem Sci 2019; 44:861-871. [PMID: 31101454 DOI: 10.1016/j.tibs.2019.04.004] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2019] [Revised: 03/28/2019] [Accepted: 04/08/2019] [Indexed: 02/07/2023]
Abstract
The important role of ligand-receptor binding kinetics in drug design and discovery is increasingly recognized by the drug research community. Over the past decade, accumulating evidence has shown that optimizing the ligand's dissociation rate constant can lead to desirable duration of in vivo target occupancy and, hence, improved pharmacodynamic properties. However, the association rate constant as a pharmacological principle remains less investigated, whereas it can play an equally important role in the selection of drug candidates. This review provides a compilation and discussion of otherwise scarce and dispersed information on this topic, bringing to light the importance of drug-target association in kinetics-directed drug design and discovery.
Collapse
Affiliation(s)
- Adriaan P IJzerman
- Division of Drug Discovery and Safety, Leiden Academic Centre for Drug Research (LACDR), Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Dong Guo
- Jiangsu Key Laboratory of New Drug Research and Clinical Pharmacy, Xuzhou Medical University, 209 Tongshan Road, Xuzhou, 221004, Jiangsu, China.
| |
Collapse
|