1
|
Prasad R, Kadam A, Padippurackal VV, Pulikuttymadom Balasubramanian A, Kumar Chandrakumaran N, Suresh Rangari K, Dnyaneshwar Khangar P, Ajith H, Natarajan K, Chandramohanadas R, Nelson-Sathi S. Discovery of small molecule entry inhibitors targeting the linoleic acid binding pocket of SARS-CoV-2 spike protein. J Biomol Struct Dyn 2024:1-15. [PMID: 38520147 DOI: 10.1080/07391102.2024.2327537] [Citation(s) in RCA: 0] [Impact Index Per Article: 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] [Received: 10/05/2023] [Accepted: 03/02/2024] [Indexed: 03/25/2024]
Abstract
Spike glycoprotein has a significant role in the entry of SARS-CoV-2 to host cells, which makes it a potential drug target. Continued accumulation of non-synonymous mutations in the receptor binding domain of spike protein poses great challenges in identifying antiviral drugs targeting this protein. This study aims to identify potential entry inhibitors of SARS-CoV-2 using virtual screening and molecular dynamics (MD) simulations from three distinct chemical libraries including Pandemic Response Box, Drugbank and DrugCentral, comprising 6971 small molecules. The molecules were screened against a binding pocket identified in the receptor-binding domain (RBD) region of the spike protein which is known as the linoleic acid binding pocket, a highly conserved motif among several SARS-CoV-2 variants. Through virtual screening and binding free energy calculations, we identified four top-scoring compounds, MMV1579787 ([2-Oxo-2-[2-(3-phenoxyphenyl)ethylamino]ethyl]phosphonic acid), Tretinoin, MMV1633963 ((2E,4E)-5-[3-(3,5-dichlorophenoxy)phenyl]penta-2,4-dienoic acid) and Polydatin, which were previously reported to have antibacterial, antifungal or antiviral properties. These molecules showed stable binding on MD simulations over 100 ns and maintained stable interactions with TYR365, PHE338, PHE342, PHE377, TYR369, PHE374 and LEU368 of the spike protein RBD that are found to be conserved among SARS-CoV-2 variants. Our findings were further validated with free energy landscape, principal component analysis and dynamic cross-correlation analysis. Our in silico analysis of binding mode and MD simulation analyses suggest that the identified compounds may impede viral entrance by interacting with the linoleic acid binding site of the spike protein of SARS-CoV-2 regardless of its variants, and they thus demand for further in vitro and in vivo research.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Roshny Prasad
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, Kerala, India
| | - Anil Kadam
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, Kerala, India
| | | | | | | | - Kartik Suresh Rangari
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, Kerala, India
| | | | - Harikrishnan Ajith
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, Kerala, India
| | - Kathiresan Natarajan
- Trans-disciplinary Biology, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, Kerala, India
| | | | - Shijulal Nelson-Sathi
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, Kerala, India
| |
Collapse
|
2
|
Prasad R, Ajith H, Kumar Chandrakumaran N, Dnyaneshwar Khangar P, Mohan A, Nelson-Sathi S. In silico study identifies peptide inhibitors that negate the effect of non-synonymous mutations in major drug targets of SARS-CoV-2 variants. J Biomol Struct Dyn 2023; 41:9551-9561. [PMID: 36377464 DOI: 10.1080/07391102.2022.2143426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Accepted: 10/29/2022] [Indexed: 11/16/2022]
Abstract
Since its advent in December 2019, SARS-CoV-2 has diverged into multiple variants with differing levels of virulence owing to the accumulation of mutations in its genome. The structural changes induced by non-synonymous mutations in major drug targets of the virus are known to alter the binding of potential antagonistic inhibitors. Here, we analyzed the effects of non-synonymous mutations in major targets of SARS-CoV-2 in response to potential peptide inhibitors. We screened 12 peptides reported to have anti-viral properties against RBD and 5 peptides against Mpro of SARS-CoV-2 variants using molecular docking and simulation approaches. The mutational landscape of RBD among SARS-CoV-2 variants had 21 non-synonymous mutations across 18 distinct sites. Among these, 14 mutations were present in the RBM region directly interacting with the hACE2 receptor. However, Only 3 non-synonymous mutations were observed in Mpro. We found that LCB1 - a de novo-synthesized peptide has the highest binding affinity to RBD despite non-synonymous mutations in variants and engages key residues of RBD-hACE2 interaction such as K417, E484, N487, and N501. Similarly, an antimicrobial peptide; 2JOS, was identified against Mpro with high binding affinity as it interacts with key residues in dimerization sites such as E166 and F140 crucial for viral replication. MD simulations affirm the stability of RBD-LCB1 and Mpro-2JOS complexes with an average RMSD of 1.902 and 2.476 respectively. We ascertain that LCB1 and 2JOS peptides are promising inhibitors to combat emerging variants of SARS-CoV-2 and thus warrant further investigations using in-vitro and in-vivo analysis.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Roshny Prasad
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, India
| | - Harikrishnan Ajith
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, India
| | | | | | - Anand Mohan
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, India
| | - Shijulal Nelson-Sathi
- Bioinformatics Laboratory, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, India
| |
Collapse
|
3
|
Sahoo RK, Manu S, Chandrakumaran NK, Vasudevan K. Nuclear and Mitochondrial Genome Assemblies of the Beetle, Zygogramma bicolorata, a Globally Important Biocontrol Agent of Invasive Weed Parthenium hysterophorus. Genome Biol Evol 2023; 15:evad188. [PMID: 37831427 PMCID: PMC10603765 DOI: 10.1093/gbe/evad188] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 10/05/2023] [Accepted: 10/09/2023] [Indexed: 10/14/2023] Open
Abstract
Implementing a genetic-based approach to achieve the full potential of classical biocontrol programs has been advocated for decades. The availability of genome-level information brings the opportunity to scrutinize biocontrol traits for their efficacy and evolvability. However, implementation of this advocacy remains limited to few instances. Biocontrol of a globally noxious weed, Parthenium hysterophorus, by the leaf-feeding beetle, Zygogramma bicolorata, has been in place for more than four decades now, with varying levels of success. As the first step in providing genetic-based improvement to this biocontrol program, we describe the nuclear and mitochondrial assemblies of Z. bicolorata. We assembled the genome from the long-read sequence data, error corrected with high-throughput short reads and checked for contaminants and sequence duplication to produce a 936 Mb nuclear genome. With 96.5% Benchmarking Universal Single-Copy Orthologs completeness and the long terminal repeat assembly index 12.91, we present a reference-quality assembly that appeared to be repeat rich at 62.7% genome-wide and consists of 29,437 protein-coding regions. We detected signature of nuclear insertion of mitochondrial fragments in 80 nuclear positions comprising 13 kb out of 17.9 kb mitochondria genome sequence. This genome, along with its annotations, provides a valuable resource to gain further insights into the biocontrol traits of Z. bicolorata for improving the control of the invasive weed P. hysterophorus.
Collapse
Affiliation(s)
- Ranjit Kumar Sahoo
- Laboratory for the Conservation of Endangered Species (LaCONES), CSIR-Centre for Cellular and Molecular Biology (CCMB), Hyderabad, India
| | - Shivakumara Manu
- Laboratory for the Conservation of Endangered Species (LaCONES), CSIR-Centre for Cellular and Molecular Biology (CCMB), Hyderabad, India
| | - Naveen Kumar Chandrakumaran
- Laboratory for the Conservation of Endangered Species (LaCONES), CSIR-Centre for Cellular and Molecular Biology (CCMB), Hyderabad, India
| | - Karthikeyan Vasudevan
- Laboratory for the Conservation of Endangered Species (LaCONES), CSIR-Centre for Cellular and Molecular Biology (CCMB), Hyderabad, India
| |
Collapse
|