Equilibrium and non-equilibrium furanose selection in the ribose isomerisation network.
Ontology highlight
ABSTRACT:
SUBMITTER: Dass AV
PROVIDER: S-EPMC8115175 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Life has most likely originated as a consequence of processes taking place in non-equilibrium conditions (e.g. in the proximity of deep-sea thermal vents) selecting states of matter that would have been otherwise unfavorable at equilibrium. Here we present a simple chemical network in which the selection of states is driven by the thermodynamic necessity of dissipating heat as rapidly as possible in the presence of a thermal gradient: states participating to faster reactions contribute the most to the dissipation rate, and are the most populated ones in non-equilibrium steady-state conditions. Building upon these results, we show that, as the complexity of the chemical network increases, the velocity of the reaction path leading to a given state determines its selection, giving rise to non-trivial localization phenomena in state space. A byproduct of our studies is that, in the presence of a temperature gradient, thermophoresis-like behavior inevitably appears depending on the transport properties of each individual state, thus hinting at a possible microscopic explanation of this intriguing yet still not fully understood phenomenon.
| S-EPMC9814615 | biostudies-literature
Project description:Unravelling underlying complex structures from limited resolution measurements is a known problem arising in many scientific disciplines. We study a stochastic dynamical model with a multiplicative noise. It consists of a stochastic differential equation living on a graph, similar to approaches used in population dynamics or directed polymers in random media. We develop a new tool for approximation of correlation functions based on spectral analysis that does not require translation invariance. This enables us to go beyond lattices and analyse general networks. We show, analytically, that this general model has different phases depending on the topology of the network. One of the main parameters which describe the network topology is the spectral dimension [Formula: see text]. We show that the correlation functions depend on the spectral dimension and that only for [Formula: see text] > 2 a dynamical phase transition occurs. We show by simulation how the system behaves for different network topologies, by defining and calculating the Lyapunov exponents on the graph. We present an application of this model in the context of Magnetic Resonance (MR) measurements of porous structure such as brain tissue. This model can also be interpreted as a KPZ equation on a graph.
| S-EPMC6156338 | biostudies-literature
Project description:Key requirements for the first cells on Earth include the ability to compartmentalize and evolve. Compartmentalization spatially localizes biomolecules from a dilute pool and an evolving cell, which, as it grows and divides, permits mixing and propagation of information to daughter cells. Complex coacervate microdroplets are excellent candidates as primordial cells with the ability to partition and concentrate molecules into their core and support primitive and complex biochemical reactions. However, the evolution of coacervate protocells by fusion, growth and fission has not yet been demonstrated. In this work, a primordial environment initiated the evolution of coacervate-based protocells. Gas bubbles inside heated rock pores perturb the coacervate protocell distribution and drive the growth, fusion, division and selection of coacervate microdroplets. Our findings provide a compelling scenario for the evolution of membrane-free coacervate microdroplets on the early Earth, induced by common gas bubbles within heated rock pores.
| S-EPMC8755537 | biostudies-literature
Project description:DNA polymerases select for the incorporation of deoxyribonucleotide triphosphates (dNTPs) using amino acid side-chains that act as a "steric-gate" to bar improper incorporation of rNTPs. An additional factor in the selection of nucleotide substrates resides in the preferred geometry for the furanose moiety of the incoming nucleotide triphosphate. We have probed the role of sugar geometry during nucleotide selection by model DNA polymerases from Sulfolobus solfataricus using fixed conformation nucleotide analogues. North-methanocarba-dATP (N-MC-dATP) locks the central ring into a RNA-type (C2'-exo, North) conformation near a C3'-endo pucker, and South-methanocarba-dATP (S-MC-dATP) locks the central ring system into a (C3'-exo, South) conformation near a C2'-endo pucker. Dpo4 preferentially inserts N-MC-dATP and in the crystal structure of Dpo4 in complex with N-MC-dAMP, the nucleotide analogue superimposes almost perfectly with Dpo4 bound to unmodified dATP. Biochemical assays indicate that the S. solfataricus B-family DNA polymerase Dpo1 can insert and extend from both N-MC-dATP and S-MC-dATP. In this respect, Dpo1 is unexpectedly more tolerant of substrate conformation than Dpo4. The crystal structure of Dpo4 bound to S-MC-dADP shows that poor incorporation of the Southern pucker by the Y-family polymerase results from a hydrogen bond between the 3'-OH group of the nucleotide analogue and the OH group of the steric gate residue, Tyr12, shifting the S-MC-dADP molecule away from the dNTP binding pocket and distorting the base pair at the primer-template junction. These results provide insights into substrate specificity of DNA polymerases, as well as molecular mechanisms that act as a barrier against insertion of rNTPs.
| S-EPMC3563936 | biostudies-literature
Project description:BackgroundBiochemical equilibria are usually modeled iteratively: given one or a few fitted models, if there is a lack of fit or over fitting, a new model with additional or fewer parameters is then fitted, and the process is repeated. The problem with this approach is that different analysts can propose and select different models and thus extract different binding parameter estimates from the same data. An alternative is to first generate a comprehensive standardized list of plausible models, and to then fit them exhaustively, or semi-exhaustively.ResultsA framework is presented in which equilibriums are modeled as pairs (g, h) where g = 0 maps total reactant concentrations (system inputs) into free reactant concentrations (system states) which h then maps into expected values of measurements (system outputs). By letting dissociation constants Kd be either freely estimated, infinity, zero, or equal to other Kd, and by letting undamaged protein fractions be either freely estimated or 1, many g models are formed. A standard space of g models for ligand-induced protein dimerization equilibria is given. Coupled to an h model, the resulting (g, h) were fitted to dTTP induced R1 dimerization data (R1 is the large subunit of ribonucleotide reductase). Models with the fewest parameters were fitted first. Thereafter, upon fitting a batch, the next batch of models (with one more parameter) was fitted only if the current batch yielded a model that was better (based on the Akaike Information Criterion) than the best model in the previous batch (with one less parameter). Within batches models were fitted in parallel. This semi-exhaustive approach yielded the same best models as an exhaustive model space fit, but in approximately one-fifth the time.ConclusionComprehensive model space based biochemical equilibrium model selection methods are realizable. Their significance to systems biology as mappings of data into mathematical models warrants their development.
| S-EPMC2268910 | biostudies-literature
Project description:Local Shannon entropy lies at the heart of modern thermodynamics, with much discussion of trajectory-dependent entropy production. When taken at both boundaries of a process in phase space, it reproduces the second law of thermodynamics over a finite time interval for small scale systems. However, given that entropy is an ensemble property, it has never been clear how one can assign such a quantity locally. Given such a fundamental omission in our knowledge, we construct a new ensemble composed of trajectories reaching an individual microstate, and show that locally defined entropy, information, and free energy are properties of the ensemble, or trajectory-independent true thermodynamic potentials. We find that the Boltzmann-Gibbs distribution and Landauer's principle can be generalized naturally as properties of the ensemble, and that trajectory-free state functions of the ensemble govern the exact mechanism of non-equilibrium relaxation.
| S-EPMC4296294 | biostudies-literature
Project description:Intra-molecular energy transport between distant functional sites plays important roles in allosterically regulating the biochemical activity of proteins. How to identify the specific intra-molecular signaling pathway from protein tertiary structure remains a challenging problem. In the present work, a non-equilibrium dynamics method based on the elastic network model (ENM) was proposed to simulate the energy propagation process and identify the specific signaling pathways within proteins. In this method, a given residue was perturbed and the propagation of energy was simulated by non-equilibrium dynamics in the normal modes space of ENM. After that, the simulation results were transformed from the normal modes space to the Cartesian coordinate space to identify the intra-protein energy transduction pathways. The proposed method was applied to myosin and the third PDZ domain (PDZ3) of PSD-95 as case studies. For myosin, two signaling pathways were identified, which mediate the energy transductions form the nucleotide binding site to the 50?kDa cleft and the converter subdomain, respectively. For PDZ3, one specific signaling pathway was identified, through which the intra-protein energy was transduced from ligand binding site to the distant opposite side of the protein. It is also found that comparing with the commonly used cross-correlation analysis method, the proposed method can identify the anisotropic energy transduction pathways more effectively.
| S-EPMC6015066 | biostudies-literature
Project description:The RNA world is one of the principal hypotheses to explain the emergence of living systems on the prebiotic Earth. It posits that RNA oligonucleotides acted as both carriers of information as well as catalytic molecules, promoting their own replication. However, it does not explain the origin of the catalytic RNA molecules. How could the transition from a pre-RNA to an RNA world occur? A starting point to answer this question is to analyze the dynamics in sequence space on the lowest level, where mononucleotide and short oligonucleotides come together and collectively evolve into larger molecules. To this end, we study the sequence-dependent self-assembly of polymers from a random initial pool of short building blocks via templated ligation. Templated ligation requires two strands that are hybridized adjacently on a third strand. The thermodynamic stability of such a configuration crucially depends on the sequence context and, therefore, significantly influences the ligation probability. However, the sequence context also has a kinetic effect, since non-complementary nucleotide pairs in the vicinity of the ligation site stall the ligation reaction. These sequence-dependent thermodynamic and kinetic effects are explicitly included in our stochastic model. Using this model, we investigate the system-level dynamics inside a non-equilibrium 'RNA reactor' enabling a fast chemical activation of the termini of interacting oligomers. Moreover, the RNA reactor subjects the oligomer pool to periodic temperature changes inducing the reshuffling of the system. The binding stability of strands typically grows with the number of complementary nucleotides forming the hybridization site. While shorter strands unbind spontaneously during the cold phase, larger complexes only disassemble during the temperature peaks. Inside the RNA reactor, strand growth is balanced by cleavage via hydrolysis, such that the oligomer pool eventually reaches a non-equilibrium stationary state characterized by its length and sequence distribution. How do motif-dependent energy and stalling parameters affect the sequence composition of the pool of long strands? As a critical factor for self-enhancing sequence selection, we identify kinetic stalling due to non-complementary base pairs at the ligation site. Kinetic stalling enables cascades of self-amplification that result in a strong reduction of occupied states in sequence space. Moreover, we discuss the significance of the symmetry breaking for the transition from a pre-RNA to an RNA world.
| S-EPMC9032526 | biostudies-literature
Project description:Island biotas emerge from the interplay between colonisation, speciation and extinction and are often the scene of spectacular adaptive radiations. A common assumption is that insular diversity is at a dynamic equilibrium, but for remote islands, such as Hawaii or Galápagos, this idea remains untested. Here, we reconstruct the temporal accumulation of terrestrial bird species of the Galápagos using a novel phylogenetic method that estimates rates of biota assembly for an entire community. We show that species richness on the archipelago is in an ascending phase and does not tend towards equilibrium. The majority of the avifauna diversifies at a slow rate, without detectable ecological limits. However, Darwin's finches form an exception: they rapidly reach a carrying capacity and subsequently follow a coalescent-like diversification process. Together, these results suggest that avian diversity of remote islands is rising, and challenge the mutual exclusivity of the non-equilibrium and equilibrium ecological paradigms.
| S-EPMC4745040 | biostudies-literature
Project description:We study mechanisms of synchronisation, coordination, and equilibrium selection in two-player coordination games on multilayer networks. We investigate three possible update rules: the replicator dynamics (RD), the best response (BR), and the unconditional imitation (UI). Players interact on a two-layer random regular network. The population on each layer plays a different game, with layer I preferring the opposite strategy to layer II. We measure the difference between the two games played on the layers by a difference in payoffs, and the inter-connectedness by a node overlap parameter. We discover a critical value of the overlap below which layers do not synchronise, i.e. they display different levels of coordination. Above this threshold both layers typically coordinate on the same strategy. Surprisingly, there is a symmetry breaking in the selection of equilibrium-for RD and UI there is a phase where only the payoff-dominant equilibrium is selected. It is not observed, however, for BR update rule. Our work is an example of previously observed differences between the update rules. Nonetheless, we took a novel approach with the game being played on two inter-connected layers. As we show, the multilayer structure enhances the abundance of the Pareto-optimal equilibrium in coordination games with imitative update rules.
| S-EPMC10362047 | biostudies-literature