Simultaneous Viral Exposure and Protection from Neanderthal Introgression.

In this issue, Enard and Petrov present intriguing results on the possibility of genetic traces left behind in our genomes from adaptation to past viral epidemics that may have been initiated by interaction with Neanderthal archaic hominins. The work highlights how powerful infectious agents can act as a selective force to shape our genetic makeup.

Discovering anti-platelet drug combinations with an integrated model of activator-inhibitor relationships, activator-activator synergies and inhibitor-inhibitor synergies.

Identifying effective therapeutic drug combinations that modulate complex signaling pathways in platelets is central to the advancement of effective anti-thrombotic therapies. However, there is no systems model of the platelet that predicts responses to different inhibitor combinations. We developed an approach which goes beyond current inhibitor-inhibitor combination screening to efficiently consider other signaling aspects that may give insights into the behaviour of the platelet as a system. We investigated combinations of platelet inhibitors and activators. We evaluated three distinct strands of information, namely: activator-inhibitor combination screens (testing a panel of inhibitors against a panel of activators); inhibitor-inhibitor synergy screens; and activator-activator synergy screens. We demonstrated how these analyses may be efficiently performed, both experimentally and computationally, to identify particular combinations of most interest. Robust tests of activator-activator synergy and of inhibitor-inhibitor synergy required combinations to show significant excesses over the double doses of each component. Modeling identified multiple effects of an inhibitor of the P2Y12 ADP receptor, and complementarity between inhibitor-inhibitor synergy effects and activator-inhibitor combination effects. This approach accelerates the mapping of combination effects of compounds to develop combinations that may be therapeutically beneficial. We integrated the three information sources into a unified model that predicted the benefits of a triple drug combination targeting ADP, thromboxane and thrombin signaling.

Alternate virtual populations elucidate the type I interferon signature predictive of the response to rituximab in rheumatoid arthritis.

BACKGROUND: Mechanistic biosimulation can be used in drug development to form testable hypotheses, develop predictions of efficacy before clinical trial results are available, and elucidate clinical response to therapy. However, there is a lack of tools to simultaneously (1) calibrate the prevalence of mechanistically distinct, large sets of virtual patients so their simulated responses statistically match phenotypic variability reported in published clinical trial outcomes, and (2) explore alternate hypotheses of those prevalence weightings to reflect underlying uncertainty in population biology. Here, we report the development of an algorithm, MAPEL (Mechanistic Axes Population Ensemble Linkage), which utilizes a mechanistically-based weighting method to match clinical trial statistics. MAPEL is the first algorithm for developing weighted virtual populations based on biosimulation results that enables the rapid development of an ensemble of alternate virtual population hypotheses, each validated by a composite goodness-of-fit criterion. RESULTS: Virtual patient cohort mechanistic biosimulation results were successfully calibrated with an acceptable composite goodness-of-fit to clinical populations across multiple therapeutic interventions. The resulting virtual populations were employed to investigate the mechanistic underpinnings of variations in the response to rituximab. A comparison between virtual populations with a strong or weak American College of Rheumatology (ACR) score in response to rituximab suggested that interferon ß (IFNß) was an important mechanistic contributor to the disease state, a signature that has previously been identified though the underlying mechanisms remain unclear. Sensitivity analysis elucidated key anti-inflammatory properties of IFNß that modulated the pathophysiologic state, consistent with the observed prognostic correlation of baseline type I interferon measurements with clinical response. Specifically, the effects of IFNß on proliferation of fibroblast-like synoviocytes and interleukin-10 synthesis in macrophages each partially counteract reductions in synovial inflammation imparted by rituximab. A multianalyte biomarker panel predictive for virtual population therapeutic responses suggested population dependencies on B cell-dependent mediators as well as additional markers implicating fibroblast-like synoviocytes. CONCLUSIONS: The results illustrate how the MAPEL algorithm can leverage knowledge of cellular and molecular function through biosimulation to propose clear mechanistic hypotheses for differences in clinical populations. Furthermore, MAPEL facilitates the development of multianalyte biomarkers prognostic of patient responses in silico.

Optimal stepwise experimental design for pairwise functional interaction studies.

MOTIVATION: Pairwise experimental perturbation is increasingly used to probe gene and protein function because these studies offer powerful insight into the activity and regulation of biological systems. Symmetric two-dimensional datasets, such as pairwise genetic interactions are amenable to an optimally designed measurement procedure because of the equivalence of cases and conditions where fewer experimental measurements may be required to extract the underlying structure. RESULTS: We show that optimal experimental design can provide improvements in efficiency when collecting data in an iterative manner. We develop a method built on a statistical clustering model for symmetric data and the Fisher information uncertainty estimates, and we also provide simple heuristic approaches that have comparable performance. Using yeast epistatic miniarrays as an example, we show that correct assignment of the major subnetworks could be achieved with <50% of the measurements in the complete dataset. Optimization is likely to become critical as pairwise functional studies extend to more complex mammalian systems where all by all experiments are currently intractable.

Discovery of small molecule inhibitors of protein-protein interactions using combined ligand and target score normalization.

Docking experiments of multiple compounds typically focus on a single protein. However, other targets provide information about relative binding efficiencies that is otherwise lacking. We developed a docking strategy that normalized results in both the ligand and target dimensions. This was applied to dock 287 approved small drugs with 35 peptide-binding proteins, including 15 true positives. The combined docking score was normalized by drug and protein and by incorporating information on contact similarity to the template protein-peptide contacts. The 20 top ranking hits included 6 true positives, and three matches with suggestive evidence in the literature: the cardiac glycoside digitoxin may inhibit WW domain interactions, the 14-3-3 zeta protein may bind negatively charged ligands, and the nuclear receptor coactivator site may bind nuclear receptor agonists. Additionally, the Bcl-2 antiapoptotic protein is predicted to bind pargyline, and the antiapoptic p53 interacting protein MDM2 is suggested to bind clofazimine. These predictions represent starting points for the experimental development of PPI inhibitors based on an existing database of approved drugs and demonstrate that two-dimensional normalization improves docking efficiency.

Universally sloppy parameter sensitivities in systems biology models.

Quantitative computational models play an increasingly important role in modern biology. Such models typically involve many free parameters, and assigning their values is often a substantial obstacle to model development. Directly measuring in vivo biochemical parameters is difficult, and collectively fitting them to other experimental data often yields large parameter uncertainties. Nevertheless, in earlier work we showed in a growth-factor-signaling model that collective fitting could yield well-constrained predictions, even when it left individual parameters very poorly constrained. We also showed that the model had a "sloppy" spectrum of parameter sensitivities, with eigenvalues roughly evenly distributed over many decades. Here we use a collection of models from the literature to test whether such sloppy spectra are common in systems biology. Strikingly, we find that every model we examine has a sloppy spectrum of sensitivities. We also test several consequences of this sloppiness for building predictive models. In particular, sloppiness suggests that collective fits to even large amounts of ideal time-series data will often leave many parameters poorly constrained. Tests over our model collection are consistent with this suggestion. This difficulty with collective fits may seem to argue for direct parameter measurements, but sloppiness also implies that such measurements must be formidably precise and complete to usefully constrain many model predictions. We confirm this implication in our growth-factor-signaling model. Our results suggest that sloppy sensitivity spectra are universal in systems biology models. The prevalence of sloppiness highlights the power of collective fits and suggests that modelers should focus on predictions rather than on parameters.

Variational method for estimating the rate of convergence of Markov-chain Monte Carlo algorithms.

We demonstrate the use of a variational method to determine a quantitative lower bound on the rate of convergence of Markov chain Monte Carlo (MCMC) algorithms as a function of the target density and proposal density. The bound relies on approximating the second largest eigenvalue in the spectrum of the MCMC operator using a variational principle and the approach is applicable to problems with continuous state spaces. We apply the method to one dimensional examples with Gaussian and quartic target densities, and we contrast the performance of the random walk Metropolis-Hastings algorithm with a "smart" variant that incorporates gradient information into the trial moves, a generalization of the Metropolis adjusted Langevin algorithm. We find that the variational method agrees quite closely with numerical simulations. We also see that the smart MCMC algorithm often fails to converge geometrically in the tails of the target density except in the simplest case we examine, and even then care must be taken to choose the appropriate scaling of the deterministic and random parts of the proposed moves. Again, this calls into question the utility of smart MCMC in more complex problems. Finally, we apply the same method to approximate the rate of convergence in multidimensional Gaussian problems with and without importance sampling. There we demonstrate the necessity of importance sampling for target densities which depend on variables with a wide range of scales.

Extracting falsifiable predictions from sloppy models.

Successful predictions are among the most compelling validations of any model. Extracting falsifiable predictions from nonlinear multiparameter models is complicated by the fact that such models are commonly sloppy, possessing sensitivities to different parameter combinations that range over many decades. Here we discuss how sloppiness affects the sorts of data that best constrain model predictions, makes linear uncertainty approximations dangerous, and introduces computational difficulties in Monte-Carlo uncertainty analysis. We also present a useful test problem and suggest refinements to the standards by which models are communicated.

Archaic inheritance: supporting high-altitude life in Tibet.

The Tibetan Plateau, often called the roof of the world, sits at an average altitude exceeding 4,500 m. Because of its extreme altitude, the Plateau is one of the harshest human-inhabited environments in the world. This, however, did not impede human colonization, and the Tibetan people have made the Tibetan Plateau their home for many generations. Many studies have quantified their markedly different physiological response to altitude and proposed that Tibetans were genetically adapted. Recently, advances in sequencing technologies led to the discovery of a set of candidate genes which harbor mutations that are likely beneficial at high altitudes in Tibetans. Since then, other studies have further characterized this impressive adaptation. Here, in this minireview, we discuss the progress made since the discovery of the genes involved in Tibetans' adaptation to high altitude with a particular emphasis on describing the series of studies that led us to conclude that archaic human DNA likely contributed to this impressive adaptation.

Web server to identify similarity of amino acid motifs to compounds (SAAMCO).

Protein-protein interactions are fundamental in mediating biological processes including metabolism, cell growth, and signaling. To be able to selectively inhibit or induce protein activity or complex formation is a key feature in controlling disease. For those situations in which protein-protein interactions derive substantial affinity from short linear peptide sequences, or motifs, we can develop search algorithms for peptidomimetic compounds that resemble the short peptide's structure but are not compromised by poor pharmacological properties. SAAMCO is a Web service ( http://bioware.ucd.ie/ approximately saamco) that facilitates the screening of motifs with known structures against bioactive compound databases. It is built on an algorithm that defines compound similarity based on the presence of appropriate amino acid side chain fragments and a favorable Root Mean Squared Deviation (RMSD) between compound and motif structure. The methodology is efficient as the available compound databases are preprocessed and fast regular expression searches filter potential matches before time-intensive 3D superposition is performed. The required input information is minimal, and the compound databases have been selected to maximize the availability of information on biological activity. "Hits" are accompanied with a visualization window and links to source database entries. Motif matching can be defined on partial or full similarity which will increase or reduce respectively the number of potential mimetic compounds. The Web server provides the functionality for rapid screening of known or putative interaction motifs against prepared compound libraries using a novel search algorithm. The tabulated results can be analyzed by linking to appropriate databases and by visualization.

Sloppy-model universality class and the Vandermonde matrix.

In a variety of contexts, physicists study complex, nonlinear models with many unknown or tunable parameters to explain experimental data. We explain why such systems so often are sloppy: the system behavior depends only on a few "stiff" combinations of the parameters and is unchanged as other "sloppy" parameter combinations vary by orders of magnitude. We observe that the eigenvalue spectra for the sensitivity of sloppy models have a striking, characteristic form with a density of logarithms of eigenvalues which is roughly constant over a large range. We suggest that the common features of sloppy models indicate that they may belong to a common universality class. In particular, we motivate focusing on a Vandermonde ensemble of multiparameter nonlinear models and show in one limit that they exhibit the universal features of sloppy models.