Search PubMed⌕ Search

Biomedical subjects

Darren R Flower

Publications and source records attributed to Darren R Flower.

At least 37 records · Page 2Linked to original sources
PubMed2005

Molecular basis of peptide recognition by the TCR: affinity differences calculated using large scale computing.

Free energy calculations of the wild-type and the variant human T cell lymphotropic virus type 1 Tax peptide presented by the MHC to the TCR have been performed using large scale massively parallel molecular dynamics simulations. The computed free energy difference (-1.86 +/- 0.44 kcal/mol) using alchemical mutation-based thermodynamic integration agrees well with experimental data (-2.9 +/- 0.2 kcal/mol). Our simulations exploit state-of-the-art hardware and codes whose algorithms have been optimized for supercomputing platforms. This enables us to simulate larger, more realistic biological systems for longer durations without the imposition of artificial constraints.

Binding Sites↗
PubMed2005

In silico identification of supertypes for class II MHCs.

The development of epitope-based vaccines, which have wide population coverage, is greatly complicated by MHC polymorphism. The grouping of alleles into supertypes, on the basis of common structural and functional features, addresses this problem directly. In the present study we applied a combined bioinformatics approach, based on analysis of both protein sequence and structure, to identify similarities in the peptide binding sites of 2225 human class II MHC molecules, and thus define supertypes and supertype fingerprints. Two chemometric techniques were used: hierarchical clustering using three-dimensional Comparative Similarity Indices Analysis fields and nonhierarchical k-means clustering using sequence-based z-descriptors. An average consensus of 84% was achieved, i.e., 1872 of 2225 class II molecules were classified in the same supertype by both techniques. Twelve class II supertypes were defined: five DRs, three DQs, and four DPs. The HLA class II supertypes and their fingerprints given in parenthesis are DR1 (Trp(9beta)), DR3 (Glu(9beta), Gln(70beta), and Gln/Arg(74beta)), DR4 (Glu(9beta), Gln/Arg(70beta), and Glu/Ala(74beta)), DR5 (Glu(9beta), Asp(70beta)), and DR9 (Lys/Gln(9beta)); DQ1 (Ala/Gly(86beta)), DQ2 (Glu(86beta), Lys(71beta)), and DQ3 (Glu(86beta), Thr/Asp(71beta)); DPw1 (Asp(84beta) and Lys(69beta)), DPw2 (Gly/Val(84beta) and Glu(69beta)), DPw4 (Gly/Val(84beta) and Lys(69beta)), and DPw6 (Asp(84beta) and Glu(69beta)). Apart from the good agreement between known binding motifs and our classification, several new supertypes, and corresponding thematic binding motifs, were also defined.

Alleles↗
PubMed2005

Towards the chemometric dissection of peptide--HLA-A*0201 binding affinity: comparison of local and global QSAR models.

The affinities of 177 nonameric peptides binding to the HLA-A*0201 molecule were measured using a FACS-based MHC stabilisation assay and analysed using chemometrics. Their structures were described by global and local descriptors, QSAR models were derived by genetic algorithm, stepwise regression and PLS. The global molecular descriptors included molecular connectivity chi indices, kappa shape indices, E-state indices, molecular properties like molecular weight and log P, and three-dimensional descriptors like polarizability, surface area and volume. The local descriptors were of two types. The first used a binary string to indicate the presence of each amino acid type at each position of the peptide. The second was also position-dependent but used five z-scales to describe the main physicochemical properties of the amino acids forming the peptides. The models were developed using a representative training set of 131 peptides and validated using an independent test set of 46 peptides. It was found that the global descriptors could not explain the variance in the training set nor predict the affinities of the test set accurately. Both types of local descriptors gave QSAR models with better explained variance and predictive ability. The results suggest that, in their interactions with the MHC molecule, the peptide acts as a complicated ensemble of multiple amino acids mutually potentiating each other.

Algorithms↗
PubMed2005

Evolutionary dynamics of hepatitis C virus envelope genes during chronic infection.

Hepatitis C virus (HCV) envelope glycoproteins E1 and E2 are important targets for the host immune response. The genes encoding these proteins exhibit a high degree of variability that gives rise to differing phenotypic traits, including alterations in receptor-binding affinity and immune recognition and escape. In order to elucidate patterns of adaptive evolution during chronic infection, a panel of full-length E1E2 clones was generated from sequential serum samples obtained from four chronically infected individuals. By using likelihood-based methods for phylogenetic inference, the evolutionary dynamics of circulating HCV quasispecies populations were assessed and a site-by-site analysis of the d(N)/d(S) ratio was performed, to identify specific codons undergoing diversifying positive selection. HCV phylogenies, coupled with the number and distribution of selected sites, differed markedly between patients, highlighting that HCV evolution during chronic infection is a patient-specific phenomenon. This analysis shows that purifying selection is the major force acting on HCV populations in chronic infection. Whilst no significant evidence for positive selection was observed in E1, a number of sites under positive selection were identified within the ectodomain of the E2 protein. All of these sites were located in regions hypothesized to be exposed to the selective environment of the host, including a number of functionally defined domains that have been reported to be involved in immune evasion and receptor binding. Dated-tip methods for estimation of underlying HCV mutation rates were also applied to the data, enabling prediction of the most recent common ancestor for each patient's quasispecies.

Amino Acid Sequence↗
PubMed2005

An update on the functional molecular immunology (FIMM) database.

Data on the major histocompatibility complex, T-cell epitopes, B-cell epitopes, antigens and diseases are heterogeneous and scattered among different databases and the literature. Since it has become increasingly difficult to obtain an integrated view of functional immune response components, we have developed and updated over several years the Functional molecular IMMunology (FIMM) database (http:// research.i2r.a-star.edu.sg/fimm/). FIMM contains integrated expert-curated data on protein antigens, and on human immunological receptors that recognise and bind them in healthy or disease states. Interfaces with multiple, intuitive query options and query reports provide immunologists with prioritised information that aids data interpretation, vaccine target discovery and immune disease research.

Allergy and Immunology↗
PubMed2004

Benchmarking B cell epitope prediction: underperformance of existing methods.

Sequence profiling is used routinely to predict the location of B-cell epitopes. In the postgenomic era, the need for reliable epitope prediction is clear. We assessed 484 amino acid propensity scales in combination with ranges of plotting parameters to examine exhaustively the correlation of peaks and epitope location within 50 proteins mapped for polyclonal responses. After examining more than 10(6) combinations, we found that even the best set of scales and parameters performed only marginally better than random. Our results confirm the null hypothesis: Single-scale amino acid propensity profiles cannot be used to predict epitope location reliably. The implication for studies using such methods is obvious.

Amino Acids↗
PubMed2004

Transporter associated with antigen processing preselection of peptides binding to the MHC: a bioinformatic evaluation.

TAP is responsible for the transit of peptides from the cytosol to the lumen of the endoplasmic reticulum. In an immunological context, this event is followed by the binding of peptides to MHC molecules before export to the cell surface and recognition by T cells. Because TAP transport precedes MHC binding, TAP preferences may make a significant contribution to epitope selection. To assess the impact of this preselection, we have developed a scoring function for TAP affinity prediction using the additive method, have used it to analyze and extend the TAP binding motif, and have evaluated how well this model acts as a preselection step in predicting MHC binding peptides. To distinguish between MHC alleles that are exclusively dependent on TAP and those exhibiting only a partial dependence on TAP, two sets of MHC binding peptides were examined: HLA-A*0201 was selected as a representative of partially TAP-dependent HLA alleles, and HLA-A*0301 represented fully TAP-dependent HLA alleles. TAP preselection has a greater impact on TAP-dependent alleles than on TAP-independent alleles. The reduction in the number of nonbinders varied from 10% (TAP-independent) to 33% (TAP-dependent), suggesting that TAP preselection is an important component in the successful in silico prediction of T cell epitopes.

ATP Binding Cassette Transporter, Subfamily B, Mem↗
PubMed2004

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide: can the alpha3 and beta2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide-MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the alpha1-alpha2 domain, the alpha3 and beta2m domains cannot be neglected.

Algorithms↗
PubMed2004

Determinants of human immunodeficiency virus type 1 escape from the primary CD8+ cytotoxic T lymphocyte response.

CD8+ cytotoxic T lymphocytes (CTLs) play an important role in containment of virus replication in primary human immunodeficiency virus (HIV) infection. HIV's ability to mutate to escape from CTL pressure is increasingly recognized; but comprehensive studies of escape from the CD8 T cell response in primary HIV infection are currently lacking. Here, we have fully characterized the primary CTL response to autologous virus Env, Gag, and Tat proteins in three patients, and investigated the extent, kinetics, and mechanisms of viral escape from epitope-specific components of the response. In all three individuals, we observed variation beginning within weeks of infection at epitope-containing sites in the viral quasispecies, which conferred escape by mechanisms including altered peptide presentation/recognition and altered antigen processing. The number of epitope-containing regions exhibiting evidence of early CTL escape ranged from 1 out of 21 in a subject who controlled viral replication effectively to 5 out of 7 in a subject who did not. Evaluation of the extent and kinetics of HIV-1 escape from >40 different epitope-specific CD8 T cell responses enabled analysis of factors determining escape and suggested that escape is restricted by costs to intrinsic viral fitness and by broad, codominant distribution of CTL-mediated pressure on viral replication.

Amino Acid Sequence↗
PubMed2004

New horizons in mouse immunoinformatics: reliable in silico prediction of mouse class I histocompatibility major complex peptide binding affinity.

Quantitative structure-activity relationship (QSAR) analysis is a main cornerstone of modern informatic disciplines. Predictive computational models, based on QSAR technology, of peptide-major histocompatibility complex (MHC) binding affinity have now become a vital component of modern day computational immunovaccinology. Historically, such approaches have been built around semi-qualitative, classification methods, but these are now giving way to quantitative regression methods. The additive method, an established immunoinformatics technique for the quantitative prediction of peptide-protein affinity, was used here to identify the sequence dependence of peptide binding specificity for three mouse class I MHC alleles: H2-D(b), H2-K(b) and H2-K(k). As we show, in terms of reliability the resulting models represent a significant advance on existing methods. They can be used for the accurate prediction of T-cell epitopes and are freely available online ( http://www.jenner.ac.uk/MHCPred).

Algorithms↗
PubMed2004

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.

Amino Acid Sequence↗
PubMed2004

Identifiying human MHC supertypes using bioinformatic methods.

Classification of MHC molecules into supertypes in terms of peptide-binding specificities is an important issue, with direct implications for the development of epitope-based vaccines with wide population coverage. In view of extremely high MHC polymorphism (948 class I and 633 class II HLA alleles) the experimental solution of this task is presently impossible. In this study, we describe a bioinformatics strategy for classifying MHC molecules into supertypes using information drawn solely from three-dimensional protein structure. Two chemometric techniques-hierarchical clustering and principal component analysis-were used independently on a set of 783 HLA class I molecules to identify supertypes based on structural similarities and molecular interaction fields calculated for the peptide binding site. Eight supertypes were defined: A2, A3, A24, B7, B27, B44, C1, and C4. The two techniques gave 77% consensus, i.e., 605 HLA class I alleles were classified in the same supertype by both methods. The proposed strategy allowed "supertype fingerprints" to be identified. Thus, the A2 supertype fingerprint is Tyr(9)/Phe(9), Arg(97), and His(114) or Tyr(116); the A3-Tyr(9)/Phe(9)/Ser(9), Ile(97)/Met(97) and Glu(114) or Asp(116); the A24-Ser(9) and Met(97); the B7-Asn(63) and Leu(81); the B27-Glu(63) and Leu(81); for B44-Ala(81); the C1-Ser(77); and the C4-Asn(77).

Alleles↗
PubMed2004

Integrative bioinformatics for functional genome annotation: trawling for G protein-coupled receptors.

G protein-coupled receptors (GPCR) are amongst the best studied and most functionally diverse types of cell-surface protein. The importance of GPCRs as mediates or cell function and organismal developmental underlies their involvement in key physiological roles and their prominence as targets for pharmacological therapeutics. In this review, we highlight the requirement for integrated protocols which underline the different perspectives offered by different sequence analysis methods. BLAST and FastA offer broad brush strokes. Motif-based search methods add the fine detail. Structural modelling offers another perspective which allows us to elucidate the physicochemical properties that underlie ligand binding. Together, these different views provide a more informative and a more detailed picture of GPCR structure and function. Many GPCRs remain orphan receptors with no identified ligand, yet as computer-driven functional genomics starts to elaborate their functions, a new understanding of their roles in cell and developmental biology will follow.

Algorithms↗
PubMed2004

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.

Amino Acid Motifs↗
PubMed2004

Quantitative online prediction of peptide binding to the major histocompatibility complex.

With its implications for vaccine discovery, the accurate prediction of T cell epitopes is one of the key aspirations of computational vaccinology. We have developed a robust multivariate statistical method, based on partial least squares, for the quantitative prediction of peptide binding to major histocompatibility complexes (MHC), the principal checkpoint on the antigen presentation pathway. As a service to the immunobiology community, we have made a Perl implementation of the method available via a World Wide Web server. We call this server MHCPred. Access to the server is freely available from the URL: http://www.jenner.ac.uk/MHCPred. We have exemplified our method with a model for peptides binding to the common human MHC molecule HLA-B*3501.

Antigen Presentation↗
PubMed2003

MHCPred: A server for quantitative prediction of peptide-MHC binding.

Accurate T-cell epitope prediction is a principal objective of computational vaccinology. As a service to the immunology and vaccinology communities at large, we have implemented, as a server on the World Wide Web, a partial least squares-based multivariate statistical approach to the quantitative prediction of peptide binding to major histocom- patibility complexes (MHC), the key checkpoint on the antigen presentation pathway within adaptive cellular immunity. MHCPred implements robust statistical models for both Class I alleles (HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3301, HLA-A*6801, HLA-A*6802 and HLA-B*3501) and Class II alleles (HLA-DRB*0401, HLA-DRB*0401 and HLA-DRB*0701). MHCPred is available from the URL: http://www.jenner.ac.uk/MHCPred.

Binding Sites↗