Gediminas Jurgis Pažėra, Thomas P. Fay, Ilia A. Solov’yov, P. J. Hore, Luca Gerhards
{"title":"Spin Dynamics of Radical Pairs Using the Stochastic Schrödinger Equation in MolSpin","authors":"Gediminas Jurgis Pažėra, Thomas P. Fay, Ilia A. Solov’yov, P. J. Hore, Luca Gerhards","doi":"10.1021/acs.jctc.4c00361","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00361","url":null,"abstract":"The chemical reactivity of radical pairs is strongly influenced by the interactions of electronic and nuclear spins. A detailed understanding of these effects requires a quantum description of the spin dynamics that considers spin-dependent reaction rates, interactions with external magnetic fields, spin–spin interactions, and the loss of spin coherence caused by coupling to a fluctuating environment. Modeling real chemical and biochemical systems, which frequently involve radicals with multinuclear spin systems, poses a severe computational challenge. Here, we implement a method based on the stochastic Schrödinger equation in the software package <i>MolSpin</i>. Large electron–nuclear spin systems can be simulated efficiently, with asymmetric spin-selective recombination reactions, anisotropic hyperfine interactions, and nonzero exchange and dipolar couplings. Spin-relaxation can be modeled using the stochastic time-dependence of spin interactions determined by molecular dynamics and quantum chemical calculations or by allowing rate coefficients to be explicitly time-dependent. The flexibility afforded by this approach opens new avenues for exploring the effects of complex molecular motions on the spin dynamics of chemical transformations.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.5,"publicationDate":"2024-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142234215","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Particle–Particle Random Phase Approximation for Predicting Correlated Excited States of Point Defects","authors":"Jiachen Li, Yu Jin, Jincheng Yu, Weitao Yang, Tianyu Zhu","doi":"10.1021/acs.jctc.4c00829","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00829","url":null,"abstract":"The particle–particle random phase approximation (ppRPA) within the hole–hole channel was recently proposed as an efficient tool for computing excitation energies of point defects in solids [<i>J. Phys. Chem. Lett.</i> 2024, 15, 2757–2764]. In this work, we investigate the application of ppRPA within the particle–particle channel for predicting correlated excited states of point defects, including the carbon-vacancy (VC) in diamond, the oxygen-vacancy (VO) in magnesium oxide (MgO), and the carbon dimer defect (C<sub>B</sub>C<sub>N</sub>) in two-dimensional hexagonal boron nitride (h-BN). Starting from a density functional theory calculation of the (<i>N</i> – 2)-electron ground state, vertical excitation energies of the <i>N</i>-electron system are obtained as the differences between the two-electron addition energies. We show that active-space ppRPA with the B3LYP functional yields accurate excitation energies, with errors mostly smaller than 0.1 eV for tested systems compared to available experimental values. We further develop a natural transition orbital scheme within ppRPA, which provides insights into the multireference character of defect states. This study, together with our previous work, establishes ppRPA as a low-cost and accurate method for investigating excited-state properties of point defect systems.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.5,"publicationDate":"2024-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142234256","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Aditya Koneru, Henry Chan, Sukriti Manna, Suvo Banik, Valeria Molinero, Subramanian K. R. S. Sankaranarayanan
{"title":"Machine Learning a Simple Interpretable Short-Range Potential for Silica","authors":"Aditya Koneru, Henry Chan, Sukriti Manna, Suvo Banik, Valeria Molinero, Subramanian K. R. S. Sankaranarayanan","doi":"10.1021/acs.jctc.4c00639","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00639","url":null,"abstract":"A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.5,"publicationDate":"2024-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142234216","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Guanxing Chen, Haohuai He, Qiujie Lv, Lu Zhao and Calvin Yu-Chian Chen*,
{"title":"MMFA-DTA: Multimodal Feature Attention Fusion Network for Drug-Target Affinity Prediction for Drug Repurposing Against SARS-CoV-2","authors":"Guanxing Chen, Haohuai He, Qiujie Lv, Lu Zhao and Calvin Yu-Chian Chen*, ","doi":"10.1021/acs.jctc.4c0066310.1021/acs.jctc.4c00663","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00663https://doi.org/10.1021/acs.jctc.4c00663","url":null,"abstract":"<p >The continuous emergence of novel infectious diseases poses a significant threat to global public health security, necessitating the development of small-molecule inhibitors that directly target pathogens. The RNA-dependent RNA polymerase (RdRp) and main protease (Mpro) of SARS-CoV-2 have been validated as potential key antiviral drug targets for the treatment of COVID-19. However, the conventional new drug R&D cycle takes 10–15 years, failing to meet the urgent needs during epidemics. Here, we propose a general multimodal deep learning framework for drug repurposing, MMFA-DTA, to enable rapid virtual screening of known drugs and significantly improve discovery efficiency. By extracting graph topological and sequence features from both small molecules and proteins, we design attention mechanisms to achieve dynamic fusion across modalities. Results demonstrate the superior performance of MMFA-DTA in drug-target affinity prediction over several state-of-the-art baseline methods on Davis and KIBA data sets, validating the benefits of heterogeneous information integration for representation learning and interaction modeling. Further fine-tuning on COVID-19-relevant bioactivity data enhances model predictions for critical SARS-CoV-2 enzymes. Case studies screening the FDA-approved drug library successfully identify etacrynic acid as the potential lead compound against both RdRp and Mpro. Molecular dynamics simulations further confirm the stability and binding affinity of etacrynic acid to these targets. This study proves the great potential and advantages of deep learning and drug repurposing strategies in supporting antiviral drug discovery. The proposed general and rapid response computational framework holds significance for preparedness against future public health events.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.7,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142309825","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"MMFA-DTA: Multimodal Feature Attention Fusion Network for Drug-Target Affinity Prediction for Drug Repurposing Against SARS-CoV-2","authors":"Guanxing Chen, Haohuai He, Qiujie Lv, Lu Zhao, Calvin Yu-Chian Chen","doi":"10.1021/acs.jctc.4c00663","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00663","url":null,"abstract":"The continuous emergence of novel infectious diseases poses a significant threat to global public health security, necessitating the development of small-molecule inhibitors that directly target pathogens. The RNA-dependent RNA polymerase (RdRp) and main protease (Mpro) of SARS-CoV-2 have been validated as potential key antiviral drug targets for the treatment of COVID-19. However, the conventional new drug R&D cycle takes 10–15 years, failing to meet the urgent needs during epidemics. Here, we propose a general multimodal deep learning framework for drug repurposing, MMFA-DTA, to enable rapid virtual screening of known drugs and significantly improve discovery efficiency. By extracting graph topological and sequence features from both small molecules and proteins, we design attention mechanisms to achieve dynamic fusion across modalities. Results demonstrate the superior performance of MMFA-DTA in drug-target affinity prediction over several state-of-the-art baseline methods on Davis and KIBA data sets, validating the benefits of heterogeneous information integration for representation learning and interaction modeling. Further fine-tuning on COVID-19-relevant bioactivity data enhances model predictions for critical SARS-CoV-2 enzymes. Case studies screening the FDA-approved drug library successfully identify etacrynic acid as the potential lead compound against both RdRp and Mpro. Molecular dynamics simulations further confirm the stability and binding affinity of etacrynic acid to these targets. This study proves the great potential and advantages of deep learning and drug repurposing strategies in supporting antiviral drug discovery. The proposed general and rapid response computational framework holds significance for preparedness against future public health events.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.5,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142198170","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Benchmark Study of Core-Ionization Energies with the Generalized Active Space-Driven Similarity Renormalization Group","authors":"Meng Huang, and , Francesco A. Evangelista*, ","doi":"10.1021/acs.jctc.4c0083510.1021/acs.jctc.4c00835","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00835https://doi.org/10.1021/acs.jctc.4c00835","url":null,"abstract":"<p >X-ray photoelectron spectroscopy (XPS) is a powerful experimental technique for probing the electronic structure of molecules and materials; however, interpreting XPS data requires accurate computational methods to model core-ionized states. This work proposes and benchmarks a new approach based on the generalized active space-driven similarity renormalization group (GAS-DSRG) for calculating core-ionization energies and treating correlation effects at the perturbative and nonperturbative levels. We tested the GAS-DSRG across three data sets. First, the vertical core-ionization energies of small molecules containing first-row elements are evaluated. GAS-DSRG achieves mean absolute errors below 0.3 eV, which is comparable to high-level coupled cluster methods. Next, the accuracy of GAS-DSRG is evaluated for larger organic molecules using the CORE65 data set, with the DSRG-MRPT3 level yielding a mean absolute error of only 0.34 eV for 65 core-ionization transitions. Insights are provided into the treatment of static and dynamic correlation, the importance of high-order perturbation theory, and notable differences from density functional theory in the predicted energy ordering of core-ionized states for specific molecules. Finally, vibrationally resolved XPS spectra of diatomic molecules (CO, N<sub>2</sub>, and O<sub>2</sub>) are simulated, showing excellent agreement with experimental data.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.7,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/epdf/10.1021/acs.jctc.4c00835","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142309820","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Benchmark Study of Core-Ionization Energies with the Generalized Active Space-Driven Similarity Renormalization Group","authors":"Meng Huang, Francesco A. Evangelista","doi":"10.1021/acs.jctc.4c00835","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00835","url":null,"abstract":"X-ray photoelectron spectroscopy (XPS) is a powerful experimental technique for probing the electronic structure of molecules and materials; however, interpreting XPS data requires accurate computational methods to model core-ionized states. This work proposes and benchmarks a new approach based on the generalized active space-driven similarity renormalization group (GAS-DSRG) for calculating core-ionization energies and treating correlation effects at the perturbative and nonperturbative levels. We tested the GAS-DSRG across three data sets. First, the vertical core-ionization energies of small molecules containing first-row elements are evaluated. GAS-DSRG achieves mean absolute errors below 0.3 eV, which is comparable to high-level coupled cluster methods. Next, the accuracy of GAS-DSRG is evaluated for larger organic molecules using the CORE65 data set, with the DSRG-MRPT3 level yielding a mean absolute error of only 0.34 eV for 65 core-ionization transitions. Insights are provided into the treatment of static and dynamic correlation, the importance of high-order perturbation theory, and notable differences from density functional theory in the predicted energy ordering of core-ionized states for specific molecules. Finally, vibrationally resolved XPS spectra of diatomic molecules (CO, N<sub>2</sub>, and O<sub>2</sub>) are simulated, showing excellent agreement with experimental data.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.5,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142198191","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Tuning Two-Photon Absorption in Rhodopsin Chromophore via Backbone Modification: The Story Told by CC2 and TD-DFT","authors":"Saruti Sirimatayanant, and , Tadeusz Andruniów*, ","doi":"10.1021/acs.jctc.4c0067510.1021/acs.jctc.4c00675","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00675https://doi.org/10.1021/acs.jctc.4c00675","url":null,"abstract":"<p >We investigate here a systematic way to tune two-photon transition strengths (δ<sup>2PA</sup>) and two-photon absorption (2PA) cross sections (σ<sup>2PA</sup>) of the rhodopsin’s chromophore 11-<i>cis</i>-retinal protonated Schiff base (RPSB) via the modulation of the methyl groups pattern along its polyene chain. Our team employed the resolution of identity, coupled cluster approximate second order (RI-CC2) method with Dunning’s aug-cc-pVDZ basis set, to determine the structural impact on δ<sup>2PA</sup>, as well as its correlation to both transition dipole moments and permanent electric dipole moments. Seven structures were probed in vacuo, including five-double-bond-conjugated model of the native chromophore, shortened by the β-ionone ring (RPSB5), and its de/methylated analogues: 9-methyl, 13-methyl, planar and twisted models of 9,10-dimethyl and 9,10,13-trimethyl. Our results demonstrate that the magnitude of δ<sup>2PA</sup> is dictated by both the position and number of methylated groups attached to its polyene chain as well as the degree of dihedral twist that is introduced due to the de/methylation. In fact, a strong correlation between δ<sup>2PA</sup> enhancement and the presence of a C13-methyl group in the planar RPSB5 species is found. Trends in δ<sup>2PA</sup> values follow the trends observed in their corresponding changes in the permanent dipole moment upon the S0–S1 excitation nearly exactly. The assessment of four DFT functionals, i.e., M11, MN15, CAM-B3LYP, and BHandHLYP, previously found most successful in predicting 2PA properties in biological chromophores, points to a long-range-corrected hybrid meta-GGA M11 as the top-performing functional, albeit still delivering underestimated δ<sup>2PA</sup> and σ<sup>2PA</sup> values by a factor of 3.3–5.3 with respect to the CC2 results. In the case of global-hybrid meta-NGA (MN15), as well as CAM-B3LYP and BHandHLYP functionals, this factor deteriorates significantly to 6.7–20.9 and is mostly related to significantly lower quality of the ground- and excited-state dipole moments.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.7,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/epdf/10.1021/acs.jctc.4c00675","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142309863","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Tuning Two-Photon Absorption in Rhodopsin Chromophore via Backbone Modification: The Story Told by CC2 and TD-DFT","authors":"Saruti Sirimatayanant, Tadeusz Andruniów","doi":"10.1021/acs.jctc.4c00675","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00675","url":null,"abstract":"We investigate here a systematic way to tune two-photon transition strengths (δ<sup>2PA</sup>) and two-photon absorption (2PA) cross sections (σ<sup>2PA</sup>) of the rhodopsin’s chromophore 11-<i>cis</i>-retinal protonated Schiff base (RPSB) via the modulation of the methyl groups pattern along its polyene chain. Our team employed the resolution of identity, coupled cluster approximate second order (RI-CC2) method with Dunning’s aug-cc-pVDZ basis set, to determine the structural impact on δ<sup>2PA</sup>, as well as its correlation to both transition dipole moments and permanent electric dipole moments. Seven structures were probed in vacuo, including five-double-bond-conjugated model of the native chromophore, shortened by the β-ionone ring (RPSB5), and its de/methylated analogues: 9-methyl, 13-methyl, planar and twisted models of 9,10-dimethyl and 9,10,13-trimethyl. Our results demonstrate that the magnitude of δ<sup>2PA</sup> is dictated by both the position and number of methylated groups attached to its polyene chain as well as the degree of dihedral twist that is introduced due to the de/methylation. In fact, a strong correlation between δ<sup>2PA</sup> enhancement and the presence of a C13-methyl group in the planar RPSB5 species is found. Trends in δ<sup>2PA</sup> values follow the trends observed in their corresponding changes in the permanent dipole moment upon the S0–S1 excitation nearly exactly. The assessment of four DFT functionals, i.e., M11, MN15, CAM-B3LYP, and BHandHLYP, previously found most successful in predicting 2PA properties in biological chromophores, points to a long-range-corrected hybrid meta-GGA M11 as the top-performing functional, albeit still delivering underestimated δ<sup>2PA</sup> and σ<sup>2PA</sup> values by a factor of 3.3–5.3 with respect to the CC2 results. In the case of global-hybrid meta-NGA (MN15), as well as CAM-B3LYP and BHandHLYP functionals, this factor deteriorates significantly to 6.7–20.9 and is mostly related to significantly lower quality of the ground- and excited-state dipole moments.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.5,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142198183","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Spectral Map for Slow Collective Variables, Markovian Dynamics, and Transition State Ensembles","authors":"Jakub Rydzewski*, ","doi":"10.1021/acs.jctc.4c0042810.1021/acs.jctc.4c00428","DOIUrl":"https://doi.org/10.1021/acs.jctc.4c00428https://doi.org/10.1021/acs.jctc.4c00428","url":null,"abstract":"<p >Understanding the behavior of complex molecular systems is a fundamental problem in physical chemistry. To describe the long-time dynamics of such systems, which is responsible for their most informative characteristics, we can identify a few slow collective variables (CVs) while treating the remaining fast variables as thermal noise. This enables us to simplify the dynamics and treat it as diffusion in a free-energy landscape spanned by slow CVs, effectively rendering the dynamics Markovian. Our recent statistical learning technique, spectral map [Rydzewski, J. <i>J. Phys. Chem. Lett.</i> <b>2023</b>, <i>14</i>(22), 5216–5220], explores this strategy to learn slow CVs by maximizing a spectral gap of a transition matrix. In this work, we introduce several advancements into our framework, using a high-dimensional reversible folding process of a protein as an example. We implement an algorithm for coarse-graining Markov transition matrices to partition the reduced space of slow CVs kinetically and use it to define a transition state ensemble. We show that slow CVs learned by spectral map closely approach the Markovian limit for an overdamped diffusion. We demonstrate that coordinate-dependent diffusion coefficients only slightly affect the constructed free-energy landscapes. Finally, we present how spectral maps can be used to quantify the importance of features and compare slow CVs with structural descriptors commonly used in protein folding. Overall, we demonstrate that a single slow CV learned by spectral map can be used as a physical reaction coordinate to capture essential characteristics of protein folding.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":null,"pages":null},"PeriodicalIF":5.7,"publicationDate":"2024-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/epdf/10.1021/acs.jctc.4c00428","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142309838","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}