Computer-generated schematic diagrams of protein structures.

Computer-generated pictures are essential for studying and comparing the structures of proteins that have been solved by x-ray crystallography. Stereoscopic pairs produced by a computer program are particularly useful in providing an intelligible portrayal of the molecular topology.

The predicted structure of immunoglobulin D1.3 and its comparison with the crystal structure.

Predictions of the structures of the antigen-binding domains of an antibody, recorded before its experimental structure determination and tested subsequently, were based on comparative analysis of known antibody structures or on conformational energy calculations. The framework, the relative positions of the hypervariable regions, and the folds of four of the hypervariable loops were predicted correctly. This portion includes all residues in contact with the antigen, in this case hen egg white lysozyme, implying that the main chain conformation of the antibody combining site does not change upon ligation. The conformations of three residues in each of the other two hypervariable loops are different in the predicted models and the experimental structure.

An atlas of serpin conformations.

The serpins are a family of proteins that inhibit chymotrypsin-like serine proteinases, with an unusual mechanism involving a large conformational change known as the stressed-->relaxed (S-->R) transition. This article is a guide to the known serpin conformations and their biological significance.

NAD-binding domains of dehydrogenases.

The nicotinamide adenine dinucleotide (NAD)-binding domains of dehydrogenases, containing a conserved double beta-alpha-beta-alpha-beta motif, are common structural feature of many enzymes that bind NAD, nicotinamide adenine dinucleotide phosphate (NADP) and related cofactors. Features of this folding pattern that create a natural binding site for such molecules have been described. The domain continues to appear in many structures, in the form of a common core with different peripheral additions or variations. Other structures that bind NAD and related molecules use entirely different topologies, although, in many, a phosphate group appears at the N terminus of an alpha helix. Ferredoxin reductase seems to show convergent evolution, containing a single beta-alpha-beta motif that is similar both in its structure and in its interactions with the ligand to a region in dehydrogenases.

Mechanisms of domain closure in proteins.

Certain enzymes respond to the binding of substrates and coenzymes by the closure of an active site that lies in a cleft between two domains. We have examined the mechanism of the domain closure in citrate synthase, for which atomic co-ordinates are available for "open" and "closed" forms. We show that the mechanism of domain closure involves small shifts and rotations of packed helices within the two domains and at their interface. Large motions of distant segments of the structure are the cumulative effect of the small relative shifts in intervening pairs of packed segments. These shifts are accommodated not by changes in packing but rather by small conformational changes in side-chains. We call this the helix interface shear mechanism of domain closure. The relative movements of packed helices follow the principles suggested by our recent study of insulin. This mechanism of domain closure is quite different from the hinge mechanisms that allow the rigid body movements of domains in immunoglobulins. The large interface between the domains of citrate synthase precludes a simple hinge mechanism for its conformational change. The helix interface shear mechanism of conformational change occurs in other enzymes that contain extensive domain-domain interfaces.

Helix movements and the reconstruction of the haem pocket during the evolution of the cytochrome c family.

Analysis of cytochromes c (tuna), c2 (Rhodospirillum rubrum), c550 (Paracoccus denitrificans) and c551 (Pseudomonas aeruginosa) shows that they contain 48 residues identifiable as homologous from superposition of the structures. The other 34 to 64 residues are in loops that vary greatly in sequence, length and conformation, or in alpha-helices that are found in only some of the structures. Of the 48 homologous residues, 17 are in three segments which pack onto the haem faces. In all four structures, these segments have the same conformations, and the same locations relative to the haem. The other 31 residues are in three alpha-helices which are in contact with each other. These form the back and one side of the haem pocket. In cytochrome c551 the positions of the three alpha-helices have shifted and rotated, in comparison with cytochromes c and c2, by up to 5 A and 25 degrees relative to the haem. These shifts, facilitated by mutations at the helix-helix interfaces, are related to the reconstruction of the propionic acid side of the haem pocket described by Almassy & Dickerson (1978). Together these effects produce alternative structures for the haem pocket. This mechanism of adaptation to mutation contrasts with that observed in the globins. In the globins, mutations also produce changes in helix interfaces and shifts of packed helices, but in the globins these shifts are coupled to conserve the structure of the haem pocket.

Canonical structures for the hypervariable regions of immunoglobulins.

We have analysed the atomic structures of Fab and VL fragments of immunoglobulins to determine the relationship between their amino acid sequences and the three-dimensional structures of their antigen binding sites. We identify the relatively few residues that, through their packing, hydrogen bonding or the ability to assume unusual phi, psi or omega conformations, are primarily responsible for the main-chain conformations of the hypervariable regions. These residues are found to occur at sites within the hypervariable regions and in the conserved beta-sheet framework. Examination of the sequences of immunoglobulins of unknown structure shows that many have hypervariable regions that are similar in size to one of the known structures and contain identical residues at the sites responsible for the observed conformation. This implies that these hypervariable regions have conformations close to those in the known structures. For five of the hypervariable regions, the repertoire of conformations appears to be limited to a relatively small number of discrete structural classes. We call the commonly occurring main-chain conformations of the hypervariable regions "canonical structures". The accuracy of the analysis is being tested and refined by the prediction of immunoglobulin structures prior to their experimental determination.

Conservation and variability in the structures of serine proteinases of the chymotrypsin family.

The chymotrypsin-like serine proteinases are a widely divergent family of enzymes appearing in animals, bacteria and viruses. All comprise two homologous domains containing six-stranded beta-barrels, with the active site between the domains. What are the structural constraints to which these proteins have been subject during their evolution, and how have the molecules explored the limits these constraints impose? We have analysed the structures of 13 widely divergent serine proteinases determined by X-ray crystallography to high resolution and well refined. We have identified the regions of the molecule that evolution has preserved both within each of the two domains and in the interdomain interface. An alignment of the sequences is presented based on structural superpositions. From this we analyse the conserved structural patterns and the interactions crucial to maintaining the structure and the active side.

Standard conformations for the canonical structures of immunoglobulins.

A comparative analysis of the main-chain conformation of the L1, L2, L3, H1 and H2 hypervariable regions in 17 immunoglobulin structures that have been accurately determined at high resolution is described. This involves 79 hypervariable regions in all. We also analysed a part of the H3 region in 12 of the 15 VH domains considered here. On the basis of the residues at key sites the 79 hypervariable regions can be assigned to one of 18 different canonical structures. We show that 71 of these hypervariable regions have a conformation that is very close to what can be defined as a "standard" conformation of each canonical structure. These standard conformations are described in detail. The other eight hypervariable regions have small deviations from the standard conformations that, in six cases, involve only the rotation of a single peptide group. Most H3 hypervariable regions have the same conformation in the part that is close to the framework and the details of this conformation are also described here.

beta-Trefoil fold. Patterns of structure and sequence in the Kunitz inhibitors interleukins-1 beta and 1 alpha and fibroblast growth factors.

Previous crystallographic analyses of the Kunitz inhibitors from soybean. Erythrina caffra and wheat, the interleukins-1 beta and 1 alpha and the acidic and basic fibroblast growth factors have shown that they contain a most unusual fold. It is formed by six two-stranded hairpins. Three of these form a barrel structure and the other three are in a triangular array that caps the barrel. The arrangement of the secondary structures gives the molecules a pseudo 3-fold axis. Although the different proteins have very similar structures, many of their sequences have no significant similarities overall. The structural determinants of this fold are described and discussed in this paper. The barrels in the different proteins have the same geometrical features: six strands tilted at 56 degrees to the barrel axis; a barrel diameter of 16 A, and the beta-sheet hydrogen bonded so that it is staggered with a shear number of 12. These features fit McLachlan's equations for ideal barrels formed by beta-sheets. The wide diameter of the barrels is filled by layers of residues that, while not identical in the different proteins, are, in almost all cases, large. The structure of the triangular array of hairpins is determined by the coiling of the strands and the packing of hairpin residues against each other and against residues from the interior of the barrel. The major sequence requirements of this fold are large or medium hydrophobic residues at 18 buried sites. In the different structures the total volume of these residues is 3000 (+/- 120) A. The polyhedron model of protein architecture is used to demonstrate that the main, and in particular the symmetrical, features of this fold arise from the ideal and equal packing of six hairpins, modified only slightly to form hydrogen bonds between the hairpins.

Determinants of a protein fold. Unique features of the globin amino acid sequences.

The three-dimensional structures of globins are known, from crystallographic analyses, to be very similar. Their amino acid sequences, however, differ greatly. Only two residues are absolutely conserved in all sequences, and the residue identities of some pairs of sequences are only 16%. We have determined the nature and exact extent of the sequence variations and the extent to which the conserved features of the globin sequences are unique to this family. The 226 globin sequences now known were aligned and analysed. Because distantly related protein sequences cannot be aligned correctly without the use of structural data, we developed a method that incorporated structural information into the alignment procedure. Analysis of the aligned sequences show that: (1) Although individual chains vary in size between 132 and 157 residues, deletions and insertions result in there being only 102 residue sites common to all globins. These sites form six separate regions. Insertions and deletions between these regions means that their separations can vary in different sequences. (2) Within the conserved regions there are 32 sites that almost always contain hydrophobic residues. In the known structures, these sites are in the protein interior. We measured the variations in the size of the residues that occur in the 226 sequences at these sites. At six sites the residues differ in size by less than 40 A3, at 11 sites they differ by 40 to 100 A3, and at 15 sites they differ by more than 100 A3. There are two other conserved buried sites: one contains the His linked to the haem iron and the other usually contains a His involved with the haem ligand. (3) Within the conserved regions there are another 32 sites that are almost always occupied by charged, polar or small non-polar (Gly or Ala) residues. In the known structures, these sites are on the protein surface. To determine the extent to which the conserved features found for the globin sequences are unique to that protein family, the following procedure was used. The six conserved regions, and the residue restrictions that occur at the 66 sites within these regions, were encoded into two "templates". One was based only on the sequences so far determined; the other was extended to include as yet unobserved substitutions that seemed plausible on the basis of size, hydrophobicity and polarity. Each of the 3286 non-globin sequences in the data bank was then examined by a computer program to see how closely it could be matched to these templates.(ABSTRACT TRUNCATED AT 400 WORDS)

Principles determining the structure of beta-sheet barrels in proteins. I. A theoretical analysis.

The major feature of many proteins is a large beta-sheet that twists and coils to form a closed structure in which the first strand is hydrogen bonded to the last: the beta-sheet barrel. McLachlan classified barrels in terms of two integral parameters: the number of strands in the beta-sheet, n, and the "shear number", S, a measure of the stagger of the strands in the beta-sheet. He showed that the mean radius of a barrel and the extent to which strands are tilted relative to its axis are determined by the values of n and S. Here we show that the (n, S) values determine all the other general structural features of regular beta-sheet barrels, in particular, optimal values of the twist and coiling angles that produce the closed beta-sheet, the hyperboloidal shape and the arrangement of residues in the barrel interior. Consideration of the residue arrangements in the interiors of different potential barrel structures, and of side-chain volumes, suggest that barrels, in which the interiors are close packed by the residues in beta-sheets with good geometries, have structures that correspond to one of only ten different combinations of n and S. In the accompanying paper, we demonstrate, by an analysis of all observed protein structures that contain beta-sheet barrels and for which atomic co-ordinates are available, the validity of these theoretical results.

Principles determining the structure of beta-sheet barrels in proteins. II. The observed structures.

In the accompanying paper we derived a set of principles that, we argue, govern the structure of beta-sheet barrels. Barrel structures are classified in terms of two integral parameters: the number of strands in the beta-sheet, n, and a measure of the stagger in the beta-sheet, S. We derived a set of equations that show how the (n, S) values of a barrel structure determine the arrangement of its strands; its general shape; the twist and coiling of the beta-sheet, and the arrangement of residues in the barrel interior. This work suggested that there are ten different combinations of n and S that form barrels with good beta-sheet geometries and interiors close packed by beta-sheet residues. In this paper we demonstrate the validity of these principles. We analyse in detail the observed structures of 39 different beta-sheet barrels. These structures include representatives of all the different barrel structures currently known and for which atomic co-ordinates are available. We show that the observed arrangement of the strands, and the extent of the twist and coiling of the beta-sheets, are very close to those calculated from the (n, S) values for the barrel. Of the 39 structures, 34 have one of the ten (n, S) values that we expect to form barrels with good beta-sheet geometries and interiors close packed by beta-sheet residues. The other five have one of two (n, S) values that give good beta-sheet geometries but radii so large the beta-sheet residues leave cavities at the centre of the barrels. In at least four of these cavities have a functional role.

Canonical structures for the hypervariable regions of T cell alphabeta receptors.

T cell alphabeta receptors have binding sites for peptide-MHC complexes formed by six hypervariable regions. Analysis of the six atomic structures known for Valpha and for Vbeta domains shows that their first and second hypervariable regions have one of three or four different main-chain conformations (canonical structures). Six of these canonical structures have the same conformation in complexes with peptide-MHC complexes, the free receptor and/or in an isolated V domain. Thus, for at least the first and second hypervariable regions in the currently known structures, the conformation of the canonical structures is well defined in the free state and is conserved on formation of complexes with peptide-MHC. We identified the key residues that are mainly responsible for the conformation of each canonical structure. The first and second hypervariable regions of Valpha and Vbeta domains are encoded by the germline V segments. Humans have 37 functional Valpha segments and 47 Vbeta segments, and mice have 20 Vbeta segments. Inspection of the size of their hypervariable regions, and of sites that contain key residues, indicates that close to 70 % of Valpha segments and 90 % of Vbeta segments have hypervariable regions with a conformation of one of the known canonical structures. The alpha and beta V gene segments in both humans and mice have only a few combinations of different canonical structure in their first and second hypervariable regions. In human Vbeta domains, the number of different sequences with these canonical structure combinations is larger than in mice, whilst for Valpha domains it is probably smaller.

Framework residue 71 is a major determinant of the position and conformation of the second hypervariable region in the VH domains of immunoglobulins.

Analysis of the immunoglobulins of known structure reveals systematic differences in the position and main-chain conformation of the second hypervariable region of the VH domain (H2). We show that the major determinant of the position of H2 is the size of the residue at site 71, a site that is in the conserved framework of the VH domain. It is likely that for about two thirds of the known VH sequences the size of the residue at this site is also a major determinant of the conformation of H2. This effect can override the predisposition of the sequence, as in the case of the H2 loop of J539, which is an exception to the rules relating sequence and conformation of short hairpin loops. Understanding the relationship between the residue at position 71 and the position and conformation of H2 has applications to the prediction and engineering of antigen-binding sites of immunoglobulins.

Correlation of co-ordinated amino acid substitutions with function in viruses related to tobacco mosaic virus.

Sequence data are available for the coat proteins of seven tobamoviruses, with homologies ranging from at least 26% to 82%, and atomic co-ordinates are known for tobacco mosaic virus (TMV) vulgare. A significant spatial relationship has been found between groups of residues with identical amino acid substitution patterns. This strongly suggest that their location is linked to a particular function, at least in viruses identical with the wild-type for these residues. The most conserved feature of TMV is the RNA binding region. Core residues are conserved in all viruses or show mutations complementary in volume. The specificity of inter-subunit contacts is achieved in different ways in the three more distantly related viruses.

Interior and surface of monomeric proteins.

The solvent-accessible surface area (As) of 46 monomeric proteins is calculated using atomic co-ordinates from high-resolution and well-refined crystal structures. The As of these proteins can be determined to within 1 to 2% and that of their individual residues to within 10 to 20%. The As values of proteins are correlated with their molecular weight (Mr) in the range 4000 to 35,000: the power law As = 6.3 M0.73 predicts protein As values to within 4% on average. The average water-accessible surface is found to be 57% non-polar, 24% polar and 19% charged, with 5% root-mean-square variations. The molecular surface buried inside the protein is 58% non-polar, 39% polar and 4% charged. The buried surface contains more uncharged polar groups (mostly peptides) than the surface that remains accessible, but many fewer charged groups. On average, 15% of residues in small proteins and 32% in larger ones may be classed as "buried residues", having less than 5% of their surface accessible to the solvent. The accessibilities of most other residues are evenly distributed in the range 5 to 50%. Although the fraction of buried residues increases with molecular weight, the amino acid compositions of the protein interior and surface show no systematic variation with molecular weight, except for small proteins that are often very rich in buried cysteines. From amino acid compositions of protein surfaces and interiors we calculate an effective coefficient of partition for each type of residue, and derive an implied set of transfer free energy values. This is compared with other sets of partition coefficients derived directly from experimental data. The extent to which groups of residues (charged, polar and non-polar) are buried within proteins correlates well with their hydrophobicity derived from amino acid transfer experiments. Within these three groups, the correlation is low.

Haemoglobin: the surface buried between the alpha 1 beta 1 and alpha 2 beta 2 dimers in the deoxy and oxy structures.

Using the newly available refined co-ordinates of deoxy and oxyhaemoglobin, we have re-examined and compared the interfaces between the dimers alpha 1 beta 1 and alpha 2 beta 2. The most extensive monomer-monomer contacts are between alpha 1 and beta 2, and, symmetrically, alpha 2 and beta 1. In oxyhaemoglobin these interfaces bury 700 A2 less protein surface than in deoxyhaemoglobin. The alpha 1 alpha 2 interface involves similar salt bridges in both forms, but in oxyhaemoglobin buries 240 A2 more surface than in deoxyhaemoglobin. There is a loosely packed beta 1 beta 2 interface burying 320 A2 of surface in oxyhaemoglobin; there is no beta 1 beta 2 interface in deoxyhaemoglobin. The greater stability of the deoxy form, in the absence of ligands, can be attributed to a combination of hydrophobic, van der Waals' and electrostatic interactions.

Conformational changes in serpins: I. The native and cleaved conformations of alpha(1)-antitrypsin.

The serpins (SERine Proteinase INhibitors) are a family of proteins with important physiological roles, including but not limited to the inhibition of chymotrypsin-like serine proteinases. The inhibitory mechan- ism involves a large conformational change known as the S-->R (stressed-->relaxed) transition. The largest structural differences occur in a region around the scissile bond called the reactive centre loop: In the native (S) state, the reactive centre is exposed, and is free to interact with proteinases. In inhibitory serpins, in the cleaved (R) state the reactive centre loop forms an additional strand within the beta-sheet. The latent state is an uncleaved state in which the intact reactive centre loop is integrated into the A sheet as in the cleaved form, to give an alternative R state. The serpin structures illustrate detailed control of conformation within a single protein. Serpins are also an unusual family of proteins in which homologues have native states with different folding topologies. Determination of the structures of inhibitory serpins in multiple conformational states permits a detailed analysis of the mechanism of the S-->R transition, and of the way in which a single sequence can form two stabilised states of different topology. Here we compare the conformations of alpha(1)-antitrypsin in native and cleaved states. Many protein conformational changes involve relative motions of large rigid subunits. We determine the rigid subunits of alpha(1)-antitrypsin and analyse the changes in their relative position and orientation. Knowing that the conformational change is initiated by cleavage at the reactive centre, we describe a mechanism of the S-->R transition as a logical sequence of mechanical effects, even though the transition likely proceeds in a concerted manner.

Conformational changes in serpins: I. The native and cleaved conformations of alpha(1)-antitrypsin.

The serpins (SERine Proteinase INhibitors) are a family of proteins with important physiological roles, including but not limited to the inhibition of chymotrypsin-like serine proteinases. The inhibitory mechan- ism involves a large conformational change known as the S-->R (stressed-->relaxed) transition. The largest structural differences occur in a region around the scissile bond called the reactive centre loop: In the native (S) state, the reactive centre is exposed, and is free to interact with proteinases. In inhibitory serpins, in the cleaved (R) state the reactive centre loop forms an additional strand within the beta-sheet. The latent state is an uncleaved state in which the intact reactive centre loop is integrated into the A sheet as in the cleaved form, to give an alternative R state. The serpin structures illustrate detailed control of conformation within a single protein. Serpins are also an unusual family of proteins in which homologues have native states with different folding topologies. Determination of the structures of inhibitory serpins in multiple conformational states permits a detailed analysis of the mechanism of the S-->R transition, and of the way in which a single sequence can form two stabilised states of different topology. Here we compare the conformations of alpha(1)-antitrypsin in native and cleaved states. Many protein conformational changes involve relative motions of large rigid subunits. We determine the rigid subunits of alpha(1)-antitrypsin and analyse the changes in their relative position and orientation. Knowing that the conformational change is initiated by cleavage at the reactive centre, we describe a mechanism of the S-->R transition as a logical sequence of mechanical effects, even though the transition likely proceeds in a concerted manner.