Enhancing systematic protein-protein docking methods using ray casting: application to ATTRACT

Systematic protein-protein docking methods need to evaluate a huge number of different probe configurations, thus leading to high computational cost. We present an efficient filter-ray casting filter (RCF)-that enables a notable speed-up of systematic protein-protein docking. The high efficiency of RCF is the outcome of the following factors: (i) extracting of pockets and protrusions on the surfaces of the proteins using visibilities; (ii) a ray casting method that finds aligned receptor pocket/probe protrusion pairs without explicit similarity computations. The RCF method enables the integration of systematic methods and local shape feature matching methods. To verify the efficiency and the accuracy of RCF, we integrated it with a systematic protein-protein docking approach (ATTRACT) based on a reduced protein representation. The test results show that the integrated docking approach is much faster. At the same time, it ranks the lowest ligand root-mean-square deviation (RMSD) (L_rms) solutions higher when docking enzyme-enzyme inhibitor complexes. Consequently, RCF not only enables much faster execution of systematic docking runs but also improves the qualities of docking predictions.     

Extinction in relation to demographic and environmental stochasticity in age-structured models

The demographic variance of an age-structured population is defined. This parameter is further split into components generated by demographic stochasticity in each vital rate. The applicability of these parameters are investigated by checking how an age-structured population process can be approximated by a diffusion with only three parameters. These are the deterministic growth rate computed from the expected projection matrix and the environmental and demographic variances. We also consider age-structured populations where the fecundity at any stage is either zero or one, and there is neither environmental stochasticity nor dependence between individual fecundity and survival. In this case the demographic variance is uniquely determined by the vital rates defining the projection matrix. The demographic variance for a long-lived bird species, the wandering albatross in the southwestern part of the Indian Ocean, is estimated. We also compute estimates of the age-specific contributions to the total demographic variance from survival, fecundity and the covariance between survival and fecundity.     

Viability of the Endangered Little Bustard <Emphasis Type="Italic">Tetrax tetrax</Emphasis> Population of Western France

Stochastic computer simulations are used to evaluate the sensitivity of Little bustard population parameters, estimating the survival probabilities of the seven endangered Little bustard populations of central-western France for which conservation actions are currentlybeing or have been implemented. Different scenarios of parameter compensation for those nuclei to establish parameter levels assuring population viability are discussed. Adult survival, productivity per female, initial population size and carrying capacity were the most sensitive parameters in a hypothetical, isolated population. Juvenile survival also affected population survival, although its sensitivity was lower. Sex ratio did not have a linear effect on population survival, but probability of extinction increased for extreme values. Productivity per female and initial population size, varied strongly among the populations studied, determining their average time of extinction and growth rate. When a metapopulation scenario was simulated, the survival probabilities of each population and the metapopulations stayed close to 1.0 if no mortality was associated to migration. When mortality during migration was included in the simulations, the metapopulation's probability of survival significantly decreased under 90%. This approach may help managers to correctly address conservation measures and design effective strategies, which should be directed mainly to improve productivity, enhance female survival, and minimise mortality during migration (e.g. promoting insect-rich nesting substrates, avoiding female killing and nest destruction at harvesting, reducing the risk of collision with powerlines, or controlling poaching).     

Bayesian inference for stochastic kinetic models using a diffusion approximation

This article is concerned with the Bayesian estimation of stochastic rate constants in the context of dynamic models of intracellular processes. The underlying discrete stochastic kinetic model is replaced by a diffusion approximation (or stochastic differential equation approach) where a white noise term models stochastic behavior and the model is identified using equispaced time course data. The estimation framework involves the introduction of m- 1 latent data points between every pair of observations. MCMC methods are then used to sample the posterior distribution of the latent process and the model parameters. The methodology is applied to the estimation of parameters in a prokaryotic autoregulatory gene network.     

Theoretical study of volume changes associated with the helix-coil transition of an alanine-rich peptide in aqueous solution

The changes in the partial molar volume (PMV) associated with the conformational transition of an alanine-rich peptide AK16 from the alpha-helix structure to various random coil structures are calculated by the three-dimensional interaction site model (3D-RISM) theory coupled with the Kirkwood-Buff theory. The volume change is analyzed by decomposing it into contributions from geometry and hydration: the changes in the van der Waals, void, thermal, and interaction volume. The total change in the PMV is positive. This is primarily due to the growth of void space within the peptide, which is canceled in part by the volume reduction resulting from the increase in the electrostatic interaction between the peptide and water molecules. The changes in the void and thermal volume of the coil structures are widely distributed and tend to compensate each other. Additionally, the relations between the hydration volume components and the surface properties are investigated. We categorize coil structures into extended coils with the PMV smaller than helix and general coils with the PMV larger than helix. The pressure therefore can both stabilize and destabilize the coil structures. The latter seems to be a more proper model of random coil structures of the peptide.     

Investigating the form-function interface in African apes: Relationships between principal moments of area and positional behaviors in femoral and humeral diaphyses

Investigations of cross-sectional geometry in nonhuman primate limb bones typically attribute shape ratios to qualitative behavioral characterizations, e.g., leaper, slow climber, brachiator, or terrestrial vs. arboreal quadruped. Quantitative positional behavioral data, however, have yet to be used in a rigorous evaluation of such shape-behavior connections. African apes represent an ideal population for such an investigation because their relatedness minimizes phylogenetic inertia, they exhibit diverse behavioral repertoires, and their locomotor behaviors are known from multiple studies. Cross-sectional data from femoral and humeral diaphyses were collected for 222 wild-shot specimens, encompassing Pan paniscus and all commonly recognized African ape subspecies. Digital representations of diaphyseal cross sections were acquired via computed tomography at three locations per diaphysis. Locomotor behaviors were pooled broadly into arboreal and terrestrial categories, then partitioned into quadrupedal walking, quadrumanous climbing, scrambling, and suspensory categories. Sex-specific taxonomic differences in ratios of principal moments of area (PMA) were statistically significant more often in the femoral diaphysis than the humeral diaphysis. While it appears difficult to relate a measure of shape (e.g., PMA ratio) to individual locomotor modes, general locomotor differences (e.g., percentage arboreal vs. terrestrial locomotion) are discerned more easily. As percentage of arboreal locomotion for a group increases, average cross sections appear more circular. Associations between PMA ratio and specific locomotor behaviors are less straightforward. Individual behaviors that integrate eccentric limb positions (e.g., arboreal scrambling) may not engender more circular cross sections than behaviors that incorporate repetitive sagittal movements (e.g., quadrupedal walking) in a straightforward manner.