The prediction of mutagenicity and pKa for pharmaceutically relevant compounds using 'quantum chemical topology' descriptors

Quantum Chemical Topology (QCT) descriptors, calculated from ab initio wave functions, have been utilised to model pKa and mutagenicity for data sets of pharmaceutically relevant compounds. The pKa of a compound is a pivotal property in both life science and chemistry since the propensity of a compound to donate or accept a proton is fundamental to understanding chemical and biological processes. The prediction of mutagenicity, specifically as determined by the Ames test, is important to aid medicinal chemists select compounds avoiding this potential pitfall in drug design. Carbocyclic and heterocyclic aromatic amines were chosen because this compounds class is synthetically very useful but also prone to positive outcomes in the battery of genotoxicity assays.The importance of pKa and genotoxic characteristics cannot be overestimated in drug design, where the multivariate optimisations of properties that influence the Absorption-Distribution-Metabolism-Excretion-Toxicity (ADMET) profiles now features very early on in the drug discovery process.Models were constructed using carboxylic acids in conjunction with the Quantum Topological Molecular Similarity (QTMS) method. The models produced Root Mean Square Error of Prediction (RMSEP) values of less than 0.5 pKa units and compared favourably to other pKa prediction methods. The ortho-substituted benzoic acids had the largest RMSEP which was significantly improved by splitting the compounds into high-correlation subsets. For these subsets, single-term equations containing one ab initio bond length were able to accurately predict pKa. The pKa prediction equations were extended to phenols and anilines.Quantitative Structure Activity Relationship (QSAR) models of acceptable quality were built based on literature data to predict the mutagenic potency (LogMP) of carbo- and heterocyclic aromatic amines using QTMS. However, these models failed to predict Ames test values for compounds screened at GSK. Contradictory internal and external data for several compounds motivated us to determine the fidelity of the Ames test for this compound class. The systematic investigation involved recrystallisation to purify compounds, analytical methods to measure the purity and finally comparative Ames testing. Unexpectedly, the Ames test results were very reproducible when 14 representative repurified molecules were tested as the freebase and the hydrochloride salt in two different solvents (water and DMSO). This work formed the basis for the analysis of Ames data at GSK and a systematic Ames testing programme for aromatic amines. So far, an unprecedentedly large list of 400 compounds has been made available to guide medicinal chemists. We constructed a model for the subset of 100 meta-/para-substituted anilines that could predict 70% of the Ames classifications. The experimental values of several of the model outliers appeared questionable after closer inspection and three of these have been retested so far. The retests lead to the reclassification of two of them and thereby to improved model accuracy of 78%. This demonstrates the power of the iterative process of model building, critical analysis of experimental data, retesting outliers and rebuilding the model.EThOS - Electronic Theses Online ServiceEPSRCGlaxoSmithKlineGBUnited Kingdo

Predictive QSPR study of the dissociation constants of diverse pharmaceutical compounds

The objective of the article was to perform a predictive analysis, based on quantitative structure-property relationships, of the dissociation constants (pKa) of different medicinal compounds (e.g., salicylic acid, salbutamol, lidocaine). Given the importance of this property in medicinal chemistry, it is of interest to develop theoretical methods for its prediction. The descriptors selection from a pool containing more than a thousand geometrical, topological, quantum-mechanical, and electronic types of descriptors was performed using the enhanced replacement method. Genetic algorithm and the replacement method (RM) techniques were used as reference points. A new methodology for the selection of the optimal number of descriptors to include in a model was presented and successfully used, showing that the best model should contain four descriptors. The best quantitative structure-property relationships linear model constructed using 62 molecular structures not previously used in this type of quantitative structure-property study showed good predictive attributes. The root mean squared error of the 26 molecules test set was 0.5600. The analysis of the quantitative structure-property relationships model suggests that the dissociation constants depend significantly on the number of acceptor atoms for H-bonds and on the number of carboxylic acids present in the molecules.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada

Predictive QSPR study of the dissociation constants of diverse pharmaceutical compounds

The objective of the article was to perform a predictive analysis, based on quantitative structure-property relationships, of the dissociation constants (pKa) of different medicinal compounds (e.g., salicylic acid, salbutamol, lidocaine). Given the importance of this property in medicinal chemistry, it is of interest to develop theoretical methods for its prediction. The descriptors selection from a pool containing more than a thousand geometrical, topological, quantum-mechanical, and electronic types of descriptors was performed using the enhanced replacement method. Genetic algorithm and the replacement method (RM) techniques were used as reference points. A new methodology for the selection of the optimal number of descriptors to include in a model was presented and successfully used, showing that the best model should contain four descriptors. The best quantitative structure-property relationships linear model constructed using 62 molecular structures not previously used in this type of quantitative structure-property study showed good predictive attributes. The root mean squared error of the 26 molecules test set was 0.5600. The analysis of the quantitative structure-property relationships model suggests that the dissociation constants depend significantly on the number of acceptor atoms for H-bonds and on the number of carboxylic acids present in the molecules.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada

Speciation of organoarsenicals in aqueous solutions by Raman spectrometry and quantum chemical calculations

.Knowledge about the existence and stability of different species of organoarsenicals in solution is of the most significant interest for fields so different as chemical, environmental, biological, toxicological and forensic. This work provides a comparative evaluation of the Raman spectra of four organoarsenicals (o-arsanilic acid, p-arsanilic acid, roxarsone and cacodylic acid) in aqueous solutions under acidic, neutral and alkaline conditions. Speciation of some of these organoarsenicals is possible by Raman spectrometry at different selected pHs. Further, we examine the proficiency of computational chemistry to obtain the theoretical Raman spectra of the four organoarsenicals compounds. To this end, we employ a computational protocol that includes explicit water molecules and conformational sampling, finding that the calculated organoarsenicals spectra agree reasonably well with those experimentally obtained in an aqueous solution in the whole pH range covered. Finally, we highlight the effectiveness of quantum chemical calculations to identify organoarsenicals in an aqueous solution.S

Palladium-Catalyzed Cross-Coupling Reaction of Diarylmethanol Derivatives with Diborylmethane

A palladium-catalyzed cross-coupling reaction of diarylmethanol derivatives with diborylmethane has been developed. The reaction proceeds chemoselectively to deliver the corresponding homobenzylic boronates in good yields.Asai K., Miura M., Hirano K.. Palladium-Catalyzed Cross-Coupling Reaction of Diarylmethanol Derivatives with Diborylmethane. Journal of Organic Chemistry. 87(11), 7436-7445, (2022), 3 June 2022; © 2022 American Chemical Society. https://doi.org/10.1021/acs.joc.2c00715

Hydrate crystal structures, radial distribution functions, and computing solubility

Solubility prediction usually refers to prediction of the intrinsic aqueous solubility, which is the concentration of an unionised molecule in a saturated aqueous solution at thermodynamic equilibrium at a given temperature. Solubility is determined by structural and energetic components emanating from solid-phase structure and packing interactions, solute–solvent interactions, and structural reorganisation in solution. An overview of the most commonly used methods for solubility prediction is given in Chapter 1. In this thesis, we investigate various approaches to solubility prediction and solvation model development, based on informatics and incorporation of empirical and experimental data. These are of a knowledge-based nature, and specifically incorporate information from the Cambridge Structural Database (CSD). A common problem for solubility prediction is the computational cost associated with accurate models. This issue is usually addressed by use of machine learning and regression models, such as the General Solubility Equation (GSE). These types of models are investigated and discussed in Chapter 3, where we evaluate the reliability of the GSE for a set of structures covering a large area of chemical space. We find that molecular descriptors relating to specific atom or functional group counts in the solute molecule almost always appear in improved regression models. In accordance with the findings of Chapter 3, in Chapter 4 we investigate whether radial distribution functions (RDFs) calculated for atoms (defined according to their immediate chemical environment) with water from organic hydrate crystal structures may give a good indication of interactions applicable to the solution phase, and justify this by comparison of our own RDFs to neutron diffraction data for water and ice. We then apply our RDFs to the theory of the Reference Interaction Site Model (RISM) in Chapter 5, and produce novel models for the calculation of Hydration Free Energies (HFEs)