Benjamin Meyer

Structure and Stability Studies of Pharmacologically Relevant S-Nitrosothiols: A Theoretical Approach

Nowadays, S-nitrosothiols (RSNOs) represent a promising class of nitric oxide (NO) donors that could be successfully used as drugs to compensate the decrease of NO production that usually arises in conjunction with cardiovascular diseases. Nevertheless, notwithstanding their pharmacological interest, the structure-stability relationship in RSNOs is still unclear, and this issue, together with the mechanism of NO donation in the physiological medium, deserves further investigation. As a first step forward in this direction, in this paper, the overall stability and structural preference of two pharmacologically relevant S-nitrosothiol molecules were studied in detail by means of computational strategies. In particular, performing calculations in implicit solvent (water) on the S-nitroso-N-acetylpenicillamine and the S-nitroso-N-acetylcysteine and analyzing the noncovalent interactions networks of their most stable conformers, we observed that the structure and the stability of these molecules can be directly related to the formation of stabilizing hydrogen-bond and chalcogen-chalcogen intramolecular interactions. The obtained results represent the starting point for further investigations to be conducted also on larger RSNOs to shed further light on the role played by intra- and intermolecular interactions and by solvation effects in stabilizing this class of molecules. The obtained insights will be hopefully helpful to design new RSNO-based drugs characterized by an enhanced pharmacological potency.

Can X-ray constrained Hartree–Fock wavefunctions retrieve electron correlation?

In this study, the X-ray constrained wavefunction approach is carefully investigated in order to assess its ability to capture the effect of electron correlation on electron density. Electron distributions obtained from highly correlated molecular wavefunctions are the benchmarks and their Fourier transforms are used to simulate X-ray intensities for the constrained wavefunction calculations.

Libraries of Extremely Localized Molecular Orbitals. 2. Comparison with the Pseudoatoms Transferability

Due to both technical and methodological difficulties, determining and analyzing charge densities of very large molecular systems represents a serious challenge that, in the crystallographers community, has been mainly tackled by observing that the so-called pseudoatoms of the electron density multipole expansions are reliably transferable from molecule to molecule. This has led to the construction of pseudoatoms databanks that have allowed successful refinements of crystallographic structures of macromolecules, while taking into account their corresponding reconstructed electron distributions. A recent alternative/complement to the previous approach is represented by techniques based on extremely localized molecular orbitals (ELMOs) that, due to their strict localization on small molecular fragments (e.g., atoms, bonds, and functional groups), are also in principle exportable from system to system. The ELMOs transferability has been already tested in detail, and, in this work, it has been compared to the one of the pseudoatoms. To accomplish this task, electron distributions obtained both through the transfer of pseudoatoms and through the transfer of extremely localized molecular orbitals have been analyzed, especially taking into account topological properties and similarity indexes. The obtained results indicate that all the considered reconstruction methods give completely reasonable and similar charge densities, and, consequently, the new ELMOs libraries will probably represent new useful tools not only for refining crystal structures but also for computing approximate electronic properties of very large molecules.

Libraries of Extremely Localized Molecular Orbitals. 1. Model Molecules Approximation and Molecular Orbitals Transferability

Despite more and more remarkable computational ab initio results are nowadays continuously obtained for large macromolecular systems, the development of new linear-scaling techniques is still an open and stimulating field of research in theoretical chemistry. In this family of methods, an important role is occupied by those strategies based on the observation that molecules are generally constituted by recurrent functional units with well-defined intrinsic features. In this context, we propose to exploit the notion of extremely localized molecular orbitals (ELMOs) that, due to their strict localization on small molecular fragments (e.g., atoms, bonds, or functional groups), are in principle transferable from one molecule to another. Accordingly, the construction of orbital libraries to almost instantaneously build up approximate wave functions and electron densities of very large systems becomes conceivable. In this work, the ELMOs transferability is further investigated in detail and, furthermore, suitable rules to construct model molecules for the computation of ELMOs to be stored in future databanks are also defined. The obtained results confirm the reliable transferability of the ELMOs and show that electron densities obtained from the transfer of extremely localized molecular orbitals are very close to the corresponding Hartree-Fock ones. These observations prompt us to construct new ELMOs databases that could represent an alternative/complement to the already popular pseudoatoms databanks both for determining electron densities and for refining crystallographic structures of very large molecules.

Exploring charge density analysis in crystals at high pressure: data collection, data analysis and advanced modelling

The possibility to determine electron-density distribution in crystals has been an enormous breakthrough, stimulated by a favourable combination of equipment for X-ray and neutron diffraction at low temperature, by the development of simplified, though accurate, electron-density models refined from the experimental data and by the progress in charge density analysis often in combination with theoretical work. Many years after the first successful charge density determination and analysis, scientists face new challenges, for example: (i) determination of the finer details of the electron-density distribution in the atomic cores, (ii) simultaneous refinement of electron charge and spin density or (iii) measuring crystals under perturbation. In this context, the possibility of obtaining experimental charge density at high pressure has recently been demonstrated [Casati et al. (2016). Nat. Commun. 7, 10901]. This paper reports on the necessities and pitfalls of this new challenge, focusing on the species syn-1,6:8,13-biscarbonyl[14]annulene. The experimental requirements, the expected data quality and data corrections are discussed in detail, including warnings about possible shortcomings. At the same time, new modelling techniques are proposed, which could enable specific information to be extracted, from the limited and less accurate observations, like the degree of localization of double bonds, which is fundamental to the scientific case under examination.

Libraries of Extremely Localized Molecular Orbitals. 3. Construction and Preliminary Assessment of the New Databanks

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.

Analyzing Fluxional Molecules Using DORI

The Density Overlap Region Indicator (DORI) is a density-based scalar field that reveals covalent bonding patterns and noncovalent interactions in the same value range. This work goes beyond the traditional static quantum chemistry use of scalar fields and illustrates the suitability of DORI for analyzing geometrical and electronic signatures in highly fluxional molecular systems. Examples include a dithiocyclophane, which possesses multiple local minima with differing extents of π-stacking interactions and a temperature dependent rotation of a molecular rotor, where the descriptor is employed to capture fingerprints of CH-π and π-π interactions. Finally, DORI serves to examine the fluctuating π-conjugation pathway of a photochromic torsional switch (PTS). Attention is also placed on postprocessing the large amount of generated data and juxtaposing DORI with a data-driven low-dimensional representation of the structural landscape.