Roman A. Laskowski

PROCHECK: a program to check the stereochemical quality of protein structures

The PROCHECK suite of programs provides a detailed check on the stereochemistry of a protein structure. Its outputs comprise a number of plots in PostScript format and a comprehensive residue-by-residue listing. These give an assessment of the overall quality of the structure as compared with well refined structures of the same resolution and also highlight regions that may need further investigation. The PROCHECK programs are useful for assessing the quality not only of protein structures in the process of being solved but also of existing structures and of those being modelled on known structures.

AQUA and PROCHECK-NMR: Programs for checking the quality of protein structures solved by NMR*

SummaryThe 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 data, 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.

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–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 chains 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 proteins 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.

The implications of alternative splicing in the ENCODE protein complement.

Alternative premessenger RNA splicing enables genes to generate more than one gene product. Splicing events that occur within protein coding regions have the potential to alter the biological function of the expressed protein and even to create new protein functions. Alternative splicing has been suggested as one explanation for the discrepancy between the number of human genes and functional complexity. Here, we carry out a detailed study of the alternatively spliced gene products annotated in the ENCODE pilot project. We find that alternative splicing in human genes is more frequent than has commonly been suggested, and we demonstrate that many of the potential alternative gene products will have markedly different structure and function from their constitutively spliced counterparts. For the vast majority of these alternative isoforms, little evidence exists to suggest they have a role as functional proteins, and it seems unlikely that the spectrum of conventional enzymatic or structural functions can be substantially extended through alternative splicing.

Protein folds and functions

BACKGROUNDnThe 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.nnnRESULTSnUsing 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.nnnCONCLUSIONSnThe 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.

BLEEP—potential of mean force describing protein–ligand interactions: I. Generating potential

We have developed BLEEP (biomolecular ligand energy evaluation protocol), an atomic level potential of mean force (PMF) describing protein–ligand interactions. The pair potentials for BLEEP have been derived from high‐resolution X‐ray structures of protein–ligand complexes in the Brookhaven Protein Data Bank (PDB), with a careful treatment of homology. The use of a broad variety of protein–ligand structures in the derivation phase gives BLEEP more general applicability than previous potentials, which have been based on limited classes of complexes, and thus represents a significant step forward. We calculate the distance distributions in protein–ligand interactions for all 820 possible pairs that can be chosen from our set of 40 different atom types, including polar hydrogen. We then use a reverse Boltzmann methodology to convert these into energy‐like pair potential functions. Two versions of BLEEP are calculated, one including and one excluding interactions between protein and water. The pair potentials are found to have the expected forms; polar and hydrogen bonding interactions show minima at short range, around 3.0 Å, whereas a typical hydrophobic interaction is repulsive at this distance, with values above 4.0 Å being preferred.u2003©1999 John Wiley & Sons, Inc.u2003J Comput Chem 20: 1165–1176, 1999

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.

Error Estimates of Protein Structure Coordinates and Deviations from Standard Geometry by Full-Matrix Refinement of γB- and βB2-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. γB-crystallin, refined at 1.49u2005A yielded average errors in atomic positions which ranged from 0.05u2005A for main-chain atoms to 0.27u2005A for unrestrained water molecules. The second structure used in this work was that of βB2-crystallin refined at 2.1u2005A resolution where the corresponding average errors were 0.08 and 0.35u2005A, 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.