RESUMEN
The comparative crystallizability and polymorphic selectivity of ritonavir, a novel protease inhibitor for the treatment of acquired immune-deficiency syndrome, as a function of solvent selection are examined through an integrated and self-consistent experimental and computational molecular modeling study. Recrystallization at high supersaturation by rapid cooling at 283.15 K is found to produce the metastable "disappeared" polymorphic form I from acetone, ethyl acetate, acetonitrile, and toluene solutions in contrast to ethanol which produces the stable form II. Concomitant crystallization of the other known solid forms is not found under these conditions. Isothermal crystallization studies using turbidometric detection based upon classical nucleation theory reveal that, for an equal induction time, the required driving force needed to initiate solution nucleation decreases with solubility in the order of ethanol, acetone, acetonitrile, ethyl acetate, and toluene consistent with the expected desolvation behavior predicted from the calculated solute solvation free energies. Molecular dynamics simulations of the molecular and intermolecular chemistry reveal the presence of conformational interplay between intramolecular and intermolecular interactions within the solution phase. These encompass the solvent-dependent formation of intramolecular O-H...O hydrogen bonding between the hydroxyl and carbamate groups coupled with differing conformations of the hydroxyl's shielding phenyl groups. These conformational preferences and their relative interaction propensities, as a function of solvent selection, may play a rate-limiting role in the crystallization behavior by not only inhibiting to different degrees the nucleation process but also restricting the assembly of the optimal intermolecular hydrogen bonding network needed for the formation of the stable form II polymorph.
Asunto(s)
Cristalización , Enlace de Hidrógeno , Simulación de Dinámica Molecular , Ritonavir , Solventes , Ritonavir/química , Solventes/química , Solubilidad , Etanol/química , Acetatos , AcetonitrilosRESUMEN
Precision measurement of the growth rate of individual single crystal facets (hkl) represents an important component in the design of industrial crystallization processes. Current approaches for crystal growth measurement using optical microscopy are labor intensive and prone to error. An automated process using state-of-the-art computer vision and machine learning to segment and measure the crystal images is presented. The accuracies and efficiencies of the new crystal sizing approach are evaluated against existing manual and semi-automatic methods, demonstrating equivalent accuracy but over a much shorter time, thereby enabling a more complete kinematic analysis of the overall crystallization process. This is applied to measure in situ the crystal growth rates and through this determining the associated kinetic mechanisms for the crystallization of ß-form l-glutamic acid from the solution phase. Growth on the {101} capping faces is consistent with a Birth and Spread mechanism, in agreement with the literature, while the growth rate of the {021} prismatic faces, previously not available in the literature, is consistent with a Burton-Cabrera-Frank screw dislocation mechanism. At a typical supersaturation of σ = 0.78, the growth rate of the {101} capping faces (3.2 × 10-8 m s-1) is found to be 17 times that of the {021} prismatic faces (1.9 × 10-9 m s-1). Both capping and prismatic faces are found to have dead zones in their growth kinetic profiles, with the capping faces (σc = 0.23) being about half that of the prismatic faces (σc = 0.46). The importance of this overall approach as an integral component of the digital design of industrial crystallization processes is highlighted.
RESUMEN
Grid-based systematic search methods are used to investigate molecule-molecule, molecule-surface, and surface-surface contributions to interparticle interactions in order to identify the crystal faces that most strongly affect particle behavior during powder blend formulation and delivery processes. The model system comprises terbutaline sulfate (TBS) as an active pharmaceutical ingredient (API) and α-form lactose monohydrate (LMH). A combination of systematic molecular modeling and X-ray computed tomography (XCT) is used to determine not only the adhesive and cohesive interparticle energies but, also the agglomeration behavior during manufacturing and de-agglomeration behavior during delivery after inhalation. This is achieved through a detailed examination of the balance between the adhesive and cohesive energies with the XCT results confirming the blend segregation tendencies, through the particle-particle de-agglomeration process. The results reveal that the cohesive interaction energies of TBS-TBS are higher than the adhesive energies between TBS and LMH, but that the cohesive energies of LMH-LMH are the smallest between molecule and molecule, molecule and surface, and surface and surface. This shows how systematic grid-search molecular modeling along with XCT can guide the digital formulation design of inhalation powders in order to achieve optimum aerosolization and efficacy for inhaled medicines. This will lead to faster pharmaceutical design with less variability, higher quality, and enhanced performance.
RESUMEN
The influence of the solution environment on the solution thermodynamics, crystallizability, and nucleation of tolfenamic acid (TFA) in five different solvents (isopropanol, ethanol, methanol, toluene, and acetonitrile) is examined using an integrated workflow encompassing both experimental studies and intermolecular modeling. The solubility of TFA in isopropanol is found to be the highest, consistent with the strongest solvent-solute interactions, and a concomitantly higher than ideal solubility. The crystallizability is found to be highly dependent on the solvent type with the overall order being isopropanol < ethanol < methanol < toluene < acetonitrile with the widest solution metastable zone width in isopropanol (24.49 to 47.41 °C) and the narrowest in acetonitrile (8.23 to 16.17 °C). Nucleation is found to occur via progressive mechanism in all the solvents studied. The calculated nucleation parameters reveal a considerably higher interfacial tension and larger critical nucleus radius in the isopropanol solutions, indicating the higher energy barrier hindering nucleation and hence lowering the nucleation rate. This is supported by diffusion coefficient measurements which are lowest in isopropanol, highlighting the lower molecular diffusion in the bulk of solution compared to the other solutions. The TFA concentration and critical supersaturation at the crystallization onset is found to be directly correlated with TFA/isopropanol solutions having the highest values of solubility and critical supersaturation. Intermolecular modeling of solute-solvent interactions supports the experimental observations of the solubility and crystallizability, highlighting the importance of understanding solvent selection and solution state structure at the molecular level in directing the solubility, solute mass transfer, crystallizability, and nucleation kinetics.
RESUMEN
Detailed atomistic interactions of 1,1,1,2-tetrafluoroethane (HFA-134a) liquid were presented in a data format, namely, DL_ANALYSER Notation for Atomic Interactions (DANAI), that annotates precisely the nature of interactions that is discoverable and searchable without having to resolve to diagrammatic illustrations. The datasets were obtained from raw atomic trajectory files of HFA-134a pure liquid models produced by using DL_POLY molecular dynamics software package. The trajectory datafiles contain expressions of atomic species in a natural chemical sense, and hence, provide localized key interactions, 'at a glance', of the liquid model on otherwise a typically disordered system consists of complex network of intermolecular interactions. The data provide insights to detailed structural behavior of molecules in liquid phase, and can be used as cheminformatics comparative investigations, linking to other molecular system models that contain similar interaction types and chemical species. This can form the foundation of investigations into the role of HFA-134a plays within different applications. For example, it can be used to compare structural and atomic interaction differences with alternative refrigerants, or as liquid propellants in pharmaceutical devices when solvating formulation ingredients.
RESUMEN
Fine chemicals produced via batch crystallization with properties dependent on the crystal size distribution require precise control of supersaturation, which drives the evolution of crystal size over time. Model predictive control (MPC) of supersaturation using a mechanistic model to represent the behavior of a crystallization process requires less experimental time and resources compared with fully empirical model-based control methods. Experimental characterization of the hexamine-ethanol crystallization system was performed in order to collect the parameters required to build a one-dimensional (1D) population balance model (PBM) in gPROMS FormulatedProducts software (Siemens-PSE Ltd.). Analysis of the metastable zone width (MSZW) and a series of seeded batch cooling crystallizations informed the suitable process conditions selected for supersaturation control experiments. The gPROMS model was integrated with the control software PharmaMV (Perceptive Engineering Ltd.) to create a digital twin of the crystallizer. Simulated batch crystallizations were used to train two statistical MPC blocks, allowing for in silico supersaturation control simulations to develop an effective control strategy. In the supersaturation set-point range of 0.012-0.036, the digital twin displayed excellent performance that would require minimal controller tuning to steady out any instabilities. The MPC strategy was implemented on a physical 500 mL crystallizer, with the simulated solution concentration replaced by in situ measurements from calibrated attenuated total reflection-Fourier transform infrared (ATR-FTIR) spectroscopy. Physical supersaturation control performance was slightly more unstable than the in silico tests, which is consistent with expected disturbances to the heat transfer, which were not specifically modeled in simulations. Overall, the level of supersaturation control in a real crystallizer was found to be accurate and precise enough to consider future adaptations to the MPC strategy for more advanced control objectives, such as the crystal size.
RESUMEN
The data presented in this article relates to the crystallisation of 8 single n-alkanes, C16H34 - C23H48 in representative diesel solvents dodecane and toluene, as well as a mixture of these 8-alkanes with a composition representative of real diesel fuel in the same solvents. For the single alkane systems, the data was collected over a range of 5 concentrations ranging from 0.09 - 0.311xi, depending upon the system, and 4 concentrations for the 8-alkane mixture, 0.1 - 0.5xi. Raw average crystallisation and dissolution points as a function of cooling rate (q) from a polythermal methodology are presented. Along with the equilibrium crystallisation and dissolution temperatures, van't Hoff fitting parameters, relative critical undercooling (uc) values as a function of q as well as the calculated values of KG and αdet.
RESUMEN
The pharmaceutical compound entacapone ((E)-2-cyano-3-(3,4-dihydroxy-5-nitrophenyl)-N,N-diethylprop-2-enamide) is important in the treatment of Parkinson's disease, exhibiting interesting polymorphic behavior upon crystallization from solution. It consistently produces its stable form A with a uniform crystal size distribution on the surface of an Au(111) template while concomitantly forming its metastable form D within the same bulk solution. Molecular modeling using empirical atomistic force-fields reveals more complex molecular and intermolecular structures for form D compared to form A, with the crystal chemistry of both polymorphs being dominated by van der Waals and π-π stacking interactions with lower contributions (ca. 20%) from hydrogen bonding and electrostatic interactions. Comparative lattice energies and convergence for the polymorphs are consistent with the observed concomitant polymorphic behavior. Synthon characterization reveals an elongated needle-like morphology for form D crystals in contrast to the more equant form A crystals with the surface chemistry of the latter exposing the molecules' cyano groups on its {010} and {011} habit faces. Density functional theory modeling of surface adsorption reveals preferential interactions between Au and the synthon GA interactions of form A on the Au surface. Molecular dynamics modeling of the entacapone/gold interface reveals the entacapone molecular structure within the first adsorbed layer to show nearly identical interaction distances, for both the molecules within form A or D with respect to the Au surface, while in the second and third layers when entacapone molecule-molecule interactions overtake the interactions between those of molecule-Au, the intermolecular structures are found to be closer to the form A structure than form D. In these layers, synthon GA (form A) could be reproduced with just two small azimuthal rotations (5° and 15°) whereas the closest alignment to a form D synthon requires larger azimuthal rotations (15° and 40°). The cyano functional group interactions with the Au template dominate interfacial interactions with these groups being aligned parallel to the Au surface and with nearest neighbor distances to Au atoms more closely matching those in form A than form D. The overall polymorph direction pathway thus encompasses consideration of molecular, crystal, and surface chemistry factors.
RESUMEN
Crystallization and dissolution are important processes to consider in drug development as well as many other industrial processes. Many current growth and dissolution models are based on bulk solution properties and do not implicitly consider concentration variation close to the crystal surface-solution interface and how this is mediated by solute diffusive mass transfer. Solution boundary layer thickness and concentration distribution, for the {120} crystal habit face of single crystals of l-alanine in saturated aqueous solutions during both growth and dissolution processes, is measured as a function of super/undersaturation using a Mach-Zehnder optical interferometer system. Further analysis allows determination of the diffusion coefficient and mass flux within the boundary layer as well as whether the processes are controlled by mass transfer or crystal interfacial kinetics. The measurement of this study revealed that the {120} face was not saturated at its surface during growth or dissolution meaning both processes were somewhat limited by their crystal interfacial kinetics. Growth was limited by crystal interfacial kinetics at all supersaturations to the same degree, whereas dissolution had a mixed dependency on crystal interfacial kinetics and mass transfer at lower undersaturations becoming more limited by mass transfer at higher undersaturations. Boundary layer thickness increased with super/undersaturation but to a lesser degree than the increase in the concentration difference between the crystal surface and bulk solution leading to a higher mass flux of solute molecules through the boundary layer. At the same relative super/undersaturation mass flux of solute molecules was faster during dissolution which was concurrent with its increased surface to bulk solution concentration difference and boundary layer thickness.
RESUMEN
Intermolecular (synthonic) modelling is used for a statistical analysis of crystal lattice energies, together with their contributing intermolecular interactions for the crystallographic structures selected from the CCDC's Drug Subset (https://doi.org/10.1016/j.xphs.2018.12.011). Analysis of this selected subset reveal similarities in packing compared to other organic crystals in the CSD with linear relationships between molecular weight and unit cell volume, void space, and packing coefficient. Crystal lattice energy calculations converge within a 30 Å intermolecular radius characterised by a mean lattice energy of ca. -36 kcal mol-1 with ca. 85% and 15% due to dispersive and electrostatic interactions, respectively. The distribution of the strongest synthons within the individual structures reveals an average strength of -5.79 kcal mol-1. The diversity of chemical space within the drug molecules is in agreement with the analysis of atom types across the selected subset with phenyl groups being found to contribute the highest mean energy of -11.28 kcal mol-1, highlighting the importance of aromatic interactions within pharmaceutical compounds. Despite an initial focus on Z'â¯=â¯1 structures, this automated approach enables rapid and consistent quantitative analysis of lattice energy, synthon strength and functional group contributions, providing solid-form informatics for pharmaceutical R&D and a helpful basis for further investigations.
Asunto(s)
Fenómenos Físicos , Cristalografía , Preparaciones FarmacéuticasRESUMEN
The data presented in this article relates to the crystallisation of the long chain hydrocarbon eicosane (C20H42), from supersaturated toluene solutions in the absence/presence of a commercially available cold-flow improver additive (IA) at different solution treat rates. Data was collected for treat rates of 0, 0.1, 0.5, 2, 3, 5 and 10 wt% IA with respect to eicosane, with each treat rate studied over four solution concentrations. Data is collected by transmission vs. solution temperature experimental investigations and is analysed through a conventional transmission analysis route (STR) and a reanalysed route that takes into account multiple phase transformation behaviour (RRT). Average crystallisation and dissolution data is provided over a range of solution concentrations and cooling rates used under a polythermal crystallisation methodology for each analysis route. Equilibrium saturation temperature, supersolubility and metastable zone width data is also presented for each treat rate, concentration and analysis route. Laser transmission as a function of solution temperature profiles are displayed for IA crystallising from toluene solutions. This data relates to the research article: Kaskiewicz, P. L., Downie, R., Dowding, P. J., George, N. & Roberts, K. J. Influence of a Polymeric Additive on the Crystallisability and Nucleation Mechanism for the Model Fuel System of Eicosane Crystallising from Supersaturated Toluene Solutions. J. Cryst. Growth 581, (2021) 126,470. https://doi.org/10.1016/j.jcrysgro.2021.126470.
RESUMEN
PURPOSE: Development of a quantitative model and associated workflow for predicting the mechanical deformation properties (plastic deformation or cleavage fracture) of organic single crystals from their crystallographic structures using molecular and crystallographic modelling. METHODS: Intermolecular synthons, hydrogen bonding, crystal morphology and surface chemistry are modelled using empirical force fields with the data integrated into the analysis of lattice deformation as computed using a statistical approach. RESULTS: The approach developed comprises three main components. Firstly, the identification of the likely direction of deformation based on lattice unit cell geometry; secondly, the identification of likely lattice planes for deformation through the calculation of the strength and stereochemistry of interplanar intermolecular interactions, surface plane rugosity and surface energy; thirdly, identification of potential crystal planes for cleavage fracture by assessing intermolecular bonding anisotropy. Pentaerythritol is predicted to fracture by brittle cleavage on the {001} lattice planes by strong in-plane hydrogen-bond interactions in the <110>, whereas pentaerythritol tetranitrate is predicted to deform by plastic deformation through the slip system {110} < 001>, with both predictions being in excellent agreement with known experimental data. CONCLUSION: A crystallographic framework and associated workflow for predicting the mechanical deformation of molecular crystals is developed through quantitative assessment of lattice energetics, crystal surface chemistry and crystal defects. The potential for the de novo prediction of the mechanical deformation of pharmaceutical materials using this approach is highlighted for its potential importance in the design of formulated drug products process as needed for manufacture by direct compression.
Asunto(s)
Tetranitrato de PentaeritritolRESUMEN
The solvent-mediated crystal morphologies of the α and ß polymorphic forms of l-glutamic acid are presented. This work applies a digital mechanistically based workflow that encompasses calculation of the crystal lattice energy and its constituent intermolecular synthons, their interaction energies, and their key role in understanding and predicting crystal morphology as well as assessing the surface chemistry, topology, and solvent binding on crystal habit growth surfaces. Through a comparison between the contrasting morphologies of the conformational polymorphs of l-glutamic acid, this approach highlights how the interfacial chemistry of organic crystalline materials and their inherent anisotropic interactions with their solvation environments direct their crystal habit with potential impact on their further downstream processing behavior.
RESUMEN
Para amino benzoic acid (PABA) has two well-characterised α- and ß-polymorphic forms and, whilst both crystallise in the monoclinic space group P21/n, they have quite different crystal chemistry and crystallisability behaviour. Previous work has shown that the molecular conformation deformation energy in the crystalline state is higher for the ß-form than for the α-form and that the lattice energy for the former converges more slowly than for the latter overall. This suggests that not only is there a higher barrier to crystallisation for the ß-form but also that low solution supersaturations might be needed for it to preferentially nucleate. Additionally, solute cluster propensity and solute solvation energetic analysis highlight the importance of an aqueous solvation environment in inhibiting the α-form's strong OHâ¯O carboxylic acid hydrogen bond (H-bond) dimer. Despite this, the detailed molecular-scale pathway from solvated molecules to 3D crystallographic structure still remains unclear, most notably regarding how the nucleation process is activated and how, as a result, this mediates the preferential formation of either of the two polymorphic forms. Molecular dynamics (MD) simulations coupled with FTIR studies and intermolecular synthon analysis address this issue through characterisation of the propensity of the incipient bulk synthons that are important in the crystallisation of the two polymorphic forms within the solution state. MD molecular trajectory analysis within crystallisation solutions reveals a greater propensity for OHâ¯O synthons (both single H-bonds and homodimers) typical of the α-form and NHâ¯O synthons found in both the α- and ß-forms when compared to aqueous solution but much lower propensities for the ß-form's "fingerprinting" OHâ¯N and π-π stacking synthons. In contrast, data from the aqueous solution environment reveals a much greater propensity for the ß-form's π-π interaction synthons. IR dilution studies in acetonitrile in the carbonyl region reveal the presence of two CîO vibrational stretching bands, whose relative intensities vary as a function of solution dilution. These were assigned to the solvated PABA monomer and a COOH dimer of PABA. Similar data in ethanol shows a main CîO stretching band with a shoulder peak suggesting a similar monomer vs. dimer speciation may exist in this solvent. The IR data is consistent with the organic solvent MD data, albeit the corresponding analysis for the aqueous solution was precluded due to the latter's strong OH vibrational mode which restricted validation in aqueous solutions.
Asunto(s)
Ácido 4-Aminobenzoico , Aminoácidos , Enlace de Hidrógeno , Conformación Molecular , Solventes/química , Agua/químicaRESUMEN
Anti-TNFα and anti-IL-23 antibodies are highly effective therapies for Crohn's disease or ulcerative colitis in a proportion of patients. V56B2 is a novel bispecific domain antibody in which a llama-derived IL-23p19-specific domain antibody, humanised and engineered for intestinal protease resistance, V900, was combined with a previously-described TNFα-specific domain antibody, V565. V56B2 contains a central protease-labile linker to create a single molecule for oral administration. Incubation of V56B2 with trypsin or human faecal supernatant resulted in a complete separation of the V565 and V900 monomers without loss of neutralising potency. Following oral administration of V900 and V565 in mice, high levels of each domain antibody were detected in the faeces, demonstrating stability in the intestinal milieu. In ex vivo cultures of colonic biopsies from IBD patients, treatment with V565 or V900 inhibited tissue phosphoprotein levels and with a combination of the two, inhibition was even greater. These results support further development of V56B2 as an oral therapy for IBD with improved safety and efficacy in a greater proportion of patients as well as greater convenience for patients compared with traditional monoclonal antibody therapies.
Asunto(s)
Anticuerpos Biespecíficos , Anticuerpos Monoclonales/administración & dosificación , Evaluación Preclínica de Medicamentos/métodos , Enfermedades Inflamatorias del Intestino/tratamiento farmacológico , Interleucina-23/inmunología , Factor de Necrosis Tumoral alfa/inmunología , Animales , Anticuerpos Biespecíficos/administración & dosificación , Anticuerpos Biespecíficos/farmacología , Anticuerpos Monoclonales/farmacología , Femenino , Humanos , Ratones , Ratones Endogámicos C57BLRESUMEN
PURPOSE: Application of multi-scale modelling workflows to characterise polymorphism in ritonavir with regard to its stability, bioavailability and processing. METHODS: Molecular conformation, polarizability and stability are examined using quantum mechanics (QM). Intermolecular synthons, hydrogen bonding, crystal morphology and surface chemistry are modelled using empirical force fields. RESULTS: The form I conformation is more stable and polarized with more efficient intermolecular packing, lower void space and higher density, however its shielded hydroxyl is only a hydrogen bond donor. In contrast, the hydroxyl in the more open but less stable and polarized form II conformation is both a donor and acceptor resulting in stronger hydrogen bonding and a more stable crystal structure but one that is less dense. Both forms have strong 1D networks of hydrogen bonds and the differences in packing energies are partially offset in form II by its conformational deformation energy difference with respect to form I. The lattice energies converge at shorter distances for form I, consistent with its preferential crystallization at high supersaturation. Both forms exhibit a needle/lath-like crystal habit with slower growing hydrophobic side and faster growing hydrophilic capping habit faces with aspect ratios increasing from polar-protic, polar-aprotic and non-polar solvents, respectively. Surface energies are higher for form II than form I and increase with solvent polarity. The higher deformation, lattice and surface energies of form II are consistent with its lower solubility and hence bioavailability. CONCLUSION: Inter-relationship between molecular, solid-state and surface structures of the polymorphic forms of ritonavir are quantified in relation to their physical-chemical properties.
Asunto(s)
Química Farmacéutica/métodos , Cristalización/métodos , Inhibidores de la Proteasa del VIH/química , Conformación Molecular , Ritonavir/química , Fenómenos Químicos , Inhibidores de la Proteasa del VIH/metabolismo , Enlace de Hidrógeno , Interacciones Hidrofóbicas e Hidrofílicas , Ritonavir/metabolismo , Solubilidad , Propiedades de SuperficieRESUMEN
V565 is an engineered TNFα-neutralising single domain antibody formulated into enteric coated mini-tablets to enable release in the intestine after oral administration as a possible oral treatment for inflammatory bowel disease (IBD). Following oral administration, ileal recovery of V565 was investigated in four patients with terminal ileostomy. Intestinal and systemic pharmacokinetics were measured in six patients with Crohn's disease and evidence of target engagement assessed in five patients with ulcerative colitis. Following oral administration, V565 was detected at micromolar concentrations in ileal fluid from the ileostomy patients and in stools of the Crohn's patients. In four of the five ulcerative colitis patients, biopsies taken after 7d dosing demonstrated V565 in the lamina propria with co-immunostaining on CD3+ T-lymphocytes and CD14+ macrophages. Phosphorylation of signalling proteins in biopsies taken after 7d oral dosing was decreased by approximately 50%. In conclusion, enteric coating of V565 mini-tablets provided protection in the stomach with gradual release in intestinal regions affected by IBD. Immunostaining revealed V565 tissue penetration and association with inflammatory cells, while decreased phosphoproteins after 7d oral dosing was consistent with V565-TNFα engagement and neutralising activity. Overall these results are encouraging for the clinical utility of V565 in the treatment of IBD.
Asunto(s)
Anticuerpos/uso terapéutico , Colitis Ulcerosa/tratamiento farmacológico , Inmunoterapia/métodos , Factor de Necrosis Tumoral alfa/antagonistas & inhibidores , Adulto , Anciano , Anticuerpos/análisis , Anticuerpos/metabolismo , Femenino , Humanos , Intestinos/química , Masculino , Persona de Mediana Edad , Proteínas Recombinantes/uso terapéutico , Factor de Necrosis Tumoral alfa/inmunologíaRESUMEN
Crystal lattice energy is a key property affecting the ease of processing pharmaceutical materials during manufacturing, as well as product performance. We present an extensive comparison of 324 force-field protocols for calculating the lattice energies of single component, organic molecular crystals (further restricted to Z' less than or equal to one), corresponding to a wide variety of force-fields (DREIDING, Universal, CVFF, PCFF, COMPASS, COMPASSII), optimization routines, and other variations, which could be implemented as part of an automated workflow using the industry standard Materials Studio software. All calculations were validated using a large new dataset (SUB-BIG), which we make publicly available. This dataset comprises public domain sublimation data, from which estimated experimental lattice energies were derived, linked to 235 molecular crystals. Analysis of pharmaceutical relevance was performed according to two distinct methods based upon (A) public and (B) proprietary data. These identified overlapping subsets of SUB-BIG comprising (A) 172 and (B) 63 crystals, of putative pharmaceutical relevance, respectively. We recommend a protocol based on the COMPASSII force field for lattice energy calculations of general organic or pharmaceutically relevant molecular crystals. This protocol was the most highly ranked prior to subsetting and was either the top ranking or amongst the top 15 protocols (top 5%) following subsetting of the dataset according to putative pharmaceutical relevance. Further analysis identified scenarios where the lattice energies calculated using the recommended force-field protocol should either be disregarded (values greater than or equal to zero and/or the messages generated by the automated workflow indicate extraneous atoms were added to the unit cell) or treated cautiously (values less than or equal to -249 kJ/mol), as they are likely to be inaccurate. Application of the recommended force-field protocol, coupled with these heuristic filtering criteria, achieved an root mean-squared error (RMSE) around 17 kJ/mol (mean absolute deviation (MAD) around 11 kJ/mol, Spearman's rank correlation coefficient of 0.88) across all 226 SUB-BIG structures retained after removing calculation failures and applying the filtering criteria. Across these 226 structures, the estimated experimental lattice energies ranged from -60 to -269 kJ/mol, with a standard deviation around 29 kJ/mol. The performance of the recommended protocol on pharmaceutically relevant crystals could be somewhat reduced, with an RMSE around 20 kJ/mol (MAD around 13 kJ/mol, Spearman's rank correlation coefficient of 0.76) obtained on 62 structures retained following filtering according to pharmaceutical relevance method B, for which the distribution of experimental values was similar. For a diverse set of 17 SUB-BIG entries, deemed pharmaceutically relevant according to method B, this recommended force-field protocol was compared to dispersion corrected density functional theory (DFT) calculations (PBE + TS). These calculations suggest that the recommended force-field protocol (RMSE around 15 kJ/mol) outperforms PBE + TS (RMSE around 37 kJ/mol), although it may not outperform more sophisticated DFT protocols and future studies should investigate this. Finally, further work is required to compare our recommended protocol to other lattice energy calculation protocols reported in the literature, as comparisons based upon previously reported smaller datasets indicated this protocol was outperformed by a number of other methods. The SUB-BIG dataset provides a basis for these future studies and could support protocol refinement.
Asunto(s)
Compuestos Orgánicos/química , Preparaciones Farmacéuticas/química , Termodinámica , Algoritmos , Cristalización , Bases de Datos Farmacéuticas , Teoría Funcional de la Densidad , Modelos Químicos , Modelos Moleculares , Programas InformáticosRESUMEN
Organic molecular crystals contain long-range dispersion interactions that can be challenging for solid-state methods such as density functional theory (DFT) to capture, and in some industrial sectors are overlooked in favor of classical methods to calculate atomistic properties. Hence, this publication addresses the critical question of whether dispersion corrected DFT calculations for organic crystals can reproduce the structural and energetic trends seen from experiment, i.e., whether the calculations can now be said to be truly "on-trend." In this work, we assess the performance of three of the latest dispersion-corrected DFT methods, in calculating the long-range, dispersion energy: the pairwise methods of D3(0) and D3(BJ) and the many-body dispersion method, MBD@rsSCS. We calculate the energetics and optimized structures of two homologous series of organic molecular crystals, namely, carboxylic acids and amino acids. We also use a classical force field method (using COMPASS II) and compare all results to experimental data where possible. The mean absolute error in lattice energies is 9.59 and 343.85 kJ/mol (COMPASS II), 10.17 and 16.23 kJ/mol (MBD@rsSCS), 10.57 and 18.76 kJ/mol [D3(0)], and 8.52 and 14.66 kJ/mol [D3(BJ)] for the carboxylic acids and amino acids, respectively. MBD@rsSCS produces structural and energetic trends that most closely match experimental trends, performing the most consistently across the two series and competing favorably with COMPASS II.