The Application of Ligand-Mapping Molecular Dynamics Simulations to the Rational Design of Peptidic Modulators of Protein-Protein Interactions

Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms

In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.

Exploring Aromatic Cage Flexibility Using Cosolvent Molecular Dynamics Simulations─An In-Silico Case Study of Tudor Domains

Cosolvent molecular dynamics (MD) simulations have proven to be powerful in silico tools to predict hotspots for binding regions on protein surfaces. In the current study, the method was adapted and applied to two Tudor domain-containing proteins, namely Spindlin1 (SPIN1) and survival motor neuron protein (SMN). Tudor domains are characterized by so-called aromatic cages that recognize methylated lysine residues of protein targets. In the study, the conformational transitions from closed to open aromatic cage conformations were investigated by performing MD simulations with cosolvents using six different probe molecules. It is shown that a trajectory clustering approach in combination with volume and atomic distance tracking allows a reasonable discrimination between open and closed aromatic cage conformations and the docking of inhibitors yields very good reproducibility with crystal structures. Cosolvent MDs are suitable to capture the flexibility of aromatic cages and thus represent a promising tool for the optimization of inhibitors.

Success probability of high-affinity DNA aptamer generation by genetic alphabet expansion

Nucleic acid aptamers as antibody alternatives bind specifically to target molecules. These aptamers are generated by isolating candidates from libraries with random sequence fragments, through an evolutionary engineering system. We recently reported a high-affinity DNA aptamer generation method that introduces unnatural bases (UBs) as a fifth letter into the library, by genetic alphabet expansion. By incorporating hydrophobic UBs, the affinities of DNA aptamers to target proteins are increased over 100-fold, as compared with those of conventional aptamers with only the natural four letters. However, there is still plenty of room for improvement of the methods for routinely generating high-affinity UB-containing DNA (UB-DNA) aptamers. The success probabilities of the high-affinity aptamer generation depend on the existence of the aptamer candidate sequences in the initial library. We estimated the success probabilities by analysing several UB-DNA aptamers that we generated, as examples. In addition, we investigated the possible improvement of conventional aptamer affinities by introducing one UB at specific positions. Our data revealed that UB-DNA aptamers adopt specific tertiary structures, in which many bases including UBs interact with target proteins for high affinity, suggesting the importance of the UB-DNA library design. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Accelerated Ligand-Mapping Molecular Dynamics Simulations for the Detection of Recalcitrant Cryptic Pockets and Occluded Binding Sites

The identification and characterization of binding sites is a critical component of structure-based drug design (SBDD). Probe-based/cosolvent molecular dynamics (MD) methods that allow for protein flexibility have been developed to predict ligand binding sites. However, cryptic pockets that appear only upon ligand binding and occluded binding sites with no access to the solvent pose significant challenges to these methods. Here, we report the development of accelerated ligand-mapping MD (aLMMD), which combines accelerated MD with LMMD, for the detection of these challenging binding sites. The method was validated on five proteins with what we term "recalcitrant" cryptic pockets, which are deeply buried pockets that require extensive movement of the protein backbone to expose, and three proteins with occluded binding sites. In all the cases, aLMMD was able to detect and sample the binding sites. Our results suggest that aLMMD could be used as a general approach for the detection of such elusive binding sites in protein targets, thus providing valuable information for SBDD.

Application of Site-Identification by Ligand Competitive Saturation in Computer-Aided Drug Design

Site Identification by Ligand Competitive Saturation (SILCS) is a molecular simulation approach that uses diverse small solutes in aqueous solution to obtain functional group affinity patterns of a protein or other macromolecule. This involves employing a combined Grand Canonical Monte Carlo (GCMC)-molecular dynamics (MD) method to sample the full 3D space of the protein, including deep binding pockets and interior cavities from which functional group free energy maps (FragMaps) are obtained. The information content in the maps, which include contributions from protein flexibilty and both protein and functional group desolvation contributions, can be used in many aspects of the drug discovery process. These include identification of novel ligand binding pockets, including allosteric sites, pharmacophore modeling, prediction of relative protein-ligand binding affinities for database screening and lead optimization efforts, evaluation of protein-protein interactions as well as in the formulation of biologics-based drugs including monoclonal antibodies. The present article summarizes the various tools developed in the context of the SILCS methodology and their utility in computer-aided drug design (CADD) applications, showing how the SILCS toolset can improve the drug-development process on a number of fronts with respect to both accuracy and throughput representing a new avenue of CADD applications.

Computational Selectivity Assessment of Protease Inhibitors against SARS-CoV-2

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a serious global health threat. Since no specific therapeutics are available, researchers around the world screened compounds to inhibit various molecular targets of SARS-CoV-2 including its main protease (Mpro) essential for viral replication. Due to the high urgency of these discovery efforts, off-target binding, which is one of the major reasons for drug-induced toxicity and safety-related drug attrition, was neglected. Here, we used molecular docking, toxicity profiling, and multiple molecular dynamics (MD) protocols to assess the selectivity of 33 reported non-covalent inhibitors of SARS-CoV-2 Mpro against eight proteases and 16 anti-targets. The panel of proteases included SARS-CoV Mpro, cathepsin G, caspase-3, ubiquitin carboxy-terminal hydrolase L1 (UCHL1), thrombin, factor Xa, chymase, and prostasin. Several of the assessed compounds presented considerable off-target binding towards the panel of proteases, as well as the selected anti-targets. Our results further suggest a high risk of off-target binding to chymase and cathepsin G. Thus, in future discovery projects, experimental selectivity assessment should be directed toward these proteases. A systematic selectivity assessment of SARS-CoV-2 Mpro inhibitors, as we report it, was not previously conducted.

Straightforward Incorporation of Multiple Ligand Types into Molecular Dynamics Simulations for Efficient Binding Site Detection and Characterization

Binding site identification and characterization is an important initial step in structure-based drug design. To account for the effects of protein flexibility and solvation, several cosolvent molecular dynamics (MD) simulation methods that incorporate small organic molecules into the protein's solvent box to probe for binding sites have been developed. However, most of these methods are highly inefficient, as they allow for the use of only one probe type at a time, which means that multiple sets of simulations have to be performed to map different types of binding sites. The high probe concentrations used in some of these methods also necessitate the use of artificial repulsive forces to prevent the probes from aggregating. Here, we present multiple-ligand-mapping MD (mLMMD), a method that incorporates multiple types of probes for simultaneous and efficient mapping of different types of binding sites without the need for introduction of artificial forces that may cause unintended mapping artifacts. We validate the method on a diverse set of 10 proteins and show that the mLMMD probes are able to reliably identify hydrophobic, hydrogen-bonding, charged, and cryptic binding sites in all of the test cases. Our results also highlight the potential utility of mLMMD for virtual screening and rational drug design.

Roles of computational modelling in understanding p53 structure, biology, and its therapeutic targeting

The transcription factor p53 plays pivotal roles in numerous biological processes, including the suppression of tumours. The rich availability of biophysical data aimed at understanding its structure–function relationships since the 1990s has enabled the application of a variety of computational modelling techniques towards the establishment of mechanistic models. Together they have provided deep insights into the structure, mechanics, energetics, and dynamics of p53. In parallel, the observation that mutations in p53 or changes in its associated pathways characterize several human cancers has resulted in a race to develop therapeutic modulators of p53, some of which have entered clinical trials. This review describes how computational modelling has played key roles in understanding structural-dynamic aspects of p53, formulating hypotheses about domains that are beyond current experimental investigations, and the development of therapeutic molecules that target the p53 pathway.

Allosteric Binding Sites On Nuclear Receptors: Focus On Drug Efficacy and Selectivity

Nuclear receptors (NRs) are highly relevant drug targets in major indications such as oncologic, metabolic, reproductive, and immunologic diseases. However, currently, marketed drugs designed towards the orthosteric binding site of NRs often suffer from resistance mechanisms and poor selectivity. The identification of two superficial allosteric sites, activation function-2 (AF-2) and binding function-3 (BF-3), as novel drug targets sparked the development of inhibitors, while selectivity concerns due to a high conservation degree remained. To determine important pharmacophores and hydration sites among AF-2 and BF-3 of eight hormonal NRs, we systematically analyzed over 10 μ s of molecular dynamics simulations including simulations in explicit water and solvent mixtures. In addition, a library of over 300 allosteric inhibitors was evaluated by molecular docking. Based on our results, we suggest the BF-3 site to offer a higher potential for drug selectivity as opposed to the AF-2 site that is more conserved among the selected receptors. Detected similarities among the AF-2 sites of various NRs urge for a broader selectivity assessment in future studies. In combination with the Supplementary Material, this work provides a foundation to improve both selectivity and potency of allosteric inhibitors in a rational manner and increase the therapeutic applicability of this promising compound class.

Simulations of mutant p53 DNA binding domains reveal a novel druggable pocket

The DNA binding domain (DBD) of the tumor suppressor p53 is the site of several oncogenic mutations. A subset of these mutations lowers the unfolding temperature of the DBD. Unfolding leads to the exposure of a hydrophobic β-strand and nucleates aggregation which results in pathologies through loss of function and dominant negative/gain of function effects. Inspired by the hypothesis that structural changes that are associated with events initiating unfolding in DBD are likely to present opportunities for inhibition, we investigate the dynamics of the wild type (WT) and some aggregating mutants through extensive all atom explicit solvent MD simulations. Simulations reveal differential conformational sampling between the WT and the mutants of a turn region (S6-S7) that is contiguous to a known aggregation-prone region (APR). The conformational properties of the S6-S7 turn appear to be modulated by a network of interacting residues. We speculate that changes that take place in this network as a result of the mutational stress result in the events that destabilize the DBD and initiate unfolding. These perturbations also result in the emergence of a novel pocket that appears to have druggable characteristics. FDA approved drugs are computationally screened against this pocket.

Efficient development of stable and highly functionalised peptides targeting the CK2α/CK2β protein-protein interaction

The discovery of new Protein-Protein Interaction (PPI) modulators is currently limited by the difficulties associated with the design and synthesis of selective small molecule inhibitors. Peptides are a potential solution for disrupting PPIs; however, they typically suffer from poor stability in vivo and limited tissue penetration hampering their wide spread use as new chemical biology tools and potential therapeutics. In this work, a combination of CuAAC chemistry, molecular modelling, X-ray crystallography, and biological validation allowed us to develop highly functionalised peptide PPI inhibitors of the protein CK2. The lead peptide, CAM7117, prevents the formation of the holoenzyme assembly in vitro, slows down proliferation, induces apoptosis in cancer cells and is stable in human serum. CAM7117 could aid the development of novel CK2 inhibitors acting at the interface and help to fully understand the intracellular pathways involving CK2. Importantly, the approach adopted herein could be applied to many PPI targets and has the potential to ease the study of PPIs by efficiently providing access to functionalised peptides.

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.

Insights into the RNA binding mechanism of human L1-ORF1p: a molecular dynamics study

The recognition and binding of nucleic acids by ORF1p, an L1 retrotransposon protein, have not yet been clearly understood due to the lack of structural knowledge. The present study attempts to identify the probable single-stranded RNA binding pathway of trimeric ORF1p using computational methods like ligand mapping methodology combined with molecular dynamics simulations. Using the ligand mapping methodology, the possible RNA interacting sites on the surface of the trimeric ORF1p were identified. The crystal structure of the ORF1p timer and an RNA molecule of 29 nucleotide bases in length were used to generate the structure of the ORF1p complex based on information on predicted binding sites as well as the functional states of the CTD. The various complexes of ORF1p-RNA were generated using polyU, polyA and L1RNA sequences and were simulated for a period of 75 ns. The observed stable interaction pattern was used to propose the possible binding pathway. Based on the binding free energy for complex formation, both polyU and L1RNA complexes were identified as stable complexes, while the complex formed with polyA was the least stable one. Furthermore, the importance of the residues in the CC domain (Lys137 and Arg141), the RRM loop (Arg206, Arg210 and Arg211) and the CTD (Arg 261 and Arg262) of all three chains in stabilizing the wrapped RNA has been highlighted in this study. The presence of several electrostatic interactions including H-bond interactions increases the affinity towards RNA and hence plays a vital role in retaining the wrapped position of RNA around ORF1p. Altogether, this study presents one of the possible RNA binding pathways of ORF1p and clearly highlights the functional state of ORF1p visited during RNA binding.

Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short-read RNA sequencing, single molecule long-read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron-containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non-conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.

Targeting YAP/TAZ-TEAD protein-protein interactions using fragment-based and computational modeling approaches

The Hippo signaling pathway, which is implicated in the regulation of organ size, has emerged as a potential target for the development of cancer therapeutics. YAP, TAZ (transcription co-activators) and TEAD (transcription factor) are the downstream transcriptional machinery and effectors of the pathway. Formation of the YAP/TAZ-TEAD complex leads to transcription of growth-promoting genes. Conversely, disrupting the interactions of the complex decreases cell proliferation. Herein, we screened a 1000-member fragment library using Thermal Shift Assay and identified a hit fragment. We confirmed its binding at the YAP/TAZ-TEAD interface by X-ray crystallography, and showed that it occupies the same hydrophobic pocket as a conserved phenylalanine of YAP/TAZ. This hit fragment serves as a scaffold for the development of compounds that have the potential to disrupt YAP/TAZ-TEAD interactions. Structure-activity relationship studies and computational modeling were also carried out to identify more potent compounds that may bind at this validated druggable binding site.

Cosolvent-Based Molecular Dynamics for Ensemble Docking: Practical Method for Generating Druggable Protein Conformations

Protein flexibility is a major hurdle in current structure-based virtual screening (VS). In spite of the recent advances in high-performance computing, protein-ligand docking methods still demand tremendous computational cost to take into account the full degree of protein flexibility. In this context, ensemble docking has proven its utility and efficiency for VS studies, but it still needs a rational and efficient method to select and/or generate multiple protein conformations. Molecular dynamics (MD) simulations are useful to produce distinct protein conformations without abundant experimental structures. In this study, we present a novel strategy that makes use of cosolvent-based molecular dynamics (CMD) simulations for ensemble docking. By mixing small organic molecules into a solvent, CMD can stimulate dynamic protein motions and induce partial conformational changes of binding pocket residues appropriate for the binding of diverse ligands. The present method has been applied to six diverse target proteins and assessed by VS experiments using many actives and decoys of DEKOIS 2.0. The simulation results have revealed that the CMD is beneficial for ensemble docking. Utilizing cosolvent simulation allows the generation of druggable protein conformations, improving the VS performance compared with the use of a single experimental structure or ensemble docking by standard MD with pure water as the solvent.

Computational design of an epitope-specific Keap1 binding antibody using hotspot residues grafting and CDR loop swapping

Therapeutic and diagnostic applications of monoclonal antibodies often require careful selection of binders that recognize specific epitopes on the target molecule to exert a desired modulation of biological function. Here we present a proof-of-concept application for the rational design of an epitope-specific antibody binding with the target protein Keap1, by grafting pre-defined structural interaction patterns from the native binding partner protein, Nrf2, onto geometrically matched positions of a set of antibody scaffolds. The designed antibodies bind to Keap1 and block the Keap1-Nrf2 interaction in an epitope-specific way. One resulting antibody is further optimised to achieve low-nanomolar binding affinity by in silico redesign of the CDRH3 sequences. An X-ray co-crystal structure of one resulting design reveals that the actual binding orientation and interface with Keap1 is very close to the design model, despite an unexpected CDRH3 tilt and V_H/V_L interface deviation, which indicates that the modelling precision may be improved by taking into account simultaneous CDR loops conformation and V_H/V_L orientation optimisation upon antibody sequence change. Our study confirms that, given a pre-existing crystal structure of the target protein-protein interaction, hotspots grafting with CDR loop swapping is an attractive route to the rational design of an antibody targeting a pre-selected epitope.

Benzene Probes in Molecular Dynamics Simulations Reveal Novel Binding Sites for Ligand Design

Protein flexibility poses a major challenge in binding site identification. Several computational pocket detection methods that utilize small-molecule probes in molecular dynamics (MD) simulations have been developed to address this issue. Although they have proven hugely successful at reproducing experimental structural data, their ability to predict new binding sites that are yet to be identified and characterized has not been demonstrated. Here, we report the use of benzenes as probe molecules in ligand-mapping MD (LMMD) simulations to predict the existence of two novel binding sites on the surface of the oncoprotein MDM2. One of them was serendipitously confirmed by biophysical assays and X-ray crystallography to be important for the binding of a new family of hydrocarbon stapled peptides that were specifically designed to target the other putative site. These results highlight the predictive power of LMMD and suggest that predictions derived from LMMD simulations can serve as a reliable basis for the identification of novel ligand binding sites in structure-based drug design.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Promiscuous proteins are commonly observed in biological systems, such as modular domains that recognize phosphopeptides during signal transduction. The use of phosphopeptides and compounds with phosphate groups as inhibitors to protein–protein interactions have attracted increasing interest for years. By using atomistic molecular dynamics simulations, we are able to perform detailed analyses of the dihedral space to explore protein fluctuation upon ligand binding to better understand promiscuous molecular recognition. Free energy calculation can further provide insights into the mechanism of binding, including both enthalpic and entropic contributions for molecular recognition, which assist in inhibitor design. Our calculation results show that pre-rigidifying a ligand is not always advantageous, suggesting the challenge in retaining optimized intermolecular interactions in pre-rigidified ligand. Instead, certain flexible ligands with multiple binding conformations can reduce entropic penalty, and therefore improves binding affinity. According to our computations, we can introduce new intermolecular interactions to flexible ligand to strengthen attractions while maintaining small entropic penalty by retaining its plasticity in the bound conformation. The study might cast light on a new general strategy for designing inhibitors targeting promiscuous modular domains and protein–protein interactions.

Driving Structure-Based Drug Discovery through Cosolvent Molecular Dynamics

Identifying binding hotspots on protein surfaces is of prime interest in structure-based drug discovery, either to assess the tractability of pursuing a protein target or to drive improved potency of lead compounds. Computational approaches to detect such regions have traditionally relied on energy minimization of probe molecules onto static protein conformations in the absence of the natural aqueous environment. Advances in high performance computing now allow us to assess hotspots using molecular dynamics (MD) simulations. MD simulations integrate protein flexibility and the complicated role of water, thereby providing a more realistic assessment of the complex kinetics and thermodynamics at play. In this review, we describe the evolution of various cosolvent-based MD techniques and highlight a myriad of potential applications for such technologies in computational drug development.

Small-Molecule and Peptide Inhibitors of the Pro-Survival Protein Mcl-1

The ability of protein-protein interactions to regulate cellular processes in both beneficial and detrimental ways has made them obvious drug targets. The Bcl-2 family of proteins undergo a series of protein-protein interactions which regulate the intrinsic cell-death pathway. The pro-survival members of the Bcl-2 family, including Bcl-2, Bcl-xL , and Mcl-1, are commonly overexpressed in a number of human cancers. Effective modulators of members of the Bcl-2 family have been developed and are undergoing clinical trials, but the efficient modulation of Mcl-1 is still not represented in the clinic. In addition, Mcl-1 is a major cause of resistance to radio- and chemotherapies, including inhibitors that target other Bcl-2 family members. Subsequently, the inhibition of Mcl-1 has become of significant interest to the scientific community. This review covers the progress made to date in modulating the activity of Mcl-1, by both stapled peptides and small molecules. The development of peptides as drug candidates, and the advancement of experimental and computational techniques used to discover small molecules are also highlighted.

Identification of Small-Molecule Inhibitors of the Antiapoptotic Protein Myeloid Cell Leukaemia-1 (Mcl-1)

Protein-protein interactions (PPIs) control many cellular processes in cancer and tumour growth. Of significant interest is the role PPIs play in regulating apoptosis. The overexpression of the antiapoptosis regulating Bcl-2 family of proteins is commonly observed in several cancers, leading to resistance towards both radiation and chemotherapies. From this family, myeloid cell leukemia-1 (Mcl-1) has proven the most difficult to target, and one of the leading causes of treatment resistance. Exploiting the selective PPI between the apoptosis-regulating protein Noxa and Mcl-1, utilising a fluorescence polarization assay, we have identified four small molecules with the ability to modulate Mcl-1. The identified compounds were computationally modelled and docked against the Mcl-1 binding interface to obtain structural information about their binding sites allowing for future analogue design. When examined for their activity towards pancreatic cell lines that overexpress Mcl-1 (MiaPaCa-2 and BxPC-3), the identified compounds demonstrated growth inhibition, suggesting effective Mcl-1 modulation.

Small Molecule Targeting of Protein-Protein Interactions through Allosteric Modulation of Dynamics

The protein–protein interaction (PPI) target class is particularly challenging, but offers potential for “first in class” therapies. Most known PPI small molecules are orthosteric inhibitors but many PPI sites may be fundamentally intractable to this approach. One potential alternative is to consider more attractive, remote small molecule pockets; however, on the whole, allostery is poorly understood and difficult to discover and develop. Here we review the literature in order to understand the basis for allostery, especially as it can apply to PPIs. We suggest that the upfront generation of sophisticated and experimentally validated dynamic models of target proteins can aid in target choice and strategy for allosteric intervention to produce the required functional effect.