Conformational Characterization of the Co-Activator Binding Site Revealed the Mechanism to Achieve the Bioactive State of FXR

In silico screening of dicoumarols as potential Mur B enzyme inhibitors in Mycobacterium tuberculosis: molecular docking, ADME, QSAR and MD simulations

UDP-N-acetylenolpyruvoylglucosamine reductase (Mur B) from Mycobacterium tuberculosis has gathered significant pharmaceutical interest as a pivotal target because of its essential role in bacterial viability. This study employed computational methods to screen and assess the inhibitory potential of dicoumarol derivatives against the Mur B protein. A diverse set of dicoumarols, sourced from PubChem and Zinc database, is subjected to molecular docking, ADME studies, and MD simulations to elucidate interacting modes and stability. A QSAR model was constructed for dicoumarol derivatives based on known inhibitor MIC values against Staphylococcus aureus and Mur B. The four best dicoumarols (CID142097979, CID54716867, CID91962283, CID54705236) aligned well with the model. Subsequently, these dicoumarols were scrutinized via 200 ns MD simulations and MM-PBSA analysis to assess their complex stability with Mur B protein. Various MD simulation parameters such as RMSD, RMSF, Rg, H-bonds, PCA, and FEL were employed. The 200 ns MD simulation analysis outcomes indicated that the Mur B-CID54705236 complex exhibited the highest stability and possessed the binding energy of -59.96 kcal/mol further verifying its stability. The post-dynamic simulation analysis showed four hydrogen bond formations with Ser70, Asn71, Leu72 and Gln137 residues at the active site of Mur B. Overall, these results underscored dicoumarol derivatives as potential Mur B inhibitors and these findings can serve as a basis for further in vitro studies against Mur B protein.

Investigate the binding of pesticides with the TLR4 receptor protein found in mammals and zebrafish using molecular docking and molecular dynamics simulations

The widespread use of pesticides poses significant threats to both environmental and human health, primarily due to their potential toxic effects. The study investigated the cardiovascular toxicity of selected pesticides, focusing on their interactions with Toll-like receptor 4 (TLR4), an important part of the innate immune system. Using computational tools such as molecular docking, molecular dynamics (MD) simulations, principal component analysis (PCA), density functional theory (DFT) calculations, and ADME analysis, this study identified C160 as having the lowest binding affinity (-8.2 kcal/mol), followed by C107 and C165 (-8.0 kcal/mol). RMSD, RMSF, Rg, and hydrogen bond metrics indicated the formation of stable complexes between specific pesticides and TLR4. PCA revealed significant structural changes upon ligand binding, affecting stability and flexibility, while DFT calculations provided information about the stability, reactivity, and polarity of the compounds. ADME studies highlighted the solubility, permeability, and metabolic stability of C107, C160, and C165, suggesting their potential for bioavailability and impact on cardiovascular toxicity. C107 and C165 exhibit higher bioactivity scores, indicating favourable absorption, metabolism, and distribution properties. C165 also violated rule where molecular weight is greater than 500 g/mol. Further, DFT and NCI analysis of post MD conformations confirmed the binding of ligands at the binding pocket. The analysis shed light on the molecular mechanisms of pesticide-induced cardiovascular toxicity, aiding in the development of strategies to mitigate their harmful effects on human health.

Molecular analysis of dUTPase of Helicobacter pylori for identification of novel inhibitors using in silico studies

The human gastric pathogen Helicobacter pylori chronically affects the gastric mucosal layer of approximately half of world's population. The emergence of resistant strains urges the need for identification of novel and selective drug against new molecular targets. A ubiquitous enzyme, Deoxyuridine 5'-triphosphate nucleotidohydrolase (dUTPase), is considered as first line of defense against uracil mis-incorporation into DNA, and essential for genome integrity. Lack of dUTPase triggers an elevated recombination frequency, DNA breaks and ultimately cell death. Hence, dUTPase can be considered as a promising target for development of novel lead inhibitor compounds in H. pylori treatment. Herein, we report the generation of three-dimensional model of the target protein using comparative modelling and its validation. To identify dUTPase inhibitors, a high throughput virtual screening approach utilizing Knowledge-based inhibitors and DrugBank database was implemented. Top ranked compounds were scrutinized based on investigations of the protein-ligand interaction fingerprints, molecular interaction maps and binding affinities and the drug potentiality. The best ligands were studied further for complex stability and intermolecular interaction profiling with respect to time under 100 ns classical molecular dynamic stimulation, establishing significant stability in dynamic states as observed from RMSD and RMSF parameters and interactions with the catalytic site residues. The binding free energy calculation computed using MM-GBSA method from the MD simulation trajectories demonstrated that our molecules possess strong binding affinity towards the Helicobacter pylori dUTPase protein. We conclude that our proposed molecules may be potential lead molecules for effective inhibition against the H. pylori dUTPase protein subject to experimental validation.Communicated by Ramaswamy H. Sarma.

7D, a small molecule inhibits dengue infection by increasing interferons and neutralizing-antibodies via CXCL4:CXCR3:p38:IRF3 and Sirt1:STAT3 axes respectively

There are a limited number of effective vaccines against dengue virus (DENV) and significant efforts are being made to develop potent anti-virals. Previously, we described that platelet-chemokine CXCL4 negatively regulates interferon (IFN)-α/β synthesis and promotes DENV2 replication. An antagonist to CXCR3 (CXCL4 receptor) reversed it and inhibited viral replication. In a concurrent search, we identified CXCR3-antagonist from our compound library, namely 7D, which inhibited all serotypes of DENV in vitro. With a half-life of ~2.85 h in plasma and no significant toxicity, 7D supplementation (8 mg/kg-body-weight) to DENV2-infected IFNα/β/γR-/-AG129 or wild-type C57BL6 mice increased synthesis of IFN-α/β and IFN-λ, and rescued disease symptoms like thrombocytopenia, leukopenia and vascular-leakage, with improved survival. 7D, having the property to inhibit Sirt-1 deacetylase, promoted acetylation and phosphorylation of STAT3, which in-turn increased plasmablast proliferation, germinal-center maturation and synthesis of neutralizing-antibodies against DENV2 in mice. A STAT3-inhibitor successfully inhibited these effects of 7D. Together, these observations identify compound 7D as a stimulator of IFN-α/β/λ synthesis via CXCL4:CXCR3:p38:IRF3 signaling, and a booster for neutralizing-antibody generation by promoting STAT3-acetylation in plasmablasts, capable of protecting dengue infection.

Illuminating ligand-induced dynamics in nuclear receptors through MD simulations

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

On the Cholesterol Raising Effect of Coffee Diterpenes Cafestol and 16-O-Methylcafestol: Interaction with Farnesoid X Receptor

The diterpene cafestol represents the most potent cholesterol-elevating compound known in the human diet, being responsible for more than 80% of the effect of coffee on serum lipids, with a mechanism still not fully clarified. In the present study, the interaction of cafestol and 16-O-methylcafestol with the stabilized ligand-binding domain (LBD) of the Farnesoid X Receptor was evaluated by fluorescence and circular dichroism. Fluorescence quenching was observed with both cafestol and 16-O-methylcafestol due to an interaction occurring in the close environment of the tryptophan W454 residue of the protein, as confirmed by docking and molecular dynamics. A conformational change of the protein was also observed by circular dichroism, particularly for cafestol. These results provide evidence at the molecular level of the interactions of FXR with the coffee diterpenes, confirming that cafestol can act as an agonist of FXR, causing an enhancement of the cholesterol level in blood serum.

Elevated glycosylation of CD36 in platelets is a risk factor for oxLDL-mediated platelet activation in type 2 diabetes

Platelet activation and related cardiovascular complications are the hallmarks of type 2 diabetes (T2D). We investigated the mechanism of platelet activation in T2D using MS-based identification of differentially expressed platelet proteins with a focus on glycosylated forms. Glycosylation is considered one of the common post-translational modifications in T2D, and N/O-linked glycosylation of glycoproteins (GPs)/integrins is known to play crucial roles in platelet activation. Our platelet proteome data revealed elevated levels of GPs GPIbα, GPIIbIIIa, GPIV (CD36), GPV and integrins in T2D patients. T2D platelets had elevated N-linked glycosylation of CD36 at asparagine (Asn)408,417 . Enrichment analysis revealed a close association of glycosylated CD36 with thrombospondin-1, fibrinogen and SERPINA1 in T2D platelets. The glycosylation of CD36 has previously been reported to increase cellular uptake of long-chain fatty acids. Our in silico molecular docking data also showed a favorable binding of cholesterol with glycosylated Asn417 CD36 compared to the non-glycosylated form. We further investigated the CD36:LDL cholesterol axis in T2D. Elevated levels of oxidized-low density lipoprotein (oxLDL) were found to cause significant platelet activation via CD36-mediated stimulation of Lyn-JNK signaling. Sulfo-N-succinimidyl oleate, an inhibitor of CD36, effectively inhibited oxLDL-mediated platelet activation and adhesion in vitro. Our study suggests increased glycosylation of CD36 in T2D platelets as a potential route for oxLDL-mediated platelet activation. The oxLDL:CD36 axis may thus be exploited as a prospective target to develop therapeutics against thrombosis in T2D.

Fangchinoline inhibits SARS-CoV-2 and MERS-CoV entry

The COVID-19 pandemic caused by SARS-CoV-2, lead to mild to severe respiratory illness and resulted in 6.9 million deaths worldwide. Although vaccines are effective in preventing COVID-19, they may not be sufficient to protect immunocompromised individuals from this respiratory illness. Moreover, novel emerging variants of SARS-CoV-2 pose a risk of new COVID-19 waves. Therefore, identification of effective antivirals is critical in controlling SARS and other coronaviruses, such as MERS-CoV. We show that Fangchinoline (Fcn), a bisbenzylisoquinoline alkaloid, inhibits replication of SARS-CoV, SARS-CoV-2, and MERS-CoV in a range of in vitro assays, by blocking entry. Therapeutic use of Fcn inhibited viral loads in the lungs, and suppressed associated airway inflammation in hACE2. Tg mice and Syrian hamster infected with SARS-CoV-2. Combination of Fcn with remdesivir (RDV) or an anti-leprosy drug, Clofazimine, exhibited synergistic antiviral activity. Compared to Fcn, its synthetic derivative, MK-04-003, more effectively inhibited SARS-CoV-2 and its variants B.1.617.2 and BA.5 in mice. Taken together these data demonstrate that Fcn is a pan beta coronavirus inhibitor, which possibly can be used to combat novel emerging coronavirus diseases.

Understanding the role of conserved proline and serine residues in the SARS-CoV-2 spike cleavage sites in the virus entry, fusion, and infectivity

The spike (S) glycoprotein of the SARS-CoV-2 virus binds to the host cell receptor and promotes the virus's entry into the target host cell. This interaction is primed by host cell proteases like furin and TMPRSS2, which act at the S1/S2 and S2´ cleavage sites, respectively. Both cleavage sites have serine or proline residues flanking either the single or polybasic region and were found to be conserved in coronaviruses. Unravelling the effects of these conserved residues on the virus entry and infectivity might facilitate the development of novel therapeutics. Here, we have investigated the role of the conserved serine and proline residues in the SARS-CoV-2 spike mediated entry, fusogenicity, and viral infectivity by using the HIV-1/spike-based pseudovirus system. A conserved serine residue mutation to alanine (S2´S-A) at the S2´ cleavage site resulted in the complete loss of spike cleavage. Exogenous treatment with trypsin or overexpression of TMPRSS2 protease could not rescue the loss of spike cleavage and biological activity. The S2´S-A mutant showed no significant responses against E-64d, TMPRSS2 or other relevant inhibitors. Taken together, serine at the S2´ site in the spike protein was indispensable for spike protein cleavage and virus infectivity. Thus, novel interventions targeting the conserved serine at the S2´ cleavage site should be explored to reduce severe disease caused by SARS-CoV-2-and novel emerging variants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03749-y.

Discovery of non-nucleoside oxindole derivatives as potent inhibitors against dengue RNA-dependent RNA polymerase

A series of thiazole linked Oxindole-5-Sulfonamide (OSA) derivatives were designed as inhibitors of RNA-dependent RNA polymerase (RdRp) activity of Dengue virus. These were synthesized and then evaluated for their efficacy in ex-vivo virus replication assay using human cell lines. Among 20 primary compounds in the series, OSA-15 was identified as a hit. A series of analogues were synthesized by replacing the difluoro benzyl group of OSA-15 with different substituted benzyl groups. The efficacy of OSA-15derivatives was less than that of the parent compound, except OSA-15-17, which has shown improved efficacy than OSA-15. The further optimization was carried out by adding dimethyl (DM) groups to both the sulfonamide and oxindole NH's to produce OSA-15-DM and OSA-15-17-DM. These two compounds were showing no detectable cytotoxicity and the latter was more efficacious. Further, both these compounds were tested for inhibition in all the serotypes of the Dengue virus using an ex-vivo assay. The EC₅₀ of OSA-15-17-DM was observed in a low micromolar range between 2.5 and 5.0 µg/ml. Computation docking and molecular dynamics simulation studies confirmed the binding of identified hits to DENV RdRp. OSA15-17-DM blocks the RNA entrance and elongation site for their biological activity with high binding affinity. Overall, the identified oxindole derivatives are novel compounds that can inhibit Dengue replication, working as non-nucleoside inhibitors (NNI) to explore as anti-viral RdRp activity.

Deciphering the Structural Determinants Critical in Attaining the FXR Partial Agonism

Elucidation of structural determinants is pivotal for structure-based drug discovery. The Farnesoid X receptor (FXR) is a proven target for NASH; however, its full agonism causes certain clinical complications. Therefore, partial agonism (PA) appears as a viable alternative for improved therapeutics. Since the agonist and PA both share the same binding site, i.e., ligand-binding pocket (LBP), which is highly dynamic and has synergy with the substrate binding site, the selective designing of PA is challenging. The identification of structural and conformational determinants is critical for PA compared with an agonist. Furthermore, the mechanism by which PA modulates the structural dynamics of FXR at the residue level, a prerequisite for PA designing, is still elusive. Here, by using ∼4.5 μs of MD simulations and residue-wise communication network analysis, we identified the structural regions which are flexible with PA but frozen with an agonist. Also, the network analysis identified the considerable changes between an agonist and PA in biologically essential zones of FXR such as helix H10/H11 and loop L:H11/H12, which lead to the modulation of synergy between LBP and the substrate binding site. Furthermore, the thermodynamic profiling suggested the methionine residues, mainly M³²⁸, M³⁶⁵, and M⁴⁵⁰, seem to be responsible for the recruitment of PA. The other residues I³⁵⁷, Y³⁶¹, L⁴⁶⁵, F³⁰⁸, Q³¹⁶, and K³²¹ are also identified, exclusively interacting with PA. This study offers novel structural and mechanistic insights that are critical for FXR targeted drug discovery for PA designing.

Characterizing (un)binding mechanism of USP7 inhibitors to unravel the cause of enhanced binding potencies at allosteric checkpoint

The ability to predict the intricate mechanistic behavior of ligands and associated structural determinants during protein-ligand (un)binding is of great practical importance in drug discovery. Ubiquitin specific protease-7 (USP7) is a newly emerging attractive cancer therapeutic target with bound allosteric inhibitors. However, none of the inhibitors have reached clinical trials, allowing opportunities to examine every aspect of allosteric modulation. The crystallographic insights reveal that these inhibitors have common properties such as chemical scaffolds, binding site and interaction fingerprinting. However, they still possess a broader range of binding potencies, ranging from 22 nM to 1,300 nM. Hence, it becomes more critical to decipher the structural determinants guiding the enhanced binding potency of the inhibitors. In this regard, we elucidated the atomic-level insights from both interacting partners, that is, protein-ligand perspective, and established the structure-activity link between USP7 inhibitors by using classical and advanced molecular dynamics simulations combined with linear interaction energy and molecular mechanics-Poisson Boltzmann surface area. We revealed the inhibitor potency differences by examining the contributions of chemical moieties and USP7 residues, the involvement of water-mediated interactions, and the thermodynamic landscape alterations. Additionally, the dissociation profiles aided in the establishment of a correlation between experimental potencies and structural determinants. Our study demonstrates the critical role of blocking loop 1 in allosteric inhibition and enhanced binding affinity. Comprehensively, our findings provide a constructive expansion of experimental outcomes and show the basis for varying binding potency using in-silico approaches. We expect this atomistic approach to be useful for effective drug design.