A systematic approach to error isolation in computerized wastewater simulation models

Activated sludge models are used extensively in the study of wastewater treatment processes. While various commercial implementations of these models are available, there are many people who need to code models themselves using the simulation packages available to them, Quality assurance of such models is difficult. While benchmarking problems have been developed and are available, the comparison of simulation data with that of commercial models leads only to the detection, not the isolation of errors. To identify the errors in the code is time-consuming. In this paper, we address the problem by developing a systematic and largely automated approach to the isolation of coding errors. There are three steps: firstly, possible errors are classified according to their place in the model structure and a feature matrix is established for each class of errors. Secondly, an observer is designed to generate residuals, such that each class of errors imposes a subspace, spanned by its feature matrix, on the residuals. Finally. localising the residuals in a subspace isolates coding errors. The algorithm proved capable of rapidly and reliably isolating a variety of single and simultaneous errors in a case study using the ASM 1 activated sludge model. In this paper a newly coded model was verified against a known implementation. The method is also applicable to simultaneous verification of any two independent implementations, hence is useful in commercial model development.

Linear models of simple cells: Correspondence to real cell responses and space spanning properties

Despite their limitations, linear filter models continue to be used to simulate the receptive field properties of cortical simple cells. For theoreticians interested in large scale models of visual cortex, a family of self-similar filters represents a convenient way in which to characterise simple cells in one basic model. This paper reviews research on the suitability of such models, and goes on to advance biologically motivated reasons for adopting a particular group of models in preference to all others. In particular, the paper describes why the Gabor model, so often used in network simulations, should be dropped in favour of a Cauchy model, both on the grounds of frequency response and mutual filter orthogonality.

Stable expression of noncytopathic Kunjin replicons simulates both ultrastructural and biochemical characteristics observed during replication of Kunjin virus

This report focuses mainly on the characterization of a Vero cell line stably expressing the flavivirus Kunjin (KUN) replicon C20SDrep (C20SDrepVero). We showed by immunofluorescence and cryoimmunoelectron microscopy that unique flavivirus-induced membrane structures, termed convoluted membranes/paracrystalline structures, were induced in the C20SDrepVero cells. These induced cytoplasmic foci were immunolabeled with KUN virus anti-NS3 antibodies and with antibodies to the cellular markers ERGIC53 (for the intermediate compartment) and protein disulfide isomerase (for the rough endoplasmic reticulum). However, in contrast to the large perinuclear inclusions observed by immunofluorescence with anti-double-stranded (ds)RNA antibodies in KUN virus-infected cells, the dsRNA in C20SDrepVero cells was localized to small isolated foci scattered throughout the cytoplasm, which were coincident with small foci dual-labeled with the trans-Golgi specific marker GaIT. importantly persistent expression of the KUN replicons in cells did not produce cytopathic effects, and the morphology of major host organelles (including Golgi, mitochondria, endoplasmic reticulum, and nucleus) was apparently unaffected. The amounts of plus- and minus-sense RNA synthesis in replicon cells were similar to those in KUN virus-infected cells until near the end of the latent period, but subsequently increases of about 10- and fourfold, respectively, occurred in infected cells. Virus-specified protein synthesis in C20SDrepVero cells was also about 10-fold greater than that in infected cells. When several KUN replicon cell lines were compared with respect to membrane induction, the relative efficiencies increased in parallel with increases in viral RNA and protein synthesis, consistent with the increases observed during the virus infectious cycle. Based on these observations, cell lines expressing less-efficient replicons may provide a useful tool to study early events in flavivirus RNA replication, which are difficult to assess in Virus infections. (C) 2001 Academic press.

The investigation of optimal bombardment parameters for transient and stable transgene expression in sorghum

This report outlines the development of optimized particle inflow gun (PIG) parameters for producing transgenic sorghum (Sorghum bicolor (L.) Moench). Both transient and stable expression were examined when determining these parameters. The uidA reporter gene (GUS) encoding beta -glucuronidase was used in transient experiments and the green fluorescent protein (GFP) used to monitor stable expression. Initially, optimization was conducted using leaf segments, as the generation of sorghum callus in sufficiently large quantities is time-consuming. Following leaf optimization, experiments were conducted using callus, identifying a high similarity between the two tissue types (r(s) = 0.83). High levels of GUS expression were observed in both leaf and callus material when most distant from the DNA expulsion point, and using a pressure greater than 1800 kPa. A higher level of expression was also observed when the aperture of the helium inlet valve was constricted. Using the optimized conditions (pressure of 2200 kPa, distance to target tissue of 15 cm from the expulsion point, and the aperture of the helium inlet valve at one full turn), three promoters (Ubiquitin, Actin1 and CaMV 35S) were evaluated over a 72-h period using GUS as the reporter gene. A significantly higher number of GUS foci were counted with the Ubiquitin construct over this period, compared to the Actin1 and CaMV 35S constructs. Stable callus sectors (on 2 mg l(-1) bialaphos) with GFP expression were visualized for as long as 6 wk post-bombardment. Using this optimized protocol, several plants were regenerated after having been bombarded with the pAHC20 construct (containing the bar gene), with molecular evidence confirming integration.

Discrimination of adsorption kinetic models for the description of hydrocarbon adsorption in activated carbon

Five kinetic models for adsorption of hydrocarbons on activated carbon are compared and investigated in this study. These models assume different mass transfer mechanisms within the porous carbon particle. They are: (a) dual pore and surface diffusion (MSD), (b) macropore, surface, and micropore diffusion (MSMD), (c) macropore, surface and finite mass exchange (FK), (d) finite mass exchange (LK), and (e) macropore, micropore diffusion (BM) models. These models are discriminated using the single component kinetic data of ethane and propane as well as the multicomponent kinetics data of their binary mixtures measured on two commercial activated carbon samples (Ajax and Norit) under various conditions. The adsorption energetic heterogeneity is considered for all models to account for the system. It is found that, in general, the models assuming diffusion flux of adsorbed phase along the particle scale give better description of the kinetic data.

Natural abundance of stable carbon and nitrogen isotopes in Cannabis sativa reflects growth conditions

Stable carbon and nitrogen isotope signatures (delta C-13 and delta N-15) of Cannabis sativa were assessed for their usefulness to trace seized Cannabis leaves to the country of origin and to source crops by determining how isotope signatures relate to plant growth conditions. The isotopic composition of Cannabis examined here covered nearly the entire range of values reported for terrestrial C-3 plants. The delta C-13 values of Cannabis from Australia, Papua New Guinea and Thailand ranged from -36 to -25 parts per thousand, and delta N-15 values ranged from -1.0 to 15.8 parts per thousand. The stable isotope content did not allow differentiation between Cannabis originating from the three countries, but delta C-13 values of plantation-grown Cannabis differed between well-watered plants (average delta C-13 of -30.0 parts per thousand) and plants that had received little irrigation (average delta C-13 of -26.4 parts per thousand). Cannabis grown under controlled conditions had delta C-13 values of -32.6 and -30.6 parts per thousand with high and low water supply, respectively. These results indicate that water availability determines leaf C-13 in plants grown under similar conditions of light, temperature and air humidity. The delta C-13 values also distinguished between indoor- and outdoor-grown Cannabis; indoor- grown plants had overall more negative delta C-13 values (average -31.8 parts per thousand) than outdoor-grown plants (average -27.9 parts per thousand). Contributing to the strong C-13-depletion of indoor- grown plants may be high relative humidity, poor ventilation and recycling of C-13-depleted respired CO2. Mineral fertilizers had mostly lower delta N-15 values (-0.2 to 2.2 parts per thousand) than manure-based fertilizers (7.6 to 22.7 parts per thousand). It was possible to link delta N-15 values of fertilizers associated with a crop site to soil and plant delta N-15 values. The strong relationship between soil, fertilizer, and plant delta N-15 suggests that Cannabis delta N-15 is determined by the isotopic composition of the nitrogen source. The distinct delta N-15 values measured in Cannabis crops make delta N-15 an excellent tool for matching seized Cannabis with a source crop. A case study is presented that demonstrates how delta C-13 and delta N-15 values can be used as a forensic tool.

Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits

Understanding the genetic architecture of quantitative traits can greatly assist the design of strategies for their manipulation in plant-breeding programs. For a number of traits, genetic variation can be the result of segregation of a few major genes and many polygenes (minor genes). The joint segregation analysis (JSA) is a maximum-likelihood approach for fitting segregation models through the simultaneous use of phenotypic information from multiple generations. Our objective in this paper was to use computer simulation to quantify the power of the JSA method for testing the mixed-inheritance model for quantitative traits when it was applied to the six basic generations: both parents (P-1 and P-2), F-1, F-2, and both backcross generations (B-1 and B-2) derived from crossing the F-1 to each parent. A total of 1968 genetic model-experiment scenarios were considered in the simulation study to quantify the power of the method. Factors that interacted to influence the power of the JSA method to correctly detect genetic models were: (1) whether there were one or two major genes in combination with polygenes, (2) the heritability of the major genes and polygenes, (3) the level of dispersion of the major genes and polygenes between the two parents, and (4) the number of individuals examined in each generation (population size). The greatest levels of power were observed for the genetic models defined with simple inheritance; e.g., the power was greater than 90% for the one major gene model, regardless of the population size and major-gene heritability. Lower levels of power were observed for the genetic models with complex inheritance (major genes and polygenes), low heritability, small population sizes and a large dispersion of favourable genes among the two parents; e.g., the power was less than 5% for the two major-gene model with a heritability value of 0.3 and population sizes of 100 individuals. The JSA methodology was then applied to a previously studied sorghum data-set to investigate the genetic control of the putative drought resistance-trait osmotic adjustment in three crosses. The previous study concluded that there were two major genes segregating for osmotic adjustment in the three crosses. Application of the JSA method resulted in a change in the proposed genetic model. The presence of the two major genes was confirmed with the addition of an unspecified number of polygenes.

Immune responses to ISCOM (R) formulations in animal and primate models

ISCOMs(R) are typically 40 nm cage-like structures comprising antigen, saponin, cholesterol and phospholipid. ISCOMs(R) have been shown to induce antibody responses and activate T helper cells and cyrolytic T lymphocytes in a number of animal species, including non-human primates. Recent clinical studies have demonstrated that ISCOMs(R) are also able to induce antibody and cellular immune responses in humans. This review describes the current understanding of the ability of ISCOMs(R) to induce immune responses and the mechanisms underlying this property. Recent progress in the characterisation and manufacture of ISCOMs(R) will also be discussed. (C) 2001 Elsevier Science Ltd. All rights reserved.

Hemoglobin-degrading, aspartic proteases of blood-feeding parasites - Substrate specificity revealed by homology models

Blood-feeding parasites, including schistosomes, hookworms, and malaria parasites, employ aspartic proteases to make initial or early cleavages in ingested host hemoglobin. To better understand the substrate affinity of these aspartic proteases, sequences were aligned with and/or three-dimensional, molecular models were constructed of the cathepsin D-like aspartic proteases of schistosomes and hookworms and of plasmepsins of Plasmodium falciparum and Plasmodium vivax, using the structure of human cathepsin D bound to the inhibitor pepstatin as the template. The catalytic subsites S5 through S4' were determined for the modeled parasite proteases. Subsequently, the crystal structure of mouse renin complexed with the nonapeptidyl inhibitor t-butyl-CO-His-Pro-Phe-His-Leu [CHOHCH2]Leu-Tyr-Tyr-Ser-NH2 (CH-66) was used to build homology models of the hemoglobin-degrading peptidases docked with a series of octapeptide substrates. The modeled octapeptides included representative sites in hemoglobin known to be cleaved by both Schistosoma japonicum cathepsin D and human cathepsin D, as well as sites cleaved by one but not the other of these enzymes. The peptidase-octapeptide substrate models revealed that differences in cleavage sites were generally attributable to the influence of a single amino acid change among the P5 to P4' residues that would either enhance or diminish the enzymatic affinity. The difference in cleavage sites appeared to be more profound than might be expected from sequence differences in the enzymes and hemoglobins. The findings support the notion that selective inhibitors of the hemoglobin-degrading peptidases of blood-feeding parasites at large could be developed as novel anti-parasitic agents.

The plausibility of long-wavelength stress correlation or stress magnitude as a mechanism for precursory seismicity: Results from two simple elastic models

Observations of accelerating seismic activity prior to large earthquakes in natural fault systems have raised hopes for intermediate-term eartquake forecasting. If this phenomena does exist, then what causes it to occur? Recent theoretical work suggests that the accelerating seismic release sequence is a symptom of increasing long-wavelength stress correlation in the fault region. A more traditional explanation, based on Reid's elastic rebound theory, argues that an accelerating sequence of seismic energy release could be a consequence of increasing stress in a fault system whose stress moment release is dominated by large events. Both of these theories are examined using two discrete models of seismicity: a Burridge-Knopoff block-slider model and an elastic continuum based model. Both models display an accelerating release of seismic energy prior to large simulated earthquakes. In both models there is a correlation between the rate of seismic energy release with the total root-mean-squared stress and the level of long-wavelength stress correlation. Furthermore, both models exhibit a systematic increase in the number of large events at high stress and high long-wavelength stress correlation levels. These results suggest that either explanation is plausible for the accelerating moment release in the models examined. A statistical model based on the Burridge-Knopoff block-slider is constructed which indicates that stress alone is sufficient to produce accelerating release of seismic energy with time prior to a large earthquake.