Surface Analysis-From Crystal Structures to Particle Properties.

Understanding the surface properties of particles is crucial for optimizing the performance of formulated products in various industries. However, acquiring this understanding often requires expensive trial-and-error studies. Here, we present advanced surface analysis tools that enable the visualization and quantification of chemical and topological information derived from crystallographic data. By employing functional group analysis, roughness calculations, and statistical interaction data, we facilitate direct comparisons of surfaces. We further demonstrate the practicality of our approach by correlating the sticking propensity of distinct ibuprofen morphologies with surface and particle descriptors calculated from a single crystal structure. Our findings support and expand upon previous work, demonstrating that the presence of a carboxylic acid group on the {011} facet leads to significant differences in particle properties and explains the higher electrostatic potential observed in the block-like morphology. While our surface analysis tools are not intended to replace the importance of chemical intuition and expertise, they provide valuable insights for formulators and particle engineers, facilitating informed, data-driven decisions to mitigate formulation risks. This research represents a significant step toward a comprehensive understanding of particle surfaces and their impact on products.

Exploring the CSD Drug Subset: An Analysis of Lattice Energies and Constituent Intermolecular Interactions for the Crystal Structures of Pharmaceuticals.

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.

Mechanistic model for production of recombinant adeno-associated virus via triple transfection of HEK293 cells.

Manufacturing of recombinant adeno-associated virus (rAAV) viral vectors remains challenging, with low yields and low full:empty capsid ratios in the harvest. To elucidate the dynamics of recombinant viral production, we develop a mechanistic model for the synthesis of rAAV viral vectors by triple plasmid transfection based on the underlying biological processes derived from wild-type AAV. The model covers major steps starting from exogenous DNA delivery to the reaction cascade that forms viral proteins and DNA, which subsequently result in filled capsids, and the complex functions of the Rep protein as a regulator of the packaging plasmid gene expression and a catalyst for viral DNA packaging. We estimate kinetic parameters using dynamic data from literature and in-house triple transient transfection experiments. Model predictions of productivity changes as a result of the varied input plasmid ratio are benchmarked against transfection data from the literature. Sensitivity analysis suggests that (1) the poorly coordinated timeline of capsid synthesis and viral DNA replication results in a low ratio of full virions in harvest, and (2) repressive function of the Rep protein could be impeding capsid production at a later phase. The analyses from the mathematical model provide testable hypotheses for evaluation and reveal potential process bottlenecks that can be investigated.

Cellular pathways of recombinant adeno-associated virus production for gene therapy.

Recombinant adeno-associated viruses (rAAVs) are among the most important vectors for in vivo gene therapies. With the rapid development of gene therapy, current rAAV manufacturing capacity faces a challenge to meet the emerging demand for these therapies in the future. To examine the bottlenecks in rAAV production during cell culture, we focus here on an analysis of cellular pathways of rAAV production, based on an overview of assembly mechanisms first in the wild-type (wt) AAV replication and then in the common methods of rAAV production. The differences analyzed between the wild-type and recombinant systems provide insights into the mechanistic differences that may correlate with viral productivity. Based on these analyses, we identify potential barriers to high productivity of rAAV and discuss future directions for improvement to meet the emerging needs set by the growth of rAAV-based therapy and the needs of patients.

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.

Analytical methods for process and product characterization of recombinant adeno-associated virus-based gene therapies.

The optimization of upstream and downstream processes for production of recombinant adeno-associated virus (rAAV) with consistent quality depends on the ability to rapidly characterize critical quality attributes (CQAs). In the context of rAAV production, the virus titer, capsid content, and aggregation are identified as potential CQAs, affecting the potency, purity, and safety of rAAV-mediated gene therapy products. Analytical methods to measure these attributes commonly suffer from long turnaround times or low throughput for process development, although rapid, high-throughput methods are beginning to be developed and commercialized. These methods are not yet well established in academic or industrial practice, and supportive data are scarce. Here, we review both established and upcoming analytical methods for the quantification of rAAV quality attributes. In assessing each method, we highlight the progress toward rapid, at-line characterization of rAAV. Furthermore, we identify that a key challenge for transitioning from traditional to newer methods is the scarcity of academic and industrial experience with the latter. This literature review serves as a guide for the selection of analytical methods targeting quality attributes for rapid, high-throughput process characterization during process development of rAAV-mediated gene therapies.

Modeling of copy number variability in Pichia pastoris.

Development of continuous biopharmaceutical manufacturing processes is an area of active research. This study considers the long-term transgene copy number stability of Pichia pastoris in continuous bioreactors. We propose a model of copy number loss that quantifies population heterogeneity. An analytical solution is derived and compared with existing experimental data. The model is then used to provide guidance for stable operating timescales. The model is extended to consider copy number dependent growth such as in the case of Zeocin supplementation. The model is also extended to analyze a continuous seeding strategy. This study is a critical step towards understanding the impact of continuous processing on the stability of Pichia pastoris and the resultant products.

Macroscopic modeling of bioreactors for recombinant protein producing Pichia pastoris in defined medium.

The methylotrophic yeast Pichia pastoris is widely used as a microbial host for recombinant protein production. Bioreactor models for P. pastoris can inform understanding of cellular metabolism and can be used to optimize bioreactor operation. This article constructs an extensive macroscopic bioreactor model for P. pastoris which describes substrates, biomass, total protein, other medium components, and off-gas components. Species and elemental balances are introduced to describe uptake and evolution rates for medium components and off-gas components. Additionally, a pH model is constructed using an overall charge balance, acid/base equilibria, and activity coefficients to describe production of recombinant protein and precipitation of medium components. The extent of run-to-run variability is modeled by distributions of a subset of the model parameters, which are estimated using the maximum likelihood method. Model prediction from the extensive macroscopic bioreactor model well describes experimental data with different operating conditions. The probability distributions of the model predictions quantified from the parameter distribution are quantifiably consistent with the run-to-run variability observed in the experimental data. The uncertainty description in this macroscopic bioreactor model identifies the model parameters that have large variability and provides guidance as to which aspects of cellular metabolism should be the focus of additional experimental studies. The model for medium components with pH and precipitation can be used for improving chemically defined medium by minimizing the amount of components needed while meeting cellular requirements.

Quality of life outcomes in tracheobronchomalacia surgery.

Tracheobronchomalacia (TBM) is an obstructive airway disease characterized by laxity and redundancy of the posterior membrane of the main airways leading to dynamic airway collapse during exhalation. The gold standard for diagnosis is dynamic computed tomography (DCT) scan and dynamic flexible bronchoscopy (DFB). Patients with complete or near-complete collapse (>90% reduction in cross-sectional area) of the airway are possible candidates for surgical management. Central airway stabilization by tracheobronchoplasty (TBP) effectively corrects malacic airways and has demonstrated significant improvement in objective functional measures, which is often but not uniformly accompanied by equal improvement in health-related quality of life (HRQOL) metrics. This article reviews HRQOL instruments used to report outcomes after TBM surgery.

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.

The CSD Drug Subset: The Changing Chemistry and Crystallography of Small Molecule Pharmaceuticals.

We report the generation and statistical analysis of the CSD drug subset: a subset of the Cambridge Structural Database (CSD) consisting of every published small-molecule crystal structure containing an approved drug molecule. By making use of InChI matching, a CSD Python API workflow to link CSD entries to the online database Drugbank.ca has been produced. This has resulted in a subset of 8632 crystal structures, representing all published solid forms of 785 unique drug molecules. We hope that this new resource will lead to improvements in targeted cheminformatics and statistical model building in a pharmaceutical setting. In addition to this, as part of the Advanced Digital Design of Pharmaceutical Therapeutics collaboration between academia and industry, we have been given the unique opportunity to run comparative analysis on the internal crystal structure databases of AstraZeneca and Pfizer, alongside comparison to the CSD as a whole.

A chiral diamine: practical implications of a three-stereoisomer cocrystallization.

A brief comparison of seven straightforward methods for molecular crystal-volume estimation revealed that their precisions are comparable. A chiral diamine, N2,N3-bis[2,6-bis(propan-2-yl)phenyl]butane-2,3-diamine, C28H44N2, has been used to illustrate the application of the methods. Three stereoisomers of the diamine cocrystallize in the centrosymmetric space group P21/c with Z' = 1.5. The molecules occupying general positions are RR and SS, whereas that residing on an inversion center is meso. This is one of only ten examples of three stereoisomers with two asymmetric atoms cocrystallizing together reported to the Cambridge Structural Database (CSD). The conformations of the SS/RR and meso molecules differ considerably and lead to statistically significantly different C(asymmetric)-C(asymmetric) bond lengths in the diastereomers. An advanced Python script-based CSD searching technique for chiral compounds is presented.

Undersized Stent Grafts for Acute Mesenteric Ischemia in Chronic Type B Dissection.

Acute ischemia in chronic type B dissections carries high rates of morbidity and mortality. A 29-year-old woman with a chronic type B dissection presented with acute abdominal pain. Imaging revealed a worsening dissection with pseudocoarctation causing near complete occlusion of the true lumen by the false lumen. We placed purposefully undersized stent grafts to treat acute mesenteric ischemia by improving true lumen flow. The patient was discharged on postoperative day 4 without adverse events. We suggest that endovascular rescue by placing undersized stent grafts can provide improved flow to the mesenteric vessels with continued false lumen flow to vital organs.

Use of the PIXEL method to investigate gas adsorption in metal-organic frameworks.

PIXEL has been used to perform calculations of adsorbate-adsorbent interaction energies between a range of metal-organic frameworks (MOFs) and simple guest molecules. Interactions have been calculated for adsorption between MOF-5 and Ar, H2, and N2; Zn2(BDC)2(TED) (BDC = 1,4-benzenedicarboxylic acid, TED = triethylenediamine) and H2; and HKUST-1 and CO2. The locations of the adsorption sites and the calculated energies, which show differences in the Coulombic or dispersion characteristic of the interaction, compare favourably to experimental data and literature energy values calculated using density functional theory.

Cardiac arrhythmias induced by glutathione oxidation can be inhibited by preventing mitochondrial depolarization.

We have previously proposed that the heterogeneous collapse of mitochondrial inner membrane potential (DeltaPsi(m)) during ischemia and reperfusion contributes to arrhythmogenesis through the formation of metabolic sinks in the myocardium, wherein clusters of myocytes with uncoupled mitochondria and high K(ATP) current levels alter electrical propagation to promote reentry. Single myocyte studies have also shown that cell-wide DeltaPsi(m) depolarization, through a reactive oxygen species (ROS)-induced ROS release mechanism, can be triggered by global depletion of the antioxidant pool with diamide, a glutathione oxidant. Here we examine whether diamide causes mitochondrial depolarization and promotes arrhythmias in normoxic isolated perfused guinea pig hearts. We also investigate whether stabilization of DeltaPsi(m) with a ligand of the mitochondrial benzodiazepine receptor (4'-chlorodiazepam; 4-ClDzp) prevents the formation of metabolic sinks and, consequently, precludes arrhythmias. Oxidation of the GSH pool was initiated by treatment with 200 microM diamide for 35 min, followed by washout. This treatment increased GSSG and decreased both total GSH and the GSH/GSSG ratio. All hearts receiving diamide transitioned from sinus rhythm into ventricular tachycardia and/or ventricular fibrillation during the diamide exposure: arrhythmia scores were 5.5+/-0.5; n=6 hearts. These arrhythmias and impaired LV function were significantly inhibited by co-administration of 4-ClDzp (64 microM): arrhythmia scores with diamide+4-ClDzp were 0.4+/-0.2 (n=5; P<0.05 vs. diamide alone). Imaging DeltaPsi(m) in intact hearts revealed the heterogeneous collapse of DeltaPsi(m) beginning 20 min into diamide, paralleling the timeframe for the onset of arrhythmias. Loss of DeltaPsi(m) was prevented by 4-ClDzp treatment, as was the increase in myocardial GSSG. These findings show that oxidative stress induced by oxidation of GSH with diamide can cause electromechanical dysfunction under normoxic conditions. Analogous to ischemia-reperfusion injury, the dysfunction depends on the mitochondrial energy state. Targeting the mitochondrial benzodiazepine receptor can prevent electrical and mechanical dysfunction in both models of oxidative stress.