Metal-mediated cross-coupling reactions offer organic chemists a wide array of stereo- and chemically-selective reactions with broad applications in fine chemical and pharmaceutical synthesis.1 Current batch-based synthesis methods are beginning to be replaced with flow chemistry strategies to take advantage of the improved consistency and process control methods offered by continuous flow systems.2,3 Most cross-coupling chemistries still encounter several issues in flow using homogeneous catalysis, including expensive catalyst recovery and air sensitivity due to the chemical nature of the catalyst ligands.1 To mitigate some of these issues, a ligand-free heterogeneous catalysis reaction was developed using palladium (Pd) loaded into a polymeric network of a silicone elastomer, poly(hydromethylsiloxane) (PHMS), that is not air sensitive and can be used with mild reaction solvents (ethanol and water).4
In this work we present a novel method of producing soft catalytic microparticles using a multiphase flow-focusing microreactor and demonstrate their application for continuous Suzuki-Miyaura cross-coupling reactions. The catalytic microparticles are produced in a coaxial glass capillary-based 3D flow-focusing microreactor. The microreactor consists of two precursors, a cross-linking catalyst in toluene and a mixture of the PHMS polymer and a divinyl cross-linker. The dispersed phase containing the polymer, cross-linker, and cross-linking catalyst is continuously mixed and then formed into microdroplets by the continuous phase of water and surfactant (sodium dodecyl sulfate) introduced in a counter-flow configuration. Elastomeric microdroplets with a diameter ranging between 50 to 300 micron are produced at 25 to 250 Hz with a size polydispersity less than 3% in single stream production. The physicochemical properties of the elastomeric microparticles such as particle swelling/softness can be tuned using the ratio of cross-linker to polymer as well as the ratio of polymer mixture to solvent during the particle formation. Swelling in toluene can be tuned up to 400% of the initial particle volume by reducing the concentration of cross-linker in the mixture and increasing the ratio of polymer to solvent during production.5
After the particles are produced and collected, they are transferred into toluene containing palladium acetate, allowing the particles to incorporate the palladium into the polymer network and then reduce the palladium to Pd0 with the Si-H functionality present on the PHMS backbones. After the reduction, the Pd-loaded particles can be washed and dried for storage or switched into an ethanol/water solution for loading into a micro-packed bed reactor (µ-PBR) for continuous organic synthesis. The in-situ reduction of Pd within the PHMS microparticles was confirmed using energy dispersive X-ray spectroscopy (EDS), X-ray photoelectron spectroscopy (XPS) and focused ion beam-SEM, and TEM techniques.
In the next step, we used the developed µ-PBR to conduct continuous organic synthesis of 4-phenyltoluene by Suzuki-Miyaura cross-coupling of 4-iodotoluene and phenylboronic acid using potassium carbonate as the base. Catalyst leaching was determined to only occur at sub ppm concentrations even at high solvent flow rates after 24 h of continuous run using inductively coupled plasma mass spectrometry (ICP-MS).
The developed µ-PBR using the elastomeric microparticles is an important initial step towards the development of highly-efficient and green continuous manufacturing technologies in the pharma industry. In addition, the developed elastomeric microparticle synthesis technique can be utilized for the development of a library of other chemically cross-linkable polymer/cross-linker pairs for applications in organic synthesis, targeted drug delivery, cell encapsulation, or biomedical imaging.
References
1. Ruiz-Castillo P, Buchwald SL. Applications of Palladium-Catalyzed C-N Cross-Coupling Reactions. Chem Rev. 2016;116(19):12564-12649.
2. Adamo A, Beingessner RL, Behnam M, et al. On-demand continuous-flow production of pharmaceuticals in a compact, reconfigurable system. Science. 2016;352(6281):61 LP-67.
3. Jensen KF. Flow Chemistry — Microreaction Technology Comes of Age. 2017;63(3).
4. Stibingerova I, Voltrova S, Kocova S, Lindale M, Srogl J. Modular Approach to Heterogenous Catalysis. Manipulation of Cross-Coupling Catalyst Activity. Org Lett. 2016;18(2):312-315.
5. Bennett JA, Kristof AJ, Vasudevan V, Genzer J, Srogl J, Abolhasani M. Microfluidic synthesis of elastomeric microparticles: A case study in catalysis of palladium-mediated cross-coupling. AIChE J. 2018;0(0):1-10.
more »
« less
Fast predictions of liquid-phase acid-catalyzed reaction rates using molecular dynamics simulations and convolutional neural networks
The rates of liquid-phase, acid-catalyzed reactions relevant to the upgrading of biomass into high-value chemicals are highly sensitive to solvent composition and identifying suitable solvent mixtures is theoretically and experimentally challenging. We show that the complex atomistic configurations of reactant–solvent environments generated by classical molecular dynamics simulations can be exploited by 3D convolutional neural networks to enable accurate predictions of Brønsted acid-catalyzed reaction rates for model biomass compounds. We develop a 3D convolutional neural network, which we call SolventNet, and train it to predict acid-catalyzed reaction rates using experimental reaction data and corresponding molecular dynamics simulation data for seven biomass-derived oxygenates in water–cosolvent mixtures. We show that SolventNet can predict reaction rates for additional reactants and solvent systems an order of magnitude faster than prior simulation methods. This combination of machine learning with molecular dynamics enables the rapid, high-throughput screening of solvent systems and identification of improved biomass conversion conditions.
more »
« less
- Award ID(s):
- 1720415
- NSF-PAR ID:
- 10208670
- Date Published:
- Journal Name:
- Chemical Science
- Volume:
- 11
- Issue:
- 46
- ISSN:
- 2041-6520
- Page Range / eLocation ID:
- 12464 to 12476
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Photochemical reactions are widely used by academic and industrial researchers to construct complex molecular architectures via mechanisms that often require harsh reaction conditions. Photodynamics simulations provide time-resolved snapshots of molecular excited-state structures required to understand and predict reactivities and chemoselectivities. Molecular excited-states are often nearly degenerate and require computationally intensive multiconfigurational quantum mechanical methods, especially at conical intersections. Non-adiabatic molecular dynamics require thousands of these computations per trajectory, which limits simulations to ∼1 picosecond for most organic photochemical reactions. Westermayr et al. recently introduced a neural-network-based method to accelerate the predictions of electronic properties and pushed the simulation limit to 1 ns for the model system, methylenimmonium cation (CH 2 NH 2 + ). We have adapted this methodology to develop the Python-based, Python Rapid Artificial Intelligence Ab Initio Molecular Dynamics (PyRAI 2 MD) software for the cis – trans isomerization of trans -hexafluoro-2-butene and the 4π-electrocyclic ring-closing of a norbornyl hexacyclodiene. We performed a 10 ns simulation for trans -hexafluoro-2-butene in just 2 days. The same simulation would take approximately 58 years with traditional multiconfigurational photodynamics simulations. We generated training data by combining Wigner sampling, geometrical interpolations, and short-time quantum chemical trajectories to adaptively sample sparse data regions along reaction coordinates. The final data set of the cis – trans isomerization and the 4π-electrocyclic ring-closing model has 6207 and 6267 data points, respectively. The training errors in energy using feedforward neural networks achieved chemical accuracy (0.023–0.032 eV). The neural network photodynamics simulations of trans -hexafluoro-2-butene agree with the quantum chemical calculations showing the formation of the cis -product and reactive carbene intermediate. The neural network trajectories of the norbornyl cyclohexadiene corroborate the low-yielding syn -product, which was absent in the quantum chemical trajectories, and revealed subsequent thermal reactions in 1 ns.more » « less
-
The ability to predict and understand complex molecular motions occurring over diverse timescales ranging from picoseconds to seconds and even hours in biological systems remains one of the largest challenges to chemical theory. Markov state models (MSMs), which provide a memoryless description of the transitions between different states of a biochemical system, have provided numerous important physically transparent insights into biological function. However, constructing these models often necessitates performing extremely long molecular simulations to converge the rates. Here, we show that by incorporating memory via the time-convolutionless generalized master equation (TCL-GME) one can build a theoretically transparent and physically intuitive memory-enriched model of biochemical processes with up to a three order of magnitude reduction in the simulation data required while also providing a higher temporal resolution. We derive the conditions under which the TCL-GME provides a more efficient means to capture slow dynamics than MSMs and rigorously prove when the two provide equally valid and efficient descriptions of the slow configurational dynamics. We further introduce a simple averaging procedure that enables our TCL-GME approach to quickly converge and accurately predict long-time dynamics even when parameterized with noisy reference data arising from short trajectories. We illustrate the advantages of the TCL-GME using alanine dipeptide, the human argonaute complex, and FiP35 WW domain.more » « less
-
Efficient transfer of halogen atoms is essential for controlling the growth of polymers in atom transfer radical polymerization (ATRP). The nature of halogens may influence the efficiency of the halogen atom transfer during the activation and deactivation processes. The effect of halogens can be associated with the C–X bond dissociation energy and the affinity of the halogens/halides to the transition metal catalyst. In this paper, we study the effect of halogens (Br vs. Cl) and reaction media in iron-catalyzed ATRP in the presence of halide anions as ligands. In Br-based initiating systems, polymerization of methacrylate monomers was well-controlled whereas Cl-based initiating systems provided limited control over the polymerization. The high affinity of the Cl atom to the iron catalyst renders it less efficient for fast deactivation of growing chains, resulting in polymers with molecular weights higher than predetermined by Δ[M]/[RX] o and with high dispersities. Conversely, Br can be exchanged with higher efficiency and hence provided good control over polymerization. Decreasing the polarity of the reaction medium improved the polymerization control. Polymerizations using ppm levels of the iron catalyst in acetonitrile (a more polar solvent) yielded polymers with larger dispersity values due to the slow rate of deactivation as opposed to the less polar solvent anisole, which afforded well-controlled polymers with dispersity <1.2.more » « less
-
more » « less
Abstract Deep learning models are seeing increased use as methods to predict mutational effects or allowed mutations in proteins. The models commonly used for these purposes include large language models (LLMs) and 3D Convolutional Neural Networks (CNNs). These two model types have very different architectures and are commonly trained on different representations of proteins. LLMs make use of the transformer architecture and are trained purely on protein sequences whereas 3D CNNs are trained on voxelized representations of local protein structure. While comparable overall prediction accuracies have been reported for both types of models, it is not known to what extent these models make comparable specific predictions and/or generalize protein biochemistry in similar ways. Here, we perform a systematic comparison of two LLMs and two structure-based models (CNNs) and show that the different model types have distinct strengths and weaknesses. The overall prediction accuracies are largely uncorrelated between the sequence- and structure-based models. Overall, the two structure-based models are better at predicting buried aliphatic and hydrophobic residues whereas the two LLMs are better at predicting solvent-exposed polar and charged amino acids. Finally, we find that a combined model that takes the individual model predictions as input can leverage these individual model strengths and results in significantly improved overall prediction accuracy.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/d0sc03261a
