Fast evolution strategies

Evolution strategies are a class of general optimisation algorithms which are applicable to functions that are multimodal, nondifferentiable, or even discontinuous. Although recombination operators have been introduced into evolution strategies, the primary search operator is still mutation. Classical evolution strategies rely on Gaussian mutations. A new mutation operator based on the Cauchy distribution is proposed in this paper. It is shown empirically that the new evolution strategy based on Cauchy mutation outperforms the classical evolution strategy on most of the 23 benchmark problems tested in this paper. The paper also shows empirically that changing the order of mutating the objective variables and mutating the strategy parameters does not alter the previous conclusion significantly, and that Cauchy mutations with different scaling parameters still outperform the Gaussian mutation with self-adaptation. However, the advantage of Cauchy mutations disappears when recombination is used in evolution strategies. It is argued that the search step size plays an important role in determining evolution strategies' performance. The large step size of recombination plays a similar role as Cauchy mutation.

Evolutionary distinctiveness and status of the endangered Lake Eacham rainbowfish (Melanotaenia eachamensis)

The Lake Eacham rainbowfish (Melanotaenia eachamensis) was declared extinct in the wild in the late 1980s after it disappeared from its only known locality, an isolated crater lake in northeast Queensland. Doubts have been raised about whether this taxon is distinct from surrounding populations of the eastern rainbowfish (Melanotaenia splendida splendida). We examined the evolutionary distinctiveness of M. eachamensis, obtained from captive stocks, relative to M. s. splendida through analysis of variation in mtDNA sequences, nuclear microsatellites, and morphometric characters Captive M. eachamensis had mtDNAs that were highly divergent from those in most populations of M. s. splendida. A broader geographic survey using RFLPs revealed some populations initially identified as M. s. splendida, that carried eachamensis mtDNA, whereas some others had mixtures of eachamensis and splendida mtDNA. The presence of eachamensis-like mtDNA in these populations could in principle be due to (1) sorting of ancestral polymorphisms, (2) introgression of M. eachamensis mtDNA into M. s. splendida, or (3) incorrect species boundaries, such that some populations currently assigned to M. s. splendida are M. eachamensis or are mixtures of the two species. These alternatives hypotheses were evaluated through comparisons of four nuclear microsatellite loci and morphometrics and meristics. In analyses of both data sets, populations of M. s. splendida with eachamensis mtDNA were more similar to captive M. eachamensis than to M. s. splendida with splendida mtDNA, supporting hypothesis 3. These results are significant for the management of M. eachamensis in several respects. First the combined molecular and morphological evidence indicates that M. eachamensis is a distinct species and a discrete evolutionarily significant unit worthy of conservation effort. Second it appears that the species boundary between M. eachamensis and M. s. splendida has been misdiagnosed such that there are extant populations on the Atherton Tableland as well as areas where both forms coexist. Accordingly we suggest that M. eachamensis be listed as vulnerable, rather than critical (or extinct in the wild). Third, the discovery of extant but genetically divergent populations of M. eachamensis on the Atherton Tableland broadens the options for future reintroductions to Lake Eacham.

Prediction of MHC class II-binding peptides using an evolutionary algorithm and artificial neural network

Motivation: Prediction methods for identifying binding peptides could minimize the number of peptides required to be synthesized and assayed, and thereby facilitate the identification of potential T-cell epitopes. We developed a bioinformatic method for the prediction of peptide binding to MHC class II molecules. Results: Experimental binding data and expert knowledge of anchor positions and binding motifs were combined with an evolutionary algorithm (EA) and an artificial neural network (ANN): binding data extraction --> peptide alignment --> ANN training and classification. This method, termed PERUN, was implemented for the prediction of peptides that bind to HLA-DR4(B1*0401). The respective positive predictive values of PERUN predictions of high-, moderate-, low- and zero-affinity binder-a were assessed as 0.8, 0.7, 0.5 and 0.8 by cross-validation, and 1.0, 0.8, 0.3 and 0.7 by experimental binding. This illustrates the synergy between experimentation and computer modeling, and its application to the identification of potential immunotheraaeutic peptides.

Cost functions and duality for stochastic technologies

Cost functions dual to stochastic production technologies are derived and their properties are discussed. These cost functions are shown to be consistent with expected-utility maximization without placing serious structural restrictions on the underlying technology.

Studies of evolutionary temperature adaptation: Muscle function and locomotor performance in Antarctic fish

1, Studies of evolutionary temperature adaptation of muscle and locomotor performance in fish are reviewed with a focus on the Antarctic fauna living at subzero temperatures. 2. Only limited data are available to compare the sustained and burst swimming kinematics and performance of Antarctic, temperate and tropical species. Available data indicate that low temperatures limit maximum swimming performance and this is especially evident in fish larvae. 3, In a recent study, muscle performance in the Antarctic rock cod Notothenia coriiceps at 0 degrees C was found to be sufficient to produce maximum velocities during burst swimming that were similar to those seen in the sculpin Myoxocephalus scorpius at 10 degrees C, indicating temperature compensation of muscle and locomotor performance in the Antarctic fish. However, at 15 degrees C, sculpin produce maximum swimming velocities greater than N, coriiceps at 0 degrees C, 4, It is recommended that strict hypothesis-driven investigations using ecologically relevant measures of performance are undertaken to study temperature adaptation in Antarctic fish, Recent detailed phylogenetic analyses of the Antarctic fish fauna and their temperate relatives will allow a stronger experimental approach by helping to separate what is due to adaptation to the cold and what is due to phylogeny alone.

Nine classes of integrable boundary conditions for the eight-state supersymmetric fermion model

Nine classes of integrable boundary conditions for the eight-state supersymmetric model of strongly correlated fermions are presented. The boundary systems are solved by using the coordinate Bethe ansatz method and the Bethe ansatz equations for all nine cases are given.

The evolutionary structure of scientific theories

David Hull's (1988c) model of science as a selection process suffers from a two-fold inability: (a) to ascertain when a lineage of theories has been established; i.e., when theories are descendants of older theories or are novelties, and what counts as a distinct lineage; and (b) to specify what the scientific analogue is of genotype and phenotype. This paper seeks to clarify these issues and to propose an abstract model of theories analogous to particulate genetic structure, in order to reconstruct relationships of descent and identity.

The dimerization and topological specificity functions of MinE reside in a structurally autonomous C-terminal domain

Correct placement of the division septum in Escherichia coli requires the co-ordinated action of three proteins, MinC, MinD and MinE. MinC and MinD interact to form a non-specific division inhibitor that blocks septation at all potential division sites. MinE is able to antagonize MinCD in a topologically sensitive manner, as it restricts MinCD activity to the unwanted division sites at the cell poles, Here, we show that the topological specificity function of MinE residues in a structurally autonomous, trypsin-resistant domain comprising residues 31-88, Nuclear magnetic resonance (NMR) and circular dichroic spectroscopy indicate that this domain includes both alpha and beta secondary structure, while analytical ultracentrifugation reveals that it also contains a region responsible for MinE homodimerization. While trypsin digestion indicates that the anti-MinCD domain of MinE (residues 1-22) does not form a tightly folded structural domain, NMR analysis of a peptide corresponding to MinE(1-22) indicates that this region forms a nascent helix in which the peptide rapidly interconverts between disordered (random coil) and alpha-helical conformations, This suggests that the N-terminal region of MinE may be poised to adopt an alpha-helical conformation when it interacts with the target of its anti-MinCD activity, presumably MinD.

Theory of modulational instability in Bragg gratings with quadratic nonlinearity

Modulational instability in optical Bragg gratings with a quadratic nonlinearity is studied. The electric field in such structures consists of forward and backward propagating components at the fundamental frequency and its second harmonic. Analytic continuous wave (CW) solutions are obtained, and the intricate complexity of their stability, due to the large number of equations and number of free parameters, is revealed. The stability boundaries are rich in structures and often cannot be described by a simple relationship. In most cases, the CW solutions are unstable. However, stable regions are found in the nonlinear Schrodinger equation limit, and also when the grating strength for the second harmonic is stronger than that of the first harmonic. Stable CW solutions usually require a low intensity. The analysis is confirmed by directly simulating the governing equations. The stable regions found have possible applications in second-harmonic generation and dark solitons, while the unstable regions maybe useful in the generation of ultrafast pulse trains at relatively low intensities. [S1063-651X(99)03005-6].

Protein fold recognition score functions: Unusual construction strategies

We describe two ways of optimizing score functions for protein sequence to structure threading. The first method adjusts parameters to improve sequence to structure alignment. The second adjusts parameters so as to improve a score function's ability to rank alignments calculated in the first score function. Unlike those functions known as knowledge-based force fields, the resulting parameter sets do not rely on Boltzmann statistics, have no claim to representing free energies and are purely constructions for recognizing protein folds. The methods give a small improvement, but suggest that functions can be profitably optimized for very specific aspects of protein fold recognition, Proteins 1999;36:454-461. (C) 1999 Wiley-Liss, Inc.

Evolutionary conservation of the membrane fusion machine

Recent structural studies of proteins mediating membrane fusion reveal intriguing similarities between diverse viral and mammalian systems. Particularly striking is the close similarity between the transmembrane envelope glycoproteins from the retrovirus HTLV-1 and the filovirus Ebola. These similarities suggest similar mechanisms of membrane fusion. The model that fits most currently available data suggests fusion activation in viral systems is driven by a symmetrical conformational change triggered by an activation event such as receptor binding or a pH change. The mammalian vesicle fusion mediated by the SNARE protein complex most likely occurs by a similar mechanism but without symmetry constraints.

Myeloproliferative disorder and leukaemia in mice induced by different classes of constitutive mutants of the human IL-3/IL-5/GM-CSF receptor common beta subunit

Several constitutively active mutant forms of the common β subunit of the human IL-3, IL-5 and GM-CSF receptors (hβc), which enable it to signal in the absence of ligand, have recently been described. Two of these, V449E and I374N, are amino acid substitutions in the transmembrane and extracellular regions of hβc, respectively. A third, FIΔ, contains a 37 amino acid duplication in the extracellular domain. We have shown previously that when expressed in primary murine haemopoietic cells, the extracellular mutants confer factor-independence on cells of the neutrophil and monocyte lineages only, whereas V449E does so on all cell types of the myeloid and erythroid compartments. To study the in vivo effects and leukaemic potential of these mutants, we have expressed all three in mice by bone marrow reconstitution using retrovirally infected donor cells. Expression of the extracellular mutants leads to an early onset, chronic myeloproliferative disorder marked by elevations in the neutrophil, monocyte, erythrocyte and platelet lineages. In contrast, expression of V449E leads to an acute leukaemia-like syndrome of anaemia, thrombocytopaenia and blast cell expansion. These data support the possibility that activating mutations in hβc are involved in haemopoietic disorders in man.

Murine cytomegalovirus homologues of cellular immunomodulatory genes

The study of 'molecular mimicry' or 'genetic piracy', with respect to the utilisation of cellular genes captured and modified during the course of virus evolution, has been an area of increasing research with the expansion in virus genome sequencing. Examples of cellular immunomodulatory genes which have been captured from hosts have been identified in a number of viruses. This review concentrates upon studies of murine cytomegalovirus (MCMV), investigating the functions of viral genes homologous to G protein-coupled receptors, MHC class I and chemokines, The study of recombinant MCMV engineered with specific disruptions of these genes has revealed their significance during virus replication and dissemination within the host, In the case of the latter two classes of genes, evidence suggests they interfere with cellular immune responses, although the detailed mechanisms underlying this interference have yet to be delineated. Copyright (C) 2000 S. Karger AG, Basel.

Computation of bifurcation boundaries for power systems: A new Delta-plane method

This paper is devoted to the problems of finding the load flow feasibility, saddle node, and Hopf bifurcation boundaries in the space of power system parameters. The first part contains a review of the existing relevant approaches including not-so-well-known contributions from Russia. The second part presents a new robust method for finding the power system load flow feasibility boundary on the plane defined by any three vectors of dependent variables (nodal voltages), called the Delta plane. The method exploits some quadratic and linear properties of the load now equations and state matrices written in rectangular coordinates. An advantage of the method is that it does not require an iterative solution of nonlinear equations (except the eigenvalue problem). In addition to benefits for visualization, the method is a useful tool for topological studies of power system multiple solution structures and stability domains. Although the power system application is developed, the method can be equally efficient for any quadratic algebraic problem.

U-q[(sl(2)over-cap vertical bar 1)] vertex operators, screen currents, and correlation functions at an arbitrary level

Bosonized q-vertex operators related to the four-dimensional evaluation modules of the quantum affine superalgebra U-q[sl((2) over cap\1)] are constructed for arbitrary level k=alpha, where alpha not equal 0,-1 is a complex parameter appearing in the four-dimensional evaluation representations. They are intertwiners among the level-alpha highest weight Fock-Wakimoto modules. Screen currents which commute with the action of U-q[sl((2) over cap/1)] up to total differences are presented. Integral formulas for N-point functions of type I and type II q-vertex operators are proposed. (C) 2000 American Institute of Physics. [S0022-2488(00)00608-3].