Hotspot site microenvironment in the deubiquitinase OTUB1 drives its stability and aggregation.

Lewy bodies (LB) are aberrant protein accumulations observed in the brain cells of individuals affected by Parkinson's disease (PD). A comprehensive analysis of LB proteome identified over a hundred proteins, many co-enriched with α-synuclein, a major constituent of LB. Within this context, OTUB1, a deubiquitinase detected in LB, exhibits amyloidogenic properties, yet the mechanisms underlying its aggregation remain elusive. In this study, we identify two critical sites in OTUB1-namely, positions 133 and 173-that significantly impact its amyloid aggregation. Substituting alanine at position 133 and lysine at position 173 enhances both thermodynamic and kinetic stability, effectively preventing amyloid aggregation. Remarkably, lysine at position 173 demonstrates the highest stability without compromising enzymatic activity. The increased stability and inhibition of amyloid aggregation are attributed mainly to the changes in the specific microenvironment at the hotspot. In our exploration of the in-vivo co-occurrence of α-synuclein and OTUB1 in LB, we observed a synergistic modulation of each other's aggregation. Collectively, our study unveils the molecular determinants influencing OTUB1 aggregation, shedding light on the role of specific residues in modulating aggregation kinetics and structural transition. These findings contribute valuable insights into the complex interplay of amino acid properties and protein aggregation, with potential implications for understanding broader aspects of protein folding and aggregation phenomena.

Polyphosphate kinase-1 regulates bacterial and host metabolic pathways involved in pathogenesis of Mycobacterium tuberculosis.

Inorganic polyphosphate (polyP) is primarily synthesized by Polyphosphate Kinase-1 (PPK-1) and regulates numerous cellular processes, including energy metabolism, stress adaptation, drug tolerance, and microbial pathogenesis. Here, we report that polyP interacts with acyl CoA carboxylases, enzymes involved in lipid biosynthesis in Mycobacterium tuberculosis. We show that deletion of ppk-1 in M. tuberculosis results in transcriptional and metabolic reprogramming. In comparison to the parental strain, the Δppk-1 mutant strain had reduced levels of virulence-associated lipids such as PDIMs and TDM. We also observed that polyP deficiency in M. tuberculosis is associated with enhanced phagosome-lysosome fusion in infected macrophages and attenuated growth in mice. Host RNA-seq analysis revealed decreased levels of transcripts encoding for proteins involved in either type I interferon signaling or formation of foamy macrophages in the lungs of Δppk-1 mutant-infected mice relative to parental strain-infected animals. Using target-based screening and molecular docking, we have identified raloxifene hydrochloride as a broad-spectrum PPK-1 inhibitor. We show that raloxifene hydrochloride significantly enhanced the activity of isoniazid, bedaquiline, and pretomanid against M. tuberculosis in macrophages. Additionally, raloxifene inhibited the growth of M. tuberculosis in mice. This is an in-depth study that provides mechanistic insights into the regulation of mycobacterial pathogenesis by polyP deficiency.

SARS-CoV-2 envelope protein attain Kac mediated dynamical interaction network to adopt 'histone mimic' at BRD4 interface.

Interface mimicry, achieved by recognition of host-pathogen interactions, is the basis by which pathogen proteins can hijack the host machinery. The envelope (E) protein of SARS-CoV-2 is reported to mimic the histones at the BRD4 surface via establishing the structural mimicry; however, the underlying mechanism of E protein mimicking the histones is still elusive. To explore the mimics at dynamic and structural residual network level an extensive docking, and MD simulations were carried out in a comparative manner between complexes of H3-, H4-, E-, and apo-BRD4. We identified that E peptide is able to attain an 'interaction network mimicry', as its acetylated lysine (Kac) achieves orientation and residual fingerprint similar to histones, including water-mediated interactions for both the Kac positions. We identified Y59 of E, playing an anchor role to escort lysine positioning inside the binding site. Furthermore, the binding site analysis confirms that E peptide needs a higher volume, similar to the H4-BRD4 where both the lysine's (Kac5 and Kac8) can accommodate nicely, however, the position of Kac8 is mimicked by two additional water molecules other than four water-mediated bridging's, strengthening the possibility that E peptide could hijack host BRD4 surface. These molecular insights seem pivotal for mechanistic understanding and BRD4-specific therapeutic intervention. KEY POINTSMolecular mimicry is reported in hijacking and then outcompeting the host counterparts so that pathogens can rewire their cellular function by overcoming the host defense mechanism.The molecular recognition process is the basis of molecular mimicry. The E peptide of SARS-CoV-2 is reported to mimic host histone at the BRD4 surface by utilizing its C-terminally placed acetylated lysine (Kac63) to mimic the N-terminally placed acetylated lysine Kac5GGKac8 histone (H4) by interaction network mimicry identified through microsecond molecular dynamics (MD) simulations and post-processing extensive analysis.There are two steps to mimic: firstly, tyrosine residues help E to anchor at the BRD4 surface to position Kac and increase the volume of the pocket. Secondary, after positioning of Kac, a common durable interaction network N140:Kac5; Kac5:W1; W1:Y97; W1:W2; W2:W3; W3:W4; W4:P82 is established between Kac5, with key residues P82, Y97, N140, and four water molecules through water mediate bridge. Furthermore, the second acetylated lysine Kac8 position and its interaction as polar contact with Kac5 were also mimicked by E peptide through interaction network P82:W5; W5:Kac63; W5:W6; W6:Kac63.The binding event at BRD4/BD1 seems an induced-fit mechanism as a bigger binding site volume was identified at H4-BRD4 on which E peptide attains its better stability than H3-BRD4.We identified the tyrosine residue Y59 of E that acts like an anchor on the BRD4 surface to position Kac inside the pocket and attain the interaction network by using aromatic residues of the BRD4 surface.Communicated by Ramaswamy H. Sarma.

Computational study of novel inhibitory molecule, 1-(4-((2S,3S)-3-amino-2-hydroxy-4-phenylbutyl)piperazin-1-yl)-3-phenylurea, with high potential to competitively block ATP binding to the RNA dependent RNA polymerase of SARS-CoV-2 virus.

For coronaviruses, RNA-dependent RNA polymerase (RdRp) is an essential enzyme that catalyses the replication from RNA template and therefore remains an attractive therapeutic target for anti-COVID drug discovery. In the present study, we performed a comprehensive in silico screening for 16,776 potential molecules from recently established drug libraries based on two important pharmacophores (3-amino-4-phenylbutan-2-ol and piperazine). Based on initial assessment, 4042 molecules were obtained suitable as drug candidates, which were following Lipinski's rule. Molecular docking implemented for the analysis of molecular interactions narrowed this number of compounds down to 19. Subsequent to screening filtering criteria and considering the critical parameters viz. docking score and MM-GBSA binding free energy, 1-(4-((2S,3S)-3-amino-2-hydroxy-4-phenylbutyl)piperazin-1-yl)-3-phenylurea (compound 1) was accomplished to score highest in comparison to the remaining 18 shortlisted drug candidates. Notably, compound 1 displayed higher docking score (-8.069 kcal/mol) and MM-GBSA binding free energy (-49.56 kcal/mol) than the control drug, remdesivir triphosphate, the active form of remdesivir as well as adenosine triphosphate. Furthermore, a molecular dynamics simulation was carried out (100 ns), which substantiated the candidacy of compound 1 as better inhibitor. Overall, our systematic in silico study predicts the potential of compound 1 to exhibit a more favourable specific activity than remdesivir triphosphate. Hence, we suggest compound 1 as a novel potential drug candidate, which should be considered for further exploration and validation of its potential against SARS-CoV-2 in wet lab experimental studies.Communicated by Ramasawamy H. Sarma.

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.

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

The main protease, Mpro/3CLpro, 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 Mpro mechanism of action. A comparative analysis of both forms of Mpro 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 Mpro 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 Mpro. 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 Mpro inhibitors.Communicated by Ramaswamy H. Sarma.

Identification of diphenyl furan derivatives via high throughput and computational studies as ArgA inhibitors of Mycobacterium tuberculosis.

Microbial amino acid biosynthetic pathways are underexploited for the development of anti-bacterial agents. N-acetyl glutamate synthase (ArgA) catalyses the first committed step in L-arginine biosynthesis and is essential for M. tuberculosis growth. Here, we have purified and optimized assay conditions for the acetylation of l-glutamine by ArgA. Using the optimized conditions, high throughput screening was performed to identify ArgA inhibitors. We identified 2,5-Bis (2-chloro-4-guanidinophenyl) furan, a dicationic diaryl furan derivatives, as ArgA inhibitor, with a MIC99 values of 1.56 µM against M. tuberculosis. The diaryl furan derivative displayed bactericidal killing against both M. bovis BCG and M. tuberculosis. Inhibition of ArgA by the lead compound resulted in transcriptional reprogramming and accumulation of reactive oxygen species. The lead compound and its derivatives showed micromolar binding with ArgA as observed in surface plasmon resonance and tryptophan quenching experiments. Computational and dynamic analysis revealed that these scaffolds share similar binding site residues with L-arginine, however, with slight variations in their interaction pattern. Partial restoration of growth upon supplementation of liquid cultures with either L-arginine or N-acetyl cysteine suggests a multi-target killing mechanism for the lead compound. Taken together, we have identified small molecule inhibitors against ArgA enzyme from M. tuberculosis.

Molecular dynamics simulation reveals the possible druggable hot-spots of USP7.

The plasticity in Ubiquitin Specific Proteases (USP7) inducing conformational changes at important areas has highlighted an intricate mechanism, by which USP7 is regulated. Given the importance of USP7 in oncogenic pathways and immune-oncology, identification of USP7 inhibitors has attracted considerable interest. Despite substantial efforts, the discovery of deubiquitinases (DUBs) inhibitors, knowledge of their binding site and understanding the possible mechanism of action has proven particularly challenging. We disclose the most likely binding site of P5091 (a potent USP7 inhibitor), which reveal a cryptic allosteric site through extensive computational studies in an inhibitor dependent and independent manner. Overall, these findings demonstrate the tractability and druggability of USP7. Through a series of molecular dynamics simulations and detailed quantitative analysis, a dynamically stable allosteric binding site near catalytic center of the inactive state of USP7 (site partially absent in active state), along with two newly identified sites have been revealed, which opens the avenue for rational structure-guided inhibitor designing in USP7 specific-manner.