Rational Computational Workflow for Structure-Guided Discovery of a Novel USP7 Inhibitor

Rationally applied, structurally guided computational methods hold the promise of identifying potent and distinct chemotypes while enabling the precise targeting of structural determinants. Here, we implemented a computational workflow combining insights from cocrystal poses and monitoring the dynamical structural determinants from our previous studies for the identification of potential candidates against USP7. We identified and tested several diverse chemical scaffolds, which underwent in vitro validation across six cancer cell lines. Among these hits, compound M15, belonging to the benzothiazole chemical class, exhibited remarkable anticancer activities, demonstrating dose-dependent reduction in cancer cell viability across all cell lines and indicating that it is a promising candidate to explore as a potent anticancer drug. Biophysical binding confirms binding of M15 on USP7. M15 also exhibited certain USP7 inhibitory activity, as observed in the enzymatic assay. A comparative cocrystal mining of reported USP7 inhibitors unveiled a distinct binding mode of M15, which nicely cross-corroborated with MD and binding-pose metadynamics simulations. Notably, M15 occupies both the determinants, i.e., BL1 and the allosteric checkpoint, which has not yet been underscored as a druggable site. In essence, our study presents a robust and multifaceted computational method for the discovery and characterization of a novel inhibitor scaffold, exemplified by the identification and mechanistic elucidation of M15 against USP7. This integrated approach not only advances our understanding of USP7 inhibition and underscores mechanistic determinants but also offers promising avenues for the discovery of target-specific therapeutic intervention.

Molecular Dynamics Simulations Reveal the Interaction Fingerprint of Remdesivir Triphosphate Pivotal in Allosteric Regulation of SARS-CoV-2 RdRp

The COVID-19 pandemic has now strengthened its hold on human health and coronavirus' lethal existence does not seem to be going away soon. In this regard, the optimization of reported information for understanding the mechanistic insights that facilitate the discovery towards new therapeutics is an unmet need. Remdesivir (RDV) is established to inhibit RNA-dependent RNA polymerase (RdRp) in distinct viral families including Ebola and SARS-CoV-2. Therefore, its derivatives have the potential to become a broad-spectrum antiviral agent effective against many other RNA viruses. In this study, we performed comparative analysis of RDV, RMP (RDV monophosphate), and RTP (RDV triphosphate) to undermine the inhibition mechanism caused by RTP as it is a metabolically active form of RDV. The MD results indicated that RTP rearranges itself from its initial RMP-pose at the catalytic site towards NTP entry site, however, RMP stays at the catalytic site. The thermodynamic profiling and free-energy analysis revealed that a stable pose of RTP at NTP entrance site seems critical to modulate the inhibition as its binding strength improved more than its initial RMP-pose obtained from docking at the catalytic site. We found that RTP not only occupies the residues K545, R553, and R555, essential to escorting NTP towards the catalytic site, but also interacts with other residues D618, P620, K621, R624, K798, and R836 that contribute significantly to its stability. From the interaction fingerprinting it is revealed that the RTP interact with basic and conserved residues that are detrimental for the RdRp activity, therefore it possibly perturbed the catalytic site and blocked the NTP entrance site considerably. Overall, we are highlighting the RTP binding pose and key residues that render the SARS-CoV-2 RdRp inactive, paving crucial insights towards the discovery of potent inhibitors.

Binding mode characterization of 13b in the monomeric and dimeric states of SARS-CoV-2 main protease using molecular dynamics simulations

The main protease, M^pro/3CL^pro, plays an essential role in processing polyproteins translated from viral RNA to produce functional viral proteins and therefore serve as an attractive target for discovering COVID-19 therapeutics. The availability of both monomer and dimer crystal bound with a common ligand, '13b' (α-ketoamide inhibitor), opened up opportunities to understand the M^pro mechanism of action. A comparative analysis of both forms of M^pro was carried out to elucidate the binding site architectural differences in the presence and absence of '13b'. Molecular dynamics simulations suggest that the presence of '13b' enhances the stability of M^pro than the unbound APO form. The N- and C- terminals of both the protomers stabilize each other, and making it's interface essential for the active form of M^pro. In comparison to monomer, the relatively high affinity of '13b' is gained in dimer pocket due to the high stability of the pocket by the interaction of S1 residue of chain B with residues F140, E166 and H172 of chain A, which is absent in monomer. The comprehensive essential dynamics, protein structure network analysis and thermodynamic profiling highlight the hot-spots, pivotal in molecular recognition process at protein-ligand and protein-protein interaction levels, cross-validated through computational alanine scanning study. A comparative description of '13b' binding mechanism in both forms illustrates valuable insights into the inhibition mechanism and the selection of critical residues suitable for the structure-based approaches for the identification of more potent M^pro inhibitors.Communicated by Ramaswamy H. Sarma.

Comparative binding analysis of noscapine and piperine with tRNA: A structural perturbation and energetic study

In this study, we have exploring the binding mechanisms of the two anticancer alkaloid noscapine (NOS) and piperine (PIP) with tRNA using different spectroscopy and computational method. Absorbance and emission spectra revealed that both the drugs show strong binding with tRNA, where NOS intercalate between the base pairs of tRNA and PIP binds in the groove of tRNA. Competitive binding study and steady state anisotropy further confirms the intercalative mode of binding between NOS and tRNA and groove binding in PIP-tRNA complex. The observed thermodynamic parameters suggested that NOS-tRNA complex formation is endothermic and entropy driven, however it was exothermic, and enthalpy driven in case of PIP-tRNA complex. CD and time resolved fluorescence studies show the structural perturbations and conformational change in tRNA structure with NOS as well as PIP. Molecular docking studies are comparable with experimental results and further confirmed that the hydrophobic interactions involved in the NOS-tRNA binding, whereas hydrogen binding and van der Waals interactions play important role in the PIP-tRNA complex formation. This study can be useful to understand the potential binding and resultant tRNA damage by alkaloids and deigned new target specific anticancer drug.

Rational design of 9-vinyl-phenyl noscapine as potent tubulin binding anticancer agent and evaluation of the effects of its combination on Docetaxel

Docetaxel (DOX) based combination therapy is a novel therapeutic strategy that attracts great interest in breast cancer treatment but its clinical utility got limited due to side effects. In contrast, noscapine, an antitussive drug showed antitumor activity against many cancers without any side effects that targets microtubules and attenuates its dynamic instability. In the quest for an increase in the anticancer activity of noscapine, we strategically designed a novel derivative, 9-vinyl phenyl noscapine (VPN), based on our in silico molecular docking and molecular dynamics simulation effort. Molecular docking of VPN and DOX onto microtubule revealed a docking score of -4.82 kcal/mol and -6.67 kcal/mol respectively, while the docking score of VPN was changed to -3.23 kcal/mol when it was docked onto the co-complex of tubulin-DOX. Further, the binding free energy (ΔG_bind,PBSA) of VPN and DOX with tubulin showed -24.04 and -18.65 kcal/mol respectively, while the binding free energy of DOX was increased further in combination with VPN (ΔG_{bind, PBSA} was reduced to -21.41 kcal/mol), denoting combination effect of both ligands. The IC₅₀ value amounted to 30.17 µM and 19.92 µM for VPN and 0.621 µM and 0.193 µM for DOX, respectively for 48 h and 72 h. The dose dependent cytotoxicity of DOX has been reduced considerably with the combination dose regimen of VPN. Further, the combine effect of both the agents improved the apoptotic cell death 28.5% compared to single agent treatment 5.71% and 10.5% for VPN and DOX, respectively. Both agents bind effectively to tubulin in single and in combination to interfere with cell cycle progression in G2/M transition. This study provides novel concept of combination treatment of DOX and VPN to amend efficiency in breast cancer treatment.Communicated by Ramaswamy H. Sarma.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Functional Analysis of γ-Tubulin Complex Proteins Indicates Specific Lateral Association via Their N-terminal Domains

Microtubules are nucleated from multiprotein complexes containing γ-tubulin and associated γ-tubulin complex proteins (GCPs). Small complexes (γTuSCs) comprise two molecules of γ-tubulin bound to the C-terminal domains of GCP2 and GCP3. γTuSCs associate laterally into helical structures, providing a structural template for microtubule nucleation. In most eukaryotes γTuSCs associate with additional GCPs (4, 5, and 6) to form the core of the so-called γ-tubulin ring complex (γTuRC). GCPs 2-6 constitute a family of homologous proteins. Previous structural analysis and modeling of GCPs suggest that all family members can potentially integrate into the helical structure. Here we provide experimental evidence for this model. Using chimeric proteins in which the N- and C-terminal domains of different GCPs are swapped, we show that the N-terminal domains define the functional identity of GCPs, whereas the C-terminal domains are exchangeable. FLIM-FRET experiments indicate that GCP4 and GCP5 associate laterally within the complex, and their interaction is mediated by their N-terminal domains as previously shown for γTuSCs. Our results suggest that all GCPs are incorporated into the helix via lateral interactions between their N-terminal domains, whereas the C-terminal domains mediate longitudinal interactions with γ-tubulin. Moreover, we show that binding to γ-tubulin is not essential for integrating into the helical complex.

Identification of the active site of human mitochondrial malonyl-coenzyme a decarboxylase: A combined computational study

Malonyl-CoA decarboxylase (MCD) can control the level of malonyl-CoA in cell through the decarboxylation of malonyl-CoA to acetyl-CoA, and plays an essential role in regulating fatty acid metabolism, thus it is a potential target for drug discovery. However, the interactions of MCD with CoA derivatives are not well understood owing to unavailable crystal structure with a complete occupancy in the active site. To identify the active site of MCD, molecular docking and molecular dynamics simulations were performed to explore the interactions of human mitochondrial MCD (HmMCD) and CoA derivatives. The findings reveal that the active site of HmMCD indeed resides in the prominent groove which resembles that of CurA. However, the binding modes are slightly different from the one observed in CurA due to the occupancy of the side chain of Lys183 from the N-terminal helical domain instead of the adenine ring of CoA. The residues 300 - 305 play an essential role in maintaining the stability of complex mainly through hydrogen bond interactions with the pyrophosphate moiety of acetyl-CoA. Principle component analysis elucidates the conformational distribution and dominant concerted motions of HmMCD. MM_PBSA calculations present the crucial residues and the major driving force responsible for the binding of acetyl-CoA. These results provide useful information for understanding the interactions of HmMCD with CoA derivatives. Proteins 2016; 84:792-802. © 2016 Wiley Periodicals, Inc.

In silico and experimental validation of protein-protein interactions between PknI and Rv2159c from Mycobacterium tuberculosis

Protein-protein interactions control the diverse and essential molecular processes inside the cell. To maintain the cellular physiology, protein kinases not only signal their substrates through reversible phosphorylation, but they also physically interact with them. PknI, a serine/threonine protein kinase of Mycobacterium tuberculosis is known to be important for cellular homoeostasis. In this study, we have identified the interacting proteins for PknI. We screened for proteins interacting with PknI using an in vitro assay, Far-western blot. This protein kinase specifically interacts with two peroxidase proteins of M. tuberculosis, Rv2159c and Rv0148. The PknI-Rv2159c interaction pair was further studied for the critical amino acid residues in Rv2159c that are responsible for the interaction. Rv2159c, a hypothetical protein is predicted to be an antioxidant with peroxidase activity. We performed homology modelling of Rv2159c protein and molecular docking using multiple docking servers such as Z-Dock and ClusPro. Further, the most favorable conformation of PknI-Rv2159c interaction was obtained using molecular dynamics simulation. The critical amino acid residues of the Rv2159c involved in interaction with PknI were identified. Mutation and docking analysis showed that the Ala1-Gly2-Trp3 residues in Rv2159c structure are responsible for the interaction. The free binding energy between the wild type and mutant complexes using MM-GBSA has provided insight about the stability of PknI-Rv2159c interaction. We propose that, PknI physically interacts with Rv2159c both in vitro and in silico studies.

Emerging microtubule targets in glioma therapy

Major advances in the genomics and epigenomics of diffuse gliomas and glioblastoma to date have not been translated into effective therapy, necessitating pursuit of alternative treatment approaches for these therapeutically challenging tumors. Current knowledge of microtubules in cancer and the development of new microtubule-based treatment strategies for high-grade gliomas are the topic in this review article. Discussed are cellular, molecular, and pharmacologic aspects of the microtubule cytoskeleton underlying mitosis and interactions with other cellular partners involved in cell cycle progression, directional cell migration, and tumor invasion. Special focus is placed on (1) the aberrant overexpression of βIII-tubulin, a survival factor associated with hypoxic tumor microenvironment and dynamic instability of microtubules; (2) the ectopic overexpression of γ-tubulin, which in addition to its conventional role as a microtubule-nucleating protein has recently emerged as a transcription factor interacting with oncogenes and kinases; (3) the microtubule-severing ATPase spastin and its emerging role in cell motility of glioblastoma cells; and (4) the modulating role of posttranslational modifications of tubulin in the context of interaction of microtubules with motor proteins. Specific antineoplastic strategies discussed include downregulation of targeted molecules aimed at achieving a sensitization effect on currently used mainstay therapies. The potential role of new classes of tubulin-binding agents and ATPase inhibitors is also examined. Understanding the cellular and molecular mechanisms underpinning the distinct behaviors of microtubules in glioma tumorigenesis and drug resistance is key to the discovery of novel molecular targets that will fundamentally change the prognostic outlook of patients with diffuse high-grade gliomas.