Populations of hydrophobic amino acids within protein globular domains: Identification of conserved “topohydrophobic” positions

The 3D structural comparison of families of divergent homologous domains revealed two main populations of hydrophobic amino acids, one with a low and the other with a significantly higher mean solvent accessibility, allowing two regions of the core of protein globular domains to be distinguished. The side chains of hydrophobic amino acids in topologically conserved positions (positions in the structural alignment where only hydrophobic amino acids are found), which we call topohydrophobic positions, are considerably less dispersed than those of the other amino acids (hydrophobic or not). Mean distances between gravity centers of amino acids in topohydrophobic positions are significantly shorter than those for non-topohydrophobic positions and show that the corresponding amino acids are almost all in direct contact in the inner core of globular domains. This study also showed that the small number of topohydrophobic positions is a characteristic of the structural differences between proteins of a family. This criterion is independent of the sequence identity between the sequences and of the root-mean-square distance between their corresponding structures. Using sensitive sequence alignment processes it will be possible, for many protein families, to identify topohydrophobic positions from sequences only. Proteins 33:329–342, 1998. © 1998 Wiley-Liss, Inc.     

Application of genetic semihomology algorithm to theoretical studies on various protein families.

Several protein families of different nature were studied for genetic relationship, correct alignment at non-homologous fragments, optimal sequence consensus construction, and confirmation of their actual relevance. A comparison of the genetic semihomology approach with statistical approaches indicates a high accuracy and cognition significance of the former. This is particularly pronounced in the study of related proteins that show a low degree of homology. The sequence multiple alignments were verified and corrected with respect to the questionable, non-homologous fragments. The verified alignments were the basis for consensus sequence formation. The frequency of six-codon amino acids occurrence versus position variability was studied and their possible role in amino acid mutational exchange at variable positions is discussed.     

Amino acid sequence of the N-terminal sixty-nine residues of heavy chain derived from a homogeneous rabbit antibody

The sequence of the N-terminal 69 residues of heavy chain from a homogeneous rabbit antibody to type III pneumococcal polysaccharide was determined. The sequence is similar to that found in heavy chains of normal pooled rabbit immunoglobulins of the same allotype Aa1. Two regions of the homogeneous heavy chain (residues 35-46 and 62-69) are very similar to corresponding regions of heavy chains from rabbit Aa2 immunoglobulin, as well as from mouse, guinea-pig and human immunoglobulins. In contrast, residues 47-62 appear to be variable. Comparison in this section with another homogeneous anti-pneumococcal antibody (Strosberg et al., 1972) of related specificity and of the same allotype indicates sequence variation in at least three positions. An antibody to group C streptococcal carbohydrate of allotype Aa2 (Fleischman, 1971) differs by five amino acids in the same region of the heavy chain. Sequence variability between these three antibodies does not occur in homologous positions within this variable section. Allotype-related sequences could not be identified in section 34-65.     

Amino acid sequences of antibody light chain variable regions of pedigreed rabbits: kappa light chain K49-501 (allotype b4 anti-streptococcal group A-variant polysaccharide antibody).

The amino acid sequence of positions 1--150 of a light chain, isolated from another monoclonal rabbit anti-streptococcal group A-variant polysaccharide antibody, was determined. The analysis was performed with 2 mumol of polypeptide chain, using a grossly modified Beckman 890B sequenator. This sequence stretch accounts for the whole variable region and a considerable part of the constant region at a total length of 218 amino acids. This allotype b4 light chain was isolated from a non-precipitating, end-group-specific antibody with a KD = 1.3 X 10(-5)M. This brings the present number of totally known rabbit VL sequences of antigen elicited antibodies to 21. A comparison of these 21 sequences reveals a building plan of ribbit VL homologous to that of human and murine VL regions. The observed variability does follow a pattern of linked amino acid substitutions, indicating that this information must be contained in the germ-line of the rabbit in the form of multiple VL region genes. This conclusion, however, does not rule out the occasional variant being due to somatic rearrangement. Finally, this comparison reveals that the joining peptide between positions 96--110 is also a separate entity in rabbit VL region sequences.     

Functional consequences of insertions and deletions in the complementarity-determining regions of human antibodies

Insertions and deletions of nucleotides in the genes encoding the variable domains of antibodies are natural components of the hypermutation process, which may expand the available repertoire of hypervariable loop lengths and conformations. Although insertion of amino acids has also been utilized in antibody engineering, little is known about the functional consequences of such modifications. To investigate this further, we have introduced single-codon insertions and deletions as well as more complex modifications in the complementarity-determining regions of human antibody fragments with different specificities. Our results demonstrate that single amino acid insertions and deletions are generally well tolerated and permit production of stably folded proteins, often with retained antigen recognition, despite the fact that the thus modified loops carry amino acids that are disallowed at key residue positions in canonical loops of the corresponding length or are of a length not associated with a known canonical structure. We have thus shown that single-codon insertions and deletions can efficiently be utilized to expand structure and sequence space of the antigen-binding site beyond what is encoded by the germline gene repertoire.     

Relationships between sequence and structure for the four-alpha-helix bundle tertiary motif in proteins.

The sequences of four-alpha-helical bundle proteins are characterized by a pattern of hydrophilic and hydrophobic amino acids which is repeated every seven residues. At each position of the heptad repeat there are specific constraints on the amino acid properties which result from the topology of the tertiary motif. These constraints give rise to patterns of amino acid distribution which are distinct from those of other proteins. The distributions in each of the heptad positions have been determined by a statistical analysis of structural and sequence data derived from seven families of aligned protein sequences. The constitution of each position is dominated by a very small number of different amino acids, with the core positions consisting overwhelmingly of Leu and Ala. The positional preferences of the individual amino acids can be generally interpreted in terms of residue properties and topological constraints. The potential for four-alpha-helix bundle folding is reflected primarily in the pattern of residue occurrence in the heptad and not in the overall amino acid composition of the protein. Possible applications of this analysis in structure predictions, sequence alignments and in the rational design and engineering of four-alpha-helical bundle proteins are discussed.     

Sequence analysis of an artificial family of RNA-binding peptides

Diverse peptide sequences recognizing the lambda boxB RNA hairpin were previously isolated from a library encoding the 22-residue lambda N peptide with random amino acids at positions 13-22 using mRNA display. We have statistically analyzed amino acid distributions in 65 unique sequences from rounds 11 and 12 of this selection and evaluated the resulting structural and functional predictions by alanine-scanning mutagenesis and circular dichroism spectrometry. This artificial sequence family has a consensus structure that continues the bent alpha helix of lambda N up to position 17 when bound to lambda boxB. A charge pair (E(14)R(15)) and hydrophobic patch (A(21)L(22) or V(21)L(22)) have important functional roles in this context. Notably, amino acid covariance reveals six specific pairs of random region positions with >95% significant linkage and strong overall helical (i+1, i+3, and i+4) couplings. The covariance analysis suggests that (1) the sequence context of every residue in each insert has been optimized, (2) selected sequences are local optima on a rugged fitness landscape, and (3) it is possible to detect more subtle structural features with artificial protein sequence families than natural homologs. Our results provide a framework for investigating the structures of in vitro selected proteins by functional minimization, reselection, and covariance analysis.     

Molecular cloning and polymorphism analysis of the prion protein gene in Tan sheep of Ningxia, China

The resistance or susceptibility of sheep to scrapie is associated with polymorphisms of the prion protein gene (PRNP), particularly, single nucleotide polymorphisms (SNPs) in amino acid positions 136, 154 and 171. The prion protein (PrP) gene sequence and the deduced amino acid alignment of prion protein in Tan sheep, a local Chinese sheep breed traditionally raised in Ningxia, northwestern China, were determined and variability of the PrP amino acids sequence was analyzed in this study. The PrP nucleic acids and amino acids sequences of 112 Tan sheep were highly homogenous, although polymorphism of the PrP gene was detected at several sites, particularly codons 106, 154, and 171. The analysis of both sequences revealed that the most predominant allele at codons 136, 154 and 171 in Tan sheep was ARQ, which was known to be associated with high susceptibility to scrapie in sheep. The result suggests that Tan sheep is potentially susceptible to scrapie. Our findings provide valuable information for future breeding projects to scrapie resistance in Tan sheep.     

Simian immunodeficiency virus mutants resistant to serum neutralization arise during persistent infection of rhesus monkeys.

We previously described the pattern of sequence variation in gp120 following persistent infection of rhesus monkeys with the pathogenic simian immunodeficiency virus SIVmac239 molecular clone (D.P.W. Burns and R.C. Desrosiers, J. Virol. 65:1843, 1991). Sequence changes were confined largely to five variable regions (V1 to V5), four of which correspond to human immunodeficiency virus type 1 (HIV-1) gp120 variable regions. Remarkably, 182 of 186 nucleotide substitutions that were documented in these variable regions resulted in amino acid changes. This is an extremely nonrandom pattern, which suggests selective pressure driving amino acid changes in discrete variable domains. In the present study, we investigated whether neutralizing-antibody responses are one selective force responsible at least in part for the observed pattern of sequence variation. Variant env sequences called 1-12 and 8-22 obtained 69 and 93 weeks after infection of a rhesus monkey with cloned SIVmac239 were recombined into the parental SIVmac239 genome, and variant viruses were generated by transfection of cultured cells with cloned DNA. The 1-12 and 8-22 recombinants differ from the parental SIVmac239 at 18 amino acid positions in gp120 and at 5 and 10 amino acid positions, respectively, in gp41. Sequential sera from the monkey infected with cloned SIVmac239 from which the 1-12 and 8-22 variants were isolated showed much higher neutralizing antibody titers to cloned SIVmac239 than to the cloned 1-12 and 8-22 variants. For example, at 55 weeks postinfection the neutralizing antibody titer against SIVmac239 was 640 while those to the variant viruses were 40 and less than 20. Two other rhesus monkeys infected with cloned SIVmac239 showed a similar pattern. Rhesus monkeys were also experimentally infected with the cloned variants so that the type-specific nature of the neutralizing antibody responses could be verified. Indeed, each of these monkeys showed neutralizing-antibody responses of much higher titer to the homologous variant used for infection. These experiments unambiguously demonstrate that SIV mutants resistant to serum neutralization arise during the course of persistent infection of rhesus monkeys.     

Conservation of intrinsic disorder in protein domains and families: I. A database of conserved predicted disordered regions

Many protein regions have been shown to be intrinsically disordered, lacking unique structure under physiological conditions. These intrinsically disordered regions are not only very common in proteomes, but also crucial to the function of many proteins, especially those involved in signaling, recognition, and regulation. The goal of this work was to identify the prevalence, characteristics, and functions of conserved disordered regions within protein domains and families. A database was created to store the amino acid sequences of nearly one million proteins and their domain matches from the InterPro database, a resource integrating eight different protein family and domain databases. Disorder prediction was performed on these protein sequences. Regions of sequence corresponding to domains were aligned using a multiple sequence alignment tool. From this initial information, regions of conserved predicted disorder were found within the domains. The methodology for this search consisted of finding regions of consecutive positions in the multiple sequence alignments in which a 90% or more of the sequences were predicted to be disordered. This procedure was constrained to find such regions of conserved disorder prediction that were at least 20 amino acids in length. The results of this work included 3,653 regions of conserved disorder prediction, found within 2,898 distinct InterPro entries. Most regions of conserved predicted disorder detected were short, with less than 10% of those found exceeding 30 residues in length.     

Evolutionarily conserved allosteric network in the Cys loop family of ligand-gated ion channels revealed by statistical covariance analyses

The Cys loop family of ligand-gated ion channels mediate fast synaptic transmission for communication between neurons. They are allosteric proteins, in which binding of a neurotransmitter to its binding site in the extracellular amino-terminal domain triggers structural changes in distant transmembrane domains to open a channel for ion flow. Although the locations of binding site and channel gating machinery are well defined, the structural basis of the activation pathway coupling binding and channel opening remains to be determined. In this paper, by analyzing amino acid covariance in a multiple sequence alignment, we have identified an energetically interconnected network in the Cys loop family of ligand-gated ion channels. Statistical coupling and correlated mutational analyses along with clustering revealed a highly coupled cluster. Mapping the positions in the cluster onto a three-dimensional structural model demonstrated that these highly coupled positions form an interconnected network linking experimentally identified binding domains through the coupling region to the gating machinery. In addition, these highly coupled positions are also condensed in the transmembrane domains, which are a recent focus for the sites of action of many allosteric modulators. Thus, our results revealed a genetically interconnected network that potentially plays an important role in the allosteric activation and modulation of the Cys loop family of ligand-gated ion channels.     

Isolation,expression, and the primary structure of HLA-Cw1 and HLA-Cw2 genes: Evolutionary aspects

The HLA-Cw1 and -Cw2 genes were identified in a genomic library and their products characterized by biochemical methods. The HLA-Cw and -Cw2 genes, upon transfection in mouse L cells, give rise to class I antigen heavy chains that associate with neither mouse nor human beta-2 microglobulin. They are indistinguishable in isoelectric point from polypeptides identified as HLA-Cw1 and -Cw2 in human cells. The nucleotide sequence of HLA-Cw1 and -Cw2 and their comparison with HLA-Cw3, the only other known HLA-C sequence, reveal a characteristic pattern of locus-specific amino acids. A comparison of 13 different human class I primary structures leads us to speculate that the most variable region in HLA class I antigens, positions 61–83, could assume an alpha helical structure of critical importance for class I antigen function. The locus specificity and the higher degree of intralocus conservation in the COON-terminal region, especially in the transmembrane and cytoplasmic domains, must reflect evolutionary ancestry rather than positive selection. In view of the pattern and types of substitutions observed for HLA-C locus products, their function as immune response gene products is questioned.     

Interplay of selective pressure and stochastic events directs evolution of the MEL172 satellite DNA library in root-knot nematodes

According to the library model, related species can have in common satellite DNA (satDNA) families amplified in differing abundances, but reasons for persistence of particular sequences in the library during long periods of time are poorly understood. In this paper, we characterize 3 related satDNAs coexisting in the form of a library in mitotic parthenogenetic root-knot nematodes of the genus Meloidogyne. Due to sequence similarity and conserved monomer length of 172 bp, this group of satDNAs is named MEL172. Analysis of sequence variability patterns among monomers of the 3 MEL172 satellites revealed 2 low-variable (LV) domains highly reluctant to sequence changes, 2 moderately variable (MV) domains characterized by limited number of mutations, and 1 highly variable (HV) domain. The latter domain is prone to rapid spread and homogenization of changes. Comparison of the 3 MEL172 consensus sequences shows that the LV domains have 6% changed nucleotide positions, the MV domains have 48%, whereas 78% divergence is concentrated in the HV domain. Conserved distribution of intersatellite variability might indicate a complex pattern of interactions in heterochromatin, which limits the range and phasing of allowed changes, implying a possible selection imposed on monomer sequences. The lack of fixed species-diagnostic mutations in each of the examined MEL172 satellites suggests that they existed in unaltered form in a common ancestor of extant species. Consequently, the evolution of these satellites seems to be driven by interplay of selective constraints and stochastic events. We propose that new satellites were derived from an ancestral progenitor sequence by nonrandom accumulation of mutations due to selective pressure on particular sequence segments. In the library of particular taxa, established satellites might be subject to differential amplification at chance due to stochastic mechanisms of concerted evolution.     

Sequence plasticity in the antigen-binding site of a therapeutic anti-HER2 antibody

We have examined the plasticity of the antigen-combining site of a high-affinity antibody. In phage-displayed Fab libraries, selected CDR positions and one FR position of the humanized anti-Her2 antibody hu4D5 were substituted with all 20 amino acids. Antigen-binding selections were used to enrich for high-affinity variants, and a large number of sequences were obtained prior to convergence of the selected pool to a small set of clones. As expected, sequence variability of the antigen-binding site is overall diminished compared to known IgG sequences; however, certain positions retain much higher variability than others. The sequence variability map of the hu4D5 binding site is compared with a map derived from previous alanine-scanning of the antibody. Affinities of soluble Fab fragments for antigen confirm that multiple variants were selected with high affinity for antigen, including one variant with a single point mutation that was about threefold improved in affinity compared to the parental hu4D5. Interestingly, this mutation is one of the most radical in terms of changing side-chain chemistry (Trp for Asp) and occurs at the most plastic site as calculated by the Wu-Kabat variability coefficient. Thus variability mapping yields information about the antibody-antigen interaction that is useful and complementary to that obtained by alanine scanning mutagenesis.     

The primary structure of the VLA-2/collagen receptor alpha 2 subunit (platelet GPIa): homology to other integrins and the presence of a possible collagen-binding domain

VLA-2 (also called gpIa/IIa on platelets) is a collagen receptor with a unique alpha subunit and a beta subunit common to other adhesion receptors in the VLA/integrin family. Multiple cDNA clones for the human VLA-2 alpha 2 subunit have been selected from a lambda gtll library by specific antibody screening. The 5,374-bp nucleotide sequence encoded for 1,181 amino acids, including a signal peptide of 29 amino acids followed by a long extracellular domain (1,103 amino acids), a transmembrane domain, and a short cytoplasmic segment (22 amino acids). Direct sequencing of purified alpha 2 protein confirmed the identity of the 15 NH2-terminal amino acids. Overall, the alpha 2 amino acid sequence was 18-25% similar to the sequences known for other integrin alpha subunits. In particular, the alpha 2 sequence matched other integrin alpha chains in (a) the positions of 17 of its 20 cysteine residues; (b) the presence of three metal-binding domains of the general structure DXDXDGXXD; and (c) the transmembrane domain sequence. In addition, the alpha 2 sequence has a 191-amino acid insert (called the I-domain), previously found only in leukocyte integrins of the beta 2 integrin family. The alpha 2 I-domain was 23-41% similar to domains in cartilage matrix protein and von Willebrand factor, which are perhaps associated with collagen binding. The NH2-terminal sequence reported here for alpha 2 does not match the previously reported alpha 2 NH2-terminal sequence (Takada, Y., J. L. Strominger, and M. E. Hemler. 1987. Proc. Natl. Acad. Sci. USA. 84:3239-3243). Resolution of this discrepancy suggests that there may be another VLA heterodimer that resembles VLA-2 in size but has a different amino acid sequence.