A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

A monolithic smoothing-gap algorithm for contact-impact based on the signed distance function

A new algorithm has been developed for smoothing the surfaces in finite element formulations of contact-impact. A key feature of this method is that the smoothing is done implicitly by constructing smooth signed distance functions for the bodies. These functions are then employed for the computation of the gap and other variables needed for implementation of contact-impact. The smoothed signed distance functions are constructed by a moving least-squares approximation with a polynomial basis. Results show that when nodes are placed on a surface, the surface can be reproduced with an error of about one per cent or less with either a quadratic or a linear basis. With a quadratic basis, the method exactly reproduces a circle or a sphere even for coarse meshes. Results are presented for contact problems involving the contact of circular bodies. Copyright (C) 2002 John Wiley Sons, Ltd.

Behavior of argon gas release from manganese oxide minerals as revealed by Ar-40/Ar-39 laser incremental heating analysis

Manganese oxides in association with paleo-weathering may provide significant insights into the multiple factors affecting the formation and evolution of weathering profiles, such as temperature, precipitation, and biodiversity. Laser probe step-heating analysis of supergene hollandite and cryptomelane samples collected from central Queensland, Australia, yield well-defined plateaus and consistent isochron ages, confirming the feasibility, dating very-fined supergene manganese oxides by Ar-40/(39) Ar technique. Two distinct structural sites hosting Ar isotopes can be identified in light of their degassing behaviors obtained by incremental heating analyses. The first site, releasing its gas fraction at the laser power 0.2-0.4 W, yields primarily Ar-40(atm), Ar-38(atm), and Ar-36(atm), (atmospheric Ar isotopes). The second sites yield predominantly Ar-40* (radiogenic Ar-40), Ar-39(K), and Ar-38(K) (nucleogenic components), at similar to0.5-1.0 W. There is no significant Ar gas released at the laser power higher than 1.0 W, indicating the breakdown of the tunnel sites hosting the radiogenic and nucleogenic components. The excellent match between the degassing behaviors of Ar-40*, Ar-39(K), and Ar-38(K) suggests that these isotopes occupy the same crystallographic sites and that Ar-39(K) loss from the tunnel site by recoil during neutron irradiation and/or bake-out procedure preceding isotopic analysis does not occur. Present investigation supports that neither the overwhelming atmospheric Ar-40 nor the very-fined nature of the supergene manganese oxides poses problems in extracting meaningful weathering geo-chronological information by analyzing supergene manganese oxides minerals.

A subset-oriented algorithm for minimizing the number of steps required for synthesis of cyclic-peptide libraries

Libraries of cyclic peptides are being synthesized using combinatorial chemistry for high throughput screening in the drug discovery process. This paper describes the min_syn_steps.cpp program (available at http://www.imb.uq.edu.au/groups/smythe/tran), which after inputting a list of cyclic peptides to be synthesized, removes cyclic redundant sequences and calculates synthetic strategies which minimize the synthetic steps as well as the reagent requirements. The synthetic steps and reagent requirements could be minimized by finding common subsets within the sequences for block synthesis. Since a brute-force approach to search for optimum synthetic strategies is impractically large, a subset-orientated approach is utilized here to limit the size of the search. (C) 2002 Elsevier Science Ltd. All rights reserved.

The tree cut and merge algorithm for estimation of network reliability

This article presents Monte Carlo techniques for estimating network reliability. For highly reliable networks, techniques based on graph evolution models provide very good performance. However, they are known to have significant simulation cost. An existing hybrid scheme (based on partitioning the time space) is available to speed up the simulations; however, there are difficulties with optimizing the important parameter associated with this scheme. To overcome these difficulties, a new hybrid scheme (based on partitioning the edge set) is proposed in this article. The proposed scheme shows orders of magnitude improvement of performance over the existing techniques in certain classes of network. It also provides reliability bounds with little overhead.

A genetic algorithm/method of moments approach to the optimization of an RF coil for MRI applications - Theoretical considerations - Abstract

A Combined Genetic Algorithm and Method of Moments design methods is presented for the design of unusual near-field antennas for use in Magnetic Resonance Imaging systems. The method is successfully applied to the design of an asymmetric coil structure for use at 190MHz and demonstrates excellent radiofrequency field homogeneity.