1
|
Dutkiewicz Z. Computational methods for calculation of protein-ligand binding affinities in structure-based drug design. PHYSICAL SCIENCES REVIEWS 2022. [DOI: 10.1515/psr-2020-0034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Abstract
Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.
Collapse
Affiliation(s)
- Zbigniew Dutkiewicz
- Department of Chemical Technology of Drugs , Poznan University of Medical Sciences , ul. Grunwaldzka 6 , 60-780 Poznań , Poznan , 60-780, Poland
| |
Collapse
|
2
|
Rocha-Santos A, Chaves EJ, Grillo IB, de Freitas AS, Araújo DAM, Rocha GB. Thermochemical and Quantum Descriptor Calculations for Gaining Insight into Ricin Toxin A (RTA) Inhibitors. ACS OMEGA 2021; 6:8764-8777. [PMID: 33842748 PMCID: PMC8027999 DOI: 10.1021/acsomega.0c02588] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Accepted: 12/30/2020] [Indexed: 05/03/2023]
Abstract
In this work, we performed a study to assess the interactions between the ricin toxin A (RTA) subunit of ricin and some of its inhibitors using modern semiempirical quantum chemistry and ONIOM quantum mechanics/molecular mechanics (QM/MM) methods. Two approaches were followed (calculation of binding enthalpies, ΔH bind, and reactivity quantum chemical descriptors) and compared with the respective half-maximal inhibitory concentration (IC50) experimental data, to gain insight into RTA inhibitors and verify which quantum chemical method would better describe RTA-ligand interactions. The geometries for all RTA-ligand complexes were obtained after running classical molecular dynamics simulations in aqueous media. We found that single-point energy calculations of ΔH bind with the PM6-DH+, PM6-D3H4, and PM7 semiempirical methods and ONIOM QM/MM presented a good correlation with the IC50 data. We also observed, however, that the correlation decreased significantly when we calculated ΔH bind after full-atom geometry optimization with all semiempirical methods. Based on the results from reactivity descriptors calculations for the cases studied, we noted that both types of interactions, molecular overlap and electrostatic interactions, play significant roles in the overall affinity of these ligands for the RTA binding pocket.
Collapse
Affiliation(s)
- Acassio Rocha-Santos
- Department
of Chemistry, Federal University of Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
| | - Elton José
Ferreira Chaves
- Department
of Biotechnology, Federal University of
Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
| | - Igor Barden Grillo
- Department
of Chemistry, Federal University of Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
| | - Amanara Souza de Freitas
- Department
of Chemical Engineering, Federal University
of Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
| | | | - Gerd Bruno Rocha
- Department
of Chemistry, Federal University of Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
- . Phone/Fax: +55-83-3216-7437
| |
Collapse
|
3
|
Pecina A, Eyrilmez SM, Köprülüoğlu C, Miriyala VM, Lepšík M, Fanfrlík J, Řezáč J, Hobza P. SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design. Chempluschem 2020; 85:2362-2371. [PMID: 32609421 DOI: 10.1002/cplu.202000120] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Revised: 05/27/2020] [Indexed: 12/17/2022]
Abstract
Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.
Collapse
Affiliation(s)
- Adam Pecina
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Saltuk M Eyrilmez
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Cemal Köprülüoğlu
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Vijay Madhav Miriyala
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| |
Collapse
|
4
|
Ryde U, Söderhjelm P. Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. Chem Rev 2016; 116:5520-66. [DOI: 10.1021/acs.chemrev.5b00630] [Citation(s) in RCA: 175] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Pär Söderhjelm
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
5
|
Tandon V, Urvashi, Yadav P, Sur S, Abbat S, Tiwari V, Hewer R, Papathanasopoulos MA, Raja R, Banerjea AC, Verma AK, Kukreti S, Bharatam PV. Design, Synthesis, and Biological Evaluation of 1,2-Dihydroisoquinolines as HIV-1 Integrase Inhibitors. ACS Med Chem Lett 2015; 6:1065-70. [PMID: 26487913 DOI: 10.1021/acsmedchemlett.5b00230] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2015] [Accepted: 08/09/2015] [Indexed: 12/30/2022] Open
Abstract
6-Endo-dig-cyclization is an efficient method for the synthesis of 1,2-dihydroisoquinolines. We have synthesized few 1,2-dihydroisoquinolines having different functionality at the C-1, C-3, C-7, and N-2 positions for evaluation against HIV-1 integrase (HIV1-IN) inhibitory activity. A direct nitro-Mannich condensation of o-alkynylaldimines and dual activation of o-alkynyl aldehydes by inexpensive cobalt chloride yielded desired compounds. Out of 24 compounds, 4m and 6c came out as potent integrase inhibitors in in vitro strand transfer (ST) assay, with IC50 value of 0.7 and 0.8 μM, respectively. Molecular docking of these compounds in integrase revealed strong interaction between metal and ligands, which stabilizes the enzyme-inhibitor complex. The ten most active compounds were subjected to antiviral assay. Out of those, 6c reduced the level of p24 viral antigen by 91%, which is comparable to RAL in antiviral assay. Interestingly, these compounds showed similar ST inhibitory activity in G140S mutant, suggesting they can act against resistant strains.
Collapse
Affiliation(s)
- Vibha Tandon
- Department
of Chemistry, University of Delhi, Delhi 110007, India
- Special
Centre for Molecular Medicine, Jawaharlal Nehru University, New Delhi 110067, India
| | - Urvashi
- Department
of Chemistry, University of Delhi, Delhi 110007, India
| | - Pooja Yadav
- Department
of Chemistry, University of Delhi, Delhi 110007, India
| | - Souvik Sur
- Department
of Chemistry, University of Delhi, Delhi 110007, India
| | - Sheenu Abbat
- National Institute of Pharmaceutical Education and Research, S. A. S Nagar, Mohali, Punjab 160062, India
| | - Vinod Tiwari
- Department
of Chemistry, University of Delhi, Delhi 110007, India
| | - Raymond Hewer
- Biomedical
Advanced Material Division, Mintek, Private Bag X3015, Randburg 2125, Johannesburg, South Africa
| | - Maria A. Papathanasopoulos
- Department
of Molecular Medicine and Haematology, University of the Witwatersrand Medical School,
Parktown 2193, Johannesburg, South Africa
| | - Rameez Raja
- Laboratory
of Virology, National Institute of Immunology, New Delhi 110067, India
| | - Akhil C. Banerjea
- Laboratory
of Virology, National Institute of Immunology, New Delhi 110067, India
| | | | - Shrikant Kukreti
- Department
of Chemistry, University of Delhi, Delhi 110007, India
| | - Prasad V. Bharatam
- National Institute of Pharmaceutical Education and Research, S. A. S Nagar, Mohali, Punjab 160062, India
| |
Collapse
|
6
|
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
|
7
|
García-Sosa AT, Maran U. Improving the use of ranking in virtual screening against HIV-1 integrase with triangular numbers and including ligand profiling with antitargets. J Chem Inf Model 2014; 54:3172-85. [PMID: 25303089 DOI: 10.1021/ci500300u] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
A delicate balance exists between a drug molecule's toxicity and its activity. Indeed, efficacy, toxicity, and side effect problems are a common cause for the termination of drug candidate compounds and development projects. To address this, an antitarget interaction profile is built and combined with virtual screening and cross docking for new inhibitors of HIV-1 integrase, in order to consider possible off-target interactions as early as possible in a drug or hit discovery program. New ranking techniques using triangular numbers improve ranking information on the compounds and recovery of known inhibitors into the top compounds using different docking programs. This improved ranking arises from using consensus of ranks between docking programs and ligand efficiencies to derive a new rank, instead of using absolute score values, or average of ranks. The triangular number rerank also allowed the objective combination of results from several protein targets or screen conditions and several programs. Triangular number reranking conserves more information than other reranking methods such as average of scores or averages of ranks. In addition, the use of triangular numbers for reranking makes possible the use of thresholds with a justified leeway based on the number of available known inhibitors, so that the majority of the compounds above the threshold in ranks compare to the compounds that have known experimentally determined biological activity. The battery of anti- or off-targets can be tailored to specific molecular or drug design challenges. In silico filters can thus be deployed in successive stages, for prefiltering, activity profiling, and for further analysis and triaging of libraries of compounds.
Collapse
|
8
|
Beesu M, Malladi SS, Fox LM, Jones CD, Dixit A, David SA. Human Toll-like receptor 8-selective agonistic activities in 1-alkyl-1H-benzimidazol-2-amines. J Med Chem 2014; 57:7325-41. [PMID: 25102141 PMCID: PMC4161153 DOI: 10.1021/jm500701q] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
![]()
Toll-like
receptor (TLR)-8 agonists strongly induce the production
of T helper 1-polarizing cytokines and may therefore serve as promising
candidate vaccine adjuvants, especially for the very young and the
elderly. Earlier structure-based ligand design led to the identification
of 3-pentyl-quinoline-2-amine as a novel, human TLR8-specific agonist.
Comprehensive structure–activity relationships in ring-contracted
1-alkyl-1H-benzimidazol-2-amines were undertaken,
and the best-in-class compound, 4-methyl-1-pentyl-1H-benzo[d]imidazol-2-amine, was found to be a pure
TLR8 agonist, evoking strong proinflammatory cytokine and Type II
interferon responses in human PBMCs, with no attendant CD69 upregulation
in natural lymphocytic subsets. The 1-alkyl-1H-benzimidazol-2-amines
represent a novel, alternate chemotype with pure TLR8-agonistic activities
and will likely prove useful not only in understanding TLR8 signaling
but also perhaps as a candidate vaccine adjuvant.
Collapse
Affiliation(s)
- Mallesh Beesu
- Department of Medicinal Chemistry, University of Kansas , Lawrence, Kansas 66047, United States
| | | | | | | | | | | |
Collapse
|
9
|
Reis M, Alves CN, Lameira J, Tuñón I, Martí S, Moliner V. The catalytic mechanism of glyceraldehyde 3-phosphate dehydrogenase from Trypanosoma cruzi elucidated via the QM/MM approach. Phys Chem Chem Phys 2013; 15:3772-85. [PMID: 23389436 DOI: 10.1039/c3cp43968b] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) has been identified as a key enzyme involved in glycolysis processes for energy production in the Trypanosoma cruzi parasite. This enzyme catalyses the oxidative phosphorylation of glyceraldehyde 3-phosphate (G3P) in the presence of inorganic phosphate (Pi) and nicotinamide adenosine dinucleotide (NAD+). The catalytic mechanism used by GAPDH has been intensively investigated. However, the individual roles of Pi and the C3 phosphate of G3P (Ps) sites, as well as some residues such as His194 in the catalytic mechanism, remain unclear. In this study, we have employed Molecular Dynamics (MD) simulations within hybrid quantum mechanical/molecular mechanical (QM/MM) potentials to obtain the Potential of Mean Force of the catalytic oxidative phosphorylation mechanism of the G3P substrate used by GAPDH. According to our results, the first stage of the reaction (oxidoreduction) takes place in the Pi site (energetically more favourable), with the formation of oxyanion thiohemiacetal and thioacylenzyme intermediates without acid-base assistance of His194. Analysis of the interaction energy by residues shows that Arg249 has an important role in the ability of the enzyme to bind the G3P substrate, which interacts with NAD+ and other important residues, such as Cys166, Glu109, Thr167, Ser247 and Thr226, in the GAPDH active site. Finally, the inhibition mechanism of the GAPDH enzyme by the 3-(p-nitrophenoxycarboxyl)-3-ethylene propyl dihydroxyphosphonate inhibitor was investigated in order to contribute to the design of new inhibitors of GAPDH from Trypanosoma cruzi.
Collapse
Affiliation(s)
- Mauro Reis
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, 66075-110, Belém, PA, Brazil
| | | | | | | | | | | |
Collapse
|
10
|
Kokatla HP, Sil D, Malladi SS, Balakrishna R, Hermanson AR, Fox LM, Wang X, Dixit A, David SA. Exquisite selectivity for human toll-like receptor 8 in substituted furo[2,3-c]quinolines. J Med Chem 2013; 56:6871-85. [PMID: 23899291 DOI: 10.1021/jm400694d] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Toll-like receptor (TLR)-8 agonists activate adaptive immune responses by inducing robust production of T helper 1-polarizing cytokines, suggesting that TLR8-active compounds may be promising candidate adjuvants. We synthesized and evaluated hitherto unexplored furo[2,3-c]quinolines and regioisomeric furo[3,2-c]quinolines derived via a tandem, one-pot Sonogashira coupling and intramolecular 5-endo-dig cyclization strategy in a panel of primary screens. We observed a pure TLR8-agonistic activity profile in select furo[2,3-c]quinolines, with maximal potency conferred by a C2-butyl group (EC50 = 1.6 μM); shorter, longer, or substituted homologues as well as compounds bearing C1 substitutions were inactive, which was rationalized by docking studies using the recently described crystal structure of human TLR8. The best-in-class compound displayed prominent proinflammatory cytokine induction (including interleukin-12 and interleukin-18), but was bereft of interferon-α inducing properties, confirming its high selectivity for human TLR8.
Collapse
Affiliation(s)
- Hari Prasad Kokatla
- Department of Medicinal Chemistry and ‡Genomics Facility, University of Kansas , Lawrence, Kansas 66047, United States
| | | | | | | | | | | | | | | | | |
Collapse
|
11
|
In silico design: extended molecular dynamic simulations of a new series of dually acting inhibitors against EGFR and HER2. J Mol Graph Model 2013; 44:220-31. [PMID: 23911931 DOI: 10.1016/j.jmgm.2013.06.004] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2013] [Revised: 06/13/2013] [Accepted: 06/17/2013] [Indexed: 01/23/2023]
Abstract
Based on the hit structures that have been identified in our previous studies against EGFR and HER2, new potential inhibitors that share the same scaffold of the hit structures are designed and screened in silico. Insights into understanding the potential inhibitory effect of the new inhibitors against both EGFR and HER2 receptors is obtained using extended molecular dynamics (MD) simulations and different scoring techniques. The binding mechanisms and dynamics are detailed with respect to two approved inhibitors against EGFR (lapatinib) and HER2 (SYR127063). The best scoring inhibitor (T9) is chosen for additional in silico investigation against both the wild-type and T790M mutant strain of EGFR and the wild-type HER2. The results reveal that certain substitution patterns increase the stability and assure stronger binding and higher H-bond occupancy of the conserved water molecule that is commonly observed with kinase crystal structures. Furthermore, the new inhibitor (T9) forms stable interactions with the mutant strain as a direct consequence of the enhanced ability to form additional hydrogen bonding interactions with binding site residues.
Collapse
|
12
|
Srivastava HK, Sastry GN. Molecular dynamics investigation on a series of HIV protease inhibitors: assessing the performance of MM-PBSA and MM-GBSA approaches. J Chem Inf Model 2012; 52:3088-98. [PMID: 23121465 DOI: 10.1021/ci300385h] [Citation(s) in RCA: 91] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
The binding free energies (ΔG(Bind)) obtained from molecular mechanics with Poisson-Boltzmann surface area (MM-PBSA) or molecular mechanics with Generalized Born surface area (MM-GBSA) calculations using molecular dynamics (MD) trajectories are the most popular procedures to measure the strength of interactions between a ligand and its receptor. Several attempts have been made to correlate the ΔG(Bind) and experimental IC(50) values in order to observe the relationship between binding strength of a ligand (with its receptor) and its inhibitory activity. The duration of MD simulations seems very important for getting acceptable correlation. Here, we are presenting a systematic study to estimate the reasonable MD simulation time for acceptable correlation between ΔG(Bind) and experimental IC(50) values. A comparison between MM-PBSA and MM-GBSA approaches is also presented at various time scales. MD simulations (10 ns) for 14 HIV protease inhibitors have been carried out by using the Amber program. MM-PBSA/GBSA based ΔG(Bind) have been calculated and correlated with experimental IC(50) values at different time scales (0-1 to 0-10 ns). This study clearly demonstrates that the MM-PBSA based ΔG(Bind) (ΔG(Bind)-PB) values provide very good correlation with experimental IC(50) values (quantitative and qualitative) when MD simulation is carried out for a longer time; however, MM-GBSA based ΔG(Bind) (ΔG(Bind)-GB) values show acceptable correlation for shorter time of simulation also. The accuracy of ΔG(Bind)-PB increases and ΔG(Bind)-GB remains almost constant with the increasing time of simulation.
Collapse
Affiliation(s)
- Hemant Kumar Srivastava
- Centre for Molecular Modelling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad 500 607, India
| | | |
Collapse
|
13
|
Lima AH, Souza PRM, Alencar N, Lameira J, Govender T, Kruger HG, Maguire GEM, Alves CN. Molecular modeling of T. rangeli, T. brucei gambiense, and T. evansi sialidases in complex with the DANA inhibitor. Chem Biol Drug Des 2012; 80:114-20. [PMID: 22416952 DOI: 10.1111/j.1747-0285.2012.01380.x] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Trypanosomal (trans-) sialidases are enzymes that catalyze the transfer of sialic acid residues between host and parasite glycoconjugates. Herein, we have used homology modeling to construct the 3D structures of sialidases from Trypanosoma brucei and Trypanosoma evansi. Hybrid quantum mechanical/molecular mechanical molecular dynamics simulations were used to determine the interaction energy between the 2-deoxy-2,3-didehydro-N-acetylneuraminic acid inhibitor and the three sialidases studied here. Our results suggest that the two constructed enzymes share the same basic fold motive of the Trypanosoma rangeli crystallographic structure. In addition, quantum mechanical/molecular mechanical molecular dynamics simulations show that the 2-deoxy-2,3-didehydro-N-acetylneuraminic acid inhibitor forms a stronger complex with Trypanosoma rangeli than with Trypanosoma brucei and Trypanosoma evansi sialidases. Finally, the interaction energy by residues shows that the arginine triad plays a decisive role to complex 2-deoxy-2,3-didehydro-N-acetylneuraminic acid with the enzyme through hydrogen bonding.
Collapse
Affiliation(s)
- Anderson H Lima
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, 66075-110 Belém, PA, Brazil
| | | | | | | | | | | | | | | |
Collapse
|
14
|
Świderek K, Martí S, Moliner V. Theoretical studies of HIV-1 reverse transcriptase inhibition. Phys Chem Chem Phys 2012; 14:12614-24. [DOI: 10.1039/c2cp40953d] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
15
|
Xue W, Liu H, Yao X. Molecular mechanism of HIV-1 integrase-vDNA interactions and strand transfer inhibitor action: A molecular modeling perspective. J Comput Chem 2011; 33:527-36. [DOI: 10.1002/jcc.22887] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2011] [Revised: 09/25/2011] [Accepted: 10/20/2011] [Indexed: 01/03/2023]
|
16
|
Carneiro AS, Lameira J, Alves CN. A theoretical study of the molecular mechanism of the GAPDH Trypanosoma cruzi enzyme involving iodoacetate inhibitor. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.08.051] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
17
|
|
18
|
Protein–ligand interaction of T. cruzi trans-sialidase inhibitors: a docking and QM/MM MD study. Struct Chem 2011. [DOI: 10.1007/s11224-011-9854-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
19
|
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
|
20
|
Silva NDF, Lameira J, Alves CN. A quantum mechanical/molecular mechanical study of the aspartic protease plasmepsin IV complexed with allophenylnorstatine-based inhibitor. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.04.085] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
21
|
Lameira J, Alves CN, Tuñón I, Martí S, Moliner V. Enzyme Molecular Mechanism as a Starting Point to Design New Inhibitors: A Theoretical Study of O-GlcNAcase. J Phys Chem B 2011; 115:6764-75. [DOI: 10.1021/jp202079e] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jerônimo Lameira
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, 66075-110 Belém, PA, Brazil
| | - Cláudio Nahum Alves
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, 66075-110 Belém, PA, Brazil
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, València, Spain
| | - Sergio Martí
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| |
Collapse
|
22
|
Computational analysis of aspartic protease plasmepsin II complexed with EH58 inhibitor: a QM/MM MD study. J Mol Model 2011; 17:2631-8. [PMID: 21264482 DOI: 10.1007/s00894-011-0963-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2010] [Accepted: 01/05/2011] [Indexed: 10/18/2022]
Abstract
Plasmepsin (PM) II is one of four enzymes in the food vacuole of Plasmodium falciparum. It has become an attractive target for combating malaria through research regarding its importance in the P. falciparum metabolism and life cycle, making it the target of choice for structure-based drug design. This paper reports the results of hybrid quantum mechanics / molecular mechanics (QM/MM) molecular dynamics (MD) simulations employed to study the details of the interactions established between PM II and N-(3-{(2-benzo[1, 3]dioxol-5-yl-ethyl)[3-(1-methyl-3-oxo-1,3-dihydro-isoindol-2-yl) propionyl]-amino}-1-benzyl-2-(hydroxyl-propyl)-4-benzyloxy-3,5dimethoxy-benzamide (EH58), a well-known potent inhibitor for this enzyme. Electrostatic binding free energy and energy terms decomposition have been computed for PM II complexed with the EH58 inhibitor. The results reveal that there is a strong interaction between Asp34, Val78, Ser79, Tyr192 and Asp214 residues and the EH58 inhibitor. In addition, we have computed the potential of the mean force (PMF) profile in order to assign the protonation state of the two catalytic aspartates in PM II-EH58 complex. The results indicate that the protonation of Asp214 favors a stable active site structure, which is consistent with our electrostatic binding free energy calculation and with previous published works.
Collapse
|
23
|
Crosby DC, Lei X, Gibbs CG, McDougall BR, Robinson WE, Reinecke MG. Design, synthesis, and biological evaluation of novel hybrid dicaffeoyltartaric/diketo acid and tetrazole-substituted L-chicoric acid analogue inhibitors of human immunodeficiency virus type 1 integrase. J Med Chem 2010; 53:8161-75. [PMID: 20977258 DOI: 10.1021/jm1010594] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Fourteen analogues of the anti-HIV-1 integrase (IN) inhibitor L-chicoric acid (L-CA) were prepared. Their IC(50) values for 3'-end processing and strand transfer against recombinant HIV-1 IN were determined in vitro, and their cell toxicities and EC(50) against HIV-1 were measured in cells (ex vivo). Compounds 1-6 are catechol/β-diketoacid hybrids, the majority of which exhibit submicromolar potency against 3'-end processing and strand transfer, though only with modest antiviral activities. Compounds 7-10 are L-CA/p-fluorobenzylpyrroloyl hybrids, several of which were more potent against strand transfer than 3'-end processing, a phenomenon previously attributed to the β-diketo acid pharmacophore. Compounds 11-14 are tetrazole bioisosteres of L-CA and its analogues, whose in vitro potencies were comparable to L-CA but with enhanced antiviral potency. The trihydroxyphenyl analogue 14 was 30-fold more potent than L-CA at relatively nontoxic concentrations. These data indicate that L-CA analogues are attractive candidates for development into clinically relevant inhibitors of HIV-1 IN.
Collapse
Affiliation(s)
- David C Crosby
- Department of Pathology and Laboratory Medicine, University of California, Irvine, California 92697-4800, USA
| | | | | | | | | | | |
Collapse
|
24
|
Lameira J, Alves CN, Moliner V, Martí S, Castillo R, Tuñón I. Quantum mechanical/molecular mechanical molecular dynamics simulation of wild-type and seven mutants of CpNagJ in complex with PUGNAc. J Phys Chem B 2010; 114:7029-36. [PMID: 20429600 DOI: 10.1021/jp9115673] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The enzyme O-glycoprotein 2-acetamino-2-deoxy-beta-d-glucopyranosidase (O-GlcNAcase) is responsible for the removal of N-acetylglucosamine moieties from 2-acetamido-2-deoxy-beta-D-glucopyranose (O-GlcNAc) residues of serine/threonine residues of modified proteins. We herein present results of hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations applied to the study of the interactions established between a bacterial Clostridium perfringens homologue (CpNagJ) and PUGNAc, a potent known inhibitor of this enzyme. Electrostatic binding free energy and energy term decomposition have been computed for the wild-type CpNagJ and several mutants: D297N, D298N, Y335F, N390A, N396A, D401A, and W490A. The theoretical results have been compared with recently experimental data based on crystallographic and mutation studies on the same system. Our results reveal that, first, there is a strong interaction between Asp401, Asp298, and Asp297 residues and the PUGNAc inhibitor; and, second, the electrostatic substrate binding free energy is higher in wild-type, N390A and W490A mutants than in D297N, D298N, Y335F, N396A, and D401A ones, in accordance with the experimental results. Finally, both our theoretical predictions and the experimental data are compatible with a substrate-assisted reaction mechanism, involving two conserved aspartate residues.
Collapse
Affiliation(s)
- Jeronimo Lameira
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, 66075-110, Belém, PA, Brazil
| | | | | | | | | | | |
Collapse
|
25
|
Lu Y, Mei Y, Zhang JZH, Zhang D. Communications: Electron polarization critically stabilizes the Mg2+ complex in the catalytic core domain of HIV-1 integrase. J Chem Phys 2010; 132:131101. [PMID: 20387913 DOI: 10.1063/1.3360769] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we present a detailed dynamics study of the catalytic core domain (CCD) of HIV-1 integrase using both polarized and nonpolarized force fields. The numerical results reveal the critical role of protein polarization in stabilizing Mg(2+) coordination complex in CCD. Specifically, when nonpolarized force field is used, a remarkable drift of the Mg(2+) complex away from its equilibrium position is observed, which causes the binding site blocked by the Mg(2+) complex. In contrast, when polarized force field is employed in MD simulation, HIV-1 integrase CCD structure is stabilized and both the position of the Mg(2+) complex and the binding site are well preserved. The detailed analysis shows the transition of alpha-helix to 3(10)-helix adjacent to the catalytic loop (residues 139-147), which correlates with the dislocation of the Mg(2+) complex. The current study demonstrates the importance of electronic polarization of protein in stabilizing the metal complex in the catalytic core domain of HIV-1 integrase.
Collapse
Affiliation(s)
- Yunpeng Lu
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371, Singapore
| | | | | | | |
Collapse
|
26
|
|
27
|
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
|
28
|
|
29
|
Lameira J, Alves CN, Moliner V, Martí S, Kanaan N, Tuñón I. A quantum mechanics/molecular mechanics study of the protein-ligand interaction of two potent inhibitors of human O-GlcNAcase: PUGNAc and NAG-thiazoline. J Phys Chem B 2008; 112:14260-6. [PMID: 18939790 DOI: 10.1021/jp804626j] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
O-glycoprotein 2-acetamino-2-deoxy-beta- d-glucopyranosidase ( O-GlcNAcase) hydrolyzes 2-acetamido-2-deoxy-beta- d-glucopyranose ( O-GlcNAc) residues of serine/threonine residues of modified proteins. O-GlcNAc is present in many intracellular proteins and appears to have a role in the etiology of several diseases including cancer, Alzheimer's disease, and type II diabetes. In this work, we have carried out molecular dynamics simulations using a hybrid quantum mechanics/molecular mechanics approach to determine the binding of two potent inhibitors, PUGNAc and NAG, with a bacterial O-GlcNAcase. The results of these simulations show that Asp-401, Asp-298, and Asp-297 residues play an important role in the protein-inhibitor interactions. These results might be useful to design compounds with more interesting inhibitory activity on the basis of its three-dimensional structure.
Collapse
Affiliation(s)
- Jeronimo Lameira
- Laboratorio de Planejamento e Desenvolvimento de Farmacos, Instituto de Ciencias Exatas e Naturais, Universidade Federal do Para, CP 11101, 66075-110, Belem, PA, Brazil
| | | | | | | | | | | |
Collapse
|
30
|
Lameira J, Alves C, Santos L, Santos A, de Almeida Santos R, Souza J, Silva C, da Silva A. A combined X-ray and theoretical study of flavonoid compounds with anti-inflammatory activity. ACTA ACUST UNITED AC 2008. [DOI: 10.1016/j.theochem.2008.04.029] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
31
|
Lodola A, Woods CJ, Mulholland AJ. Applications and Advances of QM/MM Methods in Computational Enzymology. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2008. [DOI: 10.1016/s1574-1400(08)00009-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
|
32
|
A quantum mechanic/molecular mechanic study of the wild-type and N155S mutant HIV-1 integrase complexed with diketo acid. Biophys J 2007; 94:2443-51. [PMID: 17981909 DOI: 10.1529/biophysj.107.107623] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. Recently, mutation studies have been reported that have shown that a certain degree of viral resistance to diketo acids (DKAs) appears when some amino acid residues of the IN active site are mutated. Mutations represent a fascinating experimental challenge, and we invite theoretical simulations for the disclosure of still unexplored features of enzyme reactions. The aim of this work is to understand the molecular mechanisms of HIV-1 IN drug resistance, which will be useful for designing anti-HIV inhibitors with unique resistance profiles. In this study, we use molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy for wild-type and N155S mutant HIV-1 IN, both complexed with a DKA. This hybrid methodology has the advantage of the inclusion of quantum effects such as ligand polarization upon binding, which can be very important when highly polarizable groups are embedded in anisotropic environments, for example in metal-containing active sites. Furthermore, an energy terms decomposition analysis was performed to determine contributions of individual residues to the enzyme-inhibitor interactions. The results reveal that there is a strong interaction between the Lys-159, Lys-156, and Asn-155 residues and Mg(2+) cation and the DKA inhibitor. Our calculations show that the binding energy is higher in wild-type than in the N155S mutant, in accordance with the experimental results. The role of the mutated residue has thus been checked as maintaining the structure of the ternary complex formed by the protein, the Mg(2+) cation, and the inhibitor. These results might be useful to design compounds with more interesting anti-HIV-1 IN activity on the basis of its three-dimensional structure.
Collapse
|