Weak Intermolecular CH···N Hydrogen Bonding: Determination of 13CH-15N Hydrogen-Bond Mediated J Couplings by Solid-State NMR Spectroscopy and First-Principles Calculations

Synthesis and crystal structures of two racemic 2-heteroaryl-3-phenyl-2,3-di-hydro-4H-pyrido[3,2-e][1,3]thia-zin-4-ones

3-Phenyl-2-(thio-phen-3-yl)-2,3-di-hydro-4H-pyrido[3,2-e][1,3]thia-zin-4-one (C17H12N2OS2, 1) and 2-(1H-indol-3-yl)-3-phenyl-2,3-di-hydro-4H-pyrido[3,2-e][1,3]thia-zin-4-one 0.438-hydrate (C21H15N3OS·0.438H2O, 2) crystallize in space groups P21/n and C2/c, respectively. The asymmetric unit in each case is comprised of two parent mol-ecules, albeit of mixed chirality in the case of 1 and of similar chirality in 2 with the enanti-omers occupying the neighboring asymmetric units. Structure 2 also has water mol-ecules (partial occupancies) that form continuous channels along the b -axis direction. The thia-zine rings in both structures exhibit an envelope conformation. Inter-molecular inter-actions in 1 are defined only by C-H⋯O and C-H⋯N hydrogen bonds between crystallographically independent mol-ecules. In 2, hydrogen bonds of the type N-H⋯O between independent mol-ecules and C-H⋯N(π) type, and π-π stacking inter-actions between the pyridine rings of symmetry-related mol-ecules are observed.

Ultrapotent influenza hemagglutinin fusion inhibitors developed through SuFEx-enabled high-throughput medicinal chemistry

Seasonal and pandemic-associated influenza strains cause highly contagious viral respiratory infections that can lead to severe illness and excess mortality. Here, we report on the optimization of our small-molecule inhibitor F0045(S) targeting the influenza hemagglutinin (HA) stem with our Sulfur-Fluoride Exchange (SuFEx) click chemistry-based high-throughput medicinal chemistry (HTMC) strategy. A combination of SuFEx- and amide-based lead molecule diversification and structure-guided design led to identification and validation of ultrapotent influenza fusion inhibitors with subnanomolar EC50 cellular antiviral activity against several influenza A group 1 strains. X-ray structures of six of these compounds with HA indicate that the appended moieties occupy additional pockets on the HA surface and increase the binding interaction, where the accumulation of several polar interactions also contributes to the improved affinity. The compounds here represent the most potent HA small-molecule inhibitors to date. Our divergent HTMC platform is therefore a powerful, rapid, and cost-effective approach to develop bioactive chemical probes and drug-like candidates against viral targets.

Key to High Performance Ion Hybrid Capacitor: Weakly Solvated Zinc Cations

Zinc ion hybrid capacitors suffer from lack of reversibility and dendrite formation. An electrolyte, based on a solution of a zinc salt in acetonitrile and tetramethylene sulfone, allows smooth zinc deposition with high coulombic efficiency in a Zn||stainless steel cell (99.6% for 2880 cycles at 1.0 mA cm-2 , 1.0 mAh cm-2 ). A Zn||Zn cell operates stably for at least 7940 h at 1.0 mA cm-2 with an area capacity of 10 mAh cm-2 , or 648 h at 90% depth of discharge and 1 mA cm-2 , 9.0 mAh cm-2 . Molecular dynamics simulations reveal the reason for the excellent reversibility: The zinc cation is only weakly solvated than in pure tetramethylene sulfone with the closest atoms at 3.3 to 3.8 Å. With this electrolyte, a zinc||activated-carbon hybrid capacitor exhibits an operating voltage of 2.0 to 2.5 V, an energy-density of 135 Wh kg-1 and a power-density of 613 W kg-1 at 0.5 A g-1 . At the very high current-density of 15 A g-1 , 29.3 Wh kg-1 and 14 250 W kg-1 are achieved with 81.2% capacity retention over 9000 cycles.

Exploiting the (-C-H···C-) Interaction to Design Cage-Functionalized Organic Superbases and Hyperbases: A Computational Study

A set of carbon center-based P-ylidesubstituting bases have been exploited computationally with pentacyclo[5.4.0.02,6.03,10.05.9]undecane (PCU) and pentacyclo [6.4.0.02,7.03,11.06,10] dodecane (PCD) scaffolds using the B3LYP-D3/6-311+G(d,p) level of theory. The proton affinities calculated in the gas phase are in the range of superbases and hyperbases. The Atomsin-Molecules and Natural Bond Orbital calculations reveal that the -C-H···C- interaction plays a substantial role in improving the basicity, and tuning the -C-H···C- interaction can enhance the basicity of such systems. The free activation energy for proton exchange for PCD and PCU scaffolds substituted with P-ylide is substantially low. The computed results reveal the strength and nature of such - C-H···C- interactions compared to the -N-H···N- hydrogen bonds. The isodesmic reactions suggest that the superbasicity achieved using these frameworks arises from a combination of several factors, such as the ring strain of the bases in their unprotonated form, steric repulsion, and the intramolecular -C-H···C- interaction.

Strong structuring arising from weak cooperative O-H···π and C-H···O hydrogen bonding in benzene-methanol solution

Weak hydrogen bonds, such as O-H···π and C-H···O, are thought to direct biochemical assembly, molecular recognition, and chemical selectivity but are seldom observed in solution. We have used neutron diffraction combined with H/D isotopic substitution to obtain a detailed spatial and orientational picture of the structure of benzene-methanol mixtures. Our analysis reveals that methanol fully solvates and surrounds each benzene molecule. The expected O-H···π interaction is highly localised and directional, with the methanol hydroxyl bond aligned normal to the aromatic plane and the hydrogen at a distance of 2.30 Å from the ring centroid. Simultaneously, the tendency of methanol to form chain and cyclic motifs in the bulk liquid is manifest in a highly templated solvation structure in the plane of the ring. The methanol molecules surround the benzene so that the O-H bonds are coplanar with the aromatic ring while the oxygens interact with C-H groups through simultaneous bifurcated hydrogen bonds. This demonstrates that weak hydrogen bonding can modulate existing stronger interactions to give rise to highly ordered cooperative structural motifs that persist in the liquid phase.

Discovery of novel diarypyrimidine derivatives bearing six-membered non-aromatic heterocycles as potent HIV-1 NNRTIs with improved anti-resistance and drug-like profiles

Taking our previously reported HIV-1 NNRTIs BH-11c and XJ-10c as lead compounds, series of novel diarypyrimidine derivatives bearing six-membered non-aromatic heterocycles were designed to improve anti-resistance and drug-like profiles. According to the three rounds of in vitro antiviral activity screening, compound 12g was the most active inhibitor against wild-type and five prevalent NNRTI-resistant HIV-1 strains with EC₅₀ values ranging from 0.024 to 0.0010 μM. This is obviously better than the lead compound BH-11c and the approved drug ETR. Detailed structure-activity relationship was investigated to provide valuable guidance for further optimization. The MD simulation study indicated that 12g could form additional interactions with residues around the binding site in HIV-1 RT, which provided reasonable explanations for its improved anti-resistance profile compared to ETR. Furthermore, 12g showed significant improvement in water solubility and other drug-like properties compared to ETR. The CYP enzymatic inhibitory assay indicated that 12g was unlikely to induce CYP-mediated drug-drug interactions. 12g pharmacokinetics parameters were investigated and it displayed a long half-life of 6.59 h in vivo. The properties of compound 12g make it a promising lead compound for the development of new generation of antiretroviral drugs.

Discovery of diarylpyrimidine derivatives bearing piperazine sulfonyl as potent HIV-1 nonnucleoside reverse transcriptase inhibitors

HIV-1 reverse transcriptase is one of the most attractive targets for the treatment of AIDS. However, the rapid emergence of drug-resistant strains and unsatisfactory drug-like properties seriously limit the clinical application of HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs). Here we show that a series of piperazine sulfonyl-bearing diarylpyrimidine-based NNRTIs were designed to improve the potency against wild-type and NNRTI-resistant strains by enhancing backbone-binding interactions. Among them, compound 18b1 demonstrates single-digit nanomolar potency against the wild-type and five mutant HIV-1 strains, which is significantly better than the approved drug etravirine. The co-crystal structure analysis and molecular dynamics simulation studies were conducted to explain the broad-spectrum inhibitory activity of 18b1 against reverse transcriptase variants. Besides, compound 18b1 demonstrates improved water solubility, cytochrome P450 liability, and other pharmacokinetic properties compared to the currently approved diarylpyrimidine (DAPY) NNRTIs. Therefore, we consider compound 18b1 a potential lead compound worthy of further study.

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the "gold standard" coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller-Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

To what extent do bond length and angle govern the 13C and 1H NMR response to weak CH⋯O hydrogen bonds? A case study of caffeine and theophylline cocrystals

Weak hydrogen bonds are important structure-directing elements in supramolecular chemistry and biochemistry. We consider here weak CH⋯O hydrogen bonds in a series of cocrystals of theophylline and caffeine and assess to what extent the CH⋯O distance and angle govern the observed ¹³C and ¹H isotropic chemical shifts. Gauge-including projector-augmented wave density functional theory (GIPAW DFT) calculations consistently predict a decrease in the ¹³C and ¹H magnetic shielding constants upon hydrogen bond formation on the order of 2-5 ppm (¹³C) and 1-2 ppm (¹H). These trends are reproduced using the machine-learning approach implemented in ShiftML. Experimental ¹³C and ¹H chemical shifts obtained for powdered samples using one-dimensional NMR spectroscopy as well as heteronuclear correlation (HETCOR) spectroscopy correlate well with the GIPAW DFT results. However, the experimental ¹³C NMR response only correlates moderately well with the hydrogen bond length and angle, while the experimental ¹H chemical shifts only show very weak correlations to these local structural elements. DFT computations on isolated imidazole-formaldehyde models show that the ¹³C and ¹H chemical shifts generally decrease with the C⋯O distance but show no clear dependence on the CH⋯O angle. These results demonstrate that the ¹³C and ¹H response to weak CH⋯O hydrogen bonding is influenced significantly by additional weak contacts within cocrystal heterodimeric units.

Solid-State Photo-NMR Study on Light-Induced Nitrosyl Linkage Isomers Uncovers Their Structural, Electronic, and Diamagnetic Nature

A light-induced linkage NO isomer (MS1) in trans-[Ru(¹⁵NO)(py)₄¹⁹F](ClO₄)₂ is detected and measured for the first time by solid-state MAS NMR. Chemical shift tensors of ¹⁵N and ¹⁹F, along with ⁿJ(¹⁵N-¹⁹F) spin-spin couplings and T₁ relaxation times of MS1, are compared with the ground state (GS) at temperatures T < 250 K. Isotropic chemical shifts (¹⁵N and ¹⁹F) are well resolved for two crystallographically independent cations (A and B) [Ru(¹⁵NO)(py)₄¹⁹F]²⁺, allowing to define separately both populations of MS1 isomers and thermal decay rates for two structural sites. The relaxation times T₁ of ¹⁹F in the case of GS (30/38.6 s for sites A/B) and MS1 (11.6/11.8 s for sites A/B) indicate that both isomers are diamagnetic, which is the first experimental evidence of diamagnetic properties of MS1 in ruthenium nitrosyl. After light irradiation (λ = 420 nm), the NO ligand rotates by nearly 180° from F-Ru-N-O to F-Ru-O-N, whereby the isotropic chemical shifts of δ_iso(¹⁵N) increase and those of δ_iso(¹⁹F) decrease. The ⁿJ(¹⁵N-¹⁹F) couplings increase from ²J(¹⁵N-Ru-¹⁹F)_GS = 71 Hz to ³J(¹⁵N-O-Ru-¹⁹F)_MS1 = 105 Hz. These results are interpreted on the basis of DFT-CASTEP calculations including Bader-, Mulliken-, and Hirshfeld-charge density distributions of both states.

Alkyl substituent-dependent systematic change in cold crystallization of azo molecules

The thermal behavior of alkyl-derivatized 1-(2,4-dimethylphenylazo)-4-naphthol and 1-(2,4-dimethylphenylazo)-2-naphthol (2,4-DM-4-Cn and 2,4-DM-2-Cn, respectively) was investigated. The change in the position of the alkyl substituent led to a variation in the thermal behavior, including the cold crystallization, which is a heat-storing phenomenon. In addition, a comprehensive study of the alkyl chain length revealed that 2,4-DM-4-Cn had better crystallinity and exhibited cold crystallization with short alkyl chains. The π–π, C–H⋯N, and C–H⋯π interactions stabilized the crystal structure of 2,4-DM-4-Cn. On the other hand, the polymorphism of 2,4-DM-2-Cn inhibited the formation of a uniform crystalline phase during cooling, which led to poor crystallinity. The only difference between the compounds, the position of the substituent, resulted in a clear variation in the cold crystallization and heat storage properties.

The change in alkyl substitution position and alkyl chain length of azo molecules led to the variation in cold crystallization behavior.

Mapping of N-C Bond Formation from a Series of Crystalline Peri-Substituted Naphthalenes by Charge Density and Solid-State NMR Methodologies

A combination of charge density studies and solid state nuclear magnetic resonance (NMR) 1 JNC coupling measurements supported by periodic density functional theory (DFT) calculations is used to characterise the transition from an n-π* interaction to bond formation between a nucleophilic nitrogen atom and an electrophilic sp2 carbon atom in a series of crystalline peri-substituted naphthalenes. As the N⋅⋅⋅C distance reduces there is a sharp decrease in the Laplacian derived from increasing charge density between the two groups at ca. N⋅⋅⋅C = 1.8 Å, with the periodic DFT calculations predicting, and heteronuclear spin-echo NMR measurements confirming, the 1 JNC couplings of ≈3-6 Hz for long C-N bonds (1.60-1.65 Å), and 1 JNC couplings of <1 Hz for N⋅⋅⋅C >2.1 Å.

Pure Isotropic Proton Solid State NMR

Resolution in proton solid state magic angle sample spinning (MAS) NMR is limited by the intrinsically imperfect nature of coherent averaging induced by either MAS or multiple pulse sequence methods. Here, we suggest that instead of optimizing and perfecting a coherent averaging scheme, we could approach the problem by parametrically mapping the error terms due to imperfect averaging in a k-space representation, in such a way that they can be removed in a multidimensional correlation leaving only the desired pure isotropic signal. We illustrate the approach here by determining pure isotropic ¹H spectra from a series of MAS spectra acquired at different spinning rates. For six different organic solids, the approach is shown to produce pure isotropic ¹H spectra that are significantly narrower than the MAS spectrum acquired at the fastest possible rate, with linewidths down to as little as 48 Hz. On average, we observe a 7-fold increase in resolution, and up to a factor of 20, as compared with spectra acquired at 100 kHz MAS. The approach is directly applicable to a range of solids, and we anticipate that the same underlying principle for removing errors introduced here can be applied to other problems in NMR spectroscopy.

Complex formation of titanocene dichloride anticancer and Al12N12 nano-cluster: A quantum chemical investigation of solvent, temperature and pressure effects

This study investigated the interaction between Al12N12 nano-cluster and titanocene dichloride anticancer drug complex using B3P86 functional in gas and solution phases. Non-covalent interaction (NCI) analysis of this complex was employed for illustration of the Cl⋯Al weak non-covalent interaction. The self-consistent reaction field theory (SCRF) based on the Polarizable Continuum Model (PCM) was applied for testing the solvent effects. The solvent effect on the interaction energy, dipole moment, frontier orbital energy, and global reactivity parameters was examined as well. The changes in the dipole moment, polarizability and electronic spatial extent (ESE) with solvent polarity were analyzed by applying different solvent polarity parameters based on Lippert-Mataga, Bakhshiev and Bilot-Kawski models. In addition, temperature and pressure effects on the thermodynamic parameters of this interaction were illustrated.

Revealing weak histidine 15N homonuclear scalar couplings using Solid-State Magic-Angle-Spinning NMR spectroscopy

The tautomeric structure and chemistry of the histidine imidazole ring play active roles in many structurally and functionally important proteins and polypeptides. While in NMR spectroscopy histidine chemical shifts (e.g. 15N, 13C, and 1H) have been commonly used to characterize the tautomeric structure, hydrogen bonding, and torsion angles, homonuclear 15N scalar couplings in histidine have rarely been reported. Here, we propose double spin-echo sequences to compare the observed signals with and without a 90° pulse between the two spin-echo periods, such that their signal ratio as a function of the echo time solely depends on homonuclear scalar couplings, allowing for measuring weak homonuclear scalar couplings without influence from transverse dephasing effects, thus capable of revealing hydrogen-bond mediated 15N-15N J-couplings that can provide direct and definitive evidence for the formation of N…H…N hydrogen-bonding associated with the imidazole ring. We used two 13C,15N labeled histidine samples recrystallized from solutions at pH 6.3 and pH 11.0 to demonstrate the feasibility of this method and reveal the existence of a weak two-bond scalar coupling between the Nδ1 and Nε2 sites in the histidine imidazole ring in three tautomeric states and the presence of a hydrogen-bond mediated scalar coupling between the Nδ1 site in the imidazole ring and the backbone Nα site in the histidine neutral τ and π states. Our results demonstrate that weak 15N homonuclear scalar couplings can be measured even when their values are less than their corresponding intrinsic natural linewidths, thus providing direct and definitive evidence for the formation of N…H…N hydrogen bonding that is associated with the histidine imidazole ring.