Antimicrobial Activity Classification of Imidazolium Derivatives Predicted by Artificial Neural Networks.

PURPOSE: This study assesses the Multilayer Perceptron (MLP) neural network, complemented by other Machine Learning techniques (CART, PCA), in predicting the antimicrobial activity of 140 newly designed imidazolium chlorides against Klebsiella pneumoniae before synthesis. Emphasis is on leveraging molecular properties for predictive analysis. METHODS: Classification and regression decision trees (CART) identified the top 200 predictive molecular descriptors. Principal Component Analysis (PCA) reduced these descriptors to 5 components, retaining 99.57% of raw data information. Antimicrobial activity, categorized as high or low, was based on experimentally proven minimal inhibitory concentration (MIC), with a cut-point at MIC = 0.856 mol/L. A 12-fold cross-validation trained the MLP (architecture 5-12-2 with 5 Principal Components). RESULTS: The MLP exhibited commendable performance, achieving almost 90% correct classifications across learning, validation, and test sets, outperforming models without PCA dimension reduction. Key metrics, including accuracy (0.907), sensitivity (0.905), specificity (0.909), and precision (0.891), were notably high. These results highlight the MLP model's efficacy with PCA as a high-quality classifier for determining antimicrobial activity. CONCLUSIONS: The study concludes that the MLP neural network, along with CART and PCA, is a robust tool for predicting the antimicrobial activity class of imidazolium chlorides against Klebsiella pneumoniae. CART and PCA, used in this study, allowed input variable reduction without significant information loss. High classification accuracy and associated metrics affirm the method's potential utility in pre-synthesis assessments, offering valuable insights for antimicrobial compound design.

Toward the General Mechanistic Model of Liquid Chromatographic Retention.

Large datasets of chromatographic retention times are relatively easy to collect. This statement is particularly true when mixtures of compounds are analyzed under a series of gradient conditions using chromatographic techniques coupled with mass spectrometry detection. Such datasets carry much information about chromatographic retention that, if extracted, can provide useful predictive information. In this work, we proposed a mechanistic model that jointly explains the relationship between pH, organic modifier type, temperature, gradient duration, and analyte retention based on liquid chromatography retention data collected for 187 small molecules. The model was built utilizing a Bayesian multilevel framework. The model assumes (i) a deterministic Neue equation that describes the relationship between retention time and analyte-specific and instrument-specific parameters, (ii) the relationship between analyte-specific descriptors (log P, pKa, and functional groups) and analyte-specific chromatographic parameters, and (iii) stochastic components of between-analyte and residual variability. The model utilizes prior knowledge about model parameters to regularize predictions which is important as there is ample information about the retention behavior of analytes in various stationary phases in the literature. The usefulness of the proposed model in providing interpretable summaries of complex data and in decision making is discussed.

Statistical analysis of isocratic chromatographic data using Bayesian modeling.

Chromatographic retention times are usually modeled considering only one analyte at a time. However, it has certain limitations as no information is shared between the analytes, and consequently the model predictions poorly generalize to out-of-sample analytes. In this work, a publicly available dataset was used to illustrate the benefits of pooling the individual data and analyzing them simultaneously utilizing Bayesian hierarchical approach. Statistical analysis was carried out using the Stan program coupled with R, which enables full Bayesian inference with Markov chain Monte Carlo sampling. This methodology allows (i) incorporating prior knowledge about the likely values of model parameters, (ii) considering the between-analyte variability and the correlation between the model parameters, (iii) explaining the between-analyte variability by available predictors, and (iv) sharing information across the analytes. The latter is especially valuable when only limited information is available in the data about certain model parameters. The results are obtained in the form of posterior probability distribution, which quantifies uncertainty about the model parameters and predictions. Posterior probability is also directly relevant for decision-making. In this work, we used the Neue model to describe the relationship between retention factor and acetonitrile content in the mobile phase for 1026 analytes. The model was parametrized in terms of retention factor in 100% water, retention factor in 100% acetonitrile, and curvature coefficient, and considered log P and pKa as predictors. From this analysis, we discovered that the analytes formed two clusters with different retention depending on the degree of analyte dissociation. The final model turned out to be well calibrated with the data. It gives insight into the behavior of analytes in the chromatographic column and can be used to make predictions for a structurally diverse set of analytes if their log P and pKa values are known.

Application of Bayesian Multilevel Modeling in the Quantitative Structure-Retention Relationship Studies of Heterogeneous Compounds.

Quantitative structure-retention relationships (QSRRs) are used in the field of chromatography to model the relationship between an analyte structure and chromatographic retention. Such models are typically difficult to build and validate for heterogeneous compounds because of their many descriptors and relatively limited analyte-specific data. In this study, a Bayesian multilevel model is proposed to characterize the isocratic retention time data collected for 1026 heterogeneous analytes. The QSRR considers the effects of the molecular mass and 100 functional groups (substituents) on analyte-specific chromatographic parameters of the Neue model (i.e., the retention factor in water, the retention factor in acetonitrile, and the curvature coefficient). A Bayesian multilevel regression model was used to smooth noisy parameter estimates with too few data and to consider the uncertainties in the model parameters. We discuss the benefits of the Bayesian multilevel model (i) to understand chromatographic data, (ii) to quantify the effect of functional groups on chromatographic retention, and (iii) to predict analyte retention based on various types of preliminary data. The uncertainty of isocratic and gradient predictions was visualized using uncertainty chromatograms and discussed in terms of usefulness in decision making. We think that this method will provide the most benefit in providing a unified scheme for analyzing large chromatographic databases and assessing the impact of functional groups and other descriptors on analyte retention.

Analysis of Isocratic-Chromatographic-Retention Data using Bayesian Multilevel Modeling.

The objective of this work was to develop a multilevel (hierarchical) model based on isocratic-reversed-phase-high-performance-chromatographic data collected in methanol and acetonitrile for 58 chemical compounds. Such a multilevel model is a regression model of the analyte-specific chromatographic measurements, in which all the regression parameters are given a probability model. It is a fundamentally different approach from the most common approach, where parameters are separately estimated for each analyte (without sharing information across analytes and different organic modifiers). The statistical analysis was done with Stan software implementing the Bayesian-statistics inference with Markov-chain Monte Carlo sampling. During the model-building process, a series of multilevel models of different complexity were obtained, such as (1) a model with no pooling (separate models were fitted for each analyte), (2) a model with partial pooling (a common distribution was used for analyte-specific parameters), and (3) a model with partial pooling as well as a regression model relating analyte-specific parameters and analyte-specific properties (QSRR equations). All the models were compared with each other using 10-fold cross-validation. The benefits of multilevel models in inference and predictions were shown. In particular the obtained models allowed us to (i) better understand the data and (ii) solve many routine analytical problems, such as obtaining well-calibrated predictions of retention factors for an analyte in acetonitrile-containing mobile phases given zero, one, or several measurements in methanol-containing mobile phases and vice versa.

pH-Gradient Liquid Chromatography: Fundamentals and Examples.

In this paper, we acquaint the readers with the fundamentals of gradient separation, followed by the latest innovations in this field. We describe the principles of organic modifier- and pH-gradient elution emphasizing the differences and similarities with isocratic separation. The double organic modifier-/pH-gradient is also thoroughly reviewed as a useful method for the simultaneous determination of logkw (substitute of logP) and the pka of analytes present in complex mixtures.

Comparative pharmacodynamic analysis of imidazoline compounds using rat model of ocular mydriasis with a test of quantitative structure-activity relationships.

Imidazol(in)e derivatives, having the chemical structure similar to clonidine, exert diverse pharmacological activities connected with their interactions with alpha2-adrenergic receptors, e.g. hypotension, bradycardia, sedation as well as antinociceptive, anxiolytic, antiarrhythmic, muscle relaxant and mydriatic effects. The mechanism of pupillary dilation observed after systemic administration of imidazol(in)es to rats, mice and cats depends on the stimulation of postsynaptic alpha2-adrenoceptors within the brain. It was proved that the central nervous system (CNS)-localized I1-imidazoline receptors are not engaged in those effects. It appeared interesting to analyze the CNS-mediated pharmacodynamics of imidazole(in)e agents in terms of their chromatographic and calculation chemistry-derived parameters. In the present study a systematic determination and comparative pharmacometric analysis of mydriatic effects in rats were performed on a series of 20 imidazol(in)e agents, composed of the well-known drugs and of the substances used in experimental pharmacology. The eye pupil dilatory activities of the compounds were assessed in anesthetized Wistar rats according to the established Koss method. Among twenty imidazol(in)e derivatives studied, 18 produced diverse dose-dependent mydriatic effects. In the quantitative structure-activity relationships (QSAR) analysis, the pharmacological data (half maximum mydriatic effect - ED50 in µmol/kg) were considered along with the structural parameters of the agents from molecular modeling. The theoretically calculated lipophilicity parameters, CLOGP, of imidazol(in)es, as well as their lipophilicity parameters from HPLC, logkw, were also considered. The attempts to derive statistically significant QSAR equations for a full series of the agents under study were unsuccessful. However, for a subgroup of eight apparently structurally related imidazol(in)es a significant relationship between log(1/ED50) and logkw values was obtained. The lack of "predictive" QSAR for the whole series of the structurally diverse agents is probably due to a complex mechanism of the ligand-alpha2-adrenergic receptor interactions, which are predominantly of a highly structurally specific polar nature. Such interactions are difficult to quantify with the established chemical structural descriptors, contrary to the less specific, molecular bulkiness-related interactions.

Quantitative structure-(chromatographic) retention relationship models for dissociating compounds.

The aim of this work was to develop mathematical models relating the hydrophobicity and dissociation constant of an analyte with its structure, which would be useful in predicting analyte retention times in reversed-phase liquid chromatography. For that purpose a large and diverse group of 115 drugs was used to build three QSRR models combining retention-related parameters (logkw-chromatographic measure of hydrophobicity, S-slope factor from Snyder-Soczewinski equation, and pKa) with structural descriptors calculated by means of molecular modeling for both dissociated and nondissociated forms of analytes. Lasso, Stepwise and PLS regressions were used to build statistical models. Moreover a simple QSRR equations based on lipophilicity and dissociation constant parameters calculated in the ACD/Labs software were proposed and compared with quantum chemistry-based QSRR equations. The obtained relationships were further used to predict chromatographic retention times. The predictive performances of the obtained models were assessed using 10-fold cross-validation and external validation. The QSRR equations developed were simple and were characterized by satisfactory predictive performance. Application of quantum chemistry-based and ACD-based descriptors leads to similar accuracy of retention times' prediction.

Simultaneous determination of hydrophobicity and dissociation constant for a large set of compounds by gradient reverse phase high performance liquid chromatography-mass spectrometry technique.

Fast and reliable methods for the determination of hydrophobicity and acidity are desired in pre-clinical drug development phases to eliminate compounds with poor pharmacokinetic properties. Reversed-phase high-performance liquid chromatography (RP HPLC) coupled with time-of-flight mass spectrometry (RP HPLC-ESI-TOF-MS) is a convenient technique for that purpose. In this work we determined the chromatographic measure of hydrophobicity (logkw) and dissociation constant (pKa) simultaneously for a large and diverse group of 161 drugs. Retention times were determined by means of RP HPLC-ESI-TOF-MS for a series of pH and organic modifier gradients. We were able to measure retention times for 140 out of 161 (87%) compounds. For those analytes logkw and pKa parameters were calculated and compared with literature and ACD Labs-calculated data. The determined chromatographic measure of hydrophobicity and dissociation constant was closely related to literature and theoretically calculated values. Applied methodology achieved the medium-throughput screening rate of 100 compounds per day and proved to be a simple, fast and reliable approach of assessing important physicochemical properties of drugs. This technique has certain limitations as it is not applicable for very hydrophilic analytes (logP<0.5) and compounds with identical molar masses.

Activation of Human Toll-like Receptor 4 (TLR4)·Myeloid Differentiation Factor 2 (MD-2) by Hypoacylated Lipopolysaccharide from a Clinical Isolate of Burkholderia cenocepacia.

Lung infection by Burkholderia species, in particular Burkholderia cenocepacia, accelerates tissue damage and increases post-lung transplant mortality in cystic fibrosis patients. Host-microbe interplay largely depends on interactions between pathogen-specific molecules and innate immune receptors such as Toll-like receptor 4 (TLR4), which recognizes the lipid A moiety of the bacterial lipopolysaccharide (LPS). The human TLR4·myeloid differentiation factor 2 (MD-2) LPS receptor complex is strongly activated by hexa-acylated lipid A and poorly activated by underacylated lipid A. Here, we report that B. cenocepacia LPS strongly activates human TLR4·MD-2 despite its lipid A having only five acyl chains. Furthermore, we show that aminoarabinose residues in lipid A contribute to TLR4-lipid A interactions, and experiments in a mouse model of LPS-induced endotoxic shock confirmed the proinflammatory potential of B. cenocepacia penta-acylated lipid A. Molecular modeling combined with mutagenesis of TLR4-MD-2 interactive surfaces suggests that longer acyl chains and the aminoarabinose residues in the B. cenocepacia lipid A allow exposure of the fifth acyl chain on the surface of MD-2 enabling interactions with TLR4 and its dimerization. Our results provide a molecular model for activation of the human TLR4·MD-2 complex by penta-acylated lipid A explaining the ability of hypoacylated B. cenocepacia LPS to promote proinflammatory responses associated with the severe pathogenicity of this opportunistic bacterium.

Maximum A Posteriori Bayesian Estimation of Chromatographic Parameters by Limited Number of Experiments.

The aim of this work was to develop a nonlinear mixed-effect chromatographic model able to describe the retention times of weak acids and bases in all possible combinations of organic modifier content and mobile-phase pH. Further, we aimed to identify the influence of basic covariates, like lipophilicity (log P), dissociation constant (pK(a)), and polar surface area (PSA), on the intercompound variability of chromatographic parameters. Lastly, we aimed to propose the optimal limited experimental design to the estimation process of parameters through a maximum a posteriori (MAP) Bayesian method to facilitate the method development process. The data set comprised retention times for two series of organic modifier content collected at different pH for a large series of acids and bases. The obtained typical parameters and their distribution were subsequently used as priors to improve the estimation process from reduced design with a variable number of preliminary experiments. The MAP Bayesian estimator was validated using two external-validation data sets. The common literature model was used to relate analyte retention time with mobile-phase pH and organic modifier content. A set of QSRR-based covariate relationships was established. It turned out that four preliminary experiments and prior information that includes analyte pK(a), log P, acid/base type, and PSA are sufficient to accurately predict analyte retention in virtually all combined changes of pH and organic modifier content. The MAP Bayesian estimator of all important chromatographic parameters controlling retention in pH/organic modifier gradient was developed. It can be used to improve parameter estimation using limited experimental design.

Least absolute shrinkage and selection operator and dimensionality reduction techniques in quantitative structure retention relationship modeling of retention in hydrophilic interaction liquid chromatography.

The objective of this study was to model the retention of nucleosides and pterins in hydrophilic interaction liquid chromatography (HILIC) via QSRR-based approach. Two home-made (Amino-P-C18, Amino-P-C10) and one commercial (IAM.PC.DD2) HILIC stationary phases were considered. Logarithm of retention factor at 5% of acetonitrile (logkACN) along with descriptors obtained for 16 nucleosides and 11 pterins were used to develop QSRR models. We used and compared the predictive performance of three regression techniques: partial least square (PLS), the least absolute shrinkage and selection operator (LASSO), and the LASSO followed by stepwise multiple linear regression. The highest predictive squared correlation coefficient (QLOOCV(2)) in PLS analysis was found for Amino-P-C10 (QLOOCV(2)=0.687) and IAM.PC.DD2 (QLOOCV(2)=0.506) and the lowest for IAM.PC.DD2 (QLOOCV(2)=-0.01). Much higher values were obtained for the LASSO model. The QLOOCV(2) equaled 0.9 for Amino-P-C10, 0.66 for IAM.PC.DD2 and 0.59 for Amino-P-C18. The combination of LASSO with stepwise regression provided models with comparable predictive performance as the LASSO, however with possibility of calculating the standard error of estimates. The use of LASSO itself and in combination with classical stepwise regression may offer greater stability of the developed models thanks to more smooth change of coefficients and reduced susceptibility towards chance correlation. Application of QSRR-based approach, along with the computational methods proposed in this work, may offer a useful approach in the modeling of retention of nucleoside and pterin compounds in HILIC.

Modulation of CD14 and TLR4·MD-2 activities by a synthetic lipid A mimetic.

Monosaccharide lipid A mimetics based on a glucosamine core linked to two fatty acid chains and bearing one or two phosphate groups have been synthesized. Compounds 1 and 2, each with one phosphate group, were practically inactive in inhibiting LPS-induced TLR4 signaling and cytokine production in HEK-blue cells and murine macrophages, but compound 3, with two phosphate groups, was found to be active in efficiently inhibiting TLR4 signal in both cell types. The direct interaction between compound 3 and the MD-2 coreceptor was investigated by NMR spectroscopy and molecular modeling/docking analysis. This compound also interacts directly with the CD14 receptor, stimulating its internalization by endocytosis. Experiments on macrophages show that the effect on CD14 reinforces the activity on MD-2·TLR4 because compound 3's activity is higher when CD14 is important for TLR4 signaling (i.e., at low LPS concentration). The dual targeting of MD-2 and CD14, accompanied by good solubility in water and lack of toxicity, suggests the use of monosaccharide 3 as a lead compound for the development of drugs directed against TLR4-related syndromes.

A new pH/organic modifier gradient RP HPLC method for convenient determination of lipophilicity and acidity of drugs as applied to established imidazoline agents.

Convenient drug candidates testing methods for lipophilicity and acidity are highly requested in modern pharmaceutical research and development strategy. Reversed-phase high-performance liquid chromatography (RP HPLC) might be particularly useful for the determination of both pK(a) and the apparent (pH-dependent) octanol-water partition coefficient, applicable in high-throughput analysis of multi-component mixtures. In this report the pH/organic modifier gradient RP HPLC is presented as a means of simultaneous determination of acidity and lipophilicity of a series of 26 imidazoline-like drugs. The previously theoretically elaborated approach has been applied consisting in retention measurements in a series of methanol gradient runs differing in pH range and duration of the gradient. The simultaneously determined lipophilicity and dissociation constants have been demonstrated to correlate to the respective parameters form calculation chemistry. The proposed approach can be applied to compound mixtures, it requires only minute amounts of substances, and pK(a) values can be determined in the range 3-10 units and lipophilicity log P parameter in the range 0-7 units.