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 865 records · Page 48Linked to original sources

PubMed2006

Pseudo-likelihood for non-reversible nucleotide substitution models with neighbour dependent rates.

In the field of molecular evolution genome substitution models with neighbour dependent substitution rates have recently received much attention. It is well-known that substitution of nucleotides does not occur independently of neighbouring nucleotides, but there has been less focus on the phenomenon that this substitution process is also not time-reversible. In this paper I construct a pseudo-likelihood type method for inference in non-reversible substitution models with neighbour dependent substitution rates. I also construct an EM-algorithm for maximising the pseudo-likelihood. For human-mouse aligned sequence data a number of different models are investigated, where I show that strand-symmetric models are appropriate, and that overlapping di-nucleotide models do not fit the data well.

Algorithms↗

PubMed1995

The use of side-chain packing methods in modeling bacteriophage repressor and cro proteins.

In recent years, it has been repeatedly demonstrated that the coordinates of the main-chain atoms alone are sufficient to determine the side-chain conformations of buried residues of compact proteins. Given a perfect backbone, the side-chain packing method can predict the side-chain conformations to an accuracy as high as 1.2 A RMS deviation (RMSD) with greater than 80% of the chi angles correct. However, similarly rigorous studies have not been conducted to determine how well these apply, if at all, to the more important problem of homology modeling per se. Specifically, if the available backbone is imperfect, as expected for practical application of homology modeling, can packing constraints alone achieve sufficiently accurate predictions to be useful? Here, by systematically applying such methods to the pairwise modeling of two repressor and two cro proteins from the closely related bacteriophages 434 and P22, we find that when the backbone RMSD is 0.8 A, the prediction on buried side chain is accurate with an RMS error of 1.8 A and approximately 70% of the chi angles correctly predicted. When the backbone RMSD is larger, in the range of 1.6-1.8 A, the prediction quality is still significantly better than random, with RMS error at 2.2 A on the buried side chains and 60% accuracy on chi angles. Together these results suggest the following rules-of-thumb for homology modeling of buried side chains. When the sequence identity between the modeled sequence and the template sequence is > 50% (or, equivalently, the expected backbone RMSD is < 1 A), side-chain packing methods work well. When sequence identity is between 30-50%, reflecting a backbone RMS error of 1-2 A, it is still valid to use side-chain packing methods to predict the buried residues, albeit with care. When sequence identity is below 30% (or backbone RMS error greater than 2 A), the backbone constraint alone is unlikely to produce useful models. Other methods, such as those involving the use of database fragments to reconstruct a template backbone, may be necessary as a complementary guide for modeling.

Amino Acid Sequence↗

PubMed2005

Development of a single-round and multiplex PCR method for the simultaneous detection of Babesia caballi and Babesia equi in horse blood.

With the aim of developing more simple diagnostic alternatives, a differential single-round and multiplex polymerase chain reaction (PCR) method was designed for the simultaneous detection of Babesia caballi and Babesia equi, by targeting 18S ribosomal RNA genes. The multiplex PCR amplified DNA fragments of 540 and 392 bp from B. caballi and B. equi, respectively, in one reaction. The PCR method evaluated on 39 blood samples collected from domestic horses in Mongolia yielded similar results to those obtained from confirmative PCR methods that had been established earlier. Thus, the single-round and multiplex PCR method offers a simple tool for the differential diagnosis of B. caballi and B. equi infections in routine diagnostic laboratory settings as well as in epidemiological studies.

Animals↗

PubMed2005

Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

Animals↗

PubMed1997

A theoretical method for evaluating the relative importance of positive selection and neutral drift from observed base changes.

To evaluate the relative importance of positive selection and neutral drift from the nucleotide base changes observed in the homologous alignment of genes, a theoretical equation of base changes is formulated by including both the influence of selection and the base substitutions due to mutations. Under the assumption that the average rate of base substitutions estimated from synonymous changes is the "true" mutation rate applicable at all positions, this method is applied to the vertebrate globin gene family, and evaluates the departures of base change rates from the "true" mutation rate at the first and second codon positions as a consequence of preferential selection for the conservation of important function. In addition to the strong effect of selection on the amino acid residues in the internal region mostly common to myoglobin and hemoglobin chains, the distinctive directions of selective parameter values are seen at sites on the globin surface, distinguishing the subunit contact residues of hemoglobins from the polar residues on the surface of myoglobins. Moreover, this effect of selection distinguishing between the myoglobin and hemoglobin chain genes becomes weaker in cold-blooded vertebrates, especially in fish, strongly suggesting the possibility that the clear distinction between these globins is a result of selection out of the changes regarded as neutral ones in an ancestor of vertebrates. Thus, the present method may also serve to investigate the homology of many other proteins from the aspect of molecular evolution, mainly focusing on the evolution of their biological functions.

Animals↗

PubMed2004

Aligned inclusion of hemicyanine dyes into silica zeolite films for second harmonic generation.

Silicalite-1 films (thickness = 400 nm) supported on both sides of glass plates (SL/G) were prepared, and hemicyanine dyes (HC-n) with different alkyl chain lengths (n, n = 3, 6, 9, 12, 15, 18, 22, and 24) were included into the silicalite-1 films by dipping SL/Gs into each methanol solution of HC-n (1 mM) for 1 d. The included numbers of HC-n per channel (N(C)) generally decreased with increasing n; that is, they were 6.4, 23.1, 15.4, 8.2, 5.7, 3.5, 0.9, and 1.2 molecules per channel, respectively. The d(33) value gradually increased with increasing n but decreased when n > 18; that is, they were 1.12, 0.50, 2.25, 3.59, 4.99, 5.30, 1.71, and 2.57 pm V(-1), respectively. However, d(33)/N(C) progressively increased with increasing n. The d(31) values were approximately 100 times smaller than the corresponding d(33) values, and the average d(33)/d(31) ratio was 109, which is higher than those of Langmuir-Blodgett (LB) films and poled polymers of nonlinear optical (NLO) dyes, by approximately 2-5 and approximately 30-50 times, respectively. The estimated average tilted angle of the dyes with respect to the channel direction was 7.7 degrees, and the calculated average order parameter was 0.97, which is approximately 480 times higher than the values observed from poled polymers. The degree of uniform alignment (DUA) generally increased with increasing n. The progressive increase of both DUA and d(33)/N(C) with n is attributed to the increase in the tendency of HC-n to enter hydrophobic silicalite-1 channels with the hydrophobic alkyl chain first. A more than 134-fold increase in DUA was observed upon increasing n from 6 to 24. The DUA of HC-24 in the silicalite-1 film reached close to 1. Although the observed d(33) values were lower than those of the LB films of NLO dyes due to very small dye densities of the silicalite films, this methodology bears a great potential to be developed into the methods for preparing practically viable NLO films.

Journal Article↗

PubMed1998

Structure and proposed amino-acid sequence of a pepsin from atlantic cod (Gadus morhua).

The crystal structure of a pepsin from the gastric mucosa of Atlantic cod has been determined to 2.16 A resolution. Data were collected on orthorhombic crystals with cell dimensions a = 35.98, b = 75.40 and c = 108.10 A, on a FAST area-detector system. The phase problem was solved by the molecular-replacement method using porcine pepsin (PDB entry 5PEP) as a search model. The structure has been refined to a crystallographic R factor of 20.8% using all reflections between 8.0 and 2.16 A, without prior knowledge of the primary sequence. The resulting crystal structure is very similar to the porcine enzyme, consisting of two domains with predominantly beta-sheet structure in the same sequential positions as the enzyme from pig. In the course of the model building, 122 residues were substituted and two residues deleted from the starting model to give a polypeptide chain of 324 amino acids and a sequence identity of 57.7% with the pig pepsin. No carbohydrate residues were located. Sequence alignment with available aspartic proteinases, indicates that the fish enzyme seems to be more related to mammalian gastric pepsins than to the mammalian gastricsins and chymosins, lysosomal cathepsin D's and a pepsin from tuna fish. The amino-acid composition of the cod enzyme, however, is more in accordance with the cathepsin D's.

Amino Acid Sequence↗

PubMed2000

Modeling of loops in protein structures.

Comparative protein structure prediction is limited mostly by the errors in alignment and loop modeling. We describe here a new automated modeling technique that significantly improves the accuracy of loop predictions in protein structures. The positions of all nonhydrogen atoms of the loop are optimized in a fixed environment with respect to a pseudo energy function. The energy is a sum of many spatial restraints that include the bond length, bond angle, and improper dihedral angle terms from the CHARMM-22 force field, statistical preferences for the main-chain and side-chain dihedral angles, and statistical preferences for nonbonded atomic contacts that depend on the two atom types, their distance through space, and separation in sequence. The energy function is optimized with the method of conjugate gradients combined with molecular dynamics and simulated annealing. Typically, the predicted loop conformation corresponds to the lowest energy conformation among 500 independent optimizations. Predictions were made for 40 loops of known structure at each length from 1 to 14 residues. The accuracy of loop predictions is evaluated as a function of thoroughness of conformational sampling, loop length, and structural properties of native loops. When accuracy is measured by local superposition of the model on the native loop, 100, 90, and 30% of 4-, 8-, and 12-residue loop predictions, respectively, had <2 A RMSD error for the mainchain N, C(alpha), C, and O atoms; the average accuracies were 0.59 +/- 0.05, 1.16 +/- 0.10, and 2.61 +/- 0.16 A, respectively. To simulate real comparative modeling problems, the method was also evaluated by predicting loops of known structure in only approximately correct environments with errors typical of comparative modeling without misalignment. When the RMSD distortion of the main-chain stem atoms is 2.5 A, the average loop prediction error increased by 180, 25, and 3% for 4-, 8-, and 12-residue loops, respectively. The accuracy of the lowest energy prediction for a given loop can be estimated from the structural variability among a number of low energy predictions. The relative value of the present method is gauged by (1) comparing it with one of the most successful previously described methods, and (2) describing its accuracy in recent blind predictions of protein structure. Finally, it is shown that the average accuracy of prediction is limited primarily by the accuracy of the energy function rather than by the extent of conformational sampling.

Models, Molecular↗

PubMed2005

Human intervertebral disc aggrecan inhibits endothelial cell adhesion and cell migration in vitro.

STUDY DESIGN: The effect of human intervertebral disc aggrecan on endothelial cell growth was examined using cell culture assays. OBJECTIVE: To determine the response of endothelial cells to human intervertebral disc aggrecan, and whether the amount and type of aggrecan present in the intervertebral disc may be implicated in disc vascularization. SUMMARY OF BACKGROUND DATA: Intervertebral disc degeneration has been associated with a loss of proteoglycan, and the ingrowth of blood vessels and nerves. Neovascularization is a common feature also of disc herniation. Intervertebral disc aggrecan is inhibitory to sensory nerve growth, but the effects of disc aggrecan on endothelial cell growth are not known. METHODS: Aggrecan monomers were isolated separately from the anulus fibrosus and nucleus pulposus of human lumbar intervertebral discs, and characterized to determine the amount and type of sulfated glycosaminoglycan side chains present. The effects of these aggrecan isolates on the cellular adhesion and migration of the human endothelial cell lines, HMEC-1 and EAhy-926, were examined in vitro. RESULTS: Homogenous substrata of disc aggrecan inhibited endothelial cell adhesion and cell spreading in a concentration dependent manner. In substrata choice assays, endothelial cells seeded onto collagen type I migrated over the collagen until they encountered substrata of disc aggrecan, where they either stopped migrating, retreated onto the collagen, or, more commonly, changed direction to align along the collagen-aggrecan border. The inhibitory effect of aggrecan on endothelial cell migration was concentration dependent, and reduced by enzymatic treatment of the aggrecan monomers with a combination of chondroitinase ABC and keratinase/keratinase II. Anulus fibrosus aggrecan was more inhibitory to endothelial cell adhesion than nucleus pulposus aggrecan. However, this difference did not relate to the extent to which the different aggrecan isolates were charged, as determined by colorimetric assay with 1,9-dimethylmethylene blue, or to marked differences in the distribution of chondroitin sulfated and keratan sulfated side chains. CONCLUSIONS: Human intervertebral disc aggrecan is inhibitory to endothelial cell migration, and this inhibitory effect appears to depend, in part, on the presence of glycosaminoglycan side chains on the aggrecan monomer.

Aggrecans↗

PubMed1997

PCR-RFLP analysis of the Cryptosporidium oocyst wall protein (COWP) gene discriminates between C. wrairi and C. parvum, and between C. parvum isolates of human and animal origin.

Cryptosporidium wrairi was isolated from guinea pigs during a spontaneous outbreak of cryptosporidiosis. Despite the morphological and antigenic similarities to C. parvum, C. wrairi displayed a different host range and site of infection and may represent a separate species or sub-species. We used the polymerase chain reaction to clone two distinct 550 bp-long DNA fragments, Wc-I and Wc-II, of the gene encoding the Cryptosporidium oocyst wall protein (COWP) of C. wrairi, which showed 98% identity to the C. parvum homologue. Within Wc-I, polymorphic Rsal restriction sites were used to develop a polymerase chain reaction-restriction fragment length polymorphism method able to distinguish C. wrairi from C. parvum and to identify two groups of C. parvum isolates differentially associated with animal and human infections.

Amino Acid Sequence↗

PubMed2005

Identification of sodium channel mutations in human head louse (Anoplura: Pediculidae) from Denmark.

The presence of the two mutations T9291-L932F in the voltage-sensitive sodium channel gene associated with permethrin resistance (kdr-like) was shown in head louse, Pediculus humanus capitis De Geer, populations in Denmark. The existence of one susceptible and one T929I-L932F haplotype based on the limited single nucleotide polymorphism (SNP) of these sequences was established. One louse had an SNP causing a G943A substitution in the trans-membrane segment five of domain II on the sodium channel, which has not been identified in other insect species. A polymerase chain reaction-restriction endonuclease method using genomic DNA to discriminate between resistant and susceptible alleles in head lice was developed and implemented. The SNP that results in the T9291 substitution also creates a cutting site for the restriction endonuclease SspI, and the presence of one or two SspI cutting sites in head lice is diagnostic for the T9291-L932F haplotype.

Amino Acid Sequence↗

PubMed1995

Cloning of the rat homologue of the von Hippel-Lindau tumor suppressor gene and its non-somatic mutation in rat renal cell carcinomas.

Recently, von Hippel-Lindau (VHL) gene mutations were detected in non-inherited, sporadic human renal cell carcinomas (RCs) at a high frequency. In order to determine whether or not the VHL gene is also a critical gene in rat RCs, we cloned and sequenced the rat homologue of human VHL gene and searched for mutations of the VHL gene in rat RCs. Mutations in the VHL gene were not detected in spontaneous RCs of the Eker rat model or in ferric nitrilotriacetate-induced rat RCs using the polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) method. These data indicate that mutation of the VHL tumor suppressor gene is not an event in rat renal carcinogenesis, at least in our present systems.

Amino Acid Sequence↗

PubMed2001

Prediction of unidentified human genes on the basis of sequence similarity to novel cDNAs from cynomolgus monkey brain.

BACKGROUND: The complete assignment of the protein-coding regions of the human genome is a major challenge for genome biology today. We have already isolated many hitherto unknown full-length cDNAs as orthologs of unidentified human genes from cDNA libraries of the cynomolgus monkey (Macaca fascicularis) brain (parietal lobe and cerebellum). In this study, we used cDNA libraries of three other parts of the brain (frontal lobe, temporal lobe and medulla oblongata) to isolate novel full-length cDNAs. RESULTS: The entire sequences of novel cDNAs of the cynomolgus monkey were determined, and the orthologous human cDNA sequences were predicted from the human genome sequence. We predicted 29 novel human genes with putative coding regions sharing an open reading frame with the cynomolgus monkey, and we confirmed the expression of 21 pairs of genes by the reverse transcription-coupled polymerase chain reaction method. The hypothetical proteins were also functionally annotated by computer analysis. CONCLUSIONS: The 29 new genes had not been discovered in recent explorations for novel genes in humans, and the ab initio method failed to predict all exons. Thus, monkey cDNA is a valuable resource for the preparation of a complete human gene catalog, which will facilitate post-genomic studies.

Animals↗

PubMed1992

An efficient and reliable method for cloning PCR-amplification products: a survey of point mutations in integrin cDNA.

A highly efficient, non-labor-intensive method for cloning DNA fragments produced by PCR amplification was used to carry out a rapid survey of potential point mutations in integrin alpha 6 cDNA from 17 different cell-type sources. The method includes glass powder purification of the PCR reaction mixture, followed by simultaneous treatment with T4 polynucleotide kinase and DNA polymerase I, and another glass powder purification. Sequences from multiple subclones of each cell type were readily generated, aligned and checked for mismatches. Several commonly used alternative procedures were compared for cloning efficiency and size-fidelity of inserted DNA fragments.

Amino Acid Sequence↗

PubMed1993

An automated method for modeling proteins on known templates using distance geometry.

We present an automated method incorporated into a software package, FOLDER, to fold a protein sequence on a given three-dimensional (3D) template. Starting with the sequence alignment of a family of homologous proteins, tertiary structures are modeled using the known 3D structure of one member of the family as a template. Homologous interatomic distances from the template are used as constraints. For nonhomologous regions in the model protein, the lower and the upper bounds for the interatomic distances are imposed by steric constraints and the globular dimensions of the template, respectively. Distance geometry is used to embed an ensemble of structures consistent with these distance bounds. Structures are selected from this ensemble based on minimal distance error criteria, after a penalty function optimization step. These structures are then refined using energy optimization methods. The method is tested by simulating the alpha-chain of horse hemoglobin using the alpha-chain of human hemoglobin as the template and by comparing the generated models with the crystal structure of the alpha-chain of horse hemoglobin. We also test the packing efficiency of this method by reconstructing the atomic positions of the interior side chains beyond C beta atoms of a protein domain from a known 3D structure. In both test cases, models retain the template constraints and any additionally imposed constraints while the packing of the interior residues is optimized with no short contacts or bond deformations. To demonstrate the use of this method in simulating structures of proteins with nonhomologous disulfides, we construct a model of murine interleukin (IL)-4 using the NMR structure of human IL-4 as the template. The resulting geometry of the nonhomologous disulfide in the model structure for murine IL-4 is consistent with standard disulfide geometry.

Algorithms↗

PubMed1990

Protein structure prediction.

Current methods developed for predicting protein structure are reviewed. The most widely used algorithms of Chou and Fasman and Garnier et al for predicting secondary structure are compared to the most recent ones including sequence similarity methods, neural network, pattern recognition or joint prediction methods. The best of these methods correctly predict 63-65% of the residues in the database with cross-validation for 3 conformations, helix, beta strand and coli with a standard deviation of 6-8% per protein. However, when a homologous protein is already in the database, the accuracy of prediction by the similarity peptide method of Levin and Garnier reaches about 90%. Some conclusions can be drawn on the mechanism of protein folding. As all the prediction methods only use the local sequence for prediction (+/- 8 residues maximum) one can infer that 65% of the conformation of a residue is dictated on average by the local sequence, the rest is brought by the folding. The best predicted proteins or peptide segments are those for which the folding has less effect on the conformation. Presently, prediction of tertiary structure is only of practical use when the structure of a homologous protein is already known. Amino acid alignment to define residues of equivalent spatial position is critical for modelling of the protein. We showed for serine proteases that secondary structure prediction can help to define a better alignment. Non-homologous segments of the polypeptide chain, such as loops, libraries of known loops and/or energy minimization with various force fields, are used without yet giving satisfactory solutions. An example of modelling by homology, aided by secondary structure prediction on 2 regulatory proteins, Fnr and FixK is presented.

Algorithms↗

PubMed2004

True molar surface energy and alignment of surface molecules.

By use of data for surface tension, for the first time a method is presented for calculating true molar surface properties of liquids: free energy, entropy, and enthalpy. These new data allow full comparison with other molar quantities, such as enthalpy and entropy of vaporization. All data are at the normal boiling point. There are differences in behavior between various classes of nonpolar compounds. Rare gases and tetrahydrides of Group 14 form a separate category. The results agree with the experimental findings that water and alcohol molecules are aligned with the -OH groups attached to the surface. The data indicate that hydrogen peroxide and 1,2-ethanediol also have one -OH group directed toward the surface and one directed out toward the vapor phase. Small straight-chain amines have some structure at the surface, but the larger ones behave like the corresponding alkanes. Very polar compounds, such as nitriles, nitro compounds, and aldehydes, have little or no increased degree of order of molecules at the surface. Except for hydrogen-bonded compounds, molecules at the surface have surroundings similar to those in the bulk liquid.

Journal Article↗

PubMed1983

The serine proteinase chain of human complement component C1s. Cyanogen bromide cleavage and N-terminal sequences of the fragments.

Human complement component C1s was purified from fresh blood by conventional methods of precipitation and chromatography. The single-chain zymogen form was activated by treatment with C1r. Reduction and carboxymethylation then allowed the light chain and heavy chain to be separated on DEAE-Sepharose CL-6B in 8 M-urea. Liquid-phase sequencing of the light chain determined 50 residues from the N-terminus. CNBr-cleavage fragments of the light chain were separated by high-pressure liquid chromatography on gel-permeation and reverse-phase columns. N-Terminal sequencing of these fragments determined the order of a further 138 residues, giving a total of 188 residues or about 75% of the light chain. Seven of these eight sequences could be readily aligned with the amino acid sequences of other serine proteinases. The typical serine proteinase active-site residues are clearly conserved in C1s, and the specificity-related side chain of the substrate-binding pocket is aspartic acid, as in trypsin, consistent with the proteolytic action of C1s on C4 at an arginine residue. Somewhat surprisingly, when the C1s sequence is compared with that of complement subcomponent C1r, the percentage difference (59%) is approximately the same as that found between the other mammalian serine proteinases (56-71%).

Amino Acid Sequence↗