1
|
La Force H, Freindorf M, Kraka E. Ligand Characterization and DNA Intercalation of Ru(II) Polypyridyl Complexes: A Local Vibrational Mode Study. J Phys Chem A 2024. [PMID: 38990174 DOI: 10.1021/acs.jpca.4c02954] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/12/2024]
Abstract
We investigated in this work ruthenium-ligand bonding across the RuN framework in 12 Ru(II) polypyridyl complexes in the gas phase and solution for both singlet and triplet states, in addition to their affinity for DNA binding through π-π stacking interactions with DNA nucleobases. As a tool to assess the intrinsic strength of the ruthenium-ligand bonds, we determined local vibrational force constants via our local vibrational mode analysis software. We introduced a novel local force constant that directly accounts for the intrinsic strength of the π-π stacking interaction between DNA and the intercalated Ru(II) complex. According to our findings, [Ru(phen)2(dppz)]2+ and [Ru(phen)2(11-CN-dppz)]2+ provide an intriguing trade-off between photoinduced complex excitation and the strength of the subsequent π-π stacking interaction with DNA. [Ru(phen)2(dppz)]2+ displays a small singlet-triplet splitting and a strong π-π stacking interaction in its singlet state, suggesting a favorable photoexcitation but potentially weaker interaction with DNA in the excited state. Conversely, [Ru(phen)2(11-CN-dppz)]2+ exhibits a larger singlet-triplet splitting and a stronger π-π stacking interaction with DNA in its triplet state, indicating a less favorable photoinduced transition but a stronger interaction with DNA postexcitation. We hope our study will inspire future experimental and computational work aimed at the design of novel Ru-polypyridyl drug candidates and that our new quantitative measure of π-π stacking interactions in DNA will find a general application in the field.
Collapse
Affiliation(s)
- Hunter La Force
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
2
|
Patel P, Chung J, Bowman MA, Ulusoy I, Wilson AK. Potential energy surfaces and dynamic properties via ab initio composite and density functional approaches. J Comput Chem 2024; 45:1352-1363. [PMID: 38376255 DOI: 10.1002/jcc.27333] [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: 09/19/2023] [Revised: 02/01/2024] [Accepted: 02/02/2024] [Indexed: 02/21/2024]
Abstract
Vibrational spectroscopy enables critical insight into the structural and dynamic properties of molecules. Presently, the majority of theoretical approaches to spectroscopy employ wavefunction-based ab initio or density functional methods that rely on the harmonic approximation. This approximation breaks down for large molecules with strongly anharmonic bonds or for molecules with large internuclear separations. An alternative to these methods involves generating molecular anharmonic potential energy surfaces (potentials) and using them to extrapolate the vibrational frequencies. This study examines the efficacy of density functional theory (DFT) and the correlation consistent Composite Approach (ccCA) in generating anharmonic frequencies from potentials of small main group molecules. Vibrational self-consistent field Theory (VSCF) and post-VSCF methods were used to calculate the fundamental frequencies of these molecules from their potentials. Functional choice, basis set selection, and mode-coupling are also examined as factors in influencing accuracy. The absolute deviations for the calculated frequencies using potentials at the ccCA level of theory were lower than the potentials at the DFT level. With DFT resulting in bending modes that are better described than those of ccCA, a multilevel DFT:ccCA approach where DFT potentials are used for single vibrational mode potentials and ccCA is used for vibrational mode-mode couplings can be utilized for larger polyatomic systems. The frequencies obtained with this multilevel approach using VCIPSI-PT2 were closer to experimental frequencies than the scaled harmonic frequencies, indicating the success of utilizing post-VSCF methods to generate more accurate representations of computed infrared spectra.
Collapse
Affiliation(s)
- Prajay Patel
- Department of Chemistry, Michigan State University, East Lansing, Michigan, USA
- Chemistry Department, University of Dallas, Irving, Texas, USA
| | - Joseph Chung
- Department of Chemistry, Michigan State University, East Lansing, Michigan, USA
| | - Max Aksel Bowman
- Department of Chemistry, Michigan State University, East Lansing, Michigan, USA
| | - Inga Ulusoy
- Department of Chemistry, Michigan State University, East Lansing, Michigan, USA
- Scientific Software Center, Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan, USA
| |
Collapse
|
3
|
Bodo F, Erba A, Kraka E, Moura RT. Chemical bonding in Uranium-based materials: A local vibrational mode case study of Cs 2 UO 2 Cl 4 and UCl 4 crystals. J Comput Chem 2024; 45:1130-1142. [PMID: 38279637 DOI: 10.1002/jcc.27311] [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: 09/29/2023] [Revised: 12/10/2023] [Accepted: 01/08/2024] [Indexed: 01/28/2024]
Abstract
The Local Vibrational Mode Analysis, initially applied to diverse molecular systems, was extended to periodic systems in 2019. This work introduces an enhanced version of the LModeA software, specifically designed for the comprehensive analysis of two and three-dimensional periodic structures. Notably, a novel interface with the Crystal package was established, enabling a seamless transition from molecules to periodic systems using a unified methodology. Two distinct sets of uranium-based systems were investigated: (i) the evolution of the Uranyl ion (UO 2 2 + ) traced from its molecular configurations to the solid state, exemplified by Cs 2 UO 2 Cl 4 and (ii) Uranium tetrachloride (UCl 4 ) in both its molecular and crystalline forms. The primary focus was on exploring the impact of crystal packing on key properties, including IR and Raman spectra, structural parameters, and an in-depth assessment of bond strength utilizing local mode perspectives. This work not only demonstrates the adaptability and versatility of LModeA for periodic systems but also highlights its potential for gaining insights into complex materials and aiding in the design of new materials through fine-tuning.
Collapse
Affiliation(s)
- Filippo Bodo
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
- Dipartimento di Chimica, Università di Torino, Torino, Italy
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, Torino, Italy
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
| | - Renaldo T Moura
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, Brazil
| |
Collapse
|
4
|
Filatov M, Mironov V, Kraka E. Unraveling the effect of aromaticity for the dynamics of excited states of single benzene fluorophores. J Comput Chem 2024; 45:1033-1045. [PMID: 38216513 DOI: 10.1002/jcc.27304] [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: 09/29/2023] [Revised: 11/23/2023] [Accepted: 12/23/2023] [Indexed: 01/14/2024]
Abstract
The photophysical properties of a series of recently synthesized single benzene fluorophores were investigated using ensemble density functional theory calculations. The energetic stability of the ground and excited state species were counterposed against the aromaticity index derived from local vibrational modes. It was found that the large Stokes shift of the fluorophores (up to ca. 5800 cm - 1 ) originates from the effect of electron donating and electron withdrawing substituents rather than π -delocalization and related (anti-)aromaticity. On the basis of nonadiabatic molecular dynamics simulations, the absence of fluorescence from one of the regioisomers was explained by the occurrence of easily accessible S 1 /S 0 conical intersections below the vertical excitation energy level. It is demonstrated in the manuscript that the analysis of local mode force constants and the related aromaticity index represent a useful tool for the characterization of π -delocalization effects in π -conjugated compounds.
Collapse
Affiliation(s)
- Michael Filatov
- Center for Multidimensional Carbon Materials, Institute for Basic Science (IBS), Ulsan, Republic of Korea
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
| | | | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
| |
Collapse
|
5
|
Chagas JCV, Milanez BD, Oliveira VP, Pinheiro M, Ferrão LFA, Aquino AJA, Lischka H, Machado FBC. A multi-descriptor analysis of substituent effects on the structure and aromaticity of benzene derivatives: π-Conjugation versus charge effects. J Comput Chem 2024; 45:863-877. [PMID: 38153839 DOI: 10.1002/jcc.27296] [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: 08/24/2023] [Revised: 11/29/2023] [Accepted: 12/10/2023] [Indexed: 12/30/2023]
Abstract
This work provides a detailed multi-component analysis of aromaticity in monosubstituted (X = CH3, CH 2 - , CH 2 + , NH2, NH-, NH+, OH, O-, and O+) and para-homodisubstituted (X = CH3, CH2, NH2, NH, OH, and O) benzene derivatives. We investigate the effects of substituents using single-reference (B3LYP/DFT) and multireference (CASSCF/MRCI) methods, focusing on structural (HOMA), vibrational (AI(vib)), topological (ELFπ), electronic (MCI), magnetic (NICS), and stability (S0-T1 splitting) properties. The findings reveal that appropriate π-electron-donating and π-electron-accepting substituents with suitable size and symmetry can interact with the π-system of the ring, significantly influencing π-electron delocalization. While the charge factor has a minimal impact on π-electron delocalization, the presence of a pz orbital capable of interacting with the π-electron delocalization is the primary factor leading to a deviation from the typical aromaticity characteristics observed in benzene.
Collapse
Affiliation(s)
- Julio C V Chagas
- Department of Chemistry, Aeronautics Institute of Technology, São José dos Campos, São Paulo, Brazil
| | - Bruno D Milanez
- Department of Chemistry, Aeronautics Institute of Technology, São José dos Campos, São Paulo, Brazil
| | - Vytor P Oliveira
- Department of Chemistry, Aeronautics Institute of Technology, São José dos Campos, São Paulo, Brazil
| | - Max Pinheiro
- Department of Chemistry, Aeronautics Institute of Technology, São José dos Campos, São Paulo, Brazil
| | - Luiz F A Ferrão
- Department of Chemistry, Aeronautics Institute of Technology, São José dos Campos, São Paulo, Brazil
| | - Adelia J A Aquino
- Department of Mechanical Engineering, Texas Tech University, Lubbock, Texas, USA
| | - Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas, USA
| | - Francisco B C Machado
- Department of Chemistry, Aeronautics Institute of Technology, São José dos Campos, São Paulo, Brazil
| |
Collapse
|
6
|
Freindorf M, Antonio JJ, Kraka E. Iron-histidine bonding in bishistidyl hemoproteins-A local vibrational mode study. J Comput Chem 2024; 45:574-588. [PMID: 38041830 DOI: 10.1002/jcc.27267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2023] [Revised: 11/09/2023] [Accepted: 11/11/2023] [Indexed: 12/04/2023]
Abstract
We investigated the intrinsic strength of distal and proximal FeN bonds for both ferric and ferrous oxidation states of bishistidyl hemoproteins from bacteria, animals, human, and plants, including two cytoglobins, ten hemoglobins, two myoglobins, six neuroglobins, and six phytoglobins. As a qualified measure of bond strength, we used local vibrational force constants ka (FeN) based on local mode theory developed in our group. All calculations were performed with a hybrid QM/MM ansatz. Starting geometries were taken from available x-ray structures. ka (FeN) values were correlated with FeN bond lengths and covalent bond character. We also investigated the stiffness of the axial NFeN bond angle. Our results highlight that protein effects are sensitively reflected in ka (FeN), allowing one to compare trends in diverse protein groups. Moreover, ka (NFeN) is a perfect tool to monitor changes in the axial heme framework caused by different protein environments as well as different Fe oxidation states.
Collapse
Affiliation(s)
- Marek Freindorf
- Chemistry Department, Southern Methodist University, Dallas, Texas, USA
| | - Juliana J Antonio
- Chemistry Department, Southern Methodist University, Dallas, Texas, USA
| | - Elfi Kraka
- Chemistry Department, Southern Methodist University, Dallas, Texas, USA
| |
Collapse
|
7
|
Benítez FJ, Gutiérrez-Oliva S, Herrera B, Toro-Labbé A. Basis Electronic Activity of Molecular Systems. A Theory of Bond Reactivity. J Phys Chem A 2024. [PMID: 38437616 DOI: 10.1021/acs.jpca.4c00359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2024]
Abstract
In this paper, we present a new finding, the basis electronic activity (BEA) of molecular systems; it corresponds to the significant, although nonreactive, vibrationally induced electronic activity that takes place in any molecular system. Although the molecule's BEA is composed of an equal number of local contributions as the vibrational degrees of freedom, our results indicate that only stretching modes contribute to it. To account for this electronic activity, a new descriptor, the bond electronic flux (BEF), is introduced. The BEF combined with the force constant of the potential well hosting the electronic activity gives rise to the effective bond reactivity index (EBR), which turns out to be the first density functional theory-based descriptor that simultaneously accounts for structural and electronic effects. Besides quantifying the bond reactivity, EBR provides a basis to compare the reactivities of bonds inserted in different chemical environments and paves the way for the exertion of selective control to enhance or inhibit their reactivities. The new concepts formulated in this paper and the associated computational tools are illustrated with characterization of the BEA of a set of representative molecules. In all cases, the BEFs follow the same linear pattern, whose slopes indicate the intensity of the electronic activity and quantify the reactivity of chemical bonds.
Collapse
Affiliation(s)
- Francisca J Benítez
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Avda. Vicuña Mackenna 4860 Macul, Santiago 7820436, Chile
| | - Soledad Gutiérrez-Oliva
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Avda. Vicuña Mackenna 4860 Macul, Santiago 7820436, Chile
| | - Bárbara Herrera
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Avda. Vicuña Mackenna 4860 Macul, Santiago 7820436, Chile
| | - Alejandro Toro-Labbé
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Avda. Vicuña Mackenna 4860 Macul, Santiago 7820436, Chile
| |
Collapse
|
8
|
Devore DP, Ellington TL, Shuford KL. Elucidating the Role of Electron-Donating Groups in Halogen Bonding. J Phys Chem A 2024; 128:1477-1490. [PMID: 38373286 DOI: 10.1021/acs.jpca.3c06894] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
Computational quantum chemical techniques were utilized to systematically examine how electron-donating groups affect the electronic and spectroscopic properties of halogen bond donors and their corresponding complexes. Unlike the majority of studies on halogen bonding, where electron-withdrawing groups are utilized, this work investigates the influence of electron-donating substituents within the halogen bond donors. Statistical analyses were performed on the descriptors of halogen bond donors in a prescribed set of archetype, halo-alkyne, halo-benzene, and halo-ethynyl benzene halogen bond systems. The σ-hole magnitude, binding and interaction energies, and the vibrational X···N local force constant (where X = Cl, Br, I, and At) were found to correlate very well in a monotonic and linear manner with all other properties studied. In addition, enhanced halogen bonds were found when the systems contained electron-donating groups that could form intramolecular hydrogen bonds with the electronegative belt of the halogen atom and adjacent linker features.
Collapse
Affiliation(s)
- Daniel P Devore
- Department of Chemistry and Biochemistry, Baylor University, One Bear Place #97348, Waco, Texas 76798-7348, United States
| | - Thomas L Ellington
- Department of Chemistry and Biochemistry, Baylor University, One Bear Place #97348, Waco, Texas 76798-7348, United States
| | - Kevin L Shuford
- Department of Chemistry and Biochemistry, Baylor University, One Bear Place #97348, Waco, Texas 76798-7348, United States
| |
Collapse
|
9
|
Streu K, Hunsberger S, Patel J, Wan X, Daly CA. Development of a universal method for vibrational analysis of the terminal alkyne C≡C stretch. J Chem Phys 2024; 160:074106. [PMID: 38364010 DOI: 10.1063/5.0185580] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 01/07/2024] [Indexed: 02/18/2024] Open
Abstract
The terminal alkyne C≡C stretch has a large Raman scattering cross section in the "silent" region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born-Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.
Collapse
Affiliation(s)
- Kristina Streu
- Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, USA
| | - Sara Hunsberger
- Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, USA
| | - Jeanette Patel
- Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, USA
| | - Xiang Wan
- Department of Mathematics and Statistics, Loyola University Chicago, 1032 W. Sheridan Rd., Chicago, Illinois 60660, USA
| | - Clyde A Daly
- Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, USA
| |
Collapse
|
10
|
Dangat Y, Freindorf M, Kraka E. Mechanistic Insights into S-Depalmitolyse Activity of Cln5 Protein Linked to Neurodegeneration and Batten Disease: A QM/MM Study. J Am Chem Soc 2024; 146:145-158. [PMID: 38055807 DOI: 10.1021/jacs.3c06397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/08/2023]
Abstract
Ceroid lipofuscinosis neuronal protein 5 (Cln5) is encoded by the CLN5 gene. The genetic variants of this gene are associated with the CLN5 form of Batten disease. Recently, the first crystal structure of Cln5 was reported. Cln5 shows cysteine palmitoyl thioesterase S-depalmitoylation activity, which was explored via fluorescent emission spectroscopy utilizing the fluorescent probe DDP-5. In this work, the mechanism of the reaction between Cln5 and DDP-5 was studied computationally by applying a QM/MM methodology at the ωB97X-D/6-31G(d,p):AMBER level. The results of our study clearly demonstrate the critical role of the catalytic triad Cys280-His166-Glu183 in S-depalmitoylation activity. This is evidenced through a comparison of the pathways catalyzed by the Cys280-His166-Glu183 triad and those with only Cys280 involved. The computed reaction barriers are in agreement with the catalytic efficiency. The calculated Gibb's free-energy profile suggests that S-depalmitoylation is a rate-limiting step compared to the preceding S-palmitoylation, with barriers of 26.1 and 25.3 kcal/mol, respectively. The energetics were complemented by monitoring the fluctuations in the electron density distribution through NBO charges and bond strength alterations via local mode stretching force constants during the catalytic pathways. This comprehensive protocol led to a more holistic picture of the reaction mechanism at the atomic level. It forms the foundation for future studies on the effects of gene mutations on both the S-palmitoylation and S-depalmitoylation steps, providing valuable data for the further development of enzyme replacement therapy, which is currently the only FDA-approved therapy for childhood neurodegenerative diseases, including Batten disease.
Collapse
Affiliation(s)
- Yuvraj Dangat
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Marek Freindorf
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
11
|
Tran HK, Berkelbach TC. Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals. J Chem Phys 2023; 159:194101. [PMID: 37965997 PMCID: PMC10653875 DOI: 10.1063/5.0172702] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Accepted: 10/23/2023] [Indexed: 11/16/2023] Open
Abstract
Vibrational heat-bath configuration interaction (VHCI)-a selected configuration interaction technique for vibrational structure theory-has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm-1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.
Collapse
Affiliation(s)
- Henry K. Tran
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
12
|
Phan Dang CT, Tam NM, Huynh TN, Trung NT. Revisiting conventional noncovalent interactions towards a complete understanding: from tetrel to pnicogen, chalcogen, and halogen bond. RSC Adv 2023; 13:31507-31517. [PMID: 37901266 PMCID: PMC10606978 DOI: 10.1039/d3ra06078k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 10/23/2023] [Indexed: 10/31/2023] Open
Abstract
Typical noncovalent interactions, including tetrel (TtB), pnicogen (PniB), chalcogen (ChalB), and halogen bonds (HalB), were systematically re-investigated by modeling the N⋯Z interactions (Z = Si, P, S, Cl) between NH3 - as a nucleophilic, and SiF4, PF3, SF2, and ClF - as electrophilic components, employing highly reliable ab initio methods. The characteristics of N⋯Z interactions when Z goes from Si to Cl, were examined through their changes in stability, vibrational spectroscopy, electron density, and natural orbital analyses. The binding energies of these complexes at CCSD(T)/CBS indicate that NH3 tends to hold tightly most with ClF (-34.7 kJ mol-1) and SiF4 (-23.7 kJ mol-1) to form N⋯Cl HalB and N⋯Si TtB, respectively. Remarkably, the interaction energies obtained from various approaches imply that the strength of these noncovalent interactions follows the order: N⋯Si TtB > N⋯Cl HalB > N⋯S ChalB > N⋯P PniB, that differs the order of their corresponding complex stability. The conventional N⋯Z noncovalent interactions are characterized by the local vibrational frequencies of 351, 126, 167, and 261 cm-1 for TtB, PniB, ChalB, and HalB, respectively. The SAPT2+(3)dMP2 calculations demonstrate that the primary force controlling their strength retains the electrostatic term. Accompanied by the stronger strength of N⋯Si TtB and N⋯Cl HalB, the AIM and NBO results state that they are partly covalent in nature with amounts of 18.57% and 27.53%, respectively. Among various analysis approaches, the force constant of the local N⋯Z stretching vibration is shown to be most accurate in describing the noncovalent interactions.
Collapse
Affiliation(s)
- Cam-Tu Phan Dang
- Faculty of Natural Sciences, Duy Tan University Da Nang 550000 Vietnam
- Institute of Research and Development, Duy Tan University Da Nang 550000 Vietnam
| | - Nguyen Minh Tam
- Faculty of Basic Sciences, University of Phan Thiet 225 Nguyen Thong Phan Thiet City Binh Thuan Vietnam
| | - Thanh-Nam Huynh
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology Eggenstein-Leopoldshafen 76344 Germany
| | - Nguyen Tien Trung
- Laboratory of Computational Chemistry and Modelling (LCCM), Quy Nhon University Quy Nhon City 590000 Vietnam
| |
Collapse
|
13
|
Freindorf M, Antonio J, Kraka E. Hydrogen Sulfide Ligation in Hemoglobin I of Lucina pectinata─A QM/MM and Local Mode Study. J Phys Chem A 2023; 127:8316-8329. [PMID: 37774120 DOI: 10.1021/acs.jpca.3c04399] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/01/2023]
Abstract
In this study, we investigated the interaction between the H2S ligand and the heme pocket of hemoglobin I (HbI) of Lucina pectinata for the wild-type protein; three known mutations where distal glutamine is replaced by hydrophobic valine (Gln64Val) and hydrophilic histidine in both protonation forms (Gln64Hisϵ and Gln64Hisδ); five known mutations of the so-called phenyl cage, replacing the hydrophobic phenylalanines Phe29 and Phe43 with tyrosine (Tyr), valine (Val), or leucine (Leu); and two additional mutations, Phe68Tyr and Phe68Val, in order to complement previous studies with new insights about the binding mechanism at the molecular level. A particular focus was on the intrinsic strengths of the chemical bonds involved, utilizing local vibrational force constants based on combined quantum mechanical-molecular mechanical calculations. Wild-type protein and mutations clustered into two distinct groups: Group 1 protein systems with a proton acceptor in the distal protein pocket, close to one of the H2S bonds, and Group 2 protein systems without a hydrogen acceptor close by in the active site of the protein. According to our results, the interactions between H2S and HbI of Lucina pectinata involve two important elements, namely, binding of H2S to Fe of the heme group, followed by the proton transfer from the HS bond to the distal residue. The distal residue is additionally stabilized by a second proton transfer from the distal residue to COO- of the propionate group in heme. We could identify the FeS bond as a key player and discovered that the strength of this bond depends on two mutual factors, namely, the strength of the HS bond involved in the proton transfer and the electrostatic field of the protein pocket qualifying the FeS bond as a sensitive probe for monitoring changes in H2S ligation upon protein mutations. We hope our study will inspire and guide future experimental studies, targeting new promising mutations such as Phe68Tyr, Phe68Val, or Phe43Tyr/Phe68Val.
Collapse
Affiliation(s)
- Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Juliana Antonio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
14
|
Santos CV, Monteiro SA, Soares ASC, Souto ICA, Moura RT. Decoding Chemical Bonds: Assessment of the Basis Set Effect on Overlap Electron Density Descriptors and Topological Properties in Comparison to QTAIM. J Phys Chem A 2023; 127:7997-8014. [PMID: 37703453 DOI: 10.1021/acs.jpca.3c04504] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2023]
Abstract
Quantum chemical bonding descriptors based on the total and overlap density can provide valuable information about chemical interactions in different systems. However, these descriptors can be sensitive to the basis set used. To address this, different numerical treatments of electron density have been proposed to reduce the basis set dependency. In this work, we introduce overlap properties (OPs) obtained through numerical treatment of the electron density and present the topology of overlap density (TOP) for the first time. We compare the basis set dependency of numerical OP and TOP descriptors with their quantum theory of atoms in molecules (QTAIM) counterparts, considering the total electron density. Three single (C-C, C-O, and C-F) bonds in ethane, methanol, and fluoromethane and two double (C═C and C═O) bonds in ethene and formaldehyde were analyzed. Diatomic molecules Li-X with X = F, Cl, and Br were also analyzed. Eight parameters, including QTAIM descriptors and OP/TOP descriptors, are used to assess the basis dependency at the ωB97X-D level of theory using 28 basis sets from three classes: Pople, Ahlrichs, and Dunning. The study revealed that the topological overlap electron density properties exhibit comparatively lesser dependence on the basis set compared to their total electron density counterparts. Remarkably, these properties retain their chemical significance even with reduced basis set dependency. Similarly, numerical OP descriptors show less basis set dependency than their QTAIM counterparts. The excess of polarization functions increases charge concentration in the interatomic region and influences both QTAIM and OP descriptors. The basis sets Def2TZVP, 6-31++G(d,p), 6-311++G(d,p), cc-pVDZ, cc-pVTZ, and cc-pVQZ demonstrate reduced variability for the tested bond classes in this study, with particular emphasis on the triple-ζ quality Ahlrichs' basis set. We recommend against using basis sets with numerous polarization functions, such as augmented Dunning's and Ahlrichs' quadruple-ζ.
Collapse
Affiliation(s)
- Carlos V Santos
- Department of Chemistry, Federal University of Paraiba, Joao Pessoa, Paraiba 58051-970, Brazil
| | - Shirlene A Monteiro
- Department of Chemistry, State University of Paraiba, Campina Grande, Paraiba 58051-970, Brazil
| | - Amanda S C Soares
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, Paraiba 58397-000, Brazil
| | - Isabeli C A Souto
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, Paraiba 58397-000, Brazil
| | - Renaldo T Moura
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, Paraiba 58397-000, Brazil
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| |
Collapse
|
15
|
Warren LR, Rowell A, McMaster P, Pulham CR, Morrison CA. Accounting for super-, plateau- and mesa-rate burning by lead and copper-based ballistic modifiers in double-base propellants: a computational study. Phys Chem Chem Phys 2023; 25:23995-24003. [PMID: 37646274 DOI: 10.1039/d3cp03197g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
Abstract
We present a first-principles computational study to understand the action of lead and copper-based ballistic modifiers in the combustion of double-base propellants (DBPs). We show that lead oxide clusters are easily broken down upon addition of small amounts of carbon and the resulting graphitic matrix, dispersed with weakly bound and exposed Pb sites, acts as a Lewis acid to bind small molecule Lewis bases such as NO2 and CH2O that form in the combustion flame. This accounts for super-rate burning, where the fuel burn rate is enhanced. We also show how carbon availability accounts for the plateau- and mesa-rate burning effects, where the fuel burn rate is suppressed. In contrast, cluster integrity on binding carbon to copper oxide is retained, and interaction with NO2 and CH2O is essentially negligible. Carbon binds more strongly to copper oxide, however, and we therefore propose that when carbon levels start to fall this results in the lead oxide clusters being starved of carbon, which leads to plateau and mesa burning. Taken together, the calculations support a general model that accounts for the super-, plateau- and mesa-rate ballistic modifier burning effects.
Collapse
Affiliation(s)
- Lisette R Warren
- EaSTCHEM School of Chemistry, University of Edinburgh, The King's Buildings, David Brewster Road, Edinburgh, EH9 3FJ, UK.
| | - Aaron Rowell
- EaSTCHEM School of Chemistry, University of Edinburgh, The King's Buildings, David Brewster Road, Edinburgh, EH9 3FJ, UK.
| | | | - Colin R Pulham
- EaSTCHEM School of Chemistry, University of Edinburgh, The King's Buildings, David Brewster Road, Edinburgh, EH9 3FJ, UK.
| | - Carole A Morrison
- EaSTCHEM School of Chemistry, University of Edinburgh, The King's Buildings, David Brewster Road, Edinburgh, EH9 3FJ, UK.
| |
Collapse
|
16
|
Antonio JJ, Kraka E. Noncovalent π Interactions in Mutated Aquomet-Myoglobin Proteins: A QM/MM and Local Vibrational Mode Study. Biochemistry 2023; 62:2325-2337. [PMID: 37458402 DOI: 10.1021/acs.biochem.3c00192] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Protein dynamics and function is strongly connected to the energy flow taking place. Myoglobin (Mb) and its mutations are ideal systems to study the process of vibrational energy transfer (VET) at the molecular level. Anti-Stokes ultraviolet resonance Raman studies using a tryptophan (Trp) probe, introduced at different Mb positions by amino acid replacement, have suggested that the amount of VET depends on the position of the Trp probe relative to the heme group. Inspired by this experimental work, we explored the strength of noncovalent π interactions, as well as covalent interactions for both the axial and distal ligands bound to iron in aquomet-Mb with the local vibrational mode analysis (LMA), originally developed by Konkoli and Cremer. Two sets of noncovalent interactions were investigated: (1) the interaction between the water ligand and Trp rings and (2) the interaction between the Trp and the porphyrin rings of the heme group. We assessed the strength of these noncovalent interactions via a special local mode force constant. Various Trp-modified water-bound ferric Mb proteins in the ground state were studied (6 in total) using gas-phase and QM/MM calculations followed by LMA. Our results disclose that VET is indeed dependent on the position of the Trp probe relative to the heme group but also on the tautomeric nature of distal histidine. They provide new guidelines on how to assess noncovalent π interactions in proteins utilizing LMA and how to use these data to explore VET, and more generally protein dynamics and function.
Collapse
Affiliation(s)
- Juliana J Antonio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas 75275-0314, United States
| |
Collapse
|
17
|
Peluzo BMTC, Makoś MZ, Moura RT, Freindorf M, Kraka E. Linear versus Bent Uranium(II) Metallocenes─A Local Vibrational Mode Study. Inorg Chem 2023. [PMID: 37478353 DOI: 10.1021/acs.inorgchem.3c01761] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2023]
Abstract
Uranium metallocenes have recently attracted attention driven by their use as catalysts in organometallic synthesis. In addition to bent U(IV) and U(III), an U(II) metallocene [(η5-C5i Pr5)2U] was synthesized with an unusual linear Cp-U-Cp angle. In this work, we investigated 22 U(II) metallocenes, (i) assessing the intrinsic strength of the U-ring interactions in these complexes with a novel bond strength measure based on our local vibrational mode analysis and (ii) systematically exploring what makes these U(II) metallocenes bent. We included relativistic effects through the NESCau Hamiltonian and complemented the local mode analysis with natural bonding orbital (NBO) and quantum theory of atoms in molecules (QTAIM) data. Our study led to the following results: (i) reduction of bulky U-ring ligand substituents does not lead to bent complexes for alkyl substituents (iPr and iBu) in contrast to SiMe3 ring substituents, which are all bent. (ii) The most bent complexes are [(η5-C5H4SiMe3)2U] (130°) and [η5-P5H5)2U] (143°). (iii) Linear complexes showed one hybridized NBO with s/d character, while bent structures were characterized by s/d/f mixing. (iv) We did not observe a correlation between the strength of the U-ring interaction and the amount of the ring-U-ring bend; the strongest interaction was found for [η5-Cp)2U] and the weakest for [η5-P5H5)2U]. In conclusion, our results provide a foundation for the design of U(II) metallocenes with specific physicochemical properties and increased reactivity.
Collapse
Affiliation(s)
- Bárbara M T C Peluzo
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Małgorzata Z Makoś
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37830, United States
| | - Renaldo T Moura
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraíba, Areia 58397-000, Paraíba, Brazil
| | - Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
18
|
Schiller C, Sieh D, Lindenmaier N, Stephan M, Junker N, Reijerse E, Granovsky AA, Burger P. Cleavage of an Aromatic C-C Bond in Ferrocene by Insertion of an Iridium Nitrido Nitrogen Atom. J Am Chem Soc 2023; 145:11392-11401. [PMID: 37172080 DOI: 10.1021/jacs.3c02781] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
The intermolecular cleavage of C-C bonds is a rare event. Herein, we report on a late transition-metal terminal nitrido complex, which upon oxidation undergoes insertion of the nitrido nitrogen atom into the aromatic C-C bond of ferrocene. This reaction path was confirmed through 15N and deuterium isotope labeling experiments of the nitrido complex and ferrocenium, respectively. Cyclic voltammetry and UV/vis spectroscopy monitoring of the reaction revealed that oxidation is the initial step, yielding the tentative radical cationic nitrido complex, which is experimentally supported by extended X and Q-band electron paramagnetic resonance (EPR) and ENDOR, UV/vis, vT 1H NMR, and vibrational spectroscopic data. Density functional theory (DFT) and multireference calculations of this highly reactive intermediate revealed an S = 1/2 ground state. The high reactivity can be traced to the increased electrophilicity in the oxidized complex. Based on high-level PNO-UCCSD(T) calculations and UV/vis kinetic measurements, it is proposed that the reaction proceeds by initial electrophilic exo attack of the nitrido nitrogen atom at the cyclopentadienyl ring and consecutive ring expansion to a pyridine ring.
Collapse
Affiliation(s)
- Carl Schiller
- Institut für Angewandte und Anorganische Chemie, Fachbereich Chemie, Universität Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Daniel Sieh
- Institut für Angewandte und Anorganische Chemie, Fachbereich Chemie, Universität Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Nils Lindenmaier
- Institut für Angewandte und Anorganische Chemie, Fachbereich Chemie, Universität Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Michel Stephan
- Institut für Angewandte und Anorganische Chemie, Fachbereich Chemie, Universität Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Natascha Junker
- Institut für Angewandte und Anorganische Chemie, Fachbereich Chemie, Universität Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Edward Reijerse
- Max-Planck-Institut für chemische Energiekonversion, EPR Research Group, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Alexander A Granovsky
- Institut für Angewandte und Anorganische Chemie, Fachbereich Chemie, Universität Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Peter Burger
- Institut für Angewandte und Anorganische Chemie, Fachbereich Chemie, Universität Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| |
Collapse
|
19
|
Li J, Wang C, Mo Y. Selectivity Rule of Cryptands for Anions: Molecular Rigidity and Bonding Site. Chemistry 2023; 29:e202203558. [PMID: 36538660 DOI: 10.1002/chem.202203558] [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/15/2022] [Revised: 12/15/2022] [Accepted: 12/20/2022] [Indexed: 12/24/2022]
Abstract
Cryptands utilize inside CH or NH groups as hydrogen bond (H-bond) donors to capture anions such as halides. In this work, the nature and selectivity of confined hydrogen bonds inside cryptands were computationally analyzed with the energy decomposition scheme based on the block-localized wavefunction method (BLW-ED), aiming at an elucidation of governing factors in the binding between cryptands and anions. It was revealed that the intrinsic strengths of inward hydrogen bonds are dominated by the electrostatic attraction, while the anion preferences (selectivity) of inner CH and NH hydrogen bonds are governed by the Pauli exchange repulsion and electrostatic interaction, respectively. Typical conformers of cages are classified into two groups, including the C3(h) -symmetrical conformers, in which all halide anions are located near the centroids of cages, and the "semi-open" conformers, which exhibit shifted bonding sites for different halide anions. Accordingly, the difference in governing factors of selectivity is attributed to either the rigidity of cages or the binding site of anions for these two groups. In details, the C3 conformers of NH cryptands can be enlarged more remarkably than the C3(h) -symmetrical conformers of CH cryptands as the size of anion (ionic radius) increases, resulting in the relaxation of the Pauli repulsion and a dramatic reduction in electrostatic attraction, which eventually rules the selectivity of NH cryptands for halide anions. By contrary, the CH cryptands are more rigid and cannot effectively reduce the Pauli repulsion, which subsequently governs the anion preference. Unlike C3 conformers whose rigidity determines the selectivity, semi-open conformers exhibit different binding sites for different anions. From F- to I- , the bonding site shifts toward the outside end of the pocket inside the semi-open NH cryptand, leading to the significant reduction of the electrostatic interaction that dominates the anion preference. Differently, binding sites are much less affected by the size of anion inside the semi-open CH cryptand, in which the Pauli exchange repulsion remains the key factor for the selectivity of inner hydrogen bonds.
Collapse
Affiliation(s)
- Jiayao Li
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi'an, 710119, China
| | - Changwei Wang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi'an, 710119, China
| | - Yirong Mo
- Department of Nanoscience, Joint School of Nanoscience & Nanoengineering, University of North Carolina at Greensboro, Greensboro, NC, 27401, USA
| |
Collapse
|
20
|
Ellington TL, Devore DP, Uvin G De Alwis WM, French KA, Shuford KL. Shedding Light on the Vibrational Signatures in Halogen-Bonded Graphitic Carbon Nitride Building Blocks. Chemphyschem 2022; 24:e202200812. [PMID: 36480235 DOI: 10.1002/cphc.202200812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Revised: 12/02/2022] [Accepted: 12/08/2022] [Indexed: 12/13/2022]
Abstract
The relative contributions of halogen and hydrogen bonding to the interaction between graphitic carbon nitride monomers and halogen bond (XB) donors containing C-X and C≡C bonds were evaluated using computational vibrational spectroscopy. Conventional probes into select vibrational stretching frequencies can often lead to disconnected results. To elucidate this behavior, local mode analyses were performed on the XB donors and complexes identified previously at the M06-2X/aVDZ-PP level of theory. Due to coupling between low and high energy C-X vibrations, the C≡C stretch is deemed a better candidate when analyzing XB complex properties or detecting XB formation. The local force constants support this conclusion, as the C≡C values correlate much better with the σ-hole magnitude than their C-X counterparts. The intermolecular local stretching force constants were also assessed, and it was found that attractive forces other than halogen bonding play a supporting role in complex formation.
Collapse
Affiliation(s)
- Thomas L Ellington
- Department of Chemistry and Biochemistry, Baylor University, One Bear Place #97348, Waco, TX, 76798-7348, USA
| | - Daniel P Devore
- Department of Chemistry and Biochemistry, Baylor University, One Bear Place #97348, Waco, TX, 76798-7348, USA
| | - W M Uvin G De Alwis
- Department of Chemistry and Biochemistry, Baylor University, One Bear Place #97348, Waco, TX, 76798-7348, USA
| | - Kirk A French
- Department of Chemistry and Biochemistry, Baylor University, One Bear Place #97348, Waco, TX, 76798-7348, USA
| | - Kevin L Shuford
- Department of Chemistry and Biochemistry, Baylor University, One Bear Place #97348, Waco, TX, 76798-7348, USA
| |
Collapse
|
21
|
Moura RT, Quintano M, Antonio JJ, Freindorf M, Kraka E. Automatic Generation of Local Vibrational Mode Parameters: From Small to Large Molecules and QM/MM Systems. J Phys Chem A 2022; 126:9313-9331. [DOI: 10.1021/acs.jpca.2c07871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Renaldo T. Moura
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas75275-0314, United States
- Department of Chemistry and Physics Center of Agrarian Sciences, Federal University of Paraiba, Areia, PB58397-000, Brazil
| | - Mateus Quintano
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas75275-0314, United States
| | - Juliana J. Antonio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas75275-0314, United States
| | - Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas75275-0314, United States
| |
Collapse
|
22
|
Kraka E, Quintano M, La Force HW, Antonio JJ, Freindorf M. The Local Vibrational Mode Theory and Its Place in the Vibrational Spectroscopy Arena. J Phys Chem A 2022; 126:8781-8798. [DOI: 10.1021/acs.jpca.2c05962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| | - Mateus Quintano
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| | - Hunter W. La Force
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| | - Juliana J. Antonio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| | - Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| |
Collapse
|
23
|
Hasil A, Beck D, Schröder D, Pillet S, Wenger E, Woike T, Klüfers P, Schaniel D. Pas de Deux of an NO Couple: Synchronous Photoswitching from a Double-Linear to a Double-Bent Ru(NO) 2 Core under Nitrosyl Charge Conservation. Angew Chem Int Ed Engl 2022; 61:e202210671. [PMID: 35983847 PMCID: PMC9826364 DOI: 10.1002/anie.202210671] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Indexed: 01/11/2023]
Abstract
The {Ru(NO)2 }10 dinitrosylruthenium complex [Ru(NO)2 (PPh3 )2 ] (1) shows photo-induced linkage isomerism (PLI) of a special kind: the two NO ligands switch, on photo-excitation, synchronously from the ground state (GS) with two almost linear RuNO functions to a metastable state (MS) which persists up to 230 K and can be populated to ≈50 %. The MS was experimentally characterised by photo-crystallography, IR spectroscopy and DS-calorimetry as a double-bent variant of the double-linear GS. The experimental results are confirmed by computation which unravels the GS/MS transition as a disrotatory synchronous 50° turn of the two nitrosyl ligands. Although 1 shows the usual redshift of the N-O stretch on bending the MNO unit, there is no increased charge transfer from Ru to NO along the GS-to-MS path. In terms of the effective-oxidation-state (EOS) method, both isomers of 1 and the transition state are Ru-II (NO+ )2 species.
Collapse
Affiliation(s)
- Asma Hasil
- Université de Lorraine, CNRS, CRM254000NancyFrance
| | - Daniel Beck
- Department Chemie der Ludwigs-Maximilians-UniversitätButenandtstraße 5–1381377MünchenGermany
| | - Daniel Schröder
- Department Chemie der Ludwigs-Maximilians-UniversitätButenandtstraße 5–1381377MünchenGermany
| | | | | | - Theo Woike
- Université de Lorraine, CNRS, CRM254000NancyFrance
| | - Peter Klüfers
- Department Chemie der Ludwigs-Maximilians-UniversitätButenandtstraße 5–1381377MünchenGermany
| | | |
Collapse
|
24
|
Madushanka A, Verma N, Freindorf M, Kraka E. Papaya Leaf Extracts as Potential Dengue Treatment: An In-Silico Study. Int J Mol Sci 2022; 23:ijms232012310. [PMID: 36293162 PMCID: PMC9610845 DOI: 10.3390/ijms232012310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 10/04/2022] [Accepted: 10/06/2022] [Indexed: 11/05/2022] Open
Abstract
Dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) cause serious public health problems, with nearly 390 million people affected and 20,000 deaths per year in tropical and subtropical countries. Despite numerous attempts, no antiviral drug or vaccine is currently available to combat the manifestation. The challenge of discovering an efficient vaccine is enhanced by the surplus presence of efficient vectors and drug resistance from the virus. For centuries, papaya (Carica papaya) extracts have been traditionally used to treat DF, DHF, and DSS. In the present study, we systematically investigated seven compounds isolated from papaya leaf extract with regard to their potential as inhibitors for non-structural (NS) proteins, NS3 and NS5, which play a crucial role in viral RNA replication. The computational tools applied stretched across classical molecular docking, molecular dynamics (MD) simulations and SwissADME used to calculate binding affinities; binding free energies; Absorption, Distribution, Metabolism, and Excretion (ADME); and drug-likeness properties, thus, identifying Kaempferol, Chlorogenic acid, and Quercetin as potential candidates, with Kaempferol and Quercetin scoring best. Therefore, for the Kaempferol and Quercetin complexes, hybrid quantum mechanical/molecular mechanical (QM/MM) geometry and frequency calculations were performed, followed by the local mode analysis developed in our group to quantify Kaempferol-NS and Quercetin-NS hydrogen bonding. Given the non-toxic nature and the wide availability of the Kaempferol and Quercetin papaya extract in almost all of the susceptible regions, and our results showing high NS3 and NS5 binding affinities and energies, strong hydrogen bonding with both NS3 and NS5, and excellent ADME properties, we suggest Kaempferol and Quercetin as a strong NS3 and NS5 inhibitor to be further investigated in vitro.
Collapse
|
25
|
Popp J, Kluefers P. A Diatomic Acceptor Ligand’s Bond Strength: a Reliable Measure of Its Antibond Occupation and Its Charge? Eur J Inorg Chem 2022. [DOI: 10.1002/ejic.202200374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Jens Popp
- LMU Munich Faculty for Chemistry and Pharmacy: Ludwig-Maximilians-Universitat Munchen Fakultat fur Chemie und Pharmazie Chemistry GERMANY
| | - Peter Kluefers
- Ludwig-Maximilians-Universitat Munchen Fakultat fur Chemie und Pharmazie Department of Chemistry Butenandtstrasse 5-13Haus D 81377 München GERMANY
| |
Collapse
|
26
|
Hasil A, Beck D, Schröder D, Pillet S, Wenger E, Woike T, Klüfers P, Schaniel D. Pas de Deux of an NO Couple: Synchronous Photoswitching from a Double‐Linear to a Double‐Bent Ru(NO)2 Core under Nitrosyl Charge Conservation. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202210671] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Asma Hasil
- Université de Lorraine, CNRS CRM2 FRANCE
| | - Daniel Beck
- Ludwig-Maximilians-Universität München: Ludwig-Maximilians-Universitat Munchen Chemie GERMANY
| | - Daniel Schröder
- Ludwig-Maximilians-Universität München: Ludwig-Maximilians-Universitat Munchen Chemie GERMANY
| | | | | | - Theo Woike
- Université de Lorraine, CNRS CRM2 FRANCE
| | - Peter Klüfers
- Ludwig-Maximilians-Universität München: Ludwig-Maximilians-Universitat Munchen Chemie GERMANY
| | - Dominik Schaniel
- Université de Lorraine, CNRS CRM2 Bvd des Aiguillettes 54000 Nancy FRANCE
| |
Collapse
|
27
|
Freindorf M, Delgado AAA, Kraka E. CO bonding in hexa‐ and pentacoordinate carboxy‐neuroglobin: A quantum mechanics/molecular mechanics and local vibrational mode study. J Comput Chem 2022. [DOI: 10.1002/jcc.26973] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Marek Freindorf
- Department of Chemistry Southern Methodist University Dallas Texas USA
| | | | - Elfi Kraka
- Department of Chemistry Southern Methodist University Dallas Texas USA
| |
Collapse
|
28
|
Mechanistic Details of the Sharpless Epoxidation of Allylic Alcohols—A Combined URVA and Local Mode Study. Catalysts 2022. [DOI: 10.3390/catal12070789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Abstract
In this work, we investigated the catalytic effects of a Sharpless dimeric titanium (IV)–tartrate–diester catalyst on the epoxidation of allylalcohol with methyl–hydroperoxide considering four different orientations of the reacting species coordinated at the titanium atom (reactions R1–R4) as well as a model for the non-catalyzed reaction (reaction R0). As major analysis tools, we applied the URVA (Unified Reaction Valley Approach) and LMA (Local Mode Analysis), both being based on vibrational spectroscopy and complemented by a QTAIM analysis of the electron density calculated at the DFT level of theory. The energetics of each reaction were recalculated at the DLPNO-CCSD(T) level of theory. The URVA curvature profiles identified the important chemical events of all five reactions as peroxide OO bond cleavage taking place before the TS (i.e., accounting for the energy barrier) and epoxide CO bond formation together with rehybridization of the carbon atoms of the targeted CC double bond after the TS. The energy decomposition into reaction phase contribution phases showed that the major effect of the catalyst is the weakening of the OO bond to be broken and replacement of OH bond breakage in the non-catalyzed reaction by an energetically more favorable TiO bond breakage. LMA performed at all stationary points rounded up the investigation (i) quantifying OO bond weakening of the oxidizing peroxide upon coordination at the metal atom, (ii) showing that a more synchronous formation of the new CO epoxide bonds correlates with smaller bond strength differences between these bonds, and (iii) elucidating the different roles of the three TiO bonds formed between catalyst and reactants and their interplay as orchestrated by the Sharpless catalyst. We hope that this article will inspire the computational community to use URVA complemented with LMA in the future as an efficient mechanistic tool for the optimization and fine-tuning of current Sharpless catalysts and for the design new of catalysts for epoxidation reactions.
Collapse
|
29
|
Zhao J, Chi CX, Meng LY, Jiang XL, Grunenberg J, HU HS, Zhou M, Li J, Schwarz W. Cis- and Trans-Binding Influences in [NUO · (N2)n]+ . J Chem Phys 2022; 157:054301. [DOI: 10.1063/5.0098068] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Uranium nitride-oxide cations [NUO]+ and their complexes with equatorial N2 ligands, [NUO·(N2) n]+ ( n=1-7), were synthesized in the gas phase. Mass-selected infrared photo-dissociation spectroscopy and quantum-chemical calculations confirm [NUO·(N2)5]+ as the sterically fully coordinated cation, with electronic singlet ground state of 1A1, linear [NUO]+ core, and C5v structure. The short N-U bond distances and high stretching modes, with slightly elongated U-O bond distances and lowered stretching modes, are rationalized as due to cooperative covalent and dative [ǀN≡U≡Oǀ]+ triple bonds. The mutual trans-interaction through the flexible electronic U-5f6d7sp valence shell, and the linearly increasing perturbation by an increasing number of equatorial dative N2 ligands are rationalized. It highlights the bonding and distinctiveness of uranium chemistry.
Collapse
Affiliation(s)
| | | | - Lu-Yan Meng
- East China University of Technology, Nanchang, China
| | - Xue-Lian Jiang
- Southern University of Science and Technology, Shenzhen, China
| | | | | | | | - Jun Li
- Tsinghua University, China
| | | |
Collapse
|
30
|
Margaret McCutcheon M, Freindorf M, Kraka E. Bonding in Nitrile Photo-dissociating Ruthenium Drug Candidates --A Local Vibrational Mode Study. J Chem Phys 2022; 157:014301. [DOI: 10.1063/5.0094567] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we investigated bonding features 15 ruthenium complexes of the type [Ru(tpy)(L)-(CH3CN)]n+, containing the tridentate tpy ligand (tpy = 2,2':6',2'--terpyridine) and various bidentate ancillary ligands, 12 compounds originally synthesized by Loftus et al. (J. Phys. Chem. C 123, 10291-10299 (2019)) complemented with three additional complexes. The main focus of our work was to relate these local features to the experimental data of Loftus et al. which assess the efficiency of nitrile release in an indirect way via observed quantum yields for ruthenium water association after nitrile release. As a tool to quantitatively assess Ru-NC and Ru-L bonding we utilized the local vibrational mode analysis complemented by the topological analysis of the electron density and the natural bond orbital analysis. Interestingly, the stronger Ru-NC bonds have the greater observed quantum yields, leading to the conclusion that the observed quantum yields are a result of a complex interplay of several processes excluding a direct relationship between QY and Ru-NC or Ru-L bond strengths. We identified the ST splitting as one of the key players and not the Ru-NC bond strength, as one may have thought. In summary, this work has presented a modern computational tool set for the investigation of bonding features applied to nitrile photo-dissociating ruthenium drug candidates forming a valuable basis for future design and fine tuning of nitrile releasing ruthenium compounds, as well as for the understanding of how local properties affect overall experimental outcomes.
Collapse
Affiliation(s)
| | | | - Elfi Kraka
- Chemistry, Southern Methodist University, United States of America
| |
Collapse
|
31
|
Martín Pendás Á, Francisco E. The role of references and the elusive nature of the chemical bond. Nat Commun 2022; 13:3327. [PMID: 35680893 PMCID: PMC9184482 DOI: 10.1038/s41467-022-31036-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Accepted: 05/30/2022] [Indexed: 11/16/2022] Open
Abstract
Chemical bonding theory is of utmost importance to chemistry, and a standard paradigm in which quantum mechanical interference drives the kinetic energy lowering of two approaching fragments has emerged. Here we report that both internal and external reference biases remain in this model, leaving plenty of unexplored territory. We show how the former biases affect the notion of wavefunction interference, which is purportedly recognized as the most basic bonding mechanism. The latter influence how bonding models are chosen. We demonstrate that the use of real space analyses are as reference-less as possible, advocating for their use. Delocalisation emerges as the reference-less equivalent to interference and the ultimate root of bonding. Atoms (or fragments) in molecules should be understood as a statistical mixture of components differing in electron number, spin, etc. The theory of chemical bonding relies on arbitrary references. Here the authors report a fundamental study on the chemical bond showing that considering the binding fragments as objects in real space enables to eliminate inherent biases.
Collapse
Affiliation(s)
- Ángel Martín Pendás
- Departamento de Química Física y Analítica, Universidad de Oviedo, 33006, Oviedo, Spain.
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Universidad de Oviedo, 33006, Oviedo, Spain
| |
Collapse
|
32
|
Costa Peluzo BMT, Kraka E. Uranium: The Nuclear Fuel Cycle and Beyond. Int J Mol Sci 2022; 23:ijms23094655. [PMID: 35563047 PMCID: PMC9101921 DOI: 10.3390/ijms23094655] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Revised: 04/13/2022] [Accepted: 04/18/2022] [Indexed: 01/01/2023] Open
Abstract
This review summarizes the recent developments regarding the use of uranium as nuclear fuel, including recycling and health aspects, elucidated from a chemical point of view, i.e., emphasizing the rich uranium coordination chemistry, which has also raised interest in using uranium compounds in synthesis and catalysis. A number of novel uranium coordination features are addressed, such the emerging number of U(II) complexes and uranium nitride complexes as a promising class of materials for more efficient and safer nuclear fuels. The current discussion about uranium triple bonds is addressed by quantum chemical investigations using local vibrational mode force constants as quantitative bond strength descriptors based on vibrational spectroscopy. The local mode analysis of selected uranium nitrides, N≡U≡N, U≡N, N≡U=NH and N≡U=O, could confirm and quantify, for the first time, that these molecules exhibit a UN triple bond as hypothesized in the literature. We hope that this review will inspire the community interested in uranium chemistry and will serve as an incubator for fruitful collaborations between theory and experimentation in exploring the wealth of uranium chemistry.
Collapse
|
33
|
BF3–Catalyzed Diels–Alder Reaction between Butadiene and Methyl Acrylate in Aqueous Solution—An URVA and Local Vibrational Mode Study. Catalysts 2022. [DOI: 10.3390/catal12040415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
In this study we investigate the Diels–Alder reaction between methyl acrylate and butadiene, which is catalyzed by BF3 Lewis acid in explicit water solution, using URVA and Local Mode Analysis as major tools complemented with NBO, electron density and ring puckering analyses. We considered four different starting orientations of methyl acrylate and butadiene, which led to 16 DA reactions in total. In order to isolate the catalytic effects of the BF3 catalyst and those of the water environment and exploring how these effects are synchronized, we systematically compared the non-catalyzed reaction in gas phase and aqueous solution with the catalyzed reaction in gas phase and aqueous solution. Gas phase studies were performed at the B3LYP/6-311+G(2d,p) level of theory and studies in aqueous solution were performed utilizing a QM/MM approach at the B3LYP/6-311+G(2d,p)/AMBER level of theory. The URVA results revealed reaction path curvature profiles with an overall similar pattern for all 16 reactions showing the same sequence of CC single bond formation for all of them. In contrast to the parent DA reaction with symmetric substrates causing a synchronous bond formation process, here, first the new CC single bond on the CH2 side of methyl acrylate is formed followed by the CC bond at the ester side. As for the parent DA reaction, both bond formation events occur after the TS, i.e., they do not contribute to the energy barrier. What determines the barrier is the preparation process for CC bond formation, including the approach diene and dienophile, CC bond length changes and, in particular, rehybridization of the carbon atoms involved in the formation of the cyclohexene ring. This process is modified by both the BF3 catalyst and the water environment, where both work in a hand-in-hand fashion leading to the lowest energy barrier of 9.06 kcal/mol found for the catalyzed reaction R1 in aqueous solution compared to the highest energy barrier of 20.68 kcal/mol found for the non-catalyzed reaction R1 in the gas phase. The major effect of the BF3 catalyst is the increased mutual polarization and the increased charge transfer between methyl acrylate and butadiene, facilitating the approach of diene and dienophile and the pyramidalization of the CC atoms involved in the ring formation, which leads to a lowering of the activation energy. The catalytic effect of water solution is threefold. The polar environment leads also to increased polarization and charge transfer between the reacting species, similar as in the case of the BF3 catalyst, although to a smaller extend. More important is the formation of hydrogen bonds with the reaction complex, which are stronger for the TS than for the reactant, thus stabilizing the TS which leads to a further reduction of the activation energy. As shown by the ring puckering analysis, the third effect of water is space confinement of the reacting partners, conserving the boat form of the six-member ring from the entrance to the exit reaction channel. In summary, URVA combined with LMA has led to a clearer picture on how both BF3 catalyst and aqueous environment in a synchronized effort lower the reaction barrier. These new insights will serve to further fine-tune the DA reaction of methyl acrylate and butadiene and DA reactions in general.
Collapse
|
34
|
Tao Y, Zou W, Nanayakkara S, Kraka E. LModeA-nano: A PyMOL Plugin for Calculating Bond Strength in Solids, Surfaces, and Molecules via Local Vibrational Mode Analysis. J Chem Theory Comput 2022; 18:1821-1837. [PMID: 35192350 DOI: 10.1021/acs.jctc.1c01269] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The analysis of chemical bonding in crystal structures and surfaces is an important research topic in theoretical chemistry. In this work, we present a PyMOL plugin, named LModeA-nano, as implementation of the local vibrational mode theory for periodic systems (Tao et al. J. Chem. Theory Comput. 2019, 15, 1761) assessing bond strength in terms of local stretching force constants in extended systems of one, two, and three dimensions. LModeA-nano can also analyze chemical bonds in isolated molecular systems thus enabling a head-to-head comparison of bond strength across systems with different dimensions in periodicity (0-3D). The new code is interfaced to the output generated by various solid-state modeling packages including VASP, CP2K, Quantum ESPRESSO, CASTEP, and CRYSTAL. LModeA-nano is cross-platform, open-source and freely available on GitHub: https://github.com/smutao/LModeA-nano.
Collapse
Affiliation(s)
- Yunwen Tao
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an, Shaanxi 710127, P. R. China
| | - Sadisha Nanayakkara
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
35
|
Liang X, Guan F, Ling Z, Wang H, Tao Y, Kraka E, Huang H, Yu C, Li D, He J, Fang H. Pivotal role of water molecules in the photodegradation of pymetrozine: New insights for developing green pesticides. JOURNAL OF HAZARDOUS MATERIALS 2022; 423:127197. [PMID: 34844343 DOI: 10.1016/j.jhazmat.2021.127197] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/13/2021] [Revised: 08/04/2021] [Accepted: 09/08/2021] [Indexed: 06/13/2023]
Abstract
Photodegradation of the insecticide pymetrozine (PYM) was studied on surface of wax films, and in aqueous and nonaqueous phase. The half-life of PYM on the wax surface was approximately 250 times longer than in water. Scavenging experiments, laser flash photolysis, and spectra analysis indicated the first singlet excited state of PYM (S1 *PYM) to be the most important photoinduced species initiating the photodegradation. Quantum chemistry calculations identified significant molecular torsion and changes in the structure C-CN-N of S1 *PYM, and the absolute charges of the CN atoms increased and the bond strength weakened. Free energy surface analysis, and O18 labeling experiments further confirmed that the mechanism was two-step photoinduced hydrolysis. The first step is the hydrolysis of S1 *PYM at CN upon reaction with 2-3 water molecules (one H2O molecule as the catalyst). The second step is an intramolecular hydrogen transfer coupled with the cleavage of C-N bond and formation of two cyclic products. During the interactions, water molecules experience catalytic activation by transferring protons, while there is a negligible solvent effect. Clarifying the detailed photodegradation mechanisms of PYM is beneficial for the development of green pesticides that are photostable and effective on leaf surfaces, and photolabile and detoxified in the aquatic environment.
Collapse
Affiliation(s)
- Ximei Liang
- College of Animal Science and Technology, Jiangxi Agricultural University, Nanchang 330045, China
| | - Fangling Guan
- Key Laboratory of Poyang Lake Basin Agricultural Resource and Ecology of Jiangxi Province, College of Land Resource and Environment, Jiangxi Agricultural University, Nanchang 330045, China
| | - Zhiyou Ling
- Key Laboratory of Poyang Lake Basin Agricultural Resource and Ecology of Jiangxi Province, College of Land Resource and Environment, Jiangxi Agricultural University, Nanchang 330045, China
| | - Honghong Wang
- State Key Joint Laboratory of Environment Simulation and Pollution Control, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
| | - Yunwen Tao
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314, USA
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314, USA
| | - Huajun Huang
- Key Laboratory of Poyang Lake Basin Agricultural Resource and Ecology of Jiangxi Province, College of Land Resource and Environment, Jiangxi Agricultural University, Nanchang 330045, China
| | - Chenglong Yu
- Key Laboratory of Poyang Lake Basin Agricultural Resource and Ecology of Jiangxi Province, College of Land Resource and Environment, Jiangxi Agricultural University, Nanchang 330045, China
| | - Danping Li
- Key Laboratory of Poyang Lake Basin Agricultural Resource and Ecology of Jiangxi Province, College of Land Resource and Environment, Jiangxi Agricultural University, Nanchang 330045, China
| | - Jinbao He
- Key Laboratory of Poyang Lake Basin Agricultural Resource and Ecology of Jiangxi Province, College of Land Resource and Environment, Jiangxi Agricultural University, Nanchang 330045, China
| | - Hansun Fang
- Key Laboratory of Poyang Lake Basin Agricultural Resource and Ecology of Jiangxi Province, College of Land Resource and Environment, Jiangxi Agricultural University, Nanchang 330045, China.
| |
Collapse
|
36
|
Santos-Jr CV, A. F. de Souza M, Kraka E, Moura Jr RT. Analysis of spectator chemical bonds in SN2@C and @Si reaction mechanisms in the gas phase. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2021.139282] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
37
|
Nanayakkara S, Tao Y, Kraka E. Capturing Individual Hydrogen Bond Strengths in Ices via Periodic Local Vibrational Mode Theory: Beyond the Lattice Energy Picture. J Chem Theory Comput 2021; 18:562-579. [PMID: 34928619 DOI: 10.1021/acs.jctc.1c00357] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Local stretching force constants derived from periodic local vibrational modes at the vdW-DF2 density functional level have been employed to quantify the intrinsic hydrogen bond strength of 16 ice polymorphs, ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI, XII, XIII, XIV, XV, XVII, and XIX, that are stable under ambient to elevated pressures. Based on this characterization on 1820 hydrogen bonds, relationships between local stretching force constants and structural parameters such as hydrogen bond length and angle were identified. Moreover, different bond strength distributions, from uniform to inhomogeneous, were observed for the 16 ices and could be explained in relation to different local structural elements within ices, that is, rings, that consist of different hydrogen bond types. In addition, criteria for the classification of hydrogen bonds as strong, intermediate, and weak were introduced. The latter was used to explore a different dimension of the water-ice phase diagram. These findings will provide important guidelines for assessing the credibility of new ice structures.
Collapse
Affiliation(s)
- Sadisha Nanayakkara
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Yunwen Tao
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
38
|
Kabir MP, Orozco-Gonzalez Y, Hastings G, Gozem S. The effect of hydrogen-bonding on flavin's infrared absorption spectrum. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 262:120110. [PMID: 34224983 DOI: 10.1016/j.saa.2021.120110] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 06/17/2021] [Accepted: 06/21/2021] [Indexed: 06/13/2023]
Abstract
Cluster and continuum solvation computational models are employed to model the effect of hydrogen bonding interactions on the vibrational modes of lumiflavin. Calculated spectra were compared to experimental Fourier-transform infrared (FTIR) spectra in the diagnostic 1450-1800 cm-1 range, where intense νC=C, νC=N, [Formula: see text] , and [Formula: see text] stretching modes of flavin's isoalloxazine ring are found. Local mode analysis is used to describe the strength of hydrogen-bonding in cluster models. The computations indicate that νC=C and νC=N mode frequencies are relatively insensitive to intermolecular interactions while the [Formula: see text] and [Formula: see text] modes are sensitive to direct (and also indirect for [Formula: see text] ) hydrogen-bonding interactions. Although flavin is neutral, basis sets without the diffuse functions provide incorrect relative frequencies and intensities. The 6-31+G* basis set is found to be adequate for this system, and there is limited benefit to considering larger basis sets. Calculated vibrational mode frequencies agree with experimentally determined frequencies in solution when cluster models with multiple water molecules are used. Accurate simulation of relative FTIR band intensities, on the other hand, requires a continuum (or possibly quantum mechanical/molecular mechanical) model that accounts for long-range electrostatic effects. Finally, an experimental peak at ca. 1624 cm-1 that is typically assigned to the [Formula: see text] vibrational stretching mode has a complicated shape that suggests multiple underlying contributions. Our calculations show that this band has contributions from both the C6-C7 and C2 = O stretching vibrations.
Collapse
Affiliation(s)
- Mohammad Pabel Kabir
- Department of Chemistry, Georgia State University, Atlanta, GA 30302, United States
| | | | - Gary Hastings
- Department of Physics and Astronomy, Georgia State University, Atlanta, GA 30302, United States; Center for Nano-Optics, Georgia State University, Atlanta, GA 30302, United States.
| | - Samer Gozem
- Department of Chemistry, Georgia State University, Atlanta, GA 30302, United States.
| |
Collapse
|
39
|
Taherivardanjani S, Elfgen R, Reckien W, Suarez E, Perlt E, Kirchner B. Benchmarking the Computational Costs and Quality of Vibrational Spectra from Ab Initio Simulations. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202100293] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Shima Taherivardanjani
- Mulliken Center for Theoretical Chemistry Institute for Physical and Theoretical Chemistry Beringstr. 4 Bonn D‐53115 Germany
| | - Roman Elfgen
- Mulliken Center for Theoretical Chemistry Institute for Physical and Theoretical Chemistry Beringstr. 4 Bonn D‐53115 Germany
| | - Werner Reckien
- Mulliken Center for Theoretical Chemistry Institute for Physical and Theoretical Chemistry Beringstr. 4 Bonn D‐53115 Germany
| | - Estela Suarez
- Institute for Advanced Simulation Jülich Supercomputing Centre, Forschungszentrum Jülich GmbH Wilhelm‐Johnen‐Straße Jülich D‐52425 Germany
| | - Eva Perlt
- Otto Schott Institute of Materials Research Faculty of Physics and Astronomy Friedrich‐Schiller‐Universität Jena Löbdergraben 32 Jena D‐07743 Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry Institute for Physical and Theoretical Chemistry Beringstr. 4 Bonn D‐53115 Germany
| |
Collapse
|
40
|
Delgado AAA, Sethio D, Matthews D, Oliveira V, Kraka E. Substituted hydrocarbon: a CCSD(T) and local vibrational mode investigation. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1970844] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
| | - Daniel Sethio
- Department of Chemistry – BMC, Uppsala University, Uppsala, Sweden
| | - Devin Matthews
- Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Vytor Oliveira
- Departamento de Química, Instituto Tecnológico de Aeronáutica (ITA), São José dos Campos, São Paulo, Brazil
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| |
Collapse
|
41
|
Oliveira VP, Marcial BL, Machado FBC, Kraka E. Relating Bond Strength and Nature to the Thermodynamic Stability of Hypervalent Togni-Type Iodine Compounds. Chempluschem 2021; 86:1199-1210. [PMID: 34437775 DOI: 10.1002/cplu.202100285] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2021] [Revised: 08/09/2021] [Indexed: 11/08/2022]
Abstract
The bond strength and nature of a set of 32 Togni-like reagents have been investigated at the M062X/def2-TZVP(D) level of theory in acetonitrile described with the SMD continuum solvent model, to rationalize the main factors responsible for their thermodynamic stability in different conformations, and trifluoromethylation capabilities. For the assessment of bond strength, we utilized local stretching force constants and associated bond strength orders, complemented with local features of the electron density to access the nature of the bonds. Bond dissociation energies varied from 31.6 to 79.9 kcal/mol depending on the polarizing power of the ligand trans to CF3 . Based on the analysis of the Laplacian of the density, we propose that the charge-shift bond character plays an important role in the stability of the molecules studied, especially for those containing I-O bonds. New insights on the trans influence and on possible ways to fine-tune the stability of these reagents are provided.
Collapse
Affiliation(s)
- Vytor Pinheiro Oliveira
- Instituto Tecnológico de Aeronáutica (ITA), Departamento de Química, São José dos Campos, 12228-900, São Paulo, Brazil
| | - Bruna Luana Marcial
- Instituto Federal Goiano (IF Goiano), Núcleo de Química, Campus Morrinhos, Goiás, Brazil
| | - Francisco B C Machado
- Instituto Tecnológico de Aeronáutica (ITA), Departamento de Química, São José dos Campos, 12228-900, São Paulo, Brazil
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas, 75275-0314, USA
| |
Collapse
|
42
|
Christopher IL, Michalchuk AAL, Pulham CR, Morrison CA. Towards Computational Screening for New Energetic Molecules: Calculation of Heat of Formation and Determination of Bond Strengths by Local Mode Analysis. Front Chem 2021; 9:726357. [PMID: 34354982 PMCID: PMC8329490 DOI: 10.3389/fchem.2021.726357] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2021] [Accepted: 07/02/2021] [Indexed: 11/13/2022] Open
Abstract
The reliable determination of gas-phase and solid-state heats of formation are important considerations in energetic materials research. Herein, the ability of PM7 to calculate the gas-phase heats of formation for CNHO-only and inorganic compounds has been critically evaluated, and for the former, comparisons drawn with isodesmic equations and atom equivalence methods. Routes to obtain solid-state heats of formation for a range of single-component molecular solids, salts, and co-crystals were also evaluated. Finally, local vibrational mode analysis has been used to calculate bond length/force constant curves for seven different chemical bonds occurring in CHNO-containing molecules, which allow for rapid identification of the weakest bond, opening up great potential to rationalise decomposition pathways. Both metrics are important tools in rationalising the design of new energetic materials through computational screening processes.
Collapse
Affiliation(s)
- Imogen L Christopher
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| | | | - Colin R Pulham
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| | - Carole A Morrison
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
43
|
Altun Z, Bleda EA, Trindle C. Focal Point Evaluation of Energies for Tautomers and Isomers for 3-Hydroxy-2-Butenamide: Evaluation of Competing Internal Hydrogen Bonds of Types -OH…O=, -OH…N, -NH…O=, and CH…X (X=O and N). Molecules 2021; 26:molecules26092623. [PMID: 33946257 PMCID: PMC8124485 DOI: 10.3390/molecules26092623] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Revised: 04/22/2021] [Accepted: 04/27/2021] [Indexed: 11/16/2022] Open
Abstract
The title compound is a small molecule with many structural variations; it can illustrate a variety of internal hydrogen bonds, among other noncovalent interactions. Here we examine structures displaying hydrogen bonding between carbonyl oxygen and hydroxyl H; between carbonyl oxygen and amino H; hydroxyl H and amino N; hydroxyl O and amino H. We also consider H-bonding in its tautomer 2-oxopropanamide. By extrapolation algorithms applied to Hartree-Fock and correlation energies as estimated in HF, MP2, and CCSD calculations using the cc-pVNZ correlation-consistent basis sets (N = 2, 3, and 4) we obtain reliable relative energies of the isomeric forms. Assuming that such energy differences may be attributed to the presence of the various types of hydrogen bonding, we attempt to infer relative strengths of types of H-bonding. The Atoms in Molecules theory of Bader and the Local Vibrational Modes analysis of Cremer and Kraka are applied to this task. Hydrogen bonds are ranked by relative strength as measured by local stretching force constants, with the stronger =O…HO- > NH…O= > -OH…N well separated from a cluster > NH…O= ≈ >NH…OH ≈ CH…O= of comparable and intermediate strength. Weaker but still significant interactions are of type CH…N which is stronger than CH…OH.
Collapse
Affiliation(s)
- Zikri Altun
- Physics Department, Marmara University, Göztepe Kampus, Istanbul 34772, Turkey; (Z.A.); (E.A.B.)
| | - Erdi Ata Bleda
- Physics Department, Marmara University, Göztepe Kampus, Istanbul 34772, Turkey; (Z.A.); (E.A.B.)
| | - Carl Trindle
- Chemistry Department, University of Virginia, Charlottesville, VA 22902, USA
- Correspondence: ; Tel.: +1-4347709197
| |
Collapse
|
44
|
Assessing the Intrinsic Strengths of Ion–Solvent and Solvent–Solvent Interactions for Hydrated Mg2+ Clusters. INORGANICS 2021. [DOI: 10.3390/inorganics9050031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of their intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters.
Collapse
|
45
|
Makoś MZ, Freindorf M, Tao Y, Kraka E. Theoretical Insights into [NHC]Au(I) Catalyzed Hydroalkoxylation of Allenes: A Unified Reaction Valley Approach Study. J Org Chem 2021; 86:5714-5726. [PMID: 33780251 DOI: 10.1021/acs.joc.1c00208] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Hydroxylation is an effective approach for the synthesis of carbon-oxygen bonds and allylic ethers. The [NHC]Au(I) catalyzed intermolecular hydroalkoxylation of allene was studied at the DFT and Coupled Cluster level of theory. Using the Unified Reaction Valley Approach (URVA), we carry out a comprehensive mechanistic analysis of [NHC]Au(I)-catalyzed and noncatalyzed reactions. The URVA study of several possible reaction pathways reveal that the [NHC]Au(I) catalyst enables the hydroalkoxylation reaction to occur via a two step mechanism based upon the Au ability to switch between π- and σ-complexation. The first step of the mechanism involves the formation of a CO bond after the transition state with no energy penalty. Following the CO bond breakage, the OH bond breaks and CH bond forms during the second step of the mechanism, as the catalyst transforms into the more stable π-Au complex. The URVA results were complemented with local vibrational mode analysis to provide measures of intrinsic bond strength for Au(I)-allene interactions of all stationary points, and NBO analysis was applied in order to observe charge transfer events along the reaction pathway. Overall, the π-Au C═C interactions of the products are stronger than those of the reactants adding to their exothermicity. Our work on the hydroxylation of allene provides new insights for the design of effective reaction pathways to produce allylic ethers and also unravels new strategies to form C-O bonds by activation of C═C bonds.
Collapse
Affiliation(s)
- Małgorzata Z Makoś
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Yunwen Tao
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
46
|
Beiranvand N, Freindorf M, Kraka E. Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study. Molecules 2021; 26:2268. [PMID: 33919989 PMCID: PMC8071019 DOI: 10.3390/molecules26082268] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Revised: 04/08/2021] [Accepted: 04/09/2021] [Indexed: 12/13/2022] Open
Abstract
In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.
Collapse
Affiliation(s)
| | | | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314, USA; (N.B.); (M.F.)
| |
Collapse
|
47
|
Abstract
We systematically investigated iodine–metal and iodine–iodine bonding in van Koten’s pincer complex and 19 modifications changing substituents and/or the transition metal with a PBE0–D3(BJ)/aug–cc–pVTZ/PP(M,I) model chemistry. As a novel tool for the quantitative assessment of the iodine–metal and iodine–iodine bond strength in these complexes we used the local mode analysis, originally introduced by Konkoli and Cremer, complemented with NBO and Bader’s QTAIM analyses. Our study reveals the major electronic effects in the catalytic activity of the M–I–I non-classical three-center bond of the pincer complex, which is involved in the oxidative addition of molecular iodine I2 to the metal center. According to our investigations the charge transfer from the metal to the σ* antibonding orbital of the I–I bond changes the 3c–4e character of the M–I–I three-center bond, which leads to weakening of the iodine I–I bond and strengthening of the metal–iodine M–I bond, facilitating in this way the oxidative addition of I2 to the metal. The charge transfer can be systematically modified by substitution at different places of the pincer complex and by different transition metals, changing the strength of both the M–I and the I2 bonds. We also modeled for the original pincer complex how solvents with different polarity influence the 3c–4e character of the M–I–I bond. Our results provide new guidelines for the design of pincer complexes with specific iodine–metal bond strengths and introduce the local vibrational mode analysis as an efficient tool to assess the bond strength in complexes.
Collapse
|
48
|
Tao Y, Zou W, Nanayakkara S, Freindorf M, Kraka E. A revised formulation of the generalized subsystem vibrational analysis (GSVA). Theor Chem Acc 2021; 140:31. [PMID: 33716564 PMCID: PMC7942689 DOI: 10.1007/s00214-021-02727-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2020] [Accepted: 02/08/2021] [Indexed: 11/30/2022]
Abstract
In this work, a simplified formulation of our recently developed generalized subsystem vibrational analysis (GSVA) for obtaining intrinsic fragmental vibrations (J Chem Theory Comput 14:2558, 2018) is presented. In contrast to the earlier implementation, which requires the explicit definition of a non-redundant set of internal coordinate parameters to be constructed for the subsystem, the new implementation circumvents this process by employing massless Eckart conditions to the subsystem fragment paired with a Gram-Schmidt orthogonalization to span the same internal vibration space indirectly. This revised version of GSVA (rev-GSVA) can be applied to equilibrium structure as well as transition state structure, and it has been incorporated into the open-source package UniMoVib (https://github.com/zorkzou/UniMoVib). Supplementary Information The online version contains supplementary material available at 10.1007/s00214-021-02727-y.
Collapse
Affiliation(s)
- Yunwen Tao
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314 USA
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an, 710127 Shaanxi People's Republic of China
| | - Sadisha Nanayakkara
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314 USA
| | - Marek Freindorf
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314 USA
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314 USA
| |
Collapse
|
49
|
Verma N, Tao Y, Kraka E. Systematic Detection and Characterization of Hydrogen Bonding in Proteins via Local Vibrational Modes. J Phys Chem B 2021; 125:2551-2565. [DOI: 10.1021/acs.jpcb.0c11392] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Niraj Verma
- Department of Chemistry, Southern Methodist University, Dallas Texas United States
| | - Yunwen Tao
- Department of Chemistry, Southern Methodist University, Dallas Texas United States
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, Dallas Texas United States
| |
Collapse
|
50
|
Complexes of Formaldehyde and α-Dicarbonyls with Hydroxylamine: FTIR Matrix Isolation and Theoretical Study. Molecules 2021; 26:molecules26041144. [PMID: 33672783 PMCID: PMC7924657 DOI: 10.3390/molecules26041144] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2021] [Revised: 02/16/2021] [Accepted: 02/18/2021] [Indexed: 11/29/2022] Open
Abstract
The interactions of formaldehyde (FA), glyoxal (Gly) and methylglyoxal (MGly) with hydroxylamine (HA) isolated in solid argon and nitrogen were studied using FTIR spectroscopy and ab initio methods. The spectra analysis indicates the formation of two types of hydrogen-bonded complexes between carbonyl and hydroxylamine in the studied matrices. The cyclic planar complexes are stabilized by O–H⋯O(C), and C–H⋯N interactions and the nonplanar complexes are stabilized by O–H⋯O(C) bond. Formaldehyde was found to form with hydroxylamine, the cyclic planar complex and methylglyoxal, the nonplanar one in both argon and nitrogen matrices. In turn, glyoxal forms with hydroxylamine the most stable nonplanar complex in solid argon, whereas in solid nitrogen, both types of the complex are formed.
Collapse
|