Measurement and modeling of metoclopramide hydrochloride (anti-emetic drug) solubility in supercritical carbon dioxide

Knowledge of drug solubility data in supercritical carbon dioxide (SC-CO₂) is a fundamental step in producing nano and microparticles through supercritical fluid technology. In this work, for the first time, the solubility of metoclopramide hydrochloride (MCP) in SC-CO₂ was measured in pressure and temperature range of 12 to 27 MPa and 308 to 338 K, respectively. The results represented a range mole fractions of 0.15 × 10^-5 to 5.56 × 10^-5. To expand the application of the obtained data, six semi-empirical models and three models based on the Peng-Robinson equation of state (PR + VDW, PR + WS + Wilson and PR + MHV1 + COSMOSAC) with different mixing rules and various ways to describe intermolecular interactions were investigated. Furthermore, total enthalpy, sublimation enthalpy and solvation enthalpy relevant to MCP solvating in SC-CO₂ were estimated.     

Development of machine learning model and analysis study of drug solubility in supercritical solvent for green technology development

The current research is focused on development of machine learning model for estimation of pharmaceutical solubility in supercritical CO₂ as the green solvent. The main aim is to assess the suitability of supercritical processing for preparation of nanomedicine. Oxaprozin was taken as model drug for the solubility measurements, and its solubility was determined at different operational conditions by variation of temperature and pressure of the process. Artificial Neural Network (ANN) model was implemented for simulation of the drug solubility, and the best model was obtained with R² greater than 0.99 for the training and validation as well. The tested model was then exploited to understand the process, and it turned out that both pressure and temperature had major and considerable influence on the solubility of Oxaprozin in supercritical carbon dioxide as solvent. However, the effect of pressure was shown to be more significant on the solubility compared to the effect of pressure, which was attributed to the effect of pressure on the density of the supercritical solvent. The developed ANN model was indicated to be robust in estimating the values of drug solubility in wide range of conditions which can save time and cost of the measurements.     

Machine learning simulation of pharmaceutical solubility in supercritical carbon dioxide: Prediction and experimental validation for busulfan drug

An artificial intelligence-based predictive model was developed using a support vector machine to investigate the solubility data of the drug Busulfan drug in supercritical carbon dioxide. The data for simulations were collected from literature. The model was trained and implemented in order to determine the correlation between the solubility values and the input parameters, namely, temperature and pressure. These parameters were used as the inputs as they are known to have a significant effect on the solubility of Busulfan in supercritical carbon dioxide. In the artificial intelligence model, a polynomial model with kernel function was applied to the data, and the model’s findings were compared with measured data for fitting. Good agreement was observed between the model’s outputs and the measured data with coefficient of determination greater than 0.99.     

A comparison between neural network method and semi empirical equations to predict the solubility of different compounds in supercritical carbon dioxide

Accuracy of seven semi empirical equations for the estimation of solubility of 30 different compounds in supercritical carbon dioxide has been compared with a new neural network method. To base this comparison on a fair basis, a unique set of experimental data was used for both optimization of semi empirical equations’ parameters and training, validation and testing of neural network. Results showed that neural network method with an average relative deviation of about 5.3% was more accurate than the best semi empirical equation with an average relative deviation of about 15.96% for same compounds. It was also found that the average relative deviation of semi empirical equations varies sharply among different compounds, while this quantity is less dependent on material type for neural network method.     

Solubilities of cinnamic acid,phenoxyacetic acid and 4-methoxyphenylacetic acid in supercritical carbon dioxide

The solid solubilities of cinnamic acid, phenoxyacetic acid and 4-methoxyphenylacetic acid in supercritical carbon dioxide were measured using a semi-flow apparatus. The experiments were taken at 308.2, 318.2 and 328.2 K. The pressure range was from 11 to 24 MPa. These data were confirmed as equilibrium solid solubilities based on a plug flow mass transfer model. The solid solubilities were further correlated using the equations of state or semi-empirical models. The correlation results are satisfactory with optimally fitted binary interaction parameters in the Peng–Robinson equation of state.     

Measurements for the solid solubilities of antipyrine, 4-aminoantipyrine and 4-dimethylaminoantipyrine in supercritical carbon dioxide

The solid solubilities of three active pharmaceutical ingredients (APIs) of antipyrine, 4-aminoantipyrine and 4-dimethylaminoantipyrine in supercritical carbon dioxide were measured by a semi-flow apparatus. The experiments were taken at 308.2, 318.2 and 328.2 K. The pressure range was from 10 to 22 MPa. These experimental results were correlated by the semi-empirical models of Mendez–Santiago–Teja and Chrastil. A solution model was also employed to fit the measured data. The average absolute relative deviation in solid solubility from semi-empirical models was 4–6%, and that from the solution model was 5–8%. The measured data satisfied the self-consistency test, and the parameters in the semi-empirical models are feasible for data extrapolation.     

Experimental study on equilibrium solubility (at low partial pressures), density, viscosity and corrosion rate of carbon dioxide in aqueous solutions of ascorbic acid

In this paper, ascorbic acid as a new carbon dioxide (CO₂) absorbent was investigated. The equilibrium solubility of CO₂ into 0.5, 1 and 1.5 mol dm⁻³ (M) aqueous ascorbic acid solutions were measured experimentally with a stirred batch reactor at total atmospheric pressure over the CO₂ partial pressure ranging from 0 to 45 kPa and temperatures between 298 and 313 K. The results of the gas solubility are presented as loading capacity (mol CO₂/mol ascorbic acid) as function of partial pressure of CO₂ for all experimental runs. Experimental results showed that solubility of CO₂ increases with increase in molar concentration of ascorbic acid solution at a given temperature and decreases with increase in temperature at a given concentration. The densities and viscosities of the ascorbic acid solutions were measured at the same conditions of the solubility measurement. Some corrosion rate tests were also performed on carbon steel at temperature of 308 K. It was observed that viscosity and corrosion rate increase when the molar concentration of ascorbic acid solution increases.     

Solvent effect on distribution ratio of Pd(II) in supercritical carbon dioxide extraction and solvent extraction using 2-methyl-8-quinolinol

The distribution ratio (D_M) of Pd(II) by the extraction with 2-methyl-8-quinolinol (HMQ) was determined using the supercritical carbon dioxide medium (SF-CO₂) and organic solvent media such as perfluoro-methylcyclohexane, heptane, cyclohexane, carbon tetrachloride and benzene. From experimental results of the slopes of log D_M versus pH plot and log D_M versus HMQ concentration plot, the extracted species both in the SF-CO₂ extraction (SFE) and the solvent extraction (SE) were determined to be Pd(MQ)₂. The distribution constant of HMQ (K_D,HMQ) in the SFE and SE systems were determined from the dependence of the distribution ratio of HMQ (D_HMQ) on the pH. A linear relationship was observed between log K_D,HMQ and the solubility parameter (δ) of the extraction medium based on the regular solution theory in both the SFE using SF-CO₂ at the pressure of 8.5–40 MPa and the SE systems. The difference in the slope of the log K_D,HMQ versus δ plot between the SFE and the SE systems is attributable to the extent of the specific interaction of the solute HMQ with the solvent molecules, i.e. CO₂ molecules and the organic solvent molecules. The D_M versus δ plot obtained under a given extraction condition using SF-CO₂ (11–40 MPa) and organic solvents showed clear linearity. The D_M obtained using SF-CO₂ at relatively low pressure range from 8.5 to 11 MPa was independent of the pressure and the δ of SF-CO₂, which coincides with the experimental fact that the solubility of Pd(MQ)₂ in the SF-CO₂ at 8.5–11 MPa was practically constant.     

Experimental and thermodynamic analysis on the solubility of valsartan in ethyl acetate + (butanone,isopropyl ether) binary solvent mixtures (from 278.15 to 323.15 K)

The solubility of valsartan in ethyl acetate + (butanone, isopropyl ether) binary solvent mixtures was measured at temperatures T = 278.15–323.15 K and pressure p = 0.1 MPa with a laser monitoring dynamic technique by a synthetic method. The experimental data were regressed by the modified Apelblat equation, the general single model and the hybrid model. The experimental data are well correlated with the above models because the mean deviations (MDs) are less than 3.79%. The mole fraction solubility of valsartan increases with increase in temperature and enrichment in butanone content, while it decreases with increased mole fraction of isopropyl ether at constant temperature. In addition, thermodynamic studies, including Gibbs energy, entropy and enthalpy, were calculated by van’t Hoff analysis. The results showed that the dissolution of valsartan in mixed solvents is endothermic, spontaneous and entropy-driven.     

2,5-Dichloro-1-(ROSO2)benzene [R = C6H5, C6F5, and CH2(CF2)4H]: Synthesis, molecular structure, and solubility in supercritical CO2

Highly crystalline phenyl 2,5-dichlorobenzenesulfonate (PDBS, T_melt = 86-87 °C) and pentafluorophenyl 2,5-dichlorobenzenesulfonate (FPDBS, T_melt = 120-122 °C) were synthesized. Single-crystal X-ray molecular structure determinations show that both compounds have similar three-dimensional molecular structures; however, PDBS crystals are thin platelets and FPDBS crystals form hexagonal tube-like structures that are predominately hollow at one end. PDBS crystals exhibit offset π-stacking of the phenoxy-rings that form complete two-dimensional layers each two molecules thick. Hydrogen-bonding interactions are calculated at ∼3.2 Å between the C6-hydrogen and the sulfonyl-oxygen of a neighboring molecule. On the other hand, for FPDBS, π-stacking of the dichloro-substituted ring as well as dipole-dipole interactions of the fluorinated-phenoxy rings appears to be the predominate intermolecular interactions. Neither structure exhibits any kind of side-on interaction of the phenyl rings. PDBS and FPDBS exhibit melting point depressions of 26 and 40 °C, respectively, in the presence of supercritical CO₂. Although both sulfonates exhibit high solubility in CO₂, much lower pressures are needed to dissolve FPDBS compared to PDBS. For example, at 100 °C FPDBS dissolves at 4750 psia and PDBS dissolves at 11,000 psia. The solubility data reinforce the observation that fluorinating a compound can significantly lower the conditions needed to dissolve that compound in CO₂.     

A temperature study on critical micellization concentration (CMC), solubility, and degree of counterion binding of α-sulfonatomyristic acid methyl ester in water by electroconductivity measurements

 For a sodium salt of α-sulfonatomyristic acid methyl ester (14SFNa), one of the α-SFMe series surfactants, critical micellization concentration (CMC), solubility and degree of counterion binding (β) were determined by means of electrocon-ductivity measurements at different temperatures (at every 5 °C) ranging from 15 to 50 °C. The phase diagram of 14SFNa in pure water was constructed from the CMC- and solubility-temperature data, in which the Krafft temperature (critical solution temperature) was found around 0 °C. The changes in the Gibbs energy, ΔG ⁰ _m, enthalpy, ΔH ⁰ _m, and entropy, ΔS ⁰ _m, upon micelle formation as a function of temperature were evaluated taking βvalues into calculation. Received: 28 August 1996 Accepted: 5 November 1996     

Cucurbit[n]urils (n=7, 8) binding of camptothecin and the effects on solubility and reactivity of the anticancer drug

The interaction between cucurbit[n]uril (n = 7, 8)(Q[n]) with two forms namely lactone modality and carboxylate modality of anticancer drug camptothecin (CPT) was studied. The results revealed that the combination of Q[n] with the lactone form of CPT was observed by electronic absorption spectroscopy, fluorescence spectroscopy and ¹H NMR technique in the acid solution (pH 2) and the total stability constants β were also obtained by Job plot with a host:guest ratio of 2:1; while in the phosphate buffer solution (pH 7.4), only Q[8] bound the carboxylate form of CPT in ratio 1:1, but no obvious interaction between Q[7] and the carboxylate form of CPT was observed. The solubility of CPT was enhanced up to about 70 and 8 times at pH 2 due to the formation of interaction complexes with Q[7] and Q[8], respectively, by using phase solubility method. The cytotoxicity tests revealed that compared with the free CPT, the complexes of Q[n] and CPT had the same cytotoxic activity on the human lung cancer cell line A549 and murine macrophage cell line P388D1 and the moderate depressed activity on human leukaemia cell line K562.     

Theoretical study on the adduct of chlorine trifluoride oxide and boron trifluoride—molecular and crystal structures,vibrational spectrum,and thermodynamic properties

To search for the most stable form and understand the relationship between the structure and property of the adduct of chlorine trifluoride oxide and boron trifluoride, various isomeric structures were constructed and optimized with the first principle methods. The most stable structure of the adduct has an ionic form. The vibrational spectra were calculated and assigned, and used to compute the standard thermodynamic properties. Molecular mechanics method was employed to predict the possible crystal packing of the adduct, and density functional theory method was used to refine the crystal structure and to calculate the energy band. The obtained band gap (3.40 eV) shows that the crystal is insensitive and stable under closed environment as was found experimentally. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Discovery of potential inhibitors targeting the kinase domain of polynucleotide kinase/phosphatase (PNKP): Homology modeling,virtual screening based on multiple conformations,and molecular dynamics simulation

In recent years, the level of interest has been increased in developing the DNA-repair inhibitors, to enhance the cytotoxic effects in the treatment of cancers. Polynucleotide kinase/phosphatase (PNKP) is a critical human DNA repair enzyme that repairs DNA strand breaks by catalyzing the restoration of 5’-phosphate and 3’-hydroxyl termini that are required for subsequent processing by DNA ligases and polymerases. PNKP is the only protein that repairs the 3′-hydroxyl group and 5′-phosphate group, which depicts PNKP as a potential therapeutic target. Besides, PNKP is the only DNA-repair enzyme that contains the 5′-kinase activity, therefore, targeting this kinase domain would motivate the development of novel PNKP-specific inhibitors. However, there are neither crystal structures of human PNKP nor the kinase inhibitors reported so far. Thus, in this present study, a sequential molecular docking-based virtual screening with multiple PNKP conformations integrating homology modeling, molecular dynamics simulation, and binding free energy calculation was developed to discover novel PNKP kinase inhibitors, and the top-scored molecule was finally submitted to molecular dynamics simulation to reveal the binding mechanism between the inhibitor and PNKP. Taken together, the current study could provide some guidance for the molecular docking based-virtual screening of novel PNKP kinase inhibitors.     

Investigation on the isoform selectivity of novel kinesin-like protein 1 (KIF11) inhibitor using chemical feature based pharmacophore,molecular docking,and quantum mechanical studies

Kinesin-like protein (KIF11) is a molecular motor protein that is essential in mitosis. Removal of KIF11 prevents centrosome migration and causes cell arrest in mitosis. KIF11 defects are linked to the disease of microcephaly, lymph edema or mental retardation. The human KIF11 protein has been actively studied for its role in mitosis and its potential as a therapeutic target for cancer treatment. Pharmacophore modeling, molecular docking and density functional theory approaches was employed to reveal the structural, chemical and electronic features essential for the development of small molecule inhibitor for KIF11. Hence we have developed chemical feature based pharmacophore models using Discovery Studio v 2.5 (DS). The best hypothesis (Hypo1) consisting of four chemical features (two hydrogen bond acceptor, one hydrophobic and one ring aromatic) has exhibited high correlation co-efficient of 0.9521, cost difference of 70.63 and low RMS value of 0.9475. This Hypo1 is cross validated by Cat Scramble method; test set and decoy set to prove its robustness, statistical significance and predictability respectively. The well validated Hypo1 was used as 3Dquery to perform virtual screening. The hits obtained from the virtual screening were subjected to various scrupulous drug-like filters such as Lipinski’s rule of five and ADMET properties. Finally, six hit compounds were identified based on the molecular interaction and its electronic properties. Our final lead compound could serve as a powerful tool for the discovery of potent inhibitor as KIF11 agonists.     

Corrosion inhibition performance of 2-ethyl phenyl-2, 5-dithiohydrazodicarbonamide on Fe (110)/Cu (111) in acidic/alkaline solutions: Synthesis,experimental, theoretical,and molecular dynamic studies

Herein, 2-ethyl phenyl-2,5-dithiohydrazodicarbonamide (2EPDCA) was synthesised and tested as a corrosion inhibitor for mild steel (MS) and copper (Cu) in 1 M HCl and 3.5% NaCl, respectively. Fourier transform infrared spectroscopy (FT-IR) and (NMR) nuclear magnetic resonance (¹H, ¹³C) were used to identify the chemical structure. Both experimental and computational approaches have been conducted to evaluate inhibitor efficiency on both metal systems. The electrochemical results showed that the 2EPDCA inhibition efficiency for MS systems was 95% at 1 × 10^?2 M, while in copper systems it was 97.5% at 1 × 10^?2 M. The Langmuir adsorption isotherm was fitted using adsorption surface coverage data, and for inhibitor in both systems, the kind of adsorption was mixed (physisorption and chemisorption). Through scanning electron microscopy (SEM), EDX, and atomic force microscopy (AFM) tests, we have confirmed the presence of the inhibitor molecules on the metal surface in both systems. Quantum chemistry simulations indicate that the superior corrosion inhibition efficacy of 2EPDCA on copper compared to mild steel surfaces is attributable to the former's greater electron donating propensity on copper. The adsorption of 2EPDCA molecules on Fe (110) and Cu (111) surfaces was further verified by molecular dynamic simulations, with the former having a greater adsorption energy. The results indicate that the corrosion inhibitor was effective even in harsh conditions, and it can be thought of as a novel corrosion inhibitor for mild steel and copper that provides good protection.     

Theoretical study of bifurcated hydrogen bonding effects on the 1J(N,H), 1hJ(N,H), 2hJ(N,N) couplings and 1H, 15N shieldings in model pyrroles

According to the density functional theory calculations, the X···H···N (X?N, O) intramolecular bifurcated (three‐centered) hydrogen bond with one hydrogen donor and two hydrogen acceptors causes a significant decrease of the ^1hJ(N,H) and ^2hJ(N,N) coupling constants across the N? H···N hydrogen bond and an increase of the ¹J(N,H) coupling constant across the N? H covalent bond in the 2,5‐disubsituted pyrroles. This occurs due to a weakening of the N? H···N hydrogen bridge resulting in a lengthening of the N···H distance and a decrease of the hydrogen bond angle at the bifurcated hydrogen bond formation. The gauge‐independent atomic orbital calculations of the shielding constants suggest that a weakening of the N? H···N hydrogen bridge in case of the three‐centered hydrogen bond yields a shielding of the bridge proton and deshielding of the acceptor nitrogen atom. The atoms‐in‐molecules analysis shows that an attenuation of the ^1hJ(N,H) and ^2hJ(N,N) couplings in the compounds with bifurcated hydrogen bond is connected with a decrease of the electron density ρ_H···N at the hydrogen bond critical point and Laplacian of this electron density ?²ρ_H···N. The natural bond orbital analysis suggests that the additional N? H···X interaction partly inhibits the charge transfer from the nitrogen lone pair to the σ*_{N? H} antibonding orbital across hydrogen bond weakening of the ^1hJ(N,H) and ^2hJ(N,N) trans‐hydrogen bond couplings through Fermi‐contact mechanism. An increase of the nitrogen s‐character percentage of the N? H bond in consequence of the bifurcated hydrogen bonding leads to an increase of the ¹J(N,H) coupling constant across the N? H covalent bond and deshielding of the hydrogen donor nitrogen atom. Copyright © 2010 John Wiley & Sons, Ltd.