Computing van der Waals energies in the context of the rotamer approximation

The rotamer approximation states that protein side-chain conformations can be described well using a finite set of rotational isomers. This approximation is often applied in the context of computational protein design and structure prediction to reduce the complexity of structural sampling. It is an effective way of reducing the structure space to the most relevant conformations. However, the appropriateness of rotamers for sampling structure space does not imply that a rotamer-based energy landscape preserves any of the properties of the true continuous energy landscape. Specifically, because the energy of a van der Waals interaction can be very sensitive to small changes in atomic separation, meaningful van der Waals energies are particularly difficult to calculate from rotamer-based structures. This presents a problem for computational protein design, where the total energy of a given structure is often represented as a sum of precalculated rigid rotamer self and pair contributions. A common way of addressing this issue is to modify the van der Waals function to reduce its sensitivity to atomic position, but excessive modification may result in a strongly nonphysical potential. Although many different van der Waals modifications have been used in protein design, little is known about which performs best, and why. In this paper, we study 10 ways of computing van der Waals energies under the rotamer approximation, representing four general classes, and compare their performance using a variety of metrics relevant to protein design and native-sequence repacking calculations. Scaling van der Waals radii by anywhere from 85 to 95% gives the best performance. Linearizing and capping the repulsive portion of the potential can give additional improvement, which comes primarily from getting rid of unrealistically large clash energies. On the other hand, continuously minimizing individual rotamer pairs prior to evaluating their interaction works acceptably in native-sequence repacking, but fails in protein design. Additionally, we show that the problem of predicting relevant van der Waals energies from rotamer-based structures is strongly nonpairwise decomposable and hence further modifications of the potential are unlikely to give significant improvement.     

Systematic comparison of catalytic mechanisms of hydrolysis and transfer reactions classified in the EzCatDB database

Catalytic mechanisms of 270 enzymes from 131 superfamilies, mainly hydrolases and transferases, were analyzed based on their enzyme structures. A method of systematic comparison and classification of the catalytic reactions was developed. Hydrolysis and transfer reactions closely resemble one another, displaying common mechanisms, single displacement, and double displacement. These displacement mechanisms might be further subclassified according to the type of catalytic factors and nucleophilic substitution involved. Several types of catalytic factors exist: nucleophile, acid, base, stabilizer, modulator, cofactors. Nucleophilic substitution might be categorized as S(N)1/S(N)2 (or dissociative/associative) reactions. The classification indicates that some mechanisms favor particular types of catalytic factors. In hydrolyses of amide bonds and phosphoric ester bonds, mechanisms with single displacement tend to use inorganic cofactors such as zinc and magnesium ions as important catalysts, whereas those with double displacement frequently do not use such cofactors. In contrast, hydrolyses of O-glycoside bond rarely use such cofactors, with one exception. The trypsin-like hydrolytic reaction, which is catalyzed by the classic catalytic triad comprising serine/histidine/aspartate, can be considered as a "super-reaction" because it is observed in at least three nonhomologous enzymes, whereas most reactions are singlets without any nonhomologous enzymes. By dividing complex reactions into several reactions, correlations between active site structures and catalytic functions can be suggested. This classification method is applicable to other reactions such as elimination and isomerization. Furthermore, it will facilitate annotation of enzyme functions from 3D patterns of enzyme active sites. The classification is available at http://mbs.cbrc.jp/EzCatDB/RLCP/index.html.     

Carbon–nitrogen stoichiometry in the tritrophic food chain willow, leaf beetle, and predatory ladybird beetle

Although plant quality can indirectly increase the performance of the third trophic level by bottom-up cascading effects, the mechanisms of this indirect effect are still unclear. In this study the carbon–nitrogen stoichiometry in a tri-trophic system consisting of the willow, a leaf beetle, and a predatory ladybird beetle were examined to determine the mechanisms of the bottom-up cascading effect. The bottom-up cascade is initiated by increasing leaf nitrogen, because of artificial cutting of willow trees. The relative growth rate (RGR) of the leaf beetle increased when fed on cut willow leaves, because of the high leaf nitrogen in the cut willows. Ladybird beetle RGR also increased when fed on leaf beetles fed on cut willow leaves. The increased RGR of the ladybirds cannot be explained by the quality of the prey, however, because leaf beetle nitrogen was not affected by host plant quality. Thus, the carbon–nitrogen stoichiometry could not be a mechanism of the bottom-up cascade through multiple trophic levels.     

Environmental decision support systems: A new approach to support the operation and maintenance of horizontal subsurface flow constructed wetlands

This paper describes the development and operation of an Environmental Decision Support System (EDSS) to improve the operation and maintenance of horizontal subsurface-flow constructed wetlands (EDSS-maintenance). Constructed wetlands (CWs) allow wastewater treatment in a sustainable manner since they involve low energy consumption, low construction and functioning costs and low environmental impact. However, operation and maintenance activities are essential to guarantee reliability in CWs performance. The definition of operation and maintenance protocols depends on several quantitative and qualitative aspects such as wastewater treatment plant configuration, CW design, influent characteristics, sensitivity of the receiving media, etc. Bearing this in mind and considering the limited technical knowledge about CWs, the need for a new tool to support CW performance is clear. In this sense, EDSSs offer a new approach because they can tackle problems of complex and uncertain systems. The EDSS-maintenance provides operation and maintenance manuals specifically defined for every CW. To achieve it, the required knowledge was implemented within a rule-based system, which forms the backbone of the EDSS. Several features presented in this paper demonstrate how the EDSS-maintenance provides a proper platform to support the necessary collaborative work in the ecological engineering problem of horizontal subsurface flow CWs operation and maintenance.     

Ecology of testate amoebae from mires in the Central Rhodope Mountains, Greece and development of a transfer function for palaeohydrological reconstruction

Testate amoebae are useful environmental indicators in ecological and palaeoecological studies from peatlands. Previous quantitative studies have focused on the peatlands of Northern and Central Europe, North America, and New Zealand and have considered a relatively restricted variety of peatland types, mostly ombrotrophic or Sphagnum-dominated while more minerotrophic fens have been less studied. Here we present the first quantitative ecological study of testate amoebae from four small mesotrophic fens (pH 5.5-8.1) in the Elatia Forest, northern Macedonia province, Greece. Relationships with the environmental data were investigated using redundancy analysis and mantel tests. Transfer function models were derived using a variety of techniques. Results demonstrate that as for Sphagnum-dominated mires hydrology is the most important control on amoebae community composition. Transfer function models should enable water tables to be predicted within 2.5 cm, when data selection is used this is reduced to less than 2 cm. pH is also an important environmental control on testate amoebae communities, a transfer function model enables pH prediction within 0.4 pH units. The hydrological transfer function is the best performing such model yet produced in terms of prediction error. This study provides new data on the ecology of testate amoebae in fens, and the transfer function models should allow quantitative palaeohydrological reconstruction.     

Testing the growth rate vs. geochemical hypothesis for latitudinal variation in plant nutrients

Two hypotheses have been proposed to explain increases in plant nitrogen (N) and phosphorus (P) concentrations with latitude: (i) geochemical limitation to P availability in the tropics and (ii) temperature driven variation in growth rate, where greater growth rates (requiring greater nutrient levels) are needed to complete growth and reproduction within shorter growing seasons in temperate than tropical climates. These two hypotheses were assessed in one forest type, intertidal mangroves, using fertilized plots at sites between latitudes 36º S and 27º N. The N and P concentrations in mangrove leaf tissue increased with latitude, but there were no trends in N : P ratios. Growth rates of trees, adjusted for average minimum temperature showed a significant increase with latitude supporting the Growth Rate Hypothesis. However, support for the Geochemical Hypothesis was also strong; both photosynthetic P use efficiency and nutrient resorption efficiency decreased with increasing latitude, indicating that P was less limiting to metabolism at the higher latitudes. Our study supports the hypothesis that historically low P availability in the tropics has been an important selective pressure shaping the evolution of plant traits.     

Disease mapping and spatial regression with count data

In this paper, we provide critical reviews of methods suggested for the analysis of aggregate count data in the context of disease mapping and spatial regression. We introduce a new method for picking prior distributions, and propose a number of refinements of previously used models. We also consider ecological bias, mutual standardization, and choice of both spatial model and prior specification. We analyze male lip cancer incidence data collected in Scotland over the period 1975-1980, and outline a number of problems with previous analyses of these data. In disease mapping studies, hierarchical models can provide robust estimation of area-level risk parameters, though care is required in the choice of covariate model, and it is important to assess the sensitivity of estimates to the spatial model chosen, and to the prior specifications on the variance parameters. Spatial ecological regression is a far more hazardous enterprise for two reasons. First, there is always the possibility of ecological bias, and this can only be alleviated by the inclusion of individual-level data. For the Scottish data, we show that the previously used mean model has limited interpretation from an individual perspective. Second, when residual spatial dependence is modeled, and if the exposure has spatial structure, then estimates of exposure association parameters will change when compared with those obtained from the independence across space model, and the data alone cannot choose the form and extent of spatial correlation that is appropriate.     

Disturbance persistence in managed grasslands: shifts in aboveground community structure and the weed seed bank

The length of time and form in which disturbances persist in systems depends on the intensity and frequency of disturbance and on the abilities of resident species to recover from such events. In grazed grasslands, trampling by large mammalian herbivores can periodically facilitate weed establishment by exposing patches of bare ground but whether an intense soil disturbance event results in a temporary increase in weed abundance or a persistent weed problem remains unclear. In May 2002, cattle trampling following heavy rain caused severe damage to nine-month old, rotationally grazed, cool-season pastures (Midwest USA). In September 2002, we compared the aboveground composition of paddocks (i.e., fenced pasture sections) that were heavily disturbed to those that received no damage. Relative to undisturbed paddocks, forage species relative cover was 17% lower in disturbed paddocks, and weed species and bare ground relative cover was 61% and 100% higher, respectively. By September 2004, paddock types did not differ in all aboveground community components. However, the abundance and species richness of weed seeds in the soil seed bank averaged respectively 82% and 30% higher in disturbed paddocks between 2003 and 2004. These findings indicate that a spatially extensive, intense soil disturbance event may soon become undetectable in components of aboveground pasture structure but can persist as an augmented weed seed bank. Because of high weed seed bank longevity, disturbances to formerly disturbed pastures would likely result in higher weed recruitment, with more species represented, than in those which lack previous disturbance. Disturbance history may thus be a useful predictor of weed community composition following subsequent disturbance. Based on empirical data supporting this proposition, we recommend that grassland managers explicitly incorporate disturbance history into dynamic management planning and do not rely exclusively on aboveground characters to evaluate the invasion status or colonization potential of an area by undesirable plants. We emphasize that the ecological legacies of past soil disturbance events cannot only influence the contemporary patterns and processes of grasslands, but importantly, affect their compositional trajectories following subsequent perturbation.     

Schistosoma japonicum translationally controlled tumour protein,which is associated with the development of female worms,as a target for control of schistosomiasis

Schistosomiasis is a globally important helminthic disease of both humans and animals, and is the second most common parasitic disease after malaria. Although praziquantel is extensively used for treatment of parasitic diseases, drug resistance has been reported. Therefore, new drugs and effective vaccines are needed for continuous control of schistosomiasis. Eggs produced by schistosomes are responsible for the occurrence and spread of schistosomiasis. Revealing the reproductive mechanism of schistosomes will help to control this disease. In this study, the proteomic profiles of single-sex infected female worms and bisexual infected mature female worms of Schistosoma japonicum at 18, 21, 23 and 25 days p.i. were identified with isobaric tags for relative quantitation-coupled liquid chromatography–tandem mass spectrometry. Differentially expressed proteins were subsequently used for bioinformatic analysis. Six highly expressed differentially expressed proteins in mature female worms were selected and long-term interference with small interfering RNA (siRNA) was conducted to determine biological functions. SiRNA against S. japonicum translationally controlled tumour protein (SjTCTP) resulted in the most significant effect on the growth and development of MF worms. Sjtctp mRNA expression gradually increased over time with a high level of expression maintained at 25–42 days p.i., while levels were significantly higher in mature female worms than male and SF worms. The subsequent animal immune protection experiments showed that recombinant SjTCTP (rSjTCTP) reduced the number of adults by 44.7% (P < 0.01), average egg burden per gram of liver by 57.94% (P < 0.01), egg hatching rate by 47.57% (P < 0.01), and oviposition of individual females by 43.16%. rSjTCTP induced higher levels of serum IgG, IL-2, and IL-10 in mice. Collectively, these results show that SjTCTP is vital to reproduction of female worms and, thus, is a candidate antigen for immune protection.     

When cooperation begets cooperation: the role of key individuals in galvanizing support

Life abounds with examples of conspecifics actively cooperating to a common end, despite conflicts of interest being expected concerning how much each individual should contribute. Mathematical models typically find that such conflict can be resolved by partial-response strategies, leading investors to contribute relatively equitably. Using a case study approach, we show that such model expectations can be contradicted in at least four disparate contexts: (i) bi-parental care; (ii) cooperative breeding; (iii) cooperative hunting; and (iv) human cooperation. We highlight that: (a) marked variation in contributions is commonplace; and (b) individuals can often respond positively rather than negatively to the contributions of others. Existing models have surprisingly limited power in explaining these phenomena. Here, we propose that, although among-individual variation in cooperative contributions will be influenced by differential costs and benefits, there is likely to be a strong genetic or epigenetic component. We then suggest that selection can maintain high investors (key individuals) when their contributions promote support by increasing the benefits and/or reducing the costs for others. Our intentions are to raise awareness in—and provide testable hypotheses of—two of the most poorly understood, yet integral, questions regarding cooperative ventures: why do individuals vary in their contributions and when does cooperation beget cooperation?