Predicting Inclusion Behaviour and Framework Structures in Organic Crystals

Engineering Lipophilic Aggregation of Adapalene and Adamantane-Based Cocrystals via van der Waals Forces and Hydrogen Bonding

Lipophilic aggregation using adamantanes is a widely exploited molecular property in medicinal and materials chemistry. Adamantanes are traditionally installed to molecular units via covalent bonds. However, the noncovalent installation of adamantanes has been relatively underexplored and presents the potential to bring properties associated with adamantanes to molecules without affecting their intrinsic properties (e.g., pharmacophores). Here, we systematically study a series of adamantanecarboxylic acids with varying substitution levels of methyl groups and their cocrystals with bipyridines. Specifically, single-crystal X-ray diffraction shows that while the directionality of single-component adamantanes is notably sensitive to changes in methyl substitution, hydrogen-bonded cocrystals with bipyridines show consistent and robust packing due to π-stacking predominance. Our observations are supported by Hirshfeld surface and energy framework analyses. The applicability of cocrystal formation of adamantanes bearing carboxylic acids was used to generate the first cocrystals of adapalene, an adamantane-bearing retinoid used for treating acne vulgaris. We envisage our study to inspire noncovalent (i.e., cocrystal) installation of adamantanes to generate lipophilic aggregation in multicomponent systems.

A data-driven and topological mapping approach for the a priori prediction of stable molecular crystalline hydrates

Predictions of the structures of stoichiometric, fractional, or nonstoichiometric hydrates of organic molecular crystals are immensely challenging due to the extensive search space of different water contents, host molecular placements throughout the crystal, and internal molecular conformations. However, the dry frameworks of these hydrates, especially for nonstoichiometric or isostructural dehydrates, can often be predicted from a standard anhydrous crystal structure prediction (CSP) protocol. Inspired by developments in the field of drug binding, we introduce an efficient data-driven and topologically aware approach for predicting organic molecular crystal hydrate structures through a mapping of water positions within the crystal structure. The method does not require a priori specification of water content and can, therefore, predict stoichiometric, fractional, and nonstoichiometric hydrate structures. This approach, which we term a mapping approach for crystal hydrates (MACH), establishes a set of rules for systematic determination of favorable positions for water insertion within predicted or experimental crystal structures based on considerations of the chemical features of local environments and void regions. The proposed approach is tested on hydrates of three pharmaceutically relevant compounds that exhibit diverse crystal packing motifs and void environments characteristic of hydrate structures. Overall, we show that our mapping approach introduces an advance in the efficient performance of hydrate CSP through generation of stable hydrate stoichiometries at low cost and should be considered an integral component for CSP workflows.

Minimizing Polymorphic Risk through Cooperative Computational and Experimental Exploration

We combine state-of-the-art computational crystal structure prediction (CSP) techniques with a wide range of experimental crystallization methods to understand and explore crystal structure in pharmaceuticals and minimize the risk of unanticipated late-appearing polymorphs. Initially, we demonstrate the power of CSP to rationalize the difficulty in obtaining polymorphs of the well-known pharmaceutical isoniazid and show that CSP provides the structure of the recently obtained, but unsolved, Form III of this drug despite there being only a single resolved form for almost 70 years. More dramatically, our blind CSP study predicts a significant risk of polymorphism for the related iproniazid. Employing a wide variety of experimental techniques, including high-pressure experiments, we experimentally obtained the first three known nonsolvated crystal forms of iproniazid, all of which were successfully predicted in the CSP procedure. We demonstrate the power of CSP methods and free energy calculations to rationalize the observed elusiveness of the third form of iproniazid, the success of high-pressure experiments in obtaining it, and the ability of our synergistic computational-experimental approach to “de-risk” solid form landscapes.

Influence of crystal packing on the thermal properties of cocrystals and cocrystal solvates of olanzapine: insights from computations

A multicomponent supramolecular host with adaptive guest accommodation abilities is observed in the cocrystal solvates of the olanzapine–hydroquinone system.

Supramolecular organisation of sulphate salt hydrates exemplified with brucine sulphate

The frequency of hydrate formation among organic sulphate salts is unravelled. Interconversion of the hydrates of brucine sulphate occurs with small changes in the relative humidity.

Computational modelling of solvent effects in a prolific solvatomorphic porous organic cage

Crystal structure prediction methods can enable the in silico design of functional molecular crystals, but solvent effects can have a major influence on relative lattice energies, sometimes thwarting predictions. This is particularly true for porous solids, where solvent included in the pores can have an important energetic contribution. We present a Monte Carlo solvent insertion procedure for predicting the solvent filling of porous structures from crystal structure prediction landscapes, tested using a highly solvatomorphic porous organic cage molecule, CC1. Using this method, we can understand why the predicted global energy minimum structure for CC1 is never observed from solvent crystallisation. We also explain the formation of three different solvatomorphs of CC1 from three structurally-similar chlorinated solvents. Calculated solvent stabilisation energies are found to correlate with experimental results from thermogravimetric analysis, suggesting a future computational framework for a priori materials design that factors in solvation effects.

Energy-Structure-Function Maps: Cartography for Materials Discovery

Some of the most successful approaches to structural design in materials chemistry have exploited strong directional bonds, whose geometric reliability lends predictability to solid-state assembly. For example, metal-organic frameworks are an important design platform in materials chemistry. By contrast, the structure of molecular crystals is defined by a balance of weaker intermolecular forces, and small changes to the molecular building blocks can lead to large changes in crystal packing. Hence, empirical rules are inherently less reliable for engineering the structures of molecular solids. Energy-structure-function (ESF) maps are a new approach for the discovery of functional organic crystals. These maps fuse crystal-structure prediction with the computation of physical properties to allow researchers to choose the most promising molecule for a given application, prior to its synthesis. ESF maps were used recently to discover a highly porous molecular crystal that has a high methane deliverable capacity and the lowest density molecular crystal reported to date (r = 0.41 g cm^-3 , SA_BET = 3425 m² g^-1 ). Progress in this field is reviewed, with emphasis on the future opportunities and challenges for a design strategy based on computed ESF maps.

From serendipity to supramolecular design: assessing the utility of computed crystal form landscapes in inferring the risks of crystal hydration in carboxylic acids

Calculated structural descriptors for predicted anhydrate polymorphs are used to assess the risks of crystal hydration in carboxylic acids.

Unraveling Complexity in the Solid Form Screening of a Pharmaceutical Salt: Why so Many Forms? Why so Few?

The solid form landscape of 5-HT2a antagonist 3-(4-(benzo[d]isoxazole-3-yl)piperazin-1-yl)-2,2-dimethylpropanoic acid hydrochloride (B5HCl) proved difficult to establish. Many crystalline materials were produced by solid form screening, but few forms readily grew high quality crystals to afford a clear picture or understanding of the solid form landscape. Careful control of crystallization conditions, a range of experimental methods, computational modeling of solvate structures, and crystal structure prediction were required to see potential arrangements of the salt in its crystal forms. Structural diversity in the solid form landscape of B5HCl was apparent in the layer structures for the anhydrate polymorphs (Forms I and II), dihydrate and a family of solvates with alcohols. The alcohol solvates, which provided a distinct packing from the neat forms and the dihydrate, form layers with conserved hydrogen bonding between B5HCl and the solvent, as well as stacking of the aromatic rings. The ability of the alcohol hydrocarbon moieties to efficiently pack between the layers accounted for the difficulty in growing some solvate crystals and the inability of other solvates to crystallize altogether. Through a combination of experiment and computation, the crystallization problems, form stability, and desolvation pathways of B5HCl have been rationalized at a molecular level.

Computationally-Guided Synthetic Control over Pore Size in Isostructural Porous Organic Cages

The physical properties of 3-D porous solids are defined by their molecular geometry. Hence, precise control of pore size, pore shape, and pore connectivity are needed to tailor them for specific applications. However, for porous molecular crystals, the modification of pore size by adding pore-blocking groups can also affect crystal packing in an unpredictable way. This precludes strategies adopted for isoreticular metal-organic frameworks, where addition of a small group, such as a methyl group, does not affect the basic framework topology. Here, we narrow the pore size of a cage molecule, CC3, in a systematic way by introducing methyl groups into the cage windows. Computational crystal structure prediction was used to anticipate the packing preferences of two homochiral methylated cages, CC14-R and CC15-R, and to assess the structure-energy landscape of a CC15-R/CC3-S cocrystal, designed such that both component cages could be directed to pack with a 3-D, interconnected pore structure. The experimental gas sorption properties of these three cage systems agree well with physical properties predicted by computational energy-structure-function maps.

Functional materials discovery using energy-structure-function maps

Molecular crystals cannot be designed in the same manner as macroscopic objects, because they do not assemble according to simple, intuitive rules. Their structures result from the balance of many weak interactions, rather than from the strong and predictable bonding patterns found in metal-organic frameworks and covalent organic frameworks. Hence, design strategies that assume a topology or other structural blueprint will often fail. Here we combine computational crystal structure prediction and property prediction to build energy-structure-function maps that describe the possible structures and properties that are available to a candidate molecule. Using these maps, we identify a highly porous solid, which has the lowest density reported for a molecular crystal so far. Both the structure of the crystal and its physical properties, such as methane storage capacity and guest-molecule selectivity, are predicted using the molecular structure as the only input. More generally, energy-structure-function maps could be used to guide the experimental discovery of materials with any target function that can be calculated from predicted crystal structures, such as electronic structure or mechanical properties.

Clathrate Structure Determination by Combining Crystal Structure Prediction with Computational and Experimental 129 Xe NMR Spectroscopy

An approach is presented for the structure determination of clathrates using NMR spectroscopy of enclathrated xenon to select from a set of predicted crystal structures. Crystal structure prediction methods have been used to generate an ensemble of putative structures of o- and m-fluorophenol, whose previously unknown clathrate structures have been studied by ¹²⁹ Xe NMR spectroscopy. The high sensitivity of the ¹²⁹ Xe chemical shift tensor to the chemical environment and shape of the crystalline cavity makes it ideal as a probe for porous materials. The experimental powder NMR spectra can be used to directly confirm or reject hypothetical crystal structures generated by computational prediction, whose chemical shift tensors have been simulated using density functional theory. For each fluorophenol isomer one predicted crystal structure was found, whose measured and computed chemical shift tensors agree within experimental and computational error margins and these are thus proposed as the true fluorophenol xenon clathrate structures.

Towards stable porous crystalline phases of molecular belts

A chemical modification of a molecular belt (M1) renders the molecule (M2) into a stable supramolecular nanotube porous crystal.

Application of computational methods to the design and characterisation of porous molecular materials

Composed from discrete units, porous molecular materials (PMMs) possess properties not observed for conventional, extended solids. Molecular simulations provide crucial understanding for the design and characterisation of these unique materials.

Why do Hydrates (Solvates) Form in Small Neutral Organic Molecules? Exploring the Crystal Form Landscapes of the Alkaloids Brucine and Strychnine

Computational methods were used to generate and explore the crystal structure landscapes of the two alkaloids strychnine and brucine. The computed structures were analyzed and rationalized by correlating the modelling results to a rich pool of available experimental data. Despite their structural similarity, the two compounds show marked differences in the formation of solid forms. For strychnine only one anhydrous form is reported in the literature and two new solvates from 1,4-dioxane were detected in the course of this work. In contrast, 22 solid forms are so far known to exist for brucine, comprising two anhydrates, four hydrates (HyA - HyC and a 5.25-hydrate), twelve solvates (alcohols and acetone) and four heterosolvates (mixed solvates with water and alcohols). For strychnine it is hard to produce any solid form other than the stable anhydrate while the formation of specific solid state forms of brucine is governed by a complex interplay between temperature and relative humidity/water activity and it is rather a challenging to avoid hydrate formation. Differences in crystal packing and the high tendency for brucine to form hydrates are not intuitive from the molecular structure alone, as both molecules have hydrogen bond acceptor groups but lack hydrogen bond donor groups. Only the evaluation of the crystal energy landscapes, in particular the close-packed crystal structures and high-energy open frameworks containing voids of molecular (water) dimensions, allowed us to unravel the diverse solid state behavior of the two alkaloids at a molecular level. In this study we demonstrate that expanding the analysis of anhydrate crystal energy landscapes to higher energy structures and calculating the solvent-accessible volume can be used to estimate non-stoichiometric or channel hydrate (solvate) formation, without explicitly computing the hydrate/solvate crystal energy landscapes.

Crystalline arrays of side chain modified bile acids derivatives. Two novel self-assemblies based on π-π and belly-to-belly interactions

Crystalline derivatives of side chain modified bile acids were efficiently prepared from the naturally occurring steroids by palladium-catalyzed cross coupling reaction as a key step. The solvent-free crystalline bile acids derivatives 2b-e are readily accessed by slow evaporation from selected solvents. A variety of steroidal scaffolds were found and elucidated by SXRD studies. The crystal packing of the title compounds are dominated by hydrogen-bonding interactions established between differently positioned acetyl protecting groups, which in the case of 2b and 2e take advantage of the facial amphiphilicity producing two novel steroidal supramolecular self-assemblies combining π-π and strong facial interactions. Thus, these crystalline arrays of side chain modified bile acids represent promising scaffolds for research and implementation in biomolecular materials or inclusion phenomena.

Using crystal structure prediction to rationalize the hydration propensities of substituted adamantane hydrochloride salts

The crystal energy landscapes of the salts of two rigid pharmaceutically active molecules reveal that the experimental structure of amantadine hydrochloride is the most stable structure with the majority of low-energy structures adopting a chain hydrogen-bond motif and packings that do not have solvent accessible voids. By contrast, memantine hydrochloride which differs in the substitution of two methyl groups on the adamantane ring has a crystal energy landscape where all structures within 10 kJ mol−1of the global minimum have solvent-accessible voids ranging from 3 to 14% of the unit-cell volume including the lattice energy minimum that was calculated after removing water from the hydrated memantine hydrochloride salt structure. The success in using crystal structure prediction (CSP) to rationalize the different hydration propensities of these substituted adamantane hydrochloride salts allowed us to extend the model to predict under blind test conditions the experimental crystal structures of the previously uncharacterized 1-(methylamino)adamantane base and its corresponding hydrochloride salt. Although the crystal structure of 1-(methylamino)adamantane was correctly predicted as the second ranked structure on the static lattice energy landscape, the crystallization of aZ′ = 3 structure of 1-(methylamino)adamantane hydrochloride reveals the limits of applying CSP when the contents of the crystallographic asymmetric unit are unknown.

Can computed crystal energy landscapes help understand pharmaceutical solids?

Computational crystal structure prediction (CSP) methods can now be applied to the smaller pharmaceutical molecules currently in drug development. We review the recent uses of computed crystal energy landscapes for pharmaceuticals, concentrating on examples where they have been used in collaboration with industrial-style experimental solid form screening. There is a strong complementarity in aiding experiment to find and characterise practically important solid forms and understanding the nature of the solid form landscape.

Computational and Experimental Characterization of Five Crystal Forms of Thymine: Packing Polymorphism, Polytypism/Disorder and Stoichiometric 0.8-Hydrate

New polymorphs of thymine emerged in an experimental search for solid forms, which was guided by the computationally generated crystal energy landscape. Three of the four anhydrates (AH) are homeoenergetic (A° - C) and their packing modes differ only in the location of oxygen and hydrogen atoms. AHs A° and B are ordered phases, whereas AH C shows disorder (X-ray diffuse scattering). Anhydrates AHs A° and B are ordered phases, whereas AH C shows disorder (X-ray diffuse scattering). Analysis of the crystal energy landscape for alternative AH C hydrogen bonded ribbon motifs identified a number of different packing modes, whose 3D structures were calculated to deviate by less than 0.24 kJ mol-1 in lattice energy. These structures provide models for stacking faults. The three anhydrates A° - C show strong similarity in their powder X-ray diffraction, thermoanalytical and spectroscopic (IR and Raman) characteristics. The already known anhydrate AH A° was identified as the thermodynamically most stable form at ambient conditions; AH B and AH C are metastable but show high kinetic stability. The hydrate of thymine is stable only at water activities (aw) > 0.95 at temperatures ≤ 25 °C. It was found to be a stoichiometric hydrate despite being a channel hydrate with an unusual water:thymine ratio of 0.8:1. Depending on the dehydration conditions, either AH C or AH D is obtained. The hydrate is the only known precursor to AH D. This study highlights the value and complementarity of simultaneous explorations of computationally and experimentally generated solid form landscapes of a small molecule anhydrate ↔ hydrate system.

Convergence Properties of Crystal Structure Prediction by Quasi-Random Sampling

Generating sets of trial structures that sample the configurational space of crystal packing possibilities is an essential step in the process of ab initio crystal structure prediction (CSP). One effective methodology for performing such a search relies on low-discrepancy, quasi-random sampling, and our implementation of such a search for molecular crystals is described in this paper. Herein we restrict ourselves to rigid organic molecules and, by considering their geometric properties, build trial crystal packings as starting points for local lattice energy minimization. We also describe a method to match instances of the same structure, which we use to measure the convergence of our packing search toward completeness. The use of these tools is demonstrated for a set of molecules with diverse molecular characteristics and as representative of areas of application where CSP has been applied. An important finding is that the lowest energy crystal structures are typically located early and frequently during a quasi-random search of phase space. It is usually the complete sampling of higher energy structures that requires extended sampling. We show how the procedure can first be refined, through targetting the volume of the generated crystal structures, and then extended across a range of space groups to make a full CSP search and locate experimentally observed and lists of hypothetical polymorphs. As the described method has also been created to lie at the base of more involved approaches to CSP, which are being developed within the Global Lattice Energy Explorer (Glee) software, a few of these extensions are briefly discussed.

Solid state forms of 4-aminoquinaldine - From void structures with and without solvent inclusion to close packing

Polymorphs of 4-aminoquinaldine (4-AQ) have been predicted in silico and experimentally identified and characterised. The two metastable forms, AH (anhydrate) II and AH III, crystallise in the trigonal space group [Formula: see text] and are less densely packed than the thermodynamically most stable phase AH I° (P21/c ). AH II can crystallise and exist both, as a solvent inclusion compound and as an unsolvated phase. The third polymorph, AH III, is exclusively obtained by desolvation of a carbon tetrachloride solvate. Theoretical calculations correctly estimated the experimental 0K stability order, confirmed that AH II can exist without solvents, gave access to the AH III structure, and identified that there exists a subtle balance between close packing and number of hydrogen bonding interactions in the solid state of anhydrous 4-AQ. Furthermore, the prevalence of void space and solvent inclusion in [Formula: see text] structures is discussed.

Insights into hydrate formation and stability of morphinanes from a combination of experimental and computational approaches

Morphine, codeine, and ethylmorphine are important drug compounds whose free bases and hydrochloride salts form stable hydrates. These compounds were used to systematically investigate the influence of the type of functional groups, the role of water molecules, and the Cl^– counterion on molecular aggregation and solid state properties. Five new crystal structures have been determined. Additionally, structure models for anhydrous ethylmorphine and morphine hydrochloride dihydrate, two phases existing only in a very limited humidity range, are proposed on the basis of computational dehydration modeling. These match the experimental powder X-ray diffraction patterns and the structural information derived from infrared spectroscopy. All 12 structurally characterized morphinane forms (including structures from the Cambridge Structural Database) crystallize in the orthorhombic space group P2₁2₁2₁. Hydrate formation results in higher dimensional hydrogen bond networks. The salt structures of the different compounds exhibit only little structural variation. Anhydrous polymorphs were detected for all compounds except ethylmorphine (one anhydrate) and its hydrochloride salt (no anhydrate). Morphine HCl forms a trihydrate and dihydrate. Differential scanning and isothermal calorimetry were employed to estimate the heat of the hydrate ↔ anhydrate phase transformations, indicating an enthalpic stabilization of the respective hydrate of 5.7 to 25.6 kJ mol^–1 relative to the most stable anhydrate. These results are in qualitative agreement with static 0 K lattice energy calculations for all systems except morphine hydrochloride, showing the need for further improvements in quantitative thermodynamic prediction of hydrates having water···water interactions. Thus, the combination of a variety of experimental techniques, covering temperature- and moisture-dependent stability, and computational modeling allowed us to generate sufficient kinetic, thermodynamic and structural information to understand the principles of hydrate formation of the model compounds. This approach also led to the detection of several new crystal forms of the investigated morphinanes.

Which conformations make stable crystal structures? Mapping crystalline molecular geometries to the conformational energy landscape

Flexible organic molecules often do not adopt their lowest energy conformer in crystal structures. We find that there is a preference for molecules to crystallise with high surface area conformers.

Predicted crystal energy landscapes of porous organic cages

Computational methods predict the crystal packing of porous organic cage molecules, allowing crystal structure and porosity to be predicted starting from the chemical diagram alone.

Tunable porous organic crystals: structural scope and adsorption properties of nanoporous steroidal ureas

Previous work has shown that certain steroidal bis-(N-phenyl)ureas, derived from cholic acid, form crystals in the P6(1) space group with unusually wide unidimensional pores. A key feature of the nanoporous steroidal urea (NPSU) structure is that groups at either end of the steroid are directed into the channels and may in principle be altered without disturbing the crystal packing. Herein we report an expanded study of this system, which increases the structural variety of NPSUs and also examines their inclusion properties. Nineteen new NPSU crystal structures are described, to add to the six which were previously reported. The materials show wide variations in channel size, shape, and chemical nature. Minimum pore diameters vary from ~0 up to 13.1 Å, while some of the interior surfaces are markedly corrugated. Several variants possess functional groups positioned in the channels with potential to interact with guest molecules. Inclusion studies were performed using a relatively accessible tris-(N-phenyl)urea. Solvent removal was possible without crystal degradation, and gas adsorption could be demonstrated. Organic molecules ranging from simple aromatics (e.g., aniline and chlorobenzene) to the much larger squalene (M(w) = 411) could be adsorbed from the liquid state, while several dyes were taken up from solutions in ether. Some dyes gave dichroic complexes, implying alignment of the chromophores in the NPSU channels. Notably, these complexes were formed by direct adsorption rather than cocrystallization, emphasizing the unusually robust nature of these organic molecular hosts.

Absorbing a Little Water: The Structural, Thermodynamic, and Kinetic Relationship between Pyrogallol and Its Tetarto-Hydrate

The anhydrate and the stoichiometric tetarto-hydrate of pyrogallol (0.25 mol water per mol pyrogallol) are both storage stable at ambient conditions, provided that they are phase pure, with the system being at equilibrium at aw (water activity) = 0.15 at 25 °C. Structures have been derived from single crystal and powder X-ray diffraction data for the anhydrate and hydrate, respectively. It is notable that the tetarto-hydrate forms a tetragonal structure with water in channels, a framework that although stabilized by water, is found as a higher energy structure on a computationally generated crystal energy landscape, which has the anhydrate crystal structure as the most stable form. Thus, a combination of slurry experiments, X-ray diffraction, spectroscopy, moisture (de)sorption, and thermo-analytical methods with the computationally generated crystal energy landscape and lattice energy calculations provides a consistent picture of the finely balanced hydration behavior of pyrogallol. In addition, two monotropically related dimethyl sulfoxide monosolvates were found in the accompanying solid form screen.

Why don't we find more polymorphs?

Crystal structure prediction (CSP) studies are not limited to being a search for the most thermodynamically stable crystal structure, but play a valuable role in understanding polymorphism, as shown by interdisciplinary studies where the crystal energy landscape has been explored experimentally and computationally. CSP usually produces more thermodynamically plausible crystal structures than known polymorphs. This article illustrates some reasons why: because (i) of approximations in the calculations, particularly the neglect of thermal effects (see §1.1); (ii) of the molecular rearrangement during nucleation and growth (see §1.2); (iii) the solid-state structures observed show dynamic or static disorder, stacking faults, other defects or are not crystalline and so represent more than one calculated structure (see §1.3); (iv) the structures are metastable relative to other molecular compositions (see §1.4); (v) the right crystallization experiment has not yet been performed (see §1.5) or (vi) cannot be performed (see §1.6) and the possibility (vii) that the polymorphs are not detected or structurally characterized (see §1.7). Thus, we can only aspire to a general predictive theory for polymorphism, as this appears to require a quantitative understanding of the kinetic factors involved in all possible multi-component crystallizations. For a specific molecule, analysis of the crystal energy landscape shows the potential complexity of its crystallization behaviour.

In silico design of supramolecules from their precursors: odd-even effects in cage-forming reactions

We synthesize a series of imine cage molecules where increasing the chain length of the alkanediamine precursor results in an odd-even alternation between [2 + 3] and [4 + 6] cage macrocycles. A computational procedure is developed to predict the thermodynamically preferred product and the lowest energy conformer, hence rationalizing the observed alternation and the 3D cage structures, based on knowledge of the precursors alone.

Successful prediction of a model pharmaceutical in the fifth blind test of crystal structure prediction

The range of target structures in the fifth international blind test of crystal structure prediction was extended to include a highly flexible molecule, (benzyl-(4-(4-methyl-5-(p-tolylsulfonyl)-1,3-thiazol-2-yl)phenyl)carbamate, as a challenge representative of modern pharmaceuticals. Two of the groups participating in the blind test independently predicted the correct structure. The methods they used are described and contrasted, and the implications of the capability to tackle molecules of this complexity are discussed.