Abstract:The formulation of active pharmaceutical ingredients involves discovering stable crystal packing arrangements, or polymorphs, each of which has distinct pharmaceutically relevant properties. Traditional experimental screening techniques utilizing various conditions are commonly supplemented with in silico crystal structure prediction (CSP) to expedite the crystallization process and mitigate risk. Most predictions are finalized using advanced classical force fields or quantum mechanical calculations that accurately depict the potential energy surface, but do not fully incorporate temperature, pressure or solution conditions during the search. This study proposes an innovative alchemical path that utilizes an advanced polarizable atomic multipole force field to predict crystal structures based on sampling from the NPT ensemble. The use of alchemical (i.e., nonphysical) intermediates, a novel Monte Carlo barostat, and an orthogonal space tempering bias combine to enhance sampling of the deposition/sublimation phase transition and promote an efficient search. The proposed algorithm was applied to 2-((4-(2-(3,4-dichlorophenyl)ethyl)phenyl)amino)benzoic acid (Cambridge Crystallography Database Centre ID: XAFPAY) as a case study to showcase the theory, features and efficiency of this algorithm. All experimentally determined polymorphs with one molecule in the asymmetric unit were successfully reproduced via 1,000 short 1 nsec simulations. Utilizing two threads from a recent Intel® CPU (a Xeon® Gold 6330 CPU at 2.00 GHz) 1 nsec of sampling using the polarizable AMOEBA force field can be acquired in 4 hours (or 336 nsec/day using all 112 threads / 56 cores of the CPU). These results demonstrate a step forward in the rigorous use of the NPT ensemble during a CSP search process and opens the door to future algorithms that incorporate solution conditions using continuum solvation methods.

Rich polymorphism in nicotinamide revealed by melt crystallization and crystal structure prediction

Nicotinamide: Seven New Polymorphic Structures Revealed by Melt Crystallization and Crystal Structure Prediction

A Robust Crystal Structure Prediction Method to Support Small Molecule Drug Development with Large Scale Validation and Prospective Studies

Unprecedented Packing Polymorphism of Oxindole: An Exploration Inspired by Crystal Structure Prediction

Direct Structure Determination from Spherulites using 3D Electron Diffraction

Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs

Insight into the Structure and Property Diversity of Four Nicorandil Oxalate Salt Polymorphs Based on Experiments and Quantum Chemistry Study

Recent Advances in Polymorph Discovery Methods of Organic Crystals

X-Ray Diffraction and Theoretical Calculation-Supported Formation of Polymorphic Cocrystals Discovered Through Thermal Methods: A Case Study

Insight into the Influential Mechanism of Polymorphic Parent Molecule with High Z on the Cocrystal Formation

Polymorphic Solid Solutions in Molecular Crystals: Tips, Tricks, and Switches

Ritonavir Revisited: Melt Crystallization Can Easily Find the Late-Appearing Polymorph II and Unexpectedly Discover a New Polymorph III.

The computational prediction of pharmaceutical crystal structures and polymorphism

Crystal Structure Prediction of Molecular Crystals from First Principles: Are we there yet?

Direct Structure Determination of Vemurafenib Polymorphism from Compact Spherulites Using 3D Electron Diffraction

Understanding Polymorphism in Organic Semiconductor Thin Films Through Nanoconfinement

Crystal Polymorph Search in the NPT Ensemble via a Deposition/Sublimation Alchemical Path

Rationalizing polymorphism with local correlation-based methods: A case study of pnictogen molecular crystals

Polymorphic ROYalty: The 14th ROY Polymorph Discovered via High-Throughput Crystallization

Controlling the crystal polymorph by exploiting the time dependence of nucleation rates