key: cord-0749537-isbbsece
authors: Siddiqa, M A.; Rao, D. S.; Suvarna, G.; Chennamachetty, V. K.; Verma, M. K.; Rao, M. V. R.
title: In-Silico Drug Designing of Spike Receptor with Its ACE2 Receptor and Nsp10/Nsp16 MTase Complex Against SARS-CoV-2
date: 2021-03-17
journal: Int J Pept Res Ther
DOI: 10.1007/s10989-021-10196-x
sha: 5d946114bc80fa1347794a224827ced5cc8c0aaa
doc_id: 749537
cord_uid: isbbsece

The realm Riboviria constitutes Coronaviruses, which led to the emergence of the pandemic COVID 19 in the twenty-first century affected millions of lives. At present, the management of COVID 19 largely depends on antiviral therapeutics along with the anti-inflammatory drug. The vaccine is under the final clinical phase, and emergency use is available. We aim at ACE2 and Nsp10/Nsp16 MTase as potential drug candidate in COVID 19 management in the present work. For drug designing, various computational simulation strategies have been employed like Swiss-Model, Hawk Dock, HDOCK, py Dock, and PockDrug for homology modeling, binding energies of the molecule with a target, simulate the conformation and binding poses, statistics of protein lock with target key and drug ability, respectively. The current in-silico screening depicts that the spike protein receptor is complementary to the target when bound to each other and forms a stable complex. The MMGBSA free energy binding property of receptor and ligand is critical. The intermolecular Statistics with the target Nsp10/Nsp16 MTase complex are plausible. We have also observed a high-affinity pocket binding site with the target. Therefore, the favorable intermolecular interactions and Physico-chemical properties emanate as a drug candidate treating COVID-19. This study has approached computational tools to analyze the conformation, binding affinity, and drug ability of receptor-ligand. Thus, the spike receptor with its ACE2 receptor with Nsp10/Nsp16 MTase complex would be a potent drug against SARS CoV-2 and can cure the infection as per consensus scoring.

In December 2019 the first COVID 19 outbreak reported in Wuhan, China as mysterious and contagious viral infection leading to viral pneumonia (Bello 2020) . The causative agent for COVID 19 was identified as novel severe acute respiratory syndrome corona virus 2 (SARS-CoV2). The corona viruses are non-segmented RNA viruses generally infect animals precisely wild animal however a few strain are also reported in human as well (Chellapandi and Saranya 2020) . Later, a new virus and plague was explored then by WHO as virus pandemic because of the severity of spread at an alarming level in March 2020. COVID-19 is now a pandemic influencing many countries globally . Coronaviruses are a broad family of viruses which may lead to illness in humans. In humans, several Corona Viruses cause respiratory diseases covering from the common cold to more stringent diseases such as Middle East Respiratory Syndrome (MERS) and Severe Acute Respiratory Syndrome (SARS) (Xu et al. 2020 ). Both, vaccines and therapeutics are under trials worldwide to prevent and treatment of COVID 19 (Choudhary et al. 2020 ). There are growing evidences suggesting viral infections often modulate host defense via immune regulation (Verma and Shakya 2020) .

SARS-CoV-2 is extracellular with protein-coat surrounded by a core containing nucleic acid. It is 50-200 nm in diameter. SARS-CoV-2 includes four structural proteins, referred to as the S (spike), E (envelope), M (membrane), and N (nucleocapsid) proteins (Ahmad et al. 2020) . The genome of the Coronavirus encodes a range of structural proteins that facilitate cellular entry and assembly of virions, of that the spike protein S turns out to be crucial for cellular entry (Alazmi and Motwalli 2020) . The spike protein regulates the virus to adhere to the host cell. S protein has within a receptor viz., Receptor Binding Domain. The RBD of spike protein S binds to Angiotensin-converting protein two (ACE2) to trigger cellular entry (Albini et al. 2020) . The spike protein reveals the cleavage sites to cellular proteases when bound to the ACE2 receptor. Subsequently, after the cleavage of spike protein by transmembrane protease serine 2(TMPRSS2) and varied cellular proteases starts fusion and endocytosis (Baig et al. 2020) . The spike protein contains addendum furin cleavage site and that will endorse it to be reactivated to gain access to intrude in the cell and vigorously infect once replication starts (Davidson et al. 2020 ).

Structure of Nsp10.The SARS-CoV-2 Nsp10 structure comprises 19-133 residues. An electropositive hydrophobic surface interacts with an electronegative hydrophobic pocket from Nsp16, which contribute the stabilization of the SAM binding site (Donoghue et al. 2020) . Structure of Nsp16. Nsp16 structure has all 298 residues. The structure of Nsp16 involves the 2-O-MTase chemical action core composed of a Rossmann-like β-sheet fold formed by α-helices(11), b-strands(7), and loops. The chemical action core requires one molecule of SAM (near the β1 and β2 strands of the Rossmann-like fold). In eukaryotes, many viruses that replicate inside the substance that has developed 2′-O-methyltransferases (2′-O-MTase) to switch their mRNAs and carry a cap-1 structure (m7GpppNm) at the 5′ end. Similarly, Nsp10 is the stimulatory entity for Nsp16 (El Hassab et al. 2020) . This is helpful in SAM-binding pocket and enhancing the substrate RNA-binding groove of nsp16. Intervention by tiny peptides between nsp10 and nsp16 of SARS-CoV may inhibit the activity of 2′-O-MTase. Moreover, the 5′ cap structures of organism mRNAs are crucial for RNA stability and protein translation (Gorbalenya et al. 2020) . Thereby aiding infectious agent (viral) replication and escaping innate immune recognition in host cells. The study was aimed to explore 3D structure of drug and target along with desirable functional properties for affinity of drug towards receptor (Guo et al. 2020) . We performed molecular modeling, docking and simulation studies for drug receptor interaction and affinity.

The structure of novel coronavirus spike receptor-binding domain complexes with its receptor ACE 2 (protein1) obtained from Protein Data Bank (www. rcsb. org) with PDB ID: 6LZG (Hamming et al. 2004 ). The target (protein2) viz, room temperature structure of sars-cov-2 Nsp10/Nsp16 methyl-transferase in a complex with sam determined by fixed-target serial crystallography with PDB ID: 6XKM. The protein sequence is retrieved from Uniprot KB (www. unipr ot. org) bearing Uniprot Id: P0DTC2 (SPIKE_SARS2) (Hoffman et al. 2020). The protein sequence of target is retrieved from Uniprot Id: P0DTD1 (R1AB_SARS2) was collected.

A computational tool which relies on the principles from known protein structure/unknown protein sequence. The model is obtained via X-Ray diffraction (Hussein et al. 2015) . The template selection, alignment and model building as well are automated by the server. It generates a reliable three-dimensional protein structure. It compares target-template giving GMQE, QMEAN (https:// swiss model. expasy. org/).

A network server to predict and analyze the structures of protein-protein complexes basing on computational docking and MM/GBSA (Molecular mechanics generalized born surface area) (Kaswa and Govender 2004) . MM/GBSA is a high-throughput virtual screening method widely used to predict binding free energies and to identify the correct binding conformations in the protein-protein complex (http:// cadd. zju. edu. cn/ hawkd ock/).

It is based on a hybrid algorithm which is template-based modeling and ab-initio docking. The prognosis is topranked performed docking server when submitted the CAPRI experiments (Krishnan et al. 2020) . It is highly integrated suite supports intrinsic scoring function with Protein/DNA-RNA docking (http:// hdock. phys. hust. edu. cn/).

A classic approach to the scoring of rigid-body docking poses with FTT based algorithm, which has been implemented in a program pyDock (Lai et al. 2020 ). The program is based on Columbic electrostatics along with distance-dependent dielectric constant, desolvation and Van Der Waal energy are the molecular dynamics integrated along with protein-protein binding (Letko et al. 2020 ) https:// life. bsc. es/ pid/ pydock/.

PockDrug-Server authorizes to predict drug ability probability for a protein pocket in the interest of target recognition aspect of drug discovery (Li et al. 2005 . The drug ability of pockets is predicted on holo and apo-proteins. The pocket validation is based on four different methods like prox4, prox5.5, fpocket, and DoGSite. http:// pockd rug. rpbs. univ-paris-dider ot. fr/ cgi-bin/ index. py? page= home.

Protein (Novel Corona-Virus Spike Receptor-Binding Domain Complexes with Its Receptor ACE2) has two chains: chain-A (Angiotensin-converting enzyme2) and Chain-B (Spike glycoprotein) with ligands 1xNAG, 1 × ZN. Models are built based on the target-template alignment using ProMod3 (Li 2016) . As shown in Fig. 1 , the homology modeling is carried out via x-ray diffraction method with a resolution value 2.50A0. This complements the GMQE score (0.99) which estimates the accuracy of the tertiary structure of the model. The QMEAN score is -0.45. It has torsion value as -0.89 and 0.47 as salvation .

SARS-CoV-2 NSP10/NSP16 Methyl-transferase in a Complex with SAM has two chains: chain A(2′-Omethyltransferase) Chain B (Non-structural protein 10) including ligands 1 × 8nk (7-methylguanosine 5′-diphosphate), 1 × gta (p1-7-methylguanosine-p3-adenosine-5′,5′triphosphate), 1 × sah(s-adenosyl-l-homocysteine), 2 × zn (zinc ion) model was built via x-ray diffraction with resolution power 2.10A0 it has QSQE score 0.91 which estimates the accuracy GMQE with sore 0.97.the QMEAN score is -1.35. Torsion and Solvation of the residue are − 0.60 and − 1.82 respectively (Maurya et al. 2020 

MM/GBSA is employed using HAWKDOCK to predict the binding free energy.MM-GBSA (Molecular mechanics-Generalized Born Model and Solvent Accessibility) we approached computational method to calculate the absolute binding affinities between ligand and receptor (Millet and Whittaker 2018) . Various features related to binding free energy like VdW: Vander Waal potentials, ELE: Electrostatic potentials, GB: Polar Solvation free energies are evaluated by the Generalized Born model, SA: Nonpolar contribution to the solvation free energy calculated by an empirical model, TOTAL: Final estimated binding free energy is 2.84 kcal/mol calculated from the terms above. (The lower binding free energy, the more critical residues) (Zhang et al. 2020; Mohan et al. 2021) . Hence, the scoring depicts the structure is critical residue and sensible structure as a target. The given Fig. 2 shows the prediction of (MMG-BSA) with binding free energy features and values. Top 5 binding affinities per residue in both receptor and ligand with their binding energies are as follows ( 

The protein-protein docking was performed on a real case from a CAPRI experiment using pydock. Finally, RMSD comparison with the real 3-D complex structure will be given to check the results of our "simulated CAPRI experiment". The result tabular form analyzes for protein-protein interactions of structure (Ni et al. 2020) . The docking poses are visualized using PYMOL (Fig. 4, Table 1 ).

Given two individual structures i.e. Spike receptor binding domain with its ACE2 receptor submitted as receptor and NSP10/NSP16 as Ligand (Panda et al. 2020) . The predicted binding mode possesses an acceptable accuracy according to the CAPRI criteria. The top performed binding sites of receptor and ligand using HDOCK are predicted (Fig. 5) . The dashboard also describes the summary of top-5 (Row-1: Rank of the top 5 docked models, Row-2: Docking score, Row-3: Ligand RMSD of modeled structures using homologous templates (Phan 2020) . Moreover, a quality report of the docking structure is also given as shown in Fig. 4 : Quality report of receptor:

LG score: 5.322(very good) Maxsub: 0.377(good) and Quality report of ligand:

LG score: 4.664 (good) Maxsub: 0.356 (good). Thusly, the census of receptor and ligand docking poses are cohesive with good quality report which is a consequence of stable complex when bound to each other with complementarity (Porter et al. 2019; Weng et al. 2020 ).

Drug ability, the target identification is major aspect of drug discovery. We have predicted the ability of protein pocket to bind with high affinity as a drug candidate Pock Drug. It is efficient for receptor and ligand as well which are estimated pockets using several thresholds (Pallara et al. 2017) . It scrutinizes various intermolecular interactions with scoring potential like volume of the whole structure, hydrophobicity, polar residues, Aromatic residues, otyr atom, and number of residue (Sharma et al. 2020; Singhal 2020) . The structures are put down in hierarchy according to their peak drug capability model. The descriptors are Volume of convex hull (Å3, Vol. hull), Hydrophobicity based properties of residues (Hydroph), Polar residues, Aromatic residues, Otyr atom (Toor et al. 2020) , Number of pocket residues (Nb Res), Drug Probability, Standard deviation. Therefore, pocket drug ability investigations represent a key step of compound clinical progression projects (Waterhouse et al. 2018) . Earlier, Verma and Shakya (2021) , shown binding affinity of receptor and Ligand are function of binding pose (Table 2) . Post translational modification is essential step in the flow of genetic information. Several proteins and enzyme remain associated with post translational modification of host RNA and vital RNA as well. Viruses mainly RNA viruses evolved with several enzymes including 2-O-MTase allow an mRNA modification to bypass could result in potential antiviral drug including against nSARS-CoV-2 (Wang et al. 2015; Krafcikova et al. 2020 ).

The fundamental role of the spike protein in infectivity indicates that it is a major target for immunogenic development. The stimulation of Nsp16 methyltransferase activity by Nsp10 may be a common mechanism for Coronaviruses. Our study results have implications for design and discover specific anticorona viral medicine to regulate the virus infection. The present study describes docking poses, MMGBSA, binding sites and drug ability that is aimed towards molecular associations certain to infection and replication. Therefore, these findings would help to control the exploiting pandemic condition. Thus the receptor and ligand form a stable complex according to the quality report of docked poses. A major try can turn out prospering medicine and vaccinations against prevailing and potential future SARS-CoV-2 infections to reduce the tremendous impact on human life.

SARS-CoV-2 RNA Dependent RNA polymerase (RdRp)-a drug repurposing study

Molecular basis for drug repurposing to study the interface of the S protein in SARS-CoV-2 and human ACE2 through docking, characterization, and molecular dynamics for natural drug candidates

The SARS-CoV-2 receptor, ACE-2, is expressed on many different cell types: implications for ACE-inhibitor-and angiotensin II receptor blocker-based cardiovascular therapies

Evidence of the COVID-19 virus targeting the CNS: tissue distribution, host-virus interaction, and proposed neurotropic mechanisms

Prediction of potential inhibitors of the dimeric SARS-CoV2 main proteinase through the MM/GBSA approach

Genomics insights of SARS-CoV-2 (COVID-19) into target-based drug discovery

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

Identification of SARS-CoV-2 cell entry inhibitors by drug repurposing using in silico structure-based virtual screening approach

Interaction of SARS-CoV-2 and other coronavirus with ACE (angiotensin-converting enzyme)-2 as their main receptor: therapeutic implications

A novel angiotensin-converting enzyme-related carboxypeptidase (ACE2) converts angiotensin I to angiotensin 1-9

Identification of a new potential SARS-COV-2 RNAdependent RNA polymerase inhibitor via combining fragmentbased drug design, docking, molecular dynamics, and MM-PBSA calculations

The species severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

The origin, transmission and clinical therapies on coronavirus disease 2019 (COVID-19) outbreak-an update on the status

Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis

SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor

PockDrug-server: a new web server for predicting pocket druggability on holo and apo proteins

Novel coronavirus pandemic: a clinical overview

Structural analysis of the SARS-CoV-2 methyltransferase complex involved in RNA cap creation bound to sinefungin

Structure-based drug designing towards the identification of potential anti-viral for COVID-19 by targeting endoribonuclease NSP15

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): the epidemic and the challenges

Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B beta-coronaviruses

Structure, function, and evolution of coronavirus spike proteins

Structure of SARS coronavirus spike receptor-binding domain complexed with receptor

Physiological and pathological regulation of ACE2, the SARS-CoV-2 receptor

Coronavirus disease pandemic (COVID-19): challenges and a global perspective

Structure-based drug designing for potential antiviral activity of selected natural products from ayurveda against SARS-CoV-2 spike glycoprotein and its cellular receptor

Physiological and molecular triggers for SARS-CoV membrane fusion and entry into host cells

SARS-CoV-2 in environmental perspective: occurrence, persistence, surveillance, inactivation and challenges

Role of angiotensin-converting enzyme 2 (ACE2) in COVID-19

pyDock scoring for the new modeling challenges in docking: protein-peptide, homo-multimers, and domain-domain interactions

Structure-based drug designing and immunoinformatics approach for SARS-CoV-2

Genetic diversity and evolution of SARS-CoV-2

What method to use for protein-protein docking?

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2): a global pandemic and treatment strategies

A review of coronavirus disease-2019 (COVID-19)

Computational drug re-purposing targeting the spike glycoprotein of SARS-CoV-2 as an effective strategy to neutralize COVID-19

Genetic variation in chemokines receptor 5 gene and course of HIV infection; review on genetics and immunological aspect

LRP-1 mediated endocytosis of EFE across the blood-brain barrier protein-protein interaction and molecular dynamics analysis

Coronavirus nsp10/nsp16 methyltransferase can be targeted by nsp10-derived peptide in vitro and in vivo to reduce replication and pathogenesis

SWISS-MODEL: homology modelling of protein structures and complexes

Hawk-Dock: a web server to predict and analyze the protein-protein complex based on computational docking and MM/GBSA

Comprehensive evaluation of fourteen docking programs on protein-peptide complexes

Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission

Biological, clinical and epidemiological features of COVID-19, SARS and MERS and AutoDock simulation of ACE2

The author would like to acknowledge Depart-