The advanced treatment of hydrogen bonding in quantum crystallography

Optimization of Quantum Nuclei Positions with the Adaptive Nuclear-Electronic Orbital Approach

The use of multicomponent methods has become increasingly popular over the last years. Under this framework, nuclei (commonly protons) are treated quantum mechanically on the same footing as the electronic structure problem. Under the use of atomic-centered orbitals, this can lead to some complications as the ideal location of the nuclear basis centers must be optimized. In this contribution, we propose a straightforward approach to determine the position of such centers within the self-consistent cycle of a multicomponent calculation, making use of individual proton charge centroids. We test the method on model systems including the water dimer, a protonated water tetramer, and a porphine system. Comparing to numerical gradient calculations, the adaptive nuclear-electronic orbital (NEO) procedure is able to converge the basis centers to within a few cents of an Ångström and with less than 0.1 kcal/mol differences in absolute energies. This is achieved in one single calculation and with a small added computational effort of up to 80% compared to a regular NEO- self-consistent field run. An example application for the human transketolase proton wire is also provided.

Comparative Conformational Analysis of Acyclic Sugar Alcohols Ribitol, Xylitol and d-Arabitol by Solution NMR and Molecular Dynamics Simulations

Ribitol (C5H12O5) is an acyclic sugar alcohol that was recently identified in O-mannose glycan on mammalian α-dystroglycan. The conformation and dynamics of acyclic sugar alcohols such as ribitol are dependent on the stereochemistry of the hydroxyl groups; however, the dynamics are not fully understood. To gain insights into the conformation and dynamics of sugar alcohols, we carried out comparative analyses of ribitol, d-arabitol and xylitol by a crystal structure database search, solution NMR analysis and molecular dynamics (MD) simulations. The crystal structures of the sugar alcohols showed a limited number of conformations, suggesting that only certain stable conformations are prevalent among all possible conformations. The three-bond scholar coupling constants and exchange rates of hydroxyl protons were measured to obtain information on the backbone torsion angle and possible hydrogen bonding of each hydroxyl group. The 100 ns MD simulations indicate that the ribitol backbone has frequent conformational transitions with torsion angles between 180∘ and ±60∘, while d-arabitol and xylitol showed fewer conformational transitions. Taking our experimental and computational data together, it can be concluded that ribitol is more flexible than d-arabitol or xylitol, and the flexibility is at least in part defined by the configuration of the OH groups, which may form intramolecular hydrogen bonds.

Influence of N-protonation on electronic properties of acridine derivatives by quantum crystallography

Applications of 9-aminoacridine (9aa) and its derivatives span fields such as chemistry, biology, and medicine, including anticancer and antimicrobial activities. Protonation of such molecules can alter their bioavailability as weakly basic drugs like aminoacridines exhibit reduced solubility at high pH levels potentially limiting their effectiveness in patients with elevated gastric pH. In this study, we analyse the influence of protonation on the electronic characteristics of the molecular organic crystals of 9-aminoacridine. The application of quantum crystallography, including aspherical atom refinement, has enriched the depiction of electron density in the studied systems and non-covalent interactions, providing more details than previous studies. Our experimental results, combined with a topological analysis of the electron density and its Laplacian, provided detailed descriptions of how protonation changes the electron density distribution around the amine group and water molecule, concurrently decreasing the electron density at bond critical points of N/O-H bonds. Protonation also alters the molecular architecture of the systems under investigation. This is reflected in different proportions of the N⋯H and O⋯H intermolecular contacts for the neutral and protonated forms. Periodic DFT calculations of the cohesive energies of the crystal lattice, as well as computed interaction energies between molecules in the crystal, confirm that protonation stabilises the crystal structure due to a positive synergy between strong halogen and hydrogen bonds. Our findings highlight the potential of quantum crystallography in predicting crystal structure properties and point to its possible applications in developing new formulations for poorly soluble drugs.

How do density functionals affect the Hirshfeld atom refinement?

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

X-ray constrained wavefunctions based on Hirshfeld atoms. I. Method and review

The X-ray constrained wavefunction (XCW) procedure for obtaining an experimentally reconstructed wavefunction from X-ray diffraction data is reviewed. The two-center probability distribution model used to perform nuclear-position averaging in the original paper [Grimwood & Jayatilaka (2001). Acta Cryst. A57, 87-100] is carefully distinguished from the newer one-center probability distribution model. In the one-center model, Hirshfeld atoms are used, and the Hirshfeld atom based X-ray constrained wavefunction (HA-XCW) procedure is described for the first time, as well as its efficient implementation. In this context, the definition of the related X-ray wavefunction refinement (XWR) method is refined. The key halting problem for the XCW method - the procedure by which one determines when overfitting has occurred - is named and work on it reviewed.

Hirshfeld atom refinement based on projector augmented wave densities with periodic boundary conditions

Hirshfeld atom refinement (HAR) is an X-ray diffraction refinement method that, in numerous publications, has been shown to give H-atom bond lengths in close agreement with neutron diffraction derived values. Presented here is a first evaluation of an approach using densities derived from projector augmented wave (PAW) densities with three-dimensional periodic boundary conditions for HAR. The results show an improvement over refinements that neglect the crystal environment or treat it classically, while being on a par with non-periodic approximations for treating the solid-state environment quantum mechanically. A suite of functionals were evaluated for this purpose, showing that the SCAN and revSCAN functionals are most suited to these types of calculation.

Modeling of the Response of Hydrogen Bond Properties on an External Electric Field: Geometry, NMR Chemical Shift, Spin-Spin Scalar Coupling

The response of the geometric and NMR properties of molecular systems to an external electric field has been studied theoretically in a wide field range. It has been shown that this adduct under field approach can be used to model the geometric and spectral changes experienced by molecular systems in polar media if the system in question has one and only one bond, the polarizability of which significantly exceeds the polarizability of other bonds. If this requirement is met, then it becomes possible to model even extreme cases, for example, proton dissociation in hydrogen halides. This requirement is fulfilled for many complexes with one hydrogen bond. For such complexes, this approach can be used to facilitate a detailed analysis of spectral changes associated with geometric changes in the hydrogen bond. For example, in hydrogen-bonded complexes of isocyanide C≡¹⁵N-¹H⋯X, ¹J(¹⁵N¹H) depends exclusively on the N-H distance, while δ(¹⁵N) is also slightly influenced by the nature of X.