The evolutionary dynamics of metabolic protocells
[['The RNA world hypothesis [1–3] describes a stage in the early evolution of life where RNA constituted both the catalytic and the informational component of living systems .', 'This scenario has a reasonable degree of plausibility accounted for by the evidence gathered from both prebiotic chemistry [4 , 5] and the analytic approach of evolutionary biology [6] to the study of the origins of life .', 'Constrained by this assumption , elucidating the evolutionary dynamics of primordial replicators becomes a fundamental question , initially addressed by the quasi-species model [7 , 8] .', 'This formalism assumed that error rates for early self-replicators must have been inherently high leading to a distinct evolutionary behavior , namely that a sequence , regardless of its fitness , never fixates in the population , but instead is always surrounded by its mutant neighborhood forming a quasispecies .', 'Furthermore , given a sharply peaked fitness landscape , if the mutation rate is increased to a critical value , known as the error threshold , the extinction of the population known as the error catastrophe will occur .', 'One possible solution would be to reduce the mutation rate , but that would almost certainly involve a more complex self-replicator , encoded by a sequence of considerable length , which is in turn limited by the error rate .', 'This evolutionary cul-de-sac , known as Eigen’s paradox , required a reassessment of the model of protobiotic systems in light of their poor evolutionary potential .', 'Trying to address this issue led to the consideration of distributing the information encoded by the sequences in a mutualistic community of replicators known as the hypercycle model [9] in order to relax the error rate .', 'However , the large amplitude of the oscillations of the stationary state of the hypercycle model could lead to its extinction upon stochastic perturbation which , together with the susceptibility of the hypercycle model to the invasion of selfish molecules , make it difficult for the community of replicators to emerge without a form of population viscosity [10] .', 'Viscosity could have been provided by limited spatial dispersal over a substrate like clay [11 , 12] or as a limiting case , the compartmentalization of the sequences in a protocell acquiring a multi-level selection [13 , 14] structure when a the vesicle survival is coupled to its sequence content .', 'The hypothetical protocell organization has been extensively modelled in the past decades , from population genetics multilevel selection models to explicit biophysical accounts of chemical reactions that contribute to the formation of the vesicle membrane [15 , 16] .', 'These two avenues have remained relatively separated in the literature and the population genetic models have not fully explored the consequences of deriving the fitness function from realistic models of catalysed reaction networks that could have occurred inside the vesicle .', 'On the other hand , the biophysically inspired models are mainly concerned with the coupling of chemical reactions to the division of the membrane but do not incorporate population genetic considerations such as mutation , recombination or variations in heredity mechanisms .', 'The population genetics multilevel selection protocell models have their origins in the package model [17 , 18] and the stochastic error corrector model [19 , 20] .', 'These models describe the evolutionary dynamics of a set of replicating sequences encapsulated in a protocell and consider the processes that such a system would have had to endure at the dawn of life such as random assortment of the sequences upon fission , mutations to non-functional sequence types and conflict between selection at the sequence level and the vesicle level .', 'All of these processes suppose that the protocell population exhibits an average fitness lower than the maximum achievable which is commonly referred as genetic load .', 'In these models , the protocell level fitness is commonly specified by a simple function , like the geometric average , of the frequencies of the sequence types contained in the vesicle .', 'The dynamics at the sequence level inside the protocell is described by a replication-mutation process .', 'Deleterious mutation is typically introduced into the system in the form of irreversible mutation to a defective type that does not directly affect the protocell fitness except by reducing the effective total number of functional sequences .', 'These protocell models account , in a non-systematic manner , for the study of additional complexities and processes to the protocell dynamics .', 'For example , the incorporation of mutation schemes exhibiting trade-offs in catalytic activities [21 , 22] , modifications in the mechanism of inheritance or the incorporation of protocell fusion [23] and recombination [24] .', 'Few of them have directly considered a model of the structure of the underlying metabolic network and those that have , did so in a limited manner by studying a linear chain of unsaturated Michaelis-Menten enzymes [25] .', 'Most of the vesicle models depict a scenario where the nature of the encapsulated sequences is RNA , corresponding to the idea of the ribocell as a solid candidate for the protocellular stage in the early evolution of life [26] .', 'The coupling between metabolism and replication would have been particularly strong in the hypothesized ribocells as the RNA sequences are supposed to have operated as both the catalysts and the carriers of the genetic information .', 'We conjecture that this must have limited the space of possible catalysed chemical reaction networks which could be encapsulated within the protocells .', 'In this work , we incorporate a simple implementation of enzymatic metabolism to the protocell model to address the problem of determining the complexity and architectural characteristics of early metabolic networks .', 'The approach of considering the organism structure has been taken in other studies of the evolution of biological networks [27 , 28] that highlight the relevance of accounting for some level of biological detail to study evolutionary dynamics .', 'Metabolic networks are modelled with different degrees of dynamic and kinetic detail .', 'With a high level of kinetic detail , metabolic networks can be modelled with systems of differential equations incorporating aspects of the reaction mechanisms and properties of their catalysts .', 'These same reactions can be studied in the stochastic limit when the number of particles considered is low .', 'A framework that connects both the stochastic and the deterministic limits of a reaction network is that of chemical reaction network theory [29] .', 'With a lower level of kinetic detail , one can use modelling techniques such as flux balance analysis to understand the general behavior of reaction fluxes at stationary states [30] .', 'As a compromise between realism and computability we modeled the protocellular metabolism as a mass action irreversible reaction network with a small number of catalysts .', 'In our model , each of the metabolic reactions is catalyzed by a sequence of a given type with a Michaelis-Menten kinetic mechanism and thus exhibiting saturation behavior .', 'The output of this metabolic network determines the fitness of the protocell .', 'We are interested in the transition to the coexistence of all sequence types at the protocell population level and how this depends on the kinetics and architecture of the network .', 'For simple protocell models , the transition to coexistence has been extensively studied analytically [31] .', 'Sequence coexistence provides a way to increase the complexity of protocellular metabolism by harboring different types of catalysts .', 'We find that saturated metabolic fluxes provide a way to increase the coexistence of sequence types in the protocell by allowing the necessary redundancy to alleviate the fitness load produced by assortment when the protocell divides .', 'We also find that the combination of parallel fluxes with a terminal bimolecular catalyst provides a way of increasing the complexity of the protocellular metabolism .'], ['We model the protocell evolutionary dynamics as a Wright-Fisher process following the approach in [32] .', 'We augment the package model with a fitness function given by the flux of a critical metabolite , ϕ ( M ) , coming from a small metabolic network comprised of a set of chemical reactions , r , catalysed by the sequences inside the protocell .', 'Each protocell contains a number of sequences from a finite set of sequence types , or genes , g , of cardinality |g| = |r| + 1 , meaning that there is a sequence type catalysing each reaction of the metabolic network plus a mutational sink , the type ω , that does not contribute to the metabolic activity .', 'The population dynamics is derived from the following protocell life-cycle .', 'At the beginning of the life-cycle , N sequences are found in the vesicle distributed amongst the set of sequence types and each sequence is equally likely to replicate .', 'The total number of sequences at the beginning of the protocell life-cycle will be herein referred as the protocell ploidy .', 'When the sequences replicate they can mutate to the defective type , ω , with a probability traditionally unpacked in terms of the sequence length , L , and the mutation probability per base , μ , so that the accuracy of the replication process is given by q = ( 1 − μ ) L .', 'Hence , the sequences have a probability ( 1 − q ) to mutate to the defective type .', 'After reaching twice the ploidy , 2N , the vesicle splits in half and the sequences randomly assort between the offspring protocells regaining the starting total number of sequences N . At the population level , once all the protocells have divided , half of them are selected with a survival probability proportional to the critical metabolite flux .', 'A simplified treatment of metabolic networks is used to define the protocell fitness function in our simulations .', 'We evaluated the consequences on the protocell evolutionary dynamics of varying the kinetic parameters and the organization of the catalysts in different network topologies .', 'Our model assumes the separation of the dynamics at three time scales .', 'The quickest time scale is that of metabolic reactions catalysed by the sequences .', 'The next time scale is that of the replication of sequences and division of the protocell .', 'The third and longest time scale is that of the dynamics of the population of protocells .', 'We make use of simplifying assumptions based on each of the time scales being much longer than the previous time scale .', 'We illustrate the effects of the variations in the kinetic and stoichiometric parameters of the network with the example of a protocell containing two catalysts ( Fig 1A ) .', 'This network is the starting point for the construction of larger networks which allow us study the effects of the network architecture on the protocell evolutionary dynamics .', 'The first catalyst , the keystone ( denoted as α ) , takes an environmental input A that diffuses inside the protocell to produce the critical metabolite , M , that is required for the protocell survival .', 'The keystone catalyst also produces certain amount , c w α , of a waste product , W that can be converted by a second catalyst , the recycler ( denoted as β ) , into c M β molecules of critical metabolite .', 'A protocell composed solely of keystone catalysts is viable and thus constitutes an absorbing boundary of the stochastic process .', 'For all of the networks studied , we assume that the environmental input is present at sufficiently large quantities such that its concentration is not affected by the protocell population growth and the keystone catalysts that use it are in the saturation regime .', 'The corresponding metabolic flux for this scenario according to our formalism is given by ( see Methods , section 4 . 1 . 2 ) :', 'Φ = k cat α [ α ] + c M β min ( c w α k cat α [ α ] , k cat β [ β ] ) ( 1 ) Here kcat denotes the catalytic constant of the catalysts .', 'With the total ploidy constraint , the fitness function defined as the critical metabolite flux of this metabolic network ( Fig 1B ) will have a single maximum with a position on the space of protocell types given by the kinetic parameters and the protocell composition , more concretely when c w α k cat α [ α ] = k cat β [ β ] .', 'The position of the maximum determines the behaviour of the observables of the population like the coexistence fraction; the fraction of the population with protocells containing the maximum number of functional sequence types .', 'This dependence can be shown by choosing the kinetic and stoichiometric parameters so that the flux is closer or further away from the absorbing boundary of protocells with no recycler sequences .', 'One simple way of achieving this is by increasing the catalytic constant , k cat β , of the recycler sequence to get closer to the boundary or increasing the waste output of the keystone sequence by changing the stoichiometry , c w α , to get further away from it .', 'The remaining stoichiometric coefficients can be adjusted to achieve the same maximum flux for comparison as the fitness difference between the maximum in the interior and the fitness of the absorbing boundary of the stochastic process ( that is when the protocell lacks any recycler sequences and just harbors the keystone ) also affects the ploidy at which the coexistence transition occurs .', 'The coexistence fraction of the stationary population for the flux with the maximum further away from the absorbing boundary is always higher than the one closer to it ( Fig 1C ) .', 'This is due to the effects of stochastic assortment when protocells divide that can bring a fraction of the population to fixate at the boundary .', 'This probability is smaller as we get further away from the boundary and it vanishes upon increasing the ploidy of the protocell leading to a stationary population of coexisting types only .', 'Hence , fluxes with a maximum further away from the absorbing boundary benefit coexistence as they are less affected by the loss introduced by the assortment .', 'In terms of the average fitness of the stationary population , the flux further away from the boundary also presents a higher average fitness for all ploidy values and consequently if the two fluxes were to be competed with no interaction between the protocells the flux further from the boundary would dominate at the infinite population size limit ( see S1 Fig ) .', 'This is because the flux that is further away from the boundary is also flatter than the other flux and thus suffers from less assortment load as the population can spread over a larger region of viable coexisting types .', 'When mutations leading to the defective type are taken into account , the decay of the coexistence fraction of the flux closer to the boundary is faster when increasing the mutation rate than that of the flux further away from it ( Fig 1D ) .', 'In terms of average fitness of the populations , both landscapes show similar values at low mutation rates but the flux further to the boundary becomes fitter as the mutation rate increases and the coexisting types are lost in the flux closer to the boundary .', 'As noted by [32] , the way mutation towards a defective sequence type is usually included in the protocell model has no other effect at the protocell level than reducing the total number of functional types in the vesicle .', 'Therefore , the impact of mutation on the coexistence fraction of the stationary population will be equivalent to that of reducing the ploidy to the effective ploidy given the mutation rate and the starting ploidy of the protocell .', 'For a given flux function and mutation rate , the effective ploidy will saturate at a constant value as the total starting ploidy is increased .', 'At high mutation rates the effective ploidy becomes sensitive to the assortment in the flux closer to the boundary and thus the flux further away from the boundary is also more robust to the effects of introducing non-functional mutants .', 'We now focus on analysing effects of network architecture by considering two different arrangements of three catalytic sequences ( Fig 2A ) .', 'We will add the third recycler catalyst in serial or parallel with respect to the second recycler sequence of the previous example .', 'In both cases studied , we assume for simplicity that the keystone catalyst produces two different waste products to have non-competing reactions .', 'In order to discern effects derived from the topology of the networks alone , the total fitness benefit a protocell containing all three catalysts and the marginal fitness benefit of the third sequence with respect to the other two given by the kinetic and stoichiometric parameters was set to be the same for each of the studied architectures .', 'According to our simplifying treatment ( see Supplement , section 7 . 1 ) , the output metabolite fluxes for these two architectures are:', 'Φ S e r i a l = k cat α [ α ] + c M β min ( k cat α [ α ] , k cat β [ β ] ) + c M γ min ( k cat α [ α ] , k cat β [ β ] , k cat γ [ γ ] ) ( 2 ) Φ P a r a l l e l = k cat α [ α ] + c M β min ( k cat α [ α ] , k cat β [ β ] ) + c M γ min ( k cat α [ α ] , k cat γ [ γ ] ) ( 3 ) The fluxes over the simplex of protocell types all have their maxima at the center for the parameters chosen removing the effect of the proximity to the absorbing boundaries .', 'The transition to coexistence for the parallel addition of the third catalyst occurs at a lower ploidy values than that of the serial addition ( Fig 2B ) .', 'The coexistence fraction also grows faster with ploidy for the parallel architecture .', 'Both observations can be explained in terms of the assortment load , as the parallel architecture has a larger region of viable protocell types surrounding the maximum making it more robust to assortment by allowing the population to spread over this region ( Fig 2C ) .', 'The serial architecture , on the other hand , has a smaller region of viable types surrounding the maximum and becomes more sensitive to the assortment .', 'The mean fitness of the population reflects also this pattern given by the assortment load ( S2 Fig ) .', 'It can be concluded that for our model the population observables are mostly determined by the features of the fitness landscape given by the metabolic flux .', 'The protocell population can be seen as a quasispecies with a cloud of vesicles around the fittest type spread by the assortment at a particular ploidy .', 'Adding mutation to the system does not change the conclusion of the parallel addition higher robustness as mutations factor in the dynamics by effectively reducing the ploidy of the system .', 'The process of increasing the complexity of the studied networks by sequentially connecting to them additional recycler catalysts given the considered motifs , raises the question about what is the best way to perform this addition in order to guarantee a transition to coexistence of the resulting networks at low ploidy values .', 'Adding a new sequence type to the protocell always comes at the cost of increasing the effects of the assortment , which depends on ratio the protocell ploidy and the total number of sequence types .', 'This means that any increase in the coding capacity of the protocell , in terms of the number of coexisting sequence types , requires in turn an increase in the vesicle ploidy .', 'It has been shown that by making the ploidy sufficiently large , the numbers of types that can coexist in the protocell can be made to increase arbitrarily [32 , 33] .', 'This is , however , just for the cases with low competition at the sequence level and no mutations .', 'The requirement of increasing the ploidy becomes a burden once mutation is considered , as higher ploidy protocells accumulate mutants at a faster rate ( S3 Fig ) .', 'If certain amount of mutants are considered deleterious as in the early versions of the package model , the mutation process will set a limit to the ploidy increase .', 'The mutation process also establishes a limit on the number of catalytically active sequences in the protocell by saturating to a constant effective ploidy regardless of the starting ploidy value .', 'This can bring the population to the absorbing boundaries if the effective ploidy is low enough for the assortment to dominate .', 'A way of keeping the ploidy moderate while adding new catalysts to the network is by increasing the marginal fitness benefit of each added sequence ( changing the stoichiometric coefficient cM ) .', 'This can substantially lower the critical ploidy for the coexistence transition and its effects will depend on both the fitness benefit and the topology of the previous architecture in place .', 'For instance , if we were to add another catalyst with a high marginal fitness benefit to our two previous examples , the parallel architecture would become unbalanced leading to the fixation of one of the branches while the serial addition would still be viable .', 'In order for the parallel architecture to accept the new sequence with a higher marginal fitness benefit without becoming unbalanced a bimolecular catalyst should be considered .', 'A bimolecular catalyst can be added to two independent keystone catalysts and allow them to coexist .', 'We consider the case of two keystone catalysts converting two independent environment metabolites into the critical metabolite and also producing two different waste products that are in turn required for the third recycler catalyst in a bimolecular reaction to produce more metabolite of interest ( Supplement , section 7 . 2 .', 'The parameters were chosen to have the same marginal fitness benefit and maximum fitness than the previously studied cases .', 'According to our simplifying the flux would be:', 'Φ B i m o l e c u l a r = c M α k cat α [ α ] + c M β k cat β [ β ] + c M γ min ( k cat α [ α ] , k cat β [ β ] , k cat γ [ γ ] ) ( 4 ) The bimolecular motif flux exhibits a large region of viable types around the maximum however , the protocell population is also able to spread to an absorbing boundary with a large neutral space as the two keystone catalyst have the same marginal fitness benefit ( Fig 3A ) .', 'These features of the fitness landscape contribute to a large assortment load and a substantial fraction of the population is close to the absorbing boundary even at relatively high ploidy values ( Fig 3C ) .', 'Therefore , the coexistence transition for this architecture occurs at larger ploidy values than the other two architectures in spite of a larger region of viable types surrounding the flux maximum ( Fig 3B ) .', 'On the other hand , the third bimolecular catalyst is able to provide a way for the coexistence of these two keystone catalysts .', 'Moreover , the bimolecular reaction gives us a way of adding a catalyst with a large fitness benefit to the parallel architecture without unbalancing it .', 'We now consider the fluxes derived from the addition of a fourth catalytic sequence to the previous networks as a unimolecular or bimolecular recycler ( Supplement , section 7 . 3 ) .', 'We again consider the production of an additional waste metabolite to have non-competing reactions .', 'When the added recycler catalyst has the same marginal fitness benefit as the other two recyclers , the parallel architecture shows the transition to coexistence at the lowest ploidy value , as in the three catalyst networks previously analysed ( Fig 4A ) .', 'The parallel architecture is followed by the two other architectures with parallel motifs showing that for the case of equal marginal fitness benefit the parallel architectures are better at supporting the coexistence of catalyst and thus maintaining metabolic complexity .', 'This is due to the larger region of viable types surrounding the flux maximum for the parallel architectures .', 'When the marginal fitness benefit of the fourth recycler catalyst is greater than the other two , the parallel architectures are no longer the ones showing coexistence transitions at lower ploidy values .', 'The bimolecular and the serial architectures will be the ones showing a markedly earlier transition to coexistence than the architectures with parallel motifs ( Fig 4B ) .', 'The addition of a recycler with larger marginal fitness benefit unbalances the parallel architectures unless it is added as a bimolecular motif .', 'Increasing the complexity of metabolic networks in the protocell model requires the balance between the limitations given assortment load , the marginal fitness benefit and the risk of unbalancing the previously existing architectures .', 'We find that combining bimolecular and parallel pathways one can achieve networks with transitions to coexistence at lower ploidy values that would avert the problems introduced by mutations at large ploidy values .', 'Our finding is interesting in the light of the prevalence of bow-tie modules in extant metabolic networks [34 , 35] as the concatenation of the parallel and bimolecular motifs could produce such a network feature .'], ['Mathematical models of early life forms constitute a powerful tool to test the feasibility of proposed hypothetical scenarios in the first stages of life .', 'They allow us to increase our understanding of the conjectured scenarios and provide intuition about how to overcome their limitations .', 'Our metabolic model of the protocell serves to analyse the effects of considering details of the kinetic mechanism of the catalysts encapsulated and their network architecture .', 'Our model provides an example of how the fitness load introduced by protocell fission assortment process and its effects on the coexistence of catalysts can be alleviated by modifications in their kinetic parameters and their network organization .', 'The requirement to withstand the assortment load limits at low ploidy values the metabolic architectures that can be effectively encoded in the protocell .', 'Although the networks analyzed here are simple and small , they already show interesting consequences of accounting for the saturation of the catalysts in different architectures .', 'We demonstrate that variations in the kinetic parameters that keep the catalyst saturated provide a way to increase the catalyst redundancy and diminish the assortment effects .', 'The saturation of the catalysts means that the reaction they control becomes limiting and additional catalytic sequences will consequently increase the output flux .', 'The other major limitation for the evolution of a complex protocellular metabolism is the interplay of the effects of the assortment and the network architecture .', 'We demonstrate that architectures with a larger region of viable types surrounding the maximum in the interior favour the coexistence transition as they allow for the quasispecies to spread .', 'The architectures that are favoured also depend on how they are assembled , more specifically , the marginal fitness benefit of the new catalyst with respect to the previous network catalysts .', 'This benefit can happen directly by a higher stoichiometric benefit or indirectly by increasing the catalytic efficiency of the enzyme which allows the limiting catalyst to increase their numbers in the vesicle virtue of the flux saturation .', 'For architectures containing parallel motifs it will be detrimental to accept the new catalyst if it has a higher marginal fitness benefit unless it is involved in a bimolecular reaction .', 'The bimolecular motif prevents the system from becoming unbalanced and becomes a fundamental building block for larger architectures .', 'Nevertheless , the assortment load together with the coupling of genes and catalysts , limit the architectures of the metabolic networks that can be maintained in the vesicle by coexistence of its catalysts .', 'This limit in the metabolic evolutionary potential of the protocell can be alleviated , as it has been previously pointed out [36] , by the evolution of the decoupling between genetic information and catalytic activity , either with the advent of proteins or DNA in a ribozyme world or by the divergence of complementary RNA sequences [37] .', 'The evolution of correlated inheritance like the chromosomal organization would be another mechanism that would allow the exploration of a larger space of network architectures and catalytic efficiency [38] as it would diminish the assortment load .', 'An example of releasing this constraint can be found in [25] which makes use of the unsaturated linear metabolic chain formulation of [39] to address the effects of chromosomal organization in the evolution toward specialization of promiscuous catalysts .', 'The model presented incorporates assumptions about the biochemistry of the hypothetical protocell and standard population dynamics that can be sophisticated or relaxed with interesting consequences .', 'For example , the incorporation of trade-off between catalytic and replication activity has shown interesting dynamical features at the population level such as sustained oscillations [40] .', 'The metabolism model could be further enriched in the direction of the whole cell model [41] by increasing the biochemical details considered .', 'The only limitation is that the nature of early living things is purely conjectural and it would be preferable to understand general features of the evolution of metabolic networks rather than a concrete instantiation mirroring extant life .', 'Together with this consideration , a better model for the development stage of the protocell and the corresponding population dynamics would allow us to decompose the summary statistic of the evolutionary dynamics , that is fitness , in terms of the whole interplay of the network flux and the replication of catalytic sequences throughout the developmental stage .', 'A more detailed account on protocell development would require a full stochastic treatment that properly accounts for the separation of metabolic and replicative dynamic time scales .', 'The effects of introducing other processes in the population dynamics like alternative division mechanisms in the metabolic protocell model as described by [42 , 43] or modification in the heredity mechanisms already studied in previous protocell models could also yield interesting results .', 'Further incorporating considerations about environment would lead us in the direction of understanding the ecology of the protocells .', 'If our intermediates that leak out accumulate significantly in the environment [44] , interesting frequency dependent relations can be established in the protocell populations and one could add protocell fusion to study the assembly of networks in this augmented model .', 'The protocell organization is considered a solid candidate for a transitional form of life and , as such , has been extensively explored computationally and is starting to be addressed in experimental settings [45] .', 'The utility of the protocell theoretical models has been recently validated in helping the development of chemical implementations of protocell structures [46] .', 'This also provides the experimental set-up to test the conclusions of the theoretical works in the protocell literature .', 'Furthermore , construction of the first synthetic autonomous protocell is becoming a reality and we believe that theoretical models can be useful to further this endeavour .'], ['We will model the metabolic activity of our protocell as a well-stirred reactor in which ribozymes , environmental inputs , and metabolites interact .', 'More specifically , the metabolic dynamics will be modelled as Michaelis-Menten enzyme dynamics where each flux is given by the expression', 'v = k cat [ E ] [ S ] K m + [ S ] ( 5 )', 'Where , [E] and [S] denote the concentrations of enzyme and substrate respectively and kcat is the turnover number ( catalytic constant ) of the enzyme and Km its Michaelis constant .', 'We will also study bimolecular reactions with mechanisms of the following form:', '( 6 )', 'The velocity will have the functional form', 'v r = k cat [ S 1 ] [ S 2 ] K m + a 1 c 1 + a 2 c 2 + c 1 c 2 ( 7 )', 'where the values of a1 and a2 can be determined from either the forward and the backward rates of the underlying reactions or from thermodynamic potentials .', 'Special values of these constants correspond to the cases of random order , definite order
['Introduction', 'Results', 'Discussion', 'Methods']
['Protocell multilevel selection models have been proposed to study the evolutionary dynamics of vesicles encapsulating a set of replicating , competing and mutating sequences .', 'The frequency of the different sequence types determines protocell survival through a fitness function .', 'One of the defining features of these models is the genetic load generated when the protocell divides and its sequences are assorted between the offspring vesicles .', 'However , these stochastic assortment effects disappear when the redundancy of each sequence type is sufficiently high .', 'The fitness dependence of the vesicle with its sequence content is usually defined without considering a realistic account on how the lower level dynamics would specify the protocell fitness .', 'Here , we present a protocell model with a fitness function determined by the output flux of a simple metabolic network with the aim of understanding how the evolution of both kinetic and topological features of metabolism would have been constrained by the particularities of the protocell evolutionary dynamics .', 'In our model , the sequences inside the vesicle are both the carriers of information and Michaelis-Menten catalysts exhibiting saturation .', 'We found that the saturation of the catalysts controlling the metabolic fluxes , achievable by modifying the kinetic or stoichiometric parameters , provides a mechanism to ameliorate the assortment load by increasing the redundancy of the catalytic sequences required to achieve the maximum flux .', 'Regarding the network architecture , we conclude that combinations of parallel network motifs and bimolecular catalysts are a robust way to increase the complexity of the metabolism enclosed by the protocell .']
['The protocell hypothesis conjectures the existence of a vesicle containing catalytic and replicating sequences as the primordial cellular organization during the early stages of the evolution of life .', 'Mathematical models of protocells traditionally consider RNA sequences being encapsulated and having both an informational and a catalytic role in the same molecule .', 'Because of this dual function , the protocell sequences are evolutionarily constrained .', 'Mathematical models have been extensively used to study the evolutionary dynamics of protocells with a focus on the processes , like mutation or stochastic sequence assortment upon division , that affect the protocell information capacity in terms of the coexistence of different sequence types .', 'Here we introduce a simple model of metabolic networks whose output determines the survival of the protocell with the aim of studying the effect of modifying the kinetic and architectural properties of the network on sequence coexistence .', 'We find that stochastic assortment and mutation limit the architectures able to be encapsulated by the protocell with a given fraction of the population harbouring all possible sequence types .']
['vesicles', 'population genetics', 'ploidy', 'metabolic networks', 'vertebrates', 'enzymology', 'animals', 'mammals', 'network analysis', 'enzyme metabolism', 'sequence motif analysis', 'population biology', 'cellular structures and organelles', 'enzyme chemistry', 'research and analysis methods', 'sequence analysis', 'computer and information sciences', 'catalysts', 'network motifs', 'bioinformatics', 'chemistry', 'biochemistry', 'eukaryota', 'cell biology', 'cats', 'database and informatics methods', 'catalysis', 'genetics', 'biology and life sciences', 'physical sciences', 'evolutionary biology', 'amniotes', 'organisms']