Christian B. Hübschle

ShelXle: a Qt graphical user interface for SHELXL

ShelXle is a user-friendly graphical user interface for SHELXL. It combines an editor with syntax highlighting for SHELXL-associated files with an interactive graphical display for visualization of a three-dimensional structure.

MolIso – a program for colour-mapped iso-surfaces

The program MolIso allows the visualization of colour-mapped iso-surfaces from XD grid or GAUSSIAN CUBE files. It is written in C++ and uses the OpenGL and GLUT libraries, and is therefore portable. Linux and Windows versions are currently available. Optically correct looking, transparent visualizations are possible. The molecule inside the iso-surface can be depicted as ball-and-stick or atomic displacement parameter ellipsoid–line models. Screenshots can easily be taken in Portable Pix Map format. Short video clips, with various codecs for presentations, are easy to create.

The invariom model and its application: refinement of d,l-serine at different temperatures and resolution

Three X-ray data sets of the same D,L-serine crystal were measured at temperatures of 298, 100 and 20 K. These data were then evaluated using invarioms and the Hansen & Coppens aspherical-atom model. Multipole populations for invarioms, which are pseudoatoms that remain approximately invariant in an intermolecular transfer, were theoretically predicted using different density functional theorem (DFT) basis sets. The invariom parameters were kept fixed and positional and thermal parameters were refined to compare the fitting against the multi-temperature data at different resolutions. The deconvolution of thermal motion and electron density with respect to data resolution was studied by application of the Hirshfeld test. Above a resolution of sin theta/lambda approximately 0.55 A-1, or d approximately 0.9 A, this test was fulfilled. When the Hirshfeld test is fulfilled, a successful modeling of the aspherical electron density with invarioms is achieved, which was proven by Fourier methods. Molecular geometry improves, especially for H atoms, when using the invariom method compared to the independent-atom model, as a comparison with neutron data shows. Based on this example, the general applicability of the invariom concept to organic molecules is proven and the aspherical density modeling of a larger biomacromolecule is within reach.

MoleCoolQt – a molecule viewer for charge-density research

MoleCoolQt is a molecule viewer designed for experimental charge density studies with a user-friendly graphical user interface.

Automation of invariom and of experimental charge density modelling of organic molecules with the preprocessor program InvariomTool

The program InvariomTool can be used to obtain high-quality X-ray structures of organic molecules by automating both invariom modelling and the modelling process employed in experimental charge density studies. InvariomTool is a preprocessor program for the XD package [Koritsanszky et al. (2003), Technical Report, Freie Universitat Berlin] and allows the analysis of a structure in terms of the local atomic bonding environment in order to assign Hansen–Coppens pseudoatoms that are transferable from one molecule to another (invarioms). It relies on a database of invariom entries, each containing the invariom name, local atomic site symmetry, coordinate system, model-compound name and theoretically predicted multipole population parameters. The information on chemical equivalence and local atomic site symmetry determines which multipole parameters are to be refined in the least-squares procedure of an experimental charge density study. InvariomTool allows the user to generate input files either for invariom refinement, where parameters are fixed and taken from the database, or for an experimental refinement of multipole parameters.

The generalized invariom database (GID)

Invarioms are aspherical atomic scattering factors that enable structure refinement of more accurate and more precise geometries than refinements with the conventional independent atom model (IAM). The use of single-crystal X-ray diffraction data of a resolution better than sin θ/λ = 0.6 Å(-1) (or d = 0.83 Å) is recommended. The invariom scattering-factor database contains transferable pseudoatom parameters of the Hansen-Coppens multipole model and associated local atomic coordinate systems. Parameters were derived from geometry optimizations of suitable model compounds, whose IUPAC names are also contained in the database. Correct scattering-factor assignment and orientation reproduces molecular electron density to a good approximation. Molecular properties can hence be derived directly from the electron-density model. Coverage of chemical environments in the invariom database has been extended from the original amino acids, proteins and nucleic acid structures to many other environments encountered in organic chemistry. With over 2750 entries it now covers a wide sample of general organic chemistry involving the elements H, C, N and O, and to a lesser extent F, Si, S, P and Cl. With respect to the earlier version of the database, the main modification concerns scattering-factor notation. Modifications improve ease of use and success rates of automatic geometry-based scattering-factor assignment, especially in condensed hetero-aromatic ring systems, making the approach well suited to replace the IAM for structures of organic molecules.

Anharmonic Motion in Experimental Charge Density Investigations

In the charge density study of 9-diphenylthiophosphinoylanthracene the thermal motion of several atoms needed an anharmonic description via Gram-Charlier coefficients even for data collected at 15 K. As several data sets at different temperatures were measured, this anharmonic model could be proved to be superior to a disorder model. Refinements against theoretical data showed the resemblance of an anharmonic model and a disorder model with two positions very close to each other (~0.2 Å), whereas these two models could be clearly distinguished if the second position is 0.5 Å apart. The refined multipole parameters were distorted when the anharmonic motion was not properly refined. Therefore, this study reveals the importance of detecting and properly handling anharmonic motion. Unrefined anharmonic motion leads to typical shashlik-like residual density patterns. Therefore, careful analysis of the residual density and the derived probability density function after the refinement of the Gram-Charlier coefficients proved to be the most useful tools to indicate the presence of anharmonic motion.

Verification of structural and electrostatic properties obtained by the use of different pseudoatom databases.

The existing pseudoatom databases (ELMAM, Invariom, UBDB and ELMAM2) enable structure refinement to be performed with the use of aspherical scattering factors computed from the transferable aspherical atom model (TAAM) as an alternative to independent atom model refinement. In addition, electrostatic properties can be estimated with the help of the databases. The quality of the structural and electrostatic properties obtained from the individual databases was tested. On the basis of a 100 K high-resolution single-crystal X-ray diffraction experiment on L-His-L-Ala dihydrate and 23 K high-resolution data for L-Ala [Destro & Marsh (1988). J. Phys. Chem. 92, 966-973], the structural properties, electron-density distributions and molecular electrostatic potentials obtained from different TAAMs were compared to each other and to reference models. Experimental multipolar models and theoretical models refined against theoretical structure factors computed from periodic density functional theory (DFT) calculations were compared to the TAAMs in order to determine which model best describes the crystal-field effect. Unperturbed wavefunctions based on the MP2 and DFT calculations and properties obtained directly from these were used as a reference to judge how properly the databases reproduce the properties of isolated molecules. For Gly-L-His dihydrate, D,L-His and the above-mentioned two crystal structures, deviations of the molecular dipole moments and Coulombic intermolecular interaction energies from the reference values were examined. Root-mean-square deviations (RMSDs) and correlation coefficients were used as a quantitative measure of the quality of the analysed properties. TAAM refinements reproduce X-H bond lengths optimized in theoretical periodic calculations. Structural properties obtained from different database models are similar to each other. The anisotropic displacement parameters from TAAMs are similar to the results of experimental multipolar refinement; differences are about 0.5 and 2.5% for high-resolution and low-resolution data, respectively. Differences in dipole-moment magnitudes calculated from database models are about 5%, and directions differ by up to 30°. The values of electrostatic interaction energies estimated from the individual TAAMs differ greatly from each other and from the reference values. RMSDs are about 9-15 and 22-33 kJ mol(-1) for UBDB and the other database models, respectively.

How to easily replace the independent atom model – the example of bergenin, a potential anti-HIV agent of traditional Asian medicine

Bergenin, which has been isolated from a variety of tropical plants, has several pharmacological applications in traditional Asian medicine. Its electron-density distribution was obtained from a room-temperature low-resolution X-ray data set measured with point detection making use of multipole populations from the invariom library. Two refinement models were considered. In a first step, positional parameters and ADPs were refined with fixed library multipoles (model E1). This model was suitable to be input into a second refinement of multipoles (model E2), which converged smoothly although based on Cu Kalpha room-temperature data. Quantitative results of a topological analysis of the electron density from both models were compared with Hartree-Fock and density-functional calculations. With respect to the independent atom model (IAM) more information can be extracted from invariom modelling, including the electrostatic potential and hydrogen-bond energies, which are highly useful, especially for biologically active compounds. The reliability of the applied invariom formalism was assessed by a comparison of bond-topological properties of sucrose, for which high-resolution multipole and invariom densities were available. Since a conventional X-ray diffraction experiment using basic equipment was combined with the easy-to-use invariom formalism, the procedure described here for bergenin illustrates how it can be routinely applied in pharmacological research.

Electrostatic properties of nine fluoroquinolone antibiotics derived directly from their crystal structure refinements

This article discusses the relevance of the similarity of the molecular electrostatic potential for rational drug design and optimization of lead structures on the basis of a comparison of nine fluoroquinolone antibiotics. Molecular electrostatic potentials and dipole moments were calculated directly from the electron-density distribution after crystal structure refinement with invarioms, theoretical non-spherical scattering factors. Deposited, newly determined and re-measured single-crystal diffraction data of varying quality were evaluated in this manner. An in silico validation procedure for invariom database entries is presented. Despite differences in their substitution pattern, molecular electrostatic potentials for molecules with the same protonation states show a striking degree of similarity for the whole group of compounds studied. The rapid calculation of electrostatic potentials directly from the invariom database makes the procedure suitable for high-throughput screening.