Uncovering underlying physical principles and driving forces of cell differentiation and reprogramming from single-cell transcriptomics.

Recent advances in single-cell sequencing technology have revolutionized our ability to acquire whole transcriptome data. However, uncovering the underlying transcriptional drivers and nonequilibrium driving forces of cell function directly from these data remains challenging. We address this by learning cell state vector fields from discrete single-cell RNA velocity to quantify the single-cell global nonequilibrium driving forces as landscape and flux. From single-cell data, we quantified the Waddington landscape, showing that optimal paths for differentiation and reprogramming deviate from the naively expected landscape gradient paths and may not pass through landscape saddles at finite fluctuations, challenging conventional transition state estimation of kinetic rate for cell fate decisions due to the presence of the flux. A key insight from our study is that stem/progenitor cells necessitate greater energy dissipation for rapid cell cycles and self-renewal, maintaining pluripotency. We predict optimal developmental pathways and elucidate the nucleation mechanism of cell fate decisions, with transition states as nucleation sites and pioneer genes as nucleation seeds. The concept of loop flux quantifies the contributions of each cycle flux to cell state transitions, facilitating the understanding of cell dynamics and thermodynamic cost, and providing insights into optimizing biological functions. We also infer cell-cell interactions and cell-type-specific gene regulatory networks, encompassing feedback mechanisms and interaction intensities, predicting genetic perturbation effects on cell fate decisions from single-cell omics data. Essentially, our methodology validates the landscape and flux theory, along with its associated quantifications, offering a framework for exploring the physical principles underlying cellular differentiation and reprogramming and broader biological processes through high-throughput single-cell sequencing experiments.

Deep reaction network exploration of glucose pyrolysis.

Resolving the reaction networks associated with biomass pyrolysis is central to understanding product selectivity and aiding catalyst design to produce more valuable products. However, even the pyrolysis network of relatively simple [Formula: see text]-D-glucose remains unresolved due to its significant complexity in terms of the depth of the network and the number of major products. Here, a transition-state-guided reaction exploration has been performed that provides complete pathways to most significant experimental pyrolysis products of [Formula: see text]-D-glucose. The resulting reaction network involves over 31,000 reactions and transition states computed at the semiempirical quantum chemistry level and approximately 7,000 kinetically relevant reactions and transition states characterized with density function theory, comprising the largest reaction network reported for biomass pyrolysis. The exploration was conducted using graph-based rules to explore the reactivities of intermediates and an adaption of the Dijkstra algorithm to identify kinetically relevant intermediates. This simple exploration policy surprisingly (re)identified pathways to most major experimental pyrolysis products, many intermediates proposed by previous computational studies, and also identified new low-barrier reaction mechanisms that resolve outstanding discrepancies between reaction pathways and yields in isotope labeling experiments. This network also provides explanatory pathways for the high yield of hydroxymethylfurfural and the reaction pathway that contributes most to the formation of hydroxyacetaldehyde during glucose pyrolysis. Due to the limited domain knowledge required to generate this network, this approach should also be transferable to other complex reaction network prediction problems in biomass pyrolysis.

TS-tools: Rapid and automated localization of transition states based on a textual reaction SMILES input.

Here, TS-tools is presented, a Python package facilitating the automated localization of transition states (TS) based on a textual reaction SMILES input. TS searches can either be performed at xTB or DFT level of theory, with the former yielding guesses at marginal computational cost, and the latter directly yielding accurate structures at greater expense. On a benchmarking dataset of mono- and bimolecular reactions, TS-tools reaches an excellent success rate of 95% already at xTB level of theory. For tri- and multimolecular reaction pathways - which are typically not benchmarked when developing new automated TS search approaches, yet are relevant for various types of reactivity, cf. solvent- and autocatalysis and enzymatic reactivity - TS-tools retains its ability to identify TS geometries, though a DFT treatment becomes essential in many cases. Throughout the presented applications, a particular emphasis is placed on solvation-induced mechanistic changes, another issue that received limited attention in the automated TS search literature so far.

"On-the-fly" Crystal : How to reliably and automatically characterize and construct potential energy surfaces.

In this work, the Crystal  code, developed previously by the authors to find "holes" as well as legitimate transition states in existing potential energy surface (PES) functions [JPC Lett. 11, 6468 (2020)], is retooled to perform on-the-fly "direct dynamics"-type PES explorations, as well as automatic construction of new PES functions. In all of these contexts, the chief advantage of Crystal  over other methods is its ability to globally map the PES, thereby determining the most relevant regions of configuration space quickly and reliably-even when the dimensionality is rather large. Here, Crystal  is used to generate a uniformly spaced grid of density functional theory (DFT) or ab initio points, truncated over the relevant regions, which can then be used to either: (a) hone in precisely on PES features such as minima and transition states, or; (b) create a new PES function automatically, via interpolation. Proof of concept is demonstrated via application to three molecular systems: water (H 2 O), (reduced-dimensional) methane (CH 4 ), and methylene imine (CH 2 NH).

Periselectivity and ambimodal transition states in cycloadditions of tetrachloro-o-benzoquinone with 6,6-dimethylfulvene.

The reaction mechanism of cycloadditions of tetrachloro-o-benzoquinone with 6,6-dimethylfulvene were systematically investigated with density functional theory calculations. It was found that conditional primary interactions stabilize the ambimodal transition states in the endo pathways. Ambimodal transition states lead to [6 + 4]/[4 + 2] adducts or [4 + 2]/[2 + 4] adducts, which interconvert through 3,3-sigmatropic shift reactions. The substituent effects on periselectivity were also investigated.

Energetics and mechanisms for decomposition of cationized amino acids and peptides explored using guided ion beam tandem mass spectrometry.

Fragmentation studies of cationized amino acids and small peptides as studied using guided ion beam tandem mass spectrometry (GIBMS) are reviewed. After a brief examination of the key attributes of the GIBMS approach, results for a variety of systems are examined, compared, and contrasted. Cationization of amino acids, diglycine, and triglycine with alkali cations generally leads to dissociations in which the intact biomolecule is lost. Exceptions include most lithiated species as well as a few examples for sodiated and one example for potassiated species. Like the lithiated species, cationization by protons leads to numerous dissociation channels. Results for protonated glycine, cysteine, asparagine, diglycine, and a series of tripeptides are reviewed, along with the thermodynamic consequences that can be gleaned. Finally, the important physiological process of the deamidation of asparagine (Asn) residues is explored by the comparison of five dipeptides in which the C-terminal partner (AsnXxx) is altered. The GIBMS thermochemistry is shown to correlate well with kinetic results from solution phase studies.

Transition State Stabilizing Effects of Oxygen and Sulfur Chalcogen Bond Interactions.

Non-covalent chalcogen bond (ChB) interactions have found utility in many fields, including catalysis, organic semiconductors, and crystal engineering. In this study, the transition stabilizing effects of ChB interactions of oxygen and sulfur were experimentally measured using a series of molecular rotors. The rotors were designed to form ChB interactions in their bond rotation transition states. This enabled the kinetic influences to be assessed by monitoring changes in the rotational barriers. Despite forming weaker ChB interactions, the smaller chalcogens were able to stabilize transition states and had measurable kinetic effects on the rotational barriers. Sulfur stabilized the bond rotation transition state by as much as -7.2 kcal/mol without electron-withdrawing groups. The key was to design a system where the sulfur σ ${\sigma }$ -hole was aligned with the lone pairs of the chalcogen bond acceptor. Oxygen rotors also could form transition state stabilizing ChB interactions but required electron-withdrawing groups. For both oxygen and sulfur ChB interactions, a strong correlation was observed between transition state stabilizing abilities and electrostatic potential (ESP) of the chalcogen, providing a useful predictive parameter for the rational design of future ChB systems.

Observing the base-by-base search for native structure along transition paths during the folding of single nucleic acid hairpins.

Biomolecular folding involves searching among myriad possibilities for the native conformation, but the elementary steps expected from theory for this search have never been detected directly. We probed the dynamics of folding at high resolution using optical tweezers, measuring individual trajectories as nucleic acid hairpins passed through the high-energy transition states that dominate kinetics and define folding mechanisms. We observed brief but ubiquitous pauses in the transition states, with a dwell time distribution that matched microscopic theories of folding quantitatively. The sequence dependence suggested that pauses were dominated by microbarriers from nonnative conformations during the search by each nucleotide residue for the native base-pairing conformation. Furthermore, the pauses were position dependent, revealing subtle local variations in energy-landscape roughness and allowing the diffusion coefficient describing the microscopic dynamics within the barrier to be found without reconstructing the shape of the energy landscape. These results show how high-resolution measurements can elucidate key microscopic events during folding to test fundamental theories of folding.

Constant advance replicas method for locating minimum energy paths and transition states.

The chain-of-states (CoS) constant advance replicas (CAR) method and its climbing image variant (CI-CAR) for locating minimum energy paths (MEPs) and transition states are reported. The CAR algorithm applies the Lagrange multiplier method for imposing holonomic constraints on a chain-of-replicas, aiming to maintain equal mass-weighted/scaled root-mean-square (RMS) distances between the adjacent replicas by removing the sliding-down displacements contributed by the potential gradients during path optimization. Two contextual regularization schemes with clear geometrical interpretations are implemented to jointly promote high convergence and numerical robustness of the CAR algorithm. We show that the constrained reaction path can be solved normally within 5 steps of Lagrange multiplier updates with remarkably high numerical precision via the CAR approach. The efficacy of the CAR methods is demonstrated by testing on multiple analytical, classical, and quantum mechanical transition paths: the Müller potential, the alanine dipeptide isomerization, the helix unwinding of the VIVITLVMLKKK 12-mer peptide, and the Baker set of reactions. We also explore the potential of applying adaptive momentum (AdaM) optimizers for locating optimal transition paths under complex conformational changes. Most importantly, we discuss extensively the differences and connections between our newly proposed CAR methods and several related methods, with focuses on the reaction path with holonomic constraints (RPCons) approach of Brokaw et al. [J. Chem. Theory Comput. 2009, 5 (8), 2050-2061] and the state-of-the-art string method (SM) of E et al. [J. Chem. Phys. 2007, 126 (16), 164103]. The CAR approach represents a latest update to the general theoretical framework of reaction path finding algorithms in the two-ended CoS regime.

An automated method for graph-based chemical space exploration and transition state finding.

Algorithms that automatically explore the chemical space have been limited to chemical systems with a low number of atoms due to expensive involved quantum calculations and the large amount of possible reaction pathways. The method described here presents a novel solution to the problem of chemical exploration by generating reaction networks with heuristics based on chemical theory. First, a second version of the reaction network is determined through molecular graph transformations acting upon functional groups of the reacting. Only transformations that break two chemical bonds and form two new ones are considered, leading to a significant performance enhancement compared to previously presented algorithm. Second, energy barriers for this reaction network are estimated through quantum chemical calculations by a growing string method, which can also identify non-octet species missed during the previous step and further define the reaction network. The proposed algorithm has been successfully applied to five different chemical reactions, in all cases identifying the most important reaction pathways.

Recent Advances in Pathology: the 2022 Annual Review Issue of The Journal of Pathology.

The 2022 Annual Review Issue of The Journal of Pathology, Recent Advances in Pathology, contains 15 invited reviews on research areas of growing importance in pathology. This year, the articles include those that focus on digital pathology, employing modern imaging techniques and software to enable improved diagnostic and research applications to study human diseases. This subject area includes the ability to identify specific genetic alterations through the morphological changes they induce, as well as integrating digital and computational pathology with 'omics technologies. Other reviews in this issue include an updated evaluation of mutational patterns (mutation signatures) in cancer, the applications of lineage tracing in human tissues, and single cell sequencing technologies to uncover tumour evolution and tumour heterogeneity. The tissue microenvironment is covered in reviews specifically dealing with proteolytic control of epidermal differentiation, cancer-associated fibroblasts, field cancerisation, and host factors that determine tumour immunity. All of the reviews contained in this issue are the work of invited experts selected to discuss the considerable recent progress in their respective fields and are freely available online (https://onlinelibrary.wiley.com/journal/10969896). © 2022 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Chiral explosives: A theoretical investigation of structure and chiroptical properties of triacetone triperoxide and hexamethylene triperoxide diamine.

Triacetone triperoxide (TATP) and hexamethylene triperoxide diamine (HMTD) are cyclic peroxides that exhibit atropisomerism resulting from restricted rotation around three peroxide bonds. As a result, one pair of enantiomers with D3 symmetry and another pair of enantiomers with C2 symmetry can be identified. Previous studies, based on mass spectrometry data and computational results, have shown that conformations of TATP with D3 and C2 symmetry can be isolated. Assuming that enantiomer samples of TATP and HMTD can be obtained with sufficient enantiopurity, we investigated their chiroptical properties, namely, optical rotatory dispersion (ORD), vibrational circular dichroism (VCD), and Raman optical activity (VROA). ORD curves and VCD spectra are seen to be very similar for D3 - and C2 -symmetric atropisomers with the same overall helicity. Predicted VROA results, however, show significant differences between D3 - and C2 -symmetric atropisomers with the same overall helicity. The D3 -symmetric atropisomer is predicted to exhibit considerably larger magnitude vibrational optical activity signals than the C2 -symmetric atropisomer.

Designing an Apparently Orbital-Symmetry-Forbidden [3s,5s]-Sigmatropic Shift through Transition-State Complexation and Stereoelectronic Effects.

The [3s,5s]-sigmatropic shift is an example of an orbital-symmetry forbidden pericyclic reaction that is outcompeted by the allowed [3s,3s]-sigmatropic shift. Density functional theory calculations are used to show that PdII -complexed systems with strategically placed substituents engaging in key stereoelectronic effects can select for the [3s,5s] process, thereby outcompeting both orbital-symmetry-allowed [3s,3s]- and [3s,5a]-shifts.

Search for transition states with external forces.

It is much more difficult to find on the potential energy surface (PES) a transition state (TS) than a local minimum (LM). We propose a new methodology which makes this task much easier. Applying external forces to nuclei in a molecule we can locate on PES not just one particular TS but a number of consecutive transition states. With our approach it is possible to move over PES from one transition state to another without involving any local minima. The latter can be located in a separate step through the reaction path calculations performed for every transition state found before. Preliminary results for the 2-fluorofuran molecule illustrate the usefulness of the proposed method.

Combined in Silico and in Vitro Approaches To Uncover the Oxidation and Schiff Base Reaction of Baicalein as an Inhibitor of Amyloid Protein Aggregation.

The oxidized form of baicalein (BA) leads to covalent binding with human amyloid proteins. Such adducts hamper the aggregation and deposition of fibrils. A novel reaction of BA with pentylamine (PA) as a model for the lysine side chain is described. This is the first study addressing the atomistic details of a Schiff base reaction with the trihydroxylated moiety of BA. Nuclear magnetic resonance and mass spectrometry approaches clearly indicate the formation of dehydrobaicalein in solution as well as its condensation with PA under aerobic conditions, yielding regioselectively C6-substituted products. The combined results suggest initial ion pair formation between BA and PA, followed by a redox chain reaction: the initiation by oxygen/air; an o-quinone-based chain involving oxidation and reduction steps; and extra off-chain formation of a doubly oxidized product. These mechanistic details support the anti-amyloid activity of BA and endorse its trihydroxyphenyl moiety as a pharmacophore for drug-design studies.

Efficient and Accurate Description of Diels-Alder Reactions Using Density Functional Theory.

Modeling chemical reactions using Quantum Chemistry is a widely used predictive strategy capable to complement experiments in order to understand the intrinsic mechanisms guiding the chemicals towards the most favorable reaction products. However, at this purpose, it is mandatory to use reliable and computationally tractable theoretical methods. In this work, we focus on six Diels-Alder reactions of increasing complexity and perform an extensive benchmark of middle- to low-cost computational approaches to predict the characteristic reactions energy barriers. We found that Density Functional Theory, using the ωB97XD, LC-ωPBE, CAM-B3LYP, M11 and MN12SX functionals, with empirical dispersion corrections coupled to an affordable 6-31G basis set, provides quality results for this class of reactions, at a small computational effort. Such efficient and reliable simulation protocol opens perspectives for hybrid QM/MM molecular dynamics simulations of Diels-Alder reactions including explicit solvation.

Algorithmic Explorations of Unimolecular and Bimolecular Reaction Spaces.

Algorithmic reaction exploration based on transition state searches has already made inroads into many niche applications, but its potential as a general-purpose tool is still largely unrealized. Computational cost and the absence of benchmark problems involving larger molecules remain obstacles to further progress. Here an ultra-low cost exploration algorithm is implemented and used to explore the reactivity of unimolecular and bimolecular reactants, comprising a total of 581 reactions involving 51 distinct reactants. The algorithm discovers all established reaction pathways, where such comparisons are possible, while also revealing a much richer reactivity landscape, including lower barrier reaction pathways and a strong dependence of reaction conformation in the apparent barriers of the reported reactions. The diversity of these benchmarks illustrate that reaction exploration algorithms are approaching general-purpose capability.

Interconnectivity between Surface Reactivity and Self-Assembly of Kemp Elimination Catalyzing Nanorods.

Understanding the fundamental facts behind dynamicity of catalytic processes has been a longstanding quest across disciplines. Herein, we report self-assembly of catalytically active gold nanorods that can be regulated by tuning its reactivity towards a proton transfer reaction at different pH. Unlike substrate-induced templating and co-operativity, the enhanced aggregation rate is due to alteration of catalytic surface charge only during reactivity as negatively charged transition state of reactant (5-nitrobenzisoxazole) is formed on positively charged nanorod while undergoing a concerted E2-pathway. Herein, enhanced diffusivity during catalytic processes might also act as an additional contributing factor. Furthermore, we have also shown that nanosized hydrophobic cavities of clustered nanorods can also efficiently accelerate the rate of an aromatic nucleophilic substitution reaction, which also demonstrates a catalytic phenomenon that can lead to cascading of other reactions where substrates and products of the starting reactions are not directly involved.

autodE: Automated Calculation of Reaction Energy Profiles- Application to Organic and Organometallic Reactions.

Calculating reaction energy profiles to aid in mechanistic elucidation has long been the domain of the expert computational chemist. Here, we introduce autodE (https://github.com/duartegroup/autodE), an open-source Python package capable of locating transition states (TSs) and minima and delivering a full reaction energy profile from 1D or 2D chemical representations. autodE is broadly applicable to study organic and organometallic reaction classes, including addition, substitution, elimination, migratory insertion, oxidative addition, and reductive elimination; it accounts for conformational sampling of both minima and TSs and is compatible with many electronic structure packages. The general applicability of autodE is demonstrated in complex multi-step reactions, including cobalt- and rhodium-catalyzed hydroformylation and an Ireland-Claisen rearrangement.

Evolution of the Diels-Alder Reaction Mechanism since the 1930s: Woodward, Houk with Woodward, and the Influence of Computational Chemistry on Understanding Cycloadditions.

This review article describes the evolution of Woodward's mechanistic thinking, beginning in the late 1930s and early 1940s with his proposal of a charge-transfer mechanism for the Diels-Alder reaction, eventually leading to the Woodward-Katz two-stage concerted mechanism in 1959, and then to its mechanistic solution in terms of orbital symmetry control. Houk's research in the Woodward labs, testing the predictions of this theory, is described. Subsequent modern calculations with quantum mechanics and molecular dynamics simulations have shown that Woodward indeed had perfectly described not only the cyclopentadiene dimerization mechanism, but a new class of transition states now known as ambimodal or bis-pericyclic transition states. In recent years, the Houk group has found that ambimodal reactions are operative in the [6+4] cycloaddition. Molecular dynamics simulations of many Diels-Alder and ambimodal cycloadditions provide a time-resolved picture of how these reactions occur. Lastly, Roald Hoffmann provides a Coda in which he describes his joy in "being taken along the journey" of the cycloaddition story from Woodward's youth to today's trajectory simulations.