1
|
Koenekoop L, Åqvist J. Computational Analysis of Heat Capacity Effects in Protein-Ligand Binding. J Chem Theory Comput 2024; 20:5708-5716. [PMID: 38870420 PMCID: PMC11238534 DOI: 10.1021/acs.jctc.4c00525] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2024]
Abstract
Heat capacity effects in protein-ligand binding as measured by calorimetric experiments have recently attracted considerable attention, particularly in the field of enzyme inhibitor design. A significant negative heat capacity change upon ligand binding implies a marked temperature dependence of the binding enthalpy, which is of high relevance for attempts to optimize protein-ligand interactions. In this work, we address the question of how well such heat capacity changes can be predicted by computer simulations. We examine a series of human thrombin inhibitors that all bind with ΔCp values of about -0.4 kcal/mol/K and calculate heat capacity changes from plain molecular dynamics simulations of the bound and free states of the enzyme and ligand. The results show that accurate ΔCp estimates within a few tenths of a kcal/mol/K of the experimental values can be obtained with this approach. This allows us to address the structural and energetic origin of the negative heat capacity changes for the thrombin inhibitors, and it is found that conformational equilibria of the free ligands in solution make a major contribution to the observed effect.
Collapse
Affiliation(s)
- Lucien Koenekoop
- Department of Cell & Molecular Biology, Uppsala University, Biomedical Center, SE-751 24 Uppsala, Sweden
| | - Johan Åqvist
- Department of Cell & Molecular Biology, Uppsala University, Biomedical Center, SE-751 24 Uppsala, Sweden
| |
Collapse
|
2
|
Kumar GS, Sobhia ME. Water network chemistry to exploit the nature of catalytic water molecules in Mtb DNA gyrase: a computational study to understand the binding mechanism of fluoroquinolones. J Biomol Struct Dyn 2024; 42:725-733. [PMID: 37121993 DOI: 10.1080/07391102.2023.2199869] [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: 09/23/2022] [Accepted: 03/17/2023] [Indexed: 05/02/2023]
Abstract
The dynamics of DNA gyrase and mutants of DNA gyrA such as G88A, A90V, S91P, D94A, D94G, D94N, D94Y; and double-point mutant (S91P-D94G), are meticulously investigated using computational approaches. Molecular dynamics (MD) and hydration thermodynamics have shed light on the fundamental, mechanistic basis of mutations on the conformational stability of Quinolone Binding Pocket (QBP) of DNA gyrase. Analysis of MD results revealed the displacement of a single crystal water molecule (HOH201) from the catalytic site of wild-type (WT) and mutants of DNA gyrA. This prompted our research group to probe the five crystal water molecules present in the QBP of the enzyme using water thermodynamics. Hydration thermodynamics analysis revealed the displacement of HOH201 due to unstable thermodynamic signatures. Further, the analysis highlighted significant changes in thermodynamic signatures and locations of five crystal water hydration sites upon mutation. Integrated MD simulations and water thermodynamics provided promising insights into the conformational changes and inaccessibility of the catalytic water molecule that can influence the design of DNA gyrase inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- G Siva Kumar
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), Mohali, Punjab, India
| | - M Elizabeth Sobhia
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), Mohali, Punjab, India
| |
Collapse
|
3
|
Szél V, Zsidó BZ, Jeszenői N, Hetényi C. Target-ligand binding affinity from single point enthalpy calculation and elemental composition. Phys Chem Chem Phys 2023; 25:31714-31725. [PMID: 37964670 DOI: 10.1039/d3cp04483a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
Reliable target-ligand binding thermodynamics data are essential for successful drug design and molecular engineering projects. Besides experimental methods, a number of theoretical approaches have been introduced for the generation of binding thermodynamics data. However, available approaches often neglect electronic effects or explicit water molecules influencing target-ligand interactions. To handle electronic effects within a reasonable time frame, we introduce a fast calculator QMH-L using a single target-ligand complex structure pre-optimized at the molecular mechanics level. QMH-L is composed of the semi-empirical quantum mechanics calculation of binding enthalpy with predicted explicit water molecules at the complex interface, and a simple descriptor based on the elemental composition of the ligand. QMH-L estimates the target-ligand binding free energy with a root mean square error (RMSE) of 0.94 kcal mol-1. The calculations also provide binding enthalpy values and they were compared with experimental binding thermodynamics data collected from the most reliable isothermal titration calorimetry studies of systems including various protein targets and challenging, large peptide ligands with a molecular weight of up to 2-3 thousand. The single point enthalpy calculations of QMH-L require modest computational resources and are based on short runs with open source and/or free software like Gromacs, Mopac, MobyWat, and Fragmenter. QMH-L can be applied for fast, automated scoring of drug candidates during a virtual screen, enthalpic engineering of new ligands or thermodynamic explanation of complex interactions.
Collapse
Affiliation(s)
- Viktor Szél
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Balázs Zoltán Zsidó
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Norbert Jeszenői
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Csaba Hetényi
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| |
Collapse
|
4
|
Çınaroğlu S, Biggin PC. Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems. J Chem Inf Model 2023; 63:6095-6108. [PMID: 37759363 PMCID: PMC10565830 DOI: 10.1021/acs.jcim.3c01041] [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: 07/10/2023] [Indexed: 09/29/2023]
Abstract
Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area.
Collapse
Affiliation(s)
| | - Philip C. Biggin
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
5
|
Hammerschmidt SJ, Maus H, Weldert AC, Gütschow M, Kersten C. Improving binding entropy by higher ligand symmetry? - A case study with human matriptase. RSC Med Chem 2023; 14:969-982. [PMID: 37252099 PMCID: PMC10211324 DOI: 10.1039/d3md00125c] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2023] [Accepted: 04/26/2023] [Indexed: 05/31/2023] Open
Abstract
Understanding different contributions to the binding entropy of ligands is of utmost interest to better predict affinity and the thermodynamic binding profiles of protein-ligand interactions and to develop new strategies for ligand optimization. To these means, the largely neglected effects of introducing higher ligand symmetry, thereby reducing the number of energetically distinguishable binding modes on binding entropy using the human matriptase as a model system, were investigated. A set of new trivalent phloroglucinol-based inhibitors that address the roughly symmetric binding site of the enzyme was designed, synthesized, and subjected to isothermal titration calorimetry. These highly symmetric ligands that can adopt multiple indistinguishable binding modes exhibited high entropy-driven affinity in line with affinity-change predictions.
Collapse
Affiliation(s)
- Stefan J Hammerschmidt
- Institute of Pharmaceutical and Biomedical Sciences, Johannes Gutenberg-University Staudingerweg 5 55128 Mainz Germany
| | - Hannah Maus
- Institute of Pharmaceutical and Biomedical Sciences, Johannes Gutenberg-University Staudingerweg 5 55128 Mainz Germany
| | - Annabelle C Weldert
- Institute of Pharmaceutical and Biomedical Sciences, Johannes Gutenberg-University Staudingerweg 5 55128 Mainz Germany
| | - Michael Gütschow
- Pharmaceutical Institute, Pharmaceutical & Medicinal Chemistry, University of Bonn An der Immenburg 4 53121 Bonn Germany
| | - Christian Kersten
- Institute of Pharmaceutical and Biomedical Sciences, Johannes Gutenberg-University Staudingerweg 5 55128 Mainz Germany
| |
Collapse
|
6
|
Omori A, Sasaki S, Kikukawa T, Shimono K, Miyauchi S. Elucidation of a Thermodynamical Feature Attributed to Substrate Binding to the Prokaryotic H +/Oligopeptide Cotransporter YdgR with Calorimetric Analysis: The Substrate Binding Driven by the Change in Entropy Implies the Release of Bound Water Molecules from the Binding Pocket. Biochemistry 2023. [PMID: 37163674 DOI: 10.1021/acs.biochem.2c00673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Here, we have elucidated the substrate recognition mechanism by a prokaryotic H+/oligopeptide cotransporter, YdgR, using isothermal titration calorimetry. Under acidic conditions (pH 6.0), the binding of a dipeptide, Val-Ala, to YdgR elicited endothermic enthalpy, which compensated for the increase in entropy due to dipeptide binding. A series of dipeptides were used in the binding titration. The dipeptides represent Val-X and X-Val, where X is Ala, Ser, Val, Tyr, or Phe. Most dipeptides revealed endothermic enthalpy, which was completely compensated by the increase in entropy due to dipeptide binding. The change in enthalpy due to binding correlated well with the change in entropy, whereas the Gibbs free energy involved in the binding of the dipeptide to YdgR remained unchanged irrespective of dipeptide sequences, implying that the binding reaction was driven by entropy, that is, the release of bound water molecules in the binding pocket. It is also important to clarify that, based on the prediction of water molecules in the ligand-binding pocket of YdgR, the release of three bound water molecules in the putative substrate binding pocket occurred through binding to YdgR. In the comparison of Val-X and X-Val dipeptides, the N-terminal region of the binding pocket might contain more bound water molecules than the C-terminal region. In light of these findings, we suggest that bound water molecules might play an important role in substrate recognition and binding by YdgR.
Collapse
Affiliation(s)
- Akiko Omori
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
| | - Shotaro Sasaki
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
| | - Takashi Kikukawa
- Faculty of Advanced Life Science, Hokkaido University, Kita 10, Nishi 8, Kita-ku, Sapporo 060-0810, Japan
| | - Kazumi Shimono
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
- Faculty of Pharmaceutical Sciences, Sojo University, 4-22-1 Ikeda Nishi-ku, Kumamoto 860-0082, Japan
| | - Seiji Miyauchi
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
| |
Collapse
|
7
|
Jiang L, Zhang D, Li Y, Chen W, Shi W, Wu H, Ma Z. Eukaryotic Expression of the Cytochrome c Oxidase Subunit I of Sitophilus zeamais and Its Interaction with Allyl Isothiocyanate. JOURNAL OF AGRICULTURAL AND FOOD CHEMISTRY 2023; 71:3497-3507. [PMID: 36757172 DOI: 10.1021/acs.jafc.2c08363] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Sitophilus zeamais Motschulsky (Coleoptera: Curculionidae) is a destructive pest of stored grains around the world. Allyl isothiocyanate (AITC) was shown to have good bioactivity in the control of S. zeamais. In this study, the interaction of AITC on cytochrome c oxidase core subunits I (COX I) and their binding mechanism were determined using spectroscopic, isothermal titration calorimetry and molecular docking techniques. The results indicate the binding constant (Ka) of AITC and COX I was 6.742 × 103 L/mol. Analysis of spectroscopic revealed that the binding of COX I to reduced Cyt c induced conformational changes of reduced Cyt c, while AITC could competitively bind and inhibit the activity of the COX I protein. Moreover, molecular docking results suggested a sulfur atom in the AITC structure could form a hydrogen bond having a length of 3.3 Å with the Gly- 27 of COX I.
Collapse
Affiliation(s)
- Linlin Jiang
- College of Plant Protection, Northwest A & F University, Yangling 712100, China
| | - Dan Zhang
- College of Plant Protection, Northwest A & F University, Yangling 712100, China
| | - Yue Li
- College of Plant Protection, Northwest A & F University, Yangling 712100, China
| | - Wei Chen
- College of Plant Protection, Northwest A & F University, Yangling 712100, China
| | - Weilin Shi
- College of Plant Protection, Northwest A & F University, Yangling 712100, China
| | - Hua Wu
- College of Plant Protection, Northwest A & F University, Yangling 712100, China
- Provincial Center for Bio-Pesticide Engineering, Yangling, Shaanxi Province 712100, China
| | - Zhiqing Ma
- College of Plant Protection, Northwest A & F University, Yangling 712100, China
- Provincial Center for Bio-Pesticide Engineering, Yangling, Shaanxi Province 712100, China
| |
Collapse
|
8
|
Kirley TL, Norman AB. Isothermal titration calorimetry determination of thermodynamics of binding of cocaine and its metabolites to humanized h2E2 anti-cocaine mAb. Biochem Biophys Rep 2022; 32:101354. [PMID: 36186732 PMCID: PMC9516381 DOI: 10.1016/j.bbrep.2022.101354] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2022] [Revised: 08/17/2022] [Accepted: 09/19/2022] [Indexed: 11/29/2022] Open
Abstract
We analyzed the thermodynamics of binding of cocaine and several cocaine metabolites to a humanized anti-cocaine mAb (h2E2), which is under development for the treatment of cocaine use disorders, using isothermal titration calorimetry. The calculated equilibrium dissociation (binding) constants were consistent with previous findings using other methods. All three ligands that display high affinity (nM) binding to the mAb (cocaine, cocaethylene, and benzoylecgonine) displayed similar enthalpically driven binding with substantial enthalpy-entropy compensation. The increased affinity of the cocaethylene metabolite compared to cocaine and benzoylecgonine is mostly attributable to a substantially less negative entropic binding component for cocaethylene, resulting in a more favorable binding energy, and thus, a higher affinity. The much lower affinity cocaine metabolites, norcocaine and ecgonine methyl ester, have much lower binding enthalpies than the high affinity ligands, and in contrast to the three high affinity ligands, have favorable (positive) entropic thermodynamic components of binding. Surprisingly, approximately 3.7 molecules of norcocaine are bound per mAb Fab site, as determined by isothermal titration calorimetry. This is in contrast to the three high affinity ligands, which bound with the expected stoichiometry of one drug molecule bound per one mAb Fab site. The results are discussed in relation to the previously published Fab:benzoylecgonine crystal structure for this h2E2 mAb, and compared to the isothermal titration calorimetry results published previously using an unrelated anti-cocaine mAb, mAb08. ITC was used to measure the thermodynamics of ligand binding to anti-cocaine h2E2 mAb. Binding of high and low affinity cocaine metabolites were compared to cocaine binding. Low affinity metabolites have favorable/positive entropic binding components. h2E2 mAb ITC results differ from those published for mAb08 anti-cocaine mAb. Surprisingly, ≈3.7 molecules of low affinity norcocaine are bound per mAb Fab site.
Collapse
|
9
|
Li Z, Chan KC, Nickels JD, Cheng X. Electrostatic Contributions to the Binding Free Energy of Nicotine to the Acetylcholine Binding Protein. J Phys Chem B 2022; 126:8669-8679. [PMID: 36260486 PMCID: PMC10056799 DOI: 10.1021/acs.jpcb.2c04641] [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] [Indexed: 01/11/2023]
Abstract
Biomolecular binding relies on specific attractive interactions between two partner molecules, including electrostatics, dispersion, hydrophobicity, and solvation. Assessing the contributions of electrostatic interactions to binding is key to the understanding of ligand binding mechanisms and the design of improved biomolecular binders. For example, nicotine is a well-known agonist of nicotinic acetylcholine receptors (nAChRs), but the molecular mechanisms for the differential action of nicotine on brain and muscle nAChRs remain elusive. In this work, we have chosen the acetylcholine binding protein (AChBP) in complex with nicotine as a model system to interrogate the electrostatic contributions to nicotine binding. Our absolute binding free energy simulations confirm that nicotine binds AChBP predominantly in its protonated (charged) form. By comparing energetic contributions from decomposed interactions for either neutral or charged nicotine, our calculations shed light on the nature of the binding of nicotine to the AChBP. The preferred binding of charged nicotine over neutral nicotine originates from its stronger electrostatic interactions with AChBP, a cation-π interaction to a tryptophan residue and a hydrogen bond between nicotine and the backbone carbonyl of the tryptophan, whereas the major force driving the binding process appears to be van der Waals interactions. The various nonelectrostatic terms can also indirectly modulate the electrostatic interactions through fine-tuning the binding pose of the ligand in the binding site, providing an explanation of why the binding specificity of nicotine to the brain versus muscle nAChRs is driven by electrostatic interaction, given that the immediate binding site residues, including the key tryptophan residue, are identical in the two receptors.
Collapse
Affiliation(s)
- Zoe Li
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy at The Ohio State University, Columbus, Ohio43210, United States
| | - Kevin C Chan
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy at The Ohio State University, Columbus, Ohio43210, United States
| | - Jonathan D Nickels
- Department of Chemical and Environmental Engineering, University of Cincinnati, Cincinnati, Ohio45221, United States
| | - Xiaolin Cheng
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy at The Ohio State University, Columbus, Ohio43210, United States
- Translational Data Analytics Institute (TDAI) at The Ohio State University, Columbus, Ohio43210, United States
| |
Collapse
|
10
|
Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization. Life (Basel) 2022; 12:life12091438. [PMID: 36143474 PMCID: PMC9504344 DOI: 10.3390/life12091438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 09/02/2022] [Accepted: 09/03/2022] [Indexed: 11/27/2022] Open
Abstract
Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying “magic bullets” amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand–target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural–energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.
Collapse
|
11
|
Åqvist J, van der Ent F. Calculation of Heat Capacity Changes in Enzyme Catalysis and Ligand Binding. J Chem Theory Comput 2022; 18:6345-6353. [PMID: 36094903 PMCID: PMC9558309 DOI: 10.1021/acs.jctc.2c00646] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
It has been suggested that heat capacity changes in enzyme
catalysis
may be the underlying reason for temperature optima that are not related
to unfolding of the enzyme. If this were to be a common phenomenon,
it would have major implications for our interpretation of enzyme
kinetics. In most cases, the support for the possible existence of
a nonzero (negative) activation heat capacity, however, only relies
on fitting such a kinetic model to experimental data. It is therefore
of fundamental interest to try to use computer simulations to address
this issue. One way is simply to calculate the temperature dependence
of the activation free energy and determine whether the relationship
is linear or not. An alternative approach is to calculate the absolute
heat capacities of the reactant and transition states from plain molecular
dynamics simulations using either the temperature derivative or fluctuation
formula for the enthalpy. Here, we examine these different approaches
for a designer enzyme with a temperature optimum that is not caused
by unfolding. Benchmark calculations for the heat capacity of liquid
water are first carried out using different thermostats. It is shown
that the derivative formula for the heat capacity is generally the
most robust and insensitive to the thermostat used and its parameters.
The enzyme calculations using this method give results in agreement
with direct calculations of activation free energies and show no sign
of a negative activation heat capacity. We also provide a simple scheme
for the calculation of binding heat capacity changes, which is of
clear interest in ligand design, and demonstrate it for substrate
binding to the designer enzyme. Neither in that case do the simulations
predict any negative heat capacity change.
Collapse
Affiliation(s)
- Johan Åqvist
- Department of Cell & Molecular Biology, Uppsala University, Biomedical Center, SE-751 24 Uppsala, Sweden
| | - Florian van der Ent
- Department of Cell & Molecular Biology, Uppsala University, Biomedical Center, SE-751 24 Uppsala, Sweden
| |
Collapse
|
12
|
Liu W, Jiang J, Lin Y, You Q, Wang L. Insight into Thermodynamic and Kinetic Profiles in Small-Molecule Optimization. J Med Chem 2022; 65:10809-10847. [PMID: 35969687 DOI: 10.1021/acs.jmedchem.2c00682] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Structure-activity relationships (SARs) and structure-property relationships (SPRs) have been considered the most important factors during the drug optimization process. For medicinal chemists, improvements in the potencies and druglike properties of small molecules are regarded as their major goals. Among them, the binding affinity and selectivity of small molecules on their targets are the most important indicators. In recent years, there has been growing interest in using thermodynamic and kinetic profiles to analyze ligand-receptor interactions, which could provide not only binding affinities but also detailed binding parameters for small-molecule optimization. In this perspective, we are trying to provide an insight into thermodynamic and kinetic profiles in small-molecule optimization. Through a highlight of strategies on the small-molecule optimization with specific cases, we aim to put forward the importance of structure-thermodynamic relationships (STRs) and structure-kinetic relationships (SKRs), which could provide more guidance to find safe and effective small-molecule drugs.
Collapse
Affiliation(s)
- Wei Liu
- State Key Laboratory of Natural Medicines and Jiangsu Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China.,Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Jingsheng Jiang
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Yating Lin
- State Key Laboratory of Natural Medicines and Jiangsu Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China.,Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Qidong You
- State Key Laboratory of Natural Medicines and Jiangsu Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China.,Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Lei Wang
- State Key Laboratory of Natural Medicines and Jiangsu Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China.,Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| |
Collapse
|
13
|
Gao Z, Fan T, Chen L, Yang M, Wai Wong VK, Chen D, Liu Z, Zhou Y, Wu W, Qiu Z, Zhang C, Li Y, Jiang Y. Design, synthesis and antitumor evaluation of novel 1H-indole-2-carboxylic acid derivatives targeting 14-3-3η protein. Eur J Med Chem 2022; 238:114402. [DOI: 10.1016/j.ejmech.2022.114402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Revised: 04/06/2022] [Accepted: 04/17/2022] [Indexed: 11/04/2022]
|
14
|
Begnini F, Geschwindner S, Johansson P, Wissler L, Lewis RJ, Danelius E, Luttens A, Matricon P, Carlsson J, Lenders S, König B, Friedel A, Sjö P, Schiesser S, Kihlberg J. Importance of Binding Site Hydration and Flexibility Revealed When Optimizing a Macrocyclic Inhibitor of the Keap1-Nrf2 Protein-Protein Interaction. J Med Chem 2022; 65:3473-3517. [PMID: 35108001 PMCID: PMC8883477 DOI: 10.1021/acs.jmedchem.1c01975] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Upregulation of the transcription factor Nrf2 by inhibition of the interaction with its negative regulator Keap1 constitutes an opportunity for the treatment of disease caused by oxidative stress. We report a structurally unique series of nanomolar Keap1 inhibitors obtained from a natural product-derived macrocyclic lead. Initial exploration of the structure-activity relationship of the lead, followed by structure-guided optimization, resulted in a 100-fold improvement in inhibitory potency. The macrocyclic core of the nanomolar inhibitors positions three pharmacophore units for productive interactions with key residues of Keap1, including R415, R483, and Y572. Ligand optimization resulted in the displacement of a coordinated water molecule from the Keap1 binding site and a significantly altered thermodynamic profile. In addition, minor reorganizations of R415 and R483 were accompanied by major differences in affinity between ligands. This study therefore indicates the importance of accounting both for the hydration and flexibility of the Keap1 binding site when designing high-affinity ligands.
Collapse
Affiliation(s)
- Fabio Begnini
- Department
of Chemistry—BMC, Uppsala University, Box 576, 75123 Uppsala, Sweden
| | - Stefan Geschwindner
- Mechanistic
and Structural Biology, Discovery Sciences, R&D, AstraZeneca, 43183 Mölndal, Sweden
| | - Patrik Johansson
- Mechanistic
and Structural Biology, Discovery Sciences, R&D, AstraZeneca, 43183 Mölndal, Sweden
| | - Lisa Wissler
- Mechanistic
and Structural Biology, Discovery Sciences, R&D, AstraZeneca, 43183 Mölndal, Sweden
| | - Richard J. Lewis
- Department
of Medicinal Chemistry, Research and Early Development, Respiratory
and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, 43183 Mölndal, Sweden
| | - Emma Danelius
- Department
of Chemistry—BMC, Uppsala University, Box 576, 75123 Uppsala, Sweden
| | - Andreas Luttens
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Box
596, 75124 Uppsala, Sweden
| | - Pierre Matricon
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Box
596, 75124 Uppsala, Sweden
| | - Jens Carlsson
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Box
596, 75124 Uppsala, Sweden
| | - Stijn Lenders
- Department
of Chemistry—BMC, Uppsala University, Box 576, 75123 Uppsala, Sweden
| | - Beate König
- Department
of Chemistry—BMC, Uppsala University, Box 576, 75123 Uppsala, Sweden
| | - Anna Friedel
- Department
of Chemistry—BMC, Uppsala University, Box 576, 75123 Uppsala, Sweden
| | - Peter Sjö
- Drugs
for Neglected Diseases Initiative (DNDi), 15 Chemin Camille-Vidart, 1202 Geneva, Switzerland
| | - Stefan Schiesser
- Department
of Medicinal Chemistry, Research and Early Development, Respiratory
and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, 43183 Mölndal, Sweden,
| | - Jan Kihlberg
- Department
of Chemistry—BMC, Uppsala University, Box 576, 75123 Uppsala, Sweden,
| |
Collapse
|
15
|
Wang H, Zhang M, Xu Y, Zong R, Xu N, Guo M. Agrobacterium fabrum atu0526-Encoding Protein Is the Only Chemoreceptor That Regulates Chemoattraction toward the Broad Antibacterial Agent Formic Acid. BIOLOGY 2021; 10:biology10121345. [PMID: 34943260 PMCID: PMC8698456 DOI: 10.3390/biology10121345] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Revised: 12/11/2021] [Accepted: 12/13/2021] [Indexed: 12/25/2022]
Abstract
Soil-born plant pathogens, especially Agrobacterium, generally navigate their way to hosts through recognition of the root exudates by chemoreceptors. However, there is still a lack of appropriate identification of chemoreceptors and their ligands in Agrobacterium. Here, Atu0526, a sCache-type chemoreceptor from Agrobacterium fabrum C58, was confirmed as the receptor of a broad antibacterial agent, formic acid. The binding of formic acid to Atu0526 was screened using a thermo shift assay and verified using isothermal titration calorimetry. Inconsistent with the previously reported antimicrobial properties, formic acid was confirmed to be a chemoattractant to A. fabrum and could promote its growth. The chemotaxis of A. fabrum C58 toward formic acid was completely lost with the knock-out of atu0526, and regained with the complementation of the gene, indicating that Atu0526 is the only chemoreceptor for formic acid in A. fabrum C58. The affinity of formic acid to Atu0526LBD significantly increased after the arginine at position 115 was replaced by alanine. However, in vivo experiments showed that the R115A mutation fully abolished the chemotaxis of A. fabrum toward formic acid. Molecular docking based on a predicted 3D structure of Atu0526 suggested that the arginine may provide "an anchorage" for formic acid to pull the minor loop, thereby forming a conformational change that generates the ligand-binding signal. Collectively, our findings will promote an understanding of sCache-type chemoreceptors and their signal transduction mechanism.
Collapse
|
16
|
Bennett JL, Nguyen GTH, Donald WA. Protein-Small Molecule Interactions in Native Mass Spectrometry. Chem Rev 2021; 122:7327-7385. [PMID: 34449207 DOI: 10.1021/acs.chemrev.1c00293] [Citation(s) in RCA: 53] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Small molecule drug discovery has been propelled by the continual development of novel scientific methodologies to occasion therapeutic advances. Although established biophysical methods can be used to obtain information regarding the molecular mechanisms underlying drug action, these approaches are often inefficient, low throughput, and ineffective in the analysis of heterogeneous systems including dynamic oligomeric assemblies and proteins that have undergone extensive post-translational modification. Native mass spectrometry can be used to probe protein-small molecule interactions with unprecedented speed and sensitivity, providing unique insights into polydisperse biomolecular systems that are commonly encountered during the drug discovery process. In this review, we describe potential and proven applications of native MS in the study of interactions between small, drug-like molecules and proteins, including large multiprotein complexes and membrane proteins. Approaches to quantify the thermodynamic and kinetic properties of ligand binding are discussed, alongside a summary of gas-phase ion activation techniques that have been used to interrogate the structure of protein-small molecule complexes. We additionally highlight some of the key areas in modern drug design for which native mass spectrometry has elicited significant advances. Future developments and applications of native mass spectrometry in drug discovery workflows are identified, including potential pathways toward studying protein-small molecule interactions on a whole-proteome scale.
Collapse
Affiliation(s)
- Jack L Bennett
- School of Chemistry, University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Giang T H Nguyen
- School of Chemistry, University of New South Wales, Sydney, New South Wales 2052, Australia
| | - William A Donald
- School of Chemistry, University of New South Wales, Sydney, New South Wales 2052, Australia
| |
Collapse
|
17
|
Uncertainty in protein-ligand binding constants: asymmetric confidence intervals versus standard errors. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2021; 50:661-670. [PMID: 33837826 DOI: 10.1007/s00249-021-01518-4] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2020] [Revised: 02/05/2021] [Accepted: 03/13/2021] [Indexed: 01/18/2023]
Abstract
Equilibrium binding constants (Kb) between chemical compounds and target proteins or between interacting proteins provide a quantitative understanding of biological interaction mechanisms. Reported uncertainties of measured experimental parameters are critical for decision-making in many scientific areas, e.g., in lead compound discovery processes and in comparing computational predictions with experimental results. Uncertainties in measured Kb values are commonly represented by a symmetric normal distribution, often quoted in terms of the experimental value plus-minus the standard deviation. However, in general, the distributions of measured Kb (and equivalent Kd) values and the corresponding free energy change ΔGb are all asymmetric to varying degree. Here, using a simulation approach, we illustrate the effect of asymmetric Kb distributions within the realm of isothermal titration calorimetry (ITC) experiments. Further we illustrate the known, but perhaps not widely appreciated, fact that when distributions of any of Kb, Kd and ΔGb are transformed into each other, their degree of asymmetry is changed. Consequently, we recommend that a more accurate way of expressing the uncertainties of Kb, Kd, and ΔGb values is to consistently report 95% confidence intervals, in line with other authors' suggestions. The ways to obtain such error ranges are discussed in detail and exemplified for a binding reaction obtained by ITC.
Collapse
|
18
|
Inhibitor development of MTH1 via high-throughput screening with fragment based library and MTH1 substrate binding cavity. Bioorg Chem 2021; 110:104813. [PMID: 33774493 DOI: 10.1016/j.bioorg.2021.104813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2021] [Revised: 02/25/2021] [Accepted: 03/06/2021] [Indexed: 11/22/2022]
Abstract
MutT Homolog 1 (MTH1) has been proven to hydrolyze oxidized nucleotide triphosphates during DNA repair. It can prevent the incorporation of wrong nucleotides during DNA replication and mitigate cell apoptosis. In a cancer cell, abundant reactive oxygen species can lead to substantial DNA damage and DNA mutations by base-pairing mismatch. MTH1 could eliminate oxidized dNTP and prevent cancer cells from entering cell death. Therefore, inhibition of MTH1 activity is considered to be an anti-cancer therapeutic target. In this study, high-throughput screening techniques were combined with a fragment-based library containing 2,313 compounds, which were used to screen for lead compounds with MTH1 inhibitor activity. Four compounds with MTH1 inhibitor ability were selected, and compound MI0639 was found to have the highest effective inhibition. To discover the selectivity and specificity of this action, several derivatives based on the MTH1 and MI0639 complex structure were synthesized. We compared 14 complex structures of MTH1 and the various compounds in combination with enzymatic inhibition and thermodynamic analysis. Nanomolar-range IC50 inhibition abilities by enzyme kinetics and Kd values by thermodynamic analysis were obtained for two compounds, named MI1020 and MI1024. Based on structural information and compound optimization, we aim to provide a strategy for the development of MTH1 inhibitors with high selectivity and specificity.
Collapse
|
19
|
Achazi K, Haag R, Ballauff M, Dernedde J, Kizhakkedathu JN, Maysinger D, Multhaup G. Understanding the Interaction of Polyelectrolyte Architectures with Proteins and Biosystems. Angew Chem Int Ed Engl 2021; 60:3882-3904. [PMID: 32589355 PMCID: PMC7894192 DOI: 10.1002/anie.202006457] [Citation(s) in RCA: 53] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Indexed: 02/06/2023]
Abstract
The counterions neutralizing the charges on polyelectrolytes such as DNA or heparin may dissociate in water and greatly influence the interaction of such polyelectrolytes with biomolecules, particularly proteins. In this Review we give an overview of studies on the interaction of proteins with polyelectrolytes and how this knowledge can be used for medical applications. Counterion release was identified as the main driving force for the binding of proteins to polyelectrolytes: Patches of positive charge become multivalent counterions of the polyelectrolyte and lead to the release of counterions from the polyelectrolyte and a concomitant increase in entropy. This is shown from investigations on the interaction of proteins with natural and synthetic polyelectrolytes. Special emphasis is paid to sulfated dendritic polyglycerols (dPGS). The Review demonstrates that we are moving to a better understanding of charge-charge interactions in systems of biological relevance. Research along these lines will aid and promote the design of synthetic polyelectrolytes for medical applications.
Collapse
Affiliation(s)
- Katharina Achazi
- Institut für Chemie und BiochemieFreie Universität BerlinTakustrasse 314195BerlinGermany
| | - Rainer Haag
- Institut für Chemie und BiochemieFreie Universität BerlinTakustrasse 314195BerlinGermany
| | - Matthias Ballauff
- Institut für Chemie und BiochemieFreie Universität BerlinTakustrasse 314195BerlinGermany
- IRIS AdlershofHumboldt Universität zu BerlinZum Grossen Windkanal 612489BerlinGermany
| | - Jens Dernedde
- Charité-Universitätsmedizin BerlinInstitute of Laboratory MedicineClinical Chemistry, and PathobiochemistryCVK Augustenburger Platz 113353BerlinGermany
| | - Jayachandran N. Kizhakkedathu
- Centre for Blood ResearchDepartment of Pathology and Laboratory MedicineLife Science InstituteDepartment of ChemistrySchool of Biomedical EngineeringUniversity of British ColumbiaVancouverV6T 1Z3Canada
| | - Dusica Maysinger
- Department of Pharmacology and TherapeuticsMcGill UniversityMontrealH3G 1Y6Canada
| | - Gerd Multhaup
- Department of Pharmacology and TherapeuticsMcGill UniversityMontrealH3G 1Y6Canada
| |
Collapse
|
20
|
Hüfner-Wulsdorf T, Klebe G. Advancing GIST-Based Solvent Functionals through Multiobjective Optimization of Solvent Enthalpy and Entropy Scoring Terms. J Chem Inf Model 2020; 60:6654-6665. [DOI: 10.1021/acs.jcim.0c01133] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
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
|
21
|
Duzgun Z, Eroglu Z. Role of cardiac drugs and flavonoids on the IRE1-JNK pathway as revealed by re-ranked molecular docking scores, MM/PBSA and umbrella sampling. J Biomol Struct Dyn 2020; 40:3428-3450. [PMID: 33251987 DOI: 10.1080/07391102.2020.1851299] [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: 10/22/2022]
Abstract
One of the important causes of cardiac dysfunction is the triggering of apoptosis through the IRE1-JNK signaling pathway due to excessive ER stress (endoplasmic reticulum stress). Although there are various studies on beneficial or harmful side effects of cardiac drugs, knowledge about the molecular mechanism of their interactions on this pathway is very limited. In this study, we investigated interactions of statins, ace inhibitors, antiarrhythmic drugs and flavonoids in IRE1, ASK1(apoptosis signal-regulating kinase 1) and JNK1 at an atomic level in comparison with their well-known inhibitors. The rank of scores obtained from four different docking algorithms (Autodock 4, Autodock Vina, iGEMDOCK and GOLD) were combined so that they could be compared with each other and evaluated together. According to combined results, the most potent compound for each compound group was selected for molecular dynamics simulations, MM/PBSA (molecular mechanics/Poisson-Boltzmann surface area) and umbrella sampling calculations. We observed that the statin group drugs had the best affinity by interacting with ASK1 and JNK1 by having a similar effect with their inhibitors, and atorvastatin and pitavastatin came to the fore. Norizalpinine from the flavonoid group had a strong binding interaction with IRE1, and amiodarone from the antiarrhythmic drug group had high binding affinities with IRE1, ASK1 and JNK1. Our study has shown that atorvastatin, pitavastatin, norizalpinine and amiodarone may have a role in preventing cardiac dysfunctions caused by ER stress and may shed light on further in vitro and in vivo research.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Zekeriya Duzgun
- Faculty of Medicine, Department of Medical Biology, Giresun University, Giresun, Turkey
| | - Zuhal Eroglu
- Faculty of Medicine, Department of Medical Biology, Ege University, Izmir, Turkey
| |
Collapse
|
22
|
Achazi K, Haag R, Ballauff M, Dernedde J, Kizhakkedathu JN, Maysinger D, Multhaup G. Wechselwirkung von Polyelektrolyt‐Architekturen mit Proteinen und Biosystemen. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202006457] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Katharina Achazi
- Institut für Chemie und Biochemie Freie Universität Berlin Takustraße 3 14195 Berlin Deutschland
| | - Rainer Haag
- Institut für Chemie und Biochemie Freie Universität Berlin Takustraße 3 14195 Berlin Deutschland
| | - Matthias Ballauff
- Institut für Chemie und Biochemie Freie Universität Berlin Takustraße 3 14195 Berlin Deutschland
- IRIS Adlershof Humboldt-Universität zu Berlin Zum Großen Windkanal 6 12489 Berlin Deutschland
| | - Jens Dernedde
- Charité-Universitätsmedizin Berlin Institut für Laboratoriumsmedizin Klinische Chemie und Pathobiochemie CVK Augustenburger Platz 1 13353 Berlin Deutschland
| | - Jayachandran N. Kizhakkedathu
- Centre for Blood Research Department of Pathology and Laboratory Medicine Life Science Institute Department of Chemistry School of Biomedical Engineering University of British Columbia Vancouver V6T 1Z3 Kanada
| | - Dusica Maysinger
- Department of Pharmacology and Therapeutics McGill University Montreal H3G 1Y6 Kanada
| | - Gerd Multhaup
- Department of Pharmacology and Therapeutics McGill University Montreal H3G 1Y6 Kanada
| |
Collapse
|
23
|
Åqvist J, Sočan J, Purg M. Hidden Conformational States and Strange Temperature Optima in Enzyme Catalysis. Biochemistry 2020; 59:3844-3855. [PMID: 32975950 PMCID: PMC7584337 DOI: 10.1021/acs.biochem.0c00705] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2020] [Revised: 09/17/2020] [Indexed: 11/29/2022]
Abstract
The existence of temperature optima in enzyme catalysis that occur before protein melting sets in can be described by different types of kinetic models. Such optima cause distinctly curved Arrhenius plots and have, for example, been observed in several cold-adapted enzymes from psychrophilic species. The two main explanations proposed for this behavior either invoke conformational equilibria with inactive substrate-bound states or postulate differences in heat capacity between the reactant and transition states. Herein, we analyze the implications of the different types of kinetic models in terms of apparent activation enthalpies, entropies, and heat capacities, using the catalytic reaction of a cold-adapted α-amylase as a prototypic example. We show that the behavior of these thermodynamic activation parameters is fundamentally different between equilibrium and heat capacity models, and in the α-amylase case, computer simulations have shown the former model to be correct. A few other enzyme-catalyzed reactions are also discussed in this context.
Collapse
Affiliation(s)
- Johan Åqvist
- Department of Cell &
Molecular Biology, Uppsala University, Biomedical Center, Box 596, SE-751 24 Uppsala, Sweden
| | - Jaka Sočan
- Department of Cell &
Molecular Biology, Uppsala University, Biomedical Center, Box 596, SE-751 24 Uppsala, Sweden
| | - Miha Purg
- Department of Cell &
Molecular Biology, Uppsala University, Biomedical Center, Box 596, SE-751 24 Uppsala, Sweden
| |
Collapse
|
24
|
Zhou HS, Hu LB, Zhang H, Shan WX, Wang Y, Li X, Liu T, Zhao J, You QD, Jiang ZY. Design, Synthesis, and Structure–Activity Relationships of Indoline-Based Kelch-like ECH-Associated Protein 1-Nuclear Factor (Erythroid-Derived 2)-Like 2 (Keap1-Nrf2) Protein–Protein Interaction Inhibitors. J Med Chem 2020; 63:11149-11168. [DOI: 10.1021/acs.jmedchem.0c01116] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Hai-Shan Zhou
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Lv-Bin Hu
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Han Zhang
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Wen-Xin Shan
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Yan Wang
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Xue Li
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Tian Liu
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Jing Zhao
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Qi-Dong You
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Zheng-Yu Jiang
- State Key Laboratory of Natural Medicines, and Jiang Su Key Laboratory of Drug Design and Optimization, China Pharmaceutical University, Nanjing 210009, China
- Department of Medicinal Chemistry, School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| |
Collapse
|
25
|
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
|
26
|
Syamala PPN, Würthner F. Modulation of the Self-Assembly of π-Amphiphiles in Water from Enthalpy- to Entropy-Driven by Enwrapping Substituents. Chemistry 2020; 26:8426-8434. [PMID: 32364616 PMCID: PMC7384034 DOI: 10.1002/chem.202000995] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2020] [Revised: 04/29/2020] [Indexed: 12/22/2022]
Abstract
Depending on the connectivity of solubilizing oligoethylene glycol (OEG) side chains to the π-cores of amphiphilic naphthalene and perylene bisimide dyes, self-assembly in water occurs either upon heating or cooling. Herein, we show that this effect originates from differences in the enwrapping capability of the π-cores by the OEG chains. Rylene bisimides bearing phenyl substituents with three OEG chains attached directly to the hydrophobic π-cores are strongly sequestered by the OEG chains. These molecules self-assemble at elevated temperatures in an entropy-driven process according to temperature- and concentration-dependent UV/Vis spectroscopy and calorimetric dilution studies. In contrast, for rylene bisimides in which phenyl substituents with three OEG chains are attached via a methylene spacer, leading to much weaker sequestration, self-assembly originates upon cooling in an enthalpy-driven process. Our explanation for this controversial behavior is that the aggregation in the latter case is dictated by the release of "high energy water" from the hydrophobic π-surfaces as well as dispersion interactions between the π-scaffolds which drive the self-assembly in an enthalpically driven process. In contrast, for the former case we suggest that in addition to the conventional explanation of a dehydration of hydrogen-bonded water molecules from OEG units it is in particular the increase in conformational entropy of back-folded OEG side chains upon aggregation that provides the pronounced gain in entropy that drives the aggregation process. Thus, our studies revealed that a subtle change in the attachment of solubilizing substituents can switch the thermodynamic signature for the self-assembly of amphiphilic dyes in water from enthalpy- to entropy-driven.
Collapse
Affiliation(s)
- Pradeep P N Syamala
- Institut für Organische Chemie, Universität Würzburg, Am Hubland, 97074, Würzburg, Germany.,Center for Nanosystems Chemistry (CNC), & Bavarian Polymer Institute (BPI), Universität Würzburg, Theodor-Boveri-Weg, 97074, Würzburg, Germany
| | - Frank Würthner
- Institut für Organische Chemie, Universität Würzburg, Am Hubland, 97074, Würzburg, Germany.,Center for Nanosystems Chemistry (CNC), & Bavarian Polymer Institute (BPI), Universität Würzburg, Theodor-Boveri-Weg, 97074, Würzburg, Germany
| |
Collapse
|
27
|
Mapping the S1 and S1' subsites of cysteine proteases with new dipeptidyl nitrile inhibitors as trypanocidal agents. PLoS Negl Trop Dis 2020; 14:e0007755. [PMID: 32163418 PMCID: PMC7067379 DOI: 10.1371/journal.pntd.0007755] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2019] [Accepted: 01/30/2020] [Indexed: 12/24/2022] Open
Abstract
The cysteine protease cruzipain is considered to be a validated target for therapeutic intervention in the treatment of Chagas disease. A series of 26 new compounds were designed, synthesized, and tested against the recombinant cruzain (Cz) to map its S1/S1´ subsites. The same series was evaluated on a panel of four human cysteine proteases (CatB, CatK, CatL, CatS) and Leishmania mexicana CPB, which is a potential target for the treatment of cutaneous leishmaniasis. The synthesized compounds are dipeptidyl nitriles designed based on the most promising combinations of different moieties in P1 (ten), P2 (six), and P3 (four different building blocks). Eight compounds exhibited a Ki smaller than 20.0 nM for Cz, whereas three compounds met these criteria for LmCPB. Three inhibitors had an EC50 value of ca. 4.0 μM, thus being equipotent to benznidazole according to the antitrypanosomal effects. Our mapping approach and the respective structure-activity relationships provide insights into the specific ligand-target interactions for therapeutically relevant cysteine proteases.
Collapse
|
28
|
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
|
29
|
Ferenczy GG, Keserű GM. Thermodynamic profiling for fragment-based lead discovery and optimization. Expert Opin Drug Discov 2019; 15:117-129. [PMID: 31741402 DOI: 10.1080/17460441.2020.1691166] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
Introduction: The enthalpic and entropic components of the ligand-protein binding free energy reflect the type and quality of the interactions and relate to the physicochemical properties of the ligands. These findings have significance in medicinal chemistry optimizations since they suggest that the thermodynamic profiling of the binding may help monitor and control the unfavorable size and hydrophobicity increase typically accompanying affinity improvements and leading to suboptimal pharmacokinetic properties.Areas covered: This review describes the ligand-protein binding event in terms of elementary steps, their associated interactions, and their enthalpic and entropic consequences. The relationships among the breaking and forming interactions, the binding thermodynamic profile, and the physicochemical properties of the ligands are also discussed.Expert opinion: Analysis of the size dependence of available affinity and favorable enthalpy highlights the limitation of the simultaneous optimization of these quantities. Indeed, moderate, rather than very high affinities can be conciliated with favorable physicochemical and pharmacokinetic profiles as it is supported by the affinity range of historical oral drugs. Although thermodynamic quantities are not suitable endpoints for medicinal chemistry optimizations owing to the complexity of the binding thermodynamics, thermodynamic profiling together with structural studies can be advantageously used to understand the details of the binding process and to optimize it.
Collapse
Affiliation(s)
- György G Ferenczy
- Medicinal Chemistry Research Group, Research Center for Natural Sciences, Budapest, Hungary
| | - György M Keserű
- Medicinal Chemistry Research Group, Research Center for Natural Sciences, Budapest, Hungary
| |
Collapse
|
30
|
Walkowiak J, Lu Y, Gradzielski M, Zauscher S, Ballauff M. Thermodynamic Analysis of the Uptake of a Protein in a Spherical Polyelectrolyte Brush. Macromol Rapid Commun 2019; 41:e1900421. [DOI: 10.1002/marc.201900421] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2019] [Revised: 09/19/2019] [Indexed: 12/14/2022]
Affiliation(s)
- Jacek Walkowiak
- Institut für Chemie und BiochemieFreie Universität Berlin Takustraße 3 14195 Berlin Germany
| | - Yan Lu
- Soft Matter and Functional MaterialsHelmholtz‐Zentrum Berlin für Materialen und Energie Hahn‐Meitner‐Platz 1 14109 Berlin Germany
- Institute of ChemistryUniversity of Potsdam 14467 Potsdam Germany
| | - Michael Gradzielski
- Stranski Laboratorium für Physikalische Chemie und Theoretische ChemieInstitut für ChemieStraße des 17. Juni 124Sekr. TC7Technische Universität Berlin D‐10623 Berlin Germany
| | - Stefan Zauscher
- Mechanical Engineering and Material ScienceDuke University Durham NC 27708 USA
| | - Matthias Ballauff
- Soft Matter and Functional MaterialsHelmholtz‐Zentrum Berlin für Materialen und Energie Hahn‐Meitner‐Platz 1 14109 Berlin Germany
| |
Collapse
|
31
|
Geschwindner S, Ulander J. The current impact of water thermodynamics for small-molecule drug discovery. Expert Opin Drug Discov 2019; 14:1221-1225. [DOI: 10.1080/17460441.2019.1664468] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Stefan Geschwindner
- Structure, Biophysics and Fragment-based Lead Generation, Discovery Sciences, R&D, AstraZeneca, Gothenburg, Sweden
| | - Johan Ulander
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca, Gothenburg, Sweden
| |
Collapse
|
32
|
Xu X, Ballauff M. Interaction of Lysozyme with a Dendritic Polyelectrolyte: Quantitative Analysis of the Free Energy of Binding and Comparison to Molecular Dynamics Simulations. J Phys Chem B 2019; 123:8222-8231. [DOI: 10.1021/acs.jpcb.9b07448] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Xiao Xu
- School of Chemical Engineering, Nanjing University of Science and Technology, 200 Xiao Ling Wei, 210094 Nanjing, P. R. China
| | - Matthias Ballauff
- Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, 14109 Berlin, Germany
- Institut für Physik, Humboldt-Universität zu Berlin, 12489 Berlin, Germany
| |
Collapse
|
33
|
Thermodynamic, kinetic, and structural parameterization of human carbonic anhydrase interactions toward enhanced inhibitor design. Q Rev Biophys 2019; 51:e10. [PMID: 30912486 DOI: 10.1017/s0033583518000082] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The aim of rational drug design is to develop small molecules using a quantitative approach to optimize affinity. This should enhance the development of chemical compounds that would specifically, selectively, reversibly, and with high affinity interact with a target protein. It is not yet possible to develop such compounds using computational (i.e., in silico) approach and instead the lead molecules are discovered in high-throughput screening searches of large compound libraries. The main reason why in silico methods are not capable to deliver is our poor understanding of the compound structure-thermodynamics and structure-kinetics correlations. There is a need for databases of intrinsic binding parameters (e.g., the change upon binding in standard Gibbs energy (ΔGint), enthalpy (ΔHint), entropy (ΔSint), volume (ΔVintr), heat capacity (ΔCp,int), association rate (ka,int), and dissociation rate (kd,int)) between a series of closely related proteins and a chemically diverse, but pharmacophoric group-guided library of compounds together with the co-crystal structures that could help explain the structure-energetics correlations and rationally design novel compounds. Assembly of these data will facilitate attempts to provide correlations and train data for modeling of compound binding. Here, we report large datasets of the intrinsic thermodynamic and kinetic data including over 400 primary sulfonamide compound binding to a family of 12 catalytically active human carbonic anhydrases (CA). Thermodynamic parameters have been determined by the fluorescent thermal shift assay, isothermal titration calorimetry, and by the stopped-flow assay of the inhibition of enzymatic activity. Kinetic measurements were performed using surface plasmon resonance. Intrinsic thermodynamic and kinetic parameters of binding were determined by dissecting the binding-linked protonation reactions of the protein and sulfonamide. The compound structure-thermodynamics and kinetics correlations reported here helped to discover compounds that exhibited picomolar affinities, hour-long residence times, and million-fold selectivities over non-target CA isoforms. Drug-lead compounds are suggested for anticancer target CA IX and CA XII, antiglaucoma CA IV, antiobesity CA VA and CA VB, and other isoforms. Together with 85 X-ray crystallographic structures of 60 compounds bound to six CA isoforms, the database should be of help to continue developing the principles of rational target-based drug design.
Collapse
|
34
|
Cianni L, Feldmann CW, Gilberg E, Gütschow M, Juliano L, Leitão A, Bajorath J, Montanari CA. Can Cysteine Protease Cross-Class Inhibitors Achieve Selectivity? J Med Chem 2019; 62:10497-10525. [DOI: 10.1021/acs.jmedchem.9b00683] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Lorenzo Cianni
- Medicinal Chemistry Group, Institute of Chemistry of São Carlos, University of São Paulo, Avenue Trabalhador Sancarlense, 400, 23566-590 São Carlos, SP, Brazil
- Pharmaceutical Institute, Pharmaceutical Chemistry I, University of Bonn, An der Immenburg 4, D-53121 Bonn, Germany
- Department of Life Science Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal Chemistry, Rheinische Friedrich-Wilhelms-Universität, Endenicher Allee 19c, D-53115 Bonn, Germany
| | - Christian Wolfgang Feldmann
- Department of Life Science Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal Chemistry, Rheinische Friedrich-Wilhelms-Universität, Endenicher Allee 19c, D-53115 Bonn, Germany
| | - Erik Gilberg
- Pharmaceutical Institute, Pharmaceutical Chemistry I, University of Bonn, An der Immenburg 4, D-53121 Bonn, Germany
- Department of Life Science Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal Chemistry, Rheinische Friedrich-Wilhelms-Universität, Endenicher Allee 19c, D-53115 Bonn, Germany
| | - Michael Gütschow
- Pharmaceutical Institute, Pharmaceutical Chemistry I, University of Bonn, An der Immenburg 4, D-53121 Bonn, Germany
| | - Luiz Juliano
- A. C. Camargo Cancer Center and São Paulo Medical School of Federal University of São Paulo, Rua Professor Antônio Prudente, 211, 01509-010 São Paulo, SP, Brazil
| | - Andrei Leitão
- Medicinal Chemistry Group, Institute of Chemistry of São Carlos, University of São Paulo, Avenue Trabalhador Sancarlense, 400, 23566-590 São Carlos, SP, Brazil
| | - Jürgen Bajorath
- Department of Life Science Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal Chemistry, Rheinische Friedrich-Wilhelms-Universität, Endenicher Allee 19c, D-53115 Bonn, Germany
| | - Carlos A. Montanari
- Medicinal Chemistry Group, Institute of Chemistry of São Carlos, University of São Paulo, Avenue Trabalhador Sancarlense, 400, 23566-590 São Carlos, SP, Brazil
| |
Collapse
|
35
|
Gao S, Holkar A, Srivastava S. Protein-Polyelectrolyte Complexes and Micellar Assemblies. Polymers (Basel) 2019; 11:E1097. [PMID: 31261765 PMCID: PMC6680422 DOI: 10.3390/polym11071097] [Citation(s) in RCA: 57] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2019] [Revised: 06/20/2019] [Accepted: 06/24/2019] [Indexed: 12/18/2022] Open
Abstract
In this review, we highlight the recent progress in our understanding of the structure, properties and applications of protein-polyelectrolyte complexes in both bulk and micellar assemblies. Protein-polyelectrolyte complexes form the basis of the genetic code, enable facile protein purification, and have emerged as enterprising candidates for simulating protocellular environments and as efficient enzymatic bioreactors. Such complexes undergo self-assembly in bulk due to a combined influence of electrostatic interactions and entropy gains from counterion release. Diversifying the self-assembly by incorporation of block polyelectrolytes has further enabled fabrication of protein-polyelectrolyte complex micelles that are multifunctional carriers for therapeutic targeted delivery of proteins such as enzymes and antibodies. We discuss research efforts focused on the structure, properties and applications of protein-polyelectrolyte complexes in both bulk and micellar assemblies, along with the influences of amphoteric nature of proteins accompanying patchy distribution of charges leading to unique phenomena including multiple complexation windows and complexation on the wrong side of the isoelectric point.
Collapse
Affiliation(s)
- Shang Gao
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, CA 90095, USA.
| | - Advait Holkar
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, CA 90095, USA.
| | - Samanvaya Srivastava
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, CA 90095, USA.
| |
Collapse
|
36
|
Umuhire Juru A, Patwardhan NN, Hargrove AE. Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions. ACS Chem Biol 2019; 14:824-838. [PMID: 31042354 DOI: 10.1021/acschembio.8b00945] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( Ka or Kd) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.
Collapse
Affiliation(s)
- Aline Umuhire Juru
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Neeraj N. Patwardhan
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Amanda E. Hargrove
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
37
|
Xu X, Angioletti-Uberti S, Lu Y, Dzubiella J, Ballauff M. Interaction of Proteins with Polyelectrolytes: Comparison of Theory to Experiment. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:5373-5391. [PMID: 30095921 DOI: 10.1021/acs.langmuir.8b01802] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We discuss recent investigations of the interaction of polyelectrolytes with proteins. In particular, we review our recent studies on the interaction of simple proteins such as human serum albumin (HSA) and lysozyme with linear polyelectrolytes, charged dendrimers, charged networks, and polyelectrolyte brushes. In all cases discussed here, we combined experimental work with molecular dynamics (MD) simulations and mean-field theories. In particular, isothermal titration calorimetry (ITC) has been employed to obtain the respective binding constants Kb and the Gibbs free energy of binding. MD simulations with explicit counterions but implicit water demonstrate that counterion release is the main driving force for the binding of proteins to strongly charged polyelectrolytes: patches of positive charges located on the surface of the protein become multivalent counterions of the polyelectrolyte, thereby releasing a number of counterions condensed on the polyelectrolyte. The binding Gibbs free energy due to counterion release is predicted to scale with the logarithm of the salt concentration in the system, which is verified by both simulations and experiment. In several cases, namely, for the interaction of proteins with linear polyelectrolytes and highly charged hydrophilic dendrimers, the binding constant could be calculated from simulations to very good approximation. This finding demonstrated that in these cases explicit hydration effects do not contribute to the Gibbs free energy of binding. The Gibbs free energy can also be used to predict the kinetics of protein uptake by microgels for a given system by applying dynamic density functional theory. The entire discussion demonstrates that the direct comparison of theory with experiments can lead to a full understanding of the interaction of proteins with charged polymers. Possible implications for applications, such as drug design, are discussed.
Collapse
Affiliation(s)
- Xiao Xu
- School of Chemical Engineering , Nanjing University of Science and Technology , 200 Xiao Ling Wei , Nanjing 210094 , P. R. China
| | - Stefano Angioletti-Uberti
- Department of Materials , Imperial College London , London SW7 2AZ - UK , U.K
- International Research Centre for Soft Matter , Beijing University of Chemical Technology , 100099 Beijing , PR China
| | - Yan Lu
- Soft Matter and Functional Materials , Helmholtz-Zentrum Berlin für Materialien und Energie GmbH , 14109 Berlin , Germany
- Institute of Chemistry , University of Potsdam , 14467 Potsdam , Germany
| | - Joachim Dzubiella
- Soft Matter and Functional Materials , Helmholtz-Zentrum Berlin für Materialien und Energie GmbH , 14109 Berlin , Germany
- Physikalisches Institut , Albert-Ludwigs-Universität , 79104 Freiburg , Germany
| | - Matthias Ballauff
- Soft Matter and Functional Materials , Helmholtz-Zentrum Berlin für Materialien und Energie GmbH , 14109 Berlin , Germany
- Institut für Physik , Humboldt-Universität zu Berlin , 12489 Berlin , Germany
| |
Collapse
|
38
|
Williams LJ, Schendt BJ, Fritz ZR, Attali Y, Lavroff RH, Yarmush ML. A protein interaction free energy model based on amino acid residue contributions: Assessment of point mutation stability of T4 lysozyme. TECHNOLOGY 2019; 7:12-39. [PMID: 32211456 PMCID: PMC7093156 DOI: 10.1142/s233954781950002x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Here we present a model to estimate the interaction free energy contribution of each amino acid residue of a given protein. Protein interaction energy is described in terms of per-residue interaction factors, μ. Multibody interactions are implicitly captured in μ through the combination of amino acid terms (γ) guided by local conformation indices (σ). The model enables construction of an interaction factor heat map for a protein in a given fold, allows prima facie assessment of the degree of residue-residue interaction, and facilitates a qualitative and quantitative evaluation of protein association properties. The model was used to compute thermal stability of T4 bacteriophage lysozyme mutants across seven sites. Qualitative assessment of mutational effects provides a straightforward rationale regarding whether a particular site primarily perturbs native or non-native states, or both. The presented model was found to be in good agreement with experimental mutational data (R 2 = 0.73) and suggests an approach by which to convert structure space into energy space.
Collapse
Affiliation(s)
- Lawrence J Williams
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, 123 Bevier Rd., Piscataway, NJ 08854, USA
| | - Brian J Schendt
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, 123 Bevier Rd., Piscataway, NJ 08854, USA
| | - Zachary R Fritz
- Department of Biomedical Engineering, Rutgers, The State University of New Jersey, 599 Taylor Road, Piscataway, NJ 08854, USA
| | - Yonatan Attali
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, 123 Bevier Rd., Piscataway, NJ 08854, USA
| | - Robert H Lavroff
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, 123 Bevier Rd., Piscataway, NJ 08854, USA
| | - Martin L Yarmush
- Department of Biomedical Engineering, Rutgers, The State University of New Jersey, 599 Taylor Road, Piscataway, NJ 08854, USA
| |
Collapse
|
39
|
Klebe G. Broad-scale analysis of thermodynamic signatures in medicinal chemistry: are enthalpy-favored binders the better development option? Drug Discov Today 2019; 24:943-948. [PMID: 30708050 DOI: 10.1016/j.drudis.2019.01.014] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Revised: 01/18/2019] [Accepted: 01/24/2019] [Indexed: 12/30/2022]
Abstract
Thermodynamic profiles of ligand binding, particularly enthalpically favored binding signatures, have been suggested as a criterion to support the decision-making process around which compounds to select for further optimization in drug development. The concept was enthusiastically taken up, but turned out to be too superficial, either because many aspects determining thermodynamic profiles are insufficiently appreciated or because it is difficult to compare such data on a global scale. The impact of water, changes in protonation states, along with buffer dependencies and incompatible measurement conditions that are far from standard conditions hamper such broad-scale comparisons. However, thermodynamic signatures can make us aware of the impact of these aspects and provide important hints for improving our understanding of the binding process and defining criteria for drug optimization.
Collapse
Affiliation(s)
- Gerhard Klebe
- Institute of Pharmaceutical Chemistry, Philipps-University Marburg, Marbacher Weg 6, 35032 Marburg, Germany.
| |
Collapse
|
40
|
Ziehm T, Buell AK, Willbold D. Role of Hydrophobicity and Charge of Amyloid-Beta Oligomer Eliminating d-Peptides in the Interaction with Amyloid-Beta Monomers. ACS Chem Neurosci 2018; 9:2679-2688. [PMID: 29893543 DOI: 10.1021/acschemneuro.8b00132] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
Inhibition of the self-assembly process of amyloid-beta and even more the removal of already existing toxic amyloid-beta assemblies represent promising therapeutic strategies against Alzheimer's disease. To approach this aim, we selected a d-enantiomeric peptide by phage-display based on the interaction with amyloid-beta monomers. This lead compound was successfully optimized by peptide microarrays with respect to its affinity and specificity to the target resulting in d-peptides with both increased hydrophobicity and net charge. Here, we present a detailed biophysical characterization of the interactions between these optimized d-peptides and amyloid-beta monomers in comparison to the original lead compound in order to obtain a more thorough understanding of the physicochemical determinants of the interactions. Kinetics and apparent stoichiometry of complex formation were studied using surface plasmon resonance. Potential modes of binding to amyloid-beta were identified, and the influences of ionic strength on complex stability, as well as on the inhibitory effect on amyloid-beta aggregation were investigated. The results reveal a very different mode of interaction of the optimized d-peptides based on a combination of electrostatic and hydrophobic interactions as compared to the mostly electrostatically driven interaction of the lead compound. These conclusions were supported by the thermodynamic profiles of the interaction between optimized d-peptides and Aβ monomers, which indicate an increase in binding entropy with respect to the lead compound.
Collapse
Affiliation(s)
- Tamar Ziehm
- Institute of Complex Systems, Structural Biochemistry (ICS-6), Forschungszentrum Jülich, 52425 Jülich, Germany
| | - Alexander K. Buell
- Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
| | - Dieter Willbold
- Institute of Complex Systems, Structural Biochemistry (ICS-6), Forschungszentrum Jülich, 52425 Jülich, Germany
- Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
41
|
Ran Q, Xu X, Dey P, Yu S, Lu Y, Dzubiella J, Haag R, Ballauff M. Interaction of human serum albumin with dendritic polyglycerol sulfate: Rationalizing the thermodynamics of binding. J Chem Phys 2018; 149:163324. [DOI: 10.1063/1.5030601] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qidi Ran
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
- Institute of Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Multifunctional Biomaterials for Medicine, Helmholtz Virtual Institute, Kantstr. 55, 14513 Teltow-Seehof, Germany
| | - Xiao Xu
- School of Chemical Engineering, Nanjing University of Science and Technology, 200 Xiao Ling Wei, 210094 Nanjing, People’s Republic of China
| | - Pradip Dey
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
| | - Shun Yu
- Institute of Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, 12489 Berlin, Germany
| | - Yan Lu
- Institute of Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
| | - Joachim Dzubiella
- Institute of Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Multifunctional Biomaterials for Medicine, Helmholtz Virtual Institute, Kantstr. 55, 14513 Teltow-Seehof, Germany
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| | - Rainer Haag
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
- Multifunctional Biomaterials for Medicine, Helmholtz Virtual Institute, Kantstr. 55, 14513 Teltow-Seehof, Germany
| | - Matthias Ballauff
- Institute of Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Multifunctional Biomaterials for Medicine, Helmholtz Virtual Institute, Kantstr. 55, 14513 Teltow-Seehof, Germany
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, 12489 Berlin, Germany
| |
Collapse
|
42
|
Menzer WM, Li C, Sun W, Xie B, Minh DDL. Simple Entropy Terms for End-Point Binding Free Energy Calculations. J Chem Theory Comput 2018; 14:6035-6049. [PMID: 30296084 DOI: 10.1021/acs.jctc.8b00418] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We introduce a number of computationally inexpensive modifications to the MM/PBSA and MM/GBSA estimators for binding free energies, which are based on average receptor-ligand interaction energies in simulations of a noncovalent complex, to improve the treatment of entropy: second- and higher-order terms in a cumulant expansion and a confining potential on ligand external degrees of freedom. We also consider a filter for snapshots where ligands have drifted from the initial binding pose. The variations were tested on six sets of systems for which binding modes and free energies have previously been experimentally determined. For some data sets, none of the tested estimators led to results significantly correlated with measured free energies. In data sets with nontrivial correlation, a ligand RMSD cutoff of 3 Å and a second-order truncation of the cumulant expansion was found to be comparable or better than the average interaction energy by several statistical metrics.
Collapse
|
43
|
Su H, Xu Y. Application of ITC-Based Characterization of Thermodynamic and Kinetic Association of Ligands With Proteins in Drug Design. Front Pharmacol 2018; 9:1133. [PMID: 30364164 PMCID: PMC6193069 DOI: 10.3389/fphar.2018.01133] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 09/18/2018] [Indexed: 01/20/2023] Open
Abstract
A comprehensive characterization of the thermodynamic and kinetic profiling of ligands binding to a given target protein is crucial for the hit selection as well as the hit-to-lead-to-drug evolution. Isothermal titration calorimetry (ITC), widely known as an invaluable tool to measure the thermodynamic data, has recently found its way to determine the binding kinetics too. The extensive application of ITC in measurement of both thermodynamic and kinetic data manifests unique roles of ITC in drug discovery and development. This mini-review concentrates on elaborating how to gain the thermodynamic and kinetic data using ITC, highlighting the importance of these data in lead discovery and optimization, and intends to provide an overview of the technical and conceptual advances that offer unprecedented access to protein–ligand recognition by ITC measurement.
Collapse
Affiliation(s)
- Haixia Su
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China.,School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
| | - Yechun Xu
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China.,School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
44
|
Henderson JA, Harris RC, Tsai CC, Shen J. How Ligand Protonation State Controls Water in Protein-Ligand Binding. J Phys Chem Lett 2018; 9:5440-5444. [PMID: 30188715 PMCID: PMC6467052 DOI: 10.1021/acs.jpclett.8b02440] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
The role of water in protein-ligand binding has been an intensely studied topic in recent years; however, how ligand protonation state change perturbs water has not been considered. Here we show that water dynamics and interactions can be controlled by the protonation state of ligand using continuous constant pH molecular dynamics simulations of two closely related model systems, β-secretase 1 and 2 (BACE1 and BACE2), in complex with a small-molecule inhibitor. Simulations revealed that, upon binding, the inhibitor pyrimidine ring remains deprotonated in BACE1 but becomes protonated in BACE2. Pyrimidine protonation results in water displacement, rigidification of the binding pocket, and shift in the ligand binding mode from water-mediated to direct hydrogen bonding. These findings not only support but also rationalize the most recent structure-selectivity data in BACE1 drug design. Binding-induced protonation state changes are likely common; our work offers a glimpse at how modeling protein-ligand binding while allowing ligand titration can further advance the understanding of water and structure-based drug design.
Collapse
Affiliation(s)
- Jack A Henderson
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Robert C Harris
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Cheng-Chieh Tsai
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Jana Shen
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| |
Collapse
|
45
|
Ran Q, Xu X, Dzubiella J, Haag R, Ballauff M. Thermodynamics of the Binding of Lysozyme to a Dendritic Polyelectrolyte: Electrostatics Versus Hydration. ACS OMEGA 2018; 3:9086-9095. [PMID: 31459043 PMCID: PMC6644519 DOI: 10.1021/acsomega.8b01493] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Accepted: 07/30/2018] [Indexed: 05/25/2023]
Abstract
The interaction between dendritic polyglycerol sulfate (dPGS) of the second generation and lysozyme was studied by isothermal titration calorimetry (ITC) at different temperatures and salt concentrations. Analysis by ITC showed that 2-3 lysozyme molecules were bound to each dPGS. The resulting binding constant K b and the Gibbs free energy ΔG o decreased markedly with increasing salt concentration but were nearly independent of temperature. The salt dependence of K b led to the conclusion that ca. 3 counterions bound to dPGS were released upon complex formation. The gain in entropy ΔG ci by this counterion-release scales logarithmically with salt concentration and is the main driving force for binding. The temperature dependence of ΔG o was analyzed by the nonlinear van't Hoff plot, taking into account a finite heat capacity change ΔC p,vH. This evaluation led to the binding enthalpy ΔH vH and the binding entropy ΔS vH. Both quantities varied strongly with temperature and even changed sign, but they compensated each other throughout the entire range of temperature. Coarse-grained computer simulations with explicit salt and implicit water were used to obtain the binding free energies that agreed with ITC results. Thus, electrostatic factors were the driving forces for binding whereas all hydration contributions leading to the strongly varying ΔH vH and ΔS vH canceled out. The calorimetric enthalpy ΔH ITC measured directly by ITC differed largely from ΔH vH. ITC measurements done in two buffer systems with different ionization enthalpies revealed that binding was linked to buffer ionization and a partial protonation of the protein.
Collapse
Affiliation(s)
- Qidi Ran
- Institute
of Chemistry and Biochemistry, Freie Universität
Berlin, Takustr. 3, 14195 Berlin, Germany
- Institute
of Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Multifunctional
Biomaterials for Medicine, Helmholtz Virtual Institute, Kantstr. 55, 14513 Teltow-Seehof, Germany
| | - Xiao Xu
- School
of Chemical Engineering, Nanjing University
of Science and Technology, 200 Xiao Ling Wei, 210094 Nanjing, P. R. China
| | - Joachim Dzubiella
- Institute
of Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Multifunctional
Biomaterials for Medicine, Helmholtz Virtual Institute, Kantstr. 55, 14513 Teltow-Seehof, Germany
- Physikalisches
Institut, Albert-Ludwigs-Universität, Hermann-Herder-Str. 3, 79104 Freiburg, Germany
| | - Rainer Haag
- Institute
of Chemistry and Biochemistry, Freie Universität
Berlin, Takustr. 3, 14195 Berlin, Germany
- Multifunctional
Biomaterials for Medicine, Helmholtz Virtual Institute, Kantstr. 55, 14513 Teltow-Seehof, Germany
| | - Matthias Ballauff
- Institute
of Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Multifunctional
Biomaterials for Medicine, Helmholtz Virtual Institute, Kantstr. 55, 14513 Teltow-Seehof, Germany
- Institut
für Physik, Humboldt-Universität
zu Berlin, Newtonstr. 15, 12489 Berlin, Germany
| |
Collapse
|
46
|
Biophysical, Biochemical, and Cell Based Approaches Used to Decipher the Role of Carbonic Anhydrases in Cancer and to Evaluate the Potency of Targeted Inhibitors. INTERNATIONAL JOURNAL OF MEDICINAL CHEMISTRY 2018; 2018:2906519. [PMID: 30112206 PMCID: PMC6077552 DOI: 10.1155/2018/2906519] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/01/2018] [Accepted: 06/25/2018] [Indexed: 12/12/2022]
Abstract
Carbonic anhydrases (CAs) are thought to be important for regulating pH in the tumor microenvironment. A few of the CA isoforms are upregulated in cancer cells, with only limited expression in normal cells. For these reasons, there is interest in developing inhibitors that target these tumor-associated CA isoforms, with increased efficacy but limited nonspecific cytotoxicity. Here we present some of the biophysical, biochemical, and cell based techniques and approaches that can be used to evaluate the potency of CA targeted inhibitors and decipher the role of CAs in tumorigenesis, cancer progression, and metastatic processes. These techniques include esterase activity assays, stop flow kinetics, and mass inlet mass spectroscopy (MIMS), all of which measure enzymatic activity of purified protein, in the presence or absence of inhibitors. Also discussed is the application of X-ray crystallography and Cryo-EM as well as other structure-based techniques and thermal shift assays to the studies of CA structure and function. Further, large-scale genomic and proteomic analytical methods, as well as cell based techniques like those that measure cell growth, apoptosis, clonogenicity, and cell migration and invasion, are discussed. We conclude by reviewing approaches that test the metastatic potential of CAs and how the aforementioned techniques have contributed to the field of CA cancer research.
Collapse
|
47
|
Johansson P, Kaspersson K, Gurrell IK, Bäck E, Eketjäll S, Scott CW, Cebers G, Thorne P, McKenzie MJ, Beaton H, Davey P, Kolmodin K, Holenz J, Duggan ME, Budd Haeberlein S, Bürli RW. Toward β-Secretase-1 Inhibitors with Improved Isoform Selectivity. J Med Chem 2018; 61:3491-3502. [DOI: 10.1021/acs.jmedchem.7b01716] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Patrik Johansson
- Discovery Sciences, IMED Biotech Unit, AstraZeneca, S-43183 Mölndal, Sweden
| | - Karin Kaspersson
- Discovery Sciences, IMED Biotech Unit, AstraZeneca, S-43183 Mölndal, Sweden
| | - Ian K. Gurrell
- Neuroscience, IMED Biotech Unit, AstraZeneca, Cambridge CB21 6GH, U.K
| | - Elisabeth Bäck
- Discovery Sciences, IMED Biotech Unit, AstraZeneca, S-43183 Mölndal, Sweden
| | - Susanna Eketjäll
- Cardiovascular and Metabolic Diseases, IMED Biotech Unit, AstraZeneca, 141 57 Huddinge, Sweden
| | - Clay W. Scott
- Discovery Safety, Drug Safety and Metabolism, IMED Biotech Unit, AstraZeneca, Waltham, Massachusetts 02451, United States
| | - Gvido Cebers
- Neuroscience, IMED Biotech Unit, AstraZeneca, Cambridge CB21 6GH, U.K
| | | | | | | | - Paul Davey
- Oncology Chemistry, IMED Biotech Unit, AstraZeneca, Cambridge CB4 0WG, U.K
| | | | - Jörg Holenz
- Neuroscience, IMED Biotech Unit, AstraZeneca, Cambridge CB21 6GH, U.K
| | - Mark E. Duggan
- Neuroscience, IMED Biotech Unit, AstraZeneca, Cambridge CB21 6GH, U.K
| | | | - Roland W. Bürli
- Neuroscience, IMED Biotech Unit, AstraZeneca, Cambridge CB21 6GH, U.K
| |
Collapse
|
48
|
Smirnov A, Zubrienė A, Manakova E, Gražulis S, Matulis D. Crystal structure correlations with the intrinsic thermodynamics of human carbonic anhydrase inhibitor binding. PeerJ 2018; 6:e4412. [PMID: 29503769 PMCID: PMC5831161 DOI: 10.7717/peerj.4412] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2017] [Accepted: 02/04/2018] [Indexed: 11/28/2022] Open
Abstract
The structure-thermodynamics correlation analysis was performed for a series of fluorine- and chlorine-substituted benzenesulfonamide inhibitors binding to several human carbonic anhydrase (CA) isoforms. The total of 24 crystal structures of 16 inhibitors bound to isoforms CA I, CA II, CA XII, and CA XIII provided the structural information of selective recognition between a compound and CA isoform. The binding thermodynamics of all structures was determined by the analysis of binding-linked protonation events, yielding the intrinsic parameters, i.e., the enthalpy, entropy, and Gibbs energy of binding. Inhibitor binding was compared within structurally similar pairs that differ by para- or meta-substituents enabling to obtain the contributing energies of ligand fragments. The pairs were divided into two groups. First, similar binders—the pairs that keep the same orientation of the benzene ring exhibited classical hydrophobic effect, a less exothermic enthalpy and a more favorable entropy upon addition of the hydrophobic fragments. Second, dissimilar binders—the pairs of binders that demonstrated altered positions of the benzene rings exhibited the non-classical hydrophobic effect, a more favorable enthalpy and variable entropy contribution. A deeper understanding of the energies contributing to the protein-ligand recognition should lead toward the eventual goal of rational drug design where chemical structures of ligands could be designed based on the target protein structure.
Collapse
Affiliation(s)
- Alexey Smirnov
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Vilnius University, Vilnius, Lithuania
| | - Asta Zubrienė
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Vilnius University, Vilnius, Lithuania
| | - Elena Manakova
- Department of Protein-DNA Interactions, Institute of Biotechnology, Vilnius University, Vilnius, Lithuania
| | - Saulius Gražulis
- Department of Protein-DNA Interactions, Institute of Biotechnology, Vilnius University, Vilnius, Lithuania
| | - Daumantas Matulis
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Vilnius University, Vilnius, Lithuania
| |
Collapse
|
49
|
Dissecting the Structure-Activity Relationship of Galectin-Ligand Interactions. Int J Mol Sci 2018; 19:ijms19020392. [PMID: 29382172 PMCID: PMC5855614 DOI: 10.3390/ijms19020392] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2017] [Revised: 01/18/2018] [Accepted: 01/24/2018] [Indexed: 02/08/2023] Open
Abstract
Galectins are β-galactoside-binding proteins. As carbohydrate-binding proteins, they participate in intracellular trafficking, cell adhesion, and cell-cell signaling. Accumulating evidence indicates that they play a pivotal role in numerous physiological and pathological activities, such as the regulation on cancer progression, inflammation, immune response, and bacterial and viral infections. Galectins have drawn much attention as targets for therapeutic interventions. Several molecules have been developed as galectin inhibitors. In particular, TD139, a thiodigalactoside derivative, is currently examined in clinical trials for the treatment of idiopathic pulmonary fibrosis. Herein, we provide an in-depth review on the development of galectin inhibitors, aiming at the dissection of the structure-activity relationship to demonstrate how inhibitors interact with galectin(s). We especially integrate the structural information established by X-ray crystallography with several biophysical methods to offer, not only in-depth understanding at the molecular level, but also insights to tackle the existing challenges.
Collapse
|
50
|
Kachooei E, Mozaffarian F, Khodagholi F, Sadeghi P, Karami L, Ghasemi A, Vahdat E, Saboury AA, Sheibani N, Moosavi-Movahedi AA. Paclitaxel inhibited lysozyme fibrillation by increasing colloidal stability through formation of "off-pathway" oligomers. Int J Biol Macromol 2018; 111:870-879. [PMID: 29352977 DOI: 10.1016/j.ijbiomac.2018.01.072] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Revised: 01/10/2018] [Accepted: 01/11/2018] [Indexed: 12/13/2022]
Abstract
Protein fibrillation is a challenging issue in medicine, causing many diseases, and an impediment to pharmaceutics and protein industry. Many chemicals, especially polyphenol compounds and aromatic small molecules, have been widely used as an effective strategy to combat protein fibril formation. Hence, understanding mechanisms of fibrillation inhibition and contributing forces in this process are significant. In this study, the inhibitory effect of paclitaxel on lysozyme fibrillation was investigated with respect to thermal and colloidal stability. Fibrillation was monitored with ThT fluorescence, circular dichroism, and AFM; paclitaxel-lysozyme interaction with isothermal titration calorimetry and docking; thermal and colloidal stability with differential scanning calorimetry and zeta-pulse, respectively. Paclitaxel inhibited lysozyme fibrillation, and interacted with lysozyme through hydrogen bonds and van der Waals' interactions. The viability of PC12 cells retrieved as a result of fibrillation inhibition by paclitaxel. Hydrophobic forces dominantly shielded the aggregation-prone region of lysozyme and suppressed the effective interactions between lysozyme monomers. Although paclitaxel did not affect lysozyme's thermal stability, it increased lysozyme's colloidal stability by either increasing the surface charge density or charge distribution on lysozyme. In conclusion, our results suggest a model for paclitaxel's inhibitory role through two complementary steps driving to "off-pathway" oligomer formation and attenuation of fibril formation.
Collapse
Affiliation(s)
- Ehsan Kachooei
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Faroogh Mozaffarian
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Fariba Khodagholi
- Neuroscience Research Center, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Payam Sadeghi
- Neuroscience Research Center, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Leila Karami
- Department of Cell and Molecular Biology, Faculty of Biological Sciences, Kharazmi University, Tehran, Iran
| | - Atiyeh Ghasemi
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Elham Vahdat
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Ali Akbar Saboury
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran; Center of Excellence in Biothermodynamics, University of Tehran, Tehran, Iran
| | - Nader Sheibani
- Departments of Ophthalmology, University of Wisconsin School of Medicine and Public Health, Madison, WI, USA
| | - Ali Akbar Moosavi-Movahedi
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran; Center of Excellence in Biothermodynamics, University of Tehran, Tehran, Iran.
| |
Collapse
|