Machine Learning Guided AQFEP: A Fast and Efficient Absolute Free Energy Perturbation Solution for Virtual Screening.

Structure-based methods in drug discovery have become an integral part of the modern drug discovery process. The power of virtual screening lies in its ability to rapidly and cost-effectively explore enormous chemical spaces to select promising ligands for further experimental investigation. Relative free energy perturbation (RFEP) and similar methods are the gold standard for binding affinity prediction in drug discovery hit-to-lead and lead optimization phases, but have high computational cost and the requirement of a structural analog with a known activity. Without a reference molecule requirement, absolute FEP (AFEP) has, in theory, better accuracy for hit ID, but in practice, the slow throughput is not compatible with VS, where fast docking and unreliable scoring functions are still the standard. Here, we present an integrated workflow to virtually screen large and diverse chemical libraries efficiently, combining active learning with a physics-based scoring function based on a fast absolute free energy perturbation method. We validated the performance of the approach in the ranking of structurally related ligands, virtual screening hit rate enrichment, and active learning chemical space exploration; disclosing the largest reported collection of free energy simulations to date.

Learning a reactive potential for silica-water through uncertainty attribution.

The reactivity of silicates in aqueous solution is relevant to various chemistries ranging from silicate minerals in geology, to the C-S-H phase in cement, nanoporous zeolite catalysts, or highly porous precipitated silica. While simulations of chemical reactions can provide insight at the molecular level, balancing accuracy and scale in reactive simulations in the condensed phase is a challenge. Here, we demonstrate how a machine-learning reactive interatomic potential trained on PaiNN architecture can accurately capture silicate-water reactivity. The model was trained on a dataset comprising 400,000 energies and forces of molecular clusters at the ωB97X-D3/def2-TZVP level. To ensure the robustness of the model, we introduce a general active learning strategy based on the attribution of the model uncertainty, that automatically isolates uncertain regions of bulk simulations to be calculated as small-sized clusters. The potential reproduces static and dynamic properties of liquid water and solid crystalline silicates, despite having been trained exclusively on cluster data. Furthermore, we utilize enhanced sampling simulations to recover the self-ionization reactivity of water accurately, and the acidity of silicate oligomers, and lastly study the silicate dimerization reaction in a water solution at neutral conditions and find that the reaction occurs through a flanking mechanism.

A Mechanism-Based Reaction-Diffusion Model for Accelerated Discovery of Thermoset Resins Frontally Polymerized by Olefin Metathesis.

Frontal ring-opening metathesis polymerization (FROMP) involves a self-perpetuating exothermic reaction, which enables the rapid and energy-efficient manufacturing of thermoset polymers and composites. Current state-of-the-art reaction-diffusion FROMP models rely on a phenomenological description of the olefin metathesis kinetics, limiting their ability to model the governing thermo-chemical FROMP processes. Furthermore, the existing models are unable to predict the variations in FROMP kinetics with changes in the resin composition and as a result are of limited utility toward accelerated discovery of new resin formulations. In this work, we formulate a chemically meaningful model grounded in the established mechanism of ring-opening metathesis polymerization (ROMP). Our study aims to validate the hypothesis that the ROMP mechanism, applicable to monomer-initiator solutions below 100 °C, remains valid under the nonideal conditions encountered in FROMP, including ambient to >200 °C temperatures, sharp temperature gradients, and neat monomer environments. Through extensive simulations, we demonstrate that our mechanism-based model accurately predicts the FROMP behavior across various resin compositions, including polymerization front velocities and thermal characteristics (e.g., Tmax). Additionally, we introduce a semi-inverse workflow that predicts FROMP behavior from a single experimental data point. Notably, the physiochemical parameters utilized in our model can be obtained through DFT calculations and minimal experiments, highlighting the model's potential for rapid screening of new FROMP chemistries in pursuit of thermoset polymers with superior thermo-chemo-mechanical properties.

Learning Collective Variables with Synthetic Data Augmentation through Physics-Inspired Geodesic Interpolation.

In molecular dynamics simulations, rare events, such as protein folding, are typically studied using enhanced sampling techniques, most of which are based on the definition of a collective variable (CV) along which acceleration occurs. Obtaining an expressive CV is crucial, but often hindered by the lack of information about the particular event, e.g., the transition from unfolded to folded conformation. We propose a simulation-free data augmentation strategy using physics-inspired metrics to generate geodesic interpolations resembling protein folding transitions, thereby improving sampling efficiency without true transition state samples. This new data can be used to improve the accuracy of classifier-based methods. Alternatively, a regression-based learning scheme for CV models can be adopted by leveraging the interpolation progress parameter.

Low-Dispersity Polymers via Free Radical Alternating Copolymerization: Effects of Charge-Transfer-Complexes.

Alternating copolymers are crucial for diverse applications. While dispersity (Ɖ, also known as molecular weight distribution, MWD) influences the properties of polymers, achieving low dispersities in alternating copolymers poses a notable challenge via free radical polymerizations (FRPs). In this work, we demonstrated an unexpected discovery that dispersities are affected by the participation of charge transfer complexes (CTCs) formed between monomer pairs during free radical alternating copolymerization, which have inspired the successful synthesis of various alternating copolymers with low dispersities (>30 examples, Ɖ = 1.13-1.39) under visible-light irradiation. The synthetic method is compatible with binary, ternary and quaternary alternating copolymerizations and is expandable for both fluorinated and non-fluorinated monomer pairs. DFT calculations combined with model experiments indicated that CTC-absent reaction exhibits higher propagation rates and affords fewer radical terminations, which could contribute to low dispersities. Based on the integration of Monte Carlo simulation and Bayesian optimization, we established the relationship map between FRP parameter space and dispersity, further suggested the correlation between low dispersities and higher propagation rates. Our research sheds light on dispersity control via FRPs and creates a novel platform to investigate polymer dispersity through machine learning.

Reactive Processing of Furan-Based Monomers via Frontal Ring-Opening Metathesis Polymerization for High Performance Materials.

Frontal ring-opening metathesis polymerization (FROMP) presents an energy-efficient approach to produce high-performance polymers, typically utilizing norbornene derivatives from Diels-Alder reactions. This study broadens the monomer repertoire for FROMP, incorporating the cycloaddition product of biosourced furan compounds and benzyne, namely 1,4-dihydro-1,4-epoxynaphthalene (HEN) derivatives. A computational screening of Diels-Alder products is conducted, selecting products with resistance to retro-Diels-Alder but also sufficient ring strain to facilitate FROMP. The experiments reveal that varying substituents both modulate the FROMP kinetics and enable the creation of thermoplastic materials characterized by different thermomechanical properties. Moreover, HEN-based crosslinkers are designed to enhance the resulting thermomechanical properties at high temperatures (>200 °C). The versatility of such materials is demonstrated through direct ink writing (DIW) to rapidly produce 3D structures without the need for printed supports. This research significantly extends the range of monomers suitable for FROMP, furthering efficient production of high-performance polymeric materials.

Atomic order of rare earth ions in a complex oxide: a path to magnetotaxial anisotropy.

Complex oxides offer rich magnetic and electronic behavior intimately tied to the composition and arrangement of cations within the structure. Rare earth iron garnet films exhibit an anisotropy along the growth direction which has long been theorized to originate from the ordering of different cations on the same crystallographic site. Here, we directly demonstrate the three-dimensional ordering of rare earth ions in pulsed laser deposited (EuxTm1-x)3Fe5O12 garnet thin films using both atomically-resolved elemental mapping to visualize cation ordering and X-ray diffraction to detect the resulting order superlattice reflection. We quantify the resulting ordering-induced 'magnetotaxial' anisotropy as a function of Eu:Tm ratio using transport measurements, showing an overwhelmingly dominant contribution from magnetotaxial anisotropy that reaches 30 kJ m-3 for garnets with x = 0.5. Control of cation ordering on inequivalent sites provides a strategy to control matter on the atomic level and to engineer the magnetic properties of complex oxides.

ZeoSyn: A Comprehensive Zeolite Synthesis Dataset Enabling Machine-Learning Rationalization of Hydrothermal Parameters.

Zeolites, nanoporous aluminosilicates with well-defined porous structures, are versatile materials with applications in catalysis, gas separation, and ion exchange. Hydrothermal synthesis is widely used for zeolite production, offering control over composition, crystallinity, and pore size. However, the intricate interplay of synthesis parameters necessitates a comprehensive understanding of synthesis-structure relationships to optimize the synthesis process. Hitherto, public zeolite synthesis databases only contain a subset of parameters and are small in scale, comprising up to a few thousand synthesis routes. We present ZeoSyn, a dataset of 23,961 zeolite hydrothermal synthesis routes, encompassing 233 zeolite topologies and 921 organic structure-directing agents (OSDAs). Each synthesis route comprises comprehensive synthesis parameters: 1) gel composition, 2) reaction conditions, 3) OSDAs, and 4) zeolite products. Using ZeoSyn, we develop a machine learning classifier to predict the resultant zeolite given a synthesis route with >70% accuracy. We employ SHapley Additive exPlanations (SHAP) to uncover key synthesis parameters for >200 zeolite frameworks. We introduce an aggregation approach to extend SHAP to all building units. We demonstrate applications of this approach to phase-selective and intergrowth synthesis. This comprehensive analysis illuminates the synthesis parameters pivotal in driving zeolite crystallization, offering the potential to guide the synthesis of desired zeolites. The dataset is available at https://github.com/eltonpan/zeosyn_dataset.

Design of Cytotoxic T Cell Epitopes by Machine Learning of Human Degrons.

Antigen processing is critical for therapeutic vaccines to generate epitopes for priming cytotoxic T cell responses against cancer and pathogens, but insufficient processing often limits the quantity of epitopes released. We address this challenge using machine learning to ascribe a proteasomal degradation score to epitope sequences. Epitopes with varying scores were translocated into cells using nontoxic anthrax proteins. Epitopes with a low score show pronounced immunogenicity due to antigen processing, but epitopes with a high score show limited immunogenicity. This work sheds light on the sequence-activity relationships between proteasomal degradation and epitope immunogenicity. We anticipate that future efforts to incorporate proteasomal degradation signals into vaccine designs will lead to enhanced cytotoxic T cell priming by these vaccines in clinical settings.

One-Pot Synthesis of CHA/ERI-Type Zeolite Intergrowth from a Single Multiselective Organic Structure-Directing Agent.

We report the one-pot synthesis of a chabazite (CHA)/erionite (ERI)-type zeolite intergrowth structure characterized by adjustable extents of intergrowth enrichment and Si/Al molar ratios. This method utilizes readily synthesizable 6-azaspiro[5.6]dodecan-6-ium as the exclusive organic structure-directing agent (OSDA) within a potassium-dominant environment. High-throughput simulations were used to accurately determine the templating energy and molecular shape, facilitating the selection of an optimally biselective OSDA from among thousands of prospective candidates. The coexistence of the crystal phases, forming a distinct structure comprising disk-like CHA regions bridged by ERI-rich pillars, was corroborated via rigorous powder X-ray diffraction and integrated differential-phase contrast scanning transmission electron microscopy (iDPC S/TEM) analyses. iDPC S/TEM imaging further revealed the presence of single offretite layers dispersed within the ERI phase. The ratio of crystal phases between CHA and ERI in this type of intergrowth could be varied systematically by changing both the OSDA/Si and K/Si ratios. Two intergrown zeolite samples with different Si/Al molar ratios were tested for the selective catalytic reduction (SCR) of NOx with NH3, showing competitive catalytic performance and hydrothermal stability compared to that of the industry-standard commercial NH3-SCR catalyst, Cu-SSZ-13, prevalent in automotive applications. Collectively, this work underscores the potential of our approach for the synthesis and optimization of adjustable intergrown zeolite structures, offering competitive alternatives for key industrial processes.

Autonomous, multiproperty-driven molecular discovery: From predictions to measurements and back.

A closed-loop, autonomous molecular discovery platform driven by integrated machine learning tools was developed to accelerate the design of molecules with desired properties. We demonstrated two case studies on dye-like molecules, targeting absorption wavelength, lipophilicity, and photooxidative stability. In the first study, the platform experimentally realized 294 unreported molecules across three automatic iterations of molecular design-make-test-analyze cycles while exploring the structure-function space of four rarely reported scaffolds. In each iteration, the property prediction models that guided exploration learned the structure-property space of diverse scaffold derivatives, which were realized with multistep syntheses and a variety of reactions. The second study exploited property models trained on the explored chemical space and previously reported molecules to discover nine top-performing molecules within a lightly explored structure-property space.

Effect of Framework Composition and NH3 on the Diffusion of Cu+ in Cu-CHA Catalysts Predicted by Machine-Learning Accelerated Molecular Dynamics.

Cu-exchanged zeolites rely on mobile solvated Cu+ cations for their catalytic activity, but the role of the framework composition in transport is not fully understood. Ab initio molecular dynamics simulations can provide quantitative atomistic insight but are too computationally expensive to explore large length and time scales or diverse compositions. We report a machine-learning interatomic potential that accurately reproduces ab initio results and effectively generalizes to allow multinanosecond simulations of large supercells and diverse chemical compositions. Biased and unbiased simulations of [Cu(NH3)2]+ mobility show that aluminum pairing in eight-membered rings accelerates local hopping and demonstrate that increased NH3 concentration enhances long-range diffusion. The probability of finding two [Cu(NH3)2]+ complexes in the same cage, which is key for SCR-NOx reaction, increases with Cu content and Al content but does not correlate with the long-range mobility of Cu+. Supporting experimental evidence was obtained from reactivity tests of Cu-CHA catalysts with a controlled chemical composition.

Examining graph neural networks for crystal structures: Limitations and opportunities for capturing periodicity.

Graph neural networks (GNNs) have recently been used to learn the representations of crystal structures through an end-to-end data-driven approach. However, a systematic top-down approach to evaluate and understand the limitations of GNNs in accurately capturing crystal structures has yet to be established. In this study, we introduce an approach using human-designed descriptors as a compendium of human knowledge to investigate the extent to which GNNs can comprehend crystal structures. Our findings reveal that current state-of-the-art GNNs fall short in accurately capturing the periodicity of crystal structures. We analyze this failure by exploring three aspects: local expressive power, long-range information processing, and readout function. To address these identified limitations, we propose a straightforward and general solution: the hybridization of descriptors with GNNs, which directly supplements the missing information to GNNs. The hybridization enhances the predictive accuracy of GNNs for specific material properties, most notably phonon internal energy and heat capacity, which heavily rely on the periodicity of materials.

A Model Ensemble Approach Enables Data-Driven Property Prediction for Chemically Deconstructable Thermosets in the Low-Data Regime.

Thermosets present sustainability challenges that could potentially be addressed through the design of deconstructable variants with tunable properties; however, the combinatorial space of possible thermoset molecular building blocks (e.g., monomers, cross-linkers, and additives) and manufacturing conditions is vast, and predictive knowledge for how combinations of these molecular components translate to bulk thermoset properties is lacking. Data science could overcome these problems, but computational methods are difficult to apply to multicomponent, amorphous, statistical copolymer materials for which little data exist. Here, leveraging a data set with 101 examples, we introduce a closed-loop experimental, machine learning (ML), and virtual screening strategy to enable predictions of the glass transition temperature (Tg) of polydicyclopentadiene (pDCPD) thermosets containing cleavable bifunctional silyl ether (BSE) comonomers and/or cross-linkers with varied compositions and loadings. Molecular features and formulation variables are used as model inputs, and uncertainty is quantified through model ensembling, which together with heavy regularization helps to avoid overfitting and ultimately achieves predictions within <15 °C for thermosets with compositionally diverse BSEs. This work offers a path to predicting the properties of thermosets based on their molecular building blocks, which may accelerate the discovery of promising plastics, rubbers, and composites with improved functionality and controlled deconstructability.

Design of Cytotoxic T Cell Epitopes by Machine Learning of Human Degrons.

Antigen processing is critical for producing epitope peptides that are presented by HLA molecules for T cell recognition. Therapeutic vaccines aim to harness these epitopes for priming cytotoxic T cell responses against cancer and pathogens, but insufficient processing often reduces vaccine efficacy through limiting the quantity of epitopes released. Here, we set out to improve antigen processing by harnessing protein degradation signals called degrons from the ubiquitin-proteasome system. We used machine learning to generate a computational model that ascribes a proteasomal degradation score between 0 and 100. Epitope peptides with varying degron activities were synthesized and translocated into cells using nontoxic anthrax proteins: protective antigen (PA) and the N-terminus of lethal factor (LFN). Immunogenicity studies revealed epitope sequences with a low score (<25) show pronounced T-cell activation but epitope sequences with a higher score (>75) provide limited activation. This work sheds light on the sequence-activity relationships between proteasomal degradation and epitope immunogenicity, through conserving the epitope region but varying the flanking sequence. We anticipate that future efforts to incorporate proteasomal degradation signals into vaccine designs will lead to enhanced cytotoxic T cell priming by vaccine therapeutics in clinical settings.

Mapping the Space of Photoswitchable Ligands and Photodruggable Proteins with Computational Modeling.

Light-activated drugs are a promising way to localize biological activity and minimize side effects. However, their development is complicated by the numerous photophysical and biological properties that must be simultaneously optimized. To accelerate the design of photoactive drugs, we describe a procedure that combines ligand-protein docking with chemical property prediction based on machine learning (ML). We apply this procedure to 58 proteins and 9000 photo-drug candidates based on azobenzene cis-trans isomerism. We find that most proteins display a preference for trans isomers over cis and that the binding affinities of nominally active/inactive pairs are in fact highly correlated. These findings have significant value for photopharmacology research, and reinforce the need for virtual screening to identify compounds with rare desirable properties. Further, we combine our procedure with quantum chemical validation to identify promising candidates for the photoactive inhibition of PARP1, an enzyme that is over-expressed in cancer cells. The top compounds are predicted to have long-lived active forms, differential bioactivity, and absorption in the near-infrared therapeutic window.

Entropy and Energy Profiles of Chemical Reactions.

The description of chemical processes at the molecular level is often facilitated by the use of reaction coordinates or collective variables (CVs). The CV measures the progress of the reaction and allows the construction of profiles that track how specific properties evolve as the reaction progresses. Whereas CVs are routinely used, especially alongside enhanced sampling techniques, the links among reaction profiles, thermodynamic state functions, and reaction rate constants are not rigorously exploited. Here, we report a unified treatment of such reaction profiles. Tractable expressions are derived for the free-energy, internal-energy, and entropy profiles as functions of only the CV. We demonstrate the ability of this treatment to extract quantitative insight from the entropy and internal-energy profiles of various real-world physicochemical processes, including intramolecular organic reactions, ionic transport in superionic electrolytes, and molecular transport in nanoporous materials.

Lamellar Ionenes with Highly Dissociative, Anionic Channels Provide Lower Barriers for Cation Transport.

Solid polymer electrolytes have the potential to enable safer and more energy-dense batteries; however, a deeper understanding of their ion conduction mechanisms, and how they can be optimized by molecular design, is needed to realize this goal. Here, we investigate the impact of anion dissociation energy on ion conduction in solid polymer electrolytes via a novel class of ionenes prepared using acyclic diene metathesis (ADMET) polymerization of highly dissociative, liquid crystalline fluorinated aryl sulfonimide-tagged ("FAST") anion monomers. These ionenes with various cations (Li+, Na+, K+, and Cs+) form well-ordered lamellae that are thermally stable up to 180 °C and feature domain spacings that correlate with cation size, providing channels lined with dissociative FAST anions. Electrochemical impedance spectroscopy (EIS) and differential scanning calorimetry (DSC) experiments, along with nudged elastic band (NEB) calculations, suggest that cation motion in these materials operates via an ion-hopping mechanism. The activation energy for Li+ conduction is 59 kJ/mol, which is among the lowest for systems that are proposed to operate via an ion conduction mechanism that is decoupled from polymer segmental motion. Moreover, the addition of a cation-coordinating solvent to these materials led to a >1000-fold increase in ionic conductivity without detectable disruption of the lamellar structure, suggesting selective solvation of the lamellar ion channels. This work demonstrates that molecular design can facilitate controlled formation of dissociative anionic channels that translate to significant enhancements in ion conduction in solid polymer electrolytes.

Approaching enzymatic catalysis with zeolites or how to select one reaction mechanism competing with others.

Approaching the level of molecular recognition of enzymes with solid catalysts is a challenging goal, achieved in this work for the competing transalkylation and disproportionation of diethylbenzene catalyzed by acid zeolites. The key diaryl intermediates for the two competing reactions only differ in the number of ethyl substituents in the aromatic rings, and therefore finding a selective zeolite able to recognize this subtle difference requires an accurate balance of the stabilization of reaction intermediates and transition states inside the zeolite microporous voids. In this work we present a computational methodology that, by combining a fast high-throughput screeening of all zeolite structures able to stabilize the key intermediates with a more computationally demanding mechanistic study only on the most promising candidates, guides the selection of the zeolite structures to be synthesized. The methodology presented is validated experimentally and allows to go beyond the conventional criteria of zeolite shape-selectivity.

Learning pair potentials using differentiable simulations.

Learning pair interactions from experimental or simulation data is of great interest for molecular simulations. We propose a general stochastic method for learning pair interactions from data using differentiable simulations (DiffSim). DiffSim defines a loss function based on structural observables, such as the radial distribution function, through molecular dynamics (MD) simulations. The interaction potentials are then learned directly by stochastic gradient descent, using backpropagation to calculate the gradient of the structural loss metric with respect to the interaction potential through the MD simulation. This gradient-based method is flexible and can be configured to simulate and optimize multiple systems simultaneously. For example, it is possible to simultaneously learn potentials for different temperatures or for different compositions. We demonstrate the approach by recovering simple pair potentials, such as Lennard-Jones systems, from radial distribution functions. We find that DiffSim can be used to probe a wider functional space of pair potentials compared with traditional methods like iterative Boltzmann inversion. We show that our methods can be used to simultaneously fit potentials for simulations at different compositions and temperatures to improve the transferability of the learned potentials.