as Mcl-1 Inhibitors: A Comparative Study of Experimental and Computational Modelling Data

: The human myeloid leukemia cell differentiation protein (hMcl-1) is an anti-apoptotic multi-partner protein, belonging to the B-cell lymphoma-2 (Bcl-2) family of proteins. Studies have linked hMcl-1 alleviated expression with resistance to hemopoietic chemotherapeutics, which makes it a key drug target in blood cancers. However, most of the developed small-to medium-sized hMcl-1 inhibitors have typical off-target activity towards other members of the Bcl-2 family. To improve the hMcl-1 inhibitor design, especially exploring a suitable scaffold with pharmacophoric features, we focused on natural hMcl-1 inhibitors. To date, seven classes of natural compounds have been isolated, which display a low micromolar afﬁnity for hMcl-1 and have limited biophysical studies. We screened hMcl-1 co-crystal structures, and identiﬁed nine co-crystal structures of hMcl-1 protein, which were later evaluated by multiple receptor conformations (which indicates that the differences between hMcl-1 in crystal structures are low (RMSD values between 0.52 and 1.13 Å, average RMSD of 0.638–0.888 Å, with a standard deviation of 0.102–0.185Å)), and multiple ligand conformations (which led to the selection of the PDB structure, 3WIX (RMSD value = 0.879 Å, standard deviation 0.116 Å), to accommodate various Mcl-1 ligands from a range of co-crystal PDB ﬁles) methods. Later, the three adopted docking methods were assessed for their ability to reproduce the conformation bound to the crystal as well as predict trends in Ki values based on calculated RMSD and docking energies. Iterative docking and clustering of the docked pose within ≤ 1.0 Å was used to evaluate the reproducibility of the adopted docking methods and compared with their experimentally determined hMcl-1 afﬁnity data.


Introduction
The human myeloid leukemia cell differentiation protein, hMcl-1, is a multi-partner, anti-apoptotic protein, belonging to the B-cell lymphoma-2 family (Bcl-2 family) of proteins [1].Several reports have established a link between the upregulation of hMcl-1 with the development of resistance to various anticancer chemotherapies [2][3][4][5][6].Furthermore, amplification of the hMcl-1 gene is one of the most commonly observed genetic aberrations in human blood cancers [7], which makes it a key drug target in hemopoietic malignancies.Additionally, its RNA-mediated knockdown and Mcl-1 gene silencing demonstrated a substantial decrease in cancer cell growth and restored the sensitivity of chemo-resistant cells [2,8,9] to chemotherapy.This Bcl-2 family comprises anti-apoptotic proteins other than hMcl-1 (especially Bcl-2 itself, Bcl-x L [10], and Bcl-w), as well as proapoptotic and apoptotic proteins (Bax, Bak, Bik, Bid, and NOXA).The interplay between their intracellular levels (anti-apoptotic proteins versus apoptotic proteins) and their protein-protein interactions decides the fate of cell survival.Therefore, altering or blocking these protein-protein interactions (as "PPIs") is one of the strategies used to develop inhibitors.However, targeting PPIs with small molecules can be challenging, as these PPIs mainly occur over relatively large and flat surfaces.Although, there are small pockets between these interfaces (hot spots), which accommodate high-affinity binding interactions of small functional groups and, generally, is the side chain of an amino acid of the other binding partner protein [11].In this context, the cellular expression of Mcl-1 is tightly regulated by its partner proteins, both proapoptotic and apoptotic, which utilize their alpha-helical peptide regions (more specifically, "BH3-helix domain") for selective binding to the pockets of hMcl-1 and inhibit its functions.In the case of hMcl-1 (Figure 1), for example, the BH3 α-helix uses hydrophobic amino acid side chains to bind to hMcl-1 via interactions with four hydrophobic pockets, P 1-4 .Therefore, BH3 helix-based ligand design has been used to develop various inhibitor types, such as oligopeptides, stapled α-helixes, α-/β-peptide foldamers, reverse BH3 (rBH3) peptide, and small molecules, as Mcl-1 inhibitors, but none have attained clinical application so far [12].Developed hMcl-1 inhibitors based on the BH3-helix have been plagued with affinities/activities for other anti-apoptotic proteins of the Bcl-2 family as hMcl-1 shares high structural homology with other proteins (especially Bcl-2 and Bcl-x L ), and therefore most discovered inhibitors against hMcl-1 are either dualor multi-targeting inhibitors [12].
Chemistry 2022, 4, FOR PEER REVIEW 2 intracellular levels (anti-apoptotic proteins versus apoptotic proteins) and their proteinprotein interactions decides the fate of cell survival.Therefore, altering or blocking these protein-protein interactions (as "PPIs") is one of the strategies used to develop inhibitors.However, targeting PPIs with small molecules can be challenging, as these PPIs mainly occur over relatively large and flat surfaces.Although, there are small pockets between these interfaces (hot spots), which accommodate high-affinity binding interactions of small functional groups and, generally, is the side chain of an amino acid of the other binding partner protein [11].In this context, the cellular expression of Mcl-1 is tightly regulated by its partner proteins, both proapoptotic and apoptotic, which utilize their alpha-helical peptide regions (more specifically, "BH3-helix domain") for selective binding to the pockets of hMcl-1 and inhibit its functions.In the case of hMcl-1 (Figure 1), for example, the BH3 α-helix uses hydrophobic amino acid side chains to bind to hMcl-1 via interactions with four hydrophobic pockets, P1-4.Therefore, BH3 helix-based ligand design has been used to develop various inhibitor types, such as oligopeptides, stapled α-helixes, α-/βpeptide foldamers, reverse BH3 (rBH3) peptide, and small molecules, as Mcl-1 inhibitors, but none have attained clinical application so far [12].Developed hMcl-1 inhibitors based on the BH3-helix have been plagued with affinities/activities for other anti-apoptotic proteins of the Bcl-2 family as hMcl-1 shares high structural homology with other proteins (especially Bcl-2 and Bcl-xL), and therefore most discovered inhibitors against hMcl-1 are either dual-or multi-targeting inhibitors [12].
In recent years, seven classes were reported as inhibitors of hMcl-1 (six from marine sources: gymnochrome-F, oxy-polyhalogenated diphenyl ethers, anacardic acids, endiandric acids, maritoclax, and cryptosphaerolide, while one natural compound was from a terrestrial plant source: meiogynin).Their structures are summarized in Section 3.1, and their inhibitory data are provided in Section 3.3.8.To date, no other biophysical studies were performed on these compounds prior to measuring their binding behavior with hMcl-1.Herein, we provide details of a modeling study to gain insight into their binding modes.This includes validation of the computational protocols and correlation of results The crystal structure of hMcl-1 protein bound with a BH3 α-helix peptide provides information on the binding event (PDB id: 2PQK).The structure shows 4 pockets (P 1 , P 2 , P 3 , P 4 ) in the active site, with each hydrophobic amino acid side chain of the BH3 helix (H1, H2, H3, H4) utilizing these pockets.
In recent years, seven classes were reported as inhibitors of hMcl-1 (six from marine sources: gymnochrome-F, oxy-polyhalogenated diphenyl ethers, anacardic acids, endiandric acids, maritoclax, and cryptosphaerolide, while one natural compound was from a terrestrial plant source: meiogynin).Their structures are summarized in Figure 2, and their inhibitory data are provided in Table 7.To date, no other biophysical studies were performed on these compounds prior to measuring their binding behavior with hMcl-1.Herein, we provide details of a modeling study to gain insight into their binding modes.This includes validation of the computational protocols and correlation of results obtained with published, experimentally determined affinity data, as well as with X-ray crystal structural data for hMcl-1 bound to small molecules [13].

Protein Structure Preparation for Modeling
All protein crystal structural coordinates were downloaded from the protein databank (www.rcsb.org).They were further prepared for modeling in CCG MOE 2019.01 by: (a) removal of solvent molecules from their crystal structures, (b) converting any protein downloaded in a dimer to the monomeric state (Mcl-1 protein inhibition/activation takes place through its BH3 domain, which is a monomer), (c) adding any missing atoms of amino acid residues in the loops, and (d) protonation to ensure appropriate ionization and tautomeric forms for amino acid residues of the hMcl-1 protein.
Our study is limited to ligand-based human Mcl-1 (hMcl-1) inhibitor design; therefore, those hMcl-1 with no co-crystallized ligand or sourced from a different organism other than human, or those with truncated or mutated protein structures, were ignored.In some cases, the NMR-based 15 N-HSQC-based mouse Mcl-1 (mMcl-1) perturbation studies were performed, which were taken into account to determine the homologous amino acids in hMcl-1.

Ligand Structure Preparation
The ligands were collected from various crystal structures, as downloaded from the RCSB website (see in Table 1).In the case of naturally occurring Mcl-1 inhibitors, their structures were built using the MOE builder tool.In the case of cryptosphaerolide, as the absolute configuration of three stereogenic centers remains unknown, all eight possible stereoisomers were constructed.The generated models of the structures were energyminimized in MOE (MMFF94x forcefield) with a gradient of 0.00001.All energy-minimized structures were compiled into a single library ligand database file, formatted into .mdbformat as required for MOE docking.In the case of AutoDock, structures in .molformat were converted to .pdbqtfiles as required for docking using this software, while .molformat was used directly for docking using VlifeDock.

Superposition and Docking Protocols and Analysis
The atoms of the protein backbones, of all nine X-ray co-crystal structures of hMcl-1, which had small molecules bound, were superposed in MOE version 2019.01, and RMSD values were calculated.In our previous study, we compared different placement methods of MOE (triangle matcher, alpha triangle, and alpha PMI methods) for their suitability for various types of scaffolds.During our studies, we compared the resulted docking energies of eight inhibitors (whose native X-ray conformation with fucosidase enzyme was known) docked on the fucosidase enzyme using three different placement methods (triangle matcher, alpha triangle, and alpha PMI methods), and we found that the correlation factor (R 2 = 0.908) for the triangle matcher was comparatively higher.Therefore, we used the triangle matcher method as a reasonable placement method and GB/VI as a reasonable scoring algorithm [14] for various types of ligands [15].Here, again, we used the triangle matcher placement method (100,000 poses with retention time 3000 s and 10,000 returned poses) and refinement using rigid structure with poses = 1000 (gradient = 0.0001 and iteration = 5000).The triangle matcher placement method generates the inhibitor-protein poses by aligning the ligand atoms in triplet on triplets of alpha spheres in a more methodical way.Using this method, the docked conformation of each ligand was compared with their crystal conformation to calculate RMSD values.Docked conformations were arranged into clusters of similar binding conformations from the 500 lowest energy poses (RMSD value ≤ 2.0 Å).The data for the lowest energy-minimized pose from the largest cluster are displayed in the tables and shown in the various figures.A similar approach was used to analyze output from docking attained from the AutoDock and VlifeDock methods.Genetic algorithm-based grid docking in AutoDock was used.Thus, hMcl-1 protein crystal and ligand structures were saved in .pdbqtformat for AutoDock 1.5.6 [16,17].The AutoGrid was mapped around the co-crystallized ligand using the following dimensions, 60 Å × 60 Å × 60 Å (representing the three coordinates X, Y, and Z of the grid), with 0.400 Å spacing.The docking parameters were set up as follows: (a) a maximum of 25 conformers were considered for each ligand, (b) population size was kept to 250 with the individual conformers initialized randomly, (c) maximum numbers of generations were 27,000, and (d) the number of evaluations was set as 2,500,000.Docking with VlifeDock [17] involved the use of the in-house placement method (a Genetic Algorithm) and the scoring was calculated by VlifeDock's GRIP function.A grid was initially generated for docking (60 Å × 60 Å × 60 Å).This docking method generates a population of initial poses, which evolves into the optimal binding mode.After docking, the interactions (H-bonding, hydrophobic, ionic, and aromatic π-π), measurements, and labeling were carried out with the Accelrys discovery studio visualizer.Studying interactions are not only limited to understand the biophysical properties of small molecules with macromolecules but also could be used for choosing better drug delivery materials for natural compounds (which have limited physiochemical properties), such as nanoencapsulation of medicinal essential oils [18].

Validation of Molecular Modeling Protocols
Nineteen structures (see Table 1) of coordinates for hMcl-1 co-crystallized with ligands from the protein databank (www.rcsb.org)were retrieved.Although these crystal structures belong to the same species (Homo sapiens), their quality varied.Validation metrics are shown in Table 1, which are independent of the ligand.There were hMcl-1 crystal structures that contained peptide ligands (Table 1, entries 10-19), and these showed typical interactions with amino acids of the cavity, including those in the pockets: P 1 , P 2 , P 3 , and P 4 (Figure 1).We did not proceed with further study for those that had peptide ligands, which interacted with P 1 -P 4 .Small-molecule inhibitors (SMIs) co-crystallized (Table 1, entries 1-9) with hMcl-1 have a wider diversity of scaffolds compared to peptides, with a greater variety of binding modes.It was decided that these nine SMI structures would be used to generate binding modes for the natural products.[23,24], and its synthetic analog, the maritoclax structure, can be divided into ring-A (magenta) and ring-B (red) for easy presentation.6. Gossypol [25,26] and hMcl-1 inhibitor molecule 1 (MIM1) [27], a high-throughput screening compound, is not a natural product but shares features with Gossypol.7. Cryptosphaerolide from a marine-derived ascomycete related to the genus Cryptosphaeria [28] and derivatives.8. Meiogynin (obtained from ethyl acetate extract of the bark of Meiogyne cylindrocarpa) and synthetic derivatives (8a-c) [29,30].[23,24], and its synthetic analog, the maritoclax structure, can be divided into ring-A (magenta) and ring-B (red) for easy presentation.6. Gossypol [25,26] and hMcl-1 inhibitor molecule 1 (MIM1) [27], a high-throughput screening compound, is not a natural product but shares features with Gossypol.7. Cryptosphaerolide from a marine-derived ascomycete related to the genus Cryptosphaeria [28] and derivatives.8. Meiogynin (obtained from ethyl acetate extract of the bark of Meiogyne cylindrocarpa) and synthetic derivatives (8a-c) [29,30].
Having selected nine co-crystal structures of hMcl-1 with SMIs, the multiple receptor conformations-based approach (MRC) [34] was implemented, using the ensemble of conformations from the crystal structures (Figure 3a).Thus, the backbone atoms of these conformers were superimposed (entries 1-9, Table 1, excluding the small molecules), with subsequent calculation of the various RMSD values (Table 2).The hMcl-1 backbone conformation was found to be conserved based on RMSD values between 0.52 and 1.13 Å, and these values were considered acceptable given the total number of amino acids (157) in the protein backbones (Table 2).The average RMSD values are also reported.Next, the multiple ligand conformations-based approach (MLC) was applied, which involved docking all SMIs from co-crystals to each of the nine hMcl-1 structures downloaded from the protein databank.The conformations of the docked structures were then compared to that in the X-ray crystal structure and RMSD values calculated (Table 3), and the average RMSD values for the various ligands ranged from 0.879 to 1.390 Å with standard deviations from 0.116 to 0.331 Å.The structure 3WIX was selected as the PDB structure on which to proceed with further docking experiments as it accommodated co-crystal ligands to other PDBs with the lowest RMSD average.Through this approach, we also obtained structural information that helped to identify the active and passive amino acids in the cavity (mapped in Figure 3b) involved in the binding of all nine ligands to the 3WIX PDB (the superimposed ligands are shown in Figure 3c).to proceed with further docking experiments as it accommodated co-crystal ligands to other PDBs with the lowest RMSD average.Through this approach, we also obtained structural information that helped to identify the active and passive amino acids in the cavity (mapped in Figure 3b) involved in the binding of all nine ligands to the 3WIX PDB (the superimposed ligands are shown in Figure 3c).
. Table 2. MRC approach: superposition of backbone of multiple hMcl-1 from various crystal structures and their calculated RMSD values (Å ).

Evaluation of the Precision and Accuracy of Docking Methods
To evaluate the precision and accuracy of docking methods, docking scores of a set of hMcl-1 inhibitors (as shown in Figure 4) were obtained and compared with their corresponding Ki values [35].These hMcl-1 inhibitors were chosen as they were prepared by the same group (Fesik's group) and the binding assay conditions were uniform for these inhibitors.Secondly, they have all been reported in co-crystal structures with hMcl-1, providing the experimental basis to validate docking methods.We employed three docking methods, with each ligand being docked into its own crystal structure (self-docking), and the results are summarized in Table 4 [13].The self-docking [17] generally showed that the use of MOE [36] afforded the lowest RMSD values for its docked poses when compared with the crystal structure-bound ligand conformation in hMcl-1, followed by AutoDock [37] and then VLifeDock [38].
inhibitors.Secondly, they have all been reported in co-crystal structures with hMcl-1, providing the experimental basis to validate docking methods.We employed three docking methods, with each ligand being docked into its own crystal structure (self-docking), and the results are summarized in Table 4 [13].The self-docking [17] generally showed that the use of MOE [36] afforded the lowest RMSD values for its docked poses when compared with the crystal structure-bound ligand conformation in hMcl-1, followed by AutoDock [37] and then VLifeDock [38].

Figure 4.
Fesik's group's ligands, which were obtained from co-crystal structures with hMcl-1 (PDB: 4HW2 and 4HW3 [39], PDB: 5FDO and 5FDR [40], PDB: 5IF4 [41], PDB: 6BW2 and 6BW8 [42]).Fesik's group's ligands, which were obtained from co-crystal structures with hMcl-1 (PDB: 4HW2 and 4HW3 [39], PDB: 5FDO and 5FDR [40], PDB: 5IF4 [41], PDB: 6BW2 and 6BW8 [42]).To further assess the accuracy of the docking methods, iterative self-docking was performed using 3WIX [43], where the co-crystal ligand acid] of 3WIX was self-docked five times into its binding site, and each time it was compared with the orientations of its own conformation in the co-crystal structure by RMSD, as shown in Table 5.The choice of using 3WIX for iterative self-docking was based on its co-crystallized ligand's higher selectivity (above 60-fold) to Mcl-1 protein compared to Bcl-x L (Mcl-1 IC 50 = 0.54 µM; Bcl-x L IC 50 > 30 µM).The docking was commenced from random conformations of the ligand, and this, as well as affecting accuracy, can affect precision [44].The precision of docking is also a concern in drug design, with precision defined by comparing RMSD values of poses obtained after repeated docking of the same ligand.The measurement of reproducibility of accuracy (NRMSE) was determined from the RMSD between the docked pose that is closest to the RMSD to that of the co-crystallized ligand of 3WIX, as shown in Table 5.The "RMSE" value is the difference between the docked pose that has the highest RMSD value and the lowest RMSD value of the docked pose and is thus a measure of the precision of the docking.These values are provided in Table 5, and the poses are shown in Figure 5.The MOE triangle matcher docking method was the most accurate (NRMSE = 86%) compared to AutoDock (genetic algorithm, 55%) and VLifeDock (GRIP method-based docking: 36%).The MOE method was also the most precise as it had the lowest RMSE value.We also observed comparatively fast, accurately bound conformation prediction from VLifeDock when compared with AutoDock.However, AutoDock was more precise in generating docking poses [45] but less accurate in reproducing the co-crystal pose as compared to VLifeDock.

Docking of Naturally Occurring hMcl-1 Inhibitors and Related Compounds
After evaluating the docking methods, screening of seven classes of naturally occurring hMcl-1 inhibitors was carried out.All these ligands had shown affinity in the low micromolar range.

Gymnochrome-F
To determine the binding conformation of Gymnochrome-F (1), we employed docking to 3WIX using MOE, AutoDock, and VlifeDock.Similar-pose ensemble clustering [51] confirmed that the attained docking poses from these three different methods were all within RMSD 2 Å .This low RMSD value attained from these methods could be expected due to the high rigidity possessed by the fused heteroaromatic rings in the ligand structure, which limits the number of binding conformations of Gymnochrome F. The docking of Gymnochrome-F showed its main interactions with His224 (π-π interaction, 3.81 Å) and H-bond donor and acceptor interactions with Asn260 (2.14 and2.23 Å ), as shown in Figure 6.

Docking of Naturally Occurring hMcl-1 Inhibitors and Related Compounds
After evaluating the docking methods, screening of seven classes of naturally occurring hMcl-1 inhibitors was carried out.All these ligands had shown affinity in the low micromolar range.

Gymnochrome-F
To determine the binding conformation of Gymnochrome-F (1), we employed docking to 3WIX using MOE, AutoDock, and VlifeDock.Similar-pose ensemble clustering [51] confirmed that the attained docking poses from these three different methods were all within RMSD 2 Å.This low RMSD value attained from these methods could be expected due to the high rigidity possessed by the fused heteroaromatic rings in the ligand structure, which limits the number of binding conformations of Gymnochrome F. The docking of Gymnochrome-F showed its main interactions with His224 (π-π interaction, 3.81 Å) and H-bond donor and acceptor interactions with Asn260 (2.14 and2.23 Å), as shown in Figure 6.

Sponge-Derived Oxy-Polyhalogenated Diphenyl Ethers 2
In case of the diphenyl ether derivatives, we observed various bound conformations, all with similar docking energies (within 0.2 kcal/mol).The smaller structure of these natural product derivatives enabled them to adopt a wide number of potential bound conformations in the larger active site of hMcl-1, as shown in Figure 7.The structures were in proximity to Phe270 but with no interaction with, for example, Arg 263, that interacted with anacardic acid's salicylic acid residue, for example. .

Sponge-Derived Oxy-Polyhalogenated Diphenyl Ethers 2
In case of the diphenyl ether derivatives, we observed various bound conformations, all with similar docking energies (within 0.2 kcal/mol).The smaller structure of these natural product derivatives enabled them to adopt a wide number of potential bound conformations in the larger active site of hMcl-1, as shown in Figure 7.The structures were in proximity to Phe270 but with no interaction with, for example, Arg 263, that interacted with anacardic acid's salicylic acid residue, for example.

Sponge-Derived Oxy-Polyhalogenated Diphenyl Ethers 2
In case of the diphenyl ether derivatives, we observed various bound conformations, all with similar docking energies (within 0.2 kcal/mol).The smaller structure of these natural product derivatives enabled them to adopt a wide number of potential bound conformations in the larger active site of hMcl-1, as shown in Figure 7.The structures were in proximity to Phe270 but with no interaction with, for example, Arg 263, that interacted with anacardic acid's salicylic acid residue, for example. .

Anacardic Acid Derivatives (3a-g)
In case of anacardic acids, 3a showed close binding poses from all three docking methods (1.17 Å RMSD).The docking predicted close interactions of polar phenolic head and carboxylic acid with Arg263 (2.36 and 2.21 Å ) and an intramolecular interaction (1.81 Å ), as shown in Figure 8.

Anacardic Acid Derivatives (3a-g)
In case of anacardic acids, 3a showed close binding poses from all three docking methods (1.17 Å RMSD).The docking predicted close interactions of polar phenolic head and carboxylic acid with Arg263 (2.36 and 2.21 Å) and an intramolecular interaction (1.81 Å), as shown in Figure 8.The extra length in the hydrocarbon chain of 3b compared to 3a projected the salicylic acid substructure of 3b more towards the NWGR domain of hMcl-1, where it interacted with Arg263 as a H-bond acceptor (2.0 and 2.14 Å ), which also supports its improved hMcl-1 affinity (IC50 = 5.8 µM), compared to that of 3a (IC50 = 17.7 µM), as shown in Figure 9.The extra length in the hydrocarbon chain of 3b compared to 3a projected the salicylic acid substructure of 3b more towards the NWGR domain of hMcl-1, where it interacted with Arg263 as a H-bond acceptor (2.0 and 2.14 Å), which also supports its improved hMcl-1 affinity (IC 50 = 5.8 µM), compared to that of 3a (IC 50 = 17.7 µM), as shown in Figure 9.The extra length in the hydrocarbon chain of 3b compared to 3a projected the salicylic acid substructure of 3b more towards the NWGR domain of hMcl-1, where it interacted with Arg263 as a H-bond acceptor (2.0 and 2.14 Å ), which also supports its improved hMcl-1 affinity (IC50 = 5.8 µM), compared to that of 3a (IC50 = 17.7 µM), as shown in Figure 9.In comparison to 3a and 3b, shortening of a hydrocarbon chain and incorporation of a phenyl ring, as in 3c, further improved the affinity for hMcl-1 (IC 50 = 0.6 µM), which could be understood with its additional π-π interaction with Phe270, which was missing in previously mentioned anacardic acid derivatives, as shown in Figure 10.Clustering ensemble of docking conformations was found within RMSD of 1.67 Å and showed a π-π interaction with Phe270 (3.68 Å) along with Arg263 (2.14, 2.38, and 2.82 Å).
Chemistry 2022, 4, FOR PEER REVIEW 13 In comparison to 3a and 3b, shortening of a hydrocarbon chain and incorporation of a phenyl ring, as in 3c, further improved the affinity for hMcl-1 (IC50 = 0.6 µM), which could be understood with its additional π-π interaction with Phe270, which was missing in previously mentioned anacardic acid derivatives, as shown in Figure 10.Clustering ensemble of docking conformations was found within RMSD of 1.67 Å and showed a ππ interaction with Phe270 (3.68 Å ) along with Arg263 (2.14, 2.38, and 2.82 Å ).The docking of 3d, 3e, and 3f showed almost identical docking interactions with Arg263, as shown in Figure 11, which is supported by their similar docking energies (−8.34, −8.47, and −8.22 kcal/mol, respectively).These anacardic acids' binding modes and binding energies are similar, which pointed out that the hydrocarbon chain length of 11-13 carbons plays a less significant role in influencing the binding of these molecules to hMcl-1.The docking of 3d, 3e, and 3f showed almost identical docking interactions with Arg263, as shown in Figure 11, which is supported by their similar docking energies (−8.34, −8.47, and −8.22 kcal/mol, respectively).These anacardic acids' binding modes and binding energies are similar, which pointed out that the hydrocarbon chain length of 11-13 carbons plays a less significant role in influencing the binding of these molecules to hMcl-1.The docking of 3d, 3e, and 3f showed almost identical docking interactions with Arg263, as shown in Figure 11, which is supported by their similar docking energies (−8.34, −8.47, and −8.22 kcal/mol, respectively).These anacardic acids' binding modes and binding energies are similar, which pointed out that the hydrocarbon chain length of 11-13 carbons plays a less significant role in influencing the binding of these molecules to hMcl-1.In case of 3g (IC50 = 1.2 µM, ΔG = −8.00kcal.mol−1 ), when compared to 3e (IC50 = 0.2 µM, ΔG = −8.47kcal.mol−1 ) and 3f (IC50 = 0.2 µM, ΔG = −8.22kcal.mol−1 ), it has a longer hydrocarbon chain which is also unsaturated, which might lead to a reduced affinity, which is most likely to be due to its steric effect (as shown in Figure 12).In case of 3g (IC 50 = 1.2 µM, ∆G = −8.00kcal.mol−1 ), when compared to 3e (IC 50 = 0.2 µM, ∆G = −8.47kcal.mol−1 ) and 3f (IC 50 = 0.2 µM, ∆G = −8.22kcal.mol−1 ), it has a longer hydrocarbon chain which is also unsaturated, which might lead to a reduced affinity, which is most likely to be due to its steric effect (as shown in Figure 12).

Endiandric Acid Analogues (4a-d)
The docking showed that 4a (fused ring with carboxylic acid) binding conformation utilizes similar regions of hMcl-1 as 3c (salicylic acid derivative), which could possibly be based on similarity in their molecular structures (4a and 3c).However, structurally, 4a has only a COOH group in its ring, while 3c has both the acidic phenolic OH and COOH groups, where these polar functional groups were found to consistently interact with Arg263 of hMcl-1.The constrained fused ring in 4a might prevent its adoption of the optimal binding conformation for interaction with Arg263, and therefore could account for the observation of a weak H-bond interaction with Arg263 (3.51 Å) (ΔG = −7.15kcal.mol−1 ) when compared to the 3c docking pose (has three H-bond acceptor interactions with Arg263: 2.14, 2.38, and 2.82 Å, ΔG = −7.98 kcal.mol−1 ) (see Figure 10).Additionally, docking of 4a showed a π-π interaction with Phe270 (3.51 Å), as shown in Figure 13.Similar interactions were also seen with 4b (Ki = 13 µ M), where it displayed an improved geometry for H-bonding, as evaluated by its closer proximity to Arg263 (3.02 Å ) and its π-π interaction with Phe270 (3.50 Å ), as shown in Figure 14.The latter explains the marginal improvement in the in vitro hMcl-1 affinity and docked binding energy, as indicated in Section 3.3.8.

Endiandric Acid Analogues (4a-d)
The docking showed that 4a (fused ring with carboxylic acid) binding conformation utilizes similar regions of hMcl-1 as 3c (salicylic acid derivative), which could possibly be based on similarity in their molecular structures (4a and 3c).However, structurally, 4a has only a COOH group in its ring, while 3c has both the acidic phenolic OH and COOH groups, where these polar functional groups were found to consistently interact with Arg263 of hMcl-1.The constrained fused ring in 4a might prevent its adoption of the optimal binding conformation for interaction with Arg263, and therefore could account for the observation of a weak H-bond interaction with Arg263 (3.51 Å) (∆G = −7.15kcal.mol−1 ) when compared to the 3c docking pose (has three H-bond acceptor interactions with Arg263: 2.14, 2.38, and 2.82 Å, ∆G = −7.98 kcal.mol−1 ) (see Figure 10).Additionally, docking of 4a showed a π-π interaction with Phe270 (3.51 Å), as shown in Figure 13.Similar interactions were also seen with 4b (Ki = 13 µM), where it displayed an improved geometry for Hbonding, as evaluated by its closer proximity to Arg263 (3.02 Å) and its π-π interaction with Phe270 (3.50 Å), as shown in Figure 14.The latter explains the marginal improvement in the in vitro hMcl-1 affinity and docked binding energy, as indicated in Table 7.The comparatively shorter analog in the endiandric acid class (4c) has an additional phenolic polar head and therefore shows an additional new H-bond donor interaction with the peptide backbone of Leu267, and it has proximity with Phe270 (π-π interaction, 3.85 Å ) when compared to longer analogs (4a and 4b), as shown in Figure 15; therefore, this explains why 4c achieved better in vitro hMcl-1 affinity than the latter two, as shown in the table of Section 3.3.8.The comparatively shorter analog in the endiandric acid class (4c) has an additional phenolic polar head and therefore shows an additional new H-bond donor interaction with the peptide backbone of Leu267, and it has proximity with Phe270 (π-π interaction, 3.85 Å ) when compared to longer analogs (4a and 4b), as shown in Figure 15; therefore, this explains why 4c achieved better in vitro hMcl-1 affinity than the latter two, as shown in the table of Section 3.3.8.The comparatively shorter analog in the endiandric acid class (4c) has an additional phenolic polar head and therefore shows an additional new H-bond donor interaction with the peptide backbone of Leu267, and it has proximity with Phe270 (π-π interaction, 3.85 Å) when compared to longer analogs (4a and 4b), as shown in Figure 15; therefore, this explains why 4c achieved better in vitro hMcl-1 affinity than the latter two, as shown in Table 7. Ligand 4d has a different stereochemistry in its fused cyclic ring substructure with respect to the other analogs of this series, which leads to its COOH functionality being more inclined towards Arg263 (3.29 and 3.90 Å ) as compared to a similar analog (such as 4a), as shown in Figure 16, which may account for its improved binding affinity (5.2 µ M) over 4a (14.3 µ M).Ligand 4d has a different stereochemistry in its fused cyclic ring substructure with respect to the other analogs of this series, which leads to its COOH functionality being more inclined towards Arg263 (3.29 and 3.90 Å) as compared to a similar analog (such as 4a), as shown in Figure 16, which may account for its improved binding affinity (5.2 µM) over 4a (14.3 µM).Ligand 4d has a different stereochemistry in its fused cyclic ring substructure with respect to the other analogs of this series, which leads to its COOH functionality being more inclined towards Arg263 (3.29 and 3.90 Å ) as compared to a similar analog (such as 4a), as shown in Figure 16, which may account for its improved binding affinity (5.2 µ M) over 4a (14.3 µ M).

Marinopyrrole Analogs
To study the interaction of marinopyrrole-A (maritoclax, 5a), Doi et al. [23] utilized the 1 H-15 N HSQC spectrum of 15 N-labeled mMcl-1 data, and docking with mouse Mcl-1 (mMcl-1).Their study revealed some of the important molecular interactions with mMcl-1 protein: (a) the chloro group in one of the pyrrole rings of maritoclax was facing in towards the binding groove of mMcl-1, (b) the phenol group of the same pyrrole ring had a H-bond with Gly308, (c) the phenol group of another pyrrole group had a H-bond with Thr247, and (d) the carbonyl group had a H-bond with Asn204.Later, we utilized comparative sequence alignment of hMcl-1 with mMcl-1 to identify the corresponding amino acids of hMcl-1 with respect to those of mMcl-1 (as shown Table 6).Furthermore, our study revealed those amino acids of hMcl-1 which had a preference to interact with marinopyrrole-A, with key species-specific differences: (a) identified two point-mutations, G222-D241 and V262-I281, and (b) revealed the hydrophobic character of the cavity of mMcl-1 compared to hMcl-1 (15 out of 25 in mMcl-1 and 13 out of 22 in hMcl-1 are hydrophobic amino acids) (see Table 6).Table 6.Comparison between mMcl-1 and hMcl-1 active site residues *.

Corresponding Amino Acids
* The enlisted amino acids of mMcl-1 showed chemical shift perturbation ( 1 H-15 N HSQC spectrum) during the addition of maritoclax (5a).The amino acids in black color were detected through changes in the 1 H-15 N HSQC spectrum of 15 N-labeled mMcl-1.The un-bolded and bolded amino acids of mMcl-1 in black color showed a 0.05 to 0.08 ppm chemical shift (δ) change, respectively.However, the amino acids shown in green were undetected due to a large signal intensity change.Furthermore, we used sequence alignment to identify the corresponding amino acids of hMcl-1.Red colored amino acids in hMcl-1 represent the species-based point mutations among rat and human Mcl-1 proteins.
Highly substituted non-symmetric N-heterocycles are difficult to synthesis and, therefore scaffolds derived from nature with such structural specificity are of prime importance [52].Maritoclax is one such example.Although maritoclax (5a) showed similar bound conformation as that reported for mMcl-1 in our docking studies with hMcl-1, key amino acids residues' interactions were found missing.To gain a better understanding, the hMcl-1 structure (PDB code: 3WIX) was superimposed with the mMcl-1 (PDB ID: 2JM6), where the RMSD value was found as 6.91 Å.This high deviation provides the basis to explain differences observed in binding conformations of these proteins with maritoclax (5a).In our observation, maritoclax (5a) is a dimer which contains two pyrroles: ring-A and ring-B (Figure 2).The docking studies of maritoclax (5a) with hMcl-1 showed that the phenyl substituent on ring-A displays a π-π interaction with His224 and phenolic OH makes a H-bond donor interaction with Thr266 (1.51 Å) (Figure 17), whereas ring-B showed an intramolecular H-bond interaction between NH-pyrrole and OH-phenyl and the carbonyl functionality showed two H-bond acceptor interactions with the NH-functional group of the guanidine moiety of Arg263 (2.06 and 2.64 Å).It has been observed that the maritoclax (5a) is in close proximity to Asp256, Met231, Phe254, Leu267, His224, and Phe228 found in the cavity, which were also reported by Doi et al. [23].
that the maritoclax (5a) is in close proximity to Asp256, Met231, Phe254, Leu267, His224, and Phe228 found in the cavity, which were also reported by Doi et al. [23].Later, Cheng et al. chemically modified maritoclax into cyclic derivatives which were evaluated for their Mcl-1 potencies.This work led to cyclic maritoclax (5b) with improved drug-like properties and a more potent hMcl-1 IC50 value (4.3 µ M) compared to the parent maritoclax (8.9 µ M) [53].The correlation of improved drug-like properties and selective potency with cyclization of inhibitor is based on the fact that, in drug design, cyclization restricts the conformational flexibility of a ligand to improve its selectivity and potency for a given physiological target.However, docking evaluated herein failed to explain the differences in binding energies between maritoclax and its cyclic derivative (Figure 18), which shows an intrinsic deficiency of docking algorithms to possess physicochemical parameters of ligand metrics.Later, Cheng et al. chemically modified maritoclax into cyclic derivatives which were evaluated for their Mcl-1 potencies.This work led to cyclic maritoclax (5b) with improved drug-like properties and a more potent hMcl-1 IC 50 value (4.3 µM) compared to the parent maritoclax (8.9 µM) [53].The correlation of improved drug-like properties and selective potency with cyclization of inhibitor is based on the fact that, in drug design, cyclization restricts the conformational flexibility of a ligand to improve its selectivity and potency for a given physiological target.However, docking evaluated herein failed to explain the differences in binding energies between maritoclax and its cyclic derivative (Figure 18), which shows an intrinsic deficiency of docking algorithms to possess physicochemical parameters of ligand metrics.that the maritoclax (5a) is in close proximity to Asp256, Met231, Phe254, Leu267, His224, and Phe228 found in the cavity, which were also reported by Doi et al. [23].Later, Cheng et al. chemically modified maritoclax into cyclic derivatives which were evaluated for their Mcl-1 potencies.This work led to cyclic maritoclax (5b) with improved drug-like properties and a more potent hMcl-1 IC50 value (4.3 µ M) compared to the parent maritoclax (8.9 µ M) [53].The correlation of improved drug-like properties and selective potency with cyclization of inhibitor is based on the fact that, in drug design, cyclization restricts the conformational flexibility of a ligand to improve its selectivity and potency for a given physiological target.However, docking evaluated herein failed to explain the differences in binding energies between maritoclax and its cyclic derivative (Figure 18), which shows an intrinsic deficiency of docking algorithms to possess physicochemical parameters of ligand metrics.

MIM1 (Mcl-1 Inhibitor Molecule 1)
Similar to maritoclax (5), NMR-assisted docking was performed to identify the binding pose of MIM1 (6a) by Cohen et al. [27].The ensemble clustering of docked poses from MOE, AutoDock, and VlifeDock had RMSD values of 1.07 Å, 1.61 Å, and 2.02 Å, respectively (see Figure 19).The attained docking pose in our study was compared with the previous reported binding of 6a, as was studied with hMcl-1 protein.In a previous report, it had been reported that the 6a had a cyclohexyl group hydrophobic interaction with Val216 and Leu213 [27], which was not seen in our case as both residues were away from cyclohexyl functionality (>5.5 Å), which is too long even for a weak hydrophobic interaction.However, other interactions were found similar to the previous study: (a) imidazole ring of His224 with a CH-π interaction (4.19 Å) with the methylated thiazole, and (b) H-bond donor interaction (2.49Å) with ortho-substituted -OH of pyrogallol functionality.On the other hand, the guanidine moiety of the side chain of Arg263 and the alcoholic side chain (-OH) of Thr266 showed a H-bond acceptor interaction (2.47 and 2.06 Å) with para-substituted and ortho-substituted -OH of pyrogallol functionality.We also observed that the Asp256 is in close proximity to Arg263 (2.63 Å), and could be better-suited to contribute to binding the pyrogallol group via H-bonding [27].Similar to maritoclax (5), NMR-assisted docking was performed to identify the binding pose of MIM1 (6a) by Cohen et al. [27].The ensemble clustering of docked poses from MOE, AutoDock, and VlifeDock had RMSD values of 1.07 Å, 1.61 Å, and 2.02 Å, respectively (see Figure 19).The attained docking pose in our study was compared with the previous reported binding of 6a, as was studied with hMcl-1 protein.In a previous report, it had been reported that the 6a had a cyclohexyl group hydrophobic interaction with Val216 and Leu213 [27], which was not seen in our case as both residues were away from cyclohexyl functionality (>5.5 Å ), which is too long even for a weak hydrophobic interaction.However, other interactions were found similar to the previous study: (a) imidazole ring of His224 with a CH-π interaction (4.19 Å ) with the methylated thiazole, and (b) Hbond donor interaction (2.49Å ) with ortho-substituted -OH of pyrogallol functionality.On the other hand, the guanidine moiety of the side chain of Arg263 and the alcoholic side chain (-OH) of Thr266 showed a H-bond acceptor interaction (2.47 and 2.06 Å ) with parasubstituted and ortho-substituted -OH of pyrogallol functionality.We also observed that the Asp256 is in close proximity to Arg263 (2.63 Å ), and could be better-suited to contribute to binding the pyrogallol group via H-bonding [27].and acyclic (7b).The absolute stereochemistry was reported for the cyclic portion but not for the acyclic substructure.Therefore, we generated all eight possible stereoisomers and docked them.These stereoisomers were later ranked based on their docking scores (see Table 7).As shown in Figure 20, docking of 7 showed that the cyclic substructure of all eight stereoisomers (coded in different colors in the figure) does not fit into the cavity and is projected out from the binding cavity.This explains why Fenical and co-workers [28] did not observe any hMcl-1 inhibition activity for 7a.

Cryptosphaerolide
Cryptosphaerolide (7, TR-FRET, Ki = 11.4 µM) contains 2 substructures: cyclic (7a) and acyclic (7b).The absolute stereochemistry was reported for the cyclic portion but not for the acyclic substructure.Therefore, we generated all eight possible stereoisomers and docked them.These stereoisomers were later ranked based on their docking scores (see Table 7).As shown in Figure 20, docking of 7 showed that the cyclic substructure of all eight stereoisomers (coded in different colors in the figure) does not fit into the cavity and is projected out from the binding cavity.This explains why Fenical and co-workers [28] did not observe any hMcl-1 inhibition activity for 7a.

Meiogynin-Derived hMcl-1 Inhibitors (8a-c)
Meiogynin-A (8)-derived compounds have shown a wide range of activity agains anti-apoptotic proteins of the Bcl-2 family [29].Docking predicted similar binding inter actions for meiogynins, except in the case of their parent molecule 8 (Ki = 5.2 µ M), which showed H-acceptor interactions with the imidazole side chain of His224.Meiogynin de rivatives showed π-π interactions with Phe270.Previous docking studies of 8, which showed proximal binding with Arg100 and Tyr195, were also observed in our docking experiments (meiogynin-A was found within 5 Å of these residues).However, the dock  Meiogynin-A (8)-derived compounds have shown a wide range of activity against antiapoptotic proteins of the Bcl-2 family [29].Docking predicted similar binding interactions for meiogynins, except in the case of their parent molecule 8 (Ki = 5.2 µM), which showed H-acceptor interactions with the imidazole side chain of His224.Meiogynin derivatives showed π-π interactions with Phe270.Previous docking studies of 8, which showed proximal binding with Arg100 and Tyr195, were also observed in our docking experiments (meiogynin-A was found within 5 Å of these residues).However, the docking generated distinctive poses compared to those previously disclosed [54,55] for meiogynins, as shown in Figures 21-24

Conclusions
A number of the small-sized hMcl-1 inhibitors were reported in recent years [12], but the majority of them have failed due to either their off-target binding or therapeutic efficacy, or both, which leads to adverse effects such as drug-drug interactions (e.g., alteration of eicosanoid metabolism [56]) and drug-tissue toxicity (such as an increase in the progression of neurodegeneration of neurons in geriatric patients).As nearly 70% of marketed drugs are derived from natural products, this certainly proves natural products as a good source of exploration of new drug molecules, and a number have been reported as inhibitors of Mcl-1 in recent years.Our studies in this paper focused on gaining insights into the binding mode of the seven natural product class inhibitors of hMcl-1 and a limited number of derivatives that have been synthesized.Selection of a crystal structure from the protein databank was followed by evaluation of three independent scoring functions and three docking methods as to their suitability to study docking to hMcl-1.This was carried out by self-docking using various ligands with experimentally measured hMcl-1 affinities (Ki values) and for which crystal structure data were available showing their binding modes.Based on these three scoring methods, we further disclosed the binding conformations for the most potent naturally derived compounds of these seven classes.However, we also used the clustering of similar docked poses within 2.0 Å and chose the largest cluster to identify the most likely binding conformation.We have also repetitively seen that the NWGR (asparagine-tryptophan-glycine-arginine) domain of hMcl-1 is consistently within 4.5 Å of the conformations of the natural products and closely related derivatives.This domain also has an active role in the binding of BH3 helix of proapoptotic proteins.The attained binding conformations were compared with the experimental hMcl-1 affinity data (in vitro IC 50 ) to validate or evaluate the reliability of the adopted biophysical methodologies [13].Some of the highlights of this manuscript are regarding how subtle structural changes alter the binding conformation of a molecule as well as respective affinities.The docking with cryptosphaerolide and its isomers, for example, showed that its polycyclic fragment has little if any interaction with the binding groove, whereas docking with maritoclax showed its cyclic framework having more interaction with the groove.Furthermore, various non-covalent interactions, such as hydrogen bonding, Van der Waal's forces, and electrostatic forces, are commonly observed in the binding of these natural products, which finds key intrinsic molecular features required for a structure to achieve low micromolar Mcl-1 inhibitory activity.Such insights with the methodologies adopted in this manuscript could be helpful and appealing to various types of chemists (such as natural product chemists, chemical biologists, medicinal chemists, or biophysicists) in choosing how to deploy resources when looking to develop natural products in drug discovery programs.

Figure 3 .
Figure 3. (a) Superposition of conformers of hMcl-1 found in crystal structures with small-molecule inhibitors.(b) H-donor and acceptor cavity in hMcl-1 identified by MLC (docking of ligands to 3WIX).(c) Superimposed binding poses of ligands obtained from co-crystal structures docked to 3WIX.

Figure 3 .
Figure 3. (a) Superposition of conformers of hMcl-1 found in crystal structures with small-molecule inhibitors.(b) H-donor and acceptor cavity in hMcl-1 identified by MLC (docking of ligands to 3WIX).(c) Superimposed binding poses of ligands obtained from co-crystal structures docked to 3WIX.

Figure 5 .
Figure 5. Self-docking of 3WIX PDB to its ligand using three methods (together), along with similarpose ensemble clustering [49].(a) The left-hand side (LHS) figure shows the ionized solvent-exposed surface mapping of hMcl-1, where ligands were fitted to the binding site.(b) The right-hand side (RHS) shows the similar-pose clustering [50,51]: showing the following overlayed structures-cocrystal-bound conformation of the ligand (in blue), MOE-based docking conformation (in brown), AutoDock-based docking conformation (in green), and GRIP-based docking conformation (in yellow).

Figure 5 .
Figure 5. Self-docking of 3WIX PDB to its ligand using three methods (together), along with similarpose ensemble clustering [49].(a) The left-hand side (LHS) figure shows the ionized solvent-exposed surface mapping of hMcl-1, where ligands were fitted to the binding site.(b) The right-hand side (RHS) shows the similar-pose clustering [50,51]: showing the following overlayed structures-cocrystal-bound conformation of the ligand (in blue), MOE-based docking conformation (in brown), AutoDock-based docking conformation (in green), and GRIP-based docking conformation (in yellow).

Chemistry 2022, 4 ,Figure 6 .
Figure 6.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of Gymnochrome-F from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity, and all are within 1.08 Å RMSD.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of Gymnochrome-F with hMcl-1.

Figure 6 .
Figure 6.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of Gymnochrome-F from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity, and all are within 1.08 Å RMSD.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of Gymnochrome-F with hMcl-1.

Chemistry 2022, 4 , 12 Figure 8 .
Figure 8.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 3a from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3a with hMcl-1.

Figure 8 .
Figure 8.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 3a from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3a with hMcl-1.
the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3a with hMcl-1.

Figure 9 .
Figure 9. (a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 3b from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3b.

Figure 9 .
Figure 9. (a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 3b from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3b.

Figure 10 .
Figure 10.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 3c from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3c.

Figure 10 .
Figure 10.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 3c from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3c.
cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3c.

Figure 12 .
Figure 12.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 3g from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3g.

Figure 12 .
Figure 12.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 3g from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 3g.

Figure 13 .
Figure 13.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4a from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4a.

Figure 14 .
Figure 14.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4b from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4b.

Figure 13 .Figure 13 .
Figure 13.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4a from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4a.

Figure 14 .
Figure 14.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4b from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4b.

Figure 14 .
Figure 14.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4b from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4b.

Figure 15 .
Figure 15.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4c from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4c.

Figure 16 .
Figure 16.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4d from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4d.

Figure 15 .
Figure 15.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4c from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4c.

Chemistry 2022, 4 ,Figure 15 .
Figure 15.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4c from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4c.

Figure 16 .
Figure 16.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4d from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4d.

Figure 16 .
Figure 16.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 4d from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 4d.

Figure 17 .
Figure 17.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of marinopyrrole from MOE (in blue), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods, MOE (in brown), AutoDock (in green), and VlifeDock (in yellow).(c) Binding mode of marinopyrrole.

Figure 18 .
Figure 18.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of cyclic marinopyrrole (5b) from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of cyclic marinopyrrole (5b).

Figure 17 .
Figure 17.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of marinopyrrole from MOE (in blue), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods, MOE (in brown), AutoDock (in green), and VlifeDock (in yellow).(c) Binding mode of marinopyrrole.

Figure 17 .
Figure 17.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of marinopyrrole from MOE (in blue), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods, MOE (in brown), AutoDock (in green), and VlifeDock (in yellow).(c) Binding mode of marinopyrrole.

Figure 18 .
Figure 18.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of cyclic marinopyrrole (5b) from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of cyclic marinopyrrole (5b).

Figure 18 .
Figure 18.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of cyclic marinopyrrole (5b) from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of cyclic marinopyrrole (5b).

Figure 19 .
Figure 19.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of MIM1 from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses attained from the three different docking methods.(c) Binding mode of MIM1.3.3.7.Cryptosphaerolide Cryptosphaerolide (7, TR-FRET, Ki = 11.4 µ M) contains 2 substructures: cyclic (7a)and acyclic (7b).The absolute stereochemistry was reported for the cyclic portion but not for the acyclic substructure.Therefore, we generated all eight possible stereoisomers and docked them.These stereoisomers were later ranked based on their docking scores (see Table7).As shown in Figure20, docking of 7 showed that the cyclic substructure of all eight stereoisomers (coded in different colors in the figure) does not fit into the cavity and is projected out from the binding cavity.This explains why Fenical and co-workers[28] did not observe any hMcl-1 inhibition activity for 7a.

Figure 19 .
Figure 19.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of MIM1 from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses attained from the three different docking methods.(c) Binding mode of MIM1.

Figure 20 .
Figure 20.Eight stereoisomeric analogs of cryptosphaerolide can be possible based on specific ste reochemistry at three chiral centers, therefore all eight isomers were docked.The cyclic substruc tures of all 8 isomers were not involved in binding (as seen in the binding mode in the figure fo any of the eight stereoisomers (b: on the left-hand side) in the figure with hMcl-1).

Figure 20 .
Figure 20.Eight stereoisomeric analogs of cryptosphaerolide can be possible based on specific stereochemistry at three chiral centers, therefore all eight isomers were docked.The cyclic substructures of all 8 isomers were not involved in binding (as seen in the binding mode in the figure for any of the eight stereoisomers (b: on the left-hand side) in the figure with hMcl-1). .

Figure 21 .
Figure 21.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8 from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8.

Figure 22 .
Figure 22.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8a from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8a.

Figure 21 .Figure 21 .
Figure 21.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8 from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8.

Figure 22 .
Figure 22.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8a from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8a.

Figure 22 .
Figure 22.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8a from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8a.

Figure 23 .
Figure 23.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8b from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8b.

Figure 24 .
Figure 24.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8c from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8c.

Figure 23 . 23 Figure 23 .
Figure 23.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8b from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8b.

Figure 24 .
Figure 24.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8c from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8c.

Figure 24 .
Figure 24.(a) The ionized solvent-exposed surface mapping of hMcl-1 along with the docked poses of 8c from MOE (in brown), AutoDock (in green), and VlifeDock (in yellow) show their fitting to the cavity.(b) Similar-pose ensemble clustering of all three poses was attained from the three different docking methods.(c) Binding mode of 8c.

Table 1 .
Parameters of co-crystal structures from the protein databank where hMcl-1 binds to a ligand [

Table 2 .
MRC approach: superposition of backbone of multiple hMcl-1 from various crystal structures and their calculated RMSD values (Å).

Table 4 .
Comparison of scores from self-docking using various methods with hMcl-1 inhibitor Ki values.
a Scores in kcal/mol.b RMSD in Å.

Table 5 .
Reproducibility of self-docking experiments, performed five times, on 3WIX and its co-crystal ligand *.

Table 7 .
Comparison of Ki's on hMcl-1 inhibition with their respective cumulative docking energies of naturally derived molecules.

IC 50 in µM Bcl-2 a /Bcl-x L b Ki/IC 50 = µM
Ki of 11.4 µM, so all isomers were listed here, d IC 50 values, e Ki values, NR not reported, f kcal/mol.