Buffer solutions in drug formulation and processing: How pKa values depend on temperature, pressure and ionic strength

ANI neural network potentials for small molecule pKa prediction

The pKa value of a molecule is of interest to chemists across a broad spectrum of fields including pharmacology, environmental chemistry and theoretical chemistry. Determination of pKa values can be accomplished through several experimental methods such as NMR techniques and titration together with computational techniques such as DFT calculations. However, all of these methods remain time consuming and computationally expensive. In this work we develop a method for the rapid calculation of pKa values of small molecules which utilises a combination of neural network potentials, low energy conformer searches and thermodynamic cycles. We show that neural network potentials trained on different phase and charge states can be employed in tandem to predict the full thermodynamic energy cycle of molecules. Focusing here on imidazolium derived carbene species, the method utilised can easily be extended to other functional groups of interest such as amines with further training.

Probing Chemical Equilibrium in Frozen Sodium Phosphate Buffer Solution by 31P Solid-State NMR

Phosphate buffers are crucial for cryopreservative stability in pharmaceuticals, food processing, biomedical sciences, and biology. However, their freeze concentrates lack quantitative characterization, especially regarding the physicochemical properties of phosphate salt species in equilibrium at subzero temperatures. This study employs 31P solid-state NMR (ssNMR) to analyze frozen sodium phosphate (NaP) solutions, providing insights into phase composition, ionic strength, and pH. For the first time, we have directly quantified phosphate species in frozen NaP buffer, including crystallized disodium phosphate dodecahydrate (Na2HPO4·12H2O) content and the concentrations of H2PO4- and HPO42- in the freeze concentrate. This enabled the calculation of the pH as well as the ionic strength in the freeze concentrate. Trehalose effectively mitigated pH shifts in buffer solutions by preventing the selective crystallization of salt, a spectroscopic phenomenon not previously observed experimentally.

Thermal Stability of Protic Ionic Liquids

While protic ionic liquids (ILs) have found great success as solvents for a broad range of applications, little is known about their degradation when exposed to temperatures above ambient for extended periods of time. Here, we report the thermal stability of six protic ILs, namely, ethylammonium nitrate, ethylammonium formate, ethylammonium acetate, ethanolammonium nitrate, ethanolammonium formate, and ethanolammonium acetate. The effect of heating each ionic liquid to 60 °C for 1 h or 1 week (sealed or open to the atmosphere) was evaluated by considering the changes to water content, pH, mass, thermal phase transitions, and molecular structure after each treatment. Heating each of the six ILs when sealed led to measurable shifts in their water content and 10 wt % pH, but there was no significant change in their mass, thermal phase transitions according to differential scanning calorimetry (DSC), or molecular structure using proton nuclear magnetic resonance (1H NMR) spectra, indicating that the samples were largely unchanged. The samples that were heated open to the atmosphere also displayed no significant changes after 1 h but displayed significant changes after 1 week.

Fluorescence Spectra of Prototropic Forms of Fluorescein and Some Derivatives and Their Potential Use for Calibration-Free pH Sensing

Fluorescence pH sensing has proven to be efficient but with the drawback that molecules photobleach, requiring frequent calibrations. Double-emission peak molecules allow ratiometric measurements and theoretically avoid calibration. However, they are often expensive and fragile and usually have very low quantum yields. Single emission peaks such as fluorescein and derivatives are inexpensive and have very high quantum yields. Because they are single emission peaks, the pH is assumed to be derived from the ratio of emitted intensities at measured pH and at high pH values, i.e., they require frequent calibration. However, the shape of their single emitted peak evolves slightly with pH. In this paper, we first demonstrate a simple method to calculate the emission spectrum shape of each prototropic form of fluorescein (and derivatives) as well as the values of the pKas. A complete model of the evolution of the emission spectrum shape with pH is then constructed. Second, we evaluate the potential of these molecules for pH sensing by fitting the experimental spectra with the complete emission model. The method is applied to fluorescein, FITC and FAM. Depending on the molecule, pH can be measured from pH 1.9 to pH 7.3 with standard deviations between 0.06 and 0.08 pH units. Estimating pH and pKas from shape instead of intensity allows calibration-free measurements even with single-emission peak molecules.

The curious case of proton migration under pressure in the malonic acid and 4,4'-bipyridine cocrystal

In the search for new active pharmaceutical ingredients, the precise control of the chemistry of cocrystals becomes essential. One crucial step within this chemistry is proton migration between cocrystal coformers to form a salt, usually anticipated by the empirical ΔpKa rule. Due to the effective role it plays in modifying intermolecular distances and interactions, pressure adds a new dimension to the ΔpKa rule. Still, this variable has been scarcely applied to induce proton-transfer reactions within these systems. In our study, high-pressure X-ray diffraction and Raman spectroscopy experiments, supported by DFT calculations, reveal modifications to the protonation states of the 4,4'-bipyridine (BIPY) and malonic acid (MA) cocrystal (BIPYMA) that allow the conversion of the cocrystal phase into ionic salt polymorphs. On compression, neutral BIPYMA and monoprotonated (BIPYH+MA-) species coexist up to 3.1 GPa, where a phase transition to a structure of P21/c symmetry occurs, induced by a double proton-transfer reaction forming BIPYH22+MA2-. The low-pressure C2/c phase is recovered at 2.4 GPa on decompression, leading to a 0.7 GPa hysteresis pressure range. This is one of a few studies on proton transfer in multicomponent crystals that shows how susceptible the interconversion between differently charged species is to even slight pressure changes, and how the proton transfer can be a triggering factor leading to changes in the crystal symmetry. These new data, coupled with information from previous reports on proton-transfer reactions between coformers, extend the applicability of the ΔpKa rule incorporating the pressure required to induce salt formation.

Impact of Sinapic Acid on Bovine Serum Albumin Thermal Stability

The thermal stability of bovine serum albumin (BSA) in Tris buffer, as well as the effect of sinapic acid (SA) on protein conformation were investigated via calorimetric (differential scanning microcalorimetry-μDSC), spectroscopic (dynamic light scattering-DLS; circular dichroism-CD), and molecular docking approaches. μDSC data revealed both the denaturation (endotherm) and aggregation (exotherm) of the protein, demonstrating the dual effect of SA on protein thermal stability. With an increase in ligand concentration, (i) protein denaturation shifts to a higher temperature (indicating native form stabilization), while (ii) the aggregation process shifts to a lower temperature (indicating enhanced reactivity of the denatured form). The stabilization effect of SA on the native structure of the protein was supported by CD results. High temperature (338 K) incubation induced protein unfolding and aggregation, and increasing the concentration of SA altered the size distribution of the protein population, as DLS measurements demonstrated. Complementary information offered by molecular docking allowed for the assessment of the ligand binding within the Sudlow's site I of the protein. The deeper insight into the SA-BSA interaction offered by the present study may serve in the clarification of ligand pharmacokinetics and pharmacodynamics, thus opening paths for future research and therapeutic applications.

Formic and acetic acid pKa values increase under nanoconfinement

Organic acids are prevalent in the environment and their acidity and the corresponding dissociation constants can change under varying environmental conditions. The impact of nanoconfinement (when acids are confined within nanometer-scale domains) on physicochemical properties of chemical species is poorly understood and is an emerging field of study. By combining infrared and Raman spectroscopies with molecular dynamics (MD) simulations, we quantified the effect of nanoconfinement in silica nanopores on one of the fundamental chemical reactions-the dissociation of organic acids. The pKa of formic and acetic acids confined within cylindrical silica nanopores with 4 nm diameters were measured. MD models were constructed to calculate the shifts in the pKa values of acetic acid nanoconfined within 1, 2, 3, and 4 nm silica slit pores. Both experiments and MD models indicate a decrease in the apparent acid dissociation constants (i.e., increase in the pKa values) when organic acids are nanoconfined. Therefore, nanoconfinement stabilizes the protonated species. We attribute this observation to (1) a decrease in the average dielectric response of nanoconfined aqueous solutions where charge screening may be decreased; or (2) an increase in proton concentration inside nanopores, which would shift the equilibrium towards the protonated form. Overall, the results of this study provide the first quantification of the pKa values for nanoconfined formic and acetic acids and pave the way for a unifying theory predicting the impact of nanoconfinement on acid-base chemistry.

Calculating the Aqueous pKa of Phenols: Predictions for Antioxidants and Cannabinoids

We aim to develop a theoretical methodology for the accurate aqueous pK_a prediction of structurally complex phenolic antioxidants and cannabinoids. In this study, five functionals (M06-2X, B3LYP, BHandHLYP, PBE0, and TPSS) and two solvent models (SMD and PCM) were combined with the 6-311++G(d,p) basis set to predict pK_a values for twenty structurally simple phenols. None of the direct calculations produced good results. However, the correlations between the calculated Gibbs energy difference of each acid and its conjugate base, ΔGaq(BA)°=ΔGaqA-°-ΔGaq(HA)°, and the experimental aqueous pK_a values had superior predictive accuracy, which was also tested relative to an independent set of ten molecules of which six were structurally complex phenols. New correlations were built with twenty-seven phenols (including the phenols with experimental pK_a values from the test set), which were used to make predictions. The best correlation equations used the PCM method and produced mean absolute errors of 0.26-0.27 pK_a units and R² values of 0.957-0.960. The average range of predictions for the potential antioxidants (cannabinoids) was 0.15 (0.25) pK_a units, which indicates good agreement between our methodologies. The new correlation equations could be used to make pK_a predictions for other phenols in water and potentially in other solvents where they might be more soluble.

Characterizing Oscillatory and Excitability Regimes in a Protein-Free Lipid Membrane

Observations of electric potential oscillations in artificial lipid bilayers near the order-disorder transition indicate the existence of a stable limit cycle and, therefore, the possibility of producing excitable signals close to the bifurcation. We present a theoretical investigation of membrane oscillatory and excitability regimes induced by an increase in ion permeability at the order-disorder transition. The model accounts for the coupled effects of state-dependent permeability, membrane charge density, and hydrogen ion adsorption. A bifurcation diagram shows a transition between fixed-point and limit cycle solutions, enabling both oscillatory and excitability responses at different values of the acid association parameter. Oscillations are identified in terms of the membrane state, electric potential difference, and ion concentration near the membrane. The emerging voltage and time scales agree with measurements. Excitability is demonstrated by applying an external electric current stimulus, and the emerging signals display a threshold response and the appearance of repetitive signals upon using a long-lasting stimulus. The approach highlights the important role of the order-disorder transition, enabling membrane excitability in the absence of specialized proteins.

A systematic review of various pKa determination techniques

The pK_a values of functional groups is crucial in determining the pharmacokinetic properties of a drug, affecting its absorption and thus bioavailability. This physicochemical property is also vital for the designing of drug excipients and vehicles. There are currently 13 known methods of determining a pK_a value, namely: potentiometric titration, spectrometry, fluorometry, NMR, HPLC, conductometry, electrophoresis, voltammetry, solubility, partition coefficient, calorimetry, computational, and surface tension. Some of these techniques are more widely utilized and well-established compared to others, with each having their inherent advantages and disadvantages. This review discusses each of the aforementioned techniques with emphasis on their pros and cons.

Two Opposing Effects of Monovalent Cations on the Stability of i-Motif Structure

At acidic pH, cytosine-rich single-stranded DNA can be folded into a tetraplex structure called i-motif (iM). In recent studies, the effect of monovalent cations on the stability of iM structure has been addressed, but a consensus about the issue has not been reached yet. Thus, we investigated the effects of various factors on the stability of iM structure using fluorescence resonance energy transfer (FRET)-based analysis for three types of iM derived from human telomere sequences. We confirmed that the protonated cytosine-cytosine (C:C⁺) base pair is destabilized as the concentration of monovalent cations (Li⁺, Na⁺, K⁺) increases and that Li⁺ has the greatest tendency of destabilization. Intriguingly, monovalent cations would play an ambivalent role in iM formation by making single-stranded DNA flexible and pliant for an iM structure. In particular, we found that Li⁺ has a notably greater flexibilizing effect than Na⁺ and K⁺. All taken together, we conclude that the stability of iM structure is controlled by the subtle balance of the two counteractive effects of monovalent cations: electrostatic screening and disruption of cytosine base pairing.

Degradation kinetics of artesunate for the development of an ex-tempore intravenous injection

Background

Artesunate is recommended by the World Health Organization (WHO) for parenteral treatment of severe Plasmodium falciparum malaria. However, artesunate is inherently unstable in an aqueous solution and hydrolyses rapidly after its preparation for injection. Therefore, the aim of the study was to evaluate the stabilizing effects of phosphate buffer and mannitol against short-term (ex-tempore) artesunate hydrolysis.

Methods

A HPLC–UV isocratic method was developed using a reversed-phase fused core column (HALO RP-C18) and a mobile phase consisting of a mixture of 45% ammonium formate 10 mM in water (pH 4.5) and 55% methanol. Artesunate was formulated as aqueous solutions using a design of experiment (DOE) to investigate the artesunate stabilizing effects of pH (8–10), phosphate buffer strength (0.3–0.5 M), and mannitol (0–0.22 mmol/mL). The solutions were incubated at predefined temperatures (5, 25, and 40 °C) with subsequent analysis. Arrhenius equation was applied to model and evaluate the stability results.

Results

The developed HPLC-based method using fused-core stationary phase allowed to selectively quantify artesunate in the presence of its main hydrolysis degradants; namely β-dihydroartemisinin (β-DHA) and α-dihydroartemisinin (α-DHA) within 10 min. By applying the Arrhenius equation, the rate of hydrolysis of the drug increased approximately by 3.4 as the temperature raised by 10 °C. Buffer strength was found to be the main factor affecting the hydrolysis rate constants at 5 and 25 °C (p < 0.05), the activation energy (p = 0.009), and the frequency factor (p = 0.045). However, the effect of the buffer was predominant on the activation energy and hydrolysis rate constants, revealing its stabilizing effect on the drug at lower buffer strength (0.3 M). Within the investigated range (pH = 8–10), pH was found to influence the activation energy, with a positive stabilizing effect in the pH range of 8–9. The addition of mannitol as stabilizing agent into artesunate aqueous formulation did not show an improved response.

Conclusion

Phosphate buffer was the main stability determining factor of artesunate in the aqueous intravenous (i.v.) formulation and was found to be more effective in stabilizing artesunate at a buffer strength of 0.3 M in pH 8–9, while mannitol lacked stabilizing effect.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04278-4.

Development of liquid chromatographic and UV-visible spectrophotometric methods for determination of pKa values of folic acid antimetabolites

In this study, the values of ionization/protonation constants (pK_a) of pemetrexed and raltitrexed, which are folic acid antimetabolites, were determined by reversed-phase liquid chromatography (RPLC) and UV-visible spectrophotometric methods at 25 °C. The pKass values were evaluated by retention time (t_R) in binary acetonitrile-water mixtures with acetonitrile (ACN) contents of 13%, 15%, 17% and 20% (v/v) for the RPLC method and 13%, 15% and 17% (v/v) for the spectrophotometric method. The aqueous pK_a values (pKaww) of the studied compounds were calculated from the pKass value using the macroscopic parameters (mole fraction and dielectric constant), which play an important role in solvent properties. In addition, the degree of ionization of pemetrexed and raltitrexed was calculated using the pKaww values.

OUP accepted manuscript

A simple, robust stability indicating RP-HPLC method was developed for simultaneous quantification of serdexmethyl phenidate and dexmethyl phenidate in a fixed capsule dosage form. This is the first method to be reported for simultaneous estimation and quantification of degradation products produced from forced stressing of the dosage form as per ICH guidelines. The chromatographic separation was attained on Waters X-terra C18 column using a mixture of trifluoro acetic acid and acetonitrile (70:30 v/v) as mobile phase with a flow rate of 1 mL, monitored at 265 nm over a run time of 10 min. Serdexmethyl phenidate and dexmethyl phenidate were eluted with retention times of 2.71 min and 7.33 min, respectively. The method displayed linear responses in the range of 4.2-63 μg/mL (0.9994) for serdexmethyl phenidate and 0.9 to 3.5 μg/mL (0.9998) for dexmethyl phenidate, respectively. The percentage recoveries of the two drugs were found within the acceptable limits. Forced degradation was conducted and showed considerable degradation in various stress conditions. It also confirms the specificity of the method as no interference peaks were observed concerning for to stress products. This method can be routinely used in quality control labs for simultaneous determination of respective drugs in marketed dosage form.

Evaluation of ionic equilibria in mixed-buffer isothermal titration calorimetry and continuously stirred tank reactors

Determination of an equilibrium pH value in complex aqueous solution and deconvolution of this equilibrium to evaluate phenomena related to mixing, dilution, or progress of reaction is increasingly important in areas ranging from water quality to pharmaceutical formulations and manufacturing. Linearization of pH problems by simple algebraic substitution enables equilibria within complex buffered aqueous solutions to be modeled as an eigenvalue problem. This formulation approach makes rigorous determination of equilibrium pH values and reactor dynamics more accessible than with previous calculation methods, even when activity coefficients and non-ideality are considered. This work demonstrates how such calculations can enable detailed modeling of enthalpic changes in an isothermal titration calorimeter. In support of this work, the acid dissociation constants for three furancarboxylic acids (2-furancarboxylic acid, FA; 5-formyl-2-furancarboxylic acid, FFA; and 2,5-furandicarboxylic acid, FDCA), two of them novel, were determined and compared with multi-wavelength ultraviolet-visible spectrophotometry. The thermodynamic pKa values were determined to be 3.1 for FA, 2.2 for FFA, and 2.1 and 3.4 for the first and second ionization steps of FDCA, respectively.

Enantioseparation of ketoconazole and miconazole by capillary electrophoresis and a study on their inclusion interactions with β-cyclodextrin and derivatives

A chiral separation method coupled with capillary electrophoresis (CE) analysis for ketoconazole and miconazole enantiomers using chiral selectors such as β-cyclodextrin (β-CD) and hydroxypropyl-β-CD (HP-β-CD) was developed in this study, which included the optimisation, validation and application of the method on the antifungal cream samples. The formation of inclusion complex between the hosts (β-CD and HP-β-CD) and guests (ketoconazole and miconazole) were compared and analysed using ultraviolet-visible spectrophotometry, nuclear magnetic resonance (NMR) spectroscopy and molecular docking methods. Results from the study showed that in a concentration that ranged between 0.25 and 50 mg L^-1 , the linear calibration curves of each enantiomer had a high coefficient of regression (R² > 0.999), low limit of detection (0.075 mg L^-1 ) and low limit of quantification (0.25 mg L^-1 ). The relative standard deviation (RSD) of the intraday and interday analyses ranged from 0.79% to 8.01% and 3.30% to 11.43%, respectively, while the recoveries ranged from 82.0% to 105.7% (RSD < 7%, n = 3). The most probable structure of the inclusion complexes was proposed based on the findings from the molecular docking studies conducted using the PatchDock server.

Determination of acidity constants for weak acids and bases by isothermal titration calorimetry

The advantage of isothermal titration calorimetry (ITC) to determine the acid dissociation constant (pK_a value) is the simultaneous determination of the binding constant and binding enthalpy, as well as being precise and easy to use. The pK_a can be calculated from the binding constant, and the temperature dependency of the pK_a can be calculated from the binding enthalpy. The use of ITC to study protonation reactions is less common compared to its more conventional use of studying macromolecules and ligands. Water will influence the equilibrium due to autoionization, meaning that both the conjugate base and acid will exist in the sample cell at the beginning of the experiment. These differences are accounted for by optimizing the theoretical model used to estimate the binding constant and binding enthalpy. Through simulations and experimental measurements, we show that ITC can be used to determine the pK_a for ibuprofen, ascorbic acid, 2-morpholin-4-ylethanesulfonic acid and paracetamol. The pK_a values were consistent with potentiometric or spectrophotometric determinations as well as literature values. Optimizing the theoretical model does not lead to an improved determination, so the "one set of sites" model is adequate for the determination of pK_a values.