SEARCH · PubMed Health

Results for “alignment chaining method”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 937 records · Page 52Linked to original sources

PubMed2002

Transitional cell metaplasia of the uterine cervix is related to human papillomavirus: molecular analysis in seven patients with cytohistologic correlation.

BACKGROUND: Transitional cell metaplasia of the uterine cervix is an under-recognized entity in cervical pathology. The underlying etiology and biologic significance remains uncertain. The thin-layer cytology findings and association with human papillomavirus (HPV) have not been studied thoroughly. METHODS: The authors retrospectively reviewed the clinical findings, thin-layer cytology and histologic features of pure transitional cell metaplasia of the uterine cervix occurring in seven perimenopausal or postmenopausal Chinese women at Pamela Youde Nethersole Eastern Hospital, Hong Kong, during the period from January, 1998 to April, 2001. Molecular techniques for HPV screening and genotyping using polymerase chain reaction and restriction fragment length polymorphism analysis were employed in the thin-layer cytology specimens and paraffin block material. RESULTS: In all seven patients, transitional cell metaplasia represented an incidental histologic finding. It occurred in the ectocervix, transformation zone, endocervix, or vagina. Histologically, it resembled urothelium of the urinary bladder and was comprised of multilayers of mitotically inactive, immature epithelial cells with vertically aligned oval nuclei, fine chromatin, indistinct nucleoli, and conspicuous longitudinal nuclear grooves. The superficial cells were oriented more horizontally and contained pale-staining cytoplasm similar to umbrella cells. Features consistent with transitional cell metaplasia were identified in two of seven preoperative thin-layer preparations. Cytologically, the affected parabasal cells recapitulated the features that were seen in histologic sections. In addition to the bland nuclear morphology and longitudinal nuclear grooves, the cell borders appeared distinct, and the appearance of a perinuclear cytoplasmic halo was common. Sometimes, the metaplastic cells assumed a spindle shape and appeared as cohesive, streaming cell clusters. Molecular study successfully demonstrated the presence of HPV in all seven patients, mostly in the liquid-based cytology samples. In general, the viral DNA load was relatively low; and, for samples in which HPV genotyping was feasible, HPV type 58 was the prevalent genotype. CONCLUSIONS: The current study demonstrates that transitional cell metaplasia of the uterine cervix is related to HPV. It also carries a distinctive cytologic appearance in thin-layer preparations. Based on the limited follow-up data from a small number of reported patients, transitional cell metaplasia seems to run an indolent clinical course. However, its peculiar association with HPV and its possible correlation, both morphologic and histogenetic, with cervical intraepithelial neoplasia need further investigation.

Aged↗

PubMedDate not supplied

Evaluation measures of multiple sequence alignments.

Multiple sequence alignments (MSAs) are frequently used in the study of families of protein sequences or DNA/RNA sequences. They are a fundamental tool for the understanding of the structure, functionality and, ultimately, the evolution of proteins. A new algorithm, the Circular Sum (CS) method, is presented for formally evaluating the quality of an MSA. It is based on the use of a solution to the Traveling Salesman Problem, which identifies a circular tour through an evolutionary tree connecting the sequences in a protein family. With this approach, the calculation of an evolutionary tree and the errors that it would introduce can be avoided altogether. The algorithm gives an upper bound, the best score that can possibly be achieved by any MSA for a given set of protein sequences. Alternatively, if presented with a specific MSA, the algorithm provides a formal score for the MSA, which serves as an absolute measure of the quality of the MSA. The CS measure yields a direct connection between an MSA and the associated evolutionary tree. The measure can be used as a tool for evaluating different methods for producing MSAs. A brief example of the last application is provided. Because it weights all evolutionary events on a tree identically, but does not require the reconstruction of a tree, the CS algorithm has advantages over the frequently used sum-of-pairs measures for scoring MSAs, which weight some evolutionary events more strongly than others. Compared to other weighted sum-of-pairs measures, it has the advantage that no evolutionary tree must be constructed, because we can find a circular tour without knowing the tree.

Algorithms↗

PubMed1993

Phylogeny of twenty Thermus isolates constructed from 16S rRNA gene sequence data.

The sequences of the 16S rRNA genes of 20 Thermus isolates were determined to a high fidelity by using automated DNA sequencing and fluorescent-dye-labelled primers. The strains tested included members of the three validly named Thermus species and representatives of major taxonomic clusters defined previously for this genus. The parsimony method was used to reconstruct the phylogeny of the strains from the aligned sequences, and a bootstrap analysis revealed a number of well-supported clades. Our results are not consistent with groupings inferred from numerical taxonomy data but support the conjecture that the genus Thermus contains more species than the three currently recognized species.

Base Sequence↗

PubMed2002

Three-dimensional quantitative structure-activity relationship of 1,4-dihydropyridines as antitubercular agents.

Three-dimensional quantitative structure-activity relationship (3D QSAR) methods, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), were applied on a series of 1,4-dihydropyridines possessing antitubercular activity. The study was performed using 33 compounds, in which 22 molecules were used for the derivation of the 3D QSAR models (training set) and 11 molecules were used to evaluate the predictive ability of the derived models (test set). Superimpositions were performed using three alignment rules: atom-based fitting, SYBYL QSAR rigid body field fit of the steric and electrostatic fields of the molecules, and flexible fitting (multifit). Both methods were analyzed in terms of their predictive abilities and produced comparable results with high internal as well as external predictivities. Steric and electrostatic fields of the inhibitors were found to be relevant descriptors for SAR. Use of lowest unoccupied molecular orbital energies or ClogP as additional descriptors in the QSAR table did not improve the significance of the 3D QSAR models. Both CoMFA and CoMSIA models based on multifit alignment showed better correlative and predictive properties than other models. A QSAR study using genetic function approximation was also performed for the same set of molecules using different types of physicochemical descriptors to deal with cell-based activity data. The QSAR models revealed the importance of spatial properties and conformational flexibility of side chains for antitubercular activity. Inclusion of fractional polar solvent accessible surface area as a descriptor in the model generation resulted in models with significant internal and external predictivities for the same test set molecules, which may support the possible mode of action of these compounds.

Antitubercular Agents↗

PubMed1999

The cDNA-derived amino acid sequence of hemoglobin I from Lucina pectinata.

The tropical clam Lucina pectinata contains a unique hemoglobin (HbI) which serves to transport H2S to autotrophic bacteria. The cDNA-derived amino acid sequence was obtained from overlapping clones containing the cDNA that codes for HbI. The reverse transcriptase-polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE) methods were employed to synthesize the cDNA fragments. An initial 354-bp cDNA clone encoding 118 amino acid residues of HbI was amplified from total RNA by RT-PCR using degenerate oligonucleotides. Gene-specific primers derived from the HbI-partial cDNA sequence were used for obtaining the 5' and 3' ends of the cDNA by RACE. The length of the HbI cDNA, estimated from sequence analysis of overlapping clones, was 1322 bp for the full-length cDNA. The coding region of the full-length cDNA codes for 143 amino acid residues. The most conserved amino acid residues in HbI from Lucina pectinata were identified by a multiple alignment with nonvertebrate globin sequences.

Amino Acid Sequence↗

PubMed1998

Cloning of a full-length insulin-like growth factor-I complementary DNA in the goldfish liver and ovary and development of a quantitative PCR method for its measurement.

Five forms of insulin-like growth factor-I (IGF-I) complementary DNA (cDNA) were isolated by PCR from goldfish liver and ovary, using primers based on common carp IGF-I sequence. In the goldfish liver, we cloned and sequenced three IGF-I forms (1,2, and 3), and elucidated the full-length cDNA sequence using the 5'-and 3'-RACE. Two IGF-I forms (1 and 2) were cloned from the goldfish ovary and were found to have differences with respect to both size and nucleotide sequence compared to liver IGF-I. The entire liver IGF-I form 1 sequence was found to be 833 nucleotides long, containing a 483-nucleotide open reading from encoding 161 amino acids. The deduced amino acid sequence of the mature peptide was compared to IGF-I sequences of other vertebrates, and found to have 97 and 93% similarity to carp and salmon IGF-I, respectively, In this study we also developed a competitive quantitative PCR method and demonstrated an increase in IGF-I expression following treatments with growth hormone or gonadotropin-releasing hormone in the goldfish liver.

Amino Acid Sequence↗

PubMed2004

A sensitive and quantitative single-tube real-time reverse transcriptase-PCR for detection of enteroviral RNA.

BACKGROUND: Enteroviruses (EVs) are significant human pathogens. Rapid and sensitive diagnostic techniques are desirable. OBJECTIVES: To develop a quantitative single-tube real-time reverse transcription-polymerase chain reaction (RT-PCR) for human enterovirus ribonucleic acid (RNA) (QPCR), with protection against amplimer contamination. STUDY DESIGN: The method was evaluated with serial dilutions of EV, 62 cerebrospinal fluid (CSF) specimens from meningitis patients, and the third and fourth European Union Concerted Action Enterovirus Proficiency Panels. A commercial EV PCR test was run in parallel. RESULTS: Optimisation included RNA extraction procedure, design and concentrations of primers and probes from the 5' non-coding region as well as recombinant Thermus thermophilus polymerase (rTth), Mn(OAc)(2) and thermolabile UNG concentrations. Of 62 CSF samples from cases of meningitis submitted for QPCR testing, 34 (76%) and 21 (47%) were positive by QPCR and a commercial EV RNA detection kit, respectively. The detection limit of QPCR was 0.001 TCID(50)/ml (50% tissue culture-infective dose per millilitre) for a coxsackievirus B2 preparation and <10 copies of a plasmid containing coxsackievirus B2 complementary deoxyribonucleic acid (cDNA). The relation between threshold cycle (C(t)) and amount of virus was linear (r = 0.99) over a range of 10(-3) to 10(4) TCID(50)/ml of coxsackievirus B2. CONCLUSIONS: The QPCR method allows a large number of samples to be screened rapidly. Its sensitivity, simplicity, and reproducibility make it a suitable tool for the routine laboratory.

5' Untranslated Regions↗

PubMed1998

Kinked-helices model of the nicotinic acetylcholine receptor ion channel and its complexes with blockers: simulation by the Monte Carlo minimization method.

A model of the nicotinic acetylcholine receptor ion channel was elaborated based on the data from electron microscopy, affinity labeling, cysteine scanning, mutagenesis studies, and channel blockade. A restrained Monte Carlo minimization method was used for the calculations. Five identical M2 segments (the sequence EKMTLSISVL10LALTVFLLVI20V) were arranged in five-helix bundles with various geometrical profiles of the pore. For each bundle, energy profiles for chlorpromazine, QX-222, pentamethonium, and other blocking drugs pulled through the pore were calculated. An optimal model obtained allows all of the blockers free access to the pore, but retards them at the rings of residues known to contribute to the corresponding binding sites. In this model, M2 helices are necessarily kinked. They come into contact with each other at the cytoplasmic end but diverge at the synaptic end, where N-termini of M1 segments may contribute to the pore. The kinks disengage alpha-helical H-bonds between Ala12 and Ser8. The uncoupled lone electron pairs of Ser8 carbonyl oxygens protrude into the pore, forming a hydrophilic ring that may be important for the permeation of cations. A split network of H-bonds provides a flexibility to the chains Val9-Ala12, the numerous conformations of which form only two or three intrasegment H-bonds. The cross-ectional dimensions of the interface between the flexible chains vary essentially at the level of Leu11. We suggest that conformational transitions in the chains Val9-Ala12 are responsible for the channel gating, whereas rotations of more stable alpha-helical parts of M2 segments may be necessary to transfer the channel in the desensitized state.

Amino Acid Sequence↗

PubMed2002

A DNA-based method for identification of krill species and its application to analysing the diet of marine vertebrate predators.

Accurate identification of species that are consumed by vertebrate predators is necessary for understanding marine food webs. Morphological methods for identifying prey components after consumption often fail to make accurate identifications of invertebrates because prey morphology becomes damaged during capture, ingestion and digestion. Another disadvantage of morphological methods for prey identification is that they often involve sampling procedures that are disruptive for the predator, such as stomach flushing or lethal collection. We have developed a DNA-based method for identifying species of krill (Crustacea: Malacostraca), an enormously abundant group of invertebrates that are directly consumed by many groups of marine vertebrates. The DNA-based approach allows identification of krill species present in samples of vertebrate stomach contents, vomit, and, more importantly, faeces. Utilizing samples of faeces from vertebrate predators minimizes the impact of dietary studies on the subject animals. We demonstrate our method first on samples of Adelie penguin (Pygoscelis adeliae) stomach contents, where DNA-based species identification can be confirmed by prey morphology. We then apply the method to faeces of Adelie penguins and to faeces of the endangered pygmy blue whale (Balaenoptera musculus brevicauda). In each of these cases, krill species consumed by the predators could be identified from their DNA present in faeces or stomach contents.

Animals↗

PubMed2005

Hidden Markov Models, grammars, and biology: a tutorial.

Biological sequences and structures have been modelled using various machine learning techniques and abstract mathematical concepts. This article surveys methods using Hidden Markov Model and functional grammars for this purpose. We provide a formal introduction to Hidden Markov Model and grammars, stressing on a comprehensive mathematical description of the methods and their natural continuity. The basic algorithms and their application to analyzing biological sequences and modelling structures of bio-molecules like proteins and nucleic acids are discussed. A comparison of the different approaches is discussed, and possible areas of work and problems are highlighted. Related databases and softwares, available on the internet, are also mentioned.

Algorithms↗

PubMed1997

The life cycle of Ceratomyxa shasta, a myxosporean parasite of salmonids, requires a freshwater polychaete as an alternate host.

The actinosporean life stage of Ceratomyxa shasta, a myxozoan parasite of salmonids, and the annelid worm that serves as its alternate host were identified in laboratory transmission experiments and their roles were confirmed using molecular techniques. Infection by the parasite occurred in susceptible fish that were either exposed to or force fed the freshwater polychaete, Manayunkia speciosa, infected with the actinosporean. These observations were confirmed using the polymerase chain reaction with primers designed from the C. shasta 18S rDNA sequence. DNA was amplified from polychaetes harboring the actinosporean that caused infection in the fish but not from uninfected polychaetes. Amplified DNA from an infected polychaete was sequenced and its homology with the 18S rDNA sequence of C. shasta spores verified the proposed life cycle. Ultrastructural examination of the actinosporean in the polychaete showed developmental stages in the epidermis rather than within the intestinal epithelium as described for other myxozoans. The methods described will be useful in identifying alternate hosts and morphologically diverse life stages in the complex life cycles of other myxosporea and in understanding the relationships between these parasites and their hosts.

Animals↗

PubMed2000

Species identification by means of the cytochrome b gene.

Species identification was carried out by nucleotide sequence analysis of the cytochrome b (cytb) gene. The aim of the study was to identify biological specimens from diverse vertebrate animals by extracting and amplifying DNA from 44 different animal species covering the 5 major vertebrate groups (i.e. mammals, birds, reptiles, amphibians and fishes). The sequences derived were used to identify the biological origin of the samples by aligning to cytb gene sequence entries in nucleotide databases using the program BLAST. All sequences were submitted to the GenBank including new species which were not observed in the databases. The applicability of this method to the forensic field is demonstrated by simulated casework conditions where different types of samples including problematic specimens such as hair, bone samples, bristles and feathers were investigated to identify the species.

Animals↗

PubMed1999

Controlled syntheses of natural and disulfide-mispaired regioisomers of alpha-conotoxin SI.

Methods are reported for the unambiguous syntheses of all three possible disulfide regioisomers with the sequence of alpha-conotoxin SI, a tridecapeptide amide from marine cone snail venom that binds selectively to the muscle subtype of nicotinic acetylcholine receptors. The naturally occurring peptide has two 'interlocking' disulfide bridges connecting Cys2-Cys7 and Cys3-Cys13 (2/7&3/13), while in the two mispaired isomers the disulfide bridges connect Cys2-Cys13 and Cys3-Cys7 (2/13 & 3/7, 'nested') and Cys2-Cys3 and Cys7-Cys13 (2/3 & 7/13, 'discrete'), respectively. Alignment of disulfide bridges was controlled at the level of orthogonal protection schemes for the linear precursors, assembled by Fmoc solid-phase peptide synthesis on acidolyzable tris(alkoxy)benzylamide (PAL) supports. Side-chain protection of cysteine was provided by suitable pairwise combination of the S-9H-xanthen-9-yl (Xan) and S-acetamidomethyl (Acm) protecting groups. The first disulfide bridge was formed from the corresponding bis(thiol) precursor obtained by selective deprotection of S-Xan, and the second disulfide bridge was formed by orthogonal co-oxidation of S-Acm groups on the remaining two Cys residues. It was possible to achieve the desired alignments with either order of loop formation (smaller loop before larger, or vice versa). The highest overall yields were obtained when both disulfides were formed in solution, while experiments where either the first or both bridges were formed while the peptide was on the solid support revealed lower overall yields and poorer selectivities towards the desired isomers.

Amino Acid Sequence↗

PubMed1982

Glutathione reductase from human erythrocytes. The sequences of the NADPH domain and of the interface domain.

1. Sequence analysis of the NADPH domain (residues 158--293) and of the interface domain (365--478) was based on 12 CNBr fragments, which were isolated using ion-exchange chromatography and paper methods. Fragments with more than 15 residues were digested further with trypsin and chymotrypsin. The isolated peptides were sequenced by automated solid-phase Edman degradation. All sequenced peptides were ordered and overlapped by computerized comparisons with a complete sequence guessed from the electron density map of the protein. In the case of short CNBr fragments, this alignment was confirmed by the sequence analysis of protein fragments resulting from incomplete CNBr cleavage. 2. In the NADPH domain, residue 197, which is involved in an induced-fit mechanism, was identified as a tyrosine. The structure of the NADPH domain is probably homologous with the NAD domain of lipoamide dehydrogenase and with the FAD domain of several proteins, but not with NADPH domains of known chain-fold in other proteins. 3. The paper completes the sequence analysis of glutathione reductase so that the enzyme is now known in atomic detail. The numbering scheme of the chemically determined sequence will be used henceforth in crystallographic studies also. As inferred from the sequence data each of the two identical chains contains 478 amino acid residues, the composition being Cys10, Asp21, Asn17, Thr31, Ser31, Glu29, Gln11, Pro24, Gly43, Ala42, Val44, Met15, Ile29, Leu34, Tyr13, Phe14, Lys34, His16. Arg17, and Trp3. From these data an Mr of 2 x 51 600 was calculated for the FAD-free apoenzyme and an Mr of 2 x 42 400 for the holoenzyme.

Amino Acid Sequence↗

PubMed1992

A simple and rapid method for HLA-DQA1 genotyping by polymerase chain reaction-single strand conformation polymorphism and restriction enzyme cleavage analysis.

A simple and rapid method for identification of alleles at the human leucocyte antigen (HLA)-DQA1 locus is described. The polymorphic second exon of the HLA-DQA1 locus was amplified by the polymerase chain reaction (PCR) method. The amplified DNA was analyzed by single-strand conformation polymorphism (SSCP) and restriction enzyme cleavage assay. Using this method, the eight known DQA1 alleles could be distinguished from each other. This paper suggests that the method can be used for quick genotyping of DQA1 alleles, but detecting point mutations at various positions in a fragment as well as new HLA-DQA1 genotypes should also be possible.

Alleles↗

PubMed1996

Tubulin secondary structure analysis, limited proteolysis sites, and homology to FtsZ.

The far-ultraviolet circular dichroism spectrum of the alpha beta-tubulin dimer analyzed by six different methods indicates an average content of approximately 33% alpha helix, 21% beta sheet, and 45% other secondary structure. Deconvolution of Fourier transform infrared spectra indicates 24% sheet, 37% (maximum) helix, and 38% (minimum) other structure. Separate alignments of 75 alpha-tubulin, 106 beta-tubulin, and 14 gamma-tubulin sequences and 12 sequences of the bacterial cell division protein FtsZ have been employed to predict their secondary structures with the multiple-sequence method PHD [Rost, B., & Sander, C. (1993a) J. Mol. Biol. 232, 584-599]. The predicted secondary structures average of 33% alpha helix, 24% beta sheet, and 43% loop for the alpha beta dimer. The predictions have been compared with sites of limited proteolysis by 12 proteases at the surfaces of the heterodimer and taxol-induced microtubules [de Pereda, J. M., & Andreu, J. M. (1996) Biochemistry 35, 14184-14202]. From 24 experimentally determined nicking sites, 18 are at predicted loops or at the extremes of secondary structure elements. Proteolysis zone A (including acetylable Lys40 and probably Lys60 in alpha-tubulin and Gly93 in beta-tubulin) and proteolysis zone B (extending between residues 167 and 183 in both chains) are accessible in microtubules. Proteolysis zone C, between residues 278 and 295, becomes partially occluded in microtubules. The alpha-tubulin nicking site Arg339-Ser340 is at a loop following a predicted alpha helix in proteolysis zone D. This site is protected in taxol microtubules; however, a new tryptic site appears which is probably located at the N-terminal end of the same helix. Zone D also contains beta-tubulin Cys354, which is accessible in microtubules. Proteolysis zone E includes the C-terminal hypervariable loops (10-20 residues) of each tubulin chain. These follow the two larger predicted helical zones (residues 372-395 and 405-432 in beta-tubulin), which also are the longer conserved part of the alpha- and beta-tubulin sequences. Through combination of this with other biochemical information, a set of surface and distance constraints is proposed for the folding of beta-tubulin. The FtsZ sequences are only 10-18% identical to the tubulin sequences. However, the predicted secondary structures show two clearly similar (85-87 and 51-78%) regions, at tubulin positions 95-175 and 305-350, corresponding to FtsZ 65-135 and 255-300, respectively. The first region is flanked by tubulin proteolysis zones A and B. It consists of a predicted loop1-helix-loop2-sheet-loop3-helix-loop4-sheet fold, which contains the motif (KR)GXXXXG (loop1), and the tubulin-FtsZ signature G-box motif (SAG)GGTG(SAT)G (loop3). A simple working model envisages loop1 and loop3 together at the nucleotide binding site, while loops 2 and 4 are at the surface of the protein, in agreement with proteolytic and antigenic accessibility results in tubulin. The model is compatible with studies of tubulin and FtsZ mutants. It is proposed that this region constitutes a common structural and evolutionary nucleus of tubulins and FtsZ which is different from typical GTPases.

Amino Acid Sequence↗

PubMed2003

SATCHMO: sequence alignment and tree construction using hidden Markov models.

MOTIVATION: Aligning multiple proteins based on sequence information alone is challenging if sequence identity is low or there is a significant degree of structural divergence. We present a novel algorithm (SATCHMO) that is designed to address this challenge. SATCHMO simultaneously constructs a tree and a set of multiple sequence alignments, one for each internal node of the tree. The alignment at a given node contains all sequences within its sub-tree, and predicts which positions in those sequences are alignable and which are not. Aligned regions therefore typically get shorter on a path from a leaf to the root as sequences diverge in structure. Current methods either regard all positions as alignable (e.g. ClustalW), or align only those positions believed to be homologous across all sequences (e.g. profile HMM methods); by contrast SATCHMO makes different predictions of alignable regions in different subgroups. SATCHMO generates profile hidden Markov models at each node; these are used to determine branching order, to align sequences and to predict structurally alignable regions. RESULTS: In experiments on the BAliBASE benchmark alignment database, SATCHMO is shown to perform comparably to ClustalW and the UCSC SAM HMM software. Results using SATCHMO to identify protein domains are demonstrated on potassium channels, with implications for the mechanism by which tumor necrosis factor alpha affects potassium current. AVAILABILITY: The software is available for download from http://www.drive5.com/lobster/index.htm

Algorithms↗

PubMed2001

Sequence analysis of the myosin regulatory light chain gene of the vestimentiferan Riftia pachyptila.

We have isolated and characterized a cDNA (DNA complementary to RNA) clone (Rf69) from the vestimentiferan Riftia pachyptila. The cDNA insert consists of 1169 base pairs. The aminoacid sequence deduced from the longest reading frame is 193 residues in length, and clearly characterized it as a myosin regulatory light chain (RLC). The RLC primary structure is described in relation to its function in muscle contraction. The comparison with other RLCs suggested that Riftia myosin is probably regulated through its RLC either by phosphorylation like the vertebrate smooth muscle myosins, and/or by Ca2+-binding like the mollusk myosins. Riftia RLC possesses a N-terminal extension lacking in all other species besides the earthworm Lumbricus terrestris. Aminoacid sequence comparisons with a number of RLCs from vertebrates and invertebrates revealed a relatively high identity score (64%) between Riftia RLC and the homologous gene from Lumbricus. The relationships between the members of the myosin RLCs were examined by two phylogenetic methods, i.e. distance matrix and maximum parsimony. The resulting trees depict the grouping of the RLCs according to their role in myosin activity regulation. In all trees, Riftia RLC groups with RLCs that depend on Ca2+-binding for myosin activity regulation.

Amino Acid Sequence↗