PubMed Health⌕ Search

Biomedical subjects

J M Thornton

Publications and source records attributed to J M Thornton.

At least 19 recordsLinked to original sources
PubMed2006

SPINE bioinformatics and data-management aspects of high-throughput structural biology.

SPINE (Structural Proteomics In Europe) was established in 2002 as an integrated research project to develop new methods and technologies for high-throughput structural biology. Development areas were broken down into workpackages and this article gives an overview of ongoing activity in the bioinformatics workpackage. Developments cover target selection, target registration, wet and dry laboratory data management and structure annotation as they pertain to high-throughput studies. Some individual projects and developments are discussed in detail, while those that are covered elsewhere in this issue are treated more briefly. In particular, this overview focuses on the infrastructure of the software that allows the experimentalist to move projects through different areas that are crucial to high-throughput studies, leading to the collation of large data sets which are managed and eventually archived and/or deposited.

Computational Biology↗
PubMed2006

Haplotypes in the CTLA4 region are associated with coeliac disease in the Irish population.

Chromosomal region 2q33 encodes the immune regulatory genes, CTLA4, ICOS and CD28, which are involved in regulation of T-cell activity and has been studied as a candidate gene locus in autoimmune diseases, including coeliac disease (CD). We have investigated whether an association exists between this region and CD in the Irish population using a comprehensive analysis for genetic variation. Using a haplotype-tagging approach, this gene cluster was investigated for disease association in a case-control study comprising 394 CD patients and 421 ethnically matched healthy controls. Several SNPs, including CTLA4_CT60, showed association with disease; however, after correction for multiple-testing, CTLA4-658C/T was the only polymorphism found to show significant association with disease when allele, genotype, or carrier status frequency were analysed (carrier status (Allele C), P = 0.0016). Haplotype analysis revealed a haplotype incorporating the CD28/CTLA4 and two 5' ICOS polymorphisms to be significantly associated with disease (patients 24.1%; controls 31.5%; P = 0.035), as was a shorter haplotype composed of the CTLA4 markers only (30.9 vs 34.9%; P = 0.042). The extended haplotype incorporating CD28/CTLA4 and 5' ICOS is more strongly associated with disease than haplotypes of individual genes. This suggests a causal variant associated with this haplotype may be associated with disease in this population.

Antigens, CD↗
PubMed2005

Amino acid architecture and the distribution of polar atoms on the surfaces of proteins.

We propose that a necessary condition for a protein to be soluble is the absence of large hydrophobic patches on its solvent-accessible surface, which can cause aggregation to occur. We note that the polar nature of the backbone of all amino acids guarantees a minimum polar content and hence can interrupt such patches. As a result, a carefully conserved detailed atomic placement of residues on the protein surface is not necessary for solubility. In order to demonstrate this, we construct a measure based on the average hydrophobicity of a simply defined patch. We use this measurement to compare surfaces that exhibit a clear difference in their solubility properties, namely, a) the solvent accessible surfaces for a set of homo-dimers and the surfaces buried in their interfaces and b) for a set of monomers the surfaces of fragments of secondary structure which are solvent accessible/inaccessible. Having demonstrated a difference in the first set of distributions, we characterize the solvent accessible surfaces of monomeric proteins. To test if cooperative behavior occurs between the atoms for these surfaces, we construct a set of randomized surfaces, which obey a very simple stereochemical constraint. We find that the observed and randomized distributions are much more similar than the previous sets we examined. This implies that while surfaces of soluble proteins must have sufficient polar content, the relative placement of atoms of one amino acid with respect to the atoms of neighboring amino acid need not be finely tuned, which provides an innate robustness for protein design and folding.

Amino Acids↗
PubMed2005

HTHquery: a method for detecting DNA-binding proteins with a helix-turn-helix structural motif.

SUMMARY: HTHquery is a web-based service to determine if a protein structure has a helix-turn-helix structural motif which could bind to DNA. It is based on a similarity with a set of structural templates, the accessibility of a putative structural motif and a positive electrostatic potential in the neighbourhood of the putative motif. A set of scores are computed, based on each template, using a linear predictor. From the training set used, the predictor has a true positive rate of 83.5% and a false positive rate of 0.8%. The emphasis for the website is on providing a straightforward interface which can be easily used by a bench-based scientist. AVAILABILITY: HTHquery is implemented using a set of Perl scripts and C program and can be accessed freely on the website http://www.ebi.ac.uk/thornton-srv/databases/HTHquery.

Amino Acid Motifs↗
PubMed2005

Chromosome 5q candidate genes in coeliac disease: genetic variation at IL4, IL5, IL9, IL13, IL17B and NR3C1.

Genetic predisposition to coeliac disease (CD) is determined primarily by alleles at the HLA-DQB locus, and evidence exists implicating other major histocompatibility complex-linked genes (6p21) and the CTLA4 locus on chromosome 2q33. In addition, extensive family studies have provided strong, reproducible evidence for a susceptibility locus on chromosome 5q (CELIAC2). However, the gene responsible has not been identified. We have assayed genetic variation at the IL4, IL5, IL9, IL13, IL17B and NR3C1 (GR) loci, all of which are present on chromosome 5q and have potential or demonstrated involvement in autoimmune and/or inflammatory disease, in a sample of 409 CD cases and 355 controls. Thirteen single nucleotide polymorphisms were chosen on the basis of functional relevance, prior disease association and, where possible, prior knowledge of the haplotype variation present in European populations. There were no statistically significant allele or haplotype frequency differences between cases and controls. Therefore, these results provide no evidence that these loci are associated with CD in this sample population.

Case-Control Studies↗
PubMed2004

An examination of the conservation of surface patch polarity for proteins.

MOTIVATION: The solubility of a protein is crucial for its function and is therefore an evolutionary constraint. As the solubility of a protein is related to the distribution of polar and hydrophobic residues on its solvent accessible surface, such a constraint should provide a valuable insight into the evolution of protein surfaces. We examine how the surfaces of proteins have evolved by considering how the average hydrophobicities of patches of surface residues vary across homologous proteins. We derive distributions for the average hydrophobicity/philicity of surface patches at a residue-based level-which we refer to as the residue hydrophobic density. This is computed for a set of 28 monomeric proteins and their homologues. The resulting distributions are compared with a set of randomized sequences, with the same residue content. RESULTS: We find that the patches, involving typically more than 10 residues, maintain a more hydrophilic surface than one would expect from a random substitution model, indicating a cooperative behaviour for these surfaces residues in terms of this single variable. SUPPLEMENTARY INFORMATION: Additional plots for all of the proteins examined in this paper can be found at: http://www.ebi.ac.uk/~shanahan/PCon/index.html

Algorithms↗
PubMed2004

Haplotype variation at the IBD5/SLC22A4 locus (5q31) in coeliac disease in the Irish population.

In addition to the well-established association of coeliac disease (CD) with HLA-DQ (6p21) and possibly CTLA4 (2q33), there is considerable evidence for a susceptibility locus on chromosome 5q, which contains many potential candidates for inflammatory disease, including a cluster of cytokine genes in 5q31. CD cases and controls were genotyped for four single-nucleotide polymorphism (SNP) markers that together characterize >90% of the haplotype variation at the IBD5 locus encoding, among others, the SLC22A4 gene. IBD5 and SLC22A4 map to 5q31 and have recently been associated with Crohn's disease and rheumatoid arthritis. Haplotype frequencies do not differ significantly between CD cases and controls in the Irish population, and therefore the chromosome 5 CD susceptibility locus most likely lies elsewhere on 5q.

Celiac Disease↗
PubMed2001

On the molecular discrimination between adenine and guanine by proteins.

The molecular recognition and discrimination of adenine and guanine ligand moieties in complexes with proteins have been studied using empirical observations on carefully selected crystal structures. The distribution of protein folds that bind these purines has been found to differ significantly from that across the whole PDB, but the most populated architectures and folds are also the most common in three genomes from the three different domains of life. The protein environments around the two nucleic acid bases were significantly different, in terms of the propensities of amino acid residues to be in the binding site, as well as their propensities to form hydrogen bonds to the bases. Plots of the distribution of protein atoms around the two purines clearly show different clustering of hydrogen bond donors and acceptors opposite complimentary acceptors and donors in the rings, with hydrophobic areas below and above the rings. However, the clustering pattern is fuzzy, reflecting the variety of ways that proteins have evolved to recognise the same molecular moiety. Furthermore, an analysis of the conservation of residues in the protein chains binding guanine shows that residues in contact with the base are in general better conserved than the rest of the chain.

Adenine↗
PubMed2001

Conservation helps to identify biologically relevant crystal contacts.

Some crystal contacts are biologically relevant, most are not. We assess the utility of combining measures of size and conservation to discriminate between biological and non-biological contacts. Conservation and size information is calculated for crystal contacts in 53 families of homodimers and 65 families of monomers. Biological contacts are shown to be usually conserved and typically the largest contact in the crystal. A range of neural networks accepting different combinations and encodings of this information is used to answer the following questions: (1) is a given crystal contact biological, and (2) given all crystal contacts in a homodimer, which is the biological one? Predictions for (1) are performed on both homodimer and monomer datasets. The best performing neural network combined size and conservation inputs. For the homodimers, it correctly classified 48 out of 53 biological contacts and 364 out of 366 non-biological contacts, giving a combined accuracy of 98.3 %. A more robust performance statistic, the phi-coefficient, which accounts for imbalances in the dataset, gave a value of 0.92. Taking all 535 non-biological contacts from the 65 monomers, this predictor made erroneous classifications only 4.3 % of the time. Predictions for (2) were performed on homodimers only. The best performing network achieved a prediction accuracy of 98.1 % using size information alone. We conclude that in answering question (1) size and conservation combined discriminate biological from non-biological contacts better than either measure alone. For answering question (2), we conclude that in our dataset size is so powerful a discriminant that conservation adds little predictive benefit.

Animals↗
PubMed2001

The evolution and structural anatomy of the small molecule metabolic pathways in Escherichia coli.

The 106 small molecule metabolic (SMM) pathways in Escherichia coli are formed by the protein products of 581 genes. We can define 722 domains, nearly all of which are homologous to proteins of known structure, that form all or part of 510 of these proteins. This information allows us to answer general questions on the structural anatomy of the SMM pathway proteins and to trace family relationships and recruitment events within and across pathways. Half the gene products contain a single domain and half are formed by combinations of between two and six domains. The 722 domains belong to one of 213 families that have between one and 51 members. Family members usually conserve their catalytic or cofactor binding properties; substrate recognition is rarely conserved. Of the 213 families, members of only a quarter occur in isolation, i.e. they form single-domain proteins. Most members of the other families combine with domains from just one or two other families and a few more versatile families can combine with several different partners. Excluding isoenzymes, more than twice as many homologues are distributed across pathways as within pathways. However, serial recruitment, with two consecutive enzymes both being recruited to another pathway, is rare and recruitment of three consecutive enzymes is not observed. Only eight of the 106 pathways have a high number of homologues. Homology between consecutive pairs of enzymes with conservation of the main substrate-binding site but change in catalytic mechanism (which would support a simple model of retrograde pathway evolution) occurs only six times in the whole set of enzymes. Most of the domains that form SMM pathways have homologues in non-SMM pathways. Taken together, these results imply a pervasive "mosaic" model for the formation of protein repertoires and pathways.

Bacterial Proteins↗
PubMed2001

Amino acid-base interactions: a three-dimensional analysis of protein-DNA interactions at an atomic level.

To assess whether there are universal rules that govern amino acid-base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein-DNA complex structures. DNA-backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base-amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein-DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein-DNA interface. In conclusion, the majority of amino acid-base interactions observed follow general principles that apply across all protein-DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are 'universal' and 'context-dependent'. An interactive Web-based atlas of side chain-base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.

Base Pairing↗
PubMed2001

Identification of higher brain centres that may encode the cardiorespiratory response to exercise in humans.

1. Positron emission tomography (PET) was used to identify the neuroanatomical correlates underlying 'central command' during imagination of exercise under hypnosis, in order to uncouple central command from peripheral feedback. 2. Three cognitive conditions were used: condition I, imagination of freewheeling downhill on a bicycle (no change in heart rate, HR, or ventilation, V(I)): condition II, imagination of exercise, cycling uphill (increased HR by 12 % and V(I) by 30 % of the actual exercise response): condition III, volitionally driven hyperventilation to match that achieved in condition II (no change in HR). 3. Subtraction methodology created contrast A (II minus I) highlighting cerebral areas involved in the imagination of exercise and contrast B (III minus I) highlighting areas activated in the direct volitional control of breathing (n = 4 for both; 8 scans per subject). End-tidal P(CO(2)) (P(ET,CO(2))) was held constant throughout PET scanning. 4. In contrast A, significant activations were seen in the right dorso-lateral prefrontal cortex, supplementary motor areas (SMA), the right premotor area (PMA), superolateral sensorimotor areas, thalamus, and bilaterally in the cerebellum. In contrast B, significant activations were present in the SMA and in lateral sensorimotor cortical areas. The SMA/PMA, dorso-lateral prefrontal cortex and the cerebellum are concerned with volitional/motor control, including that of the respiratory muscles. 5. The neuroanatomical areas activated suggest that a significant component of the respiratory response to 'exercise', in the absence of both movement feedback and an increase in CO(2) production, can be generated by what appears to be a behavioural response.

Adult↗
PubMed2001

Evolution of function in protein superfamilies, from a structural perspective.

The recent growth in protein databases has revealed the functional diversity of many protein superfamilies. We have assessed the functional variation of homologous enzyme superfamilies containing two or more enzymes, as defined by the CATH protein structure classification, by way of the Enzyme Commission (EC) scheme. Combining sequence and structure information to identify relatives, the majority of superfamilies display variation in enzyme function, with 25 % of superfamilies in the PDB having members of different enzyme types. We determined the extent of functional similarity at different levels of sequence identity for 486,000 homologous pairs (enzyme/enzyme and enzyme/non-enzyme), with structural and sequence relatives included. For single and multi-domain proteins, variation in EC number is rare above 40 % sequence identity, and above 30 %, the first three digits may be predicted with an accuracy of at least 90 %. For more distantly related proteins sharing less than 30 % sequence identity, functional variation is significant, and below this threshold, structural data are essential for understanding the molecular basis of observed functional differences. To explore the mechanisms for generating functional diversity during evolution, we have studied in detail 31 diverse structural enzyme superfamilies for which structural data are available. A large number of variations and peculiarities are observed, at the atomic level through to gross structural rearrangements. Almost all superfamilies exhibit functional diversity generated by local sequence variation and domain shuffling. Commonly, substrate specificity is diverse across a superfamily, whilst the reaction chemistry is maintained. In many superfamilies, the position of catalytic residues may vary despite playing equivalent functional roles in related proteins. The implications of functional diversity within supefamilies for the structural genomics projects are discussed. More detailed information on these superfamilies is available at http://www.biochem.ucl.ac.uk/bsm/FAM-EC/.

Binding Sites↗
PubMed2001

Protein-RNA interactions: a structural analysis.

A detailed computational analysis of 32 protein-RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein-double-stranded DNA and protein-single-stranded DNA complexes. The interface properties of the protein-RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein-RNA and protein-DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein-RNA complexes, backbone contacts were more dominant in the protein-DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.

Base Pairing↗
PubMed2001

A rapid classification protocol for the CATH Domain Database to support structural genomics.

In order to support the structural genomic initiatives, both by rapidly classifying newly determined structures and by suggesting suitable targets for structure determination, we have recently developed several new protocols for classifying structures in the CATH domain database (http://www.biochem.ucl.ac.uk/bsm/cath). These aim to increase the speed of classification of new structures using fast algorithms for structure comparison (GRATH) and to improve the sensitivity in recognising distant structural relatives by incorporating sequence information from relatives in the genomes (DomainFinder). In order to ensure the integrity of the database given the expected increase in data, the CATH Protein Family Database (CATH-PFDB), which currently includes 25,320 structural domains and a further 160,000 sequence relatives has now been installed in a relational ORACLE database. This was essential for developing more rigorous validation procedures and for allowing efficient querying of the database, particularly for genome analysis. The associated Dictionary of Homologous Superfamilies [Bray,J.E., Todd,A.E., Pearl,F.M.G., Thornton,J.M. and Orengo,C.A. (2000) Protein Eng., 13, 153-165], which provides multiple structural alignments and functional information to assist in assigning new relatives, has also been expanded recently and now includes information for 903 homologous superfamilies. In order to improve coverage of known structures, preliminary classification levels are now provided for new structures at interim stages in the classification protocol. Since a large proportion of new structures can be rapidly classified using profile-based sequence analysis [e.g. PSI-BLAST: Altschul,S.F., Madden,T.L., Schaffer,A.A., Zhang,J., Zhang,Z., Miller,W. and Lipman,D.J. (1997) Nucleic Acids Res., 25, 3389-3402], this provides preliminary classification for easily recognisable homologues, which in the latest release of CATH (version 1.7) represented nearly three-quarters of the non-identical structures.

Computational Biology↗
PubMed2001

Protein-protein interfaces: analysis of amino acid conservation in homodimers.

Evolutionary information derived from the large number of available protein sequences and structures could powerfully guide both analysis and prediction of protein-protein interfaces. To test the relevance of this information, we assess the conservation of residues at protein-protein interfaces compared with other residues on the protein surface. Six homodimer families are analyzed: alkaline phosphatase, enolase, glutathione S-transferase, copper-zinc superoxide dismutase, Streptomyces subtilisin inhibitor, and triose phosphate isomerase. For each family, random simulation is used to calculate the probability (P value) that the level of conservation observed at the interface occurred by chance. The results show that interface conservation is higher than expected by chance and usually statistically significant at the 5% level or better. The effect on the P values of using different definitions of the interface and of excluding active site residues is discussed.

Amino Acid Sequence↗