DevelopingWinXPRO: a software for determination of the multipole-model-based properties of crystals

Noncovalent Interactions of Silatranes and Germatranes with the Surface of Stoichiometric and Hydroxylated 2D Silica

Analysis of noncovalent interactions formed by the surface of a 2D silica bilayer and atrane molecules, as precursors of functional layers immobilized on a surface of silicatene, was performed. For this purpose, the systems of substituted silatranes and germatranes adsorbed on silicatene surfaces with different amounts of hydroxyl groups (0, 2, 4, and 30 per cell) were simulated by using quantum chemical modeling with periodic boundary conditions and full-electron basis sets. The observation results for interaction energy changes in the systems "atrane molecule-silicatene surface" with increasing silanol number of the silicatene surface can be used to predict the optimal degree of silicatene hydroxylation in order to control the effective progress of atrane deprotection on activated 2D silica materials. In addition to the typical hydrogen bonds, the ability of atranes to form noncovalent O···O and O···Hal interactions was discovered. In these bonds, the oxygen or halogen atoms of atranes act as electron-donor centers in relation to the silicatene oxygen atoms. The observed weakening of the Ge-O covalent bonds in germatranes, on which further deprotection reaction depends, is more manifested than for the Si-O bonds in adsorbed silatranes.

Valence tautomeric interconversion of bis-dioxolene cobalt complex with imino-pyridine functionalized by TEMPO moiety in solid solutions with isostructural nickel analogue: phase transitions and monocrystal destruction

Valence tautomeric complexes (VT) are promising systems for creating molecular devices. From this viewpoint, valence tautomeric complexes with a hysteresis loop on the magnetic curve are of special interest as potential memory elements. The hysteresis loop is a consequence of retarded structural rearrangements which investigation is an actual problem. Recently, we have described a new VT transition taking place in a bis-dioxolene cobalt complex with imino-pyridine having a TEMPO substituent (A. A. Zolotukhin, et al., Inorg. Chem., 2017, 56, 14751-14754). Valence tautomeric transformation occurs with a hysteresis loop and is accompanied by a phase transition. The phase transition taking place during cooling is accompanied by crystal destruction. This fact makes it impossible to monitor the structural changes responsible for the hysteresis loop. The current research attempts to resolve this problem. A nickel compound of the same composition (TEMPO-imino-pyridine)Ni(3,6-DBSQ)2 was synthesized and characterized. It was established to be isostructural with the cobalt complex. It was used as an inert matrix for the dilution of the VT cobalt complex. The number of solid solutions with Co/Ni ratios of 1 : 1, 1 : 2, 1 : 4, and 1 : 8 was obtained. Variable temperature magnetic susceptibility measurements show that VT transformation with a hysteresis loop takes place in all solid solutions. The hysteresis loop is shifted to low temperatures primarily due to the shifting of its low-temperature boundary with dilution. The hysteresis width does not change significantly with dilution. DSC detected that transformations are accompanied by phase transitions at different temperatures at cooling and heating. The phase transition at the first cooling occurs at slightly lower temperatures compared with subsequent cycles. These temperatures correspond to the transition temperatures detected from the magnetic curves. The phase transition during the first cooling is accompanied by crystal destruction. Physical destruction takes place in the crystals of all solid solutions. X-ray diffraction powder patterns confirm that phase transition is accompanied by considerable reorganization of the crystal structure typical for the first order transitions. The unit cell volume of solid solutions is larger than that of pure complexes. Especially calculated crystal invariom indicated that the "lattice energy" in a solid solution is the lowest compared with that in "pure" nickel and cobalt complexes.

Insight into crystal structures and identification of potential styrylthieno[2,3-b]pyridine-2-carboxamidederivatives against COVID-19 Mpro through structure-guided modeling and simulation approach

Anti-SARS-CoV-2 drugs are urgently needed to prevent the pandemic and for immunization. Their protease inhibitor treatment for COVID-19 has been used in clinical trials. In Calu-3 and THP1 cells, 3CL SARS-CoV-2 Mpro protease is required for viral expression, replication, and the activation of the cytokines IL-1, IL-6, and TNF-. The Mpro structure was chosen for this investigation because of its activity as a chymotrypsin-like enzyme and the presence of a cysteine-containing catalytic domain. Thienopyridine derivatives increase the release of nitric oxide from coronary endothelial cells, which is an important cell signaling molecule with antibacterial activity against bacteria, protozoa, and some viruses. Using DFT calculations, global descriptors are computed from HOMO-LUMO orbitals; the molecular reactivity sites are analyzed from an electrostatic potential map. NLO properties are calculated, and topological analysis is also part of the QTAIM studies. Both compounds 1 and 2 were designed from the precursor molecule pyrimidine and exhibited binding energies (-14.6708 kcal/mol and -16.4521 kcal/mol). The binding mechanisms of molecule 1 towards SARS-COV-2 3CL Mpro exhibited strong hydrogen bonding as well as Vdw interaction. In contrast, derivative 2 was bound to the active site protein's active studied that several residues and positions, including (His41, Cys44, Asp48, Met49, Pro52, Tyr54, Phe140, Leu141, Ser144, His163, Ser144, Cys145, His164, Met165, Glu166, Leu167, Asp187, Gln189, Thr190, and GLn192) are critical for the maintenance of inhibitors inside the active pocket. Molecular docking and 100 ns MD simulation analysis revealed that Both compounds 1 and 2 with higher binding affinity and stability toward the SARS-COV-2 3CL Mpro protein. Binding free energy calculations and other MD parameters support the finding.Communicated by Ramaswamy H. Sarma.

Metamagnetism and canted antiferromagnetic ordering in two monomeric CoII complexes with 1-(2-pyrimidyl)piperazine. Hirshfeld surface analysis and theoretical studies

Here we present the magnetic properties of two cobalt complexes formulated as: [Co(SCN)2(L)2] (1) and (H2L)2[Co(SCN)4]·H2O (2) (L = 1-(2-pyrimidyl)piperazine). The two compounds contain isolated tetrahedral CoII complexes with important intermolecular interactions that lead to the presence of a canted antiferromagnetic order below 11.5 and 10.0 K, with coercive fields at 2 K of 38 and 68 mT, respectively. Theoretical calculations have been used to explain this behaviour. Hirshfeld surface analysis shows the presence of strong intermolecular interactions in both compounds. The crystal geometries were used for geometry optimization using the DFT method. From the topological properties, electrostatic potential maps and molecular orbital analysis, information about the noncovalent interaction and chemical reactivity was obtained.

A biophysical approach of tyrphostin AG879 binding information in: bovine serum albumin, human ErbB2, c-RAF1 kinase, SARS-CoV-2 main protease and angiotensin-converting enzyme 2

Viral infections cause significant health problems all over the world, and it is critical to develop treatments for these problems. Antivirals that target viral genome-encoded proteins frequently cause the virus to become more resistant to treatment. Because viruses rely on several cellular proteins and phosphorylation processes that are essential to their life cycle, drugs targeting host-based targets could be a viable treatment option. To reduce costs and improve efficiency, existing kinase inhibitors could be repurposed as antiviral medications; however, this method rarely works, and specific biophysical approaches are required in the field. Because of the widespread use of FDA-approved kinase inhibitors, it is now possible to better understand how host kinases contribute to viral infection. The purpose of this article is to investigate the tyrphostin AG879 (Tyrosine kinase inhibitor) binding information in Bovine Serum Albumin (BSA), human ErbB2 (HER2), C-RAF1 Kinase (c-RAF), SARS-CoV-2 main protease (COVID 19), and Angiotensin-converting enzyme 2 (ACE-2).Communicated by Ramaswamy H. Sarma.

A Structural and In Silico Investigation of Potential CDC7 Kinase Enzyme Inhibitors

A crucial role in the regulation of DNA replication is played by the highly conserved CDC kinase. The CDC7 kinase could serve as a target for therapeutic intervention in cancer. The primary heterocyclic substance is pyrazole, and its derivatives offer great potential as treatments for cancer cell lines. Here, we synthesized the two pyrazole derivatives: 4-(2-(4-chlorophenyl)hydrazinyl)-5-methyl-2-tosyl-1H-pyrazol-3(2H)-one (PYRA-1) and 4-(2-(2,4-difluorophenyl)hydrazinyl)-5-methyl-2-tosyl-1H-pyrazol-3(2H)-one (PYRA-2). The structural confirmation of both the compounds at the three-dimensional level is characterized using single crystal X-ray diffraction and density functional theory. Furthermore, the in silico chemical biological properties were derived using molecular docking and molecular dynamics (MD) simulations. PYRA-1 and PYRA-2 crystallize in the P-1 (a = 8.184(9), b = 14.251(13), c = 15.601(15), α = 91.57(8), β = 97.48(9), 92.67(9), V = 1801.1(3) 3, and Z = 2) and P21/n (a = 14.8648(8), b = 8.5998(4), c = 15.5586(8), β = 116.47(7), V = 1780.4(19) 3, and Z = 4), space groups, respectively. In both PYRA-1 and PYRA-2 compounds, C-H···O intermolecular connections are common to stabilize the crystal structure. In addition, short intermolecular interactions stabilizes with C-H···π and π-π stacking. Crystal packing analysis was quantified using Hirshfeld surface analysis resulting in C···H, O···H, and H···H contacts in PYRA-1 exhibiting more contribution than in PYRA-2. The conformational stabilities of the molecules are same in the gas and liquid phases (water and DMSO). The docking scores measured for PYRA-1 and PYRA-2 with CDC7 kinase complexes are -5.421 and -5.884 kcal/mol, respectively. The MD simulations show that PYRA-2 is a more potential inhibitor than PYRA-1 against CDC7 kinase.

Computational Binding Analysis of Ethyl 3,3,5,5-Tetracyano-2-Hydroxy-2-Methyl-4,6-Diphenylcyclohexane-1-Carboxylate in Calf Thymus DNA

In the present paper, several computational binding analyses were performed on ethyl 3,3,5,5-tetracyano-2-hydroxy-2-methyl-4,6-diphenylcyclohexane-1-carboxylate which was newly synthesized by three-component condensation of benzaldehyde with ethyl acetoacetate and malononitrile in the presence of trichloroacetic acid, and the structure was finally proved by X-ray analysis. The visualization of molecular interaction was carried out through Hirshfeld surface analysis and ESP. The atomic charges, HOMO, LUMO, and electrostatic potential were also studied to explore the insight of the molecule deeper, and then, natural bonding orbitals (NBO) and non-linear optical properties (NLO) were calculated to reveal the interactions that happen to be between the filled and vacant orbitals. Afterwards, molecular docking studies predicted the compound binding mode fits in the minor groove of DNA and remained interacts via stable bonding as validated by molecular dynamics simulations. The binding energy estimation also affirmed domination van der Waals and electrostatic energies. Lastly, the compound was found as good drug-like molecule and had good pharmacokinetic profile with exception of toxic moieties.

Applicability of transferable multipole pseudo-atoms for restoring inner-crystal electronic force density fields. Chemical bonding and binding features in the crystal and dimer of 1,3-bis(2-hydroxyethyl)-6-methyluracil

We considered it timely to test the applicability of transferable multipole pseudo-atoms for restoring inner-crystal electronic force density fields. The procedure was carried out on the crystal of 1,3-bis(2-hydroxyethyl)-6-methyluracil, and some derived properties of the scalar potential and vector force fields were compared with those obtained from the experimental multipole model and from the aspherical pseudo-atom model with parameters fitted to the calculated structure factors. The procedure was shown to accurately replicate the general vector-field behavior, the peculiarities of the quantum potentials and the characteristics of the force-field pseudoatoms, such as charge, shape and volume, as well as to reproduce the relative arrangement of atomic and pseudoatomic zero-flux surfaces along internuclear regions. It was found that, in addition to the quantum-topological atoms, the force-field pseudoatoms are spatially reproduced within a single structural fragment and similar environment. In addition, the classical and nonclassical hydrogen bonds in the uracil derivative crystal, as well as the H...O, N...O and N...C interactions in the free π-stacked dimer of the uracil derivative molecules, were studied using the potential and force fields within the concepts of interatomic charge transfer and electron lone pair donation-acceptance. Remarkably, the nitrogen atoms in the N...O and N...C interactions behave rather like a Lewis base and an electron contributor. At the same time, the hydrogen atom in the H...O interaction, being a Lewis acid, also participates in the interatomic electron transfer by acting as a contributor. Thus, it has been argued that, when describing polar interatomic interactions within orbital-free considerations, it makes more physical sense to identify electronegative (electron occupier) and electropositive (electron contributor) atoms or subatomic fragments rather than nucleophilic and electrophilic sites.

1,4,9,9-tetramethyloctahydro-4,7-(epoxymethano)azulen-5(1H)-one, a natural product as a potential inhibitor of COVID-19: Extraction, crystal structure, and virtual screening approach

In the present work, we describe the extraction of a natural product namely 1,4,9,9-tetramethyloctahydro-4,7-(epoxymethano)azulen-5(1H)-one, and its structure was confirmed by single crystal X-ray diffraction analysis. The conformations of the 5-, 6-, and 7-membered rings in the title compound, C₁₅H₂₄O₂, have been probed by a Cremer-Pople puckering analysis. C-H···O hydrogen bonds generate chains in the crystal that stretch along the c-axis direction. The Hirshfeld surface analysis method was used to stabilize the crystal packing of the natural compound. Accompanied by experimental studies, quantum chemical calculations were also performed to compare the structural elucidation and the results of these geometrical parameters exhibited excellent agreement. The compound was also docked with several drug targets of the SARS-CoV-2 virus and found to show the best binding with the main protease enzyme, having a binding energy of -12.31 kcal/mol and interacting with His41 and Cys145 residues. The dynamic stability deciphered the complex to be stable with an average RMSD of 3.8 Å. The compound dynamics with the enzyme showed the compound conformation to be highly stable. The intermolecular binding free energy determined the compound-main protease enzyme to show high interaction energy of < 40 kcal/mol. Together, these studies demonstrate the compound to be a lead structure against SARS-CoV-2.

Probing the intermolecular interactions, binding affinity, charge density distribution and dynamics of silibinin in dual targets AChE and BACE1: QTAIM and molecular dynamics perspective

Alzheimer's disease (AD) is the grievous neurodegenerative disorder. Reportedly, many enzymes are responsible for this disease, in which notably, acetylcholinesterase (AChE) and β-secretase (BACE1) are largely involved for AD. An experimental study reports that silibinin molecule inhibits both AChE and BACE1 enzymes. Present study aims to understand the dual binding mechanism of silibinin in the active site of AChE and BACE1 from the intermolecular interactions, conformational flexibility, charge density distribution, binding energy and the stability of molecule. To obtain the above information, the molecular docking, molecular dynamics (MD) and QTAIM (quantum theory of atoms in molecules) calculations have been performed. The molecular docking reveals that silibinin molecule is forming strong and weak intermolecular interactions with the catalytic site of both enzymes. The QTAIM analysis for the binding pockets of both complexes shows the charge density distribution of intermolecular interactions. The electrostatic potential map displays the electronegative/positive regions at the interaction zone of silibinin with AChE and BACE1 complexes. The MD simulation confirms that the silibinin molecule is stable in the active site of AChE and BACE1 enzymes. The binding free energies of silibinin with both enzymes are more favorable to have the interactions.Communicated by Ramaswamy H. Sarma.

Charge density studies of single and transient (single to double) boron-oxygen bonds in (NH4)2B4O5(OH)4·2H2O

A H₄B₄O₉^2- ion which makes up the (NH₄)₂B₄O₅(OH)₄·2H₂O crystal structure has two types of boron-oxygen bonds, i.e. single B-O bonds and an intermediate between single and double BO bonds. Differences between these two bond types are visible not only because they differ by their lengths but also a topology of electron density distribution differs. This also gives a hint as to how to distinguish between these two bond types. Experimental results based on multipole model refinement gave excellent agreement with theoretical calculations and literature data. Calculations at bond critical points for B-O and BO (electron density, the Laplacian of electron density and the localized-orbital locator function) suggest us how boron-oxygen bonds should be categorised with respect to compounds previously reported in the literature. Additionally, a novel synthesis method for the investigated compound has been developed, which involves crystallization from an aqueous solution of BH₃NH₃ dissolved in a mixture of tetrahydrofuran and water.

Self-assembly of new cobalt complexes based on [Co (SCN)4], synthesis, empirical, antioxidant activity, and quantum theory investigations

The cobalt (II) complexes have been synthesized from the reaction of the cationic entities (3,4-dimethylaniline (1) and histamine (2)) with metallic salt CoCl₂⋅6H₂O and thiocyanate ion (SCN⁻) as a ligand in H₂O/ethanolic solution and processing by the evaporation crystal growth method at room temperature to get crystals. The synthesized complex has been fully characterized by single-crystal X-ray diffraction. UV–Visible, FTIR spectroscopy, TGA analysis, and DFT circulations were also performed. The crystal structural analysis reveals that the solid (1) {[Co(SCN)₄] (C₈H₁₂N)₃}·Cl crystallizes in the monoclinic system with the space group P2₁/n and the solid (2) {[Co(SCN)₄](C₅H₁₁N₃)₂}·2Cl crystallizes in the monoclinic space group P2₁/m. Metal cations are joined into corrugated chains parallel to the b-axis direction in (1) and (2) by four thiocyanate anions. The crystal structures of (1) and (2) were calculated using XRPD data, indicating that they are closely connected to the DRX mono-crystal results. Different interactions pack the system into a ring formed by N–H⋯Cl and N–H⋯S hydrogen bonds. C–H⋯π and the π⋯π stacking of anilinuim ring for (1) and N–H⋯S intermolecular interactions for (1) and (2) increase the crystals' robustness. Hirshfeld surface analysis cum 2D fingerprint plots visualize the main intermolecular interactions with their contributions in the solid-state phase. The molecular geometries of both complexes obtained from the crystal structure were used for quantum chemical calculation. Here, frontier orbital analysis and electrostatic potential illustrate the chemical reactivities of metal–organic complexes. QTAIM and NCI analysis reveal the strength of interactions at the electronic level.

Synthesis, virtual screening and computational approach of a quinoxaline derivative as potent anti-HIV agent targeting the reverse transcriptase enzyme

Infection by the human immunodeficiency virus still represents a continuous serious concern and a global threat to human health. Due to the appearance of multi-resistant virus strains and the serious adverse side effects of the antiretroviral therapy administered, there is an urgent need for the development of new treatment agents that are more active, less toxic, and with increased tolerability to mutations. Quinoxaline derivatives are a class of heterocyclic compounds with a wide range of organic and remedial applications. In addition, they are known to significantly inhibit HIV reverse transcriptase (RT) and HIV replication in cell cultures. For these reasons, we are investigating the synthesis and computational studies of quinoxaline derivatives with a focus on their effects on the HIV RT enzyme, and we present here the structure of one such molecule, methyl 2-[(2E)-3-oxo-1,2,3,4-tetrahydroquinoalin-2-ylidene] acetate, which was confirmed by X-ray diffraction studies. In the crystal, N-H···O and C-H···O hydrogen bonds form ribbons whose mean planes are inclined to (111) by 25.69(8)°. The ribbons are formed into stacks by C-H···π(ring) interactions and π-stacking interactions between carbonyl groups. The Hirshfeld surface map allows us to understand the nature of interactions in the contribution to crystal packing. A density functional theory (DFT) calculation was performed to optimize the geometrical parameters and then they were compared with the solid-state phase. The molecular electrostatic potential map displays reactive sites on the surface, which are responsible for intermolecular interaction in the chemical species. Computational molecular docking, in addition to molecular dynamics simulations and MMGB/PBSA binding energy techniques, was used to assess the affinity of the molecule for the HIV reverse transcriptase enzyme. The new quinoxaline derivative is more powerful in terms of binding affinity and binding conformation stability with the HIV reverse transcriptase enzyme, which suggests the molecule is a good candidate for further biological optimization.Communicated by Ramaswamy H. Sarma.

WinXPRO, 3DPlot and TrajPlot computer software: new options for orbital-free quantum crystallography studies

A multipurpose computer software package for orbital-free quantum crystallography has been developed. Based on multipole experimental electron-density parameters, this software allows users to extract the chemical-bonding information that was previously unavailable for X-ray single-crystal diffraction analysis.

Thermal-induced transformation of glutamic acid to pyroglutamic acid and self-cocrystallization: a charge–density analysis

Thermal-induced transformation of glutamic acid to pyroglutamic acid is well known. However, confusion remains over the exact temperature at which this happens. Moreover, no diffraction data are available to support the transition. In this article, we make a systematic investigation involving thermal analysis, hot-stage microscopy and single-crystal X-ray diffraction to study a one-pot thermal transition of glutamic acid to pyroglutamic acid and subsequent self-cocrystallization between the product (hydrated pyroglutamic acid) and the unreacted precursor (glutamic acid). The melt upon cooling gave a robust cocrystal, namely, glutamic acid–pyroglutamic acid–water (1/1/1), C5H7NO3·C5H9NO4·H2O, whose structure has been elucidated from single-crystal X-ray diffraction data collected at room temperature. A three-dimensional network of strong hydrogen bonds has been found. A Hirshfeld surface analysis was carried out to make a quantitative estimation of the intermolecular interactions. In order to gain insight into the strength and stability of the cocrystal, the transferability principle was utilized to make a topological analysis and to study the electron-density-derived properties. The transferred model has been found to be superior to the classical independent atom model (IAM). The experimental results have been compared with results from a multipolar refinement carried out using theoretical structure factors generated from density functional theory (DFT) calculations. Very strong classical hydrogen bonds drive the cocrystallization and lend stability to the resulting cocrystal. Important conclusions have been drawn about this transition.

Crystal engineering of a co-crystal of antipyrine and 2-chlorobenzoic acid: relative energetic contributions based on multipolar refinement

The growth and stability of a new 1 : 1 antipyrene–dichlorobenzoic acid cocrystal system has been analyzed in terms of electron density analysis and electrostatic interaction energy contributions.

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.

Orbital-free quantum crystallography: view on forces in crystals

Quantum theory of atoms in molecules and the orbital-free density functional theory (DFT) are combined in this work to study the spatial distribution of electrostatic and quantum electronic forces acting in stable crystals. The electron distribution is determined by electrostatic electron mutual repulsion corrected for exchange and correlation, their attraction to nuclei and by electron kinetic energy. The latter defines the spread of permissible variations in the electron momentum resulting from the de Broglie relationship and uncertainty principle, as far as the limitations of Pauli principle and the presence of atomic nuclei and other electrons allow. All forces are expressed via kinetic and DFT potentials and then defined in terms of the experimental electron density and its derivatives; hence, this approach may be considered as orbital-free quantum crystallography. The net force acting on an electron in a crystal at equilibrium is zero everywhere, presenting a balance of the kinetic F kin( r ) and potential forces F ( r ). The critical points of both potentials are analyzed and they are recognized as the points at which forces F kin( r ) and F ( r ) individually are zero (the Lagrange points). The positions of these points in a crystal are described according to Wyckoff notations, while their types depend on the considered scalar field. It was found that F ( r ) force pushes electrons to the atomic nuclei, while the kinetic force F kin( r ) draws electrons from nuclei. This favors formation of electron concentration bridges between some of the nearest atoms. However, in a crystal at equilibrium, only kinetic potential v kin( r ) and corresponding force exhibit the electronic shells and atomic-like zero-flux basins around the nuclear attractors. The force-field approach and quantum topological theory of atoms in molecules are compared and their distinctions are clarified.

Variations of quantum electronic pressure under the external compression in crystals with halogen bonds assembled in Cl3-, Br3-, I3-synthons

The inner-crystal quantum electronic pressure was estimated for unstrained C6Cl6, C6Br6, and C6I6 crystals and for those under external compression simulated from 1 to 20 GPa. The changes in its distribution were analyzed for the main structural elements in considered crystals: for triangles of the typical halogen bonds assembled in Hal3-synthons, where Hal = Cl, Br, I; for Hal...Hal stacking interactions, as well as for covalent bonds. Under simulated external compression, the quantum electronic pressure in the intermolecular space reduces as the electron density increases, indicating spatial areas of relatively less crystal resistance to external compression. The most compliant C6Cl6 crystal shows the largest changes of quantum electronic pressure in the centre of Cl3-synthon while the deformation of rigid I3-synthon under external compression depends only on the features of I...I halogen bonds.

Exploring the different environments effect of piperine via combined crystallographic, QM/MM and molecular dynamics simulation study

Piperine is a pungent alkaloid, largely present in the skin of pepper. It is the most active component of pepper and being used as a medicine in many Asian countries. The effect of piperine on memory impairment and neurodegeneration in Alzheimer's disease model has been investigated. In the present study, we aim to investigate the effect of piperine molecule in different environments (crystal and active site of proteins) from crystallography, molecular docking, QM/MM based charge density analysis and molecular dynamic simulation. The crystal structure of piperine has been used to determine the topological electron density of intermolecular interactions. The O-atoms of piperine is forming C-H⋅⋅⋅O interactions with the neighboring molecules in the crystal, these interactions also confirmed from the Hirshfeld surface. Further, to understand the nature of interactions and the conformational flexibility of piperine in the active site of recombinant human acetylcholinesterase (rhAChE), molecular docking analysis has been performed. The selected docked complex suggests favorable hydrogen bonding and hydrophobic interactions with rhAChE enzyme; notably, the O3 atom of piperine molecule forms strong hydrogen bonding interaction with Glu202 at 1.8 Å. To determine the charge density distribution and the electrostatic properties of piperine molecule in the active site of rhAChE, the piperine-rhAChE complex was minimized at QM/MM energy level; in which, the binding pocket with piperine was considered as QM region. The charge density analysis of piperine and the interacting amino acid groups have been carried out. The topological analysis of O3⋯H-O/Glu202 hydrogen bonding interaction exhibits strong interactions and the electron density ρ_cp(r): 0.242 eÅ^-3 and the Laplacian ∇²ρ_cp(r): 3.176 eÅ^-5 respectively. These results were compared with the corresponding molecule present in the crystal and gas phase environments of piperine. The comparison of active site structure with the corresponding crystal phase and gas phase structures reveal that piperine exhibits large conformational modification in the active site. The molecular dynamics simulation and binding free energy calculations were performed, this gives the stability, binding affinity of the molecule in the active site of rhAChE. The O3⋯H-O/Glu202 interaction shows the high stability (89.2%), this was confirmed from the stability of hydrogen bond analysis. The binding free energy was used to measure the rate of inhibition of enzyme in the presence of ligand molecule. The comparative study allows to understand the nature of piperine molecule in the gas and crystal phases, and amino acids environment.

Experimental and theoretical charge-density analysis of hippuric acid: insight into its binding with human serum albumin

In order to comprehend the binding of an important metabolite, hippuric acid, with human serum albumin and to understand its chemical and electronic nature, an experimental charge-density analysis has been carried out using high-resolution diffraction data collected under cryogenic conditions, and all the results have been compared with theoretical findings using the B3LYP/6-311++g(2d,2p) level of theory. The structure displays very strong classical hydrogen bonds as well as other noncovalent interactions, which have been fully characterized using Hirshfeld surface analysis and Bader's quantum theory of atoms in molecules. Contact analysis on the Hirshfeld surfaces shows that the O...H, C...H and C...N intermolecular interactions are enriched and gives their relative strengths. Topological analysis of the electron density shows the charge concentration/depletion of hippuric acid bonds in the crystal structure. Electrostatic parameters such as atomic charges and dipole moments were calculated. The mapping of atomic basins and the calculation of respective charges show the atomic volumes of each atom as well as their charge contributions in the hippuric acid crystal structure. The dipole-moment calculations show that the molecule is very polar in nature. Calculations of the electrostatic potential show that the chain part of the molecule has a higher concentration of negative charge than the ring, which might be instrumental in its strong binding with the polar residues of site II of human serum albumin.

Quantum pressure focusing in solids: a reconstruction from experimental electron density

Here an approach is presented for reconstructing the distribution of electronic internal quantum pressure in the electronic continuum of solids from the experimental electron density. Using the formalism of the density functional theory, the spatial inner-crystal map of the quantum pressure is obtained. The results are visualized via the indicator of quantum pressure focusing (IQPF) which reveals the regions where the pressure is concentrated or depleted due to quantum effects. IQPF contains all quantum electron-shell structure-forming contributions resulting from kinetic, exchange and correlation effects, and presents a clear picture of the chemical bond features in crystals with different type of bonding mechanisms.

Electronic criterion for categorizing the chalcogen and halogen bonds: sulfur–iodine interactions in crystals

Diversity of mutual orientations of Y–S and I–X and covalent bonds in molecular crystals complicate categorizing noncovalent chalcogen and halogen bonds. Here, the different types of S...I interactions with short interatomic distances are analysed. The selection of S...I interactions for the categorization of the chalcogen and halogen bonds has been made using angles that determine the mutual orientation of electron lone pairs and σ-holes interacted S and I atoms. In complicated cases of noncovalent interactions with `hole-to-hole' of S and I orientations, distinguishing the chalcogen and halogen bonds is only possible if the atom is uniquely determined, which also provides the electrophilic site. The electronic criterion for chalcogen/halogen bonds categorizing that is based on analysis of dispositions of electron density and electrostatic potential minima along the interatomic lines has been suggested and its effectiveness has been demonstrated.

Aspherical scattering factors for SHELXL - model, implementation and application

A new aspherical scattering factor formalism has been implemented in the crystallographic least-squares refinement program SHELXL. The formalism relies on Gaussian functions and can optionally complement the independent atom model to take into account the deformation of electron-density distribution due to chemical bonding and lone pairs. Asphericity contributions were derived from the electron density obtained from quantum-chemical density functional theory computations of suitable model compounds that contain particular chemical environments, as defined by the invariom formalism. Thanks to a new algorithm, invariom assignment for refinement in SHELXL is automated. A suitable parameterization for each chemical environment within the new model was achieved by metaheuristics. Figures of merit, precision and accuracy of crystallographic least-squares refinements improve significantly upon using the new model.

Ground-State Charge-Density Distribution in a Crystal of the Luminescent ortho-Phenylenediboronic Acid Complex with 8-Hydroxyquinoline

This contribution is devoted to the first electron density studies of a luminescent oxyquinolinato boron complex in the solid state. ortho-Phenylenediboronic acid mixed with 8-hydroxyquinoline in dioxane forms high-quality single crystals via slow solvent evaporation, which allows successful high resolution data collection (sin θ/λ = 1.2 Å^-1) and charge density distribution modeling. Particular attention has been paid to the boron-oxygen fragment connecting the two parts of the complex, and to the solvent species exhibiting anharmonic thermal motion. The experiment and theory compared rather well in terms of atomic charges and volumes, except for the boron centers. Boron atoms, as expected, constitute the most electron-deficient species in the complex molecule, whereas the neighboring oxygen and carbon atoms are the most significantly negatively charged ones. This part of the molecule appears to be very much involved in the charge transfer occurring between the acid fragment and oxyquinoline moiety leading to the observed fluorescence, as supported by the time-dependent density functional theory (TDDFT) results and the generated transition density maps. TDDFT calculations indicated that p-type atomic orbitals contributing to the HOMO-1, HOMO, and LUMO play the major role in the lowest energy transitions, and enabled further comparison with the charge density features, which is discussed in details. Furthermore, the results confirmed the known fact the Q ligand character is most important for the spectroscopic properties of this class of complexes.

Testing the tools for revealing and characterizing the iodine-iodine halogen bond in crystals

To understand what tools are really suitable to identify and classify the iodine-iodine non-covalent interactions in solid organic polyiodides, we have examined the anisotropy of the electron density within the iodine atomic basin along and across the iodine-iodine halogen bond using the Laplacian of electron density, one-electron potential and electron localization function produced by Kohn-Sham calculations with periodic boundary conditions. The Laplacian of electron density exhibits the smallest anisotropy and yields a vague picture of the outermost electronic shells. The one-electron potential does not show such a deficiency and reveals that the valence electron shell for the halogen-bond acceptor iodine is always wider than that for the halogen-bond donor iodine along its σ-hole direction. We have concluded that the one-electron potential is the most suitable for classification of the iodine-iodine bonds and interactions in complicated cases, while the electron localization function allows to distinguish the diiodine molecule bonded with the monoiodide anion from the typical triiodide anion.

Theoretical study of intermolecular interactions in crystalline arene–perhaloarene adducts in terms of the electron density

The intermolecular interactions in C₆X₆–arene crystals (X = F, Cl) and the halogen substitution effect can be quantified by the electron density.