Strict rules determine arrangements of strands in sandwich proteins.

From a computer analysis of the spatial organization of the secondary structures of beta-sandwich proteins, we find certain sets of consecutive strands that are connected by hydrogen bonds, which we call "strandons." The analysis of the arrangements of strandons in 491 protein structures that come from 69 different superfamilies reveals strict regularities in the arrangements of strandons and the formation of what we call "canonical supermotifs." Six such supermotifs account for approximately 90% of all observed structures. Simple geometric rules are described that dictate the formation of these supermotifs.

A geometric construction determines all permissible strand arrangements of sandwich proteins.

For a large class of proteins called sandwich-like proteins (SPs), the secondary structures consist of two beta-sheets packed face-to-face, with each beta-sheet consisting typically of three to five beta-strands. An important step in the prediction of the three-dimensional structure of a SP is the prediction of its supersecondary structure, namely the prediction of the arrangement of the beta-strands in the two beta-sheets. Recently, significant progress in this direction was made, where it was shown that 91% of observed SPs form what we here call "canonical motifs." Here, we show that all canonical motifs can be constructed in a simple manner that is based on thermodynamic considerations and uses certain geometric structures. The number of these structures is much smaller than the number of possible strand arrangements. For instance, whereas for SPs consisting of six strands there exist a priori 900 possible strand arrangements, there exist only five geometric structures. Furthermore, the few motifs that are noncanonial can be constructed from canonical motifs by a simple procedure.

Prediction of the structural motifs of sandwich proteins.

We investigate the supersecondary structure of a large group of proteins, the so-called sandwich proteins. The analysis of a large number of such proteins has led us to propose a set of rules that can be used to predict the possible arrangements of strands in the two beta-sheets forming a given sandwich structure. These rules imply the existence of certain invariant supersecondary substructures common to all sandwich proteins. Furthermore, they dramatically restrict the number of permissible arrangements. For example, whereas for proteins consisting of three strands in each beta-sheet 180 possible strand arrangements exist a priori, our rules imply that only 15 of them are permissible. Five of these predicted arrangements describe all currently known sandwich proteins with six strands.

Consensus and comprehensive linkage maps of bovine chromosome 24.

This study describes development of a consensus genetic linkage map of bovine chromosome 24 (BTA24). Eight participating laboratories contributed data for 58 unique markers including a total of 25 409 meioses. Eighteen markers, which were typed in more than one reference population, were used as potential anchors to generate a consensus framework map. The framework map contained 16 loci ordered with odds greater than 1000:1 and spanned 79.3 cM. Remaining markers were included in a comprehensive map relative to these anchors. The resulting BTA24 comprehensive map was 98.3 cM in length. Average marker intervals were 6.1 and 2.5 cM for framework and comprehensive maps, respectively. Marker order was generally consistent with previously reported BTA24 linkage maps. Only one discrepancy was found when comparing the comprehensive map with the published USDA-MARC linkage map. Integration of genetic information from different maps provides a high-resolution BTA24 linkage map.

Determining the roles of different chain fragments in recognition of immunoglobulin fold.

We examine sequence-to-structure specificity of beta-structural fragments of immunoglobulin domains. The structure specificity of separate chain fragments is estimated by computing the Z-score values in recognition of the native structure in gapless threading tests. To improve the accuracy of our calculations we use energy averaging over diverse homologs of immunoglobulin domains. We show that the interactions between residues of beta-structure are more determinant in recognition of the native structure than the interactions within the whole chain molecule. This result distinguishes immunoglobulins from more typical proteins where the interactions between residues of the whole chain normally recognize the native fold more accurately than interactions between the residues of the secondary structure residues alone [Reva,B. and Topiol,S. (2000) BIOCOMPUTING: Proceedings of the Pacific Symposium. World Scientific Publishing Co., pp. 168-178]. We also find that the predominant contributions of the secondary structure are produced by the four central beta-strands that form the core of the molecule. The results of this study allow us through quantitative means to understand the architecture of immunoglobulin molecules. Comparing the fold recognition data for different chain fragments one can say that beta-strands form a rigid frame for immunoglobulin molecules, whereas loops, with no structural role, can develop a broad variety of binding specificities. It is well known that protein function is determined by specific portions of a protein chain. This study suggests that the whole protein structure can be predominantly determined by a few fragments of chain which form the structural framework of the molecule. This idea may help in better understanding the mechanisms of protein evolution: strengthening a protein structure in the key framework-forming regions allows mutations and flexibility in other chain regions.

The sequence determinants of cadherin molecules.

The sequence and structural analysis of cadherins allow us to find sequence determinants-a few positions in sequences whose residues are characteristic and specific for the structures of a given family. Comparison of the five extracellular domains of classic cadherins showed that they share the same sequence determinants despite only a nonsignificant sequence similarity between the N-terminal domain and other extracellular domains. This allowed us to predict secondary structures and propose three-dimensional structures for these domains that have not been structurally analyzed previously. A new method of assigning a sequence to its proper protein family is suggested: analysis of sequence determinants. The main advantage of this method is that it is not necessary to know all or almost all residues in a sequence as required for other traditional classification tools such as BLAST, FASTA, and HMM. Using the key positions only, that is, residues that serve as the sequence determinants, we found that all members of the classic cadherin family were unequivocally selected from among 80,000 examined proteins. In addition, we proposed a model for the secondary structure of the cytoplasmic domain of cadherins based on the principal relations between sequences and secondary structure multialignments. The patterns of the secondary structure of this domain can serve as the distinguishing characteristics of cadherins.

Comprehensive linkage map of bovine chromosome 27.

The results of genotypic data contributed to the International Society for Animal Genetics (ISAG) Bovine Chromosome 27 Workshop are presented. Eight laboratories contributed 23 261 informative meioses from 44 loci. Eighteen loci were typed by at least two laboratories and were used to construct a consensus linkage map. Twenty-one loci were subsequently incorporated into a comprehensive map. The sex-averaged consensus map covered 66.9 cM. The sex-averaged comprehensive map was 75.5 cM, while the female and male maps were 73.1 and 63.7 cM, respectively. Five loci were excluded from the analysis because of ambiguous position in the linkage group and a low LOD score (less than 2.0). Average distance between loci in the comprehensive map was 1.98 cM.

Comprehensive linkage map of bovine chromosome 11.

The results of genotypic data contributed to the International Society of Animal Genetics (ISAG) Bovine Chromosome 11 (BTA11) Workshop are presented. Six laboratories contributed a total of 26 199 informative meioses from 80 loci. Thirty-six loci were typed by at least two independent laboratories and were used to construct a consensus linkage map of the chromosome. The remaining loci were subsequently incorporated into a comprehensive map. The sex-averaged consensus map covered 128.9 cM. The female consensus map was 101.2 cM, while the male consensus map was 129.8 cM. The comprehensive sex-averaged map was 134.2 cM and the average genetic distance between loci was 1.72 cM.

Intravenous gammaglobulins in refractory polymyositis: lower dose for maintenance treatment is effective.

OBJECTIVES: To test lower dose immunoglobulins as a maintenance treatment in a patient with refractory polymyositis. METHODS: In a patient with longstanding refractory polymyositis, intravenous (IV) immunoglobulin treatment was introduced at a standard dose (2 g/kg monthly). After a few treatment courses, doses were reduced to 0.8 g/kg monthly, allowing perfusion over one single day. RESULTS: Although response to the standard dose was only partial, reduction of subsequent doses did not alter the evolution. On the contrary, the evolution was marked by further improvement, which has been sustained over the following year. CONCLUSION: Lower dose IV immunoglobulins as a maintenance treatment were used with excellent results in a case of refractory polymyositis allowing considerable reduction in treatment costs. Further trials should be undertaken to evaluate this interesting alternative.

Geometric invariant core for the CL and CH1 domains of immunoglobulin molecules.

A previously developed algorithmic method for identifying a geometric invariant of protein structures, termed geometrical core, is extended to the C(L) and C(H1) domains of immunoglobulin molecules. The method uses the matrix of C(alpha) - C(alpha) distances and does not require the usual superposition of structures. The result of applying the algorithm to 53 Immunoglobulin structures led to the identification of two geometrical core sets of C(alpha) atom positions for the C(L) and C(H1) domains.

Class-defining characteristics in the mouse heavy chains of variable domains.

Analysis of residue correlation in over 2700 mouse heavy chains of the V(H) domains was carried out on three hierarchical levels. At the 'position' level, statistical analysis revealed 45 positions that conserve similar residues in almost all chains. At the 'fragment' level, the focus of investigation shifted to the study of combinations of amino acids in strands and loops. It was found that no more than 10 patterns were sufficient for describing strands and loops in the chains. At the 'sequence' level, we determined all possible combinations of these patterns and classified the mouse heavy chains. Comparison of the sequences in the eight classes revealed residues at the class-determining positions that were unique to each class. Because a strong correlation of residues was found, one only needs several residues to classify a sequence. It follows that no all residue alignment procedure is necessary to divide sequences into classes. An important corollary of our approach is the possibility of predicting residues in an incomplete sequence from a small sequence fragment. On the basis of our analysis of mouse heavy chains we hypothesize about the presently unknown mouse V(H) germline repertoire.

Algorithmic determination of core positions in the VL and VH domains of immunoglobulin molecules.

We introduce a new algorithmic method for identifying the geometrical core of proteins that does not require the usual superposition of structures. A geometrical core is defined as the set of residues such that the C alpha (I) - C alpha (J) atom distances are identical in all structures of the protein family under study, where I and J are secondary structure positions in the structural units--strands, loops, or parts of them. The result of applying the algorithm to 53 Ig structures leads to the identification of two geometrical core sets of C alpha atom positions for the VL and VH domains. Applications of the core sets are described.

Predicting amino acid sequences of the antibody human VH chains from its first several residues.

A new method for classification of Ig sequences is suggested. The defining characteristic of a class is presence of particular residues at several class-determining positions. Sequences within a class follow the same amino acid pattern, i.e., residues at identical positions are, in an overwhelming majority of sequences of that class, identical or chemically related. Thus, once the class of a sequence is determined, one can predict the residue(s) at almost any position in the sequence. In this paper, results of analysis of 1,172 human heavy chains are presented. It was shown that a sequence can be assigned to one of six classes depending on which residues are found at its positions 1, 3, 5, 6, 7, 9, 10, 12, and 13. It is important to note that it is possible to achieve same six-class classification of the human heavy chains on the basis of a different set of positions found not at the beginning but near the end of the sequence (around position 80). For every class, an amino acid pattern of an entire sequence (complementarity determining regions excepting) has been determined. Our approach allowed us to reconstruct the incomplete human heavy chains in which residues at certain positions at the beginning or end of the chain are known. We developed a software tool for analysis, classification, and prediction of residues in sequences of the Ig family.

Structural determinants in the sequences of immunoglobulin variable domain.

To determine the general relation between the sequence and structure of variable domains of immunoglobulins we have carried out an analysis of their atomic structures, some 5300 different expressed sequences and the human germline gene segments. Variable domains are formed by two beta-sheets, packed face to face, and the inter-strand turns. Comparison of the different known structures shows that they have a core of 76 residues which has the same main-chain conformation in all structures. This common core contains almost all of the beta-sheet structure and three inter-strand turns. The regions that differ in conformation are the three hypervariable regions, three other inter-strand turns and a few adjacent residues. The 5300 expressed sequences currently known for variable domains were examined to determine the residues that occur at the 76 sites. Ignoring site conservations that occur for functional reasons, there are eight sites that have the same residue in almost all sequences; 12 that have one of a small group of very similar residues, and 52 where the chemical character of the residues is strongly conserved but not their volume. The role of residues at each site in the core was determined from the examination of their accessible surface areas, contacts, packing and buried side-chain hydrogen bonds. The most strongly conserved sites form the "deep" structure of the domain at the centre of the interface between the beta-sheets. It includes eight invariant sites and 11 sites that have one of a set of very similar residues. Around the deep structure there are buried hydrophobic residues that, in different variable domains, can differ greatly in volume. These differences in volume are accommodated by conformational changes in turn regions that are outside the common core. On the surface nearly all residues not involved in function or turn conformations strongly conserve hydrophilic or neutral residues. The implications of these results for the general relations between the sequence and structure of proteins are discussed.

Geometric invariant core for the V(L) and V(H) domains of immunoglobulin molecules.

A new algorithmic method for identifying a geometric invariant of protein structures, termed geometrical core, is developed. The method used the matrix of C(alpha)-C(alpha) distances and does not require the usual superposition of structures. The result of applying the algorithm to 53 immunoglobulin structures led to the identification of two geometrical core sets of C(alpha) atoms positions for the V(L) and V(H) domains. Based on these geometric invariants a preferred coordinate system for the immunoglobulin family is constructed which serves as a basis for structural prediction. The X-ray atom coordinates for all available immunoglobulin structures are transformed to the preferred coordinate system. An affine symmetry between the V(L) and V(H) domains is defined and computed for each of the 53 immunoglobulin structures.

A very limited number of keywords (main patterns) describes all sequences of the human variable heavy (VH) and kappa (Vkappa) domains.

Sequences of the variable heavy (VH) and kappa (Vkappa) domains of Ig structures were divided into 21 fragments that correspond to strands, loops, or parts of these structural units of the variable domains. Amino acid sequences of fragments (termed "words") were collected from the 1,172 human heavy and 668 human kappa chains available in the Kabat database. Statistical analysis of words of 17 fragments was performed (fragments that comprise the complementary determining regions' fragments will not be discussed in this paper). The number of different words (those with different residues in at least one position) ranged, for various fragments, from 11 to 75 in the kappa chains, and from 23 to 189 in the heavy chains. The main result of this study is that very few keywords, or main patterns of words, were necessary to describe over 90% of the sequences (no more than two keywords per fragment in the kappa and no more than five per fragment in the heavy chains). No identical keywords were found for different fragments of the variable domains. Keywords of aligned fragments of the VH and Vkappa domains were different in all but two instances. Thus, knowing the keywords, one can determine whether any given small part of a sequence belongs to a heavy or kappa chain and predict its precise localization in the sequence. In addition, by using all of the keywords obtained through analysis of the Kabat database, it was possible to describe completely the sequences of the human VH and Vkappa germ-line segments.

Adult seronegative arthritis with antinuclear antibodies: a distinct group of patients with a different immunogenetic pattern from seropositive rheumatoid arthritis and a good outcome.

To determine whether patients with rheumatoid factor (RF)-negative, antinuclear antibody (ANA)-positive oligo/polyarthritis are clinically and immunogenetically distinct from RF-positive rheumatoid arthritis (RA) and whether this subset of patients is the adult counterpart of early-onset pauciarticular juvenile chronic arthritis (EOPA JCA), we retrospectively studied 20 adult patients with RF-negative, ANA-positive arthritis. After a median duration of 3.25 years, only half of our patients had active synovitis. Seventy-five per cent were completely or reasonably self-sufficient according to the HAQ index. In contrast to the results in a group of 30 RF-positive RA patients, the percentage of patients having at least one of the susceptibility alleles (HLA DR1 or HLA-DR4) was not significantly higher in patients with RF-negative, ANA-positive arthritis than in controls. Furthermore, none of our RF-negative, ANA-positive patients had two susceptibility alleles, whereas 16.5% of RF-positive RA patients had both DR1 and DR4 or DR4 homozygosity. In conclusion, our results show that patients with RF-negative, ANA-positive oligo/polyarthritis are immunogenetically distinct from RF-positive RA and tend to have a better articular prognosis. The absence of typical ocular features and of the characteristic HLA-DR markers suggests that these patients cannot be considered as the adult counterpart of EOPA JCA.

The invariant system of coordinates of antibody molecules: prediction of the "standard" C alpha framework of VL and VH domains.

A new approach of comparing protein structures that does not involve the procedure of superposition is suggested. An invariant system of coordinates for immunoglobulin molecules that is based on the geometrical symmetry inherent to the variable domain light-chain (VL)-heavy-chain (VH) complex is described. The coordinates of the Calpha atoms in 22 immunoglobulin structures are calculated in the invariant system of coordinates. We found that 76 identical positions in this Calpha framework are symmetrical about the twofold axis. Comparison of the identical positions in these molecules allows us to select 96 positions in the light chains and 87 positions in the heavy chains whose Calpha atom coordinates are approximately the same. To check whether the average coordinates of Calpha atoms in these positions complies with the stereochemical requirements, we calculated Calpha-Calpha distances. Seventy-three positions of the light chains and 72 positions of the heavy chains satisfy the Calpha-Calpha distance criterion. The Calpha atoms in these positions are used for constructing the "standard" Calpha framework of VL and VH complexes. The average coordinates of Calpha atoms are presented.

The IgG binding site of human FcgammaRIIIB receptor involves CC' and FG loops of the membrane-proximal domain.

Fc gamma receptors for the Fc part of IgG are the mediators for antibody effector functions. FcgammaRIII and FcgammaRII are low affinity receptors that, through the interaction with immune complexes, initiate a variety of immunological responses, such as phagocytosis, antibody-dependent cellular cytotoxicity, and release of inflammatory mediators. We set out to define the IgG binding site on human FcgammaRIII. We assumed that potential beta-turns in Ig-like domains are the most probable determinants for ligand binding, and chimeric FcgammaRIIIB/FcepsilonRI receptors as well as single residue mutants were constructed in these regions of FcgammaRIIIB. Substitution of four amino acids in the membrane-proximal domain (Gln126, Arg156, Lys162, Val164) resulted in decreased binding of human IgG1. Lys162 and Val164 were found also to be crucial for the interaction with the IgG-binding inhibitory monoclonal antibody 3G8. In a putative three-dimensional model constructed in this study, these residues map on the CC loop (Gln126), on F beta-sheet (Arg156), and on the FG loop (Lys162, Val164). Our data are consistent with the study about human FcgammaRII (Hulett, M. D., Witort, E., Brinkworth, R. I., McKenzie, I. F. C., and Hogarth, P. M. (1994) J. Biol. Chem. 269, 15287-15293), suggesting that common structural determinants, i.e. FG loop or the GFC surface of the membrane-proximal domain, can be involved in interactions with IgG by both low affinity receptor classes FcgammaRII and FcgammaRIII.

Analysis of the relation between the sequence and secondary and three-dimensional structures of immunoglobulin molecules.

Methods of structural and statistical analysis of the relation between the sequence and secondary and three-dimensional structures are developed. About 5000 secondary structures of immunoglobulin molecules from the Kabat data base were predicted. Two statistical analyses of amino acids reveal 47 universal positions in strands and loops. Eight universally conservative positions out of the 47 are singled out because they contain the same amino acid in > 90% of all chains. The remaining 39 positions, which we term universally alternative positions, were divided into five groups: hydrophobic, charged and polar, aromatic, hydrophilic, and Gly-Ala, corresponding to the residues that occupied them in almost all chains. The analysis of residue-residue contacts shows that the 47 universal positions can be distinguished by the number and types of contacts. The calculations of contact maps in the 29 antibody structures revealed that residues in 24 of these 47 positions have contacts only with residues of antiparallel beta-strands in the same beta-sheet and residues in the remaining 23 positions always have far-away contacts with residues from other beta-sheets as well. In addition, residues in 6 of the 47 universal positions are also involved in interactions with residues of the other variable or constant domains.