101
|
Khatun M, Paul S, Roy S, Dey S, Anoop A. Performance of Density Functionals and Semiempirical 3c Methods for Small Gold-Thiolate Clusters. J Phys Chem A 2023; 127:2242-2257. [PMID: 36877153 DOI: 10.1021/acs.jpca.2c07561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
Abstract
In light of the recent surge in computational studies of gold thiolate clusters, we present a comparison of popular density functionals (DFAs) and three-part corrected methods (3c-methods) on their performance by taking a data set named as AuSR18 consisting of 18 isomers of Aun(SCH3)m (m ≤ n = 1-3). We have compared the efficiency and accuracy of the DFAs and 3c-methods in geometry optimization with RI-SCS-MP2 as the reference method. Similarly, the performance for accurate and efficient energy evaluation was compared with DLPNO-CCSD(T) as the reference method. The lowest energy structure among the isomers of the largest stoichiometry from our data set, AuSR18, i.e., Au3(SCH3)3, is considered to evaluate the computational time for SCF and gradient evaluations. Alongside this, the numbers of optimization steps to locate the most stable minima of Au3(SCH3)3 are compared to assess the efficiency of the methods. A comparison of relevant bond lengths with the reference geometries was made to estimate the accuracy in geometry optimization. Some methods, such as LC-BLYP, ωB97M-D3BJ, M06-2X, and PBEh-3c, could not locate many of the minima found by most of the other methods; thus, the versatility in locating various minima is also an important criterion in choosing a method for the given project. To determine the accuracy of the methods, we compared the relative energies of the isomers in each stoichiometry and the interaction energy of the gold core with the ligands. The dependence of basis set size and relativistic effects on energies are also compared. The following are some of the highlights. TPSS has shown accuracy, while mPWPW shows comparable speed and accuracy. For the relative energies of the clusters, the hybrid range-separated DFAs are the best option. CAM-B3LYP excels, whereas B3LYP performs poorly. Overall, LC-BLYP is a balanced performer considering both the geometry and relative stability of the structures, but it lacks diversity. The 3c-methods, although fast, are less impressive in relative stability.
Collapse
Affiliation(s)
- Maya Khatun
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Sayan Paul
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Saikat Roy
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Subhasis Dey
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Anakuthil Anoop
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| |
Collapse
|
102
|
Billion A, Schorpp M, Feser R, Schmitt M, Eisele L, Scherer H, Sonoda T, Kawa H, Butschke B, Krossing I. The perfluoroadamantoxy aluminate as an ideal weakly coordinating anion? - synthesis and first applications. Dalton Trans 2023; 52:4355-4370. [PMID: 36924178 DOI: 10.1039/d3dt00199g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023]
Abstract
Weakly coordinating anions (WCAs) facilitate the stabilization and isolation of highly reactive and almost "naked" cations. Alkoxyaluminate-based WCAs such as [Al(OC(CF3)3)4]- ([pf]-) are widely used due to their synthetic accessibility and their high stability. However, small cations are still able to coordinate the oxygen atoms of the [pf]- anion or even to abstract an alkoxy ligand. The novel WCA [Al(OC10F15)4]- ([pfAd]-; OC10F15 = perfluoro-1-adamantoxy) is characterized by a very rigid core framework, thus indicating a higher stability towards fluoride-ion abstraction (DFT calculations) and providing hope to generate less disordered crystal structures. The [pfAd]- anion was generated by the reaction of the highly acidic alcohol perfluoro-1-adamantanol C10F15OH with LiAlH4 in o-DFB. Li[pfAd] could not be synthesized free of impurities (and still contains unreacted alcohol). Yet, starting from contaminated Li[pfAd], the very useful pure salts Ag[pfAd], [Ph3C][pfAd] and [H(OEt2)2][pfAd] could be synthesized. The salts were characterized by NMR spectroscopy, single-crystal X-ray diffraction and IR spectroscopy. Additionally, [NO][pfAd] could be synthesized containing alcohol impurities but nonetheless enabled the synthesis of the salt P9+[pfAd]-. The synthesis of Tl[pfAd] in a mixture of H2O/acetone/o-DFB demonstrated the water stability of the [pfAd]- anion.
Collapse
Affiliation(s)
- Andreas Billion
- Institut für Anorganische und Analytische Chemie and Freiburger Materialforschungszentrum FMF, Albert-Ludwigs-Universität Freiburg, Albertstr. 21, 79104 Freiburg, Germany.
| | - Marcel Schorpp
- Institut für Anorganische und Analytische Chemie and Freiburger Materialforschungszentrum FMF, Albert-Ludwigs-Universität Freiburg, Albertstr. 21, 79104 Freiburg, Germany.
| | - Rebecca Feser
- Institut für Anorganische und Analytische Chemie and Freiburger Materialforschungszentrum FMF, Albert-Ludwigs-Universität Freiburg, Albertstr. 21, 79104 Freiburg, Germany.
| | - Manuel Schmitt
- Institut für Anorganische und Analytische Chemie and Freiburger Materialforschungszentrum FMF, Albert-Ludwigs-Universität Freiburg, Albertstr. 21, 79104 Freiburg, Germany.
| | - Lea Eisele
- Institut für Anorganische und Analytische Chemie and Freiburger Materialforschungszentrum FMF, Albert-Ludwigs-Universität Freiburg, Albertstr. 21, 79104 Freiburg, Germany.
| | - Harald Scherer
- Institut für Anorganische und Analytische Chemie and Freiburger Materialforschungszentrum FMF, Albert-Ludwigs-Universität Freiburg, Albertstr. 21, 79104 Freiburg, Germany.
| | - Takaaki Sonoda
- Institute for Materials Chemistry and Engineering, Kyushu University, 6-1 Kasuga koen, Kasuga-shi, Fukuoka, 816-8580, Japan
| | - Hajimu Kawa
- Exfluor Research Corporation, 2350 Double Creek Drive, Round Rock, Texas 78664, USA
| | - Burkhard Butschke
- Institut für Anorganische und Analytische Chemie and Freiburger Materialforschungszentrum FMF, Albert-Ludwigs-Universität Freiburg, Albertstr. 21, 79104 Freiburg, Germany.
| | - Ingo Krossing
- Institut für Anorganische und Analytische Chemie and Freiburger Materialforschungszentrum FMF, Albert-Ludwigs-Universität Freiburg, Albertstr. 21, 79104 Freiburg, Germany.
| |
Collapse
|
103
|
Shaya J, Aloum L, Lu CS, Corridon PR, Aoudi A, Shunnar A, Alefishat E, Petroianu G. Theoretical Study of Hydroxylation of α- and β-Pinene by a Cytochrome P450 Monooxygenase Model. Int J Mol Sci 2023; 24:ijms24065150. [PMID: 36982225 PMCID: PMC10048887 DOI: 10.3390/ijms24065150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Revised: 02/09/2023] [Accepted: 02/13/2023] [Indexed: 03/30/2023] Open
Abstract
Previous studies on biocatalytic transformations of pinenes by cytochrome P450 (CYP) enzymes reveal the formation of different oxygenated products from a single substrate due to the multistate reactivity of CYP and the many reactive sites in the pinene scaffold. Up until now, the detailed mechanism of these biocatalytic transformations of pinenes have not been reported. Hereby, we report a systematic theoretical study of the plausible hydrogen abstraction and hydroxylation reactions of α- and β-pinenes by CYP using the density functional theory (DFT) method. All DFT calculations in this study were based on B3LYP/LAN computational methodology using the Gaussian09 software. We used the B3LYP functional with corrections for dispersive forces, BSSE, and anharmonicity to study the mechanism and thermodynamic properties of these reactions using a bare model (without CYP) and a pinene-CYP model. According to the potential energy surface and Boltzmann distribution for radical conformers, the major reaction products of CYP-catalyzed hydrogen abstraction from β-pinene are the doublet trans (53.4%) and doublet cis (46.1%) radical conformer at delta site. The formation of doublet cis/trans hydroxylated products released a total Gibbs free energy of about 48 kcal/mol. As for alpha pinene, the most stable radicals were trans-doublet (86.4%) and cis-doublet (13.6%) at epsilon sites, and their hydroxylation products released a total of ~50 kcal/mol Gibbs free energy. Our results highlight the likely C-H abstraction and oxygen rebounding sites accounting for the multi-state of CYP (doublet, quartet, and sextet spin states) and the formation of different conformers due to the presence of cis/trans allylic hydrogen in α-pinene and β-pinene molecules.
Collapse
Affiliation(s)
- Janah Shaya
- Department of Chemistry, College of Arts and Sciences, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
| | - Lujain Aloum
- Department of Pharmacology, College of Medicine and Health Sciences, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
| | - Chung-Shin Lu
- Department of General Education, National Taichung University of Science and Technology, Taichung 404, Taiwan, China
| | - Peter R Corridon
- Department of Immunology and Physiology, College of Medicine and Health Sciences, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
- Biomedical Engineering and Healthcare Engineering Innovation Center, Khalifa University, Abu Dhabi 127788, United Arab Emirates
- Center for Biotechnology, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
| | - Abdulrahman Aoudi
- Department of Chemistry, College of Arts and Sciences, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
| | - Abeer Shunnar
- Department of Chemistry, College of Arts and Sciences, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
| | - Eman Alefishat
- Department of Pharmacology, College of Medicine and Health Sciences, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
- Center for Biotechnology, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
- Department of Biopharmaceutics and Clinical Pharmacy, School of Pharmacy, The University of Jordan, Amman 11972, Jordan
| | - Georg Petroianu
- Department of Pharmacology, College of Medicine and Health Sciences, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
| |
Collapse
|
104
|
Villot C, Lao KU. Electronic structure theory on modeling short-range noncovalent interactions between amino acids. J Chem Phys 2023; 158:094301. [PMID: 36889981 DOI: 10.1063/5.0138032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the "gold standard" coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller-Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.
Collapse
Affiliation(s)
- Corentin Villot
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| |
Collapse
|
105
|
Iyer J, Brunsteiner M, Ray A, Davis A, Saraf I, Paudel A. Theoretical and Experimental Investigation of Autoxidation Propensity of Selected Drugs in Solution State. Mol Pharm 2023; 20:1768-1778. [PMID: 36757102 DOI: 10.1021/acs.molpharmaceut.2c00967] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
The C-H bond dissociation energy (BDE) of drug molecules is often used to estimate their relative propensities to undergo autoxidation. BDE calculations based on electronic structures provide a convenient means to estimate the risk for a given compound to degrade via autoxidation. This study aimed to verify the utility of calculated C-H BDEs of a range of drug molecules in predicting their autoxidation propensities, in the solution state. For the autoxidation study, 2,2'-azobis (2-methylpropionitrile) was employed as the solution state stressor, and the experimental reaction rate constants were determined employing ultraperformance liquid chromatographic (UPLC) methods. Reaction rates in the solution state were compared to the calculated C-H BDE values of the respective compounds. The results indicated a poor correlation for compounds in the solution state, and their relative stabilities could not be explained with C-H BDE. On the other hand, a favorable relationship was observed between the relative extent of ionization and the autoxidation rates of the selected compounds. In the solution state, factors such as the type and extent of drug ionization, degree and type of solvation have been shown to contribute to differences in reactivity. By applying the computational method involving the effect of H-atom abstraction and potential ionization sites in the molecule, the calculated C-H BDE should relate better to the experimental autoxidation rates.
Collapse
Affiliation(s)
- Jayant Iyer
- Research Center Pharmaceutical Engineering GmbH (RCPE), Graz 8010, Austria
| | | | - Andrew Ray
- Dew Modalities Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Macclesfield SK10 2NA, United Kingdom
| | - Adrian Davis
- Pfizer Worldwide Research and Development, Sandwich, Kent CT13 9NJ, United Kingdom
| | - Isha Saraf
- Research Center Pharmaceutical Engineering GmbH (RCPE), Graz 8010, Austria
| | - Amrit Paudel
- Research Center Pharmaceutical Engineering GmbH (RCPE), Graz 8010, Austria.,Graz University of Technology, Institute of Process and Particle Engineering, Graz 8010, Austria
| |
Collapse
|
106
|
Creutzberg J, Hedegård ED. A method to capture the large relativistic and solvent effects on the UV-vis spectra of photo-activated metal complexes. Phys Chem Chem Phys 2023; 25:6153-6163. [PMID: 36752122 DOI: 10.1039/d2cp04937f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
We have recently developed a method based on relativistic time-dependent density functional theory (TD-DFT) that allows the calculation of electronic spectra in solution (Creutzberg, Hedegård, J. Chem. Theory Comput.18, 2022, 3671). This method treats the solvent explicitly with a classical, polarizable embedding (PE) description. Furthermore, it employs the complex polarization propagator (CPP) formalism which allows calculations on complexes with a dense population of electronic states (such complexes are known to be problematic for conventional TD-DFT). Here, we employ this method to investigate both the dynamic and electronic effects of the solvent for the excited electronic states of trans-trans-trans-[Pt(N3)2(OH)2(NH3)2] in aqueous solution. This complex decomposes into species harmful to cancer cells under light irradiation. Thus, understanding its photo-physical properties may lead to a more efficient method to battle cancer. We quantify the effect of the underlying structure and dynamics by classical molecular mechanics simulations, refined with a subsequent DFT or semi-empirical optimization on a cluster. Moreover, we quantify the effect of employing different methods to set up the solvated system, e.g., how sensitive the results are to the method used for the refinement, and how large a solvent shell that is required. The electronic solvent effect is always included through a PE potential.
Collapse
Affiliation(s)
- Joel Creutzberg
- Division of Theoretical Chemistry, Lund University, Lund, Sweden.
| | - Erik Donovan Hedegård
- Division of Theoretical Chemistry, Lund University, Lund, Sweden. .,Department of Physics, Chemistry and Pharmacy, Campusvej 55, 5230 Odense, Denmark.
| |
Collapse
|
107
|
Weser O, Hein-Janke B, Mata RA. Automated handling of complex chemical structures in Z-matrix coordinates-The chemcoord library. J Comput Chem 2023; 44:710-726. [PMID: 36541725 DOI: 10.1002/jcc.27029] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 10/05/2022] [Accepted: 10/09/2022] [Indexed: 12/24/2022]
Abstract
In this work, we present a fully automated method for the construction of chemically meaningful sets of hierarchical nonredundant internal coordinates (ICs; also commonly denoted as Z-matrices) from the Cartesian coordinates of a molecular system. Particular focus is placed on avoiding ill-definitions of angles and dihedrals due to linear arrangements of atoms, to consistently guarantee a well-defined transformation to Cartesian coordinates, even after structural changes. The representations thus obtained are particularly well suited for pathway construction in double-ended methods for transition state search and optimizations with nonlinear constraints. Analytical gradients for the transformation between the coordinate systems were derived for analytical geometry optimizations purely in Z-matrix coordinates. The geometry optimization was coupled with a Symbolic Algebra package to support arbitrary nonlinear constraints in Z-matrix coordinates, while retaining analytical energy gradient conversion. The difference to the commonly used nonhierarchical IC transformations is discussed. Sample applications are provided for a number of common chemical reactions and illustrative examples.
Collapse
Affiliation(s)
- Oskar Weser
- Electronic Structure Theory Department, Max-Planck-Institute for Solid State Research, Stuttgart, Germany.,Institute of Physical Chemistry, University of Goettingen, Goettingen, Germany
| | - Björn Hein-Janke
- Institute of Physical Chemistry, University of Goettingen, Goettingen, Germany
| | - Ricardo A Mata
- Institute of Physical Chemistry, University of Goettingen, Goettingen, Germany
| |
Collapse
|
108
|
Kee CW. Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities. Molecules 2023; 28:molecules28041715. [PMID: 36838703 PMCID: PMC9966076 DOI: 10.3390/molecules28041715] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2022] [Revised: 02/03/2023] [Accepted: 02/05/2023] [Indexed: 02/25/2023] Open
Abstract
Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.
Collapse
Affiliation(s)
- Choon Wee Kee
- Institute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR), 1 Pesek Road, Jurong Island, Singapore 627833, Republic of Singapore
| |
Collapse
|
109
|
Sargent CT, Metcalf DP, Glick ZL, Borca CH, Sherrill CD. Benchmarking two-body contributions to crystal lattice energies and a range-dependent assessment of approximate methods. J Chem Phys 2023; 158:054112. [PMID: 36754814 DOI: 10.1063/5.0141872] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Using the many-body expansion to predict crystal lattice energies (CLEs), a pleasantly parallel process, allows for flexibility in the choice of theoretical methods. Benchmark-level two-body contributions to CLEs of 23 molecular crystals have been computed using interaction energies of dimers with minimum inter-monomer separations (i.e., closest contact distances) up to 30 Å. In a search for ways to reduce the computational expense of calculating accurate CLEs, we have computed these two-body contributions with 15 different quantum chemical levels of theory and compared these energies to those computed with coupled-cluster in the complete basis set (CBS) limit. Interaction energies of the more distant dimers are easier to compute accurately and several of the methods tested are suitable as replacements for coupled-cluster through perturbative triples for all but the closest dimers. For our dataset, sub-kJ mol-1 accuracy can be obtained when calculating two-body interaction energies of dimers with separations shorter than 4 Å with coupled-cluster with single, double, and perturbative triple excitations/CBS and dimers with separations longer than 4 Å with MP2.5/aug-cc-pVDZ, among other schemes, reducing the number of dimers to be computed with coupled-cluster by as much as 98%.
Collapse
Affiliation(s)
- Caroline T Sargent
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Derek P Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Carlos H Borca
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
110
|
Kutta RJ, Großkopf J, van Staalduinen N, Seitz A, Pracht P, Breitenlechner S, Bannwarth C, Nuernberger P, Bach T. Multifaceted View on the Mechanism of a Photochemical Deracemization Reaction. J Am Chem Soc 2023; 145:2354-2363. [PMID: 36660908 DOI: 10.1021/jacs.2c11265] [Citation(s) in RCA: 14] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Upon irradiation in the presence of a chiral benzophenone catalyst (5 mol %), a racemic mixture of a given chiral imidazolidine-2,4-dione (hydantoin) can be converted almost quantitatively into the same compound with high enantiomeric excess (80-99% ee). The mechanism of this photochemical deracemization reaction was elucidated by a suite of mechanistic experiments. It was corroborated by nuclear magnetic resonance titration that the catalyst binds the two enantiomers by two-point hydrogen bonding. In one of the diastereomeric complexes, the hydrogen atom at the stereogenic carbon atom is ideally positioned for hydrogen atom transfer (HAT) to the photoexcited benzophenone. Detection of the protonated ketyl radical by transient absorption revealed hydrogen abstraction to occur from only one but not from the other hydantoin enantiomer. Quantum chemical calculations allowed us to visualize the HAT within this complex and, more importantly, showed that the back HAT does not occur to the carbon atom of the hydantoin radical but to its oxygen atom. The achiral enol formed in this process could be directly monitored by its characteristic transient absorption signal at λ ≅ 330 nm. Subsequent tautomerization leads to both hydantoin enantiomers, but only one of them returns to the catalytic cycle, thus leading to an enrichment of the other enantiomer. The data are fully consistent with deuterium labeling experiments and deliver a detailed picture of a synthetically useful photochemical deracemization reaction.
Collapse
Affiliation(s)
- Roger Jan Kutta
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Universitätsstr. 31, RegensburgD-93053, Germany
| | - Johannes Großkopf
- Department Chemie and Catalysis Research Center (CRC), School of Natural Sciences, Technische Universität München, D-85747Garching, Germany
| | - Nils van Staalduinen
- Institut für Physikalische Chemie, RWTH Aachen University, D-52074Aachen, Germany
| | - Antonia Seitz
- Department Chemie and Catalysis Research Center (CRC), School of Natural Sciences, Technische Universität München, D-85747Garching, Germany
| | - Philipp Pracht
- Institut für Physikalische Chemie, RWTH Aachen University, D-52074Aachen, Germany.,Yusuf Hamied Department of Chemistry, University of Cambridge, CambridgeCB2 1EW, United Kingdom
| | - Stefan Breitenlechner
- Department Chemie and Catalysis Research Center (CRC), School of Natural Sciences, Technische Universität München, D-85747Garching, Germany
| | - Christoph Bannwarth
- Institut für Physikalische Chemie, RWTH Aachen University, D-52074Aachen, Germany
| | - Patrick Nuernberger
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Universitätsstr. 31, RegensburgD-93053, Germany
| | - Thorsten Bach
- Department Chemie and Catalysis Research Center (CRC), School of Natural Sciences, Technische Universität München, D-85747Garching, Germany
| |
Collapse
|
111
|
Inversion Theory Leveling as a New Methodological Approach to Antioxidant Thermodynamics: A Case Study on Phenol. Antioxidants (Basel) 2023; 12:antiox12020282. [PMID: 36829841 PMCID: PMC9952401 DOI: 10.3390/antiox12020282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Revised: 01/13/2023] [Accepted: 01/18/2023] [Indexed: 01/31/2023] Open
Abstract
Antioxidants are various types of compounds that represent a link between biology and chemistry. With the development of theoretical and computational methods, antioxidants are now being studied theoretically. Here, a novel method is presented that aims to reduce the estimated wall times for DFT calculations that result in the same or higher degree of accuracy in the second derivatives over energy than is the case with the regular computational route (i.e., optimizing the reaction system at a lower model and then recalculating the energies at a higher level of theory) by applying the inversion of theory level to the universal chemical scavenger model, i.e., phenol. The resulting accuracy and wall time obtained with such a methodological setup strongly suggest that this methodology could be generally applied to antioxidant thermodynamics for some costly DFT methods with relative absolute deviation.
Collapse
|
112
|
Wilming FM, Marazzi B, Debes PP, Becker J, Schreiner PR. Probing the Size Limit of Dispersion Energy Donors with a Bifluorenylidene Balance: Magic Cyclohexyl. J Org Chem 2023; 88:1024-1035. [PMID: 36576961 DOI: 10.1021/acs.joc.2c02444] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We report the synthesis of 14 2,2'-disubstituted 9,9'-bifluorenylidenes as molecular balances for the quantification of London dispersion interactions between various dispersion energy donors. For all balances, we measured ΔGZ/E at 333 K using 1H NMR in seven organic solvents. For various alkyl and aryl substituents, we generally observe a preference for the "folded" Z-isomer due to attractive London dispersion interactions. The cyclohexyl-substituted system shows the largest Z-preference in this study with ΔGZ/E = -0.6 ± 0.05 kcal mol-1 in all solvents, owing to the rotational freedom of cyclohexyl groups paired with their large polarizability that maximizes London dispersion interactions. On the other hand, rigid and sterically more demanding substituents like tert-butyl unexpectedly favor the unfolded E-isomer. This is a result of the close relative position in which the functional groups are positioned in this molecular balance. This close proximity is the reason for the increase of Pauli repulsion in the Z-isomers with large rigid substituents (tert-butyl, adamantyl, and diamantyl) which leads to an equilibrium shift toward the unfolded E-form. While we were able to reproduce most of our experimental trends qualitatively using contemporary computational chemistry methods, quantitative accuracy of the employed methods still needs further improvement.
Collapse
Affiliation(s)
- Finn M Wilming
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany.,Center for Materials Research (ZfM), Justus Liebig University, Heinrich-Buff-Ring 16, 35392 Giessen, Germany
| | - Benito Marazzi
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany.,Center for Materials Research (ZfM), Justus Liebig University, Heinrich-Buff-Ring 16, 35392 Giessen, Germany
| | - Paul P Debes
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany.,Center for Materials Research (ZfM), Justus Liebig University, Heinrich-Buff-Ring 16, 35392 Giessen, Germany
| | - Jonathan Becker
- Institute of Inorganic and Analytical Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| | - Peter R Schreiner
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany.,Center for Materials Research (ZfM), Justus Liebig University, Heinrich-Buff-Ring 16, 35392 Giessen, Germany
| |
Collapse
|
113
|
de Oliveira MT, Alves JMA, Vrech NL, Braga AAC, Barboza CA. A comprehensive benchmark investigation of quantum chemical methods for carbocations. Phys Chem Chem Phys 2023; 25:1903-1922. [PMID: 36541431 DOI: 10.1039/d2cp04603b] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The application of various density functional approximations (DFAs) and an emphasis on popular methods without any consensus have prevailed in computational studies dedicated to carbocations. More importantly, an extensive and rigorous benchmark investigation on density functionals for the class is still lacking. To close this gap, we present a comprehensive benchmark study of quantum chemical methods on a series of classical and nonclassical carbocations, the CARBO33 dataset. We evaluate a total of 107 DFT methods from all rungs giving particular attention to double hybrid density functionals as the potential of the class has been largely undermined in the context of carbocations. To support our findings, DLPNO-CCSD(T) at the complete basis set (CBS) limit and W1-F12 are used as reference methods. Our results indicate that the composite CBS-QB3 method performs poorly and should not be adopted for target energies. Oftentimes, the tested DFAs of a lower rung perform better than several DFAs in a higher rung of Perdew's "Jacob's ladder". Nonetheless, double hybrids DSD-PBEP86-NL and ωB97X-2-D3(BJ) stand out by showing the overall best performance. Among the hybrids evaluated, about half of them show mean absolute deviation (MAD) below 1.1 kcal mol-1, including the popular hybrids M06-2X and mPW1PW91. In this family, MN15-D3(BJ) performs particularly well (MAD = 0.77 kcal mol-1) displaying reliable results across various tests. Highly popular B3LYP exhibited one of the worst performances (MAD = 4.74 kcal mol-1), and we do not recommend its application to carbocations. We also assess the 24 general-purpose basis sets of single- up to quadruple-ζ quality. The best compromise between accuracy and computational cost is achieved with cc-pVTZ followed by def2-TZVP. Computations on larger structures of general interest, including terpene carbocations, are also presented for selected DFT methods confirming general trends in the results.
Collapse
Affiliation(s)
- Marcelo T de Oliveira
- Department of Chemistry and Physics, La Trobe Institute of Molecular Sciences, La Trobe University, Melbourne, Victoria 3086, Australia. .,Chemistry Institute of São Carlos, University of São Paulo, Av. Trabalhador São Carlense 400, 13566-590, São Carlos, SP, Brazil
| | - Júlia M A Alves
- Chemistry Institute of São Carlos, University of São Paulo, Av. Trabalhador São Carlense 400, 13566-590, São Carlos, SP, Brazil
| | - Natália L Vrech
- Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-000, São Paulo, SP, Brazil
| | - Ataualpa A C Braga
- Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-000, São Paulo, SP, Brazil
| | - Cristina A Barboza
- Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland.,Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, Warsaw, Poland
| |
Collapse
|
114
|
Zaccaria M, Dawson W, Russel Kish D, Reverberi M, Bonaccorsi di Patti MC, Domin M, Cristiglio V, Chan B, Dellafiora L, Gabel F, Nakajima T, Genovese L, Momeni B. Experimental-theoretical study of laccase as a detoxifier of aflatoxins. Sci Rep 2023; 13:860. [PMID: 36650163 PMCID: PMC9845376 DOI: 10.1038/s41598-023-27519-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Accepted: 01/03/2023] [Indexed: 01/19/2023] Open
Abstract
We investigate laccase-mediated detoxification of aflatoxins, fungal carcinogenic food contaminants. Our experimental comparison between two aflatoxins with similar structures (AFB1 and AFG2) shows significant differences in laccase-mediated detoxification. A multi-scale modeling approach (Docking, Molecular Dynamics, and Density Functional Theory) identifies the highly substrate-specific changes required to improve laccase detoxifying performance. We employ a large-scale density functional theory-based approach, involving more than 7000 atoms, to identify the amino acid residues that determine the affinity of laccase for aflatoxins. From this study we conclude: (1) AFB1 is more challenging to degrade, to the point of complete degradation stalling; (2) AFG2 is easier to degrade by laccase due to its lack of side products and favorable binding dynamics; and (3) ample opportunities to optimize laccase for aflatoxin degradation exist, especially via mutations leading to π-π stacking. This study identifies a way to optimize laccase for aflatoxin bioremediation and, more generally, contributes to the research efforts aimed at rational enzyme optimization.
Collapse
Affiliation(s)
- Marco Zaccaria
- Department of Biology, Boston College, Chestnut Hill, MA, 02467, USA
| | - William Dawson
- RIKEN Center for Computational Science, Kobe, 6500047, Japan
| | | | - Massimo Reverberi
- Department of Environmental and Evolutionary Biology, "Sapienza" University of Rome, 00185, Rome, Italy
| | | | - Marek Domin
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | | | - Bun Chan
- RIKEN Center for Computational Science, Kobe, 6500047, Japan.,Graduate School of Engineering, Nagasaki University, Nagasaki, 8528521, Japan
| | - Luca Dellafiora
- Department of Food and Drug, University of Parma, 43124, Parma, Italy
| | - Frank Gabel
- CEA/CNRS/IBS, University Grenoble Alpes, 38044, Grenoble, France
| | | | - Luigi Genovese
- CEA/INAC-MEM/L-Sim, University Grenoble Alpes, 38044, Grenoble, France
| | - Babak Momeni
- Department of Biology, Boston College, Chestnut Hill, MA, 02467, USA.
| |
Collapse
|
115
|
Li TE, Hammes-Schiffer S. QM/MM Modeling of Vibrational Polariton Induced Energy Transfer and Chemical Dynamics. J Am Chem Soc 2023; 145:377-384. [PMID: 36574620 DOI: 10.1021/jacs.2c10170] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Vibrational strong coupling (VSC) provides a novel means to modify chemical reactions and energy transfer pathways. To efficiently model chemical dynamics under VSC in the collective regime, herein a hybrid quantum mechanical/molecular mechanical (QM/MM) cavity molecular dynamics (CavMD) scheme is developed and applied to an experimentally studied chemical system. This approach can achieve linear scaling with respect to the number of molecules for a dilute solution under VSC by assuming that each QM solute molecule is surrounded by an independent MM solvent bath. Application of this approach to a dilute solution of Fe(CO)5 in n-dodecane under VSC demonstrates polariton dephasing to the dark modes and polariton-enhanced molecular nonlinear absorption. These simulations predict that strongly exciting the lower polariton may provide an energy transfer pathway that selectively excites the equatorial CO vibrations rather than the axial CO vibrations. Moreover, these simulations also directly probe the cavity effect on the dynamics of the Fe(CO)5 Berry pseudorotation reaction for comparison to recent two-dimensional infrared spectroscopy experiments. This theoretical approach is applicable to a wide range of other polaritonic systems and provides a tool for exploring the use of VSC for selective infrared photochemistry.
Collapse
Affiliation(s)
- Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
116
|
Müller M, Hansen A, Grimme S. ωB97X-3c: A composite range-separated hybrid DFT method with a molecule-optimized polarized valence double-ζ basis set. J Chem Phys 2023; 158:014103. [PMID: 36610980 DOI: 10.1063/5.0133026] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
A new composite density functional theory (DFT) method is presented. It is based on ωB97X-V as one of the best-performing density functionals for the GMTKN55 thermochemistry database and completes the family of "3c" methods toward range-separated hybrid DFT. This method is consistently available for all elements up to Rn (Z = 1-86). Its further key ingredients are a polarized valence double-ζ (vDZP) Gaussian basis set, which was fully optimized in molecular DFT calculations, in combination with large-core effective core potentials and a specially adapted D4 dispersion correction. Unlike most existing double-ζ atomic orbital sets, vDZP shows only small basis set superposition errors (BSSEs) and can compete with standard sets of triple-ζ quality. Small residual BSSE effects are efficiently absorbed by the D4 damping scheme, which overall eliminates the need for an explicit treatment or empirical corrections for BSSE. Thorough tests on a variety of thermochemistry benchmark sets show that the new composite method, dubbed ωB97X-3c, is on par with or even outperforms standard hybrid DFT methods in a quadruple-zeta basis set at a small fraction of the computational cost. Particular strengths of this method are the description of non-covalent interactions and barrier heights, for which it is among the best-performing density functionals overall.
Collapse
Affiliation(s)
- Marcel Müller
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
117
|
Dulov DA, Bogdanov AV, Dorofeev SG, Magdesieva TV. N, N'-Diaryldihydrophenazines as a Sustainable and Cost-Effective Alternative to Precious Metal Complexes in the Photoredox-Catalyzed Alkylation of Aryl Alkyl Ketones. Molecules 2022; 28:molecules28010221. [PMID: 36615415 PMCID: PMC9822323 DOI: 10.3390/molecules28010221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Revised: 12/20/2022] [Accepted: 12/21/2022] [Indexed: 12/28/2022] Open
Abstract
An inexpensive and highly efficient metal-free alternative to commonly used Ru- and Ir-based catalysts was proposed. It was shown that the new 2,7-di-tert-butyl-5,10-bis(4-trifluoromethylphenyl)-5,10-dihydrophenazine outcompeted the iridium phenylpyridyl complex in photoredox activity in the alkylation of silyl enol ethers yielding aryl alkyl ketones. The reaction occurred under visible light irradiation at room temperature and was also applicable to drug derivatives (ibuprofen and naproxen). In-depth photophysical, electrochemical, and quantum chemical studies showed that the aforementioned N,N-diaryldihydrophenazine exhibited enhanced properties that were essential for the photoredox catalysis (a long-lived triplet excited state, strong reducing ability, high stability of the radical cations formed in single-electron-transfer event, and chemical inertness of the catalyst with respect to reactants). Importantly, the substituted N,N'-diaryldihydrophenazines could be obtained directly from diaryl amines; a facile, easily handled and scaled-up one-pot synthetic procedure was elaborated.
Collapse
|
118
|
Lin J, Coles NT, Dettling L, Steiner L, Felix Witte J, Paulus B, Müller C. Phospholenes from Phosphabenzenes by Selective Ring Contraction. Chemistry 2022; 28:e202203406. [PMID: 36383093 PMCID: PMC10107248 DOI: 10.1002/chem.202203406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Revised: 11/15/2022] [Accepted: 11/15/2022] [Indexed: 11/17/2022]
Abstract
A 3-amino-functionalized phosphabenzene (phosphinine) has been synthesized and structurally characterized. The pyramidalized nitrogen atom of the dimethylamino substituent indicates only a weak interaction between the lone pair of the nitrogen atom and the aromatic phosphorus heterocycle, resulting in somewhat basic character. It turned out that the amino group can indeed be protonated by HCl. In contrast to pyridines, however, the phosphabenzene-ammonium salt undergoes a selective ring contraction to form a hydroxylphospholene oxide in the presence of additional water. Based on deuterium labeling experiments and quantum chemical calculations, a rational mechanism for this hitherto unknown conversion is proposed.
Collapse
Affiliation(s)
- Jinxiong Lin
- Freie Universität Berlin, Institut für Chemie und Biochemie, Fabeckstrasse 34/36, 14195, Berlin, Germany
| | - Nathan T Coles
- Freie Universität Berlin, Institut für Chemie und Biochemie, Fabeckstrasse 34/36, 14195, Berlin, Germany.,School of Chemistry, University of Nottingham University Park, Nottingham, NG7 2RD, UK
| | - Lea Dettling
- Freie Universität Berlin, Institut für Chemie und Biochemie, Fabeckstrasse 34/36, 14195, Berlin, Germany
| | - Luca Steiner
- Freie Universität Berlin, Institut für Chemie und Biochemie, Arnimallee 22, 14195, Berlin, Germany
| | - J Felix Witte
- Freie Universität Berlin, Institut für Chemie und Biochemie, Arnimallee 22, 14195, Berlin, Germany
| | - Beate Paulus
- Freie Universität Berlin, Institut für Chemie und Biochemie, Arnimallee 22, 14195, Berlin, Germany
| | - Christian Müller
- Freie Universität Berlin, Institut für Chemie und Biochemie, Fabeckstrasse 34/36, 14195, Berlin, Germany
| |
Collapse
|
119
|
Leach IF, Sorbelli D, Belpassi L, Belanzoni P, Havenith RWA, Klein JEMN. How reduced are nucleophilic gold complexes? Dalton Trans 2022; 52:11-15. [PMID: 35877065 PMCID: PMC9764324 DOI: 10.1039/d2dt01694j] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Nucleophilic formal gold(-I) and gold(I) complexes are investigated via Intrinsic Bond Orbital analysis and Energy Decomposition Analysis, based on density functional theory calculations. The results indicate gold(0) centres engaging in electron-sharing bonding with Al- and B- based ligands. Multiconfigurational (CASSCF) calculations corroborate the findings, highlighting the gap between the electonic structures and the oxidation state formalism.
Collapse
Affiliation(s)
- Isaac F. Leach
- Molecular Inorganic Chemistry, Stratingh Institute for Chemistry, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands,Zernike Institute for Advanced Materials, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands
| | - Diego Sorbelli
- Department of Chemistry, Biology and Biotechnology, University of PerugiaVia Elce di Sotto806123 PerugiaItaly,CNR Institute of Chemical Science and Technologies “Giulio Natta” (CNR-SCITEC)Via Elce di Sotto806123 PerugiaItaly
| | - Leonardo Belpassi
- CNR Institute of Chemical Science and Technologies “Giulio Natta” (CNR-SCITEC)Via Elce di Sotto806123 PerugiaItaly
| | - Paola Belanzoni
- Department of Chemistry, Biology and Biotechnology, University of PerugiaVia Elce di Sotto806123 PerugiaItaly,CNR Institute of Chemical Science and Technologies “Giulio Natta” (CNR-SCITEC)Via Elce di Sotto806123 PerugiaItaly
| | - Remco W. A. Havenith
- Molecular Inorganic Chemistry, Stratingh Institute for Chemistry, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands,Zernike Institute for Advanced Materials, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands,Ghent Quantum Chemistry Group, Department of Chemistry, Ghent University9000 GentBelgium
| | - Johannes E. M. N. Klein
- Molecular Inorganic Chemistry, Stratingh Institute for Chemistry, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands
| |
Collapse
|
120
|
Jensen AB, Kubečka J, Schmitz G, Christiansen O, Elm J. Massive Assessment of the Binding Energies of Atmospheric Molecular Clusters. J Chem Theory Comput 2022; 18:7373-7383. [PMID: 36417753 DOI: 10.1021/acs.jctc.2c00825] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Quantum chemical studies of the formation and growth of atmospheric molecular clusters are important for understanding aerosol particle formation. However, the search for the lowest free-energy cluster configuration is extremely time consuming. This makes high-level benchmark data sets extremely valuable in the quest for the global minimum as it allows the identification of cost-efficient computational methodologies, as well as the development of high-level machine learning (ML) models. Herein, we present a highly versatile quantum chemical data set comprising a total of 11 749 (acid)1-2(base)1-2 cluster configurations, containing up to 44 atoms. Utilizing the LUMI supercomputer, we calculated highly accurate PNO-CCSD(F12*)(T)/cc-pVDZ-F12 binding energies of the full set of cluster configurations leading to an unprecedented data set both in regard to sheer size and with respect to the level of theory. We employ the constructed benchmark set to assess the performance of various semiempirical and density functional theory methods. In particular, we find that the r2-SCAN-3c method shows excellent performance across the data set related to both accuracy and CPU time, making it a promising method to employ during cluster configurational sampling. Furthermore, applying the data sets, we construct ML models based on Δ-learning and provide recommendations for future application of ML in cluster configurational sampling.
Collapse
Affiliation(s)
| | - Jakub Kubečka
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Gunnar Schmitz
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Jonas Elm
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark.,Department of Chemistry, iClimate, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| |
Collapse
|
121
|
Froitzheim T, Grimme S, Mewes JM. Either Accurate Singlet-Triplet Gaps or Excited-State Structures: Testing and Understanding the Performance of TD-DFT for TADF Emitters. J Chem Theory Comput 2022; 18:7702-7713. [PMID: 36409831 DOI: 10.1021/acs.jctc.2c00905] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The energy gap between the lowest singlet and triplet excited states (ΔEST) is a key property of thermally activated delayed fluorescence (TADF) emitters, where these states are dominated by charge-transfer (CT) character. Despite its well-known shortcomings concerning CT states, time-dependent density functional theory (TD-DFT) is widely used to predict this gap and study TADF. Moreover, polar CT states exhibit a strong interaction with their molecular environment, which further complicates their computational description. Addressing these two major challenges, this work studies the performance of Tamm-Dancoff-approximated TD-DFT (TDA-DFT) on the recent STGABS27 benchmark set,1 exploring different strategies to include orbital and structural relaxation, as well as dielectric embedding. The results show that the best-performing strategy is to calculate ΔEST at the ground-state structure using functionals with a surprisingly small amount of Fock exchange of ≈10% and without a (complete) solvent model. However, as this approach heavily relies on error cancellation to mimic dielectric relaxation, it is not robust and exhibits large systematic deviations in excited state energies, state characters, and structures. More rigorous approaches, including state-specific solvation, do not share these systematic deviations, but their predicted ΔEST values exhibit larger statistical errors. We thus conclude that for the description of CT states in dielectric environments, none of the tested TDA-DFT methods is competitive with the recently presented ROKS/PCM approach regarding robustness, accuracy, and computational efficiency.
Collapse
Affiliation(s)
- Thomas Froitzheim
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115Bonn, Germany
| | - Jan-Michael Mewes
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115Bonn, Germany
| |
Collapse
|
122
|
Mráziková K, Kruse H, Mlýnský V, Auffinger P, Šponer J. Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions. J Chem Inf Model 2022; 62:6182-6200. [PMID: 36454943 DOI: 10.1021/acs.jcim.2c01064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.
Collapse
Affiliation(s)
- Klaudia Mráziková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| | - Pascal Auffinger
- Architecture et Réactivité de l'ARN, Université de Strasbourg, Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg67084, France
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| |
Collapse
|
123
|
Gorges J, Bädorf B, Hansen A, Grimme S. Efficient Computation of the Interaction Energies of Very Large Non-covalently Bound Complexes. Synlett 2022. [DOI: 10.1055/s-0042-1753141] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
AbstractWe present a new benchmark set consisting of 16 large non-covalently bound systems (LNCI16) ranging from 380 up to 1988 atoms and featuring diverse interaction motives. Gas-phase interaction energies are calculated with various composite DFT, semi-empirical quantum mechanical (SQM), and force field (FF) methods and are evaluated using accurate DFT reference values. Of the employed QM methods, PBEh-3c proves to be the most robust for large systems with a relative mean absolute deviation (relMAD) of 8.5% with respect to the reference interaction energies. r2SCAN-3c yields an even smaller relMAD, at least for the subset of complexes for which the calculation could be converged, but is less robust for systems with smaller HOMO–LUMO gaps. The inclusion of Fock-exchange is therefore important for the description of very large non-covalent interaction (NCI) complexes in the gas phase. GFN2-xTB was found to be the best performer of the SQM methods with an excellent result of only 11.1% deviation. From the assessed force fields, GFN-FF and GAFF achieve the best accuracy. Considering their low computational costs, both can be recommended for routine calculations of very large NCI complexes, with GFN-FF being clearly superior in terms of general applicability. Hence, GFN-FF may be routinely applied in supramolecular synthesis planning.1 Introduction2 The LNCI16 Benchmark Set3 Computational Details4 Generation of Reference Values5 Results and Discussion6 Conclusions
Collapse
|
124
|
Choudhari M, Xu J, McKay AI, Guerrin C, Forsyth C, Ma HZ, Goerigk L, O'Hair RAJ, Bonnefont A, Ruhlmann L, Aloise S, Ritchie C. A photo-switchable molecular capsule: sequential photoinduced processes. Chem Sci 2022; 13:13732-13740. [PMID: 36544719 PMCID: PMC9710229 DOI: 10.1039/d2sc04613j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Accepted: 10/21/2022] [Indexed: 12/24/2022] Open
Abstract
The metastable trilacunary heteropolyoxomolybdate [PMo9O31(py)3]3- - {PMo9}; py = pyridine) and the ditopic pyridyl bearing diarylethene (DAE) (C25H16N2F6S2) self-assemble via a facile ligand replacement methodology to yield the photo-active molecular capsule [(PMo9O31)2(DAE)3]6-. The spatial arrangement and conformation of the three DAE ligands are directed by the surface chemistry of the molecular metal oxide precursor with exclusive ligation of the photo-active antiparallel rotamer to the polyoxometalate (POM) while the integrity of the assembly in solution has been verified by a suite of spectroscopic techniques. Electrocyclisation of the three DAEs occurs sequentially and has been investigated using a combination of steady-state and time-resolved spectroscopies with the discovery of a photochemical cascade whereby rapid photoinduced ring closure is followed by electron transfer from the ring-closed DAE to the POM in the latent donor-acceptor system on subsequent excitation. This interpretation is also supported by computational and detailed spectroelectrochemical analysis. Ring-closing quantum yields were also determined using a custom quantum yield determination setup (QYDS), providing insight into the impact of POM coordination on these processes.
Collapse
Affiliation(s)
| | - Jingjing Xu
- School of Chemistry, Monash UniversityClaytonVictoria3800Australia,College of Chemistry and Materials Science, Shanghai Normal UniversityShanghai 200234China
| | | | - Clément Guerrin
- Université de Lille, CNRS, UMR 8516 – LASIRE – Laboratoire de Spectroscopie pour les Interactions, la Réactivité et l’EnvironnementLille F-59000France
| | - Craig Forsyth
- School of Chemistry, Monash UniversityClaytonVictoria3800Australia
| | - Howard Z. Ma
- School of Chemistry, University of MelbourneMelbourne3010VictoriaAustralia
| | - Lars Goerigk
- School of Chemistry, University of MelbourneMelbourne3010VictoriaAustralia
| | | | - Antoine Bonnefont
- Laboratoire d’Electrochimie et Chimie Physique du Corps Solide, Institut de Chimie de Strasbourg, UMR 7177 CNRS-Université de Strasbourg4, rue Blaise PascalStrasbourg 67000France
| | - Laurent Ruhlmann
- Laboratoire d’Electrochimie et Chimie Physique du Corps Solide, Institut de Chimie de Strasbourg, UMR 7177 CNRS-Université de Strasbourg4, rue Blaise PascalStrasbourg 67000France
| | - Stephane Aloise
- Université de Lille, CNRS, UMR 8516 – LASIRE – Laboratoire de Spectroscopie pour les Interactions, la Réactivité et l’EnvironnementLille F-59000France
| | - Chris Ritchie
- School of Chemistry, Monash UniversityClaytonVictoria3800Australia
| |
Collapse
|
125
|
Rosero-Mafla MA, Zapata-Rivera J, Gimeno MC, Visbal R. Steric and Electronic Effects in N-Heterocyclic Carbene Gold(III) Complexes: An Experimental and Computational Study. Molecules 2022; 27:molecules27238289. [PMID: 36500397 PMCID: PMC9740751 DOI: 10.3390/molecules27238289] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Revised: 11/20/2022] [Accepted: 11/22/2022] [Indexed: 11/29/2022] Open
Abstract
A series of neutral acridine-based gold(III)-NHC complexes containing the pentafluorophenyl (-C6F5) group were synthesized. All of the complexes were fully characterized by analytical techniques. The square planar geometry around the gold center was confirmed by X-ray diffraction analysis for complexes 1 (Trichloro [1-methyl-3-(9-acridine)imidazol-2-ylidene]gold(III)) and 2 (Chloro-bis(pentafluorophenyl)[1-methyl-3-(9-acridine)imidazol-2-ylidene]gold(III)). In both cases, the acridine rings play a key role in the crystal packing of the solid structures by mean of π-π stacking interactions, with centroid-centroid and interplanar distances being similar to those found in other previously reported acridine-based Au(I)-NHC complexes. A different reactivity when using a bulkier N-heterocyclic carbene ligand such as 1,3-bis-(2,6-diisopropylphenyl)-2-imidazolidinylidene (SIPr) was observed. While the use of the acridine-based NHC ligand led to the expected organometallic gold(III) species, the steric hindrance of the bulky SIPr ligand led to the formation of the corresponding imidazolinium cation stabilized by the tetrakis(pentafluorophenyl)aurate(III) [Au(C6F5)4]- anion. Computational experiments were carried out in order to figure out the ground state electronic structure and the binding formation energy of the complexes and, therefore, to explain the observed reactivity.
Collapse
Affiliation(s)
- Miguel A. Rosero-Mafla
- Departamento de Química, Facultad de Ciencias Naturales y Exactas, Universidad del Valle, A.A. 25360, Cali 760042, Colombia
| | - Jhon Zapata-Rivera
- Departamento de Química, Facultad de Ciencias, Universidad de los Andes, Cra 1 No 18A—12, Bogotá 111711, Colombia
| | - M. Concepción Gimeno
- Departamento de Química Inorgánica, Instituto de Síntesis Química y Catálisis Homogénea (ISQCH), CSIC-Universidad de Zaragoza, 50009 Zaragoza, Spain
- Correspondence: (M.C.G.); (R.V.)
| | - Renso Visbal
- Departamento de Química, Facultad de Ciencias Naturales y Exactas, Universidad del Valle, A.A. 25360, Cali 760042, Colombia
- Centro de Excelencia en Nuevos Materiales (CENM), Universidad del Valle, A.A. 25360, Cali 760031, Colombia
- Correspondence: (M.C.G.); (R.V.)
| |
Collapse
|
126
|
Seeber P, Seidenath S, Steinmetzer J, Gräfe S. Growing Spicy
ONIOMs
: Extending and generalizing concepts of
ONIOM
and many body expansions. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1644] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Affiliation(s)
- Phillip Seeber
- Institute of Physical Chemistry Friedrich Schiller University Jena Jena Germany
| | - Sebastian Seidenath
- Institute of Physical Chemistry Friedrich Schiller University Jena Jena Germany
| | | | - Stefanie Gräfe
- Institute of Physical Chemistry Friedrich Schiller University Jena Jena Germany
| |
Collapse
|
127
|
Cseh K, Geisler H, Stanojkovska K, Westermayr J, Brunmayr P, Wenisch D, Gajic N, Hejl M, Schaier M, Koellensperger G, Jakupec MA, Marquetand P, Kandioller W. Arene Variation of Highly Cytotoxic Tridentate Naphthoquinone-Based Ruthenium(II) Complexes and In-Depth In Vitro Studies. Pharmaceutics 2022; 14:2466. [PMID: 36432656 PMCID: PMC9699003 DOI: 10.3390/pharmaceutics14112466] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 11/02/2022] [Accepted: 11/05/2022] [Indexed: 11/18/2022] Open
Abstract
The main purpose of this study was to synthesize a new set of naphthoquinone-based ruthenium(II) arene complexes and to develop an understanding of their mode of action. This study systematically reviews the steps of synthesis, aiming to provide a simplified approach using microwave irradiation. The chemical structures and the physicochemical properties of this novel group of compounds were examined by 1H-NMR and 13C-NMR spectroscopy, X-ray diffractometry, HPLC-MS and supporting DFT calculations. Several aspects of the biological activity were investigated in vitro, including short- and long-term cytotoxicity tests, cellular accumulation studies, detection of reactive oxygen species generation, apoptosis induction and NAD(P)H:quinone oxidoreductase 1 (NQO1) activity as well as cell cycle analysis in A549, CH1/PA-1, and SW480 cancer cells. Furthermore, the DNA interaction ability was studied in a cell-free assay. A positive correlation was found between cytotoxicity, lipophilicity and cellular accumulation of the tested complexes, and the results offer some important insights into the effects of the arene. The most obvious finding to emerge from this study is that the usually very chemosensitive CH1/PA-1 teratocarcinoma cells showed resistance to these phthiocol-based organometallics in comparison to the usually less chemosensitive SW480 colon carcinoma cells, which pilot experiments suggest as being related to NQO1 activity.
Collapse
Affiliation(s)
- Klaudia Cseh
- Institute of Inorganic Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
| | - Heiko Geisler
- Institute of Inorganic Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
| | - Kristina Stanojkovska
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 17, 1090 Vienna, Austria
| | - Julia Westermayr
- Wilhelm-Ostwald-Institute for Physical and Theoretical Chemistry, Faculty of Chemistry and Mineralogy, University of Leipzig, Linnéstr. 2, 04103 Leipzig, Germany
| | - Philipp Brunmayr
- Institute of Inorganic Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
| | - Dominik Wenisch
- Institute of Inorganic Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
| | - Natalie Gajic
- Institute of Inorganic Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
| | - Michaela Hejl
- Institute of Inorganic Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
| | - Martin Schaier
- Institute of Analytical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 38, 1090 Vienna, Austria
- Vienna Doctoral School in Chemistry (DoSChem), University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
| | - Gunda Koellensperger
- Institute of Analytical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 38, 1090 Vienna, Austria
| | - Michael A. Jakupec
- Institute of Inorganic Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
- Research Cluster “Translational Cancer Therapy Research”, University of Vienna, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, 1090 Vienna, Austria
| | - Wolfgang Kandioller
- Institute of Inorganic Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 42, 1090 Vienna, Austria
- Research Cluster “Translational Cancer Therapy Research”, University of Vienna, 1090 Vienna, Austria
| |
Collapse
|
128
|
Göller AH. Reliable gas-phase tautomer equilibria of drug-like molecule scaffolds and the issue of continuum solvation. J Comput Aided Mol Des 2022; 36:805-824. [DOI: 10.1007/s10822-022-00480-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Accepted: 09/29/2022] [Indexed: 11/06/2022]
|
129
|
Mercurophilic Interactions in Heterometallic Ru-Hg carbonyl clusters. Inorganica Chim Acta 2022. [DOI: 10.1016/j.ica.2022.121281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
130
|
Paired electrosynthesis in the CH3CN/CHal4 system: a one-pot procedure for diarylamines and N,N’-diarylbenzidines halogenation. Electrochim Acta 2022. [DOI: 10.1016/j.electacta.2022.141217] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
131
|
Phatake RS, Averdunk A, Würtele C, Gellrich U. Piers’ Borane-Catalyzed Dimerization of Arylallenes via Transborylation: A Synthetic and Mechanistic Study. ACS Catal 2022. [DOI: 10.1021/acscatal.2c04605] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ravindra S. Phatake
- Institut für Organische Chemie, Justus-Liebig-Universität Gießen, Heinrich-Buff-Ring 17, 35392Gießen, Germany
| | - Arthur Averdunk
- Institut für Organische Chemie, Justus-Liebig-Universität Gießen, Heinrich-Buff-Ring 17, 35392Gießen, Germany
| | - Christian Würtele
- Institut für Anorganische und Analytische Chemie, Justus-Liebig-Universität Gießen, Heinrich-Buff-Ring 17, 35392Gießen, Germany
| | - Urs Gellrich
- Institut für Organische Chemie, Justus-Liebig-Universität Gießen, Heinrich-Buff-Ring 17, 35392Gießen, Germany
| |
Collapse
|
132
|
Kang CS, Zhang S, Wang H, Liu Y, Ren S, Chen Y, Li J, Bandara N, Rogachev AY, Rogers BE, Chong HS. Novel Chelating Agents for Zirconium-89-Positron Emission Tomography (PET) Imaging: Synthesis, DFT Calculation, Radiolabeling, and In Vitro and In Vivo Complex Stability. ACS OMEGA 2022; 7:37229-37236. [PMID: 36312419 PMCID: PMC9607665 DOI: 10.1021/acsomega.2c03478] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/03/2022] [Accepted: 08/19/2022] [Indexed: 06/16/2023]
Abstract
We report the synthesis and evaluation of novel chelating agents for zirconium-89 (89Zr) with positron emission tomography (PET) imaging applications. New chelating agents NODHA, NOTHA, and NODHA-PY were constructed on 1,4,7-triazacyclononane (TACN) and possess hydroxamic acid or a pyridine ring as an acyclic binding moiety. The new chelating agents were theoretically studied for complexation with Zr(IV). Structures of Zr(IV)-NODHA, Zr(IV)-NOTHA, and Zr(IV)-NODHA-PY were predicted using density functional methods. NODHA was found to form stronger bonds with Zr(IV) when compared to NOTHA and NODHA-PY. The new chelating agents were evaluated for radiolabeling efficiency in binding 89Zr. The corresponding [89Zr]Zr-labeled chelators were evaluated for complex stability in human serum. All new chelating agents rapidly bound to 89Zr in excellent radiolabeling efficiency at room temperature. Among the new [89Zr]Zr-labeled chelators evaluated, [89Zr]Zr-NODHA showed the highest stability in human serum without losing 89Zr, and [89Zr]Zr-NODHA-PY released a considerable amount of 89Zr in human serum. [89Zr]Zr-NODHA, [89Zr]Zr-NODHA-PY, and [89Zr]Zr-DFO were comparatively evaluated for in vivo complex stability by performing biodistribution studies using normal mice. [89Zr]Zr-DFO had the lowest bone uptake at all time points, while [89Zr]Zr-NODHA-PY showed poor stability in mice as evidenced by high bone accumulation at the 24 h time point. [89Zr]Zr-NODHA exhibited better renal clearance but higher bone uptake than [89Zr]Zr-DFO.
Collapse
Affiliation(s)
- Chi Soo Kang
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| | - Shuyuan Zhang
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| | - Haixing Wang
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| | - Yujie Liu
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| | - Siyuan Ren
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| | - Yanda Chen
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| | - Jingbai Li
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| | - Nilantha Bandara
- Department
of Radiation Oncology, Washington University
School of Medicine, St. Louis, Missouri 63100-1010, United States
| | - Andrey Yu Rogachev
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| | - Buck E. Rogers
- Department
of Radiation Oncology, Washington University
School of Medicine, St. Louis, Missouri 63100-1010, United States
| | - Hyun-Soon Chong
- Department
of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, United States
| |
Collapse
|
133
|
Stan A, Esch BVD, Ochsenfeld C. Fully Automated Generation of Prebiotically Relevant Reaction Networks from Optimized Nanoreactor Simulations. J Chem Theory Comput 2022; 18:6700-6712. [PMID: 36270030 DOI: 10.1021/acs.jctc.2c00754] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The nanoreactor approach first introduced by the group of Martı́nez [Wang et al. Nat. Chem. 2014, 6, 1044-1048] has recently attracted much attention because of its ability to accelerate the discovery of reaction pathways. Here, we provide a comprehensive study of various simulation parameters and present an alternative implementation for the reactivity-enhancing spherical constraint function, as well as for the detection of reaction events. In this context, a fully automated postsimulation evaluation procedure based on RDKit and NetworkX analysis is introduced. The chemical and physical robustness of the procedure is examined by investigating the reactivity of selected homogeneous systems. The optimized procedure is applied at the GFN2-xTB level of theory to a system composed of HCN molecules and argon atoms, acting as a buffer, yielding prebiotically plausible primary and secondary precursors for the synthesis of RNA. Furthermore, the formose reaction network is explored leading to numerous sugar precursors. The discovered compounds reflect experimental findings; however, new synthetic routes and a large collection of exotic, highly reactive molecules are observed, highlighting the predictive power of the nanoreactor approach for unraveling the reactive manifold.
Collapse
Affiliation(s)
- Alexandra Stan
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Beatriz von der Esch
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany.,Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
134
|
Abstract
![]()
Oxidation of a guanine nucleotide in DNA yields an 8-oxoguanine
nucleotide (oxoG) and is a mutagenic event in the genome.
Due to different arrangements of hydrogen-bond donors and acceptors, oxoG can affect the secondary structure of nucleic acids. We
have investigated base pairing preferences of oxoG in the
core of a tetrahelical G-quadruplex structure, adopted by analogues
of d(TG4T). Using spectroscopic methods, we have shown
that G-quartets can be fully substituted with oxoG nucleobases
to form an oxoG-quartet with a revamped hydrogen-bonding
scheme. While an oxoG-quartet can be incorporated into
the G-quadruplex core without distorting the phosphodiester backbone,
larger dimensions of the central cavity change the cation localization
and exchange properties.
Collapse
Affiliation(s)
- Simon Aleksič
- Slovenian NMR Centre, National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia.,Faculty of Chemistry and Chemical Technology, University of Ljubljana, Večna pot 113, 1000 Ljubljana, Slovenia
| | - Peter Podbevšek
- Slovenian NMR Centre, National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
| | - Janez Plavec
- Slovenian NMR Centre, National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia.,Faculty of Chemistry and Chemical Technology, University of Ljubljana, Večna pot 113, 1000 Ljubljana, Slovenia.,EN-FIST Centre of Excellence, Trg OF 13, 1000 Ljubljana, Slovenia
| |
Collapse
|
135
|
Bursch M, Mewes J, Hansen A, Grimme S. Best-Practice DFT Protocols for Basic Molecular Computational Chemistry. Angew Chem Int Ed Engl 2022; 61:e202205735. [PMID: 36103607 PMCID: PMC9826355 DOI: 10.1002/anie.202205735] [Citation(s) in RCA: 137] [Impact Index Per Article: 68.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2022] [Indexed: 01/11/2023]
Abstract
Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.
Collapse
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für KohlenforschungKaiser-Wilhelm-Platz 145470Mülheim an der RuhrGermany
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| | - Andreas Hansen
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| | - Stefan Grimme
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| |
Collapse
|
136
|
Jo T, Taschinski S, Leach IF, Bauer C, Hashmi ASK, Klein JEMN. On the Role of Noncovalent Ligand-Substrate Interactions in Au(I) Catalysis: An Experimental and Computational Study of Protodeauration. ACS Catal 2022; 12:13158-13163. [DOI: 10.1021/acscatal.2c03384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Revised: 09/25/2022] [Indexed: 11/29/2022]
Affiliation(s)
- Taegeun Jo
- Molecular Inorganic Chemistry, Stratingh Institute for Chemistry, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Svenja Taschinski
- Molecular Inorganic Chemistry, Stratingh Institute for Chemistry, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
- Organisch-Chemisches Institut, Heidelberg University, Im Neuenheimer Feld 270, 69120 Heidelberg, Germany
| | - Isaac F. Leach
- Molecular Inorganic Chemistry, Stratingh Institute for Chemistry, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Christina Bauer
- Molecular Inorganic Chemistry, Stratingh Institute for Chemistry, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
- Organisch-Chemisches Institut, Heidelberg University, Im Neuenheimer Feld 270, 69120 Heidelberg, Germany
| | - A. Stephen K. Hashmi
- Organisch-Chemisches Institut, Heidelberg University, Im Neuenheimer Feld 270, 69120 Heidelberg, Germany
| | - Johannes E. M. N. Klein
- Molecular Inorganic Chemistry, Stratingh Institute for Chemistry, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
137
|
Atahan-Evrenk S. Theoretical Study of the Structure and Binding Energies of Dimers of Zn(II)-Porphyrin Derivatives. J Phys Chem A 2022; 126:7102-7109. [PMID: 36194887 PMCID: PMC9574925 DOI: 10.1021/acs.jpca.2c03692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Zinc-complexed porphyrin and chlorophyll derivatives form functional aggregates with remarkable photophysical and optoelectronic properties. Understanding the type and strength of intermolecular interactions between these molecules is essential for designing new materials with desired morphology and functionality. The dimer interactions of a molecular set composed of porphyrin derivatives obtained by substitutional changes starting from free-base porphyrin is studied. It is found that the B97M-rV/def2-TZVP level of theory provides a good compromise between the accuracy and cost to get the dimer geometries and interaction energies (IEs). The neglect of the relaxation energy due to the change in the monomer configurations upon complex formation causes a more significant error than the basis set superposition error. The metal complexation increases the binding energy by about -6 to -8 kcal/mol, and the introduction of keto and hydroxy groups further stabilizes the dimers by about -20 kcal/mol. Although the saturation of one of the pyrrol double bonds does not change the IE, the addition of R groups increases it.
Collapse
Affiliation(s)
- Sule Atahan-Evrenk
- Faculty of Medicine, TOBB University of Economics and Technology, Sogutozu Cd No. 43 Sogutozu, Ankara06560, Turkey
| |
Collapse
|
138
|
Gillett AJ, Pershin A, Pandya R, Feldmann S, Sneyd AJ, Alvertis AM, Evans EW, Thomas TH, Cui LS, Drummond BH, Scholes GD, Olivier Y, Rao A, Friend RH, Beljonne D. Dielectric control of reverse intersystem crossing in thermally activated delayed fluorescence emitters. NATURE MATERIALS 2022; 21:1150-1157. [PMID: 35927434 PMCID: PMC7613666 DOI: 10.1038/s41563-022-01321-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Accepted: 06/13/2022] [Indexed: 05/22/2023]
Abstract
Thermally activated delayed fluorescence enables organic semiconductors with charge transfer-type excitons to convert dark triplet states into bright singlets via reverse intersystem crossing. However, thus far, the contribution from the dielectric environment has received insufficient attention. Here we study the role of the dielectric environment in a range of thermally activated delayed fluorescence materials with varying changes in dipole moment upon optical excitation. In dipolar emitters, we observe how environmental reorganization after excitation triggers the full charge transfer exciton formation, minimizing the singlet-triplet energy gap, with the emergence of two (reactant-inactive) modes acting as a vibrational fingerprint of the charge transfer product. In contrast, the dielectric environment plays a smaller role in less dipolar materials. The analysis of energy-time trajectories and their free-energy functions reveals that the dielectric environment substantially reduces the activation energy for reverse intersystem crossing in dipolar thermally activated delayed fluorescence emitters, increasing the reverse intersystem crossing rate by three orders of magnitude versus the isolated molecule.
Collapse
Affiliation(s)
| | - Anton Pershin
- Laboratory for Chemistry of Novel Materials, Université de Mons, Mons, Belgium
- Wigner Research Centre for Physics, Budapest, Hungary
| | - Raj Pandya
- Cavendish Laboratory, University of Cambridge, Cambridge, UK
| | - Sascha Feldmann
- Cavendish Laboratory, University of Cambridge, Cambridge, UK
| | | | | | - Emrys W Evans
- Cavendish Laboratory, University of Cambridge, Cambridge, UK
- Department of Chemistry, Swansea University, Swansea, UK
| | - Tudor H Thomas
- Cavendish Laboratory, University of Cambridge, Cambridge, UK
| | - Lin-Song Cui
- CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, China
| | | | | | - Yoann Olivier
- Unité de Chimie Physique Théorique et Structurale & Laboratoire de Physique du Solide, Namur Institute of Structured Matter, Université de Namur, Namur, Belgium
| | - Akshay Rao
- Cavendish Laboratory, University of Cambridge, Cambridge, UK
| | | | - David Beljonne
- Laboratory for Chemistry of Novel Materials, Université de Mons, Mons, Belgium.
| |
Collapse
|
139
|
Macías-Hernández CE, Romero-Chávez MM, Mojica-Sánchez JP, Pineda-Urbina K, Martínez MTS, Jimenez-Ruiz EI, Via LD, Ramos-Organillo Á. Synthesis and characterization of new monothiooxalamides containing pyridine nuclei with promising antiproliferative and antioxidant activity. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.133360] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
140
|
Bresciani G, Boni S, Zacchini S, Pampaloni G, Bortoluzzi M, Marchetti F. Alkyne-alkenyl coupling at a diruthenium complex. Dalton Trans 2022; 51:15703-15715. [PMID: 36177843 DOI: 10.1039/d2dt02866b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Dimetallic complexes are suitable platforms for the assembly of small molecular units, and the reactivity of bridging alkenyl ligands has been widely investigated to model C-C bond forming processes. Here, we report the unusual coupling of an alkenyl ligand, bridging coordinated on a diruthenium scaffold, with a series of alkynes, revealing two possible outcomes. The diruthenium complex [Ru2Cp2(Cl)(CO)(μ-CO){μ-η1:η2-C(Ph)CH(Ph)}], 2, was prepared in two steps from [Ru2Cp2(CO)2(μ-CO){μ-η1:η2-C(Ph)CH(Ph)}]BF4, [1]BF4, in 69% yield. Then, the reaction of 2 with C2(CO2Me)2, promoted by AgCF3SO3 in dichloromethane, afforded in 51% yield the complex [Ru2Cp2(CO)2{μ-η3:η2-C(Ph)CH(Ph)C(CO2Me)C(CO2Me)}]CF3SO3, [3]CF3SO3, containing a ruthenacyclopentene-based hydrocarbyl ligand. On the other hand, 2 reacted with other alkynes and AgX salts to give the butadienyl complexes [Ru2Cp2(CO)2{μ-η3:η2-C(R)CH(R')C(Ph)C(Ph)}]X (R = R' = H, [4]BF4; R = R' = Me, [5]CF3SO3; R = R' = Ph, [6]CF3SO3; R = Ph, R' = H, [7]CF3SO3), in 42-56% yields. All products were characterized by IR and NMR spectroscopy, and by single crystal X-ray diffraction in the cases of 2, [3]CF3SO3 and [6]BF4. DFT calculations highlighted the higher stability of [4-7]+-like structures with respect to the corresponding [3]+-like isomers. It is presumable that [3]+-like isomers initially form as kinetic intermediates, then undergo H-migration which is disfavoured in the presence of carboxylato substituents on the alkyne. Such hypothesis was supported by the computational optimization of the transition states for H-migration in the cases of R = R' = H and R = R' = CO2Me.
Collapse
Affiliation(s)
- Giulio Bresciani
- University of Pisa, Department of Chemistry and Industrial Chemistry, Via G. Moruzzi 13, I-56124 Pisa, Italy.
| | - Serena Boni
- University of Pisa, Department of Chemistry and Industrial Chemistry, Via G. Moruzzi 13, I-56124 Pisa, Italy.
| | - Stefano Zacchini
- University of Bologna, Department of Industrial Chemistry "Toso Montanari", Viale Risorgimento 4, I-40136 Bologna, Italy
| | - Guido Pampaloni
- University of Pisa, Department of Chemistry and Industrial Chemistry, Via G. Moruzzi 13, I-56124 Pisa, Italy.
| | - Marco Bortoluzzi
- University of Venezia "Ca' Foscari", Department of Molecular Science and Nanosystems, Via Torino 155, I-30170 Mestre (VE), Italy.
| | - Fabio Marchetti
- University of Pisa, Department of Chemistry and Industrial Chemistry, Via G. Moruzzi 13, I-56124 Pisa, Italy.
| |
Collapse
|
141
|
Viviano-Posadas AO, Romero-Mendoza U, Bazany-Rodríguez IJ, Velázquez-Castillo RV, Martínez-Otero D, Bautista-Renedo JM, González-Rivas N, Galindo-Murillo R, Salomón-Flores MK, Dorazco-González A. Efficient fluorescent recognition of ATP/GTP by a water-soluble bisquinolinium pyridine-2,6-dicarboxamide compound. Crystal structures, spectroscopic studies and interaction mode with DNA. RSC Adv 2022; 12:27826-27838. [PMID: 36320280 PMCID: PMC9520314 DOI: 10.1039/d2ra05040d] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2022] [Accepted: 09/21/2022] [Indexed: 11/04/2023] Open
Abstract
The new dicationic pyridine-2,6-dicarboxamide-based compound 1 bearing two N-alkylquinolinium units was synthesized, structurally determined by single-crystal X-ray diffraction, and studied in-depth as a fluorescent receptor for nucleotides and inorganic phosphorylated anions in pure water. The addition of nucleotides to 1 at pH = 7.0 quenches its blue emission with a selective affinity towards adenosine 5'-triphosphate (ATP) and guanosine 5'-tripohosphate (GTP) over other nucleotides such CTP, UTP, ADP, AMP, dicarboxylates and inorganic anions. On the basis of the spectroscopic tools (1H, 31P NMR, UV-vis, fluorescence), MS measurements and DFT calculations, receptor 1 binds ATP with high affinity (log K = 5.04) through the simultaneous formation of strong hydrogen bonds and π-π interactions between the adenosine fragment and quinolinium ring with binding energy calculated in 8.7 kcal mol-1. High affinity for ATP/GTP is attributed to the high acidity of amides and preorganized rigid structure of 1. Receptor 1 is an order of magnitude more selective for ATP than GTP. An efficient photoinduced electron transfer quenching mechanism with simultaneous receptor-ATP complexation in both the excited and ground states is proposed. Additionally, multiple spectroscopic studies and molecular dynamics simulations showed that 1 can intercalate into DNA base pairs.
Collapse
Affiliation(s)
- Alejandro O Viviano-Posadas
- Institute of Chemistry, National Autonomous University of Mexico Circuito Exterior, Ciudad Universitaria México, D.F., 04510 Mexico +52-55-56224514
| | - Ulises Romero-Mendoza
- Institute of Chemistry, National Autonomous University of Mexico Circuito Exterior, Ciudad Universitaria México, D.F., 04510 Mexico +52-55-56224514
| | - Iván J Bazany-Rodríguez
- Institute of Chemistry, National Autonomous University of Mexico Circuito Exterior, Ciudad Universitaria México, D.F., 04510 Mexico +52-55-56224514
| | - Rocío V Velázquez-Castillo
- Institute of Chemistry, National Autonomous University of Mexico Circuito Exterior, Ciudad Universitaria México, D.F., 04510 Mexico +52-55-56224514
| | - Diego Martínez-Otero
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM km 14.5 Carrera Toluca-Atlacomulco, Campus UAEMex "El Rosedal" San Cayetano-Toluca Toluca de Lerdo 50200 Estado de México Mexico
| | - Joanatan M Bautista-Renedo
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM km 14.5 Carrera Toluca-Atlacomulco, Campus UAEMex "El Rosedal" San Cayetano-Toluca Toluca de Lerdo 50200 Estado de México Mexico
| | - Nelly González-Rivas
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM km 14.5 Carrera Toluca-Atlacomulco, Campus UAEMex "El Rosedal" San Cayetano-Toluca Toluca de Lerdo 50200 Estado de México Mexico
| | - Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah 2000 East 30 South Skaggs 306 Salt Lake City 84112 Utah USA
| | - María K Salomón-Flores
- Institute of Chemistry, National Autonomous University of Mexico Circuito Exterior, Ciudad Universitaria México, D.F., 04510 Mexico +52-55-56224514
| | - Alejandro Dorazco-González
- Institute of Chemistry, National Autonomous University of Mexico Circuito Exterior, Ciudad Universitaria México, D.F., 04510 Mexico +52-55-56224514
| |
Collapse
|
142
|
Rummel L, König HF, Hausmann H, Schreiner PR. Silyl Groups Are Strong Dispersion Energy Donors. J Org Chem 2022; 87:13168-13177. [PMID: 36166406 DOI: 10.1021/acs.joc.2c01633] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present an experimental and computational study to investigate noncovalent interactions between silyl groups that are often employed as "innocent" protecting groups. We chose an extended cyclooctatetraene (COT)-based molecular balance comprising unfolded (1,4-disubstituted) and folded (1,6-disubstituted) valance bond isomers that typically display remote and close silyl group contacts, respectively. The thermodynamic equilibria were determined using nuclear magnetic resonance measurements. Additionally, we utilized Boltzmann weighted symmetry-adapted perturbation theory (SAPT) at the sSAPT0/aug-cc-pVDZ level of theory to dissect and quantify noncovalent interactions. Apart from the extremely bulky tris(trimethylsilyl)silyl "supersilyl" group, there is a preference for the folded 1,6-COT valence isomer, with London dispersion interactions being the main stabilizing factor. This makes silyl groups excellent dispersion energy donors, a finding that needs to be taken into account in synthesis planning.
Collapse
Affiliation(s)
- Lars Rummel
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| | - Henrik F König
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| | - Heike Hausmann
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| | - Peter R Schreiner
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| |
Collapse
|
143
|
Ethyl 5-Oxo-5-(((12-oxoindolo[2,1-b]quinazolin-6(12H)-ylidene)amino)oxy)pentanoate. MOLBANK 2022. [DOI: 10.3390/m1451] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Indolo[2,1-b]quinazolin-6,12-dione (tryptanthrin) derivatives present important types of nitrogen-containing heterocyclic compounds which are useful intermediate products in organic synthesis and have potential pharmaceutical applications. The new ethyl 5-oxo-5-(((12-oxoindolo[2,1-b]quinazolin-6(12H)-ylidene)amino)oxy)pentanoate (Compound 2) was synthesized. Compound 2 is the first example of a tryptanthrin derivative containing a dicarboxylic acid residue in the side chain. The Z,E-isomerism of Compound 2 was investigated by DFT calculations. Bioavailability was evaluated in silico using ADME predictions. According to the ADME results, Compound 2 is potentially highly bioavailable and has the prospective to be used as the main component for the development of anti-inflammatory drugs.
Collapse
|
144
|
Galleni L, Sajjadian FS, Conard T, Escudero D, Pourtois G, van Setten MJ. Modeling X-ray Photoelectron Spectroscopy of Macromolecules Using GW. J Phys Chem Lett 2022; 13:8666-8672. [PMID: 36084286 DOI: 10.1021/acs.jpclett.2c01935] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
We propose a simple additive approach to simulate X-ray photoelectron spectra (XPS) of macromolecules based on the GW method. Single-shot GW (G0W0) is a promising technique to compute accurate core-electron binding energies (BEs). However, its application to large molecules is still unfeasible. To circumvent the computational cost of G0W0, we break the macromolecule into tractable building blocks, such as isolated monomers, and sum up the theoretical spectra of each component, weighted by their molar ratio. In this work, we provide a first proof of concept by applying the method to four test polymers and one copolymer and show that it leads to an excellent agreement with experiments. The method could be used to retrieve the composition of unknown materials and study chemical reactions, by comparing the simulated spectra with experimental ones.
Collapse
Affiliation(s)
- Laura Galleni
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
- Imec, Kapeldreef 75, 3001 Leuven, Belgium
| | - Faegheh S Sajjadian
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
- Imec, Kapeldreef 75, 3001 Leuven, Belgium
| | | | - Daniel Escudero
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | | | | |
Collapse
|
145
|
Laqua H, Dietschreit JCB, Kussmann J, Ochsenfeld C. Accelerating Hybrid Density Functional Theory Molecular Dynamics Simulations by Seminumerical Integration, Resolution-of-the-Identity Approximation, and Graphics Processing Units. J Chem Theory Comput 2022; 18:6010-6020. [PMID: 36136665 DOI: 10.1021/acs.jctc.2c00509] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The computationally very demanding evaluation of the 4-center-2-electron (4c2e) integrals and their respective integral derivatives typically represents the major bottleneck within hybrid Kohn-Sham density functional theory molecular dynamics simulations. Building upon our previous works on seminumerical exact-exchange (sn-LinK) [Laqua, H., Thompsons, T. H., Kussmann, J., Ochsenfeld, C., J. Chem. Theory Comput. 2020, 16, 1465] and resolution-of-the-identity Coulomb (RI-J) [Kussmann, J., Laqua, H., Ochsenfeld, C., J. Chem. Theory Comput. 2021, 17, 1512], the expensive 4c2e integral evaluation can be avoided entirely, resulting in a highly efficient electronic structure theory method, allowing for fast ab initio molecular dynamics (AIMD) simulations even with large basis sets. Moreover, we propose to combine the final self-consistent field (SCF) step with the subsequent nuclear forces evaluation, providing the forces at virtually no additional cost after a converged SCF calculation, reducing the total runtime of an AIMD simulation by about another 25%. In addition, multiple independent MD trajectories can be computed concurrently on a single node, leading to a greatly increased utilization of the available hardware─especially when combined with graphics processing unit acceleration─improving the overall throughput by up to another 5 times in this way. With all of those optimizations combined, our proposed method provides nearly 3 orders of magnitude faster execution times than traditional 4c2e integral-based methods. To demonstrate the practical utility of the approach, quantum-mechanical/molecular-mechanical dynamics simulations on double-stranded DNA were performed, investigating the relative hydrogen bond strength between adenine-thymine and guanine-cytosine base pairs. In addition, this illustrative application also contains a general accuracy assessment of the introduced approximations (integration grids, resolution-of-the-identity) within AIMD simulations, serving as a protocol on how to apply these new methods to practical problems.
Collapse
Affiliation(s)
- Henryk Laqua
- Department of Chemistry, Chair of Theoretical Chemistry, University of Munich (LMU), D-81377 München, Germany
| | - Johannes C B Dietschreit
- Department of Chemistry, Chair of Theoretical Chemistry, University of Munich (LMU), D-81377 München, Germany
| | - Jörg Kussmann
- Department of Chemistry, Chair of Theoretical Chemistry, University of Munich (LMU), D-81377 München, Germany
| | - Christian Ochsenfeld
- Department of Chemistry, Chair of Theoretical Chemistry, University of Munich (LMU), D-81377 München, Germany.,Max Planck Institute for Solid State Research, 70569 Stuttgart, Germany
| |
Collapse
|
146
|
Bursch M, Mewes J, Hansen A, Grimme S. Best‐Practice DFT Protocols for Basic Molecular Computational Chemistry**. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202205735] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 45470 Mülheim an der Ruhr Germany
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| |
Collapse
|
147
|
Cesari C, Bortoluzzi M, Forti F, Gubbels L, Femoni C, Iapalucci MC, Zacchini S. 2-D Molecular Alloy Ru-M (M = Cu, Ag, and Au) Carbonyl Clusters: Synthesis, Molecular Structure, Catalysis, and Computational Studies. Inorg Chem 2022; 61:14726-14741. [PMID: 36069711 PMCID: PMC9490753 DOI: 10.1021/acs.inorgchem.2c02099] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The reactions of
[HRu3(CO)11]− (1) with M(I) (M = Cu, Ag, and Au) compounds such as
[Cu(CH3CN)4][BF4], AgNO3, and Au(Et2S)Cl afford the 2-D molecular alloy clusters
[CuRu6(CO)22]− (2), [AgRu6(CO)22]− (3), and [AuRu5(CO)19]− (4), respectively. The reactions of 2–4 with PPh3 result in mixtures of products, among which
[Cu2Ru8(CO)26]2– (5), Ru4(CO)12(CuPPh3)4 (6), Ru4(CO)12(AgPPh3)4 (7), Ru(CO)3(PPh3)2 (8), and HRu3(OH)(CO)7(PPh3)3 (9) have been isolated
and characterized. The molecular structures of 2–6 and 9 have been determined by single-crystal X-ray
diffraction. The metal–metal bonding within 2–5 has been computationally investigated by density functional theory
methods. In addition, the [NEt4]+ salts of 2–4 have been tested as catalyst precursors for transfer
hydrogenation on the model substrate 4-fluoroacetophenone using iPrOH as a solvent and a hydrogen source. Heterometallic Ru−M (M = Cu, Ag,
and Au) carbonyl
clusters possessing a 2-D metal core have been investigated by structural
and computational methods and tested as catalyst precursors for transfer
hydrogenation.
Collapse
Affiliation(s)
- Cristiana Cesari
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale Risorgimento 4, 40136 Bologna, Italy.,Center for Chemical Catalysis─C3, University of Bologna, Viale Risorgimento 4, 40136 Bologna, Italy
| | - Marco Bortoluzzi
- Dipartimento di Scienze Molecolari e Nanosistemi, Ca' Foscari University of Venice, Via Torino 155, 30175 Mestre (Ve), Italy
| | - Francesca Forti
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale Risorgimento 4, 40136 Bologna, Italy.,Center for Chemical Catalysis─C3, University of Bologna, Viale Risorgimento 4, 40136 Bologna, Italy
| | - Lisa Gubbels
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale Risorgimento 4, 40136 Bologna, Italy
| | - Cristina Femoni
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale Risorgimento 4, 40136 Bologna, Italy
| | - Maria Carmela Iapalucci
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale Risorgimento 4, 40136 Bologna, Italy
| | - Stefano Zacchini
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale Risorgimento 4, 40136 Bologna, Italy.,Center for Chemical Catalysis─C3, University of Bologna, Viale Risorgimento 4, 40136 Bologna, Italy
| |
Collapse
|
148
|
The Role of the Bridge in Single-Ion Magnet Behaviour: Reinvestigation of Cobalt(II) Succinate and Fumarate Coordination Polymers with Nicotinamide. INORGANICS 2022. [DOI: 10.3390/inorganics10090128] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Two previously synthesized cobalt(II) coordination polymers; {[Co(μ2-suc)(nia)2(H2O)2]·2H2O}n (suc = succinate(2−), nia = nicotinamide) and [Co(μ2-fum)(nia)2(H2O)2]n (fum = fumarate(2−)) were prepared and thoroughly characterized. Both complexes form 1D coordination chains by bonding of Co(nia)2(H2O)2 units through succinate or fumarate ligands while these chains are further linked through hydrogen bonds to 3D supramolecular networks. The intermolecular interactions of both complexes are quantified using Hirshfeld surface analysis and their infrared spectra, electronic spectra and static magnetic properties are confronted with DFT and state-of-the-art ab-initio calculations. Dynamic magnetic measurements show that both complexes exhibit single-ion magnet behaviour induced by a magnetic field. Since they possess very similar chemical structure, differing only in the rigidity of the bridge between the magnetic centres, this chemical feature is put into context with changes in their magnetic relaxation.
Collapse
|
149
|
Cheng F, Wu S, Zheng W, Su S, Nakanishi T, Xu W, Sadhukhan P, Sejima H, Ikenaga S, Yamamoto K, Gao K, Kanegawa S, Sato O. Macroscopic Polarization Change of Mononuclear Valence Tautomeric Cobalt Complexes Through the Use of Enantiopure Ligand. Chemistry 2022; 28:e202202161. [DOI: 10.1002/chem.202202161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Indexed: 11/08/2022]
Affiliation(s)
- Feng Cheng
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Shuqi Wu
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Wenwei Zheng
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Shengqun Su
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Takumi Nakanishi
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Wenhuang Xu
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Pritam Sadhukhan
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Hibiki Sejima
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Shimon Ikenaga
- Department of Physics Okayama University of Science Okayama Japan
| | - Kaoru Yamamoto
- Department of Physics Okayama University of Science Okayama Japan
| | - Kaige Gao
- College of Physical Science and Technology Yangzhou University Yangzhou Jiangsu P. R. China
| | - Shinji Kanegawa
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| | - Osamu Sato
- Institute for Materials Chemistry and Engineering & IRCCS Kyushu University 744 Motooka Nishi-ku Fukuoka 819-0395 Japan
| |
Collapse
|
150
|
Chan B. High-Level Quantum Chemistry Reference Heats of Formation for a Large Set of C, H, N, and O Species in the NIST Chemistry Webbook and the Identification and Validation of Reliable Protocols for Their Rapid Computation. J Phys Chem A 2022; 126:4981-4990. [PMID: 35878062 DOI: 10.1021/acs.jpca.2c03846] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
A recent study has examined the accuracy of NIST heats of formation for a set of C, H, and O-containing species with a proposed low-cost quantum chemistry approach. In the present study, we have used high-level methods such as W1X-2 to obtain these data more rigorously, which we have then used to assess the NIST and the previously computed values. We find that many of these NIST data that are as suggested to be unreliable by the previous study are indeed inconsistent with our high-level reference values. However, we also find substantial deviations for the previously computed values from our benchmark. Thus, we have assessed the performance of alternative low-cost methods. In our assessment, we have additionally examined C, H, N, and O-containing species for which heats of formation are available from the NIST database. We find the ωB97M-V/ma-def2-TZVP, DSD-PBEP86/ma-def2-TZVP, and CCSD(T)-F12b/aug'-cc-pVDZ methods to be adequate for obtaining heats of formation with the atomization approach, once their atomic energies are optimized with our benchmark. Notably, the low-cost ωB97M-V method yields values that agree to be within 10 kJ mol-1 for more than 90% of the (∼1500) species. A higher 20 kJ mol-1 threshold captures 98% of the data. The outlier species typically contain many electron-withdrawing (nitro) groups. In these cases, the use of isodesmic-type reactions rather than the atomization approach is more reliable. Our assessment has also identified significant outliers from the NIST database, for which experimental re-determination of the heats of formation would be desirable.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| |
Collapse
|