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.

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.

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.

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.

Probing the Electronic Properties and Interaction Landscapes in a Series of N-(Chlorophenyl)pyridinecarboxamides.

A 3 × 3 isomer grid of nine N-(chlorophenyl)pyridinecarboxamides (NxxCl) is reported with physicochemical studies and single crystal structures (Nx = pyridinoyl moiety; xCl = aminochlorobenzene ring; x = para-/meta-/ortho-), as synthesized by the reaction of the substituted p-/m-/o-pyridinecarbonyl chlorides (Nx) with p-/m-/o-aminochlorobenzenes (xCl). Several of the nine NxxCl crystal structures display structural similarities with their halogenated NxxX and methylated NxxM relatives (x = p-/m-/o-substitutions; X = F, Br; M = methyl). Indeed, five of the nine NxxCl crystal structures are isomorphous with their NxxBr analogues as the NpmCl/Br, NpoCl/Br, NmoCl/NmoBr, NopCl/Br, and NooCl/Br pairs. In the NxxCl series, the favored hydrogen bonding mode is aggregation by N-H···Npyridine interactions, though amide···amide intermolecular interactions are noted in NpoCl and NmoCl. For the NoxCl triad, intramolecular N-H···Npyridine interactions influence molecular planarity, whereas NppCl·H2O (as a monohydrate) exhibits O-H···O, N-H···O1W, and O1W-H···N interactions as the primary hydrogen bonding. Analysis of chlorine-containing compounds on the CSD is noted for comparisons. The interaction environments are probed using Hirshfeld surface analysis and contact enrichment studies. The melting temperatures (T m) depend on both the lattice energy and molecular symmetry (Carnelley's rule), and the melting points can be well predicted from a linear regression of the two variables. The relationships of the T m values with the total energy, the electrostatic component, and the strongest hydrogen bond components have been analyzed.

Biochemical and Structural Characterization of Chi-Class Glutathione Transferases: A Snapshot on the Glutathione Transferase Encoded by sll0067 Gene in the Cyanobacterium Synechocystis sp. Strain PCC 6803.

Glutathione transferases (GSTs) constitute a widespread superfamily of enzymes notably involved in detoxification processes and/or in specialized metabolism. In the cyanobacterium Synechocsytis sp. PCC 6803, SynGSTC1, a chi-class GST (GSTC), is thought to participate in the detoxification process of methylglyoxal, a toxic by-product of cellular metabolism. A comparative genomic analysis showed that GSTCs were present in all orders of cyanobacteria with the exception of the basal order Gloeobacterales. These enzymes were also detected in some marine and freshwater noncyanobacterial bacteria, probably as a result of horizontal gene transfer events. GSTCs were shorter of about 30 residues compared to most cytosolic GSTs and had a well-conserved SRAS motif in the active site (10SRAS13 in SynGSTC1). The crystal structure of SynGSTC1 in complex with glutathione adopted the canonical GST fold with a very open active site because the α4 and α5 helices were exceptionally short. A transferred multipolar electron-density analysis allowed a fine description of the solved structure. Unexpectedly, Ser10 did not have an electrostatic influence on glutathione as usually observed in serinyl-GSTs. The S10A variant was only slightly less efficient than the wild-type and molecular dynamics simulations suggested that S10 was a stabilizer of the protein backbone rather than an anchor site for glutathione.

Topological analysis of hydrogen bonds and weak interactions in protein helices via transferred experimental charge density parameters.

Helices represent the most abundant secondary structure motif in proteins and are often involved in various functional roles. They are stabilized by a network of hydrogen bonds between main chain carbonyl and amide groups. Several surveys scrutinized these hydrogen bonds, mostly based on the geometry of the interaction. Alternatively, the topological analysis of the electron density provides a powerful technique to investigate hydrogen bonds. For the first time, transferred experimental charge density parameters (ELMAM database) were used to carry out a topological analysis of the electron density in protein helices. New parameters for the description of the hydrogen bond geometry are proposed. Bonding contacts between the amide N and carbonyl O atoms (N···O) of helices, poorly addressed in the literature so far, were characterized from topological, geometrical, and local energetic analyses. Particularly, a geometrical criterion allowing for the discrimination between hydrogen bonds and N···O contacts is proposed.

Experimental and theoretical investigation of topological and energetic characteristics of Sb complexes reversibly binding molecular oxygen.

The experimental distribution of electron density in Ph(3)(4,5-OMe-3,6-Bu(t)-Cat)Sb·MeCN (1*) and Ph(3)(4,5-N(2)C(4)H(6)-3,6-Bu(t)-Cat)Sb·MeOH (2*) complexes was studied. According to atoms in molecules theory, the Sb-C(Ph), Sb-O(catecholate), and Sb···N(O) bonds are intermediate, whereas the O-C and C-C bonds are covalent, respectively. The energy of the Sb···N(MeCN) and Sb···O(MeOH) bonds are 7.0 and 11.3 kcal/mol according to the Espinosa equation. Density functional theory and Hartree-Fock calculations were carried out for a series of catecholate and amidophenolate complexes of antimony(V). It was shown that such calculations reliably reproduce geometrical and topological parameters and therefore can be used for a criterion search of dioxygen reversible binding by the catecholate and amidophenolate complexes of antimony(V). It was found that the "critical" value of the HOMO energy vary in the range from -5.197 to -5.061 eV for reversible binding of dioxygen complexes. This can serve as a thermodynamic criterion to predict the possibility of the dioxygen reversible binding by the catecholate and amidophenolate complexes of Sb(V). The HOMO energies correlate with the conversion of the catecholate and amidophenolate complexes in corresponding spiroendoperoxide derivatives as well. The contribution of the atom orbitals of the carbon atoms in the five-membered metallocycle to HOMO in complexes with different substitutes in the 4- and 5-positions of the catecholate ligand allows predicting the place of dioxygen addition.

Structural analysis and multipole modelling of quercetin monohydrate--a quantitative and comparative study.

The multipolar atom model, constructed by transferring the charge-density parameters from an experimental or theoretical database, is considered to be an easy replacement of the widely used independent atom model. The present study on a new crystal structure of quercetin monohydrate [2-(3,4-dihydroxyphenyl)-3,5,7-trihydroxy-4H-chromen-4-one monohydrate], a plant flavonoid, determined by X-ray diffraction, demonstrates that the transferred multipolar atom model approach greatly improves several factors: the accuracy of atomic positions and the magnitudes of atomic displacement parameters, the residual electron densities and the crystallographic figures of merit. The charge-density features, topological analysis and electrostatic interaction energies obtained from the multipole models based on experimental database transfer and periodic quantum mechanical calculations are found to compare well. This quantitative and comparative study shows that in the absence of high-resolution diffraction data, the database transfer approach can be applied to the multipolar electron density features very accurately.

2,2'-(ethane-1,2-diyl)bis[2-(5-bromothiophen-2-yl)-1,3-dioxolane] at 100 K refined using a multipolar atom model.

The title compound, C(16)H(16)Br(2)O(4)S(2), which is a precursor for the synthesis of oligothiophenes and their substituted homologues, was synthesized and its X-ray crystal structure determined at 100 K. The experimental electron-density parameters for the available atom types were transferred from the ELMAM2 database. The compound lies about an inversion centre, which coincides with the mid-point of a C-C bond. The molecules in the crystal are linked by several types of weak interactions; the largest contact surfaces are for H...H and H...Br.

A rush to explore protein-ligand electrostatic interaction energy with Charger.

The mutual penetration of electron densities between two interacting molecules complicates the computation of an accurate electrostatic interaction energy based on a pseudo-atom representation of electron densities. The numerical exact potential and multipole moment (nEP/MM) method is time-consuming since it performs a 3D integration to obtain the electrostatic energy at short interaction distances. Nguyen et al. [(2018), Acta Cryst. A74, 524-536] recently reported a fully analytical computation of the electrostatic interaction energy (aEP/MM). This method performs much faster than nEP/MM (up to two orders of magnitude) and remains highly accurate. A new program library, Charger, contains an implementation of the aEP/MM method. Charger has been incorporated into the MoProViewer software. Benchmark tests on a series of small molecules containing only C, H, N and O atoms show the efficiency of Charger in terms of execution time and accuracy. Charger is also powerful in a study of electrostatic symbiosis between a protein and a ligand. It determines reliable protein-ligand interaction energies even when both contain S atoms. It easily estimates the individual contribution of every residue to the total protein-ligand electrostatic binding energy. Glutathione transferase (GST) in complex with a benzophenone ligand was studied due to the availability of both structural and thermodynamic data. The resulting analysis highlights not only the residues that stabilize the ligand but also those that hinder ligand binding from an electrostatic point of view. This offers new perspectives in the search for mutations to improve the interaction between the two partners. A proposed mutation would improve ligand binding to GST by removing an electrostatic obstacle, rather than by the traditional increase in the number of favourable contacts.

Charge density and electrostatic interactions of fidarestat, an inhibitor of human aldose reductase.

The charge density and the topological features of fidarestat, an inhibitor of human aldose reductase, have been determined from ultra high-resolution X-ray diffraction data at 100 K. The modeled electron density was used to calculate the electrostatic interaction energy of fidarestat and its (2R,4S) stereoisomer with the human aldose reductase by using the ELMAM database as coded in the MoPro program. Such calculation may be extended to other protein complexes for which accurate high resolution X-ray data are available. The paper also discusses the hydrogen bonds in the fidarestat crystal. There are notably two hydrogen bonds with a pi system as an acceptor. All the chemical bonds and the intermolecular interactions, especially these two pi...H bonds, have been quantitatively studied by topological analysis. The three-dimensional electrostatic potential calculated on the molecular surface emphasizes the preferential polar binding sites of fidarestat. Theses interacting features in the molecule are crucial for drug-receptor recognition. The interactions between chemical groups in the crystal are also analyzed by computing the electrostatic energy using the latest advancements of the MoPro crystallographic software. The complexes of fidarestat and its (2R,4S) stereoisomer with human aldose reductase were modeled with a multipolar atom model transferred from our experimental electron density database. Accurate estimation of electrostatic interaction energy between inhibitors and the main residues of the protein active site is derived from this high detail level of the electron density.

Elucidation of the phosphate binding mode of DING proteins revealed by subangstrom X-ray crystallography.

PfluDING is a bacterial protein isolated from Pseudomonas fluorescens that belongs to the DING protein family, which is ubiquitous in eukaryotes and extends to prokaryotes. DING proteins and PfluDING have very similar topologies to phosphate Solute Binding Proteins (SBPs). The three-dimensional structure of PfluDING was obtained at subangstrom resolution (0.88 and 0.98 A) at two different pH's (4.5 and 8.5), allowing us to discuss the hydrogen bond network that sequesters the phosphate ion in the binding site. From this high resolution data, we experimentally elucidated the molecular basis of phosphate binding in phosphate SBPs. The phosphate ion is tightly bound to the protein via 12 hydrogen bonds between phosphate oxygen atoms and OH and NH groups of the protein. The proton on one oxygen atom of the phosphate dianion forms a 2.5 A low barrier hydrogen bond with an aspartate, with the energy released by forming this strong bond ensuring the specificity for the dianion even at pH 4.5. In particular, contrary to previous theories on phosphate SBPs, accurate electrostatic potential calculations show that the binding cleft is positively charged. PfluDING structures reveal that only dibasic phosphate binds to the protein at both acidic and basic phosphate, suggesting that the protein binding site environment stabilizes the HPO(4)(2-) form of phosphate.

Charge density and electrostatic potential analyses in paracetamol.

The electron density of monoclinic paracetamol was derived from high-resolution X-ray diffraction at 100 K. The Hansen-Coppens multipole model was used to refine the experimental electron density. The topologies of the electron density and the electrostatic potential were carefully analyzed. Numerical and analytical procedures were used to derive the charges integrated over the atomic basins. The highest charge magnitude (-1.2 e) was found for the N atom of the paracetamol molecule, which is in agreement with the observed nucleophilic attack occurring in the biological media. The electric field generated by the paracetamol molecule was used to calculate the atomic charges using the divergence theorem. This was simultaneously applied to estimate the total electrostatic force exerted on each atom of the molecule by using the Maxwell stress tensor. The interaction electrostatic energy of dimers of paracetamol in the crystal lattice was also estimated.

4-Benzoyl-3,4-dihydro-2H-1,4-benzoxazine-2-carbonitrile: refinement using a multipolar atom model.

The structural model for the title compound, C(16)H(12)N(2)O(2), was refined using a multipolar atom model transferred from an experimental electron-density database. The refinement showed some improvements of crystallographic statistical indices when compared with a conventional spherical neutral-atom refinement. The title compound adopts a half-chair conformation. The amide N atom lies almost in the plane defined by the three neighbouring C atoms. In the crystal structure, molecules are linked by weak intermolecular C-H...O and C-H...pi hydrogen bonds.

DORI Reveals the Influence of Noncovalent Interactions on Covalent Bonding Patterns in Molecular Crystals Under Pressure.

The study of organic molecular crystals under high pressure provides fundamental insight into crystal packing distortions and reveals mechanisms of phase transitions and the crystallization of polymorphs. These solid-state transformations can be monitored directly by analyzing electron charge densities that are experimentally obtained at high pressure. However, restricting the analysis to the featureless electron density does not reveal the chemical bonding nature and the existence of intermolecular interactions. This shortcoming can be resolved by the use of the DORI (density overlap region indicator) descriptor, which is capable of simultaneously detecting both covalent patterns and noncovalent interactions from electron density and its derivatives. Using the biscarbonyl[14]annulene crystal under pressure as an example, we demonstrate how DORI can be exploited on experimental electron densities to reveal and monitor changes in electronic structure patterns resulting from molecular compression. A novel approach based on a flood-fill-type algorithm is proposed for analyzing the topology of the DORI isosurface. This approach avoids the arbitrary selection of DORI isovalues and provides an intuitive way to assess how compression packing affects covalent bonding in organic solids.

A method to estimate statistical errors of properties derived from charge-density modelling.

Estimating uncertainties of property values derived from a charge-density model is not straightforward. A methodology, based on calculation of sample standard deviations (SSD) of properties using randomly deviating charge-density models, is proposed with the MoPro software. The parameter shifts applied in the deviating models are generated in order to respect the variance-covariance matrix issued from the least-squares refinement. This `SSD methodology' procedure can be applied to estimate uncertainties of any property related to a charge-density model obtained by least-squares fitting. This includes topological properties such as critical point coordinates, electron density, Laplacian and ellipticity at critical points and charges integrated over atomic basins. Errors on electrostatic potentials and interaction energies are also available now through this procedure. The method is exemplified with the charge density of compound (E)-5-phenylpent-1-enylboronic acid, refined at 0.45 Šresolution. The procedure is implemented in the freely available MoPro program dedicated to charge-density refinement and modelling.

What is the relationship between the bioaccumulation of chemical contaminants in the variegated scallop Mimachlamys varia and its health status? A study carried out on the French Atlantic coast using the Path ComDim model.

Increasing activity along the French Atlantic coast has led to chronic pollution with, in particular, mixtures of contaminants such as hydrocarbons, phytosanitary products, PCBs and heavy metals. Based on previous research, pollution biomarkers were used in this study as they can indicate health status when monitoring the impact of pollutants on coastal species such as the marine bivalve Mimachlamys varia. Mollusc bivalves were sampled in March 2016, in open and semi-open areas (a harbour zone), from thirteen sites which differed in terms of their level of pollution, and were located along the Atlantic coast from Brittany down to the Nouvelle-Aquitaine region. First, analyses of heavy metals and organic contaminants (e.g. pesticides, polycyclic aromatic hydrocarbons, polychlorobiphenyl) in the digestive gland of bivalves were performed. Second, biochemical assays were used to study defence biomarkers: oxidative stress with Superoxide Dismutase (SOD), detoxification of organic compounds with Glutathione-S Transferase (GST), lipid peroxidation with Malondialdehyde (MDA), and immune processes with Laccase. In addition to the biochemical assays, a genetic approach was used to measure genetic diversity (haplotype and nucleotide diversity) at each site. Biomarker assays and genetic diversity were correlated with the chemical contaminants in bivalves using the Path-ComDim statistical model. Our results showed specific correlations between biochemical assays in the digestive glands with heavy metal contaminants, and between genetic diversity and organic pollution. Blocks of responses were analysed for correlations in order to develop standardized tools and guidelines that could improve our understanding of the short-term and long-term impact of contaminants on physiological parameters.

On the application of an experimental multipolar pseudo-atom library for accurate refinement of small-molecule and protein crystal structures.

With an increasing number of biomacromolecular crystal structures being measured to ultra-high resolution, it has become possible to extend to large systems experimental charge-density methods that are usually applied to small molecules. A library has been built of average multipole populations describing the electron density of chemical groups in all 20 amino acids found in proteins. The library uses the Hansen & Coppens multipolar pseudo-atom model to derive molecular electron density and electrostatic potential distributions. The library values are obtained from several small peptide or amino acid crystal structures refined against ultra-high-resolution X-ray diffraction data. The library transfer is applied automatically in the MoPro software suite to peptide and protein structures measured at atomic resolution. The transferred multipolar parameters are kept fixed while the positional and thermal parameters are refined. This enables a proper deconvolution of thermal motion and valence-electron-density redistributions, even when the diffraction data do not extend to subatomic resolution. The use of the experimental library multipolar atom model (ELMAM) also has a major impact on crystallographic structure modelling in the case of small-molecule crystals at atomic resolution. Compared to a spherical-atom model, the library transfer results in a more accurate crystal structure, notably in terms of thermal displacement parameters and bond distances involving H atoms. Upon transfer, crystallographic statistics of fit are improved, particularly free R factors, and residual electron-density maps are cleaner.