Ligand-bound Complexes Research Articles

Enhancing Binding Affinity to CGG/CGG Triad: Optimizing Naphthyridine Carbamate Dimer Derivatives with Varied Linker Lengths.

This study examines the binding properties of six naphthyridine carbamate dimer (NCD) derivatives with varying linker lengths to the CGG/CGG triad, a non-canonical DNA structure linked to repeat expansion disorders. By altering the linker length from 2 to 4 methylene groups, we found changes in thermal stability of the ligand-bound complexes while maintaining a consistent 2 : 1 binding stoichiometry. Among the derivatives, CC23 showed superior binding affinity compared to the parent molecule CC33 (NCD). Spectroscopic analyses revealed that linker length influences the conformational equilibrium of NCD derivatives. Thermal melting temperature measurements demonstrated CC23's enhanced thermal stability over CC33. These findings underscore the potential of optimized NCD derivatives, like CC23, as tools to modulate CGG repeat structures, offering insights for therapeutic strategies targeting repeat expansion disorders.

Measuring Ligand-bound Protein Complexes with Proximity Labeling: A Practical Guide.

Identifying the drug-target interactome of small molecule therapeutics is essential for understanding the full pharmacological effects of a compound. These therapies often induce changes within the cellular proteome, leading to unexpected consequences such as changes in the targets complexation state or off-target interactions between the compound and additional proteins. Currently, unbiased target-ID approaches are being used to embark on this task. Here we provide an overview of the strengths and limitations of these methods, and a practical step-by-step protocol for using the BioTAC system to assist with drug target and interactome ID.

Crystal Structures of Lacticaseibacillus 4-Deoxy-L-threo-5-hexosulose-uronate Ketol-isomerase KduI in Complex with Substrate Analogs.

Some probiotics including lactobacilli, colonize host animal cells by targeting glycosaminoglycans (GAGs), such as heparin, located in the extracellular matrix. Recent studies have shown that several lactic acid bacteria degrade GAGs. Here we show the structure/function relationship of Lacticaseibacillus rhamnosus 4-deoxy-L-threo-5-hexosulose-uronate ketol-isomerase (KduI) crucial for metabolism of unsaturated glucuronic acid produced through degradation of GAGs. Crystal structures of ligand-free and bound KduIs were determined by X-ray crystallography and the enzyme was found to consist of six identical subunits and adopt a β-helix as a basic scaffold. Ligands structurally similar to the substrate were bound to the cleft of each enzyme subunit. Several residues located in the cleft interacted with ligands through hydrogen bonds and/or C-C contacts. In addition to substrate analogs, a metal ion coordinated to four residues, His198, His200, Glu205, and His248, in the cleft, and the enzyme activity was significantly inhibited by a chelator, ethylenediaminetetraacetic acid. Site-directed mutants in Arg163, Ile165, Thr184, Thr194, His200, Arg203, Tyr207, Met262, and Tyr269 in the cleft exhibited little enzyme activity, indicating that these residues and the metal ion constituted an active site in the cleft. This is the first report on the active site structure of KduI based on the ligand-bound complex.

Deciphering specificity and cross-reactivity in tachykinin NK1 and NK2 receptors

The tachykinin receptors neurokinin 1 (NK1R) and neurokinin 2 (NK2R) are G protein-coupled receptors that bind preferentially to the natural peptide ligands substance P and neurokinin A, respectively, and have been targets for drug development. Despite sharing a common C-terminal sequence of Phe-X-Gly-Leu-Met-NH2 that helps direct biological function, the peptide ligands exhibit some degree of cross-reactivity toward each other’s non-natural receptor. Here, we investigate the detailed structure-activity relationships of the ligand-bound receptor complexes that underlie both potent activation by the natural ligand and cross-reactivity. We find that the specificity and cross-reactivity of the peptide ligands can be explained by the interactions between the amino acids preceding the FxGLM consensus motif of the bound peptide ligand and two regions of the receptor: the β-hairpin of the extracellular loop 2 (ECL2) and a N-terminal segment leading into transmembrane helix 1. Positively charged sidechains of the ECL2 (R177 of NK1R and K180 of NK2R) are seen to play a vital role in the interaction. The N-terminal positions 1-3 of the peptide ligand are entirely dispensable. Mutated and chimeric receptor and ligand constructs neatly swap around ligand specificity as expected, validating the structure-activity hypotheses presented. These findings will help in developing improved agonists or antagonists for NK1R and NK2R.

How Ligands Interact with the Kinase Hinge.

ATP-competitive kinase inhibitors form hydrogen bond interactions with the kinase hinge region at the adenine binding site. Thus, it is crucial to explore hinge-ligand recognition as part of a rational drug design strategy. Here, harnessing known ligand-bound kinase structures and experimental assay resources, we first created a kinase structure-assay database (KSAD) containing 2705 nM ligand-bound kinase complexes. Then, using KSAD, we systematically investigate hinge-ligand binding patterns using interaction fingerprints, thereby delineating 15 different hydrogen-bond interaction modes. We believe these results will be valuable for de novo drug design and/or scaffold hopping of kinase-targeted drugs.

Bioactivity-Guided Synthesis: In Silico and In Vitro Studies of β-Glucosidase Inhibitors to Cope with Hepatic Cytotoxicity.

The major cause of hyperglycemia can generally be attributed to β-glucosidase as per its involvement in non-alcoholic fatty liver disease. This clinical condition leads to liver carcinoma (HepG2 cancer). The phthalimides and phthalamic acid classes possess inhibitory potential against glucosidase, forming the basis for designing new phthalimide and phthalamic acid analogs to test their ability as potent inhibitors of β-glucosidase. The study also covers in silico (molecular docking and MD simulations) and in vitro (β-glucosidase and HepG2 cancer cell line assays) analyses. The phthalimide and phthalamic acid derivatives were synthesized, followed by spectroscopic characterization. The mechanistic complexities associated with β-glucosidase inhibition were identified via the docking of the synthesized compounds inside the active site of the protein, and the results were analyzed in terms of the best binding energy and appropriate docking pose. The top-ranked compounds were subjected to extensive MD simulation studies to understand the mode of interaction of the synthesized compounds and binding energies, as well as the contribution of individual residues towards binding affinities. Lower RMSD/RMSF values were observed for 2c and 3c, respectively, in the active site, confirming more stabilized, ligand-bound complexes when compared to the free state. An anisotropic network model was used to unravel the role of loop fluctuation in the context of ligand binding and the dynamics that are distinct to the bound and free states, supported by a 3D surface plot. An in vitro study revealed that 1c (IC50 = 1.26 µM) is far better than standard acarbose (2.15 µM), confirming the potential of this compound against the target protein. Given the appreciable potential of the candidate compounds against β-glucosidase, the synthesized compounds were further tested for their cytotoxic activity against hepatic carcinoma on HepG2 cancer cell lines. The cytotoxicity profile of the synthesized compounds was performed against HepG2 cancer cell lines. The resultant IC50 value (0.048 µM) for 3c is better than the standard (thalidomide: IC50 0.053 µM). The results promise the hypothesis that the synthesized compounds might become potential drug candidates, given the fact that the β-glucosidase inhibition of 1c is 40% better than the standard, whereas compound 3c holds more anti-tumor activity (greater than 9%) against the HepG2 cell line than the known drug.

Computational exploration of natural compounds targeting Staphylococcus aureus: inhibiting AgrA promoter binding for antimicrobial intervention

Staphylococcus aureus is a highly virulent nosocomial pathogen that poses a significant threat to individuals exposed to healthcare settings. Due to its sophisticated machinery for producing virulence factors, S. aureus can cause severe and potentially fatal infections in humans. This study focuses on the response regulator AgrA, which plays a crucial role in regulating the production of virulence factors in S. aureus. The objective is to identify natural compounds that can inhibit the binding of AgrA to its promoter site, thus inhibiting the expression of virulence genes. To achieve this, a pharmacophore model was generated using known drugs and applied to screen the ZINC natural product database. The resulting compounds were subjected to molecular docking-based virtual screening against the C-terminal DNA binding domain of AgrA. Three compounds, namely ZINC000077269178, ZINC000051012304, and ZINC000004266026, were shortlisted based on their strong affinity for key residues involved in DNA binding and transcription initiation. Subsequently, the unbound and ligand-bound complexes were subjected to a 200 ns molecular dynamics simulation to assess their conformational stability. Various analyses, including RMSD, RMSF, Rg, SASA, Principal Component Analysis, and Gibbs free energy landscape, were conducted on the simulation trajectory. The RMSD profile indicated similar fluctuations in both bound and unbound structures, while the Rg profile demonstrated the compactness of the protein without any unfolding during the simulation. Furthermore, Principal component analysis revealed that ligand binding reduced the overall atomic motion of the protein whereas free energy landscape suggested the energy variations obtained in complexes. Communicated by Ramaswamy H. Sarma

Gas-phase stability and thermodynamics of ligand-bound, binary complexes of chloramphenicol acetyltransferase reveal negative cooperativity.

The biological role of the bacterial chloramphenicol (Chl)-resistance enzyme, chloramphenicol acetyltransferase (CAT), has seen renewed interest due to the resurgent use of Chl against multi-drug-resistant microbes. This looming threat calls for more rationally designed antibiotic derivatives that have improved antimicrobial properties and reduced toxicity in humans. Herein, we utilize native ion mobility spectrometry-mass spectrometry (IMS-MS) to investigate the gas-phase structure and thermodynamic stability of the type I variant of CAT from Escherichia coli (EcCATI) and several EcCATI:ligand-bound complexes. EcCATI readily binds multiple Chl without incurring significant changes to its gas-phase structure or stability. A non-hydrolyzable acetyl-CoA derivative (S-ethyl-CoA, S-Et-CoA) was used to kinetically trap EcCATI and Chl in a ternary, ligand-bound state (EcCATI:S-Et-CoA:Chl). Using collision-induced unfolding (CIU)-IMS-MS, we find that Chl dissociates from EcCATI:S-Et-CoA:Chl complexes at low collision energies, while S-Et-CoA remains bound to EcCATI even as protein unfolding occurs. Gas-phase binding constants further suggest that EcCATI binds S-Et-CoA more tightly than Chl. Both ligands exhibit negative cooperativity of subsequent ligand binding in their respective binary complexes. While we observe no significant change in structure or stability to EcCATI when bound to either or both ligands, we have elucidated novel gas-phase unfolding and dissociation behavior and provided a foundation for further characterization of alternative substrates and/or inhibitors of EcCATI.

In Vitro and in silico studies to explore potent antidiabetic inhibitor against human pancreatic alpha-amylase from the methanolic extract of the green microalga Chlorella vulgaris

Today’s era and lifestyle have led to a quick rise in cases of diabetes. Diabetes mellitus (DM) has risen to the top of the list of serious diseases and stems from different health disorders. Human pancreatic alpha-amylase (HPA) enzyme plays a critical role in the digestion of carbohydrates, and inhibitors of alpha-amylase have been investigated as a way to slow the absorption of carbohydrates and reduce postprandial (after meal) hyperglycemia in patients with diabetes. Recently algal derivatives have been studied for their potential as a new drug against diabetes and other diseases. The study is aimed to find active biochemical compounds from the methanolic extract of Chlorella vulgaris. The in vitro studies were carried out and the results revealed that methanolic extract from C. vulgaris showed abundant inhibition efficacy of the α-amylase (IC50 of about 2.66 µg/mL) compared to acarbose (IC50 of about 2.85 µg/mL), a standard, commercial inhibitor. All the bioactive compounds from the methanolic extract were identified from the GCMS study and considered for in silico evaluation. Out of 14 bioactive compounds from GCMS, compound C3 showed higher docking energy (-8.3 kcal/mol) compared to other compounds. Subsequently, the comparative molecular dynamic simulation of apo and ligand-bound (compound C3 and acarbose) α-amylase complexes showed overall structural stability for compound C3 at the active site of α-amylase from various MD analyses. Hence, we believe, the bioactive compounds identified from GCMS may assist in diabetic therapeutics. Moreover, the compound C3 identified in this study could be a potential antidiabetic therapeutic inhibitor. Communicated by Ramaswamy H. Sarma

Finding inhibitors and deciphering inhibitor-induced conformational plasticity in the Janus kinase via multiscale simulations

ABSTRACT The Janus kinase (JAK) is a master regulator of the JAK/STAT pathway. Dysregulation of this signalling cascade causes neuroinflammation and autoimmune disorders. Therefore, JAKs have been characterized as an attractive target for developing anti-inflammatory drugs. Nowadays, designing efficient, effective, and specific targeted therapeutics without being cytotoxic has gained interest. We performed the virtual screening of natural products in combination with pharmacological analyses. Subsequently, we performed molecular dynamics simulations to study the stability of the ligand-bound complexes and ligand-induced inactive conformations. Notably, inactive kinases display remarkable conformational plasticity; however, ligand-induced molecular mechanisms of these conformations are still poorly understood. Herein, we performed a free energy landscape analysis to explore the conformational plasticity of the JAK1 kinase. Leonurine, STOCK1N-68642, STOCK1N-82656, and STOCK1N-85809 bound JAK1 exhibited a smooth transition from an active (αC-in) to a completely inactive conformation (αC-out). Ligand binding induces disorders in the αC-helix. Molecular mechanics Poisson Boltzmann surface area (MM/PBSA) calculation suggested three phytochemicals, namely STOCK1N-68642, Epicatechin, and STOCK1N-98615, have higher binding affinity compared to other ligand molecules. The ligand-induced conformational plasticity revealed by our simulations differs significantly from the available crystal structures, which might help in designing allosteric drugs.

In Silico Insights on the Pro-Inflammatory Potential of Polycyclic Aromatic Hydrocarbons and the Prospective Anti-Inflammatory Capacity of Andrographis paniculata Phytocompounds.

Inflammation linked to various diseases is the biological response to certain stimuli. The pro-inflammatory potential of Polycyclic Aromatic Hydrocarbons (PAHs) as potential inducers of inflammation bound to the Toll-like Receptor 4 (TLR4) and the anti-inflammatory capacity of A. paniculata (AP) phytocompounds as prospective inhibitors of the Nuclear Factor Kappa B (NF-κB) p50 transcription factor are investigated via in silico techniques. The molecular docking of the PAHs and AP phytocompounds is performed in AutoDock Vina by calculating their binding energies. The molecular dynamics simulations (MDS) of the apo and ligand-bound complex of the top binding ligands were performed in CABS-flex. The agonists, which included the PAHs indeno(1,2,3-cd)pyrene (IP), and dibenz(a,h)anthracene (DahA), had the highest binding energies of −10 kcal/mol and −9.2 kcal/mol, respectively. The most stable antagonists in the binding site with binding energies to the NF-κB p50 were the AP phytocompounds with −5.6 kcal/mol for ergosterol peroxide and −5.3 kcal/mol for 14-deoxy-14,15-dehydroandrographolide. The MDS of the apo human TLR4 and PAH-bound TLR4, and the apo p50 and the AP phytocompound-bound NF-κB p50 showed minimal fluctuations. These results reveal that IP and DahA are significant inducers of inflammation, whereas ergosterol peroxide and 14-deoxy-14,15-dehydroandrographolide are inhibitors of the NF-κB pathway. Furthermore, the study theorizes that any inflammatory activity induced by PAH can be potentially inhibited by A. paniculata phytocompounds.

Evaluation of xanthene-appended quinoline hybrids as potential leads against antimalarial drug targets.

A series of fused heterocycle xanthene-appended quinoline 6a-n was successfully synthesized with regioselectivity and characterized using IR, 1H NMR, 13C NMR, and mass spectral data. Molecular docking was performed to find the binding efficacy of all these newly synthesized compounds towards thirteen antimalarial drug targets. Molecular dynamics simulation was carried out to predict the stability of the ligand-bound complex in a solvent medium. Blind and site-directed docking with compounds 6a-n against 13 drug targets revealed most of the ligands to have a good binding affinity with the targets. Analysis on the basis of binding energy, binding modalities of the ligands, intermolecular interactions, and pharmacophore, we identified only one of the ligand-receptor complexes to provide better results. Molecular dynamic simulation of the selected receptor-ligand complex revealed that the synthesized compound had a better binding affinity with the receptor than the native ligand complex. Further analysis of the synthesized ligand in the laboratory may prove promising results in the search for potential antimalarial drugs.

Ligand-Dependent Modulation of the Dynamics of Intracellular Loops Dictates Functional Selectivity of 5-HT2AR

The serotonin 2A receptor (5-HT2AR) subtype of the G protein-coupled receptor (GPCR) family is involved in a plethora of neuromodulatory functions (e.g., neurogenesis, sleep, and cognitive processes). 5-HT2AR is the target of pharmacologically distinct classes of ligands, binding of which either activate or inactivate the receptor. Although high-resolution structures of 5-HT2AR as well as several other 5-HT GPCRs provided snapshots of both active and inactive conformational states, these structures, representing a truncated form of the receptor, cannot fully explain the mechanism of conformational transitions during their function. Importantly, biochemical studies have suggested the importance of intracellular loops in receptor functions. In our previous study, a model of the ligand-free form of 5-HT2AR with the third intracellular loop (ICL3) has been meticulously built. Here, we have investigated the functional regulation of 5-HT2AR with intact intracellular loops in ligand-free and five distinct ligand-bound configurations using unbiased atomistic molecular dynamics (MD) simulations. The selected ligands belong to either of the full, partial, or inverse agonist classes, which exert distinct pharmacological responses. We have observed significant structural, dynamic, and thermodynamic differences within ligand-bound complexes. Our results revealed, for the first time, that either activation or inactivation of the receptor upon specific ligand binding is primarily achieved through conformational transitions of its second and third intracellular loops (ICL2 and ICL3). A remarkable allosteric cross-talk between the ligand-binding site and the distal intracellular parts of the receptor, where binding of a specific ligand thermodynamically controls (either stabilizes or destabilizes) the intracellular region, consisting of crucial dynamic elements ICL2 and ICL3, and differential conformational transitions of these loops determine ligand-dependent functional selectivity.

Molecular encapsulation by eosin yellow-β-cyclodextrin conjugate: Differential binding to quadruplex and duplex DNA

G-Quadruplexes exist in cells and are involved in several cancer-related processes. Thus, G-quadruplex binding is an effective way of aiming at cancer cells. In this paper, the binding of a new β-cyclodextrin conjugate, eosin yellow-CD, with quadruplex DNAs viz., kit22, myc22, telo24 and duplex calf-thymus DNA, is studied by UV–vis and fluorescence spectroscopic methods. Berberine, a known G-quadruplex binder, is loaded in the modified cyclodextrin. The binding constant of the host: guest complex is 1.9 × 106 mol−1 dm3. The proximity of the protons of the conjugate and Berberine with those of the internal cavity of β-cyclodextrin in the conjugate is confirmed by two-dimensional rotating-frame Overhauser effect spectroscopy. The conjugate displays a quenching of fluorescence selectively to the quadruplex kit22 and telo24 that contrasts the spectral characteristics of the ligand-bound complexes of myc22 and ctDNA. The Stern-Volmer quenching constants are 1.08 × 106 mol−1 dm3 and 1.6 × 106 mol−1 dm3, respectively for binding to kit22 and telo24. The duplex ctDNA exhibits an enhanced fluorescence in the absence and the presence of Berberine complex, whereas the quadrulex DNAs exhibit dissimilar emission profiles when the eosin yellow-CD encapsulates Berberine. We report the differences in the binding constants and the mode of binding of free- and Berberine-loaded eosin yellow-CD to the duplex and quadruplex DNAs.

In-silico study of seaweed secondary metabolites as AXL kinase inhibitors

AXL kinase is an attractive cancer target for drug design and it is involved in different cancers. A set of molecule databases with 1072 secondary metabolites from seaweeds were screened against the AXL kinase active site and eight molecules were shortlisted for further studies. From the docking analysis of the complexes, four molecules GA011, BE005, BC010, and BC005 are showing prominent binging. From the 100 ns of molecular dynamics simulations and ligand-bound complex MM-PBSA free energy analysis, two molecules BC010 (ΔG = −135.38 kJ/mol) and BE005 (ΔG = −141.72 kJ/mol) are showing molecule stability in the active site also showing very strong binding free energies. It suggests these molecules could be the potent molecules for AXL kinase.

Mechanistic Insights into SARS-CoV-2 Main Protease Inhibition Reveals Hotspot Residues.

The main protease (Mpro) is a key enzyme responsible for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication that causes the spread of the global pandemic novel coronavirus (nCOVID-19) infection. In the present study, multiple computational approaches such as docking, long-range molecular dynamics (MD) simulations, and binding free-energy (BFE) estimation techniques were employed to investigate the mechanistic basis of the high-affinity inhibitors─GC-376, Calpain XII, and Calpain II (hereafter Calpain as Cal) from the literature─binding to Mpro. Redocking GC-376 and docking Cal XII and Cal II inhibitors to Mpro were able to reproduce all crucial interactions like the X-ray conformation. Subsequently, the apo (ligand-free) and three holo (ligand-bound) complexes were subjected to extensive MD simulations, which revealed that the ligand binding did not alter the overall Mpro structural features, whereas the heatmap analysis showed that the residues located in subsites S1 and S2, the catalytic dyad, and the 45TSEDMLN51 loop in Mpro exhibit a conformational deviation. Moreover, the BFE estimation method was used to elucidate the crucial thermodynamic properties, which revealed that Coulomb, solvation surface accessibility (Solv_SA), and lipophilic components contributed significant energies for complex formation. The decomposition of the total BFE to per-residue showed that H41, H163, M165, Q166, and Q189 residues contributed maximum energies. The overall results from the current investigation might be valuable for designing novel anti-Mpro inhibitors.

A structural intermediate pre-organizes the add adenine riboswitch for ligand recognition.

Riboswitches are RNA sequences that regulate gene expression by undergoing structural changes upon the specific binding of cellular metabolites. Crystal structures of purine-sensing riboswitches have revealed an intricate network of interactions surrounding the ligand in the bound complex. The mechanistic details about how the aptamer folding pathway is involved in the formation of the metabolite binding site have been previously shown to be highly important for the riboswitch regulatory activity. Here, a combination of single-molecule FRET and SHAPE assays have been used to characterize the folding pathway of the adenine riboswitch from Vibrio vulnificus. Experimental evidences suggest a folding process characterized by the presence of a structural intermediate involved in ligand recognition. This intermediate state acts as an open conformation to ensure ligand accessibility to the aptamer and folds into a structure nearly identical to the ligand-bound complex through a series of structural changes. This study demonstrates that the add riboswitch relies on the folding of a structural intermediate that pre-organizes the aptamer global structure and the ligand binding site to allow efficient metabolite sensing and riboswitch genetic regulation.

Molecular dynamics and in silico mutagenesis on the reversible inhibitor-bound SARS-CoV-2 main protease complexes reveal the role of lateral pocket in enhancing the ligand affinity

The 2019 novel coronavirus pandemic caused by SARS-CoV-2 remains a serious health threat to humans and there is an urgent need to develop therapeutics against this deadly virus. Recent scientific evidences have suggested that the main protease (Mpro) enzyme in SARS-CoV-2 can be an ideal drug target due to its crucial role in the viral replication and transcription processes. Therefore, there are ongoing research efforts to identify drug candidates against SARS-CoV-2 Mpro that resulted in hundreds of X-ray crystal structures of ligand-bound Mpro complexes in the Protein Data Bank (PDB) describing the interactions of different fragment chemotypes within different sites of the Mpro. In this work, we performed rigorous molecular dynamics (MD) simulation of 62 reversible ligand–Mpro complexes in the PDB to gain mechanistic insights about their interactions at the atomic level. Using a total of over 3 µs long MD trajectories, we characterized different pockets in the apo Mpro structure, and analyzed the dynamic interactions and binding affinity of ligands within those pockets. Our results identified the key residues that stabilize the ligands in the catalytic sites and other pockets of Mpro. Our analyses unraveled the role of a lateral pocket in the catalytic site in Mpro that is critical for enhancing the ligand binding to the enzyme. We also highlighted the important contribution from HIS163 in the lateral pocket towards ligand binding and affinity against Mpro through computational mutation analyses. Further, we revealed the effects of explicit water molecules and Mpro dimerization in the ligand association with the target. Thus, comprehensive molecular-level insights gained from this work can be useful to identify or design potent small molecule inhibitors against SARS-CoV-2 Mpro.

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.

Screening of phytochemicals as potent inhibitor of 3-chymotrypsin and papain-like proteases of SARS-CoV2: an in silico approach to combat COVID-19

COVID-19 and its causative organism SARS-CoV2 that emerged from Wuhan city, China have paralyzed the world. With no clinically approved drugs, the global health system is struggling to find an effective treatment measure. At this crucial juncture, screening of plant-derived compounds may be an effective strategy to combat COVID-19. The present study investigated the binding affinity of phytocompounds with 3-Chymotrypsin-like (3CLpro) and Papain-like proteases (PLpro) of SARS-CoV2 using in-silico techniques. A total of 32 anti-protease phytocompounds were investigated for the binding affinity to the proteins. Docking was performed in Autodock Vina. Pharmacophore descriptors of best ligands were studied using LigandScout. Molecular dynamics (MD) simulation of apo-protein and ligand-bound complexes was carried out in YASARA software. The druglikeness properties of phytocompounds were studied using ADMETlab. Out of 32 phytochemicals, amentoflavone and gallocatechin gallate showed the best binding affinity to 3CLpro (–9.4 kcal/mol) and PLpro (–8.8 kcal/mol). Phytochemicals such as savinin, theaflavin-3,3-digallate, and kazinol-A also showed strong affinity. MD simulation revealed ligand-induced conformational changes in the protein with decreased surface area and higher stability. The RMSD/F of proteins and ligands showed stability of the protein suggesting the effective binding of the ligand in both the proteins. Both amentoflavone and gallocatechin gallate possess promising druglikeness property. The present study thus suggests that Amentoflavone and Gallocatechin gallate may be potential inhibitors of 3CLpro and PLpro proteins and effective drug candidates for SARS-CoV2. However, the findings of in silico study need to be supported by in vivo studies to establish the exact mode of action. Communicated by Ramaswamy H. Sarma