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.

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.

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.

PDBsum: summaries and analyses of PDB structures.

PDBsum is a web-based database providing a largely pictorial summary of the key information on each macromolecular structure deposited at the Protein Data Bank (PDB). It includes images of the structure, annotated plots of each protein chain's secondary structure, detailed structural analyses generated by the PROMOTIF program, summary PROCHECK results and schematic diagrams of protein-ligand and protein-DNA interactions. RasMol scripts highlight key aspects of the structure, such as the protein's domains, PROSITE patterns and protein-ligand interactions, for interactive viewing in 3D. Numerous links take the user to related sites. PDBsum is updated whenever any new structures are released by the PDB and is freely accessible via http://www.biochem.ucl.ac.uk/bsm/pdbsum.

Rfree and the rfree ratio. II. Calculation Of the expected values and variances of cross-validation statistics in macromolecular least-squares refinement.

The free R factor is used routinely as a cross-validation tool in macromolecular crystallography. However, without any means of deriving quantitative estimates of its expected value and variance, its application has been rather subjective and its usefulness therefore somewhat limited. In the first part of this series, estimates of the expected value of the ratio of the free R factor to the standard R factor at the convergence of the structure refinement were given. Here, estimates of the variance of this ratio are given and are compared with the observed deviations from the expected values for a selection of refined structures. It is discussed how errors in the functional form of the structure-factor model as well as other types of errors might influence this ratio.

Sequences annotated by structure: a tool to facilitate the use of structural information in sequence analysis.

With the aim of bridging the gap between protein sequence and structural analyses, we have developed a tool to aid the identification of new protein sequences by recognizing distant homologues using structural information. The tool generates sequence annotated by structure (SAS) files, applying structural information derived from structural analyses to a given protein sequence. A World Wide Web interface allows a given sequence to be submitted either for structural annotation or, where its structure is unknown, for search and alignment against sequences of known structure. In both cases, SAS will colour residues in the sequence of known structure according to a selection of properties, including secondary structure, interatomic contacts and active site information. SAS can also be used to inspect properties of a single structure.

Validation of protein models derived from experiment.

The growing number of protein structures solved at atomic resolution holds the promise of further improvements in geometry-based validation parameters. Additionally, the estimated standard uncertainties of the atomic coordinates have been computed for a number of X-ray structures, providing a measure of the coordinate precision. In NMR spectroscopy, a measure analogous to the crystallographic R-factor has been developed.

Rfree and the rfree ratio. I. Derivation of expected values of cross-validation residuals used in macromolecular least-squares refinement.

The last five years have seen a large increase in the use of cross validation in the refinement of macromolecular structures using X-ray data. In this technique a test set of reflections is set aside from the working set and the progress of the refinement is monitored by the calculation of a free R factor which is based only on the excluded reflections. This paper gives estimates for the ratio of the free R factor to the R factor calculated from the working set for both unrestrained and restrained refinement. It is assumed that both the X-ray and restraint observations have been weighted correctly and that there is no correlation of errors between the test and working sets. It is also shown that the least-squares weights that minimize the variances of the refined parameters, also approximately minimize the free R factor. The estimated free R-factor ratios are compared with those reported for structures in the Protein Data Bank.

Error estimates of protein structure coordinates and deviations from standard geometry by full-matrix refinement of gammaB- and betaB2-crystallin.

Faster workstations with larger memories are making error estimation from full-matrix least-squares refinement a more practicable technique in protein crystallography. Using minimum variance weighting, estimated standard deviations of atomic positions have been calculated for two eye lens proteins from the inverse of a least-squares normal matrix which was full with respect to the coordinate parameters. gammaB-crystallin, refined at 1.49 A yielded average errors in atomic positions which ranged from 0.05 A for main-chain atoms to 0.27 A for unrestrained water molecules. The second structure used in this work was that of betaB2-crystallin refined at 2.1 A resolution where the corresponding average errors were 0.08 and 0.35 A, respectively. The relative errors in atomic positions are dependent on the number and kinds of restraints used in the refinements. It is also shown that minimum variance weighting leads to mean-square deviations from target geometry in the refined structures which are smaller than the variances used in the distance weighting.

Dihydrofolate reductase: a potential drug target in trypanosomes and leishmania.

Dihydrofolate reductase has successfully been used as a drug target in the area of anti-cancer, anti-bacterial and anti-malarial chemotherapy. Little has been done to evaluate it as a drug target for treatment of the trypanosomiases and leishmaniasis. A crystal structure of Leishmania major dihydrofolate reductase has been published. In this paper, we describe the modelling of Trypanosoma cruzi and Trypanosoma brucei dihydrofolate reductases based on this crystal structure. These structures and models have been used in the comparison of protozoan, bacterial and human enzymes in order to highlight the different features that can be used in the design of selective anti-protozoan agents. Comparison has been made between residues present in the active site, the accessibility of these residues, charge distribution in the active site, and the shape and size of the active sites. Whilst there is a high degree of similarity between protozoan, human and bacterial dihydrofolate reductase active sites, there are differences that provide potential for selective drug design. In particular, we have identified a set of residues which may be important for selective drug design and identified a larger binding pocket in the protozoan than the human and bacterial enzymes.

Protein folds and functions.

BACKGROUND: The recent rapid increase in the number of available three-dimensional protein structures has further highlighted the necessity to understand the relationship between biological function and structure. Using structural classification schemes such as SCOP, CATH and DALI, it is now possible to explore global relationships between protein fold and function, something which was previously impractical. RESULTS: Using a relational database of CATH data we have generated fold distributions for arbitrary selections of proteins automatically. These distributions have been examined in the light of protein function and bound ligand. Different enzyme classes are not clearly reflected in distributions of protein class and architecture, whereas the type of bound ligand has a much more dramatic effect. CONCLUSIONS: The availability of structural classification data has enabled this novel overview analysis. We conclude that function at the top level of the EC number enzyme classification is not related to fold, as only a very few specific residues are actually responsible for enzyme activity. Conversely, the fold is much more closely related to ligand type.

New tools and resources for analysing protein structures and their interactions.

The determination of protein structures has furthered our understanding of how various proteins perform their functions. With the large number of structures currently available in the PDB, it is necessary to be able to easily study these proteins in detail. Here new software tools are presented which aim to facilitate this analysis; these include the PDBsum WWW site which provides a summary description of all PDB entries, the programs TOPS and NUCPLOT to plot schematic diagrams representing protein topology and DNA-binding interactions, SAS a WWW-based sequence-analysis tool incorporating structural data, and WWW servers for the analysis of protein-protein interfaces and analyses of over 300 haem-binding proteins.

Non-randomness in side-chain packing: the distribution of interplanar angles.

We analyze the distributions of interplanar angles between interacting side chains with well-defined planar regions, to see whether these distributions correspond to random packing or alternatively show orientational preferences. We use a non-homologous set of 79 high-resolution protein chain structures to show that the observed distributions are significantly different from the sinusoidal one expected for random packing. Overall, we see a relative excess of small angles and a paucity of large interplanar angles; the difference between the expected and observed distributions can be described as a shift of 5% of the interplanar angles from large (> or = 60 degrees) to small (< 30 degrees) values. By grouping the residue pairs into categories based on chemical similarity, we find that some categories have very non-sinusoidal interplanar angle distributions, whereas other categories have distributions that are close to sinusoidal. For a few categories, observed deviations from a sinusoidal distribution can be explained by the electrostatic anisotropy of the isolated pair potential energy. In other cases, the observed distributions reflect the longer range effects of different possible interaction geometries. In particular, geometries that disrupt external hydrogen bonding are disfavored.

NUCPLOT: a program to generate schematic diagrams of protein-nucleic acid interactions.

Proteins that bind to DNA are found in all areas of genetic activity within the cell. To help understand how these proteins perform their various functions, it is useful to analyse which residues are involved in binding to the DNA and how they interact with the bases and sugar-phosphate backbone of nucleic acids. Here we describe a program called NUCPLOT which can automatically identify these interactions from the 3D atomic coordinates of the complex from a PDB file and generate a plot that shows all the interactions in a schematic manner. The program produces a PostScript output file representing hydrogen, van der Waals and covalent bonds between the protein and the DNA. The resulting diagram is both clear and simple and allows immediate identification of important interactions within the structure. It also facilitates comparison of binding found in different structures. NUCPLOT is a completely automatic program, which can be used for any protein-DNA complex and will also work for certain protein-RNA structures.

Protein clefts in molecular recognition and function.

One of the primary factors determining how proteins interact with other molecules is the size of clefts in the protein's surface. In enzymes, for example, the active site is often characterized by a particularly large and deep cleft, while interactions between the molecules of a protein dimer tend to involve approximately planar surfaces. Here we present an analysis of how cleft volumes in proteins relate to their molecular interactions and functions. Three separate datasets are used, representing enzyme-ligand binding, protein-protein dimerization and antibody-antigen complexes. We find that, in single-chain enzymes, the ligand is bound in the largest cleft in over 83% of the proteins. Usually the largest cleft is considerably larger than the others, suggesting that size is a functional requirement. Thus, in many cases, the likely active sites of an enzyme can be identified using purely geometrical criteria alone. In other cases, where there is no predominantly large cleft, chemical interactions are required for pinpointing the correct location. In antibody-antigen interactions the antibody usually presents a large cleft for antigen binding. In contrast, protein-protein interactions in homodimers are characterized by approximately planar interfaces with several clefts involved. However, the largest cleft in each subunit still tends to be involved.

AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR.

The AQUA and PROCHECK-NMR programs provide a means of validating the geometry and restraint violations of an ensemble of protein structures solved by solution NMR. The outputs include a detailed breakdown of the restraint violations, a number of plots in PostScript format and summary statistics. These various analyses indicate both the degree of agreement of the model structures with the experimental dat, and the quality of their geometrical properties. They are intended to be of use both to support ongoing NMR structure determination and in the validation of the final results.

Derivation of 3D coordinate templates for searching structural databases: application to Ser-His-Asp catalytic triads in the serine proteinases and lipases.

It is well established that sequence templates (e.g., PROSITE) and databases are powerful tools for identifying biological function and tertiary structure for an unknown protein sequence. Here we describe a method for automatically deriving 3D templates from the protein structures deposited in the Brookhaven Protein Data Bank. As an example, we describe a template derived for the Ser-His-Asp catalytic triad found in the serine proteases and triacylglycerol lipases. We find that the resultant template provides a highly selective tool for automatically differentiating between catalytic and noncatalytic Ser-His-Asp associations. When applied to nonproteolytic proteins, the template picks out two "non-esterase" catalytic triads that may be of biological relevance. This suggests that the development of databases of 3D templates, such as those that currently exist for protein sequence templates, will help identify the functions of new protein structures as they are determined and pinpoint their functionally important regions.

X-SITE: use of empirically derived atomic packing preferences to identify favourable interaction regions in the binding sites of proteins.

A new empirically based method for predicting favourable interaction regions within the binding sites of proteins is presented. The method uses spatial distributions of atomic contact preferences derived from a non-homologous dataset of 83 high-resolution protein structures. The contact preferences are obtained for 26 different atom types relative to 163 different types of three-atom fragments. Each fragment consists of a triplet of bonded atoms, 1-2-3, which defines a reference frame for the three-dimensional distributions. In this way, directional, as well as distance, information is retained. Once derived, the distribution can be applied in a predictive manner. Given a protein's binding site, each distribution is transformed on to the three-atom fragments of the constituent residues and, when combined, can identify the favourable interaction regions for each different atom type. These predicted regions can then form the basis either for the modification of known inhibitors or for the search and design of new ones. Five known protein-ligand complexes are used to demonstrate the validity and usefulness of the approach. The results show that the method provides a powerful tool both in understanding how a given ligand exploits the interactions available to it in an active site and in helping to design improved, or novel, protein ligands.

Conformational analysis of pentapeptide sequences matching a proposed recognition motif for lysosomal degradation.

A selective pathway for the degradation of specific long-lived cytosolic proteins is activated in response to starvation in vivo or to serum withdrawal from cultured cells. It involves recognition of a targeting motif by a member of the hsp70 family. A 5-residue targeting motif has been proposed on the basis of sequence comparisons. We investigate whether there is any structural basis for this motif being the true recognition signal. We examine the conformations of four motif peptides in proteins that are either known to be serum regulated or are from related vertebrate species, and two equivalent peptides in bacterial proteins that closely resemble other regulated proteins. Our studies show that all the motif sequences are located near the ends of surface helices with one or more of the residues buried in the structure, yet it is known that members of the hsp70 family tend to interact with extended peptide chains. Furthermore, recognition by these proteins generally requires a specific ordering of key residues, yet the motif implies a largely order-independent sequence characterized by residue type only. We conclude that the proposed motif is unlikely to be the true targeting signal for lysosomal degradation unless additional factors apply.