1
|
Chen C, Zhu T, Liu X, Zhu D, Zhang Y, Wu S, Han C, Zhang H, Luo J, Kong L. Identification of a novel PHGDH covalent inhibitor by chemical proteomics and phenotypic profiling. Acta Pharm Sin B 2022; 12:246-261. [PMID: 35127383 PMCID: PMC8799887 DOI: 10.1016/j.apsb.2021.06.008] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Revised: 04/03/2021] [Accepted: 05/21/2021] [Indexed: 12/26/2022] Open
Abstract
The first rate-limiting enzyme of the serine synthesis pathway (SSP), phosphoglycerate dehydrogenase (PHGDH), is hyperactive in multiple tumors, which leads to the activation of SSP and promotes tumorigenesis. However, only a few inhibitors of PHGDH have been discovered to date, especially the covalent inhibitors of PHGDH. Here, we identified withangulatin A (WA), a natural small molecule, as a novel covalent inhibitor of PHGDH. Affinity-based protein profiling identified that WA could directly bind to PHGDH and inactivate the enzyme activity of PHGDH. Biolayer interferometry and LC-MS/MS analysis further demonstrated the selective covalent binding of WA to the cysteine 295 residue (Cys295) of PHGDH. With the covalent modification of Cys295, WA blocked the substrate-binding domain (SBD) of PHGDH and exerted an allosteric effect to induce PHGDH inactivation. Further studies revealed that with the inhibition of PHGDH mediated by WA, the glutathione synthesis was decreased and intracellular levels of reactive oxygen species (ROS) were elevated, leading to the inhibition of tumor proliferation. This study indicates WA as a novel PHGDH covalent inhibitor, which identifies Cys295 as a novel allosteric regulatory site of PHGDH and holds great potential in developing anti-tumor agents for targeting PHGDH.
Collapse
Key Words
- 3-PG, 3-phosphoglycerate
- 3-PHP, 3-phosphohydroxypyruvate
- ABPP, affinity-based protein profiling
- BLI, biolayer interferometry assay
- CETSA, cellular thermal shift assay
- Chemical proteomics
- Colon cancer
- Covalent inhibitor
- CuAAC, copper-catalyzed alkyne–azide cycloaddition
- DARTS, drug affinity responsive target stability
- GSH, glutathione
- MD, molecular dynamics
- NADPH, nicotinamide adenine dinucleotide phosphate
- Oxidative stress
- PHGDH, phosphoglycerate dehydrogenase
- PSAT, phosphoserine aminotransferase
- Phosphoglycerate dehydrogenase
- RMSD, root mean square deviation
- RMSF, root mean square fluctuations
- ROS, reactive oxygen species
- SBD, substrate-binding domain
- SSP, serine synthesis pathway
- Serine synthesis pathway
- TBTA, tris[(1-benzyl-1H-1,2,3-triazol-4-yl)methyl]amine
- TCEP, tris(2-carboxyethyl) phosphine
- Withangulatin A
- Withanolides
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | - Jianguang Luo
- Corresponding authors. Tel./fax: +86 25 83271405, +86 25 83271402.
| | - Lingyi Kong
- Corresponding authors. Tel./fax: +86 25 83271405, +86 25 83271402.
| |
Collapse
|
2
|
Kumar N, Singh A, Gulati HK, Bhagat K, Kaur K, Kaur J, Dudhal S, Duggal A, Gulati P, Singh H, Singh JV, Bedi PMS. Phytoconstituents from ten natural herbs as potent inhibitors of main protease enzyme of SARS-COV-2: In silico study. Phytomed Plus 2021; 1:100083. [PMID: 35403086 PMCID: PMC8180089 DOI: 10.1016/j.phyplu.2021.100083] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2021] [Revised: 05/31/2021] [Accepted: 06/01/2021] [Indexed: 04/13/2023]
Abstract
BACKGROUND Lack of treatment of novel Coronavirus disease led to the search of specific antivirals that are capable to inhibit the replication of the virus. The plant kingdom has demonstrated to be an important source of new molecules with antiviral potential. PURPOSE The present study aims to utilize various computational tools to identify the most eligible drug candidate that have capabilities to halt the replication of SARS-COV-2 virus by inhibiting Main protease (Mpro) enzyme. METHODS We have selected plants whose extracts have inhibitory potential against previously discovered coronaviruses. Their phytoconstituents were surveyed and a library of 100 molecules was prepared. Then, computational tools such as molecular docking, ADMET and molecular dynamic simulations were utilized to screen the compounds and evaluate them against Mpro enzyme. RESULTS All the phytoconstituents showed good binding affinities towards Mpro enzyme. Among them laurolitsine possesses the highest binding affinity i.e. -294.1533 kcal/mol. On ADMET analysis of best three ligands were simulated for 1.2 ns, then the stable ligand among them was further simulated for 20 ns. Results revealed that no conformational changes were observed in the laurolitsine w.r.t. protein residues and low RMSD value suggested that the Laurolitsine-protein complex was stable for 20 ns. CONCLUSION Laurolitsine, an active constituent of roots of Lindera aggregata, was found to be having good ADMET profile and have capabilities to halt the activity of the enzyme. Therefore, this makes laurolitsine a good drug candidate for the treatment of COVID-19.
Collapse
Key Words
- ACE-2, Angiotensin converting enzyme- 2
- ADMET
- ADMET, absorption, Distribution, metabolism, excretion and toxicity
- Ala, Alanine
- Approx., approximately
- Arg, arginine
- Asn, Asparagine
- Asp, Aspartic acid
- CADD, Computer Aided Drug Design
- CHARMM, Chemistry at Harvard Macromolecular Mechanics
- COV, coronavirus
- COVID, Novel corona-virus disease
- Covid-19
- Cys, cysteine
- DSBDS, Dassault's Systems Biovia's Discovery studio
- Gln, Glutamine
- Glu, glutamate
- Gly, Glycine
- His, histidine
- Ile, isoleucine
- K, Kelvin
- Kcal/mol, kilo calories per mol
- Leu, Leucine
- Leu, leucine
- Lys, Lysine
- MD, Molecular Dynamics
- Met, Methionine
- MoISA, Molecular Surface Area
- Molecular dynamic simulations
- Mpro protein
- Mpro, Main protease enzyme
- N protein, nucleocapsid protein
- NI, N-(4-methylpyridin-3-yl) acetamide inhibitor
- NPT, amount of substance (N), pressure (P) and temperature (T)
- NVT, amount of substance (N), volume (V) and temperature (T)
- Natural Antiviral herbs
- PDB, protein data bank
- PPB, plasma protein binding
- PSA, Polar Surface Area
- Phi, Phenylalanine
- Pro, Proline
- RCSB, Research Collaboratory for Structural Bioinformatics
- RMS, Root Mean Square
- RMSD, Root Mean Square Deviation
- RMSF, root mean square fluctuations
- RNA, Ribonucleic acid
- SAR-COV-2, severe acute respiratory syndrome coronavirus 2
- SDF, structure data format
- Ser, serine
- T, Temperature
- Thr, Threonine
- Trp, Tryptophan
- Tyr, Tyrosine
- Val, Valine
- kDa, kilo Dalton
- nCOV-19, Novel Coronavirus 2019
- ns/nsec, nano seconds
- ps, pentoseconds
- rGyr, Radius of gyration
- w.r.t., with respect to
- Å, angstrom
- α, alpha
- β, beta
Collapse
Affiliation(s)
- Nitish Kumar
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
- Drug and Pollution testing Lab, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Atamjit Singh
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Harmandeep Kaur Gulati
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Kavita Bhagat
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Komalpreet Kaur
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Jaspreet Kaur
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Shilpa Dudhal
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Amit Duggal
- Drugs Control Wing, Sector 16, Chandigarh, India, 160015
| | - Puja Gulati
- School of Pharmacy, Desh Bhagat University, Mandi Gobindgarh, Punjab, India, 147301
| | - Harbinder Singh
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Jatinder Vir Singh
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | | |
Collapse
|
3
|
Kumar N, Singh A, Gulati HK, Bhagat K, Kaur K, Kaur J, Dudhal S, Duggal A, Gulati P, Singh H, Singh JV, Bedi PMS. Phytoconstituents from ten natural herbs as potent inhibitors of main protease enzyme of SARS-COV-2: In silico study. Phytomed Plus 2021. [PMID: 35403086 DOI: 10.1016/j.phyplu.2021.100139] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
BACKGROUND Lack of treatment of novel Coronavirus disease led to the search of specific antivirals that are capable to inhibit the replication of the virus. The plant kingdom has demonstrated to be an important source of new molecules with antiviral potential. PURPOSE The present study aims to utilize various computational tools to identify the most eligible drug candidate that have capabilities to halt the replication of SARS-COV-2 virus by inhibiting Main protease (Mpro) enzyme. METHODS We have selected plants whose extracts have inhibitory potential against previously discovered coronaviruses. Their phytoconstituents were surveyed and a library of 100 molecules was prepared. Then, computational tools such as molecular docking, ADMET and molecular dynamic simulations were utilized to screen the compounds and evaluate them against Mpro enzyme. RESULTS All the phytoconstituents showed good binding affinities towards Mpro enzyme. Among them laurolitsine possesses the highest binding affinity i.e. -294.1533 kcal/mol. On ADMET analysis of best three ligands were simulated for 1.2 ns, then the stable ligand among them was further simulated for 20 ns. Results revealed that no conformational changes were observed in the laurolitsine w.r.t. protein residues and low RMSD value suggested that the Laurolitsine-protein complex was stable for 20 ns. CONCLUSION Laurolitsine, an active constituent of roots of Lindera aggregata, was found to be having good ADMET profile and have capabilities to halt the activity of the enzyme. Therefore, this makes laurolitsine a good drug candidate for the treatment of COVID-19.
Collapse
Key Words
- ACE-2, Angiotensin converting enzyme- 2
- ADMET
- ADMET, absorption, Distribution, metabolism, excretion and toxicity
- Ala, Alanine
- Approx., approximately
- Arg, arginine
- Asn, Asparagine
- Asp, Aspartic acid
- CADD, Computer Aided Drug Design
- CHARMM, Chemistry at Harvard Macromolecular Mechanics
- COV, coronavirus
- COVID, Novel corona-virus disease
- Covid-19
- Cys, cysteine
- DSBDS, Dassault's Systems Biovia's Discovery studio
- Gln, Glutamine
- Glu, glutamate
- Gly, Glycine
- His, histidine
- Ile, isoleucine
- K, Kelvin
- Kcal/mol, kilo calories per mol
- Leu, Leucine
- Leu, leucine
- Lys, Lysine
- MD, Molecular Dynamics
- Met, Methionine
- MoISA, Molecular Surface Area
- Molecular dynamic simulations
- Mpro protein
- Mpro, Main protease enzyme
- N protein, nucleocapsid protein
- NI, N-(4-methylpyridin-3-yl) acetamide inhibitor
- NPT, amount of substance (N), pressure (P) and temperature (T)
- NVT, amount of substance (N), volume (V) and temperature (T)
- Natural Antiviral herbs
- PDB, protein data bank
- PPB, plasma protein binding
- PSA, Polar Surface Area
- Phi, Phenylalanine
- Pro, Proline
- RCSB, Research Collaboratory for Structural Bioinformatics
- RMS, Root Mean Square
- RMSD, Root Mean Square Deviation
- RMSF, root mean square fluctuations
- RNA, Ribonucleic acid
- SAR-COV-2, severe acute respiratory syndrome coronavirus 2
- SDF, structure data format
- Ser, serine
- T, Temperature
- Thr, Threonine
- Trp, Tryptophan
- Tyr, Tyrosine
- Val, Valine
- kDa, kilo Dalton
- nCOV-19, Novel Coronavirus 2019
- ns/nsec, nano seconds
- ps, pentoseconds
- rGyr, Radius of gyration
- w.r.t., with respect to
- Å, angstrom
- α, alpha
- β, beta
Collapse
Affiliation(s)
- Nitish Kumar
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
- Drug and Pollution testing Lab, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Atamjit Singh
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Harmandeep Kaur Gulati
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Kavita Bhagat
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Komalpreet Kaur
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Jaspreet Kaur
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Shilpa Dudhal
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Amit Duggal
- Drugs Control Wing, Sector 16, Chandigarh, India, 160015
| | - Puja Gulati
- School of Pharmacy, Desh Bhagat University, Mandi Gobindgarh, Punjab, India, 147301
| | - Harbinder Singh
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | - Jatinder Vir Singh
- Department of Pharmaceutical Sciences, Guru Nanak Dev University, Amritsar, Punjab, India, 143005
| | | |
Collapse
|
4
|
Ahmad Mir S, Firoz A, Alaidarous M, Alshehri B, Aziz Bin Dukhyil A, Banawas S, Alsagaby SA, Alturaiki W, Ahmad Bhat G, Kashoo F, Abdel-Hadi AM. Identification of SARS-CoV-2 RNA-dependent RNA polymerase inhibitors from the major phytochemicals of Nigella sativa: An in silico approach. Saudi J Biol Sci 2021; 29:394-401. [PMID: 34518755 PMCID: PMC8426002 DOI: 10.1016/j.sjbs.2021.09.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2021] [Revised: 08/29/2021] [Accepted: 09/01/2021] [Indexed: 01/18/2023] Open
Abstract
The coronavirus disease 2019 (COVID-19), which emerged in December 2019, continues to be a serious health concern worldwide. There is an urgent need to develop effective drugs and vaccines to control the spread of this disease. In the current study, the main phytochemical compounds of Nigella sativa were screened for their binding affinity for the active site of the RNA-dependent RNA polymerase (RdRp) enzyme of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The binding affinity was investigated using molecular docking methods, and the interaction of phytochemicals with the RdRp active site was analyzed and visualized using suitable software. Out of the nine phytochemicals of N. sativa screened in this study, a significant docking score was observed for four compounds, namely α-hederin, dithymoquinone, nigellicine, and nigellidine. Based on the findings of our study, we report that α-hederin, which was found to possess the lowest binding energy (–8.6 kcal/mol) and hence the best binding affinity, is the best inhibitor of RdRp of SARS-CoV-2, among all the compounds screened here. Our results prove that the top four potential phytochemical molecules of N. sativa, especially α-hederin, could be considered for ongoing drug development strategies against SARS-CoV-2. However, further in vitro and in vivo testing are required to confirm the findings of this study.
Collapse
Affiliation(s)
- Shabir Ahmad Mir
- Department of Medical Laboratory Sciences, College of Applied Medical Science, Majmaah University, Al Majmaah 11952, Saudi Arabia
| | - Ahmad Firoz
- Department of Biological Sciences, Faculty of Science, King Abdulaziz University, Jeddah, Ssaudi Arabia
| | - Mohammed Alaidarous
- Department of Medical Laboratory Sciences, College of Applied Medical Science, Majmaah University, Al Majmaah 11952, Saudi Arabia.,Health and Basic Sciences Research Center, Majmaah University, Al Majmaah 11952, Saudi Arabia
| | - Bader Alshehri
- Department of Medical Laboratory Sciences, College of Applied Medical Science, Majmaah University, Al Majmaah 11952, Saudi Arabia
| | - Abdul Aziz Bin Dukhyil
- Department of Medical Laboratory Sciences, College of Applied Medical Science, Majmaah University, Al Majmaah 11952, Saudi Arabia
| | - Saeed Banawas
- Department of Medical Laboratory Sciences, College of Applied Medical Science, Majmaah University, Al Majmaah 11952, Saudi Arabia.,Health and Basic Sciences Research Center, Majmaah University, Al Majmaah 11952, Saudi Arabia.,Department of Biomedical Sciences, Oregon State University, Corvallis, OR 97331, USA
| | - Suliman A Alsagaby
- Department of Medical Laboratory Sciences, College of Applied Medical Science, Majmaah University, Al Majmaah 11952, Saudi Arabia
| | - Wael Alturaiki
- Department of Medical Laboratory Sciences, College of Applied Medical Science, Majmaah University, Al Majmaah 11952, Saudi Arabia
| | - Gulzar Ahmad Bhat
- Department of Clinical Biochemistry, Sher-i-Kashmir Institute of Medical Science, Srinagar, India
| | - Faizan Kashoo
- Department of Physical Therapy and Health Rehabilitation, College of Applied Medical Sciences, Majmaah University, Al Majmaah-11952, Saudi Arabia
| | - Ahmad M Abdel-Hadi
- Department of Medical Laboratory Sciences, College of Applied Medical Science, Majmaah University, Al Majmaah 11952, Saudi Arabia
| |
Collapse
|
5
|
Cruz JN, Costa JFS, Khayat AS, Kuca K, Barros CAL, Neto AMJC. Molecular dynamics simulation and binding free energy studies of novel leads belonging to the benzofuran class inhibitors of Mycobacterium tuberculosis Polyketide Synthase 13. J Biomol Struct Dyn 2018; 37:1616-1627. [PMID: 29633908 DOI: 10.1080/07391102.2018.1462734] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Abstract
In this work, the binding mechanism of new Polyketide Synthase 13 (Pks13) inhibitors has been studied through molecular dynamics simulation and free energy calculations. The drug Tam1 and its analogs, belonging to the benzofuran class, were submitted to 100 ns simulations, and according to the results obtained for root mean square deviation, all the simulations converged from approximately 30 ns. For the analysis of backbone flotation, the root mean square fluctuations were plotted for the Cα atoms; analysis revealed that the greatest fluctuation occurred in the residues that are part of the protein lid domain. The binding free energy value (ΔGbind) obtained for the Tam16 lead molecule was of -51.43 kcal/mol. When comparing this result with the ΔGbind values for the remaining analogs, the drug Tam16 was found to be the highest ranked: this result is in agreement with the experimental results obtained by Aggarwal and collaborators, where it was verified that the IC50 for Tam16 is the smallest necessary to inhibit the Pks13 (IC50 = 0.19 μM). The energy decomposition analysis suggested that the residues which most interact with inhibitors are: Ser1636, Tyr1637, Asn1640, Ala1667, Phe1670, and Tyr1674, from which the greatest energy contribution to Phe1670 was particularly notable. For the lead molecule Tam16, a hydrogen bond with the hydroxyl of the phenol not observed in the other analogs induced a more stable molecular structure. Aggarwal and colleagues reported this hydrogen bonding as being responsible for the stability of the molecule, optimizing its physic-chemical, toxicological, and pharmacokinetic properties.
Collapse
Key Words
- CNPq, National Council for Scientific and Technological Development
- CoA, coenzyme A
- FAS, fatty acid synthase
- GAFF, general amber force field
- GB, generalized born
- HB, hydrogen bonds
- INH, isoniazid
- KatG, catalase-peroxidase
- MD, molecular dynamics
- MDR, multi-drug resistant
- MM/GBSA, molecular mechanics/generalized-born surface area
- NAD, nicotinamide adenine dinucleotide
- NS, nanoseconds
- PCA, acyl carrier protein
- Pks13
- Pks13, polyketide synthase 13
- RESP, restrained electrostatic potential
- RMSD, root mean square deviation
- RMSF, root mean square fluctuations
- TB, tuberculosis
- TE, C-terminal thioesterase
- XDR, extensively drug resistant
- benzofuran
- free energy
- inhibitors
- molecular dynamics
- Δ internal energy
- Δ, Van Der Waals contributions
- Δ, electrostatic contribution
- Δ, electrostatic contributions
- Δ, energy of desolvation
- Δ, energy of the molecular mechanics
- Δ, non-polar contributions
- Δ, polar contributions
- Δ, polar solvation contribution
Collapse
Affiliation(s)
- Jorddy N Cruz
- a Laboratory of Preparation and Computation of Nanomaterials , Federal University of Pará , CP 479, 66075-110 Belém , PA , Brazil
| | - José F S Costa
- a Laboratory of Preparation and Computation of Nanomaterials , Federal University of Pará , CP 479, 66075-110 Belém , PA , Brazil
| | - André S Khayat
- b Oncology Research Center , Federal University of Pará , CP 479, 6675-105 Belém , PA , Brazil
| | - Kamil Kuca
- c Biomedical Research Center , University Hospital Hradec Kralove , Sokolska 581, 500 05 Hradec Kralove , Czech Republic
| | - Carlos A L Barros
- d Department of Pharmacy , Federal University of Pará , CP 479, 66050-160 Belém , PA , Brazil
| | - A M J C Neto
- a Laboratory of Preparation and Computation of Nanomaterials , Federal University of Pará , CP 479, 66075-110 Belém , PA , Brazil
| |
Collapse
|
6
|
Kumar R, Jade D, Gupta D. A novel identification approach for discovery of 5-HydroxyTriptamine 2A antagonists: combination of 2D/3D similarity screening, molecular docking and molecular dynamics. J Biomol Struct Dyn 2018; 37:931-943. [PMID: 29468945 DOI: 10.1080/07391102.2018.1444509] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
5-HydroxyTriptamine 2A antagonists are potential targets for treatment of various cerebrovascular and cardiovascular disorders. In this study, we have developed and performed a unique screening pipeline for filtering ZINC database compounds on the basis of similarities to known antagonists to determine novel small molecule antagonists of 5-HydroxyTriptamine 2A. The screening pipeline is based on 2D similarity, 3D dissimilarity and a combination of 2D/3D similarity. The shortlisted compounds were docked to a 5-HydroxyTriptamine 2A homology-based model, and complexes with low binding energies (287 complexes) were selected for molecular dynamics (MD) simulations in a lipid bilayer. The MD simulations of the shortlisted compounds in complex with 5-HydroxyTriptamine 2A confirmed the stability of the complexes and revealed novel interaction insights. The receptor residues S239, N343, S242, S159, Y370 and D155 predominantly participate in hydrogen bonding. π-π stacking is observed in F339, F340, F234, W151 and W336, whereas hydrophobic interactions are observed amongst V156, F339, F234, V362, V366, F340, V235, I152 and W151. The known and potential antagonists shortlisted by us have similar overlapping molecular interaction patterns. The 287 potential 5-HydroxyTriptamine 2A antagonists may be experimentally verified.
Collapse
Key Words
- , tanimoto coefficient
- 2D similarity
- 2D, two-dimensional space
- 2D/3D screening
- 3D similarity
- 3D, three-dimensional space
- 5HT
- 5HT, 5-HydroxyTryptamine
- ADHD, attention deficit hyperactivity disorders
- BLAST, basic local alignment search tool
- CNS, central nervous system
- Cl ions, chloride ions
- DOPE, discrete optimized protein energy
- G-protein coupled receptor
- GPCRs, G protein-coupled receptors
- HB, hydrogen bond
- HBA, hydrogen bond acceptors
- HBD, hydrogen bond donors
- JC virus, John Cunningham virus
- Ki, equilibrium dissociation constant for the ligand
- LBVS, ligand-based virtual screening
- MD, molecular dynamic
- MSD, mean square displacement
- MW, molecular weight
- NHB, number of hydrogen bonds
- OCD, obsessive compulsive disorder
- P5/P95, percentile calculation
- PAINS, Pan assay interference compounds
- PDB, protein data bank
- PLIP, protein–ligand interaction profiler
- PME, Particle Mesh Ewald
- PNS, peripheral nervous system
- POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine
- RMSD, root mean square deviation
- RMSF, root mean square fluctuations
- Rg, radius of gyration
- SASA, solvent accessible surface area
- SCA, stochastic clustering algorithm
- SD, steepest descent
- SDF, structure data file
- SPC, single point charge
- SPD, simple point charge
- SSE, secondary structure elements
- Sn-1/sn-2, Stereospecific number
- TM, Transmembrane
- TPSA, topological polar surface area
- drug discovery
- fs, femtosecond
- kJ/mol, kilo Joule per mol
- kcal/mol, kilocalorie per mole sn-1
- ligand-based virtual screening
- nm, nanomolar
- ns, nanosecond
- Å Ångström
- β2-AR, β2 adrenergic receptor
- μM, micromolar
Collapse
Affiliation(s)
- Rakesh Kumar
- a Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology (ICGEB) , Aruna Asaf Ali Marg, New Delhi 110067 , India
| | - Dhananjay Jade
- a Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology (ICGEB) , Aruna Asaf Ali Marg, New Delhi 110067 , India
| | - Dinesh Gupta
- a Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology (ICGEB) , Aruna Asaf Ali Marg, New Delhi 110067 , India
| |
Collapse
|