Trifloxystrobin blocks the growth of Theileria parasites and is a promising drug to treat Buparvaquone resistance.

Theileria parasites are responsible for devastating cattle diseases, causing major economic losses across Africa and Asia. Theileria spp. stand apart from other apicomplexa parasites by their ability to transform host leukocytes into immortalized, hyperproliferating, invasive cells that rapidly kill infected animals. The emergence of resistance to the theilericidal drug Buparvaquone raises the need for new anti-Theileria drugs. We developed a microscopy-based screen to reposition drugs from the open-access Medicines for Malaria Venture (MMV) Pathogen Box. We show that Trifloxystrobin (MMV688754) selectively kills lymphocytes or macrophages infected with Theileria annulata or Theileria parva parasites. Trifloxystrobin treatment reduced parasite load in vitro as effectively as Buparvaquone, with similar effects on host gene expression, cell proliferation and cell cycle. Trifloxystrobin also inhibited parasite differentiation to merozoites (merogony). Trifloxystrobin inhibition of parasite survival is independent of the parasite TaPin1 prolyl isomerase pathway. Furthermore, modeling studies predicted that Trifloxystrobin and Buparvaquone could interact distinctly with parasite Cytochrome B and we show that Trifloxystrobin was still effective against Buparvaquone-resistant cells harboring TaCytB mutations. Our study suggests that Trifloxystrobin could provide an effective alternative to Buparvaquone treatment and represents a promising candidate for future drug development against Theileria spp.

Structural Impacts of Drug-Resistance Mutations Appearing in HIV-2 Protease.

The use of antiretroviral drugs is accompanied by the emergence of HIV-2 resistances. Thus, it is important to elucidate the mechanisms of resistance to antiretroviral drugs. Here, we propose a structural analysis of 31 drug-resistant mutants of HIV-2 protease (PR2) that is an important target against HIV-2 infection. First, we modeled the structures of each mutant. We then located structural shifts putatively induced by mutations. Finally, we compared wild-type and mutant inhibitor-binding pockets and interfaces to explore the impacts of these induced structural deformations on these two regions. Our results showed that one mutation could induce large structural rearrangements in side-chain and backbone atoms of mutated residue, in its vicinity or further. Structural deformations observed in side-chain atoms are frequent and of greater magnitude, that confirms that to fight drug resistance, interactions with backbone atoms should be favored. We showed that these observed structural deformations modify the conformation, volume, and hydrophobicity of the binding pocket and the composition and size of the PR2 interface. These results suggest that resistance mutations could alter ligand binding by modifying pocket properties and PR2 stability by impacting its interface. Our results reinforce the understanding of the effects of mutations that occurred in PR2 and the different mechanisms of PR2 resistance.

Impacts of drug resistance mutations on the structural asymmetry of the HIV-2 protease.

BACKGROUND: Drug resistance is a severe problem in HIV treatment. HIV protease is a common target for the design of new drugs for treating HIV infection. Previous studies have shown that the crystallographic structures of the HIV-2 protease (PR2) in bound and unbound forms exhibit structural asymmetry that is important for ligand recognition and binding. Here, we investigated the effects of resistance mutations on the structural asymmetry of PR2. Due to the lack of structural data on PR2 mutants, the 3D structures of 30 PR2 mutants of interest have been modeled using an in silico protocol. Structural asymmetry analysis was carried out with an in-house structural-alphabet-based approach. RESULTS: The systematic comparison of the asymmetry of the wild-type structure and a large number of mutants highlighted crucial residues for PR2 structure and function. In addition, our results revealed structural changes induced by PR2 flexibility or resistance mutations. The analysis of the highlighted structural changes showed that some mutations alter protein stability or inhibitor binding. CONCLUSIONS: This work consists of a structural analysis of the impact of a large number of PR2 resistant mutants based on modeled structures. It suggests three possible resistance mechanisms of PR2, in which structural changes induced by resistance mutations lead to modifications in the dimerization interface, ligand recognition or inhibitor binding.

Exploration of the effects of sequence variations between HIV-1 and HIV-2 proteases on their three-dimensional structures.

HIV protease inhibitors (PIs) approved by the FDA (US Food and Drug Administration) are a major class of antiretroviral. HIV-2 protease (PR2) is naturally resistant to most of them as PIs were designed for HIV-1 protease (PR1). In this study, we explored the impact of amino-acid substitutions between PR1 and PR2 on the structure of protease (PR) by comparing the structural variability of 13 regions using 24 PR1 and PR2 structures complexed with diverse ligands. Our analyses confirmed structural rigidity of the catalytic region and highlighted the important role of three regions in the conservation of the catalytic region conformation. Surprisingly, we showed that the flap region, corresponding to a flexible region, exhibits similar conformations in PR1 and PR2. Furthermore, we identified regions exhibiting different conformations in PR1 and PR2, which could be explained by the intrinsic flexibility of these regions, by crystal packing, or by PR1 and PR2 substitutions. Some substitutions induce structural changes in the R2 and R4 regions that could have an impact on the properties of PI-binding site and could thus modify PI binding mode. Substitutions involved in structural changes in the elbow region could alter the flexibility of the PR2 flap regions relative to PR1, and thus play a role in the transition from the semi-open form to the closed form, and have an impact on ligand binding. These results improve the understanding of the impact of sequence variations between PR1 and PR2 on the natural resistance of HIV-2 to commercially available PIs.Communicated by Ramaswamy H. Sarma.

High Impact: The Role of Promiscuous Binding Sites in Polypharmacology.

The literature focuses on drug promiscuity, which is a drug's ability to bind to several targets, because it plays an essential role in polypharmacology. However, little work has been completed regarding binding site promiscuity, even though its properties are now recognized among the key factors that impact drug promiscuity. Here, we quantified and characterized the promiscuity of druggable binding sites from protein-ligand complexes in the high quality Mother Of All Databases while using statistical methods. Most of the sites (80%) exhibited promiscuity, irrespective of the protein class. Nearly half were highly promiscuous and able to interact with various types of ligands. The corresponding pockets were rather large and hydrophobic, with high sulfur atom and aliphatic residue frequencies, but few side chain atoms. Consequently, their interacting ligands can be large, rigid, and weakly hydrophilic. The selective sites that interacted with one ligand type presented less favorable pocket properties for establishing ligand contacts. Thus, their ligands were highly adaptable, small, and hydrophilic. In the dataset, the promiscuity of the site rather than the drug mainly explains the multiple interactions between the drug and target, as most ligand types are dedicated to one site. This underlines the essential contribution of binding site promiscuity to drug promiscuity between different protein classes.

HIV-1 protease, Gag and gp41 baseline substitutions associated with virological response to a PI-based regimen.

OBJECTIVES: To assess, at ART initiation, the impact of baseline substitutions in protease, Gag and gp41 regions on the virological response to a first-line PI-based regimen. PATIENTS AND METHODS: One hundred and fifty-four HIV-infected ART-naive patients initiating a PI-based regimen including darunavir (n = 129) or atazanavir (n = 25) were assessed, including 36 experiencing virological failure (VF). Whole pol, gag and gp41 genes were sequenced at ART baseline by ultra-deep sequencing (UDS) using Illumina® technology. Supervised data-mining analyses were performed to identify mutations associated with virological response. Structural analyses were performed to assess the impact of mutations on protease conformation. RESULTS: UDS was successful in 127, 138 and 134 samples for protease, Gag and gp41, respectively (31% subtype B and 38% CRF02_AG). Overall, T4A and S37T mutations in protease were identified as being associated with VF (P = 0.02 and P = 0.005, respectively). Among CRF02_AG sequences, I72M and E21D mutations were associated with VF (P = 0.03 for both). They all induced some conformational changes of some protease side-chain residues located near mutated residues. In Gag, mutations associated with VF were G62D, N315H and Y441S (P = 0.005, P = 0.007 and P = 0.0003, respectively). All were localized outside Gag cleavage sites (G62D, matrix; N315H, capsid; and Y441S, p1). In gp41, the I270T mutation, localized in the cytoplasmic tail, was associated with VF (P = 0.003), and the I4L mutation, in the fusion peptide, was associated with virological success (P = 0.004). CONCLUSIONS: In this study, new baseline substitutions in Gag, protease and g41, potentially impacting PI-based regimen outcome, were evidenced. Phenotypic analyses are required to confirm their role in the PI-resistance mechanism.

Characterizing the structural variability of HIV-2 protease upon the binding of diverse ligands using a structural alphabet approach.

The HIV-2 protease (PR2) is an important target for designing new drugs against the HIV-2 infection. In this study, we explored the structural backbone variability of all available PR2 structures complexed with various inhibitors using a structural alphabet approach. 77% of PR2 positions are structurally variable, meaning they exhibit different local conformations in PR2 structures. This variability was observed all along the structure, particularly in the elbow and flap regions. A part of these backbone changes observed between the 18 PR2 is induced by intrinsic flexibility, and ligand binding putatively induces others occurring in the binding pocket. These latter changes could be important for PR2 adaptation to diverse ligands and are accompanied by changes outside the binding pocket. In addition, the study of the link between structural variability of the pocket and PR2-ligand interactions allowed us to localize pocket regions important for ligand binding and catalytic function, regions important for ligand recognition that adjust their backbone in response to ligand binding and regions important for the pocket opening and closing that have large intrinsic flexibility. Finally, we suggested that differences in ligand effectiveness for PR2 could be partially explained by different backbone deformations induced by these ligands. To conclude, this study is the first characterization of the PR2 structural variability considering ligand diversity. It provides information about the recognition of PR2 to various ligands and its mechanisms to adapt its local conformation to bound ligands that could help understand the resistance of PR2 to its inhibitors, a major antiretroviral class. Communicated by Ramaswamy H. Sarma.

SAFlex: A structural alphabet extension to integrate protein structural flexibility and missing data information.

In this paper, we describe SAFlex (Structural Alphabet Flexibility), an extension of an existing structural alphabet (HMM-SA), to better explore increasing protein three dimensional structure information by encoding conformations of proteins in case of missing residues or uncertainties. An SA aims to reduce three dimensional conformations of proteins as well as their analysis and comparison complexity by simplifying any conformation in a series of structural letters. Our methodology presents several novelties. Firstly, it can account for the encoding uncertainty by providing a wide range of encoding options: the maximum a posteriori, the marginal posterior distribution, and the effective number of letters at each given position. Secondly, our new algorithm deals with the missing data in the protein structure files (concerning more than 75% of the proteins from the Protein Data Bank) in a rigorous probabilistic framework. Thirdly, SAFlex is able to encode and to build a consensus encoding from different replicates of a single protein such as several homomer chains. This allows localizing structural differences between different chains and detecting structural variability, which is essential for protein flexibility identification. These improvements are illustrated on different proteins, such as the crystal structure of an eukaryotic small heat shock protein. They are promising to explore increasing protein redundancy data and obtain useful quantification of their flexibility.

Exploration of the effect of sequence variations located inside the binding pocket of HIV-1 and HIV-2 proteases.

HIV-2 protease (PR2) is naturally resistant to most FDA (Food and Drug Administration)-approved HIV-1 protease inhibitors (PIs), a major antiretroviral class. In this study, we compared the PR1 and PR2 binding pockets extracted from structures complexed with 12 ligands. The comparison of PR1 and PR2 pocket properties showed that bound PR2 pockets were more hydrophobic with more oxygen atoms and fewer nitrogen atoms than PR1 pockets. The structural comparison of PR1 and PR2 pockets highlighted structural changes induced by their sequence variations and that were consistent with these property changes. Specifically, substitutions at residues 31, 46, and 82 induced structural changes in their main-chain atoms that could affect PI binding in PR2. In addition, the modelling of PR1 mutant structures containing V32I and L76M substitutions revealed a cooperative mechanism leading to structural deformation of flap-residue 45 that could modify PR2 flexibility. Our results suggest that substitutions in the PR1 and PR2 pockets can modify PI binding and flap flexibility, which could underlie PR2 resistance against PIs. These results provide new insights concerning the structural changes induced by PR1 and PR2 pocket variation changes, improving the understanding of the atomic mechanism of PR2 resistance to PIs.

Analysis of the HIV-2 protease's adaptation to various ligands: characterization of backbone asymmetry using a structural alphabet.

The HIV-2 protease (PR2) is a homodimer of 99 residues with asymmetric assembly and binding various ligands. We propose an exhaustive study of the local structural asymmetry between the two monomers of all available PR2 structures complexed with various inhibitors using a structural alphabet approach. On average, PR2 exhibits asymmetry in 31% of its positions-i.e., exhibiting different backbone local conformations in the two monomers. This asymmetry was observed all along its structure, particularly in the elbow and flap regions. We first differentiated structural asymmetry conserved in most PR2 structures from the one specific to some PR2. Then, we explored the origin of the detected asymmetry in PR2. We localized asymmetry that could be induced by PR2's flexibility, allowing transition from the semi-open to closed conformations and the asymmetry potentially induced by ligand binding. This latter could be important for the PR2's adaptation to diverse ligands. Our results highlighted some differences between asymmetry of PR2 bound to darunavir and amprenavir that could explain their differences of affinity. This knowledge is critical for a better description of PR2's recognition and adaptation to various ligands and for a better understanding of the resistance of PR2 to most PR2 inhibitors, a major antiretroviral class.

Exploring the potential of a structural alphabet-based tool for mining multiple target conformations and target flexibility insight.

Protein flexibility is often implied in binding with different partners and is essential for protein function. The growing number of macromolecular structures in the Protein Data Bank entries and their redundancy has become a major source of structural knowledge of the protein universe. The analysis of structural variability through available redundant structures of a target, called multiple target conformations (MTC), obtained using experimental or modeling methods and under different biological conditions or different sources is one way to explore protein flexibility. This analysis is essential to improve the understanding of various mechanisms associated with protein target function and flexibility. In this study, we explored structural variability of three biological targets by analyzing different MTC sets associated with these targets. To facilitate the study of these MTC sets, we have developed an efficient tool, SA-conf, dedicated to capturing and linking the amino acid and local structure variability and analyzing the target structural variability space. The advantage of SA-conf is that it could be applied to divers sets composed of MTCs available in the PDB obtained using NMR and crystallography or homology models. This tool could also be applied to analyze MTC sets obtained by dynamics approaches. Our results showed that SA-conf tool is effective to quantify the structural variability of a MTC set and to localize the structural variable positions and regions of the target. By selecting adapted MTC subsets and comparing their variability detected by SA-conf, we highlighted different sources of target flexibility such as induced by binding partner, by mutation and intrinsic flexibility. Our results support the interest to mine available structures associated with a target using to offer valuable insight into target flexibility and interaction mechanisms. The SA-conf executable script, with a set of pre-compiled binaries are available at http://www.mti.univ-paris-diderot.fr/recherche/plateformes/logiciels.

Statistical Profiling of One Promiscuous Protein Binding Site: Illustrated by Urokinase Catalytic Domain.

While recent literature focuses on drug promiscuity, the characterization of promiscuous binding sites (ability to bind several ligands) remains to be explored. Here, we present a proteochemometric modeling approach to analyze diverse ligands and corresponding multiple binding sub-pockets associated with one promiscuous binding site to characterize protein-ligand recognition. We analyze both geometrical and physicochemical profile correspondences. This approach was applied to examine the well-studied druggable urokinase catalytic domain inhibitor binding site, which results in a large number of complex structures bound to various ligands. This approach emphasizes the importance of jointly characterizing pocket and ligand spaces to explore the impact of ligand diversity on sub-pocket properties and to establish their main profile correspondences. This work supports an interest in mining available 3D holo structures associated with a promiscuous binding site to explore its main protein-ligand recognition tendency.

Cavity Versus Ligand Shape Descriptors: Application to Urokinase Binding Pockets.

We analyzed 78 binding pockets of the human urokinase plasminogen activator (uPA) catalytic domain extracted from a data set of crystallized uPA-ligand complexes. These binding pockets were computed with an original geometric method that does NOT involve any arbitrary parameter, such as cutoff distances, angles, and so on. We measured the deviation from convexity of each pocket shape with the pocket convexity index (PCI). We defined a new pocket descriptor called distributional sphericity coefficient (DISC), which indicates to which extent the protein atoms of a given pocket lie on the surface of a sphere. The DISC values were computed with the freeware PCI. The pocket descriptors and their high correspondences with ligand descriptors are crucial for polypharmacology prediction. We found that the protein heavy atoms lining the urokinases binding pockets are either located on the surface of their convex hull or lie close to this surface. We also found that the radii of the urokinases binding pockets and the radii of their ligands are highly correlated (r = 0.9).

Investigating the Importance of the Pocket-estimation Method in Pocket-based Approaches: An Illustration Using Pocket-ligand Classification.

Small molecules interact with their protein target on surface cavities known as binding pockets. Pocket-based approaches are very useful in all of the phases of drug design. Their first step is estimating the binding pocket based on protein structure. The available pocket-estimation methods produce different pockets for the same target. The aim of this work is to investigate the effects of different pocket-estimation methods on the results of pocket-based approaches. We focused on the effect of three pocket-estimation methods on a pocket-ligand (PL) classification. This pocket-based approach is useful for understanding the correspondence between the pocket and ligand spaces and to develop pharmacological profiling models. We found pocket-estimation methods yield different binding pockets in terms of boundaries and properties. These differences are responsible for the variation in the PL classification results that can have an impact on the detected correspondence between pocket and ligand profiles. Thus, we highlighted the importance of the pocket-estimation method choice in pocket-based approaches.

Structural Isosteres of Phosphate Groups in the Protein Data Bank.

We developed a computational workflow to mine the Protein Data Bank for isosteric replacements that exist in different binding site environments but have not necessarily been identified and exploited in compound design. Taking phosphate groups as examples, the workflow was used to construct 157 data sets, each composed of a reference protein complexed with AMP, ADP, ATP, or pyrophosphate as well other ligands. Phosphate binding sites appear to have a high hydration content and large size, resulting in U-shaped bioactive conformations recurrently found across unrelated protein families. A total of 16 413 replacements were extracted, filtered for a significant structural overlap on phosphate groups, and sorted according to their SMILES codes. In addition to the classical isosteres of phosphate, such as carboxylate, sulfone, or sulfonamide, unexpected replacements that do not conserve charge or polarity, such as aryl, aliphatic, or positively charged groups, were found.

PockDrug-Server: a new web server for predicting pocket druggability on holo and apo proteins.

Predicting protein pocket's ability to bind drug-like molecules with high affinity, i.e. druggability, is of major interest in the target identification phase of drug discovery. Therefore, pocket druggability investigations represent a key step of compound clinical progression projects. Currently computational druggability prediction models are attached to one unique pocket estimation method despite pocket estimation uncertainties. In this paper, we propose 'PockDrug-Server' to predict pocket druggability, efficient on both (i) estimated pockets guided by the ligand proximity (extracted by proximity to a ligand from a holo protein structure) and (ii) estimated pockets based solely on protein structure information (based on amino atoms that form the surface of potential binding cavities). PockDrug-Server provides consistent druggability results using different pocket estimation methods. It is robust with respect to pocket boundary and estimation uncertainties, thus efficient using apo pockets that are challenging to estimate. It clearly distinguishes druggable from less druggable pockets using different estimation methods and outperformed recent druggability models for apo pockets. It can be carried out from one or a set of apo/holo proteins using different pocket estimation methods proposed by our web server or from any pocket previously estimated by the user. PockDrug-Server is publicly available at: http://pockdrug.rpbs.univ-paris-diderot.fr.

PockDrug: A Model for Predicting Pocket Druggability That Overcomes Pocket Estimation Uncertainties.

Predicting protein druggability is a key interest in the target identification phase of drug discovery. Here, we assess the pocket estimation methods' influence on druggability predictions by comparing statistical models constructed from pockets estimated using different pocket estimation methods: a proximity of either 4 or 5.5 Å to a cocrystallized ligand or DoGSite and fpocket estimation methods. We developed PockDrug, a robust pocket druggability model that copes with uncertainties in pocket boundaries. It is based on a linear discriminant analysis from a pool of 52 descriptors combined with a selection of the most stable and efficient models using different pocket estimation methods. PockDrug retains the best combinations of three pocket properties which impact druggability: geometry, hydrophobicity, and aromaticity. It results in an average accuracy of 87.9% ± 4.7% using a test set and exhibits higher accuracy (∼5-10%) than previous studies that used an identical apo set. In conclusion, this study confirms the influence of pocket estimation on pocket druggability prediction and proposes PockDrug as a new model that overcomes pocket estimation variability.

Structure of the prolyl-acyl carrier protein oxidase involved in the biosynthesis of the cyanotoxin anatoxin-a.

Anatoxin-a and homoanatoxin-a are two potent cyanobacterial neurotoxins biosynthesized from L-proline by a short pathway involving polyketide synthases. Proline is first loaded onto AnaD, an acyl carrier protein, and prolyl-AnaD is then oxidized to 1-pyrroline-5-carboxyl-AnaD by a flavoprotein, AnaB. Three polyketide synthases then transform this imine into anatoxin-a or homoanatoxin-a. AnaB was crystallized in its holo form and its three-dimensional structure was determined by X-ray diffraction at 2.8 Šresolution. AnaB is a homotetramer and its fold is very similar to that of the acyl-CoA dehydrogenases (ACADs). The active-site base of AnaB, Glu244, superimposed very well with that of human isovaleryl-CoA dehydrogenase, confirming previous site-directed mutagenesis experiments and mechanistic proposals. The substrate-binding site of AnaB is small and is likely to be fitted for the pyrrolidine ring of proline. However, in contrast to ACADs, which use an electron-transport protein, AnaB uses molecular oxygen as the electron acceptor, as in acyl-CoA oxidases. Calculation of the solvent-accessible surface area around the FAD in AnaB and in several homologues showed that it is significantly larger in AnaB than in its homologues. A protonated histidine near the FAD in AnaB is likely to participate in oxygen activation. Furthermore, an array of water molecules detected in the AnaB structure suggests a possible path for molecular oxygen towards FAD. This is consistent with AnaB being an oxidase rather than a dehydrogenase. The structure of AnaB is the first to be described for a prolyl-ACP oxidase and it will contribute to defining the structural basis responsible for oxygen reactivity in flavoenzymes.

Insights into an original pocket-ligand pair classification: a promising tool for ligand profile prediction.

Pockets are today at the cornerstones of modern drug discovery projects and at the crossroad of several research fields, from structural biology to mathematical modeling. Being able to predict if a small molecule could bind to one or more protein targets or if a protein could bind to some given ligands is very useful for drug discovery endeavors, anticipation of binding to off- and anti-targets. To date, several studies explore such questions from chemogenomic approach to reverse docking methods. Most of these studies have been performed either from the viewpoint of ligands or targets. However it seems valuable to use information from both ligands and target binding pockets. Hence, we present a multivariate approach relating ligand properties with protein pocket properties from the analysis of known ligand-protein interactions. We explored and optimized the pocket-ligand pair space by combining pocket and ligand descriptors using Principal Component Analysis and developed a classification engine on this paired space, revealing five main clusters of pocket-ligand pairs sharing specific and similar structural or physico-chemical properties. These pocket-ligand pair clusters highlight correspondences between pocket and ligand topological and physico-chemical properties and capture relevant information with respect to protein-ligand interactions. Based on these pocket-ligand correspondences, a protocol of prediction of clusters sharing similarity in terms of recognition characteristics is developed for a given pocket-ligand complex and gives high performances. It is then extended to cluster prediction for a given pocket in order to acquire knowledge about its expected ligand profile or to cluster prediction for a given ligand in order to acquire knowledge about its expected pocket profile. This prediction approach shows promising results and could contribute to predict some ligand properties critical for binding to a given pocket, and conversely, some key pocket properties for ligand binding.

Characterization of Arabidopsis calcium-dependent protein kinases: activated or not by calcium?

CDPKs (calcium-dependent protein kinases), which contain both calmodulin-like calcium binding and serine/threonine protein kinase domains, are only present in plants and some protozoans. Upon activation by a stimulus, they transduce the signal through phosphorylation cascades to induce downstream responses, including transcriptional regulation. To understand the functional specificities of CDPKs, 14 Arabidopsis CPKs (CDPKs in plants) representative of the three main subgroups were characterized at the biochemical level, using HA (haemagglutinin)-tagged CPKs expressed in planta. Most of them were partially or mainly associated with membranes, in agreement with acylation predictions. Importantly, CPKs displayed highly variable calcium-dependences for their kinase activities: seven CPKs from subgroups 1 and 2 were clearly sensitive to calcium with different intensities, whereas six CPKs from subgroup 3 exhibited low or no calcium sensitivity to two generic substrates. Interestingly, this apparent calcium-independence correlated with significant alterations in the predicted EF-hands of these kinases, although they all bound calcium. The noticeable exception, CPK25, was calcium-independent owing to the absence of functional EF-hands. Taken together, the results of the present study suggest that calcium binding differentially affects CDPK isoforms that may be activated by distinct molecular mechanisms.