Biomedical subjects

E C Uberbacher

Publications and source records attributed to E C Uberbacher.

16 recordsLinked to original sources

PubMed1995

Correcting sequencing errors in DNA coding regions using a dynamic programming approach.

This paper presents an algorithm for detecting and 'correcting' sequencing errors that occur in DNA coding regions. The types of sequencing errors addressed are insertions and deletions (indels) of DNA bases. The goal is to provide a capability which makes single-pass or low-redundancy sequence data more informative, reducing the need for high-redundancy sequencing for gene identification and characterization purposes. This would permit improved sequencing efficiency and reduce genome sequencing costs. The algorithm detects sequencing errors by discovering changes in the statistically preferred reading frame within a putative coding region and then inserts a number of 'neutral' bases at a perceived reading frame transition point to make the putative exon candidate frame consistent. We have implemented the algorithm as a front-end subsystem of the GRAIL DNA sequence analysis system to construct a version which is very error tolerant and also intend to use this as a testbed for further development of sequencing error-correction technology. Preliminary test results have shown the usefulness of this algorithm and also exhibited some of its weakness, providing possible directions for further improvement. On a test set consisting of 68 human DNA sequences with 1% randomly generated indels in coding regions, the algorithm detected and corrected 76% of the indels. The average distance between the position of an indel and the predicted one was 9.4 bases. With this subsystem in place, GRAIL correctly predicted 89% of the coding messages with 10% false message on the 'corrected' sequences, compared to 69% correctly predicted coding messages and 11% falsely predicted messages on the 'corrupted' sequences using standard GRAIL II method (version 1.2).(ABSTRACT TRUNCATED AT 250 WORDS)

Algorithms↗

PubMed1995

Predicting protein folding classes without overly relying on homology.

An important open problem in molecular biology is how to use computational methods to understand the structure and function of proteins given only their primary sequences. We describe and evaluate an original machine-learning approach to classifying protein sequences according to their structural folding class. Our work is novel in several respects: we use a set of protein classes that previously have not been used for classifying primary sequences, and we use a unique set of attributes to represent protein sequences to the learners. We evaluate our approach by measuring its ability to correctly classify proteins that were not in its training set. We compare our input representation to a commonly used input representation--amino acid composition--and show that our approach more accurately classifies proteins that have very limited homology to the sequences on which the systems are trained.

Algorithms↗

PubMed1994

An improved system for exon recognition and gene modeling in human DNA sequences.

A new version of the GRAIL system (Uberbacher and Mural, 1991; Mural et al., 1992; Uberbacher et al., 1993), called GRAIL II, has recently been developed (Xu et al., 1994). GRAIL II is a hybrid AI system that supports a number of DNA sequence analysis tools including protein-coding region recognition, PolyA site and transcription promoter recognition, gene model construction, translation to protein, and DNA/protein database searching capabilities. This paper presents the core of GRAIL II, the coding exon recognition and gene model construction algorithms. The exon recognition algorithm recognizes coding exons by combining coding feature analysis and edge signal (acceptor/donor/translation-start sites) detection. Unlike the original GRAIL system (Uberbacher and Mural, 1991; Mural et al., 1992), this algorithm uses variable-length windows tailored to each potential exon candidate, making its performance almost exon length-independent. In this algorithm, the recognition process is divided into four steps. Initially a large number of possible coding exon candidates are generated. Then a rule-based prescreening algorithm eliminates the majority of the improbable candidates. As the kernel of the recognition algorithm, three neural networks are trained to evaluate the remaining candidates. The outputs of the neural networks are then divided into clusters of candidates, corresponding to presumed exons. The algorithm makes its final prediction by picking the best canadidate from each cluster. The gene construction algorithm (Xu, Mural and Uberbacher, 1994) uses a dynamic programming approach to build gene models by using as input the clusters predicted by the exon recognition algorithm. Extensive testing has been done on these two algorithms.(ABSTRACT TRUNCATED AT 250 WORDS)

Algorithms↗

PubMed1991

Locating protein-coding regions in human DNA sequences by a multiple sensor-neural network approach.

Genes in higher eukaryotes may span tens or hundreds of kilobases with the protein-coding regions accounting for only a few percent of the total sequence. Identifying genes within large regions of uncharacterized DNA is a difficult undertaking and is currently the focus of many research efforts. We describe a reliable computational approach for locating protein-coding portions of genes in anonymous DNA sequence. Using a concept suggested by robotic environmental sensing, our method combines a set of sensor algorithms and a neural network to localize the coding regions. Several algorithms that report local characteristics of the DNA sequence, and therefore act as sensors, are also described. In its current configuration the "coding recognition module" identifies 90% of coding exons of length 100 bases or greater with less than one false positive coding exon indicated per five coding exons indicated. This is a significantly lower false positive rate than any method of which we are aware. This module demonstrates a method with general applicability to sequence-pattern recognition problems and is available for current research efforts.

Animals↗

PubMed1989

Structure of the nucleosome core particle at 8 A resolution.

The x-ray crystallographic structure of the nucleosome core particle has been determined using 8 A resolution diffraction data. The particle has a mean diameter of 106 A and a maximum thickness of 65 A in the superhelical axis direction. The longest chord through the histone core measures 85 A and is in a non-axial direction. The 1.87 turn superhelix consists of B-DNA with about 78 base pairs or 7.6 helical repeats per superhelical turn. The mean DNA helical repeat contains 10.2 +/- 0.05 base pairs and spans 35 A, slightly more than standard B-DNA. The superhelix varies several Angstroms in radius and pitch, and has three distinct domains of curvature (with radii of curvature of 60, 45 and 51 A). These regions are separated by localized sharper bends +/- 10 and +/- 40 base pairs from the center of the particle, resulting in an overall radius of curvature about 43 A. Compression of superhelical DNA grooves on the inner surface and expansion on the outer surface can be seen throughout the DNA electron density. This density has been fit with a double helical ribbon model providing groove width estimates of 12 +/- 1 A inside vs. 19 +/- 1 A outside for the major groove, and 8 +/- 1 A inside vs. 13 +/- 1 A outside for the minor groove. The histone core is primarily contained within the bounds defined by the superhelical DNA, contacting the DNA where the phosphate backbone faces in toward the core. Possible extensions of density between the gyres have been located, but these are below the significance level of the electron density map. In cross-section, a tripartite organization of the histone octamer is apparent, with the tetramer occupying the central region and the dimers at the extremes. Several extensions of histone density are present which form contacts between nucleosomes in the crystal, perhaps representing flexible or "tail" histone regions. The radius of gyration of the histone portion of the electron density is calculated to be 30.4 A (in reasonable agreement with solution scattering values), and the histone core volume in the map is 93% of its theoretical volume.

Animals↗

PubMed1988

Rapid, large-scale purification and characterization of 'Ada protein' (O6 methylguanine-DNA methyltransferase) of E. coli.

The E. coli Ada protein (O6-methylguanine-DNA methyltransferase) has been purified using a high-level expression vector with a yield of about 3 mg per liter of E. coli culture. The 39-kDa protein has an extinction coefficient (E280 nm (1%)) of 5.3. Its isoelectric point of 7.1 is lower than that predicted from the amino acid content. The homogeneous Ada protein is fully active as a methyl acceptor from O6-methylguanine in DNA. Its reaction with O6-methylguanine in a synthetic DNA has a second-order rate constant of 1.1 x 10(9) M-1 min-1 at O degree C. Both the native form and the protein methylated at Cys-69 are monomeric. The CD spectrum suggests a low alpha-helical content and the radius of gyration of 23 A indicates a compact, globular shape. The middle region of the protein is sensitive to a variety of proteases, including an endogenous activity in E. coli, suggesting that the protein is composed of N-terminal and C-terminal domains connected by a hinge region. E. coli B has a higher level of this protease than does K12.

Circular Dichroism↗

PubMed1988

Domain interaction in rabbit muscle pyruvate kinase. II. Small angle neutron scattering and computer simulation.

The effects of ligands on the structure of rabbit muscle pyruvate kinase were studied by small angle neutron scattering. The radius of gyration, RG, decreases by about 1 A in the presence of the substrate phosphoenolpyruvate, but increases by about the same magnitude in the presence of the allosteric inhibitor phenylalanine. With increasing pH or in the absence of Mg2+ and K+, the RG of pyruvate kinase increases. Hence, there is a 2-A difference in RG between two alternative conformations. Length distribution analysis indicates that, under all experimental conditions which increase the radius of gyration, there is a pronounced increase observed in the probability for interatomic distance between 80 and 110 A. These small angle neutron scattering results indicate a "contraction" and "expansion" of the enzyme when it transforms between its active and inactive forms. Using the alpha-carbon coordinates of crystalline cat muscle pyruvate kinase, a length distribution profile was calculated, and it matches the scattering profile of the inactive form. These observations are expected since the crystals were grown in the absence of divalent cations (Stuart, D. I., Levine, M., Muirhead, H., and Stammers, D. K. (1979) J. Mol. Biol. 134, 109-142). Hence, results from neutron scattering, x-ray crystallographic, and sedimentation studies (Oberfelder, R. W., Lee, L. L.-Y., and Lee, J.C. (1984) Biochemistry 23, 3813-3821) are totally consistent with each other. With the aid of computer modeling, the crystal structure has been manipulated in order to effect changes that are consistent with the conformational change described by the solution scattering data. The structural manipulation involves the rotation of the B domain relative to the A domain, leading to the closure of the cleft between these domains. These manipulations resulted in the generation of new sets of atomic (C-alpha) coordinates, which were utilized in calculations, the result of which compared favorably with the solution data.

Algorithms↗

PubMed1988

DNA sequence patterns in precisely positioned nucleosomes.

Several investigators have recognized the importance of non-periodic DNA sequence information in determining the translational position of precisely positioned nucleosomes. The purpose of this study is to determine the extent of such information, in addition to the character of periodic information present. This is accomplished by examining the half-nucleosome DNA sequences of a considerable number of precisely positioned nucleosomes, and determining the probability of occurrence of each dinucleotide type as a function of position from the nucleosome center to the terminus (positions 0 to 72). By the nature of this procedure, no assumptions of periodicity are made. The results show the importance of several DNA sequence periodicities including 6-7, 10, and 21 base pairs, in addition to significant nonperiodic information. The results demonstrate that each dinucleotide type is unique in terms of its positional preference in precisely positioned nucleosomes (for example AA not equal to TT). The probabilities of occurrence for the dinucleotide types can be used to predict the translational positions of a number of observed nucleosomes.

Base Composition↗

PubMed1987

Small-angle neutron-scattering and electron microscope studies of the chicken liver fatty acid synthase.

A structural model for the chicken liver fatty acid synthase is proposed based on electron microscope and small-angle neutron-scattering studies of the enzyme. The model has the overall appearance of two side by side cylinders with dimensions of 160 X 146 X 73 A, with each subunit 160 A in length and 73 A in diameter. The model was constructed by dividing each cylinder into three domains having lengths of 32, 82, and 46 A, with the domain structures in the two subunits being related to each other by a dyad axis. The model is consistent with chemical cross-linking studies which indicated that the subunits are arranged in a head to tail fashion. The cross-linking studies further showed that the beta-ketoacyl synthase active site contains a cysteine and a pantetheine residue from adjacent subunits. It is proposed that the domains which catalyze the addition of C2 units from malonate to the growing fatty acid chain lie in the crevice between the two subunits and that the two independent sets of fatty acid-synthesizing centers lie on the major axis of the model on opposite ends of the molecular dyad.

Animals↗

PubMed1986

Shape analysis of the histone octamer in solution.

The conformation of the histone octamer is shown to depend upon the specific salt used to solubilize it. In 2M sodium chloride the octamer is similar in size and shape to the histone component of crystallized core nucleosomes. In contrast, in 3.5M ammonium sulfate the octamer is elongated, resembling an ellipsoid with approximate dimensions of 114 by 62 by 62 angstroms. These results indicate that the elongated conformation seen in the 3.3 angstroms electron density map of the histone octamer crystallized in ammonium sulfate is due to the particular salt conditions used.

Ammonium Sulfate↗

PubMed1986

Formation of actin dimers as studied by small angle neutron scattering.

Small angle neutron scattering has been used to study the dimensions of G-actin and the formation of low molecular weight actin oligomers under conditions where rapid polymerization does not take place. In the presence of 200 microM Ca2+, actin in solution consists of a single component with a radius of gyration (Rg) of 19.9 +/- 0.4 A, consistent with the known molecular dimensions of the G-actin molecule. In the presence of 50 microM Mg2+, however, formation of an actin species with a larger Rg occurs over a 4-h period. Multicomponent fits were tried and the data were best fit assuming two components, the monomer and a species with an Rg of 29 +/- 1 A. This latter value is consistent with the dimensions expected for certain actin dimers. The apparent dissociation constant for dimer formation is approximately 150 microM with forward and reverse rate constants of 6.0 X 10(-7) microM-1 s-1 and 8.8 X 10(-5) s-1, respectively. Kinetic fluorescence experiments show that the dimer formed in the presence of low levels of Mg2+ is a nonproductive complex which does not participate in the polymerization process. However, the addition of cytochalasin D to actin in the presence of 50 microM Mg2+ rapidly induces the formation of dimers, presumably related to cytochalasin's ability to nucleate actin polymerization.

Actins↗

PubMed1985

X-ray structure of the nucleosome core particle.

Two monoclinic crystal forms (P2(1),C2) of chicken erythrocyte nucleosomes have been under study in this laboratory. The x-ray structure of the P2(1) crystal form has been solved to 15 A resolution. The B-DNA superhelix has a relatively uniform curvature, with only several local distortions observed in the superhelix. The individual histone domains have been localized and specific contacts between each histone and the DNA can be observed. Histone contacts to the inner surface of the DNA superhelix occur predominantly at the minor groove sites. Most of the histone core is contained within the inner surface of the superhelical DNA, except for part of H2A which extends between the DNA gyres near the terminus of the DNA. No part of H2A blocks the DNA terminus or would prevent a smooth exit of the DNA into the linker region. A similar extension of a portion of histone H4 between the DNA gyres occurs close to the dyad axis. Both unique nucleosomes in the P2(1) asymmetric unit demonstrate good dyad symmetry and are similar to each other throughout the histone core and DNA regions.

Animals↗

PubMed1983

Neutron scattering studies of nucleosome structure at low ionic strength.

Ionic strength studies using homogeneous preparations of chicken erythrocyte nucleosomes containing either 146 or 175 base pairs of DNA show a single unfolding transition at about 1.5 mM ionic strength as determined by small-angle neutron scattering. The transition seen by some investigators at between 2.9 and 7.5 mM ionic strength is not observed by small-angle neutron scattering in either type of nucleosome particle. The two contrasts measured (H2O and D2O) indicate that only small conformational changes occur in the protein core, but the DNA is partially unfolded below the transition point. Patterson inversion of the data and analysis of models indicate that the DNA in both types of particle is unwinding from the ends, leaving about one turn of supercoiled DNA bound to the histone core in approximately its normal (compact) conformation. The mechanism of unfolding appears to be similar for both types of particles and in both cases occurs at the same ionic strength. The unfolding observed for nucleosomes in this study is in definite disagreement with extended superhelical models for the DNA and also disagrees with models incorporating an unfolded histone core.

Animals↗

PubMed1982

Conformation of the HMG 14 nucleosome core complex from flow birefringence.

Flow birefringence and extinction angles have been measured for HMG 14 complexes with nucleosome core particles from chicken erythrocytes under cooperative "tight" binding conditions, and for the uncomplexed core particles used in the preparations. Results are interpreted using optical models for the observed DNA anisotropy, and are compared to recent small angle neutron scattering results. (19) The studies effectively rule out highly distorted DNA conformations and configurations in which DNA ends are unwound and extended. It is concluded that the most likely conformation of the complex is one in which the DNA superhelix is radially increased, either uniformly or bilaterally, with the DNA ends remaining tightly bound to the particle. This conformation does not require large changes in spatial relationships between the DNA ends compared to the uncomplexed core as would accompany, for example, significant unwinding of the ends. However, it may lead to more subtle but possibly highly significant differences in the angles at which the DNA exits the core particle.

Animals↗

PubMed1982

Neutron scattering studies and modeling of high mobility group 14 core nucleosome complex.

Considerable evidence relates the nonhistone proteins high mobility group (HMG) 14 and HMG 17 with the structure of active or potentially active chromatin. In this study, bulk nucleosome core particles prepared from chicken erythrocytes and the complex formed by binding two HMG 14 molecules per nucleosome core were studied by use of small-angle neutron scattering techniques. By varying the H2O/2H2O ratio, and hence the contrast between the solvent and the particles, it was possible to determine the radius of gyration of the protein and of the DNA independently and as a function of HMG 14 binding. The results show an increase of 0.9 +/- 0.6 A (mean +/- SEM) in the protein radius of gyration and of 2.7 +/- 0.6 A in the DNA radius of gyration upon binding of HMG 14 to the nucleosome. These changes are considered in the light of several postulated modes for the unfolding or perturbation of the nucleosome structure. Modeling calculations demonstrate that the observed changes in radius of gyration for the DNA and for the protein are too small to be consistent with an overall unfolding or opening of the core particle upon HMG 14 binding. However, the observed changes are consistent with several models that involve only minor changes in the structure. It is postulated that the differences observed may be an indication of the type of conformational change occurring in active nucleosomes.

Animals↗

PubMedDate not supplied

An artificial intelligence approach to DNA sequence feature recognition.

The ultimate goal of the Human Genome project is to extract the biologically relevant information recorded in the estimated 100,000 genes encoded by the 3 x 10(9) bases of the human genome. This necessitates development of reliable computer-based methods capable of analysing and correctly identifying genes in the vast amounts of DNA-sequence data generated. Such tools may save time and labour by simplifying, for example, screening of cDNA libraries. They may also facilitate the localization of human disease genes by identifying candidate genes in promising regions of anonymous DNA sequence.

Artificial Intelligence↗