PubMed HealthSearch

Biomedical subjects

R A Laskowski

Publications and source records attributed to R A Laskowski.

At least 19 recordsLinked to original sources

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.

Amino Acid Sequence

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.

Binding Sites

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.

Algorithms

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.

Animals

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.

Crystallography, X-Ray

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.

Amino Acid Sequence

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.

Amino Acid Sequence

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.

Computer Graphics

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.

Amides

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.

Algorithms

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.

Amino Acid Sequence

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.

Amino Acid Isomerases

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.

Animals

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.

Amino Acid Sequence

SURFNET: a program for visualizing molecular surfaces, cavities, and intermolecular interactions.

The SURFNET program generates molecular surfaces and gaps between surfaces from 3D coordinates supplied in a PDB-format file. The gap regions can correspond to the voids between two or more molecules, or to the internal cavities and surface grooves within a single molecule. The program is particularly useful in clearly delineating the regions of the active site of a protein. It can also generate 3D contour surfaces of the density distributions of any set of 3D data points. All output surfaces can be viewed interactively, along with the molecules or data points in question, using some of the best-known molecular modeling packages. In addition, PostScript output is available, and the generated surfaces can be rendered using various other graphics packages.

Amino Acid Sequence

LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions.

The LIGPLOT program automatically generates schematic 2-D representations of protein-ligand complexes from standard Protein Data Bank file input. The output is a colour, or black-and-white, PostScript file giving a simple and informative representation of the intermolecular interactions and their strengths, including hydrogen bonds, hydrophobic interactions and atom accessibilities. The program is completely general for any ligand and can also be used to show other types of interaction in proteins and nucleic acids. It was designed to facilitate the rapid inspection of many enzyme complexes, but has found many other applications.

Algorithms