Magnitude and Mechanism of Charge Enhancement of CH··O Hydrogen Bonds

Repurposing FIASMAs against Acid Sphingomyelinase for COVID-19: A Computational Molecular Docking and Dynamic Simulation Approach

Over the past few years, COVID-19 has caused widespread suffering worldwide. There is great research potential in this domain and it is also necessary. The main objective of this study was to identify potential inhibitors against acid sphingomyelinase (ASM) in order to prevent coronavirus infection. Experimental studies revealed that SARS-CoV-2 causes activation of the acid sphingomyelinase/ceramide pathway, which in turn facilitates the viral entry into the cells. The objective was to inhibit acid sphingomyelinase activity in order to prevent the cells from SARS-CoV-2 infection. Previous studies have reported functional inhibitors against ASM (FIASMAs). These inhibitors can be exploited to block the entry of SARS-CoV-2 into the cells. To achieve our objective, a drug library containing 257 functional inhibitors of ASM was constructed. Computational molecular docking was applied to dock the library against the target protein (PDB: 5I81). The potential binding site of the target protein was identified through structural alignment with the known binding pocket of a protein with a similar function. AutoDock Vina was used to carry out the docking steps. The docking results were analyzed and the inhibitors were screened based on their binding affinity scores and ADME properties. Among the 257 functional inhibitors, Dutasteride, Cepharanthine, and Zafirlukast presented the lowest binding affinity scores of −9.7, −9.6, and −9.5 kcal/mol, respectively. Furthermore, computational ADME analysis of these results revealed Cepharanthine and Zafirlukast to have non-toxic properties. To further validate these findings, the top two inhibitors in complex with the target protein were subjected to molecular dynamic simulations at 100 ns. The molecular interactions and stability of these compounds revealed that these inhibitors could be a promising tool for inhibiting SARS-CoV-2 infection.

Cycloaddition Reactivities Analyzed by Energy Decomposition Analyses and the Frontier Molecular Orbital Model

This Account describes our quest to understand and predict organic reactivity, a principal goal of physical and theoretical organic chemistry. The focus is on the development and testing of models for the prediction of cycloaddition reactivities and selectivities. We describe the involvement of the Houk group, and other groups, in the evolution of theoretical models that can achieve ever greater accuracy as well as provide practical heuristic models for understanding and prediction.Is the venerable frontier molecular orbital (FMO) model, the basis of Kenichi Fukui's 1981 Nobel Prize, still useful, or must it be replaced with more advanced models? In particular, models such as Conceptual Density Functional, the Pauli Exclusion Model, and the recent popularity of Electrostatic Potential Plots and Dispersion Energies have not only added to our understanding, but they have also created uncertainty about whether the simple FMO heuristic model has a place in 21st century discussions. This Account addresses this issue and asserts the value of the FMO model.Beginning with brief descriptions of selected models for cycloaddition reactivity starting with early donor-acceptor (nucleophile-electrophile) charge-transfer concepts, this Account reviews Fukui's frontier molecular orbital model, Salem and Klopman's orbital, electrostatic and Pauli repulsion model, the conceptual DFT model by Parr and later by Domingo and others, the recent Houk and Bickelhaupt Distortion/Interaction Activation Strain model, and the Bickelhaupt-Hamlin's Pauli-repulsion lowering model.Computations and analyses of four well-studied Diels-Alder cycloadditions, both normal and inverse electron-demand types, are presented. Most were studied earlier in our published work but are presented here with new insights from calculations with modern methods. Depending on the types of substrates (cycloaddends), the dominant factors controlling reactivity can be orbital interactions, electrostatics and polarization, or Pauli repulsion and dispersion effects, or a combination of all of these.By comparing orbital interactions, especially the frontier molecular orbital interactions, with the other factors that influence reactivity, we show why the FMO model is such a powerful─and theoretically meaningful─heuristic for understanding and predicting reactivity. We also present a method to understand Pauli repulsion effects on activation barriers, ρ(1.1). The use of a new reaction coordinate, the extent of Pauli repulsion along the reaction path, is advocated to emphasize the role of repulsive occupied orbital interactions on reactivity.Fukui's frontier molecular orbital model is effective because FMO interactions parallel all the quantities that influence reactivity. The FMO model continues to provide a practical model to understand and guide experiments.

Unusually Large Effects of Charge-assisted C-H⋅⋅⋅F Hydrogen Bonds to Anionic Fluorine in Organic Solvents: Computational Study of 19 F NMR Shifts versus Thermochemistry

A comparison of computed ¹⁹ F NMR chemical shifts and experiment provides evidence for large specific solvent effects for fluoride-type anions interacting with the σ*(C-H) orbitals in organic solvents like MeCN or CH₂ Cl₂ . We show this for systems ranging from the fluoride ion and the bifluoride ion [FHF]^- to polyhalogen anions [ClF_x ]^- . Discrepancies between computed and experimental shifts when using continuum solvent models like COSMO or force-field-based descriptions like the 3D-RISM-SCF model show specific orbital interactions that require a quantum-mechanical treatment of the solvent molecules. This is confirmed by orbital analyses of the shielding constants, while less negatively charged fluorine atoms (e. g., in [EF₄ ]^- ) do not require such quantum-mechanical treatments to achieve reasonable accuracy. The larger ¹⁹ F solvent shift of fluoride in MeCN compared to water is due to the larger coordination number in the former. These observations are due to unusually strong charge-assisted C-H⋅⋅⋅F^- hydrogen bonds, which manifest beyond some threshold negative natural charge on fluorine of ca. < -0.6 e. The interactions are accompanied by sizable free energies of solvation, in the order F^- ≫[FHF]^- >[ClF₂ ]^- >[ClF₄ ]^- . COSMO-RS solvation free energies tend to moderately underestimate those from the micro-solvated cluster treatment. Red-shifted and intense vibrational C-H stretching bands, potentially accessible in bulk solution, are further spectroscopic finger prints.

Molecular structure, NBO analysis of the hydrogen-bonded interactions, spectroscopic (FT-IR, FT-Raman), drug likeness and molecular docking of the novel anti COVID-2 molecule (2E)-N-methyl-2-[(4-oxo-4H-chromen-3-yl)methylidene]-hydrazinecarbothioamide (Dimer) - quantum chemical approach

Prospective antiviral molecule (2E)-N-methyl-2-[(4-oxo-4H-chromen-3-yl)methylidene]-hydrazinecarbothioamide has been probed using Fourier transform infrared (FTIR), FT-Raman and quantum chemical computations. The geometry equilibrium and natural bond orbital analysis have been carried out with density functional theory employing Becke, 3-parameter, Lee-Yang-Parr method with the 6-311G++(d,p) basis set. The vibrational assignments pertaining to different modes of vibrations have been augmented by normal coordinate analysis, force constant and potential energy distributions. Drug likeness and oral activity have been carried out based on Lipinski's rule of five. The inhibiting potency of 2(2E)-methyl-2-[(4-oxo-4H-chromen-3-yl)methylidene]-hydrazinecarbothioamide has been investigated by docking simulation against SARS-CoV-2 protein. The optimized geometry shows a planar structure between the chromone and the side chain. Differences in the geometries due to the substitution of the electronegative atom and intermolecular contacts due to the chromone and hydrazinecarbothioamide were analyzed. NBO analysis confirms the presence of two strong stable hydrogen bonded NH⋯O intermolecular interactions and two weak hydrogen bonded CH⋯O interactions. The red shift in NH stretching frequency exposed from IR substantiates the formation of NH⋯O intermolecular hydrogen bond and the blue shift in CH stretching frequency substantiates the formation of CH⋯O intermolecular hydrogen bond. Drug likeness, absorption, distribution, metabolism, excretion and toxicity property gives an idea about the pharmacokinetic properties of the title molecule. The binding energy of the nonbonding interaction with Histidine 41 and Cysteine 145, present a clear view that 2(2E)-methyl-2-[(4-oxo-4H-chromen-3-yl)methylidene]-hydrazinecarbothioamide can irreversibly interact with SARS-CoV-2 protease.

Inspiration for revival of old drugs: improving solubility and avoiding hygroscopicity of pipemidic acid by forming two pharmaceutical salts based on charge-assisted hydrogen bond recognitions

Improving solubility and avoiding hygroscopicity of pipemidic acid by forming pharmaceutical salts based on CAHBs.

The Distinct Conformational Landscapes of 4S-Substituted Prolines That Promote an endo Ring Pucker

4-Substitution on proline directly impacts protein main chain conformational preferences. The structural effects of N-acyl substitution and of 4-substitution were examined by NMR spectroscopy and X-ray crystallography on minimal molecules with a proline 4S-nitrobenzoate. The effects of N-acyl substitution on conformation were attenuated in the 4S-nitrobenzoate context, due to the minimal role of the n→π* interaction in stabilizing extended conformations. By X-ray crystallography, an extended conformation was observed for most molecules. The formyl derivative adopted a δ conformation that is observed at the i+2 position of β-turns. Computational analysis indicated that the structures observed crystallographically represent the inherent conformational preferences of 4S-substituted prolines with electron-withdrawing 4-position substituents. The divergent conformational preferences of 4R- and 4S-substituted prolines suggest their wider structure-specific application in molecular design. In particular, the proline endo ring pucker favored by 4S-substituted prolines uniquely promotes the δ conformation [(ϕ, ψ) ≈(-80°, 0°)] found in β-turns. In contrast to other acyl capping groups, the pivaloyl group strongly promoted trans amide bond and polyproline II helix conformation, with a close n→π* interaction in the crystalline state, despite the endo ring pucker, suggesting its special capabilities in promoting compact conformations in ϕ due to its strongly electron-donating character.

A New Microbial Pathway for Organophosphonate Degradation Catalyzed by Two Previously Misannotated Non-Heme-Iron Oxygenases

The assignment of biochemical functions to hypothetical proteins is challenged by functional diversification within many protein structural superfamilies. This diversification, which is particularly common for metalloenzymes, renders functional annotations that are founded solely on sequence and domain similarities unreliable and often erroneous. Definitive biochemical characterization to delineate functional subgroups within these superfamilies will aid in improving bioinformatic approaches for functional annotation. We describe here the structural and functional characterization of two non-heme-iron oxygenases, TmpA and TmpB, which are encoded by a genomically clustered pair of genes found in more than 350 species of bacteria. TmpA and TmpB are functional homologues of a pair of enzymes (PhnY and PhnZ) that degrade 2-aminoethylphosphonate but instead act on its naturally occurring, quaternary ammonium analogue, 2-(trimethylammonio)ethylphosphonate (TMAEP). TmpA, an iron(II)- and 2-(oxo)glutarate-dependent oxygenase misannotated as a γ-butyrobetaine (γbb) hydroxylase, shows no activity toward γbb but efficiently hydroxylates TMAEP. The product, ( R)-1-hydroxy-2-(trimethylammonio)ethylphosphonate [( R)-OH-TMAEP], then serves as the substrate for the second enzyme, TmpB. By contrast to its purported phosphohydrolytic activity, TmpB is an HD-domain oxygenase that uses a mixed-valent diiron cofactor to enact oxidative cleavage of the C-P bond of its substrate, yielding glycine betaine and phosphate. The high specificities of TmpA and TmpB for their N-trimethylated substrates suggest that they have evolved specifically to degrade TMAEP, which was not previously known to be subject to microbial catabolism. This study thus adds to the growing list of known pathways through which microbes break down organophosphonates to harvest phosphorus, carbon, and nitrogen in nutrient-limited niches.

N+-C-H···O Hydrogen bonds in protein-ligand complexes

In the context of drug design, C-H···O hydrogen bonds have received little attention so far, mostly because they are considered weak relative to other noncovalent interactions such as O-H···O hydrogen bonds, π/π interactions, and van der Waals interactions. Herein, we demonstrate the significance of hydrogen bonds between C-H groups adjacent to an ammonium cation and an oxygen atom (N⁺-C-H···O hydrogen bonds) in protein-ligand complexes. Quantum chemical calculations revealed details on the strength and geometrical requirements of these N⁺-C-H···O hydrogen bonds, and a subsequent survey of the Protein Data Bank (PDB) based on these criteria suggested that numerous protein-ligand complexes contain such N⁺-C-H···O hydrogen bonds. An ensuing experimental investigation into the G9a-like protein (GLP)-inhibitor complex demonstrated that N⁺-C-H···O hydrogen bonds affect the activity of the inhibitors against the target enzyme. These results should provide the basis for the use of N⁺-C-H···O hydrogen bonds in drug discovery.

Crystallographic and Computational Characterization of Methyl Tetrel Bonding in S-Adenosylmethionine-Dependent Methyltransferases

Tetrel bonds represent a category of non-bonding interaction wherein an electronegative atom donates a lone pair of electrons into the sigma antibonding orbital of an atom in the carbon group of the periodic table. Prior computational studies have implicated tetrel bonding in the stabilization of a preliminary state that precedes the transition state in S_N2 reactions, including methyl transfer. Notably, the angles between the tetrel bond donor and acceptor atoms coincide with the prerequisite geometry for the S_N2 reaction. Prompted by these findings, we surveyed crystal structures of methyltransferases in the Protein Data Bank and discovered multiple instances of carbon tetrel bonding between the methyl group of the substrate S-adenosylmethionine (AdoMet) and electronegative atoms of small molecule inhibitors, ions, and solvent molecules. The majority of these interactions involve oxygen atoms as the Lewis base, with the exception of one structure in which a chlorine atom of an inhibitor functions as the electron donor. Quantum mechanical analyses of a representative subset of the methyltransferase structures from the survey revealed that the calculated interaction energies and spectral properties are consistent with the values for bona fide carbon tetrel bonds. The discovery of methyl tetrel bonding offers new insights into the mechanism underlying the S_N2 reaction catalyzed by AdoMet-dependent methyltransferases. These findings highlight the potential of exploiting these interactions in developing new methyltransferase inhibitors.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

Opposite-charge repulsive cation and anion pair cooperative organocatalysis in ring-opening polymerization

A new strained ion pair catalysis was proposed in ring-opening polymerization.

New tricks for the glycyl radical enzyme family

Glycyl radical enzymes (GREs) are important biological catalysts in both strict and facultative anaerobes, playing key roles both in the human microbiota and in the environment. GREs contain a backbone glycyl radical that is post-translationally installed, enabling radical-based mechanisms. GREs function in several metabolic pathways including mixed acid fermentation, ribonucleotide reduction and the anaerobic breakdown of the nutrient choline and the pollutant toluene. By generating a substrate-based radical species within the active site, GREs enable C-C, C-O and C-N bond breaking and formation steps that are otherwise challenging for nonradical enzymes. Identification of previously unknown family members from genomic data and the determination of structures of well-characterized GREs have expanded the scope of GRE-catalyzed reactions as well as defined key features that enable radical catalysis. Here, we review the structures and mechanisms of characterized GREs, classifying members into five categories. We consider the open questions about each of the five GRE classes and evaluate the tools available to interrogate uncharacterized GREs.

Quantum Chemical Modeling of Hydrogen Bonding in Ionic Liquids

Hydrogen bonding (H-bonding) is an important and very general phenomenon. H-bonding is part of the basis of life in DNA, key in controlling the properties of water and ice, and critical to modern applications such as crystal engineering, catalysis applications, pharmaceutical and agrochemical development. H-bonding also plays a significant role for many ionic liquids (IL), determining the secondary structuring and affecting key physical parameters. ILs exhibit a particularly diverse and wide range of traditional as well as non-standard forms of H-bonding, in particular the doubly ionic H-bond is important. Understanding the fundamental nature of the H-bonds that form within ILs is critical, and one way of accessing this information, that cannot be recovered by any other computational method, is through quantum chemical electronic structure calculations. However, an appropriate method and basis set must be employed, and a robust procedure for determining key structures is essential. Modern generalised solvation models have recently been extended to ILs, bringing both advantages and disadvantages. QC can provide a range of information on geometry, IR and Raman spectra, NMR spectra and at a more fundamental level through analysis of the electronic structure.

Torsional and Electronic Factors Control the C-H⋅⋅⋅O Interaction

The precise role of non‐conventional hydrogen bonds such as the C−H⋅⋅⋅O interaction in influencing the conformation of small molecules remains unresolved. Here we survey a series of β‐turn mimetics using X‐ray crystallography and NMR spectroscopy in conjunction with quantum calculation, and conclude that favourable torsional and electronic effects are important for the population of states with conformationally influential C−H⋅⋅⋅O interactions. Our results also highlight the challenge in attempting to deconvolute a myriad of interdependent noncovalent interactions in order to focus on the contribution of a single one. Within a small molecule that is designed to resemble the complexity of the environment within peptides and proteins, the interplay of different steric burdens, hydrogen‐acceptor/‐donor properties and rotational profiles illustrate why unambiguous conclusions based solely on NMR chemical shift data are extremely challenging to rationalize.

Molecular Basis of C-N Bond Cleavage by the Glycyl Radical Enzyme Choline Trimethylamine-Lyase

Deamination of choline catalyzed by the glycyl radical enzyme choline trimethylamine-lyase (CutC) has emerged as an important route for the production of trimethylamine, a microbial metabolite associated with both human disease and biological methane production. Here, we have determined five high-resolution X-ray structures of wild-type CutC and mechanistically informative mutants in the presence of choline. Within an unexpectedly polar active site, CutC orients choline through hydrogen bonding with a putative general base, and through close interactions between phenolic and carboxylate oxygen atoms of the protein scaffold and the polarized methyl groups of the trimethylammonium moiety. These structural data, along with biochemical analysis of active site mutants, support a mechanism that involves direct elimination of trimethylamine. This work broadens our understanding of radical-based enzyme catalysis and will aid in the rational design of inhibitors of bacterial trimethylamine production.

Sulfur-Oxygen Chalcogen Bonding Mediates AdoMet Recognition in the Lysine Methyltransferase SET7/9

Recent studies have demonstrated that carbon-oxygen (CH···O) hydrogen bonds have important roles in S-adenosylmethionine (AdoMet) recognition and catalysis in methyltransferases. Here, we investigate noncovalent interactions that occur between the AdoMet sulfur cation and oxygen atoms in methyltransferase active sites. These interactions represent sulfur-oxygen (S···O) chalcogen bonds in which the oxygen atom donates a lone pair of electrons to the σ antibonding orbital of the AdoMet sulfur atom. Structural, biochemical, and computational analyses of an asparagine mutation in the lysine methyltransferase SET7/9 that abolishes AdoMet S···O chalcogen bonding reveal that this interaction enhances substrate binding affinity relative to the product S-adenosylhomocysteine. Corroborative quantum mechanical calculations demonstrate that sulfonium systems form strong S···O chalcogen bonds relative to their neutral thioether counterparts. An inspection of high-resolution crystal structures reveals the presence of AdoMet S···O chalcogen bonding in different classes of methyltransferases, illustrating that these interactions are not limited to SET domain methyltransferases. Together, these results demonstrate that S···O chalcogen bonds contribute to AdoMet recognition and can enable methyltransferases to distinguish between substrate and product.

Doubly ionic hydrogen bond interactions within the choline chloride–urea deep eutectic solvent

Computational analysis indicates flexibility and diversity in the hydrogen bonding, but limited charge delocalisation, within the choline chloride–urea eutectic.

Synthesis of Bi2Te3and (BixSb1−x)2Te3nanoparticles using the novel IL [C4mim]3[Bi3I12]

[C₄mim]₃[Bi₃I₁₂] is a promising Bi-source for the ionothermal synthesis of binary (Bi₂Te₃) and ternary tetradymite-type nanoparticles (Bi_xSb_1−x)₂Te₃(x= 0.25, 0.5, 0.75) in ionic liquid.

Quantification of Conventional and Nonconventional Charge-Assisted Hydrogen Bonds in the Condensed and Gas Phases

Charge-assisted hydrogen bonds (CAHBs) play critical roles in many systems from biology through to materials. In none of these areas has the role and function of CAHBs been explored satisfactorily because of the lack of data on the energy of CAHBs in the condensed phases. We have, for the first time, quantified three types of CAHBs in both the condensed and gas phases for 1-(2'-hydroxylethyl)-3-methylimidazolium acetate ([C2OHmim][OAc]). The energy of conventional OH···[OAc](-) CAHBs is ∼10 kcal·mol(-1), whereas nonconventional C(sp2)H···[OAc](-) and C(sp3)H···[OAc](-) CAHBs are weaker by ∼5-7 kcal·mol(-1). In the gas phase, the strength of the nonconventional CAHBs is doubled, whereas the conventional CAHBs are strengthened by <20%. The influence of cooperativity effects on the ability of the [OAc](-) anion to deprotonate the imidazolium cation is evaluated. The ability to quantify CAHBs in the condensed phase on the basis of easier accessible gas-phase estimates is highlighted.

Hydrogen bonding in ionic liquids

Ionic liquids (IL) and hydrogen bonding (H-bonding) are two diverse fields for which there is a developing recognition of significant overlap. Doubly ionic H-bonds occur when a H-bond forms between a cation and anion, and are a key feature of ILs. Doubly ionic H-bonds represent a wide area of H-bonding which has yet to be fully recognised, characterised or explored. H-bonds in ILs (both protic and aprotic) are bifurcated and chelating, and unlike many molecular liquids a significant variety of distinct H-bonds are formed between different types and numbers of donor and acceptor sites within a given IL. Traditional more neutral H-bonds can also be formed in functionalised ILs, adding a further level of complexity. Ab initio computed parameters; association energies, partial charges, density descriptors as encompassed by the QTAIM methodology (ρBCP), qualitative molecular orbital theory and NBO analysis provide established and robust mechanisms for understanding and interpreting traditional neutral and ionic H-bonds. In this review the applicability and extension of these parameters to describe and quantify the doubly ionic H-bond has been explored. Estimating the H-bonding energy is difficult because at a fundamental level the H-bond and ionic interaction are coupled. The NBO and QTAIM methodologies, unlike the total energy, are local descriptors and therefore can be used to directly compare neutral, ionic and doubly ionic H-bonds. The charged nature of the ions influences the ionic characteristics of the H-bond and vice versa, in addition the close association of the ions leads to enhanced orbital overlap and covalent contributions. The charge on the ions raises the energy of the Ylp and lowers the energy of the X-H σ* NBOs resulting in greater charge transfer, strengthening the H-bond. Using this range of parameters and comparing doubly ionic H-bonds to more traditional neutral and ionic H-bonds it is clear that doubly ionic H-bonds cover the full range of weak through to very strong H-bonds.

Manipulating unconventional CH-based hydrogen bonding in a methyltransferase via noncanonical amino acid mutagenesis

Recent studies have demonstrated that the active sites of S-adenosylmethionine (AdoMet)-dependent methyltransferases form strong carbon-oxygen (CH···O) hydrogen bonds with the substrate's sulfonium group that are important in AdoMet binding and catalysis. To probe these interactions, we substituted the noncanonical amino acid p-aminophenylalanine (pAF) for the active site tyrosine in the lysine methyltransferase SET7/9, which forms multiple CH···O hydrogen bonds to AdoMet and is invariant in SET domain enzymes. Using quantum chemistry calculations to predict the mutation's effects, coupled with biochemical and structural studies, we observed that pAF forms a strong CH···N hydrogen bond to AdoMet that is offset by an energetically unfavorable amine group rotamer within the SET7/9 active site that hinders AdoMet binding and activity. Together, these results illustrate that the invariant tyrosine in SET domain methyltransferases functions as an essential hydrogen bonding hub and cannot be readily substituted by residues bearing other hydrogen bond acceptors.

C–H/O interactions of nucleic bases with a water molecule: a crystallographic and quantum chemical study

The C–H/O interactions of nucleic bases are substantially stronger than the C–H/O interactions of benzene and pyridine. These results can be very important for molecular recognition of DNA and RNA.