In this study, we present a novel methodology to infer indel parameters from multiple sequence alignments (MSAs) based on simulations. Our algorithm searches for the set of evolutionary parameters describing indel dynamics which best fits a given input MSA. In each step of the search, we use parametric bootstraps and the Mahalanobis distance to estimate how well a proposed set of parameters fits input data. Using simulations, we demonstrate that our methodology can accurately infer the indel parameters for a large variety of plausible settings. Moreover, using our methodology, we show that indel parameters substantially vary between three genomic data sets: Mammals, bacteria, and retroviruses. Finally, we demonstrate how our methodology can be used to simulate MSAs based on indel parameters inferred from real data sets.

Although parasitic organisms are found worldwide, the relative importance of host specificity and geographic isolation for parasite speciation has been explored in only a few systems. Here, we study Plasmodium parasites known to infect Asian nonhuman primates, a monophyletic group that includes the lineage leading to the human parasite Plasmodium vivax and several species used as laboratory models in malaria research. We analyze the available data together with new samples from three sympatric primate species from Borneo: The Bornean orangutan and the long-tailed and the pig-tailed macaques. We find several species of malaria parasites, including three putatively new species in this biodiversity hotspot. Among those newly discovered lineages, we report two sympatric parasites in orangutans. We find no differences in the sets of malaria species infecting each macaque species indicating that these species show no host specificity. Finally, phylogenetic analysis of these data suggests that the malaria parasites infecting Southeast Asian macaques and their relatives are speciating three to four times more rapidly than those with other mammalian hosts such as lemurs and African apes. We estimate that these events took place in approximately a 3–4-Ma period. Based on the genetic and phenotypic diversity of the macaque malarias, we hypothesize that the diversification of this group of parasites has been facilitated by the diversity, geographic distributions, and demographic histories of their primate hosts.

Mathematical models of infectious diseases are a valuable tool in understanding the mechanisms and patterns of disease transmission. It is, however, a difficult subject to teach, requiring both mathematical expertise and extensive subject-matter knowledge of a variety of disease systems. In this article, we explore several uses of zombie epidemics to make mathematical modeling and infectious disease epidemiology more accessible to public health professionals, students, and the general public. We further introduce a web-based simulation, White Zed (http://cartwrig.ht/apps/whitezed/), that can be deployed in classrooms to allow students to explore models before implementing them. In our experience, zombie epidemics are familiar, approachable, flexible, and an ideal way to introduce basic concepts of infectious disease epidemiology.

Network reconstruction is a fundamental problem for understanding many complex systems with unknown interaction structures. In many complex systems, there are indirect interactions between two individuals without immediate connection but with common neighbors. Despite recent advances in network reconstruction, we continue to lack an approach for reconstructing complex networks with indirect interactions. Here we introduce a two-step strategy to resolve the reconstruction problem, where in the first step, we recover both direct and indirect interactions by employing the Lasso to solve a sparse signal reconstruction problem, and in the second step, we use matrix transformation and optimization to distinguish between direct and indirect interactions. The network structure corresponding to direct interactions can be fully uncovered. We exploit the public goods game occurring on complex networks as a paradigm for characterizing indirect interactions and test our reconstruction approach. We find that high reconstruction accuracy can be achieved for both homogeneous and heterogeneous networks, and a number of empirical networks in spite of insufficient data measurement contaminated by noise. Although a general framework for reconstructing complex networks with arbitrary types of indirect interactions is yet lacking, our approach opens new routes to separate direct and indirect interactions in a representative complex system.

Time-resolved fluorescence spectroscopy was used to explore the pathway and kinetics of energy transfer in photosynthetic membrane vesicles (chromatophores) isolated from Rhodobacter (Rba.) sphaeroides cells harvested 2, 4, 6 or 24 hours after a transition from growth in high to low level illumination. As previously observed, this light intensity transition initiates the remodeling of the photosynthetic apparatus and an increase in the number of light harvesting 2 (LH2) complexes relative to light harvesting 1 (LH1) and reaction center (RC) complexes. It has generally been thought that the increase in LH2 complexes served the purpose of increasing the overall energy transmission to the RC. However, fluorescence lifetime measurements and analysis in terms of energy transfer within LH2 and between LH2 and LH1 indicate that, during the remodeling time period measured, only a portion of the additional LH2 generated are well connected to LH1 and the reaction center. The majority of the additional LH2 fluorescence decays with a lifetime comparable to that of free, unconnected LH2 complexes. The presence of large LH2-only domains has been observed by atomic force microscopy in Rba. sphaeroides chromatophores (Bahatyrova et al., Nature, 2004, 430, 1058), providing structural support for the existence of pools of partially connected LH2 complexes. These LH2-only domains represent the light-responsive antenna complement formed after a switch in growth conditions from high to low illumination, while the remaining LH2 complexes occupy membrane regions containing mixtures of LH2 and LH1–RC core complexes. The current study utilized a multi-parameter approach to explore the fluorescence spectroscopic properties related to the remodeling process, shedding light on the structure-function relationship of the photosynthetic assembles. Possible reasons for the accumulation of these largely disconnected LH2-only pools are discussed.

Our ability to uncover complex network structure and dynamics from data is fundamental to understanding and controlling collective dynamics in complex systems. Despite recent progress in this area, reconstructing networks with stochastic dynamical processes from limited time series remains to be an outstanding problem. Here we develop a framework based on compressed sensing to reconstruct complex networks on which stochastic spreading dynamics take place. We apply the methodology to a large number of model and real networks, finding that a full reconstruction of inhomogeneous interactions can be achieved from small amounts of polarized (binary) data, a virtue of compressed sensing. Further, we demonstrate that a hidden source that triggers the spreading process but is externally inaccessible can be ascertained and located with high confidence in the absence of direct routes of propagation from it. Our approach thus establishes a paradigm for tracing and controlling epidemic invasion and information diffusion in complex networked systems.

Two classes of scaling behaviours, namely the super-linear scaling of links or activities, and the sub-linear scaling of area, diversity, or time elapsed with respect to size have been found to prevail in the growth of complex networked systems. Despite some pioneering modelling approaches proposed for specific systems, whether there exists some general mechanisms that account for the origins of such scaling behaviours in different contexts, especially in socioeconomic systems, remains an open question. We address this problem by introducing a geometric network model without free parameter, finding that both super-linear and sub-linear scaling behaviours can be simultaneously reproduced and that the scaling exponents are exclusively determined by the dimension of the Euclidean space in which the network is embedded. We implement some realistic extensions to the basic model to offer more accurate predictions for cities of various scaling behaviours and the Zipf distribution reported in the literature and observed in our empirical studies. All of the empirical results can be precisely recovered by our model with analytical predictions of all major properties. By virtue of these general findings concerning scaling behaviour, our models with simple mechanisms gain new insights into the evolution and development of complex networked systems.

Background: Chemistry and particularly enzymology at surfaces is a topic of rapidly growing interest, both in terms of its role in biological systems and its application in biocatalysis. Existing protein immobilization approaches, including noncovalent or covalent attachments to solid supports, have difficulties in controlling protein orientation, reducing nonspecific absorption and preventing protein denaturation. New strategies for enzyme immobilization are needed that allow the precise control over orientation and position and thereby provide optimized activity.

Methodology/Principal Findings: A method is presented for utilizing peptide ligands to immobilize enzymes on surfaces with improved enzyme activity and stability. The appropriate peptide ligands have been rapidly selected from high-density arrays and when desirable, the peptide sequences were further optimized by single-point variant screening to enhance both the affinity and activity of the bound enzyme. For proof of concept, the peptides that bound to β-galactosidase and optimized its activity were covalently attached to surfaces for the purpose of capturing target enzymes. Compared to conventional methods, enzymes immobilized on peptide-modified surfaces exhibited higher specific activity and stability, as well as controlled protein orientation.

Conclusions/Significance: A simple method for immobilizing enzymes through specific interactions with peptides anchored on surfaces has been developed. This approach will be applicable to the immobilization of a wide variety of enzymes on surfaces with optimized orientation, location and performance, and provides a potential mechanism for the patterned self-assembly of multiple enzymes on surfaces.