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.

Structural analysis and molecular model of a self-incompatibility RNase from wild tomato

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S-3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S-3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S-3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S-3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

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.

Calcium is essential for the structural integrity of the cysteine-rich, ligand-binding repeat of the low-density lipoprotein receptor

Seven cysteine-rich repeats form the ligand-binding region of the low-density lipoprotein (LDL) receptor. Each of these repeats is assumed to bind a calcium ion, which is needed for association of the receptor with its ligands, LDL and beta-VLDL. The effects of metal ions on the folding of the reduced N-terminal cysteine-rich repeat have been examined by using reverse-phase high-performance liquid chromatography to follow the formation of fully oxidized isomers with different disulfide connectivities. in the absence of calcium many of the 15 possible isomers formed on oxidation, whereas in its presence the predominant product at equilibrium had the native disulfide bond connectivities. Other metals were far less effective at directing disulfide bond formation: Mn2+ partly mimicked the action of Ca2+, but Ba2+, Sr2+, and Mg2+ had little effect. This metal-ion specificity was also observed in two-dimensional H-1 NMR spectral studies: only Ca2+ induced the native three-dimensional fold. The two paramagnetic ions, Gd3+ and Mn2+, and Cd2+ did not promote adoption of a well-defined structure, and the two paramagnetic ions did not displace calcium ions. The location of calcium ion binding sites in the repeat was also explored by NMR spectroscopy. The absence of chemical shift changes for the side chain proton resonances of Asp26, Asp36, and Glu37 from pH 3.9 to 6.8 in the presence of calcium ions and their proximal location in the NMR structures implicated these side chains as calcium ligands. Deuterium exchange NMR experiments also revealed a network of hydrogen bonds that stabilizes the putative calcium-binding loop.

Structural and functional characterisation of human sulfotransferases

The human aryl sulfotransferases HAST4 and HAST4v vary by only two amino acids but exhibit markedly different affinity towards the sulfonate acceptor p-nitrophenol and the sulfonate donor 3'-phosphoadenosine-5'-phosphosulfate (PAPS). To determine the importance of each of these amino acid differences, chimeric constructs were made of HAST4 and HAST4v. By attaching the last 120 amino acids of HAST-4v to HAST4 (changing Thr235 to Asn235) we have been able to produce a protein that has a K-m for PAPS similar to HAST4v. The reverse construct, HAST4v/4 produces a protein with a K-m for PAPS similar to HAST4. These data suggests that the COOH-terminal of sulfotransferases is involved in co-factor binding. (C) 1998 Elsevier Science Ireland Ltd. All rights reserved.

Functional and structural characterization of isolated perfused stingray liver including effects of ischaemia/reperfusion

The morphological and functional characteristics of stingray liver were studied, including the effect of ischaemia/reperfusion. With an isolated perfused model, it was shown that the stingray liver was more resistant than the rat liver to ischaemia/reperfusion injury; this was consistent with the differing partial oxygen tensions usually present in the two species. This study confirmed that whereas stingray hepatocytes form tubules with central bile canaliculi as in other fish, the stingray liver has portal triads and a lobular architecture as in mammals. Apoptosis of hepatocytes, demonstrated in the normal liver, was only marginally enhanced by ischaemia/reperfusion. Resulting apoptotic bodies were phagocytized by macrophage-like cells in hepatocyte tubules. In contrast to rat liver, the stingray liver showed no necrosis after ischaemia-reperfusion. (C) 1998 W.B. Saunders Company Limited.

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.

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.

Structure determination of the three disulfide bond isomers of alpha-conotoxin GI: A model for the role of disulfide bonds in structural stability

The three possible disulfide bonded isomers of alpha-conotoxin GI have been selectively synthesised and their structures determined by H-1 NMR spectroscopy. alpha-Conotoxin GI derives from the venom of Conus geographus and is a useful neuropharmacological tool as it selectively binds to the nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel involved in nerve signal transmission. The peptide has the sequence ECCNPACGRHYSC-NH2, and the three disulfide bonded isomers are referred to as GI(2-7;3-13), GI(2-13;3-7) and GI(2-3;7-13). The NMR structure for the native isomer GI(2-7;3-13) is of excellent quality, with a backbone pairwise RMSD of 0.16 Angstrom for a family of 35 structures, and comprises primarily a distorted 3(10),, helix between residues 5 to 11. The two non-native isomers exhibit multiple conformers in solution, with the major populated forms being different in structure both from each other and from the native form. Structure-activity relationships for the native GI(2-7;3-13) as well as the role of the disulfide bonds on folding and stability of the three isomers are examined. It is concluded that the disulfide bonds in alpha-conotoxin GI play a crucial part in determining both the structure and stability of the peptide. A trend for increased conformational heterogeneity was observed in the order of GI(2-7;3-13) < GI(2-13;3-7) < GI(2-3;7-13). It was found that the peptide bond joining Cys2 to Cys3 in GI(2-3;7-13) is predominantly trans, rather than cis as theoretically predicted. These structural data are used to interpret the varying nAChR binding of the non-native forms. A model for the binding of native GI(2-7;3-13) to the mammalian nAChR is proposed, with an alpha-subunit binding face made up of Cys2, Asn4, Pro5, Ala6 and Cys7 and a selectivity face, comprised of Arg9 and His10. These two faces orient the molecule between the alpha and delta subunits of the receptor. The structure of the CCNPAC sequence of the native GI(2-7;3-13) is compared to the structure of the identical sequence from the toxic domain of heat-stable enterotoxins, which forms part of the receptor binding region of the enterotoxins, but which has a different disulfide connectivity. (C) 1998 Academic Press Limited.

Stage-specific proteophosphoglycan from Leishmania mexicana amastigotes - Structural characterization of novel mono-, di-, and triphosphorylated phosphodiester-linked oligosaccharides

Intracellular amastigotes of the protozoan parasite Leishmania mexicana secrete a macromolecular proteophosphoglycan (aPPG) into the phagolysosome of their host cell, the mammalian macrophage. The structures of aPPG glycans were analyzed by a combination of high pH anion exchange high pressure liquid chromatography, gas chromatography-mass spectrometry, enzymatic digestions, electrospray-mass spectrometry as well as H-1 and P-31 NMR spectroscopy. Some glycans are identical to oligosaccharides known from Leishmania mexicana promastigote lipophosphoglycan and secreted acid phosphatase, However, the majority of the aPPG glycans represent amastigote stage-specific and novel structures. These include neutral glycans ([Glc beta(1-3)](1-2)Gal beta 1-4Man, Gal beta 1-3Gal beta 1-4Man, Gal beta 1-3Glc beta 1-3Gal beta 1-4Man), several monophosphorylated glycans containing the conserved phosphodisaccharide backbone (R-3-[PO4-6-Gal]beta 1-4Man) but carrying stage-specific modifications (R = Gal beta 1-, [Glc beta 1-3](1-2)Glc beta 1-), and monophosphorylated aPPG tri- and tetrasaccharides that are uniquely phosphorylated on the terminal hexose (PO4-6-Glc beta 1-3Gal beta 1-4Man, PO4-6-Glc beta 1-3Glc beta 1-3Gal beta 1-4Man, PO4-6-Gal beta 1-3Glc beta 1-3Gal beta 1-4Man), In addition aPPG contains highly unusual di- and triphosphorylated glycans whose major species are PO4-6-Glc beta 1-3Glc beta 1-3[PO4-6-Gal]beta 1-4Man, PO4-6-Gal beta 1-3Glc beta 1-3 [PO4-6-Gal]beta 1-4Man, PO4-6-GaL beta 1-3Glc beta 1-3Glc beta 1-3[PO4-6-Gal]beta 1-4Man, PO4-6-Glc beta 1-3[PO4-6-Glc]beta 1-3[PO4-6-Gal]beta 1-4Man, PO4-6Gal beta 1-3[PO4-6-Glc]beta 1-3Glc beta 1-3[PO4-6-Gal]beta 1-4Man, and PO4-6-Glc beta 1-3[PO4-6-Glc]beta 1-3Glc beta 1-3[PO4-6-Gal]beta 1-4Man. These glycans are linked together by the conserved phosphodiester R-Man alpha 1-PO4-6-Gal-R or the novel phosphodiester R-Man alpha 1-PO4-6-Glc-R and are connected to Ser(P) of the protein backbone most likely via the linkage R-Man alpha 1-PO4-Ser. The variety of stage-specific glycan structures in Leishmania mexicana aPPG suggests the presence of developmentally regulated amastigote glycosyltransferases which may be potential anti-parasite drug targets.

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.

A compact superconducting magnet for magnetic resonance microscopy

Magnetic resonance microscopy (MRM) depends on the use of high field, superconducting magnet systems for its operation. The magnets that are conventionally used are those that were initially designed for chemical structural analysis work. A novel, compact magnet designed specifically for MRM is presented here, and while preserving high field, high homogeneity conditions, has a length less than one-third that of conventional systems. This enables much better access to samples, an important consideration in many MRM experiments. As the homogeneity of a magnet is strongly dependent on its length, novel geometries and optimization techniques are required to meet the requirements of MRM in a compact system. An important outcome of the stochastic optimization performed in this work, is that the use used of a thin superconducting solenoid surrounded by counterwound disk windings provides a mechanism for drastic length reductions over conventional magnet designs. (C) 1998 American Institute of Physics.

Iron(III) and iron(II) complexes of 1-thia-4,7-diazacyclononane ([9]aneN(2)S) and 1,4-dithia-7-azacyclononane ([9]aneNS(2)). X-Ray structural analyses, magnetic susceptibility, Mossbauer, EPR and electronic spectroscopy

The complexes [Fe([9]aneN(2)S)(2)][ClO4](2), [Fe([9]aneN(2)S)(2)][ClO4](3) and [Fe([9]aneNS(2))(2)][ClO4](2) ([9]aneN(2)S = 1-thia-4. 7-diazacyclononane and [9]aneNS(2) = 1,4-dithia-7-azacyclononane) have been prepared and the latter two characterised by X-ray crystallography. The Mossbauer spectra (isomer shift/mm s(-1), quadrupole splitting/mm s(-1), 4.2 K) for [Fe([9]aneN(2)S)(2)][ClO4](2) (0.52, 0.57), [Fe([9]aneN(2)S)(2)][ClO4](3) (0.25, 2.72) and [Fe([9]aneNS(2))(2)][ClO4](2) (0.43, 0.28) are typical for iron(II) and iron(III) complexes. Variable-temperature susceptibility measurements for [Fe([9]aneN(2)S)(2)][ClO4](2) (2-300 K) revealed temperature-dependent behaviour in both the solid state [2.95 mu(B) (300 K)-0.5 mu(B) (4.2 K)] and solution (Delta H degrees 20-22 kJ mol(-1), Delta S degrees 53-60 J mol(-1) K-1). For [Fe([9]aneN(2)S)(2)][ClO4](3) in the solid state [2.3 mu(B) (300 K)-1.9 mu(B) (4.2 K)] the magnetic data were fit to a simple model (H = -lambda L . S + mu L-z) to give the spin-orbit coupling constant (lambda) of -260 +/- 10 cm(-1). The solid-state X-band EPR spectrum of [Fe([9]aneN(2)S)(2)][ClO4](3) revealed axial symmetry (g(perpendicular to) = 2.607, g(parallel to) = 1.599). Resolution of g(perpendicular to) into two components at Q-band frequencies indicated a rhombic distortion. The low-temperature single-crystal absorption spectra of [Fe([9]aneN(2)S)(2)][ClO4](2) and [Fe([9]aneNS(2))(2)][ClO4](2) exhibited additional bands which resembled pseudotetragonal low-symmetry splitting of the parent octahedral (1)A(1g) --> T-1(2g) and (1)A(1g) ---> T-1(1g) transitions. However, the magnitude of these splittings was too large, requiring 10Dq for the thioether donors to be significantly larger than for the amine donors. Instead, these bands were tentatively assigned to weak, low-energy S --> Fe-II charge-transfer transitions. Above 200 K, thermal occupation of the high-spin T-5(2g) ground state resulted in observation of the T-5(2g) --> E-5(g) transition in the crystal spectrum of [Fe([9]aneN(2)S)(2)][ClO4](2). From a temperature-dependence study, the separation of the low-spin (1)A(1g) and high-spin T-5(2g) ground states was approximately 1700 cm(-1). The spectrum of the iron(III) complex [Fe([9]aneN(2)S)(2)][ClO4](3) is consistent with a low-spin d(5) configuration.

Folding, calcium binding, and structural characterization of a concatemer of the first and second ligand-binding modules of the low-density lipoprotein receptor

The ligand-binding domain of the low-density lipoprotein (LDL) receptor is comprised of seven tandemly repeated ligand-binding modules, each being approximately 40 amino acids long and containing six conserved cysteine residues. We have expressed and characterized a concatemer of the first two modules (LB1 and LB2) of the human LDL receptor. Oxidative folding of the recombinant concatemer (rLB(1-2)), in the presence of calcium ions, gave a single dominant isomer with six disulfide bonds. Peptic cleavage of the short Linker region that connects the last cysteine residue of LB1 and the first cysteine residue of LB2 yielded two discrete fragments, thus excluding the presence of intermodule disulfide bonds. The N-terminal module, LB1, reacted with a conformation-specific monoclonal antibody (IgG-C7) made to LB1 in the native LDL receptor. From this, we concluded that the first module was correctly folded, with the same set of disulfide bonds as LB1 of the LDL receptor. The disulfide bond connections of LB2 were identified from mass spectral analysis of fragments formed by digestion of the C-terminal peptic fragment with elastase. These data showed that the disulfide bonds of LB2 connected Cys(I) and Cys(III), Cys(II) and Cys(V), and Cys(IV) and Cys(VI). This pattern is identical to that found for recombinant LB1 and LB2. The concatemer has two high-affinity calcium-binding sites, one per module. An analysis of the secondary chemical shifts of C alpha protons shows that the conformations of LB1 and LB2 in the concatemer are very similar to those of the individual modules, with no evidence for strong interactions between the two modules.

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.