Towards an unified hydrogen-bond theory

Mechanistic and Structural Insights on Difluoromethyl-1,3,4-oxadiazole Inhibitors of HDAC6

Histone deacetylase 6 (HDAC6) is increasingly recognized for its potential in targeted disease therapy. This study delves into the mechanistic and structural nuances of HDAC6 inhibition by difluoromethyl-1,3,4-oxadiazole (DFMO) derivatives, a class of non-hydroxamic inhibitors with remarkable selectivity and potency. Employing a combination of nuclear magnetic resonance (NMR) spectroscopy and liquid chromatography-mass spectrometry (LC-MS) kinetic experiments, comprehensive enzymatic characterizations, and X-ray crystallography, we dissect the intricate details of the DFMO-HDAC6 interaction dynamics. More specifically, we find that the chemical structure of a DMFO and the binding mode of its difluoroacetylhydrazide derivative are crucial in determining the predominant hydrolysis mechanism. Our findings provide additional insights into two different mechanisms of DFMO hydrolysis, thus contributing to a better understanding of the HDAC6 inhibition by oxadiazoles in disease modulation and therapeutic intervention.

Thermoswitchable catalysis to inhibit and promote plastic flow in vitrimers

Acid-base catalysis is a common strategy to induce covalent bond exchanges in dynamic polymer networks. Strong acids or strong bases can promote rapid network rearrangements, and are simultaneously preferred catalysts for chemical reactions where maximum efficiency at the lowest possible temperature is aimed for. However, within the context of dynamic polymer networks, the incorporation of highly active catalysts can negatively affect the longer term application potential. Network dynamicity can diminish through catalyst ageing or quenching and highly active catalysts may prematurely activate bond exchanges, leading to dimensional instability and thus low creep resistance of the polymer networks. Herein, we present several examples where we explicitly explored weak acids (carboxylic acids) as catalysts for dynamic bond exchanges, using vinylogous urethanes (VU) as a well-understood protic acid catalysed vitrimer chemistry. Surprisingly, we have found that the sought-after long-term stability offered by a weak acid does not necessarily bring lower activity at high temperature. In fact, the weak acids show a remarkable thermoswitchable catalytic behaviour, going from an inactive hydrogen bonded state to an active state where the polymer matrix is protonated, with a profound impact on the network reactivity and rheology. Carboxylic acids with different electronic or steric environments show clear reactivity trends and their fine-tuning resulted in the most thermally responsive VU vitrimers studied to date. Our findings point out that catalyst choice and design for vitrimers is only poorly informed by catalyst performance in more traditional chemical reactions (in solvent), and that a more tailored catalyst design holds great promise for the field of vitrimers.

Effects of substitution and conjugation on photophysical properties of ESIPT-based fluorophores with the core of 4-aminophthalimide

4-Aminophthalimide is a highly fluorescent signaling unit with excellent photophysical properties and wide application foregrounds. Based on this, a range of theoretical investigations are conducted on the fluorescent probe (E)-5-((2-hydroxybenzylidene) amino) isoindoline-1, 3-dione (HID) with the core of 4-aminophthalimide using density functional theory (DFT) and time-containing density functional theory (TD-DFT) methods in this paper. The optimized configurations, vertical excitation and emission energies, electronic characteristics and excited-state intramolecular proton transfer (ESIPT) behaviors of the probe HID are discussed in detail. Furthermore, to enhance the luminescent properties of HID, five novel compounds have been designed based on the structure of HID by introducing amino, methoxy and naphthalene groups (-NH2, -OMe and C10H8). Our work thoroughly explores how the property and position of substituents and conjugation affect photophysical characteristics and ESIPT processes. We find that the ESIPT dynamics can be modulated by the substitution and conjugation effects. Specifically, the introduction of amino and methoxy groups at the ortho-position and the introduction of the naphthalene group promote the ESIPT behavior of HID1, whereas the introduction of amino and methoxy groups at the meta-position exhibits the contrary impact. Therefore, we boldly infer that the introduction of electron-donating groups in the ortho-position and the introduction of the conjugated group make the ESIPT process more effortless to occur, whereas the introduction of substituents with opposing natures in the meta-position makes the ESIPT process more difficult to occur. In addition, the ionization potentials (IP), electron affinities (EA) and reorganization energies (λh and λe) of molecules are calculated to assess their potential as luminescent materials. Our work not only reveals the luminescence and ESIPT mechanism of the probe HID1, but also proposes to modulate the ESIPT process through the substitution and conjugation effects. In particular, the designed molecules have better photoelectric properties as a result of their red-shifted absorption and fluorescence spectra, smaller energy gaps, larger transferred charges and greater charge transferred distances, which offers some valuable ideas for the experimental development of more efficient organic luminescent materials with ESIPT properties.

2,2'-Dithiobispyrazine: about the disulfide bond

X-ray diffraction studies reveal that pyrazine-2-thiol undergoes condensation to 2,2'-dithiobispyrazine [systematic name: 2-(pyrazin-2-yldisulfanyl)pyrazine], C8H6N4S2 (I), under aerial conditions. In the molecule of I, the pyrazine rings are arranged in an almost perpendicular manner, with an absolute value of the C-S-S-C torsion angle of -91.45 (6)°. A search in the Cambridge Structural Database confirmed that such a conformation is typical for disulfide compounds. Three different rotamers of disulfide I were studied using quantum theoretical studies. The rotamer of lowest energy was observed in the crystalline state in the structure stabilized by hydrogen-bond, chalcogen-bond and stacking interactions. Further quantum chemical computations confirm that 2,2'-dithiobispyrazine can react according to the SN2 mechanism.

Experimental charge density and topological analysis of tetraaquabis(hydrogenmaleato)nickel(II): a comparison with Hirshfeld atom refinement

Experimental charge density analysis is conducted on the coordination compound tetraaquabis(hydrogenmaleato)nickel(II), which exhibits a short intramolecular hydrogen bond. Through topological analysis, the nature of Ni-O bonds is concluded to be intermediate between ionic and covalent, but mainly presenting an ionic character, while the short hydrogen bond is classified as covalent in nature. The compound was also analysed after Hirshfeld atom refinement performed using NoSpherA2. A topological analysis was conducted on the molecular wavefunction and the results are compared with those obtained from experiment. In general, there is good agreement between the refinements, and the chemical bonds involving H atoms are in better agreement with what is expected from neutron data after HAR than they are after multipole refinement.

Identification of a Novel SSTR3 Full Agonist for the Treatment of Nonfunctioning Pituitary Adenomas

Somatostatin receptor (SSTR) agonists have been extensively used for treating neuroendocrine tumors. Synthetic therapeutic agonists showing selectivity for SSTR2 (Octreotide) or for SSTR2 and SSTR5 (Pasireotide) have been approved for the treatment of patients with acromegaly and Cushing's syndrome, as their pituitary tumors highly express SSTR2 or SSTR2/SSTR5, respectively. Nonfunctioning pituitary adenomas (NFPAs), which express high levels of SSTR3 and show only modest response to currently available SSTR agonists, are often invasive and cannot be completely resected, and therefore easily recur. The aim of the present study was the evaluation of ITF2984, a somatostatin analog and full SSTR3 agonist, as a new potential treatment for NFPAs. ITF2984 shows a 10-fold improved affinity for SSTR3 compared to Octreotide or Pasireotide. Molecular modeling and NMR studies indicated that the higher affinity for SSTR3 correlates with a higher stability of a distorted β-I turn in the cyclic peptide backbone. ITF2984 induces receptor internalization and phosphorylation, and triggers G-protein signaling at pharmacologically relevant concentrations. Furthermore, ITF2984 displays antitumor activity that is dependent on SSTR3 expression levels in the MENX (homozygous mutant) NFPA rat model, which closely recapitulates human disease. Therefore, ITF2984 may represent a novel therapeutic option for patients affected by NFPA.

Role of Non-Covalent Interactions in Carbonic Anhydrase I-Topiramate Complex Based on QM/MM Approach

Carbonic anhydrase (CA) I with a Topiramate (TPM) complex was investigated on the basis of a Quantum Mechanics/Molecular Mechanics (QM/MM) approach. The QM part was treated using Density Functional Theory (DFT) while the MM was simulated using Amberff14SB and GAFF force fields. In addition, the TIP3P model was applied to reproduce the polar environment influence on the studied complex. Next, three snapshots (after 5 ps, 10 ps, and 15 ps of the simulation time) were taken from the obtained trajectory to provide an insight into the non-covalent interactions present between the ligand and binding pocket of the protein. Our special attention was devoted to the binding site rearrangement, which is known in the literature concerning the complex. This part of the computations was performed using ωB97X functional with Grimme D3 dispersion corrections as well as a Becke-Johnson damping function (D3-BJ). Two basis sets were applied: def2-SVP (for larger models) and def2-TZVPD (for smaller models), respectively. In order to detect and describe non-covalent interactions between amino acids of the binding pocket and the ligand, Independent Gradient Model based on Hirshfeld partitioning (IGMH), Interaction Region Indicator (IRI), Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Orbitals (NBO) methods were employed. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was applied for energy decomposition between the ligand and protein. It was found that during the simulation time, the ligand position in the binding site was preserved. Nonetheless, amino acids interacting with TPM were exchanging during the simulation, thus showing the binding site reorganization. The energy partitioning revealed that dispersion and electrostatics are decisive factors that are responsible for the complex stability.

Pyridinium-Inspired Organocatalysts for Carbon Dioxide Fixation: A Density Functional Theory Inspection

The current project aims to apply the virtues of minimalism to examine the catalytic ability of commercially organic compounds of small chemical structures to catalyze the coupling reaction between carbon dioxide and propylene oxide (PO) under mild conditions. The proposed catalysts are pyridinium iodide (A), 2-hydroxypyridinium iodide (B), and piperidinium iodide (C), where their structure is based on cooperative acidic and nucleophilic motifs. The quantum chemistry model, M062X-D3/def2-TZVP//M062X-D3/def2-SVPP, was used to understand the reaction mechanism and the catalytic performance. Since the coupling reaction was performed under excess PO, we proposed that PO serves as a reactant and solvent. Therefore, calculations were performed in gas and liquid phases for comparison. The findings indicated that the rate-determining step depends on the chemical structure of the catalyst and whether the phase is a gas or liquid phase. In general, modeling in the liquid phase produces potential energy surfaces of lower energy barriers. The noncovalent interactions reflect the role of hydrogen bonding in controlling the kinetic behavior of the coupling reaction. Based on the finding, catalyst A is the best candidate for transforming CO₂ into cyclic carbonates.

Various hydrogen bonds make different fates of pharmaceutical contaminants on oxygen-rich nanomaterials

Various hydrogen bonds, especially charge-assisted hydrogen bond (CAHB), is considered as one of vital mechanisms affecting the environmental behavior and risk of pharmaceutical contaminants (PCs). Herein the sorption/desorption of three PCs including clofibric acid (CA), acetaminophen (ACT), and sulfamerazine (SMZ) on three Oxygen-rich (O-rich) nanoparticles (nano-silica: Nano-SiO₂, nano-alumina: Nano-Al₂O₃, and oxidized carbon nanotubes: O-CNTs) were investigated to explore the effect of various hydrogen bonds with different strengths on environmental behaviors of PCs. The results indicated that although solvent-assisted CAHB, solvent-uninvolved CAHB, and ordinary hydrogen bond (OHB) all played a crucial role in sorption of PCs on three O-rich nanomaterials, they showed significantly different effects on the desorption behaviors of PCs from three sorbents. Compared with OHB (hysteresis index ≤0.0766), the stronger CAHB (hysteresis index ≥0.1981) between PCs and O-rich nanoparticles having comparable pK_a with PCs, caused obvious desorption hysteresis of PCs, resulting in their better immobilization on O-rich nanomaterials. The FTIR characterization found that both solvent-assisted and solvent-uninvolved CAHB formation resulted in a new characteristic peak appeared in the high frequency (3660 cm^-1 for Nano-SiO₂, 3730 cm^-1 for Nano-Al₂O₃, and 3780 cm^-1 for O-CNTs). Also, density functional theory (DFT) calculation verified that the smaller |ΔpK_a| between PCs and O-rich sorbents, the shorter bond length, and the larger bond angle resulted in the stronger hydrogen bond formed, thereby leading to the greater immobilization of PCs. These results provide in-depth understanding of the environmental behavior and risk of PCs, and light new idea for designed materials to control PCs pollution in the environment.

Exploring the Dynamical Nature of Intermolecular Hydrogen Bonds in Benzamide, Quinoline and Benzoic Acid Derivatives

The hydrogen bonds properties of 2,6-difluorobenzamide, 5-hydroxyquinoline and 4-hydroxybenzoic acid were investigated by Car-Parrinello and path integral molecular dynamics (CPMD and PIMD), respectively. The computations were carried out in vacuo and in the crystalline phase. The studied complexes possess diverse networks of intermolecular hydrogen bonds (N-H…O, O-H…N and O-H…O). The time evolution of hydrogen bridges gave a deeper insight into bonds dynamics, showing that bridged protons are mostly localized on the donor side; however, the proton transfer phenomenon was registered as well. The vibrational features associated with O-H and N-H stretching were analyzed on the basis of the Fourier transform of the atomic velocity autocorrelation function. The spectroscopic effects of hydrogen bond formation were studied. The PIMD revealed quantum effects influencing the hydrogen bridges providing more accurate free energy sampling. It was found that the N…O or O…O interatomic distances decreased (reducing the length of the hydrogen bridge), while the O-H or N-H covalent bond was elongated, which led to the increase in the proton sharing. Furthermore, Quantum Theory of Atoms in Molecules (QTAIM) was used to give insight into electronic structure parameters. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was employed to estimate the energy contributions to the interaction energy of the selected dimers.

Visible-Light Organophotoredox-Mediated [3 + 2] Cycloaddition of Arylcyclopropylamine with Structurally Diverse Olefins for the Construction of Cyclopentylamines and Spiro[4.n] Skeletons

We developed a visible-light-mediated [3 + 2] cycloaddition of arylcyclopropylamine with structurally diverse olefins using QXPT-NPh as a highly efficient organic photoredox catalyst. We first achieved the use of various alkyl-substituted alkenes in intermolecular [3 + 2] cycloadditions with cyclopropylamine. We also developed a general and efficient strategy for the construction of structurally diverse cyclopentane-based spiro[4.n] skeletons with 1,3-difunctional groups, which broadly exist in natural products and synthetic molecules. Furthermore, we proposed a hydrogen-bond mode between the arylcyclopropylamine and the photocatalyst QXPT-NPh.

Effects of substitution and conjugation on ESIPT behavior of Schiff base derivatives

The excited-state proton transfer (ESIPT) behavior of organic fluorophores has been of great interest due to their unique photophysical properties. In this work, we have focused on the excited state kinetic behavior of four Schiff base organic molecules (i.e. CPMP, CPMMP, CPMDP, and CPMN) in acetonitrile solvents. The electron-donating of substituents and conjugation effects on the photophysical properties and ESIPT process of the Schiff base derivatives are investigated by theoretical methods. The results show that the hydrogen bonds are all enhanced in the excited states, which could provide the impetus for the ESIPT process. To further reveal the reaction process of ESIPT, we have scanned the potential energy curves of the ESIPT process and compared the potential barriers. It is found that the stronger the substituents give electrons and the conjugation effects the more favorable the excited state proton transfer (ESIPT). In the meantime, this study paves the way for the development of new Schiff base materials based on ESIPT.

Effect of pH-Dependent Homo/Heteronuclear CAHB on Adsorption and Desorption Behaviors of Ionizable Organic Compounds on Carbonaceous Materials

Herein, the adsorption/desorption behaviors of benzoic acid (BA) and phthalic acid (PA) on three functionalized carbon nanotubes (CNTs) at various pH were investigated, and the charge-assisted H-bond (CAHB) was verified by DFT and FTIR analyses to play a key role. The results indicated that the adsorption order of BA and PA on CNTs was different from K_ow of that at pH 2.0, 4.0, and 7.0 caused by the CAHB interaction. The strength of homonuclear CAHB (≥78.96 kJ·mol^-1) formed by BA/PA on oxidized CNTs is stronger than that of heteronuclear CAHB formed between BA/PA and amino-functionalized CNTs (≤51.66 kJ·mol^-1). Compared with the heteronuclear CAHB (Hysteresis index, HI ≥ 1.47), the stronger homonuclear CAHB leads to clearly desorption hysteresis (HI ≥ 3.51). Additionally, the contribution of homonuclear CAHB (≥52.70%) was also greater than that of heteronuclear CAHB (≤45.79%) at pH 7.0. These conclusions were further confirmed by FTIR and DFT calculation, and the crucial evidence of CAHB formation in FTIR was found. The highlight of this work is the identification of the importance and difference of pH-dependent homonuclear/heteronuclear CAHB on the adsorption and desorption behaviors of ionizable organic compounds on carbonaceous materials, which can provide a deeper understanding for the removal of ionizable organic compounds by designed carbonaceous materials.

Core fucose-specific Pholiota squarrosa lectin (PhoSL) as a potent broad-spectrum inhibitor of SARS-CoV-2 infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein (S protein) is highly N-glycosylated, and a "glycan shield" is formed to limit the access of other molecules; however, a small open area coincides with the interface to the host's receptor and also neutralising antibodies. Most of the variants of concern have mutations in this area, which could reduce the efficacy of existing antibodies. In contrast, N-glycosylation sites are relatively invariant, and some are essential for infection. Here, we observed that the S proteins of the ancestral (Wuhan) and Omicron strains bind with Pholiota squarrosa lectin (PhoSL), a 40-amino-acid chemically synthesised peptide specific to core-fucosylated N-glycans. The affinities were at a low nanomolar level, which were ~ 1000-fold stronger than those between PhoSL and the core-fucosylated N-glycans at the micromolar level. We demonstrated that PhoSL inhibited infection by both strains at similar submicromolar levels, suggesting its broad-spectrum effect on SARS-CoV-2 variants. Cryogenic electron microscopy revealed that PhoSL caused an aggregation of the S protein, which was likely due to the multivalence of both the trimeric PhoSL and S protein. This characteristic is likely relevant to the inhibitory mechanism. Structural modelling of the PhoSL-S protein complex indicated that PhoSL was in contact with the amino acids of the S protein, which explains the enhanced affinity with S protein and also indicates the significant potential for developing specific binders by the engineering of PhoSL.

Adsorption Mechanism of High-Concentration Ammonium by Chinese Natural Zeolite with Experimental Optimization and Theoretical Computation

Natural zeolite, as an abundant aluminosilicate mineral with a hierarchically porous structure, has a strong affinity to ammonium in solutions. Adsorption mechanism of high-concentration ammonium (1000~4000 mg-N/L) in an aqueous solution without pH adjustment onto Chinese natural zeolite with the dosage of 5 g/L was revealed by the strategy of experimental optimization integrated with Molecular Dynamics (MD) simulation, and found the maximum ammonium adsorption capacity was 26.94 mg/g. The adsorption kinetics and isotherm analysis showed that this adsorption process fitted better with descriptions of the pseudo-second-order kinetics and Freundlich model. The theoretical calculations and infrared-spectrum characterization results verified the existence of hydrogen bonds and chemisorption. Therefore, the adsorption mechanism by natural zeolites of high-concentration NH4+ is defined as a process under the joint influence of multiple effects, which is mainly promoted by the synergy of the ion exchange process, electrostatic attraction, and chemisorption. Meanwhile, the hydrogen bond also plays an auxiliary role in this efficient adsorption. This study presents important theoretical significance for enriching the mechanism of zeolites adsorbing NH4+ from water, and provides reference and theoretical guidance for further exploring the potential application of natural zeolites.

Hepatitis C virus NS3/4A inhibitors and other drug-like compounds as covalent binders of SARS-CoV-2 main protease

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2), which causes coronavirus disease 2019 (COVID-19), threatens global public health. The world needs rapid development of new antivirals and vaccines to control the current pandemic and to control the spread of the variants. Among the proteins synthesized by the SARS-CoV-2 genome, main protease (M^pro also known as 3CL^pro) is a primary drug target, due to its essential role in maturation of the viral polyproteins. In this study, we provide crystallographic evidence, along with some binding assay data, that three clinically approved anti hepatitis C virus drugs and two other drug-like compounds covalently bind to the M^pro Cys145 catalytic residue in the active site. Also, molecular docking studies can provide additional insight for the design of new antiviral inhibitors for SARS-CoV-2 using these drugs as lead compounds. One might consider derivatives of these lead compounds with higher affinity to the M^pro as potential COVID-19 therapeutics for further testing and possibly clinical trials.

Direct Spectroscopic Evidence for Charge-Assisted Hydrogen-Bond Formation between Ionizable Organic Chemicals and Carbonaceous Materials

The direct evidence for the formation of charge-assisted hydrogen bond (CAHB) between the charged groups of ionizable organic chemicals (IOCs) and carbonaceous materials with similar proton affinity remains elusive. We therefore selected three pharmaceutical contaminants (PCs) as representative IOCs to provide the direct evidence of CAHB formation between IOCs and functionalized carbon nanotubes (CNTs) and its intensity/contribution to PC sorption on CNTs by NMR, FTIR, and DFT analyses. Sorption of PCs on functionalized CNTs resulted in the FTIR characteristic peak that appeared at a higher frequency (3787 cm^-1) and the ¹H NMR characteristic peak that emerged at an extremely low-field region (<18.0 ppm), which can be used as the direct spectroscopic evidence for CAHB formation. Both homonuclear CAHB (HM-CAHB, e.g., [O-H···O]^-) and heteronuclear CAHB (HT-CAHB, e.g., [N⁺-H···O^-]/[O-H···N]⁺) exhibited a much higher sorption energy (|E_ads| ≥ 56.24 kJ/mol) than ordinary hydrogen bond (OHB, |E_ads| ≤ 6.136 kJ/mol), leading to a greater sorption contribution (HM-/HT-CAHB ≥ 42.3%, OHB ≤ 36.5%) and irreversibility (hysteresis index: HM-/HT-CAHB ≥ 1.69, OHB ≤ 0.43) of PCs on CNTs. This work presents the direct evidence for CAHB formation between IOCs and CNTs, which is significant for understanding and predicting the environmental fate and risk of IOCs, thus providing new insights for controlling their pollution using specifically designed carbonaceous materials.

Small Molecules Promote Selective Denaturation and Degradation of Tubulin Heterodimers through a Low-Barrier Hydrogen Bond

Here, we report a novel mechanism to selectively degrade target proteins. 3-(3-Phenoxybenzyl)amino-β-carboline (PAC), a tubulin inhibitor, promotes selective degradation of αβ-tubulin heterodimers. Biochemical studies have revealed that PAC specifically denatures tubulin, making it prone to aggregation that predisposes it to ubiquitinylation and then degradation. The degradation is mediated by a single hydrogen bond formed between the pyridine nitrogen of PAC and βGlu198, which is identified as a low-barrier hydrogen bond (LBHB). In contrast, another two tubulin inhibitors that only form normal hydrogen bonds with βGlu198 exhibit no degradation effect. Thus, the LBHB accounts for the degradation. We then screened for compounds capable of forming an LBHB with βGlu198 and demonstrated that BML284, a Wnt signaling activator, also promotes tubulin heterodimer degradation through the LBHB. Our study provided a unique example of LBHB function and identified a novel approach to obtain tubulin degraders.

The Chemistry of the Passivation Mechanism of Perovskite Films with Ionic Liquids

Passivation of perovskite films by ionic liquids (ILs) improves the performance (efficiency and stability) of perovskite solar cells (PSCs). However, the role of ILs in the passivation of perovskite films is not fully understood. Here, we report the reactions of commonly used ILs with the components of perovskites. The reaction of ILs with perovskite precursors (PbI₂ and methylammonium iodide or formamidinium iodide) in a 1:1:1 molar ratio affords one-dimensional (1D) salts composed of the IL cation interspersed along infinite 1D polymeric [PbI₃]^-_n chains. If the IL is applied in excess, the resulting crystal is composed of six cations surrounding a discrete [Pb₃I₁₂]^6- cluster. All the isolated salts were unambiguously characterized by single-crystal X-ray diffraction analysis, which also reveals extensive hydrogen-bonding interactions.

Solvent polarity dependent ESIPT behavior for the novel flavonoid-based solvatofluorochromic chemosensors

In this work, we explore the excited-state intramolecular proton transfer (ESIPT) mechanisms and relative solvent effects for three novel 3-hydroxylflavone derivatives (i.e., HOF, SHOF, and NSHOF) in acetonitrile, dichloromethane, and toluene solvents. Through calculations, we optimize the structures of HOF, SHOF, and NSHOF. Through the analysis of a series of structural parameters related to hydrogen bonding interactions, it could be found that the hydrogen bonds of the three derivatives are all enhanced in the S₁ state, and more importantly, the excited-state hydrogen bonds of HOF are stronger than those of SHOF and NSHOF. In order to explore the effects of solvent polarity, we analyze the core-valence bifurcation (CVB) index, infrared (IR) vibration spectrum, and the potential energy curves. We find that for HOF, SHOF, and NSHOF, the strength of the excited-state hydrogen bonds increases as the solvent polarity decreases. The solvent polarity dependent ESIPT mechanisms pave the way for further designing novel flavonoid-based solvatofluorochromic probes in future.

Effects of azole rings with different chalcogen atoms on ESIPT behavior for benzochalcogenazolyl-substituted hydroxyfluorenes

ESIPT behavior has attracted a lot of eyes of researchers in recent years because of its unique optical properties. Due to its large Stokes shift and double emission fluorescence, white light can be generated in the fluorophore based on the excited state intramolecular proton transfer (ESIPT) principle. The excited state proton transfer behavior of hydroxylated benzoxazole (BO-OH), benzothiazole (BS-OH) and benzoselenazole (BSe-OH) have been investigated in heptane, chloroform and DMF solvents. By comparing the infrared vibration spectra and the variation of bond parameters from the S₀ to S₁ states, and analyzing the frontier molecular orbitals, the influence of hydrogen bond dynamics, the solvent polarity, charge redistribution and the effects of different proton acceptors on proton transfer were observed. The only structural difference among the three substituted hydroxyfluorenes is the heteroatom in the azole ring (oxygen, sulfur and selenium, respectively). We have scanned the potential energy curve of the ESIPT process, and compared the potential barrier, it is found that the heavier chalcogen atoms are more favorable for proton transfer. At the same time, the potential application of changing heteroatoms in the azole ring by walking down the chalcogenic group in crystal luminescence color regulation is also discussed.

Co-crystals synthesis of 2- and 4-mercaptopyridine with thiourea and its analogue, trithiocyanuric acid

This article describes the synthesis of four new co-crystals obtained via cross-crystallisation of 2- and 4-mercaptopyridine with thiourea and trithiocyanuric acid, and examines their crystal structures. The resulting structures are...

Resonance-Assisted Hydrogen Bond-Revisiting the Original Concept in the Context of Its Criticism in the Literature

In the presented research, we address the original concept of resonance-assisted hydrogen bonding (RAHB) by means of the many-body interaction approach and electron density delocalization analysis. The investigated molecular patterns of RAHBs are open chains consisting of two to six molecules in which the intermolecular hydrogen bond stabilizes the complex. Non-RAHB counterparts are considered to be reference systems. The results show the influence of the neighbour monomers on the unsaturated chains in terms of the many-body interaction energy contribution. Exploring the relation between the energy parameters and the growing number of molecules in the chain, we give an explicit extrapolation of the interaction energy and its components in the infinite chain. Electron delocalization within chain motifs has been analysed from three different points of view: three-body delocalization between C=C-C, two-body hydrogen bond delocalization indices and also between fragments (monomers). A many-body contribution to the interaction energy as well as electron density helps to establish the assistance of resonance in the strength of hydrogen bonds upon the formation of the present molecular chains. The direct relation between interaction energy and delocalization supports the original concept, and refutes some of the criticisms of the RAHB idea.

Hydrogen bonded trimesic acid networks on Cu(111) reveal how basic chemical properties are imprinted in HR-AFM images

High resolution non-contact atomic force microscopy measurements characterize assemblies of trimesic acid molecules on Cu(111) and the link group interactions, providing the first fingerprints utilizing CO-based probes for this widely studied paradigm for hydrogen bond driven molecular self assembly. The enhanced submolecular resolution offered by this technique uniquely reveals key aspects of the competing interactions. Accurate comparison between full-density-based modeled images and experiment allows to identify key structural elements in the assembly in terms of the electron-withdrawing character of the carboxylic groups, interactions of those groups with Cu atoms in the surface, and the valence electron density in the intermolecular region of the hydrogen bonds. This study of trimesic acid assemblies on Cu(111) combining high resolution atomic force microscopy measurements with theory and simulation forges clear connections between fundamental chemical properties of molecules and key features imprinted in force images with submolecular resolution.

The temperature-dependent conformational ensemble of SARS-CoV-2 main protease (Mpro)

The COVID-19 pandemic, instigated by the SARS-CoV-2 coronavirus, continues to plague the globe. The SARS-CoV-2 main protease, or Mpro, is a promising target for development of novel antiviral therapeutics. Previous X-ray crystal structures of Mpro were obtained at cryogenic temperature or room temperature only. Here we report a series of high-resolution crystal structures of unliganded Mpro across multiple temperatures from cryogenic to physiological, and another at high humidity. We interrogate these datasets with parsimonious multiconformer models, multi-copy ensemble models, and isomorphous difference density maps. Our analysis reveals a temperature-dependent conformational landscape for Mpro, including mobile solvent interleaved between the catalytic dyad, mercurial conformational heterogeneity in a key substrate-binding loop, and a far-reaching intramolecular network bridging the active site and dimer interface. Our results may inspire new strategies for antiviral drug development to counter-punch COVID-19 and combat future coronavirus pandemics.

Halogen Bonding in N-Alkyl-3-halogenopyridinium Salts

We performed a structural study of N-alkylated halogenopyridinium cations to examine whether choice of the N-substituent has any considerable effect on the halogen bonding capability of the cations. For that purpose, we prepared a series of N-ethyl-3-halopyridinium iodides and compared them with their N-methyl-3-halopyridinium analogues. Structural analysis revealed that N-ethylated halogenopyridinium cations form slightly shorter C−X⋯I− halogen bonds with iodide anion. We have also attempted synthesis of ditopic symmetric bis-(3-iodopyridinium) dications. Although successful in only one case, the syntheses have afforded two novel ditopic asymmetric monocations with an iodine atom bonded to the pyridine ring and another on the aliphatic N-substituent. Here, the C−I⋯I− halogen bond lengths involving pyridine iodine atom were notably shorter than those involving an aliphatic iodine atom as a halogen bond donor. This trend in halogen bond lengths is in line with the charge distribution on the Hirshfeld surfaces of the cations—the positive charge is predominantly located in the pyridine ring making the pyridine iodine atom σ-hole more positive than the one on the alkyl chan.

Naphthazarin Derivatives in the Light of Intra- and Intermolecular Forces

Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller-Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers' stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm-1-1700 cm-1 and 2300 cm-1-3400 cm-1 in the gas phase and 600 cm-1-1800 cm-1 and 2200 cm-1-3400 cm-1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm-1-1700 cm-1 and 2300 cm-1-3300 cm-1 for the gas phase and one broad absorption region in the solid state between 700 cm-1 and 3100 cm-1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.

A quantum crystallographic approach to short hydrogen bonds

In this work we use high-resolution synchrotron X-ray diffraction for electron density mapping, in conjunction with ab initio modelling, to study short O-H⋯O and O+-H⋯O- hydrogen bonds whose behaviour is known to be tuneable by temperature. The short hydrogen bonds have donor-acceptor distances in the region of 2.45 Å and are formed in substituted urea and organic acid molecular complexes of N,N'-dimethylurea oxalic acid 2 : 1 (1), N,N-dimethylurea 2,4-dinitrobenzoate 1 : 1 (2) and N,N-dimethylurea 3,5-dinitrobenzoic acid 2 : 2 (3). From the combined analyses, these complexes are found to fall within the salt-cocrystal continuum and exhibit short hydrogen bonds that can be characterised as both strong and electrostatic (1, 3) or very strong with a significant covalent contribution (2). An additional charge assisted component is found to be important in distinguishing the relatively uncommon O-H⋯O pseudo-covalent interaction from a typical strong hydrogen bond. The electron density is found to be sensitive to the extent of static proton transfer, presenting it as a useful parameter in the study of the salt-cocrystal continuum. From complementary calculated hydrogen atom potentials, we attribute changes in proton position to the molecular environment. Calculated potentials also show zero barrier to proton migration, forming an 'energy slide' between the donor and acceptor atoms. The better fundamental understanding of the short hydrogen bond in the 'zone of fluctuation' presented in a salt-cocrystal continuum, enabled by studies like this, provide greater insight into their related properties and can have implications in the regulation of pharmaceutical materials.

Entropy-Enthalpy Compensations Fold Proteins in Precise Ways

Exploring the protein-folding problem has been a longstanding challenge in molecular biology and biophysics. Intramolecular hydrogen (H)-bonds play an extremely important role in stabilizing protein structures. To form these intramolecular H-bonds, nascent unfolded polypeptide chains need to escape from hydrogen bonding with surrounding polar water molecules under the solution conditions that require entropy-enthalpy compensations, according to the Gibbs free energy equation and the change in enthalpy. Here, by analyzing the spatial layout of the side-chains of amino acid residues in experimentally determined protein structures, we reveal a protein-folding mechanism based on the entropy-enthalpy compensations that initially driven by laterally hydrophobic collapse among the side-chains of adjacent residues in the sequences of unfolded protein chains. This hydrophobic collapse promotes the formation of the H-bonds within the polypeptide backbone structures through the entropy-enthalpy compensation mechanism, enabling secondary structures and tertiary structures to fold reproducibly following explicit physical folding codes and forces. The temperature dependence of protein folding is thus attributed to the environment dependence of the conformational Gibbs free energy equation. The folding codes and forces in the amino acid sequence that dictate the formation of β-strands and α-helices can be deciphered with great accuracy through evaluation of the hydrophobic interactions among neighboring side-chains of an unfolded polypeptide from a β-strand-like thermodynamic metastable state. The folding of protein quaternary structures is found to be guided by the entropy-enthalpy compensations in between the docking sites of protein subunits according to the Gibbs free energy equation that is verified by bioinformatics analyses of a dozen structures of dimers. Protein folding is therefore guided by multistage entropy-enthalpy compensations of the system of polypeptide chains and water molecules under the solution conditions.

Effects of Intermolecular Hydrogen Bonding and Solvation on Enol-Keto Tautomerism and Photophysics of Azomethine-BODIPY Dyads

Boron-dipyrromethene derivatives (BODIPYs) are a category of molecules with excellent photophysical properties and can be applied to various fields. This work investigates the fluorescent properties of two azomethine-BODIPY dyads in different solvents based on the time-dependent density functional theory (TD-DFT) method. The potential energy curves (PECs) show that the polar protic solvent and the enhanced π-conjugation effect can lower the proton-transfer (PT) barriers, causing the main configuration of NA-BODIPY in methanol to be the keto form, while the main configuration of NA-BODIPY in toluene and SA-BODIPY in methanol and toluene is the enol form. The keto forms of the two compounds possess the twisted intramolecular charge transfer (TICT) decay pathway in the excited state identified by the optimized twisted configurations and the appropriate barriers of the TICT process, whereas the twisted configurations of the enol forms are nonexistent. TICT successfully competes with excited-state proton transfer (ESIPT) of the keto form, which leads to the fluorescence quenching of NA-BODIPY in methanol. This work provides new ideas for the influence of enol-keto tautomerism and the competitiveness of TICT and ESIPT on the photophysical properties of BODIPYs and is expected to provide guidance for the design of new BODIPY functional molecules.

Hydrogen Bonds and n → π* Interactions in the Acetylation of Propranolol Catalyzed by Candida antarctica Lipase B: A QTAIM Study

Enzyme-substrate interactions play a crucial role in enzymatic catalysis. Quantum theory of atoms in molecules (QTAIM) calculations are extremely useful in computational studies of these interactions because they provide very detailed information about the strengths and types of molecular interactions. QTAIM also provides information about the intramolecular changes that occur in the catalytic reaction. Here, we analyze the enzyme-substrate interactions and the topological properties of the electron density in the enantioselective step of the acylation of (R,S)-propranolol, an aminoalcohol with therapeutic applications, catalyzed by Candida antarctica lipase B. Eight reaction paths (four for each enantiomer) are investigated and the energies, atomic charges, hydrogen bonds, and n → π* interactions of propranolol, the catalytic triad (composed of D187, H224, and S105), and the oxyanion hole are analyzed. It is found that D187 acts as an electron density reservoir for H224, and H224 acts as an electron density reservoir for the active site of the protein. It releases electron density when the tetrahedral intermediate is formed from the Michaelis complex and receives it when the enzyme-product complex is formed. Hydrogen bonds can be grouped into noncovalent and covalent hydrogen bonds. The latter are stronger and more important for the reaction than the former. We also found weak n → π* interactions, which are characterized by QTAIM and the natural bond orbital (NBO) analysis.

Dynamics of proton, ion, molecule, and crystal lattice in functional molecular assemblies

Dynamic molecular processes, such as short- or long-range proton (H⁺) and ion (M⁺) motions, and molecular rotations in electrical conducting and magnetic molecular assemblies enable the fabrication of electron-H⁺ (or M⁺) coupling systems, while crystal lattice dynamics and molecular conformation changes in hydrogen-bonded molecular crystals have been utilised in external stimuli responsive reversible gas-induced gate opening and molecular adsorption/desorption behavior. These dynamics of the polar structural units are responsible for the dielectric measurements. The H⁺ dynamics are formed from ferroelectrics and H⁺ conductors, while the dynamic M⁺ motions of Li⁺ and Na⁺ involve ionic conductors and coupling to the conduction electrons. In n-type organic semiconductors, the crystal lattices are modulated by replacing M⁺ cations, with cations such as Li⁺, Na⁺, K⁺, Rb⁺, and Cs⁺. The use of polar rotator or inversion structures such as alkyl amides, m-fluoroanilinium cations, and bowl-shaped trithiasumanene π-cores enables the formation of ferroelectric molecular assemblies. The host-guest molecular systems of ESIPT fluorescent chromic molecules showed interesting molecular sensing properties using various bases, where the dynamic transformation of the crystal lattice and the molecular conformational change were coupled to each other.

Covalence and π-electron delocalization influence on hydrogen bonds in proton transfer process of o-hydroxy aryl Schiff bases: A combined NMR and QTAIM analysis

Within the framework of the density functional theory approach, we studied the relationship between the chemical nature of intramolecular hydrogen bonds (HBs) and nuclear magnetic resonance (NMR) parameters, J-couplings and ¹H-chemical shifts [δ(¹H)], of the atoms involved in such bonds in o-hydroxyaryl Schiff bases during the proton transfer process. For the first time, the shape of the dependence of the degree of covalence in HBs on ¹J(N-H), ¹J(O-H), ^2hJ(O-N), and δ(¹H) during the proton transfer process in o-hydroxyaryl Schiff bases was analyzed. Parameters obtained from Bader's theory of atoms in molecules were used to assess the dependence of covalent character in HBs with both the NMR properties. The influence of π-electronic delocalization on ^2hJ(N-O) under the proton transfer process was investigated. ^2hJ(O-N) in a Mannich base was also studied in order to compare the results with an unsaturated system. In addition, substituent effects on the phenolic ring were investigated. Our results indicate that the covalent character of HBs on both sides of the transition state undergoes a smooth exponential increase as the δ(¹H) moves downfield. The degree of covalence of the N⋯H (O⋯H) bond increases linearly as ¹J(N-H) (¹J(O-H)) becomes more negative, even after reaching the transition state. Non-vanishing values of spin dipolar (SD) and paramagnetic spin orbital terms of ^2hJ(O-N) show that π-electronic delocalization has a non-negligible effect on tautomeric equilibrium and gives evidence of the presence of the resonance assisted HB.Variation of the SD term of ^2hJ(O-N) follows a similar pattern as the change in the para-delocalization aromaticity index of the chelate ring.

Protein Dynamics and Substrate Protonation States Mediate the Catalytic Action of trans-4-Hydroxy-l-Proline Dehydratase

The enzyme trans-4-hydroxy-l-proline (Hyp) dehydratase (HypD) is among the most abundant glycyl radical enzymes (GREs) in the healthy human gut microbiome and is considered a promising antibiotic target for the prominent antibiotic-resistant pathogen Clostridium difficile. Although an enzymatic mechanism has been proposed, the role of the greater HypD protein environment in mediating radical reactivity is not well understood. To fill this gap in understanding, we investigate HypD across multiple time- and length-scales using electronic structure modeling and classical molecular dynamics. We observe that the Hyp substrate protonation state significantly alters both its enzyme-free reactivity and its dynamics within the enzyme active site. Accurate coupled-cluster modeling suggests the deprotonated form of Hyp to be the most reactive protonation state for C5-H_pro-S activation. In the protein environment, hydrophobic interactions modulate the positioning of the Cys434 radical to enhance the reactivity of C5-H_pro-S abstraction. Long-time dynamics reveal that changing Hyp protonation states triggers the switching of a Leu643-gated water tunnel, a functional feature that has not yet been observed for members of the GRE superfamily.