Dynamics of allosteric modulation of lymphocyte function associated antigen-1 closure-open switch: unveiling the structural mechanisms associated with outside-in signaling activation

Battling BTK mutants with noncovalent inhibitors that overcome Cys481 and Thr474 mutations in Waldenström macroglobulinemia therapy: structural mechanistic insights on the role of fenebrutinib

Recently, the non-covalent Bruton tyrosine kinase (BTK) inhibitor fenebrutinib was presented as a therapeutic option with strong inhibitory efficacy against a single (C481S) and double (T474S/C481S) BTK variant in the treatment of Waldenström macroglobulinemia (WM). However, the molecular events surrounding its inhibition mechanism towards this variant remain unresolved. Herein, we employed in silico methods such as molecular dynamic simulation coupled with binding free energy estimations to explore the mechanistic activity of the fenebrutinib on (C481S) and (T474S/C481S) BTK variant, at a molecular level. Our investigations reveal that amino acid arginine contributed immensely to the total binding energy, this establishing the cruciality of amino acid residues, Arg132 and Arg156 in (C481S) and Arg99, Arg137, and Arg132 in (T474S/C481S) in the binding of fenebrutinib towards both BTK variants. The structural orientations of fenebrutinib within the respective hydrophobic pockets allowed favorable interactions with binding site residues, accounting for its superior binding affinity by 24.5% and relative high hydrogen bond formation towards (T474S/C481S) when compared with (C481S) BTK variants. Structurally, fenebrutinib impacted the stability, flexibility, and solvent accessible surface area of both BTK variants, characterized by various alterations observed in the bound and unbound structures, which proved enough to disrupt their biological function. Findings from this study, therefore, provide insights into the inhibitory mechanism of fenebrutinib at the atomistic level and reveal its high selectivity towards BTK variants. These insights could be key in designing and developing BTK mutants' inhibitors to treat Waldenström macroglobulinemia (WM).

Probing Alterations in MDM2 Catalytic Core Structure Effect of Garcinia Mangostana Derivatives: Insight from Molecular Dynamics Simulations

The MDM2-p53 protein-protein interaction is a promising model for researchers to design, study, and discover new anticancer drugs. The design of therapeutically active compounds that can maintain or restore the binding of MDM2 to p53 has been found to limit the oncogenic activities of both. This led to the current development of a group of xanthone-core and cis-imidazoline analogs compounds, among which γ-Mangostin (GM), α-Mangostin (AM), and Nutlin exhibited their MDM2-p53 interaction inhibitory effects. Therefore, in this study, we seek to determine the mechanisms by which these compounds elicit MDM2-p53 interaction targeting. Unique to the binding of GM, AM, and Nutlin, from our findings, they share the same three active site residues Val76, Tyr50, and Gly41, which represent the top active side residues that contribute to high electrostatic energy. Consequently, the free binding energy contributed enormously to the binding of these compounds, which culminated in the high binding affinities of GM, AM, and Nutlin with high values. Furthermore, GM, AM, and Nutlin commonly interrupted the stable and compact conformation of MDM2 coupled with its active site, where Cα deviations were relatively high. We believe that our findings would assist in the design of more potent active anticancer drugs.

Co-Binding of JQ1 and Venetoclax Exhibited Synergetic Inhibitory Effect for Cancer Therapy; Potential Line of Treatment for the Waldenström Macroglobulinemia Lymphoma

In recent times, the development of combination therapy has been a focal point in drug discovery. This article explores the potential synergistic effect of co-administration of Bcl2 inhibitor Venetoclax and BET inhibitor JQ1. We envisioned that the 'dual-site'-binding of Bcl2 has significant prospects and paves the way for the next round of rational design of potent Waldenström macroglobulinemia (WM) therapy. The preferential binding mechanisms of the multi-catalytic sites of the Bcl2 enzyme have been a subject of debate in the literature. This study conducted a systematic procedure to explore the preferred binding modes and the structural effects of co-binding at each catalytic active site. Interestingly, a mutual enhanced binding effect was observed - Venetoclax increased the binding affinity of JQ1 by 11.5 %, while JQ1 boosted the binding affinity of Venetoclax by 16.3 % when compared with individual inhibition of each drug. This synergistic binding effect has significantly increased protein stability, with substantial correlated movements and multiple van der Waals interactions. The structural and thermodynamic insights unveiled in this report would assist the future design of improved combined therapy against WM.

Multi-catalytic Sites Inhibition of Bcl2 Induces Expanding of Hydrophobic Groove: A New Avenue Towards Waldenström Macroglobulinemia Therapy

B-cell lymphoma 2 (Bcl2) is a key protein regulator of apoptosis. The hydrophobic groove in Bcl2 is a unique structural feature to this class of enzymes and found to have a profound impact on protein overall structure, function, and dynamics. Dynamics of the hydrophobic groove is an essential determinant of the catalytic activity of Bcl2, an implicated protein in Waldenström macroglobulinemia (WM). The mobility of α3-α4 helices around the catalytic site of the protein remains crucial to its activity. The preferential binding mechanisms of the multi-catalytic sites of the Bcl2 enzyme have been a subject of debate in the literature. In addition to our previous report on the same protein, herein, we further investigate the preferential binding modes and the conformational implications of Venetoclax-JQ1 dual drug binding at both catalytic active sites of Bcl2. Structural analysis revealed asymmetric α3-α4 helices movement with the expansion of the distance between the α3 and α4 helix in Venetoclax-JQ1 dual inhibition by 15.2% and 26.3%, respectively when compared to JQ1 and Venetoclax individual drug inhibition-resulting in remarkable widening of hydrophobic groove. Moreso, a reciprocal enhanced binding effect was observed: Venetoclax increased the binding affinity of JQ1 by 11.5%, while the JQ1 fostered the binding affinity of Venetoclax by 16.3% compared with individual inhibition of each drug. This divergence has also resulted in higher protein stability, and prominent correlated motions were observed with the least fluctuations and multiple van der Waals interactions. Findings offer vital conformational dynamics and structural mechanisms of enzyme-single ligand and enzyme-dual ligand interactions, which could potentially shift the current therapeutic protocol of Waldenström macroglobulinemia.

Inhibitory potential of repurposed drugs against the SARS-CoV-2 main protease: a computational-aided approach

The exponential increase in cases and mortality of coronavirus disease (COVID-19) has called for a need to develop drugs to treat this infection. Using in silico and molecular docking approaches, this study investigated the inhibitory effects of Pradimicin A, Lamivudine, Plerixafor and Lopinavir against SARS-CoV-2 M^pro. ADME/Tox of the ligands, pharmacophore hypothesis of the co-crystalized ligand and the receptor, and docking studies were carried out on different modules of Schrodinger (2019-4) Maestro v12.2. Among the ligands subjected to ADME/Tox by QikProp, Lamivudine demonstrated drug-like physico-chemical properties. A total of five pharmacophore binding sites (A3, A4, R9, R10, and R11) were predicted from the co-crystalized ligand and the binding cavity of the SARS-CoV-2 M^pro. The docking result showed that Lopinavir and Lamivudine bind with a higher affinity and lower free energy than the standard ligand having a glide score of -9.2 kcal/mol and -5.3 kcal/mol, respectively. Plerixafor and Pradimicin A have a glide score of -3.7 kcal/mol and -2.4 kcal/mol, respectively, which is lower than the co-crystallized ligand with a glide score of -5.3 kcal/mol. Molecular dynamics confirmed that the ligands maintained their interaction with the protein with lower RMSD fluctuations over the trajectory period of 100 nsecs and that GLU166 residue is pivotal for binding. On the whole, present study specifies the repurposing aptitude of these molecules as inhibitors of SARS-CoV-2 M^pro with higher binding scores and forms energetically stable complexes with M^pro.Communicated by Ramaswamy H. Sarma.

The Dual-Targeting Activity of the Metabolite Substrate of Para-amino Salicyclic Acid in the Mycobacterial Folate Pathway: Atomistic and Structural Perspectives

Therapeutic targeting of folate biosynthetic pathway has recently been explored as a viable strategy in the treatment of tuberculosis. The bioactive metabolite substrate of Para-amino salicyclic acid (PAS-M) reportedly dual-targets dihydrofolate reductase (DHFR) and flavin-dependent thymidylate synthase (FDTS), two essential enzymes in folate biosynthetic pathway. However, the molecular mechanisms and structural dynamics of this dual inhibitory activity of the PAS-M remain elusive. Molecular dynamics simulations revealed that binding of PAS-M towards DHFR is characterized by a recurrence of strong conventional hydrogen bond interactions between a peculiar DHFR binding site residue (Asp27) and the 2-amino-decahydropteridin-4-ol group of PAS-M. Similarly, the binding of PAS-M towards FDTS also involved consistent strong conventional hydrogen bond interactions between some specific residues (Tyr101, Arg172, Thr4, Gln103, Arg87 and Gln106) and, the 2-amino-decahydropteridin-4-ol group, thus establishing the cruciality of the group. Structural dynamics of the bound complexes of both enzymes revealed that, upon binding, PAS-M is anchored at the entrance of hydrophobic pockets by strong hydrogen bond interactions while the rest of the structure gains access to deeper hydrophobic residues to engage in favorable interactions. Further analysis of atomistic changes of both enzymes showed increased C-α atom deviations as well as an increase C-α atoms radius of gyration consistent with structural disorientations. These conformational changes possibly interfered with the biological functions of the enzymes and hence their inhibition as experimentally reported. Structural Insights provided could open up a novel paradigm of structure-based design of multi-targeting inhibitors of biological targets in the folate biosynthetic pathway toward tuberculosis therapy.

'Piperazining' the catalytic gatekeepers: unraveling the pan-inhibition of SRC kinases; LYN, FYN and BLK by masitinib

Aim: Blocking oncogenic signaling of B-cell receptor (BCR) has been explored as a viable strategy in the treatment of diffuse large B-cell lymphoma. Masitinib is shown to multitarget LYN, FYN and BLK kinases that propagate BCR signals to downstream effectors. However, the molecular mechanisms of its selectivity and pan-inhibition remain elusive. Materials & methods: This study therefore employed molecular dynamics simulations coupled with advanced post-molecular dynamics simulation techniques to unravel the structural mechanisms that inform the reported multitargeting ability of masitinib. Results: Molecular dynamics simulations revealed initial selective targeting of catalytic residues (Asp334/Glu335 - LYN; Asp130/Asp148/Glu54 - FYN; Asp89 - BLK) by masitinib, with high-affinity interactions via its piperazine ring at the entrance of the ATP-binding pockets, before systematic access into the hydrophobic deep pocket grooves. Conclusion: Identification of these 'gatekeeper' residues could open up a novel paradigm of structure-based design of highly selective pan-inhibitors of BCR signaling in the treatment of diffuse large B-cell lymphoma.

Solving the riddle: Unraveling the mechanisms of blocking the binding of leukotoxin by therapeutic antagonists in periodontal diseases

Aggregatibacter actinomycetemcomitans is a Gram-negative bacteria that has gained wide recognition for its causative role in the development of various immune diseases, which includes localized aggressive periodontitis. Its ability to evade host defense mechanisms is mediated by the secretion of leukotoxin (LtxA), which induces death of white blood cells (leukocytes) by specific binding to their surface-expressed leukocyte function-associated receptor (LFA-1) in its active state. Therapeutic compounds that interfere with this pathogenic process and abrogate A. actinomycetemcomitans virulence have been reported in literature. These include doxycycline, and more recently phytochemical compounds such as hamamelitanin, resveratrol, naringin, and quercetin. However, the question remains how do they work? Therefore, with the aid of computational tools, we explore the molecular mechanisms by which they possibly elicit their therapeutic functions. Molecular mechanics Poisson/Boltzmann surface area analyses revealed that these compounds bind favorably to active LFA-1 with high affinity and considerable stability, indicative of their ability to occupy the LtxA binding site (LBS) and prevent LtxA binding. The conformational transition of open LFA-1 to its closed state further describe the mechanistic activity of these compounds. In addition to notable reductions in structural mobility and flexibility, the burial of surface-exposed interactive side chains at the LBS was observed, an occurrence that could alter the complementary binding of LtxA. It is also important to mention that these occurrences were induced more prominently by the phytochemicals. We believe that these findings will enhance the scope of drug design and discovery for potent LtxA antagonists with improved activities and therapeutic efficacies in the treatment of virulent A. actinomycetemcomitans diseases.

An overview of recent molecular dynamics applications as medicinal chemistry tools for the undruggable site challenge

Molecular dynamics (MD) has become increasingly popular due to the development of hardware and software solutions and the improvement in algorithms, which allowed researchers to scale up calculations in order to speed them up. MD simulations are usually used to address protein folding issues or protein-ligand complex stability through energy profile analysis over time. In recent years, the development of new tools able to deeply explore a potential energy surface (PES) has allowed researchers to focus on the dynamic nature of the binding recognition process and binding-induced protein conformational changes. Moreover, modern approaches have been demonstrated to be effective and reliable in calculating some kinetic and thermodynamic parameters behind the host-guest recognition process. Starting from all of these considerations, several efforts have been made in order to integrate MD within the virtual screening process in drug discovery. Knowledge retrieved from MD can, in fact, be exploited as a starting point to build pharmacophores or docking constraints in the early stage of the screening campaign as well as to define key features, in order to unravel hidden binding modes and help the optimisation of the molecular structure of a lead compound. Based on these outcomes, researchers are nowadays using MD as an invaluable tool to discover and target previously considered undruggable binding sites, including protein-protein interactions and allosteric sites on a protein surface. As a matter of fact, the use of MD has been recognised as vital to the discovery of selective protein-protein interaction modulators. The use of a dynamic overview on how the host-guest recognition occurs and of the relative conformational modifications induced allows researchers to optimise small molecules and small peptides capable of tightly interacting within the cleft between two proteins. In this review, we aim to present the most recent applications of MD as an integrated tool to be used in the rational design of small molecules or small peptides able to modulate undruggable targets, such as allosteric sites and protein-protein interactions.

Across the blood-brain barrier: Neurotherapeutic screening and characterization of naringenin as a novel CRMP-2 inhibitor in the treatment of Alzheimer's disease using bioinformatics and computational tools

The discovery and developmental processes of CNS drugs have been limited by the inability of potential drug molecules to pass through the blood-brain barrier (BBB). This presents a significant setback in the treatment of neurodegenerative disorders such as Alzheimer's disease (AD), hence the need for compounds that can adhere strictly to the selective criteria of suitable CNS drugs. Collapsin response mediator protein-2 (CRMP-2) has been recently identified as a viable target in neurotherapeutics due to its involvement in the etiology of AD. As shown in previous studies, Naringenin (NAR), a small molecule derivative of Drynaria rhizome (DR) extract, specifically binds CRMP-2 and reduces its phosphorylation. This was shown to facilitate axonal regrowth, with improvement in cognition and learning. Herein, we report the first account of the use of cheminformatics techniques to define the CNS drug-suitability of NAR using selective criteria, coupled with the prediction of possible biological activities and toxicities. Also, we evaluated the mechanistic activity of NAR by modeling its molecular interaction with human CRMP-2 (hCRMP-2). Physicochemical analyses revealed the suitability of NAR as a CNS drug and its ability to transverse the BBB. Possible neurogenic, anti-carcinogenic and cardioprotective activities were also predicted. NAR exhibited favorable binding to CRMP-2 and formed strong bonds with active site residues, which accounts for its stabilization and affinity. Moreover, NAR induced notable conformational changes in CRMP-2, an occurrence that could possibly disrupt kinase-mediated phosphorylation. These findings will aid in the optimization of NAR and improve its neurotherapeutic activities in the treatment of AD.