Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Reaction coordinate of an enzymatic reaction revealed by transition path sampling.

ABSTRACT: 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.

SUBMITTER: Quaytman SL

PROVIDER: S-EPMC1941458 | biostudies-literature | 2007 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Reaction coordinate of an enzymatic reaction revealed by transition path sampling.

Quaytman Sara L SL Schwartz Steven D SD

Proceedings of the National Academy of Sciences of the United States of America 20070717 30

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 sub ...[more]

PMID: 17640885

Similar Datasets

Transition Path Sampling Study of the Reaction Catalyzed by Purine Nucleoside Phosphorylase.

Project description:The Transition Path Sampling (TPS) method is a powerful technique for studying rare events in complex systems, that allows description of reactive events in atomic detail without prior knowledge of reaction coordinates and transition states. We have applied TPS in combination with a hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) method to study the enzyme human purine nucleoside phosphorylase (hPNP). This enzyme catalyzes the reversible phosphorolysis of 6-oxypurine (deoxy)nucleosides to generate the corresponding purine base and (deoxy)ribose 1-phosphate. Hundreds of reactive trajectories were generated. Analysis of this transition path ensembles provides insight into the detailed mechanistic dynamics of reaction in the enzyme. Our studies have indicated a reaction mechanism involving the cleavage of the N-ribosidic bond to form transition states with substantial ribooxacarbenium ion character, that is then followed by conformational changes in the enzyme and the ribosyl group leading to migration of the anomeric carbon of the ribosyl group toward phosphate to form the product ribose 1-phosphate. This latter process is crucial in PNP, because several strong H-bonds form between active site residues in order to capture and align the phosphate nucleophile. Calculations of the commitment probability along reactive paths demonstrated the presence of a broad energy barrier at the transition state. Analysis of these transition state structures showed that bond-breaking and bond-forming distances are not a good choice for the reaction coordinate, but that the pseudorotational phase of the ribose ring is also a significant variable.

| S-EPMC2909608 | biostudies-literature

Association of Cl with C2H2 by unified variable-reaction-coordinate and reaction-path variational transition-state theory.

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

Mechanistic Insights on Human Phosphoglucomutase Revealed by Transition Path Sampling and Molecular Dynamics Calculations.

Project description:Human α-phosphoglucomutase 1 (α-PGM) catalyzes the isomerization of glucose-1-phosphate into glucose-6-phosphate (G6P) through two sequential phosphoryl transfer steps with a glucose-1,6-bisphosphate (G16P) intermediate. Given that the release of G6P in the gluconeogenesis raises the glucose output levels, α-PGM represents a tempting pharmacological target for type 2 diabetes. Here, we provide the first theoretical study of the catalytic mechanism of human α-PGM. We performed transition-path sampling simulations to unveil the atomic details of the two catalytic chemical steps, which could be key for developing transition state (TS) analogue molecules with inhibitory properties. Our calculations revealed that both steps proceed through a concerted SN 2-like mechanism, with a loose metaphosphate-like TS. Even though experimental data suggests that the two steps are identical, we observed noticeable differences: 1) the transition state ensemble has a well-defined TS region and a late TS for the second step, and 2) larger coordinated protein motions are required to reach the TS of the second step. We have identified key residues (Arg23, Ser117, His118, Lys389), and the Mg2+ ion that contribute in different ways to the reaction coordinate. Accelerated molecular dynamics simulations suggest that the G16P intermediate may reorient without leaving the enzymatic binding pocket, through significant conformational rearrangements of the G16P and of specific loop regions of the human α-PGM.

| S-EPMC5801119 | biostudies-literature

Transition-path sampling of beta-hairpin folding.

Project description:We examine the dynamical folding pathways of the C-terminal beta-hairpin of protein G-B1 in explicit solvent at room temperature by means of a transition-path sampling algorithm. In agreement with previous free-energy calculations, the resulting path ensembles reveal a folding mechanism in which the hydrophobic residues collapse first followed by backbone hydrogen-bond formation, starting with the hydrogen bonds inside the hydrophobic core. In addition, the path ensembles contain information on the folding kinetics, including solvent motion. Using the recently developed transition interface sampling technique, we calculate the rate constant for unfolding of the protein fragment and find it to be in reasonable agreement with experiments. The results support the validation of using all-atom force fields to study protein folding.

| S-EPMC218724 | biostudies-literature

Evolution alters the enzymatic reaction coordinate of dihydrofolate reductase.

Project description:How evolution has affected enzyme function is a topic of great interest in the field of biophysical chemistry. Evolutionary changes from Escherichia coli dihydrofolate reductase (ecDHFR) to human dihydrofolate reductase (hsDHFR) have resulted in increased catalytic efficiency and an altered dynamic landscape in the human enzyme. Here, we show that a subpicosecond protein motion is dynamically coupled to hydride transfer catalyzed by hsDHFR but not ecDHFR. This motion propagates through residues that correspond to mutational events along the evolutionary path from ecDHFR to hsDHFR. We observe an increase in the variability of the transition states, reactive conformations, and times of barrier crossing in the human system. In the hsDHFR active site, we detect structural changes that have enabled the coupling of fast protein dynamics to the reaction coordinate. These results indicate a shift in the DHFR family to a form of catalysis that incorporates rapid protein dynamics and a concomitant shift to a more flexible path through reactive phase space.

| S-EPMC4306496 | biostudies-literature

Enzymatic Kinetic Isotope Effects from First-Principles Path Sampling Calculations.

Project description:In this study, we develop and test a method to determine the rate of particle transfer and kinetic isotope effects in enzymatic reactions, specifically yeast alcohol dehydrogenase (YADH), from first-principles. Transition path sampling (TPS) and normal mode centroid dynamics (CMD) are used to simulate these enzymatic reactions without knowledge of their reaction coordinates and with the inclusion of quantum effects, such as zero-point energy and tunneling, on the transferring particle. Though previous studies have used TPS to calculate reaction rate constants in various model and real systems, it has not been applied to a system as large as YADH. The calculated primary H/D kinetic isotope effect agrees with previously reported experimental results, within experimental error. The kinetic isotope effects calculated with this method correspond to the kinetic isotope effect of the transfer event itself. The results reported here show that the kinetic isotope effects calculated from first-principles, purely for barrier passage, can be used to predict experimental kinetic isotope effects in enzymatic systems.

| S-EPMC4829465 | biostudies-literature

The enzymatic reaction catalyzed by lactate dehydrogenase exhibits one dominant reaction path.

Project description:Enzymes are the most efficient chemical catalysts known, but the exact nature of chemical barrier crossing in enzymes is not fully understood. Application of transition state theory to enzymatic reactions indicates that the rates of all possible reaction paths, weighted by their relative probabilities, must be considered in order to achieve an accurate calculation of the overall rate. Previous studies in our group have shown a single mechanism for enzymatic barrier passage in human heart lactate dehydrogenase (LDH). To ensure that this result was not due to our methodology insufficiently sampling reactive phase space, we implement high-perturbation transition path sampling in both microcanonical and canonical regimes for the reaction catalyzed by human heart LDH. We find that, although multiple, distinct paths through reactive phase space are possible for this enzymatic reaction, one specific reaction path is dominant. Since the frequency of these paths in a canonical ensemble is inversely proportional to the free energy barriers separating them from other regions of phase space, we conclude that the rarer reaction paths are likely to have a negligible contribution. Furthermore, the non-dominate reaction paths correspond to altered reactive conformations and only occur after multiple steps of high perturbation, suggesting that these paths may be the result of non-biologically significant changes to the structure of the enzymatic active site.

| S-EPMC4215548 | biostudies-literature

Catalytic Mechanism of Processive GlfT2: Transition Path Sampling Investigation of Substrate Translocation.

Project description:We applied the transition path sampling (TPS) method to study the translocation step of the catalytic mechanism of galactofuranosyl transferase 2 (GlfT2). Using TPS in the field of enzymatic reactions is still relatively rare, and we show its effectiveness on this enzymatic system. We decipher an unknown mechanism of the translocation step and, thus, provide a complete understanding of the catalytic mechanism of GlfT2 at the atomistic level. The GlfT2 enzyme is involved in the formation of the mycobacterial cell wall and transfers galactofuranose (Galf) from UDP-Galf onto a growing acceptor Galf chain. The biosynthesis of the galactan chain is accomplished in a processive manner, with the growing acceptor substrate remaining bound to GlfT2. The glycosidic bond formed by GlfT2 between the two Galf residues alternates between ?-(1-6) and ?-(1-5) linkages. The translocation of the growing galactan between individual additions of Galf residues is crucial for the function of GlfT2. Analysis of unbiased trajectory ensembles revealed that the translocation proceeds differently depending on the glycosidic linkage between the last two Galf residues. We also showed that the protonation state of the catalytic residue Asp372 significantly influences the translocation. Approximate transition state structures and potential energy reaction barriers of the translocation process were determined. The calculated potential reaction barriers in the range of 6-14 kcal/mol show that the translocation process is not the rate-limiting step in galactan biosynthesis.

| S-EPMC7469130 | biostudies-literature

Past-future information bottleneck for sampling molecular reaction coordinate simultaneously with thermodynamics and kinetics.

Project description:The ability to rapidly learn from high-dimensional data to make reliable bets about the future is crucial in many contexts. This could be a fly avoiding predators, or the retina processing gigabytes of data to guide human actions. In this work we draw parallels between these and the efficient sampling of biomolecules with hundreds of thousands of atoms. For this we use the Predictive Information Bottleneck framework used for the first two problems, and re-formulate it for the sampling of biomolecules, especially when plagued with rare events. Our method uses a deep neural network to learn the minimally complex yet most predictive aspects of a given biomolecular trajectory. This information is used to perform iteratively biased simulations that enhance the sampling and directly obtain associated thermodynamic and kinetic information. We demonstrate the method on two test-pieces, studying processes slower than milliseconds, calculating free energies, kinetics and critical mutations.

| S-EPMC6687748 | biostudies-literature

Elucidating the locking mechanism of peptides onto growing amyloid fibrils through transition path sampling.

Project description:We investigate the molecular mechanism of monomer addition to a growing amyloid fibril composed of the main amyloidogenic region from the insulin peptide hormone, the LVEALYL heptapeptide. Applying transition path sampling in combination with reaction coordinate analysis reveals that the transition from a docked peptide to a locked, fully incorporated peptide can occur in two ways. Both routes involve the formation of backbone hydrogen bonds between the three central amino acids of the attaching peptide and the fibril, as well as a reorientation of the central Glu side chain of the locking peptide toward the interface between two β-sheets forming the fibril. The mechanisms differ in the sequence of events. We also conclude that proper docking is important for correct alignment of the peptide with the fibril, as alternative pathways result in misfolding.

| S-EPMC3446680 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data