Search PubMed⌕ Search

Biomedical subjects

J M Thornton

Publications and source records attributed to J M Thornton.

At least 91 records · Page 5Linked to original sources
PubMed1997

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↗
PubMed1997

Structures of N-termini of helices in proteins.

We have surveyed 393 N-termini of alpha-helices and 156 N-termini of 3(10)-helices in 85 high resolution, non-homologous protein crystal structures for N-cap side-chain rotamer preferences, hydrogen bonding patterns, and solvent accessibilities. We find very strong rotamer preferences that are unique to N-cap sites. The following rules are generally observed for N-capping in alpha-helices: Thr and Ser N-cap side chains adopt the gauche - rotamer, hydrogen bond to the N3 NH and have psi restricted to 164 +/- 8 degrees. Asp and Asn N-cap side chains either adopt the gauche - rotamer and hydrogen bond to the N3 NH with psi = 172 +/- 10 degrees, or adopt the trans rotamer and hydrogen bond to both the N2 and N3 NH groups with psi = 1-7 +/- 19 degrees. With all other N-caps, the side chain is found in the gauche + rotamer so that the side chain does not interact unfavorably with the N-terminus by blocking solvation and psi is unrestricted. An i, i + 3 hydrogen bond from N3 NH to the N-cap backbone C = O in more likely to form at the N-terminus when an unfavorable N-cap is present. In the 3(10)-helix Asn and Asp remain favorable N-caps as they can hydrogen bond to the N2 NH while in the trans rotamer; in contrast, Ser and Thr are disfavored as their preferred hydrogen bonding partner (N3 NH) is inaccessible. This suggests that Ser is the optimum choice of N-cap when alpha-helix formation is to be encouraged while 3(10)-helix formation discouraged. The strong energetic and structural preferences found for N-caps, which differ greatly from positions within helix interiors, suggest that N-caps should be treated explicitly in any consideration of helical structure in peptides or proteins.

Hydrogen Bonding↗
PubMed1997

TESS: a geometric hashing algorithm for deriving 3D coordinate templates for searching structural databases. Application to enzyme active sites.

It is well established that sequence templates such as those in the PROSITE and PRINTS databases are powerful tools for predicting the biological function and tertiary structure for newly derived protein sequences. The number of X-ray and NMR protein structures is increasing rapidly and it is apparent that a 3D equivalent of the sequence templates is needed. Here, we describe an algorithm called TESS that automatically derives 3D templates from structures deposited in the Brookhaven Protein Data Bank. While a new sequence can be searched for sequence patterns, a new structure can be scanned against these 3D templates to identify functional sites. As examples, 3D templates are derived for enzymes with an O-His-O "catalytic triad" and for the ribonucleases and lysozymes. When these 3D templates are applied to a large data set of nonidentical proteins, several interesting hits are located. This suggests that the development of a 3D template database may help to identify the function of new protein structures, if unknown, as well as to design proteins with specific functions.

Algorithms↗
PubMed1997

Modelling protein unfolding: hen egg-white lysozyme.

A novel modelling procedure, which rapidly unfolds a protein by enhancing solvent penetration of its core, was used to investigate the unfolding pathway of hen egg-white lysozyme. Early on the unfolding pathway there is a dramatic disruption of the tertiary contacts within the protein, which decouples its domains. Subsequently, the helical domain slowly loses its compactness and the helices fluctuate rapidly. The protein then adopts a 'molten globule-like' structure in which the native beta-sheet is essentially intact. The modelled structures have properties similar to those of lysozyme's experimentally characterized partially folded states and provide insight into its complex (un)folding process. The sequence of unfolding events shows how the unfolding pathway of a multidomain protein may be most similar to its fastest, but not necessarily its dominant, folding pathway.

Animals↗
PubMed1996

Deviations from planarity of the peptide bond in peptides and proteins.

The work described here is the result of a survey of the peptide omega angles in the Cambridge Structural Database of small molecules, which was carried out to establish "ideal" or "target" values for their distribution. We have shown that substantial deviations from planarity can be tolerated with a standard deviation in the angle of up to 6 degrees about a mean value for the trans peptide that is less than 180 degrees . The distortion can arise from pyramidalization at the amino nitrogen atom as well as simple twist about the peptide bond. We include an analysis of omega angles in the existing database of protein structure (PDB) and show that their distributions can depend on the refinement method used, but no correlation with resolution is evident. A surprising finding was a systematic variation of omega in phi,psi space in proteins as well as in the linear and cyclic peptides. This is particularly manifest as a consistent difference between the mean omega values in chains of left and right-hand chirality. This dichotomy is observed for all the standard amino acids and is especially striking in the absence of secondary structure. The phenomenon is discussed in the context of theoretical work on peptide analogues, and the implications for protein conformation and structure are briefly considered.

Amino Acids↗
PubMed1996

Multiple solution conformations of the integrin-binding cyclic pentapeptide cyclo(-Ser-D-Leu-Asp-Val-Pro-). Analysis of the (phi, psi) space available to cyclic pentapeptides.

The aqueous solution structure of the cyclic pentapeptide cyclo(-Ser-D-Leu-Asp-Val-Pro-) has been determined by two-dimensional 1H-NMR spectroscopy, combined with a conformational search and distance-geometry calculations. As many as five conformers in slow exchange were observed, and the rate of interconversion between components was measured from the build-up rates of exchange peaks. NMR data allowed the structures of the two predominant conformers to be determined. The major component (66%) contained a cis-proline as part of a type-VIa2 beta-turn encompassing residues Asp-Val-cis-Pro-Ser. The second component (16%) contained only trans-amide bonds, and a type-VIII beta-turn formed by residues Val-Pro-Ser-D-Leu. These structures are discussed in relation to the (phi, psi), space available to the cyclic pentapeptide, determined by a conformational search, and in relation to previously published cyclic-pentapeptide structures. The molecule exhibits activity in a scintillation-proximity assay for the inhibition of the interaction between the integrin very-late antigen-4 (VLA-4; alpha 4 beta 1) and vascular-cell-adhesion molecule-1 (VCAM-1). The structure/activity relationship of the LDV sequence is discussed and related to the recently published X-ray structure of VCAM-1. The relevance of the work to the design of anti-inflammatory drugs is discussed.

Amino Acid Sequence↗
PubMed1996

Structural families in loops of homologous proteins: automatic classification, modelling and application to antibodies.

Loop regions of polypeptide in homologous proteins may be classified into structural families. A method is described by which this classification may be performed automatically and "key residue" templates, which may be responsible for the loop adopting a given conformation, are defined. The technique has been applied to the hypervariable loops of antibodies and the results are compared with the previous definition of canonical classes. We have extended these definitions and provide complete sets of structurally determining residues (SDRs) for the observed clusters including the first set of key residues for seven-residue CDR-H3 loops.

Antibodies↗
PubMed1996

Protein recognition of adenylate: an example of a fuzzy recognition template.

The interaction between protein and adenylate in a non-homologous dataset of 18 high-resolution protein/nucleotide crystal structures is analysed. We find that each constituent of adenylate, adenine, ribose and phosphate, is substantially buried. Adenine has a largely hydrophobic protein interface, while phosphate interacts primarily with hydrophilic residues; ribose is intermediate. A detailed study of hydrogen bonding in these complexes shows hydrogen bonds between protein and adenine to be surprisingly scarce. There does not seem to be a conserved hydrogen-bonding pattern for adenine recognition. The hydrogen bonds that are seen have geometries close to energy minima found in our Distributed Multipole Analysis based model calculations. The experimental hydrogen-bonded geometries have a characteristic signature in our model energy calculations, with a dominant attractive electrostatic term. For stacked interactions, however, the dispersion energy dominates. Finally, we present the concept of a fuzzy recognition template, as a useful means of describing the protein/adenylate interactions presented here, which will also be a valuable concept for characterising other protein/ligand interactions.

Adenine↗
PubMed1996

Antibody-antigen interactions: contact analysis and binding site topography.

We have analysed antigen-contacting residues and combining site shape in the antibody Fv and Fab crystal structures now available from the Protein Data Bank. Antigen-contacting propensities are presented for each antibody residue, allowing a new definition for the complementarity determining regions (CDRs) to be proposed based on observed antigen contacts. Contacts are more common at CDR residues which are located centrally within the combining site; some less central CDR residues are only contacted by large antigens. Non-contacting residues within the CDRs coincide with residues identified by Chothia and co-workers as important in defining "canonical" conformations. An objective means of classifying protein surfaces by gross topography has been developed and applied to the antibody combining site surfaces. The surfaces have been clustered into four topographic classes: concave and moderately concave (mostly hapten binders), ridged (mostly peptide binders) and planar (mostly protein binders). We have determined the topographic classes for ten pairs of complexed and uncomplexed antibody-antigen crystal structures; four change topographic class on complexation. The results will be of use in antibody engineering, antigen docking and in clinical immunology. To demonstrate one application, we show how the data can be used to locate the antigen binding pocket on antibody structures.

Antibodies↗
PubMed1996

Analysis of domain structural class using an automated class assignment protocol.

The extent to which the contemporary dataset of protein structures can be segregated into four structural "classes" as originally defined by Levitt & Chothia in 1976 is examined and a simple method presented for the assignment of protein domains into these classes. Assignments are based on known three-dimensional structures, and for successful assignment it was found that helix/sheet content, contacts between secondary structures and their sequential order had to be used. The procedure attempts to maximise the automatic separation into classes for a dataset of 197 manually classified, non-homologous domains. It was found that approximately 90% of the structures were classified automatically; the remainder were borderline and were left for manual inspection. The method was then applied to a test set of 43 protein domains with similar results. The data support the concept of distinct classes of protein structure, although a few intermediate structures are found, demonstrating that it is possible to define relatively simple parameters complying with commonly accepted nomenclature that automatically define 90% of protein domains with essentially 100% accuracy. However, re-examination of the data also suggested that the previously separate alpha/beta and alpha + beta classes show considerable overlap and are more naturally represented as a single alpha beta class. This large alpha beta class can then be most easily subdivided by consideration of whether the sheets are mainly parallel, antiparallel or mixed. The correlation between structural class and function is discussed, together with the conservation of class within a sequence superfamily. This represents the first step in an automated phenetic description of protein structure complementing the usual phylogenetic approach to protein structure classification.

Algorithms↗
PubMed1996

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↗
PubMed1996

Analysis of main chain torsion angles in proteins: prediction of NMR coupling constants for native and random coil conformations.

Using a data base of 85 high resolution protein crystal structures the distributions of main chain torsion angles, both in secondary structure and in coil regions where no secondary structure is present, have been analysed. These torsion angle distributions have been used to predict NMR homonuclear and heteronuclear coupling constants for residues in secondary structure using known Karplus relationships. For alpha helices, 3(10) helices and beta strands mean predicted 3JHN alpha coupling constants are 4.8, 5.6 and 8.5 Hz, respectively. These values differ significantly from those expected for the ideal phi angles (3.9, 3.0 and 8.9 Hz; phi = -57 degrees, -49 degrees, -139 degrees for alpha and 3(10) helices and beta strands (antiparallel), respectively) in regular secondary structure, but agree well with available experimental NMR data for nine proteins. The crystallographic data set has also been used to provide a basis for interpreting coupling constants measured for peptides and denatured proteins. Using a model for a random coil, in which all residues adopt distributions of phi, psi angles equivalent to those seen for residues in the coil regions of native folded proteins, predicted 3JHN alpha values for different residue types have been found to range from 5.9 Hz and 6.1 Hz for glycine and alanine, respectively, to 7.7 Hz for valine. A good correlation has been found between the predicted 3JHN alpha coupling constants for this model and experimental values for a set of peptides that other evidence suggest are highly unstructured. For other peptides, however, deviations from the predictions of the model are clear and provide evidence for additional interactions within otherwise disordered states. The values of homonuclear and heteronuclear coupling constants derived from the protein data base listed here therefore provide a basis not only for analysing the secondary structure of native proteins in solution but for assessing and interpreting the extent of structure present in peptides and non-native states of proteins.

Databases, Factual↗
PubMed1996

Principles of protein-protein interactions.

This review examines protein complexes in the Brookhaven Protein Databank to gain a better understanding of the principles governing the interactions involved in protein-protein recognition. The factors that influence the formation of protein-protein complexes are explored in four different types of protein-protein complexes--homodimeric proteins, heterodimeric proteins, enzyme-inhibitor complexes, and antibody-protein complexes. The comparison between the complexes highlights differences that reflect their biological roles.

Animals↗
PubMed1996

Towards meeting the Paracelsus Challenge: The design, synthesis, and characterization of paracelsin-43, an alpha-helical protein with over 50% sequence identity to an all-beta protein.

In response to the Paracelsus Challenge (Rose and Creamer, Proteins, 19:1-3, 1994), we present here the design, synthesis, and characterization of a helical protein, whose sequence is 50% identical to that of an all-beta protein. The new sequence was derived by applying an inverse protein folding approach, in which the sequence was optimized to "fit" the new helical structure, but constrained to retain 50% of the original amino acid residues. The program utilizes a genetic algorithm to optimize the sequence, together with empirical potentials of mean force to evaluate the sequence-structure compatibility. Although the designed sequence has little ordered (secondary) structure in water, circular dichroism and nuclear magnetic resonance data show clear evidence for significant helical content in water/ethylene glycol and in water/methanol mixtures at low temperatures, as well as melting behavior indicative of cooperative folding. We believe that this represents a significant step toward meeting the Paracelsus Challenge.

Algorithms↗
PubMed1996

PROMOTIF--a program to identify and analyze structural motifs in proteins.

We describe a suite of programs, PROMOTIF, that analyzes a protein coordinate file and provides details about the structural motifs in the protein. The program currently analyzes the following structural features: secondary structure; beta-and gamma-turns; helical geometry and interactions; beta-strands and beta-sheet topology; beta-bulges; beta-hairpins; beta-alpha-beta units and psi-loops; disulphide bridges; and main-chain hydrogen bonding patterns. PROMOTIF creates postscript files showing the examples of each type of motif in the protein, and a summary page. The program can also be used to compare motifs in a group of related structures, such as an ensemble of NMR structures.

Amino Acid Sequence↗
PubMed1996

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↗
PubMed1996

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↗