Motivation: Recent large-scale omics initiatives have catalogued the somatic alterations of cancer cell line panels along with their pharmacological response to hundreds of compounds. In this study, we have explored these data to advance computational approaches that enable more effective and targeted use of current and future anticancer therapeutics. Results: We modelled the 50% growth inhibition bioassay end-point (GI50) of 17 142 compounds screened against 59 cancer cell lines from the NCI60 panel (941 831 data-points, matrix 93.08% complete) by integrating the chemical and biological (cell line) information. We determine that the protein, gene transcript and miRNA abundance provide the highest predictive signal when modelling the GI50 endpoint, which significantly outperformed the DNA copy-number variation or exome sequencing data (Tukey’s Honestly Significant Difference, P <0.05). We demonstrate that, within the limits of the data, our approach exhibits the ability to both interpolate and extrapolate compound bioactivities to new cell lines and tissues and, although to a lesser extent, to dissimilar compounds. Moreover, our approach outperforms previous models generated on the GDSC dataset. Finally, we determine that in the cases investigated in more detail, the predicted drug-pathway associations and growth inhibition patterns are mostly consistent with the experimental data, which also suggests the possibility of identifying genomic markers of drug sensitivity for novel compounds on novel cell lines. Contact:terez@pasteur.fr; ab454@ac.cam.uk Supplementary information: Supplementary data are available at Bioinformatics online.
The ability of pathogens to cause disease depends on their aptitude to escape the immune system. Type IV pili are extracellular filamentous virulence factors composed of pilin monomers and frequently expressed by bacterial pathogens. As such they are major targets for the host immune system. In the human pathogen Neisseria meningitidis, strains expressing class I pilins contain a genetic recombination system that promotes variation of the pilin sequence and is thought to aid immune escape. However, numerous hypervirulent clinical isolates express class II pilins that lack this property. This raises the question of how they evade immunity targeting type IV pili. As glycosylation is a possible source of antigenic variation it was investigated using top-down mass spectrometry to provide the highest molecular precision on the modified proteins. Unlike class I pilins that carry a single glycan, we found that class II pilins display up to 5 glycosylation sites per monomer on the pilus surface. Swapping of pilin class and genetic background shows that the pilin primary structure determines multisite glycosylation while the genetic background determines the nature of the glycans. Absence of glycosylation in class II pilins affects pilus biogenesis or enhances pilus-dependent aggregation in a strain specific fashion highlighting the extensive functional impact of multisite glycosylation. Finally, molecular modeling shows that glycans cover the surface of class II pilins and strongly decrease antibody access to the polypeptide chain. This strongly supports a model where strains expressing class II pilins evade the immune system by changing their sugar structure rather than pilin primary structure. Overall these results show that sequence invariable class II pilins are cloaked in glycans with extensive functional and immunological consequences. During infection pathogens and their host engage in a series of measures and counter-measures to promote their own survival: pathogens express virulence factors, the immune system targets these surface structures and pathogens modify them to evade detection. Like numerous bacterial pathogens, Neisseria meningitidis express type IV pili, long filamentous adhesive structures composed of pilins. Intriguingly the amino acid sequences of pilins from most hypervirulent strains do not vary, raising the question of how they evade the immune system. This study shows that the pilus structure is completely coated with sugars thus limiting access of antibodies to the pilin polypeptide chain. We propose that multisite glycosylation and thus variation in the type of sugar mediates immune evasion in these strains.
Motivation: Scoring functions provided by the docking software are still a major limiting factor in virtual screening (VS) process to classify compounds. Score analysis of the docking is not able to find out all active compounds. This is due to a bad estimation of the ligand binding energies. Making the assumption that active compounds should have specific contacts with their target to display activity, it would be possible to discriminate active compounds from inactive ones with careful analysis of interatomic contacts between the molecule and the target. However, compounds clustering is very tedious due to the large number of contacts extracted from the different conformations proposed by docking experiments. Results: Structural analysis of docked structures is processed in three steps: (i) a Kohonen self-organizing map (SOM) training phase using drug–protein contact descriptors followed by (ii) an unsupervised cluster analysis and (iii) a Newick file generation for results visualization as a tree. The docking poses are then analysed and classified quickly and automatically by AuPosSOM (Automatic analysis of Poses using SOM). AuPosSOM can be integrated into strategies for VS currently employed. We demonstrate that it is possible to discriminate active compounds from inactive ones using only mean protein contacts' footprints calculation from the multiple conformations given by the docking software. Chemical structure of the compound and key binding residues information are not necessary to find out active molecules. Thus, contact–activity relationship can be employed as a new VS process. Availability: AuPosSOM is available at http://www.aupossom.com. Contact:contact@aupossom.com; gildas.bertho@parisdescartes.fr Supplementary information: Supplementary data are available at Bioinformatics online.
The determination of protein structures satisfying distance constraints is an important problem in structural biology. Whereas the most common method currently employed is simulated annealing, there have been other methods previously proposed in the literature. Most of them, however, are designed to find one solution only. In order to explore exhaustively the feasible conformational space, we propose here an interval Branch-and-Prune algorithm (iBP) to solve the Distance Geometry Problem (DGP) associated to protein structure determination. This algorithm is based on a discretization of the problem obtained by recursively constructing a search space having the structure of a tree, and by verifying whether the generated atomic positions are feasible or not by making use of pruning devices. The pruning devices used here are directly related to features of protein conformations. We described the new algorithm iBP to generate protein conformations satisfying distance constraints, that would potentially allows a systematic exploration of the conformational space. The algorithm iBP has been applied on three α-helical peptides.