Revealing Non-covalent Interactions in Molecular Crystals through Their Experimental Electron Densities

Formation of left-handed helices by C2'-fluorinated nucleic acids under physiological salt conditions

Recent findings in cell biology have rekindled interest in Z-DNA, the left-handed helical form of DNA. We report here that two minimally modified nucleosides, 2'F-araC and 2'F-riboG, induce the formation of the Z-form under low ionic strength. We show that oligomers entirely made of these two nucleosides exclusively produce left-handed duplexes that bind to the Zα domain of ADAR1. The effect of the two nucleotides is so dramatic that Z-form duplexes are the only species observed in 10 mM sodium phosphate buffer and neutral pH, and no B-form is observed at any temperature. Hence, in contrast to other studies reporting formation of Z/B-form equilibria by a preference for purine glycosidic angles in syn, our NMR and computational work revealed that sequential 2'F…H2N and intramolecular 3'H…N3' interactions stabilize the left-handed helix. The equilibrium between B- and Z- forms is slow in the 19F NMR time scale (≥ms), and each conformation exhibited unprecedented chemical shift differences in the 19F signals. This observation led to a reliable estimation of the relative population of B and Z species and enabled us to monitor B-Z transitions under different conditions. The unique features of 2'F-modified DNA should thus be a valuable addition to existing techniques for specific detection of new Z-binding proteins and ligands.

Exploring the effects of mono-bromination on hole-electron transport and distribution in dibenzofuran and dibenzothiophene isomers: a first-principles study

CONTEXT

This study delves into hole-electron transport and distribution properties inherent in mono-brominated dibenzofuran (DBF) and dibenzothiophene (DBT) isomers. As determined by frontier molecular orbitals, all brominated structures have narrower bandgaps than their primary structures. The TD-DFT calculation showed that 2BDBT had the highest absorption wavelength of all molecules at 315.35 nm. Notably, the study unveils remarkably low electron and hole reorganization energies due to bromine substitution in DBF and DBT molecules. Specifically, the 4BDBF has the lowest hole reorganization energy of all DBF configurations, 0.229 eV. In addition, 3BDBF has 0.226 eV less electron reorganization energy than all other molecules. Compared to DBT, 3BDBT has the lowest electron reorganization energy of 0.254 eV. Overall, this research sheds significant light on the fundamental electronic and hole transport characteristics of bromine-substituted DBF and DBT isomers, highlighting their promising role in polymer design as donors/acceptors for advanced organic electronic applications.

METHODS

Molecular structures were optimized using Density Functional Theory (DFT) B3LYP/6-311 +  + G (d, p) level of theory, and the study further elucidates these molecules' energy levels and absorption spectra through Time-Dependent Density Functional Theory TD-DFT; these calculations were performed using Gaussian 09W software package. The key parameters such as reorganization energies, Electron Localization Function map, Laplacian Bond Order, and NCI-RDG were meticulously examined for the molecules with the results of DFT calculations were analyzed and displayed by utilizing the software packages VMD 1.9.4 and Multiwfn 3.8, aiming to comprehend their charge transport and distribution properties.

Exploring the Sensing Potential of g-C3N4 versus Li/g-C3N4 Nanoflakes toward Hazardous Organic Volatiles: A DFT Simulation Study

Ab initio calculations were performed to determine the sensing behavior of g-C3N4 and Li metal-doped g-C3N4 (Li/g-C3N4) quantum dots toward toxic compounds acetamide (AA), benzamide (BA), and their thio-analogues, namely, thioacetamide (TAA) and thiobenzamide (TAA). For optimization and interaction energies, the ωB97XD/6-31G(d,p) level of theory was used. Interaction energies (Eint) illustrate the high thermodynamic stabilities of the designed complexes due to the presence of the noncovalent interactions. The presence of electrostatic forces in some complexes is also observed. The observed trend of Eint in g-C3N4 complexes was BA > TAA > AA > TBA, while in Li/g-C3N4, the trend was BA > AA > TBA > TAA. The electronic properties were studied by frontier molecular orbital (FMO) and natural bond orbital analyses. According to FMO, lithium metal doping greatly enhanced the conductivity of the complexes by generating new HOMOs near the Fermi level. A significant amount of charge transfer was also observed in complexes, reflecting the increase in charge conductivity. NCI and QTAIM analyses evidenced the presence of significant noncovalent dispersion and electrostatic forces in Li/g-C3N4 and respective complexes. Charge decomposition analysis gave an idea of the transfer of charge density between quantum dots and analytes. Finally, TD-DFT explained the optical behavior of the reported complexes. The findings of this study suggested that both bare g-C3N4 and Li/g-C3N4 can effectively be used as atmospheric sensors having excellent adsorbing properties toward toxic analytes.

Revealing the Role of Noncovalent Interactions on the Conformation of the Methyl Group in Tricyclic Orthoamide

Tricyclic orthoamides are valuable molecules with wide-ranging applications, including organic synthesis and molecular recognition. Their structural properties make them intriguing, particularly the eclipsed all-trans conformer, which is typically less stable than the alternated conformation and is a rare phenomenon in organic chemistry. However, it gains stability in crystalline and hydrated settings, challenging the existing theoretical explanations. This study investigates which factors make eclipsed conformers more stable using experimentally reported anhydrous (ATO) and hydrated (HTO) crystal structures. Employing the quantum theory of atoms in molecules, noncovalent interaction index, and pairwise energy decomposition analysis, we delve into the noncovalent interaction environment surrounding the molecule of interest. In ATO, dispersive interactions dominate, whereas in HTO, both dispersive and electrostatic contributions are observed due to the presence of water molecules. Anchored to the lone pairs of the nitrogen atom in the orthoamide tricycle, water molecules prompt the methyl group's eclipsing through intermolecular and intramolecular interactions. This work resolves the long-standing conflict behind why tricyclic orthoamide has an eclipsed conformation by establishing the stabilization factors. These insights have implications for crystal engineering and design, enhancing our understanding of structural behavior in both crystalline and hydrated environments.

Double Centrosymmetric Si···π Tetrel Bonds as New Synthons─Evidence from Crystal Structures and DFT Calculations

The crystal structure of bis((μ2-ethynylsilyloxo)-dichloro-aluminum), BEDCA, and a few related structures are characterized by the occurrence of tetrel bonds that link molecules. Particularly, centosymmetric dimers in such structures occur that are connected by two equivalent Si···π tetrel bonds. The dimer of BEDCA and dimers of other model species that similarly are linked by two equivalent Si···π tetrel bonds are analyzed theoretically. Some of the complexes calculated here are also characterized by the occurrence of triel bonds. Thus, ωB97XD/aug-cc-pVTZ calculations are performed and these DFT results are further supported by calculations with the use of other theoretical approaches: the quantum theory of atoms in molecules, QTAIM; the natural bond orbital, NBO; the energy decomposition analysis, EDA; and the noncovalent interactions method, NCI. The results show that the tetrel bonds analyzed here are rather weak, and they are not detected by the QTAIM approach; however, they are detected by other approaches, like NBO, for example. On the other hand, the triel bonds that occur in a few complexes discussed here are very strong and possess characteristics of covalent bonds.

Combined Experimental and Computational Study of V-Substituted Lindqvist Polyoxotungstate: Screening by Docking for Potential Antidiabetic Activity

In the current work, a novel vanadotungstate compound, (C6H9N2)4[V2W4O19]·2H2O (1), is isolated by a simple stepwise synthesis method and characterized by a combined experimental and computational study. Molecular docking is conducted for the first time for this kind of substituted Lindqvist polyoxometalates to elucidate for potential antidiabetic activity. Hence, the modeling results revealed a significant docking score of the reported compound to bind to the active sites of α-glucosidase with the lowest binding energy of -5.7 kcal/mol, where the standard drug acarbose (ACB) had -4.6 kcal/mol binding energy. The stability of binding was enhanced by strong H-bonding, van der Waals, and electrostatic interactions occurring in the three-dimensional (3D) supramolecular network of polyanionic vanadotungstate subunits templated with organic moieties as shown by X-ray diffraction and Hirshfeld analyses. Furthermore, density functional theory (DFT) calculations supported with photophysical measurements are also discussed to predict the most chemical and biological reactivity. In this view, the complete description of electronic and biological features of (1) is enhanced by determination of the highest occupied molecular orbital (HOMO)/least unoccupied molecular orbital (LUMO) energy, electronic density, ionization potential, electron affinity, etc. These chemical descriptors, intermolecular interactions, docking score, and binding free energy estimation are essential in understanding the reactivity of this bioactive compound offering potential inhibition of the α-glucosidase enzyme.

Dissolution mechanism of cellulose in a benzyltriethylammonium/urea deep eutectic solvent (DES): DFT-quantum modeling, molecular dynamics and experimental investigation

In this paper, a benzyltriethylammonium/urea DES was investigated as a new green and eco-friendly medium for the progress of organic chemical reactions, particularly the dissolution and the functionalization of cellulose. In this regard, the viscosity-average molecular weight of cellulose (M̄w) during the dissolution/regeneration process was investigated, showing no significant degradation of the polymer chains. Moreover, X-ray diffraction patterns indicated that the cellulose dissolution process in the BTEAB/urea DES decreased the crystallinity index from 87% to 75%, and there was no effect on type I cellulose polymorphism. However, a drastic impact of the cosolvents (water and DMSO) on the melting point of the DES was observed. Besides, to understand the evolution of cellulose-DES interactions, the formation mechanism of the system was studied in terms of H-bond density and radial distribution function (RDF) using molecular dynamics modeling. Furthermore, density functional theory (DFT) was used to evaluate the topological characteristics of the polymeric system such as potential energy density (PED), laplacian electron density (LED), energy density, and kinetic energy density (KED) at bond critical points (BCPs) between the cellulose and the DES. The quantum theory of atoms in molecules (AIM), Bader's quantum theory (BQT), and reduced density gradient (RDG) scatter plots have been exploited to estimate and locate non-covalent interactions (NCIs). The results revealed that the dissolution process is attributed to the physical interactions, mainly the strong H-bond interactions.

The NCIWEB Server: A Novel Implementation of the Noncovalent Interactions Index for Biomolecular Systems

It is well-known that the activity and function of proteins is strictly correlated with their secondary, tertiary, and quaternary structures. Their biological role is regulated by their conformational flexibility and global fold, which, in turn, is largely governed by complex noncovalent interaction networks. Because of the large size of proteins, the analysis of their noncovalent interaction networks is challenging, but can provide insights into the energetics of conformational changes or protein-protein and protein-ligand interactions. The noncovalent interaction (NCI) index, based on the reduced density gradient, is a well-established tool for the detection of weak contacts in biological systems. In this work, we present a web-based application to expand the use of this index to proteins, which only requires a molecular structure as input and provides a mapping of the number, type, and strength of noncovalent interactions. Structure preparation is automated and allows direct importing from the PDB database, making this server (https://nciweb.dsi.upmc.fr) accessible to scientists with limited experience in bioinformatics. A quick overview of this tool and concise instructions are presented, together with an illustrative application.

Design, synthesis, and computational studies of novel imidazo[1,2-a]pyrimidine derivatives as potential dual inhibitors of hACE2 and spike protein for blocking SARS-CoV-2 cell entry

In the present work, a new series of imidazo[1,2-a]pyrimidine Schiff base derivatives have been obtained using an easy and conventional synthetic route. The synthesized compounds were spectroscopically characterized using ¹H, ¹³C NMR, LC-MS(ESI), and FT-IR techniques. Green metric calculations indicate adherence to several green chemistry principles. The energy of Frontier Molecular Orbitals (FMO), Molecular Electrostatic Potential (MEP), Quantum Theory of Atoms in Molecules (QTAIM), and Reduced Density Gradient (RDG) were determined by the Density Functional Theory (DFT) method at B3LYP/6-31 G (d, p) as the basis set. Moreover, molecular docking studies targeting the human ACE2 and the spike, key entrance proteins of the severe acute respiratory syndrome coronavirus-2 were carried out along with hACE2 natural ligand Angiotensin II, the MLN-4760 inhibitor as well as the Cannabidiolic Acid CBDA which has been demonstrated to bind to the spike protein and block cell entry. The molecular modeling results showed auspicious results in terms of binding affinity as the top-scoring compound exhibited a remarkable affinity (-9.1 and -7.3 kcal/mol) to the ACE2 and spike protein respectively compared to CBDA (-5.7 kcal/mol), the MLN-4760 inhibitor (-7.3 kcal/mol), and angiotensin II (-9.2 kcal/mol). These findings suggest that the synthesized compounds may potentially act as effective entrance inhibitors, preventing the SARS-CoV-2 infection of human cells. Furthermore, in silico, ADMET, and drug-likeness prediction expressed promising drug-like characteristics.

Charge equilibration model with shielded long-range Coulomb for reactive molecular dynamics simulations

Atomic description of electrochemical systems requires reactive interaction potential to explicitly describe the chemistry between atoms and molecules and the evolving charge distribution and polarization effects. Calculating Coulomb electrostatic interactions and polarization effects requires a better estimate of the partial charge distribution in molecular systems. However, models such as reactive force fields and charge equilibration (QEq) include Coulomb interactions up to a short-distance cutoff for better computational speeds. Ignoring long-distance electrostatic interaction affects the ability to describe electrochemistry in large systems. We studied the long-range Coulomb effects among charged particles and extended the QEq method to include long-range effects. By this extension, we anticipate a proper account of Coulomb interactions in reactive molecular dynamics simulations. We validate the approach by computing charges on a series of metal-organic frameworks and some simple systems. Results are compared to regular QEq and quantum mechanics calculations. The study shows slightly overestimated charge values in the regular QEq approach. Moreover, our method was combined with Ewald summation to compute forces and evaluate the long-range effects of simple capacitor configurations. There were noticeable differences between the calculated charges with/without long-range Coulomb interactions. The difference, which may have originated from the long-range influence on the capacitor ions, makes the Ewald method a better descriptor of Coulomb electrostatics for charged electrodes. The approach explored in this study enabled the atomic description of electrochemical systems with realistic electrolyte thickness while accounting for the electrostatic effects of charged electrodes throughout the dielectric layer in devices like batteries and emerging solid-state memory.

Interactions between Paracetamol and Formaldehyde: Theoretical Investigation and Topological Analysis

In this work, noncovalent interactions including hydrogen bonds, C···C, N···O, and van der Waals forces between paracetamol and formaldehyde were investigated using the second-order perturbation theory MP2 in conjunction with the correlation consistent basis sets (aug-cc-pVDZ and aug-cc-pVTZ). Two molecular conformations of paracetamol were considered. Seven equilibrium geometries of dimers were found from the result of the interactions with formaldehyde for each conformation of paracetamol. Interaction energies of complexes with both ZPE and BSSE corrections range from -7.0 to -21.7 kJ mol^-1. Topological parameters (such as electron density, its Laplacian, and local electron energy density at the bond critical points) of the bonds from atoms in molecules theory were analyzed in detail. The natural bond orbital analysis showed that the stability of complexes was controlled by noncovalent interactions including O-H···O, N-H···O, C-H···O, C-H···N, C-H···H-C, C···C, and N···O. The red- and blue-shifted hydrogen bonds could both be observed in these complexes. The properties of these interactions were also further examined in water using a polarized continuum model. In water, the stability of the complex was slightly reduced as compared to that in the gas phase.

Exploring the Non-Covalent Bonding in Water Clusters

QTAIM and source function analysis were used to explore the non-covalent bonding in twelve different water clusters (H2O)n obtained by considering n = 2–7 and various geometrical arrangements. A total of seventy-seven O−H⋯O hydrogen bonds (HBs) were identified in the systems under consideration, and the examination of the electron density at the bond critical point (BCP) of these HBs revealed the existence of a great diversity of O−H⋯O interactions. Furthermore, the analysis of quantities, such as |V(r)|/G(r) and H(r), allowed a further description of the nature of analogous O−H⋯O interactions within each cluster. In the case of 2-D cyclic clusters, the HBs are nearly equivalent between them. However, significant differences among the O−H⋯O interactions were observed in 3-D clusters. The assessment of the source function (SF) confirmed these findings. Finally, the ability of SF to decompose the electron density (ρ) into atomic contributions allowed the evaluation of the localized or delocalized character of these contributions to ρ at the BCP associated to the different HBs, revealing that weak O−H⋯O interactions have a significant spread of the atomic contributions, whereas strong interactions have more localized atomic contributions. These observations suggest that the nature of the O−H⋯O hydrogen bond in water clusters is determined by the inductive effects originated by the different spatial arrangements of the water molecules in the studied clusters.

Vibrational Spectroscopies, Global Reactivity, Molecular Docking, Thermodynamic Properties and Linear and Nonlinear Optical Parameters of Monohydrate Arsonate Salt of 4-Aminopyridine

In this work, a structural and electronic properties of a novel organic arsenate template by 4-aminopyridine, with the general formula (C₅H₇N₂)(C₅H₈N₂)[AsO₄]·H₂O ((4-APH)(4-APH₂)[AsO₄]·H₂O) have been presented. The density functional theory (DFT) along with B3LYP hybrid functional is employed. The optimized structure was found to be in well consistent with the X-ray diffraction geometry. The examination of the vibrational spectrum was correlated by DFT calculation using the unit cell parameters obtained from the experiment data. Besides, the thermodynamic functions (heat capacity, entropy, enthalpy) from spectroscopic data by statistical methods were obtained for the range of temperature 100–1000 K. In addition, the molecular orbital calculations such as Natural Bond Orbitals (NBOs), AIM approach, HOMO–LUMO energy gap, NLO characteristic and Hirshfeld surface analysis were also performed with the same level of DFT. Electronic stability of the compound arising from hyper conjugative interactions and charge delocalization were also investigated based on the natural bond orbital (NBO) analysis. Molecular docking studies were also conducted as part of this study. The theoretical results showed an excellent agreement with the experimental values.

Extracting Quantitative Information at Quantum Mechanical Level from Noncovalent Interaction Index Analyses

The noncovalent interaction (NCI) index is nowadays a well-known strategy to detect NCIs in molecular systems. Even though it initially provided only qualitative descriptions, the technique has been recently extended to also extract quantitative information. To accomplish this task, integrals of powers of the electron distribution were considered, with the requirement that the overall electron density can be clearly decomposed as sum of distinct fragment contributions to enable the definition of the (noncovalent) integration region. So far, this was done by only exploiting approximate promolecular electron densities, which are given by the sum of spherically averaged atomic electron distributions and thus represent too crude approximations. Therefore, to obtain more quantum mechanically (QM) rigorous results from NCI index analyses, in this work, we propose to use electron densities obtained through the transfer of extremely localized molecular orbitals (ELMOs) or through the recently developed QM/ELMO embedding technique. Although still approximate, the electron distributions resulting from the abovementioned methods are fully QM and, above all, are again partitionable into subunit contributions, which makes them completely suitable for the NCI integral approach. Therefore, we benchmarked the integrals resulting from NCI index analyses (both those based on the promolecular densities and those based on ELMO electron distributions) against interaction energies computed at a high quantum chemical level (in particular, at the coupled cluster level). The performed test calculations have indicated that the NCI integrals based on ELMO electron densities outperform the promolecular ones. Furthermore, it was observed that the novel quantitative NCI-(QM/)ELMO approach can be also profitably exploited both to characterize and evaluate the strength of specific interactions between ligand subunits and protein residues in protein-ligand complexes and to follow the evolution of NCIs along trajectories of molecular dynamics simulations. Although further methodological improvements are still possible, the new quantitative ELMO-based technique could be already exploited in situations in which fast and reliable assessments of NCIs are crucial, such as in computational high-throughput screenings for drug discovery.

Quantum chemical studies on hydrogen bonds in helical secondary structures

We present a brief review of our recent computational studies of hydrogen bonds (H-bonds) in helical secondary structures of proteins, α-helix and 3₁₀-helix, using a Negative Fragmentation Approach with density functional theory. We found that the depolarized electronic structures of the carbonyl oxygen of the ith residue and the amide hydrogen of the (i + 4)th residue cause weaker H-bond in an α-helix than in an isolated H-bond. Our calculations showed that the H-bond energies in the 3₁₀-helix were also weaker than those of the isolated H-bonds. In the 3₁₀-helices, the adjacent N-H group at the (i + 1)th residue was closer to the C=O group of the H-bond pair than the adjacent C=O group in the 3₁₀-helices, whereas the adjacent C=O group at the (i + 1)th residue was close to the H-bond acceptor in α-helices. Therefore, the destabilization of the H-bond is attributed to the depolarization caused by the adjacent residue of the helical backbone connecting the H-bond donor and acceptor. The differences in the change in electron density revealed that such depolarizations were caused by the local electronic interactions in their neighborhood inside the helical structure and redistributed the electron density. We also present the improvements in the force field of classical molecular simulation, based on our findings.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01034-5.

Modeling the structure and reactivity landscapes of a pyrazole-ammonium ionic derivative using wavefunction-dependent characteristics and screening for potential anti-inflammatory activity

Spectroscopic investigations of 1-phenyl -2,3-dimethyl-5-oxo-1,2-dihydro-1H-pyrazol-4-ammonium 2[(2-carboxyphenyl) disulfanyl]benzoate (PACB) reported experimentally and theoretically. NH-O interaction is observed and there is a very large downshift for NH-O stretching frequency. Reactive sites are identified from the chemical and electronic properties. For PACB the maximum repulsion was around H33, H55 and H57 atom. LOL shows red regions between C-C and blue around C atoms are surrounded by a delocalized electron cloud. The red ring is a hallmark of electron density depletion from the NCI plot due to electrostatic repulsion and its existences suggests that coordination sphere for PACB is minimally strained around the central ion. Atomic contact energy values and high score of the docking results obtained propose that, PACB may have inhibitory properties and have a significant function in pharmacological chemistry. Molecular dynamics simulation was performed to validate the stability of the title compound with the Bovine thrombin-activatable fibrinolysis inhibitor protein.Communicated by Ramaswamy H. Sarma.

Supramolecular Diversity, Theoretical Investigation and Antibacterial Activity of Cu, Co and Cd Complexes Based on the Tridentate N,N,O-Schiff Base Ligand Formed In Situ

The four new complexes, [Cu(HL¹)(L²)Cl] (1), [Cu(HL¹)(L¹)]∙Cl∙2H₂O (2), [Co(L¹)₂]∙Cl (3) and [Cd(HL¹)I₂]∙dmso (4), have been prepared by one-pot reactions of the respective chloride or iodide metal salt with a non-aqueous solution of the polydentate Schiff base, HL¹, resulted from in situ condensation of benzhydrazide and 2-pyridinecarboxaldehyde, while a ligand HL², in case of 1, has been formed due to the oxidation of 2-pyridinecarboxaldehyde under reaction conditions. The crystallographic analysis revealed that the molecular building units in 1-4 are linked together into complex structures by hydrogen bonding, resulting in 1D, 2D and 3D supramolecular architectures for 1, 2 and 4, respectively, and the supramolecular trimer for 3. The electronic structures of 1-4 were investigated by the DFT theoretical calculations. The non-covalent interactions in the crystal structures of 1-4 were studied by means of the Hirshfeld surface analysis and the QTAIM theory with a special focus on the C-H⋯Cl bonding. From the DFT/DLPNO-CCSD(T) calculations, using a series of charged model {R₃C-H}⁰⋯Cl^- assemblies, we propose linear regressions for assessment of the interaction enthalpy (ΔH, kcal mol^-1) and the binding energy (BE, kcal mol^-1) between {R₃C-H}⁰ and Cl^- sites starting from the electron density at the bond critical point (ρ(r_BCP), a.u.): ΔH = -678 × ρ(r) + 3 and BE = -726 × ρ(r) + 4. It was also has been found that compounds 1, 3 and 4 during in vitro screening showed an antibacterial activity toward the nine bacteria species, comprising both Gram-positive and Gram-negative, with MIC values ranging from 156.2 to 625 mg/L. The best results have been obtained against Acinetobacter baumannii MβL.

Depolarizing Effects in Hydrogen Bond Energy in 310-Helices Revealed by Quantum Chemical Analysis

Hydrogen-bond (H-bond) energies in 3₁₀-helices of short alanine peptides were systematically examined by precise DFT calculations with the negative fragmentation approach (NFA), a modified method based on the molecular tailoring approach. The contribution of each H-bond was evaluated in detail from the 3₁₀-helical conformation of total energies (whole helical model, WH_3-10 model), and the results were compared with the property of H-bond in α-helix from our previous study. The H-bond energies of the WH_3-10 model exhibited tendencies different from those exhibited by the α-helix in that they depended on the helical position of the relevant H-bond pair. H-bond pairs adjacent to the terminal H-bond pairs were observed to be strongly destabilized. The analysis of electronic structures indicated that structural characteristics cause the destabilization of the H-bond in 3₁₀-helices. We also found that the longer the helix length, the more stable the H-bond in the terminal pairs of the WH_3-10 model, suggesting the action of H-bond cooperativity.

Ratiometric pH Sensing, Photophysics, and Cell Imaging of Nonaromatic Light-Emitting Polymers

Here, four nontraditional fluorescent polymers (NTFPs) of varying N,N-dimethyl-2-propenamide (DMPA) and butyl prop-2-enoate (BPE) mole ratios, i.e., 2:1 (NTFP1), 4:1 (NTFP2), 8:1 (NTFP3), and 16:1 (NTFP4), are prepared via random polymerization in water. The maximum fluorescence enhancement of NTFP3 makes it suitable for ratiometric pH sensing, Cu(II) sensing, and pH-dependent cell imaging of Madin-Darby canine kidney (MDCK) cells. The oxygen donor functionalities of NTFP3 involved in binding and sensing with Cu(II) ions are studied by absorption, emission, electron paramagnetic resonance, Fourier transform infrared (FTIR), and O1s/Cu2p X-ray photoelectron spectroscopies (XPS). The spectral responses of the ratiometric pH sensor within 1.5-11.5 confirm 22 and 44 nm red shifts in absorption and ratiometric emission, respectively. The striking color changes from blue (436 nm) to green (480 nm) via an increase in pH are thought to be the stabilization of the charged canonical form of tertiary amide, i.e., -C(O^-)═N⁺(CH₃)₂, realized from the changes in the absorption/fluorescence spectra and XPS/FTIR analyses. The through-space n-π* interactions in the NTFP3 aggregate, N-branching-associated rigidity, and nonconventional intramolecular hydrogen bondings of adjacent NTFP3 moieties in the NTFP3 aggregate contribute to aggregation-enhanced emissions (AEEs). Here, structures of NTFP3, NTFP3 aggregate, and Cu(II)-NTFP3; absorption; n-π* interactions; hydrogen bondings; AEEs; and binding with Cu(II) are ascertained by density functional theory, time-dependent density functional theory, and reduced density gradient calculations. The excellent limits of detection and Stern-Volmer constants of NTFP3 are 2.24 nM/0.14234 ppb and 4.26 × 10³ M^-1 at pH = 6.5 and 0.95 nM/0.06037 ppb and 4.90 × 10³ M^-1 at pH = 8.0, respectively. Additionally, the Stokes shift and binding energy of NTFP3 are 13,636 cm^-1/1.69 eV and -4.64 eV, respectively. The pH-dependent MDCK cell imaging ability of noncytotoxic NTFP3 is supported via fluorescence imaging and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay.

Si-H···Se Chalcogen-Hydride Bond Quantified by Diffraction and Topological Analyses

The Si-H···Se contact in 1-mesitylselanyl-8-(dimethylsilyl)naphthalene (1), which exhibits the spatial arrangement of a δ-agostic interaction from geometric considerations, was investigated. Is this just enforced by close 1,8-proximity or is this a favorable interaction? Charge density studies are best suited to investigate the exact origin of the interaction and to quantify the properties. Hence, they are most elucidating. High-resolution X-ray diffraction data of 1 were collected, and a multipole refinement followed by a topological analysis using Bader's quantum theory of atoms in molecules was employed. The resulting bond properties were set in relation to high-level computational parameters. The comparison to Si-H···[M] agostics, hydride bonding, chalcogen bonds, and charge-inverted hydrogen bonds qualified the Si-H···Se noncovalent interaction to be best classified as a chalcogen-hydride bond.

Chalcogen Bonding in the Molecular Dimers of WCh2 (Ch = S, Se, Te): On the Basic Understanding of the Local Interfacial and Interlayer Bonding Environment in 2D Layered Tungsten Dichalcogenides

Layered two-dimensional transition metal dichalcogenides and their heterostructures are of current interest, owing to the diversity of their applications in many areas of materials nanoscience and technologies. With this in mind, we have examined the three molecular dimers of the tungsten dichalcogenide series, (WCh₂)₂ (Ch = S, Se, Te), using density functional theory to provide insight into which interactions, and their specific characteristics, are responsible for the interfacial/interlayer region in the room temperature 2H phase of WCh₂ crystals. Our calculations at various levels of theory suggested that the Te···Te chalcogen bonding in (WTe₂)₂ is weak, whereas the Se···Se and S···S bonding interactions in (WSe₂)₂ and (WS₂)₂, respectively, are of the van der Waals type. The presence and character of Ch···Ch chalcogen bonding interactions in the dimers of (WCh₂)₂ are examined with a number of theoretical approaches and discussed, including charge-density-based approaches, such as the quantum theory of atoms in molecules, interaction region indicator, independent gradient model, and reduced density gradient non-covalent index approaches. The charge-density-based topological features are shown to be concordant with the results that originate from the extrema of potential on the electrostatic surfaces of WCh₂ monomers. A natural bond orbital analysis has enabled us to suggest a number of weak hyperconjugative charge transfer interactions between the interacting monomers that are responsible for the geometry of the (WCh₂)₂ dimers at equilibrium. In addition to other features, we demonstrate that there is no so-called van der Waals gap between the monolayers in two-dimensional layered transition metal tungsten dichalcogenides, which are gapless, and that the (WCh₂)₂ dimers may be prototypes for a basic understanding of the physical chemistry of the chemical bonding environments associated with the local interfacial/interlayer regions in layered 2H-WCh₂ nanoscale systems.

Structural, surface, and computational analysis of two vitamin-B1 crystals with sulfonimide-based anions

Two crystals incorporating the thiamine·HCl cation and the fluorinated anion 1,3-disulfonylhexafluoropropyleneimide have been characterized via single-crystal X-ray diffraction. The host-guest interactions of thiamine with the anions are analyzed and characterized using Hirshfeld surface analysis. The cations in both structures form a dimer in the solid-state via reciprocal hydrogen bonding through the amine and hydroxyl moieties. Additional investigation into the interactions responsible for dimer formation found that the sulfur atom in the thiazolium ring interacting with several hydrogen atoms to form stabilizing interactions. These interactions in the dimer are further analyzed using reduced density gradient analysis and the results are correlated to the fingerprint plots derived from the Hirshfeld surfaces. Moreover, specific interactions are observed from the cyclical anions, with both the fluorine and sulfonyl oxygen atoms participating in bridging interactions, displaying the diverse host-guest properties of thiamine.

Orbital-Free Quantum Crystallographic View on Noncovalent Bonding: Insights into Hydrogen Bonds, π⋅⋅⋅π and Reverse Electron Lone Pair⋅⋅⋅π Interactions

A detailed analysis of a complete set of the local potentials that appear in the Euler equation for electron density is carried out for noncovalent interactions in the crystal of a uracil derivative using experimental X-ray charge density. The interplay between the quantum theory of atoms in molecules and crystals and the local potentials and corresponding inner-crystal electronic forces of electrostatic and kinetic origin is explored. Partitioning of crystal space into atomic basins and atomic-like potential basins led us to the definite description of interatomic interaction and charge transfer. Novel physically grounded bonding descriptors derived within the orbital-free quantum crystallography provided the detailed examination of π-stacking and intricate C=O⋅⋅⋅π interactions and nonclassical hydrogen bonds present in the crystal. The donor-acceptor character of these interactions is revealed by analysis of Pauli and von Weizsäcker potentials together with well-known functions, e. g., deformation electron density and electron localization function. In this way, our analysis throws light on aspects of these closed-shell interactions hitherto hidden from the description.

Multi-stage screening to predict the specific anticancer activity of Ni(II) mixed-ligand complex on gastric cancer cells; biological activity, FTIR spectrum, DNA binding behavior and simulation studies

The anticancer activity of a transition metal complex with [Ni(L¹)₂L²]H₂O (where L¹ and L² were acetylacetonato (acac) and 2-aminopyridine (2-ampy), respectively) was evaluated in MKN45 cell line. Methyl thiazolyl tetrazolium (MTT) assay was performed to assess the antitumor capacity of the Ni(II) complex against gastric cancer cell line MKN45. The complexexhibited high in vitro antitumor activity against MKN45 cells with IC₅₀values of 1.99 μM in 48 hrs. The alterations in the structure of cellular biomolecules (proteins, lipids, carbohydrates, and especially DNA) by the Ni(II) complex were confirmed by bio spectroscopic studies. Fourier Transformed Infrared (FTIR) spectroscopy analysis revealed significant differences between untreated and treated MKN45 cell line in the region of glycogen, nucleic acid, amide I and amide II bands (1000, 1100, ~1650, and ~1577 cm^-1). The absorption bands 1150 cm^-1 and 1020-1025 cm^-1 can be assigned to the CO bond of glycogen and other carbohydrates and are significantly overlapped by DNA. The interaction of calf thymus (CT) DNA with Ni(II) complex was explored using absorption spectral method. The UV-visible studies demonstrated that this complex was able to bind with DNA via groove, non-covalent, and electrostatic interactions, and binding constant (K_b) was found to be 3 * 10⁴. Docking simulation and Non Covalent Interaction (NCI) topological analysis were conducted to provide insights into the nature of DNA/complex interactions. The binding affinity and binding stability of complex was validated by 400-ns MD simulations.

Strength of the [Z-I···Hal]- and [Z-Hal···I]- Halogen Bonds: Electron Density Properties and Halogen Bond Length as Estimators of Interaction Energy

Bond energy is the main characteristic of chemical bonds in general and of non-covalent interactions in particular. Simple methods of express estimates of the interaction energy, Eint, using relationships between Eint and a property which is easily accessible from experiment is of great importance for the characterization of non-covalent interactions. In this work, practically important relationships between Eint and electron density, its Laplacian, curvature, potential, kinetic, and total energy densities at the bond critical point as well as bond length were derived for the structures of the [Z-I···Hal]- and [Z-Hal···I]- types bearing halogen bonds and involving iodine as interacting atom(s) (totally 412 structures). The mean absolute deviations for the correlations found were 2.06-4.76 kcal/mol.

Synthesis and characterization of N-heterocyclic carbene-MOEt2 complexes (M = Cu, Ag, Au). Analysis of solvated auxiliary-ligand free [(NHC)M]+ species

We report the synthesis, characterization and computational analysis of coinage metal-ether complexes supported by N-heterocyclic carbenes (NHC), SIPr and Et2CAAC. The related water adducts are also included. The [(NHC)M]+(M = Cu, Ag, Au) species show the noteworthy ability to bind Et2O and H2O. This interaction towards Et2O and H2O is partly ascribed to a σ-hole bonding with an almost linear disposition, taking advantage of the enhanced σ-hole potential evaluated for such [(NHC)M]+ species. This enhanced ability is larger than those found for non-covalent interactions involving main group species.

A Step toward the Quantification of Noncovalent Interactions in Large Biological Systems: The Independent Gradient Model-Extremely Localized Molecular Orbital Approach

The independent gradient model (IGM) is a recent electron density-based computational method that enables to detect and quantify covalent and noncovalent interactions. When applied to large systems, the original version of the technique still relies on promolecular electron densities given by the sum of spherically averaged atomic electron distributions, which leads to approximate evaluations of the inter- and intramolecular interactions occurring in systems of biological interest. To overcome this drawback and perform IGM analyses based on quantum mechanically rigorous electron densities also for macromolecular systems, we coupled the IGM approach with the recently constructed libraries of extremely localized molecular orbitals (ELMOs) that allow fast and reliable reconstructions of polypeptide and protein electron densities. The validation tests performed on small polypeptides and peptide dimers have shown that the novel IGM-ELMO strategy provides results that are systematically closer to the fully quantum mechanical ones and outperforms the IGM method based on the crude promolecular approximation, but still keeping a quite low computational cost. The results of the test calculations carried out on proteins have also confirmed the trends observed for the IGM analyses conducted on small systems. This makes us envisage the future application of the novel IGM-ELMO approach to unravel complicated noncovalent interaction networks (e.g., in protein-protein contacts) or to rationally design new drugs through molecular docking calculations and virtual high-throughput screenings.

Coinage-metal pillarplexes hosts. Insights into host-guest interaction nature and luminescence quenching effects

Host-guest chemistry is a relevant issue in materials science, which encourages further development of versatile host structures. Here the particular features of coinage-metal pillarplexes are evaluated towards formation of host-guest aggregates by the inclusion of 1,8-diaminooctane, as characterized for [M8(LMe)2]4+ (M = Ag, and, Au). The obtained results denotes the main contribution from van der Waals type interaction (50%), followed by a contribution from orbital polarization and electrostatic nature (20% and 30%), involving both orbitalary and electrostatic terms. Throughout the different coinage-metal based hosts (M = Cu, Ag, and Au), a similar interaction energy is found given by the large contribution of the π-surface from the organic ligand backbone to both van de Waals and electrostatic interactions. This suggests that a similar host structure can be obtained for the lighter copper counterpart, retaining similar how-guest features. Moreoves, the [Au8(LMe)2]4+ host exhibits inherent luminescent properties, involving the shortening of Au(i)-Au(i) contacts at the excited state, which is partially avoided when the guest is incorporated, accounting for the observed quenching from titration experiments. This results encourages further exploration of coinage metal hosts in the formation of inclusion complexes.

Exotic Two-Dimensional Structure: The First Case of Hexagonal NaCl

NaCl is one of the simplest compounds and was thought to be well-understood, and yet, unexpected complexities related to it were uncovered at high pressure and in low-dimensional states. Here, exotic hexagonal NaCl thin films on the (110) diamond surface were crystallized in the experiment following a theoretical prediction based on ab initio evolutionary algorithm USPEX. State-of-the-art calculations and experiments showed the existence of a hexagonal NaCl thin film, which is due to the strong chemical interaction of the NaCl film with the diamond substrate.

The Volumetric Source Function: Looking Inside van der Waals Interactions

The study of van der Waals interactions plays a central role in the understanding of bonding across a range of biological, chemical and physical phenomena. The presence of van der Waals interactions can be identified through analysis of the reduced density gradient, a fundamental parameter at the core of Density Functional Theory. An extension of Bader’s Quantum Theory of Atoms in Molecules is developed here through combination with the analysis of the reduced density gradient. Through this development, a new quantum chemical topological tool is presented: the volumetric source function. This technique allows insight into the atomic composition of van der Waals interactions, offering the first route towards applying the highly successful source function to these disperse interactions. A new algorithm has been implemented in the open-source code, CRITIC2, and tested on acetone, adipic and maleic acids molecular crystals, each stabilized by van der Waals interactions. This novel technique for studying van der Waals interactions at an atomic level offers unprecedented opportunities in the fundamental study of intermolecular interactions and molecular design for crystal engineering, drug design and bio-macromolecular processes.