Noncovalent Carbon-Bonding Interactions in Proteins

Cooperativity between Intermolecular Hydrogen and Carbon Bonds in ZY···CH3CN/CH3NC···HX Trimers (ZY = H2O, H2S, HF, HCl, HBr, NH3, and H2CO; HX = HF, HCl, and HBr)

Hydrogen-bonding and carbon-bonding interactions are widespread in nature. We studied the cooperativity between these interactions in 42 trimeric complexes ZY···CH3CN/CH3NC···HX, where ZY molecules are H2O, H2S, HF, HCl, HBr, NH3, and H2CO, and HX molecules are HF, HCl, and HBr. Acetonitrile (CH3CN) and isoacetonitrile (CH3NC) act as hydrogen bond acceptors as well as carbon bond donors in these trimers. Various theoretical methods, such as electronic structure calculations, quantum theory of atoms in molecule (QTAIM), natural bond orbital (NBO), and reduced density gradient analysis, are employed to study these trimers, and the results are compared with the corresponding ZY···CH3CN/CH3NC and CH3CN/CH3NC···HX dimers. Electronic structure calculations are performed at the second-order Mo̷ller-Plesset perturbation theory using the 6-311++G(2d,2p) basis set. We show that both the interactions act synergistically in these trimers leading to an increase in their bond strength as compared to the strength in the individual dimers. The cooperative energies for these trimers are in the range of 0.69 to 3.22 kJ/mol. It is seen that the carbon bonds benefit more from the cooperativity than the hydrogen bonds. The trends of cooperativity and correlations of interaction energies and cooperative energies with relevant QTAIM and NBO parameters are reported.

Metal-Induced Enhancement of Tetrel Bonding. The Case of C⋅⋅⋅X-IrIII (X=Cl, Br) Tetrel Bond Involving a Methyl Group

In X-ray structures of the isomorphic mer-[IrX3(THT)(CNXyl)2] (X=Cl 1, Br 2; THT=tetrahydrothiophene; Xyl=2,6-Me2C6H3-) complexes, we revealed short intermolecular contacts between the C-atom of an isocyanide methyl group and halide ligands of another molecule. Geometrical consideration of the X-ray data and analysis of appropriate DFT studies allowed the attribution of these contacts to CMe⋅⋅⋅X-IrIII (X=Cl, Br) tetrel bond. Specifically, through the application of DFT calculations and various theoretical models, the presence of tetrel bonding interactions was validated, and the contribution of the CMe⋅⋅⋅X-IrIII interaction was assessed. The reinforcement of the tetrel bond upon the isocyanide coordination to iridium(III) is substantiated by molecular electrostatic potential (MEP) surface calculations. To distinguish the tetrel bonding characteristics of CMe⋅⋅⋅X-IrIII (X=Cl, Br) interactions from conventional hydrogen bonding, we employed multiple computational methodologies, including Natural Bond Orbital (NBO) analysis and Electron Localization Function (ELF) analysis. Additionally, Energy Decomposition Analysis (EDA) was applied to selected model systems to explore the underlying physical nature of these interactions.

Triel Bonds between BH3/C5H4BX and M(MDA)2 (X = H, CN, F, CH3, NH2; M = Ni, Pd, Pt, MDA = Enolated Malondialdehyde) and Group 10 Transition Metal Electron Donors

A systematic theoretical study was conducted on the triel bonds (TrB) within the BH3∙∙∙M(MDA)2 and C5H4BX∙∙∙M(MDA)2 (M = Ni, Pd, Pt, X = H, CN, F, CH3, NH2, MDA = enolated malondialdehyde) complexes, with BH3 and C5H4BX acting as the electron acceptors and the square-coordinated M(MDA)2 acting as the electron donor. The interaction energies of these systems range between -4.71 and -33.18 kcal/mol. The larger the transition metal center M, the greater the enhancement of the TrB, with σ-hole TrBs found to be stronger than π-hole TrBs. In the σ-hole TrB complex, an electron-withdrawing substituent on the C opposite to the B atom enhances the TrB, while an electron-donating substituent has little effect on the strength of TrB in the Pd and Pt complexes but enhances the TrB in the Ni-containing complexes. The van der Waals interaction plays an important role in stabilizing these binary systems, and its contribution diminishes with increasing M size. The orbital effect within these systems is largely due to charge transfer from the dz2 orbital of M into the empty pz orbital of B.

Bioinspired Supramolecular Hydrogel from Design to Applications

Nature offers a wealth of opportunities to solve scientific and technological issues based on its unique structures and function. The dynamic non-covalent interaction is considered to be the main base of living functions of creatures including humans, animals, and plants. Supramolecular hydrogels formed by non-covalent bonding interactions has become a unique platform for constructing promising materials for medicine, energy, electronic, and biological substitute. In this review, the self-assemble principle of supramolecular hydrogels is summarized. Next, the stimulation of external environment that triggers the assembly or disassembly of supramolecular hydrogels are recapitulated, including temperature, mechanics, light, pH, ions, etc. The main applications of bioinspired supramolecular hydrogels in terms of bionic objects including humans, animals, and plants are also described. Although so many efforts are done for revealing the synergized mechanism of the function and non-covalent interactions on the supramolecular hydrogel, the complexity and variability between stimulus and non-covalent bonding in the supramolecular system still require impeccable theories. As an outlook, the bioinspired supramolecular hydrogel is just beginning to exhibit its great potential in human life, offering significant opportunities in drug delivery and screening, implantable devices and substitutions, tissue engineering, micro-fluidic devices, and biosensors.

Spodium bonding with noble gas atoms

The nature of the bonding between a neutral group 12 member (Zn3, Cd3 and Hg3) ring and a noble gas atom was explored using quantum chemical simulations. Natural bond orbital, quantum theory of atoms in molecules, symmetry-adapted perturbation theory, and molecular electrostatic potential surface analysis were also used to investigate the type of interaction between the noble gas atom and the metal rings (Zn3, Cd3 and Hg3). The Zn3, Cd3 and Hg3 rings are bonded to the noble gas through non-covalent interactions, which was revealed by the non-covalent interaction index. Additionally, energy decomposition analysis reveals that dispersion energy is the key factor in stabilizing these systems.

Search for Osme Bonds with π Systems as Electron Donors

The Osme bond is defined as pairing a Group 8 metal atom as an electron acceptor in a noncovalent interaction with a nucleophile. DFT calculations with the ωB97XD functional consider MO4 (M = Ru, Os) as the Lewis acid, paired with a series of π electron donors C2H2, C2H4, C6H6, C4H5N, C4H4O, and C4H4S. The calculations establish interaction energies in the range between 9.5 and 26.4 kJ/mol. Os engages in stronger interactions than does Ru, and those involving more extensive π-systems within the aromatic rings form stronger bonds than do the smaller ethylene and acetylene. Extensive analysis questions the existence of a true Osme bond, as the bonding chiefly involves interactions with the three O atoms of MO4 that lie closest to the π-system, via π(C-C)→σ*(M-O) transfers. These interactions are supplemented by back donation from M-O bonds to the π*(CC) antibonding orbitals of the π-systems. Dispersion makes a large contribution to these interactions, higher than electrostatics and much greater than induction.

Cooperativity between H-bonds and tetrel bonds. Transformation of a noncovalent C⋯N tetrel bond to a covalent bond

The dimers and trimers formed by imidazole (IM) and F2TO (T = C, Si, Ge) are studied by ab initio calculations. IM can engage in either a NH⋯O H-bond with F2TO or a T⋯N tetrel bond (TB) with the π-hole above the T atom. The latter is a true noncovalent TB for T = C but is a much shorter and stronger covalent bond with F2SiO or F2GeO. When a second IM is added, the cooperativity emerging from its H-bond with the first IM makes it a stronger nucleophile, leading to two minima with F2CO. The first structure contains a long noncovalent C⋯O TB and there is a much shorter covalent bond in the other, with a small energy barrier separating them. The same sort of double minimum occurs when the two IM units are situated parallel to one another in a stacked geometry.

The Tetrel Bonds of Hypervalent Halogen Compounds

The tetrel bond between PhXF2Y(TF3) (T = C and Si; X = Cl, Br, and I; Y = F and Cl) and the electron donor MCN (M = Li and Na) was investigated at the M06-2X/aug-cc-pVDZ level of theory. As the electronegativity of the halogen atom X increases, the strength of the tetrel bond also increases, but as the electronegativity of the halogen atom Y increases, the strength of the tetrel bond decreases. The magnitude of the interaction energy in most -CF3 complexes was found to be less than 10 kcal/mol, but to exceed 11 kcal/mol for PhClF2Cl(CF3)⋯NCNa. The tetrel bond is greatly enhanced when the -SiF3 group interacts with LiCN or NaCN, with the largest interaction energy approaching 100 kcal/mol and displaying a covalent Si⋯N interaction. Along with this enhancement, the Si⋯N distance was found to be less than the X-Si bond length, the -SiF3 group to be closer to the N atom, and in most -SiF3 systems, the X-Si-F angle to be less than 90°; the -SiF3 group therefore undergoes inversion and complete transfer in some systems.

Topology of electrostatic potential and electron density reveals a covalent to non-covalent carbon-carbon bond continuum

The covalent and non-covalent nature of carbon-carbon (CC) interactions in a wide range of molecular systems can be characterized using various methods, including the analysis of molecular electrostatic potential (MESP), represented as V(r), and the molecular electron density (MED), represented as ρ(r). These techniques provide valuable insights into the bonding between carbon atoms in different molecular environments. By uncovering a fundamental exponential relationship between the distance of the CC bond and the highest eigenvalue (λv1) of V(r) at the bond critical point (BCP), this study establishes the continuum model for all types of CC interactions, including transition states. The continuum model is further delineated into three distinct regions, namely covalent, borderline cases, and non-covalent, based on the gradient, , with the bond distance of the CC interaction. For covalent interactions, this parameter exhibits a more negative value than -5.0 a.u. Å-1, while for non-covalent interactions, it is less negative than -1.0 a.u. Å-1. Borderline cases, which encompass transition state structures, fall within the range of -1.0 to -5.0 a.u. Å-1. Furthermore, this study expands upon Popelier's analysis of the Laplacian of the MED, denoted as ∇2ρ, to encompass the entire spectrum of covalent, non-covalent, and borderline cases of CC interactions. Therefore, the present study presents compelling evidence supporting the concept of a continuum model for CC bonds in chemistry. Additionally, this continuum model is further explored within the context of C-N, C-O, C-S, N-N, O-O, and S-S interactions, albeit with a limited dataset.

C∙∙∙O and Si∙∙∙O Tetrel Bonds: Substituent Effects and Transfer of the SiF3 Group

The tetrel bond (TB) between 1,2-benzisothiazol-3-one-2-TF3-1,1-dioxide (T = C, Si) and the O atom of pyridine-1-oxide (PO) and its derivatives (PO-X, X = H, NO2, CN, F, CH3, OH, OCH3, NH2, and Li) is examined by quantum chemical means. The Si∙∙∙O TB is quite strong, with interaction energies approaching a maximum of nearly 70 kcal/mol, while the C∙∙∙O TB is an order of magnitude weaker, with interaction energies between 2.0 and 2.6 kcal/mol. An electron-withdrawing substituent on the Lewis base weakens this TB, while an electron-donating group has the opposite effect. The SiF3 group transfers roughly halfway between the N of the acid and the O of the base without the aid of cooperative effects from a third entity.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Computational Study of Driving Forces in ATSP, PDIQ, and P53 Peptide Binding: C═O···C═O Tetrel Bonding Interactions at Work

Understanding the molecular interactions that drive peptide folding is crucial to chemistry and biology. In this study, we analyzed the role of CO···CO tetrel bonding (TtB) interactions in the folding mechanism of three different peptides (ATSP, pDIQ, and p53), which exhibit a different propensity to fold in an α helix motif. To achieve this goal, we used both a recently developed Bayesian inference approach (MELDxMD) and Quantum Mechanics (QM) calculations at the RI-MP2/def2-TZVP level of theory. These techniques allowed us to study the folding process and to evaluate the strength of the CO···CO TtBs as well as the synergies between TtBs and hydrogen-bonding (HB) interactions. We believe that the results derived from our study will be helpful for those scientists working in computational biology, peptide chemistry, and structural biology.

Tetrel bonds involving a CF3 group participate in protein-drug recognition: a combined crystallographic and computational study

In this study, the ability of CF₃ groups to bind to the electron-rich side chains and backbone groups of proteins has been investigated by combining a Protein Data Bank (PDB) survey and ab initio quantum mechanics calculations. More precisely, an inspection of the PDB involving organic ligands containing a CF₃ group and electron-rich atoms (A = N, O and S) in the vicinity revealed 419 X-ray structures exhibiting CF₃⋯A tetrel bonds (TtBs). In a posterior stage, those hits that exhibited the most relevant features in terms of directionality and intermolecular distance were selected for theoretical calculations at the RI-MP2/def2-TZVPD level of theory. Also, Hammett's regression plots of several TtB complexes involving meta- and para-substituted benzene derivatives were computed to shed light on the substituent effects. Moreover, the TtBs were characterized through several state-of-the-art computational techniques, such as the Quantum Theory of Atoms in Molecules (QTAIM) and Noncovalent Interactions plot (NCIplot) methodologies. We believe that the results gathered from our study will be useful for rational drug design and biological communities as well as for further expanding the role of this interaction to biomedical applications.

Substituent Effects in Tetrel Bonds Involving Aromatic Silane Derivatives: An ab initio Study

In this manuscript substituent effects in several silicon tetrel bonding (TtB) complexes were investigated at the RI-MP2/def2-TZVP level of theory. Particularly, we have analysed how the interaction energy is influenced by the electronic nature of the substituent in both donor and acceptor moieties. To achieve that, several tetrafluorophenyl silane derivatives have been substituted at the meta and para positions by several electron donating and electron withdrawing groups (EDG and EWG, respectively), such as -NH₂, -OCH₃, -CH₃, -H, -CF₃ and -CN substituents. As electron donor molecules, we have used a series of hydrogen cyanide derivatives using the same EDGs and EWGs. We have obtained the Hammett's plots for different combinations of donors and acceptors and in all cases we have obtained good regression plots (interaction energies vs. Hammet's σ parameter). In addition, we have used the electrostatic potential (ESP) surface analysis as well as the Bader's theory of atoms in molecules (AIM) and noncovalent interaction plot (NCI plot) techniques to further characterize the TtBs studied herein. Finally, a Cambridge Structural Database (CSD) inspection was carried out, retrieving several structures where halogenated aromatic silanes participate in tetrel bonding interactions, being an additional stabilization force of their supramolecular architectures.

Geminal Charge-Assisted Tetrel Bonds in Bis-Pyridinium Methylene Salts

C(sp³) atoms are known to act as electrophilic sites in self-assembly processes, and in all cases reported till now, they form only one interaction with nucleophiles; that is, they function as monodentate tetrel bond donors. This manuscript reports experimental (X-ray structural analysis) and theoretical evidence (DFT calculations), proving that the methylene carbon in bis-pyridinium methylene salts establishes two short and directional C(sp³)···anion interactions; that is, they function as bidentate tetrel bond donors.

Beryllium bonding with noble gas atoms

Quantum chemical calculations were carried out to investigate the nature of the bonding between a neutral Be₃ ring and noble gas atom. Electronic structure calculation for these complexes was carried out at different computational levels in association with natural bond orbital, quantum theory of atoms in molecules, electron localization function, symmetry adapted perturbation theory, and molecular electrostatic potential surface analysis of Be₃ complexes. The Be atoms in the Be₃ moiety are chemically bonded to one another, with the BeBe bond dissociation energy being ~125 kJ mol^-1 . The Be₃ ring interacts with the noble gases through non-covalent interactions. The binding energies of the noble gas atoms with the Be₃ ring increases with increase in their atomic number. The non-covalent interaction index, density overlap region indicator and independent gradient model analyses reveal the presence of non-covalent inter-fragment interactions in the complexes. Energy decomposition analysis reveals that dispersion plays the major role towards stabilizing these systems.

Non-covalent interactions reveal the protein chain δ conformation in a flexible single-residue model

The δ conformation is a local secondary structure in proteins that implicates a π_amide N-H⋯N interaction between a backbone N atom and the NH of the following residue. Small-molecule models thereof have been limited so far to rigid proline-type compounds. We show here that in derivatives of a cyclic amino acid with a sulphur atom in the γ-position, specific side-chain/backbone N-H⋯S interactions stabilize the δ conformation sufficiently to allow it to compete with classical C5 and C7 H-bonded conformers.

Intermolecular noncovalent interactions with carbon in solution

One of the most familiar carbon-centered noncovalent interactions (NCIs) involving an antibonding π*-orbital situated at the Bürgi-Dunitz angle from the electron donor, mostly lone pairs of electrons, is known as n → π* interactions, and if it involves a σ* orbital in a linear fashion, then it is known as the carbon bond. These NCIs can be intra- or inter-molecular and are usually weak in strength but have a paramount effect on the structure and function of small-molecular crystals and proteins. Surprisingly, the experimental evidence of such interactions in the solution phase is scarce. It is even difficult to determine the interaction energy in the solution. Using NMR spectroscopy aided with molecular dynamics (MD) simulation and high-level quantum mechanical calculations, herein we provide the experimental evidence of intermolecular carbon-centered NCIs in solution. The challenge was to find appropriate heterodimers that could sustain room temperature thermal energy and collisions from the solvent molecules. However, after several trial model compounds, the pyridine-N-oxide:dimethyltetracyanocyclopropane (PNO-DMTCCP) complex was found to be a good candidate for the investigation. NBO analyses show that the PNO:DMTCCP complex is stabilized mainly by intermolecular n → π* interaction when a weaker carbon bond gives extra stability to the complex. From the NMR study, it is observed that the NCIs between DMTCCP and PNO are enthalpy driven with an enthalpy change of -28.12 kJ mol-1 and dimerization energy of ∼-38 kJ mol-1 is comparable to the binding energies of a conventional hydrogen-bonded dimer. This study opens up a new strategy to investigate weak intermolecular interactions such as n → π* interaction and carbon bonds in the solution phase.

Anion···Anion [AuI4]-···[AuI2]- Complex Trapped in the Solid State by Tetramethylammonium Cations

A discrete π-hole···σ-hole dimer is synthesized and X-ray characterized. It presents a perfect thumbtack geometry where the σ-hole of the linear [AuI₂]^- anion points to the π-hole located above the central Au-atom of the [AuI₄]^- anion. Such discrete π-hole···σ-hole dimers are unprecedented in literature, since all mixed-valence gold(I/III) iodide compounds reported to date form infinite ···([AuI₄]^-···[AuI₂]^-) _n ·· chains in the solid state. If an excess of iodine is used for the synthesis, triiodide [I₃]^- ions are partially incorporated into the [AuI₂]^- sites, forming infinite chains. The nature of the anion···anion interaction has been studied considering two possibilities: (i) a π-hole coinage bond or (ii) σ-hole halogen bond using high-level density functional theory calculations, the quantum theory of atoms in molecules, and the noncovalent interaction plot index.

Unique Dispersion-Induced Tetrel Bond with Co-operative σ-hole-Induced Pnicogen Bond in the POCl3-Acetone Heterodimer: Experimental Confirmation at Low Temperatures

Both tetrel and pnicogen bonds are known to be induced through σ-/π-holes. This work reports computational and experimental evidence of the carbonyl carbon of acetone hosting a tetrel bond by dispersion rather electrostatic forces, for the first time, while phosphorus of POCl₃ sustains pnicogen bonding via the σ-hole. Heterodimers of POCl₃ with acetone (CH₃COCH₃) have been isolated within inert gas matrixes of Ar and N₂ at 12 K. Characteristic vibrational bands at P═O stretching of POCl₃ and C═O stretching of CH₃COCH₃ have been obtained in support of the computations. The potential energy surface has been traced computationally using ab initio and density functional methods. CH₃COCH₃ harboring such a tetrel bond, in itself, is quite intriguing. The interplay of these interactions has been comprehended by the quantum theory of atoms in molecules, natural bond orbital, energy decomposition, electrostatic potential mapping, and noncovalent interaction analyses.

The Relevance of Experimental Charge Density Analysis in Unraveling Noncovalent Interactions in Molecular Crystals

The work carried out by our research group over the last couple of decades in the context of quantitative crystal engineering involves the analysis of intermolecular interactions such as carbon (tetrel) bonding, pnicogen bonding, chalcogen bonding, and halogen bonding using experimental charge density methodology is reviewed. The focus is to extract electron density distribution in the intermolecular space and to obtain guidelines to evaluate the strength and directionality of such interactions towards the design of molecular crystals with desired properties. Following the early studies on halogen bonding interactions, several "sigma-hole" interaction types with similar electrostatic origins have been explored in recent times for their strength, origin, and structural consequences. These include interactions such as carbon (tetrel) bonding, pnicogen bonding, chalcogen bonding, and halogen bonding. Experimental X-ray charge density analysis has proved to be a powerful tool in unraveling the strength and electronic origin of such interactions, providing insights beyond the theoretical estimates from gas-phase molecular dimer calculations. In this mini-review, we outline some selected contributions from the X-ray charge density studies to the field of non-covalent interactions (NCIs) involving elements of the groups 14-17 of the periodic table. Quantitative insights into the nature of these interactions obtained from the experimental electron density distribution and subsequent topological analysis by the quantum theory of atoms in molecules (QTAIM) have been discussed. A few notable examples of weak interactions have been presented in terms of their experimental charge density features. These examples reveal not only the strength and beauty of X-ray charge density multipole modeling as an advanced structural chemistry tool but also its utility in providing experimental benchmarks for the theoretical studies of weak interactions in crystals.

Noncovalent interactions in proteins and nucleic acids: beyond hydrogen bonding and π-stacking

Understanding the noncovalent interactions (NCIs) among the residues of proteins and nucleic acids, and between drugs and proteins/nucleic acids, etc., has extraordinary relevance in biomolecular structure and function. It helps in interpreting the dynamics of complex biological systems and enzymatic activity, which is esential for new drug design and efficient drug delivery. NCIs like hydrogen bonding (H-bonding) and π-stacking have been researchers' delight for a long time. Prominent among the recently discovered NCIs are halogen, chalcogen, pnictogen, tetrel, carbo-hydrogen, and spodium bonding, and n → π* interaction. These NCIs have caught the imaginations of various research groups in recent years while explaining several chemical and biological processes. At this stage, a holistic view of these new ideas and findings lying scattered can undoubtedly trigger our minds to explore more. The present review attempts to address NCIs beyond H-bonding and π-stacking, which are mainly n → σ*, n → π* and σ → σ* type interactions. Five of the seven NCIs mentioned earlier are linked to five non-inert end groups of the modern periodic table. Halogen (group-17) bonding is one of the oldest and most explored NCIs, which finds its relevance in biomolecules due to the phase correction and inhibitory properties of halogens. Chalcogen (group 16) bonding serves as a redox-active functional group of different active sites of enzymes and acts as a nucleophile in proteases and phosphates. Pnictogen (group 15), tetrel (group 14), triel (group 13) and spodium (group 12) bonding does exist in biomolecules. The n → π* interactions are linked to backbone carbonyl groups and protein side chains. Thus, they are crucial in determining the conformational stability of the secondary structures in proteins. In addition, a more recently discovered to and fro σ → σ* type interaction, namely carbo-hydrogen bonding, is also present in protein-ligand systems. This review summarizes these grand epiphanies routinely used to elucidate the structure and dynamics of biomolecules, their enzymatic activities, and their application in drug discovery. It also briefs about the future perspectives and challenges posed to the spectroscopists and theoreticians.

Carbon-Centered Hydrogen Bonds in Proteins

Hydrogen bonding (H-bonding) without lone pair(s) of electrons and π-electrons is a concept developed 2-3 years ago. H-bonds involving less electronegative tetrahedral carbon are beyond the classical concept of H-bonds. Herein, we present the first report on H-bonds with tetravalent carbons in proteins. A special bonding arrangement is needed to increase the negative charge density around the sp³-hybridized carbon atom. Therefore, less electronegative elements such as As and Mg, when bonded to sp³-C, enable the C-atoms as H-bond acceptors. Careful protein structure analysis aided by several quantum chemical calculations suggests that these H-bonds are weak to moderate in strength. We developed an empirical equation to estimate the C-H···C H-bond energy in proteins from the distances between the C- and H-atoms. In proteins, the binding energies range from -5.4 to -14.0 kJ/mol. The C-H···C H-bonds assist the substrate binding in proteins. We also explored the potential role of these carbon-centered H-bonds in C-H bond activation through σ-bond metathesis. To our surprise, contribution from these H-bonds is almost of similar magnitude as that from C-H···π H-bonds for C-H bond activation.

Structural topologies involving energetically significant antiparallel π-stacking and unconventional N(nitrile)⋯π(fumarate) contacts in dinuclear Zn(ii) and polymeric Mn(ii) compounds: antiproliferative evaluation and theoretical studies

Anticancer activities considering cell viability, apoptosis and molecular docking have been explored in dinuclear Zn(ii) and polymeric Mn(ii) compounds involving energetically significant unconventional N(nitrile)⋯π(fum) contacts.

Enhancement of tetrel bond involving tetrazole-TtR3 (Tt = C, Si; R = H, F). Promotion of SiR3 transfer by a triel bond

The combination of a CR3 (R = H, F) with a tetrazole can result in a moderate carbon bond, which can be further strengthened by a triel bond. On the other hand, SiR3 group is half transferred between the two N atoms in these conditions.

DeepNCI: DFT Noncovalent Interaction Correction with Transferable Multimodal Three-Dimensional Convolutional Neural Networks

A multimodal deep learning model, DeepNCI, is proposed for improving noncovalent interactions (NCIs) calculated via density functional theory (DFT). DeepNCI is composed of a three-dimensional convolutional neural network (3D CNN) for abstracting critical and comprehensive features from 3D electron density, and a neural network for modeling one-dimensional quantum chemical properties. By merging features from two networks, DeepNCI is able to reduce the root-mean-square error of DFT-calculated NCI from 1.19 kcal/mol to ∼0.2 kcal/mol for a NCI molecular database (>1000 molecules). The representativeness of the joint features can be visualized by t-distributed stochastic neighbor embedding (t-SNE), where they can distinguish categorized NCI systems quite well. Therefore, the fused model performs better than its component networks. In addition, the 3D CNN takes electron density as inputs that are in the same range, despite the size of molecular systems, so it can promote model applicability and transferability. To clarify the applicability of DeepNCI, an application domain (AD) has been defined with merged features using the K-nearest-neighbor method. The calculations for external test sets are shown that AD can properly monitor the reliability for a prediction. The model transferability is tested with a small database of homolysis bond dissociation energy including only dozens of samples. With NCI database pretrained parameters, the same or better performance than the reported results is achieved by transfer learning. This suggests that the DeepNCI model is transferable and it may transfer to other relative tasks, which possibly can resolve some small sampling problems. The source code of DeepNCI can be freely accessed at https://github.com/wenzelee/DeepNCI.

Theoretical and Experimental Insights into the Possible Interfacial Interactions between β-Glucan and Fat Molecules in Aqueous Media

Excessive body fat and high cholesterol are one of the leading reasons for triggering cardiovascular risk factors, obesity, and type 2 diabetes. Beta-glucan (BG)-based dietary fibers are found to be effective for lowering fat digestion in the gastrointestinal tract. However, the fat capturing mechanism of BG in aqueous medium is still elusive. In this report, we studied the dietary effect of barley-extracted BG on docosahexaenoic acid (DHA, a model fat molecule) uptake and the impact of the aqueous medium on their interactions using computational modeling and experimental parameters. The possible microscale and macroscale molecular interactions between BG and DHA in an aqueous medium were analyzed through density functional theory (DFT), Monte-Carlo (MC), and molecular dynamics (MD) simulations. DFT analysis revealed that the BG polymer extends hydrogen bonding and nonbonding interactions with DHA. Bulk simulation with multiple DHA molecules on a long-chain BG showed that a viscous colloidal system is formed upon increasing DHA loading. Experimental size and zeta potential measurements also confirmed the electrostatic interaction between BG-DHA systems. Furthermore, simulated and experimental diffusion and viscosity measurements showed excellent agreement. These simulated and experimental results revealed the mechanistic pathway of how BG fibers form colloidal systems with fat molecules, which is probably responsible for BG-induced delayed fat digestion and further halting of fatty molecule absorption in the GI tract.

Enhancement of the Tetrel Bond by the Effects of Substituents, Cooperativity, and Electric Field: Transition from Noncovalent to Covalent Bond

The T⋅⋅⋅N tetrel bond (TB) formed between TX₃ OH (T=C, Si, Ge; X=H, F) and the Lewis base N≡CM (M=H, Li, Na) is studied by ab initio calculations at the MP2/aug-cc-pVTZ level. Complexes involving TH₃ OH contain a conventional TB with interaction energy less than 10 kcal/mol. This bond is substantially strengthened, approaching 35 kcal/mol and covalent character, when fluorosubstituted TF₃ OH is combined with NCLi or NCNa. Along with this enhanced binding comes a near equalization of the TB T⋅⋅⋅N and the internal T-O bond lengths, and the associated structure acquires a trigonal bipyramidal shape, despite a high internal deformation energy. This structural transformation becomes more complete, and the TB is further strengthened upon adding an electron acceptor BeCl₂ to the Lewis acid and a base to the NCM unit. This same TB strengthening can be accomplished also by imposition of an external electric field.

N-H⋯X interactions stabilize intra-residue C5 hydrogen bonded conformations in heterocyclic α-amino acid derivatives

Nature makes extensive and elaborate use of hydrogen bonding to assemble and stabilize biomolecular structures. The shapes of peptides and proteins rely significantly on N–H⋯O Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> C interactions, which are the linchpins of turns, sheets and helices. The C5 H-bond, in which a single residue provides both donor and acceptor, is generally considered too weak to force the backbone to adopt extended structures. Exploiting the synergy between gas phase (experimental and quantum chemistry) and solution spectroscopies to decipher IR spectroscopic data, this work demonstrates that the extended C5-based conformation in 4-membered ring heterocyclic α-amino acid derivatives is significantly stabilized by the formation of an N–H⋯X H-bond. In this synergic system the strength of the C5 interaction remains constant while the N–H⋯X H-bond strength, and thereby the support provided by it, varies with the heteroatom.

In 4-membered ring heterocyclic α-amino acid derivatives, extended conformations based on intraresidue C5 H-bonds can be stabilized by N–H⋯X H-bonds, making the combined C5–C6γ structures prominent in both gas phase and in weakly polar solutions.

Se⋅⋅⋅O/S and S⋅⋅⋅O Chalcogen Bonds in Small Molecules and Proteins: A Combined CSD and PDB Study

The importance of selenium-centered noncovalent chalcogen bonds represented as Se⋅⋅⋅A (A=O/S) has been explored for short directional contacts in small molecules and proteins. In addition, S⋅⋅⋅O centered contacts have been analyzed. Computational analyses involving the quantitative assessment of the associated energetics, the molecular electrostatic potentials (MEP), and electron density derived topological parameters, namely, quantum theory of atom in molecules (QTAIM) analyses, and NBO (natural bond orbital) based calculations, have been performed to unequivocally establish the strength, stability, and attractive role of chalcogen bonds in the solid-state. This investigation has been performed in molecules from both the Cambridge Structural Database (CSD) and Protein Data Bank (PDB). Thus futuristic materials may be designed keeping in mind the significance of these interactions, including their relevance in biology.

Computational Evaluation of Me2 TCCP as Lewis Acid

Supramolecular adducts between dimethyl-2,2,3,3-tetracyanocyclopropane (Me2 TCCP) with 21 small (polar) molecules and 10 anions were computed with DFT (B3LYP-D3/def2-TZVP). Their optimized geometries were used to obtain interaction energies, and perform energy decomposition and 'atoms-in-molecules' analyses. A set of 38 other adducts were also evaluated for comparison purposes. Selected examples were further scrutinized by inspection of the molecular electrostatic potential maps, Noncovalent Interaction index plots, the Laplacian, the orbital interactions, and by estimating the Gibbs free energy of complexation in hexane solution. These calculations divulge the thermodynamic feasibility of Me2 TCCP adducts and show that complexation is typically driven by dispersion with less polarized partners, but by orbital interactions when more polarized or anionic guests are deployed. Most Me2 TCCP adducts are more stable than simple hydrogen bonding with water, but less stable than traditional Lewis adducts involving Me3 B, or a strong halogen bond such as with Br2 . Several bonding analyses showed that the locus of interaction is found near the electron poor sp3 -hydridized (NC)2 C-C(CN)2 carbon atoms. An empty hybrid σ*/π* orbital on Me2 TCCP was identified that can be held responsible for the stability of the most stable adducts due to donor-acceptor interactions.

Molecular Electrostatic Potential and Noncovalent Interactions in Derivatives of Group 8 Elements

This communication reports experimental and theoretical evidences of σ‐hole interactions in adducts between nitrogen or oxygen nucleophiles and tetroxides of osmium or other group 8 elements. Cocrystals between pyridine or pyridine N‐oxide derivatives and osmium tetroxide are characterized through various techniques and rationalized as σ‐hole interactions using DFT calculations and several other computational tools. We propose the term “osme bond” (OmB, Om=Fe, Ru, Os, (Hs)) for naming the noncovalent interactions wherein group 8 elements have the role of the electrophile. The word osme is the transcription of ὀσμή, the ancient Greek word for smell that was used to name the heaviest group 8 element in relation to the smoky odor of its tetroxide.

Methyl groups as widespread Lewis bases in noncovalent interactions

It is well known that, under certain conditions, C(sp³) atoms behave, via their σ-hole, as Lewis acids in tetrel bonding. Here, we show that methyl groups, when bound to atoms less electronegative than carbon, can counterintuitively participate in noncovalent interactions as electron density donors. Thousands of experimental structures are found in which methyl groups behave as Lewis bases to establish alkaline, alkaline earth, triel, tetrel, pnictogen, chalcogen and halogen bonds. Theoretical calculations confirm the high directionality and significant strength of the interactions that arise from a common pattern based on the electron density holes model. Moreover, despite the absence of lone pairs, methyl groups are able to transfer charge from σ bonding orbitals into empty orbitals of the electrophile to reinforce the attractive interaction.

Azide···Oxygen Interaction: A Crystal Engineering Tool for Conformational Locking

We have designed, synthesized and crystallized 36 compounds, each containing an azide group and an oxygen atom separated by three bonds. Crystal structure analysis revealed that each of these molecules adopts a conformation in which the azide and oxygen groups orient syn to each other with a short O ··· N b  contact. Geometry-optimized structures [using  M06-2X/6-311G(d,p) level of theory ] also showed the syn conformation in all 36 of these cases, suggesting that this not merely a crystal packing effect. Quantum topological analysis using Bader's Atoms in Molecules (AIM) theory revealed bond paths and bond critical points (BCP) in these structures suggesting its nature and energetics to be similar to weak hydrogen bonding. The NCI-RDG plot clearly revealed the attractive interaction consisting of electrostatic or dispersive components in all the 36 systems. NBO analysis suggested a weak orbital-relaxation (charge-transfer) contribution of energy for a few  (sp2) O-donor systems. Natural population analysis (NPA) and molecular electrostatic potential mapping (MESP) of these crystal structures further revealed the existence of favorable azide-oxygen interaction.   A CSD search indicated the frequent and consistent occurrence of this interaction and its role dictating the syn conformation of azide and oxygen in molecules where these groups are separated by 2-4 bonds.

Quantification of the electric field inside protein active sites and fullerenes

While electrostatic interactions are exceedingly accountable for biological functions, no simple method exists to directly estimate or measure the electrostatic field in protein active sites. The electrostatic field inside the protein is generally inferred from the shift in the vibrational stretching frequencies of nitrile and thionitrile probes at the active sites through several painstaking and time-consuming experiments like vibrational Stark effect spectroscopy (VSS). Here we present a simple, fast, and reliable methodology, which can efficiently predict the vibrational Stark tuning rates (VSRs) of a large variety of probes within 10% error of the reported experimental data. Our methodology is based on geometry optimization and frequency calculations in the presence of an external electric field to predict the accurate VSR of newly designed nitrile/thionitrile probes. A priori information of VSRs is useful for difficult experiments such as catalytic/enzymatic study and in structural biology. We also applied our methodology successfully to estimate the electric field inside fullerenes and nano-onions, which is encouraging for researchers to adopt it for further applications in materials science and supramolecular chemistry.

Gram-Scale Synthesis of 1,8-Naphthyridines in Water: The Friedlander Reaction Revisited

The products of the Friedlander reaction, i.e., 1,8-naphthyridines, have far-reaching impacts in materials science, chemical biology, and medicine. The reported synthetic methodologies elegantly orchestrate the diverse synthetic routes of naphthyridines but require harsh reaction conditions, organic solvents, and expensive metal catalysts. Here, we introduce gram-scale synthesis of 1,8-naphthyridines in water using an inexpensive and biocompatible ionic liquid (IL) as a catalyst. This is the first-ever report on the synthesis of naphthyridines in water. This is a one-step reaction, and the product separation is relatively easy. The choline hydroxide (ChOH) is used as a metal-free, nontoxic, and water-soluble catalyst. In comparison to other catalysts reported in the literature, ChOH has the advantage of forming an additional hydrogen bond with the reactants, which is the vital step for the reaction to happen in water. Density functional theory (DFT) and noncovalent interaction (NCI) plot index analysis provide the plausible reaction mechanism for the catalytic cycle and confirm that hydrogen bonds with the IL catalyst are pivotal to facilitate the reaction. Molecular docking and molecular dynamics (MD) simulations are also performed to demonstrate the potentialities of the newly synthesized products as drugs. Through MD simulations, it was established that the tetrahydropyrido derivative of naphthyridine (10j) binds to the active sites of the ts3 human serotonin transporter (hSERT) (PDB ID: 6AWO) without perturbing the secondary structure, suggesting that 10j can be a potential preclinical drug candidate for hSERT inhibition and depression treatment.

Anion⋅⋅⋅Anion Coinage Bonds: The Case of Tetrachloridoaurate

Interactions in crystalline tetrachloridoaurates of acetylcholine and dimethylpropiothetine are characterized by Au⋅⋅⋅Cl and Au⋅⋅⋅O short contacts. The former interactions assemble the AuCl4 - units into supramolecular anionic polymers, while the latter interactions append the acetylcholine and propiothetine units to the polymer. The distorted octahedral geometry of the bonding pattern around the gold center is rationalized on the basis of the anisotropic distribution of the electron density, which enables gold to behave as an electrophile (π-hole coinage-bond donor). Computational studies prove that gold atoms in negatively charged species can function as acceptors of electron density. The attractive nature of the Au⋅⋅⋅Cl/O interactions described here complement the known aurophilic bonds involved in gold-centered interactions.

Supramolecular Antiparallel β-Sheet Formation by Tetrapeptides Based on Amyloid Sequence

Self-assembly of short peptides has emerged as an interesting research field for a wide range of applications. Recently, several truncated fragments of long-chain peptides or proteins responsible for different neurodegenerative diseases were studied to understand whether they can mimic the property and function of native peptides or not. It was reported that such a kind of peptide adopts a β-sheet structure in the disease state. It was observed that aromatic amino acid-rich peptide fragments possess a high tendency to adopt a β-sheet conformation. In this article, we are first time reporting the crystal structure of two tetrapeptides: Boc-GAII-OMe (Peptide 1) and Boc-GGVV-OMe (Peptide 2), composed of aliphatic amino acids, and the sequences are similar to the Aβ-peptide fragments Aβ_29-32 and Aβ_37-40 , respectively. In the solid-state, they are self-assembled in an antiparallel β-sheet fashion. The peptide units are connected by the strong amide hydrogen-bonding (N-H···O) interactions. Apart from that, other noncovalent interactions are also present, which help to stabilize the cross-β-sheet arrangement. Interestingly, in the crystal structure of Peptide 1, noncovalent C···C interaction between the electron-deficient carbonyl carbon, and the electron-rich sp³-carbon atom is observed, which is quite rare in the literature. The calculated torsion angles for these peptides are lying in the β-sheet region of the Ramachandran plot. FT-IR studies also indicate the formation of an antiparallel β-sheet structure in the solid-state. Circular dichroism of the peptides in the aqueous solution also suggests the presence of predominantly β-sheet-like conformation in the aqueous solution. Under cross-polarized light, Congo Red stained both peptides showed green-gold color due to birefringence indicating their amyloidogenic nature. This result indicates that the short peptide composed of aliphatic amino acid is capable of forming a β-sheet structure in the absence of aromatic amino acid and also can mimic the function of the native amyloid peptide.

Relative Strengths of a Pnicogen and a Tetrel Bond and Their Mutual Effects upon One Another

The ability of the T and Z atoms of TR₃ZR₂ to engage in a noncovalent interaction with NH₃ is assessed by DFT calculations, where the T atom refers to C, Si, and Ge; Z = As, Sb, and P; and substituents R = H and F. In most instances, the tetrel bond (TB) is both stronger and shorter than the pnicogen bond (ZB). These two bond strengths can be equalized, or preference shifted to the ZB, if F substituents are placed on the Z and H on the T atoms. Employing C as the T atom results in a very weak TB, with the ZB clearly favored energetically. The simultaneous formation of both TB and ZB weakens both, particularly the latter, but both bonds survive intact. Geometric and spectroscopic perturbations of the subunits reflect the two types of noncovalent bonds.

Hydrogen versus tetrel bonds in complexes of 3-oxetanone with water and formaldehyde

The ability and preference of 3-oxetanone to form hydrogen or tetrel bonds have been investigated in its complexes with water and formaldehyde by using Fourier transform microwave spectroscopy complemented with quantum chemical calculations. Different types of interactions and internal dynamics have been observed in the targeted complexes. With water, the ether oxygen of 3-oxetanone is the favoured interaction site forming a classical O-HO hydrogen bond. Quite differently, the carbonyl group of 3-oxetanone plays the dual role as a tetrel donor and a proton acceptor in the 3-oxetanone-formaldehyde complex, featuring the CO tetrel bond and C-HO weak hydrogen bond interactions. Splittings originated from the internal rotation of formaldehyde around its C₂ axis were also observed. The V₂ barrier was estimated to be 375(10) cm^-1 based on Meyer's one-dimensional flexible model. The changes in geometries and electronic densities upon complexation would shed light on the impact of archetype solvent and organic substrate molecules on the reactivity of 3-oxetanone.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

On the importance of RH3C⋯N tetrel bonding interactions in the solid state of a dinuclear zinc complex with a tetradentate Schiff base ligand

The tetrel bonding and π-stacking interactions in a new dinuclear zinc complex using a tetradentate N₂O₂ donor Schiff base have been analysed energetically using DFT calculations and several computational tools.

Origins and properties of the tetrel bond

The tetrel bond (TB) recruits an element drawn from the C, Si, Ge, Sn, Pb family as electron acceptor in an interaction with a partner Lewis base. The underlying principles that explain this attractive interaction are described in terms of occupied and vacant orbitals, total electron density, and electrostatic potential. These principles facilitate a delineation of the factors that feed into a strong TB. The geometric deformation that occurs within the tetrel-bearing Lewis acid monomer is a particularly important issue, with both primary and secondary effects. As a first-row atom of low polarizability, C is a reluctant participant in TBs, but its preponderance in organic and biochemistry make it extremely important that its potential in this regard be thoroughly understood. The IR and NMR manifestations of tetrel bonding are explored as spectroscopy offers a bridge to experimental examination of this phenomenon. In addition to the most common σ-hole type TBs, discussion is provided of π-hole interactions which are a result of a common alternate covalent bonding pattern of tetrel atoms.

Towards developing a criterion to characterize non-covalent bonds: a quantum mechanical study

Chemical bonds are central to chemistry, biology, and allied fields, but still, the criterion to characterize an interaction as a non-covalent bond has not been studied rigorously.