1
|
Spinello A, Borišek J, Pavlin M, Janoš P, Magistrato A. Computing Metal-Binding Proteins for Therapeutic Benefit. ChemMedChem 2021; 16:2034-2049. [PMID: 33740297 DOI: 10.1002/cmdc.202100109] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Indexed: 01/18/2023]
Abstract
Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.
Collapse
Affiliation(s)
- Angelo Spinello
- National Research Council of Italy (CNR)-, Institute of Materials (IOM) c/o International School for Advanced Studies (SISSA), via Bonomea 265, 34136, Trieste, Italy
| | - Jure Borišek
- National Institute of Chemistry Institution Hajdrihova ulica 19, 1000, Ljubljana, Slovenia
| | - Matic Pavlin
- Laboratory of Microsensor Structures and Electronics Faculty of Electrical Engineering, University of Ljubljana Tržaška cesta 25, 1000, Ljubljana, Slovenia
| | - Pavel Janoš
- National Research Council of Italy (CNR)-, Institute of Materials (IOM) c/o International School for Advanced Studies (SISSA), via Bonomea 265, 34136, Trieste, Italy
| | - Alessandra Magistrato
- National Research Council of Italy (CNR)-, Institute of Materials (IOM) c/o International School for Advanced Studies (SISSA), via Bonomea 265, 34136, Trieste, Italy
| |
Collapse
|
2
|
Crosby DC, Lei X, Gibbs CG, Reinecke MG, Robinson WE. Mutagenesis of Lysines 156 and 159 in Human Immunodeficiency Virus Type 1 Integrase (IN) Reveals Differential Interactions between these Residues and Different IN Inhibitors. Nat Prod Commun 2015. [DOI: 10.1177/1934578x1501000129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Human immunodeficiency virus (HIV) type 1 integrase (IN) active site, and viral DNA-binding residues K156 and K159 are predicted to interact both with strand transfer-selective IN inhibitors (STI), e.g. L-731,988, Elvitegravir (EVG), and the FDA-approved IN inhibitor, Raltegravir (RGV), and strand transfer non-selective inhibitors, e.g. dicaffeoyltartaric acids (DCTAs), e.g. L-chicoric acid (L-CA). To test posited roles for these two lysine residues in inhibitor action we assayed the potency of L-CA and several STI against a panel of K156 and K159 mutants. Mutagenesis of K156 conferred resistance to L-CA and mutagenesis of either K156 or K159 conferred resistance to STI indicating that the cationic charge at these two viral DNA-binding residues is important for inhibitor potency. IN K156N, a reported polymorphism associated with resistance to RGV, conferred resistance to L-CA and STI as well. To investigate the apparent preference L-CA exhibits for interactions with K156, we assayed the potency of several hybrid inhibitors containing combinations of DCTA and STI pharmacophores against recombinant IN K156A or K159A. Although K156A conferred resistance to diketo acid-branched bis-catechol hybrid inhibitors, neither K156A nor K159A conferred resistance to their monocatechol counterparts, suggesting that bis-catechol moieties direct DCTAs toward K156. In contrast, STI were more promiscuous in their interaction with K156 and K159. Taken together, the results of this study indicate that DCTAs interact with IN in a manner different than that of STI and suggest that DCTAs are an attractive candidate chemotype for development into drugs potent against STI-resistant IN.
Collapse
Affiliation(s)
- David C. Crosby
- Department of Pathology and Laboratory Medicine, University of California, Irvine, CA 92697-4800, USA
| | - Xiangyang Lei
- Department of Chemistry, Texas Christian University, Fort Worth, TX 76129, USA
| | - Charles G. Gibbs
- Department of Chemistry, Texas Christian University, Fort Worth, TX 76129, USA
| | - Manfred G. Reinecke
- Department of Chemistry, Texas Christian University, Fort Worth, TX 76129, USA
| | - W. Edward Robinson
- Department of Pathology and Laboratory Medicine, University of California, Irvine, CA 92697-4800, USA
- Department of Microbiology and Molecular Genetics, University of California, Irvine, CA 92697-4800, USA
| |
Collapse
|
3
|
|
4
|
Huang M, Grant GH, Richards WG. Binding modes of diketo-acid inhibitors of HIV-1 integrase: a comparative molecular dynamics simulation study. J Mol Graph Model 2011; 29:956-64. [PMID: 21531158 PMCID: PMC3101338 DOI: 10.1016/j.jmgm.2011.04.002] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2010] [Revised: 03/30/2011] [Accepted: 04/01/2011] [Indexed: 11/18/2022]
Abstract
HIV-1 integrase (IN) has become an attractive target since drug resistance against HIV-1 reverse transcriptase (RT) and protease (PR) has appeared. Diketo acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN: however the action mechanism is not well understood. Here, to study the inhibition mechanism of DKAs we performed 10 ns comparative molecular dynamics simulations on HIV-1 IN bound with three most representative DKA inhibitors: Shionogi inhibitor, S-1360 and two Merck inhibitors L-731,988 and L-708,906. Our simulations show that the acidic part of S-1360 formed salt bridge and cation-π interactions with Lys159. In addition, the catalytic Glu152 in S-1360 was pushed away from the active site to form an ion-pair interaction with Arg199. The Merck inhibitors can maintain either one or both of these ion-pair interaction features. The difference in potencies of the DKA inhibitors is thus attributed to the different binding modes at the catalytic site. Such structural information at atomic level, not only demonstrates the action modes of DKA inhibitors but also provides a novel starting point for structural-based design of HIV-1 IN inhibitors.
Collapse
Affiliation(s)
- Meilan Huang
- School of Chemistry and Chemical Engineering, David Keir Building, Queens University Belfast, Stranmillis Road, Belfast BT95AG, UK.
| | | | | |
Collapse
|
5
|
Abstract
Integration of the HIV-1 viral DNA generated by reverse transcription of the RNA genome into the host cell chromosomes is a key step of viral replication, catalyzed by the viral integrase. In October 2007, the first integrase inhibitor, raltegravir, was approved for clinical use under the name of Isentress™. The results of the various clinical trials that have evaluated raltegravir have been very encouraging with regard to the immunological and virological efficacy and tolerance. However, as observed for other anti-retrovirals, specific resistance mutations have been identified in patients failing to respond to treatment with raltegravir. Although knowledge of the integrase structural biology remains fragmentary, the structures and modeling data available might provide relevant clues on the origin of the emergence of these resistance mutations. In this review, we describe the mechanism of action of this drug and the main data relating to its use in vivo, together with recent structural data important to our understanding of the origin of viral resistance.
Collapse
Affiliation(s)
- Jean-Francois Mouscadet
- LBPA, CNRS UMR8113, Ecole Normale Superieure de Cachan, 61 avenue du President Wilson, 94235 Cachan Cedex, France.
| | | |
Collapse
|
6
|
Mouscadet JF, Arora R, André J, Lambry JC, Delelis O, Malet I, Marcelin AG, Calvez V, Tchertanov L. HIV-1 IN alternative molecular recognition of DNA induced by raltegravir resistance mutations. J Mol Recognit 2010; 22:480-94. [PMID: 19623602 DOI: 10.1002/jmr.970] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Virologic failure during treatment with raltegravir, the first effective drug targeting HIV integrase, is associated with two exclusive pathways involving either Q148H/R/K, G140S/A or N155H mutations. We carried out a detailed analysis of the molecular and structural effects of these mutations. We observed no topological change in the integrase core domain, with conservation of a newly identified Omega-shaped hairpin containing the Q148 residue, in particular. In contrast, the mutations greatly altered the specificity of DNA recognition by integrase. The native residues displayed a clear preference for adenine, whereas the mutant residues strongly favored pyrimidines. Raltegravir may bind to N155 and/or Q148 residues as an adenine bioisoster. This may account for the selected mutations impairing raltegravir binding while allowing alternative DNA recognition by integrase. This study opens up new opportunities for the design of integrase inhibitors active against raltegravir-resistant viruses.
Collapse
Affiliation(s)
- Jean-François Mouscadet
- LBPA, CNRS, Ecole Normale Supérieure de Cachan, 61 Avenue du Président Wilson, 94235 Cachan, France
| | | | | | | | | | | | | | | | | |
Collapse
|
7
|
Ranaghan KE, Mulholland AJ. Investigations of enzyme-catalysed reactions with combined quantum mechanics/molecular mechanics (QM/MM) methods. INT REV PHYS CHEM 2010. [DOI: 10.1080/01442350903495417] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
8
|
Rungrotmongkol T, Decha P, Sompornpisut P, Malaisree M, Intharathep P, Nunthaboot N, Udommaneethanakit T, Aruksakunwong O, Hannongbua S. Combined QM/MM mechanistic study of the acylation process in furin complexed with the H5N1 avian influenza virus hemagglutinin's cleavage site. Proteins 2009; 76:62-71. [PMID: 19089976 DOI: 10.1002/prot.22318] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Combined quantum mechanical/molecular mechanical (QM/MM) techniques have been applied to investigate the detailed reaction mechanism of the first step of the acylation process by furin in which the cleavage site of the highly pathogenic avian influenza virus subtype H5N1 (HPH5) acts as its substrate. The energy profile shows a simultaneous mechanism, known as a concerted reaction, of the two subprocesses: the proton transfer from Ser368 to His194 and the nucleophilic attack on the carbonyl carbon of the scissile peptide of the HPH5 cleavage site with a formation of tetrahedral intermediate (INT). The calculated energy barrier for this reaction is 16.2 kcal.mol(-1) at QM/MM B3LYP/6-31+G*//PM3-CHARMM22 level of theory. Once the reaction proceeds, the ordering of the electrostatic stabilization by protein environment is of the enzyme-substrate < transition state < INT complexes. Asp153 was found to play the most important role in the enzymatic reaction by providing the highest degree of intermediate complex stabilization. In addition, the negatively charged carbonyl oxygen of INT is well stabilized by the oxyanion hole constructed by Asn295's carboxamide and Ser368's backbone.
Collapse
Affiliation(s)
- Thanyada Rungrotmongkol
- Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
| | | | | | | | | | | | | | | | | |
Collapse
|
9
|
Abstract
Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.
Collapse
Affiliation(s)
- Hans Martin Senn
- Department of Chemistry, WestCHEM and University of Glasgow, Glasgow G12 8QQ, UK.
| | | |
Collapse
|
10
|
|