The protein c-Myc is a transcription factor that remains largely intrinsically disordered and is known to be involved in various biological processes and is overexpressed in various cancers, making it an attractive drug target. However, intrinsically disordered proteins such as c-Myc do not show funnel-like basins in their free-energy landscapes; this makes their druggability a challenge. For the first time, we propose a heterodimer model of c-Myc/Max in full length in this work. We used Gaussian-accelerated molecular dynamics (GaMD) simulations to explore the behavior of c-Myc and its various regions, including the transactivation domain (TAD) and the basic helix-loop-helix-leucine-zipper (bHLH-Zipper) motif in three different conformational states: (a) monomeric c-Myc, (b) c-Myc when bound to its partner protein, Max, and (c) when Max was removed after binding. We analyzed the GaMD trajectories using root-mean-square deviation (RMSD), radius of gyration, root-mean-square fluctuation, and free-energy landscape (FEL) calculations to elaborate the behaviors of these regions. The results showed that the monomeric c-Myc structure showed a higher RMSD fluctuation as compared with the c-Myc/Max heterodimer in the bHLH-Zipper motif. This indicated that the bHLH-Zipper motif of c-Myc is more stable when it is bound to Max. The TAD region in both monomeric and Max-bound states showed similar plasticity in terms of RMSD. We also conducted residue decomposition calculations and showed that the c-Myc and Max interaction could be driven mainly by electrostatic interactions and the residues Arg299, Ile403, and Leu420 seemed to play important roles in the interaction. Our work provides insights into the behavior of c-Myc and its regions that could support the development of drugs that target c-Myc and other intrinsically disordered proteins.
Polo-like kinase 1 (PLK1) and p38γ mitogen-activated protein kinase (p38γ) play important roles in cancer pathogenesis by controlling cell cycle progression and are therefore attractive cancer targets. The design of multitarget inhibitors may offer synergistic inhibition of distinct targets and reduce the risk of drug-drug interactions to improve the balance between therapeutic efficacy and safety. We combined deep-learning-based quantitative structure-activity relationship (QSAR) modeling and hybrid-based consensus scoring to screen for inhibitors with potential activity against the targeted proteins. Using this combination strategy, we identified a potent PLK1 inhibitor (compound 4) that inhibited PLK1 activity and liver cancer cell growth in the nanomolar range. Next, we deployed both our QSAR models for PLK1 and p38γ on the Enamine compound library to identify dual-targeting inhibitors against PLK1 and p38γ. Likewise, the identified hits were subsequently subjected to hybrid-based consensus scoring. Using this method, we identified a promising compound (compound 14) that could inhibit both PLK1 and p38γ activities. At nanomolar concentrations, compound 14 inhibited the growth of human hepatocellular carcinoma and hepatoblastoma cells in vitro. This study demonstrates the combined screening strategy to identify novel potential inhibitors for existing targets.
Lipoteichoic acid synthase (LtaS) is a key enzyme for the cell wall biosynthesis of Gram-positive bacteria. Gram-positive bacteria that lack lipoteichoic acid (LTA) exhibit impaired cell division and growth defects. Thus, LtaS appears to be an attractive antimicrobial target. The pharmacology around LtaS remains largely unexplored with only two small-molecule LtaS inhibitors reported, namely "compound 1771" and the Congo red dye. Structure-based drug discovery efforts against LtaS remain unattempted due to the lack of an inhibitor-bound structure of LtaS. To address this, we combined the use of a molecular docking technique with molecular dynamics (MD) simulations to model a plausible binding mode of compound 1771 to the extracellular catalytic domain of LtaS (eLtaS). The model was validated using alanine mutagenesis studies combined with isothermal titration calorimetry. Additionally, lead optimization driven by our computational model resulted in an improved version of compound 1771, namely, compound 4 which showed greater affinity for binding to eLtaS than compound 1771 in biophysical assays. Compound 4 reduced LTA production in S. aureus dose-dependently, induced aberrant morphology as seen for LTA-deficient bacteria, and significantly reduced bacteria titers in the lung of mice infected with S. aureus. Analysis of our MD simulation trajectories revealed the possible formation of a transient cryptic pocket in eLtaS. Virtual screening (VS) against the cryptic pocket led to the identification of a new class of inhibitors that could potentiate β-lactams against methicillin-resistant S. aureus. Our overall workflow and data should encourage further drug design campaign against LtaS. Finally, our work reinforces the importance of considering protein conformational flexibility to a successful VS endeavor.