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.

A comparison of low-storage strategies for spectral analysis in dissipative systems: filter diagonalisation in the Lanczos representation and harmonic inversion of the Chebychev-order-domain autocorrelation function

We compare the performance of two different low-storage filter diagonalisation (LSFD) strategies in the calculation of complex resonance energies of the HO2, radical. The first is carried out within a complex-symmetric Lanczos subspace representation [H. Zhang, S.C. Smith, Phys. Chem. Chem. Phys. 3 (2001) 2281]. The second involves harmonic inversion of a real autocorrelation function obtained via a damped Chebychev recursion [V.A. Mandelshtam, H.S. Taylor, J. Chem. Phys. 107 (1997) 6756]. We find that while the Chebychev approach has the advantage of utilizing real algebra in the time-consuming process of generating the vector recursion, the Lanczos, method (using complex vectors) requires fewer iterations, especially for low-energy part of the spectrum. The overall efficiency in calculating resonances for these two methods is comparable for this challenging system. (C) 2001 Elsevier Science B.V. All rights reserved.

Simulation of binary mixture adsorption of methane and CO2 at supercritical conditions in carbons

Knowledge of the adsorption behavior of coal-bed gases, mainly under supercritical high-pressure conditions, is important for optimum design of production processes to recover coal-bed methane and to sequester CO2 in coal-beds. Here, we compare the two most rigorous adsorption methods based on the statistical mechanics approach, which are Density Functional Theory (DFT) and Grand Canonical Monte Carlo (GCMC) simulation, for single and binary mixtures of methane and carbon dioxide in slit-shaped pores ranging from around 0.75 to 7.5 nm in width, for pressure up to 300 bar, and temperature range of 308-348 K, as a preliminary study for the CO2 sequestration problem. For single component adsorption, the isotherms generated by DFT, especially for CO2, do not match well with GCMC calculation, and simulation is subsequently pursued here to investigate the binary mixture adsorption. For binary adsorption, upon increase of pressure, the selectivity of carbon dioxide relative to methane in a binary mixture initially increases to a maximum value, and subsequently drops before attaining a constant value at pressures higher than 300 bar. While the selectivity increases with temperature in the initial pressure-sensitive region, the constant high-pressure value is also temperature independent. Optimum selectivity at any temperature is attained at a pressure of 90-100 bar at low bulk mole fraction of CO2, decreasing to approximately 35 bar at high bulk mole fractions. (c) 2005 American Institute of Chemical Engineers.

On the nonlinear Neumann problem with critical and supercritical nonlinearities

We investigate the solvability of the Neumann problem (1.1) involving a critical Sobolev exponent. In the first part of this work it is assumed that the coeffcients Q and h are at least continuous. Moreover Q is positive on overline Omega and lambda > 0 is a parameter. We examine the common effect of the mean curvature and the shape of the graphs of the coeffcients Q and h on the existence of low energy solutions. In the second part of this work we consider the same problem with Q replaced by - Q. In this case the problem can be supercritical and the existence results depend on integrability conditions on Q and h.

Protein fold recognition without Boltzmann statistics or explicit physical basis

We present a fast method for finding optimal parameters for a low-resolution (threading) force field intended to distinguish correct from incorrect folds for a given protein sequence. In contrast to other methods, the parameterization uses information from >10(7) misfolded structures as well as a set of native sequence-structure pairs. In addition to testing the resulting force field's performance on the protein sequence threading problem, results are shown that characterize the number of parameters necessary for effective structure recognition.

Adjustment to a low-control situation: Reexamining the role of coping responses

The aim of the study was to test a revised conceptualization of the role of coping in adjustment to a low-control stressor-women's adjustment to a failed in vitro fertilization (NF) attempt. Data were collected prior to the IVF treatment (Time I) and twice after the failed NF attempt (1 to 2 weeks after finding out the results, n = 171, and fi weeks later, n = 139). Initial adjustment was assessed at Time I, whereas measures of coping and both self-report and partner ratings of adjustment were obtained at Times ? and 3. As predicted, escapist strategies and problem-management strategies (mainly at Time ?) were associated with poor adjustment, whereas problem-appraisal coping was associated with better adjustment., There was also support for the proposed positive relationship between adjustment and emotional approach coping (on self-report adjustment).

Adsorption of binary hydrocarbon mixtures in carbon slit pores: A density functional theory study

Adsorption of binary hydrocarbon mixtures involving methane in carbon slit pores is theoretically studied here from the viewpoints of separation and of the effect of impurities on methane storage. It is seen that even small amounts of ethane, propane, or butane can significantly reduce the methane capacity of carbons. Optimal pore sizes and pressures, depending on impurity concentration, are noted in the present work, suggesting that careful adsorbent and process design can lead to enhanced separation. These results are consistent with earlier literature studies for the infinite dilution limit. For methane storage applications a carbon micropore width of 11.4 Angstrom (based on distance between centers of carbon atoms on opposing walls) is found to be the most suitable from the point of view of lower impurity uptake during high-pressure adsorption and greater impurity retention during low-pressure delivery. The results also theoretically confirm unusual recently reported observations of enhanced methane adsorption in the presence of a small amount of heavier hydrocarbon impurity.

Low symbiont diversity in southern Great Barrier Reef corals, relative to those of the Caribbean

The specific identity of endosymbiotic dinoflagellates (Symbiodinium spp.) from most zooxanthellate corals is unknown. In a survey of symbiotic cnidarians from the southern Great Barrier Reef (GBR), 23 symbiont types were identified from 86 host species representing 40 genera. A majority (>85%) of these symbionts belong to a single phylogenetic clade or subgenus (C) composed of closely related (as assessed by sequence data from the internal transcribed spacer region and the ribosomal large subunit gene), yet ecologically and physiologically distinct, types. A few prevalent symbiont types, or generalists, dominate the coral community of the southern GBR, whereas many rare and/or specific symbionts, or specialists, are found uniquely within certain host taxa. The comparison of symbiont diversity between southern GBR and Caribbean reefs shows an inverse relationship between coral diversity and symbiont diversity, perhaps as a consequence of more-rapid diversification of Caribbean symbionts. Among clade C types, generalists C1 and C3 are common to both Caribbean and southern GBR symbiont assemblages, whereas the rest are regionally endemic. Possibly because of environmental changes in the Caribbean after geographic isolation through the Quaternary period, a high proportion of Caribbean fauna associate with symbiont taxa from two other distantly related Symbiodinium clades (A and B) that rarely occur in Pacific hosts. The resilience of Porites spp. and the resistance of Montipora digitata to thermal stress and bleaching are partially explained by their association with a thermally tolerant symbiont type, whereas the indiscriminant widespread bleaching and death among certain Pacific corals, during El Nino Southern Oscillation events, are influenced by associations with symbionts possessing higher sensitivity to thermal stress.

Binary games with many players

We examine a problem with n players each facing the same binary choice. One choice is superior to the other. The simple assumption of competition - that an individual's payoff falls with a rise in the number of players making the same choice, guarantees the existence of a unique symmetric equilibrium (involving mixed strategies). As n increases, there are two opposing effects. First, events in the middle of the distribution - where a player finds itself having made the same choice as many others - become more likely, but the payoffs in these events fall. In opposition, events in the tails of the distribution - where a player finds itself having made the same choice as few others - become less likely, but the payoffs in these events remain high. We provide a sufficient condition (strong competition) under which an increase in the number of players leads to a reduction in the equilibrium probability that the superior choice is made.

Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

Low intraspecific variation in the rRNA internal transcribed spacer 2 (ITS2) of the Australian paralysis tick, Ixodes holocyclus

Ixodes holocyclus has a narrow, discontinuous distribution along the east coast of Australia. We studied ticks from 17 localities throughout the geographic range of this tick. The ITS2 of I. holocyclus is 793 bp long. We found nucleotide variation at eight of the 588 nucleotide positions (1.4%) that were compared for all ticks. There were eight different nucleotide sequences. Most sequences were not restricted to a particular geographic region. However, sequences F, G and H, which had an adenine at position 197, were found only in the far north of Queensland - all other ticks had a guanine at this position. The low level of intraspecific variation in this tick (0.7%) contrasts with the sequence divergence between L holocyclus and its close relative, I. cornuatus (13.1 %). These data indicate that L holocyclus does not contain cryptic species despite possible geographic isolation of some populations. We conclude that variation in the ITS2 is likely to be informative about the phylogeny of the group.

Application of a simplified method based on regular regime approach to determine the effective moisture diffusivity of mixture of low molecular weight sugars and maltodextrin during desorption

The binary diffusivities of water in low molecular weight sugars; fructose, sucrose and a high molecular weight carbohydrate; maltodextrin (DE 11) and the effective diffusivities of water in mixtures of these sugars (sucrose, glucose, fructose) and maltodextrin (DE 11) were determined using a simplified procedure based on the Regular Regime Approach. The effective diffusivity of these mixtures exhibited both the concentration and molecular weight dependence. Surface stickiness was observed in all samples during desorption, with fructose exhibiting the highest and maltodextrin the lowest. (C) 2002 Elsevier Science Ltd. All rights reserved.

Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

Short-term autocorrelation in Australian equities

This paper examines the statistical and economic significance of short-term autocorrelation in Australian equities. We document large negative first-order autocorrelation in individual stock returns. Preliminary results suggest this autocorrelation is economically significant, as two simple trading strategies based on the autocorrelation structure appear to yield large risk-adjusted returns. Further analysis, however, shows that these results are driven by the inclusion of small-capitalisation and low-priced stocks which are vulnerable to a number of market-microstructure-related problems. After revising the dataset to mitigate these problems, little evidence of economic significance remains.