Toward the prediction of class I and II mouse major histocompatibility complex-peptide-binding affinity:in silico bioinformatic step-by-step guide using quantitative structure-activity relationships

Quantitative structure-activity relationship (QSAR) analysis is a cornerstone of modern informatics. Predictive computational models of peptide-major histocompatibility complex (MHC)-binding affinity based on QSAR technology have now become important components of modern computational immunovaccinology. Historically, such approaches have been built around semiqualitative, classification methods, but these are now giving way to quantitative regression methods. We review three methods--a 2D-QSAR additive-partial least squares (PLS) and a 3D-QSAR comparative molecular similarity index analysis (CoMSIA) method--which can identify the sequence dependence of peptide-binding specificity for various class I MHC alleles from the reported binding affinities (IC50) of peptide sets. The third method is an iterative self-consistent (ISC) PLS-based additive method, which is a recently developed extension to the additive method for the affinity prediction of class II peptides. The QSAR methods presented here have established themselves as immunoinformatic techniques complementary to existing methodology, useful in the quantitative prediction of binding affinity: current methods for the in silico identification of T-cell epitopes (which form the basis of many vaccines, diagnostics, and reagents) rely on the accurate computational prediction of peptide-MHC affinity. We have reviewed various human and mouse class I and class II allele models. Studied alleles comprise HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3101, HLA-A*6801, HLA-A*6802, HLA-B*3501, H2-K(k), H2-K(b), H2-D(b) HLA-DRB1*0101, HLA-DRB1*0401, HLA-DRB1*0701, I-A(b), I-A(d), I-A(k), I-A(S), I-E(d), and I-E(k). In this chapter we show a step-by-step guide into predicting the reliability and the resulting models to represent an advance on existing methods. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made are freely available online at the URL http://www.jenner.ac.uk/MHCPred.

Quantitative structure-activity relationships and the prediction of MHC supermotifs

The underlying assumption in quantitative structure–activity relationship (QSAR) methodology is that related chemical structures exhibit related biological activities. We review here two QSAR methods in terms of their applicability for human MHC supermotif definition. Supermotifs are motifs that characterise binding to more than one allele. Supermotif definition is the initial in silico step of epitope-based vaccine design. The first QSAR method we review here—the additive method—is based on the assumption that the binding affinity of a peptide depends on contributions from both amino acids and the interactions between them. The second method is a 3D-QSAR method: comparative molecular similarity indices analysis (CoMSIA). Both methods were applied to 771 peptides binding to 9 HLA alleles. Five of the alleles (A*0201, A* 0202, A*0203, A*0206 and A*6802) belong to the HLA-A2 superfamily and the other four (A*0301, A*1101, A*3101 and A*6801) to the HLA-A3 superfamily. For each superfamily, supermotifs defined by the two QSAR methods agree closely and are supported by many experimental data.

Automated synthesis and evaluation of potential new anti-microbial agents

A series of N1-benzylideneheteroarylcarboxamidrazones was prepared in an automated fashion, and tested against Mycobacterium fortuitum in a rapid screen for antimycobacterial activity. Many of the compounds from this series were also tested against Mycobacterium tuberculosis, and the usefulness as M.fortuitum as a rapid, initial screen for anti-tubercular activity evaluated. Various deletions were made to the N1-benzylideneheteroarylcarboxamidrazone structure in order to establish the minimum structural requirements for activity. The N1-benzylideneheteroarylcarbox-amidrazones were then subjected to molecular modelling studies and their activities against M.fortuitum and M.tuberculosis were analysed using quantitative structure-analysis relationship (QSAR) techniques in the computational package TSAR (Oxford Molecular Ltd.). A set of equations predictive of antimycobacterial activity was hereby obtained. The series of N1-benzylidenehetero-arylcarboxamidrazones was also tested against a multidrug-resistant strain of Staphylococcus aureus (MRSA), followed by a panel of Gram-positive and Gram-negative bacteria, if activity was observed for MRSA. A set of antimycobacterial N1-benzylideneheteroarylcarboxamidrazones was hereby discovered, the best of which had MICs against m. fortuitum in the range 4-8μgml-1 and displayed 94% inhibition of M.tuberculosis at a concentration of 6.25μgml-1. The antimycobacterial activity of these compounds appeared to be specific, since the same compounds were shown to be inactive against other classes of organisms. Compounds which were found to be sufficiently active in any screen were also tested for their toxicity against human mononuclear leucocytes. Polyethylene glycol (PEG) was used as a soluble polymeric support for the synthesis of some fatty acid derivatives, containing an isoxazoline group, which may inhibit mycolic acid synthesis in mycobacteria. Both the PEG-bound products and the cleaved, isolated products themselves were tested against M.fortuitum and some low levels of antimycobacterial activity were observed, which may serve as lead compounds for further studies.

Empirical prediction of peptide octanol-water partition coefficients

Peptides are of great therapeutic potential as vaccines and drugs. Knowledge of physicochemical descriptors, including the partition coefficient P (commonly expressed in logarithm form: logP), is useful for screening out unsuitable molecules and also for the development of predictive Quantitative Structure-Activity Relationships (QSARs). In this paper we develop a new approach to the prediction of LogP values for peptides based on an empirical relationship between global molecular properties and measured physical properties. Our method was successful in terms of peptide prediction (total r2 = 0.641). The final model consisted of 5 physicochemical descriptors (molecular weight, number of single bonds, 2D-VDW volume, 2D-VSA hydrophobic and 2D-VSA polar). The approach is peptide specific and its predictive accuracy was high. Overall, 67% of the peptides were able to be predicted within +/-0.5 log units from the experimental values. Our method thus represents a novel prediction method with proven predictive ability.

On the hydrophobicity of peptides:comparing empirical predictions of peptide log P values

Peptides are of great therapeutic potential as vaccines and drugs. Knowledge of physicochemical descriptors, including the partition coefficient logP, is useful for the development of predictive Quantitative Structure-Activity Relationships (QSARs). We have investigated the accuracy of available programs for the prediction of logP values for peptides with known experimental values obtained from the literature. Eight prediction programs were tested, of which seven programs were fragment-based methods: XLogP, LogKow, PLogP, ACDLogP, AlogP, Interactive Analysis's LogP and MlogP; and one program used a whole molecule approach: QikProp. The predictive accuracy of the programs was assessed using r(2) values, with ALogP being the most effective (r( 2) = 0.822) and MLogP the least (r(2) = 0.090). We also examined three distinct types of peptide structure: blocked, unblocked, and cyclic. For each study (all peptides, blocked, unblocked and cyclic peptides) the performance of programs rated from best to worse is as follows: all peptides - ALogP, QikProp, PLogP, XLogP, IALogP, LogKow, ACDLogP, and MlogP; blocked peptides - PLogP, XLogP, ACDLogP, IALogP, LogKow, QikProp, ALogP, and MLogP; unblocked peptides - QikProp, IALogP, ALogP, ACDLogP, MLogP, XLogP, LogKow and PLogP; cyclic peptides - LogKow, ALogP, XLogP, MLogP, QikProp, ACDLogP, IALogP. In summary, all programs gave better predictions for blocked peptides, while, in general, logP values for cyclic peptides were under-predicted and those of unblocked peptides were over-predicted.

Coupling in silico and in vitro analysis of peptide-MHC binding:a bioinformatic approach enabling prediction of superbinding peptides and anchorless epitopes

The ability to define and manipulate the interaction of peptides with MHC molecules has immense immunological utility, with applications in epitope identification, vaccine design, and immunomodulation. However, the methods currently available for prediction of peptide-MHC binding are far from ideal. We recently described the application of a bioinformatic prediction method based on quantitative structure-affinity relationship methods to peptide-MHC binding. In this study we demonstrate the predictivity and utility of this approach. We determined the binding affinities of a set of 90 nonamer peptides for the MHC class I allele HLA-A*0201 using an in-house, FACS-based, MHC stabilization assay, and from these data we derived an additive quantitative structure-affinity relationship model for peptide interaction with the HLA-A*0201 molecule. Using this model we then designed a series of high affinity HLA-A2-binding peptides. Experimental analysis revealed that all these peptides showed high binding affinities to the HLA-A*0201 molecule, significantly higher than the highest previously recorded. In addition, by the use of systematic substitution at principal anchor positions 2 and 9, we showed that high binding peptides are tolerant to a wide range of nonpreferred amino acids. Our results support a model in which the affinity of peptide binding to MHC is determined by the interactions of amino acids at multiple positions with the MHC molecule and may be enhanced by enthalpic cooperativity between these component interactions.

The use of computational QSAR analysis in the toxicological evaluation of a series of 2-pyridylcarboxamidrazone candidate anti-tuberculosis compounds

A series of N1-benzylidene pyridine-2-carboxamidrazone anti-tuberculosis compounds has been evaluated for their cytotoxicity using human mononuclear leucocytes (MNL) as target cells. All eight compounds were significantly more toxic than dimethyl sulphoxide control and isoniazid (INH) with the exception of a 4-methoxy-3-(2-phenylethyloxy) derivative, which was not significantly different in toxicity compared with INH. The most toxic agent was an ethoxy derivative, followed by 3-nitro, 4-methoxy, dimethylpropyl, 4-methylbenzyloxy, 3-methoxy-4-(-2-phenylethyloxy) and 4-benzyloxy in rank order. In comparison with the effect of selected carboxamidrazone agents on cells alone, the presence of either N-acetyl cysteine (NAC) or glutathione caused a significant reduction in the toxicity of INH, as well as on the 4-benzyloxy derivative, although both increased the toxicity of a 4-N,N-dimethylamino-1-naphthylidene and a 2-t-butylthio derivative. The derivatives from this and three previous studies were subjected to computational analysis in order to derive equations designed to establish quantitative structure activity relationships for these agents. Twenty-five compounds were thus resolved into two groups (1 and 2), which on analysis yielded equations with r2 values in the range 0.65-0.92. Group 1 shares a common mode of toxicity related to hydrophobicity, where cytotoxicity peaked at logP of 3.2, while Group 2 toxicity was strongly related to ionisation potential. The presence of thiols such as NAC and GSH both promoted and attenuated toxicity in selected compounds from Group 1, suggesting that secondary mechanisms of toxicity were operating. These studies will facilitate the design of future low toxicity high activity anti-tubercular carboxamidrazone agents. © 2003 Elsevier Science B.V. All rights reserved.

Structure and regulation of type 2 transglutaminase in relation to its physiological functions and pathological roles

This chapter contains sections titled: Introduction Structure and Regulation Physiologic Functions of TG2 Disruption of TG2 Functions in Pathologic Conditions Perspectives for Pharmacologic Interventions Concluding Comments Acknowledgements References

Infer user interests via link structure regularization

Learning user interests from online social networks helps to better understand user behaviors and provides useful guidance to design user-centric applications. Apart from analyzing users' online content, it is also important to consider users' social connections in the social Web. Graph regularization methods have been widely used in various text mining tasks, which can leverage the graph structure information extracted from data. Previously, graph regularization methods operate under the cluster assumption that nearby nodes are more similar and nodes on the same structure (typically referred to as a cluster or a manifold) are likely to be similar. We argue that learning user interests from complex, sparse, and dynamic social networks should be based on the link structure assumption under which node similarities are evaluated based on the local link structures instead of explicit links between two nodes. We propose a regularization framework based on the relation bipartite graph, which can be constructed from any type of relations. Using Twitter as our case study, we evaluate our proposed framework from social networks built from retweet relations. Both quantitative and qualitative experiments show that our proposed method outperforms a few competitive baselines in learning user interests over a set of predefined topics. It also gives superior results compared to the baselines on retweet prediction and topical authority identification. © 2014 ACM.

Exploring the relation between the qualitative and quantitative uses of the determiner some

This article attempts to repair the neglect of the qualitative uses of some and to suggest an explanation which could cover the full range of usage with this determiner - both quantitative and qualitative - showing how a single underlying meaning, modulated by contextual and pragmatic factors, can give rise to the wide variety of messages expressed by some in actual usage. Both the treatment of some as an existential quantifier and the scalar model which views some as evoking a less-than-expected quantity on a pragmatic scale are shown to be incapable of handling the qualitative uses of this determiner. An original analysis of some and the interaction of its meaning with the defining features of the qualitative uses is proposed, extending the discussion as well to the role of focus and the adverbial modifier quite. The crucial semantic feature of some for the explanation of its capacity to express qualitative readings is argued to be non-identification of a referent assumed to be particular. Under the appropriate conditions, this notion can give rise to qualitative denigration (implying it is not even worth the bother to identify the referent) or qualitative appreciation (implying the referent to be so outstanding that it defies identification). The explanation put forward is also shown to cover some's use as an approximator, thereby enhancing its plausibility even further. © Cambridge University Press 2012.

The structure and dynamics of saxicolous lichen communities

The saxicolous lichen vegetation on Ordovician slate rock at the mouth of the River Dovey, South Merionethshire, Wales was described in relation to several environmental variables which include aspect, slope angle, light intensity, rock porosity, rock microtopography and rock stability. Each of the measured environmental variables was shown to influence the lichen vegetation. A number of groups of species which were characteristic of certain environments were described. The data from the saxicolous lichen communities were analysed using multivariate analysis. Qualitative and quantitative data were ordinated, the qualitative data being easier to interpret ecologically, and site number (which reflects distance from the sea and altitude), rock porosity and light intensity were shown to be important environmental variables. A classification of the data was also carried out. The results of the ordination and classification were combined together and a model constructed which describes saxicolous lichen vegetation. A method which uses the model as an aid to the design and interpretation of field experiments is described. The model is applied to an experiment which investigates the effect on growth of transplanting four saxicolous lichens to different aspects. Growth was inhibited in Physcia orbicularis and Parmelia conspersa on rock surfaces of northwest aspect compared with growth on rock surfaces of southeast aspect. Growth was inhibited in Parmelia glabratula ssp. fuliginosa on rock surfaces of southeast aspect compared with rock surfaces of northwesr aspect. The growth of Parmelia saxatilis was similar at both southeast and northwesr aspects. Growth inhibition or stimulation in thalli of Physcia orbicularis, Parmelia conspersa and Parmelia glabratula ssp. fuliginosa after transplantation was consistent with the predictions of the model while the results for Parmelia saxatilis were not as expected. There was evidence that the frequency of Parmelia conspersa and Parmelia glabratula at a site is related to an effect of the environment on the growth of the thalli. There was also evidence that the frequency of Physcia orbicularis at a site is related to an effect of the environment on the establishment phase of the thalli and for the competitive exclusion of Parmelia saxatilis thalli from southeast facing rock surfaces. The distribution of lichens in relation to height on nine rock surfaces was investigated. It was suggested that the distribution of the lichens was influenced by microclimatic factors which are related to height on the rock, environmental variables which are associated with the rock substratum (e.g. rock porosity and rock microtopography) and by historical factors. The pattern of one crustose and one foliose lichen on four rock surfaces of different aspect and slope was investigated. On the vertically inclined surface the density of small thalli of Buellia aethalea and Parmelia glabratula ssp fuliginosa was correlated with the microtopography of the surface in transects horizontally across the rock surface but not in transects vertically down the rock surface. there were consitent differences in the scale and intensity of pattern horizontally and vertically and also a decrease in the intensity of pattern vertically as the slope of the rock surface decreased. These results were consistent with the suggestion of a gradient of microclimatic factors up the rock. The differences in the scale and intensity of pattern in different size classes in the population were consistent with the changes in pattern with time which have been shown to occur during succession in sand dune and salt marsh vegetation. The relationship between thallus size and height on a rock surface and between the radial growth rate and location of a thallus on a rock surface were investigated. Thalli of Parmelia glabratula ssp. fuliginosa were larger at the top of the rock surface than at the bottom and the data were consistent with the suggestion that the colonisation of the rock surface began at the top and, in time, spread downwards. The radial growth rate of the thalli could not be related to variation in slope, porosity, microtopography or directly to height on the rock but could be related to the horizontal location of the thalli on the rock. These results were consistent with the suggestion that here is a gradient of microclimatic factors across the rock surface which is also modified by height on the rock surface. The succession of lichen communities was described by relating the vegetation to rock porosity, rock microtopography, species diversity and rock stability. An initial stage dominated by crustose lichens leads to communities dominated by crustose, foliose and fruticose species. In the late stages of the succession on some rock surfaces crustose species again become dominant. The occurrence of the climax state and cyclic vegetation change in lichen communities are discussed. A mthod of estimating the age structure of a lichen population by relating thallus size to growth rate is described. The sources of error in the method are discussed in some detail and several refinements suggested to increase the accuracy of the method. The population dynamics of Parmelia glabratula ssp. fuliginosa was investigated by applying life tables to the age structures of eight different populations. The data were consistent with a period of relatively constant recruitment of thalli into the populations. Mortality in lichen populations was divided into deaths which occur after fragmentation of the thallus and deaths which occur after catastrophic environmental events. THe data suggest that the rate of fragmenting death is dependent on the age of the thallus while the rate of catastrophic death is dependent on the number of thalli established in an age class. A comparison of the numbers of thalli in each age class in the eight populations suggested that population density is controlled firstly, by climate and secondly, by variables related to the local rock surface environment. The rate of fragmenting death is related to the diversity of the community and the influence of diversity together with environmental variables in fluctuating or cyclic changes in population number.

Ownership structure and operating performance: evidence from the European maritime industry

In this paper we examine the relation between ownership structure and operating performance for European maritime firms. Using a sample of 266 firm-year observations, during the period 2002–2004, we provide evidence that operating performance is positively related with foreign held shares and investment corporation held shares, indicating better investor protection from managerial opportunism. We also find no relation between operating performance and employee held shares, suggesting no relation between employee commitment and firms’ economic performance. Furthermore, we find no relation between operating performance and government held shares, indicating that government may not adequately protect shareholders’ interests from managerial opportunism. Finally, we do find a positive relation between operating performance and portfolio held shares for code law maritime firms but not for common law maritime firms. Results are robust after adjusting for various firm and country risk characteristics. Overall, our results on the importance of the ownership structure are new to this setting and add to a large body of evidence linking ownership characteristics to corporate performance.

A quantitative study of the pathological changes in the cerebellum in 15 cases of variant Creutzfeldt–Jakob disease (vCJD)

Aims: To determine in the cerebellum in variant Creutzfeldt–Jakob disease (vCJD): (i) whether the pathology affected all laminae; (ii) the spatial topography of the pathology along the folia; (iii) spatial correlations between the pathological changes; and (iv) whether the pathology was similar to that of the common methionine/methionine Type 1 subtype of sporadic CJD. Methods: Sequential cerebellar sections of 15 cases of vCJD were stained with haematoxylin and eosin, or immunolabelled with monoclonal antibody 12F10 against prion protein (PrP) and studied using spatial pattern analysis. Results: Loss of Purkinje cells was evident compared with control cases. Densities of the vacuolation and the protease-resistant form of prion protein (PrPSc) (diffuse and florid plaques) were greater in the granule cell layer (GL) than the molecular layer (ML). In the ML, vacuoles and PrPSc plaques occurred in clusters regularly distributed along the folia with larger clusters of vacuoles and diffuse plaques in the GL. There was a negative spatial correlation between the vacuoles and the surviving Purkinje cells in the ML. There was a positive spatial correlation between the vacuoles and diffuse PrPSc plaques in the ML and GL. Conclusions: (i) all laminae were affected by the pathology, the GL more severely than the ML; (ii) the pathology was topographically distributed along the folia especially in the Purkinje cell layer and ML; (iii) pathological spread may occur in relation to the loop of anatomical connections involving the cerebellum, thalamus, cerebral cortex and pons; and (iv) there were pathological differences compared with methionine/methionine Type 1 sporadic CJD.

Statistical study of talk spurts and silence lengths with relation to various languages and speakers for advanced silence suppression in Next Generation Wireless Communication Devices

As wireless network technologies evolve towards an All-IP framework, Next Generation Wireless Communication Devices demand better use of spectral resources by employing advanced techniques of silence suppression. This paper presents an analysis of VoIP call data and compares the statistical results based on observed patterns of talk spurts and silence lengths to those achieved by a modified on-off voice model for silence suppression in wireless networks. As talk spurts and silence lengths are sensitive to varying word lengths, temporal structure and other prosodic aspects of speech, the impact of the use of various languages, dialects and gender of speakers on these results is also assessed.