Substrate specificity in HIV-1 protease by a biased sequence search method

Drug resistance in HIV-1 protease can also occasionally confer a change in the substrate specificity. Through the use of computational techniques, a relationship can be determined between the substrate sequence and three-dimensional structure of HIV-1 protease, and be utilized to predict substrate specificity. In this study, we introduce a biased sequence search threading (BSST) methodology to analyze the preferences of substrate positions and correlations between them that might also identify which positions within known substrates can likely tolerate sequence variability and which cannot. The potential sequence space was efficiently explored using a low-resolution knowledge-based scoring function. The low-energy substrate sequences generated by the biased search are correlated with the natural substrates. Octameric sequences were predicted using the probabilities of residue positions in the sequences generated by BSST in three ways: considering each position in the substrate independently, considering pairwise interdependency, and considering triple-wise interdependency. The prediction of octameric sequences using the triple-wise conditional probabilities produces the most accurate results, reproducing most of the sequences for five of the nine natural substrates and implying that there is a complex interdependence between the different substrate residue positions. This likely reflects that HIV-1 protease recognizes the overall shape of the substrate more than its specific sequence.     

GB virus B and hepatitis C virus NS3 serine proteases share substrate specificity.

GB virus B (GBV-B) is a recently discovered virus responsible for hepatitis in tamarins (Saguinus species). GBV-B belongs to the Flaviviridae family and is closely related to the human pathogen hepatitis C virus (HCV). Nonstructural protein 3 (NS3) of HCV has been shown to encompass a serine protease domain required for viral maturation. GBV-B and HCV share only about 30% of the amino acid sequence within the NS3 protease domain. The catalytic triad is conserved, and the residue Phe-154, presumed to be a crucial amino acid for determining the S1 specificity pocket of the HCV NS3 protease, is also conserved. We have expressed a synthetic gene encoding the GBV-B NS3 protease domain in Escherichia coli and have characterized the purified recombinant protein for its activity on HCV substrates. We have shown that the NS3 region of the GBV-B genome actually encodes a serine protease that, despite the low sequence homology, shares substrate specificity with the HCV NS3 protease.     

Virus-specific cofactor requirement and chimeric hepatitis C virus/GB virus B nonstructural protein 3

GB virus B (GBV-B) is closely related to hepatitis C virus (HCV) and causes acute hepatitis in tamarins (Saguinus species), making it an attractive surrogate virus for in vivo testing of anti-HCV inhibitors in a small monkey model. It has been reported that the nonstructural protein 3 (NS3) serine protease of GBV-B shares similar substrate specificity with its counterpart in HCV. Authentic proteolytic processing of the HCV polyprotein junctions (NS4A/4B, NS4B/5A, and NS5A/5B) can be accomplished by the GBV-B NS3 protease in an HCV NS4A cofactor-independent fashion. We further characterized the protease activity of a full-length GBV-B NS3 protein and its cofactor requirement using in vitro-translated GBV-B substrates. Cleavages at the NS4A/4B and NS5A/5B junctions were readily detectable only in the presence of a cofactor peptide derived from the central region of GBV-B NS4A. Interestingly, the GBV-B substrates could also be cleaved by the HCV NS3 protease in an HCV NS4A cofactor-dependent manner, supporting the notion that HCV and GBV-B share similar NS3 protease specificity while retaining a virus-specific cofactor requirement. This finding of a strict virus-specific cofactor requirement is consistent with the lack of sequence homology in the NS4A cofactor regions of HCV and GBV-B. The minimum cofactor region that supported GBV-B protease activity was mapped to a central region of GBV-B NS4A (between amino acids Phe22 and Val36) which overlapped with the cofactor region of HCV. Alanine substitution analysis demonstrated that two amino acids, Val27 and Trp31, were essential for the cofactor activity, a finding reminiscent of the two critical residues in the HCV NS4A cofactor, Ile25 and Ile29. A model for the GBV-B NS3 protease domain and NS4A cofactor complex revealed that GBV-B might have developed a similar structural strategy in the activation and regulation of its NS3 protease activity. Finally, a chimeric HCV/GBV-B bifunctional NS3, consisting of an N-terminal HCV protease domain and a C-terminal GBV-B RNA helicase domain, was engineered. Both enzymatic activities were retained by the chimeric protein, which could lead to the development of a chimeric GBV-B virus that depends on HCV protease function.     

Sequence analysis of NS3 protease gene in clinical strains of hepatitis C virus

The amino terminal region of the non structural gene 3 (NS3) of hepatitis C virus (HCV) is a chymotripsinlike serine-protease responsible for cleavage of the non structural proteins of Hepatitis C virus (HCV). In order to investigate the genetic variation of this region, we developed a nested PCR to obtain NS3 protease sequences from 54 patients chronically infected with HCV genotypes 1a, 1b and 3, respectively. Comparison of nucleotide and amino acids sequences of NS3 protease domain with consensus sequence obtained within the same genotype, showed 3.73% nucleotide divergence and 1.64% amino acid divergence in isolates of genotype 3a, whereas isolates 1a exhibited 4.45% nucleotide and 4% amino acid change, respectively. Finally, NS3 sequence from 1b isolates revealed 6.47% nucleotide and 3.5 % aa changes. Comparison of consensus amino acid sequences derived from isolates 1a, 1b and 3, with the HCV prototypes showed a low amino acid sequence diversity. However, the consensus sequence of HCV genotype 3 isolates showed an amino acid changed from the prototype, that was located within a region important for enzyme structure and activity. These results indicated that the NS3 protease gene is highly conserved within the same HCV genotype. The domains involved in enzyme function were highly conserved in 1a and 1b strains, whereas consensus sequence of isolates 3a showed that the majority of these strains were not perfectly conserved in one of such regions. These findings altogether suggested that the NS3 protease enzyme of HCV may constitute an important target for antiviral therapy, but the NS3 protease variability of isolates 3 within a region that is a potential target for antiviral therapy could pose a problem for structure based drug development.     

Probing the substrate specificity of the dengue virus type 2 NS3 serine protease by using internally quenched fluorescent peptides

The NS3 (dengue virus non-structural protein 3) serine protease of dengue virus is an essential component for virus maturation, thus representing an attractive target for the development of antiviral drugs directed at the inhibition of polyprotein processing. In the present study, we have investigated determinants of substrate specificity of the dengue virus NS3 protease by using internally quenched fluorogenic peptides containing Abz (o-aminobenzoic acid; synonymous to anthranilic acid) and 3-nitrotyrosine (nY) representing both native and chimaeric polyprotein cleavage site sequences. By using this combinatorial approach, we were able to describe the substrate preferences and determinants of specificity for the dengue virus NS2B(H)-NS3pro protease. Kinetic parameters (kcat/K(m)) for the hydrolysis of peptide substrates with systematic truncations at the prime and non-prime side revealed a length preference for peptides spanning the P4-P3' residues, and the peptide Abz-RRRRSAGnY-amide based on the dengue virus capsid protein processing site was discovered as a novel and efficient substrate of the NS3 protease (kcat/K(m)=11087 M(-1) x s(-1)). Thus, while having confirmed the exclusive preference of the NS3 protease for basic residues at the P1 and P2 positions, we have also shown that the presence of basic amino acids at the P3 and P4 positions is a major specificity-determining feature of the dengue virus NS3 protease. Investigation of the substrate peptide Abz-KKQRAGVLnY-amide based on the NS2B/NS3 polyprotein cleavage site demonstrated an unexpected high degree of cleavage efficiency. Chimaeric peptides with combinations of prime and non-prime sequences spanning the P4-P4' positions of all five native polyprotein cleavage sites revealed a preponderant effect of non-prime side residues on the K(m) values, whereas variations at the prime side sequences had higher impact on kcat.     

Performance of random forest when SNPs are in linkage disequilibrium

Background

Proteases of human pathogens are becoming increasingly important drug targets, hence it is necessary to understand their substrate specificity and to interpret this knowledge in practically useful ways. New methods are being developed that produce large amounts of cleavage information for individual proteases and some have been applied to extract cleavage rules from data. However, the hitherto proposed methods for extracting rules have been neither easy to understand nor very accurate. To be practically useful, cleavage rules should be accurate, compact, and expressed in an easily understandable way.

Results

A new method is presented for producing cleavage rules for viral proteases with seemingly complex cleavage profiles. The method is based on orthogonal search-based rule extraction (OSRE) combined with spectral clustering. It is demonstrated on substrate data sets for human immunodeficiency virus type 1 (HIV-1) protease and hepatitis C (HCV) NS3/4A protease, showing excellent prediction performance for both HIV-1 cleavage and HCV NS3/4A cleavage, agreeing with observed HCV genotype differences. New cleavage rules (consensus sequences) are suggested for HIV-1 and HCV NS3/4A cleavages. The practical usability of the method is also demonstrated by using it to predict the location of an internal cleavage site in the HCV NS3 protease and to correct the location of a previously reported internal cleavage site in the HCV NS3 protease. The method is fast to converge and yields accurate rules, on par with previous results for HIV-1 protease and better than previous state-of-the-art for HCV NS3/4A protease. Moreover, the rules are fewer and simpler than previously obtained with rule extraction methods.

Conclusion

A rule extraction methodology by searching for multivariate low-order predicates yields results that significantly outperform existing rule bases on out-of-sample data, but are more transparent to expert users. The approach yields rules that are easy to use and useful for interpreting experimental data.     

Identification of Peptide Leads to Inhibit Hepatitis C Virus: Inhibitory Effect of Plectasin Peptide Against Hepatitis C Serine Protease

The emerging of hepatitis C virus (HCV) resistant strains has been considered as a main drawback of the available drugs. Since HCV has a large inactive surface, we would like to hypothesis that the mutation occur in HCV is minimal and causing less resistance against inhibition. In this study, a short peptide inhibitor of HCV namely plectasin was identified by HCV NS3-4A serine protease assay. Plectasin peptide showed considerable inhibition against HCV NS3-4A serine protease. Enzymatic activity of the recombinant NS3-4Apro was analysed by fluorescence release from several fluorogenic peptide substrates which resembling the dibasic cleavage site sequences of the flavivirus polyprotein precursor. Of all amc-labelled peptides, Pyr-RTKR-amc was the most efficiently cleaved substrate with the lowest Km value of 20 µM. The kinetic assay showed that plectasin peptide inhibited NS3-4Apro activity with an IC₅₀ value of 4.3 μM compared to the aprotinin as a standard proteases inhibitor with an IC₅₀ of 6.1 μM. From the results, plectasin peptide also demonstrated a dose-dependent inhibition of HCV replication with a considerable reduction in RLuc activity at 15 µM using HCV replicon- containing Huh-7 cells. Our study has identified a unique natural peptide that can be used to highlight novel structures for the development of drug derivatives with high efficacy of HCV NS3-4A protease inhibitors.     

Functional profiling of recombinant NS3 proteases from all four serotypes of dengue virus using tetrapeptide and octapeptide substrate libraries

Regulated proteolysis by the two-component NS2B/NS3 protease of dengue virus is essential for virus replication and the maturation of infectious virions. The functional similarity between the NS2B/NS3 proteases from the four genetically and antigenically distinct serotypes was addressed by characterizing the differences in their substrate specificity using tetrapeptide and octapeptide libraries in a positional scanning format, each containing 130,321 substrates. The proteases from different serotypes were shown to be functionally homologous based on the similarity of their substrate cleavage preferences. A strong preference for basic amino acid residues (Arg/Lys) at the P1 positions was observed, whereas the preferences for the P2-4 sites were in the order of Arg > Thr > Gln/Asn/Lys for P2, Lys > Arg > Asn for P3, and Nle > Leu > Lys > Xaa for P4. The prime site substrate specificity was for small and polar amino acids in P1' and P3'. In contrast, the P2' and P4' substrate positions showed minimal activity. The influence of the P2 and P3 amino acids on ground state binding and the P4 position for transition state stabilization was identified through single substrate kinetics with optimal and suboptimal substrate sequences. The specificities observed for dengue NS2B/NS3 have features in common with the physiological cleavage sites in the dengue polyprotein; however, all sites reveal previously unrecognized suboptimal sequences.     

HCV RNA-dependent RNA polymerase replicates in vitro the 3' terminal region of the minus-strand viral RNA more efficiently than the 3' terminal region of the plus RNA.

The NS5B protein, or RNA-dependent RNA polymerase of the hepatitis virus type C, catalyzes the replication of the viral genomic RNA. Little is known about the recognition domains of the viral genome by the NS5B. To better understand the initiation of RNA synthesis on HCV genomic RNA, we used in vitro transcribed RNAs as templates for in vitro RNA synthesis catalyzed by the HCV NS5B. These RNA templates contained different regions of the 3' end of either the plus or the minus RNA strands. Large differences were obtained depending on the template. A few products shorter than the template were synthesized by using the 3' UTR of the (+) strand RNA. In contrast the 341 nucleotides at the 3' end of the HCV minus-strand RNA were efficiently copied by the purified HCV NS5B in vitro. At least three elements were found to be involved in the high efficiency of the RNA synthesis directed by the HCV NS5B with templates derived from the 3' end of the minus-strand RNA: (a) the presence of a C residue as the 3' terminal nucleotide; (b) one or two G residues at positions +2 and +3; (c) other sequences and/or structures inside the following 42-nucleotide stretch. These results indicate that the 3' end of the minus-strand RNA of HCV possesses some sequences and structure elements well recognized by the purified NS5B.     

Effect of naturally occurring active site mutations on hepatitis C virus NS3 protease specificity

A comparison of the DNA sequences from all available genotypes of HCV indicate that the active site residues of the NS3 protease are strictly conserved with the exception of positions 123 and 168, which border the S(4) subsite. In genotype 3, the canonic arginine and aspartic acid have been replaced with threonine and glutamine, respectively. To determine if these differences contribute to an altered specificity, we characterized single-chain NS3 proteases from strains 1a, 1b, and 3a with peptide substrates and product inhibitors on the basis of the natural cleavage junction sequences, in addition to polyprotein substrates derived from the 1a strain. No statistically significant differences in specificity were observed. To demonstrate that the active sites were actually different, we generated and evaluated peptide substrates with unnatural extended side-chains. These studies confirmed that there are measurable differences between the NS3 proteases of genotypes 1 and 3. Specifically, a 5-fold difference in K(i) was observed between the proteases from genotypes 1 and 3 when a D-Glu occupied P(5), and a 30-fold difference was seen when this position contained a D-homoglutamate. The contribution of residues 123 and 168 toward the altered specificity was then evaluated individually by site-directed mutagenesis. These mutants showed that potency differences within this series could be attributed to the residue that occupied position 123 of the protease. Modeling these unnatural substrate/mutant protease interactions, on the basis of cocrystal structures of enzyme-substrate complexes, provides a structural basis for these observations. Proteins 2001;43:82-88.     

Selection of 3'-template bases and initiating nucleotides by hepatitis C virus NS5B RNA-dependent RNA polymerase

De novo RNA synthesis by hepatitis C virus (HCV) nonstructural protein 5B (NS5B) RNA-dependent RNA polymerase has been investigated using short RNA templates. Various templates including those derived from the HCV genome were evaluated by examining the early steps of de novo RNA synthesis. NS5B was shown to be able to produce an initiation dinucleotide product from templates as short as 4-mer and from the 3'-terminal sequences of both plus and minus strands of the HCV RNA genome. GMP, GDP, and guanosine were able to act as an initiating nucleotide in de novo RNA synthesis, indicating that the triphosphate moiety is not absolutely required by an initiating nucleotide. Significant amounts of the initiation product accumulated in de novo synthesis, and elongation from the dinucleotide was observed when large amounts of dinucleotide were available. This result suggests that NS5B, a template, and incoming nucleotides are able to form an initiation complex that aborts frequently by releasing the dinucleotide product before transition to an elongation complex. The transition is rate limiting. Furthermore, we discovered that the secondary structure of a template was not essential for de novo initiation and that 3'-terminal bases of a template conferred specificity in selection of an initiation site. Initiation can occur at the +1, +2, or +3 position numbered from the 3' end of a template depending on base composition. Pyrimidine bases at any of the three positions are able to serve as an initiation site, while purine bases at the +2 and +3 positions do not support initiation. This result implies that HCV possesses an intrinsic ability to ensure that de novo synthesis is initiated from the +1 position and to maintain the integrity of the 3' end of its genome. This assay system should be an important tool for investigating the detailed mechanism of de novo initiation by HCV NS5B as well as other viral RNA polymerases.     

Assessment of substrate specificity of hepatitis G virus NS3 protease by a genetic method

The RNA genome of hepatitis G virus (HGV) encodes a large polyprotein that is processed to mature proteins by viral-encoded proteases. The HGV NS3 protease is responsible for the cleavage of the HGV polyprotein at four different locations. No conserved sequence motif has been identified for the cleavage sites of the NS3 protease. To determine the substrate specificity of the NS3 protease, amino acid sequences cleaved by the NS3 protease were obtained from randomized sequence libraries by using a screening method referred to as GASP (Genetic Assay for Site-specific Proteolysis). Based on statistical analyses of the obtained cleavable sequences, a consensus substrate sequence was deduced: Gln-Glu-Thr-Leu-Val downward arrow Ser, with the scissile bond located between Val and Ser. The relevance of this peptide as a cleavable substrate was further supported by molecular modeling of the NS3 protease. Our result would provide an insight on the molecular activity of the NS3 protease and may be useful for the design of substrate-based inhibitors.     

丙型肝炎病毒丝氨酸蛋白酶特异性结合肽的筛选和鉴定

丙型肝炎病毒丝氨酸蛋白酶在病毒复制和包装中的重要作用使其成为特异性抗病毒药物研究的首选靶标。根据丝氨酸蛋白酶晶体结构特点,用柔性连接子连接NS3丝氨酸蛋白酶结构域和NS4A的核心序列,构建成单链丝氨酸蛋白酶基因并且在大肠杆菌中获得高水平的可溶性表达,纯化后的目的蛋白能够切割重组蛋白底物NS5ab。随后,以单链丝氨酸蛋白酶为靶分子对噬菌体展示的随机十二肽库进行了三轮淘筛,挑选的44个克隆中有37个克隆能够特异性地结合丝氨酸蛋白酶,并且这种结合作用为竞争性ELISA试验结果所支持。对13个克隆进行序列测定,得到6种序列,它们在氨基酸组成上存在明显偏性,富含组氨酸和色氨酸,缺乏酸性氨基酸;6种序列存在一个共有序列。     

Construction, expression, and characterization of a novel fully activated recombinant single-chain hepatitis C virus protease.

Efficient proteolytic processing of essential junctions of the hepatitis C virus (HCV) polyprotein requires a heterodimeric complex of the NS3 bifunctional protease/helicase and the NS4A accessory protein. A single-chain recombinant form of the protease has been constructed in which NS4A residues 21-32 (GSVVIVGRIILS) were fused in frame to the amino terminus of the NS3 protease domain (residues 3-181) through a tetrapeptide linker. The single-chain recombinant protease has been overexpressed as a soluble protein in E. coli and purified to homogeneity by a combination of metal chelate and size-exclusion chromatography. The single-chain recombinant protease domain shows full proteolytic activity cleaving the NS5A-5B synthetic peptide substrate, DTEDVVCCSMSYTWTGK with a Km and k(cat) of 20.0 +/- 2.0 microM and 9.6 +/- 2.0 min(-1), respectively; parameters identical to those of the authentic NS3(1-631)/NS4A(1-54) protein complex generated in eukaryotic cells (Sali DL et al., 1998, Biochemistry 37:3392-3401).     

Hepatitis C virus nonstructural 5A protein and interferon resistance: a new model for testing the reliability of mutational analyses

Presumably due to the capability of the hepatitis C virus (HCV) to evade the antiviral effects of alpha interferon, treatment is ineffective in more than half of chronically genotype HCV type 1 (HCV-1)-infected patients. Previous approaches to correlate the number of amino acid mutations within regions of HCV nonstructural (NS)-5A protein with virologic treatment response provided conflicting results. In the present study, we developed a new mathematical model to investigate NS5A sequences of HCV-1-infected patients. The mean number of all mutations within the complete NS5A protein was significantly higher in virologic responders compared to nonresponders (P = 0.008 and P = 0.0001 for amino acid residues predicted on the surface of the NS5A protein). Differences did not achieve statistical significance for NS5A regions that are currently assumed to be functionally relevant (e.g., the interferon sensitivity-determining region, the RNA-activated protein kinase-binding domain, etc.). Analyses of smoothed mutational frequencies showed that the number of mutations in other NS5A regions correlated with virologic response. Such a correlation was observed for both genuine and randomly generated NS5A sequences. The existence of local accumulations of mutations within genuine NS5A isolates that truly correlated with treatment response was defined by a refined test procedure. Upon considering the predicted residue accessibility, we identified the main focus of mutations correlating with treatment response to be the sequence from amino acids 2350 to 2370. Thus, evaluation of NS5A mutations in correlation with treatment response is improved by consideration of functional and predicted conformational amino acid properties. As shown by simulations with randomly generated sequences, multiple analyses of simple counts of local NS5A amino acid mutations and correlation with treatment response are insufficient. For improvement of mutational analysis, a refined specific functional data test procedure is proposed.     

HCV NS3/4A蛋白酶与Faldaprevir类似物的分子识别机制及其复合物运动模式

Faldaprevir类似物(Faldaprevir analogue molecule,FAM)能有效抑制HCV NS3/4A蛋白酶的催化活性,是一种潜在抗HCV先导化合物。通过生物信息学统计分析了已报道的HCV NS3/4A蛋白酶晶体结构,得到了FAM-HCV NS3/4A蛋白酶晶体结构。对FAM-HCV NS3/4A蛋白酶复合物进行了20.4 ns的分子动力学模拟,重点从氢键和结合自由能两个角度分析了二者分子识别中的关键残基及结合驱动力。氢键和范德华力是促使FAM特异性结合到蛋白V132?S139、F154?D168、D79?D81和V55的活性口袋中的主要驱动力,这与实验数据较为吻合。耐药性突变实验分析了R155K、D168E/V和V170T定点突变对FAM分子识别的影响,为可能存在的FAM耐药性提供了分子依据。最后,用自由能曲面和构象聚类两个方法探讨了体系的构象变化,给出体系的4种优势构象,为后续的基于HCV NS3/4A蛋白酶结构的Faldaprevir类似物抑制剂分子设计提供一定的理论帮助。     

Chimeric Sindbis viruses dependent on the NS3 protease of hepatitis C virus.

The hepatitis C virus (HCV) NS3 protease cleaves the viral polyprotein at specific sites to release the putative components of the HCV replication machinery. Selective inhibition of this enzyme is predicted to block virus replication, and NS3 is thus considered an attractive candidate for development of anti-HCV therapeutics. To set up a system for analysis of NS3 protease activity in cultured cells, we constructed a family of chimeric Sindbis viruses which carry sequences coding for NS3 and its activator, NS4A, in their genomes. HCV sequences were fused to the gene coding for the Sindbis virus structural polyprotein via an NS3-specific cleavage site, with the expectation that processing of the chimeric polyprotein, nucleocapsid assembly, and generation of viable viral particles would occur only upon NS3-dependent proteolysis. Indeed, the chimeric genomes encoding an active NS3 protease produced infectious viruses in mammalian cells, while those encoding NS3 inactivated by alanine substitution of the catalytic serine did not. However, in infected cells chimeric genomes recombined, splicing out HCV sequences and reverting to pseudo-wild-type Sindbis virus. To force retention of HCV sequences, we modified one of the initial chimeras by introducing a second NS3 cleavage site in the Sindbis virus portion of the recombinant polyprotein, anticipating that revertants not encoding an active NS3 protease would not be viable. The resulting chimera produced infectious viruses which replicated at a lower rate than the parental construct and displayed a marked temperature dependence in the formation of lysis plaques yet stably expressed NS3.     

Comparative molecular model building of two serine proteinases from cytotoxic T lymphocytes

Two genes that are expressed when precursor cytotoxic T lymphocytes are transformed to T killer cells have been cloned and sequenced. The derived amino acid sequences, coding for cytotoxic cell protease 1 (CCP1) and Hannuka factor (HF) are highly homologous to members of the serine proteinase family. Comparative molecular model building using the known three-dimensional structures and the derived amino acid sequences of the lymphocyte enzymes has provided useful structural information, especially in predicting the conformations of the substrate binding sites. In applying this modelling procedure, we used the X-ray structures of four serine proteinases to provide a structurally based sequence alignment: alpha-chymotrypsin (CHT), bovine trypsin (BT), Streptomyces griseus trypsin (SGT), and rat mast cell protease 2 (RMCP2). The root mean square differences in alpha-carbon atom positions among these four structures when compared in a pairwise fashion range from 0.79 to 0.97 A for structurally equivalent residues. The sequences of the two lymphocyte enzymes were then aligned to these proteinases using chemical criteria and the superimposed X-ray structures as guides. The alignment showed that the sequence of CCP1 was most similar to RMCP2, whereas HF has regions of homology with both RMCP2 and BT. With RMCP2 as a template for CCP1 and the two enzymes RMCP2 and BT as templates for HF, the molecular models were constructed. Intramolecular steric clashes that resulted from the replacement of amino acid side chains of the templates by the aligned residues of CCP1 and HF were relieved by adjustment of the side chain conformational angles in an interactive computer graphics device. This process was followed by energy minimization of the enzyme model to optimize the stereochemical geometry and to relieve any remaining unacceptably close nonbonded contacts. The resulting model of CCP1 has an arginine residue at position 226 in the specificity pocket, thereby predicting a substrate preference for P1 aspartate or glutamate residues. The model also predicts favorable binding for a small hydrophobic residue at the P2 position of the substrate. The primary specificity pocket of HF resembles that of BT and therefore predicts a lysine or arginine preference for the P1 residue. The arginine at position 99 in the model of HF suggests a preference for aspartate or glutamate side chains in the P2 position of the substrate. Both CCP1 and HF have a free cysteine in the segment of polypeptide 88 to 93.(ABSTRACT TRUNCATED AT 400 WORDS)