Statistical significance of protein structure prediction by threading

In this study, we estimate the statistical significance of structure prediction by threading. We introduce a single parameter ɛ that serves as a universal measure determining the probability that the best alignment is indeed a native-like analog. Parameter ɛ takes into account both length and composition of the query sequence and the number of decoys in threading simulation. It can be computed directly from the query sequence and potential of interactions, eliminating the need for sequence reshuffling and realignment. Although our theoretical analysis is general, here we compare its predictions with the results of gapless threading. Finally we estimate the number of decoys from which the native structure can be found by existing potentials of interactions. We discuss how this analysis can be extended to determine the optimal gap penalties for any sequence-structure alignment (threading) method, thus optimizing it to maximum possible performance.

Local densities orthogonal to beta-sheet amide planes: patterns of packing in globular proteins.

We have investigated the efficiency of packing by calculating intramolecular packing density above and below peptide planes of internal beta-pleated sheet residues in five globular proteins. The orientation of interest was chosen to allow study of regions that are approximately perpendicular to the faces of beta-pleated sheets. In these locations, nonbonded van der Waals packing interactions predominate over hydrogen bonding and solvent interactions. We observed considerable variability in packing densities within these regions, confirming that the interior packing of a protein does not result in uniform occupation of the available space. Patterns of fluctuation in packing density suggest that the regular backbone-to-backbone network of hydrogen bonds is not likely to be interrupted to maximize van der Waals interactions. However, high-density packing tends to occur toward the ends of beta-structure strands where hydrogen bonds are more likely to involve nonpolar side-chain groups or solvent molecules. These features result in internal protein folding with a central low-density core surrounded by a higher-density subsurface shell, consistent with our previous calculations regarding overall protein packing density.

Extraordinarily high spider densities on islands: flow of energy from the marine to terrestrial food webs and the absence of predation.

Some islands in the Gulf of California support very high densities of spiders. Spider density is negatively correlated with island size; many small islands support 50-200 spiders per m3 of cactus. Energy for these spiders comes primarily from the ocean and not from in situ productivity by land plants. We explicitly connect the marine and terrestrial systems to show that insular food webs represent one endpoint of the marine web. We describe two conduits for marine energy entering these islands: shore drift and seabird colonies. Both conduits are related to island area, having a much stronger effect on smaller islands. This asymmetric effect helps to explain the exceptionally high spider densities on small islands. Although productivity sets the maximal potential densities, predation (by scorpions) limits realized spider abundance. Thus, prey availability and predation act in concert to set insular spider abundance.