Protein Docking by the Underestimation of Free Energy Funnels in the Space of Encounter Complexes

Similarly to protein folding, the association of two proteins is driven by a free energy funnel, determined by favorable interactions in some neighborhood of the native state. We describe a docking method based on stochastic global minimization of funnel-shaped energy functions in the space of rigid body motions (SE(3)) while accounting for flexibility of the interface side chains. The method, called semi-definite programming-based underestimation (SDU), employs a general quadratic function to underestimate a set of local energy minima and uses the resulting underestimator to bias further sampling. While SDU effectively minimizes functions with funnel-shaped basins, its application to docking in the rotational and translational space SE(3) is not straightforward due to the geometry of that space. We introduce a strategy that uses separate independent variables for side-chain optimization, center-to-center distance of the two proteins, and five angular descriptors of the relative orientations of the molecules. The removal of the center-to-center distance turns out to vastly improve the efficiency of the search, because the five-dimensional space now exhibits a well-behaved energy surface suitable for underestimation. This algorithm explores the free energy surface spanned by encounter complexes that correspond to local free energy minima and shows similarity to the model of macromolecular association that proceeds through a series of collisions. Results for standard protein docking benchmarks establish that in this space the free energy landscape is a funnel in a reasonably broad neighborhood of the native state and that the SDU strategy can generate docking predictions with less than 5 � ligand interface Ca root-mean-square deviation while achieving an approximately 20-fold efficiency gain compared to Monte Carlo methods.

Quantitative Analysis of Single Nucleotide Polymorphisms within Copy Number Variation

Background Single nucleotide polymorphisms (SNPs) have been used extensively in genetics and epidemiology studies. Traditionally, SNPs that did not pass the Hardy-Weinberg equilibrium (HWE) test were excluded from these analyses. Many investigators have addressed possible causes for departure from HWE, including genotyping errors, population admixture and segmental duplication. Recent large-scale surveys have revealed abundant structural variations in the human genome, including copy number variations (CNVs). This suggests that a significant number of SNPs must be within these regions, which may cause deviation from HWE. Results We performed a Bayesian analysis on the potential effect of copy number variation, segmental duplication and genotyping errors on the behavior of SNPs. Our results suggest that copy number variation is a major factor of HWE violation for SNPs with a small minor allele frequency, when the sample size is large and the genotyping error rate is 0~1%. Conclusions Our study provides the posterior probability that a SNP falls in a CNV or a segmental duplication, given the observed allele frequency of the SNP, sample size and the significance level of HWE testing.

Nutrition or Detoxification: Why Bats Visit Mineral Licks of the Amazonian Rainforest

Many animals in the tropics of Africa, Asia and South America regularly visit so-called salt or mineral licks to consume clay or drink clay-saturated water. Whether this behavior is used to supplement diets with locally limited nutrients or to buffer the effects of toxic secondary plant compounds remains unclear. In the Amazonian rainforest, pregnant and lactating bats are frequently observed and captured at mineral licks. We measured the nitrogen isotope ratio in wing tissue of omnivorous short-tailed fruit bats, Carollia perspicillata, and in an obligate fruit-eating bat, Artibeus obscurus, captured at mineral licks and at control sites in the rainforest. Carollia perspicillata with a plant-dominated diet were more often captured at mineral licks than individuals with an insect-dominated diet, although insects were more mineral depleted than fruits. In contrast, nitrogen isotope ratios of A. obscurus did not differ between individuals captured at mineral lick versus control sites. We conclude that pregnant and lactating fruit-eating bats do not visit mineral licks principally for minerals, but instead to buffer the effects of secondary plant compounds that they ingest in large quantities during periods of high energy demand. These findings have potential implications for the role of mineral licks for mammals in general, including humans.

Patient Retention in Antiretroviral Therapy Programs in Sub-Saharan Africa: A Systematic Review

Paper published in PLoS Medicine in 2007.