Parallel implementation of stochastic simulation for large scale cellular processes

Experimental and theoretical studies have shown the importance of stochastic processes in genetic regulatory networks and cellular processes. Cellular networks and genetic circuits often involve small numbers of key proteins such as transcriptional factors and signaling proteins. In recent years stochastic models have been used successfully for studying noise in biological pathways, and stochastic modelling of biological systems has become a very important research field in computational biology. One of the challenge problems in this field is the reduction of the huge computing time in stochastic simulations. Based on the system of the mitogen-activated protein kinase cascade that is activated by epidermal growth factor, this work give a parallel implementation by using OpenMP and parallelism across the simulation. Special attention is paid to the independence of the generated random numbers in parallel computing, that is a key criterion for the success of stochastic simulations. Numerical results indicate that parallel computers can be used as an efficient tool for simulating the dynamics of large-scale genetic regulatory networks and cellular processes

A first-principles density-functional calculation of the electronic and vibrational structure of the key melanin monomers

We report first-principles density-functional calculations for hydroquinone (HQ), indolequinone (IQ), and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of biomacromolecules with important biological functions (including photoprotection) and with the potential for certain bioengineering applications. We have used the difference of self-consistent fields method to study the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, HL. We show that HL is similar in IQ and SQ, but approximately twice as large in HQ. This may have important implications for our understanding of the observed broadband optical absorption of the eumelanins. The possibility of using this difference in HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules have significantly different infrared and Raman signatures, and so these spectra could be used in situ to nondestructively identify the monomeric content of macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure of eumelanin macromolecules and hence in helping to determine the structure-property-function relationships that control the behavior of the eumelanins.

Employee adjustment to large-scale organisational change

The present study was designed to test the utility of a stress-coping model of employee adjustment to organisational change. Specifically, it was proposed that employee adjustment to this type of work stress would be influenced by the characteristics of the change situation, employees' appraisals of the situation, their coping strategies, and the extent of their personal resources. Data were collected from 140 middle managers and supervisors involved in a large-scale public sector integration. The results of the research provided some support for the proposed model: high levels of psychological distress were related to a reliance on informal sources of information, high appraised stress, low appraised certainty, and the use of avoidant rather than problem-focused strategies, whereas poor social functioning was associated with low self-esteem, high levels or disruption across the period of change, a reliance on informal sources of information, and the use of avoidant coping strategies. There was no evidence that coping strategies mediated the effects of the event characteristics, situational appraisals, and personal resources on adjustment; however, there was some evidence linking these variables to coping strategies, in particular, problem-focused coping. There was also some evidence to indicate that the experience of organisational change was different for managers and supervisors: levels of threat were higher for the managers than the supervisors, but there was no difference between the groups of employees in terms of adjustment.

Autoinhibition by an internal nuclear localization signal revealed by the crystal structure of mammalian importin alpha

Importin alpha is the nuclear import receptor that recognizes classical monopartite and bipartite nuclear localization signals (NLSs). The structure of mouse importin alpha has been determined at 2.5 Angstrom resolution. The structure shows a large C-terminal domain containing armadillo repeats, and a less structured N-terminal importin beta-binding domain containing an internal NLS bound to the NLS-binding site. The structure explains the regulatory switch between the cytoplasmic, high-affinity form, and the nuclear, low-affinity form for NLS binding of the nuclear import receptor predicted by the current models of nuclear import. Importin beta conceivably converts the low- to high-affinity form by binding to a site overlapping the autoinhibitory sequence. The structure also has implications for understanding NLS recognition, and the structures of armadillo and HEAT repeats.

Effects of habitat modification in mangroves on the structure of mollusc and crab assemblages

The abundance and species richness of mollusc and crab assemblages were examined in a subtropical mangrove forest in Moreton Bay, Queensland, Australia, which has been disturbed and damaged by the construction of a wooden boardwalk and a path. Sections of the forest immediately adjacent to the boardwalk and path were compared with reference areas to determine whether changes to the small-scale structural complexity within the forest affected the benthic fauna. The disturbed area was characterised by having 65-80% fewer pneumatophores, significantly fewer species and individuals of molluscs, but significantly more species and individuals of crabs than the reference areas. The abundance of mangrove pneumatophores and the attached epiphytic algae were manipulated at two sites to determine whether observed differences in these features could account for the differences in the assemblage of molluscs in the disturbed area of the forest compared with reference areas. Five experimental treatments were used: undisturbed controls, pneumatophore removals (abundance reduced by ca. 65%), epiphytic algal removals (algae removed from ca. 65% of pneumatophores), pneumatophore disturbance controls and algal disturbance controls. The experimental reduction of the abundance of mangrove pneumatophores and the associated epiphytic algae led to significant declines (by as much as 83%) in the number of molluscs utilising the substratum in the modified plots. There was no significant difference in the abundance of molluscs in the pneumatophore and algal removal plots suggesting any effect was primarily related to removal of the epiphytic algae from the surface of the pneumatophores. The responses by the biota to the changes in the physical environment demonstrate that even relatively small-scale modifications to the physical structure of subtropical mangrove forests can lead to significant effects on the diversity and abundance of macrobenthic organisms in these habitats. Such modifications have the potential to cause cascading effects at higher trophic levels with a deterioration in the value of these habitats as nursery and feeding grounds. Future efforts at conservation of these estuarine environments must focus on the prevention or reduction of modifications to the physical structure and integrity of the system, rather than just on the prevention of loss of entire patches of habitat. (C) 2000 Elsevier Science B.V. All rights reserved.

Regional-scale investigation of the spatial distribution and origin of soluble salts in central north Queensland

A method for regional assessment of the distribution of saline outbreaks is demonstrated for a large area (68 000 km(2)) in north Queensland, Australia. Soil samples were used in conjunction with a digital elevation model and a map of potentially saline discharge zones to examine the landscape distribution of soluble salts in the region. The hypothesis of atmospheric accession of salt was tested for the topographically defined catchment regions feeding into each potentially saline discharge area. Most catchments showed a salt distribution consistent with this hypothesis, i.e. %TSS was large near the discharge areas and decreased rapidly with distance uphill from the discharge areas. In some catchments, however, local saline outbreaks were apparent at significant distances uphill from discharge areas. The possibility of geological sources of this salt was examined by comparing random point distributions with the location of saline points with distance downhill from geological units (excluding points near discharge zones). The distribution of some saline outbreaks was consistent with the occurrence of Cambro-Ordovician metasediments, Devonian limestone, Upper Devonian-Lower Carboniferous volcanics, and Triassic sediments. Copyright (C) 2000 John Wiley & Sons, Ltd.

The Fornax spectroscopic survey I. Survey strategy and preliminary results on the redshift distribution of a complete sample of stars and galaxies

The Fornax Spectroscopic Survey will use the Two degree Field spectrograph (2dF) of the Angle-Australian Telescope to obtain spectra for a complete sample of all 14000 objects with 16.5 less than or equal to b(j) less than or equal to 19.7 in a 12 square degree area centred on the Fornax Cluster. The aims of this project include the study of dwarf galaxies in the cluster (both known low surface brightness objects and putative normal surface brightness dwarfs) and a comparison sample of background field galaxies. We will also measure quasars and other active galaxies, any previously unrecognised compact galaxies and a large sample of Galactic stars. By selecting all objects-both stars and galaxies-independent of morphology, we cover a much larger range of surface brightness and scale size than previous surveys. In this paper we first describe the design of the survey. Our targets are selected from UK Schmidt Telescope sky survey plates digitised by the Automated Plate Measuring (APM) facility. We then describe the photometric and astrometric calibration of these data and show that the APM astrometry is accurate enough for use with the 2dF. We also describe a general approach to object identification using cross-correlations which allows us to identify and classify both stellar and galaxy spectra. We present results from the first 2dF field. Redshift distributions and velocity structures are shown for all observed objects in the direction of Fornax, including Galactic stars? galaxies in and around the Fornax Cluster, and for the background galaxy population. The velocity data for the stars show the contributions from the different Galactic components, plus a small tail to high velocities. We find no galaxies in the foreground to the cluster in our 2dF field. The Fornax Cluster is clearly defined kinematically. The mean velocity from the 26 cluster members having reliable redshifts is 1560 +/- 80 km s(-1). They show a velocity dispersion of 380 +/- 50 km s(-1). Large-scale structure can be traced behind the cluster to a redshift beyond z = 0.3. Background compact galaxies and low surface brightness galaxies are found to follow the general galaxy distribution.

Three-dimensional structure of RK-1: A novel alpha-defensin peptide

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of RK-1, an antimicrobial peptide from rabbit kidney recently discovered from homology screening based on the distinctive physicochemical properties of the corticostatins/defensins. RK-1 consists of 32 residues, including six cysteines arranged into three disulfide bonds. It exhibits antimicrobial activity against Escherichia coli and activates Ca2+ channels in vitro. Through its physicochemical similarity, identical cysteine spacing, and linkage to the corticostatins/defensins, it was presumed to be a member of this family. However, RK-1 lacks both a large number of arginines in the primary sequence and a high overall positive charge, which are characteristic of this family of peptides. The three-dimensional solution structure, determined by NMR, consists of a triple-stranded antiparallel beta -sheet and a series of turns and is similar to the known structures of other alpha -defensins. This has enabled the definitive classification of RK-1 as a member of this family of antimicrobial peptides. Ultracentrifuge measurements confirmed that like rabbit neutrophil defensins, RK-1 is monomeric in solution, in contrast to human neutrophil defensins, which are dimeric.

Identifying rate-limiting nodes in large-scale cortical networks for visuospatial processing: an illustration using fMRI

With the advent of functional neuroimaging techniques, in particular functional magnetic resonance imaging (fMRI), we have gained greater insight into the neural correlates of visuospatial function. However, it may not always be easy to identify the cerebral regions most specifically associated with performance on a given task. One approach is to examine the quantitative relationships between regional activation and behavioral performance measures. In the present study, we investigated the functional neuroanatomy of two different visuospatial processing tasks, judgement of line orientation and mental rotation. Twenty-four normal participants were scanned with fMRI using blocked periodic designs for experimental task presentation. Accuracy and reaction time (RT) to each trial of both activation and baseline conditions in each experiment was recorded. Both experiments activated dorsal and ventral visual cortical areas as well as dorsolateral prefrontal cortex. More regionally specific associations with task performance were identified by estimating the association between (sinusoidal) power of functional response and mean RT to the activation condition; a permutation test based on spatial statistics was used for inference. There was significant behavioral-physiological association in right ventral extrastriate cortex for the line orientation task and in bilateral (predominantly right) superior parietal lobule for the mental rotation task. Comparable associations were not found between power of response and RT to the baseline conditions of the tasks. These data suggest that one region in a neurocognitive network may be most strongly associated with behavioral performance and this may be regarded as the computationally least efficient or rate-limiting node of the network.

Community structure of Quaternary coral reefs compared with Recent life and death assemblages

This paper assesses the reliability with which fossil reefs record the diversity and community structure of adjacent Recent reefs. The diversity and taxonomic composition of Holocene raised fossil reefs was compared with those of modern reef coral life and death assemblages in adjacent moderate and low-energy shallow reef habitats Of Madang Lagoon, Papua New Guinea. Species richness per sample area and Shannon-Wiener diversity (H') were highest in the fossil reefs, intermediate in the life assemblages, and lowest in the death assemblages. The taxonomic composition of the fossil reefs was most similar to the combination of the life and death assemblages from the modern reefs adjacent to the two fossil reefs. Depth zonation was recorded accurately in the fossil reefs. The Madang fossil reefs represent time-averaged composites of the combined life and death assemblages as they existed at the time the reef was uplifted. Because fossil reefs include overlapping cohorts from the life and death assemblages, lagoonal facies of fossil reefs are dominated by the dominant sediment-producing taxa, which are not necessarily the most abundant in the life assemblage. Rare or slow-growing taxa accumulate more slowly than the encasing sediments and are underrepresented in fossil reef lagoons. Time-averaging dilutes the contribution of rare taxa, rather than concentrating their contribution. Consequently, fidelity indices developed for mollusks in sediments yield low values in coral reef death and fossil assemblages. Branching corals dominate lagoonal facies of fossil reefs because they are abundant, they grow and produce sediment rapidly, and most of the sediment they produce is not exported. Fossil reefs distinguished kilometer-scale variations in community structure more clearly than did the modern life assemblages. This difference implies that fossil,reefs may provide a better long-term record of community structure than modern reefs. This difference also suggests that modern kilometer-scale variation in coral reef community structure may have been reduced by anthropogenic degradation, even in the relatively unimpacted reefs of Madang Lagoon. Holocene and Pleistocene fossil reefs provide a time-integrated historical record of community composition and may be used as long-term benchmarks for comparison with modern, degraded, nearshore reefs. Comparisons between fossil reefs and degraded modern reefs display gross changes in community structure more effectively than they demonstrate local extinction of rare taxa.

A catalog of HI-selected galaxies from the south celestial cap region of sky

The first deep catalog of the H I Parkes All Sky Survey (HIPASS) is presented, covering the south celestial cap (SCC) region. The SCC area is similar to2400 deg(2) and covers delta < -62&DEG;. The average rms noise for the survey is 13 mJy beam(-1). Five hundred thirty-six galaxies have been cataloged according to their neutral hydrogen content, including 114 galaxies that have no previous cataloged optical counterpart. This is the largest sample of galaxies from a blind H I survey to date. Most galaxies in optically unobscured regions of sky have a visible optical counterpart; however, there is a small population of low-velocity H I clouds without visible optical counterparts whose origins and significance are unclear. The rms accuracy of the HIPASS positions is found to be 1.'9. The H I mass range of galaxies detected is from &SIM;10(6) to &SIM;10(11) M-.. There are a large number of late-type spiral galaxies in the SCC sample (66%), compared with 30% for optically selected galaxies from the same region in the NASA Extragalactic Database. The average ratio of H I mass to B luminosity of the sample increases according to optical type, from 1.8 M-./L-. for early types to 3.2 M-./L-. for late-type galaxies. The H I-detected galaxies tend to follow the large-scale structure traced by galaxies found in optical surveys. From the number of galaxies detected in this region of sky, we predict the full HIPASS catalog will contain &SIM;5000 galaxies, to a peak flux density limit of &SIM;39 mJy (3 σ), although this may be a conservative estimate as two large voids are present in the region. The H I mass function for this catalog is presented in a subsequent paper.

Mechanism of ribonuclease inhibition by ribonuclease inhibitor protein based on the crystal structure of its complex with ribonuclease A

We describe the mechanism of ribonuclease inhibition by ribonuclease inhibitor, a protein built of leucine-rich repeats, based on the crystal structure of the complex between the inhibitor and ribonuclease A. The structure was determined by molecular replacement and refined to an R(cryst) of 19.4% at 2.5 Angstrom resolution. Ribonuclease A binds to the concave region of the inhibitor protein comprising its parallel beta-sheet and loops. The inhibitor covers the ribonuclease active site and directly contacts several active-site residues. The inhibitor only partially mimics the RNase-nucleotide interaction and does not utilize the pi phosphate-binding pocket of ribonuclease A, where a sulfate ion remains bound. The 2550 Angstrom(2) of accessible surface area buried upon complex formation may be one of the major contributors to the extremely tight association (K-i = 5.9 x 10(-14) M). The interaction is predominantly electrostatic; there is a high chemical complementarity with 18 putative hydrogen bonds and salt links, but the shape complementarity is lower than in most other protein-protein complexes. Ribonuclease inhibitor changes its conformation upon complex formation; the conformational change is unusual in that it is a plastic reorganization of the entire structure without any obvious hinge and reflects the conformational flexibility of the structure of the inhibitor. There is a good agreement between the crystal structure and other biochemical studies of the interaction. The structure suggests that the conformational flexibility of RI and an unusually large contact area that compensates for a lower degree of complementarity may be the principal reasons for the ability of RI to potently inhibit diverse ribonucleases. However, the inhibition is lost with amphibian ribonucleases that have substituted most residues corresponding to inhibitor-binding residues in RNase A, and with bovine seminal ribonuclease that prevents inhibitor binding by forming a dimer. (C) 1996 Academic Press Limited

Quasi-stationarity in populations that are subject to large-scale mortality or emigration

We shall examine a model, first studied by Brockwell et al. [Adv Appl Probab 14 (1982) 709.], which can be used to describe the longterm behaviour of populations that are subject to catastrophic mortality or emigration events. Populations can suffer dramatic declines when disease, such as an introduced virus, affects the population, or when food shortages occur, due to overgrazing or fluctuations in rainfall. However, perhaps surprisingly, such populations can survive for long periods and, although they may eventually become extinct, they can exhibit an apparently stationary regime. It is useful to be able to model this behaviour. This is particularly true of the ecological examples that motivated the present study, since, in order to properly manage these populations, it is necessary to be able to predict persistence times and to estimate the conditional probability distribution of population size. We shall see that although our model predicts eventual extinction, the time till extinction can be long and the stationary exhibited by these populations over any reasonable time scale can be explained using a quasistationary distribution. (C) 2001 Elsevier Science Ltd. All rights reserved.