An official website of the United States government
Ferrocene derivatives have a wide range of applications, including as ligands in asymmetric catalysis, due to their chemical stability, rigid backbone, steric bulk, and ability to encode stereochemical information via planar chirality. Unfortunately, few of the available molecular mechanics force fields incorporate parameters for the accurate study of this important building block. Here, we present a MM3* force field for ferrocenyl ligands, which was generated using the quantum-guided molecular mechanics (Q2MM) method. Detailed validation by comparison to DFT calculations and crystal structures demonstrates the accuracy of the parameters and uncovers the physical origin of deviations through excess energy analysis. Combining the ferrocene force field with a force field for Pd–allyl complexes and comparing the crystal structures shows the compatibility with previously developed MM3* force fields. Finally, the ferrocene force field was combined with a previously published transition-state force field to predict the stereochemical outcomes of the aminations of Pd–allyl complexes with different amines and different chiral ferrocenyl ligands, with an R2 of ∼0.91 over 10 examples.
more »
« less
Transition State Force Field for the Asymmetric Redox-Relay Heck Reaction
A transition state force field (TSFF) was developed using the quantum-guided molecular mechanics (Q2MM) method to describe the stereodetermining migratory insertion step of the enantioselective redox-relay Heck reaction for a range of multisubstituted alkenes. We show that the TSFF is highly predictive through an external validation of the TSFF against 151 experimentally determined stereoselectivities resulting in an R2 of 0.89 and MUE of 1.8 kJ/mol. In addition, limitations in the underlying force field were identified by comparison of the TSFF results to DFT level calculations. A novel application of the TSFF was demonstrated for 31 cases where the enantiomer predicted by the TSFF differed from the originally published values. Experimental determination of the absolute configuration demonstrated that the computational predictions were accurate, suggesting that TSFFs can be used for the rapid prediction of the absolute stereochemistry for a class of reactions. Finally, a virtual ligand screen was conducted utilizing both the TSFF and a simple molecular correlation method. Both methods were similarly predictive, but the TSFF was able to show greater utility through transferability, speed, and interpretability.
more »
« less
- Award ID(s):
- 1855908
- PAR ID:
- 10179623
- Date Published:
- Journal Name:
- Journal of the American Chemical Society
- Volume:
- 142
- Issue:
- 21
- ISSN:
- 0002-7863
- Page Range / eLocation ID:
- 9700-9707
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Molecular dynamics at the atomistic scale is increasingly being used to predict material properties and speed up the materials design and development process. However, the accuracy of molecular dynamics predictions is sensitively dependent on the force fields. In the traditional force field calibration process, a specific property, predicted by the model, is compared with the experimental observation and the force field parameters are adjusted to minimize the difference. This leads to the issue that the calibrated force fields are not generic and robust enough to predict different properties. Here, a new calibration method based on multi-objective Bayesian optimization is developed to speed up the development of molecular dynamics force fields that are capable of predicting multiple properties accurately. This is achieved by reducing the number of simulation runs to generate the Pareto front with an efficient sequential sampling strategy. The methodology is demonstrated by generating a new coarse-grained force field for polycaprolactone, where the force field can predict mechanical properties and water diffusion in the polymer.more » « less
-
Gaussian process (GP) emulator has been used as a surrogate model for predicting force field and molecular potential, to overcome the computational bottleneck of ab initio molecular dynamics simulation. Integrating both atomic force and energy in predictions was found to be more accurate than using energy alone, yet it requires O(( NM) 3 ) computational operations for computing the likelihood function and making predictions, where N is the number of atoms and M is the number of simulated configurations in the training sample due to the inversion of a large covariance matrix. The high computational cost limits its applications to the simulation of small molecules. The computational challenge of using both gradient information and function values in GPs was recently noticed in machine learning communities, whereas conventional approximation methods may not work well. Here, we introduce a new approach, the atomized force field model, that integrates both force and energy in the emulator with many fewer computational operations. The drastic reduction in computation is achieved by utilizing the naturally sparse covariance structure that satisfies the constraints of the energy conservation and permutation symmetry of atoms. The efficient machine learning algorithm extends the limits of its applications on larger molecules under the same computational budget, with nearly no loss of predictive accuracy. Furthermore, our approach contains an uncertainty assessment of predictions of atomic forces and energies, useful for developing a sequential design over the chemical input space.more » « less
-
null (Ed.)We developed a new coarse-grained (CG) molecular dynamics force field for polyacrylamide (PAM) polymer based on fitting to the quantum mechanics (QM) equation of state (EOS). In this method, all nonbond interactions between representative beads are parameterized using a series of QM-EOS, which significantly improves the accuracy in comparison to common CG methods derived from atomistic molecular dynamics. This CG force-field has both higher accuracy and improved computational efficiency with respect to the OPLS atomistic force field. The nonbond components of the EOS were obtained from cold-compression curves on PAM crystals with rigid chains, while the covalent terms that contribute to the EOS were obtained using relaxed chains. For describing PAM gels we developed water–PAM interaction parameters using the same method. We demonstrate that the new CG-PAM force field reproduces the EOS of PAM crystals, isolated PAM chains, and water–PAM systems, while successfully predicting such experimental quantities as density, specific heat capacity, thermal conductivity and melting point.more » « less
-
Despite the enormous potential for the use of stereospecific cross-coupling reactions to rationally manipulate the three-dimensional structure of organic molecules, the factors that control the transfer of stereochemistry in these reactions remain poorly understood. Here we report a mechanistic and synthetic investigation into the use of enantioenriched alkylboron nucleophiles in stereospecific Pd-catalyzed Suzuki cross-coupling reactions. By developing a suite of molecular descriptors of phosphine ligands, we could apply predictive statistical models to select or design distinct ligands that respectively promoted stereoinvertive and stereoretentive cross-coupling reactions. Stereodefined branched structures were thereby accessed through the predictable manipulation of absolute stereochemistry, and a general model for the mechanism of alkylboron transmetallation was proposed.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/jacs.0c01979
