Small-molecule Crystal Structures Research Articles

Molecular replacement for small-molecule crystal structure determination from X-ray and electron diffraction data with reduced resolution.

The resolution of 3D electron diffraction (ED) data of small-molecule crystals is often relatively poor, due to either electron-beam radiation damage during data collection or poor crystallinity of the material. Direct methods, used as standard for crystal structure determination, are not applicable when the data resolution falls below the commonly accepted limit of 1.2 Å. Therefore an evaluation was carried out of the performance of molecular replacement (MR) procedures, regularly used for protein structure determination, for structure analysis of small-molecule crystal structures from 3D ED data. In the course of this study, two crystal structures of Bi-3812, a highly potent inhibitor of the oncogenic transcription factor BCL6, were determined: the structure of α-Bi-3812 was determined from single-crystal X-ray data, the structure of β-Bi-3812 from 3D ED data, using direct methods in both cases. These data were subsequently used for MR with different data types, varying the data resolution limit (1, 1.5 and 2 Å) and by using search models consisting of connected or disconnected fragments of BI-3812. MR was successful with 3D ED data at 2 Å resolution using a search model that represented 74% of the complete molecule.

Reducing errors and increasing accuracy of small-molecule crystal structures

Reducing errors and increasing accuracy of small-molecule crystal structures

The Cambridge structural database (CSD): important resources for teaching concepts in structural chemistry and intermolecular interactions

Abstract The Cambridge Structural Database (CSD) is a repository of all published organic and metal-organic crystal structures of small molecules. These compounds have been crystallized under different conditions and have variable bond parameters and molecular landscapes. Entries in the database therefore serve as real models that can be used to illustrate structural properties such as bond angles, bond distances, torsion angles and other intra- and intermolecular interactions. This paper illustrates how the CSD programs ConQuest and Mercury can be used to search the database for 3D molecular structures to teach concepts such as molecular geometry, symmetry and group theory, organometallic chemistry and intermolecular interactions involving hydrogen bonding and Full Interaction Maps (FIMs) to explore molecular landscapes for halogen bond interactions. Results obtained from these studies are beneficial to understand and predict crystal properties as well as the structural properties of molecules and ions in crystal environments.

Proline C-H Bonds as Loci for Proline Assembly via C-H/O Interactions.

Proline residues within proteins lack a traditional hydrogen bond donor. However, the hydrogens of the proline ring are all sterically accessible, with polarized C-H bonds at Hα and Hδ that exhibit greater partial positive character and can be utilized as alternative sites for molecular recognition. C-H/O interactions, between proline C-H bonds and oxygen lone pairs, have been previously identified as modes of recognition within protein structures and for higher-order assembly of protein structures. In order to better understand intermolecular recognition of proline residues, a series of proline derivatives was synthesized, including 4R-hydroxyproline nitrobenzoate methyl ester, acylated on the proline nitrogen with bromoacetyl and glycolyl groups, and Boc-4S-(4-iodophenyl)hydroxyproline methyl amide. All three derivatives exhibited multiple close intermolecular C-H/O interactions in the crystallographic state, with H⋅⋅⋅O distances as close as 2.3 Å. These observed distances are well below the 2.72 Å sum of the van der Waals radii of H and O, and suggest that these interactions are particularly favorable. In order to generalize these results, we further analyzed the role of C-H/O interactions in all previously crystallized derivatives of these amino acids, and found that all 26 structures exhibited close intermolecular C-H/O interactions. Finally, we analyzed all proline residues in the Cambridge Structural Database of small-molecule crystal structures. We found that the majority of these structures exhibited intermolecular C-H/O interactions at proline C-H bonds, suggesting that C-H/O interactions are an inherent and important mode for recognition of and higher-order assembly at proline residues. Due to steric accessibility and multiple polarized C-H bonds, proline residues are uniquely positioned as sites for binding and recognition via C-H/O interactions.

Multipolar Atom Types from Theory and Statistical Clustering (MATTS) Data Bank: Restructurization and Extension of UBDB.

A fast and accurate operational model of electron density is crucial in many scientific disciplines including crystallography, molecular biology, pharmaceutical, and structural chemistry. In quantum crystallography, the aspherical refinement of crystal structures is becoming increasingly popular because of its accurate description in terms of physically meaningful properties. The transferable aspherical atom model (TAAM) is quick and precise, though it requires a robust algorithm for atom typing and coverage of the most popular atom types present in small organic molecules. Thus, the University at Buffalo Databank (UBDB) has been renamed to the Multipolar Atom Types from Theory and Statistical clustering (MATTS) data bank, broadened, restructured, and implemented into the software DiSCaMB with 651 atom types obtained from 2316 small-molecule crystal structures containing C, H, N, O, P, S, F, Cl, and Br atoms. MATTS2021 data bank now covers most of the small molecules, peptides, RNA, DNA, and some frequently occurring cations and anions in biological, pharmaceutical, and organic materials, including the majority of known crystal structures composed of the above elements. The multipole model parameters (Pval, κ, κ′, Plm) obtained for different atom types were greatly influenced by neighboring atom types, hybridization, geometrical strain in the ring system, and charges on the molecule. Contrary to previous findings, the atoms showing variable oxidation states and ions deviate from the linear dependence of monopole-derived charges on the expansion–contraction κ parameter.

Refinement of RNA Structures Using Amber Force Fields.

Atomic models for nucleic acids derived from X-ray diffraction data at low resolution provide much useful information, but the observed scattering intensities can be fit with models that can differ in structural detail. Tradtional geometric restraints favor models that have bond length and angle terms derived from small molecule crystal structures. Here we explore replacing these restraints with energy gradients derived from force fields, including recently developed integral equation models to account for the effects of water molecules and ions that are not part of the explicit model. We compare conventional and force-field based refinements for 22 RNA crystals, ranging in resolution from 1.1 to 3.6 Å. As expected, it can be important to account for solvent screening of charge–charge interactions, especially in the crowded environment of a nucleic acid crystal. The newly refined models can show improvements in torsion angles and hydrogen-bonding interactions, and can significantly reduce unfavorable atomic clashes, while maintaining or improving agreement with observed scattering intensities.

Development of a DNA damage response (DDR) therapeutics platform for oncology.

e15036 Background: Cancer cells respond to increases in DNA damage by deploying their DNA damage response (DDR) pathways. We are building a platform for the discovery and development of target-specific DDR therapeutics, including small molecule inhibitors and targeted protein degradation warheads, founded on fragment- and structure-based drug discovery. Methods: Our DDR platform, which includes hit-to-lead, lead optimization and candidate selection, starts with hit generation from a new technology that uses high-throughput protein X-ray crystallography to directly screen compound libraries. Our hit generation process produces empirical evidence of direct target engagement. The elucidation of high-quality ligand-bound 3D structures reveals the location and pose of the ligand and details of the protein-ligand interactions. Thus we can predict the structure-function consequences of the hit molecule engagement, which sets the stage for rapid assessment of synthetic tractability and intellectual property. After hit identification, we apply a multi-pronged approach in hit-to-lead conversion and lead optimization using iterative biophysical and biochemical assays, coupled with crystallography. We are applying our approach to several new targets in DDR and will present some early progress in this space. Results: Apurinic/apyrimidinic endonuclease 1 (APE1) is the major repair enzyme for abasic sites in DNA and contributes to DNA strand break processing. Many studies have associated increased APE1 levels with enhanced growth, migration, and drug resistance in human tumor cells, and with decreased patient survival. APE1 has been implicated in over 20 human cancers, including glioblastoma, making the protein an attractive target for the development of anticancer therapeutics. Despite intensive effort, there are no clinical endonuclease inhibitors of APE1. We have identified 25 diverse small molecule fragments that bind to APE1 at two distinct sites, including the endonuclease site. Pol eta (or PolH) is a DNA polymerase implicated, among other things, in the development of cisplatin resistance in a subset of ovarian cancers. In our quest to develop PolH inhibitors, we have identified 5 diverse fragments that bind to two distinct sites in the polymerase including a new potential allosteric site. Our results on APE1 and PolH represent the first known cases of crystal structures of small molecules bound to these proteins. Flap endonuclease (FEN1) is implicated in several cancers including for example ER/tamoxifen-resistant breast cancer. We are developing a targeted protein degradation approach using PROTACs (Proteolysis-Targeting Chimeras) toward the development of novel therapeutics against FEN1. Conclusions: Our results will help us develop small molecule inhibitors and targeted protein degradation against DDR targets that may be effective as single therapies or be used to make existing therapies more effective.

Decisive Influence of Environment on Aromatic/Aromatic Interaction Geometries. Comparison of Aromatic/Aromatic Interactions in Crystal Structures of Small Molecules and in Protein Structures

Geometries of aromatic/aromatic interactions in crystal structures of small molecules from the Cambridge Structural Database (CSD) (benzene/benzene, toluene/toluene, and p-phenol/p-phenol interacti...

First global analysis of the GSK database of small molecule crystal structures

Analysis of the molecular and structural features of the GSK crystal structure database and Cambridge Structural Database leads to improved reliability in hydrogen bond propensity models for pharmaceutical polymorphs.

Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. III. Application to crystal structures via the Ewald and direct summation methods

The previously reported exact potential and multipole moment (EP/MM) method for fast and accurate evaluation of the intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density [Volkov, Koritsanszky & Coppens (2004). Chem. Phys. Lett. 391, 170-175; Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536; Nguyen & Volkov (2019). Acta Cryst. A75, 448-464] is extended to the calculation of electrostatic interaction energies in molecular crystals using two newly developed implementations: (i) the Ewald summation (ES), which includes interactions up to the hexadecapolar level and the EP correction to account for short-range electron-density penetration effects, and (ii) the enhanced EP/MM-based direct summation (DS), which at sufficiently large intermolecular separations replaces the atomic multipole moment approximation to the electrostatic energy with that based on the molecular multipole moments. As in the previous study [Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536], the EP electron repulsion integral is evaluated analytically using the Löwdin α-function approach. The resulting techniques, incorporated in the XDPROP module of the software package XD2016, have been tested on several small-molecule crystal systems (benzene, L-dopa, paracetamol, amino acids etc.) and the crystal structure of a 181-atom decapeptide molecule (Z = 4) using electron densities constructed via the University at Buffalo Aspherical Pseudoatom Databank [Volkov, Li, Koritsanszky & Coppens (2004). J. Phys. Chem. A, 108, 4283-4300]. Using a 2015 2.8 GHz Intel Xeon E3-1505M v5 computer processor, a 64-bit implementation of the Löwdin α-function and one of the higher optimization levels in the GNU Fortran compiler, the ES method evaluates the electrostatic interaction energy with a numerical precision of at least 10-5 kJ mol-1 in under 6 s for any of the tested small-molecule crystal structures, and in 48.5 s for the decapeptide structure. The DS approach is competitive in terms of precision and speed with the ES technique only for crystal structures of small molecules that do not carry a large molecular dipole moment. The electron-density penetration effects, correctly accounted for by the two described methods, contribute 28-64% to the total electrostatic interaction energy in the examined systems, and thus cannot be neglected.

Frequency and hydrogen bonding of nucleobase homopairs in small molecule crystals.

We used the high resolution and accuracy of the Cambridge Structural Database (CSD) to provide detailed information regarding base pairing interactions of selected nucleobases. We searched for base pairs in which nucleobases interact with each other through two or more hydrogen bonds and form more or less planar structures. The investigated compounds were either free forms or derivatives of adenine, guanine, hypoxanthine, thymine, uracil and cytosine. We divided our findings into categories including types of pairs, protonation patterns and whether they are formed by free bases or substituted ones. We found base pair types that are exclusive to small molecule crystal structures, some that can be found only in RNA containing crystal structures and many that are native to both environments. With a few exceptions, nucleobase protonation generally followed a standard pattern governed by pKa values. The lengths of hydrogen bonds did not depend on whether the nucleobases forming a base pair were charged or not. The reasons why particular nucleobases formed base pairs in a certain way varied significantly.

New Insights into 4-Anilinoquinazolines as Inhibitors of Cardiac Troponin I–Interacting Kinase (TNNi3K)

We report the synthesis of several related 4-anilinoquinazolines as inhibitors of cardiac troponin I–interacting kinase (TNNi3K). These close structural analogs of 3-((6,7-dimethoxyquinazolin-4-yl)amino)-4-(dimethylamino)-N-methylbenzenesulfonamide (GSK114) provide new understanding of structure–activity relationships between the 4-anilinoquinazoline scaffold and TNNi3K inhibition. Through a small focused library of inhibitors, we observed that the N-methylbenzenesulfonamide was driving the potency in addition to the more traditional quinazoline hinge-binding motif. We also identified a compound devoid of TNNi3K kinase activity due to the addition of a methyl group in the hinge binding region. This compound could serve as a negative control in the study of TNNi3K biology. Small molecule crystal structures of several quinazolines have been solved, supporting observations made about overall conformation and TNNi3K inhibition.

The impact of the Cambridge Structural Database and the small molecule crystal structures it contains: a bibliographic and literature study

A bibliographic and literature-based analysis of the impact of the Cambridge Structural Database (CSD) and the papers associated with crystal structures in the CSD has been undertaken.

Comprehensive Assessment of Torsional Strain in Crystal Structures of Small Molecules and Protein-Ligand Complexes using ab Initio Calculations.

The energetics of rotation around single bonds (torsions) is a key determinant of the three-dimensional shape that druglike molecules adopt in solution, the solid state, and in different biological environments, which in turn defines their unique physical and pharmacological properties. Therefore, accurate characterization of torsion angle preference and energetics is essential for the success of computational drug discovery and design. Here, we analyze torsional strain in crystal structures of druglike molecules in Cambridge structure database (CSD) and bioactive ligand conformations in protein data bank (PDB), expressing the total strain energy as a sum of strain energy from constituent rotatable bonds. We utilized cloud computing to generate torsion scan profiles of a very large collection of chemically diverse neutral fragments at DFT(B3LYP)/6-31G*//6-31G** or DFT(B3LYP)/6-31+G*//6-31+G** (for sulfur-containing molecule). With the data generated from these ab initio calculations, we performed rigorous analysis of strain due to deviation of observed torsion angles relative to their ideal gas-phase geometries. Contrary to the previous studies based on molecular mechanics, we find that in the crystalline state, molecules generally adopt low-strain conformations, with median per-torsion strain energy in CSD and PDB under one-tenth and one-third of a kcal/mol, respectively. However, for a small fraction (<5%) of motifs, external effects such as steric hindrance and hydrogen bonds result in strain penalty exceeding 2.5 kcal/mol. We find that due to poor quality of PDB structures in general, bioactive structures tend to have higher torsional strain compared to small-molecule crystal conformations. However, in the absence of structural fitting artifacts in PDB structures, protein-induced strain in bioactive conformations is quantitatively similar to those due to the packing forces in small-molecule crystal structures. This analysis allows us to establish strain energy thresholds to help identify biologically relevant conformers in a given ensemble. The work presented here is the most comprehensive study to date that demonstrates the utility and feasibility of gas-phase quantum mechanics (QM) calculations to study conformational preference and energetics of drug-size molecules. Potential applications of this study in computational lead discovery and structure-based design are discussed.

Bridging Crystal Engineering and Drug Discovery by Utilizing Intermolecular Interactions and Molecular Shapes in Crystals.

Most structure-based drug discovery methods utilize crystal structures of receptor proteins. Crystal engineering, on the other hand, utilizes the wealth of chemical information inherent in small-molecule crystal structures in the Cambridge Structural Database (CSD). We show that the interaction surfaces and shapes of molecules in experimentally determined small-molecule crystal structures can serve as effective tools in drug discovery. Our description of the shape and interaction propensities of molecules in their crystal structures can be used to screen them for specific binding compatibility with protein targets, as demonstrated through the high-throughput profiling of around 138 000 small-molecule structures in the CSD and a series of drug-protein crystal structures. Electron-density-based intermolecular boundary surfaces in small-molecule crystal structures and in target-protein pockets are utilized to identify potential ligand molecules from the CSD based on 3D shape and intermolecular interaction matching.

Development of Potent and Selective Pyrazolopyrimidine IRAK4 Inhibitors.

A series of pyrazolopyrimidine inhibitors of IRAK4 were developed from a high-throughput screen (HTS). Modification of an HTS hit led to a series of bicyclic heterocycles with improved potency and kinase selectivity but lacking sufficient solubility to progress in vivo. Structure-based drug design, informed by cocrystal structures with the protein and small-molecule crystal structures, yielded a series of dihydrobenzofurans. This semisaturated bicycle provided superior druglike properties while maintaining excellent potency and selectivity. Improved physicochemical properties allowed for progression into in vivo experiments, where lead molecules exhibited low clearance and showed target-based inhibition of IRAK4 signaling in an inflammation-mediated PK/PD mouse model.

Hunting down hydrogen: applying neutron macromolecular crystallography to galectins

Galectins are a medically important family of small, soluble carbohydrate-binding proteins that have been implicated in a wide variety of diseases, and are therefore targets for the development of new drugs. Many X-ray crystal structures of small molecules in complex with galectins have been determined and have given invaluable information on how natural ligands and synthetic inhibitors are recognized. However, until now the hydrogen-bonding patterns involved in ligand recognition have been a matter of informed guesswork, as X-ray crystallography can't reliably tell us where hydrogen atoms are positioned. We describe how neutron macromolecular crystallography has been applied to galectins for the first time to shed light on this question.

Accurate geometrical restraints for Watson-Crick base pairs.

Geometrical restraints provide key structural information for the determination of biomolecular structures at lower resolution by experimental methods such as crystallography or cryo-electron microscopy. In this work, restraint targets for nucleic acids bases are derived from three different sources and compared: small-molecule crystal structures in the Cambridge Structural Database (CSD), ultrahigh-resolution structures in the Protein Data Bank (PDB) and quantum-mechanical (QM) calculations. The best parameters are those based on CSD structures. After over two decades, the standard library of Parkinson et al. [(1996), Acta Cryst. D52, 57-64] is still valid, but improvements are possible with the use of the current CSD database. The CSD-derived geometry is fully compatible with Watson-Crick base pairs, as comparisons with QM results for isolated and paired bases clearly show that the CSD targets closely correspond to proper base pairing. While the QM results are capable of distinguishing between single and paired bases, their level of accuracy is, on average, nearly two times lower than for the CSD-derived targets when gauged by root-mean-square deviations from ultrahigh-resolution structures in the PDB. Nevertheless, the accuracy of QM results appears sufficient to provide stereochemical targets for synthetic base pairs where no reliable experimental structural information is available. To enable future tests for this approach, QM calculations are provided for isocytosine, isoguanine and the iCiG base pair.

Pushing the Limits of Molecular Crystal Structure Determination From Powder Diffraction Data in High-Throughput Chemical Environments

Crystal structure determination from powder diffraction data (SDPD) using the DASH software package is evaluated for data recorded using transmission capillary, transmission flat plate, and reflection flat plate geometries on a selection of pharmaceutical compounds. We show that transmission capillary geometry remains the best option when crystal structure determination is the primary consideration and, as expected, reflection flat plate geometry is not recommended for SDPD because of preferred orientation effects. However, the quality of crystal structures obtained from transmission plate instruments can be excellent, and the convenience factor for sample preparation, throughput, and retrieval is higher than that of transmission capillary instruments. Indeed, it is possible to solve crystal structures within an hour of a polycrystalline sample arriving in the laboratory, which has clear implications for making small-molecule crystal structures more routinely available to the practicing laboratory medicinal chemist. With appropriate modifications to crystal structure determination software, it can be imagined that SDPD could become a rapid turn-around walk-up analytical service in high-throughput chemical environments.

Using more than 801 296 small-molecule crystal structures to aid in protein structure refinement and analysis. Addendum.

An addendum to the Introduction of Cole et al. [(2017), Acta Cryst. D73, 234-239] is made to recognize the work of Bricogne, Smart and others in the development of methods to make use of Cambridge Structural Database data in protein structure solution.