Invariom refinement of a new monoclinic solvate of thiostrepton at 0.64 Å resolution

The Relevance of Experimental Charge Density Analysis in Unraveling Noncovalent Interactions in Molecular Crystals

The work carried out by our research group over the last couple of decades in the context of quantitative crystal engineering involves the analysis of intermolecular interactions such as carbon (tetrel) bonding, pnicogen bonding, chalcogen bonding, and halogen bonding using experimental charge density methodology is reviewed. The focus is to extract electron density distribution in the intermolecular space and to obtain guidelines to evaluate the strength and directionality of such interactions towards the design of molecular crystals with desired properties. Following the early studies on halogen bonding interactions, several "sigma-hole" interaction types with similar electrostatic origins have been explored in recent times for their strength, origin, and structural consequences. These include interactions such as carbon (tetrel) bonding, pnicogen bonding, chalcogen bonding, and halogen bonding. Experimental X-ray charge density analysis has proved to be a powerful tool in unraveling the strength and electronic origin of such interactions, providing insights beyond the theoretical estimates from gas-phase molecular dimer calculations. In this mini-review, we outline some selected contributions from the X-ray charge density studies to the field of non-covalent interactions (NCIs) involving elements of the groups 14-17 of the periodic table. Quantitative insights into the nature of these interactions obtained from the experimental electron density distribution and subsequent topological analysis by the quantum theory of atoms in molecules (QTAIM) have been discussed. A few notable examples of weak interactions have been presented in terms of their experimental charge density features. These examples reveal not only the strength and beauty of X-ray charge density multipole modeling as an advanced structural chemistry tool but also its utility in providing experimental benchmarks for the theoretical studies of weak interactions in crystals.

The bacterial thiopeptide thiostrepton. An update of its mode of action, pharmacological properties and applications

The bacterial thiopeptide thiostrepton (TS) is used as a veterinary medicine to treat bacterial infections. TS is a protein translation inhibitor, essentially active against Gram-positive bacteria and some Gram-negative bacteria. In procaryotes, TS abrogates binding of GTPase elongation factors to the 70S ribosome, by altering the structure of rRNA-L11 protein complexes. TS exerts also antimalarial effects by disrupting protein synthesis in the apicoplast genome of Plasmodium falciparum. Interestingly, the drug targets both the infectious pathogen (bacteria or parasite) and host cell, by inducing endoplasmic reticulum stress-mediated autophagy which contributes to enhance the host cell defense. In addition, TS has been characterized as a potent chemical inhibitor of the oncogenic transcription factor FoxM1, frequently overexpressed in cancers or other diseases. The capacity of TS to crosslink FoxM1, and a few other proteins such as peroxiredoxin 3 (PRX3) and the 19S proteasome, contributes to the anticancer effects of the thiopeptide. The anticancer activities of TS evidenced using diverse tumor cell lines, in vivo models and drug combinations are reviewed here, together with the implicated targets and mechanisms. The difficulty to formulate TS is a drag on the pharmaceutical development of the natural product. However, the design of hemisynthetic analogues and the use of micellar drug delivery systems should facilitate a broader utilization of the compound in human and veterinary medicines. This review shed light on the many pharmacological properties of TS, with the objective to promote its use as a pharmacological tool and medicinal product.

Fast energy minimization of the CCDC drug-subset structures by molecule-in-cluster computations allows independent structure validation and model completion

Optimizing structures with computations on clusters of molecules permits generation of structure-specific restraints for refinement and structure validation.

Fast and Accurate Quantum Crystallography: From Small to Large, from Light to Heavy

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.

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.

Amino acid modifications for conformationally constraining naturally occurring and engineered peptide backbones: Insights from the Protein Data Bank

Extensive efforts invested in understanding the rules of protein folding are now being applied, with good effect, in de novo design of proteins/peptides. For proteins containing standard α-amino acids alone, knowledge derived from experimentally determined three-dimensional (3D) structures of proteins and biologically active peptides are available from the Protein Data Bank (PDB), and the Cambridge Structural Database (CSD). These help predict and design protein structures, with reasonable confidence. However, our knowledge of 3D structures of biomolecules containing backbone modified amino acids is still evolving. A major challenge in de novo protein/peptide design concerns the engineering of conformationally constrained molecules with specific structural elements and chemical groups appropriately positioned for biological activity. This review explores four classes of amino acid modifications that constrain protein/peptide backbone structure. Systematic analysis of peptidic molecule structures (eg, bioactive peptides, inhibitors, antibiotics, and designed molecules), containing these backbone-modified amino acids, found in the PDB and CSD are discussed. The review aims to provide structure-function insights that will guide future design of proteins/peptides.

DiSCaMB: a software library for aspherical atom model X-ray scattering factor calculations with CPUs and GPUs

It has been recently established that the accuracy of structural parameters from X-ray refinement of crystal structures can be improved by using a bank of aspherical pseudoatoms instead of the classical spherical model of atomic form factors. This comes, however, at the cost of increased complexity of the underlying calculations. In order to facilitate the adoption of this more advanced electron density model by the broader community of crystallographers, a new software implementation called DiSCaMB, 'densities in structural chemistry and molecular biology', has been developed. It addresses the challenge of providing for high performance on modern computing architectures. With parallelization options for both multi-core processors and graphics processing units (using CUDA), the library features calculation of X-ray scattering factors and their derivatives with respect to structural parameters, gives access to intermediate steps of the scattering factor calculations (thus allowing for experimentation with modifications of the underlying electron density model), and provides tools for basic structural crystallographic operations. Permissively (MIT) licensed, DiSCaMB is an open-source C++ library that can be embedded in both academic and commercial tools for X-ray structure refinement.

Invariom modeling of disordered structures: case studies on a dipeptide, an amino acid, and cefaclor, a cephalosporin antibiotic

Routines to facilitate the treatment of disorder in invariom modeling have been implemented in the open-source program MolecoolQt, a visualization program for charge-density work, and InvariomTool, a pre-processor program. Two published structures of an amino acid and a dipeptide and the new structure of cefaclor, a cephalosporin antibiotic, provide examples with increasing amounts of disorder, which can now be successfully modeled with invarioms. Like for ordered structures, these non-spherical scattering factors predicted by density functional theory significantly improve the structural model (figures of merit and standard deviations) also in these cases. Furthermore, they allow rapid calculation and comparison of the electrostatic potential and the molecular dipole moment for the different conformers present in the crystal structures.

Invariom based electron density studies on the C/Si analogues haloperidol/sila-haloperidol and venlafaxine/sila-venlafaxine

Hirshfeld surfaces of haloperidol hydropicrate (left) and sila-haloperidol hydrochloride (right) show comparable sites of ED concentrations due to comparable intermolecular environments.

Contributions of charge-density research to medicinal chemistry

This article reviews efforts in accurate experimental charge-density studies with relevance to medicinal chemistry. Initially, classical charge-density studies that measure electron density distribution via least-squares refinement of aspherical-atom population parameters are summarized. Next, interaction density is discussed as an idealized situation resembling drug-receptor interactions. Scattering-factor databases play an increasing role in charge-density research, and they can be applied both to small-molecule and macromolecular structures in refinement and analysis; software development facilitates their use. Therefore combining both of these complementary branches of X-ray crystallography is recommended, and examples are given where such a combination already proved useful. On the side of the experiment, new pixel detectors are allowing rapid measurements, thereby enabling both high-throughput small-molecule studies and macromolecular structure determination to higher resolutions. Currently, the most ambitious studies compute intermolecular interaction energies of drug-receptor complexes, and it is recommended that future studies benefit from recent method developments. Selected new developments in theoretical charge-density studies are discussed with emphasis on its symbiotic relation to crystallography.

On the temperature dependence of H-U(iso) in the riding hydrogen model

The temperature dependence of hydrogen U_iso and parent U_eq in the riding hydrogen model is investigated by neutron diffraction, aspherical-atom refinements and QM/MM and MO/MO cluster calculations. Fixed values of 1.2 or 1.5 appear to be underestimated, especially at temperatures below 100 K.

The temperature dependence of H-U_iso in N-acetyl-l-4-hydroxyproline monohydrate is investigated. Imposing a constant temperature-independent multiplier of 1.2 or 1.5 for the riding hydrogen model is found to be inaccurate, and severely underestimates H-U_iso below 100 K. Neutron diffraction data at temperatures of 9, 150, 200 and 250 K provide benchmark results for this study. X-ray diffraction data to high resolution, collected at temperatures of 9, 30, 50, 75, 100, 150, 200 and 250 K (synchrotron and home source), reproduce neutron results only when evaluated by aspherical-atom refinement models, since these take into account bonding and lone-pair electron density; both invariom and Hirshfeld-atom refinement models enable a more precise determination of the magnitude of H-atom displacements than independent-atom model refinements. Experimental efforts are complemented by computing displacement parameters following the TLS+ONIOM approach. A satisfactory agreement between all approaches is found.