Structure-based profiling of putative therapeutics against monkeypox virus VP39 using pharmacophore modelling and molecular dynamics simulation studies
{"title":"Structure-based profiling of putative therapeutics against monkeypox virus VP39 using pharmacophore modelling and molecular dynamics simulation studies","authors":"Gideon Ampoma Gyebi, Saheed Sabiu","doi":"10.1016/j.compbiolchem.2025.108458","DOIUrl":null,"url":null,"abstract":"<div><div>The growing global health threat of the monkeypox virus (MPXV) underscores the critical need for effective antiviral agents, since there are currently no therapeutics. The MPXV VP39, a methyltransferase, is essential for viral replication, hence a potential target for anti-MPXV drug candidates. Herein, a structure-based pharmacophore modelling and molecular docking approach was employed to screen natural compounds (NCs: 581,426) from the COCONUT database for potential inhibitors of MPXV VP39. After ranking of the docking scores, an ensemble-based docking of the top-ranked 20 NCs against multiple conformations obtained from ttcluster analysis of the molecular dynamics simulation trajectory of unbound MPXV VP39 further identified five leads with favourable interaction profiles, drug-likeness, ADMET properties, and synthetic features when compared to the reference standard (sinefungin). Further analysis of the thermodynamic stability of the resulting complexes of the leads over a 100-ns MD simulation period revealed varying degrees of thermodynamic stability while maintaining the structural integrity of MPXV VP39. Furthermore, the thermodynamic binding free energy calculation, while corroborating the docking analysis, identified CNP0297833 (-39.07 kcal/mol), CNP0371756 (-25.76 kcal/mol), and CNP0402319 (-19.26 kcal/mol) as the most promising candidates, with better modulatory effect against MPXV VP39 relative to sinefungin (-3.68 kcal/mol). These leads were stabilised with hydrophobic (Phe115, Val139, and Val116) and electrostatic (Glu46 and Asp138) interactions in different conformational clusters. In addition to the observed consistent interaction patterns, favourable binding energies, pharmacokinetics, ADMET, thermodynamic stability, and molecular orbital energies of these leads, the potential for optimisation for enhanced binding features for the active site of MPXV VP39 was elucidated. Further <em>in vitro</em> investigation to validate these findings is suggested to establish the putative leads as therapeutics targeting the replication phase of MPXV.</div></div>","PeriodicalId":10616,"journal":{"name":"Computational Biology and Chemistry","volume":"118 ","pages":"Article 108458"},"PeriodicalIF":2.6000,"publicationDate":"2025-04-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computational Biology and Chemistry","FirstCategoryId":"99","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S1476927125001185","RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"BIOLOGY","Score":null,"Total":0}
引用次数: 0
Abstract
The growing global health threat of the monkeypox virus (MPXV) underscores the critical need for effective antiviral agents, since there are currently no therapeutics. The MPXV VP39, a methyltransferase, is essential for viral replication, hence a potential target for anti-MPXV drug candidates. Herein, a structure-based pharmacophore modelling and molecular docking approach was employed to screen natural compounds (NCs: 581,426) from the COCONUT database for potential inhibitors of MPXV VP39. After ranking of the docking scores, an ensemble-based docking of the top-ranked 20 NCs against multiple conformations obtained from ttcluster analysis of the molecular dynamics simulation trajectory of unbound MPXV VP39 further identified five leads with favourable interaction profiles, drug-likeness, ADMET properties, and synthetic features when compared to the reference standard (sinefungin). Further analysis of the thermodynamic stability of the resulting complexes of the leads over a 100-ns MD simulation period revealed varying degrees of thermodynamic stability while maintaining the structural integrity of MPXV VP39. Furthermore, the thermodynamic binding free energy calculation, while corroborating the docking analysis, identified CNP0297833 (-39.07 kcal/mol), CNP0371756 (-25.76 kcal/mol), and CNP0402319 (-19.26 kcal/mol) as the most promising candidates, with better modulatory effect against MPXV VP39 relative to sinefungin (-3.68 kcal/mol). These leads were stabilised with hydrophobic (Phe115, Val139, and Val116) and electrostatic (Glu46 and Asp138) interactions in different conformational clusters. In addition to the observed consistent interaction patterns, favourable binding energies, pharmacokinetics, ADMET, thermodynamic stability, and molecular orbital energies of these leads, the potential for optimisation for enhanced binding features for the active site of MPXV VP39 was elucidated. Further in vitro investigation to validate these findings is suggested to establish the putative leads as therapeutics targeting the replication phase of MPXV.
期刊介绍:
Computational Biology and Chemistry publishes original research papers and review articles in all areas of computational life sciences. High quality research contributions with a major computational component in the areas of nucleic acid and protein sequence research, molecular evolution, molecular genetics (functional genomics and proteomics), theory and practice of either biology-specific or chemical-biology-specific modeling, and structural biology of nucleic acids and proteins are particularly welcome. Exceptionally high quality research work in bioinformatics, systems biology, ecology, computational pharmacology, metabolism, biomedical engineering, epidemiology, and statistical genetics will also be considered.
Given their inherent uncertainty, protein modeling and molecular docking studies should be thoroughly validated. In the absence of experimental results for validation, the use of molecular dynamics simulations along with detailed free energy calculations, for example, should be used as complementary techniques to support the major conclusions. Submissions of premature modeling exercises without additional biological insights will not be considered.
Review articles will generally be commissioned by the editors and should not be submitted to the journal without explicit invitation. However prospective authors are welcome to send a brief (one to three pages) synopsis, which will be evaluated by the editors.