Network ecology: Tie fitness in social context(s).

Social relations are embedded in material, cultural, and institutional settings that affect network dynamics and the resulting topologies. For example, romantic entanglements are subject to social and cultural norms, interfirm alliances are constrained by country-specific legislation, and adolescent friendships are conditioned by classroom settings and neighborhood effects. In short, social contexts shape social relations and the networks they give rise to. However, how and when they do so remain to be established. This paper presents network ecology as a general framework for identifying how the proximal environment shapes social networks by focusing interactions and social relations, and how these interactions and relations in turn shape the environment in which social networks form. Tie fitness is introduced as a metric that quantifies how well particular dyadic social relations would align with the setting. Using longitudinal networks collected on two cohorts each in 18 North American schools, i.e., 36 settings, we develop five generalizable observations about the time-varying fitness of adolescent friendship. Across all 252 analyzed networks, tie fitness predicted new tie formation, tie longevity, and tie survival. Dormant fit ties cluster in relational niches, thereby establishing a resource base for social identities competing for increased representation in the relational system.

Overlapping timescales obscure early warning signals of the second COVID-19 wave.

Early warning indicators based on critical slowing down have been suggested as a model-independent and low-cost tool to anticipate the (re)emergence of infectious diseases. We studied whether such indicators could reliably have anticipated the second COVID-19 wave in European countries. Contrary to theoretical predictions, we found that characteristic early warning indicators generally decreased rather than increased prior to the second wave. A model explains this unexpected finding as a result of transient dynamics and the multiple timescales of relaxation during a non-stationary epidemic. Particularly, if an epidemic that seems initially contained after a first wave does not fully settle to its new quasi-equilibrium prior to changing circumstances or conditions that force a second wave, then indicators will show a decreasing rather than an increasing trend as a result of the persistent transient trajectory of the first wave. Our simulations show that this lack of timescale separation was to be expected during the second European epidemic wave of COVID-19. Overall, our results emphasize that the theory of critical slowing down applies only when the external forcing of the system across a critical point is slow relative to the internal system dynamics.

Demographic noise can reverse the direction of deterministic selection.

Deterministic evolutionary theory robustly predicts that populations displaying altruistic behaviors will be driven to extinction by mutant cheats that absorb common benefits but do not themselves contribute. Here we show that when demographic stochasticity is accounted for, selection can in fact act in the reverse direction to that predicted deterministically, instead favoring cooperative behaviors that appreciably increase the carrying capacity of the population. Populations that exist in larger numbers experience a selective advantage by being more stochastically robust to invasions than smaller populations, and this advantage can persist even in the presence of reproductive costs. We investigate this general effect in the specific context of public goods production and find conditions for stochastic selection reversal leading to the success of public good producers. This insight, developed here analytically, is missed by the deterministic analysis as well as by standard game theoretic models that enforce a fixed population size. The effect is found to be amplified by space; in this scenario we find that selection reversal occurs within biologically reasonable parameter regimes for microbial populations. Beyond the public good problem, we formulate a general mathematical framework for models that may exhibit stochastic selection reversal. In this context, we describe a stochastic analog to [Formula: see text] theory, by which small populations can evolve to higher densities in the absence of disturbance.

Spectral gap optimization of order parameters for sampling complex molecular systems.

In modern-day simulations of many-body systems, much of the computational complexity is shifted to the identification of slowly changing molecular order parameters called collective variables (CVs) or reaction coordinates. A vast array of enhanced-sampling methods are based on the identification and biasing of these low-dimensional order parameters, whose fluctuations are important in driving rare events of interest. Here, we describe a new algorithm for finding optimal low-dimensional CVs for use in enhanced-sampling biasing methods like umbrella sampling, metadynamics, and related methods, when limited prior static and dynamic information is known about the system, and a much larger set of candidate CVs is specified. The algorithm involves estimating the best combination of these candidate CVs, as quantified by a maximum path entropy estimate of the spectral gap for dynamics viewed as a function of that CV. The algorithm is called spectral gap optimization of order parameters (SGOOP). Through multiple practical examples, we show how this postprocessing procedure can lead to optimization of CV and several orders of magnitude improvement in the convergence of the free energy calculated through metadynamics, essentially giving the ability to extract useful information even from unsuccessful metadynamics runs.

A Kinetic Analysis of Coupled (or Auxiliary) Enzyme Reactions.

As a case study, we consider a coupled (or auxiliary) enzyme assay of two reactions obeying the Michaelis-Menten mechanism. The coupled reaction consists of a single-substrate, single-enzyme non-observable reaction followed by another single-substrate, single-enzyme observable reaction (indicator reaction). In this assay, the product of the non-observable reaction is the substrate of the indicator reaction. A mathematical analysis of the reaction kinetics is performed, and it is found that after an initial fast transient, the coupled reaction is described by a pair of interacting Michaelis-Menten equations. Moreover, we show that when the indicator reaction is fast, the quasi-steady-state dynamics are governed by three fast variables and one slow variable. Timescales that approximate the respective lengths of the indicator and non-observable reactions, as well as conditions for the validity of the Michaelis-Menten equations, are derived. The theory can be extended to deal with more complex sequences of enzyme-catalyzed reactions.

Frequency-dependent growth in class-structured populations: continuous dynamics in the limit of weak selection.

In this paper we consider class-structured populations in discrete time in the limit of weak selection and with the inverse of the intensity of selection as unit of time. The aim is to establish a continuous model that approximates the discrete model. More precisely, we study frequency-dependent growth in an infinite haploid population structured into a finite number of classes such that individuals in each class contribute to a given subset of classes from one time step to the next. These contributions take the form of generalized fecundity parameters with perturbations of order 1 / N that depends on the class frequencies of each type and the type frequencies. Moreover, they satisfy some mild conditions that ensure mixing in the long run. The dynamics in the limit as [Formula: see text] with N time steps as unit of time is considered first in the case of a single type, and second in the case of multiple types. The main result is that the type frequencies as [Formula: see text] obey the replicator equation with instantaneous growth rates for the different types that depend only on instantaneous equilibrium class frequencies and reproductive values. An application to evolutionary game theory complemented by simulation results is presented.

Resilience of the slow component in timescale-separated synchronized oscillators.

Physiological networks are usually made of a large number of biological oscillators evolving on a multitude of different timescales. Phase oscillators are particularly useful in the modelling of the synchronization dynamics of such systems. If the coupling is strong enough compared to the heterogeneity of the internal parameters, synchronized states might emerge where phase oscillators start to behave coherently. Here, we focus on the case where synchronized oscillators are divided into a fast and a slow component so that the two subsets evolve on separated timescales. We assess the resilience of the slow component by, first, reducing the dynamics of the fast one using Mori-Zwanzig formalism. Second, we evaluate the variance of the phase deviations when the oscillators in the two components are subject to noise with possibly distinct correlation times. From the general expression for the variance, we consider specific network structures and show how the noise transmission between the fast and slow components is affected. Interestingly, we find that oscillators that are among the most robust when there is only a single timescale, might become the most vulnerable when the system undergoes a timescale separation. We also find that layered networks seem to be insensitive to such timescale separations.

Quasi-robust control of biochemical reaction networks via stochastic morphing.

One of the main objectives of synthetic biology is the development of molecular controllers that can manipulate the dynamics of a given biochemical network that is at most partially known. When integrated into smaller compartments, such as living or synthetic cells, controllers have to be calibrated to factor in the intrinsic noise. In this context, biochemical controllers put forward in the literature have focused on manipulating the mean (first moment) and reducing the variance (second moment) of the target molecular species. However, many critical biochemical processes are realized via higher-order moments, particularly the number and configuration of the probability distribution modes (maxima). To bridge the gap, we put forward the stochastic morpher controller that can, under suitable timescale separations, morph the probability distribution of the target molecular species into a predefined form. The morphing can be performed at a lower-resolution, allowing one to achieve desired multi-modality/multi-stability, and at a higher-resolution, allowing one to achieve arbitrary probability distributions. Properties of the controller, such as robustness and convergence, are rigorously established, and demonstrated on various examples. Also proposed is a blueprint for an experimental implementation of stochastic morpher.

A structure preserving model order reduction method for calcium homeostatic system.

Calcium Homeostasis is a complex physiological process. Its mathematical model results in high order differential equation. In this paper, a model order reduction technique, based on time scale separation is proposed for a 27th order Calcium Homeostasis and Bone Remodeling (CHBR) system. The original state-space model after linearization has been decoupled into three reduced order subsystems: "Very-Slow", "Slow" and "Fast", at the same time preserving the structure of the system. The time and frequency response of individual reduced order model has been compared with the response of the original system. Furthermore, the effect of administering a therapeutic daily moderate dose of PTH has been studied with the help of reduced order models.

A theory of reactant-stationary kinetics for a mechanism of zymogen activation.

A theoretical analysis is performed on the nonlinear ordinary differential equations that govern the dynamics of a reaction mechanism of zymogen activation. The reaction consists of a primary non-observable zymogen activation reaction that it is coupled to an indicator (observable) reaction. The product of the first reaction is the enzyme of the indicator reaction, and both reactions are governed by the Michaelis-Menten reaction mechanism. Using singular perturbation methods, we derive asymptotic solutions that are valid under the quasi-steady-state and reactant-stationary assumptions. In particular, we obtain closed form solutions that are analogous to the Schnell-Mendoza equation for Michaelis-Menten type reactions. These closed-form solutions approximate the evolution of the observable reaction and provide the mathematical link necessary to measure the enzyme activity of the non-observable reaction. Conditions for the validity of the asymptotic solutions are also derived, and we demonstrate that these asymptotic expressions are applicable under reactant-stationary kinetics.

Approaches for the estimation of timescales in nonlinear dynamical systems: Timescale separation in enzyme kinetics as a case study.

The derivation of timescales is frequently introduced as an art form in papers and textbooks. The best scaling techniques require the application of physical intuition to identify dimensionless variables that are one unit order of magnitude and small parameters, which can simplify nonlinear differential equations. However, physical intuition requires prior knowledge of the solution to the dynamical systems under investigation. There are problems where the application of physical intuition is not straightforward. Therefore, it is necessary to apply mathematical techniques to estimate scales for the separation of timescales and simplification. In this review, we present three mathematical techniques - determination of pairwise balances, principle of minimum simplification and scaling by inverse rates - to scale dynamical systems with limited prior knowledge of model behavior. We illustrate the application of these techniques with the Michaelis-Menten reaction, which is widely studied to introduce scaling and simplification techniques in textbooks. We show that the pairwise balance approach, though commonly introduced as a method for nondimensionalization, can fail to derive a separation between timescales. The other techniques we review here can be applied to a number of dynamical systems, where the separation of timescales can lead to the simplification of a complex nonlinear problem.

Clusters of reaction rates and concentrations in protein networks such as the phosphotransferase system.

To understand the functioning of living cells, it is often helpful or even necessary to exploit inherent timescale disparities and focus on long-term dynamic behaviour. In the present study, we explore this type of behaviour for the biochemical network of the phosphotransferase system. We show that, during the slow phase that follows a fast initial transient, the network reaction rates are partitioned into clusters corresponding to connected parts of the reaction network. Rates within any of these clusters assume essentially the same value: differences within each cluster are vastly smaller than that from one cluster to another. This rate clustering induces an analogous clustering of the reactive compounds: only the molecular concentrations on the interface between these clusters are produced and consumed at substantially different rates and hence change considerably during the slow phase. The remaining concentrations essentially assume their steady-state values already by the end of the transient phase. Further, we find that this clustering phenomenon occurs for a large number of parameter values and also for models with different topologies; to each of these models, there corresponds a particular network partitioning. Our results show that, in spite of its complexity, the phosphotransferase system tends to behave in a rather simple (yet versatile) way. The persistence of clustering for the perturbed models we examined suggests that it is likely to be encountered in various environmental conditions, as well as in other signal transduction pathways with network structures similar to that of the phosphotransferase system.