The role of deer as vehicles to move ticks, Ixodes ricinus, between contrasting habitats

In Europe the most important hosts maintaining Ixodes ricinus tick populations are deer. Therefore, excluding deer by fencing or culling are potential tick management tools. Here we test the hypothesis that deer act as vehicles for moving ticks between two distinct habitats: forest and open heather moorland. We utilised an ideal “natural experiment” whereby forests were either fenced or unfenced to prevent or allow deer to move between habitats. We aimed to test the hypothesis that deer cause a net movement of ticks from high tick density areas, i.e. forests, to low tick density areas, i.e. open moorland. We recorded I. ricinus and host abundance in 10 unfenced and seven fenced forests and their respective surrounding heather moorland. We found that fenced forests had fewer deer and fewer I. ricinus nymphs than unfenced forests. However, we found no evidence that fencing forests reduced I. ricinus abundance on adjacent heather moorland. Thus there was insufficient evidence for our hypothesis that deer cause a net movement of ticks from forest onto adjacent moorland. However, we found that deer abundance generally correlates with I. ricinus abundance. We conclude that fencing can be used as a tool to reduce ticks and disease risk in forests, but that fencing forests is unlikely to reduce ticks or disease risk on adjacent moorland. Instead, reducing deer numbers could be a potential tool to reduce tick abundance with implications for disease mitigation.     

Differences in beech (Fagus crenata) regeneration between two types of Japanese beech forest and along a snow gradient

Differences in beech (Fagus crenata) regeneration were quantitatively investigated using power function analysis for the size–class (diameter at breast height, DBH) distribution and juvenile-to-canopy tree (J/C) ratio along a snow gradient throughout Japan. In snowy areas, all species combined, as well as F. crenata alone, showed constant regeneration, with parameter b≈−1.6 for the power function y=ax ^b (x=DBH, y=density), which is related to the DBH–class distribution. The good fit of the data to the function suggests that beech regenerates constantly with self-thinning patch dynamics. Parameter a, which indicates the abundance of small trunks, was high. Furthermore, the mean J/C ratio was ≈8, i.e., each parent beech tree produced eight juveniles. These results suggest that beech regenerates constantly with gap dynamics in snowy beech forests on the Japan Sea side of Japan (snowy). However, the fit of F. crenata was lower and nonsignificant in some forests in less snowy areas, despite the high fit of all species combined. In these areas, the mean of a was low, and b was often near zero for F. crenata regressions. These results suggest that the abundance of beech was low, and self-thinning was not evident because of the initial low abundance. Moreover, the mean J/C ratio was <1.0, suggesting that juvenile density was lower than that of canopy trees. Thus, the regeneration of F. crenata on the Pacific Ocean side of Japan (less snowy) is rather sporadic. Less snowy conditions may promote seed desiccation, predation of beechnuts and seedlings, and water stress. Lower F. crenata density may also reduce predator satiation and wind pollination.     

Estimation of amino acid residue substitution rates at local spatial regions and application in protein function inference: a Bayesian Monte Carlo approach

The amino acid sequences of proteins provide rich information for inferring distant phylogenetic relationships and for predicting protein functions. Estimating the rate matrix of residue substitutions from amino acid sequences is also important because the rate matrix can be used to develop scoring matrices for sequence alignment. Here we use a continuous time Markov process to model the substitution rates of residues and develop a Bayesian Markov chain Monte Carlo method for rate estimation. We validate our method using simulated artificial protein sequences. Because different local regions such as binding surfaces and the protein interior core experience different selection pressures due to functional or stability constraints, we use our method to estimate the substitution rates of local regions. Our results show that the substitution rates are very different for residues in the buried core and residues on the solvent-exposed surfaces. In addition, the rest of the proteins on the binding surfaces also have very different substitution rates from residues. Based on these findings, we further develop a method for protein function prediction by surface matching using scoring matrices derived from estimated substitution rates for residues located on the binding surfaces. We show with examples that our method is effective in identifying functionally related proteins that have overall low sequence identity, a task known to be very challenging.     

Role of disulfide bonds in modulating internal motions of proteins to tune their function: molecular dynamics simulation of scorpion toxin Lqh III

A series of 1-ns MD simulations were performed on the scorpion toxin Lqh III in native and disulfide bond broken states. The removal of disulfide bonds has caused hydrogen bond network alteration in the five-residue turn, the long loop, the alpha-helix, the loop connecting strands II and III, and the C-terminal region. In addition and more importantly, it has influenced the amplitude of the fluctuations of five-residue turn, loops, and C-terminal region with a minor effect on the fluctuations of the cysteines in the broken bond sites. These findings suggest that disulfide bonds are not the most important factors in rigidifying their own locations, while they have more important effects at a global scale. Furthermore, our results reveal that disulfide bonds have considerable influence on the functionally important essential modes of motions and the correlations between the motions of the binding site residues. Therefore, we can conclude that disulfide bonds have a crucial role in modulating the function via adjusting the dynamics of scorpion toxin molecules. Although this conclusion cannot be generalized to all peptides and proteins, it demonstrates the importance of more investigations on this aspect of disulfide bond efficacy.     

Differentiation-dependent expression of hypothetical proteins in the neuroblastoma cell line N1E-115

Several protein cascades, including signaling, cytoskeletal, chaperones, metabolic, and antioxidant proteins, have been shown to be involved in the process of neuronal differentiation (ND) of neuroblastoma cell lines. No systematic approach to detect hitherto unknown and unnamed proteins or structures that have been predicted upon nucleic acid sequences in ND has been published so far. We therefore decided to screen hypothetical protein (HP) expression by protein profiling. Two-dimensional gel electrophoresis with subsequent matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF/TOF) identification was used for expression analysis of undifferentiated and dimethylsulfoxide-induced neuronally differentiated N1E-115 cells. We unambiguously identified six HPs: Q8C520, Q99LF4, Q9CXS1, Q9DAF8, Q91WT0, and Q8C5G2. A prefoldin domain in Q91WT0, a t-SNARE domain in Q9CXS1, and a bromodomain were observed in Q8C5G2. For the three remaining proteins, no putative function using Pfam, BLOCKS, PROSITE, PRINTS, InterPro, Superfamily, CoPS, and ExPASy could be assigned. While two proteins were present in both cell lines, Q9CXS1 was switched off (i.e., undetectably low) in differentiated cells only, and Q9DAF8, Q91WT0, and Q8C5G2 were switched on in differentiated cells exclusively. Herein, using a proteomic approach suitable for screening and identification of HP, we present HP structures that have been only predicted so far based upon nucleic acid sequences. The four differentially regulated HPs may play a putative role in the process of ND.     

Assessing protein disorder and induced folding

Intrinsically disordered proteins (IDPs) defy the structure-function paradigm as they fulfill essential biological functions while lacking well-defined secondary and tertiary structures. Conformational and spectroscopic analyses showed that IDPs do not constitute a uniform family, and can be divided into subfamilies as a function of their residual structure content. Residual intramolecular interactions are thought to facilitate binding to a partner and then induced folding. Comprehensive information about experimental approaches to investigate structural disorder and induced folding is still scarce. We herein provide hints to readily recognize features typical of intrinsic disorder and review the principal techniques to assess structural disorder and induced folding. We describe their theoretical principles and discuss their respective advantages and limitations. Finally, we point out the necessity of using different approaches and show how information can be broadened by the use of multiples techniques.     

On the nature of cavities on protein surfaces: application to the identification of drug-binding sites

In this article we introduce a new method for the identification and the accurate characterization of protein surface cavities. The method is encoded in the program SCREEN (Surface Cavity REcognition and EvaluatioN). As a first test of the utility of our approach we used SCREEN to locate and analyze the surface cavities of a nonredundant set of 99 proteins cocrystallized with drugs. We find that this set of proteins has on average about 14 distinct cavities per protein. In all cases, a drug is bound at one (and sometimes more than one) of these cavities. Using cavity size alone as a criterion for predicting drug-binding sites yields a high balanced error rate of 15.7%, with only 71.7% coverage. Here we characterize each surface cavity by computing a comprehensive set of 408 physicochemical, structural, and geometric attributes. By applying modern machine learning techniques (Random Forests) we were able to develop a classifier that can identify drug-binding cavities with a balanced error rate of 7.2% and coverage of 88.9%. Only 18 of the 408 cavity attributes had a statistically significant role in the prediction. Of these 18 important attributes, almost all involved size and shape rather than physicochemical properties of the surface cavity. The implications of these results are discussed. A SCREEN Web server is available at http://interface.bioc.columbia.edu/screen.