Prediction of Secondary Structure,Spatial Organization and Distribution of Antigenic Determinants for Hepatitis A Virus Proteins

Abstract

On the basis of the secondary structure calculations from the known amino acid sequence we came to the conclusion that hepatitis A virus capsid proteins have the typical antiparallel beta-sheet bilayer structure.

The predicted secondary structure of the HAV proteins can be well aligned with those of the poliovirus (type 1 Mahoney) and human rhinovirus (type 14). It enabled us to use the X-ray structure of the PV-1M and HRV-14 proteins as a template and then, firstly, to localize the positions of alpha and beta regions in the architecture of the HAV protein molecules and, secondly, to discover the amino acid homologies of the secondary structure regions aligned. The obtained model of the three-dimensional structure for HAV proteins helped us to indicate the exposed regions of the polypeptide chains and to pinpoint the potential neutralizing antigenic sites.     

Detection of a genome-linked protein (VPg) of hepatitis A virus and its comparison with other picornaviral VPgs.

The nucleotide sequence corresponding to the P3 region of the hepatitis A virus (HAV) polyprotein genome was determined from cloned cDNA and translated into an amino acid sequence. Comparison of the amino acid sequences of the genome-linked proteins (VPgs) of other picornaviruses with the predicted amino acid sequence of HAV was used to locate the primary structure of a putative VPg within the genome of HAV. The sequence of HAV VPg, like those of other picornaviral VPg molecules, contains a tyrosine residue as a potential binding site for HAV RNA in position 3 from its N terminus. The potential cleavage sites to generate VPg from a putative HAV polyprotein are between glutamic acid and glycine at the N terminus and glutamic acid and serine or glutamine and serine at the C terminus. A synthetic peptide corresponding to 10 amino acids of the predicted C terminus of HAV VPg induced anti-peptide antibodies in rabbits when it was conjugated to thyroglobulin as a carrier. These antibodies were specific for the peptide and precipitated VPg, linked to HAV RNA, from purified HAV and from lysates of HAV-infected cells. The precipitation reaction was blocked by the synthetic peptide (free in solution or coupled to carrier proteins) and prevented by pretreatment of VPg RNA with protease. Thus, our predicted amino acid sequence is colinear with the nucleotide sequence of the VPg gene in the HAV genome. From our results we concluded that HAV has the typical organization of picornavirus genes in this part of its genome. Similarity among hydrophobicity patterns of amino acid sequences of different picornaviral VPgs was revealed in hydropathy plots. Thus, the VPg of HAV appears to be closely related to VPg1 and VPg2 of foot-and-mouth disease virus. In contrast, HAV VPg has a unique isoelectric point (pI = 7.15) among the picornavirus VPgs.     

Occurrence,conformational features and amino acid propensities for the pi-helix

The most abundant helix type in proteins is the alpha-helix, accounting for about 31% of amino acid secondary structure states, while the 3(10)-helix accounts for about 4%. The pi-helix appears to be extremely rare and is considered to be unstable. Existing secondary structure definition methods find very few within the Protein Data Bank. Using an improved pi-helix definition algorithm to search a non-redundant subset of high-resolution and well-refined protein structures, we found that almost every tenth protein contained a pi-helix. This enabled us to show for the first time that the pi-helix has structural parameters that are different from the hypothesized model values. It also has distinctive amino acid preferences and it is conserved within functionally related proteins. Features that may contribute to the stability of the pi-helical structure have also been identified. In addition to hydrogen bonds, several other factors contribute to the stability of pi-helices. The pi-helix may have some functional advantages over other helical structures. Thus, we describe cases where the side chains of functionally important residues at every fourth position within a pi-helix could be aligned and brought close together in a way that would not be allowed by any other helix type.     

Characterization of a simian hepatitis A virus (HAV): antigenic and genetic comparison with human HAV.

PA21, a strain of hepatitis A virus (HAV) recovered from a naturally infected captive owl monkey, is indistinguishable from human HAV in polyclonal radioimmunoassays and cross-neutralization studies. However, cDNA-RNA hybridization has suggested a significant difference at the genomic level between PA21 and a reference human virus, HM175. Further characterization of this unique HAV was undertaken in an effort to determine the extent of genetic divergence from human HAV and its relation to the conserved antigenic structure of the virus. The close similarity between PA21 and HM175 antigens was confirmed with an extended panel of 18 neutralizing murine monoclonal antibodies: a reproducible difference in binding to the two viruses was detected with only one antibody (B5-B3). The nucleotide sequence of the P1 region of the PA21 genome had only 83.2% identity with HM175 virus, a difference approximately twice as great as that found between any two human strains. Most nucleotide changes were in third base positions, and the amino acid sequences of the capsid proteins were largely conserved. Amino acid replacements were clustered in the carboxy terminus of VP1 and the amino-terminal regions of VP2 and VP1. These data indicate that PA21 virus represents a unique genotype of HAV and suggest the existence of an ecologically isolated niche for HAV among feral owl monkeys.     

The complete primary structure of ribosomal proteins L1, L14, L15, L23, L24 and L29 from Bacillus stearothermophilus

The amino acid sequences of ribosomal proteins L1, L14, L15, L23, L24 and L29 from Bacillus stearothermophilus have been completely determined. This has been achieved by sequence analyses of peptides derived from enzymatic digestions of the proteins with trypsin, chymotrypsin, pepsin, Staphylococcus aureus protease, and Armillaria mellea protease as well as by chemical cleavage with hydroxylamine and cyanogen bromide. Based on the primary structures of the six proteins, their secondary structures were predicted using four different computer prediction programs. A comparison of the amino acid sequences of the studied proteins from B. stearothermophilus with the homologous proteins from Escherichia coli revealed that in four proteins (L1, L15, L24 and L29) between 40-50% of the residue in the sequences are identical, whereas this value is significantly higher (69%) for L14 and lower (28%) for L23. The distribution of those amino acid residues which are identical in the corresponding proteins from the two bacteria is not random along the protein chain: some regions are highly conserved whereas others are not. This finding indicates that the regions which are conserved during evolution are important for the spatial structure and/or function of the protein.     

Fibroblast growth factor receptors contain a conserved HAV region common to cadherins and influenza strain A hemagglutinins: a role in protein-protein interactions?

The first extracellular domain of the cadherins has been shown to exhibit extensive sequence homology with the amino termini of the HA1 chains of influenza strain A hemagglutinins. These regions of homology are known to be functionally important in both the cadherins and the hemagglutinins. The homologous regions harbor the tripeptide HAV, which has been identified as being the cadherin cell adhesion recognition sequence. Here we report that members of the rapidly expanding family of fibroblast growth factor receptors also possess HAV-containing regions. These regions are homologous to the HAV-containing regions present within both the hemagglutinins and the cadherins and appear to be involved in regulating the function of the fibroblast growth factor receptors. We speculate that the HAV motif may represent an evolutionarily conserved amino acid sequence that will prove to be functionally important in a wide variety of proteins.     

Complete nucleotide sequence of wild-type hepatitis A virus: comparison with different strains of hepatitis A virus and other picornaviruses.

The complete nucleotide sequence of wild-type hepatitis A virus (HAV) HM-175 was determined. The sequence was compared with that of a cell culture-adapted HAV strain (R. Najarian, D. Caput, W. Gee, S.J. Potter, A. Renard, J. Merryweather, G.V. Nest, and D. Dina, Proc. Natl. Acad. Sci. USA 82:2627-2631, 1985). Both strains have a genome length of 7,478 nucleotides followed by a poly(A) tail, and both encode a polyprotein of 2,227 amino acids. Sequence comparison showed 624 nucleotide differences (91.7% identity) but only 34 amino acid differences (98.5% identity). All of the dipeptide cleavage sites mapped in this study were conserved between the two strains. The sequences of these two HAV strains were compared with the partial sequences of three other HAV strains. Most amino acid differences were located in the capsid region, especially in VP1. Whereas changes in amino acids were localized to certain portions of the genome, nucleotide differences occurred randomly throughout the genome. The most extensive nucleotide homology between the strains was in the 5' noncoding region (96% identity for cell culture-adapted strains versus wild type; greater than 99% identity among cell culture-adapted strains). HAV proteins are less homologous with those of any other picornavirus than the latter proteins are when compared with each other. When the sequences of wild-type and cell culture-adapted HAV strains are compared, the nucleotide differences in the 5' noncoding region and the amino acid differences in the capsid region suggest areas that may contain markers for cell culture adaptation and for attenuation.     

Predicting reliable regions in protein alignments from sequence profiles

For applications such as comparative modelling one major issue is the reliability of sequence alignments. Reliable regions in alignments can be predicted using sub-optimal alignments of the same pair of sequences. Here we show that reliable regions in alignments can also be predicted from multiple sequence profile information alone.Alignments were created for a set of remotely related pairs of proteins using five different test methods. Structural alignments were used to assess the quality of the alignments and the aligned positions were scored using information from the observed frequencies of amino acid residues in sequence profiles pre-generated for each template structure. High-scoring regions of these profile-derived alignment scores were a good predictor of reliably aligned regions.These profile-derived alignment scores are easy to obtain and are applicable to any alignment method. They can be used to detect those regions of alignments that are reliably aligned and to help predict the quality of an alignment. For those residues within secondary structure elements, the regions predicted as reliably aligned agreed with the structural alignments for between 92% and 97.4% of the residues. In loop regions just under 92% of the residues predicted to be reliable agreed with the structural alignments. The percentage of residues predicted as reliable ranged from 32.1% for helix residues to 52.8% for strand residues.This information could also be used to help predict conserved binding sites from sequence alignments. Residues in the template that were identified as binding sites, that aligned to an identical amino acid residue and where the sequence alignment agreed with the structural alignment were in highly conserved, high scoring regions over 80% of the time. This suggests that many binding sites that are present in both target and template sequences are in sequence-conserved regions and that there is the possibility of translating reliability to binding site prediction.     

Structure of the hepatitis A virion: peptide mapping of the capsid region.

Milligram amounts of highly purified hepatitis A virus (HAV) were obtained from persistently infected cell cultures. The HAV polypeptides were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and transferred to nitrocellulose for detection by an enzyme-linked immunotransfer blot procedure. The HAV nucleotide-derived amino acid sequence was subjected to computer analysis to identify potential immunogenic regions within the HAV capsid polypeptides. Synthetic peptides corresponding to selected regions of each of the larger putative capsid polypeptides were coupled to keyhole limpet hemocyanin and used to immunize rabbits. Four of six anti-HAV peptide sera were strongly reactive. Antipeptide serum generated against amino acids (a.a.) 75 through 82 reacted with the 27,000-molecular-weight (MW) polypeptide; serum against a.a. 279 through 285 reacted with the 29,000-MW HAV polypeptide; and sera against a.a. 591 through 602 and 606 through 618 reacted with the 33,000-MW HAV polypeptide. These reactions enabled the identification of the gene order of the larger HAV P1 region gene products. Our data indicate the following molecular weights: HAV VP2 or 1B, 27,000; HAV VP3 or 1C, 29,000; and HAV VP1 or 1D, 33,000.     

Environment-specific amino acid substitution tables: tertiary templates and prediction of protein folds.

The local environment of an amino acid in a folded protein determines the acceptability of mutations at that position. In order to characterize and quantify these structural constraints, we have made a comparative analysis of families of homologous proteins. Residues in each structure are classified according to amino acid type, secondary structure, accessibility of the side chain, and existence of hydrogen bonds from the side chains. Analysis of the pattern of observed substitutions as a function of local environment shows that there are distinct patterns, especially for buried polar residues. The substitution data tables are available on diskette with Protein Science. Given the fold of a protein, one is able to predict sequences compatible with the fold (profiles or templates) and potentially to discriminate between a correctly folded and misfolded protein. Conversely, analysis of residue variation across a family of aligned sequences in terms of substitution profiles can allow prediction of secondary structure or tertiary environment.     

Properties of polyproline II, a secondary structure element implicated in protein-protein interactions

The polyproline II (PPII) conformation of protein backbone is an important secondary structure type. It is unusual in that, due to steric constraints, its main-chain hydrogen-bond donors and acceptors cannot easily be satisfied. It is unable to make local hydrogen bonds, in a manner similar to that of alpha-helices, and it cannot easily satisfy the hydrogen-bonding potential of neighboring residues in polyproline conformation in a manner analogous to beta-strands. Here we describe an analysis of polyproline conformations using the HOMSTRAD database of structurally aligned proteins. This allows us not only to determine amino acid propensities from a much larger database than previously but also to investigate conservation of amino acids in polyproline conformations, and the conservation of the conformation itself. Although proline is common in polyproline helices, helices without proline represent 46% of the total. No other amino acid appears to be greatly preferred; glycine and aromatic amino acids have low propensities for PPII. Accordingly, the hydrogen-bonding potential of PPII main-chain is mainly satisfied by water molecules and by other parts of the main-chain. Side-chain to main-chain interactions are mostly nonlocal. Interestingly, the increased number of nonsatisfied H-bond donors and acceptors (as compared with alpha-helices and beta-strands) makes PPII conformers well suited to take part in protein-protein interactions.     

Computational methods for predicting structure of membrane proteins using amino acid sequences

The structure of membrane proteins specifies their functional properties, which are important for medicine and pharmacology and, therefore, is of significant interest. The repetition of transmembrane regions that consist of hydrophobic amino acids is a characteristic and organic feature of polytopic membrane proteins. The ordered repetition (periodicity) can be detected by the Fourier method applied to a digital image of the symbolic amino acid sequence of a protein. In the present work, this investigation was carried out for 24 transmembrane proteins (successfully for 14 of them). If the repetition of transmembrane regions is aperiodic, it can be revealed by another method, that is, the method of the reiterated (four to five times) averaging of the protein hydrophobicity function in a window within the limits of 9–11 amino acids that moves along the sequence. This novel method was applied to the 24 transmembrane proteins (successfully for 19 of them) and demonstrated higher suitability than the Fourier method for predicting the secondary structure of these proteins and the corresponding functional properties.     

Reinvestigating the codon and amino acid usage of S. cerevisiae genome: a new insight from protein secondary structure analysis

Biased usage of synonymous codons has been elucidated under the perspective of cellular tRNA abundance for quite a long time now. Taking advantage of publicly available gene expression data for Saccharomyces cerevisiae, a systematic analysis of the codon and amino acid usages in two different coding regions corresponding to the regular (helix and strand) as well as the irregular (coil) protein secondary structures, have been performed. Our analyses suggest that apart from tRNA abundance, mRNA folding stability is another major evolutionary force in shaping the codon and amino acid usage differences between the highly and lowly expressed genes in S. cerevisiae genome and surprisingly it depends on the coding regions corresponding to the secondary structures of the encoded proteins. This is obviously a new paradigm in understanding the codon usage in S. cerevisiae. Differential amino acid usage between highly and lowly expressed genes in the regions coding for the irregular protein secondary structure in S. cerevisiae is expounded by the stability of the mRNA folded structure. Irrespective of the protein secondary structural type, the highly expressed genes always tend to encode cheaper amino acids in order to reduce the overall biosynthetic cost of production of the corresponding protein. This study supports the hypothesis that the tRNA abundance is a consequence of and not a reason for the biased usage of amino acid between highly and lowly expressed genes.     

Sequence alignment approach to pick up conformationally similar protein fragments.

Crystal structure data of globular proteins were used to prepare (phi, psi) probability maps of 20 proteinous amino acids. These maps were compared grid-wise with each other and a conformational similarity index was calculated for each pair of amino acids. A weight matrix, called Conformational Similarity Weight (CSW) matrix, was prepared using the conformational similarity index. This weight matrix was used to align sequences of 21 pairs of proteins whose crystal structures are known. The aligned regions with more than seven contiguous amino acids were further analysed by plotting average weight (W) values of overlapping hepatapeptides in these regions and carrying out curve fitting by Fourier series having TEN harmonics. The protein fragments corresponding to the half-linewidth of peaks were predicted as fragments having similar conformation in the protein pair under consideration. Such an approach allows us to pick up conformationally similar protein fragments with more than 67% accuracy.     

Structural adaptation of enzymes to low temperatures

A systematic comparative analysis of 21 psychrophilic enzymes belonging to different structural families from prokaryotic and eukaryotic organisms is reported. The sequences of these enzymes were multiply aligned to 427 homologous proteins from mesophiles and thermophiles. The net flux of amino acid exchanges from meso/thermophilic to psychrophilic enzymes was measured. To assign the observed preferred exchanges to different structural environments, such as secondary structure, solvent accessibility and subunit interfaces, homology modeling was utilized to predict the secondary structure and accessibility of amino acid residues for the psychrophilic enzymes for which no experimental three-dimensional structure is available. Our results show a clear tendency for the charged residues Arg and Glu to be replaced at exposed sites on alpha-helices by Lys and Ala, respectively, in the direction from 'hot' to 'cold' enzymes. Val is replaced by Ala at buried regions in alpha-helices. Compositional analysis of psychrophilic enzymes shows a significant increase in Ala and Asn and a decrease in Arg at exposed sites. Buried sites in beta-strands tend to be depleted of VAL: Possible implications of the observed structural variations for protein stability and engineering are discussed.     

Approaching a complete classification of protein secondary structure

A complete classification of types of the protein secondary structure is developed on the basis of computer analysis of the crystallographic structural data deposited in the protein Data Bank. The majority of amino acid residues fall into five conformation types. A conclusion is drawn that the number of sequence variants of torsion angles phi, psi in globular proteins is limited and is essentially less than the number of possible amino acid sequences for this chain length. Along with alpha-helix and beta-structure, the distribution analysis assigning every maximum of distribution of amino acid conformations on Ramachandran map to a certain type of the secondary structure exposed a third type of the secondary structure that was previously neglected. This type of the structure is extended left-handed helical conformation, designated as mobile (M-) conformation. A full set of M-conformation fragments that seems to play a major role in protein globule dynamics has been obtained, a small radius of correlation for the polypeptide chain in M-conformation is demonstrated. It explains a prevalence of short segments of mobile conformation revealed in globular proteins. For secondary structure types, the frequency of occurrence of amino acid residues has been computed.     

Analysis of the amino acid sequence of homologous proteins on the personal computer "Iskra-226"

The presented software package allows: 1) to input and edit amino acid sequences; 2) to list aligned amino acid sequences of homologous proteins; 3) pairwise comparison of homologous sequences; 4) construction of phylogenetic trees; 5) comparison of two groups of protein sequences from the same family of homologous proteins; 6) graphic identification of conservative and variable regions of homologous sequences. The stepwise application of the programs allows to study the process of amino acid replacement accumulation during certain intervals of species evolution.     

Exploring the sequence patterns in the alpha-helices of proteins

This paper reports an extensive sequence analysis of the alpha-helices of proteins. alpha-Helices were extracted from the Protein Data Bank (PDB) and were divided into groups according to their sizes. It was found that some amino acids had differential propensity values for adopting helical conformation in short, medium and long alpha-helices. Pro and Trp had a significantly higher propensity for helical conformation in short helices than in medium and long helices. Trp was the strongest helix conformer in short helices. Sequence patterns favoring helical conformation were derived from a neighbor-dependent sequence analysis of proteins, which calculated the effect of neighboring amino acid type on the propensity of residues for adopting a particular secondary structure in proteins. This method produced an enhanced statistical significance scale that allowed us to explore the positional preference of amino acids for alpha-helical conformations. It was shown that the amino acid pair preference for alpha-helix had a unique pattern and this pattern was not always predictable by assuming proportional contributions from the individual propensity values of the amino acids. Our analysis also yielded a series of amino acid dyads that showed preference for alpha-helix conformation. The data presented in this study, along with our previous study on loop sequences of proteins, should prove useful for developing potential 'codes' for recognizing sequence patterns that are favorable for specific secondary structural elements in proteins.     

Nucleotide and amino acid sequence comparisons of preprosomatostatins

The amino acid and nucleotide sequences have been aligned for five preprosomatostatins: two from catfish, two from anglerfish, and one from rat. The physical characteristics of the aligned residues as well as substitutions in the nucleotide sequences suggest considerable mutability in one of the catfish and anglerfish precursors while three of the preprohormones which give rise to somatostatins-14 are closely related. Although there is a considerable number of amino acid substitutions within the aligned somatostatin-14 precursors, there is a high degree of structural relatedness among them. These results suggest that additional physiological roles for the amino acid sequences outside of the hormone region are likely. The predicted secondary structures of the three somatostatin-14 precursors are dominated by helical and turn configurations which are correlated with observed co- and post-translational processing of the preprohormones.     

Identification and characterization of functionally important elements in the multidrug resistance protein 1 COOH-terminal region

The ATP binding cassette (ABC) transporter, multidrug resistance protein 1 (MRP1/ABCC1), transports a broad spectrum of conjugated and unconjugated compounds, including natural product chemotherapeutic agents. In this study, we have investigated the importance of the COOH-terminal region of MRP1 for transport activity and basolateral plasma membrane trafficking. The COOH-terminal regions of some ABCC proteins have been implicated in protein trafficking, but the function of this region of MRP1 has not been defined. In contrast to results obtained with other ABCC proteins, we found that the COOH-proximal 30 amino acids of MRP1 can be removed without affecting trafficking to basolateral membranes. However, the truncated protein is inactive. Furthermore, removal of as few as 4 COOH-terminal amino acids profoundly decreases transport activity. Although amino acid sequence conservation of the COOH-terminal regions of ABC proteins is low, secondary structure predictions indicate that they consist of a broadly conserved helix-sheet-sheet-helix-helix structure. Consistent with a conservation of secondary and tertiary structure, MRP1 hybrids containing the COOH-terminal regions of either the homologous MRP2 or the distantly related P-glycoprotein were fully active and trafficked normally. Using mutated proteins, we have identified structural elements containing five conserved hydrophobic amino acids that are required for activity. We show that these are important for binding and hydrolysis of ATP by nucleotide binding domain 2. Based on crystal structures of several ABC proteins, we suggest that the conserved amino acids may stabilize a helical bundle formed by the COOH-terminal three helices and may contribute to interactions between the COOH-terminal region and the protein's two nucleotide binding domains.