Reactive Vortexes in a Naturally Activated Process: Non-Diffusive Rotational Fluxes at Transition State Uncovered by Persistent Homology

Abstract

The dynamics of reaction coordinates during barrier-crossing are key to understanding activated processes in complex systems such as proteins. The default assumption from Kramers' physical intuition is that of a diffusion process. However, the dynamics of barrier-crossing in natural complex molecules are largely unexplored. Here we investigate the transition dynamics of alanine dipeptide isomerization, the simplest complex system with a large number of non-reaction coordinates that can serve as an adequate thermal bath feeding energy into the reaction coordinates. We separate conformations along the time axis and construct the dynamic probability surface of reaction. We quantify its topological structure and rotational flux using persistent homology and differential form. Our results uncovered a region with a strong reactive vortex in the configuration-time space, where the highest probability peak and the transition state ensemble are located. This reactive region contains strong rotational fluxes: Most reactive trajectories swirl multiple times around this region in the subspace of the two most important reaction coordinates. Furthermore, the rotational fluxes result from cooperative movement along the isocommitter surfaces and orthogonal barrier-crossing. Overall, our findings offer a first glimpse into the reactive vortex regions that characterize the non-diffusive dynamics of barrier-crossing of a naturally occurring activation process.

Memory Unlocks the Future of Biomolecular Dynamics: Transformative Tools to Uncover Physical Insights Accurately and Efficiently

“…Thus, investigations of memory kernels may provide new insights into biomolecular mechanisms. For example, recent work has leveraged persistent homology analysis to identify the existence of flux vortices during the conformational transition of alanine dipeptide. , We expect that such vortices will leave an imprint on the memory kernel. Yet, the extent to which such features can be simply gleaned and interpreted from the structure of the non-Markovian generatorswhether these are memory kernels, their integrated moments, or their time-local analogsremains an open question.…”

Gentlest Ascent Dynamics on Manifolds Defined by Adaptively Sampled Point-Clouds

“…One interesting aspect of diffusion maps is its relation to the infinitesimal generator of a diffusion process, 25 which in turn has connections to the committor function, which is an optimal reaction coordinate. 26−28 We conclude this discussion on collective variables by pointing out that a priori knowledge of good reaction coordinates (some recent examples can be found in 29,30 ) can often be put in one-to-one correspondence with diffusion map coordinates. 31,32 Mean Force.…”

Biomolecular Electrostatic Phenomena: An Evergreen Field

“…A more theoretical contribution focused on the characterization of the nondiffusive dynamics of barrier-crossing of naturally occurring activation processes involving complex systems such as biomolecules. 42 In another contribution, Gaussian accelerated molecular dynamics, a method to accelerate the sampling of the conformational space, was implemented in the OpenMM simulation package, showing remarkable capabilities in reproducing energy profiles and folding pathways of model systems. 43 On the applicative side, the acetylcholine binding protein in a complex with nicotine was chosen as a model system to investigate the electrostatic components of nicotine binding, 44 and a simple predictor for the affinity of binding between the Receptor Binding Domain (RBD) of the spike protein of the infamous SARS-CoV-2 virus and its binding target, the angiotensin-converting enzyme 2, was obtained based on the total charge of the RBD.…”

“…The problem of free energy calculation is the focus of the last four works described in this synopsis. A more theoretical contribution focused on the characterization of the nondiffusive dynamics of barrier-crossing of naturally occurring activation processes involving complex systems such as biomolecules . In another contribution, Gaussian accelerated molecular dynamics, a method to accelerate the sampling of the conformational space, was implemented in the OpenMM simulation package, showing remarkable capabilities in reproducing energy profiles and folding pathways of model systems .…”