Algorithms for sequence analysis via mutagenesis

Despite many successes of conventional DNA sequencing methods, some DNAs remain difficult or impossible to sequence. Unsequenceable regions occur in the genomes of many biologically important organisms, including the human genome. Such regions range in length from tens to millions of bases, and may contain valuable information such as the sequences of important genes. The authors have recently developed a technique that renders a wide range of problematic DNAs amenable to sequencing. The technique is known as sequence analysis via mutagenesis (SAM). This paper presents a number of algorithms for analysing and interpreting data generated by this technique.

Comparison of BR and QR Eigenvalue Algorithms for Power System Small Signal Stability Analysis

The BR algorithm is a novel and efficient method to find all eigenvalues of upper Hessenberg matrices and has never been applied to eigenanalysis for power system small signal stability. This paper analyzes differences between the BR and the QR algorithms with performance comparison in terms of CPU time based on stopping criteria and storage requirement. The BR algorithm utilizes accelerating strategies to improve its performance when computing eigenvalues of narrowly banded, nearly tridiagonal upper Hessenberg matrices. These strategies significantly reduce the computation time at a reasonable level of precision. Compared with the QR algorithm, the BR algorithm requires fewer iteration steps and less storage space without depriving of appropriate precision in solving eigenvalue problems of large-scale power systems. Numerical examples demonstrate the efficiency of the BR algorithm in pursuing eigenanalysis tasks of 39-, 68-, 115-, 300-, and 600-bus systems. Experiment results suggest that the BR algorithm is a more efficient algorithm for large-scale power system small signal stability eigenanalysis.

Identification of two novel mutations in the hydroxymethylbilane synthase gene in three patients from two unrelated families with acute intermittent porphyria

We have screened the hydroxymethylbilane synthase cDNAs of 3 patients from 2 families suffering from acute intermittent porphyria (AIP) from Scotland and South Africa using heteroduplex and chemical cleavage of mismatch analyses, Direct sequencing was used to characterise the mutations, The two novel mutations identified were a missense mutation at nucleotide position 64 in exon 3 (R22C) and a single base-pair deletion in exon 15, These mutations are predicted to affect the normal function of the enzyme and, therefore, are expected to be the primary cause of disease in these patients.

A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

Extended GCD and Hermite normal form algorithms via lattice basis reduction

Extended gcd calculation has a long history and plays an important role in computational number theory and linear algebra. Recent results have shown that finding optimal multipliers in extended gcd calculations is difficult. We present an algorithm which uses lattice basis reduction to produce small integer multipliers x(1), ..., x(m) for the equation s = gcd (s(1), ..., s(m)) = x(1)s(1) + ... + x(m)s(m), where s1, ... , s(m) are given integers. The method generalises to produce small unimodular transformation matrices for computing the Hermite normal form of an integer matrix.

Spectroscopic identification of a dinuclear metal centre in manganese(II)-activated aminopeptidase P from Escherichia coli: implications for human prolidase

Electron paramagnetic resonance (EPR) spectra and X-ray absorption (EXAFS and XANES) data have been recorded for the manganese enzyme aminopeptidase P (AMPP, PepP protein) from Escherichia coli. The biological function of the protein, a tetramer of 50-kDa subunits, is the hydrolysis of N-terminal Xaa-Pro peptide bonds. Activity assays confirm that the enzyme is activated by treatment with Mn2+. The EPR spectrum of Mn2+-activated AMPP at liquid-He temperature is characteristic of an exchange-coupled dinuclear Mn(II) site, the Mn-Mn separation calculated from the zero-field splitting D of the quintet state being 3.5 (+/- 0.1) Angstrom. In the X-ray absorption spectrum of Mn2+-activated AMPP at the Mn K edge, the near-edge features are consistent with octahedrally coordinated Mn atoms in oxidation state +2. EXAFS data, limited to k less than or equal to 12 Angstrom(-1) by traces of Fe in the protein, are consistent with a single coordination shell occupied predominantly by O donor atoms at an average Mn-ligand distance of 2.15 Angstrom, but the possibility of a mixture of O and N donor atoms is not excluded. The Mn-Mn interaction at 3.5 Angstrom, is not detected in the EXAFS, probably due to destructive interference from light outer-shell atoms. The biological function, amino acid sequence and metal-ion dependence of E. coli AMPP are closely related to those of human prolidase, an enzyme that specifically cleaves Xaa-Pro dipeptides. Mutations that lead to human prolidase deficiency and clinical symptoms have been identified. Several known inhibitors of prolidase also inhibit AMPP. When these inhibitors are added to Mn2+-activated AMPP, the EPR spectrum and EXAFS remain unchanged. It can be inferred that the inhibitors either do not bind directly to the Mn centres, or substitute for existing Mn ligands without a significant change in donor atoms or coordination geometry. The conclusions from the spectroscopic measurements on AMPP have been verified by, and complement, a recent crystal structure analysis.

Identification in Australia of the quarantine pathogen of sunflower Phoma macdonaldii (Teleomorph : Leptosphaeria lindquistii)

Isolations from black stem lesions of sunflower growing in south-eastern Queensland yielded fungi putatively identified as species of Phoma. Pathogenicity assays showed that these isolates were capable of killing sunflower plants under glasshouse conditions. The isolates were compared with authentic cultures of Phoma macdonaldii and other isolates of Phoma taken from sunflower from around the world. Random amplified polymorphic DNA analysis showed that all the Australian isolates examined were very similar to the holotype culture of Phoma macdonaldii from Canada. Sequencing of the internal transcribed spacer regions also revealed the relatedness of the Australian isolates to the holotype. This is the first official record of P. macdonaldii in Australia.

Effects of data characteristics on the results of reserve selection algorithms

We tested the effects of four data characteristics on the results of reserve selection algorithms. The data characteristics were nestedness of features (land types in this case), rarity of features, size variation of sites (potential reserves) and size of data sets (numbers of sites and features). We manipulated data sets to produce three levels, with replication, of each of these data characteristics while holding the other three characteristics constant. We then used an optimizing algorithm and three heuristic algorithms to select sites to solve several reservation problems. We measured efficiency as the number or total area of selected sites, indicating the relative cost of a reserve system. Higher nestedness increased the efficiency of all algorithms (reduced the total cost of new reserves). Higher rarity reduced the efficiency of all algorithms (increased the total cost of new reserves). More variation in site size increased the efficiency of all algorithms expressed in terms of total area of selected sites. We measured the suboptimality of heuristic algorithms as the percentage increase of their results over optimal (minimum possible) results. Suboptimality is a measure of the reliability of heuristics as indicative costing analyses. Higher rarity reduced the suboptimality of heuristics (increased their reliability) and there is some evidence that more size variation did the same for the total area of selected sites. We discuss the implications of these results for the use of reserve selection algorithms as indicative and real-world planning tools.

Identification and characterisation of a cytochrome P450 gene and processed pseudogene from an arachnid: the cattle tick, Boophilus microplus

We isolated and sequenced the first known cytochrome P450 gene and pseudogene from an arachnid, the cattle tick, Boophilus microplus. Bath the gene and pseudogene belong to the family CYP4, but a new subfamily, CYP4W, had to be created for these genes because they are substantially different to other CYP4 genes. The gene, CPP4W1, has greatest homology with CYP4C1 from a cockroach, Blaberus discoidalis. The predicted molecular weight of the protein encoded by CYP4W1 (63 KDa) is greater than that of the other CYP4 genes. The pseudogene, CYP4W1P, is probably a processed pseudogene derived from the functional gene CYP4W1. This is only the third CYP processed pseudogene to be identified. The pseudogene is 98% identical to the functional gene, CYP4W1, therefore we hypothesise that this pseudogene evolved recently from the functional gene. The CYP4 genes from arthropods have diverged from each other more than those of mammals; consequently the phylogeny of the arthropod genes could not be resolved. (C) 1999 Elsevier Science Ltd. All rights reserved.

Identification of a new ligand binding domain in the al subunit of the inhibitory glycine receptor

Four discontinuous extracellular sequence domains have been proposed to form the ligand binding sites of the ligand-gated ion channel receptor superfamily. In this study, we investigated the role of 12 contiguous residues of the inhibitory glycine receptor that define the proposed loop A ligand binding domain; Using the techniques of site-directed mutagenesis and patch-clamp electrophysiology, four of the 12 residues were shown to have impaired ligand binding. Three mutants, I93A, A101H, and N102A, resulted in significant (17-44-fold) increases in the agonist EC50 values as compared with the wild-type glycine receptor, whereas Hill coefficients, I-max values, and antagonist affinity remained largely unaffected. Consideration of receptor efficacy values indicates that these residues are involved in ligand binding rather than channel activation. A fourth mutant, W94A, failed to give rise to any glycine-activated currents, although cell-surface expression was observed, suggesting that this residue may also be involved in agonist binding. These data provide the most extensive characterization of the loop A ligand binding domain available to date and define two new residue locations, Ile(93) and Asn(102), as contributing to the four-loop model of ligand binding.

Identification of residues involved in binding of IL5 to beta(com) using beta(IL3) and beta(com) chimeras

In mice there are two forms of the beta chain used in the IL3 receptor system, beta(com) and beta(IL3). beta(com) is used by the IL3, ILS and GM-CSF receptors whereas Pns is only used in the IL3 receptor. In this work an assay was developed to identify residues of beta(IL3) that restrict IL5 activity. It was found that such residues reside within the 2nd CRM of the molecule. Furthermore, when residues in the beta(IL3) B'-C' loop were replaced with beta(com) sequence a form of beta(IL3) was produced that was able to respond to IL5. This region is also responsible for IL3 binding to beta(IL3) in the absence of alpha chain. It is therefore an important structural motif of beta(com) and beta(IL3) responsible for both ligand interaction and specificity. (C) 1999 Federation of European Biochemical Societies.

Identification of initiation sites for T4 lysozyme folding using CD and NMR spectroscopy of peptide fragments

Using CD and 2D H-1 NMR spectroscopy, we have identified potential initiation sites for the folding of T4 lysozyme by examining the conformational preferences of peptide fragments corresponding to regions of secondary structure. CD spectropolarimetry showed most peptides were unstructured in water, but adopted partial helical conformations in TFE and SDS solution. This was also consistent with the H-1 NMR data which showed that the peptides were predominantly disordered in water, although in some cases, nascent or small populations of partially folded conformations could be detected. NOE patterns, coupling constants, and deviations from random coil Her chemical shift values complemented the CD data and confirmed that many of the peptides were helical in TFE and SDS micelles. In particular, the peptide corresponding to helix E in the native enzyme formed a well-defined helix in both TFE and SDS, indicating that helix E potentially forms an initiation site for T4 lysozyme folding. The data for the other peptides indicated that helices D, F, G, and H are dependent on tertiary interactions for their folding and/or stability. Overall, the results from this study, and those of our earlier studies, are in agreement with modeling and IID-deuterium exchange experiments, and support an hierarchical model of folding for T4 lysozyme.