Leveraging the Fragment Molecular Orbital Method to Explore the PLK1 Kinase Binding Site and Polo-Box Domain for Potent Small-Molecule Drug Design

Polo-like kinase 1 (PLK1) plays a pivotal role in cell division regulation and emerges as a promising therapeutic target for cancer treatment. Consequently, the development of small-molecule inhibitors targeting PLK1 has become a focal point in contemporary research. The adenosine triphosphate (ATP)-binding site and the polo-box domain in PLK1 present crucial interaction sites for these inhibitors, aiming to disrupt the protein's function. However, designing potent and selective small-molecule inhibitors can be challenging, requiring a deep understanding of protein-ligand interaction mechanisms at these binding sites. In this context, our study leverages the fragment molecular orbital (FMO) method to explore these site-specific interactions in-depth. Through the FMO approach, we used the FMO method to elucidate the molecular mechanisms of small-molecule drug binding to these sites to design PLK1 inhibitors that are both potent and selective. Our investigation further entailed a comparative analysis of various PLK1 inhibitors, each characterized by distinct structural attributes. This comparison enhanced our understanding of the structure-activity relationships of these inhibitors, underscoring the efficacy of the FMO method in identifying critical binding features and predicting binding modes for small-molecule ligands. Furthermore, Our research spotlighted "hot spot" residues instrumental for selective and robust PLK1 binding. Our findings offer profound insights, priming the rational design of innovative potential PLK1 inhibitors with significant implications for developing anticancer therapeutics.

Keywords:

Subject:

Biology and Life Sciences - Biology and Biotechnology

Supplementary Material

supplementary.docx (1.31MB )

1. Introduction:

Polo-like kinases (PLKs) form a family of serine/tyrosine kinase proteins with wide distribution in eukaryotic cells and play crucial roles in various cell-cycle phases. Currently, the PLK protein family comprises five members: PLK1, PLK2, PLK3, PLK4, and PLK5 (Figure 1(A)). Among these members, PLK1 has undergone the most comprehensive research to understand the regulatory mechanisms influencing its functions and potential as a target for drug design [1]. Generally, PLK1 exhibits a structure akin to other kinase family members, comprising an N-terminal serine/threonine kinase domain (KD) and a C-terminal repetition of the polobox domain (PBD) [2]. The enzymatic activity of PLK1 is directly influenced by its phosphorylation [3]. In addition, PLK1 contains two distinct druggable domains that can be used for the inhibitory design. Moreover, the expression pattern of PLK1 is linked to mitotic progression. [4]. In particular, PLK1 is typically expressed at low levels during interphase, and its expression gradually increases as the cell enters the G2/M phase, reaching its peak at this stage [5]. Following mitosis, PLK1 undergoes significant degradation, rapidly decreasing its protein levels [6]. During the cell cycle, PLK1 is crucial in regulating various processes. These include checkpoint recovery [7], the timing of mitotic entry [8], centrosome maturation, bipolar spindle assembly [9], microtubule-kinetochore attachment stabilization [10], and proper chromosomal segregation during anaphase [11].

Cell-cycle disruption is a prominent characteristic of various cancer types [12]. In malignant cells, there is an upregulation of PLK1 expression, leading to multiple defects in mitosis and cytokinesis. This heightened PLK1 activity increases chromosomal instability [13], often associated with a high tumor grade and unfavorable patient prognosis [12]. Beyond its role in cell-cycle dysregulation, PLK1 significantly promotes cancer progression through metabolic reprogramming, including enhancing the flux of the pentose phosphate pathway and directing glucose to pathways involved in the synthesis of macromolecules [14]. Furthermore, the downregulation of PLK1 expression typically reduces the proliferation of diverse cancer cells [15]. Thus, PLK1 is widely regarded as a potent proto-oncogene and a promising target for cancer therapy.

The PLK1 protein can be divided into two primary regions: the N-terminal kinase catalytic domain, which is highly conserved and spans approximately 252 amino acids to the 303 amino acids position, and the unique C-terminal PBD comprising two PBDs, each containing 60 to 70 amino acids [16]. The N-terminal catalytic KD (amino acid residues 53 to 303) possesses characteristic residues found in conserved serine/threonine kinases, including a T-loop and an adenosine triphosphate (ATP)-binding cassette. This domain contains crucial residues, such as Lys82, Cys133, and Asp194, which are vital for ATP binding. Additionally, residues 134 to 136 form a functional bipartite nuclear localization signal sequence within this domain [17]. During the G2 phase of the cell cycle, the conserved residue 210 in the T-loop of PLK1 undergoes phosphorylation. This phosphorylation event is conducted by the upstream kinase Aurora A and its cofactor Bora, activating PLK1 [18]. The ATP-binding pocket of PLK1 consists of specific residues: Phe183 (at the bottom of the binding site), Leu59, Cys67, and Ala80 (at the top), effectively sandwiching the surface of the adenine ring system. The gatekeeper residues (Leu130 and Val129) are at the back of the adenine-binding pocket. Additionally, the PBD is elucidated in crystalline conjugates with phosphopeptides [19].

The C-terminal region of PLK1 (residues 345 to 603) primarily consists of two distinctive PBDs. The polocap located at the N-terminus of polo-box (PB) 1 can fold around PB 2, restricting its movement. Although PB 1 and 2 exhibit limited homology, their three-dimensional (3D) structures display remarkable similarities [16]. These proteins are pivotal in regulating PLK1 catalytic activity, localization, and substrate binding [20]. By binding to the PLK1 KD, these proteins act together to inhibit PLK1 kinase activity by inducing conformational changes in the protein. The interaction of the PLK1 PBD with Map205 stabilizes the autoinhibited state of PLK1 and sequesters it from its substrates [19]. Collectively, the flexible KD with its hinge region, along with the T-loop and the mentioned proteins, forms a switch that facilitates intricate spatiotemporal regulation of PLK1 activity [21]. The conformation of the PBD is crucial to its interaction with a phosphopeptide, enabling the main chains to interact and form a short antiparallel β-sheet between the peptide and PBD strand. Four specific residues, Trp414, Leu490, His538, and Lys540, significantly mediate other interactions. Particularly noteworthy are His538 and Lys540, which interact with the threonine phosphate group [22].

Comprising two domains, the KD and PBD, PLK1 offers inherent targets for PLK1 inhibitor development. Recently, researchers have made significant progress in creating numerous novel small-molecule PLK1 inhibitors, encompassing both KD and PBD inhibitors [23]. The ATP-competitive inhibitors target the deep groove in the kinase ATP-binding domain. At present, over 10 commercially available PLK1-specific inhibitors exist, with at least four of them (BI2536, BI6727 (Volasertib), GSK461364, and NMS-1286937 (Onvansertib)) having undergone evaluation in clinical trials (Table 1) [24]. All four PLK1 inhibitors function as ATP-competitive inhibitors and exhibit a similar mode of action. Among them, BI2536, a specific PLK1 inhibitor, has been extensively studied as a cytotoxic drug for treating various cancer types [25]. Nevertheless, BI2536 restricts antitumor activity during clinical trials and exhibits dose-limiting side effects [26]. In a recent study, the reduced effectiveness of BI2536 against progressive hepatocellular carcinoma was linked to low intratumoral drug levels [27]. Also known as Volasertib, BI6727 has exhibited great promise as a PLK1 inhibitor. Several preclinical experiments have demonstrated its effectiveness in inducing tumor regression [28]. Onvansertib, an orally administered and remarkably potent PLK1 inhibitor, exhibits high specificity for PLK1, excellent bioavailability, and a short half-life. Notably, it can potentially mitigate the toxicity of previous nonspecific PLK1 inhibitors [29].

Despite the promise of ATP-competitive inhibitors, resistance often emerges in cancer patients due to the high conservation of ATP-binding domains among various kinases and frequent mutations in ATP-binding sites [30]. Furthermore, ATP-competitive inhibitors may exhibit activity against other kinases, leading to a lack of specificity for PLK1.

In contrast, the PBD is unique to PLKs, making it a more viable and specific target for developing selective PLK1 inhibitors. Three selective PLK1 inhibitors (poloxin, thymoquinone, and purpurogallin) target the PBD of PLK1 [31]. Both Poloxin and thymoquinone can obstruct the correct orientation of PLK1, effectively impeding the mitosis of cancer cells [32]. A recent study revealed that poloxin-2, an optimized analog of poloxin, exhibits substantially enhanced potency and selectivity compared to poloxin. This improvement efficiently induces mitotic arrest and apoptosis in cultured human cancer cells [33].

Although no small-molecule PLK1 inhibitor has been approved for marketing yet, at least 10 drugs have entered clinical trials (Table 1). These drugs have undergone evaluation in various cancer treatments, revealing promising progress. However, most are still in early clinical trials, and certain caveats regarding their potential use have been reported.

In contrast, the PBD is exclusive to PLKs, suggesting it could be a more appropriate target for developing selective PLK1 inhibitors. Poloxin, thymoquinone, and purpurogallin are specific PLK1 inhibitors that target the PBD of PLK1 [31]. Poloxin and thymoquinone can disrupt the proper orientation of PLK1, impeding the mitosis of cancer cells [32]. A recent study found that poloxin-2, an optimized analog of poloxin, demonstrates substantially enhanced potency and selectivity compared to poloxin when inducing mitotic arrest and apoptosis in cultured human cancer cells [32]. As mentioned previously, PLK1 inhibitors have dose-limiting toxicities. Therefore, various approaches such as bivalent inhibitors were introduced. Andrej and his colleagues successfully identified the initial bifunctional inhibitors of PLK1, which are bridge kinase inhibitors and PBD peptides [34].

In 1999, Kitaura et al. [35] introduced the Fragment Molecular Orbital (FMO) method, partitioning a target molecule into smaller fragments for molecular orbital calculations [35,36]. The FMO method has recently gained popularity in new drug designs due to its accurate energy calculations, particularly in analyzing protein-ligand and protein-protein interactions (PPIs) [37,38]. By employing the FMO method, we can calculate fragment-fragment interaction energies, providing valuable insight into protein-ligand interactions. The energies of fragment-fragment pairs are called interfragment interaction energy. Nakano et al. developed packages in ABINIT-MP that implement the FMO method [39]. The software suite incorporates pair interaction energy (PIE) decomposition analysis (PIEDA) to estimate inter-fragmental interactions between functional group units based on their distinct contributions from electrostatic and dispersion forces.

Recently, the FMO method has gained recognition as a powerful tool for analyzing biological interactions between ligands and proteins in various contexts, such as estrogen receptors [40], human immunodeficiency virus proteases [41], influenza neuraminidases [42], G protein-coupled receptors, serine/threonine protein kinase Pim-1 [43], and YAP binding transcription factor TEAD [44]. This study involved an FMO/3D scattered PIE (SPIE) analysis to explore interaction systems at the ATP-binding site. Additionally, to investigate PPIs, the FMO/3D-SPIE analysis was introduced to identify crucial interactions in previous work. This analysis efficiently correlated the results with experimental site-directed mutagenesis findings [45].

In this study, we employed the FMO/3D-SPIE analysis to investigate PPIs between the substrate peptide and the PBD of PLK1 at the quantum mechanical level. We utilized the FMO method to analyze interactions within the ATP-binding site, as well as with PBDs, as part of our PLK1 hot-spot analysis. This research investigates the ATP and substrate peptide binding sites of the KD and PBD of PLK1, respectively. We explored the interactions of inhibitors BI2536, Onvansertib, and GSK461364, shedding light on the key residues responsible for selectivity. Our FMO analysis furnished insights into their respective binding modes. Furthermore, our study emphasized the significance of particular pockets in PBD binding. Molecular dynamic (MD) simulations and solvation analysis further substantiated our findings, pinpointing potential avenues for enhancing inhibitor efficacy. Taken together, this research provides invaluable insights for drug design and deepens our understanding of the binding mechanisms of PLK1.

2. Result:

2.1. Hot-spot Analysis for the Kinase Domain in PLK1

In this first stage of this work, we investigated hot spots of PLK1 using FMO analysis. We focused on two major druggable binding regions of PLK1: the kinase- and PBD-binding sites (Figure 1(A)). The complexes in each domain were derived from the protein data bank (PDB). To investigate the KD in PLK1, we compared the binding pose for known inhibitors (Figure 2(A)). Two classes of binding pose tendency were observed among the inhibitors. The A group was expanded to the solvent exposure site with the para position of the benzene group, whereas the B group was expanded to the solvent exposure site with the meta position of the benzene group, which can directly interact with Glu140. The conserved core moieties among the inhibitors are illustrated in Figure 2(B). Thus, we used three crystal structures of the most potent inhibitors, including ATP (Figure 2(A)), BI2536 (Figure 2(B), left), and Onvansertib binding structures (Figure 2(B), right).

We also analyzed the binding structures of known compounds, which are not revealed as the crystal structure, GSK461364. We performed molecular docking and MD simulations to predict the noncrystal structure binding pose. We considered any interaction with an absolute ΔE^int ≥3.0 kcal/mol significant and treated the residues involved in this interaction as hot spots. In addition, ΔE^int is represented as a PIE plot, and decompositions of intermolecular interactions for ΔPIE, which is divided into five energy terms: electrostatic (∆E^es), exchange-repulsion (∆E^ex), charge transfer with a higher-order mixed term (∆E^ct+mix), dispersion (∆E^di), and solvation energy (∆G^sol). Furthermore, we implemented an advanced analysis technique involving an MD simulation and solvation analysis to deepen our understanding of the binding mode of KD inhibitors. This approach was pursued while considering the dynamic behavior of the protein-ligand complex and the concomitant movement of water molecules. This advanced methodological approach offers a comprehensive perspective on the interplay between the inhibitor and the KD, underlining the significance of water molecules in influencing this interaction.

2.1.1. Fragment molecular orbital analysis for kinase domain: ATP–binding site

To investigate the KD of PLK1, first, we analyzed the ATP-binding structure using the FMO method. The crystal structure of the human PLK1 KD with ATP (PDB ID: 2OU7) was subjected to an FMO analysis (Figure 2(A)). The ATP-binding site is classified into the adenine, ribose, and phosphate pockets and solvent channels (Figure 3(A)) [46]. The FMO result of the holo ATP-binding structure is presented in Figure 4 and Table S1. In the adenine pocket, the adenine group interacts with the backbone of the Cys133, which has -10.425 kcal/mol PIE (Table S1). Moreover, Phe183 can interact with aromatic rings from adenine moiety and has significant binding energy. Although there is no direct interaction during visual inspection in the ribose pocket, Gly62, Arg136, and Cys67 contribute essential energy to the ribose part. Lys82 has the lowest PIE in the phosphate pocket with ATP, which binds with the phosphate part. Many groups have already demonstrated that when Lys82 mutates to Met, the kinase activity loses function, and there is no effect on S-phase progression [47]. These results suggest specific interactions between the ATP and KD of PLK1, with specific amino acids playing critical roles in stabilizing the binding.

2.1.2. Fragment molecular orbital analysis for kinase domain: BI2536 and Onvansertib

We investigated the x-ray crystal structures of known inhibitors with nanomolar potency: BI2536 has an efficacy of IC₅₀ 0.83 nM, and Onvansertib has an efficacy of IC₅₀ 2 nM (PDB ID: 2RKU, 2YAC). Although both compounds exhibit a similar nanomolar range potency, only Onvansertib demonstrates selectivity for PLK1 over PLK2 and PLK3[48]. Additionally, BI2536 and Onvansertib display markedly different binding modes, particularly in the solvent channel region (Figure 3(B)). Additionally, BI2536 occupies the solvent Channel 1 position with N-methylpiperidine, whereas Onvansertib occupies the solvent Channel 2 position with N-methylpiperazine (Figure 2(B)). The results of the FMO analyses of BI2536 and Onvansertib with PLK1 are depicted in Figure 5(A, B). Further, BI2536 exhibits significant interaction energy with 10 amino acids of PLK1: Arg57, Leu59, Gly60, Gly62, Cys67, Glu69, Leu132, Cys133, Arg136, and Phe183 (Table S2). Onvansertib also displays significant interaction energy with 14 amino acids in PLK1: Arg57, Leu59, Gly60, Cys67, Lys82, Leu132, Cys133, Arg134, Arg136, Ser137, Glu140, Asn181, Phe183, and Asp194 (Table S3). These two inhibitors share common core interactions with PLK1 involving Arg57, Leu59, Gly60, Cys67, Leu132, Cys133, Arg136, and Phe183. Particularly, Cys133 is a key residue for hinge-binding interactions in the backbone atoms forming two hydrogen bonds with the inhibitors. Onvansertib features unique interactions, distinguishable from BI2536, through the N-methylpiperazine moiety with robust electrostatic interactions with Glu140 in PLK1 and the amide moiety, forming strong hydrogen bonds with Lys82 and Asp194 in PLK1.

2.1.3. Fragment molecular orbital analysis for kinase domain: GSK461364

Like Onvansertib, GSK461364 demonstrates the same hinge interactions in the benzimidazole moiety with Cys133 in PLK1, significantly strong electrostatic interactions in the N-methylpiperazine moiety with Glu140 in PLK1, and strong hydrogen bonds in the amide moiety with Lys82 and Asp194 (Figure 6(A)). Although the drug structures of GSK and Onvansertib differ, the molecular binding tendency toward PLK1, through the FMO analysis, was consistent in both cases (Figure 6 and Table S4). Onvansertib and the GSK461364 drug possess selectivity toward PLK1, among other PLK series[49]. Compared to drugs that lack selectivity for PLK1, such as BI2536, these two drugs have been confirmed to bind to Glu140 of PLK1 specifically. Similar to Onvansertib, GSK461364 demonstrates the same hinge interactions with Cys133 in PLK1 via the benzimidazole moiety, along with significantly strong electrostatic interactions with Glu140 in PLK1 via the N-methylpiperazine moiety. A notable hydrogen bond interaction with Lys82 and Asp194 in PLK1 via the amide moiety also exists (Figure 6), implying that, despite the distinct structural differences between the GSK and Onvansertib drugs, their molecular binding tendencies toward PLK1 are consistent, as deduced from the FMO analysis. Onvansertib and GSK461364 exhibit selectivity toward PLK1 over other PLK series [49]. The superposition of the crystal structures of PLK1–3 and the sequence alignment are presented in Figure 7. The N-methylpiperazine moiety might play a vital role in imparting selectivity to the compound toward PLK1 in contrast to PLK2–3, as it establishes a polar interaction with the side chain of Glu140. Nevertheless, in PLK2 and PLK3, where Glu140 is substituted with histidine, the formation of the polar interaction is impeded (Figure 7).

2.1.4. Molecular dynamic simulation and solvation analysis for the kinase domain

We conducted a comprehensive MD simulation analysis of the holo-form of the KD by performing 12 crystal structure determinations (Table S7) and a docking structure of GSK461364, each for a duration of 100 ns. This approach ensured that the initial input docking structure did not bias the results. Subsequently, we analyzed the trajectories of each complex from the simulations by assessing ligand-protein contacts. The percentage in the findings represents the frequency of interactions occurring during the simulation. We depicted the summary of these results in a heatmap, as illustrated in Figure 8. Remarkably, all 13 ligands displayed a consistently high percentage of interaction with the hinge residue Cys133. In the phosphate pocket (Figure 3(A)), ligands from Group A interacted with Lys82 and Asp194 via a water bridge. In contrast, Group B exhibited direct and robust interactions with the residues in the phosphate pocket. Only ligands from Group B demonstrated an interaction with Glu140, a selective residue among PLK1, PLK2, and PLK 3. Based on the patterns observed in the heatmap, GSK461364 likely belongs to Group B.

Using WaterMap calculations, we performed a solvation analysis that effectively identified specific hydration sites localized near the ATP-binding site, providing valuable insight into their thermodynamic energy profiles [50]. Hydration sites with a predicted free energy (ΔG) greater than 0 kcal/mol upon water addition were classified as unfavorable sites, whereas hydration sites with a predicted free energy ΔG of less than 0 kcal/mol were considered favorable sites [51]. A comprehensive thermodynamic analysis of the hydration sites has significant importance in drug design because removing high-energy hydration sites from the protein binding site is a critical factor influencing binding affinity [52]. Comparing the WaterMap results for BI2536 and Onvansertib, we observed an unfavorable hydration site at the carbonyl position of Onvansertib, forming a hydrogen bond with Lys82 in the phosphate pocket (Figure 9(A)), which was absent in the WaterMap analysis of Onvansertib (Figure 9(C)). Furthermore, two additional unfavorable hydration sites were localized at the N-methylpiperazine position, which did not coincide with the corresponding positions in the WaterMap of Onvansertib (Figure 9(B, D)). These results suggest that Onvansertib obtains selectivity by functional groups by replacing high-energy hydration sites found in the WaterMap of BI2536.

2.2. Hot-spot Analysis for the Polo-Box Domain in PLK1

In the subsequent phase of the investigation, attention was primarily focused on the PBD of PLK1. Through the application of FMO calculations, we scrutinized its PPIs with the substrate peptide of polo-box interaction protein 1 (PBIP1), two peptide mimic inhibitors (4j, 4a), and a small molecule (KBJK557). Primarily, we evaluated the PPI hot-spot region of PBD and the substrate peptide of PBIP1 in depth. Subsequently, we performed the FMO analysis on the most potent peptide ligands and small-molecule inhibitors. Furthermore, we also performed the MD simulation and solvation analysis to gain insight into the binding mode of PBD ligands, considering the movement of the protein-ligand complex and water molecules. As a result of the comprehensive investigation, we acquired crucial insight suggesting that enhancing inhibitor potency could be effectively achieved by focusing on specific key points discovered in this study.

2.2.1. Fragment molecular orbital analysis for the polo-box domain: PBIP1-PBD

As mentioned, the PBD substrate binding pocket can be divided into phosphate, pyrrolidine, and Tyr-rich pockets [53] (Figure 10). The crystal structure (PDB ID: 3P37) of the PBD complex with FDPPLHSpTA phosphopeptide from PBIP1[54] was calculated using the FMO method to investigate the hot spot of PBIP1-PBD interactions. The FMO results of the substrate peptide interaction with the PBD are presented in graphs (Figure 11). The substrate peptide was divided into nine fragments: Ala1, p-Thr2, Ser3, His4, Leu5, Pro6, Pro7, and Phe9 (Figure 11(A)). These nine fragments have noticeable contact with the PBD at 12 residues, with PIE values lower than -3 kcal/mol. The p-Thr2 primarily interacts with Lys540 and His538 at -71.283 kcal/mol and -67.246 kcal/mol, respectively. The backbone of Ser3 interacts with the backbone of Trp414 from the PBD through two hydrogen bonds (Figure 11(B)). These results agree with previous studies that the immunoprecipitation of wild-type PBIP1 strongly interacts with the PBD, but the active signal disappears with one of the serine77 (Ser3) and threonine78 (p-Thr2) residues [55]. Further, Ala1, p-Thr2, and Ser3 occupy phosphate and pyrrolidine pockets, which maintain high energy with Leu490, Leu491, His538, Lys540, and Trp414 (Figure 11(B), right table). The remaining seven residues contribute interaction energy to the Tyr-rich pocket, displayed as a light-pink surface (Figure 10) and a light-pink stick (Figure 11(B)). In particular, Tyr417 from the Tyr-rich pocket interacts simultaneously with Leu5, Pro6, Pro7, Asp8, and Phe9 from the substrate peptide. Sharma et al. demonstrated that Tyr-rich pocket mutation selectivity impairs PBD binding to the substrate protein PBIP1[56]. The mutagenesis data are related to the FMO results in this study.

2.2.2. Fragment molecular orbital analysis for the polo-box domain: peptide ligands

Ligand 4j is the most potent peptide mimetic ligand developed from the PLHSpT sequence and has an alkyl-phenyl group with a PBD-binding IC₅₀ value of 0.12 μM. However, the PBD binding with a 4a IC₅₀ value is more than 200 μM, as the alkyl-phenyl group is absent [57]. We calculated PIE values with the crystal structure of the 4j binding in the PBD and the docking model 4a in the PBD using the FMO method. Two peptide ligands were fragmented based on peptide bonds (Figure 12(A) and Figure 13(A)). In the phosphate and pyrrolidine pocket, the p-Thr2 and Ser3 fragments from 4j have more contribution energy residues than 4a (Figure 12 and Figure 13). In particular, the PIE between 4j with His538 and Lys540 is much lower than 4a (Figure 12 and Figure 13). We found more significant energy differences between the two peptide ligands in the Tyr-rich pocket. Seven essential key residues have significant PIE values with F-Akyl from 4j, but only four with F-CH3 from 4a (Figure S1). Moreover, Tyr417 has the most frequent interactions with substrate peptides and has significant interaction energy with 4j but not with 4a. Thus, the FMO results of the two peptides are related to those from the structure analysis relationship mentioned above [57].

2.2.3. Fragment molecular orbital analysis for the polo-box domain: small molecules

In addition, KBJK557 was identified as a potent PBD-binding small molecule, confirming that a four-carbon alkyl chain attached to pyrazole nitrogen might increase the potency through structure-activity relationship (SAR) studies [58]. For negative control, we also chose KBJK-4a, which does not have the four-carbon alkyl chain, to perform the in silico analysis. First, we performed the docking generation of KBJK557 in the PBD to determine the best pose with the lowest Glide Emodel score suitable for pose selection. In this best pose, the barbituric acid interaction with Lys540, His538, and Trp414 formed hydrogen bonds and a salt bridge, and the pyrazole part exhibited π-π stacking interaction with Trp414 (Figure S2(A)). Furthermore, the phenyl group in the four-carbon alkyl tail engaged in a π-π stacking interaction with Tyr417 (Figure S2(A)). Second, we ran the 50ns MD simulation of the docking structure to optimize it. We chose the optimized structure from Frame 489 of the MD simulation trajectory for the next step (Figure S2(B) and Figure 14(A)). From the analysis of the MD simulation, the interactions found with the docking structure continued during the simulation (Figure S2(C, D)). We also docked the KBJK-4a in the protein structure from Frame 489, and those interactions in KBJK557 almost exist except for the π-π stacking interaction with Tyr417, particularly (Figure 14(D)).

Finally, we performed the FMO study for two small-molecule complex structures (Figure 14). The KBJK557 interaction residues in the PBD were His538, Lys540, Trp414, Asp416, Tyr485, Tyr417, Leu490, Val415, Leu491, and Phe535 (Figure 14(A) and Table S5). Some of these residues have been mentioned in previous hot-spot analyses for PBIP1-PBD PPIs and peptide ligands. Moreover, His538 and Lys540 interacted most with KBJK557, with PIE values of -25 kcal/mol and -23.556 kcal/mol. Electrostatic terms primarily drove these interactions. The Trp414 contacted KBJK557 with a PIE of -21.232 kcal/mol, driven by electrostatic and dispersion terms. The KBJK557 also engaged in hydrophobic interactions in the Tyr-rich pocket with Val415, Tyr417, and Tyr485 (Figure 14(B, C)). However, the interaction with Tyr417 disappeared in KBJK-4a (Figure 14(E, F) and Table S6). These results are related to the previously mentioned analysis of peptide ligands 4j and 4a (Figure 12, Figure 13, and Figure S1). Moreover, these results align with previous SAR studies [58].

To gain insight into increasing the inhibition potency of the PBD in PLK1, we summarized all FMO results into one heatmap (Figure 15). We listed the amino acids of the PBD selected within 5 Å of the ligand and highlighted the critical residues in each pocket. The PIE was summed for each residue, especially for the substrate peptide and peptide mimetic ligands. Exception KBJK-4a, no dramatic energy changes occurred between the interaction with His538, Lys540, and Trp414 from the phosphate and pyrrolidine pockets (Figure 15). However, in the Tyr-rich pocket, peptide ligand 4a, and small molecule KBJK-4a, which have no or low potency in the experiment [58], the PIE became unstable with Tyr417, Tyr481, Phe482, and Tyr485 (Figure 15). The PIE of KBJK557 with Tyr481 and Phe482 was also reduced because the phenyl group in the four-carbon alkyl chain was short of occupying the entire Tyr-rich pocket (Figure S3). Taken together, the heatmap of the FMO results of five structures illustrates that the Tyr-rich pocket, in particular, Tyr417, Tyr481, Phe482, and Tyr485, may include crucial residues to contribute to the binding affinity of the PBD. These results relate to the mutagenesis experiments from the Sharma groups mentioned above [56].

2.2.4. Molecular dynamics simulation and solvation analysis for the polo-box domain

We assessed complex structures, including substrate peptide PBIP1, two peptide ligands, and KBJK557 of the PBD, for 100 ns. In KBJK557, we used the selected snapshot structure of KBJK557 as an initial structure for MD simulation. Then, the trajectory of each complex for each run was analyzed using ligand-protein contacts. The percentage represents interactions that occur during the simulation time. The heatmap of the protein-ligand contact from the MD simulation is presented in Figure 16, and ligands represent substrate peptides, peptide ligands, and small molecules that bind in the PBD. The Trp414 from the pyrrolidine pocket maintains a high percentage with all ligands, which were also displaced in the FMO analysis (Figure 15 and Figure 16). This outcome indicates that Trp414 is an essential residue in the PBD-binding site, which is in the central position. Peptide inhibitor 4j has a contact map similar to the substrate peptide, but 4a loses the interaction contact in the Tyr-rich pocket. Small molecules, KBJK557 and KBJK-4a, demonstrate mediate interaction contact with His538 and Lys540 compared to the substrate and peptide ligands—residues belonging to the phosphate-binding pocket (Figure 16), which might be why the potency of small molecules is lower than that of the peptide ligands.

In the solvation analysis for the PBD, we might discover why the potency is different between 4j and 4a. In the WaterMap of 4a, three unstable hydration sites occupy the al-kyl-phenyl moiety position of 4j, forming a π-π stacking interaction with Phe482 and Tyr417 (Figure 17(A, B)). However, these hydration sites are not exhibited in the WaterMap of 4j (Figure 17(C)). Four unfavorable solvation sites were also found at the same position in the WaterMap of KBJK557, around the Tyr-rich pocket (Figure 17(D)). These results suggest that efforts to improve the inhibition potency of the PBD should focus on extending the moieties into the Tyr-rich pocket.

3. Discussion:

A group of PLKs, consisting of five serine/threonine kinases, can be found in different eukaryotic organisms[59]. These kinases play a crucial role in regulating cell proliferation, particularly in controlling the progression of the cell cycle [2]. The PLK1 protein consists of an N-terminal serine/threonine kinase domain (KD) and a C-terminal repeat of the Polo-box domain (PBD), with the phosphorylation of the latter directly affecting the enzymatic activity of PLK1 [53]. The PBD recognizes phosphorylated serine/threonine protein substrates to regulate PLK KD's phosphorylation activity [20]. Each domain comprises druggable binding sites, including an ATP-binding site in the KD and a substrate peptide binding site in the PBD.

Numerous ATP-competitive inhibitors have been developed, and some have even advanced to clinical trials (as indicated in Table 1). However, these ATP-competitive inhibitors have shown limitations in antitumor activity during clinical trials, often exhibiting dose-dependent side effects. Consequently, a recent focus has been on inhibitors that bind to the PBD instead [46]. In this study, we focus on analyzing the ATP binding site and investigating the binding site in PBD to design an improved inhibitor. We analyzed using the FMO method, MD simulations, and solvation analysis.

In the former part of our study, we analyzed the ATP-binding site in the kinase domain. It is widely recognized that the adenine group of ATP interacts with the backbone of the hinge residue in kinase inhibitors [60]. Our FMO study shows that PLK1 inhibitors strongly interact with Cys133 via their core scaffold, as reflected in the high absolute PIE value. These findings underscore the pivotal role of hinge interactions in achieving the high potency of ATP-competitive inhibitors. However, inhibitors of PLK1 that are ATP-competitive have exhibited a consistent problem with selectivity towards PLK2 and PLK3. To address this challenge, we have highlighted the significance of Glu140 in PLK1. In contrast, the corresponding positions in PLK2 and PLK3 are His169 and His149, respectively (Figure 7). This insight arose from the analysis of two well-known inhibitors, BI2536 and Onvansertib, both of which exhibited nanomolar potency against PLK1. Notably, only Onvansertib demonstrated selectivity for PLK1 over PLK2 and PLK3. While both inhibitors share common core interactions, Onvansertib exhibits additional interactions facilitated by specific moieties, such as N-methylpiperidine, contributing to its unique selectivity and binding properties. The N-methylpiperazine moiety in Onvansertib plays a crucial role in conferring selectivity towards PLK1, primarily due to its polar interaction with the side chain of Glu140. This insight deepens our understanding of the SAR of these inhibitors and offers guidance for the development of more potent and selective ATP-competitive PLK1 inhibitors in the future.

In the latter part of our study, we delved into identifying the binding site for the substrate peptide within the PBD. Over the years, a range of ATP-competitive inhibitors have been crafted, with several advancing to clinical trial stages. However, their effectiveness has often been overshadowed by dose-dependent adverse reactions. As a result, there has been a growing interest in inhibitors targeting the PBD. Yet, the design of nanomolar small molecule inhibitors specific to the PBD remains a formidable task. To navigate this, we probed the PBD substrate binding pocket, aiming to elucidate strategies to enhance inhibitor binding potency. Our FMO analyses reasonably pinpointed the Tyr-rich pocket as an essential component in enhancing PBD inhibitor efficacy. In our quest, we evaluated both peptide ligands and small molecules for which SAR data was accessible. Our focus was on contrasting ligands equipped with a binding moiety for the Tyr-rich pocket against those devoid of it. Notably, ligands lacking this specific moiety yielded considerably diminished absolute PIE values in FMO calculations or manifested heightened energy solvation sites. This underscores the potential of tailoring ligand moieties to fit the Tyr-rich pocket better, positing it as a potent strategy to optimize PBD inhibition. Our discoveries offer profound insights, enriching our understanding of binding dynamics and paving the way to design more potent PBD inhibitors.

Collectively, our results compellingly advocate for targeting both KD and PBD pockets with bivalent inhibitors as a potent approach to inhibit PLK1. Bivalent kinase inhibitors, characterized as a novel subset of small molecule compounds, are designed to engage two binding sites on kinase enzymes concurrently. These inhibitors stand out, offering a distinct edge over traditional inhibitors by simultaneously targeting both the ATP-binding (active) site and another site, typically termed as an "allosteric" or "regulatory" site [61,62]. Paralley, recent studies from other research groups have showcased the efficacy of bivalent inhibitors that bridge KD and PBD inhibitors for PLK1 [34,63]. Such groundbreaking strategies, epitomized by bivalent inhibitors, hold the potential to overcome the clinical limitations associated with ATP-competitive inhibitors of PLK1. The conception and development of these bivalent inhibitors demand a synergistic fusion of structural understanding, advanced computational modeling, and nuanced medicinal chemistry. Anchored by the findings of our study, we believe we stand ready to provide indispensable insights for the next wave of bivalent PLK1 inhibitor designs.

4. Conclusion

In conclusion, our study integrated the FMO methodology, MD simulations, and solvation analyses to investigate the KD and PBD domains of the PLK1 protein. The insights shed light on the intricate interactions and binding preferences of various PLK1 inhibitors. Our findings underscore the consistent binding patterns of ATP-competitive inhibitors to PLK1, emphasizing their selectivity. Furthermore, our investigation into the PBD's substrate peptide binding site pinpointed crucial residues vital for interactions, paving the way to bolster inhibitor potency. Collectively, these revelations considerably elevate our comprehension of the mechanisms dictating binding, selectivity, and the prospective therapeutic avenues for PLK1-targeting inhibitors. Grounded in our research, we are poised to contribute invaluable perspectives for the next generation of bivalent PLK1 inhibitor designs. This wealth of knowledge will indubitably influence the blueprint of upcoming pharmaceuticals, aiming to craft more efficacious and discerning inhibitors that modulate cell growth and orchestrate the cell cycle.

5. Method

5.1. Protein Structure Preparation

Fourteen crystal structures were retrieved from the PDB (Table S7). The structure was prepared using the Protein Preparation Wizard tool. All missing loops were filled using Prime implemented in the Maestro program (v. 2022; Schrödinger, LLC, New York, NY, 2022). Hydrogen atoms were added to the complex structure at a pH of 7.4, and their positions were optimized using PROPKA implemented in the Maestro program. Restrained energy minimization was performed with an OPLS4 force field within a 0.3 Å root mean square deviation.

5.2. Molecular Docking

Molecular docking was conducted to generate the docking pose of GSK461364, KBJK557, and KBJK-4a using Glide [64] (Schrödinger, LLC, New York, NY, 2022). In addition, GSK461364 was docked in the KD of PLK1, and KBJK557 and KBJK-4a were docked in the PBD of PLK1.

5.3. Fragment Molecular Orbital Calculations

All FMO computations were computed using the GAMESS program [65]. The two-body FMO method was employed in all such calculations, with each residue in the protein and ligand being classified as a fragment. Per the hybrid orbital projection scheme fragmentation [66], an in-house code was used to prepare all input files.

The two-body FMO method comprises four stages: (1) fragmentation, (2) single-fragment self-consistent field (SCF) calculation, (3) two fragment SCF computation, and (4) total property evaluation, the specifics of which have been previously detailed [67]. During the fragmentation phase (1), each residue in the protein, ligand, and water molecules can be designated as a fragment. Unlike regular peptide bond separation in proteins, all residues were divided at the sp3 bond between the alpha carbon and carbonyl carbon atoms in the backbone structure per the hybrid orbital projection scheme [66]. This method significantly decreases computational expenses and corrects inaccuracies from artificial fragmentation through a projection operator [66].

In the second and third stages (2, 3), all molecular orbitals (MOs) on a fragment were optimized using the SCF theory in the overall electrostatic field, and all electron densities were self-consistently resolved through self-consistent charge iterations [36,67]. The difference between these stages lies in the Hamiltonian operators [36]. The second stage optimized all MOs in a fragment, including the electrostatic potential from N-1 fragments, whereas the third stage optimized all MOs in two fragments and included the potential from N-2 fragments.

All MO results obtained in Stages 2 and 3 were combined to assess the system’s total properties (4) [67]. The PIE values between two fragments were computed, and the energy decomposition of the PIE values was conducted to comprehend the contributions of the five energy terms. These PIE values between fragments in the FMO computations were divided into five energy terms, as defined in Equation 1: electrostatic (∆E^es), exchange-repulsion (∆E^ex), charge transfer with a higher order mixed term (∆E^ct+mix), dispersion (∆E^di), and solvation energy (∆G^sol) derived from the polarizable continuum model (PCM):

∆E^int = ∆E^es + ∆E^ex + ∆E^ct+mix + ∆E^di + ∆G^sol

(1)

To scrutinize the crucial interactions of several protein-peptide complex or protein-ligand structures, we applied the second-order Møller–Plesset perturbation theory (MP2) [68] and PCM [69,70] with a 6–31G** basis set (FMO2-MP2/6–31G**/PCM level) to the protein-ligand or protein-peptide complex. We also calculated the crystal structure of the PLK1 complex with ATP using the FMO2/DFTB3 method [71]. Subsequently, we carried out significant PPIs in the complex by selecting PIE values that were more stable than -3.0 kcal/mol and had a distance of less than 5.4 Å between two fragments { Bodkin.2016}.

5.4. Molecular Dynamics Simulation

The MD simulation was performed using Desmond in Maestro (v. 2022-4; Schrödinger, LLC, New York, NY, 2022). The 14 protein-ligand crystal complexes (Table S7) and docking structure of GSK461364, 4a, KBJK557, and KBJK-4a were inserted into an orthorhombic box filled with explicit water molecules (TIP3P model) and with a buffer distance of 10 Å. The MD simulation was studied, and the OPLS4 force field was used. Ions (Na⁺ and Cl^-) were added to simulate a physiological concentration of monovalent ions (0.15 M). A constant number of particles at a constant temperature of 300K with a pressure of 1.01325 bar was used as an ensemble class. The particle-mesh Ewald method [72] was used to calculate long-range electrostatic interactions, with a cutoff for van der Waals and short-range electrostatic interactions of 9 Å. Nose-Hoover thermostats [73] were employed to maintain a constant simulation temperature, whereas the Martina–Tobias–Klein [74] method controlled the pressure. A RESPA integrator [75] integrated the equation of motion with an inner time step of 2.0 fs for bonded and non-bonded interactions within the short-range cutoff. The default protocol in Desmond was applied to reach system equilibration, and a 100ns simulation was performed for each complex [76]. The analysis plots and Figures were sketched using the Desmond simulation interaction diagram panel of the Maestro program.

5.5. WaterMap Calculations

WaterMap, which has been extensively explained in other sources, employs a molecular dynamics methodology to forecast the thermodynamic characteristics of water molecules within a protein setting, encompassing entropy, enthalpy, and free energy for each water site [77]. The calculations for WaterMap (v. 2018-4; Schrödinger, LLC, New York, NY, 2018) were conducted using default parameters, involving a 2ns simulation time, with water molecules analyzed within a 10Å radius of the ligand and using the prepared cocrystal structures (Table S7) and docking structures mentioned in this paper with each ligand present.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org.

References

Raab, C.A.; Raab, M.; Becker, S.; Strebhardt, K. Non-Mitotic Functions of Polo-like Kinases in Cancer Cells. Biochimica et Biophysica Acta (BBA) - Reviews on Cancer 2021, 1875, 188467. [Google Scholar] [CrossRef]
Barr, F.A.; Silljé, H.H.W.; Nigg, E.A. Polo-like Kinases and the Orchestration of Cell Division. Nature Reviews Molecular Cell Biology 2004, 5, 429–441. [Google Scholar] [CrossRef] [PubMed]
De Cárcer, G.; Manning, G.; Malumbres, M. From Plk1 to Plk5. Cell Cycle 2011, 10, 2255–2262. [Google Scholar] [CrossRef]
Clay, F.J.; McEwen, S.J.; Bertoncello, I.; Wilks, A.F.; Dunn, A.R. Identification and Cloning of a Protein Kinase-Encoding Mouse Gene, Plk, Related to the Polo Gene of Drosophila. Proc. Natl. Acad. Sci. 1993, 90, 4882–4886. [Google Scholar] [CrossRef] [PubMed]
Petronczki, M.; Lénárt, P.; Peters, J.-M. Polo on the Rise—From Mitotic Entry to Cytokinesis with Plk1. Dev. Cell 2008, 14, 646–659. [Google Scholar] [CrossRef] [PubMed]
Liu, Z.; Sun, Q.; Wang, X. PLK1, A Potential Target for Cancer Therapy. Translational Oncology 2017, 10, 22–32. [Google Scholar] [CrossRef] [PubMed]
Wang, L.; Guo, Q.; Fisher, L.A.; Liu, D.; Peng, A. Regulation of Polo-like Kinase 1 by DNA Damage and PP2A/B55α. Cell Cycle 2015, 14, 157–166. [Google Scholar] [CrossRef] [PubMed]
Gheghiani, L.; Loew, D.; Lombard, B.; Mansfeld, J.; Gavet, O. PLK1 Activation in Late G2 Sets Up Commitment to Mitosis. Cell Rep. 2017, 19, 2060–2073. [Google Scholar] [CrossRef]
Lee, K.; Rhee, K. PLK1 Phosphorylation of Pericentrin Initiates Centrosome Maturation at the Onset of Mitosis. J. Cell Biol. 2011, 195, 1093–1101. [Google Scholar] [CrossRef]
Kim, J.H.; Shim, J.; Ji, M.-J.; Jung, Y.; Bong, S.M.; Jang, Y.-J.; Yoon, E.-K.; Lee, S.-J.; Kim, K.G.; Kim, Y.H.; et al. The Condensin Component NCAPG2 Regulates Microtubule–Kinetochore Attachment through Recruitment of Polo-like Kinase 1 to Kinetochores. Nat. Commun. 2014, 5, 4588. [Google Scholar] [CrossRef]
Cunningham, C.E.; MacAuley, M.J.; Vizeacoumar, F.S.; Abuhussein, O.; Freywald, A.; Vizeacoumar, F.J. The CINs of Polo-Like Kinase 1 in Cancer. Cancers 2020, 12, 2953. [Google Scholar] [CrossRef]
Yu, Z.; Deng, P.; Chen, Y.; Liu, S.; Chen, J.; Yang, Z.; Chen, J.; Fan, X.; Wang, P.; Cai, Z.; et al. Inhibition of the PLK1-Coupled Cell Cycle Machinery Overcomes Resistance to Oxaliplatin in Colorectal Cancer. Adv. Sci. 2021, 8, 2100759. [Google Scholar] [CrossRef]
Gheghiani, L.; Wang, L.; Zhang, Y.; Moore, X.T.R.; Zhang, J.; Smith, S.C.; Tian, Y.; Wang, L.; Turner, K.; Jackson-Cook, C.K.; et al. PLK1 Induces Chromosomal Instability and Overrides Cell-Cycle Checkpoints to Drive Tumorigenesis. Cancer Res. 2021, 81, 1293–1307. [Google Scholar] [CrossRef] [PubMed]
Ma, X.; Wang, L.; Huang, D.; Li, Y.; Yang, D.; Li, T.; Li, F.; Sun, L.; Wei, H.; He, K.; et al. Polo-like Kinase 1 Coordinates Biosynthesis during Cell Cycle Progression by Directly Activating Pentose Phosphate Pathway. Nat. Commun. 2017, 8, 1506. [Google Scholar] [CrossRef]
Spänkuch-Schmitt, B.; Wolf, G.; Solbach, C.; Loibl, S.; Knecht, R.; Stegmüller, M.; von Minckwitz, G.; Kaufmann, M.; Strebhardt, K. Downregulation of Human Polo-like Kinase Activity by Antisense Oligonucleotides Induces Growth Inhibition in Cancer Cells. Oncogene 2002, 21, 3162–3171. [Google Scholar] [CrossRef]
Strebhardt, K.; Ullrich, A. Targeting Polo-like Kinase 1 for Cancer Therapy. Nat. Rev. Cancer 2006, 6, 321–330. [Google Scholar] [CrossRef]
Taniguchi, E.; Toyoshima-Morimoto, F.; Nishida, E. Nuclear Translocation of Plk1 Mediated by Its Bipartite Nuclear Localization Signal*. J. Biol. Chem. 2002, 277, 48884–48888. [Google Scholar] [CrossRef] [PubMed]
Bruinsma, W.; Macůrek, L.; Freire, R.; Lindqvist, A.; Medema, R.H. Bora and Aurora-A Continue to Activate Plk1 in Mitosis. J. Cell Sci. 2013, 127, 801–811. [Google Scholar] [CrossRef]
Xu, J.; Shen, C.; Wang, T.; Quan, J. Structural Basis for the Inhibition of Polo-like Kinase 1. Nat. Struct. Mol. Biol. 2013, 20, 1047–1053. [Google Scholar] [CrossRef] [PubMed]
Elia, A.E.H.; Rellos, P.; Haire, L.F.; Chao, J.W.; Ivins, F.J.; Hoepker, K.; Mohammad, D.; Cantley, L.C.; Smerdon, S.J.; Yaffe, M.B. The Molecular Basis for Phosphodependent Substrate Targeting and Regulation of Plks by the Polo-Box Domain. Cell 2003, 115, 83–95. [Google Scholar] [CrossRef]
Zitouni, S.; Nabais, C.; Jana, S.C.; Guerrero, A.; Bettencourt-Dias, M. Polo-like Kinases: Structural Variations Lead to Multiple Functions. Nature Reviews Molecular Cell Biology 2014, 15, 433–452. [Google Scholar] [CrossRef] [PubMed]
Cheng, K.; Lowe, E.D.; Sinclair, J.; Nigg, E.A.; Johnson, L.N. The Crystal Structure of the Human Polo-like Kinase-1 Polo Box Domain and Its Phospho-peptide Complex. EMBO J. 2003, 22, 5757–5768. [Google Scholar] [CrossRef] [PubMed]
Liao, C.; Yao, R. Diversity Evolution and Jump of Polo-like Kinase 1 Inhibitors. Sci. China Chem. 2013, 56, 1392–1401. [Google Scholar] [CrossRef]
Zeidan, A.M.; Ridinger, M.; Lin, T.L.; Becker, P.S.; Schiller, G.J.; Patel, P.A.; Spira, A.I.; Tsai, M.L.; Samuëlsz, E.; Silberman, S.L.; et al. A Phase Ib Study of Onvansertib, a Novel Oral PLK1 Inhibitor, in Combination Therapy for Patients with Relapsed or Refractory Acute Myeloid Leukemia. Clinical Cancer Research 2020, 26, 6132–6140. [Google Scholar] [CrossRef]
Gutteridge, R.E.A.; Ndiaye, M.A.; Liu, X.; Ahmad, N. Plk1 Inhibitors in Cancer Therapy: From Laboratory to Clinics. Molecular Cancer Therapeutics 2016, 15, 1427–1435. [Google Scholar] [CrossRef]
Awad, M.M.; Chu, Q.S.-C.; Gandhi, L.; Stephenson, J.J.; Govindan, R.; Bradford, D.S.; Bonomi, P.D.; Ellison, D.M.; Eaton, K.D.; Fritsch, H.; et al. An Open-Label, Phase II Study of the Polo-like Kinase-1 (Plk-1) Inhibitor, BI 2536, in Patients with Relapsed Small Cell Lung Cancer (SCLC). Lung Cancer 2017, 104, 126–130. [Google Scholar] [CrossRef] [PubMed]
Haupenthal, J.; Bihrer, V.; Korkusuz, H.; Kollmar, O.; Schmithals, C.; Kriener, S.; Engels, K.; Pleli, T.; Benz, A.; Canamero, M.; et al. Reduced Efficacy of the Plk1 Inhibitor BI 2536 on the Progression of Hepatocellular Carcinoma Due to Low Intratumoral Drug Levels. Neoplasia 2012, 14, 410–IN10. [Google Scholar] [CrossRef] [PubMed]
Gorlick, R.; Kolb, E.A.; Keir, S.T.; Maris, J.M.; Reynolds, C.P.; Kang, M.H.; Carol, H.; Lock, R.; Billups, C.A.; Kurmasheva, R.T.; et al. Initial Testing (Stage 1) of the Polo-like Kinase Inhibitor Volasertib (BI 6727), by the Pediatric Preclinical Testing Program. Pediatr. Blood Cancer 2014, 61, 158–164. [Google Scholar] [CrossRef]
Sero, V.; Tavanti, E.; Vella, S.; Hattinger, C.M.; Fanelli, M.; Michelacci, F.; Versteeg, R.; Valsasina, B.; Gudeman, B.; Picci, P.; et al. Targeting Polo-like Kinase 1 by NMS-P937 in Osteosarcoma Cell Lines Inhibits Tumor Cell Growth and Partially Overcomes Drug Resistance. Investig. N. Drugs 2014, 32, 1167–1180. [Google Scholar] [CrossRef]
Daub, H.; Specht, K.; Ullrich, A. Strategies to Overcome Resistance to Targeted Protein Kinase Inhibitors. Nat. Rev. Drug Discov. 2004, 3, 1001–1010. [Google Scholar] [CrossRef]
Liao, C.; Park, J.-E.; Bang, J.K.; Nicklaus, M.C.; Lee, K.S. Probing Binding Modes of Small Molecule Inhibitors to the Polo-Box Domain of Human Polo-like Kinase 1. ACS Med. Chem. Lett. 2010, 1, 110–114. [Google Scholar] [CrossRef]
Reindl, W.; Yuan, J.; Krämer, A.; Strebhardt, K.; Berg, T. Inhibition of Polo-like Kinase 1 by Blocking Polo-Box Domain-Dependent Protein-Protein Interactions. Chem. Biol. 2008, 15, 459–466. [Google Scholar] [CrossRef]
Scharow, A.; Raab, M.; Saxena, K.; Sreeramulu, S.; Kudlinzki, D.; Gande, S.; Dötsch, C.; Kurunci-Csacsko, E.; Klaeger, S.; Kuster, B.; et al. Optimized Plk1 PBD Inhibitors Based on Poloxin Induce Mitotic Arrest and Apoptosis in Tumor Cells. ACS Chem. Biol. 2015, 10, 2570–2579. [Google Scholar] [CrossRef] [PubMed]
Scharow, A.; Knappe, D.; Reindl, W.; Hoffmann, R.; Berg, T. Development of Bifunctional Inhibitors of Polo-Like Kinase 1 with Low-Nanomolar Activities Against the Polo-Box Domain. ChemBioChem 2016, 17, 759–767. [Google Scholar] [CrossRef]
Kitaura, K.; Ikeo, E.; Asada, T.; Nakano, T.; Uebayasi, M. Fragment Molecular Orbital Method: An Approximate Computational Method for Large Molecules. Chem. Phys. Lett. 1999, 313, 701–706. [Google Scholar] [CrossRef]
Tanaka, S.; Mochizuki, Y.; Komeiji, Y.; Okiyama, Y.; Fukuzawa, K. Electron-Correlated Fragment-Molecular-Orbital Calculations for Biomolecular and Nano Systems. Phys. Chem. Chem. Phys. 2014, 16, 10310–10344. [Google Scholar] [CrossRef]
Lim, H.; Baek, A.; Kim, J.; Kim, M.S.; Liu, J.; Nam, K.-Y.; Yoon, J.; No, K.T. Hot Spot Profiles of SARS-CoV-2 and Human ACE2 Receptor Protein Protein Interaction Obtained by Density Functional Tight Binding Fragment Molecular Orbital Method. Sci. Rep. 2020, 10, 16862. [Google Scholar] [CrossRef] [PubMed]
Tanaka, S.; Watanabe, C.; Honma, T.; Fukuzawa, K.; Ohishi, K.; Maruyama, T. Identification of Correlated Inter-Residue Interactions in Protein Complex Based on the Fragment Molecular Orbital Method. J. Mol. Graph. Model. 2020, 100, 107650. [Google Scholar] [CrossRef]
Nakano, T.; Kaminuma, T.; Sato, T.; Fukuzawa, K.; Akiyama, Y.; Uebayasi, M.; Kitaura, K. Fragment Molecular Orbital Method: Use of Approximate Electrostatic Potential. Chem. Phys. Lett. 2002, 351, 475–480. [Google Scholar] [CrossRef]
Watanabe, C.; Fukuzawa, K.; Okiyama, Y.; Tsukamoto, T.; Kato, A.; Tanaka, S.; Mochizuki, Y.; Nakano, T. Three- and Four-Body Corrected Fragment Molecular Orbital Calculations with a Novel Subdividing Fragmentation Method Applicable to Structure-Based Drug Design. J. Mol. Graph. Model. 2013, 41, 31–42. [Google Scholar] [CrossRef]
Yoshida, T.; Fujita, T.; Chuman, H. Novel Quantitative Structure-Activity Studies of HIV-1 Protease Inhibitors of the Cyclic Urea Type Using Descriptors Derived from Molecular Dynamics and Molecular Orbital Calculations. Curr. Comput. Aided-Drug Des. 2009, 5, 38–55. [Google Scholar] [CrossRef]
Hitaoka, S.; Matoba, H.; Harada, M.; Yoshida, T.; Tsuji, D.; Hirokawa, T.; Itoh, K.; Chuman, H. Correlation Analyses on Binding Affinity of Sialic Acid Analogues and Anti-Influenza Drugs with Human Neuraminidase Using Ab Initio MO Calculations on Their Complex Structures – LERE-QSAR Analysis (IV). J. Chem. Inf. Model. 2011, 51, 2706–2716. [Google Scholar] [CrossRef]
Heifetz, A.; Chudyk, E.I.; Gleave, L.; Aldeghi, M.; Cherezov, V.; Fedorov, D.G.; Biggin, P.C.; Bodkin, M.J. The Fragment Molecular Orbital Method Reveals New Insight into the Chemical Nature of GPCR–Ligand Interactions. J. Chem. Inf. Model. 2016, 56, 159–172. [Google Scholar] [CrossRef] [PubMed]
Kim, J.; Lim, H.; Moon, S.; Cho, S.Y.; Kim, M.; Park, J.H.; Park, H.W.; No, K.T. Hot Spot Analysis of YAP-TEAD Protein-Protein Interaction Using the Fragment Molecular Orbital Method and Its Application for Inhibitor Discovery. Cancers 2021, 13, 4246. [Google Scholar] [CrossRef] [PubMed]
Lim, H.; Chun, J.; Jin, X.; Kim, J.; Yoon, J.; No, K.T. Investigation of Protein-Protein Interactions and Hot Spot Region between PD-1 and PD-L1 by Fragment Molecular Orbital Method. Sci. Rep. 2019, 9, 16727. [Google Scholar] [CrossRef]
Zhang, J.; Zhang, L.; Wang, J.; Ouyang, L.; Wang, Y. Polo-like Kinase 1 Inhibitors in Human Cancer Therapy: Development and Therapeutic Potential. Journal of Medicinal Chemistry 2022, 65, 10133–10160. [Google Scholar] [CrossRef]
Shen, M.; Cai, Y.; Yang, Y.; Yan, X.; Liu, X.; Zhou, T. Centrosomal Protein FOR20 Is Essential for S-Phase Progression by Recruiting Plk1 to Centrosomes. Cell Res. 2013, 23, 1284–1295. [Google Scholar] [CrossRef]
Beria, I.; Bossi, R.T.; Brasca, M.G.; Caruso, M.; Ceccarelli, W.; Fachin, G.; Fasolini, M.; Forte, B.; Fiorentini, F.; Pesenti, E.; et al. NMS-P937, a 4,5-Dihydro-1H-Pyrazolo[4,3-h]Quinazoline Derivative as Potent and Selective Polo-like Kinase 1 Inhibitor. Bioorganic & Medicinal Chemistry Letters 2011, 21, 2969–2974. [Google Scholar] [CrossRef]
Gilmartin, A.G.; Bleam, M.R.; Richter, M.C.; Erskine, S.G.; Kruger, R.G.; Madden, L.; Hassler, D.F.; Smith, G.K.; Gontarek, R.R.; Courtney, M.P.; et al. Distinct Concentration-Dependent Effects of the Polo-like Kinase 1–Specific Inhibitor GSK461364A, Including Differential Effect on Apoptosis. Cancer Research 2009, 69, 6969–6977. [Google Scholar] [CrossRef]
Biswal, J.; Jayaprakash, P.; Rayala, S.K.; Venkatraman, G.; Rangaswamy, R.; Jeyaraman, J. WaterMap and Molecular Dynamic Simulation-Guided Discovery of Potential PAK1 Inhibitors Using Repurposing Approaches. ACS Omega 2021, 6, 26829–26845. [Google Scholar] [CrossRef]
Leger, P.R.; Hu, D.X.; Biannic, B.; Bui, M.; Han, X.; Karbarz, E.; Maung, J.; Okano, A.; Osipov, M.; Shibuya, G.M.; et al. Discovery of Potent, Selective, and Orally Bioavailable Inhibitors of USP7 with In Vivo Antitumor Activity. J. Med. Chem. 2020, 63, 5398–5420. [Google Scholar] [CrossRef]
Beuming, T.; Farid, R.; Sherman, W. High-energy Water Sites Determine Peptide Binding Affinity and Specificity of PDZ Domains. Protein Sci. 2009, 18, 1609–1619. [Google Scholar] [CrossRef]
Lee, K.S.; Burke, T.R.; Park, J.-E.; Bang, J.K.; Lee, E. Recent Advances and New Strategies in Targeting Plk1 for Anticancer Therapy. Trends in Pharmacological Sciences 2015, 36, 858–877. [Google Scholar] [CrossRef]
Śledź, P.; Stubbs, C.J.; Lang, S.; Yang, Y.; McKenzie, G.J.; Venkitaraman, A.R.; Hyvönen, M.; Abell, C. From Crystal Packing to Molecular Recognition: Prediction and Discovery of a Binding Site on the Surface of Polo-Like Kinase 1. Angewandte Chemie International Edition 2011, 50, 4003–4006. [Google Scholar] [CrossRef] [PubMed]
Kang, Y.H.; Park, J.-E.; Yu, L.-R.; Soung, N.-K.; Yun, S.-M.; Bang, J.K.; Seong, Y.-S.; Yu, H.; Garfield, S.; Veenstra, T.D.; et al. Self-Regulated Plk1 Recruitment to Kinetochores by the Plk1-PBIP1 Interaction Is Critical for Proper Chromosome Segregation. Molecular Cell 2006, 24, 409–422. [Google Scholar] [CrossRef] [PubMed]
Sharma, P.; Mahen, R.; Rossmann, M.; Stokes, J.E.; Hardwick, B.; Huggins, D.J.; Emery, A.; Kunciw, D.L.; Hyvönen, M.; Spring, D.R.; et al. A Cryptic Hydrophobic Pocket in the Polo-Box Domain of the Polo-like Kinase PLK1 Regulates Substrate Recognition and Mitotic Chromosome Segregation. Scientific Reports 2019, 9, 15930. [Google Scholar] [CrossRef]
Liu, F.; Park, J.-E.; Qian, W.-J.; Lim, D.; Gräber, M.; Berg, T.; Yaffe, M.B.; Lee, K.S.; Burke, T.R. Serendipitous Alkylation of a Plk1 Ligand Uncovers a New Binding Channel. Nature Chemical Biology 2011, 7, 595–601. [Google Scholar] [CrossRef]
Gunasekaran, P.; Yim, M.S.; Ahn, M.; Soung, N.-K.; Park, J.-E.; Kim, J.; Bang, G.; Shin, S.C.; Choi, J.; Kim, M.; et al. Development of a Polo-like Kinase-1 Polo-Box Domain Inhibitor as a Tumor Growth Suppressor in Mice Models. Journal of Medicinal Chemistry 2020, 63, 14905–14920. [Google Scholar] [CrossRef] [PubMed]
Iliaki, S.; Beyaert, R.; Afonina, I.S. Polo-like Kinase 1 (PLK1) Signaling in Cancer and Beyond. Biochemical Pharmacology 2021, 193, 114747. [Google Scholar] [CrossRef] [PubMed]
Xing, L.; Klug-Mcleod, J.; Rai, B.; Lunney, E.A. Kinase Hinge Binding Scaffolds and Their Hydrogen Bond Patterns. Bioorg. Med. Chem. 2015, 23, 6520–6527. [Google Scholar] [CrossRef]
Gower, C.M.; Chang, M.E.K.; Maly, D.J. Bivalent Inhibitors of Protein Kinases. Crit. Rev. Biochem. Mol. Biol. 2014, 49, 102–115. [Google Scholar] [CrossRef]
Goebel, G.L.; Qiu, X.; Wu, P. Kinase-Targeting Small-Molecule Inhibitors and Emerging Bifunctional Molecules. Trends Pharmacol. Sci. 2022, 43, 866–881. [Google Scholar] [CrossRef]
Tsuji, K.; Hymel, D.; Ma, B.; Tamamura, H.; Nussinov, R.; Burke, T.R. Development of Ultra-High Affinity Bivalent Ligands Targeting the Polo-like Kinase 1. RSC Chem. Biol. 2022, 3, 1111–1120. [Google Scholar] [CrossRef]
Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K.; et al. Glide: A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Accuracy. J. Med. Chem. 2004, 47, 1739–1749. [Google Scholar] [CrossRef] [PubMed]
Alexeev, Y.; Mazanetz, M.P.; Ichihara, O.; Fedorov, D.G. GAMESS As a Free Quantum-Mechanical Platform for Drug Research. Curr. Top. Med. Chem. 2012, 12, 2013–2033. [Google Scholar] [CrossRef] [PubMed]
Nakano, T.; Kaminuma, T.; Sato, T.; Akiyama, Y.; Uebayasi, M.; Kitaura, K. Fragment Molecular Orbital Method: Application to Polypeptides. Chem. Phys. Lett. 2000, 318, 614–618. [Google Scholar] [CrossRef]
Fedorov, D.G.; Nagata, T.; Kitaura, K. Exploring Chemistry with the Fragment Molecular Orbital Method. Phys. Chem. Chem. Phys. 2012, 14, 7562–7577. [Google Scholar] [CrossRef]
Fedorov, D.G.; Kitaura, K. Second Order Møller-Plesset Perturbation Theory Based upon the Fragment Molecular Orbital Method. J. Chem. Phys. 2004, 121, 2483–2490. [Google Scholar] [CrossRef]
Fedorov, D.G.; Kitaura, K.; Li, H.; Jensen, J.H.; Gordon, M.S. The Polarizable Continuum Model (PCM) Interfaced with the Fragment Molecular Orbital Method (FMO). J. Comput. Chem. 2006, 27, 976–985. [Google Scholar] [CrossRef]
Fedorov, D.G. Solvent Screening in Zwitterions Analyzed with the Fragment Molecular Orbital Method. J. Chem. Theory Comput. 2019, 15, 5404–5416. [Google Scholar] [CrossRef]
Gaus, M.; Goez, A.; Elstner, M. Parametrization and Benchmark of DFTB3 for Organic Molecules. J. Chem. Theory Comput. 2013, 9, 338–354. [Google Scholar] [CrossRef]
Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577–8593. [Google Scholar] [CrossRef]
Hoover, W.G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [PubMed]
Martyna, G.J.; Tobias, D.J.; Klein, M.L. Constant Pressure Molecular Dynamics Algorithms. J. Chem. Phys. 1994, 101, 4177–4189. [Google Scholar] [CrossRef]
Humphreys, D.D.; Friesner, R.A.; Berne, B.J. A Multiple-Time-Step Molecular Dynamics Algorithm for Macromolecules. J. Phys. Chem. 1994, 98, 6885–6892. [Google Scholar] [CrossRef]
Lim, H.; Jin, X.; Kim, J.; Hwang, S.; Shin, K.B.; Choi, J.; Nam, K.-Y.; No, K.T. Investigation of Hot Spot Region in XIAP Inhibitor Binding Site by Fragment Molecular Orbital Method. Comput. Struct. Biotechnol. J. 2019, 17, 1217–1225. [Google Scholar] [CrossRef]
Abel, R.; Young, T.; Farid, R.; Berne, B.J.; Friesner, R.A. Role of the Active-Site Solvent in the Thermodynamics of Factor Xa Ligand Binding. J. Am. Chem. Soc. 2008, 130, 2817–2831. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Domain information of the five human polo-like kinase (PLK) families and full-length structure of the PLK1 predicted by Alphafold2. (A) Numbers represent the amino acid sequence numbers. (B) Illustration of the kinase (blue) and polo-box (orange) domains.

Figure 2. Structure-based selectivity pattern considerations. (A) Superposition of 13 crystal structures in the holo-form of human polo-like kinase 1 (PLK1). The basis on para and meta, PLK1 kinase domain inhibitors can be divided into two groups. The inhibitors from the B group are PLK1 selective inhibitors toward PLK2 and 3. They all interact with Glu140 in PLK1. (B) Two-dimensional chemical structure of inhibitors from two groups. Hinge-binding regions are highlighted in light yellow, solvent exposure regions in light blue, and phosphate regions in light pink.

Figure 3. Characteristics of the ATP-binding pocket of the kinase domain in polo-like kinase 1 (PLK1). (A) Crystal structure of the human PLK1 kinase domain with ATP (protein data bank ID: 2OU7). The ATP-binding site is divided into adenine, ribose, and phosphate pockets and solvent Channels 1 and 2. (B) Comparison of the conformation change of the two crystal structures of PLK1 complex with two inhibitors, BI2536 and Onvansertib. Nonselective inhibitor BI2536 only occupies solvent Channel 1, but potent PLK1 selective inhibitor Onvansertib binds solvent Channels 1 and 2.

Figure 4. Fragment molecular orbital (FMO) analysis of ATP. (A) FMO results of the crystal structure of the polo-like kinase 1 (PLK1) complex with the ATP analog. The ligand is blue; key residues are green. (B) Pair interaction energy (PIE) values of the significant residues in the ATP-binding site. (C) PIE decomposition analysis of these critical interactions.

Figure 5. Fragment molecular orbital (FMO) analysis of BI2536 and Onvansertib. (A) FMO results of the crystal structure of the polo-like kinase 1 (PLK1) complex with BI2536. The ligand is blue; the key residues are green. (B) FMO results of the crystal structure of the PLK1 complex with Onvansertib. The ligand is light pink, key residues are green, and nitrogen and oxygen atoms are blue and red, respectively. (A, B) Middle bar plots describe the pair interaction energy (PIE) values of the significant residues in the ATP-binding site, whereas bottom bar plots describe the PIE decomposition analysis of these critical interactions.

Figure 6. Docking structure and fragment molecular orbital (FMO) analysis of GSK461364. (A) Docking structure of GSK461364. The ligand is purple; key residues are green. (B) Docking structure of GSK461364 overlayed with the crystal structure of Onvansertib (light pink). (C) Bar plot describing the pair interaction energy (PIE) values of the significant residues in the ATP-binding site and GSK461364. (D) The bar plot describes the PIE decomposition analysis of these critical interactions.

Figure 7. Superposition of the kinase domain of polo-like kinase (PLK) 1–3. Superposition of the crystal structures of the PLK1-3 and sequence alignment. Four residues occur around the ATP-binding site. Glu140 of PLK1 has a negative charge, except His169 and His149 of PLK2, and three have a positive charge.

Figure 8. Heatmap of the protein-ligand contact from the molecular dynamic simulation of 13 ligands. Percentages represent interactions occurring during the simulation. Direct interaction <30% is yellow, 30% to 60% is green, and >60% is red. The water bridge is light blue.

Figure 9. WaterMap analysis of the ATP-binding pocket of the kinase domain. Hydration sites are represented as spheres, with colors reflecting the predicted associated free energies. Green spheres signify favorable free energies, whereas red spheres indicate unfavorable free energies. (A, B) Onvansertib (light pink) is overlayed in the WaterMap of the crystal structure of BI2536, focusing on the hydration site of the phosphate-binding site and solvent Channel 2. (C, D) WaterMap of the crystal structure of Onvansertib.

Figure 10. Characteristics of the polo-box domain (PBD). Pockets of the PBD are divided into the phosphate (blue), pyrrolidine (orange), and Tyr-rich (magenta) pockets.

Figure 11. Hot-spot analysis of substrate peptide complex with the polo-box domain (PBD) of polo-like kinase 1 (PLK1). (A) The substrate peptide was divided into nine fragments: Ala01, p-Thr2, Ser3, His4, Leu5, Pro6, Pro7, Asp8, and Phe9 (green and yellow sticks). (B) PLK1 PBD hot-spot residues are blue, orange, and light pink sticks. Pair interaction energy (PIE) values are described in the right table.

Figure 12. Hot-spot analysis of substrate peptide complex with the polo-box domain (PBD) of polo-like kinase 1 (PLK1). (A) Hot-spot residues of substrate peptide: Ala01, P-Thr2, Ser3, His4, Leu5, Pro6, Pro7, Asp8, and Phe9 (green and yellow sticks). (B) PLK1 PBD hot-spot residues are blue, orange, and light pink sticks. Pair interaction energy (PIE) values are described in the right table.

Figure 13. Hot-spot analysis of the 4j complex with the polo-box domain (PBD) of polo-like kinase 1 (PLK1). (A) Fragments of 4j: p-Thr, Ser, F-Akyl, Leu, and Pro (green and yellow sticks). (B) Significant interaction residues in the PBD pocket are blue, orange, and light pink sticks. Pair interaction energy (PIE) values are described in the right table.

Figure 14. Fragment molecular orbital (FMO) analysis of KBJK557 and KBJK-4a. (A) FMO results of the structure from Frame 489 of the molecular dynamic simulation. The ligand is blue; key protein residues are blue, orange, and light pink sticks. (B) Bar plot describing the pair interaction energy (PIE) values of the significant residues in the polo-box domain (PBD). (C) Bar plot describing the PIE decomposition analysis (PIEDA) of these critical interactions. (D) FMO results of the docking structure of the KBJK-4a. The ligand is blue; key protein residues are blue, orange, and light pink sticks. (E) Bar plot describing the PIE values of the significant residues in the PBD. (F) Bar plot describing the PIEDA of these critical interactions.

Figure 15. Heatmap of the fragment molecular orbital (FMO) results of six structures. The FMO results of the complex structure of the polo-box domain (PBD) of polo-like kinase 1 (PLK1) complex with substrate peptides 4j and 4a, KBJK557, and KBJK-4a. The critical residues are highlighted blue for the phosphate pocket, orange for the pyrrolidine pocket, and pink for the Tyr-rich pocket. The pair interaction energy was summed for each residue. Darker red indicates a lower energy value.

Figure 17. WaterMap analysis of the binding pocket of the polo-box domain (PBD). Hydration sites are represented as spheres, with colors reflecting their predicted free energies. Green spheres signify favorable free energies; red spheres indicate unfavorable free energies. (A) WaterMap analysis of the 4a (lime) binding structure in the PBD. (B) Focusing on the Tyr-rich pocket in the 4a (cyan surface) water map analysis with 4j superposition (orange stick). (C) WaterMap analysis of the crystal structure of the 4j binding in the PBD. (D) WaterMap focusing on KBJK (surfaced with a partial charge) and hydration sites in the Tyr-rich pocket.

Table 1. PLK1 Kinase Domain Inhibitors in Clinical Trials.

Drugs	Potency	Selectivity to PLK2 and PLK3	Clinical Phase
BI2536	PLK1 IC₅₀ = 0.83nM	PLK2 IC50 = 3.5nM	II
BI2536	PLK1 IC₅₀ = 0.83nM	PLK3 IC₅₀ = 9.0nM	II
Volasertib	PLK1 IC₅₀ = 0.87nM	PLK2 IC₅₀ = 5nM	I/II/III
Volasertib	PLK1 IC₅₀ = 0.87nM	PLK3 IC₅₀ = 56nM	I/II/III
Onvansertib	PLK1 IC₅₀ = 2nM	PLK2 IC₅₀ > 10 000nM	II/III
Onvansertib	PLK1 IC₅₀ = 2nM	PLK3 IC₅₀ > 10 000nM	II/III
TAK-960	PLK1 IC₅₀ = 0.8nM	PLK2 IC₅₀ = 16.9nM	I
TAK-960	PLK1 IC₅₀ = 0.8nM	PLK3 IC₅₀ = 50.2nM	I
GSK461364	PLK1 Kiapp<0.5 nM	PLK2 Kiapp = 860 nM	I
GSK461364	PLK1 Kiapp<0.5 nM	PLK3 Kiapp = 1000 nM	I

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Alerts

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer