Protein tyrosine phosphatases (PTP) play important roles in the pathogenesis of many diseases. The fact that no PTP inhibitors have reached the market so far has raised many questions about their druggability. In this study, the active sites of 17 PTPs were characterized and assessed for its ability to bind drug-like molecules. Consequently, PTPs were classified according to their druggability scores into four main categories. Only four members showed intermediate to very druggable pocket; interestingly, the rest of them exhibited poor druggability.

Particularly focusing on PTP1B, we also demonstrated the influence of several factors on the druggability of PTP active site. For instance, the open conformation showed better druggability than the closed conformation, while the tight-bound water molecules appeared to have minimal effect on the PTP1B druggability.

Finally, the allosteric site of PTP1B was found to exhibit superior druggability compared to the catalytic pocket. This analysis can prove useful in the discovery of new PTP inhibitors by assisting researchers in predicting hit rates from high throughput or virtual screening and saving unnecessary cost, time, and efforts via prioritizing PTP targets according to their predicted druggability. Introduction Protein phosphorylation and dephosphorylation are key regulatory mechanisms in living organisms; aberrant protein phosphorylation is highly correlated with many human diseases ranging from diabetes to many forms of cancers. These two reactions are catalyzed by two major enzymes categories – protein kinases that catalyze the phosphorylation process and phosphatases that are responsible for the dephosphorylation process. Although the two processes are of even importance as molecular switches, the scientific world paid more attention toward the kinases; it was not until recently when protein phosphatases started to gain proper attention in the medicinal chemistry field., Protein phosphatases have been grouped into two major categories – serine/threonine phosphatases and tyrosine phosphatases.

The latter group, that is protein tyrosine phosphatases (PTPs), has gained considerable amount of attention for their pivotal non-redundant role in many cellular pathways and in the pathogenesis of many diseases. Thus, the growing number of human diseases associated with aberrant PTP function has brought them in the spotlight as a key drug target. The first bona fide piece of evidence came with the PTP1B knockout mouse which clearly demonstrated that PTP1B acts as a negative regulator of insulin signaling., Since then, a large list of key PTPs associated with different human abnormalities has emerged as possible drug targets. The strong relationship between PTP and different human diseases made them hot targets for pharmaceutical companies trying to design suitable inhibitors. The design of PTP inhibitors has been rendered by the difficulty in identifying cell permeable compound with sufficient selectivity; so far, no PTP inhibitors have reached the market.– displays examples of different PTP inhibitors that have inhibitory potency ranging between low and medium nanomolar levels. Obviously, the shown compounds carry multiple charges on their structures, so as to become capable of binding with the PTP positively charged pocket; this can partially explain the cell permeability issues generally accompanied with such inhibitors.

Examples of potent PTP inhibitors. SiteMap has been able to help researchers in assessing the active site of important targets, for example, the bromodomain acetyl-lysine proteins were found to have significantly different druggability, with a median druggability score (Dscore) ranging from 0.38 to 1.08. In another study, enoyl-acyl carrier protein reductases were all predicted to be druggable by SiteMap, with average Dscore ranging from 1.05 to 1.30; nonetheless, the size of their active site varied significantly, with the average pocket size ranging from 101 to 284. In this study, we assessed the active site of PTPs, studying their druggability and other physical characteristics using SiteMap, and then analyzed/classified them according to the resultant scores. In this work, we also studied the effect of different structural aspects on the PTP druggability, that is, loop flexibility, solvation effect, and allosteric site, which can be proven valuable for designing new PTP inhibitors. Selection and preparation of PTP crystal structures The PTPs considered in this study are known to be involved in many diseases., All PTP crystal structures were obtained from the Protein Data Bank.

All PTP structures considered in the study were solved by X-crystallography and have no mutated residues in their active site; otherwise, they were discarded. PTP structures with covalently bound inhibitors were not considered in the final list. Our resultant PTP crystal structure library finally contained a total number of 116 crystal structures that belong to 17 PTP enzymes. Preparation of PTP crystal structures The crystal structures obtained from PDB were categorized according to the PTP member they belong to and then subcategorized according to whether or not they have a ligand molecule in their active site. Subsequently, each crystal structure was processed manually to keep one protein chain.

All water molecules were removed (where applicable); all other solvent molecules and co-crystallized ligands were also removed. Afterward, each crystal structure was processed automatically via the protein preparation wizard in the MOE software. This wizard prepares a PTP crystal structure by checking the protein part for any missing atoms, residues, or loops; accordingly, it conducts all required corrections. Next, the MOE protein preparation wizard assigns the protonation state of each ionizable side chain (using the Protonate 3D algorithm). (2) where n is the number of site points, e is the enclosure score, and p is the hydrophilic score. As shown in the equations, both scoring functions use the same terms but with different coefficients so that highly polar active sites are more penalized by the Dscore function than SiteScore.

Thus, a polar pocket may score low in Dscore, although a high score can be achieved when the other algorithm, SiteScore, is employed. Such a pocket is more likely to encourage the binding of highly polar ligands that do not display druglikeness. Another parameter calculated by SiteMap is the size of the active site; the number of points that constitute an examined pocket is a measure of the size of the site. Approximately, two to three site points correspond to each atom of the bound ligand, including hydrogen atoms. Analyzing scores obtained by SiteMap SiteMap has been validated by a set of 538 complexes and was able to accurately identify the ligand binding site in 86% of the cases.

Pockets with a SiteScore value of ≥0.8 are identified as a possible binding pocket by SiteMap and those with a SiteScore of ≥1.0 are classified as binding sites with particular importance, tight-binding sites. A slightly different classification system used by Vidler et al will be used in this work where the Dscore function is employed to assess the PTP binding pocket for their druggability. This system categorizes proteins into four main classes: very druggable (Dscore ≥1.0), druggable (Dscore =0.8–1.0), intermediate (Dscore =0.7–0.8), and difficult (Dscore ≤0.7). Druggability assessment and binding site analysis The term “druggable” was defined previously by Hopkins and Groom as a protein that is able, or predicted to be able, to bind drug-like molecules. This concept is considered in this study to assess the catalytic pocket of different members from the PTP family.

The fact that no PTP inhibitors have reached the market, so far, highlights the difficulty that medicinal chemists face to design new PTP inhibitors with adequate cell permeability. In order to attain the full picture, the catalytic pocket of PTP enzymes was studied and assessed for druggability. Searching the protein data bank resulted in obtaining 116 crystal structures that belong to 17 members of the PTP family, all of which were assessed for their druggability using SiteMap. Interestingly, the results show that PTPs had a very wide range of median Dscores (0.25–1.20; ). The mycobacterial phosphatase, MPtpB, showed the highest Dscore value among all tested PTPs, while JSP-1 exhibited the poorest druggability.

As shown in, PTPs are not well distributed into these four classes, as only one member was found under the “very druggable” category (MPtpB: Dscore =1.20) and one member under the “drug-gable” category (GLEPP-1: Dscore 0.89). Only two PTPs were found in the next category “intermediate” (ie, SHP-1 and PTPε), with Dscore values ranging between 0.78 and 0.71. The rest of the PTP members, including the mostly targeted enzyme PTP1B, were classified as “difficult”. Median values of druggability score (Dscore) and other SiteMap factors obtained by the 17 pathogenic PTPs studied in this work Looking at a representative example of each class, differences in the active site of the PTPs can be readily detected.

As shown in and, MPtpB showed a large (pocket size 196 spheres) and well-defined cavity (enclosure score 0.83) with a low philic score (0.7). A similar pocket was also shown by GLEPP-1 but with a smaller pocket size (140 spheres), a less well-defined cavity (enclosure score 0.73), and higher philic score (1.5), in comparison with MPtpB. Although the “moderately druggable” class representative PTPε showed a less hydrophilic pocket (philic score 1.3) than GLEPP-1, its active site was observed to have a less well-defined cavity (enclosure 0.67) with a significantly smaller pocket size (59 spheres). The least druggable member among the PTP family, JSP-1, exhibited a very shallow cavity (enclosure 0.55) with a very small pocket size (=24 spheres), in addition to high hydrophilic nature (philic score =1.6), as shown in. The binding site of four PTP classes exemplified by the surfaces of representative example.

It is worth noting that some of the PTP members were classified as difficult targets, although they demonstrated reasonably well-defined cavities in their active pocket. For instance, PTP1B and CDC25B have a medium pocket size (60–70 spheres) with a comparable enclosure score (0.7) to the intermediate druggable member PTPε , which led all of them to score well by SiteScore (0.81–0.9). However, the former PTPs seem to have 31% more polar pockets (philic factor =1.7) than PTPε, which prevented them from obtaining good druggability scores (median Dscore 0.6). In other words, such PTPs have binding pockets that are good enough to bind ligand molecules but are not good enough to particularly bind ligands with high druglikeness.

Classifying PTP1B via SiteMap as difficult to target by drug-like molecules may further explain why PTP1B inhibitors usually suffer from cell permeability problems. Since Dscore is calculated based on size, enclosure, and hydrophilicity of a tested pocket, it would be interesting to see which of these three factors has the greatest correlation with the obtained Dscore values of the assessed PTPs. Shows that the pocket size factor has a direct and a proportional correlation with the PTP druggability scores. In contrast, pocket enclosure demonstrates a proportional correlation with Dscore; however, this correlation looks less significant compared to pocket size.

Environment

Illustrates a strong correlation between the hydrophilicity and druggability factors; it can be readily detected that when the pocket hydrophilicity increases, the druggability decreases. Hence it can be concluded that the most influential factors on the PTP druggability are the size and the hydrophilicity of their catalytic pocket. The effect of pocket conformation on PTP druggability Proteins are known to be flexible enough to change their binding site conformation according to their pocket status – the apo form or ligand-bound form. Moreover, their binding pocket can change its conformation in accordance with the shape and size of bound ligands, so ligands with different physical and chemical characteristics can still bind.

Thus, it was just sensible to separate those structures in the apo form from those having ligands in their binding pocket, and then reassess them for their druggability. As shown in, the apo and ligand-bound forms were only seen in five PTPs, four of which exhibited different Dscore and pocket size values when the two forms are compared. Three of these PTPs (ie, PTP1B, YopH, and MPtpA) have their apo form more druggable than the ligand-bound form, and only CDC25B showed the opposite. Median values for Dscore and pocket size (number of spheres) obtained by the apo form and ligand-bound form of the 17 PTPs studied here shows examples of CDC25B apo protein and ligand-bound forms aligned on each other.

Moe Molecular Operating Environment

The apo form of CDC25B showed as low Dscore as 0.33, while the ligand-bound form has more than double druggability score (median Dscore =0.76; ). The main reason behind having such druggability variation in the CDC25B active site seems to be due to the side chain flexibility of two arginine residues, Arg479 and Arg482. In the apo protein form, the two residues adopt certain conformations making the binding pocket really small (=40 spheres), preventing it from having an adequate Dscore value. However, when a ligand binds to the enzyme, the two arginine residues rotate away, channeling the active site with an adjacent pocket so that to get large enough (pocket size =83 spheres) to bind drug-like molecules. Consequently, the ligand-bound form of CDC25B may be individually classified as “intermediate” rather than “difficult” (where the apo form is listed in).

The active site of two crystal structures of CDC25B aligned on each other. The other PTP that had druggability differences between its apo form and ligand-bound form is the YopH. In contrast to CDC25B, the YopH apo structure showed a better druggability score (Dscore =0.68) when compared to its ligand-bound form (median Dscore =0.45). Shows all YopH crystal structures aligned on each other.

It can be observed that when YopH is in the unbound form, the WPD-loop adopts an open conformation (pocket size =94 spheres). However, binding of a ligand causes the WPD-loop to swing down to adopt a closed conformation, making the YopH binding pocket significantly smaller in size (median pocket size =37 spheres) and less able to achieve good druggability score. This means YopH would have even less chance to bind drug-like molecules if the catalytic pocket deviated from the apo conformation, giving emphasis to its druggability issue. Active site conformational effect on the PTP1B druggability. Lists the results obtained by SiteMap for the open and closed conformations of PTP1B.

It can be observed that the open conformation was able to score a greater druggability score (Dscore =0.79), that almost touched the lower border of the “druggable” class, when compared to the closed conformation (Dscore =0.61, classified as “difficult”). However, the open and closed conformations of PTP1B scored well by SiteScore (0.9). This means both conformations are equally capable of binding small organic molecules, but the closed conformation lacks the same capability to bind drug-like molecules (the open conformation looks superior in this context). The open structure showed a 16% increase in the pocket size and a 15% decrease in the philic score compared to the closed pocket , although the latter exhibited a 7% greater enclosure than the open conformation. Hence, the pocket size and hydrophilicity seem to be the main responsible factors for the PTP1B open conformation superiority over the closed one. Median values of Dscore and other SiteMap factors obtained by the open and closed structures of PTP1B In fact, druggability data obtained here come in line with previous results obtained from our seeding experiments on PTP1B. In that study, the open conformation was able to retrieve a larger number of known PTP1B inhibitors (that obey Lipinski’s rule of five), in comparison with the closed conformation.

That was explained at the time by the fact that the PTP1B open conformation has a more spacious binding pocket than the closed conformation, and therefore, a greater range of ligands can be readily accommodated in the active site of the open conformation. It can also be added now that, as mentioned previously, the open conformation demonstrates a relatively less polar surface area than the closed conformation, which should, theoretically, assist and facilitate the binding of less polar ligands (that exhibit greater druglikeness) to the open pocket. Practically, so far, most of the open-conformation-bound inhibitors still possess a significant amount of polarity that prevents them from having a good cell membrane permeability profile.

More studies are thus needed in this area to prove the druggability of the PTP1B open conformation. The effect of water on PTP druggability Different-ordered water molecules buried just behind the WPD loop have been observed in the closed PTP1B conformation.– Out of the 82 PTP1B structures considered in this study, the highly ordered water molecules were seen in as much as 43 structures, acting as a hydrogen bonding bridge between the co-crystallized ligand and the corresponding amino acid. Therefore, the ordered water molecules represent an important aspect in the PTP1B drug discovery, and hence it was important to study whether these solvent molecules have an influence on the druggability characteristics of the PTP1B active site. Shows the Dscore values of the solvated and unsolvated forms of the PTP1B active site. It was found that the solvated PTP1B structures have less druggable catalytic pocket by 7% than the unsolvated structure. This appeared to result mainly from the pocket size and enclosure factors as the unsolvated form showed a bigger pocket and a greater enclosure score by 14% and 6%, respectively, when compared to the solvated form.

Surprisingly, although the solvated pocket contains more polar sites (ie, water molecules) than the unsolvated form, the former pocket scored a 7% lower philic score (=1.69; ) when compared to the unsolvated one (philic score =1.82). It seems that these conserved water molecules are covering an even more polar spots in the active site, making the PTP1B pocket less hydrophilic in nature. Overall, it can be concluded that water does not appear to have a great influence on the PTP1B druggability as both forms were not able to escape from the “difficult” class. It is noteworthy that, some of the known PTP1B inhibitors, such as isothiazolidinone, demonstrated high ability to displace all ordered solvent molecules from the PTP1B active site., Thus, water should be handled with caution when designing new PTP inhibitors. The location of catalytic and allosteric pockets on the PTP1B structure along with their co-crystallized ligands. Shows the data obtained from SiteMap for both PTP1B pockets.

Obviously, the allosteric site is much more druggable than the catalytic pocket, scoring as high Dscore as 0.97, which is 52% greater than the catalytic pocket Dscore (Dscore =0.64). Although both pockets had comparable SiteScore (which means that both of them can bind organic molecules), they were not able to score similarly in Dscore. As shown in, the active pocket scored greater philic score (=1.7) by more than threefold compared to the allosteric pocket (philic score =0.5), making polarity the supreme contributor to druggability variations between the two PTP1B pockets. Comparison of druggability data obtained by the catalytic pocket allosteric pockets of PTP1B It is worth noting that this pocket seems to be not the only allosteric site on the PTP1B structure.

A recently discovered non-competitive inhibitor, MSI-1436 , showed an inhibition activity in the low micromolar level with a unique binding to PTP1B. This binding was found to partially involve the aforementioned allosteric site, in addition to a distinct pocket located within the 20 amino acids at the C-terminus of the catalytic domain. MSI-1436 is a promising PTP1B inhibitor as it successfully passed Phase I clinical trials where it showed good pharmacodynamic and pharmacokinetic profiles with minimal side effect.

However, MSI-1436 was not advanced further to Phase II trials due to limited financial resources. In summary, combining the successful story of MSI-1436 along with the druggability findings obtained in this study, allosteric pocket can be considered as a viable target for medicinal chemists. Conclusion Druggability assessment of 17 PTPs has revealed that only two of the tested PTPs seem to be druggable. MPtpB and GLEPP-1 were found to be the only druggable among the PTP family as they showed large hydrophobic active site.

The rest of the PTPs were ranging from intermediate to undruggable because of their hydrophilic and/or shallow catalytic pocket. Gintama 2 subtitle. Focusing on PTP1B, several factors affecting the active site druggability were studied.

Our results showed that the open conformation of the PTP active site has a better druggability than the closed form, while the tight bound water molecule has negligible effect on the PTP1B druggability. Finally, the PTP1B allosteric site demonstrated a significantly improved druggability compared to the catalytic pocket, which suggests that allosteric modulation may become an important pillar for the future of PTP drug discovery. The outcomes of these findings can shed more light in the search for new PTP inhibitors and can assist researchers to avoid extra coast and effort.