We present a comprehensive investigation on the different role of CO in carboxyneuroglobin (1) as ligand of the heme group in the active site forming a bond with the heme iron and (2) dissociated from the heme group but still trapped inside the active site, focusing on two specific orientations, one with CO perpendicular to the plane defined by the distal histidine of the enzyme (form A) and one with CO located parallel to that plane (form B). Our study includes wild type carboxy-neuroglobin and nine known protein mutations. Considering that the distal histidine interacting with the heme group can adapt two different tautomeric forms and the two possible orientations of the dissociated CO, a total of 36 protein systems were analyzed in this study. Fully optimized geometries and vibrational frequencies were calculated at the QM/MM level, followed by the local mode analysis, to decode CO bond properties.The intrinsic bond strengths derived from the local mode analysis, complemented with NBO and QTAIM data, reveal that the strength of the CO bond, in the hexacoordinate (where CO is a ligand of the heme group) and pentacoordinate (where CO is dissociated from the heme group) scenarios, is dominated by through bond and through space charge transfer between CO and Fe, fine-tuned by electrostatic and dispersion interactions with the side chain amino acids in the distal heme pocket.Suggestions are made as to advise on how protein modifications can influence the molecular properties of the coordinated or dissociated CO, which could serve the fine-tuning of existing and the design of new neuroglobin models with specific FeC and CO bond strengths.

For decades one has strived to synthesize a compound with the longest covalent C-C bond applying predominantly steric hindrance and/or strain to achieve this goal. On the other hand electronic effects have been added to the repertoire, such as realized in the electron deficient ethane radical cation in its D3d form. Recently, negative hyperconjugation effects occurring in diamino-o-carborane analogs such as di-N,N-dimethylamino-o-carborane have been held responsible for their long C-C bonds. In this work we systematically analyzed CC bonding in a diverse set of 53 molecules including clamped bonds, highly sterically strained complexes such as diamondoid dimers, electron deficient species, and di-N,N-dimethylamino-o-carborane to cover the whole spectrum of possibilities for elongating a covalent C-C bond to the limit. As a quantitative intrinsic bond strength measure, we utilized local vibrational CC stretching force constants ka(CC) and related bond strength orders BSO n(CC), computed at the ωB97X-D/aug-cc-pVTZ level of theory. Our systematic study quantifies for the first time that whereas steric hindrance and/or strain definitely elongate a C-C bond, electronic effects can lead to even longer and weaker C-C bonds. Within our set of molecules the electron deficient ethane radical cation, in D3d symmetry, acquires the longest C-C bond with a length of 1.935 Å followed by di-N,N-dimethylamino-o-carborane with a bond length of 1.930 Å. However, the C-C bond in di-N,N-dimethylamino-o-carborane is the weakest with a BSO n value of 0.209 compared to 0.286 for the ethane radical cation; another example that the longer bond is not always the weaker bond. Based on our findings we provide new guidelines for the general characterization of CC bonds based on local vibrational CC stretching force constants and for future design of compounds with long C-C bonds.

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 Ca2+ 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(H2O)n2+ (n = 1–8) clusters correlates with the electron density. The application of a second hydration shell to Ca(H2O)n2+ (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 Ca2+, 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.

Two-dimensional infrared spectroscopy has reported highly delocalized in-plane base vibrations in the fingerprint region of nucleotide monophosphates, suggesting the involvement of base pair C=O and C=C ring bonds and considerable interaction between C=O bond stretches. The high delocalization results in congested vibrational spectra, which complicates the assignment of the peaks. This congestion also extends to Watson-Crick base pairs. We applied in this work the characterization of normal mode procedure, a special feature of our local mode analysis, and could for the first time identify the C=O and C=C bonds being engaged in base pare coupling and quantify their contribution to each of the delocalized fingerprint vibration. In addition, a detailed and quantitative description of the hydrogen bonds involved in the Watson-Crick base pairs was provided. Based on the results of this study, we developed a new protocol to elucidate on the assignment of bands in the vibrational spectra of nucleic acids by probing the vibrational space for specific interactions between functional groups prior to and upon base pairing. This protocol will aid to fill the gap between DNA structural information and vibrational spectroscopy experiments by facilitating the interpretation of spectra on a quantitative basis.

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.