1
|
Olech B, Brázda P, Palatinus L, Dominiak PM. Dynamical refinement with multipolar electron scattering factors. IUCRJ 2024; 11:309-324. [PMID: 38512772 PMCID: PMC11067749 DOI: 10.1107/s2052252524001763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 02/22/2024] [Indexed: 03/23/2024]
Abstract
Dynamical refinement is a well established method for refining crystal structures against 3D electron diffraction (ED) data and its benefits have been discussed in the literature [Palatinus, Petříček & Corrêa, (2015). Acta Cryst. A71, 235-244; Palatinus, Corrêa et al. (2015). Acta Cryst. B71, 740-751]. However, until now, dynamical refinements have only been conducted using the independent atom model (IAM). Recent research has shown that a more accurate description can be achieved by applying the transferable aspherical atom model (TAAM), but this has been limited only to kinematical refinements [Gruza et al. (2020). Acta Cryst. A76, 92-109; Jha et al. (2021). J. Appl. Cryst. 54, 1234-1243]. In this study, we combine dynamical refinement with TAAM for the crystal structure of 1-methyluracil, using data from precession ED. Our results show that this approach improves the residual Fourier electrostatic potential and refinement figures of merit. Furthermore, it leads to systematic changes in the atomic displacement parameters of all atoms and the positions of hydrogen atoms. We found that the refinement results are sensitive to the parameters used in the TAAM modelling process. Though our results show that TAAM offers superior performance compared with IAM in all cases, they also show that TAAM parameters obtained by periodic DFT calculations on the refined structure are superior to the TAAM parameters from the UBDB/MATTS database. It appears that multipolar parameters transferred from the database may not be sufficiently accurate to provide a satisfactory description of all details of the electrostatic potential probed by the 3D ED experiment.
Collapse
Affiliation(s)
- Barbara Olech
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Warsaw, Poland
| | - Petr Brázda
- Institute of Physics of the Czech Academy of Sciences, Na Slovance 1999/2, 182 00 Prague, Czechia
| | - Lukas Palatinus
- Institute of Physics of the Czech Academy of Sciences, Na Slovance 1999/2, 182 00 Prague, Czechia
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Warsaw, Poland
| |
Collapse
|
2
|
Chodkiewicz M, Patrikeev L, Pawlędzio S, Woźniak K. Transferable Hirshfeld atom model for rapid evaluation of aspherical atomic form factors. IUCRJ 2024; 11:249-259. [PMID: 38446457 PMCID: PMC10916294 DOI: 10.1107/s2052252524001507] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 02/14/2024] [Indexed: 03/07/2024]
Abstract
Form factors based on aspherical models of atomic electron density have brought great improvement in the accuracies of hydrogen atom parameters derived from X-ray crystal structure refinement. Today, two main groups of such models are available, the banks of transferable atomic densities parametrized using the Hansen-Coppens multipole model which allows for rapid evaluation of atomic form factors and Hirshfeld atom refinement (HAR)-related methods which are usually more accurate but also slower. In this work, a model that combines the ideas utilized in the two approaches is tested. It uses atomic electron densities based on Hirshfeld partitions of electron densities, which are precalculated and stored in a databank. This model was also applied during the refinement of the structures of five small molecules. A comparison of the resulting hydrogen atom parameters with those derived from neutron diffraction data indicates that they are more accurate than those obtained with the Hansen-Coppens based databank, and only slightly less accurate than those obtained with a version of HAR that neglects the crystal environment. The advantage of using HAR becomes more noticeable when the effects of the environment are included. To speed up calculations, atomic densities were represented by multipole expansion with spherical harmonics up to l = 7, which used numerical radial functions (a different approach to that applied in the Hansen-Coppens model). Calculations of atomic form factors for the small protein crambin (at 0.73 Å resolution) took only 68 s using 12 CPU cores.
Collapse
Affiliation(s)
- Michał Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Leonid Patrikeev
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Sylwia Pawlędzio
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| |
Collapse
|
3
|
Saifina AF, Kartashov SV, Saifina LF, Fayzullin RR. 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. IUCRJ 2023; 10:584-602. [PMID: 37668216 PMCID: PMC10478519 DOI: 10.1107/s2052252523007108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Accepted: 08/11/2023] [Indexed: 09/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Alina F. Saifina
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| | - Sergey V. Kartashov
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| | - Liliya F. Saifina
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| | - Robert R. Fayzullin
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| |
Collapse
|
4
|
Jha KK, Kleemiss F, Chodkiewicz ML, Dominiak PM. Aspherical atom refinements on X-ray data of diverse structures including disordered and covalent organic framework systems: a time-accuracy trade-off. J Appl Crystallogr 2023; 56:116-127. [PMID: 36777135 PMCID: PMC9901929 DOI: 10.1107/s1600576722010883] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Accepted: 11/13/2022] [Indexed: 12/24/2022] Open
Abstract
Aspherical atom refinement is the key to achieving accurate structure models, displacement parameters, hydrogen-bond lengths and analysis of weak interactions, amongst other examples. There are various quantum crystallographic methods to perform aspherical atom refinement, including Hirshfeld atom refinement (HAR) and transferable aspherical atom model (TAAM) refinement. Both HAR and TAAM have their limitations and advantages, the former being more accurate and the latter being faster. With the advent of non-spherical atoms in Olex2 (NoSpherA2), it is now possible to overcome some limitations, like treating disorder, twinning and network structures, in aspherical refinements using HAR, TAAM or both together. TAAM refinement in NoSpherA2 showed significant improvement in refinement statistics compared with independent atom model (IAM) refinements on a diverse set of X-ray diffraction data. The sensitivity of TAAM towards poor data quality and disorder was observed in terms of higher refinement statistics for such structures. A comparison of IAM with TAAM and HAR in NoSpherA2 indicated that the time taken by TAAM refinements was of the same order of magnitude as that taken by IAM, while in HAR the time taken using a minimal basis set was 50 times higher than for IAM and rapidly increased with increasing size of the basis sets used. The displacement parameters for hydrogen and non-hydrogen atoms were very similar in both HAR and TAAM refinements. The hydrogen-bond lengths were slightly closer to neutron reference values in the case of HAR with higher basis sets than in TAAM. To benefit from the advantages of each method, a new hybrid refinement approach has been introduced, allowing a combination of IAM, HAR and TAAM in one structure refinement. Refinement of coordination complexes involving metal-organic compounds and network structures such as covalent organic frameworks and metal-organic frameworks is now possible in a hybrid mode such as IAM-TAAM or HAR-TAAM, where the metal atoms are treated via either the IAM or HAR method and the organic part via TAAM, thus reducing the computational costs without compromising the accuracy. Formal charges on the metal and ligand can also be introduced in hybrid-mode refinement.
Collapse
Affiliation(s)
- Kunal Kumar Jha
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Florian Kleemiss
- Fakultät für Chemie und Pharmazie, Universität Regensburg, Universitätstrasse 31, Regensburg, Bayern 93053, Germany
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| |
Collapse
|
5
|
Budniak U, Karolak NK, Kulik M, Młynarczyk K, Górna MW, Dominiak PM. The Role of Electrostatic Interactions in IFIT5-RNA Complexes Predicted by the UBDB+EPMM Method. J Phys Chem B 2022; 126:9152-9167. [PMID: 36326196 PMCID: PMC9677429 DOI: 10.1021/acs.jpcb.2c04519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Electrostatic energy has a significant contribution to intermolecular interaction energy, especially in biological systems. Unfortunately, precise quantum mechanics calculations are not feasible for large biological systems; hence, simpler calculation methods are required. We propose a method called UBDB+EPMM (University at Buffalo Pseudoatom DataBank + Exact Potential Multipole Moments), which shortens computational time without losing accuracy. Here, we characterize electrostatic interactions in selected complexes of IFIT proteins with RNA. IFIT proteins are effectors of the innate immune system, and by binding foreign RNA, they prevent the synthesis of viral proteins in human host cells; hence, they block the propagation of viruses. We show that by using the UBDB+EPMM method it is possible to describe protein-RNA interactions not only qualitatively but also quantitatively. Looking at the charge penetration contribution to electrostatic interactions, we find all amino acid residues with strong local interactions. Moreover, we confirm that electrostatic interaction of IFIT5 with pppRNA does not depend on the sequence of the RNA.
Collapse
Affiliation(s)
- Urszula
Anna Budniak
- Biological
and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Natalia Katarzyna Karolak
- Biological
and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland,Nencki
Institute of Experimental Biology, Polish
Academy of Sciences, ul. Ludwika Pasteura 3, 02-093 Warszawa, Poland
| | - Marta Kulik
- Biological
and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Krzysztof Młynarczyk
- Biological
and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Maria Wiktoria Górna
- Biological
and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland,E-mail:
| | - Paulina Maria Dominiak
- Biological
and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland,E-mail:
| |
Collapse
|
6
|
Electron density is not spherical: the many applications of the transferable aspherical atom model. Comput Struct Biotechnol J 2022; 20:6237-6243. [DOI: 10.1016/j.csbj.2022.10.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2022] [Revised: 10/12/2022] [Accepted: 10/13/2022] [Indexed: 11/20/2022] Open
|
7
|
Jha KK, Gruza B, Sypko A, Kumar P, Chodkiewicz ML, Dominiak PM. Multipolar Atom Types from Theory and Statistical Clustering (MATTS) Data Bank: Restructurization and Extension of UBDB. J Chem Inf Model 2022; 62:3752-3765. [PMID: 35943747 PMCID: PMC9400107 DOI: 10.1021/acs.jcim.2c00144] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A fast and accurate operational model of electron density
is crucial
in many scientific disciplines including crystallography, molecular
biology, pharmaceutical, and structural chemistry. In quantum crystallography,
the aspherical refinement of crystal structures is becoming increasingly
popular because of its accurate description in terms of physically
meaningful properties. The transferable aspherical atom model (TAAM)
is quick and precise, though it requires a robust algorithm for atom
typing and coverage of the most popular atom types present in small
organic molecules. Thus, the University at Buffalo Databank (UBDB)
has been renamed to the Multipolar Atom Types from Theory and Statistical
clustering (MATTS) data bank, broadened, restructured, and implemented
into the software DiSCaMB with 651 atom types obtained from 2316 small-molecule
crystal structures containing C, H, N, O, P, S, F, Cl, and Br atoms.
MATTS2021 data bank now covers most of the small molecules, peptides,
RNA, DNA, and some frequently occurring cations and anions in biological,
pharmaceutical, and organic materials, including the majority of known
crystal structures composed of the above elements. The multipole model
parameters (Pval, κ, κ′, Plm) obtained for different
atom types were greatly influenced by neighboring atom types, hybridization,
geometrical strain in the ring system, and charges on the molecule.
Contrary to previous findings, the atoms showing variable oxidation
states and ions deviate from the linear dependence of monopole-derived
charges on the expansion–contraction κ parameter.
Collapse
Affiliation(s)
- Kunal Kumar Jha
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Barbara Gruza
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Aleksandra Sypko
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Prashant Kumar
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|
8
|
Rybicka PM, Kulik M, Chodkiewicz ML, Dominiak PM. Multipolar Atom Types from Theory and Statistical Clustering (MATTS) Data Bank: Impact of Surrounding Atoms on Electron Density from Cluster Analysis. J Chem Inf Model 2022; 62:3766-3783. [PMID: 35943739 PMCID: PMC9400106 DOI: 10.1021/acs.jcim.2c00145] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The multipole model (MM) uses an aspherical approach
to describe
electron density and can be used to interpret data from X-ray diffraction
in a more accurate manner than using the spherical approximation.
The MATTS (multipolar atom types from theory and statistical clustering)
data bank gathers MM parameters specific for atom types in proteins,
nucleic acids, and organic molecules. However, it was not fully understood
how the electron density of particular atoms responds to their surroundings
and which factors describe the electron density in molecules within
the MM. In this work, by applying clustering using descriptors available
in the MATTS data bank, that is, topology and multipole parameters,
we found the topology features with the biggest impact on the multipole
parameters: the element of the central atom, the number of first neighbors,
and planarity of the group. The similarities in the spatial distribution
of electron density between and within atom type classes revealed
distinct and unique atom types. The quality of existing types can
be improved by adding better parametrization, definitions, and local
coordinate systems. Future development of the MATTS data bank should
lead to a wider range of atom types necessary to construct the electron
density of any molecule.
Collapse
Affiliation(s)
- Paulina Maria Rybicka
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Marta Kulik
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|
9
|
Kulik M, Chodkiewicz ML, Dominiak PM. Theoretical 3D electron diffraction electrostatic potential maps of proteins modeled with a multipolar pseudoatom data bank. Acta Crystallogr D Struct Biol 2022; 78:1010-1020. [PMID: 35916225 PMCID: PMC9344478 DOI: 10.1107/s2059798322005836] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2021] [Accepted: 05/31/2022] [Indexed: 11/10/2022] Open
Abstract
The availability of atomic resolution experimental maps of electrostatic potential from 3D electron diffraction (3D ED) extends the possibility of investigating the electrostatic potential beyond the determination of non-H-atom positions. However, accurate tools to calculate this potential for macromolecules, without the use of expensive quantum calculations, are lacking. The University at Buffalo Data Bank (UBDB) gathers atom types that can be used to calculate accurate electrostatic potential maps via structure-factor calculations. Here, the transferable aspherical atom model (TAAM) is applied with UBDB to investigate theoretically obtained electrostatic potential maps of lysozyme and proteinase K, and compare them with experimental maps from 3D ED. UBDB better reproduces the molecular electrostatic potential of molecules within their entire volume compared with the neutral spherical models used in the popular independent atom model (IAM). Additionally, the theoretical electron-density maps of the studied proteins are shown and compared with the electrostatic potential maps. The atomic displacement parameters (B factors) may affect the electrostatic potential maps in a different way than in the case of electron-density maps. The computational method presented in this study could potentially facilitate the interpretation of the less resolved regions of cryo-electron microscopy density maps and pave the way for distinguishing between different ions/water molecules in the active sites of macromolecules in high-resolution structures, which is of interest for drug-design purposes.
Collapse
Affiliation(s)
- Marta Kulik
- Biological and Chemical Research Centre, Faculty of Chemistry, University of Warsaw, Zwirki i Wigury 101, 02-089 Warsaw, Poland
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Faculty of Chemistry, University of Warsaw, Zwirki i Wigury 101, 02-089 Warsaw, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Faculty of Chemistry, University of Warsaw, Zwirki i Wigury 101, 02-089 Warsaw, Poland
| |
Collapse
|
10
|
Chodkiewicz M, Pawlędzio S, Woińska M, Woźniak K. Fragmentation and transferability in Hirshfeld atom refinement. IUCRJ 2022; 9:298-315. [PMID: 35371499 PMCID: PMC8895009 DOI: 10.1107/s2052252522000690] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2021] [Accepted: 01/19/2022] [Indexed: 05/06/2023]
Abstract
Hirshfeld atom refinement (HAR) is one of the most effective methods for obtaining accurate structural parameters for hydrogen atoms from X-ray diffraction data. Unfortunately, it is also relatively computationally expensive, especially for larger molecules due to wavefunction calculations. Here, a fragmentation approach has been tested as a remedy for this problem. It gives an order of magnitude improvement in computation time for larger organic systems and is a few times faster for metal-organic systems at the cost of only minor differences in the calculated structural parameters when compared with the original HAR calculations. Fragmentation was also applied to polymeric and disordered systems where it provides a natural solution to problems that arise when HAR is applied. The concept of fragmentation is closely related to the transferable aspherical atom model (TAAM) and allows insight into possible ways to improve TAAM. Hybrid approaches combining fragmentation with the transfer of atomic densities between chemically similar atoms have been tested. An efficient handling of intermolecular interactions was also introduced for calculations involving fragmentation. When applied in fragHAR (a fragmentation approach for polypeptides) as a replacement for the original approach, it allowed for more efficient calculations. All of the calculations were performed with a locally modified version of Olex2 combined with a development version of discamb2tsc and ORCA. Care was taken to efficiently use the power of multicore processors by simple implementation of load-balancing, which was found to be very important for lowering computational time.
Collapse
Affiliation(s)
- Michał Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
- Correspondence e-mail: ,
| | - Sylwia Pawlędzio
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Magdalena Woińska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
- Correspondence e-mail: ,
| |
Collapse
|
11
|
Vuković V, Leduc T, Jelić-Matošević Z, Didierjean C, Favier F, Guillot B, Jelsch C. A rush to explore protein-ligand electrostatic interaction energy with Charger. Acta Crystallogr D Struct Biol 2021; 77:1292-1304. [PMID: 34605432 DOI: 10.1107/s2059798321008433] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 08/13/2021] [Indexed: 11/10/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Vedran Vuković
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
| | - Theo Leduc
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
| | - Zoe Jelić-Matošević
- Department of Chemistry, Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia
| | | | | | - Benoît Guillot
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
| | | |
Collapse
|
12
|
Weatherly J, Macchi P, Volkov A. On the calculation of the electrostatic potential, electric field and electric field gradient from the aspherical pseudoatom model. II. Evaluation of the properties in an infinite crystal. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2021; 77:399-419. [PMID: 34473095 DOI: 10.1107/s2053273321005532] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Accepted: 05/28/2021] [Indexed: 11/10/2022]
Abstract
The previously reported exact potential and multipole moment (EP/MM) method for fast and precise evaluation of the intermolecular electrostatic interaction energies in molecular crystals using the pseudoatom representation of the electron density [Nguyen, Macchi & Volkov (2020), Acta Cryst. A76, 630-651] has been extended to the calculation of the electrostatic potential (ESP), electric field (EF) and electric field gradient (EFG) in an infinite crystal. The presented approach combines an efficient Ewald-type summation (ES) of atomic multipoles up to the hexadecapolar level in direct and reciprocal spaces with corrections for (i) the net polarization of the sample (the `surface term') due to a net dipole moment of the crystallographic unit cell (if present) and (ii) the short-range electron-density penetration effects. The rederived and reported closed-form expressions for all terms in the ES algorithm have been augmented by the expressions for the surface term available in the literature [Stenhammar, Trulsson & Linse (2011), J. Chem. Phys. 134, 224104] and the exact potential expressions reported in a previous study [Volkov, King, Coppens & Farrugia (2006), Acta Cryst. A62, 400-408]. The resulting algorithm, coded using Fortran in the XDPROP module of the software package XD, was tested on several small molecular crystal systems (formamide, benzene, L-dopa, paracetamol, amino acids etc.) and compared with a series of EP/MM-based direct-space summations (DS) performed within a certain number of unit cells generated along both the positive and negative crystallographic directions. The EP/MM-based ES technique allows for a noticeably more precise determination of the EF and EFG and significantly better precision of the evaluated ESP when compared with the DS calculations, even when the latter include contributions from an array of symmetry-equivalent atoms generated within four additional unit cells along each crystallographic direction. In terms of computational performance, the ES/EP/MM method is significantly faster than the DS calculations performed within the extended unit-cell limits but trails the DS calculations within the reduced summation ranges. Nonetheless, the described EP/MM-based ES algorithm is superior to the direct-space summations as it does not require the user to monitor continuously the convergence of the evaluated properties as a function of the summation limits and offers a better precision-performance balance.
Collapse
Affiliation(s)
- Jessie Weatherly
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| | - Piero Macchi
- Department of Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Via Mancinelli 7, Milano, 20131, Italy
| | - Anatoliy Volkov
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| |
Collapse
|
13
|
Jha KK, Gruza B, Chodkiewicz ML, Jelsch C, Dominiak PM. Refinements on electron diffraction data of β-glycine in MoPro: a quest for an improved structure model. J Appl Crystallogr 2021. [DOI: 10.1107/s160057672100580x] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
The advancement in 3D electron diffraction (3D ED) techniques that lead to a revolution in molecular structure determination using nano-sized crystals is now achieving atomic resolution. The structures can be obtained from 3D ED data with tools similar to those used for X-ray structure determination. In this context, the MoPro software, originally designed for structure and charge density refinements using X-ray diffraction data, has been adapted. Structure refinement on 3D ED data was achieved via implementation of electron scattering factors available in the literature and by application of the Mott–Bethe equation to X-ray scattering factors computed from the multipolar atom model. The multipolar model was parametrized using the transferable pseudoatom databanks ELMAM2 and UBDB. Applying the independent atom model (IAM), i.e. spherical neutral atom refinement, to 3D ED data on β-glycine in MoPro resulted in structure and refinement statistics comparable to those obtained from other well known software. Use of the transferred aspherical atom model (TAAM) led to improvement of the refinement statistics and a better fit of the model to the 3D ED data as compared with the spherical atom refinement. The anisotropic displacement parameters of non-H atoms appear underestimated by typically 0.003 Å2 for the non-H atoms in IAM refinement compared with TAAM. Thus, MoPro is shown to be an effective tool for crystal structure refinement on 3D ED data and allows use of a spherical or a multipolar atom model. Electron density databases can be readily transferred with no further modification needed when the Mott–Bethe equation is applied.
Collapse
|
14
|
Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease. Molecules 2021; 26:molecules26133872. [PMID: 34202892 PMCID: PMC8270314 DOI: 10.3390/molecules26133872] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Revised: 06/19/2021] [Accepted: 06/22/2021] [Indexed: 11/23/2022] Open
Abstract
Computational analysis of protein–ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions, such as in docking programs, all rely on equations that sum each type of protein–ligand interactions in order to predict the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions constitute one of the most important part of total interactions between macromolecules. Unlike dispersion forces, they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this study, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and darunavir) were analyzed by using charge densities from the transferable aspherical-atom University at Buffalo Databank (UBDB). Moreover, we analyzed the electrostatic interaction energy for an ensemble of structures, using molecular dynamic simulations to highlight the main features of electrostatic interactions important for binding affinity.
Collapse
|
15
|
Shteingolts SA, Saifina AF, Saifina LF, Semenov VE, Fukin GK, Fayzullin RR. X-ray charge density study of the 6-methyluracil derivative in the crystal: Revealing, consequences, and multipole refinement of minor static disorder. J Mol Struct 2021. [DOI: 10.1016/j.molstruc.2020.129724] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
16
|
Ikemoto K, Isobe H. Geodesic Phenine Frameworks. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2021. [DOI: 10.1246/bcsj.20200284] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Koki Ikemoto
- Department of Chemistry, The University of Tokyo, Hongo, Tokyo 113-0033, Japan
| | - Hiroyuki Isobe
- Department of Chemistry, The University of Tokyo, Hongo, Tokyo 113-0033, Japan
| |
Collapse
|
17
|
Affiliation(s)
- Piero Macchi
- Department, Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Milano, Italy
- Center for Nano Science and Technology CNST@polimi, Italian Institute of Technology, Milano, Italy
| |
Collapse
|
18
|
Chodkiewicz ML, Woińska M, Woźniak K. Hirshfeld atom like refinement with alternative electron density partitions. IUCRJ 2020; 7:1199-1215. [PMID: 33209330 PMCID: PMC7642787 DOI: 10.1107/s2052252520013603] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Accepted: 10/12/2020] [Indexed: 05/03/2023]
Abstract
Hirshfeld atom refinement is one of the most successful methods for the accurate determination of structural parameters for hydrogen atoms from X-ray diffraction data. This work introduces a generalization of the method [generalized atom refinement (GAR)], consisting of the application of various methods of partitioning electron density into atomic contributions. These were tested on three organic structures using the following partitions: Hirshfeld, iterative Hirshfeld, iterative stockholder, minimal basis iterative stockholder and Becke. The effects of partition choice were also compared with those caused by other factors such as quantum chemical methodology, basis set, representation of the crystal field and a combination of these factors. The differences between the partitions were small in terms of R factor (e.g. much smaller than for refinements with different quantum chemistry methods, i.e. Hartree-Fock and coupled cluster) and therefore no single partition was clearly the best in terms of experimental data reconstruction. In the case of structural parameters the differences between the partitions are comparable to those related to the choice of other factors. We have observed the systematic effects of the partition choice on bond lengths and ADP values of polar hydrogen atoms. The bond lengths were also systematically influenced by the choice of electron density calculation methodology. This suggests that GAR-derived structural parameters could be systematically improved by selecting an optimal combination of the partition and quantum chemistry method. The results of the refinements were compared with those of neutron diffraction experiments. This allowed a selection of the most promising partition methods for further optimization of GAR settings, namely the Hirshfeld, iterative stockholder and minimal basis iterative stockholder.
Collapse
Affiliation(s)
- Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa, 02-089 Warszawa, Poland
| | - Magdalena Woińska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa, 02-089 Warszawa, Poland
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa, 02-089 Warszawa, Poland
| |
Collapse
|
19
|
Nguyen D, Macchi P, Volkov A. Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. III. Application to crystal structures via the Ewald and direct summation methods. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2020; 76:630-651. [PMID: 33125348 DOI: 10.1107/s2053273320009584] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2020] [Accepted: 07/13/2020] [Indexed: 12/28/2022]
Abstract
The previously reported exact potential and multipole moment (EP/MM) method for fast and accurate evaluation of the intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density [Volkov, Koritsanszky & Coppens (2004). Chem. Phys. Lett. 391, 170-175; Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536; Nguyen & Volkov (2019). Acta Cryst. A75, 448-464] is extended to the calculation of electrostatic interaction energies in molecular crystals using two newly developed implementations: (i) the Ewald summation (ES), which includes interactions up to the hexadecapolar level and the EP correction to account for short-range electron-density penetration effects, and (ii) the enhanced EP/MM-based direct summation (DS), which at sufficiently large intermolecular separations replaces the atomic multipole moment approximation to the electrostatic energy with that based on the molecular multipole moments. As in the previous study [Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536], the EP electron repulsion integral is evaluated analytically using the Löwdin α-function approach. The resulting techniques, incorporated in the XDPROP module of the software package XD2016, have been tested on several small-molecule crystal systems (benzene, L-dopa, paracetamol, amino acids etc.) and the crystal structure of a 181-atom decapeptide molecule (Z = 4) using electron densities constructed via the University at Buffalo Aspherical Pseudoatom Databank [Volkov, Li, Koritsanszky & Coppens (2004). J. Phys. Chem. A, 108, 4283-4300]. Using a 2015 2.8 GHz Intel Xeon E3-1505M v5 computer processor, a 64-bit implementation of the Löwdin α-function and one of the higher optimization levels in the GNU Fortran compiler, the ES method evaluates the electrostatic interaction energy with a numerical precision of at least 10-5 kJ mol-1 in under 6 s for any of the tested small-molecule crystal structures, and in 48.5 s for the decapeptide structure. The DS approach is competitive in terms of precision and speed with the ES technique only for crystal structures of small molecules that do not carry a large molecular dipole moment. The electron-density penetration effects, correctly accounted for by the two described methods, contribute 28-64% to the total electrostatic interaction energy in the examined systems, and thus cannot be neglected.
Collapse
Affiliation(s)
- Daniel Nguyen
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| | - Piero Macchi
- Department of Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Via Mancinelli 7, Milano 20131, Italy
| | - Anatoliy Volkov
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| |
Collapse
|
20
|
|
21
|
Sanjuan-Szklarz WF, Woińska M, Domagała S, Dominiak PM, Grabowsky S, Jayatilaka D, Gutmann M, Woźniak K. On the accuracy and precision of X-ray and neutron diffraction results as a function of resolution and the electron density model. IUCRJ 2020; 7:920-933. [PMID: 32939284 PMCID: PMC7467170 DOI: 10.1107/s2052252520010441] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Accepted: 07/28/2020] [Indexed: 05/20/2023]
Abstract
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.
Collapse
Affiliation(s)
- W. Fabiola Sanjuan-Szklarz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| | - Magdalena Woińska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| | - Sławomir Domagała
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| | - Paulina M. Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| | - Simon Grabowsky
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland
| | - Dylan Jayatilaka
- School of Molecular Sciences, University of Western Australia, 35 Stirling Highway, Perth WA 6009, Australia
| | - Matthias Gutmann
- Rutherford Appleton Laboratory, ISIS Facility, Chilton, Didcot, Oxfordshire OX11 OQX, United Kingdom
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| |
Collapse
|
22
|
Faroque MU, Mehmood A, Noureen S, Ahmed M. Crystal engineering and electrostatic properties of co-crystals of pyrimethamine with benzoic acid and gallic acid. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2020.128183] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
|
23
|
Vénosová B, Koziskova J, Kožíšek J, Herich P, Lušpai K, Petricek V, Hartung J, Müller M, Hübschle CB, van Smaalen S, Bucinsky L. Charge density of 4-methyl-3-[(tetrahydro-2H-pyran-2-yl)oxy]thiazole-2(3H)-thione. A comprehensive multipole refinement, maximum entropy method and density functional theory study. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2020; 76:450-468. [PMID: 32831263 DOI: 10.1107/s2052520620005533] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Accepted: 04/20/2020] [Indexed: 06/11/2023]
Abstract
The structure of 4-methyl-3-[(tetrahydro-2H-pyran-2-yl)oxy]thiazole-2(3H)-thione (MTTOTHP) was investigated using X-ray diffraction and computational chemistry methods for determining properties of the nitrogen-oxygen bond, which is the least stable entity upon photochemical excitation. Experimentally measured structure factors have been used to determine and characterize charge density via the multipole model (MM) and the maximum entropy method (MEM). Theoretical investigation of the electron density and the electronic structure has been performed in the finite basis set density functional theory (DFT) framework. Quantum Theory of Atoms In Molecules (QTAIM), deformation densities and Laplacians maps have been used to compare theoretical and experimental results. MM experimental results and predictions from theory differ with respect to the sign and/or magnitude of the Laplacian at the N-O bond critical point (BCP), depending on the treatment of n values of the MM radial functions. Such Laplacian differences in the N-O bond case are discussed with respect to a lack of flexibility in the MM radial functions also reported by Rykounov et al. [Acta Cryst. (2011), B67, 425-436]. BCP Hessian eigenvalues show qualitatively matching results between MM and DFT. In addition, the theoretical analysis used domain-averaged fermi holes (DAFH), natural bond orbital (NBO) analysis and localized (LOC) orbitals to characterize the N-O bond as a single σ bond with marginal π character. Hirshfeld atom refinement (HAR) has been employed to compare to the MM refinement results and/or neutron dataset C-H bond lengths and to crystal or single molecule geometry optimizations, including considerations of anisotropy of H atoms. Our findings help to understand properties of molecules like MTTOTHP as progenitors of free oxygen radicals.
Collapse
Affiliation(s)
- Barbora Vénosová
- Institute of Physical Chemistry and Chemical Physics, Faculty of Chemical and Food Technology, Slovak Technical University of Technology in Bratislava, Radlinského 9, Bratislava, SK-81237, Slovak Republic
| | - Julia Koziskova
- Institute of Physical Chemistry and Chemical Physics, Faculty of Chemical and Food Technology, Slovak Technical University of Technology in Bratislava, Radlinského 9, Bratislava, SK-81237, Slovak Republic
| | - Jozef Kožíšek
- Institute of Physical Chemistry and Chemical Physics, Faculty of Chemical and Food Technology, Slovak Technical University of Technology in Bratislava, Radlinského 9, Bratislava, SK-81237, Slovak Republic
| | - Peter Herich
- Institute of Physical Chemistry and Chemical Physics, Faculty of Chemical and Food Technology, Slovak Technical University of Technology in Bratislava, Radlinského 9, Bratislava, SK-81237, Slovak Republic
| | - Karol Lušpai
- Institute of Physical Chemistry and Chemical Physics, Faculty of Chemical and Food Technology, Slovak Technical University of Technology in Bratislava, Radlinského 9, Bratislava, SK-81237, Slovak Republic
| | - Vaclav Petricek
- Institute of Physics, Czech Academy of Sciences, Na Slovance 1999/2, Praha 8, 182 21, Czech Republic
| | - Jens Hartung
- Fachbereich Chemie, Organische Chemie, Technische Universität Kaiserslautern, Erwin-Schrödinger-Straße, Kaiserslautern, D-67663, Germany
| | - Mike Müller
- Fachbereich Chemie, Organische Chemie, Technische Universität Kaiserslautern, Erwin-Schrödinger-Straße, Kaiserslautern, D-67663, Germany
| | - Christian B Hübschle
- Laboratory of Crystallography, University of Bayreuth, Universitätsstrasse 30, Bayreuth, 95447, Germany
| | - Sander van Smaalen
- Laboratory of Crystallography, University of Bayreuth, Universitätsstrasse 30, Bayreuth, 95447, Germany
| | - Lukas Bucinsky
- Institute of Physical Chemistry and Chemical Physics, Faculty of Chemical and Food Technology, Slovak Technical University of Technology in Bratislava, Radlinského 9, Bratislava, SK-81237, Slovak Republic
| |
Collapse
|
24
|
Jha KK, Gruza B, Kumar P, Chodkiewicz ML, Dominiak PM. TAAM: a reliable and user friendly tool for hydrogen-atom location using routine X-ray diffraction data. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2020; 76:296-306. [PMID: 32831250 DOI: 10.1107/s2052520620002917] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2019] [Accepted: 03/02/2020] [Indexed: 06/11/2023]
Abstract
Hydrogen is present in almost all of the molecules in living things. It is very reactive and forms bonds with most of the elements, terminating their valences and enhancing their chemistry. X-ray diffraction is the most common method for structure determination. It depends on scattering of X-rays from electron density, which means the single electron of hydrogen is difficult to detect. Generally, neutron diffraction data are used to determine the accurate position of hydrogen atoms. However, the requirement for good quality single crystals, costly maintenance and the limited number of neutron diffraction facilities means that these kind of results are rarely available. Here it is shown that the use of Transferable Aspherical Atom Model (TAAM) instead of Independent Atom Model (IAM) in routine structure refinement with X-ray data is another possible solution which largely improves the precision and accuracy of X-H bond lengths and makes them comparable to averaged neutron bond lengths. TAAM, built from a pseudoatom databank, was used to determine the X-H bond lengths on 75 data sets for organic molecule crystals. TAAM parametrizations available in the modified University of Buffalo Databank (UBDB) of pseudoatoms applied through the DiSCaMB software library were used. The averaged bond lengths determined by TAAM refinements with X-ray diffraction data of atomic resolution (dmin ≤ 0.83 Å) showed very good agreement with neutron data, mostly within one single sample standard deviation, much like Hirshfeld atom refinement (HAR). Atomic displacements for both hydrogen and non-hydrogen atoms obtained from the refinements systematically differed from IAM results. Overall TAAM gave better fits to experimental data of standard resolution compared to IAM. The research was accompanied with development of software aimed at providing user-friendly tools to use aspherical atom models in refinement of organic molecules at speeds comparable to routine refinements based on spherical atom model.
Collapse
Affiliation(s)
- Kunal Kumar Jha
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| | - Barbara Gruza
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| | - Prashant Kumar
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| | - Michal Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| |
Collapse
|
25
|
Post-Hartree-Fock methods for Hirshfeld atom refinement: are they necessary? Investigation of a strongly hydrogen-bonded molecular crystal. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2020.127934] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
26
|
A nitrogen-doped nanotube molecule with atom vacancy defects. Nat Commun 2020; 11:1807. [PMID: 32286324 PMCID: PMC7156684 DOI: 10.1038/s41467-020-15662-6] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2019] [Accepted: 03/18/2020] [Indexed: 12/02/2022] Open
Abstract
Nitrogen-doped carbon nanotubes have attracted attention in various fields, but lack of congeners with discrete molecular structures has hampered developments based on in-depth, chemical understandings. In this study, a nanotube molecule doped periodically with multiple nitrogen atoms has been synthesized by combining eight 2,4,6-trisubstituted pyridine units with thirty-two 1,3,5-trisubstituted benzene units. A synthetic strategy involving geodesic phenine frameworks is sufficiently versatile to tolerate pyridine units without requiring synthetic detours. Crystallographic analyses adopting aspherical multipole atom models reveal the presence of axially rotated structures as a minor disordered structure, which also provides detailed molecular and electronic structures. The nitrogen atoms on the nanotube serve as chemically distinct sites covered with negatively charged surfaces, and they increase the chance of electron injections by lowering the energy levels of the unoccupied orbitals that should serve as electron acceptors. Replacing carbon atoms in nanocarbons with heteroatoms alters their intrinsic properties, and nitrogen-doped nanocarbons attract much attention in various fields. Here, the authors synthesize a discrete nitrogen-doped nanotube molecule and clarify its structure to reveal unique features of nitrogen dopants.
Collapse
|
27
|
Dos Santos LH. Applications of charge-density analysis to the rational design of molecular materials: A mini review on how to engineer optical or magnetic crystals. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2019.127431] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
28
|
Gruza B, Chodkiewicz ML, Krzeszczakowska J, Dominiak PM. Refinement of organic crystal structures with multipolar electron scattering factors. Acta Crystallogr A Found Adv 2020; 76:92-109. [PMID: 31908353 PMCID: PMC8127334 DOI: 10.1107/s2053273319015304] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Accepted: 11/13/2019] [Indexed: 12/14/2022] Open
Abstract
A revolution in resolution is occurring now in electron microscopy arising from the development of methods for imaging single particles at cryogenic temperatures and obtaining electron diffraction data from nanocrystals of small organic molecules or macromolecules. Near-atomic or even atomic resolution of molecular structures can be achieved. The basis of these methods is the scattering of an electron beam due to the electrostatic potential of the sample. To analyse these high-quality experimental data, it is necessary to use appropriate atomic scattering factors. The independent atom model (IAM) is commonly used although various more advanced models, already known from X-ray diffraction, can also be applied to enhance the analysis. In this study a comparison is presented of IAM and TAAM (transferable aspherical atom model), the latter with the parameters of the Hansen-Coppens multipole model transferred from the University at Buffalo Databank (UBDB). By this method, TAAM takes into account the fact that atoms in molecules are partially charged and are not spherical. Structure refinements were performed on a carbamazepine crystal using electron structure-factor amplitudes determined experimentally [Jones et al. (2018). ACS Cent. Sci. 4, 1587-1592] or modelled with theoretical quantum-mechanical methods. The results show the possibilities and limitations of the TAAM method when applied to electron diffraction. Among others, the method clearly improves model fitting statistics, when compared with IAM, and allows for reliable refinement of atomic thermal parameters. The improvements are more pronounced with poorer-resolution diffraction data.
Collapse
Affiliation(s)
- Barbara Gruza
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Joanna Krzeszczakowska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| |
Collapse
|
29
|
Grabowsky S, Genoni A, Thomas SP, Jayatilaka D. The Advent of Quantum Crystallography: Form and Structure Factors from Quantum Mechanics for Advanced Structure Refinement and Wavefunction Fitting. 21ST CENTURY CHALLENGES IN CHEMICAL CRYSTALLOGRAPHY II 2020. [DOI: 10.1007/430_2020_62] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
30
|
Malaspina LA, Wieduwilt EK, Bergmann J, Kleemiss F, Meyer B, Ruiz-López MF, Pal R, Hupf E, Beckmann J, Piltz RO, Edwards AJ, Grabowsky S, Genoni A. Fast and Accurate Quantum Crystallography: From Small to Large, from Light to Heavy. J Phys Chem Lett 2019; 10:6973-6982. [PMID: 31633355 DOI: 10.1021/acs.jpclett.9b02646] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The coupling of the crystallographic refinement technique Hirshfeld atom refinement (HAR) with the recently constructed libraries of extremely localized molecular orbitals (ELMOs) gives rise to the new quantum-crystallographic method HAR-ELMO. This method is significantly faster than HAR but as accurate and precise, especially concerning the free refinement of hydrogen atoms from X-ray diffraction data, so that the first fully quantum-crystallographic refinement of a protein is presented here. However, the promise of HAR-ELMO exceeds large molecules and protein crystallography. In fact, it also renders possible electron-density investigations of heavy elements in small molecules and facilitates the detection and isolation of systematic errors from physical effects.
Collapse
Affiliation(s)
- Lorraine A Malaspina
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2 - Biologie/Chemie , Universität Bremen , Leobener Straße 3 und 7 , 28359 Bremen , Germany
| | - Erna K Wieduwilt
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2 - Biologie/Chemie , Universität Bremen , Leobener Straße 3 und 7 , 28359 Bremen , Germany
- Université de Lorraine , CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT) , 1 Boulevard Arago , 57078 Metz , France
| | - Justin Bergmann
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2 - Biologie/Chemie , Universität Bremen , Leobener Straße 3 und 7 , 28359 Bremen , Germany
| | - Florian Kleemiss
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2 - Biologie/Chemie , Universität Bremen , Leobener Straße 3 und 7 , 28359 Bremen , Germany
- Departement für Chemie und Biochemie , Universität Bern , Freiestrasse 3 , 3012 Bern , Switzerland
| | - Benjamin Meyer
- Université de Lorraine , CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT) , 1 Boulevard Arago , 57078 Metz , France
| | - Manuel F Ruiz-López
- Université de Lorraine , CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT) , 1 Boulevard Arago , 57078 Metz , France
| | - Rumpa Pal
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2 - Biologie/Chemie , Universität Bremen , Leobener Straße 3 und 7 , 28359 Bremen , Germany
| | - Emanuel Hupf
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2 - Biologie/Chemie , Universität Bremen , Leobener Straße 3 und 7 , 28359 Bremen , Germany
| | - Jens Beckmann
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2 - Biologie/Chemie , Universität Bremen , Leobener Straße 3 und 7 , 28359 Bremen , Germany
| | - Ross O Piltz
- Australian Nuclear Science and Technology Organisation , Australian Centre for Neutron Scattering , New Illawarra Road , Lucas Heights , NSW 2234 , Australia
| | - Alison J Edwards
- Australian Nuclear Science and Technology Organisation , Australian Centre for Neutron Scattering , New Illawarra Road , Lucas Heights , NSW 2234 , Australia
| | - Simon Grabowsky
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2 - Biologie/Chemie , Universität Bremen , Leobener Straße 3 und 7 , 28359 Bremen , Germany
- Departement für Chemie und Biochemie , Universität Bern , Freiestrasse 3 , 3012 Bern , Switzerland
| | - Alessandro Genoni
- Université de Lorraine , CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT) , 1 Boulevard Arago , 57078 Metz , France
| |
Collapse
|
31
|
Luger P, Dittrich B. Electron densities of two cyclononapeptides from invariom application. ZEITSCHRIFT FUR NATURFORSCHUNG SECTION B-A JOURNAL OF CHEMICAL SCIENCES 2019. [DOI: 10.1515/znb-2019-0154] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Abstract
The nonapeptides cyclo(Val-Leu-Pro-Ile-Leu-Leu-Leu-Val-Leu) (I) and cyclo(Val-Leu-Pro-Ala-Leu-Leu-Leu-Val-Leu) (II) were identified as promising candidates for the development as potential anti-cancer drugs. We report a re-refinement of deposited single-crystal X-ray diffraction data with aspherical scattering factors from the invariom database. A subsequent evaluation of the molecular electron density distribution and of the differences in their molecular electrostatic potentials provides insight in their activities. The sequences differ only in residue 4, Ile in (I) and Ala in (II). Since the anti-tumor potency is reduced for the Ala peptide (II), the causes for the differences seen in activity between (I) and (II) were examined from a structural and from an electron density (ED) point of view. The exchange at residue 4 does not lead to significant changes in molecular geometry. Molecular Hirshfeld surfaces and electrostatic potential (ESP) isosurfaces show accumulations of intermolecular interactions in regions adjacent to the Ile/Ala residues indicating preferred interactions with a potential receptor in these regions. The concentrations of intermolecular interactions were localized on the Hirshfeld surfaces through an extended basin of ED concentration close to the Ile/Ala residues. Differences in the electrostatic potentials (ESPs) between (I) and (II) were only found at the Ile/Ala site and were very close to zero otherwise.
Collapse
Affiliation(s)
- Peter Luger
- Institut für Chemie und Biochemie, Anorganische Chemie , Freie Universität Berlin , Fabeckstr. 36a, D-14195 Berlin , Germany
| | - Birger Dittrich
- Institut für Anorganische Chemie und Strukturchemie II , Heinrich-Heine-Universität Düsseldorf , Universitätsstr. 1, D-40225 Düsseldorf , Germany
| |
Collapse
|
32
|
Fabrizio A, Grisafi A, Meyer B, Ceriotti M, Corminboeuf C. Electron density learning of non-covalent systems. Chem Sci 2019; 10:9424-9432. [PMID: 32055318 PMCID: PMC6991182 DOI: 10.1039/c9sc02696g] [Citation(s) in RCA: 67] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2019] [Accepted: 09/08/2019] [Indexed: 11/21/2022] Open
Abstract
Chemists continuously harvest the power of non-covalent interactions to control phenomena in both the micro- and macroscopic worlds. From the quantum chemical perspective, the strategies essentially rely upon an in-depth understanding of the physical origin of these interactions, the quantification of their magnitude and their visualization in real-space. The total electron density ρ(r) represents the simplest yet most comprehensive piece of information available for fully characterizing bonding patterns and non-covalent interactions. The charge density of a molecule can be computed by solving the Schrödinger equation, but this approach becomes rapidly demanding if the electron density has to be evaluated for thousands of different molecules or very large chemical systems, such as peptides and proteins. Here we present a transferable and scalable machine-learning model capable of predicting the total electron density directly from the atomic coordinates. The regression model is used to access qualitative and quantitative insights beyond the underlying ρ(r) in a diverse ensemble of sidechain–sidechain dimers extracted from the BioFragment database (BFDb). The transferability of the model to more complex chemical systems is demonstrated by predicting and analyzing the electron density of a collection of 8 polypeptides. Machine learning model of the electron densities for analyzing non-covalent interaction patterns in peptides.![]()
Collapse
Affiliation(s)
- Alberto Fabrizio
- Laboratory for Computational Molecular Design , Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne , CH-1015 Lausanne , Switzerland . .,National Centre for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne , 1015 Lausanne , Switzerland
| | - Andrea Grisafi
- Laboratory of Computational Science and Modeling , IMX , École Polytechnique Fédérale de Lausanne , 1015 Lausanne , Switzerland.,National Centre for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne , 1015 Lausanne , Switzerland
| | - Benjamin Meyer
- Laboratory for Computational Molecular Design , Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne , CH-1015 Lausanne , Switzerland . .,National Centre for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne , 1015 Lausanne , Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling , IMX , École Polytechnique Fédérale de Lausanne , 1015 Lausanne , Switzerland.,National Centre for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne , 1015 Lausanne , Switzerland
| | - Clemence Corminboeuf
- Laboratory for Computational Molecular Design , Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne , CH-1015 Lausanne , Switzerland . .,National Centre for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne , 1015 Lausanne , Switzerland
| |
Collapse
|
33
|
Leduc T, Aubert E, Espinosa E, Jelsch C, Iordache C, Guillot B. Polarization of Electron Density Databases of Transferable Multipolar Atoms. J Phys Chem A 2019; 123:7156-7170. [PMID: 31294565 DOI: 10.1021/acs.jpca.9b05051] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Théo Leduc
- Université de Lorraine, CNRS, CRM2 , F-54000 Nancy , France
| | | | | | | | | | - Benoît Guillot
- Université de Lorraine, CNRS, CRM2 , F-54000 Nancy , France
| |
Collapse
|
34
|
Hasil A, Mehmood A, Ahmed M. Experimental and theoretical charge-density analysis of hippuric acid: insight into its binding with human serum albumin. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2019; 75:750-762. [PMID: 32830731 DOI: 10.1107/s2052520619007911] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2019] [Accepted: 06/02/2019] [Indexed: 11/10/2022]
Abstract
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.
Collapse
Affiliation(s)
- Asma Hasil
- Materials Chemistry Laboratory, Department of Chemistry, The Islamia University of Bahawalpur, Baghdad-ul-Jadeed Campus 63100, Pakistan
| | - Arshad Mehmood
- Department of Chemistry and Biochemistry, Texas Christian University, Fort Worth, Texas 76129, USA
| | - Maqsood Ahmed
- Materials Chemistry Laboratory, Department of Chemistry, The Islamia University of Bahawalpur, Baghdad-ul-Jadeed Campus 63100, Pakistan
| |
Collapse
|
35
|
Urzhumtsev AG, Lunin VY. Introduction to crystallographic refinement of macromolecular atomic models. CRYSTALLOGR REV 2019. [DOI: 10.1080/0889311x.2019.1631817] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexandre G. Urzhumtsev
- Centre for Integrative Biology, IGBMC, CNRS–INSERM–UdS, Illkirch, France
- Département de Physique, Faculté des Sciences et des Technologies, Université de Lorraine, Vandoeuvre-lès-Nancy, France
| | - Vladimir Y. Lunin
- Institute of Mathematical Problems of Biology RAS, Keldysh Institute of Applied Mathematics of Russian Academy of Sciences, Pushchino, Russia
| |
Collapse
|
36
|
Korlyukov AA, Nelyubina YV. Quantum chemical methods in charge density studies from X-ray diffraction data. RUSSIAN CHEMICAL REVIEWS 2019. [DOI: 10.1070/rcr4866] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
37
|
Kofoed PM, Hoser AA, Diness F, Capelli SC, Madsen AØ. X-ray diffraction data as a source of the vibrational free-energy contribution in polymorphic systems. IUCRJ 2019; 6:558-571. [PMID: 31316801 PMCID: PMC6608639 DOI: 10.1107/s2052252519003014] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/07/2018] [Accepted: 02/27/2019] [Indexed: 06/10/2023]
Abstract
In this contribution we attempt to answer a general question: can X-ray diffraction data combined with theoretical computations be a source of information about the thermodynamic properties of a given system? Newly collected sets of high-quality multi-temperature single-crystal X-ray diffraction data and complementary periodic DFT calculations of vibrational frequencies and normal mode vectors at the Γ point on the yellow and white polymorphs of di-methyl 3,6-di-chloro-2,5-di-hydroxy-terephthalate are combined using two different approaches, aiming to obtain thermodynamic properties for the two compounds. The first approach uses low-frequency normal modes extracted from multi-temperature X-ray diffraction data (normal coordinate analysis), while the other uses DFT-calculated low-frequency normal mode in the refinement of the same data (normal mode refinement). Thermodynamic data from the literature [Yang et al. (1989), Acta Cryst. B45, 312-323] and new periodic ab initio DFT supercell calculations are used as a reference point. Both approaches tested in this work capture the most essential features of the systems: the polymorphs are enantiotropically related, with the yellow form being the thermodynamically stable system at low temperature, and the white form at higher temperatures. However, the inferred phase transition temperature varies between different approaches. Thanks to the application of unconventional methods of X-ray data refinement and analysis, it was additionally found that, in the case of the yellow polymorph, anharmonicity is an important issue. By discussing contributions from low- and high-frequency modes to the vibrational entropy and enthalpy, the importance of high-frequency modes is highlighted. The analysis shows that larger anisotropic displacement parameters are not always related to the polymorph with the higher vibrational entropy contribution.
Collapse
Affiliation(s)
| | - Anna A. Hoser
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Frederik Diness
- Department of Chemistry, University of Copenhagen, Copenhagen Denmark
| | - Silvia C. Capelli
- ISIS Neutrons and Muons Facility, Science and Technical Facility Council, Rutherford Appleton Laboratory, Harwell Science Campus, Didcot OX11 OQX, UK
| | - Anders Østergaard Madsen
- Department of Pharmacy, University of Copenhagen, Universitetsparken 2, Copenhagen 2100, Denmark
| |
Collapse
|
38
|
Nguyen D, Volkov A. Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. II. The Fourier transform method. Acta Crystallogr A Found Adv 2019; 75:448-464. [PMID: 31041901 DOI: 10.1107/s2053273319002535] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2018] [Accepted: 02/18/2019] [Indexed: 11/10/2022] Open
Abstract
The Fourier transform method for analytical determination of the two-center Coulomb integrals needed for evaluation of the electrostatic interaction energies between pseudoatom-based charge distributions is presented, and its Fortran-based implementation using the 128-bit floating-point arithmetic in the XDPROP module of the XD software is described. In combination with mathematical libraries included in the Lahey/Fujitsu LF64 Linux compiler, the new implementation outperforms the previously reported Löwdin α-function technique [Nguyen et al. (2018). Acta Cryst. A74, 524-536] in terms of precision of the determined individual Coulomb integrals regardless of whether the latter uses the 64-, 80- or 128-bit precision floating-point format, all the while being only marginally slower. When the Löwdin α-function or Fourier transform method is combined with a multipole moment approximation for large interatomic separations (such a hybrid scheme is called the analytical exact potential and multipole moment method, aEP/MM) the resulting electrostatic interaction energies are evaluated with a precision of ≤5 × 10-5 kJ mol-1 for the current set of benchmark systems composed of H, C, N and O atoms and ranging in size from water-water to dodecapeptide-dodecapeptide dimers. Using a 2012 4.0 GHz AMD FX-8350 computer processor, the two recommended aEP/MM implementations, the 80-bit precision Löwdin α-function and 128-bit precision Fourier transform methods, evaluate the total electrostatic interaction energy between two 225-atom monomers of the benchmark dodecapeptide molecule in 6.0 and 7.9 s, respectively, versus 3.1 s for the previously reported 64-bit Löwdin α-function approach.
Collapse
Affiliation(s)
- Daniel Nguyen
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| | - Anatoliy Volkov
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| |
Collapse
|
39
|
Kumar P, Gruza B, Bojarowski SA, Dominiak PM. Extension of the transferable aspherical pseudoatom data bank for the comparison of molecular electrostatic potentials in structure-activity studies. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2019; 75:398-408. [PMID: 30821272 DOI: 10.1107/s2053273319000482] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Accepted: 01/09/2019] [Indexed: 12/31/2022]
Abstract
The transferable aspherical pseudoatom data bank, UBDB2018, is extended with over 130 new atom types present in small and biological molecules of great importance in biology and chemistry. UBDB2018 can be applied either as a source of aspherical atomic scattering factors in a standard X-ray experiment (dmin ≃ 0.8 Å) instead of the independent atom model (IAM), and can therefore enhance the final crystal structure geometry and refinement parameters; or as a tool to reconstruct the molecular charge-density distribution and derive the electrostatic properties of chemical systems for which 3D structural data are available. The extended data bank has been extensively tested, with the focus being on the accuracy of the molecular electrostatic potential computed for important drug-like molecules, namely the HIV-1 protease inhibitors. The UBDB allows the reconstruction of the reference B3LYP/6-31G** potentials, with a root-mean-squared error of 0.015 e bohr-1 computed for entire potential grids which span values from ca 200 e bohr-1 to ca -0.1 e bohr-1 and encompass both the inside and outside regions of a molecule. UBDB2018 is shown to be applicable to enhancing the physical meaning of the molecular electrostatic potential descriptors used to construct predictive quantitative structure-activity relationship/quantitative structure-property relationship (QSAR/QSPR) models for drug discovery studies. In addition, it is suggested that electron structure factors computed from UBDB2018 may significantly improve the interpretation of electrostatic potential maps measured experimentally by means of electron diffraction or single-particle cryo-EM methods.
Collapse
Affiliation(s)
- Prashant Kumar
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Barbara Gruza
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Sławomir Antoni Bojarowski
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa 02-089, Poland
| |
Collapse
|
40
|
Grisafi A, Fabrizio A, Meyer B, Wilkins DM, Corminboeuf C, Ceriotti M. Transferable Machine-Learning Model of the Electron Density. ACS CENTRAL SCIENCE 2019; 5:57-64. [PMID: 30693325 PMCID: PMC6346381 DOI: 10.1021/acscentsci.8b00551] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/10/2018] [Indexed: 05/24/2023]
Abstract
The electronic charge density plays a central role in determining the behavior of matter at the atomic scale, but its computational evaluation requires demanding electronic-structure calculations. We introduce an atom-centered, symmetry-adapted framework to machine-learn the valence charge density based on a small number of reference calculations. The model is highly transferable, meaning it can be trained on electronic-structure data of small molecules and used to predict the charge density of larger compounds with low, linear-scaling cost. Applications are shown for various hydrocarbon molecules of increasing complexity and flexibility, and demonstrate the accuracy of the model when predicting the density on octane and octatetraene after training exclusively on butane and butadiene. This transferable, data-driven model can be used to interpret experiments, accelerate electronic structure calculations, and compute electrostatic interactions in molecules and condensed-phase systems.
Collapse
Affiliation(s)
- Andrea Grisafi
- Laboratory
of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| | - Alberto Fabrizio
- Laboratory
for Computational Molecular Design, Institute of Chemical Sciences
and Engineering, École Polytechnique
Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| | - Benjamin Meyer
- Laboratory
for Computational Molecular Design, Institute of Chemical Sciences
and Engineering, École Polytechnique
Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| | - David M. Wilkins
- Laboratory
of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Clemence Corminboeuf
- Laboratory
for Computational Molecular Design, Institute of Chemical Sciences
and Engineering, École Polytechnique
Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| | - Michele Ceriotti
- Laboratory
of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
41
|
Faroque MU, Noureen S, Mirza SH, Tahir MN, Ahmed M. Structure and electrostatic properties of the pyrimethamine-3,5-dihydroxybenzoic acid cocrystal in water solvent studied using transferred electron-density parameters. Acta Crystallogr C Struct Chem 2019; 75:46-53. [PMID: 30601130 DOI: 10.1107/s2053229618017060] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2018] [Accepted: 11/30/2018] [Indexed: 11/10/2022] Open
Abstract
Pyrimethamine is an antimalarial drug. The cocrystal salt form of pyrimethamine with 3,5-dihydroxybenzoic acid in water solvent has been synthesized, namely 2,4-diamino-5-(4-chlorophenyl)-6-ethylpyrimidin-1-ium 3,5-dihydroxybenzoate hemihydrate, C12H14ClN4+·C7H5O4-·0.5H2O. X-ray diffraction data were collected at room temperature. Refinement of the crystal structure was carried out using the classical Independent Atom Model (IAM), while the electrostatic properties were studied by transferring electron-density parameters from an electron-density database. The Cl atom was refined anharmonically. The results of both refinement methods were compared. Topological analyses were carried out using Bader's theory of Atoms in Molecules (AIM). The three-dimensional Hirshfeld surface analysis and the two-dimensional fingerprint maps of individual molecules revealed that the crystal structures are dominated by H...O/O...H and H...H contacts. Other close contacts are also present, including weak C...H/H...C contacts. Charge transfer between the pyrimethamine and 3,5-dihydroxybenzoic acid molecules results in a molecular assembly based on strong intermolecular hydrogen bonds. This is further validated by the calculation of the electrostatic potential based on transferred electron-density parameters. The current work proves the significance of the transferability principle in studying the electron-density-derived properties of molecules in cases where high-resolution diffraction data at low temperature are not available.
Collapse
Affiliation(s)
- Muhammad Umer Faroque
- Materials Chemistry Laboratory, Department of Chemistry, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
| | - Sajida Noureen
- Materials Chemistry Laboratory, Department of Chemistry, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
| | | | | | - Maqsood Ahmed
- Materials Chemistry Laboratory, Department of Chemistry, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
| |
Collapse
|
42
|
Lübben J, Wandtke CM, Hübschle CB, Ruf M, Sheldrick GM, Dittrich B. Aspherical scattering factors for SHELXL - model, implementation and application. Acta Crystallogr A Found Adv 2019; 75:50-62. [PMID: 30575583 PMCID: PMC6302932 DOI: 10.1107/s2053273318013840] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2018] [Accepted: 09/29/2018] [Indexed: 11/23/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jens Lübben
- Institut für Anorganische Chemie der Universität Göttingen, Tammannstrasse 4, Göttingen, D-37077, Germany
- Bruker AXS Inc., 5465 E. Cheryl Parkway, Madison, WI 53711, USA
| | - Claudia M. Wandtke
- Institut für Anorganische Chemie der Universität Göttingen, Tammannstrasse 4, Göttingen, D-37077, Germany
| | | | - Michael Ruf
- Bruker AXS Inc., 5465 E. Cheryl Parkway, Madison, WI 53711, USA
| | - George M. Sheldrick
- Institut für Anorganische Chemie der Universität Göttingen, Tammannstrasse 4, Göttingen, D-37077, Germany
| | - Birger Dittrich
- Heinrich-Heine Universität Düsseldorf, Institut für Anorganische Chemie und Strukturchemie, Material- und Strukturforschung, Gebäude: 26.42, Universitätsstrasse 1, 40225 Düsseldorf, Germany
| |
Collapse
|
43
|
Bojarowski SA, Kumar P, Wandtke CM, Dittrich B, Dominiak PM. Universal Method for Electrostatic Interaction Energies Estimation with Charge Penetration and Easily Attainable Point Charges. J Chem Theory Comput 2018; 14:6336-6345. [PMID: 30359528 DOI: 10.1021/acs.jctc.8b00781] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sławomir A Bojarowski
- Biological and Chemical Research Centre, Department of Chemistry , University of Warsaw , ul. Żwirki i Wigury 101 , 02-089 Warszawa , Poland
| | - Prashant Kumar
- Biological and Chemical Research Centre, Department of Chemistry , University of Warsaw , ul. Żwirki i Wigury 101 , 02-089 Warszawa , Poland
| | - Claudia M Wandtke
- Institut fur Anorganische Chemie , Georg-August-Universitat , Tammannstrasse 4 , 37077 Göttingen , Germany
| | - Birger Dittrich
- Anorganische Chemie und Strukturchemie , Heinrich-Heine Universitat Düsseldorf , Universitatsstrasse 1, Gëbaude 26.42.01.21 , 40225 Düsseldorf , Germany
| | - Paulina M Dominiak
- Biological and Chemical Research Centre, Department of Chemistry , University of Warsaw , ul. Żwirki i Wigury 101 , 02-089 Warszawa , Poland
| |
Collapse
|
44
|
Meyer B, Genoni A. Libraries of Extremely Localized Molecular Orbitals. 3. Construction and Preliminary Assessment of the New Databanks. J Phys Chem A 2018; 122:8965-8981. [PMID: 30339393 DOI: 10.1021/acs.jpca.8b09056] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The fast and reliable determination of wave functions and electron densities of macromolecules has been one of the goals of theoretical chemistry for a long time, and in this context, several linear scaling techniques have been successfully devised over the years. Different approaches have been adopted to tackle this problem, and one of them exploits the fact that, according to the traditional chemical perception, molecules can be seen as constituted of recurring units (e.g., functional groups) with well-defined chemical features. This has led to the development of methods in which the global wave functions or electron densities of macromolecules are obtained by simply transferring density matrices or fuzzy electron densities associated with molecular fragments. In this context, we propose an alternative strategy that aims at quickly reconstructing wave functions and electron densities of proteins through the transfer of extremely localized molecular orbitals (ELMOs), which are orbitals strictly localized on small molecular units and, for this reason, easily transferable from molecule to molecule. To accomplish this task we have constructed original libraries of ELMOs that cover all the possible elementary fragments of the 20 natural amino acids in all their possible protonation states and forms. Our preliminary test calculations have shown that, compared to more traditional methods of quantum chemistry, the transfers from the novel ELMO databanks allow to obtain wave function and electron densities of large polypeptides and proteins at a significantly reduced computational cost. Furthermore, notwithstanding expected discrepancies, the obtained electron distributions and electrostatic potentials are in very good agreement with those obtained at Hartree-Fock and density functional theory (DFT) levels. Therefore, the results encourage to use the new libraries as alternatives to the popular pseudoatom-databases of crystallography in the refinement of crystallographic structures of macromolecules. In particular, in this context, we have already envisaged the coupling of the ELMO databanks with the promising Hirshfeld atom refinement technique to extend the applicability of the latter to very large systems.
Collapse
Affiliation(s)
- Benjamin Meyer
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019 , 1 Boulevard Arago , F-57078 Metz , France
| | - Alessandro Genoni
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019 , 1 Boulevard Arago , F-57078 Metz , France
| |
Collapse
|
45
|
Nguyen D, Kisiel Z, Volkov A. Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. I. The Löwdin α-function method. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2018; 74:524-536. [PMID: 30182939 DOI: 10.1107/s2053273318008690] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/06/2018] [Accepted: 06/13/2018] [Indexed: 11/11/2022]
Abstract
The previously reported [Volkov et al. (2004). Chem. Phys. Lett. 391, 170-175] exact potential and multipole moment (EP/MM) method for evaluation of intermolecular electrostatic interaction energies using the nuclei-centered pseudoatom representation of electron densities is significantly improved in terms of both speed and accuracy by replacing the numerical quadrature integration of the exact potential with a fully analytical technique. The resulting approach, incorporated in the XDPROP module of the software package XD, has been tested on several molecular systems ranging in size from water-water to dodecapeptide-dodecapeptide dimers using electron densities constructed via the University at Buffalo Aspherical Atom Databank. The improved hybrid method provides electrostatic interaction energies within the uncertainty of ≤0.2 kJ mol-1 for all benchmark systems. The running time for a dimer of a sizable, 225-atom dodecapeptide is under 4 s on a 2012 central processing unit (2.8 GHz AMD Opteron 6348) and under 3 s on a relatively modern processor (2.8 GHz Intel Xeon E3-1505M v5).
Collapse
Affiliation(s)
- Daniel Nguyen
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| | - Zbigniew Kisiel
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warszawa, Poland
| | - Anatoliy Volkov
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| |
Collapse
|
46
|
Sun CL, Liu LP, Tian F, Ding F, Wang LW. Charge-patching method for the calculation of electronic structure of polypeptides. Phys Chem Chem Phys 2018; 20:23301-23310. [DOI: 10.1039/c8cp01803k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Based on the CPM method, the charge densities of polypeptides can be generated and their electronic structure can be further calculated.
Collapse
Affiliation(s)
- Chang-Liang Sun
- Center of Physical Chemistry Test
- Shenyang University of Chemical Technology
- Shenyang 110142
- People's Republic of China
- Materials Science Division
| | - Li-Ping Liu
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- School of Physics
| | - Fubo Tian
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- College of Physics
| | - Fu Ding
- Center of Physical Chemistry Test
- Shenyang University of Chemical Technology
- Shenyang 110142
- People's Republic of China
| | - Lin-Wang Wang
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
| |
Collapse
|
47
|
Faroque MU, Noureen S, Ahmed M, Tahir MN. Electrostatic properties of the pyrimethamine–2,4-dihydroxybenzoic acid cocrystal in methanol studied using transferred electron-density parameters. ACTA CRYSTALLOGRAPHICA SECTION C-STRUCTURAL CHEMISTRY 2018; 74:100-107. [DOI: 10.1107/s2053229617017788] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2017] [Accepted: 12/12/2017] [Indexed: 11/11/2022]
Abstract
The crystal structure of the cocrystal salt form of the antimalarial drug pyrimethamine with 2,4-dihydroxybenzoic acid in methanol [systematic name: 2,4-diamino-5-(4-chlorophenyl)-6-ethylpyrimidin-1-ium 2,4-dihydroxybenzoate methanol monosolvate, C12H14ClN4
+·C7H5O4
−·CH3OH] has been studied using X-ray diffraction data collected at room temperature. The crystal structure was refined using the classical Independent Atom Model (IAM) and the Multipolar Atom Model by transferring electron-density parameters from the ELMAM2 database. The Cl atom was refined anharmonically. The results of both refinement methods have been compared. The intermolecular interactions have been characterized on the basis of Hirshfeld surface analysis and topological analysis using Bader's theory of Atoms in Molecules. The results show that the molecular assembly is built primarily on the basis of charge transfer between 2,4-dihydroxybenzoic acid and pyrimethamine, which results in strong intermolecular hydrogen bonds. This fact is further validated by the calculation of the electrostatic potential based on transferred electron-density parameters.
Collapse
|
48
|
Kovalenko AA, Nelyubina YV, Korlyukov AA, Lyssenko KA, Ananyev IV. The truth is out there: the metal-π interactions in crystal of Cr(CO)3(pcp) as revealed by the study of vibrational smearing of electron density. ACTA ACUST UNITED AC 2017. [DOI: 10.1515/zkri-2017-2085] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Abstract
The vibrational smearing of electron density was studied in the crystal of complex of Cr(CO)3 with [2.2]paracyclophane. The combination of theoretical and experimental methods, including periodic calculations and screening of DFT calculated and multipole-decomposed electron densities, was utilized to reveal the vibrational smearing of electron density and its influence on the multipole-constructed electron density. The multipole model, commonly used to treat the high-resolution X-ray diffraction data, was shown to be rather inaccurate in description of electron density and its vibrational smearing in metal-π complex where the interchange between diatomic interactions can occur. Namely, some bond critical points can be hidden while analyzing multipole-decomposed electron density with proved effects of vibrational smearing even if the deconvolution problem is overcome by using the invariom approach. On the contrary, the recently proposed “clouds of critical point variation” (CCPV) approach is demonstrated as the route to gather all reasonable bonding trends and to reconstruct static electron density pattern in metal-π complexes.
Collapse
Affiliation(s)
- Anna A. Kovalenko
- M.V. Lomonosov Moscow State University , GSP-1, Leninskie Gory Str. , 119991, Moscow , Russia
- A.N. Nesmeyanov Institute of Organoelement Compounds , Russian Academy of Sciences , Vavilova Str., 28 , 119991, Moscow , Russia
- Pirogov Russian National Research Medical University , Ostrovitianov Str., 1 , 117997, Moscow , Russia
- N.D. Zelinsky Institute of Organic Chemistry RAS , Leninsky Prospect, 47 , 119991, Moscow , Russian Federation
| | - Yulia V. Nelyubina
- M.V. Lomonosov Moscow State University , GSP-1, Leninskie Gory Str. , 119991, Moscow , Russia
- A.N. Nesmeyanov Institute of Organoelement Compounds , Russian Academy of Sciences , Vavilova Str., 28 , 119991, Moscow , Russia
- Pirogov Russian National Research Medical University , Ostrovitianov Str., 1 , 117997, Moscow , Russia
- N.D. Zelinsky Institute of Organic Chemistry RAS , Leninsky Prospect, 47 , 119991, Moscow , Russian Federation
| | - Alexander A. Korlyukov
- M.V. Lomonosov Moscow State University , GSP-1, Leninskie Gory Str. , 119991, Moscow , Russia
- A.N. Nesmeyanov Institute of Organoelement Compounds , Russian Academy of Sciences , Vavilova Str., 28 , 119991, Moscow , Russia
- Pirogov Russian National Research Medical University , Ostrovitianov Str., 1 , 117997, Moscow , Russia
- N.D. Zelinsky Institute of Organic Chemistry RAS , Leninsky Prospect, 47 , 119991, Moscow , Russian Federation
| | - Konstantin A. Lyssenko
- M.V. Lomonosov Moscow State University , GSP-1, Leninskie Gory Str. , 119991, Moscow , Russia
- A.N. Nesmeyanov Institute of Organoelement Compounds , Russian Academy of Sciences , Vavilova Str., 28 , 119991, Moscow , Russia
- Pirogov Russian National Research Medical University , Ostrovitianov Str., 1 , 117997, Moscow , Russia
- N.D. Zelinsky Institute of Organic Chemistry RAS , Leninsky Prospect, 47 , 119991, Moscow , Russian Federation
| | - Ivan V. Ananyev
- A.N. Nesmeyanov Institute of Organoelement Compounds , Russian Academy of Sciences , Vavilova Str., 28 , 119991, Moscow , Russia
- N.D. Zelinsky Institute of Organic Chemistry RAS , Leninsky Prospect, 47 , 119991, Moscow , Russian Federation
| |
Collapse
|
49
|
Wandtke CM, Weil M, Simpson J, Dittrich B. Using invariom modelling to distinguish correct and incorrect central atoms in `duplicate structures' with neighbouring 3d elements. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2017; 73:794-804. [PMID: 28980983 PMCID: PMC5628397 DOI: 10.1107/s2052520617010745] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/10/2017] [Accepted: 07/20/2017] [Indexed: 06/07/2023]
Abstract
Modelling coordination compounds has been shown to be feasible using the invariom method; for the best fit to a given set of diffraction data, additional steps other than using lookup tables of scattering factors need to be carried out. Here such procedures are applied to a number of `duplicate structures', where structures of two or more supposedly different coordination complexes with identical ligand environments, but with different 3d metal ions, were published. However, only one metal atom can be plausibly correct in these structures, and other spectroscopic data are unavailable. Using aspherical scattering factors, a structure can be identified as correct from the deposited Bragg intensities alone and modelling only the ligand environment often suffices to make this distinction. This is not possible in classical refinements using the independent atom model. Quantum-chemical computations of the better model obtained after aspherical-atom refinement further confirm the assignment of the element in the respective figures of merit.
Collapse
Affiliation(s)
- Claudia M. Wandtke
- Institut für Anorganische Chemie der Universität Göttingen, Tammannstrasse 4, Göttingen D-37077, Germany
| | - Matthias Weil
- Technische Universität Wien, Getreidemarkt 9/164-SC Stg 1, A-1060 Wien, Austria
| | - Jim Simpson
- University of Otago, PO Box 56, Dunedin, New Zealand
| | - Birger Dittrich
- Heinrich-Heine Universität Düsseldorf, Institut für Anorganische Chemie und Strukturchemie, Material- und Strukturforschung, Gebäude: 26.42, Universitätsstrasse 1, 40225 Düsseldorf, Germany
| |
Collapse
|
50
|
Bojarowski SA, Kumar P, Dominiak PM. Interplay of point multipole moments and charge penetration for intermolecular electrostatic interaction energies from the University at Buffalo pseudoatom databank model of electron density. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2017; 73:598-609. [PMID: 28762970 DOI: 10.1107/s2052520617005510] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/07/2017] [Accepted: 04/12/2017] [Indexed: 06/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sławomir A Bojarowski
- Biological and Chemical Research Center, Department of Chemistry, University of Warsaw, ul. Zwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Prashant Kumar
- Biological and Chemical Research Center, Department of Chemistry, University of Warsaw, ul. Zwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Paulina M Dominiak
- Biological and Chemical Research Center, Department of Chemistry, University of Warsaw, ul. Zwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|