1
|
Jokinen EM, Niemeläinen M, Kurkinen ST, Lehtonen JV, Lätti S, Postila PA, Pentikäinen OT, Niinivehmas SP. Virtual Screening Strategy to Identify Retinoic Acid-Related Orphan Receptor γt Modulators. Molecules 2023; 28:molecules28083420. [PMID: 37110655 PMCID: PMC10145393 DOI: 10.3390/molecules28083420] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 04/06/2023] [Accepted: 04/12/2023] [Indexed: 04/29/2023] Open
Abstract
Molecular docking is a key method used in virtual screening (VS) campaigns to identify small-molecule ligands for drug discovery targets. While docking provides a tangible way to understand and predict the protein-ligand complex formation, the docking algorithms are often unable to separate active ligands from inactive molecules in practical VS usage. Here, a novel docking and shape-focused pharmacophore VS protocol is demonstrated for facilitating effective hit discovery using retinoic acid receptor-related orphan receptor gamma t (RORγt) as a case study. RORγt is a prospective target for treating inflammatory diseases such as psoriasis and multiple sclerosis. First, a commercial molecular database was flexibly docked. Second, the alternative docking poses were rescored against the shape/electrostatic potential of negative image-based (NIB) models that mirror the target's binding cavity. The compositions of the NIB models were optimized via iterative trimming and benchmarking using a greedy search-driven algorithm or brute force NIB optimization. Third, a pharmacophore point-based filtering was performed to focus the hit identification on the known RORγt activity hotspots. Fourth, free energy binding affinity evaluation was performed on the remaining molecules. Finally, twenty-eight compounds were selected for in vitro testing and eight compounds were determined to be low μM range RORγt inhibitors, thereby showing that the introduced VS protocol generated an effective hit rate of ~29%.
Collapse
Affiliation(s)
- Elmeri M Jokinen
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship Center, University of Turku, FI-20014 Turku, Finland
| | - Miika Niemeläinen
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
| | - Sami T Kurkinen
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship Center, University of Turku, FI-20014 Turku, Finland
| | - Jukka V Lehtonen
- Structural Bioinformatics Laboratory, Biochemistry, Faculty of Science and Engineering, Åbo Akademi University, FI-20500 Turku, Finland
- InFLAMES Research Flagship Center, Åbo Akademi University, FI-20500 Turku, Finland
| | - Sakari Lätti
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship Center, University of Turku, FI-20014 Turku, Finland
| | - Pekka A Postila
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship Center, University of Turku, FI-20014 Turku, Finland
| | - Olli T Pentikäinen
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship Center, University of Turku, FI-20014 Turku, Finland
| | - Sanna P Niinivehmas
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship Center, University of Turku, FI-20014 Turku, Finland
| |
Collapse
|
2
|
Islam MA, Barshetty MM, Srinivasan S, Dudekula DB, Rallabandi VPS, Mohammed S, Natarajan S, Park J. Identification of Novel Ribonucleotide Reductase Inhibitors for Therapeutic Application in Bile Tract Cancer: An Advanced Pharmacoinformatics Study. Biomolecules 2022; 12:biom12091279. [PMID: 36139117 PMCID: PMC9496582 DOI: 10.3390/biom12091279] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 09/05/2022] [Accepted: 09/08/2022] [Indexed: 11/28/2022] Open
Abstract
Biliary tract cancer (BTC) is constituted by a heterogeneous group of malignant tumors that may develop in the biliary tract, and it is the second most common liver cancer. Human ribonucleotide reductase M1 (hRRM1) has already been proven to be a potential BTC target. In the current study, a de novo design approach was used to generate novel and effective chemical therapeutics for BTC. A set of comprehensive pharmacoinformatics approaches was implemented and, finally, seventeen potential molecules were found to be effective for the modulation of hRRM1 activity. Molecular docking, negative image-based ShaEP scoring, absolute binding free energy, in silico pharmacokinetics, and toxicity assessments corroborated the potentiality of the selected molecules. Almost all molecules showed higher affinity in comparison to gemcitabine and naphthyl salicylic acyl hydrazone (NSAH). On binding interaction analysis, a number of critical amino acids was found to hold the molecules at the active site cavity. The molecular dynamics (MD) simulation study also indicated the stability between protein and ligands. High negative MM-GBSA (molecular mechanics generalized Born and surface area) binding free energy indicated the potentiality of the molecules. Therefore, the proposed molecules might have the potential to be effective therapeutics for the management of BTC.
Collapse
Affiliation(s)
- Md Ataul Islam
- 3BIGS Omicscore Private Limited, 909 Lavelle Building, Richmond Circle, Bangalore 560025, India
| | | | - Sridhar Srinivasan
- 3BIGS Omicscore Private Limited, 909 Lavelle Building, Richmond Circle, Bangalore 560025, India
| | - Dawood Babu Dudekula
- 3BIGS Omicscore Private Limited, 909 Lavelle Building, Richmond Circle, Bangalore 560025, India
| | | | - Sameer Mohammed
- 3BIGS Omicscore Private Limited, 909 Lavelle Building, Richmond Circle, Bangalore 560025, India
| | | | - Junhyung Park
- 3BIGS Co., Ltd., B-831, Geumgang Penterium IX Tower, Hwaseong 18469, Korea
- Correspondence:
| |
Collapse
|
3
|
Identification of Potential Cytochrome P450 3A5 Inhibitors: An Extensive Virtual Screening through Molecular Docking, Negative Image-Based Screening, Machine Learning and Molecular Dynamics Simulation Studies. Int J Mol Sci 2022; 23:ijms23169374. [PMID: 36012627 PMCID: PMC9409045 DOI: 10.3390/ijms23169374] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 08/11/2022] [Accepted: 08/15/2022] [Indexed: 11/17/2022] Open
Abstract
Cytochrome P450 3A5 (CYP3A5) is one of the crucial CYP family members and has already proven to be an important drug target for cardiovascular diseases. In the current study, the PubChem database was screened through molecular docking and high-affinity molecules were adopted for further assessment. A negative image-based (NIB) model was used for a similarity search by considering the complementary shape and electrostatics of the target and small molecules. Further, the molecules were segregated into active and inactive groups through six machine learning (ML) matrices. The active molecules found in each ML model were used for in silico pharmacokinetics and toxicity assessments. A total of five molecules followed the acceptable pharmacokinetics and toxicity profiles. Several potential binding interactions between the proposed molecules and CYP3A5 were observed. The dynamic behavior of the selected molecules in the CYP3A5 was explored through a molecular dynamics (MD) simulation study. Several parameters obtained from the MD simulation trajectory explained the stability of the protein–ligand complexes in dynamic states. The high binding affinity of each molecule was revealed by the binding free energy calculation through the MM-GBSA methods. Therefore, it can be concluded that the proposed molecules might be potential CYP3A5 molecules for therapeutic application in cardiovascular diseases subjected to in vitro/in vivo validations.
Collapse
|
4
|
Ligand-Enhanced Negative Images Optimized for Docking Rescoring. Int J Mol Sci 2022; 23:ijms23147871. [PMID: 35887220 PMCID: PMC9323918 DOI: 10.3390/ijms23147871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Revised: 07/14/2022] [Accepted: 07/15/2022] [Indexed: 12/04/2022] Open
Abstract
Despite the pivotal role of molecular docking in modern drug discovery, the default docking scoring functions often fail to recognize active ligands in virtual screening campaigns. Negative image-based rescoring improves docking enrichment by comparing the shape/electrostatic potential (ESP) of the flexible docking poses against the target protein’s inverted cavity volume. By optimizing these negative image-based (NIB) models using a greedy search, the docking rescoring yield can be improved massively and consistently. Here, a fundamental modification is implemented to this shape-focused pharmacophore modelling approach—actual ligand 3D coordinates are incorporated into the NIB models for the optimization. This hybrid approach, labelled as ligand-enhanced brute-force negative image-based optimization (LBR-NiB), takes the best from both worlds, i.e., the all-roundedness of the NIB models and the difficult to emulate atomic arrangements of actual protein-bound small-molecule ligands. Thorough benchmarking, focused on proinflammatory targets, shows that the LBR-NiB routinely improves the docking enrichment over prior iterations of the R-NiB methodology. This boost can be massive, if the added ligand information provides truly essential binding information that was lacking or completely missing from the cavity-based NIB model. On a practical level, the results indicate that the LBR-NiB typically works well when the added ligand 3D data originates from a high-quality source, such as X-ray crystallography, and, yet, the NIB model compositions can also sometimes be improved by fusing into them, for example, with flexibly docked solvent molecules. In short, the study demonstrates that the protein-bound ligands can be used to improve the shape/ESP features of the negative images for effective docking rescoring use in virtual screening.
Collapse
|
5
|
Bhowmick S, AlFaris NA, Zaidan ALTamimi J, ALOthman ZA, Patil PC, Aldayel TS, Wabaidur SM, Saha A. Identification of bio-active food compounds as potential SARS-CoV-2 PLpro inhibitors-modulators via negative image-based screening and computational simulations. Comput Biol Med 2022; 145:105474. [PMID: 35395517 PMCID: PMC8973019 DOI: 10.1016/j.compbiomed.2022.105474] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Revised: 03/10/2022] [Accepted: 03/29/2022] [Indexed: 12/16/2022]
Abstract
Despite significant studies on the COVID-19 pandemic, scientists around the world are still battling to find a definitive therapy against the ongoing severe global health crisis. In this study, advanced computational approaches have been employed to identify bioactive food constituents as potential SARS-CoV-2 PLpro inhibitors-modulators. As a validated antiviral drug target, PLpro has gained tremendous attention for therapeutics developments. Therefore, targeting the multifunctional SARS-CoV-2 PLpro protein, ∼1039 bioactive dietary compounds have been screened extensively through novel techniques like negative image-based (NIB) screening and molecular docking approaches. In particular, the three different models of NIB screening have been generated and used to re-score the dietary compounds based on the negative image which is created by reversing the shape and electrostatics features of PLpro protein's ligand-binding cavity. Further, 100 ns molecular dynamics simulation has been performed and MM-GBSA based binding free energies have been estimated for the final proposed four dietary compounds (PC000550, PC000361, PC000558, and PC000573) as potential inhibitors/modulators of SARS-CoV-2 PLpro protein. Employed computational study outcome also has been compared with respect to the earlier experimentally investigated compound GRL0617 against SARS-CoV-2 PLpro protein, which suggests much greater interaction potential in terms of binding affinity and other energetic contributions for the proposed dietary compounds. Hence, the present study suggests that proposed dietary compounds can be suitable chemical entities for modulating the activity of PLpro protein or can be further utilized for optimizing or screening of novel chemical surrogates, however also needs experimental evaluation for entry in clinical studies for better assessment.
Collapse
Affiliation(s)
- Shovonlal Bhowmick
- Department of Chemical Technology, University of Calcutta, 92 A.P.C. Road, Kolkata, India
| | - Nora Abdullah AlFaris
- Nutrition and Food Science, Department of Physical Sport Science, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia,Corresponding author
| | - Jozaa Zaidan ALTamimi
- Nutrition and Food Science, Department of Physical Sport Science, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia
| | - Zeid A. ALOthman
- Chemistry Department, College of Science, King Saud University, Riyadh, Saudi Arabia
| | - Pritee Chunarkar Patil
- Department of Bioinformatics, Rajiv Gandhi Institute of IT and Biotechnology, Bharati Vidyapeeth Deemed University, Pune-Satara Road, Pune, India
| | - Tahany Saleh Aldayel
- Nutrition and Food Science, Department of Physical Sport Science, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia
| | | | - Achintya Saha
- Department of Chemical Technology, University of Calcutta, 92 A.P.C. Road, Kolkata, India,Corresponding author
| |
Collapse
|
6
|
Jokinen EM, Gopinath K, Kurkinen ST, Pentikäinen OT. Detection of Binding Sites on SARS-CoV-2 Spike Protein Receptor-Binding Domain by Molecular Dynamics Simulations in Mixed Solvents. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2021; 18:1281-1289. [PMID: 33914685 PMCID: PMC8791430 DOI: 10.1109/tcbb.2021.3076259] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Revised: 04/13/2021] [Accepted: 04/20/2021] [Indexed: 06/12/2023]
Abstract
The novel SARS-CoV-2 uses ACE2 (Angiotensin-Converting Enzyme 2) receptor as an entry point. Insights on S protein receptor-binding domain (RBD) interaction with ACE2 receptor and drug repurposing has accelerated drug discovery for the novel SARS-CoV-2 infection. Finding small molecule binding sites in S protein and ACE2 interface is crucial in search of effective drugs to prevent viral entry. In this study, we employed molecular dynamics simulations in mixed solvents together with virtual screening to identify small molecules that could be potential inhibitors of S protein -ACE2 interaction. Observation of organic probe molecule localization during the simulations revealed multiple sites at the S protein surface related to small molecule, antibody, and ACE2 binding. In addition, a novel conformation of the S protein was discovered that could be stabilized by small molecules to inhibit attachment to ACE2. The most promising binding site on RBD-ACE2 interface was targeted with virtual screening and top-ranked compounds (DB08248, DB02651, DB03714, and DB14826) are suggested for experimental testing. The protocol described here offers an extremely fast method for characterizing key proteins of a novel pathogen and for the identification of compounds that could inhibit or accelerate spreading of the disease.
Collapse
|
7
|
Negative Image-Based Screening: Rigid Docking Using Cavity Information. Methods Mol Biol 2021; 2266:125-140. [PMID: 33759124 DOI: 10.1007/978-1-0716-1209-5_7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Rational drug discovery relies heavily on molecular docking-based virtual screening, which samples flexibly the ligand binding poses against the target protein's structure. The upside of flexible docking is that the geometries of the generated docking poses are adjusted to match the residue alignment inside the target protein's ligand-binding pocket. The downside is that the flexible docking requires plenty of computing resources and, regardless, acquiring a decent level of enrichment typically demands further rescoring or post-processing. Negative image-based screening is a rigid docking technique that is ultrafast and computationally light but also effective as proven by vast benchmarking and screening experiments. In the NIB screening, the target protein cavity's shape/electrostatics is aligned and compared against ab initio-generated ligand 3D conformers. In this chapter, the NIB methodology is explained at the practical level and both its weaknesses and strengths are discussed candidly.
Collapse
|
8
|
Gopinath K, Jokinen EM, Kurkinen ST, Pentikäinen OT. Screening of Natural Products Targeting SARS-CoV-2-ACE2 Receptor Interface - A MixMD Based HTVS Pipeline. Front Chem 2020; 8:589769. [PMID: 33330376 PMCID: PMC7717977 DOI: 10.3389/fchem.2020.589769] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2020] [Accepted: 10/20/2020] [Indexed: 12/11/2022] Open
Abstract
The COVID-19 pandemic, caused by novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a severe global health crisis now. SARS-CoV-2 utilizes its Spike protein receptor-binding domain (S-protein) to invade human cell through binding to Angiotensin-Converting Enzyme 2 receptor (ACE2). S-protein is the key target for many therapeutics and vaccines. Potential S-protein-ACE2 fusion inhibitor is expected to block the virus entry into the host cell. In many countries, traditional practices, based on natural products (NPs) have been in use to slow down COVID-19 infection. In this study, a protocol was applied that combines mixed solvent molecular dynamics simulations (MixMD) with high-throughput virtual screening (HTVS) to search NPs to block SARS-CoV-2 entry into the human cell. MixMD simulations were employed to discover the most promising stable binding conformations of drug-like probes in the S-protein-ACE2 interface. Detected stable sites were used for HTVs of 612093 NPs to identify molecules that could interfere with the S-protein-ACE2 interaction. In total, 19 NPs were selected with rescoring model. These top-ranked NP-S-protein complexes were subjected to classical MD simulations for 300 ns (3 replicates of 100 ns) to estimate the stability and affinity of binding. Three compounds, ZINC000002128789, ZINC000002159944 and SN00059335, showed better stability in all MD runs, of which ZINC000002128789 was predicted to have the highest binding affinity, suggesting that it could be effective modulator in RBD-ACE2 interface to prevent SARS-CoV-2 infection. Our results support that NPs may provide tools to fight COVID-19.
Collapse
Affiliation(s)
| | | | | | - Olli T. Pentikäinen
- Faculty of Medicine, Integrative Physiology and Pharmacology, Institute of Biomedicine, University of Turku, Turku, Finland
| |
Collapse
|
9
|
A Perspective: Active Role of Lipids in Neurotransmitter Dynamics. Mol Neurobiol 2019; 57:910-925. [PMID: 31595461 PMCID: PMC7031182 DOI: 10.1007/s12035-019-01775-7] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Accepted: 09/01/2019] [Indexed: 12/30/2022]
Abstract
Synaptic neurotransmission is generally considered as a function of membrane-embedded receptors and ion channels in response to the neurotransmitter (NT) release and binding. This perspective aims to widen the protein-centric view by including another vital component—the synaptic membrane—in the discussion. A vast set of atomistic molecular dynamics simulations and biophysical experiments indicate that NTs are divided into membrane-binding and membrane-nonbinding categories. The binary choice takes place at the water-membrane interface and follows closely the positioning of the receptors’ binding sites in relation to the membrane. Accordingly, when a lipophilic NT is on route to a membrane-buried binding site, it adheres on the membrane and, then, travels along its plane towards the receptor. In contrast, lipophobic NTs, which are destined to bind into receptors with extracellular binding sites, prefer the water phase. This membrane-based sorting splits the neurotransmission into membrane-independent and membrane-dependent mechanisms and should make the NT binding into the receptors more efficient than random diffusion would allow. The potential implications and notable exceptions to the mechanisms are discussed here. Importantly, maintaining specific membrane lipid compositions (MLCs) at the synapses, especially regarding anionic lipids, affect the level of NT-membrane association. These effects provide a plausible link between the MLC imbalances and neurological diseases such as depression or Parkinson’s disease. Moreover, the membrane plays a vital role in other phases of the NT life cycle, including storage and release from the synaptic vesicles, transport from the synaptic cleft, as well as their synthesis and degradation.
Collapse
|