Structure optimization combining soft-core interaction functions, the diffusion equation method, and molecular dynamics

Smoothing the potential energy surface for structure optimization is a general and commonly applied strategy. We propose a combination of soft-core potential energy functions and a variation of the diffusion equation method to smooth potential energy surfaces, which is applicable to complex systems such as protein structures; The performance of the method was demonstrated by comparison with simulated annealing using the refinement of the undecapeptide Cyclosporin A as a test case. Simulations were repeated many times using different initial conditions and structures since the methods are heuristic and results are only meaningful in a statistical sense.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

Rapid computation of static fields produced by thick circular solenoids

A straightforward method is proposed for computing the magnetic field produced by a circular coil that contains a large number of turns wound onto a solenoid of rectangular cross section. The coil is thus approximated by a circular ring containing a continuous constant current density, which is very close to the real situation when sire of rectangular cross section is used. All that is required is to evaluate two functions, which are defined as integrals of periodic quantities; this is done accurately and efficiently using trapezoidal-rule quadrature. The solution can be obtained so rapidly that this procedure is ideally suited for use in stochastic optimization, An example is given, in which this approach is combined with a simulated annealing routine to optimize shielded profile coils for NMR.

Robust parameter settings of evolutionary algorithms for the optimisation of agricultural systems models

Numerical optimisation methods are being more commonly applied to agricultural systems models, to identify the most profitable management strategies. The available optimisation algorithms are reviewed and compared, with literature and our studies identifying evolutionary algorithms (including genetic algorithms) as superior in this regard to simulated annealing, tabu search, hill-climbing, and direct-search methods. Results of a complex beef property optimisation, using a real-value genetic algorithm, are presented. The relative contributions of the range of operational options and parameters of this method are discussed, and general recommendations listed to assist practitioners applying evolutionary algorithms to the solution of agricultural systems. (C) 2001 Elsevier Science Ltd. All rights reserved.

An ESR study of the gamma radiolysis of aromatic polyesters containing isomeric naphthalene links

Six polyesters were synthesised from 4,4 ' -oxy-bis(benzoyl chloride) and 1,4-, 1,5-, 1,6-, 2,3-, 2,6-, and 2,7-naphthalenediol isomers. The structures of the polyesters were characterised by means of IR, inherent viscosities in tetrachloroethane (TCE), solutions at 303 K and thermal analysis. The glass transition temperatures were in the range of 425-494 K by DSC thermal analysis. All of the polyesters were irradiated in an AECL Gammacell 220 unit at a dose rate of approximately 6.7 kGy/h to doses in the range of 0-15 kGy at 77 and 300 K. ESR spectroscopy was used to examine the radicals formed during radiolysis and to measure their yields. The G-values for radical formation in the polyesters were found to be in the range 0.18-1.41 at 77 K and 0.19-0.78 at 300 K. At 77 K, up to 15% of the radicals formed on radiolysis were found to be photo-bleachable anion radicals. Annealing experiments were carried out in order to identify the neutral radicals, which were assigned to naphthyl- or phenyl- and phenoxyl-type radicals. (C) 2001 Elsevier Science Ltd. All rights reserved.

An ESR study on gamma-irradiated poly(vinyl alcohol)

The formation of radicals in poly(vinyl alcohol), PVA, powder irradiated at 77 K by gamma -rays and the transformations of these radicals during photolysis with visible wavelengths and on thermal annealing have been studied. After irradiation a four-line ESR spectrum was observed. It was assigned to a triplet of the C-alpha-radical (38%), with a splitting of 3.27 mT, superimposed on a doublet (62%) with a splitting of 2.7 mT. The doublet appears to be composed of two radicals, one of which is photo-bleachable (58%) and the other which is not photo-bleachable (42%). This suggests that the latter radical is a neutral radical. The photo-bleachable component of the doublet has been assigned to a carbonyl anion radical. but the second doublet due to a neutral radical is unassigned. The total G-value for formation of radicals at 77 K was found to be 2.41 +/- 0.03. Upon illumination with visible light, the anion radicals were removed and the doublet components or the spectrum diminished in intensity, while the three-line spectrum of the C-alpha-radical became more clearly visible. This transition was due to the photo-detachment of electrons from traps which were proposed to be located on carbonyl groups in the polymer resulting from incomplete hydrolysis of the vinyl acetate. The photo-decay of the anion radicals could be satisfactorily described by a two-stage process. The first stage comprised the decay of approximately 80% of the anion radicals present, while the second stage was associated with the decay of the remaining 20%. Subsequent thermal annealing of a photolysed sample to 290 K led to a change in the shape of the spectrum to form a more clearly defined triplet, As the doublet of the neutral radical decays on thermal annealing between 150 and 250K, the C-alpha-radical is formed. (C) 2001 Elsevier Science Ltd. All rights reserved.

Three-dimensional structure of RTD-1, a cyclic antimicrobial defensin from rhesus macaque leukocytes

Most mammalian defensins are cationic peptides of 29-42 amino acids long, stabilized by three disulfide bonds. However, recently Tang et al. (1999, Science 286, 498-502) reported the isolation of a new defensin type found in the leukocytes of rhesus macaques. In contrast to all the other defensins found so far, rhesus theta defensin-1 (RTD-1) is composed of just 18 amino acids with the backbone cyclized through peptide bonds. Antibacterial activities of both the native cyclic peptide and a linear form were examined, showing that the cyclic form was 3-fold more active than the open chain analogue [Tang et al. (1999) Science 286, 498-502]. To elucidate the three-dimensional structure of RTD-1 and its open chain analogue, both peptides were synthesized using solid-phase peptide synthesis and tert-butyloxycarbonyl chemistry. The structures of both peptides in aqueous solution were determined from two-dimensional H-1 NMR data recorded at 500 and 750 MHz. Structural constraints consisting of interproton distances and dihedral angles were used as input for simulated-annealing calculations and water refinement with the program CNS. RTD-1 and its open chain analogue oRTD-1 adopt very similar structures in water. Both comprise an extended beta -hairpin structure with turns at one or both ends. The turns are well defined within themselves and seem to be flexible with respect to the extended regions of the molecules. Although the two strands of the beta -sheet are connected by three disulfide bonds, this region displays a degree of flexibility. The structural similarity of RTD-1 and its open chain analogue oRTD-1, as well as their comparable degree of flexibility, support the theory that the additional charges at the termini of the open chain analogue rather than overall differences in structure or flexibility are the cause for oRTD-1's lower antimicrobial activity. In contrast to numerous other antimicrobial peptides, RTD-1 does not display any amphiphilic character, even though surface models of RTD-1 exhibit a certain clustering of positive charges. Some amide protons of RTD-1 that should be solvent-exposed in monomeric beta -sheet structures show low-temperature coefficients, suggesting the possible presence of weak intermolecular hydrogen bonds.

Glass transition in Australian honeys

Differential scanning calorimetry (DSC) was used to study the glass transition in 10 Australian honeys by scanning at 10 degrees Cmin(-1) from -130 to 50 degreesC after annealing at -50 degreesC. The honeys had moisture contents 14.9 to 18.0%, seven were from Eucalyptus species. The glass transition temperatures (T-g) ranged from -46 degrees to -38 degreesC and were significantly (p

The stochastic design of force-minimized compact magnets for high-field magnetic resonance imaging applications

New designs for force-minimized compact high-field clinical MRI magnets are described. The design method is a modified simulated annealing (SA) procedure which includes Maxwell forces in the error function to be minimized. This permits an automated force reduction in the magnet designs while controlling the overall dimensions of the system. As SA optimization requires many iterations to achieve a final design, it is important that each iteration in the procedure is rapid. We have therefore developed a rapid force calculation algorithm. Novel designs for short 3- and 4-T clinical MRI systems are presented in which force reduction has been invoked. The final designs provide large homogeneous regions and reduced stray fields in remarkable short magnets. A shielded 4-T design that is approximately 30% shorter than current designs is presented. This novel magnet generates a full 50-cm diameter homogeneous region.

Mathematical methods for spatially cohesive reserve design

The problem of designing spatially cohesive nature reserve systems that meet biodiversity objectives is formulated as a nonlinear integer programming problem. The multiobjective function minimises a combination of boundary length, area and failed representation of the biological attributes we are trying to conserve. The task is to reserve a subset of sites that best meet this objective. We use data on the distribution of habitats in the Northern Territory, Australia, to show how simulated annealing and a greedy heuristic algorithm can be used to generate good solutions to such large reserve design problems, and to compare the effectiveness of these methods.

The radiation chemistry of ultem at 77 K as revealed by ESR spectroscopy

Radical formation in ultem following gamma-radiolysis has been reassessed, and the G(R*) values at different temperatures have been determined by ESR spectroscopy. The radical assignment and radical reactivity have been re-examined by photobleaching and thermal annealing studies. Photobleachable radical anions were found to comprise approximate to40% of the total number of radicals formed on radiolysis at 77 K. Spectral subtraction methods, ESR spectral simulations, measurement of g-values and the hyperfine splitting constants were used to identify the other radical intermediates. The principal chain scission radicals are formed due to scission of the main-chain at (i) the ether linkage, (ii) the isopropylidene group and (iii) the imide ring in the main chain. The side chain methyl groups of the isopropylidine units also lose hydrogen to form methylene radicals. The five-line spectrum observed to decay in the temperature range 370-430 K, which has not been assigned previously, has been identified as being characteristic of a di-substituted benzyl radical. (C) 2002 Elsevier Science Ltd. All rights reserved.

The Three-dimensional Solution Structure of NaD1, a New Floral Defensin from Nicotiana alata and its Application to a Homology Model of the Crop Defense Protein alfAFP

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of NaD1, a novel antifungal and insecticidal protein isolated from the flowers of Nicotiana alata. NaD1 is a basic, cysteine-rich protein of 47 residues and is the first example of a plant defensin from flowers to be characterized structurally. Its three-dimensional structure consists of an a-helix and a triple-stranded anti-parallel beta-sheet that are stabilized by four intramolecular disulfide bonds. NaD1 features all the characteristics of the cysteine-stabilized up motif that has been described for a variety of proteins of differing functions ranging from antibacterial insect defensins and ion channel-perturbing scorpion toxins to an elicitor of the sweet taste response. The protein is biologically active against insect pests, which makes it a potential candidate for use in crop protection. NaD1 shares 31% sequence identity with alfAFP, an antifungal protein from alfalfa that confers resistance to a fungal pathogen in transgenic potatoes. The structure of NaD1 was used to obtain a homology model of alfAFP, since NaD1 has the highest level of sequence identity with alfAFP of any structurally characterized antifungal defensin. The structures of NaD1 and alfAFP were used in conjunction with structure - activity data for the radish defensin Rs-AFP2 to provide an insight into structure-function relationships. In particular, a putative effector site was identified in the structure of NaD1 and in the corresponding homology model of alfAFP. (C) 2002 Elsevier Science Ltd. All rights reserved.

Using siting algorithms in the design of marine reserve networks

Using benthic habitat data from the Florida Keys (USA), we demonstrate how siting algorithms can help identify potential networks of marine reserves that comprehensively represent target habitat types. We applied a flexible optimization tool-simulated annealing-to represent a fixed proportion of different marine habitat types within a geographic area. We investigated the relative influence of spatial information, planning-unit size, detail of habitat classification, and magnitude of the overall conservation goal on the resulting network scenarios. With this method, we were able to identify many adequate reserve systems that met the conservation goals, e.g., representing at least 20% of each conservation target (i.e., habitat type) while fulfilling the overall aim of minimizing the system area and perimeter. One of the most useful types of information provided by this siting algorithm comes from an irreplaceability analysis, which is a count of the number of, times unique planning units were included in reserve system scenarios. This analysis indicated that many different combinations of sites produced networks that met the conservation goals. While individual 1-km(2) areas were fairly interchangeable, the irreplaceability analysis highlighted larger areas within the planning region that were chosen consistently to meet the goals incorporated into the algorithm. Additionally, we found that reserve systems designed with a high degree of spatial clustering tended to have considerably less perimeter and larger overall areas in reserve-a configuration that may be preferable particularly for sociopolitical reasons. This exercise illustrates the value of using the simulated annealing algorithm to help site marine reserves: the approach makes efficient use of;available resources, can be used interactively by conservation decision makers, and offers biologically suitable alternative networks from which an effective system of marine reserves can be crafted.

Evidence for a QTL on chromosome 19 influencing LDL cholesterol levels in the general population

The genetic basis of cardiovascular disease (CVD) with its complex etiology is still largely elusive. Plasma levels of lipids and apolipoproteins are among the major quantitative risk factors for CVD and are well-established intermediate traits that may be more accessible to genetic dissection than clinical CVD end points. Chromosome 19 harbors multiple genes that have been suggested to play a role in lipid metabolism and previous studies indicated the presence of a quantitative trait locus (QTL) for cholesterol levels in genetic isolates. To establish the relevance of genetic variation at chromosome 19 for plasma levels of lipids and apolipoproteins in the general, out-bred Caucasian population, we performed a linkage study in four independent samples, including adolescent Dutch twins and adult Dutch, Swedish and Australian twins totaling 493 dizygotic twin pairs. The average spacing of short-tandem-repeat markers was 6 - 8 cM. In the three adult twin samples, we found consistent evidence for linkage of chromosome 19 with LDL cholesterol levels ( maximum LOD scores of 4.5, 1.7 and 2.1 in the Dutch, Swedish and Australian sample, respectively); no indication for linkage was observed in the adolescent Dutch twin sample. The QTL effects in the three adult samples were not significantly different and a simultaneous analysis of the samples increased the maximum LOD score to 5.7 at 60 cM pter. Bivariate analyses indicated that the putative LDL-C QTL also contributed to the variance in ApoB levels, consistent with the high genetic correlation between these phenotypes. Our study provides strong evidence for the presence of a QTL on chromosome 19 with a major effect on LDL-C plasma levels in outbred Caucasian populations.

A study of the radiation chemistry of poly(chlorotrifluoroethylene) by ESR spectroscopy

The ESR spectra of poly(chlorotrifluoroethylene) were recorded following gamma-radiolysis under vacuum at room temperature and 77 K. The very broad spectrum at 77 K revealed little fine structure with which to identity the radicals formed upon irradiation, but subsequent photobleaching and annealing studies, together with radiolytic studies at higher temperatures, afforded scope for making radical assignments. Both main-chain radicals and a range of chain-end radicals have been identified. The G-values for radical formation were 1.55, 0.36 and 0.32 at 77 K, 273 K and room temperature, respectively. (C) 2003 Elsevier Science Ltd. All rights reserved.