A structural model for the retroviral proteases

Active foamy virus proteinase is essential for virus infectivity but not for formation of a Pol polyprotein

To analyze proteolytic processing of foamy (spuma) retroviruses, two mutations were generated in the presumed active-site triplet Asp-Ser-Gly in the predicted proteinase (PR) region of the human foamy virus (HSRV). The mutations changed either the presumed catalytic aspartic acid residue to a catalytically incompetent alanine or the adjacent serine to a threonine found in most cellular and retroviral proteases at this position. Both mutations were cloned into the full-length infectious HSRV DNA clone. Wild-type and S/T mutant genomes directed the synthesis of particles with similar infectious titers, while the HSRV D/A PR mutant was noninfectious. Immunoblot analysis of transfected cells revealed identical patterns for the wild-type and for the S/T PR mutant. HSRV D/A mutant-transfected cells expressed only a single Gag polyprotein of 78 kDa instead of the 78-kDa-74-kDa doublet found in HSRV-infected or wild-type-transfected cells. Analysis with pol-specific antisera yielded a protein of approximately 120 kDa reactive with antisera against pol- but not gag-specific domains. No Gag-Pol polyprotein was detected in this study. Electron microscopy analysis of transfected cells showed heterogeneous particle morphology in the case of the D/A mutant, with particles of normal appearance and particles of aberrant size and shape. These results indicate that foamy viruses have an aspartic PR that is essential for infectivity but not for formation of the 120-kDa Pol polyprotein.

Trans-dominant inhibitory human immunodeficiency virus type 1 protease monomers prevent protease activation and virion maturation

Production of infectious human immunodeficiency virus (HIV) requires proper polyprotein processing by the dimeric viral protease. The trans-dominant inhibitory activity of a defective protease monomer with the active site Asp-25 changed to Asn was measured by transient transfection. A proviral plasmid that included the drug-selectable Escherichia coli gpt gene was used to deliver the wild-type (wt) or mutant proteases to cultured cells. Coexpression of the wt proviral DNA (HIV-gpt) with increasing amounts of the mutant proviral DNA (HIV-gpt D25N) results in a concomitant decrease in proteolytic activity monitored by in vivo viral polyprotein processing. The viral particles resulting from inactivation of the protease were mostly immature, consisting predominantly of unprocessed p55gag and p160gag-pol polyproteins. In the presence of HIV-1 gp160 env, the number of secreted noninfectious particles correlated with the presence of increasing amounts of the defective protease. Greater than 97% reduction in infectivity was observed at a 1:6 ratio of wt to defective protease DNA. This provides an estimate of the level of inhibition required for effectively preventing virion processing. Stable expression of the defective protease in monkey cells reduced the yield of infectious particles from these cells by 90% upon transfection with the wt proviral DNA. These results show that defective subunits of the viral protease exert a trans-dominant inhibitory effect resulting from the formation of catalytically compromised heterodimers in vivo, ultimately yielding noninfectious viral particles.

Sigma factor and sporulation genes in Clostridium

The genus Clostridium, represented by Gram-positive, anaerobic, spore-forming bacteria, is well known for its clinical importance and considerable biotechnological potential. Recently, evidence for a functional role of the transcription factors sigma A, sigma E, sigma G, and sigma K in this genus was provided by cloning and sequencing these genes from C. acetobutylicum. In C. kluyveri, a partially sequenced open reading frame was found to encode the N terminus of the putative sigma factor L with significant similarity to members of the sigma 54 family. The identification of sequences with high similarity to the Bacillus sigma F (C. acetobutylicum), sigma H (several clostridial species), and sigma D (C. thermocellum)-controlled consensus promoters renders the existence of these transcription factors in clostridia very likely. These data are in agreement with information obtained by RNA transcript mapping (sigma A, sigma H), heterologous DNA hybridization (sigma D, sigma H), and immuno characterization of purified proteins (sigma A) from various clostridial species. Thus, the picture emerges that a fundamental similarity exists at the genetic level between the regulation of various cellular responses, in particular sporulation, in the genera Bacillus and Clostridium. The different induction patterns of sporulation in Bacillus spp. (nutrient starvation) and many clostridial species (cessation of growth or exposure to oxygen in the presence of excess nutrients) are most interestingly not reflected in the general regulatory features of this developmental process.

HIV protease inhibitor HOE/BAY 793, structure-activity relationships in a series of C2-symmetric diols

A detailed structure-activity relationship of C2-symmetric diol inhibitors of HIV-1 protease leads to inhibitor 6 (HOE/BAY 793) which is outstanding in the inhibition of the enzyme and in the inhibition of viral replication in HIV infected cell culture (IC50: 0.3 nM; EC50: 3 nM). There are well defined steric requirements for the design of the side chains P1-P3 of the inhibitors. In addition, all three side chains need to be lipophilic. While the enzyme tolerates hydrophilic substituents in some cases, drastic reductions in anti-HIV activity are observed in cell culture, most likely due to insufficient cell penetration.

Connective tissue and repair in the heart. Potential regulatory mechanisms

The heart is composed of highly differentiated cardiac myocytes, which constitute parenchyma, and stroma or connective tissue. Fibrillar collagen turnover in the heart and its valve leaflets, in particular, is dynamic and essential to tissue repair. Emerging evidence further suggests connective tissue is a metabolically active entity, where peptide hormones are generated and degraded and, in turn, these peptides regulate collagen turnover. This concept arose from quantitative in vitro autoradiography using an iodinated derivative of lisinopril (125I-351A) as ligand to localize angiotensin converting enzyme (ACE) binding density within the heart. A heterogeneous distribution was found: low-density ACE binding within atria and ventricles; high ACE binding density at sites of high collagen turnover, such as valve leaflets, adventitia, and fibrous tissue of diverse etiologic origins. ACE-producing cells at these latter sites were identified by monoclonal ACE antibody. They included valvular interstitial cells (VIC) and fibroblast-like cells each of which also contained alpha-smooth muscle actin and the transcript for type I collagen (in situ hybridization). Substrate utilization in cultured VIC was found to include angiotensin I and bradykinin. Angiotensin II and bradykinin receptor-ligand binding was observed in VIC and at fibrous tissue sites. Connective tissue ACE is independent of circulating angiotensin II. In vivo, fibrous tissue formation is attenuated by ACE inhibition or antagonism of AT1 receptor. Angiotensin II and bradykinin are stimulatory and inhibitory, respectively, to cultured adult cardiac fibroblast collagen synthesis suggesting a paradigm of reciprocal regulation to fibroblast collagen turnover. Stroma and its cellular constituents represent a dynamic metabolic entity that regulates its own peptide hormone composition and turnover of fibrillar collagen. These findings may provide insights that could be used to advantage to either promote or forestall fibrous tissue formation depending on the nature of cardiovascular disease.

Determination of L-735 524, an human immunodeficiency virus protease inhibitor, in human plasma and urine via high-performance liquid chromatography with column switching

A method for the determination of an HIV protease inhibitor, L-735 524, in human plasma and urine is described. Isolation of the analyte and the internal standard from the matrices was achieved via multiple liquid-liquid extractions with methyl tert.-butyl ether. The analyte lacks significant UV absorption at wavelengths greater than 220 nm, hence a column switching system using a cyano and C18 column was used to further purify the extracts prior to UV detection at 210 nm. The assay has been found to be linear and has been validated over the concentration range of 5 to 500 ng/ml, when 1-ml aliquots of plasma or urine were extracted. The assay has been utilized to support human pharmacokinetic studies.

Peptides Containing 2-, 3- or 4-Fluorophenylalanine or 2,2- or 3,3-Difluorophenylalanine as Potential Inhibitors of HIV Protease

A series of peptide mimetic substrates of the HIV-1 protease, corresponding substantially to the amino acid sequence of the Tyr-Pro cleavage site contained in the gag polyprotein, in which the tyrosine amino acid of the scissile bond was replaced by 2-, 3- or 4-fluorophenylalanine or 2,2- or 3,3-difluorophenylalanine, has been developed. The synthesis of these compounds and the results of inhibition studies are described.

Selection and analysis of human immunodeficiency virus type 1 variants with increased resistance to ABT-538, a novel protease inhibitor

Inhibitors of the human immunodeficiency virus protease represent a promising new class of antiretroviral drugs for the treatment of AIDS. We now report the in vitro selection of viral variants with decreased sensitivity to a symmetry-based protease inhibitor, ABT-538, currently being tested in clinical trials. Molecular characterization of the variants shows that an isoleucine-to-valine substitution at position 84 results in a substantial decrease in sensitivity to the drug. Moreover, an additional mutation at position 82, valine to phenylalanine, further decreases viral susceptibility to ABT-538. Three-dimensional analysis of the protease-drug complex provides a structural explanation for the relative drug resistance induced by these two mutations. These findings emphasize the importance of closely monitoring patients receiving ABT-538 for the emergence of viral resistance and provide information that may prove useful in designing the next generation of protease inhibitors.

Structure-assisted design of nonpeptide human immunodeficiency virus-1 protease inhibitors

The discovery of therapeutic agents has, in the past, started primarily with random screening efforts. These screens, although effective, are time-consuming and expensive. Attempts are now being made to design more efficient methods of screening that take into account available information about the three-dimensional structure of a target receptor or enzyme. In the case of acquired immunodeficiency syndrome, the structure of a proteolytic enzyme, the human immunodeficiency virus (HIV) protease, an aspartyl protease that plays a crucial role in the viral life cycle, has been determined and extensively characterized. Using the protease structure and the computer program DOCK, the active site of the protease was mapped and its shape used to screen a subset of the Cambridge Structural Database. Among the molecules whose shape was complementary to the active site was the antipsychotic agent haloperidol. This molecule and several chemically modified derivatives were shown to bind competitively with micromolar affinity to the HIV protease but not to cellular aspartyl proteases. X-ray structures of the HIV protease complexed with haloperidol derivatives show the molecules binding in the predicted position at the active site. In an attempt to overcome the problems associated with low-affinity competitive inhibitors, reactive groups that enable the molecule to serve as a specific irreversible inhibitor of the HIV protease were introduced onto the haloperidol scaffold. These inhibitors demonstrated an ability to block viral polyprotein processing in a tissue culture model of HIV-1 infection, although their cytotoxicity is pronounced.(ABSTRACT TRUNCATED AT 250 WORDS)

Tas, a retrotransposon from the parasitic nematode Ascaris lumbricoides

The cloned retrotransposon Tas OE3 from the genome of the parasitic nematode Ascaris lumbricoides was completely sequenced. The element is flanked by long terminal repeats (LTR) and contains three distinct regions encoding putative proteins typical for retroid elements. The first region, ORF1, encodes a putative Gag protein including a 'Leu zipper', a nucleic acid binding motif, as well as an aspartic protease domain. The second region contains an incomplete ORF (ORF2) with sequence similarities to known retroviral reverse transcriptases (RT), ribonucleases H and integrases. A third ORF, which is located adjacent to the 3' LTR, might encode an env-like protein. Based on amino-acid sequence analysis of the RT domain, Tas falls into a new subgroup of LTR-containing retrotransposons.

Sequence of MKT1, needed for propagation of M2 satellite dsRNA of the L-A virus of Saccharomyces cerevisiae

MKT1 is required for maintenance of K2 above 30 degrees C in strains with the L-A-HN variant of the L-A double-stranded RNA virus of Saccharomyces cerevisiae. We report that MKT1 encodes a 92,979 Da protein with serine-rich regions and the retroviral protease signature, DTG, but with no substantial homology to proteins presently in the databases.

Alpha-ketoamide Phe-Pro isostere as a new core structure for the inhibition of HIV protease

Studies on the inhibition of HIV-1 protease utilizing a core isostere with replacement of the scissle bond for an alpha-amino-ketone have resulted in the development of an alpha-keto-amide isosteric replacement of the Phe-Pro scissle amide bond. The simple dipeptide isostere was shown to be a promising new core structure for the development of the enzyme inhibitors. The Ki of this core structure was determined to be 6 microM, compared to 230 microM and > 50 microM for the corresponding phosphinic acid and hydroxyethylamine isosteres.

Sequence conservation and correlation measures in protein structure prediction

The rapid elucidation of protein sequences has allowed multiple sequence alignments to be calculated for a wide variety of proteins. Such alignments reveal positions that exhibit amino acid conservation--either of specific chemical groups in active and binding sites or of the more chemically inert hydrophobic residues that contribute to the protein core. The latter can provide constraints on the position of the protein chain and any local periodicity can suggest the type of secondary structure. Conservation measures, however, cannot provide specific pairwise packing information (each conserved hydrophobic position might pack against any other). However, if correlated changes between positions were observed then specific pairs of residue could be identified as interacting and therefore probably spatially adjacent. Most 'observations' of correlated changes have been anecdotal and of the few systematic studies that have been made, most have mistakenly incorporated a strong bias towards selecting conserved positions. When the conservation effect is separated (as best as possible) then little correlation signal remains to help identify adjacent positions.

Characterization of a monoclonal antibody produced in an attempt to mimic the active site of HIV aspartyl protease using haptens based on inhibitor models

The high binding affinity and specificity of antibodies for a great variety of ligands has been widely exploited in structure-activity relationship studies of biomolecules and more recently in the development of new catalysts for several chemical reactions. It is assumed that antibodies generated against haptenic protease inhibitors would recognize both these haptens and the substrate of the model proteolytic reaction. We have produced antibodies against HIV PRp12 aspartyl protease substrate analogues, chemically modified at the scissile bond, Phe-Pro. Identical chemical modifications have been reported for related HIV protease inhibitors. We finally selected an anti-hapten monoclonal antibody that specifically recognized the substrate and those haptens with both the phenylalanyl side chain and the prolyl pyrrolidine ring. This selectivity of recognition suggests that such an antibody might mimic the catalytic site of the model protease.

Potency and selectivity of inhibition of human immunodeficiency virus protease by a small nonpeptide cyclic urea, DMP 323

DMP 323 is a potent inhibitor of the protease of human immunodeficiency virus (HIV), with antiviral activity against both HIV type 1 and HIV type 2. This compound is representative of a class of small, novel, nonpeptide cyclic urea inhibitors of HIV protease that were designed on the basis of three-dimensional structural information and three-dimensional database searching. We report here studies of the kinetics of DMP 323 inhibition of the cleavage of peptide and HIV-1 gag polyprotein substrates. DMP 323 acts as a rapidly binding, competitive inhibitor of HIV protease. DMP 323 is as potent against both peptide and viral polyprotein substrates as A-80987, Q8024, and Ro-31-8959, which are among the most potent inhibitors of HIV protease described in the literature to date. Incubation with human plasma or serum did not decrease the effective potency of DMP 323 for HIV protease, suggesting that plasma protein binding is of a low affinity relative to that of HIV protease. DMP 323 was also assessed for its ability to inhibit the mammalian proteases renin, pepsin, cathepsin D, cathepsin G, and chymotrypsin. No inhibition of greater than 12% was observed for any of these enzymes at concentrations of DMP 323 that were 350 to 40,000 times higher than that required to inhibit the viral protease 50%.

Protein structure modelling from remote sequence similarity

Many methods exist for taking a sequence that exhibits similarity to another of known structure and building a molecular model. However, when the sequence similarity is very remote and fragmentary, this 'modelling-by-homology' approach is less reliable. Current methods that tackle this problem are reviewed below, taking as an example the construction of a predicted model for the retroviral protease. This earlier work, which was only partially automatic, identified many of the outstanding difficulties that have subsequently been automated in computer programs, developed both by the author and many others. Because of the rapid proliferation of methods and their variants, an exhaustive review of the literature has not been possible and the following survey concentrates on the developments of the author and colleagues to explain the basic methods.

The recognition of protein structure and function from sequence: adding value to genome data

The explosion of DNA sequence data from genome projects presents many challenges. For instance, we must extend our current knowledge of protein structure and function so that it can be applied to these new sequences. The derivation of rules for the relationships between sequence and structure allow us to recognize a common fold by the use of tertiary templates. New techniques enable us to begin to meet the challenge of rule-based modelling of distantly related proteins. This paper describes an integrated and knowledge-based approach to the prediction of protein structure and function which can maximize the value of sequence information.

Isolation and characterization of a gene, pmrD, from Salmonella typhimurium that confers resistance to polymyxin when expressed in multiple copies

We have isolated from Salmonella typhimurium a gene, designated pmrD, that confers resistance to the membrane-damaging drug, polymyxin B when expressed from the medium-copy-number plasmid pHSG576. The gene maps to 46 min on the standard genetic map, near the menB gene, and is therefore distinct from the previously described pmrA locus. We have mapped the polymyxin resistance activity to a 1.3-kb ClaI-PvuII fragment which contains a small open reading frame that could encode an 85-amino-acid peptide. When an omega-Tet insertion was made into the putative pmrD open reading frame (pmrD2::omega-Tet), the resulting plasmid no longer conferred polymyxin resistance, whereas an omega-Tet insertion into vector sequences had no effect. Maxicell analysis confirmed that a protein of the expected size is made in vivo. The PmrD protein shows no significant homology to any known protein, but it does show limited homology across the active site of the p15 acid protease from Rous sarcoma virus, indicating that the protein may have proteolytic activity. However, changing the aspartic acid residue at the putative active site to alanine reduced but did not eliminate polymyxin resistance. When pmrD2::omega-Tet replaced the chromosomal copy of pmrD, the resulting strain showed wild-type sensitivity to polymyxin and could be complemented to resistance by a plasmid that carried pmrD. The pmrA505 allele confers resistance to polymyxin when present in single copy on the chromosome or when present on a plasmid in pmrA+ pmrD+ cells. In combination with the pmrD(2)::-Tet mutation, the effect o the pmrA505 allele on polymyxin resistance was reduced, whether pmrA505 was present in the chromosome or on a plasmid. Conversely, a strain carrying an insertion in pmrA could be complemented to polymyxin resistance by a plasmid carrying the pmrA505 allele but not by a plasmid carrying pmrD. On the basis of these results, we suggest that polymyxin resistance is mediated by an interaction between PmrA or a PmrA-regulated gene product and PmrD.

Viral gene products and replication of the human immunodeficiency type 1 virus

The acquired immunodeficiency syndrome (AIDS) epidemic represents a modern-day plague that has not only resulted in a tragic loss of people from a wide spectrum of society but has reshaped our viewpoints regarding health care, the treatment of infectious diseases, and social issues regarding sexual behavior. There is little doubt now that the cause of the disease AIDS is a virus known as the human immunodeficiency virus (HIV). The HIV virus is a member of a large family of viruses termed retroviruses, which have as a hallmark the capacity to convert their RNA genome into a DNA form that then undergoes a process of integration into the host cell chromosome, followed by the expression of the viral genome and translation of viral proteins in the infected cell. This review describes the organization of the HIV-1 viral genome, the expression of viral proteins, as well as the functions of the accessory viral proteins in HIV replication. The replication of the viral genome is divided into two phases, the early phase and the late phase. The early phase consists of the interaction of the virus with the cell surface receptor (CD4 molecule in most cases), the uncoating and conversion of the viral RNA genome into a DNA form, and the integration into the host cell chromosome. The late phase consists of the expression of the viral proteins from the integrated viral genome, the translation of viral proteins, and the assembly and release of the virus. Points in the HIV-1 life cycle that are targets for therapeutic intervention are also discussed.

L-735,524: an orally bioavailable human immunodeficiency virus type 1 protease inhibitor

To date, numerous inhibitors of the human immunodeficiency virus type 1 protease have been reported, but few have been studied extensively in humans, primarily as a consequence of poor oral bioavailability in animal models. L-735,524 represents a class of human immunodeficiency virus type 1 protease inhibitors, termed hydroxyaminopentane amides, that incorporate a basic amine into the hydroxyethylene inhibitor backbone. L-735,524 is a potent inhibitor of virus replication in cell culture and inhibits the protease-mediated cleavage of the viral precursor polyproteins that results in the production of noninfectious progeny viral particles. The compound is effective against viruses resistant to reverse transcriptase inhibitors and is synergistically active when used in combination with reverse transcriptase inhibitors. Most importantly, L-735,524 exhibits good oral bioavailability and plasma pharmacokinetic profiles in two species of laboratory animals by using clinically acceptable formulations. Accordingly, the compound was selected for evaluation of safety and pharmacokinetic studies in humans.

Structure-based identification and clustering of protein families and superfamilies

We describe an approach to protein structure comparison designed to detect distantly related proteins of similar fold, where the procedure must be sufficiently flexible to take into account the elasticity of protein folds without losing specificity. Protein structures are represented as a series of secondary structure elements, where for each element a local environment describes its relations with the elements that surround it. Secondary structures are then aligned by comparing their features and local environments. The procedure is illustrated with searches of a database of 468 protein structures in order to identify proteins of similar topology to porcine pepsin, porphobilinogen deaminase and serum amyloid P-component. In all cases the searches correctly identify protein structures of similar fold as the search proteins. Multiple cross-comparisons of protein structures allow the clustering of proteins of similar fold. This is exemplified with a clustering of alpha/beta- and beta-class protein structures. We discuss applications of the comparison and clustering of three-dimensional protein structures to comparative modelling and structure-based protein design.

Structural features of mag, a gypsy-like retrotransposon of Bombyx mori, with unusual short terminal repeats

Mag is a retrotransposon found as an insert in the Sericin 2 gene. It is present in a few copies--4 to 15--dispersed in the genome of different strains of Bombyx mori as well as in Bombyx mandarina. Flanked by a 5 bp target sequence with no sequence specificity, it is bordered by direct repeats of 77 nucleotides. Despite their unusual short size, these terminal repeats and their immediately adjacent sequences present all the signals necessary for transcription into genomic RNA and for reverse transcription. Mag contains two overlapping open reading frames which are organized as the gag and pol genes of retroviruses and encode putative nucleic acid binding peptide, protease, reverse transcriptase, RNase H and endonuclease in this order. Sequence comparison of these proteins places mag within the gypsy group of LTR retrotransposons next to the echinoderm element SURL.

Knowledge-based protein modeling

Knowledge, both from the three-dimensional structures of homologous proteins and from the general analysis of protein structure, is of value in modeling a protein of known sequence but unknown structure. While many models are still constructed at least in part by manual methods on graphics devices, automated procedures have come into greater use. These procedures include those that assemble fragments of structure from other known structures and those that derive coordinates for the model from the satisfaction of restraints placed on atomic positions.

Proteolytic activity of novel human immunodeficiency virus type 1 proteinase proteins from a precursor with a blocking mutation at the N terminus of the PR domain

The mature human immunodeficiency virus type 1 proteinase (PR; 11 kDa) can cleave all interdomain junctions in the Gag and Gag-Pol polyprotein precursors. To determine the activity of the enzyme in its precursor form, we blocked release of mature PR from a truncated Gag-Pol polyprotein by introducing mutations into the N-terminal Phe-Pro cleavage site of the PR domain. The mutant precursor autoprocessed efficiently upon expression in Escherichia coli. No detectable mature PR was released; however, several PR-related products ranging in size from approximately 14 to 18 kDa accumulated. Products of the same size were generated when mutant precursors were digested with wild-type PR. Thus, PR can utilize cleavage sites in the region upstream of the PR domain, resulting in the formation of extended PR species. On the basis of active-site titration, the PR species generated from mutated precursor exhibited wild-type activity on peptide substrates. However, the proteolytic activity of these extended enzymes on polyprotein substrates provided exogenously was low when equimolar amounts of extended and wild-type PR proteins were compared. Mammalian cells expressing the mutated precursor produced predominantly precursor and considerably reduced amounts of mature products. Released particles consisted mostly of uncleaved or partially cleaved polyproteins. Our results suggest that precursor forms of PR can autoprocess but are less efficient in processing of the Gag precursor for formation of mature virus particles.

Probing structure-function relationships in human immunodeficiency virus type 1 protease via molecular dynamics simulation

This chapter has focused on the application of molecular dynamics computer simulations and related molecular modeling techniques to the study of HIV protease structure and structure-function relationships. The abundance of crystallographic data provides ample experimental quantities (average structures, temperature factors, and hydrogen bond topography) to validate the computational techniques employed. Furthermore, these studies provide insight into the structure and functional energetics of HIV-1 protease that would be difficult or impossible to study experimentally. This chapter covers studies that investigate correlated motion between and within subunits of the protease, mutants of the protease that disrupt the tertiary structure and dimer formation, and studies of HIV-1 protease-inhibitor complexes that rationalize both the protonation state of the active site and the observed binding strength of these complexes. These studies demonstrate that MD is capable of contributing to our understanding of structure-function relationships and may aid in the design of potential therapeutics.

Motif-biased protein sequence alignment

A method was developed for pairwise protein sequence alignment to emulate the effect of structural knowledge or multiple sequences. Runs of matches of the preferred length were emphasized through the use of a product-bias allowing short motifs to influence the alignment to a degree that was a realistic reflection of their infrequency of occurrence. This gave motifs a locally high scoring match, making their alignment relatively less sensitive to the value of the gap penalty. This property should be a great advantage when a large number of sequence comparisons are made with a fixed set of parameter values, as typically occurs in the scan of a sequence databank with a probe or in the development of a multiple alignment.

Studies on the specificity of HIV protease: an application of Markov chain theory

A sequence-coupled (Markov chain) model is proposed to predict the cleavage sites in proteins by proteases with extended specificity subsites. In addition to the probability of an amino acid occurring at each of these subsites as observed from a training set of oligopeptides known cleavable by HIV protease, the conditional probabilities as reflected by the neighbor-coupled effect along the subsite sequence are also taken into account. These conditional probabilities are derived from an expanded training set consisting of sufficiently large peptide sequences generated by the Monte Carlo sampling process. Very high accuracy was obtained in predicting protein cleavage sites by both HIV-1 and HIV-2 proteases. The new method provides a rapid and accurate means for analyzing the specificity of HIV protease, and hence can be used to help find effective inhibitors of HIV protease as potential drugs against AIDS. The principle of this method can also be used to study the specificity of any multisubsite enzyme.

Expression of virus-encoded proteinases: functional and structural similarities with cellular enzymes

Many viruses express their genome, or part of their genome, initially as a polyprotein precursor that undergoes proteolytic processing. Molecular genetic analyses of viral gene expression have revealed that many of these processing events are mediated by virus-encoded proteinases. Biochemical activity studies and structural analyses of these viral enzymes reveal that they have remarkable similarities to cellular proteinases. However, the viral proteinases have evolved unique features that permit them to function in a cellular environment. In this article, the current status of plant and animal virus proteinases is described along with their role in the viral replication cycle. The reactions catalyzed by viral proteinases are not simple enzyme-substrate interactions; rather, the processing steps are highly regulated, are coordinated with other viral processes, and frequently involve the participation of other factors.

Knowledge-based model building of proteins: concepts and examples

We describe how to build protein models from structural templates. Methods to identify structural similarities between proteins in cases of significant, moderate to low, or virtually absent sequence similarity are discussed. The detection and evaluation of structural relationships is emphasized as a central aspect of protein modeling, distinct from the more technical aspects of model building. Computational techniques to generate and complement comparative protein models are also reviewed. Two examples, P-selectin and gp39, are presented to illustrate the derivation of protein model structures and their use in experimental studies.

Protein fold recognition

An important, yet seemingly unattainable, goal in structural molecular biology is to be able to predict the native three-dimensional structure of a protein entirely from its amino acid sequence. Prediction methods based on rigorous energy calculations have not yet been successful, and best results have been obtained from homology modelling and statistical secondary structure prediction. Homology modelling is limited to cases where significant sequence similarity is shared between a protein of known structure and the unknown. Secondary structure prediction methods are not only unreliable, but also do not offer any obvious route to the full tertiary structure. Recently, methods have been developed whereby entire protein folds are recognized from sequence, even where little or no sequence similarity is shared between the proteins under consideration. In this paper we review the current methods, including our own, and in particular offer a historical background to their development. In addition, we also discuss the future of these methods and outline the developments under investigation in our laboratory.