1
|
Diamataris IG, Peristeras LD, Papavasileiou KD, Melissas VS, Boulougouris GC. Statistical Inference of Rate Constants in Chemical and Biochemical Reaction Networks Using an "Inverse" Event-Driven Kinetic Monte Carlo Method. J Phys Chem B 2023; 127:9132-9143. [PMID: 37823789 DOI: 10.1021/acs.jpcb.3c03649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/13/2023]
Abstract
The use of rate models for networks of stochastic reactions is frequently used to comprehend the macroscopically observed dynamic properties of finite size reactive systems as well as their relationship to the underlying molecular events. Τhis particular approach usually stumbles on parameter derivation associated with stochastic kinetics, a quite demanding procedure. The present study incorporates a novel algorithm, which infers kinetic parameters from the system's time evolution, manifested as changes in molecular species populations. The proposed methodology reconstructs distributions required to infer kinetic parameters of a stochastic process pertaining to either a simulation or experimental data. The suggested approach accurately replicates rate constants of the stochastic reaction networks, which have evolved over time by event-driven Monte Carlo (MC) simulations using the Gillespie algorithm. Furthermore, our approach has been successfully used to estimate rate constants of association and dissociation events between molecular species developing during molecular dynamics (MD) simulations. We certainly believe that our method will be remarkably helpful for considering the macroscopic characteristic molecular roots related to stochastic physical and biological processes.
Collapse
Affiliation(s)
- Ioannis G Diamataris
- Laboratory of Computational Physical-Chemistry, Department of Molecular Biology and Genetics, University of Thrace, Alexandroupoulis GR-681 00, Greece
| | - Loukas D Peristeras
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modelling of Materials Laboratory, National Center for Scientific Research "Demokritos", Attikis, Agia Paraskevi GR-153 10, Greece
| | - Konstantinos D Papavasileiou
- Department of ChemoInformatics, NovaMechanics Ltd., Nicosia CY-1070, Cyprus
- Division of Data Driven Innovation, Entelos Institute, Larnaca CY-6059, Cyprus
- Department of ChemoInformatics, NovaMechanics MIKE., Piraeus GR-185 45, Greece
| | | | - Georgios C Boulougouris
- Laboratory of Computational Physical-Chemistry, Department of Molecular Biology and Genetics, University of Thrace, Alexandroupoulis GR-681 00, Greece
| |
Collapse
|
2
|
Frechette L, Stratt RM. The inherent dynamics of isotropic- and nematic-phase liquid crystals. J Chem Phys 2017; 144:234505. [PMID: 27334177 DOI: 10.1063/1.4953618] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
The geodesic (shortest) pathways through the potential energy landscape of a liquid can be thought of as defining what its dynamics would be if thermal noise were removed, revealing what we have called the "inherent dynamics" of the liquid. We show how these inherent paths can be located for a model liquid crystal former, showing, in the process, how the molecular mechanisms of translation and reorientation compare in the isotropic and nematic phases of these systems. These mechanisms turn out to favor the preservation of local orientational order even under macroscopically isotropic conditions (a finding consistent with the experimental observation of pseudonematic domains in these cases), but disfavor the maintenance of macroscopic orientational order, even in the nematic phase. While the most efficient nematic pathways that maintain nematic order are indeed shorter than those that do not, it is apparently difficult for the system to locate these paths, suggesting that molecular motion in liquid-crystal formers is dynamically frustrated, and reinforcing the sense that there are strong analogies between liquid crystals and supercooled liquids.
Collapse
Affiliation(s)
- Layne Frechette
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Richard M Stratt
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
3
|
Lempesis N, Boulougouris GC, Theodorou DN. Temporal disconnectivity of the energy landscape in glassy systems. J Chem Phys 2013; 138:12A545. [PMID: 23556796 DOI: 10.1063/1.4792363] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
An alternative graphical representation of the potential energy landscape (PEL) has been developed and applied to a binary Lennard-Jones glassy system, providing insight into the unique topology of the system's potential energy hypersurface. With the help of this representation one is able to monitor the different explored basins of the PEL, as well as how--and mainly when--subsets of basins communicate with each other via transitions in such a way that details of the prior temporal history have been erased, i.e., local equilibration between the basins in each subset has been achieved. In this way, apart from detailed information about the structure of the PEL, the system's temporal evolution on the PEL is described. In order to gather all necessary information about the identities of two or more basins that are connected with each other, we consider two different approaches. The first one is based on consideration of the time needed for two basins to mutually equilibrate their populations according to the transition rate between them, in the absence of any effect induced by the rest of the landscape. The second approach is based on an analytical solution of the master equation that explicitly takes into account the entire explored landscape. It is shown that both approaches lead to the same result concerning the topology of the PEL and dynamical evolution on it. Moreover, a "temporal disconnectivity graph" is introduced to represent a lumped system stemming from the initial one. The lumped system is obtained via a specially designed algorithm [N. Lempesis, D. G. Tsalikis, G. C. Boulougouris, and D. N. Theodorou, J. Chem. Phys. 135, 204507 (2011)]. The temporal disconnectivity graph provides useful information about both the lumped and the initial systems, including the definition of "metabasins" as collections of basins that communicate with each other via transitions that are fast relative to the observation time. Finally, the two examined approaches are compared to an "on the fly" molecular dynamics-based algorithm [D. G. Tsalikis, N. Lempesis, G. C. Boulougouris, and D. N. Theodorou, J. Chem. Theory Comput. 6, 1307 (2010)].
Collapse
Affiliation(s)
- Nikolaos Lempesis
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece
| | | | | |
Collapse
|
4
|
Lempesis N, Tsalikis DG, Boulougouris GC, Theodorou DN. Lumping analysis for the prediction of long-time dynamics: From monomolecular reaction systems to inherent structure dynamics of glassy materials. J Chem Phys 2011; 135:204507. [DOI: 10.1063/1.3663207] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
|
5
|
Mosayebi M, Del Gado E, Ilg P, Ottinger HC. Probing a critical length scale at the glass transition. PHYSICAL REVIEW LETTERS 2010; 104:205704. [PMID: 20867040 DOI: 10.1103/physrevlett.104.205704] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2010] [Indexed: 05/29/2023]
Abstract
We give evidence of a clear structural signature of the glass transition, in terms of a static correlation length with the same dependence on the system size, which is typical of critical phenomena. Our approach is to introduce an external, static perturbation to extract the structural information from the system's response. In particular, we consider the transformation behavior of the local minima of the underlying potential energy landscape (inherent structures), under a static deformation. The finite-size scaling analysis of our numerical results indicate that the correlation length diverges at a temperature Tc, below the temperatures where the system can be equilibrated. Our numerical results are consistent with random first order theory, which predicts such a divergence with a critical exponent ν=2/3 at the Kauzmann temperature, where the extrapolated configurational entropy vanishes.
Collapse
Affiliation(s)
- Majid Mosayebi
- Polymer Physics, ETH Zürich, Department of Materials, CH-8093 Zürich, Switzerland
| | | | | | | |
Collapse
|
6
|
Tsalikis DG, Lempesis N, Boulougouris GC, Theodorou DN. Temperature Accelerated Dynamics in Glass-Forming Materials. J Phys Chem B 2010; 114:7844-53. [DOI: 10.1021/jp908975d] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Dimitrios G. Tsalikis
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, Engineering Informatics and Telecommunications, University of Western Macedonia, Konstantinou Karamanli 55,GR-50100 Kozani, Greece, Department of Chemical Engineering, University of Patras, GR-26500 Patras, Greece, and Scienomics SARL, 17, Square Edouard VII, 75009 Paris
| | - Nikolaos Lempesis
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, Engineering Informatics and Telecommunications, University of Western Macedonia, Konstantinou Karamanli 55,GR-50100 Kozani, Greece, Department of Chemical Engineering, University of Patras, GR-26500 Patras, Greece, and Scienomics SARL, 17, Square Edouard VII, 75009 Paris
| | - Georgios C. Boulougouris
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, Engineering Informatics and Telecommunications, University of Western Macedonia, Konstantinou Karamanli 55,GR-50100 Kozani, Greece, Department of Chemical Engineering, University of Patras, GR-26500 Patras, Greece, and Scienomics SARL, 17, Square Edouard VII, 75009 Paris
| | - Doros N. Theodorou
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, Engineering Informatics and Telecommunications, University of Western Macedonia, Konstantinou Karamanli 55,GR-50100 Kozani, Greece, Department of Chemical Engineering, University of Patras, GR-26500 Patras, Greece, and Scienomics SARL, 17, Square Edouard VII, 75009 Paris
| |
Collapse
|
7
|
Tsalikis DG, Lempesis N, Boulougouris GC, Theodorou DN. Efficient Parallel Decomposition of Dynamical Sampling in Glass-Forming Materials Based on an “On the Fly” Definition of Metabasins. J Chem Theory Comput 2010. [DOI: 10.1021/ct9004245] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Dimitrios G. Tsalikis
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, Engineering Informatics and Telecommunications, University of Western Macedonia, Konstantinou Karamanli 55, GR-50100 Kozani, Greece, Department of Chemical Engineering, University of Patras, GR-26500 Patras, Greece, and Scienomics SARL, 17, Square Edouard VII, 75009 Paris
| | - Nikolaos Lempesis
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, Engineering Informatics and Telecommunications, University of Western Macedonia, Konstantinou Karamanli 55, GR-50100 Kozani, Greece, Department of Chemical Engineering, University of Patras, GR-26500 Patras, Greece, and Scienomics SARL, 17, Square Edouard VII, 75009 Paris
| | - Georgios C. Boulougouris
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, Engineering Informatics and Telecommunications, University of Western Macedonia, Konstantinou Karamanli 55, GR-50100 Kozani, Greece, Department of Chemical Engineering, University of Patras, GR-26500 Patras, Greece, and Scienomics SARL, 17, Square Edouard VII, 75009 Paris
| | - Doros N. Theodorou
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, Engineering Informatics and Telecommunications, University of Western Macedonia, Konstantinou Karamanli 55, GR-50100 Kozani, Greece, Department of Chemical Engineering, University of Patras, GR-26500 Patras, Greece, and Scienomics SARL, 17, Square Edouard VII, 75009 Paris
| |
Collapse
|
8
|
Tsalikis DG, Lempesis N, Boulougouris GC, Theodorou DN. On the Role of Inherent Structures in Glass-Forming Materials: I. The Vitrification Process. J Phys Chem B 2008; 112:10619-27. [DOI: 10.1021/jp801296k] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dimitrios G. Tsalikis
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, and Scienomics, 17, Square Eduard VII, 75009 Paris, France
| | - Nikolaos Lempesis
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, and Scienomics, 17, Square Eduard VII, 75009 Paris, France
| | - Georgios C. Boulougouris
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, and Scienomics, 17, Square Eduard VII, 75009 Paris, France
| | - Doros N. Theodorou
- School of Chemical Engineering, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece, and Scienomics, 17, Square Eduard VII, 75009 Paris, France
| |
Collapse
|