Methodological artefacts cause counter-intuitive evolutionary conclusions in a simulation study.

In their simulation study, Garcia-Costoya et al. (2023) conclude that evolutionary constraints might aid populations facing climate change. However, we are concerned that this conclusion is largely a consequence of the simulated temperature variation being too small, and, most importantly, that uneven limitations to standing variation disadvantage unconstrained populations.

A time course analysis through diapause reveals dynamic temporal patterns of microRNAs associated with endocrine regulation in the butterfly Pieris napi.

Organisms inhabiting highly seasonal environments must cope with a wide range of environmentally induced challenges. Many seasonal challenges require extensive physiological modification to survive. In winter, to survive extreme cold and limited resources, insects commonly enter diapause, which is an endogenously derived dormant state associated with minimized cellular processes and low energetic expenditure. Due to the high degree of complexity involved in diapause, substantial cellular regulation is required, of which our understanding primarily derives from the transcriptome via messenger RNA expression dynamics. Here we aim to advance our understanding of diapause by investigating microRNA (miRNA) expression in diapausing and direct developing pupae of the butterfly Pieris napi. We identified coordinated patterns of miRNA expression throughout diapause in both head and abdomen tissues of pupae, and via miRNA target identification, found several expression patterns to be enriched for relevant diapause-related physiological processes. We also identified two candidate miRNAs, miR-14-5p and miR-2a-3p, that are likely involved in diapause progression through their activity in the ecdysone pathway, a critical regulator of diapause termination. miR-14-5p targets phantom, a gene in the ecdysone synthesis pathway, and is upregulated early in diapause. miR-2a-3p has been found to be expressed in response to ecdysone, and is upregulated during diapause termination. Together, the expression patterns of these two miRNAs match our current understanding of the timing of hormonal regulation of diapause in P. napi and provide interesting candidates to further explore the mechanistic role of microRNAs in diapause regulation.

Influence of Solvent Selection on the Crystallizability and Polymorphic Selectivity Associated with the Formation of the "Disappeared" Form I Polymorph of Ritonavir.

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.

Genetic constraints in genes exhibiting splicing plasticity in facultative diapause.

Phenotypic plasticity is produced and maintained by processes regulating the transcriptome. While differential gene expression is among the most important of these processes, relatively little is known about other sources of transcriptional variation. Previous work suggests that alternative splicing plays an extensive and functionally unique role in transcriptional plasticity, though plastically spliced genes may be more constrained than the remainder of expressed genes. In this study, we explore the relationship between expression and splicing plasticity, along with the genetic diversity in those genes, in an ecologically consequential polyphenism: facultative diapause. Using 96 samples spread over two tissues and 10 timepoints, we compare the extent of differential splicing and expression between diapausing and direct developing pupae of the butterfly Pieris napi. Splicing differs strongly between diapausing and direct developing trajectories but alters a smaller and functionally unique set of genes compared to differential expression. We further test the hypothesis that among these expressed loci, plastically spliced genes are likely to experience the strongest purifying selection to maintain seasonally plastic phenotypes. Genes with unique transcriptional changes through diapause consistently had the lowest nucleotide diversity, and this effect was consistently stronger among genes that were differentially spliced compared to those with just differential expression through diapause. Further, the strength of negative selection was higher in the population expressing diapause every generation. Our results suggest that maintenance of the molecular mechanisms involved in diapause progression, including post-transcriptional modifications, are highly conserved and likely to experience genetic constraints, especially in northern populations of P. napi.

Bumblebee thermoregulation at increasing temperatures is affected by behavioral state.

Over the past decades, increasing environmental temperatures have been identified as one of the causes of major insect population declines and biodiversity loss. However, it is unclear how these rising temperatures affect endoheterothermic insects, like bumblebees, that have evolved thermoregulatory capacities to exploit cold and temperate habitats. To investigate this, we measured head, thoracic, and abdominal temperature of bumblebee (Bombus terrestris) workers across a range of temperatures (24 °C-32 °C) during three distinct behaviors. In resting bumblebees, the head, abdomen, and thorax conformed to the environmental temperature. In pre-flight bumblebees, the head and abdominal temperatures were elevated with respect to the environmental temperature, while the thoracic temperature was maintained, indicating a pre-flight muscle warming stage. In post-flight bumblebees, abdominal temperature increased at the same rate as environmental temperature, but the head and the thoracic temperature did not. By calculating the excess temperature ratio, we show that thermoregulation in bumblebees during flight is partially achieved by the active transfer of heat produced in the thorax to the abdomen, where it can more easily be dissipated. These results provide the first indication that the thermoregulatory abilities of bumblebees are plastic and behavior dependent. We also show that the flight speed and number of workers foraging increase with increasing temperature, suggesting that bees do not avoid flying at these temperatures despite its impact on behavioral performance.

Genome-wide DNA methylation analysis implicates enrichment of interferon pathway in African American patients with Systemic Lupus Erythematosus and European Americans with lupus nephritis.

Systemic Lupus Erythematosus (SLE) is a chronic, multisystem, inflammatory autoimmune disease that disproportionately affects women. Trends in SLE prevalence and clinical course differ by ancestry, with those of African American ancestry presenting with more active, severe and rapidly progressive disease than European Americans. Previous research established altered epigenetic signatures in SLE patients compared to controls. However, the contribution of aberrant DNA methylation (DNAm) to the risk of SLE by ancestry and differences among patients with SLE-associated Lupus Nephritis (LN) has not been well described. We evaluated the DNA methylomes of 87 individuals including 41 SLE patients, with and without LN, and 46 controls enrolled in an ancestry diverse, well-characterized cohort study of established SLE (41 SLE patients [20 SLE-LN+, 21 SLE-LN-] and 46 sex-, race- and age-matched controls; 55% African American, 45% European American). Participants were genotyped using the Infinium Global Diversity Array (GDA), and genetic ancestry was estimated using principal components. Genome-wide DNA methylation was initially measured using the Illumina MethylationEPIC 850K Beadchip array followed by methylation-specific qPCR to validate the methylation status at putative loci. Differentially Methylated Positions (DMP) were identified using a case-control approach adjusted for ancestry. We identified a total of 51 DMPs in CpGs among SLE patients compared to controls. Genes proximal to these CpGs were highly enriched for involvement in type I interferon signaling. DMPs among European American SLE patients with LN were similar to African American SLE patients with and without LN. Our findings were validated using an orthogonal, methyl-specific PCR for three SLE-associated DMPs near or proximal to MX1, USP18, and IFITM1. Our study confirms previous reports that DMPs in CpGs associated with SLE are enriched in type I interferon genes. However, we show that European American SLE patients with LN have similar DNAm patterns to African American SLE patients irrespective of LN, suggesting that aberrant DNAm alters activity of type I interferon pathway leading to more severe disease independent of ancestry.

Predicting the Strength of Cohesive and Adhesive Interparticle Interactions for Dry Powder Inhalation Blends of Terbutaline Sulfate with α-Lactose Monohydrate.

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.

Seasonal energetics: are insects constrained by energy during dormancy?

In seasonal environments, many animals, including insects, enter dormancy, where they are limited to a fixed energy budget. The inability to replenish energetic stores during these periods suggests insects should be constrained by pre-dormancy energy stores. Over the last century, the community of researchers working on survival during dormancy has operated under the strong assumption that energy limitation is a key fitness trait driving the evolution of seasonal strategies. That is, energy use has to be minimized during dormancy because insects otherwise run out of energy and die during dormancy, or are left with too little energy to complete development, reproductive maturation or other costly post-dormancy processes such as dispersal or nest building. But if energy is so strongly constrained during dormancy, how can some insects - even within the same species and population - be dormant in very warm environments or show prolonged dormancy for many successive years? In this Commentary, we discuss major assumptions regarding dormancy energetics and outline cases where insects appear to align with our assumptions and where they do not. We then highlight several research directions that could help link organismal energy use with landscape-level changes. Overall, the optimal energetic strategy during dormancy might not be to simply minimize metabolic rate, but instead to maintain a level that matches the demands of the specific life-history strategy. Given the influence of temperature on energy use rates of insects in winter, understanding dormancy energetic strategies is critical in order to determine the potential impacts of climate change on insects in seasonal environments.

Landscape heterogeneity provides co-benefits to predator and prey.

Predator populations are imperiled globally, due in part to changing habitat and trophic interactions. Theoretical and laboratory studies suggest that heterogeneous landscapes containing prey refuges acting as source habitats can benefit both predator and prey populations, although the importance of heterogeneity in natural systems is uncertain. Here, we tested the hypothesis that landscape heterogeneity mediates predator-prey interactions between the California spotted owl (Strix occidentalis occidentalis)-a mature forest species-and one of its principal prey, the dusky-footed woodrat (Neotoma fuscipes)-a younger forest species-to the benefit of both. We did so by combining estimates of woodrat density and survival from live trapping and very high frequency tracking with direct observations of prey deliveries to dependent young by owls in both heterogeneous and homogeneous home ranges. Woodrat abundance was ~2.5 times higher in owl home ranges (14.12 km2 ) featuring greater heterogeneity in vegetation types (1805.0 ± 50.2 SE) compared to those dominated by mature forest (727.3 ± 51.9 SE), in large part because of high densities in young forests appearing to act as sources promoting woodrat densities in nearby mature forests. Woodrat mortality rates were low across vegetation types and did not differ between heterogeneous and homogeneous home ranges, yet all observed predation by owls occurred within mature forests, suggesting young forests may act as woodrat refuges. Owls exhibited a type 1 functional response, consuming ~2.5 times more woodrats in heterogeneous (31.1/month ± 5.2 SE) versus homogeneous (12.7/month ± 3.7 SE) home ranges. While consumption of smaller-bodied alternative prey partially compensated for lower woodrat consumption in homogeneous home ranges, owls nevertheless consumed 30% more biomass in heterogeneous home ranges-approximately equivalent to the energetic needs of producing one additional offspring. Thus, a mosaic of vegetation types including young forest patches increased woodrat abundance and availability that, in turn, provided energetic and potentially reproductive benefits to mature forest-associated spotted owls. More broadly, our findings provide strong empirical evidence that heterogeneous landscapes containing prey refuges can benefit both predator and prey populations. As anthropogenic activities continue to homogenize landscapes globally, promoting heterogeneous systems with prey refuges may benefit imperiled predators.

Transcriptome profiling of a synergistic volumetric muscle loss repair strategy.

Volumetric muscle loss overwhelms skeletal muscle's ordinarily capable regenerative machinery, resulting in severe functional deficits that have defied clinical repair strategies. In this manuscript we pair the early in vivo functional response induced by differing volumetric muscle loss tissue engineering repair strategies that are broadly representative of those explored by the field (scaffold alone, cells alone, or scaffold + cells) to the transcriptomic response induced by each intervention. We demonstrate that an implant strategy comprising allogeneic decellularized skeletal muscle scaffolds seeded with autologous minced muscle cellular paste (scaffold + cells) mediates a pattern of increased expression for several genes known to play roles in axon guidance and peripheral neuroregeneration, as well as several other key genes related to inflammation, phagocytosis, and extracellular matrix regulation. The upregulation of several key genes in the presence of both implant components suggests a unique synergy between scaffolding and cells in the early period following intervention that is not seen when either scaffolds or cells are used in isolation; a finding that invites further exploration of the interactions that could have a positive impact on the treatment of volumetric muscle loss.

Genetic associations with age at dementia onset in the PSEN1 E280A Colombian kindred.

INTRODUCTION: Genetic associations with Alzheimer's disease (AD) age at onset (AAO) could reveal genetic variants with therapeutic applications. We present a large Colombian kindred with autosomal dominant AD (ADAD) as a unique opportunity to discover AAO genetic associations. METHODS: A genetic association study was conducted to examine ADAD AAO in 340 individuals with the PSEN1 E280A mutation via TOPMed array imputation. Replication was assessed in two ADAD cohorts, one sporadic early-onset AD study and four late-onset AD studies. RESULTS: 13 variants had p<1×10-7 or p<1×10-5 with replication including three independent loci with candidate associations with clusterin including near CLU. Other suggestive associations were identified in or near HS3ST1, HSPG2, ACE, LRP1B, TSPAN10, and TSPAN14. DISCUSSION: Variants with suggestive associations with AAO were associated with biological processes including clusterin, heparin sulfate, and amyloid processing. The detection of these effects in the presence of a strong mutation for ADAD reinforces their potentially impactful role.

The impact of metabolic plasticity on winter energy use models.

Understanding the energetic consequences of climate change is critical to identifying organismal vulnerabilities, particularly for dormant organisms relying on finite energy budgets. Ecophysiological energy use models predict long-term energy use from metabolic rate, but we do not know the degree to which plasticity in metabolism impacts estimates. We quantified metabolic rate-temperature relationships of dormant willow leaf beetles (Chrysomela aeneicollis) monthly from February to May under constant and variable acclimation treatments. Metabolic rate increased as diapause progressed, and acclimation to variable conditions altered both metabolic intensity and thermal sensitivity. However, incorporating these two types of metabolic plasticity into energy use models did not improve energy use estimates, validated by empirical measurements of energy stores. While metabolic rate-temperature relationships are plastic during winter, the magnitude of inter-individual variability in energy stores overshadows the effects of incorporating plasticity into energy use models, highlighting the importance of within-population variation in energy reserves.

The structural pathway from its solvated molecular state to the solution crystallisation of the α- and ß-polymorphic forms of para amino benzoic acid.

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.

Prediction of the Mechanical Deformation Properties of Organic Crystals Based upon their Crystallographic Structures: Case Studies of Pentaerythritol and Pentaerythritol Tetranitrate.

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.

Delayed vaginal SHIV infection in VRC01 and anti-α4ß7 treated rhesus macaques.

VRC01 protects macaques from vaginal SHIV infection after a single high-dose challenge. Infusion of a simianized anti-α4ß7 mAb (Rh-α4ß7) just prior to, and during repeated vaginal exposures to SIVmac251 partially protected macaques from vaginal SIV infection and rescued CD4+ T cells. To investigate the impact of combining VRC01 and Rh-α4ß7 on SHIV infection, 3 groups of macaques were treated with a suboptimal dosing of VRC01 alone or in combination with Rh-α4ß7 or with control antibodies prior to the initiation of weekly vaginal exposures to a high dose (1000 TCID50) of SHIVAD8-EO. The combination Rh-α4ß7-VRC01 significantly delayed SHIVAD8-EO vaginal infection. Following infection, VRC01-Rh-α4ß7-treated macaques maintained higher CD4+ T cell counts and exhibited lower rectal SIV-DNA loads compared to controls. Interestingly, VRC01-Rh-α4ß7-treated macaques had fewer IL-17-producing cells in the blood and the gut during the acute phase of infection. Moreover, higher T cell responses to the V2-loop of the SHIVAD8-EO envelope in the VRC01-Rh-α4ß7 group inversely correlated with set point viremia. The combination of suboptimal amounts of VRC01 and Rh-α4ß7 delayed infection, altered antiviral immune responses and minimized CD4+ T cell loss. Further exploration of the effect of combining bNAbs with Rh-α4ß7 on SIV/HIV infection and antiviral immune responses is warranted and may lead to novel preventive and therapeutic strategies.

Snow modulates winter energy use and cold exposure across an elevation gradient in a montane ectotherm.

Snow insulates the soil from air temperature, decreasing winter cold stress and altering energy use for organisms that overwinter in the soil. As climate change alters snowpack and air temperatures, it is critical to account for the role of snow in modulating vulnerability to winter climate change. Along elevational gradients in snowy mountains, snow cover increases but air temperature decreases, and it is unknown how these opposing gradients impact performance and fitness of organisms overwintering in the soil. We developed experimentally validated ecophysiological models of cold and energy stress over the past decade for the montane leaf beetle Chrysomela aeneicollis, along five replicated elevational transects in the Sierra Nevada mountains in California. Cold stress peaks at mid-elevations, while high elevations are buffered by persistent snow cover, even in dry years. While protective against cold, snow increases energy stress for overwintering beetles, particularly at low elevations, potentially leading to mortality or energetic tradeoffs. Declining snowpack will predominantly impact mid-elevation populations by increasing cold exposure, while high elevation habitats may provide refugia as drier winters become more common.

Molecular, Solid-State and Surface Structures of the Conformational Polymorphic Forms of Ritonavir in Relation to their Physicochemical Properties.

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.

Additive manufacturing of platinum group element (PGE) reference materials with a silica matrix.

RATIONALE: The microanalytical community has an outstanding need for platinum group element (PGE) reference materials, particularly for trace element analysis by laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS). National Institute of Standards and Technology (NIST) glasses contain Rh, Pd, and Pt, but lack Ru, Os, and Ir. Synthesis of silicate PGE standards has proven difficult due the tendency of PGEs to form metallic nuggets. METHODS: Additive manufacturing methods were used to produce PGE standards with a silica matrix. Monodispersed submicron PGE-doped Stöber particles were used as feedstock materials for electrophoretic deposition (EPD). Two-cm-sized samples produced by EPD were subsequently densified by thermal processing. The homogeneity of PGEs was tested using LA-ICPMS and concentrations were measured by laser ablation and solution ICPMS. RESULTS: The PGE concentrations ranged from 0.5 to 3 µg/g. The inhomogeneity was at the 3% RSD level for Ru, Rh, Ir, and Os throughout and 5% for Pt and Pd in the interior of the samples. Based on LA-ICPMS analyses, the interiors of the two samples have near identical concentrations in PGEs. CONCLUSIONS: The samples fabricated in this study represent the most complete and homogeneous PGE standards produced with a silicate matrix. The ability to produce multiple samples with the same composition provides opportunities for validating methods, monitoring long-term reproducibility, and facilitating interlaboratory comparisons.

Evaluation of Force-Field Calculations of Lattice Energies on a Large Public Dataset, Assessment of Pharmaceutical Relevance, and Comparison to Density Functional Theory.

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.

Off-the-shelf DFT-DISPersion methods: Are they now "on-trend" for organic molecular crystals?

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.