Effects of data characteristics on the results of reserve selection algorithms

We tested the effects of four data characteristics on the results of reserve selection algorithms. The data characteristics were nestedness of features (land types in this case), rarity of features, size variation of sites (potential reserves) and size of data sets (numbers of sites and features). We manipulated data sets to produce three levels, with replication, of each of these data characteristics while holding the other three characteristics constant. We then used an optimizing algorithm and three heuristic algorithms to select sites to solve several reservation problems. We measured efficiency as the number or total area of selected sites, indicating the relative cost of a reserve system. Higher nestedness increased the efficiency of all algorithms (reduced the total cost of new reserves). Higher rarity reduced the efficiency of all algorithms (increased the total cost of new reserves). More variation in site size increased the efficiency of all algorithms expressed in terms of total area of selected sites. We measured the suboptimality of heuristic algorithms as the percentage increase of their results over optimal (minimum possible) results. Suboptimality is a measure of the reliability of heuristics as indicative costing analyses. Higher rarity reduced the suboptimality of heuristics (increased their reliability) and there is some evidence that more size variation did the same for the total area of selected sites. We discuss the implications of these results for the use of reserve selection algorithms as indicative and real-world planning tools.

The crack tip strain field of AISI 4340 - Part I - Measurement technique

This is the first paper in a study on the influence of the environment on the crack tip strain field for AISI 4340. A stressing stage for the environmental scanning electron microscope (ESEM) was constructed which was capable of applying loads up to 60 kN to fracture-mechanics samples. The measurement of the crack tip strain field required preparation (by electron lithography or chemical etching) of a system of reference points spaced at similar to 5 mu m intervals on the sample surface, loading the sample inside an electron microscope, image processing procedures to measure the displacement at each reference point and calculation of the strain field. Two algorithms to calculate strain were evaluated. Possible sources of errors were calculation errors due to the algorithm, errors inherent in the image processing procedure and errors due to the limited precision of the displacement measurements. Estimation of the contribution of each source of error was performed. The technique allows measurement of the crack tip strain field over an area of 50 x 40 mu m with a strain precision better than +/- 0.02 at distances larger than 5 mu m from the crack tip. (C) 1999 Kluwer Academic Publishers.

On modifications to the long-term survival mixture model in the presence of competing risks

A mixture model for long-term survivors has been adopted in various fields such as biostatistics and criminology where some individuals may never experience the type of failure under study. It is directly applicable in situations where the only information available from follow-up on individuals who will never experience this type of failure is in the form of censored observations. In this paper, we consider a modification to the model so that it still applies in the case where during the follow-up period it becomes known that an individual will never experience failure from the cause of interest. Unless a model allows for this additional information, a consistent survival analysis will not be obtained. A partial maximum likelihood (ML) approach is proposed that preserves the simplicity of the long-term survival mixture model and provides consistent estimators of the quantities of interest. Some simulation experiments are performed to assess the efficiency of the partial ML approach relative to the full ML approach for survival in the presence of competing risks.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Chemical synthesis, antibacterial activity and conformation of diptericin, an 82-mer peptide originally isolated from insects

The small amounts of antibacterial peptides that can be isolated from insects do not allow detailed studies of their range of activity, side-chain sugar requirements, or their conformation, factors that frequently play roles in the mode of action. In this paper, we report the solid-phase step-by-step synthesis of diptericin, an 82-mer peptide, originally isolated from Phormia terranovae. The unglycosylated peptide was purified to homogeneity by conventional reversed-phase high performance liquid chromatography, and its activity spectrum was compared to that Of synthetic unglycosylated drosocin, which shares strong sequence homology with diptericin's N-terminal domain. Diptericin appeared to have antibacterial activity:for only a limited number of Gram-negative bacteria. Diptericin's submicromolar potency against Escherichia coli strains indicated that, in a manner similar to drosocin, the presence of the carbohydrate side chain is not,necessary to kill bacteria. Neither the N-terminal, drosocin-analog fragment, nor the C-terminal, glycine-rich attacin-analog region was active against any of the bacterial strains studied, regardless of whether the Gal-GalNAc disaccharide units were attached. This suggested that the active site of diptericin fell outside the drosocin or attacin homology domains. In addition, the conformation of diptericin did not seem to play a role in the antibacterial activity, as was demonstrated by the complete lack of ordered structure by two-dimensional nuclear magnetic resonance spectroscopy and circular dichroism. Diptericin completely killed bacteria within I h, considerably faster than drosocin and the attacins; unlike some other, fast-acting antibacterial peptides, diptericin did not lyse normal mammalian cells. Taken together, these data suggest diptericin does not belong to any known class of antibacterial peptides.

Knowledge discovery and data mining in biological databases

The new technologies for Knowledge Discovery from Databases (KDD) and data mining promise to bring new insights into a voluminous growing amount of biological data. KDD technology is complementary to laboratory experimentation and helps speed up biological research. This article contains an introduction to KDD, a review of data mining tools, and their biological applications. We discuss the domain concepts related to biological data and databases, as well as current KDD and data mining developments in biology.

Activation of the human diaphragm during a repetitive postural task

1. The co-ordination between respiratory and postural functions of the diaphragm was investigated during repetitive upper Limb movement. It was hypothesised that diaphragm activity would occur either tonically or phasically in association with the forces from each movement and that this activity would combine with phasic respiratory activity. 2. Movements of the upper limb and ribcage were measured while standing subjects performed repetitive upper limb movements 'as fast as possible'. Electromyographic (EMG) recordings of the costal diaphragm were made using intramuscular electrodes in four subjects. Surface electrodes were placed over the deltoid and erector spinae muscles. 3. In contrast to standing at rest, diaphragm activity was present throughout expiration at 78 +/- 17% (mean +/- S.D.) of its peak inspiratory magnitude during repeated upper limb movement. 4. Bursts of deltoid and erector spinae EMG activity occurred at the Limb movement frequency (similar to 2.9 Hz). Although the majority of diaphragm EMG power was at the respiratory frequency (similar to 0.4 Hz), a peak was also present at the movement frequency. This finding was corroborated by averaged EMG activity triggered from upper limb movement. In addition, diaphragm EMG activity was coherent with ribcage motion at the respiratory frequency and with upper limb movement at the movement frequency. 5. The diaphragm response was similar when movement was performed while sitting. In addition, when subjects moved with increasing frequency the peak upper limb acceleration correlated with diaphragm EMG amplitude. These findings support the argument that diaphragm contraction is related to trunk control. 6. The results indicate that activity of human phrenic motoneurones is organised such that it contributes to both posture and respiration during a task which repetitively challenges trunk posture.

Modeling of hepatic elimination and organ distribution kinetics with the extended convection-dispersion model

The conventional convection-dispersion (also called axial dispersion) model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. An extended form of the convection-dispersion model has been developed to adequately describe the outflow concentration-time profiles for vascular markers at both short and long times after bolus injections into perfused livers. The model, based on flux concentration and a convolution of catheters and large vessels, assumes that solute elimination in hepatocytes follows either fast distribution into or radial diffusion in hepatocytes. The model includes a secondary vascular compartment, postulated to be interconnecting sinusoids. Analysis of the mean hepatic transit time (MTT) and normalized variance (CV2) of solutes with extraction showed that the discrepancy between the predictions of MTT and CV2 for the extended and conventional models are essentially identical irrespective of the magnitude of rate constants representing permeability, volume, and clearance parameters, providing that there is significant hepatic extraction. In conclusion, the application of a newly developed extended convection-dispersion model has shown that the unweighted conventional convection-dispersion model can be used to describe the disposition of extracted solutes and, in particular, to estimate hepatic availability and clearance in booth experimental and clinical situations.

Design of the multicenter Australian study of epidural anesthesia and analgesia in major surgery: The MASTER Trial

The Multicenter Australian Study of Epidural Anesthesia and Analgesia in Major Surgery (The MASTER Trial) was designed to evaluate the possible benefit of epidural block in improving outcome in high-risk patients. The trial began in 1995 and is scheduled to reach the planned sample size of 900 during 2001. This paper describes the trial design and presents data comparing 455 patients randomized in 21 institutions in Australia, Hong Kong, and Malaysia, with 237 patients from the same hospitals who were eligible but not randomized. Nine categories of high-risk patients were defined as entry criteria for the trial. Protocols for ethical review, informed consent, randomization, clinical anesthesia and analgesia, and perioperative management were determined following extensive consultation with anesthesiologists throughout Australia. Clinical and research information was collected in participating hospitals by research staff who may not have been blind to allocation. Decisions about the presence or absence of endpoints were made primarily by a computer algorithm, supplemented by blinded clinical experts. Without unblinding the trial, comparison of eligibility criteria and incidence of endpoints between randomized and nonrandomized patients showed only small differences. We conclude that there is no strong evidence of important demographic or clinical differences between randomized and nonrandomized patients eligible for the MASTER Trial. Thus, the trial results are likely to be broadly generalizable. Control Clin Trials 2000;21:244-256 (C) Elsevier Science Inc. 2000.

Muscle fibre types and size distribution in sub-antarctic notothenioid fishes

The presumptive tonic muscles fibres of Cottoperca gobio, Champsocephalus esox, Harpagifer bispinis, Eleginops maclovinus, Patagonothen tessellata, P. cornucola and Paranotothenia magellanica stained weakly or were unstained for glycogen, lipid, succinic dehydrogenase (SDHase) and myosin ATPase (mATPase) activity. Slow, intermediate and fast twitch muscle fibres, distinguished on the basis of the pH stability of their mATPases, showed intense, moderate and low staining activity for SDHase, respectively. Slow fibres were the major component of the pectoral fin adductor profundis muscle. The proportion of different muscle fibre types varied from the proximal to distal end of the muscle, but showed relatively little variation between species. The myotomes contained a lateral superficial strip of red muscle composed of presumptive tonic, slow twitch and intermediate fibres, thickening to a major wedge at the horizontal septum. All species also had characteristic secondary dorsal and ventral wedges of red muscle. The relative abundance and localization of muscle fibre types in the red muscle varied between species and with body size in the protandric hermaphrodite E. maclovinus. The frequency distribution of diameters for fast twitch muscle fibres, the major component of deep white muscle, was determined in fish of a range of body sizes. The absence of fibres <20 mu m diameter was used as a criterion for the cessation of muscle fibre recruitment. Fibre recruitment had stopped in P, tessellata of 13.8 cm L-T and E, maclovinus of 32.8 cm L-T, equivalent to 49 and 36.5% of their recorded maximum sizes respectively. As a result in 20-cm P. tessellata, the maximum fibre diameter was 300 mu m and 36% of fibres were in excess of 200 mu m The unusually large maximum fibre diameter, the general arrangement of the red muscle layer and the extreme pH lability of the mATPase of fast twitch fibres are all common characters of the sub-Antarctic and Antarctic Notothenioids, including Cottoperca gobio, the suggested sister group to the Notothenidae. (C) 2000 The Fisheries Society of the British Isles.

Computational chemistry on Fujitsu vector-parallel processors: Development and performance of applications software

In this and a preceding paper, we provide an introduction to the Fujitsu VPP range of vector-parallel supercomputers and to some of the computational chemistry software available for the VPP. Here, we consider the implementation and performance of seven popular chemistry application packages. The codes discussed range from classical molecular dynamics to semiempirical and ab initio quantum chemistry. All have evolved from sequential codes, and have typically been parallelised using a replicated data approach. As such they are well suited to the large-memory/fast-processor architecture of the VPP. For one code, CASTEP, a distributed-memory data-driven parallelisation scheme is presented. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

Effects of Mg2+ on Ca2+ release from sarcoplasmic reticulum of skeletal muscle fibres from yabby (crustacean) and rat

1. The role of myoplasmic [Mg2+] on Ca2+ release from the sarcoplasmic reticulum (SR) was examined in the two major types of crustacean muscle fibres, the tonic, long sarcomere fibres and the phasic, short sarcomere fibres of the fresh mater decapod crustacean Cherax: destructor (yabby) and in the fast-twitch rat muscle fibres using the mechanically skinned muscle fibre preparation. 2. A robust Ca2+-induced Ca2+-release (CICR) mechanism was present in both long and short sarcomere fibres and 1 mM Mg2+ exerted a strong inhibitory action on the XR Ca2+ release in both fibre types. 3. The XR displayed different properties with respect to Ca2+ loading in the long and the short sarcomere fibres and marked functional differences were identified with respect to Mg2+ inhibition between the two crustacean fibre types. Thus, in long sarcomere fibres, the submaximally loaded XR was able to release Ca2+ when [Mg2+] was lowered from 1 to 0.01 mw in the presence of 8 mM ATP(total) and in the virtual absence of Ca2+ (< 5 nM) even when the CICR was suppressed. In contrast, negligible Ca2+ was released from the submaximally loaded SR of short sarcomere yabby fibres when [Mg2+] was lowered from 1. to 0.01 mM under the same conditions as for the long sarcomere fibres. Nevertheless, the rate of XR Ca2+ release in short sarcomere fibres increased markedly when [Mg2+] was lowered in the presence of [Ca2+] approaching the normal resting levels (50-100 nM). 4. Rat fibres were able to release SR Ca2+ at a faster rate than the long sarcomere yabby fibres when [Mg2+] was lowered from 1 to 0.01 mM in the virtual absence of Ca2+ but, unlike with yabby fibres, the net rate of Ca2+ release was actually increased for conditions that were considerably less favourable to CICR. 5. In summary it is concluded that crustacean skeletal muscles have more that one functional type of Ca2+-release channels, that these channels display properties that are intermediate between those of mammalian skeletal and cardiac isoforms, that the inhibition exerted by Mg2+ at rest on the crustacean SR Ca2+-release channels must be removed during excitation-contraction coupling and that, unlike in crustacean fibres, CICR cannot play the major role in the activation of XR Ca2+-release channels in the rat skeletal muscle.

X-ray absorption and phase contrast imaging to study the interplay between plant roots and soil structure

Plant performance is, at least partly, linked to the location of roots with respect to soil structure features and the micro-environment surrounding roots. Measurements of root distributions from intact samples, using optical microscopy and field tracings have been partially successful but are imprecise and labour-intensive. Theoretically, X-ray computed micro-tomography represents an ideal solution for non-invasive imaging of plant roots and soil structure. However, before it becomes fast enough and affordable or easily accessible, there is still a need for a diagnostic tool to investigate root/soil interplay. Here, a method for detection of undisturbed plant roots and their immediate physical environment is presented. X-ray absorption and phase contrast imaging are combined to produce projection images of soil sections from which root distributions and soil structure can be analyzed. The clarity of roots on the X-ray film is sufficient to allow manual tracing on an acetate sheet fixed over the film. In its current version, the method suffers limitations mainly related to (i) the degree of subjectivity associated with manual tracing and (ii) the difficulty of separating live and dead roots. The method represents a simple and relatively inexpensive way to detect and quantify roots from intact samples and has scope for further improvements. In this paper, the main steps of the method, sampling, image acquisition and image processing are documented. The potential use of the method in an agronomic perspective is illustrated using surface and sub-surface soil samples from a controlled wheat trial. Quantitative characterization of root attributes, e.g. radius, length density, branching intensity and the complex interplay between roots and soil structure, is presented and discussed.

Robust Mixture Modelling Using the t Distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.

Nitrogen dynamics and the physiological basis of stay-green in sorghum

Sorghum [Sorghum bicolor (L,) Moench] hybrids containing the stay-green trait retain more photosynthetically active leaves under drought than do hybrids that do not contain this trait. Since the Longevity and photosynthetic capacity of a leaf are related to its N status, it is important to clarify the role of N in extending leaf greenness in stay-green hybrids. Field studies were conducted in northeastern Australia to examine the effect of three water regimes and nine hybrids on N uptake and partitioning among organs. Nine hybrids varying in the B35 and KS19 sources of stay-green were grown under a fully irrigated control, post-flowering water deficit, and terminal water deficit. For hybrids grown under terminal water deficit, stay-green was viewed as a consequence of the balance between N demand by the grain and N supply during gain filling. On the demand side, grain numbers were 16% higher in the four stay-green than in the five senescent hybrids. On the supply side, age-related senescence provided an average of 34 and 42 kg N ha(-1) for stay-green and senescent hybrids, respectively. In addition, N uptake during grain filling averaged 116 and 82 kg ha(-1) in stay-green and senescent hybrids. Matching the N supply from these two sources with grain N demand found that the shortfall in N supply for grain filling in the stay-green and senescent hybrids averaged 32 and 41 kg N ha(-1) resulting in more accelerated leaf senescence in the senescent hybrids. Genotypic differences in delayed onset and reduced rate of leaf senescence were explained by differences in specific leaf nitrogen and N uptake during grain filling. Leaf nitrogen concentration at anthesis was correlated with onset (r = 0.751**, n = 27) and rate (r = -0.783**, n = 27) of leaf senescence ender terminal water deficit.