The Hydrophobic Stabilization of Pseudomonas aeruginosa Bacteriophage F8 and the Influence of Modified Bacteriophage Preparation on Biofilm Degradation.

The bacteriophage F8 belongs to the Myoviridae group of phages and is a pathogen of Pseudomonas aeruginosa. Since Pseudomonas aeruginosa is a multidrug-resistant opportunistic bacterium and can cause serious challenges for health services, studying the potential use of phages against them is a promising approach. Pseudomonas aeruginosa can be found on medical devices because bacteria can attach to surfaces and develop biofilms, which are difficult to eradicate because of their high resistance to environmental conditions and antimicrobial therapeutics. Phage therapy is becoming promising as an alternative for the treatment of antibiotic-resistant infections, but there is still a lack of standardized protocols approved by health organizations for possible use in the clinic. In our research, we focused on the potential use of 1-octanol, which was previously used by our team to develop a method for phage purification from bacterial lysate. 1-octanol is a fatty alcohol that is mostly used in the cosmetics industry, and its advantage is that it is approved by the FDA as a food additive. In this paper, we studied the protective properties of the addition of 1-octanol for storing phage liquid preparations. We demonstrated the stabilization effect of 1-octanol addition on F8 bacteriophage preparation during storage under various conditions. Interestingly, more effective biofilm reduction was observed after treatment with the purified bacteriophage and with 1-octanol addition compared to crude lysate.

Inside out Approach to Rotator State in Hydrogen-Bonded System-Experimental and Theoretical Cross-Examination in n-Octanol.

The experimental and theoretical description of premelting behavior is one of the most challenging tasks in contemporary material science. In this paper, n-octanol was studied using a multi-method approach to investigate it at macroscopic and molecular levels. The experimental infrared (IR) spectra were collected in the solid state and liquid phase at temperature range from -84∘C to -15 ∘C to detect temperature-related indicators of pretransitional phenomena. Next, the nonlinear dielectric effect (NDE) was measured at various temperatures (from -30 ∘C to -15 ∘C) to provide insight into macroscopic effects of premelting. As a result, a two-step mechanism of premelting in n-octanol was established based on experimental data. It was postulated that it consists of a rotator state formation followed by the surface premelting. In order to shed light onto molecular-level processes, classical molecular dynamics (MD) was performed to investigate the time evolution of the changes in metric parameters as a function of simulation temperature. The applied protocol enabled simulations in the solid state as well as in the liquid (the collapse of the ordered crystal structure). The exact molecular motions contributing to the rotator state formation were obtained, revealing an enabling of the rotational freedom of the terminal parts of the chains. The Car-Parrinello molecular dynamics (CPMD) was applied to support and interpret experimental spectroscopic findings. The vibrational properties of the stretching of OH within the intermolecular hydrogen bond were studied using Fourier transformation of the autocorrelation function of both dipole moments and atomic velocity. Finally, path integral molecular dynamics (PIMD) was carried out to analyze the quantum effect's influence on the bridged proton position in the hydrogen bridge. On the basis of the combined experimental and theoretical conclusions, a novel mechanism of the bridged protons dynamics has been postulated-the interlamellar hydrogen bonding pattern, resulting in an additional OH stretching band, visible in the solid-state experimental IR spectra.

Experimental Examination of Solubility and Lipophilicity as Pharmaceutically Relevant Points of Novel Bioactive Hybrid Compounds.

The important physicochemical properties of three novel bioactive hybrid compounds with different groups (-CH3, -F and -Cl) were studied, including kinetic and thermodynamic solubility in pharmaceutically relevant solvents (buffer solutions and 1-octanol) as well as partition coefficient in system 1-octanol/buffer pH 7.4. The aqueous solubility of these chemicals is poor and ranged from 0.67 × 10-4 to 1.98 × 10-3 mol·L-1. The compounds studied are more soluble in the buffer pH 2.0, simulating the gastrointestinal tract environment (by an order of magnitude) than in the buffer pH 7.4 modelling plasma of blood. The solubility in 1-octanol is significantly higher; that is because of the specific interactions of the compounds with the solvent. The prediction solubility behaviour of the hybrid compounds using Hansen's three-parameter approach showed acceptable results. The experimental solubility of potential drugs was successfully correlated by means of two commonly known equations: modified Apelblat and van't Hoff. The temperature dependencies of partition coefficients of new hybrids in the model system 1-octanol/buffer pH 7.4 as a surrogate lipophilicity were measured by the shake flask method. It was found that compounds demonstrated a lipophilic nature and have optimal values of partition coefficients for oral absorption. Bioactive assay manifested that prepared compounds showed antifungal activities equal to or greater than fluconazole. In addition, the thermodynamic aspects of dissolution and partition processes have been examined. Bioactive assay manifested that prepared compounds showed antifungal activities equal to or greater than the reference drug.

Predicting octanol/water partition coefficients using molecular simulation for the SAMPL7 challenge: comparing the use of neat and water saturated 1-octanol.

Blind predictions of octanol/water partition coefficients at 298.15 K for 22 drug-like compounds were made for the SAMPL7 challenge. The octanol/water partition coefficients were predicted using solvation free energies computed using molecular dynamics simulations, wherein we considered the use of both pure and water-saturated 1-octanol to model the octanol-rich phase. Water and 1-octanol were modeled using TIP4P and TrAPPE-UA, respectively, which have been shown to well reproduce the experimental mutual solubility, and the solutes were modeled using GAFF. After the close of the SAMPL7 challenge, we additionally made predictions using TIP4P/2005 water. We found that the predictions were sensitive to the choice of water force field. However, the effect of water in the octanol-rich phase was found to be even more significant and non-negligible. The effect of inclusion of water was additionally sensitive to the chemical structure of the solute.

Multiple linear regression models for predicting the n­octanol/water partition coefficients in the SAMPL7 blind challenge.

A multiple linear regression model called MLR-3 is used for predicting the experimental n-octanol/water partition coefficient (log PN) of 22 N-sulfonamides proposed by the organizers of the SAMPL7 blind challenge. The MLR-3 method was trained with 82 molecules including drug-like sulfonamides and small organic molecules, which resembled the main functional groups present in the challenge dataset. Our model, submitted as "TFE-MLR", presented a root-mean-square error of 0.58 and mean absolute error of 0.41 in log P units, accomplishing the highest accuracy, among empirical methods and also in all submissions based on the ranked ones. Overall, the results support the appropriateness of multiple linear regression approach MLR-3 for computing the n-octanol/water partition coefficient in sulfonamide-bearing compounds. In this context, the outstanding performance of empirical methodologies, where 75% of the ranked submissions achieved root-mean-square errors < 1 log P units, support the suitability of these strategies for obtaining accurate and fast predictions of physicochemical properties as partition coefficients of bioorganic compounds.

SAMPL7 physical property prediction from EC-RISM theory.

Inspired by the successful application of the embedded cluster reference interaction site model (EC-RISM), a combination of quantum-mechanical calculations with three-dimensional RISM theory to predict Gibbs energies of species in solution within the SAMPL6.1 (acidity constants, pKa) and SAMPL6.2 (octanol-water partition coefficients, log P) the methodology was applied to the recent SAMPL7 physical property challenge on aqueous pKa and octanol-water log P values. Not part of the challenge but provided by the organizers, we also computed distribution coefficients log D7.4 from predicted pKa and log P data. While macroscopic pKa predictions compared very favorably with experimental data (root mean square error, RMSE 0.72 pK units), the performance of the log P model (RMSE 1.84) fell behind expectations from the SAMPL6.2 challenge, leading to reasonable log D7.4 predictions (RMSE 1.69) from combining the independent calculations. In the post-submission phase, conformations generated by different methodology yielded results that did not significantly improve the original predictions. While overall satisfactory compared to previous log D challenges, the predicted data suggest that further effort is needed for optimizing the robustness of the partition coefficient model within EC-RISM calculations and for shaping the agreement between experimental conditions and the corresponding model description.

COSMO-RS predictions of logP in the SAMPL7 blind challenge.

We applied the COSMO-RS method to predict the partition coefficient logP between water and 1-octanol for 22 small drug like molecules within the framework of the SAMPL7 blind challenge. We carefully collected a set of thermodynamically meaningful microstates, including tautomeric forms of the neutral species, and calculated the logP using the current COSMOtherm implementation on the most accurate level. With this approach, COSMO-RS was ranked as the 6st most accurate method (Measured by the mean absolute error (MAE) of 0.57) over all 17 ranked submissions. We achieved a root mean square deviation (RMSD) of 0.78. The largest deviations from experimental values are exhibited by five SAMPL molecules (SM), which seem to be shifted in most SAMPL7 contributions. In context with previous SAMPL challenges, COSMO-RS demonstrates a wide range of applicability and one of the best in class reliability and accuracy among the physical methods.

SAMPL7 blind challenge: quantum-mechanical prediction of partition coefficients and acid dissociation constants for small drug-like molecules.

The physicochemical properties of a drug molecule determine the therapeutic effectiveness of the drug. Thus, the development of fast and accurate theoretical approaches for the prediction of such properties is inevitable. The participation to the SAMPL7 challenge is based on the estimation of logP coefficients and pKa values of small drug-like sulfonamide derivatives. Thereby, quantum mechanical calculations were carried out in order to calculate the free energy of solvation and the transfer energy of 22 drug-like compounds in different environments (water and n-octanol) by employing the SMD solvation model. For logP calculations, we studied eleven different methodologies to calculate the transfer free energies, the lowest RMSE value was obtained for the M06L/def2-TZVP//M06L/def2-SVP level of theory. On the other hand, we employed an isodesmic reaction scheme within the macro pKa framework; this was based on selecting reference molecules similar to the SAMPL7 challenge molecules. Consequently, highly well correlated pKa values were obtained with the M062X/6-311+G(2df,2p)//M052X/6-31+G(d,p) level of theory.

Energy-entropy prediction of octanol-water logP of SAMPL7 N-acyl sulfonamide bioisosters.

Partition coefficients quantify a molecule's distribution between two immiscible liquid phases. While there are many methods to compute them, there is not yet a method based on the free energy of each system in terms of energy and entropy, where entropy depends on the probability distribution of all quantum states of the system. Here we test a method in this class called Energy Entropy Multiscale Cell Correlation (EE-MCC) for the calculation of octanol-water logP values for 22 N-acyl sulfonamides in the SAMPL7 Physical Properties Challenge (Statistical Assessment of the Modelling of Proteins and Ligands). EE-MCC logP values have a mean error of 1.8 logP units versus experiment and a standard error of the mean of 1.0 logP units for three separate calculations. These errors are primarily due to getting sufficiently converged energies to give accurate differences of large numbers, particularly for the large-molecule solvent octanol. However, this is also an issue for entropy, and approximations in the force field and MCC theory also contribute to the error. Unique to MCC is that it explains the entropy contributions over all the degrees of freedom of all molecules in the system. A gain in orientational entropy of water is the main favourable entropic contribution, supported by small gains in solute vibrational and orientational entropy but offset by unfavourable changes in the orientational entropy of octanol, the vibrational entropy of both solvents, and the positional and conformational entropy of the solute.

An easy-to-use and general nuclear magnetic resonance method for the determination of partition coefficients of drugs and natural products.

The lipophilicity of a drug is an important parameter for its eventual development by the pharmaceutical industry. It is usually measured by HPLC following partition of the compound between water and 1-octanol. We present here an alternative, simple, sensitive and quantitative 1 H nuclear magnetic resonance (NMR) method for the experimental measurement of partition coefficients of natural compounds and pharmaceutical drugs. It is based on measuring concentrations in the water phase, before and after partitioning and equilibration between water and octanol, using the ERETIC (Electronic Reference To Access In Vivo Concentration) technique. The signal to noise ratio is improved by a Water Suppression by Excitation Sculpting sequence. Quantification is based on an electronic reference signal and does not need addition of a reference compound. The log P values of 22 natural metabolites and four pharmaceutical drugs were determined and the experimental results are in excellent agreement with literature data. The experiments were run on ~2 mg material. This technique proved to be robust, reproducible and suitable for log P values between -2 and +2.

Isavuconazole: Thermodynamic Evaluation of Processes Sublimation, Dissolution and Partition in Pharmaceutically Relevant Media.

A temperature dependence of saturated vapor pressure of isavuconazole (IVZ), an antimycotic drug, was found by using the method of inert gas-carrier transfer and the thermodynamic functions of sublimation were calculated at a temperature of 298.15 K. The value of the compound standard molar enthalpy of sublimation was found to be 138.1 ± 0.5 kJ·mol-1. The IVZ thermophysical properties-melting point and enthalpy-equaled 302.7 K and 29.9 kJ mol-1, respectively. The isothermal saturation method was used to determine the drug solubility in seven pharmaceutically relevant solvents within the temperature range from 293.15 to 313.15 K. The IVZ solubility in the studied solvents increased in the following order: buffer pH 7.4, buffer pH 2.0, buffer pH 1.2, hexane, 1-octanol, 1-propanol, ethanol. Depending on the solvent chemical nature, the compound solubility varied from 6.7 × 10-6 to 0.3 mol·L-1. The Hansen s approach was used for evaluating and analyzing the solubility data of drug. The results show that this model well-described intermolecular interactions in the solutions studied. It was established that in comparison with the van't Hoff model, the modified Apelblat one ensured the best correlation with the experimental solubility data of the studied drug. The activity coefficients at infinite dilution and dissolution excess thermodynamic functions of IVZ were calculated in each of the solvents. Temperature dependences of the compound partition coefficients were obtained in a binary 1-octanol/buffer pH 7.4 system and the transfer thermodynamic functions were calculated. The drug distribution from the aqueous solution to the organic medium was found to be spontaneous and entropy-driven.

A blind SAMPL6 challenge: insight into the octanol-water partition coefficients of drug-like molecules via a DFT approach.

In this study quantum mechanical methods were used to predict the solvation energies of a series of drug-like molecules both in water and in octanol, in the context of the SAMPL6 n-octanol/water partition coefficient challenge. In pharmaceutical design, n-octanol/water partition coefficient, LogP, describes the drug's hydrophobicity and membrane permeability, thus, a well-established theoretical method that rapidly determines the hydrophobicity of a drug, enables the progress of the drug design. In this study, the solvation free energies were obtained via six different methodologies (B3LYP, M06-2X and ωB97XD functionals with 6-311+G** and 6-31G* basis sets) by taking into account the environment implicitly; the methodology chosen (B3LYP/6-311+G**) was used later to evaluate ΔGsolv by using explicit water as solvent. We optimized each conformer in different solvents separately, our calculations have shown that the stability of the conformers is highly dependent on the solvent environment. We have compared implicitly and explicitly solvated systems, the interaction of one explicit water with drug-molecules at the proper location leads to the prediction of more accurate LogP values.

Predicting partition coefficients of drug-like molecules in the SAMPL6 challenge with Drude polarizable force fields.

The water-octanol partition coefficient is an important physicochemical property for small molecule drug design. Here, we report our participation in the SAMPL6 logP prediction challenge with free energy perturbation (FEP) calculations in the water phase and in the 1-octanol phase using Drude polarizable force fields. Root mean square error (RMSE) and mean absolute error (MAE) of our prediction are equal to 1.85 and 1.25 logP units. The errors are not evenly distributed. Out of eleven SAMPL6 solutes, FEP/Drude performed very badly on three molecules (deviations all larger than 2 logP units) but good on the remaining eight (deviations all less than 1 logP unit). We find while FEP converges well within one nanosecond in water, simulations in 1-octanol need much longer simulation time and possibly more independent runs for sampling. We also find out that 1-octanol, albeit being a non-polar solvent, still polarizes solute molecules and forms stable hydrogen bonds with them. At the end, we attempt to reweight FEP trajectories with QM/Drude calculations and discuss possible caveats in our simulation setup.

The SAMPL6 challenge on predicting octanol-water partition coefficients from EC-RISM theory.

Results are reported for octanol-water partition coefficients (log P) of the neutral states of drug-like molecules provided during the SAMPL6 (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenge from applying the "embedded cluster reference interaction site model" (EC-RISM) as a solvation model for quantum-chemical calculations. Following the strategy outlined during earlier SAMPL challenges we first train 1- and 2-parameter water-free ("dry") and water-saturated ("wet") models for n-octanol solvation Gibbs energies with respect to experimental values from the "Minnesota Solvation Database" (MNSOL), yielding a root mean square error (RMSE) of 1.5 kcal mol-1 for the best-performing 2-parameter wet model, while the optimal water model developed for the pKa part of the SAMPL6 challenge is kept unchanged (RMSE 1.6 kcal mol-1 for neutral compounds from a model trained on both neutral and ionic species). Applying these models to the blind prediction set yields a log P RMSE of less than 0.5 for our best model (2-parameters, wet). Further analysis of our results reveals that a single compound is responsible for most of the error, SM15, without which the RMSE drops to 0.2. Since this is the only compound in the challenge dataset with a hydroxyl group we investigate other alcohols for which Gibbs energy of solvation data for both water and n-octanol are available in the MNSOL database to demonstrate a systematic cause of error and to discuss strategies for improvement.

Affinity of Antifungal Isoxazolo[3,4-b]pyridine-3(1H)-Ones to Phospholipids in Immobilized Artificial Membrane (IAM) Chromatography.

Currently, rapid evaluation of the physicochemical parameters of drug candidates, such as lipophilicity, is in high demand owing to it enabling the approximation of the processes of absorption, distribution, metabolism, and elimination. Although the lipophilicity of drug candidates is determined using the shake flash method (n-octanol/water system) or reversed phase liquid chromatography (RP-LC), more biosimilar alternatives to classical lipophilicity measurement are currently available. One of the alternatives is immobilized artificial membrane (IAM) chromatography. The present study is a continuation of our research focused on physiochemical characterization of biologically active derivatives of isoxazolo[3,4-b]pyridine-3(1H)-ones. The main goal of this study was to assess the affinity of isoxazolones to phospholipids using IAM chromatography and compare it with the lipophilicity parameters established by reversed phase chromatography. Quantitative structure-retention relationship (QSRR) modeling of IAM retention using differential evolution coupled with partial least squares (DE-PLS) regression was performed. The results indicate that in the studied group of structurally related isoxazolone derivatives, discrepancies occur between the retention under IAM and RP-LC conditions. Although some correlation between these two chromatographic methods can be found, lipophilicity does not fully explain the affinities of the investigated molecules to phospholipids. QSRR analysis also shows common factors that contribute to retention under IAM and RP-LC conditions. In this context, the significant influences of WHIM and GETAWAY descriptors in all the obtained models should be highlighted.

Quantification of controlled release leuprolide and triptorelin in rabbit plasma using electromembrane extraction coupled with HPLC-UV.

An electromembrane extraction followed by HPLC-UV technique was developed and validated for quantification of leuprolide and triptorelin in rabbit plasma. The influencing parameters on the extraction efficiency were optimized using experimental design methodology. The optimized conditions were found to be; supported liquid membrane: a mixture of 1-octanol and 2-ethyl hexanol (1:1) containing 10% v/v di(2-ethylhexyl) phosphate, applied voltage: 5 V, extraction time: 5 min, pH of the donor phase: 4.5 and pH of the acceptor phase: 1.0. The optimized method was validated for linearity, intraday and interday precision, and accuracy in rabbit plasma. The range of quantification for both peptides was 0.5-1000 ng/mL with regression coefficients higher than 0.994. Relative recoveries of leuprolide and triptorelin were found to be 80.3 and 75.5%, respectively. Limits of quantification and detection for both peptides were found to be 0.5 and 0.15 ng/mL, respectively. The validated method was successfully applied to pharmacokinetic study of the 1-month depot formulations of each peptide after subcutaneous administration to rabbits.

Uncovering abnormal changes in logP after fluorination using molecular dynamics simulations.

The fluorination-induced changes in the logP (1-octanol/water partition coefficient) of ligands were examined by molecular dynamics simulations. The protocol and force field parameters were first evaluated by calculating the logP values for n-alkanes, and their monofluorinated and monochlorinated analogs. Then, the logP values of several test sets (1-butanol, 3-propyl-1H-indole, and analogs fluorinated at the terminal methyl group) were calculated. The calculated results agree well with experiment, and the root mean square error values are 0.61 and 0.68 log units for the GAFF and GAFF2 force fields, respectively. Finally, the logP estimation was extended to a drug molecule, TAK-438, for which fluorination-induced abnormal logP reduction has been observed experimentally. This abnormal change was qualitatively reproduced by the molecular dynamics simulations. We found that the abnormal logP reduction can be mainly attributed to the effect of fluorination-induced dipole change. Our results suggest that molecular simulation is a useful strategy to predict the fluorination-induced change in logP for drug discovery applications.

Electromembrane extraction of verapamil and riluzole from urine and wastewater samples using a mixture of organic solvents as a supported liquid membrane: Study on electric current variations.

In this study, the application of a mixture of organic solvents as a supported liquid membrane for improving the efficiency of the electromembrane extraction procedure was investigated. The extraction process was followed by high-performance liquid chromatography analysis of two model drugs (verapamil and riluzole). In this research, four organic solvents, including 1-heptanol, 1-octanol, 2-nitrophenyl octyl ether, and 2-ethyl hexanol, were selected as model solvents and different binary mixtures (v/v 2:1, 1:1 and 1:2) were used as the supported liquid membrane. The mixture of 2-ethyl hexanol and 1-otanol (v/v, 2:1) improved the extraction efficiency of model drugs by 1.5 to 12 times. It was found that extraction efficiency is greatly influenced by the level of electric current. In this study, for various mixtures of organic solvents, the electric current fluctuated between 50 and 2500 µA, and the highest extraction efficiencies were obtained with low and stable electric currents. Finally, the optimized extraction condition was validated and applied for the determination of model drugs in urine and wastewater samples.

QSAR characterization of new synthesized hydantoins with antiproliferative activity.

Hydantois have been identified as constituents of a number of pharmacologically active molecules. In the present study, we have examined in vitro antiproliferative activity against human colon cancer cell lines HCT-116 of three series of 3-(4-substituted benzyl)-hydantoins with various substituent attached in position 5 of the hydantoin ring. Since the investigated compounds have recently been synthesized and show antiproliferative activity, a good understanding of the properties of the potential drug responsible for their pharmacokinetics is an important goal for their further development. One of the important properties is lipophilicity. Lipophilicity has been assessed by reversed-phase liquid chromatography (high-performance thin-layer chromatography and high-pressure liquid chromatography) by means of direct and indirect (using calibration curve) methods. Chromatographic lipophilicity indices in addition to calculated logP values were compared by hierarchical cluster analysis. The linear solvation energy relationship approach was used to understand and compare the types and relative strength of the molecular interactions that occur in the chromatographic as well as in the n-octanol-water partitioning systems. Finally, correlation between in silico pharmacokinetic predictors and antiproliferative activity was examined. Preliminary quantitative structure-activity relationship modeling indicates that pharmacokinetic predictors capture only one-quarter of all chemical features that are important for antiproliferative activity itself. Among selected descriptors are chromatographic lipophilicity indices.

Acute and subacute effects of drugs in embryos of Danio rerio. QSAR grouping and modelling.

The final fate of many drugs is release into the natural aquatic environment. It is necessary to assess the toxicity caused by this situation and the associated concerns for human beings. Zebrafish (Danio rerio) is a common biomodel used to assess toxicity in aquatic environments. The zebrafish embryo toxicity test was selected to evaluate the acute toxicity of several drugs (diphenhydramine, gentamicin, tobramycin, enalapril and lidocaine) due to the lack of such information. Lethal and sublethal effects were detected, and the LC50 values of the drugs ranged from 11.0 mg/L to 422·102 mg/L. For all of the drugs tested, these values were higher than the concentrations found in the natural environment. Therefore, there was a low environmental toxicological risk. Nevertheless, teratogenic effects were also recorded when embryos of zebrafish were exposed to caffeine (control drug), diphenhydramine and lidocaine at lower concentrations than the respective LC50 values. Quantitative structure-activity relationship analysis was also performed to analyse these drugs and other chemicals with pharmaceutical uses as well as previous toxicological data in this vertebrate after 48 h of exposure. It is estimated that the partition coefficient, log P, is the main physicochemical property related to the ecotoxicological data and can be used for the development of a mathematical model.