1
|
Antipas GSE, Reul R, Voges K, Kyeremateng SO, Ntallis NA, Karalis KT, Miroslaw L. System-agnostic prediction of pharmaceutical excipient miscibility via computing-as-a-service and experimental validation. Sci Rep 2024; 14:15106. [PMID: 38956156 PMCID: PMC11219749 DOI: 10.1038/s41598-024-65978-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2024] [Accepted: 06/26/2024] [Indexed: 07/04/2024] Open
Abstract
We applied computing-as-a-service to the unattended system-agnostic miscibility prediction of the pharmaceutical surfactants, Vitamin E TPGS and Tween 80, with Copovidone VA64 polymer at temperature relevant for the pharmaceutical hot melt extrusion process. The computations were performed in lieu of running exhaustive hot melt extrusion experiments to identify surfactant-polymer miscibility limits. The computing scheme involved a massively parallelized architecture for molecular dynamics and free energy perturbation from which binodal, spinodal, and mechanical mixture critical points were detected on molar Gibbs free energy profiles at 180 °C. We established tight agreement between the computed stability (miscibility) limits of 9.0 and 10.0 wt% vs. the experimental 7 and 9 wt% for the Vitamin E TPGS and Tween 80 systems, respectively, and identified different destabilizing mechanisms applicable to each system. This paradigm supports that computational stability prediction may serve as a physically meaningful, resource-efficient, and operationally sensible digital twin to experimental screening tests of pharmaceutical systems. This approach is also relevant to amorphous solid dispersion drug delivery systems, as it can identify critical stability points of active pharmaceutical ingredient/excipient mixtures.
Collapse
Affiliation(s)
| | - Regina Reul
- AbbVie Deutschland GmbH & Co. KG, Development Sciences, 67061, Ludwigshafen, Germany
| | - Kristin Voges
- AbbVie Deutschland GmbH & Co. KG, Development Sciences, 67061, Ludwigshafen, Germany
| | - Samuel O Kyeremateng
- AbbVie Deutschland GmbH & Co. KG, Development Sciences, 67061, Ludwigshafen, Germany.
| | | | | | - Lukasz Miroslaw
- Azure High Performance Computing and Artificial Intelligence, Microsoft Switzerland, The Circle 02, 8058, Zurich, Switzerland
| |
Collapse
|
2
|
Karalis K, Karalis N, Karkalos N, Ntallis Ν, Antipas GSE, Xenidis A. Three-dimensional computational fluid dynamics analysis of an electric submerged arc furnace. Sci Rep 2021; 11:17637. [PMID: 34480043 PMCID: PMC8417227 DOI: 10.1038/s41598-021-96085-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2021] [Accepted: 06/30/2021] [Indexed: 11/09/2022] Open
Abstract
A computational fluid dynamics (CFD) method is proposed to analyze the operation of a submerged electric arc furnace (SAF) used in ferronickel production. A three-dimensional mathematical model was used for the time-dependent solution of the fluid flow, heat transfer and electromagnetic phenomena. The slag's physical properties, which play a crucial role in the SAF operation, were previously determined using classical molecular dynamics simulations and empirical relationships. The analysis revealed that the main slag properties affecting SAF operation are density, viscosity and electrical conductivity-the latter two being mutually dependent. The high electrical conductivity values of the slag favor melting via the high Joule heat produced within the slag region. Calculation of the dimensionless Péclet and Reynolds numbers revealed that the slag velocities play a decisive role in heat transfer and further indicate that the slag flow is laminar. The average slag velocity calculated 0.0001 m/s with maxima in the vicinity of the electrodes.
Collapse
Affiliation(s)
- K Karalis
- Institute of Geological Sciences, University of Bern, 3012, Bern, Switzerland.
| | - N Karalis
- Mechanical Engineer, Mining and Metallurgical Engineer, National Technical University of Athens, Athens, Greece
| | - N Karkalos
- School of Mechanical Engineering, National Technical University of Athens, Zografou Campus, 15780, Athens, Greece
| | - Ν Ntallis
- Molecular Modelling Laboratory, Park Innovaare, 5234, Villigen, Switzerland
| | - G S E Antipas
- Molecular Modelling Laboratory, Park Innovaare, 5234, Villigen, Switzerland
| | - A Xenidis
- Mechanical Engineer, Mining and Metallurgical Engineer, National Technical University of Athens, Athens, Greece
| |
Collapse
|
3
|
Antipas GS, Karalis KT. Direct determination of amorphous number density from the reduced pair distribution function. MethodsX 2019; 6:601-605. [PMID: 30984568 PMCID: PMC6446046 DOI: 10.1016/j.mex.2019.03.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Accepted: 03/05/2019] [Indexed: 11/16/2022] Open
Abstract
The inference of amorphous bulk density, while straightforward for nonporous, soluble materials, may present a formidable challenge in some of the most important classes of industrial applications, involving melts, porous solids, and non-soluble organic pharmaceuticals, with varied implications depending on the material's level of technological interest. Within nanotechnology and the life sciences in particular, accurate determination of amorphous true density is a frequent requirement and a regular puzzle, when, e.g., neither the Archimedean principle nor gas pycnometry may be applied, the former being only applicable to insoluble compounds, while the latter yielding skeletal density - an overestimate of true density to the extent of blind pores - and its efficiency is affected by the choice of the gas medium. In these cases, it is feasible to infer amorphous density from diffraction experiments through the use of the reduced Pair Distribution Function (PDF). Although an estimate of crystalline density has been known to be possible via the PDF shape, here we outline a new method extending this facility to include the estimation of amorphous density. •Amorphous density may be inferred from the position of a local minimum of the reduced PDF profile, the latter extracted via a Fourier transformation of collected diffraction intensity.•The PDF minimum is located within the PDF range bounded by rmin = 2π/Qmax and the position of the first coordination peak, where Qmax is the maximum length of the scattering vector achieved in the diffraction experiment.•Amorphous density is calculated as the ratio of the value of the reduced PDF at the local minimum, divided by the term 4πr, where r is the real space coordinate of the local minimum.
Collapse
Affiliation(s)
| | - Konstantinos T. Karalis
- Laboratory for Scientific Computing and Modelling, Paul Scherrer Institute, 5232 Villigen PSI, Switzerland
| |
Collapse
|
4
|
Karalis K, Karkalos N, Antipas GSE, Xenidis A. Pragmatic analysis of the electric submerged arc furnace continuum. ROYAL SOCIETY OPEN SCIENCE 2017; 4:170313. [PMID: 28989738 PMCID: PMC5627078 DOI: 10.1098/rsos.170313] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/06/2017] [Accepted: 08/02/2017] [Indexed: 06/07/2023]
Abstract
A transient mathematical model was developed for the description of fluid flow, heat transfer and electromagnetic phenomena involved in the production of ferronickel in electric arc furnaces. The key operating variables considered were the thermal and electrical conductivity of the slag and the shape, immersion depth and applied electric potential of the electrodes. It was established that the principal stimuli of the velocities in the slag bath were the electric potential and immersion depth of the electrodes and the thermal and electrical conductivities of the slag. Additionally, it was determined that, under the set of operating conditions examined, the maximum slag temperature ranged between 1756 and 1825 K, which is in accordance with industrial measurements. Moreover, it was affirmed that contributions to slag stirring due to Lorentz forces and momentum forces due to the release of carbon monoxide bubbles from the electrode surface were negligible.
Collapse
Affiliation(s)
- K. Karalis
- School of Mining Engineering and Metallurgy, National Technical University of Athens, Zografou Campus, Athens 15780, Greece
| | | | | | | |
Collapse
|
5
|
Bona-fide method for the determination of short range order and transport properties in a ferro-aluminosilicate slag. Sci Rep 2016; 6:30216. [PMID: 27455915 PMCID: PMC4960545 DOI: 10.1038/srep30216] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2015] [Accepted: 06/30/2016] [Indexed: 11/09/2022] Open
Abstract
The thermodynamics, structural and transport properties (density, melting point, heat capacity, thermal expansion coefficient, viscosity and electrical conductivity) of a ferro-aluminosilicate slag have been studied in the solid and liquid state (1273-2273 K) using molecular dynamics. The simulations were based on a Buckingham-type potential, which was extended here, to account for the presence of Cr and Cu. The potential was optimized by fitting pair distribution function partials to values determined by Reverse Monte Carlo modelling of X-ray and neutron diffraction experiments. The resulting short range order features and ring statistics were in tight agreement with experimental data and created consensus for the accurate prediction of transport properties. Accordingly, calculations yielded rational values both for the average heat capacity, equal to 1668.58 J/(kg·K), and for the viscosity, in the range of 4.09-87.64 cP. The potential was consistent in predicting accurate values for mass density (i.e. 2961.50 kg/m(3) vs. an experimental value of 2940 kg/m(3)) and for electrical conductivity (5.3-233 S/m within a temperature range of 1273.15-2273.15 K).
Collapse
|
6
|
Antipas GSE, Germenis AE. Atomic Coordination Reflects Peptide Immunogenicity. Front Mol Biosci 2016; 2:77. [PMID: 26793714 PMCID: PMC4707223 DOI: 10.3389/fmolb.2015.00077] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2015] [Accepted: 12/19/2015] [Indexed: 11/13/2022] Open
Abstract
We demonstrated that the immunological identity of variant peptides may be accurately predicted on the basis of atomic coordination of both unprotonated and protonated tertiary structures, provided that the structure of the native peptide (index) is known. The metric which was discovered to account for this discrimination is the coordination difference between the variant and the index; we also showed that increasing coordination difference in respect to the index was correlated to a correspondingly weakening immunological outcome of the variant. Additionally, we established that this metric quickly seizes to operate beyond the peptide scale, e.g., within a coordination shell inclusive of atoms up to a distance of 7 Å away from the peptide or over the entire pMHC-TCR complex. Analysis of molecular orbital interactions for a range of formal charges further revealed that the N-terminus of the agonists was always able to sustain a stable ammonium (NH3+) group which was consistently absent in antagonists. We deem that the presence of NH3+ constitutes a secondary observable with a biological consequence, signifying a change in T cell activation. While our analysis of protonated structures relied on the quantum chemical relaxation of the H species, the results were consistent across a wide range of peptide charge and spin polarization conditions.
Collapse
Affiliation(s)
- Georgios S E Antipas
- Division of Materials Technology, National Technical University of AthensAthens, Greece; Department of Molecular Medicine, Hellenic Pasteur InstituteAthens, Greece
| | - Anastasios E Germenis
- Department of Molecular Medicine, Hellenic Pasteur InstituteAthens, Greece; Department of Immunology and Histocompatibility, School of Medicine, University of ThessalyLarissa, Greece
| |
Collapse
|
7
|
Antipas GSE, Germenis AE. The quantum chemical causality of pMHC-TCR biological avidity: Peptide atomic coordination data and the electronic state of agonist N termini. Data Brief 2015. [PMID: 26217741 PMCID: PMC4510142 DOI: 10.1016/j.dib.2015.02.021] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
The quantum state of functional avidity of the synapse formed between a peptide-Major Histocompatibility Complex (pMHC) and a T cell receptor (TCR) is a subject not previously touched upon. Here we present atomic pair correlation meta-data based on crystalized tertiary structures of the Tax (HTLV-1) peptide along with three artificially altered variants, all of which were presented by the (Class I) HLA-A201 protein in complexation with the human (CD8(+)) A6TCR. The meta-data reveal the existence of a direct relationship between pMHC-TCR functional avidity (agonist/antagonist) and peptide pair distribution function (PDF). In this context, antagonist peptides are consistently under-coordinated in respect to Tax. Moreover, Density Functional Theory (DFT) datasets in the BLYP/TZ2P level of theory resulting from relaxation of the H species on peptide tertiary structures reveal that the coordination requirement of agonist peptides is also expressed as a physical observable of the protonation state of their N termini: agonistic peptides are always found to retain a stable ammonium (NH3 (+)) terminal group while antagonist peptides are not.
Collapse
Affiliation(s)
- Georgios S E Antipas
- Division of Materials Technology, National Technical University of Athens, Zografou Campus, Athens 15780, Greece
| | - Anastasios E Germenis
- Department of Immunology & Histocompatibility, School of Medicine, University of Thessaly, Biopolis, Larissa 41110, Greece
| |
Collapse
|
8
|
Antipas GSE, Germenis AE. Quantum chemical calculations predict biological function: the case of T cell receptor interaction with a peptide/MHC class I. Front Chem 2015; 3:9. [PMID: 25713797 PMCID: PMC4322848 DOI: 10.3389/fchem.2015.00009] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Accepted: 01/22/2015] [Indexed: 11/13/2022] Open
Abstract
A combination of atomic correlation statistics and quantum chemical calculations are shown to predict biological function. In the present study, various antigenic peptide-Major Histocompatibility Complex (pMHC) ligands with near-identical stereochemistries, in complexation with the same T cell receptor (TCR), were found to consistently induce distinctly different quantum chemical behavior, directly dependent on the peptide's electron spin density and intrinsically expressed by the protonation state of the peptide's N-terminus. Furthermore, the cumulative coordination difference of any variant in respect to the native peptide was found to accurately reflect peptide biological function and immerges as the physical observable which is directly related to the immunological end-effect of pMHC-TCR interaction.
Collapse
Affiliation(s)
- Georgios S E Antipas
- Division of Materials Technology, National Technical University of Athens Athens, Greece
| | - Anastasios E Germenis
- Department of Immunology and Histocompatibility, School of Medicine, University of Thessaly Biopolis, Larissa, Greece
| |
Collapse
|
9
|
Antipas GSE, Temleitner L, Karalis K, Pusztai L, Xenidis A. Atomic order and cluster energetics of a 17 wt% Si-based glass versus the liquid phase. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2013; 25:454206. [PMID: 24141046 DOI: 10.1088/0953-8984/25/45/454206] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
Aerodynamic levitation of a multicomponent 17 wt% Si glass formed by rapid quenching of the melt phase was studied by high resolution x-ray diffraction (XRD) and reverse Monte Carlo (RMC) modelling. The main local atomic order features comprised interactions between Si, Fe and Mg polyhedra, the stereochemistry of which was on a par with the literature. Both the glass and the liquid state appeared to consist of the same fundamental Si-O, Fe-O and Mg-O clusters, with only the relative number of each varying between the two. Transition from liquid to glass involved a three-fold decrease in uncoordinated O (to within the first minimum of the total g(r)) and a marked increase of Fe-Si-Mg polyhedra bridging O. Octahedral Fe coordination was not suggested by the RMC data. All-electron open-shell density functional theory (DFT) calculations of the most prominent clusters suggested independence between the Fe oxidation state and its polyhedra O-coordination. Of secondary thermodynamic importance were indications of network-forming Fe(2+) and Fe(3+) distorted trigonal and tetrahedral polyhedra. In all occasions, the Fe ferrous and ferric states involved comparable binding energies within similar clusters which indicate a dynamic equilibrium between the two.
Collapse
Affiliation(s)
- G S E Antipas
- School of Mining Engineering and Metallurgy, National Technical University of Athens, Zografou Campus, Athens 15780, Greece
| | | | | | | | | |
Collapse
|