Strong N−H···π Hydrogen Bonding in Amide−Benzene Interactions

Boosting the Modeling of Infrared and Raman Spectra of Bulk Phase Chromophores with Machine Learning

In the field of vibrational spectroscopy simulation, hybrid approximations to Kohn-Sham density-functional theory (KS-DFT) are often considered computationally prohibitive due to the significant effort required to evaluate the exchange-correlation potential in planewave codes. In this Letter, we show that by leveraging the porting of KS-DFT on GPU and incorporating machine-learning techniques, simulating IR and Raman spectra of real-life chromophores in bulk aqueous solution becomes a routine application at this level of theory.

Recent Developments and Challenges in Solid-Contact Ion-Selective Electrodes

Solid-contact ion-selective electrodes (SC-ISEs) have the advantages of easy miniaturization, even chip integration, easy carrying, strong stability, and more favorable detection in complex environments. They have been widely used in conjunction with portable, wearable, and intelligent detection devices, as well as in on-site analysis and timely monitoring in the fields of environment, industry, and medicine. This article provides a comprehensive review of the composition of sensors based on redox capacitive and double-layer capacitive SC-ISEs, as well as the ion-electron transduction mechanisms in the solid-contact (SC) layer, particularly focusing on strategies proposed in the past three years (since 2021) for optimizing the performance of SC-ISEs. These strategies include the construction of ion-selective membranes, SC layer, and conductive substrates. Finally, the future research direction and possibilities in this field are discussed and prospected.

A Micromechanical Study of Interactions of Cyanate Ester Monomer with Graphene or Boron Nitride Monolayer

Polymer composites, hailed for their ultra-strength and lightweight attributes, stand out as promising materials for the upcoming era of space vehicles. The selection of the polymer matrix plays a pivotal role in material design, given its significant impact on bulk-level properties through the reinforcement/polymer interface. To aid in the systematic design of such composite systems, molecular-level calculations are employed to establish the relationship between interfacial characteristics and mechanical response, specifically stiffness. This study focuses on the interaction of fluorinated and non-fluorinated cyanate ester monomers with graphene or a BN monolayer, representing non-polymerized ester composites. Utilizing micromechanics and the density functional theory method to analyze interaction energy, charge density, and stiffness, our findings reveal that the fluorinated cyanate-ester monomer demonstrates lower interaction energy, reduced pull-apart force, and a higher separation point compared to the non-fluorinated counterpart. This behavior is attributed to the steric hindrance caused by fluorine atoms. Furthermore, the BN monolayer exhibits enhanced transverse stiffness due to increased interfacial strength, stemming from the polar nature of B-N bonds on the surface, as opposed to the C-C bonds of graphene. These molecular-level results are intended to inform the design of next-generation composites incorporating cyanate esters, specifically for structural applications.

E- or Z-Isomers Arising from the Geometries of Ligands in the Mercury Complex of 2-(Anthracen-9-ylmethylene)-N-phenylhydrazine Carbothioamide

An anionic mercury(II) complex of 2-(anthracen-9-ylmethylene)-N-phenylhydrazine carbothioamide (HATU) and two isomers of a neutral mercury(II) complex of the anion of the same ligand (ATU) were reported. The anionic complex [Hg(HATU)2Cl2]·CH2Cl2 had a monodentate HATU ligand (a neutral form of the ligand) and chloride ligands. The two conformational isomers were of the neutral mercury(II) complex Hg(ATU)2·2DMF. The two isomers were from the E or Z geometry of the ligands across the conjugated C=N-N=C-N scaffold of the coordinated ligand. The two isomers of the complex were independently prepared and characterized. The spectroscopic properties of the isomers in solution were studied by 1H NMR as well as fluorescence spectroscopy. Facile conversion of the E-isomer to the Z-isomer in solution was observed. Density functional theory (DFT) calculations revealed that the Z-isomer of the complex was stable compared to the E-isomer by an energy of 14.35 kJ/mol; whereas, E isomer of the ligand was more stable than Z isomer by 8.37 KJ/mol. The activation barrier for the conversion of the E-isomer to the Z-isomer of the ligand was 167.37 kJ/mol. The role of the mercury ion in the conversion of the E-form to the Z-form was discussed. The mercury complex [Hg(HATU)2Cl2]·CH2Cl2 had the E-form of the ligand. Distinct photophysical features of these mercury complexes were presented.

Modulation of n → π* Interaction in the Complexes of p-Substituted Pyridines with Aldehydes: A Theoretical Study

n → π* interaction is analogous to the hydrogen bond in terms of the delocalization of the electron density between the two orbitals. Studies on the intermolecular complexes stabilized by the n → π* interaction are scarce in the literature. Herein, we have studied intermolecular N···C═O n → π* interactions in the complexes of p-substituted pyridines (p-R-Py) with formaldehyde (HCHO), formyl chloride (HCOCl), and acetaldehyde (CH₃CHO) using quantum chemistry calculations. We have shown that the strength of the n → π* interaction can be modulated by varying the electronic substituents at the donor and acceptor sites in the complexes. Variation of the substituents at the para position of the pyridine ring from the electron-withdrawing groups (EWGs) to the electron-donating groups (EDGs) results in a systematic increase in the strength of the n → π* interaction. The strength of this interaction is also modulated by tuning the electron density toward the carbonyl bond by substituting the hydrogen atom of HCHO with the methyl and chloro groups. The modulation of this interaction due to the electronic substitutions at the n → π* donor and acceptor sites in the complexes is monitored by probing the relevant geometrical parameters, binding energies, C═O frequency redshift, NBO energies, and electron density for this interaction derived from QTAIM and NCI index analyses. Energy decomposition analysis reveals that the electrostatic interaction dominates the binding energies of these complexes, while the charge transfer interaction, which is representative of the n → π* interaction, also has a significant contribution to these.

Effects of phenyl acids on different degradation phases during thermophilic anaerobic digestion

Aromatic compounds like phenyl acids (PA) can accumulate during anaerobic digestion (AD) of organic wastes due to an increased entry of lignocellulose, secondary plant metabolites or proteins, and thermodynamic challenges in degrading the benzene ring. The effects of aromatic compounds can be various - from being highly toxic to be stimulating for methanogenesis - depending on many parameters like inoculum or molecular characteristics of the aromatic compound. To contribute to a better understanding of the consequences of PA exposure during AD, the aim was to evaluate the effects of 10 mM PA on microbial communities degrading different, degradation phase-specific substrates in thermophilic batch reactors within 28  days: Microcrystalline cellulose (MCC, promoting hydrolytic to methanogenic microorganisms), butyrate or propionate (promoting syntrophic volatile fatty acid (VFA) oxidisers to methanogens), or acetate (promoting syntrophic acetate oxidisers to methanogens). Methane production, VFA concentrations and pH were evaluated, and microbial communities and extracellular polymeric substances (EPS) were assessed. The toxicity of PA depended on the type of substrate which in turn determined the (i) microbial diversity and composition and (ii) EPS quantity and quality. Compared with the respective controls, methane production in MCC reactors was less impaired by PA than in butyrate, propionate and acetate reactors which showed reductions in methane production of up to 93%. In contrast to the controls, acetate concentrations were high in all PA reactors at the end of incubation thus acetate was a bottle-neck intermediate in those reactors. Considerable differences in EPS quantity and quality could be found among substrates but not among PA variants of each substrate. Methanosarcina spp. was the dominant methanogen in VFA reactors without PA exposure and was inhibited when PA were present. VFA oxidisers and Methanothermobacter spp. were abundant in VFA assays with PA exposure as well as in all MCC reactors. As MCC assays showed higher methane yields, a higher microbial diversity and a higher EPS quantity and quality than VFA reactors when exposed to PA, we conclude that EPS in MCC reactors might have been beneficial for absorbing/neutralising phenyl acids and keeping (more susceptible) microorganisms shielded in granules or biofilms.

All-solid-state paper-based potentiometric combined sensor modified with reduced graphene oxide (rGO) and molecularly imprinted polymer for monitoring losartan drug in pharmaceuticals and biological samples

A cost-effective, highly selective and sensitive paper-based potentiometric combined sensor for losartan potassium drug (LOS) is fabricated, characterized and used for the drug monitoring. The sensor consists of 2 strips of filter paper (20 × 5 mm each) as platform, each imprinted with 4 mm diameter circular spot of carbon. One carbon spot is covered by a reduced graphene oxide (rGO) for use as a substrate for the recognition sensor and the other without rGO is used for the reference electrode. LOS molecularly imprinted drug polymer (MIP) is applied onto the graphene oxide containing strip to act as a drug recognition sensing material and a solid-state polyvinyl butyral (PVB) is applied onto the second carbon spot to act as a reference electrode. Performance characteristics of the combined sensor are examined with chronopotentiometry (CP) and electrochemical impedance spectroscopy (EIS). Increase effect of rGO on the interfacial double-layer capacitance of the sensing membrane and consequently on the potential stability is confirmed. The developed combined sensor (strip cell) displays a Nernstian slope of -58.2 ± 0.3 mV/decade (R² = 0.9994) over the linear range 8.5 × 10^-7 - 6.9 × 10^-2 M with a detection limit of 2.7 ± 0.3 × 10^-7 M. The sensor shows remarkable selectivity toward various related compounds especially those commonly used by the COVID-19 patients such as paracetamol, ascorbic acid and dextromethorphan. The assay method is validated and proved to be satisfactory for direct potentiometric determination of LOS-K in some pharmaceutical formulations and in spiked human urine samples. An average recovery of 96.3 ± 0.3-98.7 ± 0.6% of the nominal or spiked concentration and a mean relative standard deviation of ±0.6% are obtained. The use of an indicating and a reference electrodes combined into a single flexible disposable paper platform enables applications to a minimum sample volume due to the close proximity of the responsive membrane and the liquid junction. The efficiency of the proposed sensor in complex urine matrix suggests its application in hospitals for rapid diagnosis of overdose patients and for quality control/quality assurance tests in pharmaceutical industry.

Microhydration of Phenyl Formate: Gas-Phase Laser Spectroscopy, Microwave Spectroscopy, and Quantum Chemistry Calculations

Herein, we have investigated the structure of phenyl formate⋅⋅⋅water (PhOF⋅⋅⋅H₂ O) dimer and various non-covalent interactions present there using gas-phase laser spectroscopy and microwave spectroscopy combined with quantum chemistry calculations. Two conformers of PhOF⋅⋅⋅H₂ O (C1 and T1), built on the two cis/trans conformers of the bare molecule, have been observed in the experiment. In cis-PhOF, there is an n_CO → π A r * ${{{\rm \pi }}_{{\rm A}{\rm r}}^{{\rm {^\ast}}}}$ interaction between the lone-pair orbital of the carbonyl oxygen atom and the π* orbital of the phenyl ring, which persists in the monohydrated C1 conformer of PhOF⋅⋅⋅H₂ O according to the NBO and NCI analyses. On the other hand, this interaction is absent in the trans-PhOF conformer as the C=O group is away from the phenyl ring. The C1 conformer is primarily stabilized by an interplay between O-H⋅⋅⋅O=C hydrogen bond and O-H⋅⋅⋅π interactions, while the stability of the T1 conformer is primarily governed by the O-H⋅⋅⋅O=C hydrogen bond. The most important finding of the present work is that the conformational preference of the PhOF monomer is retained in its monohydrated complex.

The retaining β-Kdo glycosyltransferase WbbB uses a double-displacement mechanism with an intermediate adduct rearrangement step

WbbB, a lipopolysaccharide O-antigen synthesis enzyme from Raoultella terrigena, contains an N-terminal glycosyltransferase domain with a highly modified architecture that adds a terminal β-Kdo (3-deoxy-D-manno-oct-2-ulosonic acid) residue to the O-antigen saccharide, with retention of stereochemistry. We show, using mass spectrometry, that WbbB forms a covalent adduct between the catalytic nucleophile, Asp232, and Kdo. We also determine X-ray structures for the CMP-β-Kdo donor complex, for Kdo-adducts with D232N and D232C WbbB variants, for a synthetic disaccharide acceptor complex, and for a ternary complex with both a Kdo-adduct and the acceptor. Together, these structures show that the enzyme-linked Asp232-Kdo adduct rotates to reposition the Kdo into a second sub-site, which then transfers Kdo to the acceptor. Retaining glycosyltransferases were thought to use only the front-side S_Ni substitution mechanism; here we show that retaining glycosyltransferases can also potentially use double-displacement mechanisms, but incorporating an additional catalytic subsite requires rearrangement of the protein's architecture.

Improving Antimicrobial Activity and Physico-Chemical Properties by Isosteric Replacement of 2-Aminothiazole with 2-Aminooxazole

Antimicrobial drug resistance is currently one of the most critical health issues. Pathogens resistant to last-resort antibiotics are increasing, and very few effective antibacterial agents have been introduced in recent years. The promising drug candidates are often discontinued in the primary stages of the drug discovery pipeline due to their unspecific reactivity (PAINS), toxicity, insufficient stability, or low water solubility. In this work, we investigated a series of substituted N-oxazolyl- and N-thiazolylcarboxamides of various pyridinecarboxylic acids. Final compounds were tested against several microbial species. In general, oxazole-containing compounds showed high activity against mycobacteria, especially Mycobacterium tuberculosis (best MIC_H37Ra = 3.13 µg/mL), including the multidrug-resistant strains. Promising activities against various bacterial and fungal strains were also observed. None of the compounds was significantly cytotoxic against the HepG2 cell line. Experimental measurement of lipophilicity parameter log k’_w and water solubility (log S) confirmed significantly (typically two orders in logarithmic scale) increased hydrophilicity/water solubility of oxazole derivatives in comparison with their thiazole isosteres. Mycobacterial β-ketoacyl-acyl carrier protein synthase III (FabH) was suggested as a probable target by molecular docking and molecular dynamics simulations.

Dibenzo-18-Crown-6-based Carbon Paste Sensors for the Nanomolar Potentiometric Determination of Daclatasvir Dihydrochloride: An Anti-HCV Drug and a Potential Candidate for Treatment of SARS-CoV-2

Daclatasvir dihydrochloride (DAC) is an anti-hepatitis C virus (HCV) drug that has recently proven to be a promising candidate for the treatment of SARS-CoV-2. Still, there is a lack of sensitive potentiometric methods for its determination. In this work, carbon paste sensors based on dibenzo-18-crown-6 (DB18C6) were fabricated and optimized for the sensitive and selective potentiometric determination of DAC in Daclavirocyrl® tablets, serum, and urine samples. The best performance was obtained by two sensors referred to as sensor I and sensor II. Both sensors exhibited a wide linear response range of 5×10⁻⁹ − 1×10⁻³ mol/L, and Nernstian slopes of 29.8 ± 1.18 and 29.5 ± 1.00 mV/decade, with limits of detection, 4.8×10⁻⁹ and 3.2×10⁻⁹ mol/L, for the sensors I and II, respectively. Sensors I and II displayed fast response times of 5–8 and 5–6 s, respectively, with great reversibility and no memory effect. Moreover, the sensors exhibited a lifetime of 16 days. For the study of sensors morphology and elucidation of the interaction mechanism, the scanning electron microscope (SEM), Fourier-transform infrared spectroscopy (FTIR), and nuclear magnetic resonance (¹H NMR) techniques were performed. A selectivity study was performed, and the proposed sensors exhibited good discrimination between DAC and potentially coexisting interferents with sensor II displaying better selectivity. Finally, sensor II was successfully applied for the determination of DAC in the above-mentioned samples, with recovery values ranging from 99.25 to 101.42%, and relative standard deviation (RSD) values ranging from 0.79 to 1.53% which reflected the high accuracy and precision.

Heat adaptation of phage T7 under an extended genetic code

While bacteriophages have previously been used as a model system to understand thermal adaptation, most adapted genomes observed to date contain very few modifications and cover a limited temperature range. Here, we set out to investigate genome adaptation to thermal stress by adapting six populations of T7 bacteriophage virions to increasingly stringent heat challenges. Further, we provided three of the phage populations’ access to a new genetic code in which Amber codons could be read as selenocysteine, potentially allowing the formation of more stable selenide-containing bonds. Phage virions responded to the thermal challenges with a greater than 10°C increase in heat tolerance and fixed highly reproducible patterns of non-synonymous substitutions and genome deletions. Most fixed mutations mapped to either the tail complex or to the three internal virion proteins that form a pore across the E. coli cell membrane during DNA injection. However, few global changes in Amber codon usage were observed, with only one natural Amber codon being lost. These results reinforce a model in which adaptation to thermal stress proceeds via the cumulative fixation of a small set of highly adaptive substitutions and that adaptation to new genetic codes proceeds only slowly, even with the possibility of potential phenotypic advantages.

In silico investigation of potential inhibitors to main protease and spike protein of SARS-CoV-2 in propolis

Docking analysis of propolis's natural compound was successfully performed against SARS-CoV-2 main protease (Mpro) and spike protein subunit 2 (S2). Initially, the propolis's protein was screened using chromatography analysis and successfully identified 22 compounds in the propolis. Four compounds were further investigated, i.e., neoblavaisoflavone, methylophiopogonone A, 3'-Methoxydaidzin, and genistin. The binding affinity of 3'-Methoxydaidzin was -7.7 kcal/mol, which is similar to nelfinavir (control), while the others were -7.6 kcal/mol. However, we found the key residue of Glu A:166 in the methylophiopogonone A and genistin, even though the predicted binding energy slightly higher than nelfinavir. In contrast, the predicted binding affinity of neoblavaisoflavone, methylophiopogonone A, 3'-Methoxydaidzin, and genistin against S2 were -8.1, -8.2, -8.3, and -8.3 kcal/mol, respectively, which is far below of the control (pravastatin, -7.3 kcal/mol). Instead of conventional hydrogen bonding, the π bonding influenced the binding affinity against S2. The results reveal that this is the first report about methylophiopogonone A, 3'-Methoxydaidzin, and genistin as candidates for anti-viral agents. Those compounds can then be further explored and used as a parent backbone molecule to develop a new supplementation for preventing SARS-CoV-2 infections during COVID-19 outbreaks.

A Photochemical Intramolecular C-N Coupling Toward the Synthesis of Benzimidazole-Fused Phenanthridines

Herein, we report a direct photochemical dehydrogenative C-N coupling of unactivated C(sp²)-H and N(sp²)-H bonds. The catalysts or additive-free transformation of 2-([1,1'-biphenyl]-2-yl)-1H-benzo[d]imidazole to benzo[4,5]imidazo[1,2-f]phenanthridine was achieved at ∼350 nm of irradiation via ε-hydrogen abstraction. DFT calculations helped to understand that the N-H···π interaction was essential for the reaction to proceed at a lower energy than expected.

Charge-Free, Stabilizing Amide-π Interactions Can Be Used to Control Collagen Triple-Helix Self-Assembly

There is a noted lack of understood, controllable interactions for directing the organization of collagen triple helices. While the field has had success using charge-pair interactions and cation-π interactions in helix design, these alone are not adequate for achieving the degree of specificity desirable for these supramolecular structures. Furthermore, because of the reliance on electrostatic interactions, designed heterotrimeric systems have been heavily charged, a property undesirable in some applications. Amide-π interactions are a comparatively understudied class of charge-free interactions, which could potentially be harnessed for triple-helix design. Herein, we propose, validate, and utilize pairwise amino acid amide-π interactions in collagen triple-helix design. Glutamine-phenylalanine pairs, when arranged in an axial geometry, are found to exhibit a moderately stabilizing effect, while in the lateral geometry, this pair is destabilizing. Together this allows glutamine-phenylalanine pairs to effectively set the register of triple helices. In contrast, interactions between asparagine and phenylalanine appear to have little effect on triple-helical stability. After deconvoluting the contributions of these amino acids to triple-helix stability, we demonstrate these new glutamine-phenylalanine interactions in the successful design of a heterotrimeric triple helix. The results of all of these analyses are used to update our collagen triple-helix thermal stability prediction algorithm, Scoring function for Collagen Emulating Peptides' Temperature of Transition (SCEPTTr).

O-H stretching frequency red shifts do not correlate with the dissociation energies in the dimethylether and dimethylsulfide complexes of phenol derivatives

In this perspective, we present a comprehensive report on the spectroscopic and computational investigations of the hydrogen bonded (H-bonded) complexes of Me2O and Me2S with seven para-substituted H-bond donor phenols. The salient finding was that although the dissociation energies, D0, of the Me2O complexes were consistently higher than those of the analogous Me2S complexes, the red-shifts in phenolic O-H frequencies, Δν(O-H), showed the exactly opposite trend. This is in contravention of the general perception that the red shift in the X-H stretching frequency in the X-HY hydrogen bonded complexes is a reliable indicator of H-bond strength (D0), a concept popularly known as the Badger-Bauer rule. This is also in contrast to the trend reported for the H-bonded complexes of H2S/H2O with several para substituted phenols of different pKa values wherein the oxygen centered hydrogen bonded (OCHB) complexes consistently showed higher Δν(O-H) and D0 compared to those of the analogous sulfur centered hydrogen bonded (SCHB) complexes. Our effort was to understand these intriguing observations based on the spectroscopic investigations of 1 : 1 complexes in combination with a variety of high level quantum chemical calculations. Ab initio calculations at the MP2 level and the DFT calculations using various dispersion corrected density functionals (including DFT-D3) were performed on counterpoise corrected surfaces to compute the dissociation energy, D0, of the H-bonded complexes. The importance of anharmonic frequency computations is underscored as they were able to correctly reproduce the observed trend in the relative OH frequency shifts unlike the harmonic frequency computations. We have attempted to find a unified correlation that would globally fit the observed red shifts in the O-H frequency with the H-bonding strength for the four bases, namely, H2S, H2O, Me2O, and Me2S, in this set of H-bond donors. It was found that the proton affinity normalized Δν(O-H) values scale very well with the H-bond strength.

A docosyl-terminated polyamine amphiphile-bonded stationary phase for multimodal separations in liquid chromatography

A convenient synthetic approach to a linear alkyl-polyamine amphiphilic chromatographic selector was proposed. Successive immobilization of the amphiphile onto silica gel afforded a multimodal stationary phase for high-performance liquid chromatography (HPLC). The as-prepared silica material was studied comparatively with a conventional octadecyl (C18) and an amide-embedded C18 stationary phase. The new uniform docosyl-triamine tandem was featured by an enhanced shape selectivity towards geometric isomers, and a low silanol activity towards alkaline solutes. The presence of multiple amino groups rendered the new adsorbent operable in different modes, such as hydrophilic interaction and ion-exchange modes. The satisfactory performance of the said stationary phase in separating different classes of analytes, including polycyclic aromatic hydrocarbons, flavonoids, tricyclic antidepressants, calcium channel blockers, aromatic acids, inorganic anions, nucleosides and estrogens, revealed its great potential and high adaptability for multipurpose LC utility.

5-[(Pyren-9-ylmethyl)amino]isophthalic acid with nitrogen containing heterocycles: stacking, N–H⋯π interactions and photoluminescence

Synthons guided the types of N–H⋯π interactions and stacking to cause quenching of emissions.

Sequence specific hydrogen bond of DNA with denaturants affects its stability: Spectroscopic and simulation studies

BACKGROUND

Several different small molecules have been used to target the DNA helix in order to treat the diseases caused by its mutation. Guanidinium(Gdm⁺) and urea based drugs have been used for the diseases related to central nervous system, also as the anti-inflammatory and chemotherapeutic agent. However, the role of Gdm⁺ and urea in the stabilization/destabilization of DNA is not well understood.

METHODS

Spectroscopic techniques along with molecular dynamics (MD) simulation have been performed on different sequences of DNA in the presence of guanidinium chloride (GdmCl) and urea to decode the binding of denaturants with DNA and the role of hydrogen bond with the different regions of DNA in its stability/destability.

RESULTS AND CONCLUSION

Our study reveals that, Gdm⁺ of GdmCl and urea both intrudes into the groove region of DNA along with the interaction with its phosphate backbone. However, interaction of Gdm⁺ and urea with the nucleobases in the groove region is different. Gdm⁺ forms the intra-strand hydrogen bond with the central region of the both sequences of DNA whereas inter-strand hydrogen bond along with water assisted hydrogen bond takes place in the case of urea. The intra-strand hydrogen bond formation capability of Gdm⁺ with the nucleobases in the minor groove of DNA decreases its groove width which probably causes the stabilization of B-DNA in GdmCl. In contrast, the propensity of the formation of inter-strand hydrogen bond of urea with the nucleobases in the groove region of DNA without affecting the groove width destabilizes B-DNA as compared to GdmCl. This study depicts that the opposite effect of GdmCl and urea on the stability is a general property of B-DNA. However, the extent of stabilization/destabilization of DNA in Gdm⁺ and urea depend on its sequence probably due to the difference in the intra/inter-strand hydrogen bonding with different bases present in both the sequences of DNA.

GENERAL SIGNIFICANCE

The information obtained from this study will be useful for the designing of Gdm⁺ based drug molecule which can target the DNA more specifically and selectively.

How Significant Are Unusual Protein-Ligand Interactions? Insights from Database Mining

We present a new approach to derive interaction propensities of protein-ligand atom pairs from mining of the Protein Data Bank. To ensure solid statistics, we use a line-of-sight contact filter and normalize the observed frequency of hits by a statistical null model based on exposed surface areas of atom types in the protein-ligand binding site. This allows us to investigate which intermolecular interactions and geometries are found more often than expected by chance in protein-ligand complexes. We focus our study on some of the unusual interactions that were postulated to be favorable, including σ-hole bonding of halogen and sulfur atoms, weak hydrogen bonding with fluorine as acceptor, and different types of dipolar interactions. Our results confirm some and challenge other common assumptions on these interactions and highlight other contact types that are yet underexplored in structure-based drug design.

Cooperative nature of the sulfur centered hydrogen bond: investigation of (H2S)n (n = 2-4) clusters using an affordable yet accurate level of theory

Existing studies have shown that appreciably high level quantum chemical calculations are required to reproduce experimental energetic and geometric features of a H₂S dimer. This condition severely restricts any practical possibility of obtaining reliable results for higher order H₂S clusters. We have shown here that the binding energies calculated at the CCSD(T)/CBS level with counterpoise corrected geometries calculated at the MP2/aug-cc-pV(Q+d)Z level of theory excellently match with the experimental results for the H₂S dimer. Subsequently, the above mentioned levels of theory were used for trimers and tetramers. (H₂S)_n (n = 2-4) clusters were found to show cooperative strengthening of S-HS hydrogen bonds, which is clearly evident from the evolution of binding energies and hydrogen bond lengths, with increasing cluster size. Localized molecular orbital energy decomposition analyses have been carried out to understand how the contributions of various energy components modulate with the size of the clusters and what are their relative contributions towards the overall stabilization of the clusters. Natural bond orbital and atoms in molecules analyses were also carried out in order to look into the evolution of the electronic charge transfer and electron density topology with cluster size.

Superstrong Noncovalent Interface between Melamine and Graphene Oxide

There have been growing academic interests in the study of strong organic molecule-graphene [or graphene oxide (GO)] systems, owing to their essential noncovalent nature and the consequent chemomechanical behavior within the interface. A more recent experimental measurement [ Chem 2018, 4, 896-910] reported that the melamine-GO interface exhibits a remarkable noncovalent binding strength up to ∼1 nN, even comparable with typical covalent bonds. But the poor understanding on the complex noncovalent nature in particular makes it challenging to unveil the mystery of this high-performance interface. Herein, we carry out first-principles calculations to investigate the atomistic origin of ultrastrong noncovalent interaction between the melamine molecule and the GO sheet, as well as the chemomechanical synergy in interfacial behavior. The anomalous O-H···N hydrogen bonding, formed between the triazine moiety of melamine and the -OH in GO, is found cooperatively enhanced by the pin-like NH₂-π interaction, which is responsible for the strong interface. Following static pulling simulations validates the 1 nN level rupture strength and the contribution of each noncovalent interaction within the interface. Moreover, our results show that the -OH hydrogen bonding will mainly augment the interfacial adhesion strength, whereas the -NH₂ group cooperating with the -OH hydrogen bonding and conjugating with the GO surface will greatly improve the interfacial shear performance. Our work deepens the understanding on the chemomechanical behaviors within the noncovalent interface, which is expected to provide new potential strategies in designing high-performance graphene-based artificial nacreous materials.

The interplay and strength of the π⋯HF, C⋯HF, F⋯HF and F⋯HC hydrogen bonds upon the formation of multimolecular complexes based on C2H2⋯HF and C2H4⋯HF small dimers

The conception of this theoretical research was idealized aiming to unveil the intermolecular structures of complexes formed by acetylene or ethylene and hydrofluoric acid. At light of computational calculations by using the B3LYP/6-311++G(d,p) method, the geometries of the C₂H₂⋯(HF), C₂H₂⋯2(HF), C₂H₂⋯4(HF), C₂H₄⋯(HF), C₂H₄⋯2(HF) and C₂H₄⋯4(HF) hydrogen-bonded complexes were fully optimized. Moreover, the Post-Hartree-Fock calculations MP2/6-311++G(d,p), MP2/aug-cc-pVTZ, MP4(SDQ)/6-311++G(d,p) and CCSD/6-311++G(d,p) also were also used. The infrared spectra were analyzed in order to identify the new vibrational modes and frequencies of the proton donors shifted to red region. Through the modeling of charge-fluxes on the basis of the Quantum Theory of Atoms In Molecules (QTAIM) and, by contradicting the expectation of the hydrofluorination mechanisms of acetylene or ethylene, C⋯HF was recognized as a new type of hydrogen bond instead of the already well known π⋯H. The calculations of the Natural Bonding Orbital (NBO) and Charges derived from the Electrostatic Potential Grid-based (ChElPG) were also applied to interpret the shifting frequencies as well as measuring of the punctual charge-transfer after the formation of the complexes. Finally, the determination of the stabilization energy was carried out through the arguments of the Fock matrix in NBO basis and through the supermolecule approach. Also it is worthwhile to notice that some algebraic formulations were used for determining the electronic cooperative effect (CE).

Ab Initio, Transition State Theory, and Kinetic Modeling Study of the HO2-Assisted Keto-Enol Tautomerism Propen-2-ol + HO2 ⇔ Acetone + HO2 under Combustion, Atmospheric, and Interstellar Conditions

Keto-enol tautomerisms are important reactions in gaseous and liquid systems with implications in different chemical environments, but their kinetics have not been widely investigated. These reactions can proceed via a unimolecular process or may be catalyzed by another molecule. This work presents a theoretical study of the HO₂-catalyzed tautomerism that converts propen-2-ol into acetone at conditions relevant to combustion, atmospheric and interstellar chemistry. We performed CCSD(T)/aug-cc-pVTZ//M06-2X/cc-pVTZ ab initio and multistructural torsional variational transition state theory calculations to compute the forward and reverse rate constants. These rate constants have not been investigated previously, and modelers approximate the kinetics by comparison to analogue reactions. Two features of the potential energy surface of the studied tautomerism are highlighted. First, the HO₂ radical exhibits a pronounced catalytic effect by inducing a double hydrogen atom transfer reaction with a much lower barrier than that of the unimolecular process. Second, a prereactive complex is formed with a strong OH···π hydrogen bond. The role of the studied reaction under combustion conditions has been assessed via chemical kinetic modeling of 2-butanol (a potential alternative fuel) oxidation. The HO₂-assisted process was found to not be competitive with the unimolecular and HCOOH-assisted tautomerisms. The rate constants for the formation of the prereactive complex were calculated with the variable reaction coordinate transition state theory, and pressure effects were estimated with the system-specific quantum Rice-Ramsperger-Kassel theory; this allowed us to investigate the role of the complex by using the canonical unified statistical model. The formation and equilibration of the prereactive complex, which is also important at low pressures, enhances the reactivity by inducing a large tunneling effect that leads to a significant increase of the rate constants at cold and ultracold temperatures. These findings may help to understand and model the fate of complex organic molecules in the interstellar medium, and suggest an alternative route for the high energy barrier keto-enol tautomerism which otherwise is not kinetically favored at low temperatures.

C-Terminal Residue of Ultrashort Peptides Impacts on Molecular Self-Assembly, Hydrogelation, and Interaction with Small-Molecule Drugs

Single molecular changes on a tripeptide can have dramatic effects on their self-assembly and hydrogelation. Herein, we explore C-terminal residue variation on two consistent ultrashort peptide backbones, i.e. acetylated-Leu-Ile-Val-Ala-Gly-Xaa and acetylated-Ile-Val-Xaa (Xaa = His, Arg, Asn). The objective of this study is to identify candidates that can form hydrogels for small-molecule drug (SMD) delivery. Haemolysis and cytotoxicity (with human adipose-derived mesenchymal stem cells) assays showed that the new soluble peptides (Xaa = His, Arg) are cytocompatible. Gelation studies showed that all but acetylated-Ile-Val-Arg could gel under physiological conditions. Longer peptidic backbones drive self-assembly more effectively as reflected in field emission scanning electron microscopy (FESEM) and circular dichroism spectroscopy studies. Rheological studies revealed that the resultant hydrogels have varying stiffness and yield stress, depending on the backbone and C-terminal residue. Visible spectroscopy-based elution studies with SMDs (naltrexone, methotrexate, doxorubicin) showed that besides the C-terminal residue, the shape of the SMD also determines the rate and extent of SMD elution. Based on the elution assays, infrared spectroscopy, and FESEM, we propose models for the peptide fibril-SMD interaction. Our findings highlight the importance of matching the molecular properties of the self-assembling peptide and SMD in order to achieve the desired SMD release profile.

Similarities and differences in the crystal packing of halogen-substituted indole derivatives

The role of different intermolecular interactions in the crystal structures of halogen-substituted indoles which are fused with six-membered or seven-membered cyclic rings is investigated here. Several crystal structures show isostructural characteristics due to the presence of similar supramolecular motifs. In the absence of any strong hydrogen bonds, the molecular packing of reported structures is primarily stabilized by the presence of non-classical N-H...π and C-H...π interactions in addition to C-H...X (X = F/Cl/Br) interactions. The nature and energetics of primary and secondary dimeric motifs are partitioned into the electrostatics, polarization, dispersion and exchange-repulsion components using the PIXEL method. Short and directional N-H...π interactions are further explored by a topological analysis of the electron density based on quantum theory of atoms in molecules.

Improving the first hyperpolarizability of anthracene through interaction with HX molecules (XF, Cl, Br): A theoretical study

The variations in nonlinear optical activity (NLO) of anthracene (C₁₄H₁₀) was investigated via intermolecular interactions between C₁₄H₁₀ and HX molecules (XF, Cl and Br) using B3LYP-D3 method at 6-311++G(d,p) basis set. The stabilization of those complexes was investigated via vibrational analysis, quantum theory of atoms in molecules, molecular electrostatic potential, natural bond orbitals and symmetry-adapted perturbation theory (SAPT) analysis. Furthermore, the optical spectra and the first hyperpolarizabilities of C₁₄H₁₀⋯HX complexes were computed. The adsorption of hydrogen halide through C₁₄H₁₀⋯HX complex formation, didn't change much the linear optical activities of C₁₄H₁₀ molecule, but the magnitude of the first hyperpolarizability of the C₁₄H₁₀⋯HX complexes to be as much as that of urea.

Novel biocatalytic systems for maintaining the nucleotide balance based on adenylate kinase immobilized on carbon nanostructures

In this study graphene oxide (GO), carbon quantum dots (CQD) and carbon nanoonions (CNO) have been characterized and applied for the first time as a matrix for recombinant adenylate kinase (AK, EC 2.7.4.3) immobilization. AK is an enzyme fulfilling a key role in metabolic processes. This phosphotransferase catalyzes the interconversion of adenine nucleotides (ATP, ADP and AMP) and thereby participates in nucleotide homeostasis, monitors a cellular energy charge as well as acts as a component of purinergic signaling system. The AK activity in all obtained biocatalytic systems was higher as compared to the free enzyme. We have found that the immobilization on carbon nanostructures increased both activity and stability of AK. Moreover, the biocatalytic systems consisting of AK immobilized on carbon nanostructures can be easily and efficiently lyophilized without risk of desorption or decrease in the catalytic activity of the investigated enzyme. The positive action of AK-GO biocatalytic system in maintaining the nucleotide balance in in vitro cell culture was proved.

Towards understanding the E. coli PNP binding mechanism and FRET absence between E. coli PNP and formycin A

The aim of this study is threefold: (1) augmentation of the knowledge of the E. coli PNP binding mechanism; (2) explanation of the previously observed 'lack of FRET' phenomenon and (3) an introduction of the correction (modified method) for FRET efficiency calculation in the PNP-FA complexes. We present fluorescence studies of the two E. coli PNP mutants (F159Y and F159A) with formycin A (FA), that indicate that the aromatic amino acid is indispensable in the nucleotide binding, additional hydroxyl group at position 159 probably enhances the strength of binding and that the amino acids pair 159-160 has a great impact on the spectroscopic properties of the enzyme. The experiments were carried out in hepes and phosphate buffers, at pH7 and 8.3. Two methods, a conventional and a modified one, that utilizes the dissociation constant, for calculations of the energy transfer efficiency (E) and the acceptor-to-donor distance (r) between FA and the Tyr (energy donor) were employed. Total difference spectra were calculated for emission spectra (λ_ex 280nm, 295nm, 305nm and 313nm) for all studied systems. Time-resolved techniques allowed to conclude the existence of a specific structure formed by amino acids at positions 159 and 160. The results showed an unexpected pattern change of FRET in the mutants, when compared to the wild type enzyme and a probable presence of a structure created between 159 and 160 residue, that might influence the binding efficiency. Additionally, we confirmed the indispensable role of the modification of the FRET efficiency (E) calculation on the fraction of enzyme saturation in PNP-FA systems.

Interplay among Electrostatic, Dispersion, and Steric Interactions: Spectroscopy and Quantum Chemical Calculations of π-Hydrogen Bonded Complexes

π-Hydrogen bonding interactions are ubiquitous in both materials and biology. Despite their relatively weak nature, great progress has been made in their investigation by experimental and theoretical methods, but this becomes significantly more complicated when secondary intermolecular interactions are present. In this study, the effect of successive methyl substitution on the supramolecular structure and interaction energy of indole⋅⋅⋅methylated benzene (ind⋅⋅⋅n-mb, n=1-6) complexes is probed through a combination of supersonic jet experiments and benchmark-quality quantum chemical calculations. It is demonstrated that additional secondary interactions introduce a subtle interplay among electrostatic and dispersion forces, as well as steric repulsion, which fine-tunes the overall structural motif. Resonant two-photon ionization and IR-UV double-resonance spectroscopy techniques are used to probe jet-cooled ind⋅⋅⋅n-mb (n=2, 3, 6) complexes, with redshifting of the N-H IR stretching frequency showing that increasing the degree of methyl substitution increases the strength of the primary N-H⋅⋅⋅π interaction. Ab initio harmonic frequency and binding energy calculations confirm this trend for all six complexes. Electronic spectra of the three dimers are broad and structureless, with quantum chemical calculations revealing that this is likely to be due to multiple tilted conformations of each dimer possessing similar stabilization energies.

Synergistic Coassembly of Two Structurally Different Molecular Gelators

Coassembly of molecules can produce materials with improved properties and functionalities. To this end, achieving a molecular level understanding of the interactions governing the coassembly is essential. In this work, two molecular gelators with significantly different structures and main intermolecular forces for assembly were coassembled. The elastic moduli of the hybrid gels are more than 1 order of magnitude higher than those of the gels formed by the individual gelators, showing an obvious synergistic effect. The interactions between the gelators were investigated with confocal microscopy and both one-dimensional and two-dimensional nuclear magnetic resonance. It was found that the two gelators coassemble to form fibers due to the nonspecific van der Waals interactions between their alkyl chains and the specific interactions between their functional groups. Switching from one gelator-dominated fiber network to the other gelator-dominated fiber network was achieved at a critical molar ratio of the gelators. The two gelators serve as additives of each other to tune the nucleation and growth of the fiber networks. The observations of this work are significant to the development of materials with improved properties by coassembly of different molecules.

Direct Spectroscopic Evidence for an n→π* Interaction

The n→π* interaction is an extremely weak but very important noncovalent interaction. Although this interaction is widely present in biomolecules and materials, its existence is counterintuitive and so has been debated extensively. Herein, we have reported direct spectroscopic evidence for an n→π* interaction for the first time by probing the carbonyl stretching frequency in phenyl formate using isolated gas-phase IR spectroscopy. This result also demonstrates that the conformational preference for the cis conformer of phenyl formate compared to the trans conformer arises due to the presence of the n→π* interaction in the former. The direct proof reported herein for this controversial but important noncovalent interaction should stimulate further experimental and theoretical investigation on this intriguing research topic.