Insights into Mcl-1 Conformational States and Allosteric Inhibition Mechanism from Molecular Dynamics Simulations, Enhanced Sampling, and Pocket Crosstalk Analysis

Targeting Myeloid Leukemia-1 in Cancer Therapy: Advances and Directions

As a tripartite cell death switch, B-cell lymphoma protein 2 (Bcl-2) family members precisely regulate the endogenous apoptosis pathway in response to various cell signal stresses through protein-protein interactions. Myeloid leukemia-1 (Mcl-1), a key anti-apoptotic Bcl-2 family member, is positioned downstream in the endogenous apoptotic pathway and plays a central role in regulating mitochondrial function. Mcl-1 is highly expressed in a variety of hematological malignancies and solid tumors, contributing to tumorigenesis, poor prognosis, and chemoresistance, making it an attractive target for cancer treatment. This Perspective aims to discuss the mechanism by which Mcl-1 regulates apoptosis and non-apoptotic functions in cancer cells and to outline the discovery and optimization process of potent Mcl-1 modulators. In addition, we summarize the structural characteristics of potent inhibitors that bind to Mcl-1 through multiple co-crystal structures and analyze the cardiotoxicity caused by current Mcl-1 inhibitors, providing prospects for rational targeting of Mcl-1.

In silico design of potential Mcl-1 peptide-based inhibitors

CONTEXT

BIM (Bcl-2 interacting mediator of apoptosis)-derived peptides that specifically target over-expressed Mcl-1 (myeloid cell leukemia-1) protein and induce apoptosis are potentially anti-cancer agents. Since the helicity of BIM-derived peptides has a crucial role in their functionality, a range of strategies have been used to increase the helicity including the introduction of unnatural residues and stapling methods that have some drawbacks such as the accumulation in the liver. To avoid these drawbacks, this study aimed to design a more helical peptide by utilizing bioinformatics algorithms and molecular dynamics simulations without exploiting unnatural residues and stapling methods. MM-PBSA results showed that the mutations of A4fE and A2eE in analogue 5 demonstrate a preference towards binding with Mcl-1. As evidenced by Circular dichroism results, the helicity increases from 18 to 34%, these findings could enhance the potential of analogue 5 as an anti-cancer agent targeting Mcl-1. The applied strategies in this research could shed light on the in silico peptide design. Moreover, analogue 5 as a drug candidate can be evaluated in vitro and in vivo studies.

METHODS

The sequence of the lead peptide was determined using the ApInAPDB database and PRALINE program. Contact finder and PDBsum web server softwares were used to determine the contact involved amino acids in complex with Mcl-1. All identified salt bridge contributing residues were unaltered to preserve the binding affinity. After proposing novel analogues, their secondary structures were predicted by Cham finder web server software and GOR, Neural Network, and Chou-Fasman algorithms. Finally, molecular dynamics simulations run for 100 ns were done using the GROMACS, version 5.0.7, with the CHARMM36 force field. MM-PBSA was used to assess binding affinity specificity in targeting Mcl-1 and Bcl-xL (B-cell lymphoma extra-large).

Molecular dynamics simulations suggest Thiosemicarbazones can bind p53 cancer mutant R175H

Cancer pathologies are associated with the unfolding and aggregation of most recurring mutations in the DNA Binding Domain (DBD) of p53 that coordinate the destabilization of protein. Substitution at the 175th codon with arginine to histidine (R175H, a mutation of large to small side-chain amino acid) destabilizes the DBD by 3 kcal/mol and triggers breasts, lung cancer, etc. Stabilizing the p53 mutant by small molecules offers an attractive drug-targeted anti-cancer therapy. The thiosemicarbazone (TSC) molecules NPC and DPT are known to act as zinc-metallochaperones to reactivate p53R175H. Here, a combination of LESMD simulations for 10 TSC conformations with a p53R175H receptor, single ligand-protein conformation MD, and ensemble docking with multiple p53R175H conformations observed during simulations is suggested to identify the potential binding site of the target protein in light of their importance for the direct TSC - p53R175H binding. NPC binds mutant R175H in the loop region L2-L3, forming pivotal hydrogen bonds with HIS175, pi‑sulfur bonds with TYR163, and pi-alkyl linkages with ARG174 and PRO190, all of which are contiguous to the zinc-binding native site on p53DBD. DPT, on the other hand, was primarily targeting alternative binding sites such as the loop-helix L1/H2 region and the S8 strand. The similar structural characteristics of TSC-bound p53R175H complexes with wild-type p53DBD are thought to be attributable to involved interactions that favour binding free energy contributions of TSC ligands. Our findings may be useful in the identification of novel pockets with druggable properties.

MCL-1 promiscuity and the structural resilience of its binding partners

The allosteric protein MCL-1 and its natural inhibitors, the BH3-only proteins PUMA, BIM, and NOXA regulate apoptosis by interacting promiscuously within an entangled binding network. Little is known about the transient processes and dynamic conformational fluctuations that are the basis for the formation and stability of the MCL-1/BH3-only complex. In this study, we designed photoswitchable versions of MCL-1/PUMA and MCL-1/NOXA, and investigated the protein response after an ultrafast photo-perturbation with transient infrared spectroscopy. We observed partial α-helical unfolding in all cases, albeit on strongly varying timescales (1.6 ns for PUMA, 9.7 ns for the previously studied BIM, and 85 ns for NOXA). These differences are interpreted as a BH3-only-specific "structural resilience" to defy the perturbation while remaining in MCL-1's binding pocket. Thus, the presented insights could help to better understand the differences between PUMA, BIM, and NOXA, the promiscuity of MCL-1, in general, and the role of the proteins in the apoptotic network.

Myeloid cell leukemia-1: a formidable barrier to anticancer therapeutics and the quest of targeting it

The antiapoptotic B cell lymphoma-2 (Bcl-2) family members are apical regulators of the intrinsic pathway of apoptosis that orchestrate mitochondrial outer membrane permeabilization (MOMP) through interactions with their proapoptotic counterparts. Overexpression of antiapoptotic Bcl-2 family proteins has been linked to therapy resistance and poor prognosis in diverse cancers. Among the antiapoptotic Bcl-2 family members, predominant overexpression of the prosurvival myeloid cell leukemia-1 (Mcl-1) has been reported in a myriad of hematological malignancies and solid tumors, contributing to therapy resistance and poor outcomes, thus making it a potential druggable target. The unique structure of Mcl-1 and its complex regulatory mechanism makes it an adaptive prosurvival switch that ensures tumor cell survival despite therapeutic intervention. This review focusses on diverse mechanisms adopted by tumor cells to maintain sustained elevated levels of Mcl-1 and how high Mcl-1 levels contribute to resistance in conventional as well as targeted therapies. Moreover, recent developments in the Mcl-1-targeted therapeutics and the underlying challenges and considerations in designing novel Mcl-1 inhibitors are also discussed.

Divide and Conquer. Pocket-Opening Mixed-Solvent Simulations in the Perspective of Docking Virtual Screening Applications for Drug Discovery

The existence of a druggable binding pocket is a prerequisite for computational drug-target interaction studies including virtual screening. Retrospective studies have shown that extended sampling methods like Markov State Modeling and mixed-solvent simulations can identify cryptic pockets relevant for drug discovery. Here, we apply a combination of mixed-solvent molecular dynamics (MD) and time-structure independent component analysis (TICA) to four retrospective case studies: NPC2, the CECR2 bromodomain, TEM-1, and MCL-1. We compare previous experimental and computational findings to our results. It is shown that the successful identification of cryptic pockets depends on the system and the cosolvent probes. We used alternative TICA internal features such as the unbiased backbone coordinates or backbone dihedrals versus biased interatomic distances. We found that in the case of NPC2, TEM-1, and MCL-1, the use of unbiased features is able to identify cryptic pockets, although in the case of the CECR2 bromodomain, more specific features are required to properly capture a pocket opening. In the perspective of virtual screening applications, it is shown how docking studies with the parent ligands depend critically on the conformational state of the targets.

Drug Repurposing: Deferasirox Inhibits the Anti-Apoptotic Activity of Mcl-1

Introduction

With the aim of repositioning commercially available drugs for the inhibition of the anti-apoptotic myeloid cell leukemia protein, Mcl-1, implied in various cancers, five molecules, highlighted from a published theoretical screening, were selected to experimentally validate their affinity toward Mcl-1.

Results

A detailed NMR study revealed that only two of the five tested drugs, Torsemide and Deferasirox, interacted with Mcl-1. NMR data analysis allowed the complete characterization of the binding mode of both drugs to Mcl-1, including the estimation of their affinity for Mcl-1. Biological assays evidenced that the biological activity of Torsemide was lower as compared to the Deferasirox, which was able to efficiently and selectively inhibit the anti-apoptotic activity of Mcl-1. Finally, docking and molecular dynamics led to a 3D model for the Deferasirox:Mcl-1 complex and revealed the positioning of the drug in the Mcl-1 P2/P3 pockets as well as almost all synthetic Mcl-1 inhibitors. Interestingly, contrary to known synthetic Mcl-1 inhibitors which interact through Arg263, Deferasirox, establishes a salt bridge with Lys234.

Conclusion

Deferasirox could be a potential candidate for drug repositioning as Mcl-1 inhibitor.

Cryptic Pockets Repository through Pocket Dynamics Tracking and Metadynamics on Essential Dynamics Space: Applications to Mcl-1

Detection of cryptic pockets (hidden protein pockets) is a hot topic in structure-based drug discovery, especially for drugging the yet undruggable proteome. The experimental detection of cryptic pockets is still considered an expensive endeavor. Thus, computational methods, such as atomistic simulations, are used instead. These simulation methods can provide a perspective on protein dynamics that overpasses the experimental X-ray structures' static and average view. Nonetheless, unbiased molecular dynamics (MD) simulations fall short to detect transient and cryptic pockets requiring the crossing of high-energy barriers. Enhanced sampling methods, such as Metadynamics, provide a solution to overcome the time-scale problem faced by unbiased MD simulations. However, these methods are still limited by the availability of collective variable space to capture the intricate parameters, leading to the opening of cryptic pockets. Unfortunately, the design of such collective variables requires a priori knowledge of the binding site, information that is by definition lacking for cryptic pockets. In this work, we evaluated the use of the Metadynamics biasing scheme on essential coordinates space as a general method for cryptic pocket detection. This approach was applied to an antiapoptotic protein: Mcl-1 as a test model. In addition to providing a broader characterization of Mcl-1's conformational space, we show the effectiveness of this method in drawing the full repository of Mcl-1's known and novel cryptic pockets in an unsupervised manner.

Myeloid cell leukemia 1 (MCL-1): Structural characteristics and application in cancer therapy

Apoptosis, a major hallmark of cancer cells, regulates cellular fate and homeostasis. BCL-2 (B-cell CLL/Lymphoma 2) protein family is popularly known to mediate the intrinsic mode of apoptosis, of which MCL-1 is a crucial member. Myeloid cell leukemia 1 (MCL-1) is an anti-apoptotic oncoprotein and one of the most investigated members of the BCL-2 family. It is commonly known to be genetically altered, aberrantly overexpressed, and primarily associated with drug resistance in various human cancers. Recent advancements in the development of selective MCL-1 inhibitors and evaluating their effectiveness in cancer treatment establish its popularity as a molecular target. The overall aim is the selective induction of apoptosis in cancer cells by using a single or combination of BCL-2 family inhibitors. Delineating the precise molecular mechanisms associated with MCL-1-mediated cancer progression will certainly improve the efficacy of clinical interventions aimed at MCL-1 and hence patient survival. This review is structured to highlight the structural characteristics of MCL-1, its specific interactions with NOXA, MCL-1-regulatory microRNAs, and at the same time focus on the emerging therapeutic strategies targeting our protein of interest (MCL-1), alone or in combination with other treatments.