Bernd Schweizer

Significant progress in predicting the crystal structures of small organic molecules – a report on the fourth blind test

We report on the organization and outcome of the fourth blind test of crystal structure prediction, an international collaborative project organized to evaluate the present state in computational methods of predicting the crystal structures of small organic molecules. There were 14 research groups which took part, using a variety of methods to generate and rank the most likely crystal structures for four target systems: three single-component crystal structures and a 1:1 cocrystal. Participants were challenged to predict the crystal structures of the four systems, given only their molecular diagrams, while the recently determined but as-yet unpublished crystal structures were withheld by an independent referee. Three predictions were allowed for each system. The results demonstrate a dramatic improvement in rates of success over previous blind tests; in total, there were 13 successful predictions and, for each of the four targets, at least two groups correctly predicted the observed crystal structure. The successes include one participating group who correctly predicted all four crystal structures as their first ranked choice, albeit at a considerable computational expense. The results reflect important improvements in modelling methods and suggest that, at least for the small and fairly rigid types of molecules included in this blind test, such calculations can be constructively applied to help understand crystallization and polymorphism of organic molecules.

A third blind test of crystal structure prediction

Following the interest generated by two previous blind tests of crystal structure prediction (CSP1999 and CSP2001), a third such collaborative project (CSP2004) was hosted by the Cambridge Crystallographic Data Centre. A range of methodologies used in searching for and ranking the likelihood of predicted crystal structures is represented amongst the 18 participating research groups, although most are based on the global minimization of the lattice energy. Initially the participants were given molecular diagrams of three molecules and asked to submit three predictions for the most likely crystal structure of each. Unlike earlier blind tests, no restriction was placed on the possible space group of the target crystal structures. Furthermore, Z = 2 structures were allowed. Part-way through the test, a partial structure report was discovered for one of the molecules, which could no longer be considered a blind test. Hence, a second molecule from the same category (small, rigid with common atom types) was offered to the participants as a replacement. Success rates within the three submitted predictions were lower than in the previous tests - there was only one successful prediction for any of the three ;blind molecules. For the ;simplest rigid molecule, this lack of success is partly due to the observed structure crystallizing with two molecules in the asymmetric unit. As in the 2001 blind test, there was no success in predicting the structure of the flexible molecule. The results highlight the necessity for better energy models, capable of simultaneously describing conformational and packing energies with high accuracy. There is also a need for improvements in search procedures for crystals with more than one independent molecule, as well as for molecules with conformational flexibility. These are necessary requirements for the prediction of possible thermodynamically favoured polymorphs. Which of these are actually realised is also influenced by as yet insufficiently understood processes of nucleation and crystal growth.

Crystal structure prediction of small organic molecules: a second blind test.

The first collaborative workshop on crystal structure prediction (CSP1999) has been followed by a second workshop (CSP2001) held at the Cambridge Crystallographic Data Centre. The 17 participants were given only the chemical diagram for three organic molecules and were invited to test their prediction programs within a range of named common space groups. Several different computer programs were used, using the methodology wherein a molecular model is used to construct theoretical crystal structures in given space groups, and prediction is usually based on the minimum calculated lattice energy. A maximum of three predictions were allowed per molecule. The results showed two correct predictions for the first molecule, four for the second molecule and none for the third molecule (which had torsional flexibility). The correct structure was often present in the sorted low-energy lists from the participants but at a ranking position greater than three. The use of non-indexed powder diffraction data was investigated in a secondary test, after completion of the ab initio submissions. Although no one method can be said to be completely reliable, this workshop gives an objective measure of the success and failure of current methodologies.

Cleft-Type Diamidinium Receptors for Dicarboxylate Binding in Protic Solvents

A series of potential cleft-type receptors for dicarboxylate substrates were prepared by attachment of two phenylamidinium ions to either naphthalene or 1,1′-binaphthalene scaffolds. Their synthesis (Schemesu20051u2009–u20094) involved the Pd0-catalyzed cross-coupling of aryl nitriles to the central scaffold, followed by transformation of the nitrile into amidinium groups using the Garigipati reaction. The 1,1′-binaphthalene derivative (±)-1 with phenylamidinium residues attached to the 6,6′-positions in the major groove was found to be a highly efficient receptor for dicarboxylate guests, such as glutarate and isophthalates, even in competing protic solvents such as CD3OD (Tableu20051). The vant Hoff analysis of variable-temperature 1H-NMR (VT-NMR) titrations (Tableu20052 and Fig.u20053) and isothermal microcalorimetry (ITC; Tableu20053 and Fig.u20054) revealed that complexation in MeOH is strongly entropically driven with an unfavorable enthalpic change, which partially compensates the entropic gain. These thermodynamic quantities are best explained by a particularly favorable solvation of the binding partners in the unbound state and the release of the MeOH molecules, which solvate the free ions into the bulk upon complexation. Receptor (±)-1 binds flexible glutarate and rigid isophthalates with similar association strength. This lack in response to guest preorganization and reduced guest selectivity is explained with the non-directionality of the coulombic charge-charge interactions in the complexes.

Deorphaning Pyrrolopyrazines as Potent Multi‐Target Antimalarial Agents

The discovery of pyrrolopyrazines as potent antimalarial agents is presented, with the most effective compounds exhibiting EC50 values in the low nanomolar range against asexual blood stages of Plasmodium falciparum in human red blood cells, and Plasmodium berghei liver schizonts, with negligible HepG2 cytotoxicity. Their potential mode of action is uncovered by predicting macromolecular targets through avant-garde computer modeling. The consensus prediction method suggested a functional resemblance between ligand binding sites in non-homologous target proteins, linking the observed parasite elimination to IspD, an enzyme from the non-mevalonate pathway of isoprenoid biosynthesis, and multi-kinase inhibition. Further computational analysis suggested essential P. falciparum kinases as likely targets of our lead compound. The results obtained validate our methodology for ligand- and structure-based target prediction, expand the bioinformatics toolbox for proteome mining, and provide unique access to deciphering polypharmacological effects of bioactive chemical agents.