Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “HMM”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25Linked to original sources

The precursor of a metalloendopeptidase from human rheumatoid synovial fibroblasts. Purification and mechanisms of activation by endopeptidases and 4-aminophenylmercuric acetate.

Two active forms (Mr 45,000 and 28,000) of a metalloendopeptidase that digest proteoglycans and other extracellular matrix components of connective tissues have previously been purified from rheumatoid synovial cells and characterized [Okada, Nagase & Harris (1986) J. Biol. Chem. 261, 14245-14255]. To study the mechanisms of activation the precursor of this metalloendopeptidase has now been purified. The final products are homogeneous on SDS/polyacrylamide-gel electrophoresis and identified as a set of zymogens of Mr 57,000 and 59,000, in which the latter form is probably the product of post-translational glycosylation of the Mr 57,000 zymogen, as it binds to concanavalin A. The zymogen can be activated by trypsin, chymotrypsin, plasma kallikrein, plasmin and thermolysin, but not by thrombin. Although the activated metalloendopeptidase is further degraded by trypsin, plasma kallikrein and thermolysin during a prolonged incubation, it is relatively stable against plasmin and chymotrypsin. Activation with 4-aminophenylmercuric acetate is dependent on its concentration. It requires the reaction with the zymogen, possibly through thiol groups, and the continued presence of the agent. During this treatment the zymogen undergoes a sequential processing; first it becomes active without changing its apparent molecular mass, and then it is processed to low-Mr species of Mr 46,000, 45,000 (HMM) and 28,000 (LMM). The rate of conversion of the precursor into an initial intermediate of Mr 46,000 follows first-order kinetics (t1/2 2.0 h with 1.5 mM-4-amino-phenylmercuric acetate at 37 degrees C) and is independent of the initial concentration of the zymogen or the presence of up to a 676-fold molar excess of substrate, whereas the generation of HMM and LMM species is affected by these parameters. These results indicate that activation of the prometalloendopeptidase by an organomercurial compound is initiated by the molecular perturbation of the zymogen that results in conversion into the 46,000-Mr intermediate by an intramolecular action; the subsequent processing of this intermediate in HMM and LMM species is a bimolecular reaction. In vivo it is probable that the precursor of this metalloendopeptidase is activated either by direct limited proteolysis by tissue or plasma endopeptidases, or, alternatively, by factors that cause certain conformational changes in the zymogen molecule.

Arthritis, Rheumatoid↗

Genome-wide analysis of the general stress response in Bacillus subtilis.

Bacteria respond to diverse growth-limiting stresses by producing a large set of general stress proteins. In Bacillus subtilis and related Gram-positive pathogens, this response is governed by the sigma(B) transcription factor. To establish the range of cellular functions associated with the general stress response, we compared the transcriptional profiles of wild and mutant strains under conditions that induce sigma(B) activity. Macroarrays representing more than 3900 annotated reading frames of the B. subtilis genome were hybridized to (33)P-labelled cDNA populations derived from (i) wild-type and sigB mutant strains that had been subjected to ethanol stress; and (ii) a strain in which sigma(B) expression was controlled by an inducible promoter. On the basis of their significant sigma(B)-dependent expression in three independent experiments, we identified 127 genes as prime candidates for members of the sigma(B) regulon. Of these genes, 30 were known previously or inferred to be sigma(B) dependent by other means. To assist in the analysis of the 97 new genes, we constructed hidden Markov models (HMM) that identified possible sigma(B) recognition sequences preceding 21 of them. To test the HMM and to provide an independent validation of the hybridization experiments, we mapped the sigma(B)-dependent messages for seven representative genes. For all seven, the 5' end of the message lay near typical sigma(B) recognition sequences, and these had been predicted correctly by the HMM for five of the seven examples. Lastly, all 127 gene products were assigned to functional groups by considering their similarity to known proteins. Notably, products with a direct protective function were in the minority. Instead, the general stress response increased relative message levels for known or predicted regulatory proteins, for transporters controlling solute influx and efflux, including potential drug efflux pumps, and for products implicated in carbon metabolism, envelope function and macromolecular turnover.

Bacillus subtilis↗

Structural alterations of alpha-crystallin during its chaperone action.

The small heat-shock protein, alpha-crystallin, has chaperone ability whereby it stabilises proteins under stress conditions. In this study, alterations in the structure of alpha-crystallin during its interaction with a variety of substrate proteins (insulin, alpha-lactalbumin, ovotransferrin and serum albumin) under stress conditions have been examined using visible absorption, 31P-NMR and 1H-NMR and fluorescence spectroscopy. The fluorescence and 31P-NMR data imply that during the chaperone action of alpha-crystallin under reducing conditions, there is a slight increase in hydrophilicity of its N-terminal region and an alteration in flexibility of its C-terminal region, but overall, alpha-crystallin does not undergo a gross structural change. The fluorescence data suggest that substrate proteins interact with alpha-crystallin in a molten globule or intermediately folded state. The same conclusion is made from 1H-NMR spectroscopic monitoring of the interaction of alpha-crystallin with substrate proteins, e.g. the insulin B chain. The stoichiometry of interaction between alpha-crystallin and the various substrate proteins reveals that steric factors are important in determining the efficiency of interaction between the two proteins, i.e. on a molar subunit basis, alpha-crystallin is a more efficient chaperone protein with smaller substrate proteins. Comparison is also made between the high-molecular-mass (HMM) complexes formed between alpha-crystallin and ovotransferrin when reduced and heat stressed. Under heating conditions, fluorescence spectroscopy indicates that the HMM complex has a greater exposure of hydrophobicity to solution than that formed by reduction. Furthermore, in interacting with heated ovotransferrin, the C-terminal extension of the alphaB-crystallin subunit preferentially loses its flexibility suggesting that it is involved in stabilising bound ovotransferrin. By contrast, this extension is only partially reduced in flexibility in the HMM complex formed after reduction of ovotransferrin. The functional role of the C-terminal extensions in the chaperone action and the overall quaternary structure of alpha-crystallin is discussed.

Anilino Naphthalenesulfonates↗

Lipopolysaccharide O-chain microheterogeneity of Salmonella serotypes Enteritidis and Typhimurium.

Variability in the lipopolysaccharide (LPS) of the two most prevalent Salmonella serotypes causing food-borne salmonellosis was assessed using gas chromatography analysis of neutral sugars from 43 Salmonella enterica serovar Enteritidis (S. Enteritidis) and 20 Salmonella enterica serovar Typhimurium (S. Typhimurium) isolates. Four substantially different types of O-chain chemotypes were detected using cluster analysis of sugar compositions; these were low-molecular-mass (LMM) LPS, glucosylated LMM LPS, high-molecular-mass (HMM) LPS and glucosylated HMM LPS. Nineteen out of 20 S. Typhimurium isolates yielded glucosylated LMM. In contrast, S. Enteritidis produced a more diverse structure, which varied according to the source and history of the isolate: 45.5% of egg isolates yielded glucosylated HMM LPS; 100% of stored strains lacked glucosylation but retained chain length in some cases; and 83.3% of fresh isolates from the naturally infected house mouse Mus musculus produced glucosylated LMM LPS. A chain length determinant (wzz) mutant of S. Enteritidis produced a structure similar to that of S. Typhimurium and was used to define what constituted significant differences in structure using cluster analysis. Fine mapping of the S. Enteritidis chromosome by means of a two-restriction enzyme-ribotyping technique suggested that mouse isolates producing glucosylated LMM LPS were closely related to orally invasive strains obtained from eggs, and that stored strains were accumulating genetic changes that correlated with suppression of LPS O-chain glucosylation. These results suggest that the determination of LPS chemotype is a useful tool for epidemiological monitoring of S. Enteritidis, which displays an unusual degree of diversity in its LPS O-chain.

Animals↗

Genomewide identification of Pseudomonas syringae pv. tomato DC3000 promoters controlled by the HrpL alternative sigma factor.

The ability of Pseudomonas syringae pv. tomato DC3000 to parasitize tomato and Arabidopsis thaliana depends on genes activated by the HrpL alternative sigma factor. To support various functional genomic analyses of DC3000, and specifically, to identify genes involved in pathogenesis, we developed a draft sequence of DC3000 and used an iterative process involving computational and gene expression techniques to identify virulence-implicated genes downstream of HrpL-responsive promoters. Hypersensitive response and pathogenicity (Hrp) promoters are known to control genes encoding the Hrp (type III protein secretion) machinery and a few type III effector proteins in DC3000. This process involved (i) identification of 9 new virulence-implicated genes in the Hrp regulon by miniTn5gus mutagenesis, (ii) development of a hidden Markov model (HMM) trained with known and transposon-identified Hrp promoter sequences, (iii) HMM identification of promoters upstream of 12 additional virulence-implicated genes, and (iv) microarray and RNA blot analyses of the HrpL-dependent expression of a representative subset of these DC3000 genes. We found that the Hrp regulon encodes candidates for 4 additional type III secretion machinery accessory factors, homologs of the effector proteins HopPsyA, AvrPpiB1 (2 copies), AvrPpiC2, AvrPphD (2 copies), AvrPphE, AvrPphF, and AvrXv3, and genes associated with the production or metabolism of virulence factors unrelated to the Hrp type III secretion system, including syringomycin synthetase (SyrE), N(epsilon)-(indole-3-acetyl)-l-lysine synthetase (IaaL), and a subsidiary regulon controlling coronatine production. Additional candidate effector genes, hopPtoA2, hopPtoB2, and an avrRps4 homolog, were preceded by Hrp promoter-like sequences, but these had HMM expectation values of relatively low significance and were not detectably activated by HrpL.

Bacterial Proteins↗

Distinct metastable atmospheric regimes despite nearly Gaussian statistics: a paradigm model.

A controversial topic in the recent climate modeling literature is the fashion in which metastable low-frequency regimes in the atmosphere occur despite nearly Gaussian statistics for these planetary waves. Here a simple 57-mode paradigm model for such metastable atmospheric regime behavior is introduced and analyzed through hidden Markov model (HMM) analysis of the time series of suitable low-frequency planetary waves. The analysis of this paradigm model elucidates how statistically significant metastable regime transitions between blocked and zonal statistical states occur despite nearly Gaussian behavior in the associated probability distribution function and without a significant role for the low-order truncated nonlinear dynamics alone; turbulent backscatter onto the three-dimensional subspace of low-frequency modes is responsible for these effects. It also is demonstrated that suitable stochastic mode reduction strategies, which include both augmented cubic nonlinearity and multiplicative noise, are also capable of capturing the metastable low-frequency regime behavior through a single stochastic differential equation compared with the full turbulent chaotic 57-mode model. This feature is attractive for issues such as long-term weather predictability. Although there have been many applications of HMM in other sciences, this work presents a previously undescribed application of HMM analysis to atmospheric low-frequency variability and points the way for further applications including their use in extended range predictability.

Journal Article↗

Myosin II folding is mediated by a molecular chaperonin.

The folding pathway of the heavy meromyosin subfragment (HMM) of a skeletal muscle myosin has been investigated by in vitro synthesis of the myosin heavy and light chains in a coupled transcription and translation assay. Analysis of the nascent translation products for folding intermediates has identified a major intermediate that contains all three myosin subunits in a complex with the eukaryotic cytosolic chaperonin. Partially folded HMM is released from this complex in an ATP-dependent manner. However, biochemical and functional assays reveal incomplete folding of the myosin motor domain. Dimerization of myosin heavy chains and association of heavy and light chains are accomplished early in the folding pathway. To test for other factors necessary for the complete folding of myosin, a cytoplasmic extract was prepared from myotubes produced by a mouse myogenic cell line. This extract dramatically enhanced the folding of HMM, suggesting a role for muscle-specific factors in the folding pathway. We conclude that the molecular assembly of myosin is mediated by the eukaryotic cytosolic chaperonin with folding of the motor domain as the slow step in the pathway.

Adenosine Triphosphate↗

ADP binding induces an asymmetry between the heads of unphosphorylated myosin.

Light chain phosphorylation is the key event that regulates smooth and non-muscle myosin II ATPase activity. Here we show that both heads of smooth muscle heavy meromyosin (HMM) bind tightly to actin in the absence of nucleotide, irrespective of the state of light chain phosphorylation. In striking contrast, only one of the two heads of unphosphorylated HMM binds to actin in the presence of ADP, and the heads have different affinities for ADP. This asymmetry suggests that phosphorylation alters the mechanical coupling between the heads of HMM. A model that incorporates strain between the two heads is proposed to explain the data, which have implications for how one head of a motor protein can gate the response of the other.

Actins↗

H(2)O(2)-induced kinetic and chemical modifications of smooth muscle myosin: correlation to effects of H(2)O(2) on airway smooth muscle.

The effect of H(2)O(2) on smooth muscle heavy meromyosin (HMM) and subfragment 1 (S1) was examined. The number of molecules that retained the ability to bind ATP and the actinactivated rate of P(i) release were measured by single-turnover kinetics. H(2)O(2) treatment caused a decrease in HMM regulation from 800- to 27-fold. For unphosphorylated and phosphorylated heavy meromyosin and for S1, approximately 50% of the molecules lost the ability to bind to ATP. H(2)O(2) treatment in the presence of EDTA protected against ATPase inactivation and against the loss of total ATP binding. Inactivation of S1 versus time correlated to a loss of reactive thiols. Treatment of H(2)O(2)-inactivated phosphorylated HMM or S1 with dithiothreitol partially reactivated the ATPase but had no effect on total ATP binding. H(2)O(2)-inactivated S1 contained a prominent cross-link between the N-terminal 65-kDa and C-terminal 26-kDa heavy chain regions. Mass spectral studies revealed that at least seven thiols in the heavy chain and the essential light chain were oxidized to cysteic acid. In thiophosphorylated porcine tracheal muscle strips at pCa 9 + 2.1 mM ATP, H(2)O(2) caused a approximately 50% decrease in the amplitude but did not alter the rate of force generation, suggesting that H(2)O(2) directly affects the force generating complex. Dithiothreitol treatment reversed the H(2)O(2) inhibition of the maximal force by approximately 50%. These data, when compared with the in vitro kinetic data, are consistent with a H(2)O(2)-induced loss of functional myosin heads in the muscle.

Actins↗

Effects of actin filament cross-linking and filament length on actin-myosin interaction.

We have used two actin-binding proteins of the intestinal brush border, TW 260/240 and villin, to examine the effects of filament cross-linking and filament length on myosin-actin interactions. TW 260/240 is a nonerythroid spectrin that is a potent cross-linker of actin filaments. In the presence of this cross-linker we observed a concentration-dependent enhancement of skeletal muscle actomyosin ATPase activity (150-560% of control; maximum enhancement at a 1:70-80 TW 260/240:actin molar ratio). TW 260/240 did not cause a similar enhancement of either acto-heavy meromyosin (HMM) ATPase or acto-myosin subfragment-one (S1) ATPase. Villin, a Ca2+-dependent filament capping and severing protein of the intestinal microvillus, was used to generate populations of actin filaments of various lengths from less than 20 nm to 2.0 microns; (villin:actin ratios of 1:2 to 1:4,000). The effect of filament length on actomyosin ATPase was biphasic. At villin:actin molar ratios of 1:2-25 actin-activated myosin ATPase activity was inhibited to 20-80% of control values, with maximum inhibition observed at the highest villin:actin ratio. The ATPase activities of acto-HMM and acto-S1 were also inhibited at these short filament lengths. At intermediate filament lengths generated at villin:actin ratios of 1:40-400 (average lengths 0.26-1.1 micron) an enhancement of actomyosin ATPase was observed (130-260% of controls), with a maximum enhancement at average filament lengths of 0.5 micron. The levels of actomyosin ATPase fell off to control values at low concentrations of villin where filament length distributions were almost those of controls. Unlike intact myosin, the actin-activated ATPase of neither HMM nor S1 showed an enhancement at these intermediate actin filament lengths.

Actin Cytoskeleton↗

A mutant heterodimeric myosin with one inactive head generates maximal displacement.

Each of the heads of the motor protein myosin II is capable of supporting motion. A previous report showed that double-headed myosin generates twice the displacement of single-headed myosin (Tyska, M.J., D.E. Dupuis, W.H. Guilford, J.B. Patlak, G.S. Waller, K.M. Trybus, D.M. Warshaw, and S. Lowey. 1999. Proc. Natl. Acad. Sci. USA. 96:4402-4407). To determine the role of the second head, we expressed a smooth muscle heterodimeric heavy meromyosin (HMM) with one wild-type head, and the other locked in a weak actin-binding state by introducing a point mutation in switch II (E470A). Homodimeric E470A HMM did not support in vitro motility, and only slowly hydrolyzed MgATP. Optical trap measurements revealed that the heterodimer generated unitary displacements of 10.4 nm, strikingly similar to wild-type HMM (10.2 nm) and approximately twice that of single-headed subfragment-1 (4.4 nm). These data show that a double-headed molecule can achieve a working stroke of approximately 10 nm with only one active head and an inactive weak-binding partner. We propose that the second head optimizes the orientation and/or stabilizes the structure of the motion-generating head, thereby resulting in maximum displacement.

Actins↗

Studies of muscle proteins in embryonic myocardial cells of cardiac lethal mutant mexican axolotls (Ambystoma mexicanum) by use of heavy meromyosin binding and sodium dodecyl sulfate polyacrylamide gel electrophoresis.

In the Mexican axolotl Ambystoma mexicanum recessive mutant gene c, by way of abnormal inductive processes from surrounding tissues, results in an absence of embryonic heart function. The lack of contractions in mutant heart cells apparently results from their inability to form normally organized myofibrils, even though a few actin-like (60-A) and myosin-like (150-A) filaments are present. Amorphous "proteinaceous" collections are often visible. In the present study, heavy meromyosin (HMM) treatment of mutant heart tissue greatly increases the number of thin filaments and decorates them in the usual fashion, confirming that they are actin. The amorphous collections disappear with the addition of HMM. In addition, an analysis of the constituent proteins of normal and mutant embryonic hearts and other tissues is made by sodium dodecyl sulfate (SDS) gel electrophoresis. These experiments are in full agreement with the morphological and HMM binding studies. The gels show distinct 42,000-dalton bands for both normal and mutant hearts, supporting the presence of normal actin. During early developmental stages (Harrison's stage 34) the cardiac tissues in normal and mutant siblings have indistinguishable banding patterns, but with increasing development several differences appear. Myosin heavy chain (200,000 daltons) increases substantially in normal hearts during development but very little in mutants. Even so the quantity of 200,000-dalton protein in mutant hearts is significantly more than in any of the nonmuscle tissues studied (i.e. gut, liver, brain). Unlike normal hearts, the mutant hearts lack a prominent 34,000-dalton band, indicating that if mutants contain muscle tropomyosin at all, it is present in drastically reduced amounts. Also, mutant hearts retain large amounts of yolk proteins at stages when the platelets have virtually disappeared from normal hearts. The morphologies and electrophoresis patterns of skeletal muscle from normal and mutant siblings are identical, confirming that gene c affects only heart muscle differentiation and not skeletal muscle. The results of the study suggest that the precardiac mesoderm in cardiac lethal mutant axolotl embryos initiates but then fails to complete its differentiation into functional muscle tissue. It appears that this single gene mutation, by way of abnormal inductive processes, affects the accumulation and organization of several different muscle proteins, including actin, myosin, and tropomyosin.

Actins↗

Heavy meromyosin labeling of intermediate filaments in cultured connective tissue cells.

Mild treatment with trypsin causes a radical change in the heavy meromyosin (HMM) binding properties of intermediate filaments in glycerinated, myosin-extracted cultured chick embryo connective tissue cells. In non-trypsin-treated cells, HMM labeling of filaments was often indistinct and variable in its distribution. By contrast, in cells treated with trypsin (under conditions which allowed most intermediate filaments to survive), virtually all filaments, including those of intermediate size, decorated with HMM to give distinct arrowhead patterns. We suggest that most intermediate filaments in such cells contain a core of F-actin masked by trypsin-labile accessory proteins.

Actins↗

Banding and polarity of actin filaments in interphase and cleaving cells.

Heavy meromyosin (HMM) decoration of actin filaments was used to detect the polarity of microfilaments in interphase and cleaving rat kangaroo (PtK2) cells. Ethanol at -20 degrees C was used to make the cells permeable to HMM followed by tannic acid-glutaraldehyde fixation for electron microscopy. Uniform polarity of actin filaments was observed at cell junctions and central attachment plaques with the HMM arrowheads always pointing away from the junction or plaque. Stress fibers were banded in appearance with their component microfilaments exhibiting both parallel and antiparallel orientation with respect to one another. Identical banding of microfilament bundles was also seen in cleavage furrows with the same variation in filament polarity as found in stress fibers. Similarly banded fibers were not seen outside the cleavage furrow in mitotic cells. By the time that a mid-body was present, the actin filaments in the cleavage furrow were no longer in banded fibers. The alternating dark and light bands of both the stress fibers and cleavage furrow fibers are approximately equal in length, each measuring approximately 0.16 micrometer. Actin filaments were present in both bands, and individual decorated filaments could sometimes be traced through four band lengths. Undecorated filaments, 10 nm in diameter, could often be seen within the light bands. A model is proposed to explain the arrangement of filaments in stress fibers and cleavage furrows based on the striations observed with tannic acid and the polarity of the actin filaments.

Actins↗

Qualitative analysis of skeletal myosin as substrate of Ca2+-activated neutral protease: comparison of filamentous and soluble, native, and L2-deficient myosin.

Ca2+ -activated neutral protease (CAF) was capable of degrading myosin over a 200-fold range of protease concentrations. CAF selected the heavy chain of myosin, although either prolonged exposure to or high concentrations of the protease degraded the L1, but not the L2 or L3, light chains of myosin. The following results indicated that during the first hour of digestion, under conditions where native myosin was the substrate, CAF selected for the "head" region of the myosin heavy chain: (a) large heavy chain fragments of identical molecular weight were produced from filamentous and from soluble myosin; (b) light meromyosin was not a substrate; (c) agents known to bind to the head of myosin (actin, MgATP, and L2) had both a qualitative and quantitative effect on degradation; and (d) similar cleavage sites could be demonstrated for myosin and for heavy meromyosin (HMM) despite the fact that HMM was a much poorer substrate than myosin. This observation is interpreted as an indication that the conformation of myosin heavy chain is altered in the preparation of HMM. The principal cleavage sites on the heavy chain of myosin were 20,000, 35,000 and 50,000 D from the N-terminus, producing large fragments with molecular weights of 180,000, 165,000, and 150,000 which comprised a "nicked" species of myosin. This nicked species retained both normal solubility properties and normal hydrolytic activities. For this reason, it is concluded that "nicked myosin" is an important pathophysiological species.

Animals↗

transition priors for protein hidden Markov models: an empirical study towards maximum discrimination.

Insertions and deletions in a profile hidden Markov model (HMM) are modeled by transition probabilities between insert, delete and match states. These are estimated by combining observed data and prior probabilities. The transition prior probabilities can be defined either ad hoc or by maximum likelihood (ML) estimation. We show that the choice of transition prior greatly affects the HMM's ability to discriminate between true and false hits. HMM discrimination was measured using the HMMER 2.2 package applied to 373 families from Pfam. We measured the discrimination between true members and noise sequences employing various ML transition priors and also systematically scanned the parameter space of ad hoc transition priors. Our results indicate that ML priors produce far from optimal discrimination, and we present an empirically derived prior that considerably decreases the number of misclassifications compared to ML. Most of the difference stems from the probabilities for exiting a delete state. The ML prior, which is unaware of noise sequences, estimates a delete-to-delete probability that is relatively high and does not penalize noise sequences enough for optimal discrimination.

Algorithms↗

Hidden Markov models for detecting remote protein homologies.

MOTIVATION: A new hidden Markov model method (SAM-T98) for finding remote homologs of protein sequences is described and evaluated. The method begins with a single target sequence and iteratively builds a hidden Markov model (HMM) from the sequence and homologs found using the HMM for database search. SAM-T98 is also used to construct model libraries automatically from sequences in structural databases. METHODS: We evaluate the SAM-T98 method with four datasets. Three of the test sets are fold-recognition tests, where the correct answers are determined by structural similarity. The fourth uses a curated database. The method is compared against WU-BLASTP and against DOUBLE-BLAST, a two-step method similar to ISS, but using BLAST instead of FASTA. RESULTS: SAM-T98 had the fewest errors in all tests-dramatically so for the fold-recognition tests. At the minimum-error point on the SCOP (Structural Classification of Proteins)-domains test, SAM-T98 got 880 true positives and 68 false positives, DOUBLE-BLAST got 533 true positives with 71 false positives, and WU-BLASTP got 353 true positives with 24 false positives. The method is optimized to recognize superfamilies, and would require parameter adjustment to be used to find family or fold relationships. One key to the performance of the HMM method is a new score-normalization technique that compares the score to the score with a reversed model rather than to a uniform null model. AVAILABILITY: A World Wide Web server, as well as information on obtaining the Sequence Alignment and Modeling (SAM) software suite, can be found at http://www.cse.ucsc.edu/research/compbi o/ CONTACT: karplus@cse.ucsc.edu; http://www.cse.ucsc.edu/karplus

Algorithms↗

Profile hidden Markov models.

The recent literature on profile hidden Markov model (profile HMM) methods and software is reviewed. Profile HMMs turn a multiple sequence alignment into a position-specific scoring system suitable for searching databases for remotely homologous sequences. Profile HMM analyses complement standard pairwise comparison methods for large-scale sequence analysis. Several software implementations and two large libraries of profile HMMs of common protein domains are available. HMM methods performed comparably to threading methods in the CASP2 structure prediction exercise.

Humans↗