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.

Seven retinal specializations in the tubular eye of the deep-sea pearleye, Scopelarchus michaelsarsi: A case study in visual optimization

The deep-sea pearleye, Scopelarchus michaelsarsi (Scopelarchidae) is a mesopelagic teleost with asymmetric or tubular eyes. The main retina subtends a large dorsal binocular field, while the accessory retina subtends a restricted monocular field of lateral visual space. Ocular specializations to increase the lateral visual field include an oblique pupil and a corneal lens pad. A detailed morphological and topographic study of the photoreceptors and retinal ganglion cells reveals seven specializations: a centronasal region of the main retina with ungrouped rod-like photoreceptors overlying a retinal tapetum; a region of high ganglion cell density (area centralis of 56.1x10(3) cells per mm(2)) in the centrolateral region of the main retina; a centrotemporal region of the main retina with grouped rod-like photoreceptors; a region (area giganto cellularis) of large (32.2+/-5.6 mu m(2)), alpha-like ganglion cells arranged in a regular array (nearest neighbour distance 53.5+/-9.3 mu m with a conformity ratio of 5.8) in the temporal main retina; an accessory retina with grouped rod-like photoreceptors; a nasotemporal band of a mixture of rod-and cone-like photoreceptors restricted to the ventral accessory retina; and a retinal diverticulum comprised of a ventral region of differentiated accessory retina located medial to the optic nerve head. Retrograde labelling from the optic nerve with DiI shows that approximately 14% of the cells in the ganglion cell layer of the main retina are displaced amacrine cells at 1.5 mm eccentricity. Cryosectioning of the tubular eye confirms Matthiessen's ratio (2.59), and calculations of the spatial resolving power suggests that the function of the area centralis (7.4 cycles per degree/8.1 minutes of are) and the cohort of temporal alpha-like ganglion cells (0.85 cycles per degree/70.6 minutes of are) in the main retina may be different. Low summation ratios in these various retinal zones suggests that each zone may mediate distinct visual tasks in a certain region of the visual field by optimizing sensitivity and/or resolving power.

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.

Small molecular probes for G-protein-coupled C5a receptors: Conformationally constrained antagonists derived from the C terminus of the human plasma protein C5a

Activation of the human complement system of plasma proteins in response to infection or injury produces a 4-helix bundle glycoprotein (74 amino acids) known as C5a. C5a binds to G-protein-coupled receptors on cell surfaces triggering receptor-ligand internalization, signal transduction, and powerful inflammatory responses. Since excessive levels of C5a are associated with autoimmune and chronic inflammatory disorders, inhibitors of receptor activation may have therapeutic potential. We now report solution structures and receptor-binding and antagonist activities for some of the first small molecule antagonists of C5a derived from its hexapeptide C terminus. The antagonist NMe-Phe-Lys-Pro-D-Cha-Trp-D-Arg-CO2H (1) surprisingly shows an unusually well-defined solution structure as determined by H-1 NMR spectroscopy. This is one of the smallest acyclic peptides found to possess a defined solution conformation, which can be explained by the constraining role of intramolecular hydrogen bonding. NOE and coupling constant data, slow deuterium exchange, and a low dependence on temperature for the chemical shift of the D-Cha-NH strongly indicate an inverse gamma turn stabilized by a D-Cha-NH ... OC-Lys hydrogen bond. Smaller conformational populations are associated with a hydrogen bond between Trp-NH ... OC-Lys, defining a type II beta turn distorted by the inverse gamma turn incorporated within it. An excellent correlation between receptor-affinity and antagonist activity is indicated for a limited set of synthetic peptides. Conversion of the C-terminal carboxylate of 1 to an amide decreases antagonist potency 5-fold, but potency is increased up to 10-fold over 1 if the amide bond is made between the C-terminal carboxylate and a Lys/Orn side chain to form a cyclic analogue. The solution structure of cycle 6 also shows gamma and beta turns; however, the latter occurs in a different position, and there are clear conformational changes in 6 vs 1 that result in enhanced activity. These results indicate that potent C5a antagonists can be developed by targeting site 2 alone of the C5a receptor and define a novel pharmacophore for developing powerful receptor probes or drug candidates.

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.

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.

Analysis of small signal stability margins using genetic optimization

Power system small signal stability analysis aims to explore different small signal stability conditions and controls, namely: (1) exploring the power system security domains and boundaries in the space of power system parameters of interest, including load flow feasibility, saddle node and Hopf bifurcation ones; (2) finding the maximum and minimum damping conditions; and (3) determining control actions to provide and increase small signal stability. These problems are presented in this paper as different modifications of a general optimization to a minimum/maximum, depending on the initial guesses of variables and numerical methods used. In the considered problems, all the extreme points are of interest. Additionally, there are difficulties with finding the derivatives of the objective functions with respect to parameters. Numerical computations of derivatives in traditional optimization procedures are time consuming. In this paper, we propose a new black-box genetic optimization technique for comprehensive small signal stability analysis, which can effectively cope with highly nonlinear objective functions with multiple minima and maxima, and derivatives that can not be expressed analytically. The optimization result can then be used to provide such important information such as system optimal control decision making, assessment of the maximum network's transmission capacity, etc. (C) 1998 Elsevier Science S.A. All rights reserved.

Simulation and optimization of heated, inviscid flows in scramjet ducts

A space-marching code for the simulation and optimization of inviscid supersonic flow in three dimensions is described. The now in a scramjet module with a relatively complex three-dimensional geometry is examined and wall-pressure estimates are compared with experimental data. Given that viscous effects are not presently included, the comparison is reasonable. The thermodynamic compromise of adding heat in a diverging combustor is also examined. The code is then used to optimize the shape of a thrust surface for a simpler (box-section) scramjet module in the presence of uniform and nonuniform heat distributions. The optimum two-dimensional profiles for the thrust surface are obtained via a perturbation procedure that requires about 30-50 now solutions. It is found that the final shapes are fairly insensitive to the details of the heat distribution.

A consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements

To translate and transfer solution data between two totally different meshes (i.e. mesh 1 and mesh 2), a consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements is presented in this paper. The proposed algorithm has the following significant advantages: (1) The use of a point-searching strategy allows a point in one mesh to be accurately related to an element (containing this point) in another mesh. Thus, to translate/transfer the solution of any particular point from mesh 2 td mesh 1, only one element in mesh 2 needs to be inversely mapped. This certainly minimizes the number of elements, to which the inverse mapping is applied. In this regard, the present algorithm is very effective and efficient. (2) Analytical solutions to the local co ordinates of any point in a four-node quadrilateral element, which are derived in a rigorous mathematical manner in the context of this paper, make it possible to carry out an inverse mapping process very effectively and efficiently. (3) The use of consistent interpolation enables the interpolated solution to be compatible with an original solution and, therefore guarantees the interpolated solution of extremely high accuracy. After the mathematical formulations of the algorithm are presented, the algorithm is tested and validated through a challenging problem. The related results from the test problem have demonstrated the generality, accuracy, effectiveness, efficiency and robustness of the proposed consistent point-searching algorithm. Copyright (C) 1999 John Wiley & Sons, Ltd.

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.

Assisting decision-making in Queensland barley production through chance constrained programming

A chance constrained programming model is developed to assist Queensland barley growers make varietal and agronomic decisions in the face of changing product demands and volatile production conditions. Unsuitable or overlooked in many risk programming applications, the chance constrained programming approach nonetheless aptly captures the single-stage decision problem faced by barley growers of whether to plant lower-yielding but potentially higher-priced malting varieties, given a particular expectation of meeting malting grade standards. Different expectations greatly affect the optimal mix of malting and feed barley activities. The analysis highlights the suitability of chance constrained programming to this specific class of farm decision problem.

Modelling life history strategies with capture-recapture data: Evolutionary demography of the water skink Eulamprus tympanum

Matrix population models, elasticity analysis and loop analysis can potentially provide powerful techniques for the analysis of life histories. Data from a capture-recapture study on a population of southern highland water skinks (Eulamprus tympanum) were used to construct a matrix population model. Errors in elasticities were calculated by using the parametric bootstrap technique. Elasticity and loop analyses were then conducted to identify the life history stages most important to fitness. The same techniques were used to investigate the relative importance of fast versus slow growth, and rapid versus delayed reproduction. Mature water skinks were long-lived, but there was high immature mortality. The most sensitive life history stage was the subadult stage. It is suggested that life history evolution in E. tympanum may be strongly affected by predation, particularly by birds. Because our population declined over the study, slow growth and delayed reproduction were the optimal life history strategies over this period. Although the techniques of evolutionary demography provide a powerful approach for the analysis of life histories, there are formidable logistical obstacles in gathering enough high-quality data for robust estimates of the critical parameters.