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.

Swapping space for time and unfair tests of ecological models

Testing ecological models for management is an increasingly important part of the maturation of ecology as an applied science. Consequently, we need to work at applying fair tests of models with adequate data. We demonstrate that a recent test of a discrete time, stochastic model was biased towards falsifying the predictions. If the model was a perfect description of reality, the test falsified the predictions 84% of the time. We introduce an alternative testing procedure for stochastic models, and show that it falsifies the predictions only 5% of the time when the model is a perfect description of reality. The example is used as a point of departure to discuss some of the philosophical aspects of model testing.

PMCA1 mRNA expression in rat aortic myocytes: A real-time RT-PCR study

The plasmalemmal Ca2+ adenosine triphosphatase (PMCA) is a key regulator of Ca2+ efflux in vascular smooth muscle. In these studies are developed a realtime reverse transcriptase-polymerase chain reaction (real-time RT-PCR) assay for assessing PMCA1 mRNA levels in rat primary cultured aortic myocytes. This assay detected fetal bovine serum-induced increases in PMCA1 mRNA (relative to 18S rRNA) 4, 8, and 24 h after stimulation. Early fetal bovine serum-induced increases in PMCA1 mRNA were insensitive to the Ca2+ channel blockers nifedipine, flunarizine, and SKF-96365. These studies demonstrate the feasibility of real-time RT-PCR to assess mRNA levels of PMCA1 and illustrate dynamic regulation of this Ca2+ pump isoform in rat primary cultured aortic myocytes, (C) 2000 Academic Press.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

Occupation, hours worked, and leisure-time physical activity

Background. International research indicates that blue-collar employees typically exhibit lower rates of leisure-time physical activity. While lack of time and work demands are commonly reported barriers to activity, the extent to which time-at-work mediates the relationship between occupation and leisure-time physical activity is unclear. This study investigated the association between occupation, time spent in paid employment, and participation in leisure-time physical activity. Methods. This was a secondary analysis of cross-sectional data from the 1995 Australian Health Survey, focusing on employed persons ages 18-64 years (n = 24,454), Occupation was coded as per the Australian Standard Classification of Occupations and collapsed into three categories (professional, white-collar, blue-collar). Hours worked was categorized into eight levels, ranging from 1-14 to more than 50 h per week. Participation in leisure-time physical activity was categorized as either insufficient or sufficient for health, consistent with recommended levels of energy expenditure (1600 METS-min/fortnight). The relationship between occupation, hours worked, and leisure-time physical activity was examined using logistic regression. Analyses were conducted separately for male and female, and the results are presented as a series of models that successively adjust for a range of potential covariates: age, living arrangement, smoking status, body mass index, and self-reported health. Results. Individuals in blue-collar occupations were approximately 50% more likely to be classified as insufficiently active. This occupational variability in leisure-time physical activity was not explained by hours worked. There was a suggested relationship between hours worked and leisure-time physical activity; however, this differed between men and women, and was difficult to interpret. Conclusions. Occupational variability in leisure-time physical activity cannot be explained by hours worked. Therefore, reports that work constitutes a barrier to participation should be explored further. Identification of the factors contributing to occupational variability in leisure-time physical activity will add to our understanding of why population subgroups differ in their health risk profiles, and assist in the development of health promotion strategies to reduce rates of sedentariness and health inequalities. (C) 2000 American Health Foundation and Academic Press.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.

Residence time distributions of gas flowing through rotating drum bioreactors

Residence time distribution studies of gas through a rotating drum bioreactor for solid-state fermentation were performed using carbon monoxide as a tracer gas. The exit concentration as a function of time differed considerably from profiles expected for plug flow, plug flow with axial dispersion, and continuous stirred tank reactor (CSTR) models. The data were then fitted by least-squares analysis to mathematical models describing a central plug flow region surrounded by either one dead region (a three-parameter model) or two dead regions (a five-parameter model). Model parameters were the dispersion coefficient in the central plug flow region, the volumes of the dead regions, and the exchange rates between the different regions. The superficial velocity of the gas through the reactor has a large effect on parameter values. Increased superficial velocity tends to decrease dead region volumes, interregion transfer rates, and axial dispersion. The significant deviation from CSTR, plug flow, and plug flow with axial dispersion of the residence time distribution of gas within small-scale reactors can lead to underestimation of the calculation of mass and heat transfer coefficients and hence has implications for reactor design and scaleup. (C) 2001 John Wiley & Sons, Inc.

Estimating HIV incidence rates from age prevalence data in epidemic situations

We present a method of estimating HIV incidence rates in epidemic situations from data on age-specific prevalence and changes in the overall prevalence over time. The method is applied to women attending antenatal clinics in Hlabisa, a rural district of KwaZulu/Natal, South Africa, where transmission of HIV is overwhelmingly through heterosexual contact. A model which gives age-specific prevalence rates in the presence of a progressing epidemic is fitted to prevalence data for 1998 using maximum likelihood methods and used to derive the age-specific incidence. Error estimates are obtained using a Monte Carlo procedure. Although the method is quite general some simplifying assumptions are made concerning the form of the risk function and sensitivity analyses are performed to explore the importance of these assumptions. The analysis shows that in 1998 the annual incidence of infection per susceptible woman increased from 5.4 per cent (3.3-8.5 per cent; here and elsewhere ranges give 95 per cent confidence limits) at age 15 years to 24.5 per cent (20.6-29.1 per cent) at age 22 years and declined to 1.3 per cent (0.5-2.9 per cent) at age 50 years; standardized to a uniform age distribution, the overall incidence per susceptible woman aged 15 to 59 was 11.4 per cent (10.0-13.1 per cent); per women in the population it was 8.4 per cent (7.3-9.5 per cent). Standardized to the age distribution of the female population the average incidence per woman was 9.6 per cent (8.4-11.0 per cent); standardized to the age distribution of women attending antenatal clinics, it was 11.3 per cent (9.8-13.3 per cent). The estimated incidence depends on the values used for the epidemic growth rate and the AIDS related mortality. To ensure that, for this population, errors in these two parameters change the age specific estimates of the annual incidence by less than the standard deviation of the estimates of the age specific incidence, the AIDS related mortality should be known to within +/-50 per cent and the epidemic growth rate to within +/-25 per cent, both of which conditions are met. In the absence of cohort studies to measure the incidence of HIV infection directly, useful estimates of the age-specific incidence can be obtained from cross-sectional, age-specific prevalence data and repeat cross-sectional data on the overall prevalence of HIV infection. Several assumptions were made because of the lack of data but sensitivity analyses show that they are unlikely to affect the overall estimates significantly. These estimates are important in assessing the magnitude of the public health problem, for designing vaccine trials and for evaluating the impact of interventions. Copyright (C) 2001 John Wiley & Sons, Ltd.

Autonomic activity during human sleep as a function of time and sleep stage

While there is a developing understanding of the influence of sleep on cardiovascular autonomic activity in humans, there remain unresolved issues. In particular, the effect of time within the sleep period, independent of sleep stage, has not been investigated. Further, the influence of sleep on central sympathetic nervous system (SNS) activity is uncertain because results using the major method applicable to humans, the low frequency (LF) component of heart rate Variability (HRV), have been contradictory, and because the method itself is open to criticism. Sleep and cardiac activity were measured in 14 young healthy subjects on three nights. Data was analysed in 2-min epochs. All epochs meeting specified criteria were identified, beginning 2 h before, until 7 h after, sleep onset. Epoch values were allocated to 30-min bins and during sleep were also classified into stage 2, slow wave sleep (SWS) and rapid eye movement (REM) sleep. The measures of cardiac activity were heart irate (HR), blood pressure (BP), high frequency (HF) and LF components of HRV and pre-ejection period (PEP). During non-rapid eye movement (NREM) sleep autonomic balance shifted from sympathetic to parasympathetic dominance, although this appeared to be more because of a shift in parasympathetic nervous system (PNS) activity. Autonomic balance during REM was in general similar to wakefulness. For BP and the HF and LF components the change occurred abruptly at sleep onset and was then constant over time within each stage of sleep, indicating that any change in autonomic balance over the sleep period is a consequence of the changing distribution of sleep stages. Two variables, HR and PEP, did show time effects reflecting a circadian influence over HR and perhaps time asleep affecting PEP. While both the LF component and PEP showed changes consistent with reduced sympathetic tone during sleep, their pattern of change over time differed.

A time-harmonic inverse methodology for the design of RF coils in MRI

An inverse methodology is described to assist in the design of radio-frequency (RF) coils for magnetic resonance imaging (MRI) applications. The time-harmonic electromagnetic Green's functions are used to calculate current on the coil and shield cylinders that will generate a specified internal magnetic field. Stream function techniques and the method of moments are then used to implement this theoretical current density into an RF coil. A novel asymmetric coil operating for a 4.5 T MRI machine was designed and constructed using this methodology and the results are presented.

Structuring real-time Object-Z specifications

This paper presents a means of structuring specifications in real-time Object-Z: an integration of Object-Z with the timed refinement calculus. Incremental modification of classes using inheritance and composition of classes to form multi-component systems are examined. Two approaches to the latter are considered: using Object-Z's notion of object instantiation and introducing a parallel composition operator similar to those found in process algebras. The parallel composition operator approach is both more concise and allows more general modelling of concurrency. Its incorporation into the existing semantics of real-time Object-Z is presented.

A computational analysis of substrate binding strength by phosphorylase kinase and protein kinase A

Protein kinases exhibit various degrees of substrate specificity. The large number of different protein kinases in the eukaryotic proteomes makes it impractical to determine the specificity of each enzyme experimentally. To test if it were possible to discriminate potential substrates from non-substrates by simple computational techniques, we analysed the binding enthalpies of modelled enzyme-substrate complexes and attempted to correlate it with experimental enzyme kinetics measurements. The crystal structures of phosphorylase kinase and cAMP-dependent protein kinase were used to generate models of the enzyme with a series of known peptide substrates and non-substrates, and the approximate enthalpy of binding assessed following energy minimization. We show that the computed enthalpies do not correlate closely with kinetic measurements, but the method can distinguish good substrates from weak substrates and non-substrates. Copyright (C) 2002 John Wiley Sons, Ltd.

Computational immunology: The coming of age

The explosive growth in biotechnology combined with major advancesin information technology has the potential to radically transformimmunology in the postgenomics era. Not only do we now have readyaccess to vast quantities of existing data, but new data with relevanceto immunology are being accumulated at an exponential rate. Resourcesfor computational immunology include biological databases and methodsfor data extraction, comparison, analysis and interpretation. Publiclyaccessible biological databases of relevance to immunologists numberin the hundreds and are growing daily. The ability to efficientlyextract and analyse information from these databases is vital forefficient immunology research. Most importantly, a new generationof computational immunology tools enables modelling of peptide transportby the transporter associated with antigen processing (TAP), modellingof antibody binding sites, identification of allergenic motifs andmodelling of T-cell receptor serial triggering.