The advanced treatment of hydrogen bonding in quantum crystallography

Although hydrogen bonding is one of the most important motifs in chemistry and biology, H-atom parameters are especially problematic to refine against X-ray diffraction data. New developments in quantum crystallography offer a remedy. This article reports how hydrogen bonds are treated in three different quantum-crystallographic methods: Hirshfeld atom refinement (HAR), HAR coupled to extremely localized molecular orbitals and X-ray wavefunction refinement. Three different compound classes that form strong intra- or intermolecular hydrogen bonds are used as test cases: hydrogen maleates, the tripeptide l-alanyl-glycyl-l-alanine co-crystallized with water, and xylitol. The differences in the quantum-mechanical electron densities underlying all the used methods are analysed, as well as how these differences impact on the refinement results.

Similarities and Differences between Crystal and Enzyme Environmental Effects on the Electron Density of Drug Molecules

The crystal interaction density is generally assumed to be a suitable measure of the polarization of a low-molecular weight ligand inside an enzyme, but this approximation has seldomly been tested and has never been quantified before. In this study, we compare the crystal interaction density and the interaction electrostatic potential for a model compound of loxistatin acid (E64c) with those inside cathepsin B, in solution, and in vacuum. We apply QM/MM calculations and experimental quantum crystallography to show that the crystal interaction density is indeed very similar to the enzyme interaction density. Less than 0.1 e are shifted between these two environments in total. However, this difference has non-negligible consequences for derived properties.

Dynamical properties of enzyme-substrate complexes disclose substrate specificity of the SARS-CoV-2 main protease as characterized by the electron density descriptors

A dynamical approach is proposed to discriminate between reactive (rES) and nonreactive (nES) enzyme-substrate complexes taking the SARS-CoV-2 main protease (Mpro) as an important example. Molecular dynamics simulations with the quantum mechanics/molecular mechanics potentials (QM(DFT)/MM-MD) followed by the electron density analysis are employed to evaluate geometry and electronic properties of the enzyme with different substrates along MD trajectories. We demonstrate that mapping the Laplacian of the electron density and the electron localization function provides easily visible images of the substrate activation that allow one to distinguish rES and nES. The computed fractions of reactive enzyme-substrate complexes along MD trajectories well correlate with the findings of recent experimental studies on the substrate specificity of Mpro. The results of our simulations demonstrate the role of the theory level used in QM subsystems for a proper description of the nucleophilic attack of the catalytic cysteine residue in Mpro. The activation of the carbonyl group of a substrate is correctly characterized with the hybrid DFT functional PBE0, whereas the use of a GGA-type PBE functional, that lacks the admixture of the Hartree-Fock exchange fails to describe activation.

Invariom-based comparative electron density studies of iso-sildenafil and sildenafil

The title compounds have raised considerable medical and broad public interest in that sildenafil is used as an agent against male erectile dysfunction; iso-sildenafil is not in clinical use. A comparison of their structural and electronic properties therefore seems of interest. The electron densities of iso-sildenafil and the cationic and neutral forms of sildenafil were examined by the application of the invariom formalism relying on diffraction data reported in the literature. The electron-density distributions obtained were subjected to topological analysis using the quantum theory of atoms in molecules (QTAIM) formalism to yield bond topological and atomic properties. Moreover, molecular Hirshfeld surfaces and electrostatic potentials (ESPs) were calculated. A number of structural and electronic differences were thus identified between sildenafil and the iso-analog. In both sildenafil structures, the phenyl ring and the pyrazolopyrimidine fragment are practically coplanar (planar conformation), whereas in the iso-analog they exhibit an angle of 44° (inclined form). Related to differences in molecular structure are completely different hydrogen bonding patterns and differences in the ESPs, the latter ones being influenced by different methylation at the pyrazolopyrimidine fragment. Iso-sildenafil is present as a hydrogen-bond dimer in the crystal, and the ESP of this dimer is dominated by a surrounding positive potential.