AlphaSpace 2.0: Representing Concave Biomolecular Surfaces using Beta-Clusters

A. Beta-cluster representation of pockets

In order to generate a pseudomolecular representation of fragment-centric pockets detected by AlphaSpace with better connectivity properties, we cluster alpha-atoms using their pairwise distances and the complete linkage method⁴⁸ in the hierarchical clustering implementation of Scipy package⁴⁹ with a distance cut-off of 1.6 Å, a typical aliphatic carbon-carbon bond distance, to eliminate unphysical and extremely close-contact alpha atoms while retaining its overall shape. The complete linkage method generates spherical clusters and each cluster centroid is represented as a beta-atom. The combination of all beta-atoms in a given fragment-centric pocket forms a beta-cluster, which is the pocket’s pseudomolecular representation, as illustrated in Fig 1B. The alpha-space value for each alpha atom in the cluster is summed to obtain the alpha-space of the beta-cluster. Other beta-cluster attributes, including the solvent accessible surface area and the volume, are calculated the same as for a real molecule with a van der Waals radius of 1.6 Å for each beta atom.

B. Beta Cluster Score (BScore).

With beta-cluster mimicking potential molecular binder of a fragment-centric pocket detected by AlphaSpace, we define beta-score, a maximum docking score of the beta-cluster, as an indicator of pocket ligandability. Based on the AutoDock Vina scoring function⁵⁰, a probe atom corresponding to various Vina atom types (C, N, O, F, P, S, Cl, Br, I) is placed at each beta atom location and AutoDock Vina score is calculated for that probe atom. The beta atom type and score are determined by the probe type with the best vina score (highest affinity). Beta atom types can also be generally classified as polar/non-polar and hydrogen bond donors/acceptors based on the strength of the Vina interaction terms. Finally, the BScore is obtained by summing the best vina score of each beta atom in the beta cluster, which provides an estimate of pocket ligandability.

C. Pocket mapping with aligned protein structures.

Given several aligned protein structures and their determined beta-atoms, consistent pocket mapping among different structures can be achieved by determining ensemble pockets. First, several proteins are aligned together, and each one is processed separately to generate corresponding beta atoms along concave surfaces. Second, all determined beta atoms from different structures are superimposed based on structural alignment and re-clustered altogether by hierarchical clustering using the average linkage method⁵¹. This places spatially close beta-atoms from different protein structures into an ensemble beta-cluster. After the clustering step, each individual beta atom in ensemble beta-clusters is re-mapped to its corresponding protein structure. Statistical properties such as the mean, median, or variance of the distribution pocket attributes can be calculated for the ensemble beta-clusters adding more information for the analysis of concave surfaces. A collection of individual pockets coming from the same ensemble beta-cluster are considered to be matched and forms an ensemble pocket. One advantage of ensemble pockets over our previous pocket matching method is the use of structural alignment of the protein to generate consistent pocket definitions that are less sensitive to minor structural changes on the protein surface. Because the alignment is done before the clustering step, it allows pocket-mapping among homologous proteins that share a similar protein fold.

In this work where we apply the ensemble pocket approach to identify the same pockets across different surface states or different homologous proteins, we used the Matchmaker⁵² program implemented in Chimera⁵³. Matchmaker allows for the superimposition of protein structures using their global sequence alignment or a user-specified alignment based on select residues from the target structures. Restricting the alignment to active site or interface residues allows for consistent matching of surface pockets that are less affected by large conformational changes observed for multidomain proteins. We employed a two-step process to superimpose the protein structures: an initial superimposition was performed using the global sequence alignment, if the global RMSD of the aligned structures fell below 1.5 Å then the structure pairs were considered sufficiently aligned for processing of the ensemble pockets; then for structure pairs that have RMSD >1.5 Å , a second superimposition step was performed using residues selected to be within 5.0 Å from the largest ligand-bound structure in the ensembles to be analyzed. Finally, visual inspection of beta atoms on their respective protein surface was performed to ensure alignment quality.

D. Beta-Cluster and Ligand Features.

Since beta-cluster is sufficiently molecule-like, it allows the direct comparison between beta-cluster and ligand. Meanwhile, a variety of features can be easily calculated for both beta-clusters and ligands, as listed in Table 1, which include geometric features like surface area and volume, as well as chemoinformatic features like those from ChemoPy⁵⁴. To show the “ligand-like” properties of the beta clusters, we also calculate shape similarity features using the distribution of moments with the USR formulation^55–57.

E. Datasets

The small-molecule dataset was constructed based on the protein-ligand complexes retrieved from the 2013 release of the scPDB data sets⁶¹. First we screened the scPDB 2013 release and removed structures that did not have any binding data according to the 2016 release of the PDBbind refined set⁶², which reduces the number of protein-ligand structures from ~10,000 to 996. The scPDB/PDBbind refined set was further reduced by aggregating all protein family structures using their UNIPROT accession number. The protein-ligand complex structure corresponding to the highest affinity ligand from PDBBind was selected as the representative structure for that UNIPROT group. This non-redundant dataset is called the small-molecule dataset and is composed of 293 structures with properties listed in SI. This data set has been used to determine the clustering parameters for beta-cluster and to investigate correlations between the properties of the ligand and the binding site as represented by the beta-cluster.

To demonstrate applicability of beta-cluster representation, we applied our beta-cluster approach to various data sets that are relevant to fragment-based drug-design (FBDD) and identification of cryptic pockets. The Ichihara Set⁶³ was used to investigate the ability of beta-clusters to make meaningful suggestions for FBDD, which consists of 21 lead/fragment pairs of structures where a lead ligand retains its fragment-binding pose. We applied BScore and pocket matching on a reduced set of the latest release of the 2P2I dataset⁶⁴ in which, only the structure of the maximum affinity small-molecule PPI inhibitor is analyzed for each PPI interface. The CryptoSite database⁶⁵ was analyzed to demonstrate the pocket matching capabilities of beta clusters across multiple states of the same protein. The CryptoSite database consists of 79 pairs of unbound (apo) and ligand bound (holo) structures of the same protein selected for their structural and functional diversity. Calmodulin and Ca²⁺ ATPase 1 were excluded due to the lack of a good structural alignment between the apo and holo pairs while Kynurenine aminotransferase II and Serum transferrin were excluded due to the small size and charged nature of their cognate ligands. This leaves 75 pairs of structures in the CryptoSite dataset for our analysis.

Beta-clusters Mimic the Shape of Binding-Site Ligands

To establish the pseudomolecular nature of beta-cluster, we determined shape similarities and property correlations of beta clusters against their cognate ligand partners for the high-quality, small-molecule dataset which consists of 293 protein-ligand complexes. The results in Fig 2. show that there are good correlations between ligand and beta cluster features, including volume, surface area and span. This clearly demonstrates the ability of beta clusters to mimic the shape of potential molecular binders. Meanwhile, one basic assumption in rational inhibitor design is that similar binding sites would interact with similar ligands. Thus we are motivated to directly examine correlations between ligand similarities and beta-cluster similarities using two similarity measures: one is USR similarity which is a shape similarity based on statistical moments of pairwise distance distribution derived from the USR features^{56, 57}; the other is OASA similarity which is based on the total and pharmacophore-atom types of the occluded surface area of the protein surface⁶⁰ encoded into a feature fingerprint. Both are alignment independent measures, which allow for structural patterns to be compared without the need for computationally exhaustive methods like clique detection or geometric hashing. The results in Figure 3 show the correlation of the similarity features with an R of 0.74 for the USR similarities and 0.88 for the OASA fingerprint. Non-contacting beta atoms and ligand atoms contribute to the noise and reduce the shape similarity for the USR features while the OASA fingerprints suffer from less noise due to common occluded surface area for both beta clusters and ligands. Fig. S1 illustrates two examples how highly similar ligands bind to pockets represented by highly similar beta clusters. Only considering overlapped beta atoms and ligand atoms greatly increases the correlation of USR similarity with an R of 0.94 while the OASA similarity fingerprint also increases to 0.97, which indicates that beta clusters can mimic important pocket-occupying portions of the ligand to a high degree.

Using Beta-cluster Score (BScore) to Estimate Pocket Ligandability

Besides quantifying geometric similarity between beta-clusters and ligands, we have analyzed their energetic similarity using the AutoDock Vina Scoring function by plotting the ligand vina score against the corresponding contact pockets’ BScore. The results in Fig. 4a show that BScore tends to be more negative than the calculated vina score for the ligand among all structures in the small-molecule dataset except for 1QK4. This is not surprising since the contact pockets tend to have unoccupied beta atoms that would otherwise contribute to favorable protein-ligand interactions. Meanwhile, since the overall shape of beta clusters is derived from the Voronoi tessellation of the concave binding site, they naturally represent a maximum possible pseudomolecule packing for that site. Thus, it would be reasonable to consider BScore as an estimate of pocket ligandability. Ligands that have almost complete beta cluster coverage tend to fall closer to the diagonal line indicating that these ligands are well optimized energetically to their binding sites. On the other hand, ligands that have large Vina Score affinity differences with their beta clusters tend to have low ligand coverage of the binding site and can possibly be further optimized by maximizing the overlap of the ligand and the beta cluster as shown in the case of 3ISS in Fig 4b. The individual energetics of the beta atoms also reveals characteristics of localized sites that can facilitate rational ligand design.

Detected Pocket Communities Represent Potential Binding Sites

We used BScore and alpha space volume to rank detected pocket communities for the small-molecule dataset. Table 2 shows the performance of the community ranking in predicting the small-molecule binding site as highest ranking communities. Top 1 ranking communities based on the alpha space volume and the BScore were predicted to be the small-molecule binding site of the cognate ligand with a success rate of >78%. We observe that all small molecule binding sites can be detected by AlphaSpace as one pocket community. Top 3 ranking communities include the small-molecule binding site for any criteria with a success rate of >95%. The last entry in Table 2 shows the success rate by selecting top ranking communities based on either the alpha-space or the BScore criteria. By using the union of top ranking communities from each criteria , the top 1 success rate improves to 87% while the top 3 success rate improves to 98%. These results further confirm that AlphaSpace can be a useful tool in detecting potential binding sites, which allows the prioritization of high-ranking communities for subsequent studies.

It should be noted that high or mid ranking communities that do not overlap with the binding site ligands may have other functional roles such as allosteric binding sites⁶⁶ or protein-protein/peptide interfaces. Cyclin-dependent kinase 2 (CDK2) provides an excellent example to illustrate how AlphaSpace and BScore can be used to identify multiple functionally important binding interfaces for a single protein. Based on the CDK2 monomer structure obtained from the CDK2/Cyclin A complex⁶⁷, in which CDK2 is bound to ATP and a short substrate peptide, we carried out AlphaSpace analysis and detected 16 different pocket communities, as shown in Fig. 5, and their properties are listed in Table 3. We can see that the ATP binding site associated community (community 1) is ranked at the top in terms of both the BScore and alpha-space, followed by one of the substrate binding communities (community 2). Among 13 remaining detected communities, by analyzing other available CDK2 and related CDK crystal structures, we find that 6 of these are indeed involved in binding other proteins, peptides, and even fragments as shown in Table 3 and Fig 6. For example, important binding residues from the p27 peptide inhibitor⁶⁸ are found to overlap with the 5^th ranking community. Adjacent to this community is a low-ranking community that binds a conserved helix-loop found in all cyclins. The 4^th and 9^th ranking communities formed the cavity in the C-lobe and A-loop region, which is found to bind loops and helices from Cyclin A⁶⁷ and Cyclin E⁶⁹ that may stabilize the open A-loop conformation of the active state of CDK2. Alignment of the CDK1/CKS2⁷⁰ complex to CDK2 shows that the 13th community belongs to the binding site for this regulatory protein. Lastly, a high-scoring fragment-binding site, the 7^th community, is found at the solvent exposed region of the C-lobe; a largely uncharacterized region of kinases that may have some functional importance as an allosteric site based on its ability to bind small molecules as shown in multiple crystal structures^71–73. In total, 9 out of the 16 detected communities are confirmed to be involved in true binding sites which indicates that the AlphaSpace community definition is sufficiently sensitive for potential binding site identification.

Beta-Cluster can Guide Ligand Optimization

The results for the retrospective analysis of Fragment to Lead pairs from the Ichihara set are listed in Table 4. The Ichihara Set⁶³ is composed of 21 pairs of ligand bound structures covering various important drug targets. For each pair, its fragment member forms the “core” scaffold of the ``lead” and both members have the similar binding mode. Here we employ this dataset to demonstrate how beta-cluster can be used to help guide ligand optimization. For 16 out of 21 systems, the highest-ranking community is the small-molecule binding site with BScore ranging from −15.5 kcal/mol for DPP to −42.8 kcal/mol for PPAR. The large affinity of PPAR is not surprising considering that this target is known to be highly druggable with several contiguous pockets in the binding site that determine the binding mode and activity of PPAR modulators. For 19 out of 21 structures the small-molecule binding site is found among the top 2 ranking communities, thus the community ranking demonstrates the ability of beta clusters to detect and enrich large surface regions of interests which can be useful for focused, high-throughput studies involving only ligandable sites. Within each binding community, the sub-pocket partitioning of AlphaSpace allows for the identification of localized interaction sites that is the preferred binding site for the initial fragment. These core pockets are akin to hot-spots found in PPIs and function as “seeds” for the fragment-based drug design process. Assigning core pockets allows for finding complimentary ligand shapes and pharmacophoric interactions that can maximize initial fragment binding which would carry into an elaborated lead molecule. Out of the 21 systems, only p38 and PKB have core binding community pockets not in the top 1 rank. Nonetheless, all 21 fragment associated core pockets are in the top 2 indicating the ability of BScore to identify the most important pockets for fragment binding. Using BScore ranking information at the global and local level, beta clusters can be used as a prospective tool for the judicious selection of ligandable communities and preferred fragment-targetable pockets within the community.

We also apply the concept of optimizability of the binding to the Ichihara dataset, i.e. %optimized, which is calculated as the ratio between the BScore of contact beta atoms within 1.8 Å of any ligand atoms and sum of the BScore of all the beta atoms of the binding site associated community. The higher the value of %optimized indicates the greater overlap between ligand and beta-cluster, i.e., the higher pocket-ligand complementarity. On the other hand, a lower %optimized value suggests that there is more room for further inhibitor optimization. Fig. 7A shows the %optimized for each fragment/lead pair in the dataset while Fig 7B shows the relationship between the change in the %optimized scores with the normalized experimental affinities in going from the fragment to the lead structure. For all systems investigated in the Ichihara set, an improvement in the %optimized score of the fragment to the lead is correlated with an improvement in the affinity as measured by the pIC50 values compiled by Ferenczy and Keseru⁷⁴.

To illustrate how AlphaSpace2 would be used for ligand optimization, we performed a more detailed retrospective analysis of the PDE and jak2 fragment/lead pairs. Fig 8 shows the overlay of the binding site communities for PDE and jak2 with the fragment and lead pairs shown along with other high-affinity ligands. For both systems, the initial fragment is observed to partially occupy several high-scoring pockets, with neighboring unoccupied pockets highlighted in orange which would be preferable sites for elaboration. For PDE, a nitrophenyl group was used to target this pocket which improves the affinity of the elaborated ligand from 82 μM to 21 nM⁷⁵. For jak2, the indazole core fragment was grown with a phenyl moiety to target this neighboring pocket, which increased the affinity from 40.9 μM to 1.6 μM⁷⁶. After this first iteration, another unoccupied pocket is still accessible for further fragment extension, which was subsequently utilized by the final lead compound with a butylsulfonamide moiety. This leads to a potent inhibitor with an affinity of 78 nM. We can see that the fragment-centric pocket representation can clearly indicate those preferred sites for further fragment extension. Considering that there remain unoccupied pockets in PDE and jak2 for the lead member in the Ichihara Set⁶³ , we have searched the PDBBind database and literature for other ligands that can possibly target remaining unoccupied pockets in PDE and jak2 to further improve potency, as shown in Fig. 8. For the PDE lead, the original authors have already performed a QSAR study⁷⁷ for the initial carboxyl pyrazole scaffold with R₁=Phenyl, R₂=−4-Chloro-phenyl, R₃=-H. This particular substitution extends into a mid-scoring pocket and has the second highest potency amongst the first-round optimizations (150 nM compared to 31 nM) but was unfortunately not further optimized. Comparison with the ligand from PDB ID: 4WCU indicates that moieties from R₂ can be extended to interact with this pocket and thus possibly increase the potency for this lead⁷⁸. Likewise, a phenyl fragment from the ligand in PDB ID: 5K1I extends into this partially occupied pocket, while a benzoic acid like fragment extends into another high-scoring pocket greatly contributing to ligand binding affinity⁷⁹. Jak2 appears to be further optimizable based on the presence of two nearby high-scoring pockets close to the E1 and gatekeeper regions of the kinase binding site. This can be observed by the overlay with the ligand of PDB ID: 4JI9, in which its methyl-piperazine moiety extends into the unoccupied E1 pocket⁸⁰. Similarly, a fluoro-hydroxy phenyl fragment of the ligand from PDB ID: 5TQ8 points towards the high-scoring pocket associated with the gatekeeper site⁸¹. Both unoccupied pockets can be easily elaborated from the aminoindazole core of the selected Jak2 lead molecule which is a shared scaffold of the ligand in 5TQ8. Finally, we also extended our community based analysis to other systems in the Ichihara set that had an excess of >5 pockets that would be ideal candidates for further elaboration of either the initial fragment or the more potent lead ligands. Fig. S2 shows the overlay of these selected systems with the potency of the leads from the Ichihara set as well as more elaborated ligands from the PDB that co-target nearby unoccupied pockets. We can see that the more advanced ligands not only occupy lead-molecule binding pockets but also display extensive overlaps to nearby, high-scoring, unoccupied pockets suggesting that AlphaSpace and beta clusters not only identify preferred sites of initial fragment binding but can also be used to guide the optimization of lower affinity fragments into more potent lead molecules.

Analysis of Protein-Protein Interactions using Beta Clusters

Pocket matching was performed for the apo, inhibitor-bound (iPPI), and holo-PPI complexes of the 2P2I set to determine how certain pocket structures are conserved in various states of the PPI. All of the iPPI associated matched ensembled pockets were included in the analysis across the different surface states of the 2P2I set. Fig 9 shows BScores of pockets for various proteins in the 2P2I. Bcl-2 and Bcl-xL have the best BScore for iPPI pockets consistent with the wide availability of high-affinity inhibitors for these systems. On the other hand, TNFa has the worst BScore for its iPPI pockets followed by H-Ras and K-Ras, all of which are targets that have been well known to be challenging for inhibitor design. The variability in BScore among different surface states indicates the effect of protein plasticity on ligandability of iPPI binding sites. For 20 out of 21 systems, BScore difference among apo, iPPI, and PPI beta clusters is less than 8 kcal/mol indicating that “hot spot” regions are conserved for various states of the PPI interface. Only Bcl-xL shows a drastic BScore difference between apo and iPPI/PPI bound states, nonetheless the apo state of Bcl-xL may be still more ligandable than some other PPIs. It should be noted that other authors have commented on the flexibility of Bcl-2 family proteins and how ligands can induce different conformations from the apo state, leading to the formation of more druggable, fragment-centric cavities^{82, 83}.

Meanwhile, we examined complementarity of the ligand partners for the iPPI and PPI associated pockets, as shown in Fig. 10. Only 2 out of 21 systems with %optimized value higher for the PPI structures than for the iPPI structures indicates that, at least for those iPPI associated pockets, small-molecule ligands have been optimized to have higher complementarity for the binding site to some extent. The better complementarity for iPPI structures is not surprising considering that small-molecules take advantage of shape moieties not available to natural amino-acids. On the other hand, shallow binding sites such as in cIAP1, XIAP, and ZIPa can only accommodate similarly flat molecular shapes thus showing similar %optimized values for iPPI and PPI structures.

Pocket Comparison for the CrypticSite Dataset

One-to-one pocket matching was performed for apo/holo pairs of the CryptoSite dataset and BScores of binding site pockets were calculated and compared as shown in Fig 11A. As expected, the apo state of the proteins tends to be less ligandable than the holo state, with the mean BScore of −12.1 kcal/mol for the apo state versus −18.8 kcal/mol for the holo state. In total, 72 proteins have improved holo state Bscores while only 3 proteins have apo states having better BScores with difference less than 3.0 kcal/mol. On the other hand, there are 17 highly perturbed apo/holo pairs, with the holo state more ligandable by BScore difference larger than 10.0 kcal/mol. To quantify the structural changes between the apo and holo states, we calculate the OASA fingerprint similarity scores between the matched contact pockets of the apo/holo pairs. The distribution of the OASA similarity scores for the apo/holo pairs is shown in Fig. 11B. We observe 21 apo/holo pairs that have high OASA similarity scores >0.80. Visualization of the top 3 most similar apo/holo pairs reveals similar surface states at the binding site with global RMSD differences all less than 0.6 Å. On the other hand, there are 13 apo/holo pairs with OASA similarities less than 0.5. Examination of three least similar apo/holo pairs reveals that the main binding surface difference between apo/holo pairs is the loss of pockets in the apo state due to occlusion by loops or large movements in the secondary structure.

As shown in Fig. 11A, some apo-structures in the CrypticSite dataset consist of ligandable binding sites with BScore more negative than −20.0 kcal/mol, which suggests that these apo-structures may already be targetable in the first place consistent with the extended analysis of the CryptoSite by Vajda et al using their hot-spot mapping method FTMAP⁸⁴.

Pocket Comparison of NMR Structures for MDM2 and MDM4

The oncogene proteins MDM2 and MDM4 play an important role in many cancers, functioning as a negative regulator of the p53 tumor suppressor protein^85–87. Both of them can bind a conserved α-helix segment of p53⁸⁸. Because of their importance in cancers, dual targeting of MDM2 and MDM4 proteins has been an active area of research, yielding small-molecule PPI disruptors, such as Nutlins^{89, 90}. Here we applied our ensemble pocket approach to help understand similarities and differences between binding interfaces of MDM2 and MDM4 on a total of 40 NMR structures from both MDM2/p73⁹¹ (PDB: 2MPS) and MDM4/p53⁹² (PDB: 2MWY) by using the first structure of MDM2 as the reference for alignment. From the results in Fig 12, we can see that TRP, PHE and LEU pockets corresponding to hot-spot residues of binding partners for MDM2/MDM4 are important local interaction sites for MDM2 and MDM4 based on BScore. Among them, the LEU pocket is slightly more favorable for binding in MDM2 compared to MDM4 which has a BScore difference of ~1.0 kcal/mol. Among the rest of the detected pockets at the interface, pockets 8 and 9 are uniquely associated with MDM4, which have favorable levels of ligandability of −3.1±1.2 and −1.7±0.9 kcal/mol respectively. On the other hand, pocket 7 of the MDM2 interface is more ligandable than its MDM4 counterpart with a score of −2.5±0.6 kcal/mol. Of the remaining pockets, the mid-scoring pocket 1 is optimally positioned near the LEU pocket and may be further explored for the design of dual MDM2/MDM4 PPI disruptors. Throughout this analysis, we leveraged the ability of our ensemble pockets to consistently assign pockets sets for the interface allowing for direct comparison of MDM2 and MDM4 despite that their PPI interface residues are not the same.

AlphaSpace2 Studies of Bcl-2 family Proteins

Finally, we applied beta AlphaSpace2 to study an important class of PPI related disease targets, Bcl-xL and Bcl-2. Both targets belong to the Bcl-2 family proteins and act as regulators of apoptotic pathways by binding and sequestering pro-apoptotic signaling ligands⁹³. Increase in expression of Bcl-2 or Bcl-xL is associated with resistance to pro-death signals, contributing to multiple diseases such as cancers and tumorigenesis⁹⁴. By applying AlphaSpace to map the Bcl-xL/BAK (PDB ID: 1BXL) PPI interface, we were able to detect and assign hot-spot associated pockets. As shown in Fig 13, twelve Bcl-xL pockets are detected at the PPI interface and are labeled P1 to P4 corresponding to the hot-spot residue annotation for Bcl-2 family proteins and pockets 1 to 8 corresponding to additional pockets on the interface. We can see that BScores are able to correctly predict and rank hot-spot residue associated pockets as the top 4 ranking pockets with scores more negative than −3.0 kcal/mol, which is very consistent with experimental ΔΔG from alanine scanning mutation experiments of the BAK peptide⁹⁵. This demonstrates the ability of BScore in correctly identifying important hot-spot interaction sites. The %-optimized score of each pocket also shows the potential to improve ligand/pocket complementarity based on the overlap of ligand and beta atoms. Among hot-spot pockets, P4 has the lowest %optimized score of 20% while other hot-spot pockets have %optimized scores >60%, suggesting a path to improve the affinity of BAK based peptides by designing mutants on the I85 position that maximize the occupation of P4. We also identified mid-scoring, auxiliary pockets 2 and 4 which are adjacent to the hot-spot, core pockets that may further improve the affinity of the BAK peptide. Finally, mid-scoring unoccupied pockets 1, 3, and 5 might be further utilized in the design of non-native peptide-based Bcl-xL PPI disruptors.

There has been significant interest to develop Bcl-2 or Bcl-xL specific inhibitors, while the main challenge comes from the close homology and sequence similarity of Bcl-2 family proteins⁹⁶. Here we applied AlphaSpace to analyze BAX bound complexes of Bcl-2⁹⁷ (PDB ID: 2XA0) and Bcl-xL⁹⁸ (PDB ID: 3PL7), as shown in Fig 14. The ensemble pocket clustering is able to assign the important P1, P2, P3, and P4 pockets on both Bcl-2 or Bcl-xL surfaces with a score of ~2.5 kcal/mol while also assigning unique pockets specific to either Bcl-2 or Bcl-xL. Among the conserved hot-spot pockets, significant BScore differences observed for P1, P3 and P4 pockets indicate that these pockets may contribute to differential binding profiles between Bcl-2 and Bcl-xL. The more favorable BScore for the P4 pocket of Bcl-xL suggests that this pocket contributes strongly to the affinity of Bcl-xL ligands while this interaction site may be dispensable for Bcl-2. Indeed, Elmore et al⁹⁹.were able to design Bcl-2 specific inhibitors that intentionally avoided the P4 site causing a drastic drop in affinity for Bcl-xL while only leading to a modest change for Bcl-2, ultimately leading to the design of the first FDA approved Bcl-2 family inhibitor, Venetoclax.