Computational analysis of vibrational frequencies and rovibrational spectroscopic constants of hydrogen sulfide dimer using MP2 and CCSD(T)

Docking and Electronic Structure of Rutin, Myricetin, and Baicalein Targeting 3CLpro

Understanding the role of 3CLpro protease for SARS-CoV-2 replication and knowing the potential of flavonoid molecules like rutin, myricetin, and baicalein against 3CLpro justify an investigation into their inhibition. This study investigates possible bonds and reactivity descriptors of rutin, myricetin, and baicalein through conformational and electronic properties. Density functional theory was used to determine possible interactions. Analyses were carried out through the molecular electrostatic potential, electron localization function, Fukui function descriptors based on frontier orbitals, and non-covalent interactions. A docking study was performed using a resolution of 1.55 Å for 3CLpro to analyze the interactions of rutin, myricetin, and baicalein. Scores of structures showed that rutin is the best ligand, followed by myricetin and baicalein. Docking studies showed that baicalein and rutin can establish effective interactions with residues of the catalytic dyad (Cys145 and His41), but just rutin forms a hydrogen bond. Myricetin, in turn, could not establish an effective interaction with Cys145. Baicalein interaction arose with active residues such as Arg188, Val186, Gln189, and Gln192. Interactions of rutin and myricetin with Arg188 and Gln189 were also found. A critical interaction was observed only for rutin with the hydroxyls of ring A with His41, and also for Cys145 with rings B and C, which is probably related to the highest score of rutin.

Ab Initio Approach to the Structure, Vibrational Properties, and Electron Binding Energies of H2S∙∙∙SO2

The present study employs high-level ab initio calculations to investigate the structure, vibrational frequencies, and electronic properties of H2S∙∙∙SO2. The analysis of vibrational frequencies reveals an intramolecular vibrational energy transfer phenomenon, where energy from the stretching modes of H2S is transferred to the ν1s mode of SO2. At the CCSD(T)/aug-cc-pVQZ level, the interaction energy between H2S and SO2 is predicted to be 2.78 kcal/mol. Electron propagator theory calculations yield a HOMO-LUMO gap of 8.24 eV for H2S∙∙∙SO2. Furthermore, by utilizing ab initio results for the adiabatic ionization energy and electron affinity, the electrophilicity of H2S∙∙∙SO2 is estimated to be 2.01 eV. This value is similar to the electrophilicity of SO2, suggesting comparable reactivity and chemical behavior. The non-covalent interaction (NCI) analysis of the H2S∙∙∙SO2 complex emphasizes the significant contribution of non-covalent van der Waals interactions in its energetic stabilization.

Formations of bimolecular barbituric acid complexes through hydrogen bonding interactions: DFT analyses of structural and electronic features

Formations of bimolecular barbituric acid (BA) complexes through hydrogen-bonding (HB) interactions were investigated in this work. BA has been known as a starting compound of pharmaceutical compounds developments, in which the molecular and atomic features of parent BA in homo-paring with another BA molecule were investigated here. The models were optimized to reach the stabilized structures and their properties were evaluated at the molecular and atomic scales. Density functional theory (DFT) calculations were performed to provide required information for achieving the goal of this work. Six dimer models were obtained finally according to examining all possible starting dimers configurations for involving in optimization calculations. N-H . . . O and C-H . . . O interactions were also involved in dimers formations besides participation of the X-center of parent BA in interaction. Molecular and atomic scales features were evaluated for characterizing the dimers formations. As a consequence, several configurations of BA dimers were obtained showing the importance of performing such structural analyses for developing further compounds from BA.

Assessing the Intrinsic Strengths of Ion–Solvent and Solvent–Solvent Interactions for Hydrated Mg2+ Clusters

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.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

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.

Halogen Bonding Involving I2 and d8 Transition-Metal Pincer Complexes

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.

Local vibrational mode analysis of ion-solvent and solvent-solvent interactions for hydrated Ca2+ clusters

Hydrated calcium ion clusters have received considerable attention due to their essential role in biological processes such as bone development, hormone regulation, blood coagulation, and neuronal signaling. To better understand the biological role of the cation, the interactions between the Ca²⁺ ions and water molecules have been frequently investigated. However, a quantitative measure for the intrinsic Ca-O (ion-solvent) and intermolecular hydrogen bond (solvent-solvent) interactions has been missing so far. Here, we report a topological electron density analysis and a natural population analysis to analyze the nature of these interactions for a set of 14 hydrated calcium clusters via local mode stretching force constants obtained at the ωB97X-D/6-311++G(d,p) level of theory. The results revealed that the strength of inner Ca-O interactions for Ca(H₂O)_n ²⁺ (n = 1-8) clusters correlates with the electron density. The application of a second hydration shell to Ca(H₂O)_n ²⁺ (n = 6-8) clusters resulted in stronger Ca-O interactions where a larger electron charge transfer between lp(O) of the first hydration shell and the lower valence of Ca prevailed. The strength of the intermolecular hydrogen bonds, formed between the first and second hydration shells, became stronger when the charge transfers between hydrogen bond (HB) donors and HB acceptors were enhanced. From the local mode stretching force constants of implicitly and explicitly solvated Ca²⁺, we found the six-coordinated cluster to possess the strongest stabilizations, and these results prove that the intrinsic bond strength measures for Ca-O and hydrogen bond interactions form new effective tools to predict the coordination number for the hydrated calcium ion clusters.

The validation of molecular interaction among dimer chitosan with urea and creatinine using density functional theory: In application for hemodyalisis membrane

The formation of chitosan dimer and its interaction with urea and creatinine have been investigated at the density functional theory (DFT) level (B3LYP-D3/6-31++G**) to study the transport phenomena in hemodialysis membrane. The interaction energy of chitosan-creatinine and chitosan-urea complexes are in range -4 kcal/mol < interaction energy <-20 kcal/mol which were classified in medium hydrogen bond interaction. The chemical reactivity parameter proved that creatinine was more electrophilic and easier to bind chitosan than urea. The energy gap of HOMO-LUMO of chitosan-creatinine complex was lower than chitosan-urea complex that indicating chitosan-creatinine complex was more reactive and easier to transport electron than chitosan-urea complex. Moreover, the natural bond orbital (NBO) analysis showed a high contribution of hydrogen bond between chitosan-creatinine and chitosan-urea. The chitosan-creatinine interaction has a stronger hydrogen bond than chitosan-urea through the interaction O18-H34....N56 with stabilizing energy = -13 kcal/mol. The quantum theory atom in molecule (QTAIM) also supported NBO data. All data presented that creatinine can make hydrogen bond interaction stronger with chitosan than urea, that indicated creatinine easier to transport in the chitosan membrane than urea during hemodialysis process.

Modeling Hydrogen Release from Water with Borane and Alane Catalysts: A Unified Reaction Valley Approach

The unified reaction valley approach combined with the local vibrational mode and ring puckering analysis is applied to investigate the hydrogen evolution from water in the presence of small hydrides such as BH₃, metal hydrides as AlH₃, and their derivatives. We studied a series of reactions involving BH₃, AlH₃, B₂H₆, Al₂H₆, and AlH₃BH₃ with one- and two-water molecules, considering multiple reaction paths. In addition, the influence of the aqueous medium was examined. A general reaction mechanism was identified for most of the reactions. Those that deviate could be associated with unusually high reaction barriers with no hydrogen release. The charge transfer along the reaction path suggests that a viable hydrogen release is achieved when the catalyst adopts the role of a charge donor during the chemical processes. The puckering analysis showed that twistboat and boat forms are the predominant configurations in the case of an intermediate six-membered ring formation, which influences the activation barrier. The local mode analysis was used as a tool to detect the H-H bond formation as well as to probe catalyst regenerability. Based on the correlation between the activation energy and the change in the charge separation for cleaving O-H and B(Al)-H bonds, two promising subsets of reactions could be identified along with prescriptions for lowering the reaction barrier individually with electron-donating/withdrawing substituents.