Solvent reaction coordinate for an SN2 reaction

Enhanced Sampling of Configuration and Path Space in a Generalized Ensemble by Shooting Point Exchange

The computer simulation of many molecular processes is complicated by long timescales caused by rare transitions between long-lived states. Here, we propose a new approach to simulate such rare events, which combines transition path sampling with enhanced exploration of configuration space. The method relies on exchange moves between configuration and trajectory space, carried out based on a generalized ensemble. This scheme substantially enhances the efficiency of the transition path sampling simulations, particularly for systems with multiple transition channels, and yields information on thermodynamics, kinetics and reaction coordinates of molecular processes without distorting their dynamics. The method is illustrated using the isomerization of proline in the KPTP tetrapeptide.

An exploration of machine learning models for the determination of reaction coordinates associated with conformational transitions

Determining collective variables (CVs) for conformational transitions is crucial to understanding their dynamics and targeting them in enhanced sampling simulations. Often, CVs are proposed based on intuition or prior knowledge of a system. However, the problem of systematically determining a proper reaction coordinate (RC) for a specific process in terms of a set of putative CVs can be achieved using committor analysis (CA). Identifying essential degrees of freedom that govern such transitions using CA remains elusive because of the high dimensionality of the conformational space. Various schemes exist to leverage the power of machine learning (ML) to extract an RC from CA. Here, we extend these studies and compare the ability of 17 different ML schemes to identify accurate RCs associated with conformational transitions. We tested these methods on an alanine dipeptide in vacuum and on a sarcosine dipeptoid in an implicit solvent. Our comparison revealed that the light gradient boosting machine method outperforms other methods. In order to extract key features from the models, we employed Shapley Additive exPlanations analysis and compared its interpretation with the "feature importance" approach. For the alanine dipeptide, our methodology identifies ϕ and θ dihedrals as essential degrees of freedom in the C7ax to C7eq transition. For the sarcosine dipeptoid system, the dihedrals ψ and ω are the most important for the cisαD to transαD transition. We further argue that analysis of the full dynamical pathway, and not just endpoint states, is essential for identifying key degrees of freedom governing transitions.

The SN2 reaction and its relationship with the Walden inversion, the Finkelstein and Menshutkin reactions together with theoretical calculations for the Finkelstein reaction

This communication gives an overview of the relationships between four reactions that although related were not always perceived as such: SN2, Walden, Finkelstein, and Menshutkin. Binary interactions (SN2 & Walden, SN2 & Menshutkin, SN2 & Finkelstein, Walden & Menshutkin, Walden & Finkelstein, Menshutkin & Finkelstein) were reported. Carbon, silicon, nitrogen, and phosphorus as central atoms and fluorides, chlorides, bromides, and iodides as lateral atoms were considered. Theoretical calculations provide Gibbs free energies that were analyzed with linear models to obtain the halide contributions. The M06-2x DFT computational method and the 6-311++G(d,p) basis set have been used for all atoms except for iodine where the effective core potential def2-TZVP basis set was used. Concerning the central atom pairs, carbon/silicon vs. nitrogen/phosphorus, we reported here for the first time that the effect of valence expansion was known for Si but not for P. Concerning the lateral halogen atoms, some empirical models including the interaction between F and I as entering and leaving groups explain the Gibbs free energies.

Transition Path Sampling Study of the Feruloyl Esterase Mechanism

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine-histidine-aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the "moving histidine" mechanism, thus calling into question the "ring flip" mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s^-1. Overall, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

Reaction Coordinate and Thermodynamics of Base Flipping in RNA

Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.