Exploiting the full potential of cryo-EM maps.

The Coulomb potential maps generated by electron microscopy (EM) experiments contain not only information about the position but also about the charge state of the atom. This feature of EM maps allows the identification of specific ions and the protonation state of amino acid side chains in the sample. Here, we summarize qualitative observations of charges in EM maps, discuss the difficulties in interpreting the charge in Coulomb potential maps with respect to distinguishing it from radiation damage, and outline considerations to implement the correct charge in fitting algorithms.

How do density functionals affect the Hirshfeld atom refinement?

In this work, the effect of mixing different amounts of Hartree-Fock (HF) exchange with hybrid density functionals applied to the Hirshfeld atom refinement (HAR) of urea and oxalic acid dihydrate is explored. Together, the influence of using different basis sets, methods (including MP2 and HF) and cluster sizes (to model bulk effects) is studied. The results show that changing the amount of HF exchange, no matter the level of theory, has an impact almost exclusively on the H atom refinement parameters. Contrary to pure quantum mechanical calculations where good geometries are obtained with intermediate HF exchange mixtures, in the HAR the best match with neutron diffraction reference values is not necessarily found for these admixtures. While the non-hydrogen covalent bond lengths are insensitive to the combination of method or basis set employed, the X-H bond lengths always increase proportionally to the HF exchange for the analysed systems. This outcome is opposite to what is normally observed from geometry optimisations, i.e., shorter bonds are obtained with greater HF exchange. Additionally, the thermal ellipsoids tend to shrink with larger HF exchange, especially for the H atoms involved in strong hydrogen bonding. Thus, it may be the case that the development of density functionals or basis sets suitable for quantum crystallography should take a different path than those fitted for quantum chemistry calculations.

Electron density is not spherical: the many applications of the transferable aspherical atom model.

In the simplest approach and at low resolution, the electron density of an atom can be approximated by a sphere. However, the resolutions currently achieved in X-ray macromolecular crystallography reach atomic resolutions, where a more sophisticated approach is inevitable. This review summarizes the available electron density and scattering factor models and their applications with emphasis on the transferable aspherical atom model (TAAM).

X-ray wavefunction refinement and comprehensive structural studies on bromo-substituted analogues of 2-deoxy-d-glucose in solid state and solution.

The structural studies on two bromo-substituted derivatives of 2-deoxy-d-glucose (2-DG), namely 2-deoxy-2-bromo-d-glucose (2-BG) and 2-deoxy-2-bromo-d-mannose (2-BM) are described. 2-DG itself is an inhibitor of hexokinase, the first enzyme in the glycolysis process, playing a vital role in both cancer cell metabolism and viral replication in host cells. Because of that, 2-DG derivatives are considered as potential anti-cancer and anti-viral drugs. An X-ray quantum crystallography approach allowed us to obtain more accurate positions of hydrogen atoms by applying Hirshfeld atom refinement, providing a better description of hydrogen bonding even in the case of data from routine X-ray experiments. Obtained structures showed that the introduction of bromine at the C2 position in the pyranose ring has a minor influence on its conformation but still, it has a noticeable effect on the crystal structure. Bromine imposes the formation of a layered supramolecular landscape containing hydrogen bonds, which involves the bromine atom. Periodic DFT calculations of cohesive and interaction energies (at the B3LYP level of theory) have supported these findings and highlighted energetic changes upon bromine substitution. Based on molecular wavefunction from the refinement, we calculated the electrostatic potential, Laplacian, and ELI-D, and applied them to charge-density studies, which confirmed the geometry of hydrogen bonding and involvement of the bromine atom with these intermolecular interactions. NMR studies in the solution show that both compounds do not display significant differences in their anomeric equilibria compared to 2-DG, and the pyranose ring puckering is similar in both aqueous and solid state.

Correction: X-ray wavefunction refinement and comprehensive structural studies on bromo-substituted analogues of 2-deoxy-d-glucose in solid state and solution.

[This corrects the article DOI: 10.1039/D1RA08312K.].

The PI3K pathway preserves metabolic health through MARCO-dependent lipid uptake by adipose tissue macrophages.

Adipose tissue macrophages (ATMs) display tremendous heterogeneity depending on signals in their local microenvironment and contribute to the pathogenesis of obesity. The phosphoinositide 3-kinase (PI3K) signalling pathway, antagonized by the phosphatase and tensin homologue (PTEN), is important for metabolic responses to obesity. We hypothesized that fluctuations in macrophage-intrinsic PI3K activity via PTEN could alter the trajectory of metabolic disease by driving distinct ATM populations. Using mice harbouring macrophage-specific PTEN deletion or bone marrow chimeras carrying additional PTEN copies, we demonstrate that sustained PI3K activity in macrophages preserves metabolic health in obesity by preventing lipotoxicity. Myeloid PI3K signalling promotes a beneficial ATM population characterized by lipid uptake, catabolism and high expression of the scavenger macrophage receptor with collagenous structure (MARCO). Dual MARCO and myeloid PTEN deficiencies prevent the generation of lipid-buffering ATMs, reversing the beneficial actions of elevated myeloid PI3K activity in metabolic disease. Thus, macrophage-intrinsic PI3K signalling boosts metabolic health by driving ATM programmes associated with MARCO-dependent lipid uptake.

On the accuracy and precision of X-ray and neutron diffraction results as a function of resolution and the electron density model.

X-ray diffraction is the main source of three-dimensional structural information. In total, more than 1.5 million crystal structures have been refined and deposited in structural databanks (PDB, CSD and ICSD) to date. Almost 99.7% of them were obtained by approximating atoms as spheres within the independent atom model (IAM) introduced over a century ago. In this study, X-ray datasets for single crystals of hydrated α-oxalic acid were refined using several alternative electron density models that abandon the crude spherical approximation: the multipole model (MM), the transferable aspherical atom model (TAAM) and the Hirshfeld atom refinement (HAR) model as a function of the resolution of X-ray data. The aspherical models (MM, TAAM, HAR) give far more accurate and precise single-crystal X-ray results than IAM, sometimes identical to results obtained from neutron diffraction and at low resolution. Hence, aspherical approaches open new routes for improving existing structural information collected over the last century.

Universal Method for Electrostatic Interaction Energies Estimation with Charge Penetration and Easily Attainable Point Charges.

Our new model of electron density augmented by point charges (aug-PROmol) provides an estimation of electrostatic interaction energies including penetration effects ( ChemPhysChem 2016, 17, 2455-2460). In this paper we prove that it can be applied using sources of point charges other than those from direct restrained fitting to electrostatic potential (RESP). We used a newly established databank of tabulated invariom point charges and a widely known semiempirical method. Both sources perform equivalently to the basic aug-PROmol method as well as to reference energies at the DFT-SAPT/aug-cc-pVTZ level of theory. This is possible due to the universal character of the penetration model included in the aug-PROmol. Aug-PROmol may become a basis for development of new nonbonded terms in force fields or a high success rate scoring function.

Quantum Crystallography: Current Developments and Future Perspectives.

Crystallography and quantum mechanics have always been tightly connected because reliable quantum mechanical models are needed to determine crystal structures. Due to this natural synergy, nowadays accurate distributions of electrons in space can be obtained from diffraction and scattering experiments. In the original definition of quantum crystallography (QCr) given by Massa, Karle and Huang, direct extraction of wavefunctions or density matrices from measured intensities of reflections or, conversely, ad hoc quantum mechanical calculations to enhance the accuracy of the crystallographic refinement are implicated. Nevertheless, many other active and emerging research areas involving quantum mechanics and scattering experiments are not covered by the original definition although they enable to observe and explain quantum phenomena as accurately and successfully as the original strategies. Therefore, we give an overview over current research that is related to a broader notion of QCr, and discuss options how QCr can evolve to become a complete and independent domain of natural sciences. The goal of this paper is to initiate discussions around QCr, but not to find a final definition of the field.

Validation of X-ray Wavefunction Refinement.

In this work, the quality of the electron density in crystals reconstructed by the multipolar model (MM) and by X-ray wavefunction refinement (XWR) is tested on a set of high-resolution X-ray diffraction data sets of four amino acids and six tripeptides. It results in the first thorough validation of XWR. Agreement statistics, figures of merit, residual- and deformation-density maps, as well as atomic displacement parameters are used to measure the quality of the reconstruction relative to the measured structure factors. Topological analysis of the reconstructed density is carried out to obtain atomic and bond-topological properties, which are subsequently compared to the values derived from benchmarking periodic DFT geometry optimizations. XWR is simultaneously in better agreement than the MM with both benchmarking theory and the measured diffraction pattern. In particular, the obvious problems with the description of polar bonds in the MM are significantly reduced by using XWR. Similarly, modeling of electron density in the vicinity of hydrogen atoms with XWR is visibly improved.

Multi-temperature study of potassium uridine-5'-monophosphate: electron density distribution and anharmonic motion modelling.

Uridine, a nucleoside formed of a uracil fragment attached to a ribose ring via a ß-N1-glycosidic bond, is one of the four basic components of ribonucleic acid. Here a new anhydrous structure and experimental charge density distribution analysis of a uridine-5'-monophosphate potassium salt, K(UMPH), is reported. The studied case constitutes the very first structure of a 5'-nucleotide potassium salt according to the Cambridge Structural Database. The excellent crystal quality allowed the collection of charge density data at various temperatures, i.e. 10, 100, 200 and 300 K on one single crystal. Crystal structure and charge density data were analysed thoroughly in the context of related literature-reported examples. Detailed analysis of the charge density distribution revealed elevated anharmonic motion of part of the uracil ring moiety relatively weakly interacting with the neighbouring species. The effect was manifested by alternate positive and negative residual density patterns observed for these atoms, which `disappear' at low temperature. It also occurred that the potassium cation, quite uniformly coordinated by seven O atoms from all molecular fragments of the UMPH- anion, including the O atom from the ribofuranose ring, can be treated as spherical in the charge density model which was supported by theoretical calculations. Apart from the predominant electrostatic interactions, four relatively strong hydrogen bond types further support the stability of the crystal structure. This results in a compact and quite uniform structure (in all directions) of the studied crystal, as opposed to similar cases with layered architecture reported in the literature.

Interplay of point multipole moments and charge penetration for intermolecular electrostatic interaction energies from the University at Buffalo pseudoatom databank model of electron density.

The strength of the University at Buffalo DataBank (UBDB) in Ees estimation is mainly due to charge overlap effects because the UBDB offers continuous representation of charge density which allows for a direct account of charge penetration in the derivation of electrostatic energies. In the UBDB model, these effects begin to play an important role at distances below twice the equilibrium distance and significantly increase as distances decrease. At equilibrium distances they are responsible for 30-50% of Ees for polar molecules and around 90% of Ees for nonpolar molecules. When the energy estimation from the UBDB is reduced to point multipoles, the results are comparable to point charges fitted to electrostatic potentials. On the other hand, particular components of energy from point multipole moments from the UBDB model are sensitive to the type of interaction and might be helpful in the characterization of interactions.

Hydrogen atoms can be located accurately and precisely by x-ray crystallography.

Precise and accurate structural information on hydrogen atoms is crucial to the study of energies of interactions important for crystal engineering, materials science, medicine, and pharmacy, and to the estimation of physical and chemical properties in solids. However, hydrogen atoms only scatter x-radiation weakly, so x-rays have not been used routinely to locate them accurately. Textbooks and teaching classes still emphasize that hydrogen atoms cannot be located with x-rays close to heavy elements; instead, neutron diffraction is needed. We show that, contrary to widespread expectation, hydrogen atoms can be located very accurately using x-ray diffraction, yielding bond lengths involving hydrogen atoms (A-H) that are in agreement with results from neutron diffraction mostly within a single standard deviation. The precision of the determination is also comparable between x-ray and neutron diffraction results. This has been achieved at resolutions as low as 0.8 Å using Hirshfeld atom refinement (HAR). We have applied HAR to 81 crystal structures of organic molecules and compared the A-H bond lengths with those from neutron measurements for A-H bonds sorted into bonds of the same class. We further show in a selection of inorganic compounds that hydrogen atoms can be located in bridging positions and close to heavy transition metals accurately and precisely. We anticipate that, in the future, conventional x-radiation sources at in-house diffractometers can be used routinely for locating hydrogen atoms in small molecules accurately instead of large-scale facilities such as spallation sources or nuclear reactors.

A Universal and Straightforward Approach to Include Penetration Effects in Electrostatic Interaction Energy Estimation.

To compensate for the lack of the explicit treatment of charge penetration in classical force fields, we propose a new charge-distribution model based on a promolecule augmented with point charges (aug-PROmol). It relies on a superposition of spherical atomic electron densities obtained for each chemical element from SCF energy optimized atomic orbitals. Atomic densities are further rescaled by partial point charges computed from fits to the molecular electrostatic potential. Aug-PROmol was tested on the S66 benchmark dataset extended to nonequilibrium geometries (J. Chem. Theory Comput., 2011, 7, 3466). The model does not need any additional parametrization other than point charges. Despite its simplicity, aug-PROmol approximates the electrostatic energy with good agreement (RMSE=0.76 kcal mol(-1) to DFT-SAPT with B3LYP/aug-cc-pVTZ).

The influence of fluorine position on the properties of fluorobenzoxaboroles.

5-Fluoro-2,1-benzoxaborol-1(3H)-ol, a potent antifungal drug also known as Tavaborole or AN2690, has been compared with its three isomers in terms of its activity against several fungi as well as pKa and multinuclear NMR characterization. The molecular and crystal structure of 6-fluoro-2,1-benzoxaborol-1(3H)-ol was determined and compared with that of AN2690.

Electrostatic interactions in aminoglycoside-RNA complexes.

Electrostatic interactions often play key roles in the recognition of small molecules by nucleic acids. An example is aminoglycoside antibiotics, which by binding to ribosomal RNA (rRNA) affect bacterial protein synthesis. These antibiotics remain one of the few valid treatments against hospital-acquired infections by Gram-negative bacteria. It is necessary to understand the amplitude of electrostatic interactions between aminoglycosides and their rRNA targets to introduce aminoglycoside modifications that would enhance their binding or to design new scaffolds. Here, we calculated the electrostatic energy of interactions and its per-ring contributions between aminoglycosides and their primary rRNA binding site. We applied either the methodology based on the exact potential multipole moment (EPMM) or classical molecular mechanics force field single-point partial charges with Coulomb formula. For EPMM, we first reconstructed the aspherical electron density of 12 aminoglycoside-RNA complexes from the atomic parameters deposited in the University at Buffalo Databank. The University at Buffalo Databank concept assumes transferability of electron density between atoms in chemically equivalent vicinities and allows reconstruction of the electron densities from experimental structural data. From the electron density, we then calculated the electrostatic energy of interaction using EPMM. Finally, we compared the two approaches. The calculated electrostatic interaction energies between various aminoglycosides and their binding sites correlate with experimentally obtained binding free energies. Based on the calculated energetic contributions of water molecules mediating the interactions between the antibiotic and rRNA, we suggest possible modifications that could enhance aminoglycoside binding affinity.

Hirshfeld atom refinement for modelling strong hydrogen bonds.

High-resolution low-temperature synchrotron X-ray diffraction data of the salt L-phenylalaninium hydrogen maleate are used to test the new automated iterative Hirshfeld atom refinement (HAR) procedure for the modelling of strong hydrogen bonds. The HAR models used present the first examples of Z' > 1 treatments in the framework of wavefunction-based refinement methods. L-Phenylalaninium hydrogen maleate exhibits several hydrogen bonds in its crystal structure, of which the shortest and the most challenging to model is the O-H...O intramolecular hydrogen bond present in the hydrogen maleate anion (O...O distance is about 2.41 Å). In particular, the reconstruction of the electron density in the hydrogen maleate moiety and the determination of hydrogen-atom properties [positions, bond distances and anisotropic displacement parameters (ADPs)] are the focus of the study. For comparison to the HAR results, different spherical (independent atom model, IAM) and aspherical (free multipole model, MM; transferable aspherical atom model, TAAM) X-ray refinement techniques as well as results from a low-temperature neutron-diffraction experiment are employed. Hydrogen-atom ADPs are furthermore compared to those derived from a TLS/rigid-body (SHADE) treatment of the X-ray structures. The reference neutron-diffraction experiment reveals a truly symmetric hydrogen bond in the hydrogen maleate anion. Only with HAR is it possible to freely refine hydrogen-atom positions and ADPs from the X-ray data, which leads to the best electron-density model and the closest agreement with the structural parameters derived from the neutron-diffraction experiment, e.g. the symmetric hydrogen position can be reproduced. The multipole-based refinement techniques (MM and TAAM) yield slightly asymmetric positions, whereas the IAM yields a significantly asymmetric position.

Sunitinib: from charge-density studies to interaction with proteins.

Protein kinases are targets for the treatment of a number of diseases. Sunitinib malate is a type I inhibitor of tyrosine kinases and was approved as a drug in 2006. This contribution constitutes the first comprehensive analysis of the crystal structures of sunitinib malate and of complexes of sunitinib with a series of protein kinases. The high-resolution single-crystal X-ray measurement and aspherical atom databank approach served as a basis for reconstruction of the charge-density distribution of sunitinib and its protein complexes. Hirshfeld surface and topological analyses revealed a similar interaction pattern in the sunitinib malate crystal structure to that in the protein binding pockets. Sunitinib forms nine preserved bond paths corresponding to hydrogen bonds and also to the C-H···O and C-H···π contacts common to the VEGRF2, CDK2, G2, KIT and IT kinases. In general, sunitinib interacts with the studied proteins with a similar electrostatic interaction energy and can adjust its conformation to fit the binding pocket in such a way as to enhance the electrostatic interactions, e.g. hydrogen bonds in ligand-kinase complexes. Such behaviour may be responsible for the broad spectrum of action of sunitinib as a kinase inhibitor.

A Comparative Study of Transferable Aspherical Pseudoatom Databank and Classical Force Fields for Predicting Electrostatic Interactions in Molecular Dimers.

Accurate and fast evaluation of electrostatic interactions in molecular systems is one of the most challenging tasks in the rapidly advancing field of macromolecular chemistry and drug design. Electrostatic interactions are of crucial importance in biological systems. They are well represented by quantum mechanical methods; however, such calculations are computationally expensive. In this study, we have evaluated the University of Buffalo Pseudoatom Databank (UBDB)1,2 approach for approximation of electrostatic properties of macromolecules and their complexes. We selected the S663 and JSCH-20054 data sets (208 molecular complexes in total) for this study. These complexes represent a wide range of chemical and biological systems for which hydrogen bonding, electrostatic, and van der Waals interactions play important roles. Reference electrostatic energies were obtained directly from wave functions at the B3LYP/aug-cc-pVTZ level of theory using the SAPT (Symmetry-Adapted Perturbation Theory) scheme for calculation of electrostatic contributions to total intermolecular interaction energies. Electrostatic energies calculated on the basis of the UBDB were compared with corresponding reference results. Results were also compared with energies computed using a point charge model from popular force fields (AM1-BCC and RESP used in AMBER and CGenFF from CHARMM family). The energy trends are quite consistent (R2 ≈ 0.98) for the UBDB method as compared to the AMBER5 and CHARMM force field methods6(R2 ≈ 0.93 on average). The RSMEs do not exceed 3.2 kcal mol-1 for the UBDB and are in the range of 3.7-7.6 kcal mol-1 for the point charge models. We also investigated the discrepancies in electrostatic potentials and magnitudes of dipole moments among the tested methods. This study shows that estimation of electrostatic interaction energies using the UBDB databank is accurate and reasonably fast when compared to other known methods, which opens potential new applications to macromolecules.

Transferability of atomic multipoles in amino acids and peptides for various density partitions.

Multipole expansion of electron density distribution is an efficient tool for evaluating the energy of interactions in molecules. In as much as atoms in macromolecules such as peptides are modeled with certain types of atoms derived from small organic molecules, investigating transferability of atomic multipoles for various partitions of molecular electron density is an important issue. In this study, multipole moments up to hexadecapoles for types of atoms present in selected amino acids, as well as di- and tripeptides composed of these amino acids, are computed using three density partitions: Hansen-Coppens aspherical pseudoatoms formalism, Hirshfeld's stockholder partition, and Bader's atoms in molecules theory. Electron density of relevant molecules is derived in a procedure including molecular wave function ab initio calculation for isolated molecules in geometry from X-ray measurements, calculation of theoretical structure factors for molecules put in a pseudocubic cell, and multipole refinement as in crystallographic data processing and computation of multipoles. The results were compared to calculations of multipole moments in AIM and in stockholder density partitions obtained directly from molecular wave functions. The presented comparison does not point unambiguously to any particular influence of multipole refinement on moments obtained from these two partitions. The advantage of stockholder partitioning in terms of transferability of atomic multipoles is affirmed. AIM and pseudoatoms provide slightly less transferable multipoles of lower order. Higher rank of multipole expansion reveals a transferability improvement in the case of AIM and meaningful deterioration for pseudoatoms.