scholarly journals Estimate of the Free Energy Using Calculations of Dynamics in the Semiclassical Holstein Model

Author(s):  
В.Д. Лахно ◽  
V.D. Lakhno

We considered two ways of finding the free energy, using thermodynamic equilibrium characteristics, which are calculated by direct computational experiments. The results of calculations are described for homogeneous nucleotide dimers AA, GG and TT. Both of these methods yield similar results for biologically relevant temperatures.

2019 ◽  
Vol 20 (1) ◽  
Author(s):  
Guangyao Zhou ◽  
Jackson Loper ◽  
Stuart Geman

Abstract Background A folding RNA molecule encounters multiple opportunities to form non-native yet energetically favorable pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between “bound” molecules, which are active as part of ribonucleoprotein (RNP) complexes, and “unbound,” with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an “ambiguity index”—a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative (“gold standard”) and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules. Results We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule—an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure. Conclusions A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.


2019 ◽  
Author(s):  
Guangyao Zhou ◽  
Jackson Loper ◽  
Stuart Geman

Abstract Background : A folding RNA molecule encounters multiple opportunities to form non-native yet energetically favorable pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between "bound" molecules, which are active as part of ribonucleoprotein (RNP) complexes, and "unbound," with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an "ambiguity index"---a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative ("gold standard") and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules. Results : We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule---an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure. Conclusions : A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.


2019 ◽  
Author(s):  
Guangyao Zhou ◽  
Jackson Loper ◽  
Stuart Geman

Abstract Background: A folding RNA molecule encounters multiple opportunities to form non-native yet energetically favorable pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between "bound" molecules, which are active as part of ribonucleoprotein (RNP) complexes, and "unbound," with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an "ambiguity index"---a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative ("gold standard") and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules. Results: We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule---an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure. Conclusions: A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.


2019 ◽  
Author(s):  
Guangyao Zhou ◽  
Jackson Loper ◽  
Stuart Geman

Abstract Background : A folding RNA molecule encounters multiple opportunities to form non-native yet energetically favorable pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between "bound" molecules, which are active as part of ribonucleoprotein (RNP) complexes, and "unbound," with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an "ambiguity index"---a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative ("gold standard") and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules. Results : We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule---an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure. Conclusions : A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.


2019 ◽  
Author(s):  
Guangyao Zhou ◽  
Jackson Loper ◽  
Stuart Geman

Abstract Background: A folding RNA molecule encounters multiple opportunities to form non-native yet energetically favorable pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between "bound" molecules, which are active as part of ribonucleoprotein (RNP) complexes, and "unbound," with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an "ambiguity index"---a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative ("gold standard") and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules. Results: We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule---an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure. Conclusions: A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.


2018 ◽  
Author(s):  
Guangyao Zhou ◽  
Jackson Loper ◽  
Stuart Geman

AbstractBackgroundA folding RNA molecule encounters multiple opportunities to form non-native yet energetically favorable pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between “bound” molecules, which are active as part of ribonucleoprotein (RNP) complexes, and “unbound,” with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an “ambiguity index”—a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative (“gold standard”) and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules.ResultsWe have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule—an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure.ConclusionsA statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.


Sci ◽  
2019 ◽  
Vol 1 (3) ◽  
pp. 58

It is proposed that the primary living cells (“probionts”) cannot emerge of organic substance simply by continuous chemical complication of prebiotic macromolecules and microsystems. The complication must be accompanied by the radical thermodynamic transformation (“jump”) of prebiotic microsystems that resulted in the acquired ability to extract free energy from the environment and export entropy. This transformation is called “the thermodynamic inversion” The inversion may occur by means of the efficient (intensified) response of the microsystems on the oscillations of physic-chemical parameters in hydrothermal environment. In this case the surplus available free energy within a microsystem, when combined with the informational modality, facilitates its conversion into a new microsystem—a living probiont. It is shown the schematic representation of an oscillating prebiotic microsystem that is transforming into a living probiont. A new kind of laboratory and computational experiments on prebiotic chemistry under oscillating conditions is offered to verify the inversion concept.


Sci ◽  
2020 ◽  
Vol 2 (1) ◽  
pp. 14
Author(s):  
Vladimir Kompanichenko

It is proposed that the primary living cells (“probionts”) cannot emerge of organic substance simply by continuous chemical complication of prebiotic macromolecules and microsystems. The complication must be accompanied by the radical thermodynamic transformation (“jump”) of prebiotic microsystems that resulted in the acquired ability to extract free energy from the environment and export entropy. This transformation is called “the thermodynamic inversion” The inversion may occur by means of the efficient (intensified) response of the microsystems on the oscillations of physic-chemical parameters in hydrothermal environment. In this case the surplus available free energy within a microsystem, when combined with the informational modality, facilitates its conversion into a new microsystem—a living probiont. It is shown the schematic representation of an oscillating prebiotic microsystem that is transforming into a living probiont. A new kind of laboratory and computational experiments on prebiotic chemistry under oscillating conditions is offered to verify the inversion concept.


2019 ◽  
Vol 18 (02) ◽  
pp. 1950012 ◽  
Author(s):  
Wojciech Plazinski ◽  
Karolina Gaweda ◽  
Anita Plazinska

The conformation of five-membered furanose rings is a crucial issue for the structural analysis of many biologically-relevant molecules, including DNA and RNA. Oxolane can be treated as a prototypical furanose, composed only of saturated unsubstituted ring. In spite of its structural simplicity, providing the accurate quantitative description of the oxolane conformational features remains a great challenge for both the experimental and theoretical techniques. Here we show the method of recovering the free-energy profiles describing the conformational equilibrium in the oxolane ring (i.e. pseudorotation) based on the experimentally-inferred NMR data ([Formula: see text] coupling constants). The results remain in agreement with the quantum-mechanical-based molecular dynamics simulations and emphasize the large contributions of all ring conformers, even those located at the free-energy barriers. This includes the significant populations of limiting 3T2/2T3 and OE/EO shapes. Our findings provide another example of a poor applicability of the two-state model, which is routinely applied to analyze the NMR data in terms of population of different ring conformers.


2021 ◽  
Author(s):  
Yingying Wang ◽  
Xi Chen ◽  
Katharina Spengler ◽  
Karoline Terberger ◽  
Marko Boehm ◽  
...  

The decarboxylation of pyruvate is a central reaction in the carbon metabolism of all organisms. Both the pyruvate:ferredoxin oxidoreductase (PFOR) and the pyruvate dehydrogenase (PDH) complex catalyze this reaction. Whereas PFOR reduces ferredoxin, the PDH complex utilizes NAD+. Anaerobes rely on PFOR, which was replaced during evolution by the PDH complex found in aerobes. Cyanobacteria possess both. Our data challenge the view that PFOR is exclusively utilized for fermentation. Instead, we show, that the cyanobacterial PFOR is stable in the presence of oxygen in vitro and is required for optimal photomixotrophic growth under aerobic conditions while the PDH complex is inactivated under the same conditions. We found that cells rely on a general shift from utilizing NAD(H)-dependent to ferredoxin-dependent enzymes under these conditions. The utilization of ferredoxins instead of NAD(H) saves a greater share of the Gibbs free energy, instead of wasting it as heat. This obviously simultaneously decelerates metabolic reactions as they operate closer to their thermodynamic equilibrium. It is common thought that during evolution, ferredoxins were replaced by NAD(P)H due to their higher stability in an oxidizing atmosphere. However, utilization of NAD(P)H could also have been favored due to a higher competitiveness because of an accelerated metabolism.


Sign in / Sign up

Export Citation Format

Share Document