@article{1077, abstract = {Viral capsids are structurally constrained by interactions among the amino acids (AAs) of their constituent proteins. Therefore, epistasis is expected to evolve among physically interacting sites and to influence the rates of substitution. To study the evolution of epistasis, we focused on the major structural protein of the fX174 phage family by first reconstructing the ancestral protein sequences of 18 species using a Bayesian statistical framework. The inferred ancestral reconstruction differed at eight AAs, for a total of 256 possible ancestral haplotypes. For each ancestral haplotype and the extant species, we estimated, in silico, the distribution of free energies and epistasis of the capsid structure. We found that free energy has not significantly increased but epistasis has. We decomposed epistasis up to fifth order and found that higher-order epistasis sometimes compensates pairwise interactions making the free energy seem additive. The dN/dS ratio is low, suggesting strong purifying selection, and that structure is under stabilizing selection. We synthesized phages carrying ancestral haplotypes of the coat protein gene and measured their fitness experimentally. Our findings indicate that stabilizing mutations can have higher fitness, and that fitness optima do not necessarily coincide with energy minima.}, author = {Fernandes Redondo, Rodrigo A and Vladar, Harold and Włodarski, Tomasz and Bollback, Jonathan P}, issn = {17425689}, journal = {Journal of the Royal Society Interface}, number = {126}, publisher = {Royal Society of London}, title = {{Evolutionary interplay between structure, energy and epistasis in the coat protein of the ϕX174 phage family}}, doi = {10.1098/rsif.2016.0139}, volume = {14}, year = {2017}, } @misc{9864, abstract = {Viral capsids are structurally constrained by interactions among the amino acids (AAs) of their constituent proteins. Therefore, epistasis is expected to evolve among physically interacting sites and to influence the rates of substitution. To study the evolution of epistasis, we focused on the major structural protein of the ϕX174 phage family by, first, reconstructing the ancestral protein sequences of 18 species using a Bayesian statistical framework. The inferred ancestral reconstruction differed at eight AAs, for a total of 256 possible ancestral haplotypes. For each ancestral haplotype and the extant species, we estimated, in silico, the distribution of free energies and epistasis of the capsid structure. We found that free energy has not significantly increased but epistasis has. We decomposed epistasis up to fifth order and found that higher-order epistasis sometimes compensates pairwise interactions making the free energy seem additive. The dN/dS ratio is low, suggesting strong purifying selection, and that structure is under stabilizing selection. We synthesized phages carrying ancestral haplotypes of the coat protein gene and measured their fitness experimentally. Our findings indicate that stabilizing mutations can have higher fitness, and that fitness optima do not necessarily coincide with energy minima.}, author = {Fernandes Redondo, Rodrigo A and de Vladar, Harold and Włodarski, Tomasz and Bollback, Jonathan P}, publisher = {The Royal Society}, title = {{Data from evolutionary interplay between structure, energy and epistasis in the coat protein of the ϕX174 phage family}}, doi = {10.6084/m9.figshare.4315652.v1}, year = {2016}, } @article{2917, abstract = {The search for extra-terrestrial intelligence (SETI) has been performed principally as a one-way survey, listening of radio frequencies across the Milky Way and other galaxies. However, scientists have engaged in an active messaging only rarely. This suggests the simple rationale that if other civilizations exist and take a similar approach to ours, namely listening but not broadcasting, the result is a silent universe. A simple game theoretical model, the prisoner's dilemma, explains this situation: each player (civilization) can passively search (defect), or actively search and broadcast (cooperate). In order to maximize the payoff (or, equivalently, minimize the risks) the best strategy is not to broadcast. In fact, the active search has been opposed on the basis that it might be dangerous to expose ourselves. However, most of these ideas have not been based on objective arguments, and ignore accounting of the possible gains and losses. Thus, the question stands: should we perform an active search? I develop a game-theoretical framework where civilizations can be of different types, and explicitly apply it to a situation where societies are either interested in establishing a two-way communication or belligerent and in urge to exploit ours. The framework gives a quantitative solution (a mixed-strategy), which is how frequent we should perform the active SETI. This frequency is roughly proportional to the inverse of the risk, and can be extremely small. However, given the immense amount of stars being scanned, it supports active SETI. The model is compared with simulations, and the possible actions are evaluated through the San Marino scale, measuring the risks of messaging.}, author = {Vladar, Harold}, journal = {International Journal of Astrobiology}, number = {1}, pages = {53 -- 62}, publisher = {Cambridge University Press}, title = {{The game of active search for extra terrestrial intelligence Breaking the Great Silence }}, doi = {10.1017/S1473550412000407}, volume = {12}, year = {2012}, } @article{3166, abstract = {There is evidence that the genetic code was established prior to the existence of proteins, when metabolism was powered by ribozymes. Also, early proto-organisms had to rely on simple anaerobic bioenergetic processes. In this work I propose that amino acid fermentation powered metabolism in the RNA world, and that this was facilitated by proto-adapters, the precursors of the tRNAs. Amino acids were used as carbon sources rather than as catalytic or structural elements. In modern bacteria, amino acid fermentation is known as the Stickland reaction. This pathway involves two amino acids: the first undergoes oxidative deamination, and the second acts as an electron acceptor through reductive deamination. This redox reaction results in two keto acids that are employed to synthesise ATP via substrate-level phosphorylation. The Stickland reaction is the basic bioenergetic pathway of some bacteria of the genus Clostridium. Two other facts support Stickland fermentation in the RNA world. First, several Stickland amino acid pairs are synthesised in abiotic amino acid synthesis. This suggests that amino acids that could be used as an energy substrate were freely available. Second, anticodons that have complementary sequences often correspond to amino acids that form Stickland pairs. The main hypothesis of this paper is that pairs of complementary proto-adapters were assigned to Stickland amino acids pairs. There are signatures of this hypothesis in the genetic code. Furthermore, it is argued that the proto-adapters formed double strands that brought amino acid pairs into proximity to facilitate their mutual redox reaction, structurally constraining the anticodon pairs that are assigned to these amino acid pairs. Significance tests which randomise the code are performed to study the extent of the variability of the energetic (ATP) yield. Random assignments can lead to a substantial yield of ATP and maintain enough variability, thus selection can act and refine the assignments into a proto-code that optimises the energetic yield. Monte Carlo simulations are performed to evaluate the establishment of these simple proto-codes, based on amino acid substitutions and codon swapping. In all cases, donor amino acids are assigned to anticodons composed of U+G, and have low redundancy (1-2 codons), whereas acceptor amino acids are assigned to the the remaining codons. These bioenergetic and structural constraints allow for a metabolic role for amino acids before their co-option as catalyst cofactors. Reviewers: this article was reviewed by Prof. William Martin, Prof. Eors Szathmary (nominated by Dr. Gaspar Jekely) and Dr. Adam Kun (nominated by Dr. Sandor Pongor)}, author = {Vladar, Harold}, journal = {Biology Direct}, publisher = {BioMed Central}, title = {{Amino acid fermentation at the origin of the genetic code}}, doi = {10.1186/1745-6150-7-6}, volume = {7}, year = {2012}, } @inbook{3277, abstract = {The problem of the origin of metazoa is becoming more urgent in the context of astrobiology. By now it is clear that clues to the understanding of this crucial transition in the evolution of life can arise in a fourth pathway besides the three possibilities in the quest for simplicity outlined by Bonner in his classical book. In other words, solar system exploration seems to be one way in the long-term to elucidate the simplicity of evolutionary development. We place these ideas in the context of different inheritance systems, namely the genotypic and phenotypic replicators with limited or unlimited heredity, and ask which of these can support multicellular development, and to which degree of complexity. However, the quest for evidence on the evolution of biotas from planets around other stars does not seem to be feasible with present technology with direct visualization of living organisms on exoplanets. But this may be attempted on the Galilean moons of Jupiter where there is a possibility of detecting reliable biomarkers in the next decade with the Europa Jupiter System Mission, in view of recent progress by landing micropenetrators on planetary, or satellite surfaces. Mars is a second possibility in the inner Solar System, in spite of the multiple difficulties faced by the fleet of past, present and future missions. We discuss a series of preliminary ideas for elucidating the origin of metazoan analogues with available instrumentation in potential payloads of feasible space missions to the Galilean moons.}, author = {de Vladar, Harold and Chela Flores, Julian}, booktitle = {Life on Earth and other planetary bodies}, pages = {387 -- 405}, publisher = {Springer}, title = {{Can the evolution of multicellularity be anticipated in the exploration of the solar system?}}, doi = {10.1007/978-94-007-4966-5_22}, volume = {24}, year = {2012}, } @article{3391, abstract = {Evolutionary biology shares many concepts with statistical physics: both deal with populations, whether of molecules or organisms, and both seek to simplify evolution in very many dimensions. Often, methodologies have undergone parallel and independent development, as with stochastic methods in population genetics. Here, we discuss aspects of population genetics that have embraced methods from physics: non-equilibrium statistical mechanics, travelling waves and Monte-Carlo methods, among others, have been used to study polygenic evolution, rates of adaptation and range expansions. These applications indicate that evolutionary biology can further benefit from interactions with other areas of statistical physics; for example, by following the distribution of paths taken by a population through time}, author = {de Vladar, Harold and Barton, Nicholas H}, journal = {Trends in Ecology and Evolution}, number = {8}, pages = {424 -- 432}, publisher = {Cell Press}, title = {{The contribution of statistical physics to evolutionary biology}}, doi = {10.1016/j.tree.2011.04.002}, volume = {26}, year = {2011}, } @article{3375, abstract = {By exploiting an analogy between population genetics and statistical mechanics, we study the evolution of a polygenic trait under stabilizing selection, mutation and genetic drift. This requires us to track only four macroscopic variables, instead of the distribution of all the allele frequencies that influence the trait. These macroscopic variables are the expectations of: the trait mean and its square, the genetic variance, and of a measure of heterozygosity, and are derived from a generating function that is in turn derived by maximizing an entropy measure. These four macroscopics are enough to accurately describe the dynamics of the trait mean and of its genetic variance (and in principle of any other quantity). Unlike previous approaches that were based on an infinite series of moments or cumulants, which had to be truncated arbitrarily, our calculations provide a well-defined approximation procedure. We apply the framework to abrupt and gradual changes in the optimum, as well as to changes in the strength of stabilizing selection. Our approximations are surprisingly accurate, even for systems with as few as five loci. We find that when the effects of drift are included, the expected genetic variance is hardly altered by directional selection, even though it fluctuates in any particular instance. We also find hysteresis, showing that even after averaging over the microscopic variables, the macroscopic trajectories retain a memory of the underlying genetic states.}, author = {de Vladar, Harold and Barton, Nicholas H}, journal = {Journal of the Royal Society Interface}, number = {58}, pages = {720 -- 739}, publisher = {The Royal Society}, title = {{The statistical mechanics of a polygenic character under stabilizing selection mutation and drift}}, doi = {10.1098/rsif.2010.0438}, volume = {8}, year = {2011}, } @phdthesis{4232, author = {Harold Vladar}, publisher = {Faculty of mathematical and natural sciences, University of Groningen}, title = {{Stochasticity and Variability in the dynamics and genetics of populations}}, doi = {3811}, year = {2009}, } @article{4234, abstract = {We study a generalised model of population growth in which the state variable is population growth rate instead of population size. Stochastic parametric perturbations, modelling phenotypic variability, lead to a Langevin system with two sources of multiplicative noise. The stationary probability distributions have two characteristic power-law scales. Numerical simulations show that noise suppresses the explosion of the growth rate which occurs in the deterministic counterpart. Instead, in different parameter regimes populations will grow with "anomalous" stochastic rates and (i) stabilise at "random carrying capacities", or (ii) go extinct in random times. Using logistic fits to reconstruct the simulated data, we find that even highly significant estimations do not recover or reflect information about the deterministic part of the process. Therefore, the logistic interpretation is not biologically meaningful. These results have implications for distinct model-aided calculations in biological situations because these kinds of estimations could lead to spurious conclusions. (c) 2006 Elsevier B.V. All rights reserved.}, author = {de Vladar, Harold and Pen, I.}, journal = {Physica A}, pages = {477 -- 485}, publisher = {Elsevier}, title = {{Determinism, noise, and spurious estimations in a generalised model of population growth}}, doi = {10.1016/j.physa.2006.06.025}, volume = {373}, year = {2007}, } @inproceedings{4233, author = {Harold Vladar}, editor = {Falcón,N. and Loyo de Sardi,Y.}, pages = {91 -- 109}, publisher = {Consejo de Desarrollo Cientifico y Tecnologico}, title = {{Alternativas prebióticas para la síntesis de amino- ácidos y otras moléculas relacionadas}}, doi = {3808}, year = {2007}, } @article{4237, abstract = {The growth function of populations is central in biomathematics. The main dogma is the existence of density-dependence mechanisms, which can be modelled with distinct functional forms that depend on the size of the Population. One important class of regulatory functions is the theta-logistic, which generalizes the logistic equation. Using this model as a motivation, this paper introduces a simple dynamical reformulation that generalizes many growth functions. The reformulation consists of two equations, one for population size, and one for the growth rate. Furthermore, the model shows that although population is density-dependent, the dynamics of the growth rate does not depend either on population size, nor on the carrying capacity. Actually, the growth equation is uncoupled from the population size equation, and the model has only two parameters, a Malthusian parameter rho and a competition coefficient theta. Distinct sign combinations of these parameters reproduce not only the family of theta-logistics, but also the van Bertalanffy, Gompertz and Potential Growth equations, among other possibilities. It is also shown that, except for two critical points, there is a general size-scaling relation that includes those appearing in the most important allometric theories, including the recently proposed Metabolic Theory of Ecology. With this model, several issues of general interest are discussed such as the growth of animal population, extinctions, cell growth and allometry, and the effect of environment over a population. (c) 2005 Elsevier Ltd. All rights reserved.}, author = {de Vladar, Harold}, journal = {Journal of Theoretical Biology}, number = {2}, pages = {245 -- 256}, publisher = {Elsevier}, title = {{Density-dependence as a size-independent regulatory mechanism}}, doi = {3802}, volume = {238}, year = {2006}, } @article{4235, author = {Harold Vladar and González,J. A}, journal = {Journal of Theoretical Biology}, pages = {91 -- 109}, publisher = {Elsevier}, title = {{Dynamic response of cancer under the influence of immunological activity and therapy}}, year = {2006}, } @inbook{4239, author = {Harold Vladar and Cipriani, Roberto and Scharifker, Benjamin and Bubis, Jose}, booktitle = {Life in the Universe From the Miller Experiment to the Search for Life on Other Worlds}, editor = {Seckbach,J. and Chela-Flores,J. and Owen,T. and Raulin,F.}, pages = {83 -- 87}, publisher = {Springer}, title = {{A Mechanism for the Prebiotic Emergence of Proteins}}, doi = {3807}, volume = {7}, year = {2004}, } @article{4238, abstract = {The dynamical basis of tumoral growth has been controversial. Many models have been proposed to explain cancer development. The descriptions employ exponential, potential, logistic or Gompertzian growth laws. Some of these models are concerned with the interaction between cancer and the immunological, system. Among other properties, these models are concerned with the microscopic behavior of tumors and the emergence of cancer. We propose a modification of a previous model by Stepanova, which describes the specific immunological response against cancer. The modification consists of the substitution of a Gompertian law for the exponential rate used for tumoral growth. This modification is motivated by the numerous works confirming that Gompertz's equation correctly describes solid tumor growth. The modified model predicts that near zero, tumors always tend to grow. Immunological contraposition never suffices to induce a complete regression of the tumor. Instead, a stable microscopic equilibrium between cancer and immunological activity can be attained. In other words, our model predicts that the theory of immune surveillance is plausible. A macroscopic equilibrium in which the system develops cancer is also possible. In this case, immunological activity is depleted. This is consistent with the phenomena of cancer tolerance. Both equilibrium points can coexist or can exist without the other. In all cases the fixed point at zero tumor size is unstable. Since immunity cannot induce a complete tumor regression, a therapy is required. We include constant-dose therapies and show that they are insufficient. Final levels of immunocompetent cells and tumoral cells are finite, thus post-treatment regrowth of the tumor is certain. We also evaluate late-intensification therapies which are successful. They induce an asymptotic regression to zero tumor size. Immune response is also suppressed by the therapy, and thus plays a negligible role in the remission. We conclude that treatment evaluation should be successful without taking into account immunological effects. (C) 2003 Elsevier Ltd. All rights reserved.}, author = {de Vladar, Harold and González, J.}, journal = {Journal of Theoretical Biology}, number = {3}, pages = {335 -- 348}, publisher = {Elsevier}, title = {{Dynamic response of cancer under the influence of immunological activity and therapy}}, doi = {3801}, volume = {227}, year = {2004}, } @inbook{4230, author = {Harold Vladar and Cipriani, Roberto and Scharifker, Benjamin and Bubis, Jose}, booktitle = {Life in the Universe From the Miller Experiment to the Search for Life on Other Worlds}, editor = {Hanslmeier,A. and Kempe,S. and Seckbach,J.}, pages = {83 -- 87}, publisher = {Springer}, title = {{A mechanism for the prebiotic emergence of proteins}}, year = {2004}, } @phdthesis{4236, author = {de Vladar, Harold}, publisher = {Centro de estudios avazados, IVIC}, title = {{Métodos no lineales y sus aplicaciones en dinámicas aleatorias de poblaciones celulares}}, doi = {3810}, year = {2004}, }