Molecular Mechanics Poisson-Boltzmann Research Articles

Prevention of SARS-CoV-2 Proliferation with a Novel and Potent Main Protease Inhibitor by Docking, ADMET, MM-PBSA, and Molecular Dynamics Simulation

COVID-19 is the last disease caused by SARS-CoV-2 associated with a severe immune response and lung damage. The main protease (Mpro) has a vital role in SARS-CoV-2 proliferation. Moreover, humans lack homologous Mpro, which makes the Mproa suitable drug target for the development of SARS-CoV-2 drugs. The purchasable L5000 library (Selleckchem Inc) includes 99,040 compounds that were used for virtual screening. After molecular docking and ADME studies, we selected a compound (WAY-604395) with a potent binding affinity to the Mproactive site and acceptable ADME properties compared to the reference drug (nelfinavir). Molecular dynamics (MD) simulation outcomes have proved that the Mpro-WAY604395 complex possesses a considerable value of flexibility, stability, compactness and binding energy. Our Molecular Mechanics Poisson–Boltzmann Surface Area (MM-PBSA) calculation demonstrates that WAY-604395 is more potent ([Formula: see text]272.19[Formula: see text]kcal mol[Formula: see text]) in comparison with nelfinavir ([Formula: see text]173.39[Formula: see text]kcal[Formula: see text]mol[Formula: see text]) against SARS-CoV-2 Mpro. In conclusion, we suggest that WAY-604395 has the potential for the treatment of SARS-CoV-2 by inhibition of the Mpro.

Natural compounds as potential inhibitors of SARS-CoV-2 main protease

Objective: To explore natural compounds as potential inhibitors against main protease (Mpro) of SARS-CoV-2. Methods: In the current study, systematic molecular docking analysis was conducted using AutoDock 4.2 to determine the binding affinities and interactions between natural compounds and Mpro. Selected natural compounds were further validated using a combination of molecular dynamic (MD) simulations and molecular mechanic Poisson-Boltzmann surface area (MM/PBSA) free energy calculations. Results: Out of twenty natural compounds, four natural metabolites namely, amentoflavone, guggulsterone, puerarin, and piperine were found to have strong interaction with Mpro of SARS-CoV-2 based on docking analysis. During MD simulations, all four natural compounds bound to Mpro at 50 ns and MM/G/P/BSA free energy calculations showed that all four shortlisted ligands had stable and favorable energies with strong binding to Mpro protein. Conclusions: Guggulsterone is a potential inhibitor of COVID- 19 main protease Mpro. Further in vitro and pre-clinical studies are needed.

Computational drug repurposing study elucidating simultaneous inhibition of entry and replication of novel corona virus by Grazoprevir

Outcomes of various clinical studies for the coronavirus disease 2019 (COVID-19) treatment indicated that the drug acts via inhibition of multiple pathways (targets) is likely to be more successful and promising. Keeping this hypothesis intact, the present study describes for the first-time, Grazoprevir, an FDA approved anti-viral drug primarily approved for Hepatitis C Virus (HCV), mediated multiple pathway control via synergistic inhibition of viral entry targeting host cell Angiotensin-Converting Enzyme 2 (ACE-2)/transmembrane serine protease 2 (TMPRSS2) and viral replication targeting RNA-dependent RNA polymerase (RdRP). Molecular modeling followed by in-depth structural analysis clearly demonstrated that Grazoprevir interacts with the key residues of these targets. Futher, Molecular Dynamics (MD) simulations showed stability and burial of key residues after the complex formation. Finally, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) analysis identified the governing force of drug-receptor interactions and stability. Thus, we believe that Grazoprevir could be an effective therapeutics for the treatment of the COVID-19 pandemic with a promise of unlikely drug resistance owing to multiple inhibitions of eukaryotic and viral proteins, thus warrants further clinical studies.

A Molecular Modeling Approach to Identify Potential Antileishmanial Compounds Against the Cell Division Cycle (cdc)-2-Related Kinase 12 (CRK12) Receptor of Leishmania donovani.

The huge burden of leishmaniasis caused by the trypanosomatid protozoan parasite Leishmania is well known. This illness was included in the list of neglected tropical diseases targeted for elimination by the World Health Organization. However, the increasing evidence of resistance to existing antimonial drugs has made the eradication of the disease difficult to achieve, thus warranting the search for new drug targets. We report here studies that used computational methods to identify inhibitors of receptors from natural products. The cell division cycle-2-related kinase 12 (CRK12) receptor is a plausible drug target against Leishmania donovani. This study modelled the 3D molecular structure of the L. donovani CRK12 (LdCRK12) and screened for small molecules with potential inhibitory activity from African flora. An integrated library of 7722 African natural product-derived compounds and known inhibitors were screened against the LdCRK12 using AutoDock Vina after performing energy minimization with GROMACS 2018. Four natural products, namely sesamin (NANPDB1649), methyl ellagic acid (NANPDB1406), stylopine (NANPDB2581), and sennecicannabine (NANPDB6446) were found to be potential LdCRK12 inhibitory molecules. The molecular docking studies revealed two compounds NANPDB1406 and NANPDB2581 with binding affinities of −9.5 and −9.2 kcal/mol, respectively, against LdCRK12 which were higher than those of the known inhibitors and drugs, including GSK3186899, amphotericin B, miltefosine, and paromomycin. All the four compounds were predicted to have inhibitory constant (Ki) values ranging from 0.108 to 0.587 μM. NANPDB2581, NANPDB1649 and NANPDB1406 were also predicted as antileishmanial with Pa and Pi values of 0.415 and 0.043, 0.391 and 0.052, and 0.351 and 0.071, respectively. Molecular dynamics simulations coupled with molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) computations reinforced their good binding mechanisms. Most compounds were observed to bind in the ATP binding pocket of the kinase domain. Lys488 was predicted as a key residue critical for ligand binding in the ATP binding pocket of the LdCRK12. The molecules were pharmacologically profiled as druglike with inconsequential toxicity. The identified molecules have scaffolds that could form the backbone for fragment-based drug design of novel leishmanicides but warrant further studies to evaluate their therapeutic potential.

9-Arylimino noscapinoids as potent tubulin binding anticancer agent: chemical synthesis and cellular evaluation against breast tumour cells

ABSTRACT A library of 9-arylimino derivatives of noscapine was developed by coupling of Schiff base containing imine groups. Virtual screening using molecular docking with tubulin revealed three molecules, 12–14 that bind with high affinity. An improved predicted free energy of binding (FEB) of −5.390, −6.506 and −6.679 kcal/mol for the molecules 12–14 was found compared to noscapine (−5.135 kcal/mol). Furthermore, molecular dynamics simulation in combination with Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) revealed robust binding free energy of −166.03, −169.75 and −170.63 kcal/mol for the molecules 12, 13 and 14, respectively. These derivatives were strategically synthesized and experimentally validated for their anticancer activity. Tubulin binding assay revealed substantial binding of molecules 12–14 with purified tubulin. Further, their anticancer activity was demonstrated using two cancer cell lines (MCF-7 and MDAMB-231) and a panel of primary breast tumour cells. All these derivatives inhibited cellular proliferation in all the cancer cells that ranged between 30.1 and 5.8 µM, which is 1.7 to 7.52 fold lower than that of noscapine. Further, these novel derivatives arrest cell cycle in the G2/M-phase followed by induction of apoptosis. Thus, 9-arylimino noscapinoids 12–14 have a great potential to be a novel therapeutic agent for breast cancers.

In silico drug designing for COVID-19: an approach of high-throughput virtual screening, molecular, and essential dynamics simulations

Severe acute respiratory syndrome-coronavirus2 (SARS-CoV2), a new coronavirus has emerged in Wuhan city of China, December 2019 causing pneumonia named Coronavirus disease-19 (COVID-19), which has spread to the entire world. By January 2021, number of confirmed cumulative cases crossed ∼104 million worldwide. Till date, no effective treatment or drug is available for this virus. Availability of X-ray structures of SARS-CoV2 main protease (Mpro) provides the potential opportunity for structure-based drug designing. Here, we have made an attempt to do computational drug design by targeting main protease of SARS-CoV2. High-throughput virtual screening of million molecules and natural compounds databases were performed followed by docking. After that, the protein-ligand complexes were optimized and rescoring of binding energies were accomplished through molecular dynamics simulation and Molecular mechanics Poisson Boltzmann surface area approaches, respectively. In addition, conformational effect of various ligands on protein was also examined through essential dynamics simulation. Three compounds namely ZINC14732869, ZINC19774413, and ZINC19774479 were finally filtered that displayed better binding affinities than N3 (known) inhibitor and formed conformationally stable complexes. Hence, the current study features the potential novel inhibitors against main protease of SARS-CoV2 which might provide an effective therapeutic strategy against COVID-19. Communicated by Ramaswamy H. Sarma

A molecular docking study of potential inhibitors and repurposed drugs against SARS-CoV-2 main protease enzyme

COVID-19 has affected millions of people. Although many drugs are in use to combat disease, there is not any sufficient treatment yet. Having critical role in propagation of the novel coronavirus (SARS-CoV-2) works Main Protease up into a significant drug target. We have performed a molecular docking study to define possible inhibitor candidates against SARS-CoV-2 Main Protease enzyme. Besides docking Remdesivir, Ribavirin, Chloroquine and 28 other antiviral inhibitors (totally 31 inhibitors) to Main Protease enzyme, we have also performed a molecular docking study of 2177 ligands, which are used against Main Protease for the first time by using molecular docking program Autodock4. All ligands were successfully docked into Main Protease enzyme binding site. Among all ligands, EY16 coded ligand which previously used as EBNA1-DNA binding blocker candidate showed the best score for Main Protease with a binding free energy of −10.83 ​kcal/mol which was also lower than re-docking score of N3 ligand (−10.72 ​kcal/mol) contained in crystal structure of Main Protease. After analyzing the docking modes and docking scores we have found that our ligands have better binding free energy values than the inhibitors in use of treatment. We believe that further studies such as molecular dynamics or Molecular Mechanic Poisson Boltzmann Surface Area studies can make contribution that is more exhaustive to the docking results.

Molecular Modelling Studies on Pyrazole Derivatives for the Design of Potent Rearranged during Transfection Kinase Inhibitors.

RET (rearranged during transfection) kinase, one of the receptor tyrosine kinases, plays a crucial role in the development of the human nervous system. It is also involved in various cell signaling networks responsible for the normal cell division, growth, migration, and survival. Previously reported clinical studies revealed that deregulation or aberrant activation of RET signaling can cause several types of human cancer. For example, medullary thyroid carcinoma (MTC) and multiple endocrine neoplasia (MEN2A, MEN2B) occur due to sporadic mutation or germline RET mutation. A number of RET kinase inhibitors have been approved by the FDA for the treatment of cancer, such as cabozantinib, vandetanib, lenvatinib, and sorafenib. However, each of these drugs is a multikinase inhibitor. Hence, RET is an important therapeutic target for cancer drug design. In this work, we have performed various molecular modelling studies, such as molecular docking and dynamics simulation for the most active compound of the pyrazole series as RET kinase inhibitors. Furthermore, molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) free energy calculation and 3-dimensional quantitative structure–activity relationship (3D-QSAR) were performed using g_mmpbsa and SYBYL-X 2.1 package. The results of this study revealed the crucial binding site residues at the active site of RET kinase and contour map analysis showed important structural characteristics for the design of new highly active inhibitors. Therefore, we have designed ten RET kinase inhibitors, which showed higher inhibitory activity than the most active compound of the series. The results of our study provide insights to design more potent and selective RET kinase inhibitors.

Theoretical Evaluation of Novel Thermolysin Inhibitors from Bacillus thermoproteolyticus. Possible Antibacterial Agents.

The search for new antibacterial agents that could decrease bacterial resistance is a subject in continuous development. Gram-negative and Gram-positive bacteria possess a group of metalloproteins belonging to the MEROPS peptidase (M4) family, which is the main virulence factor of these bacteria. In this work, we used the previous results of a computational biochemistry protocol of a series of ligands designed in silico using thermolysin as a model for the search of antihypertensive agents. Here, thermolysin from Bacillus thermoproteolyticus, a metalloprotein of the M4 family, was used to determine the most promising candidate as an antibacterial agent. Our results from docking, molecular dynamics simulation, molecular mechanics Poisson–Boltzmann (MM-PBSA) method, ligand efficiency, and ADME-Tox properties (Absorption, Distribution, Metabolism, Excretion, and Toxicity) indicate that the designed ligands were adequately oriented in the thermolysin active site. The Lig783, Lig2177, and Lig3444 compounds showed the best dynamic behavior; however, from the ADME-Tox calculated properties, Lig783 was selected as the unique antibacterial agent candidate amongst the designed ligands.

Interactions mechanism of commonly used drugs for the treatment of Covid-19

In this study conformation analysis of seven drugs commonly used in the treatment of COVID-19 was performed. The most stable conformers of the drug molecules were used as initial data for docking analysis. Using the Cavityplus program, the probable most active binding sites of both apo and holo forms of COVID-19 main protease enzyme (Mpro) and spike glycoprotein of SARSCoV-2 receptors were determined. The interaction mechanisms of the 7 FDA approved drugs (arbidol, colchicine, dexamethasone, favipiravir, galidesivir, hydroxychloroquine, remdesivir) were examined using the AutoDock Vina program. The six of the seven drugs were found to be more stable in binding to apo form of COVID-19 Mpro and spike glycoprotein. Moreover, a set of molecular mechanics (MM) Poisson-Boltzmann (PB) surface area (SA) calculations on the investigated drugs-protein systems were performed and the estimated binding free energy of remdesivir and the apo form of Mpro system was found to be the best. The interaction results of FDA drugs with the apo form of COVID-19 Mpro and spike glycoprotein can play an important role for the treatment of COVID-19.
 
 KEY WORDS: COVID-19, Drugs, Molecular modelling, Conformational analysis, Molecular docking
 
 Bull. Chem. Soc. Ethiop. 2020, 34(3), 613-623.
 DOI: https://dx.doi.org/10.4314/bcse.v34i3.16

Elucidating Interactions Between SARS-CoV-2 Trimeric Spike Protein and ACE2 Using Homology Modeling and Molecular Dynamics Simulations.

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19). As of October 21, 2020, more than 41.4 million confirmed cases and 1.1 million deaths have been reported. Thus, it is immensely important to develop drugs and vaccines to combat COVID-19. The spike protein present on the outer surface of the virion plays a major role in viral infection by binding to receptor proteins present on the outer membrane of host cells, triggering membrane fusion and internalization, which enables release of viral ssRNA into the host cell. Understanding the interactions between the SARS-CoV-2 trimeric spike protein and its host cell receptor protein, angiotensin converting enzyme 2 (ACE2), is important for developing drugs and vaccines to prevent and treat COVID-19. Several crystal structures of partial and mutant SARS-CoV-2 spike proteins have been reported; however, an atomistic structure of the wild-type SARS-CoV-2 trimeric spike protein complexed with ACE2 is not yet available. Therefore, in our study, homology modeling was used to build the trimeric form of the spike protein complexed with human ACE2, followed by all-atom molecular dynamics simulations to elucidate interactions at the interface between the spike protein and ACE2. Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) and in silico alanine scanning were employed to characterize the interacting residues at the interface. Twenty interacting residues in the spike protein were identified that are likely to be responsible for tightly binding to ACE2, of which five residues (Val445, Thr478, Gly485, Phe490, and Ser494) were not reported in the crystal structure of the truncated spike protein receptor binding domain (RBD) complexed with ACE2. These data indicate that the interactions between ACE2 and the tertiary structure of the full-length spike protein trimer are different from those between ACE2 and the truncated monomer of the spike protein RBD. These findings could facilitate the development of drugs and vaccines to prevent SARS-CoV-2 infection and combat COVID-19.

Scrutiny of the mechanism of β-amyloid protein captures HSV-2 to protect the brain infection

Alzheimer’s disease (AD) is an age-related neurodegenerative disorder. β-amyloid protein (Aβ) is the key protein which involved in AD. But the physiological function of Aβ is needed to be investigated. Many experimental studies have shown that Aβ could bind to glycoproteins D (gD) on the surface of the herpes virus. However the mechanism is still unclear. In the present study, we elucidate the molecular mechanism of the interaction between Aβ and gD of herpes simplex virus type 2 (HSV-2) by molecular docking and molecular dynamics simulation. Molecular dynamics simulations displayed that Aβ could stably bind to the HSV-2 gD owing to the presence of several interactions. Analysis binding free energy by molecular mechanics Poisson–Boltzmann surface area (MM–PBSA) method revealed that hot residues including Glu3, Glu11, Glu22 and Ala42 of Aβ1-42 were involved in binding with HSV-2 gD. Thus, the HSV-2 gD can be entrapped by Aβ which will be utilized for prevent and therapy of AD in future.

Multi-step molecular docking and dynamics simulation-based screening of large antiviral specific chemical libraries for identification of Nipah virus glycoprotein inhibitors

Nipah virus (NiV) infections are highly contagious and can cause severe febrile encephalitis. An outbreak of NiV infection has reported high mortality rates in Southeast Asian countries including Bangladesh, East Timor, Malaysia, Papua New Guinea, Vietnam, Cambodia, Indonesia, Madagascar, Philippines, Thailand and India. Considering the high risk for an epidemic outbreak, the World Health Organization (WHO) declared NiV as an emerging priority pathogen. However, there are no effective therapeutics or any FDA approved drugs available for the treatment of this infection. Among the known nine proteins of NiV, glycoprotein plays an important role in initiating the entry of viruses and attaching to the host cell receptors. Herein, three antiviral databases consisting of 79,892 chemical entities have been computationally screened against NiV glycoprotein (NiV-G). Particularly, multi-step molecular docking followed by extensive molecular binding interactions analyses, binding free energy estimation, in silico pharmacokinetics, synthetic accessibility and toxicity profile evaluations have been carried out for initial identification of potential NiV-G inhibitors. Further, molecular dynamics (MD) simulation has been performed to understand the dynamic properties of NiV-G protein-bound with proposed five inhibitors (G1-G5) and their interactions behavior, and any conformational changes in NiV-G protein during simulations. Moreover, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) based binding free energies (∆G) has been calculated from all MD simulation trajectories to understand the energy contribution of each proposed compound in maintaining and stabilizing the complex binding interactions with NiV-G protein. Proposed compounds showed high negative ∆G values ranging from −166.246 to −226.652 kJ/mol indicating a strong affinity towards the NiV-G protein.

Identification of FDA approved drugs and nucleoside analogues as potential SARS-CoV-2 A1pp domain inhibitor: An in silico study

Coronaviruses are known to infect respiratory tract and intestine. These viruses possess highly conserved viral macro domain A1pp having adenosine diphosphate (ADP)-ribose binding and phosphatase activity sites. A1pp inhibits adenosine diphosphate (ADP)-ribosylation in the host and promotes viral infection and pathogenesis. We performed in silico screening of FDA approved drugs and nucleoside analogue library against the recently reported crystal structure of SARS-CoV-2 A1pp domain. Docking scores and interaction profile analyses exhibited strong binding affinity of eleven FDA approved drugs and five nucleoside analogues NA1 (−13.84), nadide (−13.65), citicholine (−13.54), NA2 (−12.42), and NA3 (−12.27). The lead compound NA1 exhibited significant hydrogen bonding and hydrophobic interaction at the natural substrate binding site. The root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent accessible surface (SASA), hydrogen bond formation, principle component analysis, and free energy landscape calculations for NA1 bound protein displayed stable complex formation in 100 ns molecular dynamics simulation, compared to unbound macro domain and natural substrate adenosine-5-diphosphoribose bound macro domain that served as a positive control. The molecular mechanics Poisson–Boltzmann surface area analysis of NA1 demonstrated binding free energy of −175.978 ± 0.401 kJ/mol in comparison to natural substrate which had binding free energy of −133.403 ± 14.103 kJ/mol. In silico analysis by modelling tool ADMET and prediction of biological activity of these compounds further validated them as putative therapeutic molecules against SARS-CoV-2. Taken together, this study offers NA1 as a lead SARS-CoV-2 A1pp domain inhibitor for future testing and development as therapeutics against human coronavirus.

Detailed Characterization of the Cooperative Binding of Piperine with Heat Shock Protein 70 by Molecular Biophysical Approaches.

In this work, for the first time, details of the complex formed by heat shock protein 70 (HSP70) independent nucleotide binding domain (NBD) and piperine were characterized through experimental and computational molecular biophysical methods. Fluorescence spectroscopy results revealed positive cooperativity between the two binding sites. Circular dichroism identified secondary conformational changes. Molecular dynamics along with molecular mechanics Poisson Boltzmann surface area (MM/PBSA) reinforced the positive cooperativity, showing that the affinity of piperine for NBD increased when piperine occupied both binding sites instead of one. The spontaneity of the complexation was demonstrated through the Gibbs free energy (∆G < 0 kJ/mol) for different temperatures obtained experimentally by van’t Hoff analysis and computationally by umbrella sampling with the potential of mean force profile. Furthermore, the mean forces which drove the complexation were disclosed by van’t Hoff and MM/PBSA as being the non-specific interactions. In conclusion, the work revealed characteristics of NBD and piperine interaction, which may support further drug discover studies.

High-throughput virtual screening of novel potent inhibitor(s) for Human Vanin-1 enzyme

Vanin-1 (VNN1) is a glycosylphosphatidylinositol (GPI)-anchored ectoenzyme which hydrolyzes pantetheine to pantothenic acid and cysteamine. It has emerged as a promising drug target for many human diseases associated with oxidative stress and inflammatory pathways. In the present study we used structure-based virtual screening approach for the identification of small molecule inhibitors of vanin-1. A chemical library consisting of natural compounds, synthetic compounds and RRV analogs were screened for drug-like molecules. The filtered molecules were subjected to molecular docking studies. Three potential hits—ZINC04073864 (Natural compound), CID227017 (synthetic compound) and CID129558381 (RRV analog)—were identified for the target enzyme. The molecules form good number of hydrogen bonds with the catalytic residues such as Glu79, Lys178 and Cys211. The apo-VNN1 and VNN1-ligand complexes were subjected to molecular dynamics (MD) simulation for 30 ns. The geometric properties such as root mean square deviation, radius of gyration, solvent accessible surface area, number of hydrogen bonds and the distance between the catalytic triad residues-Glu79, Lys178 and Cys211 were altered upon binding of the compounds. Essential dynamics and entropic studies further confirmed that the fluctuations in VNN1 decrease upon binding of the compounds. The lead molecules were stable throughout the simulation time period. Molecular Mechanics Poisson–Boltzmann Surface Area (MM/PBSA) studies showed that Van der Waals interaction energy contributes significantly to the total binding free energy. Thus, our study reveals three lead molecules—ZINC04073864, CID227017 and CID129558381 as potential inhibitors of Vanin-1 which can be validated through further studies. Communicated by Ramaswamy H. Sarma

Ligand-Based Virtual Screening, Molecular Docking, Molecular Dynamics, and MM-PBSA Calculations towards the Identification of Potential Novel Ricin Inhibitors.

Ricin is a toxin found in the castor seeds and listed as a chemical weapon by the Chemical Weapons Convention (CWC) due to its high toxicity combined with the easiness of obtention and lack of available antidotes. The relatively frequent episodes of usage or attempting to use ricin in terrorist attacks reinforce the urge to develop an antidote for this toxin. In this sense, we selected in this work the current RTA (ricin catalytic subunit) inhibitor with the best experimental performance, as a reference molecule for virtual screening in the PubChem database. The selected molecules were then evaluated through docking studies, followed by drug-likeness investigation, molecular dynamics simulations and Molecular Mechanics Poisson–Boltzmann Surface Area (MM-PBSA) calculations. In every step, the selection of molecules was mainly based on their ability to occupy both the active and secondary sites of RTA, which are located right next to each other, but are not simultaneously occupied by the current RTA inhibitors. Results show that the three PubChem compounds 18309602, 18498053, and 136023163 presented better overall results than the reference molecule itself, showing up as new hits for the RTA inhibition, and encouraging further experimental evaluation.

Computational Investigation of Structural Dynamics of SARS-CoV-2 Methyltransferase-Stimulatory Factor Heterodimer nsp16/nsp10 Bound to the Cofactor SAM.

Recently, a highly contagious novel coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, has emerged, posing a global threat to public health. Identifying a potential target and developing vaccines or antiviral drugs is an urgent demand in the absence of approved therapeutic agents. The 5′-capping mechanism of eukaryotic mRNA and some viruses such as coronaviruses (CoVs) are essential for maintaining the RNA stability and protein translation in the virus. SARS-CoV-2 encodes S-adenosyl-L-methionine (SAM) dependent methyltransferase (MTase) enzyme characterized by nsp16 (2′-O-MTase) for generating the capped structure. The present study highlights the binding mechanism of nsp16 and nsp10 to identify the role of nsp10 in MTase activity. Furthermore, we investigated the conformational dynamics and energetics behind the binding of SAM to nsp16 and nsp16/nsp10 heterodimer by employing molecular dynamics simulations in conjunction with the Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) method. We observed from our simulations that the presence of nsp10 increases the favorable van der Waals and electrostatic interactions between SAM and nsp16. Thus, nsp10 acts as a stimulator for the strong binding of SAM to nsp16. The hydrophobic interactions were predominately identified for the nsp16-nsp10 interactions. Also, the stable hydrogen bonds between Ala83 (nsp16) and Tyr96 (nsp10), and between Gln87 (nsp16) and Leu45 (nsp10) play a vital role in the dimerization of nsp16 and nsp10. Besides, Computational Alanine Scanning (CAS) mutagenesis was performed, which revealed hotspot mutants, namely I40A, V104A, and R86A for the dimer association. Hence, the dimer interface of nsp16/nsp10 could also be a potential target in retarding the 2′-O-MTase activity in SARS-CoV-2. Overall, our study provides a comprehensive understanding of the dynamic and thermodynamic process of binding nsp16 and nsp10 that will contribute to the novel design of peptide inhibitors based on nsp16.

Revealing the Inhibition Mechanism of RNA-Dependent RNA Polymerase (RdRp) of SARS-CoV-2 by Remdesivir and Nucleotide Analogues: A Molecular Dynamics Simulation Study.

Antiviral drug therapy against SARS-CoV-2 is not yet established and posing a serious global health issue. Remdesivir is the first antiviral compound approved by the US FDA for the SARS-CoV-2 treatment for emergency use, targeting RNA-dependent RNA polymerase (RdRp) enzyme. In this work, we have examined the action of remdesivir and other two ligands screened from the library of nucleotide analogues using docking and molecular dynamics (MD) simulation studies. The MD simulations have been performed for all the ligand-bound RdRp complexes for the 30 ns time scale. This is one of the earlier reports to perform the MD simulations studies using the SARS-CoV-2 RdRp crystal structure (PDB ID 7BTF). The MD trajectories were analyzed and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations were performed to calculate the binding free energy. The binding energy data reveal that compound-17 (-59.6 kcal/mol) binds more strongly as compared to compound-8 (-46.3 kcal/mol) and remdesivir (-29.7 kcal/mol) with RdRp. The detailed analysis of trajectories shows that the remdesivir binds in the catalytic site and forms a hydrogen bond with the catalytic residues from 0 to 0.46 ns. Compound-8 binds in the catalytic site but does not form direct hydrogen bonds with catalytic residues. Compound-17 showed the formation of hydrogen bonds with catalytic residues throughout the simulation process. The MD simulation results such as hydrogen bonding, the center of mass distance analysis, snapshots at a different time interval, and binding energy suggest that compound-17 binds strongly with RdRp of SARS-CoV-2 and has the potential to develop as a new antiviral against COVID-19. Further, the frontier molecular orbital analysis and molecular electrostatic potential (MESP) iso-surface analysis using DFT calculations shed light on the superior binding of compound-17 with RdRp compared to remdesivir and compound-8. The computed as well as the experimentally reported pharmacokinetics and toxicity parameters of compound-17 is encouraging and therefore can be one of the potential candidates for the treatment of COVID-19.

Molecular docking and molecular dynamics simulation of Bacillus thuringiensis dehalogenase against haloacids, haloacetates and chlorpyrifos

The high dependency and surplus use of agrochemical products have liberated enormous quantities of toxic halogenated pollutants into the environment and threaten the well-being of humankind. Herein, this study performed molecular docking, molecular dynamic (MD) simulations, molecular mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculations on the DehH2 from Bacillus thuringiensis, to identify the order of which the enzyme degrades different substrates, haloacids, haloacetate and chlorpyrifos. The study discovered that the DehH2 favored the degradation of haloacids and haloacetates (−3.3 − 4.6 kcal/mol) and formed three hydrogen bonds with Asp125, Arg201 and Lys202. Despite the inconclusive molecular docking result, chlorpyrifos was consistently shown to be the least favored substrate of the DehH2 in MD simulations and MM-PBSA calculations. Results of MD simulations revealed the DehH2-haloacid- (RMSD 0.15 − 0.25 nm) and DehH2-haloacetates (RMSF 0.05 − 0.25 nm) were more stable, with the DehH2-L-2CP complex being the most stable while the least was the DehH2-chlorpyrifos (RMSD 0.295 nm; RMSF 0.05 − 0.59 nm). The Molecular Mechanics Poisson-Boltzmann Surface Area calculations showed the DehH2-L-2CP complex (−24.27 kcal/mol) having the lowest binding energy followed by DehH2-MCA (−22.78 kcal/mol), DehH2-D-2CP (−21.82 kcal/mol), DehH2-3CP (−21.11 kcal/mol), DehH2-2,2-DCP (−18.34 kcal/mol), DehH2-2,3-DCP (−8.34 kcal/mol), DehH2-TCA (−7.62 kcal/mol), while chlorpyrifos was unable to spontaneously bind to DehH2 (+127.16 kcal/mol). In a nutshell, the findings of this study offer valuable insights into the rational tailoring of the DehH2 for expanding its substrate specificity and catalytic activity in the near future. Communicated by Ramaswamy H. Sarma