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.

Structure-activity relationship of calcitonin gene-related peptide and its receptor

DUE TO COPYRIGHT RESTRICTIONS ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY AND INFORMATION SERVICES WITH PRIOR ARRANGEMENT

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.

Novel peptide antagonists of adrenomedullin and calcitonin gene-related peptide receptors:identification, pharmacological characterization, and interactions with position 74 in receptor activity-modifying protein 1/3

Human adrenomedullin (AM) is a 52-amino acid peptide belonging to the calcitonin peptide family, which also includes calcitonin gene-related peptide (CGRP) and AM2. The two AM receptors, AM(1) and AM(2), are calcitonin receptor-like receptor (CL)/receptor activity-modifying protein (RAMP) (RAMP2 and RAMP3, respectively) heterodimers. CGRP receptors comprise CL/RAMP1. The only human AM receptor antagonist (AM(22-52)) is a truncated form of AM; it has low affinity and is only weakly selective for AM(1) over AM(2) receptors. To develop novel AM receptor antagonists, we explored the importance of different regions of AM in interactions with AM(1), AM(2), and CGRP receptors. AM(22-52) was the framework for generating further AM fragments (AM(26-52) and AM(30-52)), novel AM/alphaCGRP chimeras (C1-C5 and C9), and AM/AM(2) chimeras (C6-C8). cAMP assays were used to screen the antagonists at all receptors to determine their affinity and selectivity. Circular dichroism spectroscopy was used to investigate the secondary structures of AM and its related peptides. The data indicate that the structures of AM, AM2, and alphaCGRP differ from one another. Our chimeric approach enabled the identification of two nonselective high-affinity antagonists of AM(1), AM(2), and CGRP receptors (C2 and C6), one high-affinity antagonist of AM(2) receptors (C7), and a weak antagonist selective for the CGRP receptor (C5). By use of receptor mutagenesis, we also determined that the C-terminal nine amino acids of AM seem to be responsible for its interaction with Glu74 of RAMP3. We provide new information on the structure-activity relationship of AM, alphaCGRP, and AM2 and how AM interacts with CGRP and AM(2) receptors.

Structure-activity relationships for α-calcitonin gene-related peptide

Calcitonin gene-related peptide (CGRP) is a member of the calcitonin (CT) family of peptides. It is a widely distributed neuropeptide implicated in conditions such as neurogenic inflammation. With other members of the CT family, it shares an N-terminal disulphide-bonded ring which is essential for biological activity, an area of potential α-helix, and a C-terminal amide. CGRP binds to the calcitonin receptor-like receptor (CLR) in complex with receptor activity-modifying protein 1 (RAMP1), a member of the family B (or secretin-like) GPCRs. It can also activate other CLR or calcitonin-receptor/RAMP complexes. This 37 amino acid peptide comprises the N-terminal ring that is required for receptor activation (residues 1-7); an α-helix (residues 8-18), a region incorporating a β-bend (residues 19-26) and the C-terminal portion (residues 27-37), that is characterized by bends between residues 28-30 and 33-34. A few residues have been identified that seem to make major contributions to receptor binding and activation, with a larger number contributing either to minor interactions (which collectively may be significant), or to maintaining the conformation of the bound peptide. It is not clear if CGRP follows the pattern of other family B GPCRs in binding largely as an α-helix. Linked Articles This article is part of a themed section on Neuropeptides. To view the other articles in this section visit http://dx.doi.org/10.1111/bph.2013.170.issue-7 © 2012 The Authors. British Journal of Pharmacology © 2012 The British Pharmacological Society.

Structure-activity relationships of the N-terminus of calcitonin gene-related peptide:key roles of alanine-5 and threonine-6 in receptor activation

Background and purpose - The N-terminus of calcitonin gene-related peptide (CGRP) is important for receptor activation, especially the disulphide-bonded ring (residues 1-7). However, the roles of individual amino acids within this region have not been examined and so the molecular determinants of agonism are unknown. This study has examined the role of residues 1, 3-6 and 8-9, excluding Cys-2 and Cys-7. Experimental approach - CGRP derivatives were substituted with either cysteine or alanine; further residues were introduced at position 6. Their affinity was measured by radioligand binding and their efficacy by measuring cAMP production in SK-N-MC cells and ß-arrestin 2 translocation in CHO-K1 cells at the CGRP receptor. Key results - Substitution of Ala-5 by cysteine reduced affinity 270-fold and reduced efficacy for production of cAMP in SK-N-MCs. Potency at ß-arrestin translocation was reduced by 9-fold. Substitution of Thr-6 by cysteine destroyed all measurable efficacy of both cAMP and ß-arrestin responses; substitution with either alanine or serine impaired potency. Substitutions at positions 1, 4, 8 and 9 resulted in approximately 10-fold reductions in potency at both responses. Similar observations were made at a second CGRP-activated receptor, the AMY1(a) receptor. Conclusions and implications - Ala-5 and Thr-6 are key determinants of agonist activity for CGRP. Ala-5 is also very important for receptor binding. Residues outside of the 1-7 ring also contribute to agonist activity.

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.

Structure-activity relations in Cs-doped heteropolyacid catalysts for biodiesel production

A series of insoluble heteropolytungstate (H3PW12O40 HPW) salts, CsxH3−xPW12O40 (x=0.9–3x=0.9–3), were synthesized and characterized using a range of bulk and surface sensitive probes including N2 porosimetry, powder XRD, FTIR, XPS, 31P MAS NMR, and NH3 calorimetry. Materials with Cs content in the range x=2.0–2.7x=2.0–2.7 were composed of dispersed crystallites with surface areas ∼100 m2 g−1 and high Brönsted acid strengths [ΔH0ads(NH3)=−150 kJmol−1], similar to the parent heteropolyacid. The number of accessible surface acid sites probed by α -pinene isomerization correlated well with those determined by NH3 adsorption calorimetry and surface area measurements. CsxH3−xPW12O40 were active toward the esterification of palmitic acid and transesterification of tributyrin, important steps in fatty acid and ester processing for biodiesel synthesis. Optimum performance occurs for Cs loadings of x=2.0–2.3x=2.0–2.3, correlating with the accessible surface acid site density. These catalysts were recoverable with no leaching of soluble HPW.

Synthesis and evaluation of N1-substituted-3-propyl-1,4-benzodiazepine-2-ones as cholecystokinin (CCK2) receptor ligands

A novel synthetic approach towards N1-alkylated 3-propyl-1,4-benzodiazepines was developed in five synthetic steps from 2-amino-4-chlorobenzophenone, in which the N-oxide 4 served as a key intermediate. The structure-activity relationship optimization of this 3-prophyl-1,4-benzodiazepine template was carried out on the N1-position by selective alkylation reactions and resulted in a ligand with an improved affinity on the cholecystokinin (CCK2) receptor. The N-allyl-3-propyl-benzodiazepine 6d displayed an affinity towards the CCK2 (CCK-B) receptor of 170 nM in a radiolabelled receptor-binding assay. The anxiolytic activity of this allyl-3-propyl-1,4-benzodiazepine 6d was subsequently determined in in-vivo psychotropic assays. This novel ligand had ED50 values of 4.7 and 5.2 mg kg-1 in the black and white box test and the x-maze, respectively, and no significant sedation/muscle relaxation was observed.

Synthesis and evaluation of 5-arylated 2(5H)-furanones and 2-arylated pyridazin-3(2H)-ones as anti-cancer agents

Bis-cyclic butenolides, 5-arylated 2(5H)-furanones 6a-c, 7a, b and the 3(2H)-pyridazones 9a-d were prepared by using the aldehyde form of muco halogen acids in electrophilic substitution reactions and in an aldol-like condensation reaction. The cytotoxicity of these simple and bis-cyclic butenolides have been evaluated in tissue culture studies on MAC 13 and MAC 16 murine colon cancer cell lines. The butyl furanone 3 displayed the highest cytotoxicity of 3 μM, as one selected example of a series of dichlorinated pseudoesters. The 5-arylated 2(5H)-furanones 6 and 7 did not show a structure-activity relationship (SAR) depending on the substitution pattern of the aromatic system. An IC50 (concentration inhibiting growth by 50%) was found within a range of 30-50 and 40-50 μM for the MAC 13 and MAC 16 cell lines, respectively. The pyridazine series 9 showed a maximum in-vitro activity for the p-methoxydrivative 9b, having an IC50 of 17 in MAC 13 and 11 μM in MAC 16 cell lines. Selected examples of each series and further novel 2(5H)-furanones such as the hydrazone 5 and the hydantoin 8 have been screened in-vivo in mice and the data are presented. For the pyridazines 9a-d, the in-vitro cytotoxicity correlated with an in-vivo inhibition of tumour growth. The ring expansion of the 5-membered 2(5H)-furanone ring system such as 6a into the 6-membered 3(2H)-pyridazone 9b led to an agent with improved antineoplastic properties. On the resistant MAC 16 cell line the pyridazone 9b displayed 52% tumour inhibition in mice at a dose of 50 mg kg-1 compared with 27% for the 5-FU standard.

Production of membrane proteins in yeast

Background Yeast is an important and versatile organism for studying membrane proteins. It is easy to cultivate and can perform higher eukaryote-like post-translational modifications. S. cerevisiae has a fully-sequenced genome and there are several collections of deletion strains available, whilst P. pastoris can produce very high cell densities (230 g/l). Results We have used both S. cerevisiae and P. pastoris to over-produce the following His6 and His10 carboxyl terminal fused membrane proteins. CD81 – 26 kDa tetraspanin protein (TAPA-1) that may play an important role in the regulation of lymphoma cell growth and may also act as the viral receptor for Hepatitis C-Virus. CD82 – 30 kDa tetraspanin protein that associates with CD4 or CD8 cells and delivers co-stimulatory signals for the TCR/CD3 pathway. MC4R – 37 kDa seven transmembrane G-protein coupled receptor, present on neurons in the hypothalamus region of the brain and predicted to have a role in the feast or fast signalling pathway. Adt2p – 34 kDa six transmembrane protein that catalyses the exchange of ADP and ATP across the yeast mitochondrial inner membrane. Conclusion We show that yeasts are flexible production organisms for a range of different membrane proteins. The yields are such that future structure-activity relationship studies can be initiated via reconstitution, crystallization for X-ray diffraction or NMR experiments.

Synthesis and cytotoxity of oxysterols:studies on the A,B ring polyoxygenated sterols

Oxysterols (OS), the polyoxygenated sterols, represent a class of potent regulatory molecules for important biological actions. Cytotoxicity of OS is one of the most important aspects in studies of OS bioactivities. However, studies, the structure-activity relationship (SAR) study in particular, have been hampered by the limited availability of structurally diverse OS in numbers and amounts. The aim of this project was to develop robust synthetic methods for the preparation of polyhydroxyl sterols, thereof, evaluate their cytotoxicity and establish structure-activity relationship. First, we found hydrophobicity of the side chain is essential for 7-HC's cytotoxicity, and a limited number of hydroxyl groups and a desired configuration on the A, B ring are required for a potent cytotoxicity of an OS, after syntheses and tests of a number of 7-HC's analogues against cancer cell lines. Then polyoxygenation of cholesterol A, B rings was explored. A preparative method for the synthesis of four diastereomerically pure cholest-4-en-3,6-diols was developed. Epoxidation on these cholest-4-en-3,6-diols showed that an allyl group exerts an auxiliary role in producing products with desired configuration in syntheses of the eight diastereomerically pure 45-epoxycholestane-3,6-diols. Reduction of the eight 45-epoxycholestane-3,6-diols produced all eight isomers of the cytotoxic 5α-acholestane 3β,5,6β-triol (CT) for the first time. Epoxide ring opening with protic or Lewis acids on the eight 45-epoxycholestane-3,6-diols are carefully studied. The results demonstrated a combination of an acid and a solvent affected the outcomes of a reaction dramatically. Acyl group participation and migration play an important role with numbers of substrates under certain conditions. All the eight 4,5-trans cholestane- 3,4,5,6-tetrols were synthesised through manipulation of acyl participation. Furthermore these reaction conditions were tested when a number of cholestane-3,4, 5,6,7-pentols and other C3-C7 oxygenated sterols were synthesised for the first time. Introduction of an oxygenated functional group through cholest-2-ene derivatives was studied. The elimination of 3-(4-toluenesulfonate) esters showed the interaction between the existing hydroxyls or acyls with the reaction centre often resulted in different products. The allyl oxidation, epoxidation and Epoxide ring opening reactions are investigated with these cholest-2-enes.