Probing the Isolobal Relation between Cp'''NiP3 and White Phosphorus by Experimental Charge Density Analysis.

An in-depth analysis of the description of bonding within Cp'''Ni-cyclo-P3 (Cp'''=1,2,4-tri-tert-butylcyclopentadienyl, [Ni]P3) employing X-ray diffraction based multipolar modeling, density functional theory (DFT) as well as an "experimental wavefunction" obtained from X-ray restrained wavefunction (XRW) fitting is presented. The results are compared to DFT calculations on white phosphorus - an isolobal analogue to [Ni]P3. A complementary bonding analysis shows insights into the reactivity of [Ni]P3. The isolobal principle is reflected in every aspect of our analysis and the employed methods seamlessly predict the differences in reactivity of [Ni]P3 and P4. Crystallographic modeling, solid-state NMR, and DFT calculations describe the dynamic behavior of the cyclo-P3 unit in the title molecule.

Machine learning-driven identification of drugs inhibiting cytochrome P450 2C9.

Cytochrome P450 2C9 (CYP2C9) is a major drug-metabolizing enzyme that represents 20% of the hepatic CYPs and is responsible for the metabolism of 15% of drugs. A general concern in drug discovery is to avoid the inhibition of CYP leading to toxic drug accumulation and adverse drug-drug interactions. However, the prediction of CYP inhibition remains challenging due to its complexity. We developed an original machine learning approach for the prediction of drug-like molecules inhibiting CYP2C9. We created new predictive models by integrating CYP2C9 protein structure and dynamics knowledge, an original selection of physicochemical properties of CYP2C9 inhibitors, and machine learning modeling. We tested the machine learning models on publicly available data and demonstrated that our models successfully predicted CYP2C9 inhibitors with an accuracy, sensitivity and specificity of approximately 80%. We experimentally validated the developed approach and provided the first identification of the drugs vatalanib, piriqualone, ticagrelor and cloperidone as strong inhibitors of CYP2C9 with IC values <18 µM and sertindole, asapiprant, duvelisib and dasatinib as moderate inhibitors with IC50 values between 40 and 85 µM. Vatalanib was identified as the strongest inhibitor with an IC50 value of 0.067 µM. Metabolism assays allowed the characterization of specific metabolites of abemaciclib, cloperidone, vatalanib and tarafenacin produced by CYP2C9. The obtained results demonstrate that such a strategy could improve the prediction of drug-drug interactions in clinical practice and could be utilized to prioritize drug candidates in drug discovery pipelines.

New Crystal Form of Human Neuropilin-1 b1 Fragment with Six Electrostatic Mutations Complexed with KDKPPR Peptide Ligand.

Neuropilin 1 (NRP1), a cell-surface co-receptor of a number of growth factors and other signaling molecules, has long been the focus of attention due to its association with the development and the progression of several types of cancer. For example, the KDKPPR peptide has recently been combined with a photosensitizer and a contrast agent to bind NRP1 for the detection and treatment by photodynamic therapy of glioblastoma, an aggressive brain cancer. The main therapeutic target is a pocket of the fragment b1 of NRP1 (NRP1-b1), in which vascular endothelial growth factors (VEGFs) bind. In the crystal packing of native human NRP1-b1, the VEGF-binding site is obstructed by a crystallographic symmetry neighbor protein, which prevents the binding of ligands. Six charged amino acids located at the protein surface were mutated to allow the protein to form a new crystal packing. The structure of the mutated fragment b1 complexed with the KDKPPR peptide was determined by X-ray crystallography. The variant crystallized in a new crystal form with the VEGF-binding cleft exposed to the solvent and, as expected, filled by the C-terminal moiety of the peptide. The atomic interactions were analyzed using new approaches based on a multipolar electron density model. Among other things, these methods indicated the role played by Asp320 and Glu348 in the electrostatic steering of the ligand in its binding site. Molecular dynamics simulations were carried out to further analyze the peptide binding and motion of the wild-type and mutant proteins. The simulations revealed that specific loops interacting with the peptide exhibited mobility in both the unbound and bound forms.

Polarization of Electron Density Databases of Transferable Multipolar Atoms.

Polarizability is a key molecular property involved in either macroscopic (i.e., dielectric constant) and microscopic properties (i.e., interaction energies). In rigid molecules, this property only depends on the ability of the electron density (ED) to acquire electrostatic moments in response to applied electric fields. Databases of transferable electron density fragments are a cheap and efficient way to access molecular EDs. This approach is rooted in the relative conservation of the atomic ED between different molecules, termed transferability principle. The present work discusses the application of this transferability principle to the polarizability, an electron density-derived property, partitioned in atomic contributions using the Quantum Theory of Atoms In Molecules topology. The energetic consequences of accounting for in situ deformation (polarization) of database multipolar atoms are investigated in detail by using a high-quality quantum chemical benchmark.

Syn vs Anti Carboxylic Acids in Hybrid Peptides: Experimental and Theoretical Charge Density and Chemical Bonding Analysis.

A comparative study of syn vs anti carboxylic acids in hybrid peptides based on experimental electron density studies and theoretical calculations shows that, in the anti form, all three bond angles surrounding Ccarboxyl of the -COOH group are close to ∼120°, as expected for a C-sp2 atom, whereas in the syn form, the ∠Cα-C(O)-Ohydroxyl angle is significantly smaller by 5-10°. The oxygen atom in the carboxyl group is more electronegative in the anti form, so the polarity of the acidic O-H bond is higher in the anti form compared to the syn form, as observed within the limitations of H atom treatment in X-ray diffraction. Consequently, the investigated anti carboxylic acid forms the strongest O-H···O hydrogen bond among all model compounds. Furthermore, according to natural bond orbital analysis, the oxygen lone pairs are clearly nonequivalent, as opposed to the general notion of hybridization of equivalent sp2 and sp3 lone pairs on carbonyl or hydroxyl oxygen atoms. The hybridization of the lone pairs is directly related to the directionality and strength of hydrogen bonds.

Bonding in Uranium(V) Hexafluoride Based on the Experimental Electron Density Distribution Measured at 20 K.

The electron density distribution of [PPh4][UF6] was obtained from high-resolution X-ray diffraction data measured at 20 K. The electron density was modeled with an augmented Hansen-Coppens multipolar formalism. Topological analysis reveals that the U-F bond is of incipient covalent nature. Theoretical calculations add further support to the bonding description gleaned from the experimental model. The impact of the uranium anomalous dispersion terms on the refinement is also discussed.

Carbohydrate-based peptidomimetics targeting neuropilin-1: Synthesis, molecular docking study and in vitro biological activities.

Neuropilin-1 (NRP-1), a transmembrane glycoprotein acting as a co-receptor of VEGF-A, is expressed by cancer and angiogenic endothelial cells and is involved in the angiogenesis process. Taking advantage of functionalities and stereodiversities of sugar derivatives, the design and the synthesis of carbohydrate based peptidomimetics are here described. One of these compounds (56) demonstrated inhibition of VEGF-A165 binding to NRP-1 (IC50=39µM) and specificity for NRP-1 over VEGF-R2. Biological evaluations were performed on human umbilical vein endothelial cells (HUVECs) through activation of downstream proteins (AKT and ERK phosphorylation), viability/proliferation assays and in vitro measurements of anti-angiogenic abilities.

Betulin Phosphonates; Synthesis, Structure, and Cytotoxic Activity.

Betulin derivatives are a widely studied group of compounds of natural origin due to their wide spectrum of biological activities. This paper describes new betulin derivatives, containing a phosphonate group. The allyl-vinyl isomerization and synthesis of acetylenic derivatives have been reported. Structural identification of products as E and Z isomers has been carried out using ¹H-, (13)C-, (31)P-NMR, and crystallographic analysis. The crystal structure in the orthorhombic space group and analysis of crystal packing contacts for 29-diethoxyphosphoryl-28-cyclopropylpropynoyloxy-lup-20E(29)-en-3ß-ol 8a are reported. All new compounds were tested in vitro for their antiproliferative activity against human T47D (breast cancer), SNB-19 (glioblastoma), and C32 (melanoma) cell lines.

Cholesterol oxidase: ultrahigh-resolution crystal structure and multipolar atom model-based analysis.

Examination of protein structure at the subatomic level is required to improve the understanding of enzymatic function. For this purpose, X-ray diffraction data have been collected at 100 K from cholesterol oxidase crystals using synchrotron radiation to an optical resolution of 0.94 Å. After refinement using the spherical atom model, nonmodelled bonding peaks were detected in the Fourier residual electron density on some of the individual bonds. Well defined bond density was observed in the peptide plane after averaging maps on the residues with the lowest thermal motion. The multipolar electron density of the protein-cofactor complex was modelled by transfer of the ELMAM2 charge-density database, and the topology of the intermolecular interactions between the protein and the flavin adenine dinucleotide (FAD) cofactor was subsequently investigated. Taking advantage of the high resolution of the structure, the stereochemistry of main-chain bond lengths and of C=O···H-N hydrogen bonds was analyzed with respect to the different secondary-structure elements.

Crystal structure of 2-amino-5-methyl-sulfanyl-1,3,4-thia-diazol-3-ium chloride monohydrate.

The title salt, C3H6N3S2 (+)·Cl(-)·H2O, crystallized with two organic cations, two chloride anions and two water mol-ecules in the asymmetric unit. The methyl C atoms deviate from their respective bound ring planes by 0.081 and 0.002 Å. In the crystal, the components are connected via N-H⋯O, N-H⋯Cl and O-H⋯Cl hydrogen bonds, forming sheets lying parallel to (100). The sheets are linked into bilayers by O-H⋯Cl hydrogen bonds involving the chloride ions and water mol-ecules. Within the bilayers there are π-π inter-actions [inter-centroid distances = 3.4654 (4) and 3.4789 (4) Å] involving inversion-related cations.

Bis(2-amino-4-methyl-6-oxo-3,6-di-hydro-pyrimidin-1-ium) sulfate monohydrate.

In the title hydrated mol-ecular salt, 2C5H8N3O(+)·SO4 (2-)·H2O, the components are linked by N-H⋯Os and Ow-H⋯Os (s = sulphate, w = water) hydrogen bonds, generating a layer by a+b+c and 2a-b translations. The cations are arranged nearly in parallel and show displaced π-π stacking centroid-centroid distance = 4.661 (2) Šbetween adjacent layers.

Poly[diaquatris(µ6-4,6-dioxo-1,4,5,6-tetra-hydro-1,3,5-triazine-2-carboxylato)tripotassium].

The asymmetric unit of the title compound, [K3(C4H2N3O4)3(H2O)2] n , contains two potassium cations (one in general position, one located on a twofold rotation axis), one and a half oxonate anions (the other half generated by twofold symmetry) and one water mol-ecule. As a result of the twofold symmetry, one H atom of the symmetric anion is statistically occupied. Both potassium cations are surrounded by eight oxygen atoms in the form of distorted polyhedra. Adjacent cations are inter-connected by oxygen bridges, generating layers parallel to (100). The aromatic ring system of the oxonate anions link these layers into a network structure. The crystal packing is stabilized by N-H⋯O, O-H⋯O and O-H⋯N hydrogen bonds, three of which are bifurcated. In addition, inter-molecular π-π stacking inter-actions exist between neighboring aromatic rings with a centroid-centroid distance of 3.241 (2) Å.

Experimental and theoretical charge density analysis of a bromoethyl sulfonium salt.

Bromoethyl sulfonium trifluoromethanesulfonate is a salt complex in which a sulfur atom makes three covalent bonds. This molecule has been proved to act as an efficient annulation reagent which results in formation of synthetically challenging and pharmaceutically important 4-, 5-, 6-, and 7-membered heterocycles in excellent yields. The charge density of the molecule was determined from both experimentally and theoretically derived diffraction data. The stereochemistry and electron density topology of the sulfonium group was analyzed. To understand the chemical reactivity of the molecule, the electrostatic potential difference between the two carbon atoms of the bromoethyl group was investigated. It has been considered that the hydrogen atoms on the carbon atom bound to sulfur are more acidic in character due to their vicinity with the triply covalently bonded positively charged sulfur atom. The electropositivity of the S-attached and Br-attached methylene groups are compared in the experimental and theoretical charge densities using topological atomic charges and electrostatic potential at the molecular surface.

The organic-inorganic hybrid material 1-cyclohexylpiperazine-1,4-diium tetrachloridozincate.

In the crystal structure of the title organic-inorganic hybrid material, (C10H22N2)[ZnCl4], the tetrachloridozincate anions and 1-cyclohexylpiperazine-1,4-diium dications are interconnected via N-H...Cl and C-H...Cl hydrogen bonds to form layers parallel to the (001) plane. The cyclohexyl groups from adjacent chains interdigitate, thus building the three-dimensional structure. The piperazinium and cyclohexyl rings exhibit regular spatial chair conformations. The title salt was also characterized by FT-IR and Raman spectroscopic analyses.

o-Phenyl-enediaminium chloride nitrate.

In the title mol-ecular salt, C6H10N2 (2+)·NO3 (-)·Cl(-), the complete cation is generated by a crystallographic mirror plane. The complete nitrate ion is also generated by reflection, with the N atom and one O atom lying on the mirror plane; the chloride ion also lies on the reflection plane. In the crystal, the components are linked by N-H⋯Cl and N-H⋯(N,O) hydrogen bonds, forming (001) layers with the benzene rings projecting into the inter-layer regions. The layers are linked by weak C-H⋯O hydrogen bonds, generating a three-dimensional network.

Deciphering the driving forces in crystal packing by analysis of electrostatic energies and contact enrichment ratios.

Hirshfeld surface analysis is a widely used tool for identifying the types of intermolecular contacts that contribute most significantly to crystal packing stabilization. One useful metric for analyzing these contacts is the contact enrichment descriptor, which indicates the types of contacts that are over- or under-represented. In this statistical study, enrichment ratios were combined with electrostatic energy (Eelec) data for a variety of compound families. To compute the electrostatic interaction energy between atoms, charge density models from the ELMAM2 database of multipolar atoms were used. As expected, strong hydrogen bonds such as O/N-H...N and O/N-H...O typically display large enrichment values and have the most negative (i.e. favorable) electrostatic energies. Conversely, contacts that are repulsive from an electrostatic perspective are usually the most under-represented. Analyzing the enrichment ratio and electrostatic energy indicators was shown to help identify which favorable contacts are the most competitive with each other. For weaker interactions, such as hydrophobic contacts, the behavior is less clear cut and can depend on other factors such as the chemical content of the molecule. The anticorrelation between contact enrichment and Eelec is generally lost for weaker contacts. However, we observed that C...C contacts are often enriched in crystal structures containing heterocycles, despite the low electrostatic attraction. For molecules with only weak hydrogen bond donors/acceptors and hydrophobic groups, the correlation between contact enrichment and Eelec is still evident for the strongest of these interactions. However, there are some exceptions where the most favorable contacts from an electrostatic perspective are not the most over-represented. This can occur in cases where the shape of the molecule is complex or elongated, favoring dispersion forces and shape complementarity in the packing.

Racemic cis-bis-[bis-(pyrimidin-2-yl)amine-κN]bis(dicyanamido-κN1)iron(II) dihydrate: synthesis, crystal structure and Hirshfeld surface analysis.

The title compound, [Fe(C2N3)2(C8H7N5)2]·2H2O, has been synthesized solvothermally and characterized by single-crystal X-ray diffraction. The octa-hedral iron coordination polyhedron contains two di(pyrimidin-2-yl)amine ligands coordinated in a bidentate fashion, and two monodentate dicyanimido ligands, each coordinated via a terminal N atom, with the latter in a cis orientation. The ligand configuration about the iron atom is chiral, although the compound crystallizes as a racemic mixture: the Fe-N distances (> 2.07 Å) are characteristic of high-spin iron(II). In the crystal, an extensive series of N-H⋯N, O-H⋯N and O-H⋯O hydrogen bonds links the independent mol-ecular components into a three-dimensional framework. The H atoms of both water mol-ecules are disordered. The structure also features some π-π and anion-π inter-actions. The inter-molecular inter-actions were investigated by Hirshfeld surface analysis and two-dimensional fingerprint plots. Comparisons are made with some related compounds.

Experimental and database-transferred electron-density analysis and evaluation of electrostatic forces in coumarin-102 dye.

The electron-density distribution of a new crystal form of coumarin-102, a laser dye, has been investigated using the Hansen-Coppens multipolar atom model. The charge density was refined versus high-resolution X-ray diffraction data collected at 100 K and was also constructed by transferring the charge density from the Experimental Library of Multipolar Atom Model (ELMAM2). The topology of the refined charge density has been analysed within the Bader `Atoms In Molecules' theory framework. Deformation electron-density peak heights and topological features indicate that the chromen-2-one ring system has a delocalized π-electron cloud in resonance with the N (amino) atom. The molecular electrostatic potential was estimated from both experimental and transferred multipolar models; it reveals an asymmetric character of the charge distribution across the molecule. This polarization effect is due to a substantial charge delocalization within the molecule. The molecular dipole moments derived from the experimental and transferred multipolar models are also compared with the liquid and gas-phase dipole moments. The substantial molecular dipole moment enhancements observed in the crystal environment originate from the crystal field and from intermolecular charge transfer induced and controlled by C-H···O and C-H···N intermolecular hydrogen bonds. The atomic forces were integrated over the atomic basins and compared for the two electron-density models.

π-Hole bonding in a new co-crystal hydrate of gallic acid and pyrazine: static and dynamic charge density analysis.

A new cocrystal hydrate of gallic acid with pyrazine (4GA, Py, 4H2O; GA4PyW4) was obtained and characterized by single crystal X-ray diffraction. In addition to structure determination, experimental charge density analysis was carried out in terms of Multipole Modelling (MP), X-ray wavefunction refinement (XWR) and maximum entropy method (MEM). As a part of XWR, the structural refinement via Hirshfeld atom refinement was carried out and resulted in O-H bond lengths close to values from neutron diffraction. A systematic comparison of molecular conformations and aromatic interactions in this new cocrystal hydrate was performed with other existing polymorphs of gallic acid. In GA4PyW4, the two symmetry-independent gallic acid molecules have a syn COOH orientation and form the common (COOH)2 dimeric synthon. The carboxyl C atom displays the characteristics of π-holes with electropositive regions above and below the molecular plane and engages in acceptor-donor interactions with oxygen atoms of acidic O-H groups and phenol groups of neighbouring gallic acid molecules. The signature of the π-hole was identified from experimental charge density analysis, both in static density maps in MP and XWR as well as dynamic density in MEM, but it cannot be pinned down to a specific atom-atom interaction. This study presents the first comparison between an XWR and a MEM experimental electron-density determination.

Theoretical search of crystal polymorphs of temozolomide.

Possible polymorphic forms of the chemotherapy drug, temozolomide were predicted from the ab initio and DFT methods. The lattice minimization via distributed multipole analysis was carried out for the hypothetical generated structures. A crystal with unit cell parameters close to the real one and of same space group was retrieved, with partly similar packing and interactions. The analysis of inter molecular interaction (through Hirshfeld surface) and electrostatic potential reveals the complementary sites in the molecule. The 26 predicted structures were analyzed with respect to two computed lattice energies and hydrogen-bond propensity. The lattice energy of the real crystal [EXP] packing ranked number 6 compared on the basis of DMACRYS software and number 3 on the basis of the total lattice energy issued from the Crystalexplorer17 software at the B3LYP/6-31G∗∗ level of theory. The molecule has two strong hydrogen bond donors and five strong acceptors. The predicted packings are stabilized by one or two strong N-H…O/N-H…N as well as weak C-H…O/C-H…N and H…π hydrogen bonds. While the real structure with Z' = 1, EXP, forms only one strong H-bond (N-H…O=C), several of the predicted packings form two strong H-bonds. Two predicted crystal packings have unit cell parameters close to the real structure, one of them shares several common intermolecular interactions.