Toward Identification of the reaction coordinate directly from the transition state ensemble using the kernel PCA method.
Ontology highlight
ABSTRACT:
SUBMITTER: Antoniou D
PROVIDER: S-EPMC3058940 | biostudies-literature | 2011 Mar
REPOSITORIES: biostudies-literature
Publications
Toward Identification of the reaction coordinate directly from the transition state ensemble using the kernel PCA method.
Antoniou Dimitri D Schwartz Steven D SD
The journal of physical chemistry. B 20110221 10
We propose a new method for analyzing an ensemble of transition states to extract components of the reaction coordinate. We use the kernel principal component analysis (kPCA), which is a generalization of the ordinary PCA that does not make a linearization approximation We applied this method to a TPS study of human LDH we had previously published [Quaytman, S.; Schwartz, S. D. Proc. Natl. Acad. Sci. U.S.A.2007, 104, 12253] and extracted a reasonable representation for the reaction coordinate. ...[more]
Similar Datasets
Project description:Barrierless unimolecular association reactions are prominent in atmospheric and combustion mechanisms but are challenging for both experiment and kinetics theory. A key datum for understanding the pressure dependence of association and dissociation reactions is the high-pressure limit, but this is often available experimentally only by extrapolation. Here we calculate the high-pressure limit for the addition of a chlorine atom to acetylene molecule (Cl + C2H2→C2H2Cl). This reaction has outer and inner transition states in series; the outer transition state is barrierless, and it is necessary to use different theoretical frameworks to treat the two kinds of transition state. Here we study the reaction in the high-pressure limit using multifaceted variable-reaction-coordinate variational transition-state theory (VRC-VTST) at the outer transition state and reaction-path variational transition state theory (RP-VTST) at the inner turning point; then we combine the results with the canonical unified statistical (CUS) theory. The calculations are based on a density functional validated against the W3X-L method, which is based on coupled cluster theory with single, double, and triple excitations and a quasiperturbative treatment of connected quadruple excitations [CCSDT(Q)], and the computed rate constants are in good agreement with some of the experimental results. The chlorovinyl (C2H2Cl) adduct has two isomers that are equilibrium structures of a double-well C≡C-H bending potential. Two procedures are used to calculate the vibrational partition function of chlorovinyl; one treats the two isomers separately and the other solves the anharmonic energy levels of the double well. We use these results to calculate the standard-state free energy and equilibrium constant of the reaction.
| S-EPMC7084150 | biostudies-literature
Project description:The transition path sampling method previously applied in our group to the reaction catalyzed by lactate dehydrogenase was used to generate a transition path ensemble for this reaction. Based on analysis of the reactive trajectories generated, important residues behind the active site were implicated in a compressional motion that brought the donor-acceptor atoms of the hydride closer together. In addition, residues behind the active site were implicated in a relaxational motion, locking the substrate in product formation. Although this suggested that the compression-relaxation motions of these residues were important to catalysis, it remained unproven. In this work, we used committor distribution analysis to show that these motions are integral components of the reaction coordinate.
| S-EPMC1941458 | biostudies-literature
Project description:BackgroundKernel methods have been proven to be a powerful tool for the integration and analysis of high-throughput technologies generated data. Kernels offer a nonlinear version of any linear algorithm solely based on dot products. The kernelized version of principal component analysis is a valid nonlinear alternative to tackle the nonlinearity of biological sample spaces. This paper proposes a novel methodology to obtain a data-driven feature importance based on the kernel PCA representation of the data.ResultsThe proposed method, kernel PCA Interpretable Gradient (KPCA-IG), provides a data-driven feature importance that is computationally fast and based solely on linear algebra calculations. It has been compared with existing methods on three benchmark datasets. The accuracy obtained using KPCA-IG selected features is equal to or greater than the other methods' average. Also, the computational complexity required demonstrates the high efficiency of the method. An exhaustive literature search has been conducted on the selected genes from a publicly available Hepatocellular carcinoma dataset to validate the retained features from a biological point of view. The results once again remark on the appropriateness of the computed ranking.ConclusionsThe black-box nature of kernel PCA needs new methods to interpret the original features. Our proposed methodology KPCA-IG proved to be a valid alternative to select influential variables in high-dimensional high-throughput datasets, potentially unravelling new biological and medical biomarkers.
| S-EPMC10337068 | biostudies-literature
Project description:Digital image forgery is becoming easier to perform because of the rapid development of various manipulation tools. Image splicing is one of the most prevalent techniques. Digital images had lost their trustability, and researches have exerted considerable effort to regain such trustability by focusing mostly on algorithms. However, most of the proposed algorithms are incapable of handling high dimensionality and redundancy in the extracted features. Moreover, existing algorithms are limited by high computational time. This study focuses on improving one of the image splicing detection algorithms, that is, the run length run number algorithm (RLRN), by applying two dimension reduction methods, namely, principal component analysis (PCA) and kernel PCA. Support vector machine is used to distinguish between authentic and spliced images. Results show that kernel PCA is a nonlinear dimension reduction method that has the best effect on R, G, B, and Y channels and gray-scale images.
| S-EPMC4177829 | biostudies-other
Project description:Linear models offer a robust, flexible, and computationally efficient set of tools for modeling quantitative structure-activity relationships (QSARs) but have been eclipsed in performance by nonlinear methods. Support vector machines (SVMs) and neural networks are currently among the most popular and accurate QSAR methods because they learn new representations of the data that greatly improve modelability. In this work, we use shallow representation learning to improve the accuracy of L1 regularized logistic regression (LASSO) and meet the performance of Tanimoto SVM. We embedded chemical fingerprints in Euclidean space using Tanimoto (a.k.a. Jaccard) similarity kernel principal component analysis (KPCA) and compared the effects on LASSO and SVM model performance for predicting the binding activities of chemical compounds against 102 virtual screening targets. We observed similar performance and patterns of improvement for LASSO and SVM. We also empirically measured model training and cross-validation times to show that KPCA used in concert with LASSO classification is significantly faster than linear SVM over a wide range of training set sizes. Our work shows that powerful linear QSAR methods can match nonlinear methods and demonstrates a modular approach to nonlinear classification that greatly enhances QSAR model prototyping facility, flexibility, and transferability.
| S-EPMC5942586 | biostudies-literature
Project description:The concept of the protein transition state ensemble (TSE), a collection of the conformations that have 50% probability to convert rapidly to the folded state and 50% chance to rapidly unfold, constitutes the basis of the modern interpretation of protein engineering experiments. It has been conjectured that conformations constituting the TSE in many proteins are the expanded and distorted forms of the native state built around a specific folding nucleus. This view has been supported by a number of on-lattice and off-lattice simulations. Here we report a direct observation and characterization of the TSE by molecular dynamic folding simulations of the C-Src SH3 domain, a small protein that has been extensively studied experimentally. Our analysis reveals a set of key interactions between residues, conserved by evolution, that must be formed to enter the kinetic basin of attraction of the native state.
| S-EPMC1302427 | biostudies-other
Project description:Along with the development of information and technology, the quality characteristics of a product cannot be monitored separately in the different types of control charts. In the past, conventional control charts were developed to monitor only one type of quality characteristic. The variable control charts are used to observe the variable or metric quality characteristics. Meanwhile, in monitoring non-metric characteristics or categorical data, attribute control charts are employed. To accommodate these two types of data, the PCA Mix control chart is suggested to simultaneously monitor these two types of data in one chart. However, some drawbacks occur when this chart is applied to monitor non-metric data which has an imbalanced proportion. Therefore, the Kernel PCA Mix control chart is created to overcome the gaps that occurred in the PCA Mix chart. Similar to the previous chart, this chart is also constructed using Hotelling's T2 statistics with Kernel Density Estimation control limit. Several simulations are used to evaluate the performance of the proposed control charts. The simulation results show that the proposed chart has a better result than the previous control chart, especially for a small mean shift with an imbalanced proportion of non-metric data. However, the PCA Mix chart has a similar performance to the proposed chart when it is applied to monitor the balanced proportion of categorical data with a large mean shift. The application with simulated data with various scenarios and the real-world case also shows that the Kernel PCA Mix chart performs better compared to the performance of the PCA Mix chart.
| S-EPMC9489795 | biostudies-literature
Project description:The atmospheric reaction of H2S with Cl has been reinvestigated to check if, as previously suggested, only explicit dynamical computations can lead to an accurate evaluation of the reaction rate because of strong recrossing effects and the breakdown of the variational extension of transition state theory. For this reason, the corresponding potential energy surface has been thoroughly investigated, thus leading to an accurate characterization of all stationary points, whose energetics has been computed at the state of the art. To this end, coupled-cluster theory including up to quadruple excitations has been employed, together with the extrapolation to the complete basis set limit and also incorporating core-valence correlation, spin-orbit, and scalar relativistic effects as well as diagonal Born-Oppenheimer corrections. This highly accurate composite scheme has also been paralleled by less expensive yet promising computational approaches. Moving to kinetics, variational transition state theory and its variable reaction coordinate extension for barrierless steps have been exploited, thus obtaining a reaction rate constant (8.16 × 10-11 cm3 molecule-1 s-1 at 300 K and 1 atm) in remarkable agreement with the experimental counterpart. Therefore, contrary to previous claims, there is no need to invoke any failure of the transition state theory, provided that sufficiently accurate quantum-chemical computations are performed. The investigation of the puzzling case of the H2S + Cl system allowed us to present a robust approach for disclosing the thermochemistry and kinetics of reactions of atmospheric and astrophysical interest.
| S-EPMC8009477 | biostudies-literature
Project description:Characterizing the conformations of protein in the transition state ensemble (TSE) is important for studying protein folding. A promising approach pioneered by Vendruscolo et al. [Nature (London) 409, 641 (2001)] to study TSE is to generate conformations that satisfy all constraints imposed by the experimentally measured φ values that provide information about the native likeness of the transition states. Faísca et al. [J. Chem. Phys. 129, 095108 (2008)] generated conformations of TSE based on the criterion that, starting from a TS conformation, the probabilities of folding and unfolding are about equal through Markov Chain Monte Carlo (MCMC) simulations. In this study, we use the technique of constrained sequential Monte Carlo method [Lin et al., J. Chem. Phys. 129, 094101 (2008); Zhang et al. Proteins 66, 61 (2007)] to generate TSE conformations of acylphosphatase of 98 residues that satisfy the φ-value constraints, as well as the criterion that each conformation has a folding probability of 0.5 by Monte Carlo simulations. We adopt a two stage process and first generate 5000 contact maps satisfying the φ-value constraints. Each contact map is then used to generate 1000 properly weighted conformations. After clustering similar conformations, we obtain a set of properly weighted samples of 4185 candidate clusters. Representative conformation of each of these cluster is then selected and 50 runs of Markov chain Monte Carlo (MCMC) simulation are carried using a regrowth move set. We then select a subset of 1501 conformations that have equal probabilities to fold and to unfold as the set of TSE. These 1501 samples characterize well the distribution of transition state ensemble conformations of acylphosphatase. Compared with previous studies, our approach can access much wider conformational space and can objectively generate conformations that satisfy the φ-value constraints and the criterion of 0.5 folding probability without bias. In contrast to previous studies, our results show that transition state conformations are very diverse and are far from nativelike when measured in cartesian root-mean-square deviation (cRMSD): the average cRMSD between TSE conformations and the native structure is 9.4 Å for this short protein, instead of 6 Å reported in previous studies. In addition, we found that the average fraction of native contacts in the TSE is 0.37, with enrichment in native-like β-sheets and a shortage of long range contacts, suggesting such contacts form at a later stage of folding. We further calculate the first passage time of folding of TSE conformations through calculation of physical time associated with the regrowth moves in MCMC simulation through mapping such moves to a Markovian state model, whose transition time was obtained by Langevin dynamics simulations. Our results indicate that despite the large structural diversity of the TSE, they are characterized by similar folding time. Our approach is general and can be used to study TSE in other macromolecules.
| S-EPMC3071304 | biostudies-other
Project description:Since its discovery in 1979, left-handed Z-DNA has evolved from an in vitro curiosity to a challenging DNA structure with crucial roles in gene expression, regulation and recombination. A fundamental question that has puzzled researchers for decades is how the transition from B-DNA, the prevalent right-handed form of DNA, to Z-DNA is accomplished. Due to the complexity of the B-Z-DNA transition, experimental and computational studies have resulted in several different, apparently contradictory models. Here, we use molecular dynamics simulations coupled with state-of-the-art enhanced sampling techniques operating through non-conventional reaction coordinates, to investigate the B-Z-DNA transition at the atomic level. Our results show a complex free energy landscape, where several phenomena such as over-stretching, unpeeling, base pair extrusion and base pair flipping are observed resulting in interconversions between different DNA conformations such as B-DNA, Z-DNA and S-DNA. In particular, different minimum free energy paths allow for the coexistence of different mechanisms (such as zipper and stretch-collapse mechanisms) that previously had been proposed as independent, disconnected models. We find that the B-Z-DNA transition--in absence of other molecular partners--can encompass more than one mechanism of comparable free energy, and is therefore better described in terms of a reaction path ensemble.
| S-EPMC3592462 | biostudies-literature