Dataset Information

Equilibrium and non-equilibrium furanose selection in the ribose isomerisation network.

ABSTRACT: The exclusive presence of β-D-ribofuranose in nucleic acids is still a conundrum in prebiotic chemistry, given that pyranose species are substantially more stable at equilibrium. However, a precise characterisation of the relative furanose/pyranose fraction at temperatures higher than about 50 °C is still lacking. Here, we employ a combination of NMR measurements and statistical mechanics modelling to predict a population inversion between furanose and pyranose at equilibrium at high temperatures. More importantly, we show that a steady temperature gradient may steer an open isomerisation network into a non-equilibrium steady state where furanose is boosted beyond the limits set by equilibrium thermodynamics. Moreover, we demonstrate that nonequilibrium selection of furanose is maximum at optimal dissipation, as gauged by the temperature gradient and energy barriers for isomerisation. The predicted optimum is compatible with temperature drops found in hydrothermal vents associated with extremely fresh lava flows on the seafloor.

SUBMITTER: Dass AV

PROVIDER: S-EPMC8115175 | biostudies-literature | 2021 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Equilibrium and non-equilibrium furanose selection in the ribose isomerisation network.

Dass Avinash Vicholous AV Georgelin Thomas T Westall Frances F Foucher Frédéric F De Los Rios Paolo P Busiello Daniel Maria DM Liang Shiling S Piazza Francesco F

Nature communications 20210512 1

The exclusive presence of β-D-ribofuranose in nucleic acids is still a conundrum in prebiotic chemistry, given that pyranose species are substantially more stable at equilibrium. However, a precise characterisation of the relative furanose/pyranose fraction at temperatures higher than about 50 °C is still lacking. Here, we employ a combination of NMR measurements and statistical mechanics modelling to predict a population inversion between furanose and pyranose at equilibrium at high temperature ...[more]

PMID: 33980850

Similar Datasets

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.

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.

Dataset Information

Equilibrium and non-equilibrium furanose selection in the ribose isomerisation network.

Publications

Equilibrium and non-equilibrium furanose selection in the ribose isomerisation network.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets