Multidimensional Free Energy and Accelerated Quantum Library Methods Provide a Gateway to Glycoenzyme Conformational, Electronic, and Reaction Mechanisms

Insights into the missing apiosylation step in flavonoid apiosides biosynthesis of Leguminosae plants

Apiose is a natural pentose containing an unusual branched-chain structure. Apiosides are bioactive natural products widely present in the plant kingdom. However, little is known on the key apiosylation reaction in the biosynthetic pathways of apiosides. In this work, we discover an apiosyltransferase GuApiGT from Glycyrrhiza uralensis. GuApiGT could efficiently catalyze 2″-O-apiosylation of flavonoid glycosides, and exhibits strict selectivity towards UDP-apiose. We further solve the crystal structure of GuApiGT, determine a key sugar-binding motif (RLGSDH) through structural analysis and theoretical calculations, and obtain mutants with altered sugar selectivity through protein engineering. Moreover, we discover 121 candidate apiosyltransferase genes from Leguminosae plants, and identify the functions of 4 enzymes. Finally, we introduce GuApiGT and its upstream genes into Nicotiana benthamiana, and complete de novo biosynthesis of a series of flavonoid apiosides. This work reports an efficient phenolic apiosyltransferase, and reveals mechanisms for its sugar donor selectivity.

Navigating phase space transport with the origin-fate map

We introduce and demonstrate the usage of the origin-fate map (OFM) as a tool for the detailed investigation of phase space transport in reactant-product-type systems. For these systems, which exhibit clearly defined start and end states, it is possible to build a comprehensive picture of the lobe dynamics by considering backward and forward integration of sets of initial conditions to index their origin and fate. We illustrate the method and its utility in the study of a two degrees of freedom caldera potential with four exits, demonstrating that the OFM not only recapitulates results from classical manifold theory but even provides more detailed information about complex lobe structures. The OFM allows the detection of dynamically significant transitions caused by the creation of new lobes and is also able to guide the prediction of the position of unstable periodic orbits (UPOs). Further, we compute the OFM on the periodic orbit dividing surface (PODS) associated with the transition state of a caldera entrance, which allows for a powerful analysis of reactive trajectories. The intersection of the manifolds corresponding to this UPO with other manifolds in the phase space results in the appearance of lobes on the PODS, which are directly classified by the OFM. This allows computations of branching ratios and the exploration of a fractal cascade of lobes as the caldera is stretched, which results in fluctuations in the branching ratio and chaotic selectivity. The OFM is found to be a simple and very useful tool with a vast range of descriptive and quantitative applications.

Molecular mechanisms from reaction coordinate graph enabled multidimensional free energies illustrated on water dimer hydrogen bonding

Computing the free energies of molecular mechanisms in multidimensional space relies on combinations of geometrically complex reaction coordinates. We show how a graph theory implementation reduces complexity, and illustrate this on the arrangements of hydrogen bonding of a water dimer. The reaction coordinates and forces are computed using graphs that define the dependencies on the atoms in the Free Energy from Adaptive Reaction Coordinate Forces (FEARCF) library. The library can be interfaced with classical molecular dynamics as well as quantum molecular dynamics packages. Multidimensional interdependent reaction coordinates are constructed to produce complex free energy hypersurfaces. The reaction coordinates are graphed from atomic and molecular components to define points, distances, vectors, angles, planes and combinations thereof. The resultant free energy surfaces that are a function of the distance, angles, planes, and so on, can represent molecular mechanisms in reduced dimensions from the component atomic Cartesian coordinate degrees of freedom. The FEARCF library can be interfaced with any molecular package. Here, we demonstrate the link to NWChem to compute a hyperdimensional DFT (aug‐cc‐pVDZ basis set and X3LYP exchange correlation functionals) free energy space of a water dimer. Analysis of the water dimer free energy hypervolume reveals that while the chain and cyclic hydrogen bonding configurations are located in stable minimum energy wells, the bifurcated hydrogen bond configuration is a gateway to instability and dimer dissociation.