Prediction of drug solubility from molecular structure using a drug-like training set

Using a training set of 191 drug-like compounds extracted from the AQUASOL database a quantitative structure-property relationship (QSPR) study was conducted employing a set of simple structural and physicochemical properties to predict aqueous solubility. The resultant regression model comprised five parameters (ClogP, molecular weight, indicator variable for aliphatic amine groups, number of rotatable bonds and number of aromatic rings) and demonstrated acceptable statistics (r ² = 0.87, s = 0.51, F = 243.6, n = 191). The model was applied to two test sets consisting of a drug-like set of compounds (r ² = 0.80, s = 0.68, n = 174) and a set of agrochemicals (r ² = 0.88, s = 0.65, n = 200). Using the established general solubility equation (GSE) on the training and drug-like test set gave poorer results than the current study. The agrochemical test set was predicted with equal accuracy using the GSE and the QSPR equation. The results of this study suggest that increasing molecular size, rigidity and lipophilicity decrease solubility whereas increasing conformational flexibility and the presence of a non-conjugated amine group increase the solubility of drug-like compounds. Indeed, the proposed structural parameters make physical sense and provide simple guidelines for modifying solubility during lead optimisation.     

CoMFA and CoMSIA studies on a new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line

Abstract  

A new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line is examined to determine the relationship between the structural properties and the biological activity of these compounds—the 3-D quantitative structure–activity relationship (3D-QSAR)—using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMFA and CoMSIA models were obtained using the atom-based alignment of 33 compounds, 22 training compounds and 11 tested compounds, and these give desirable statistics; those for the CoMFA standard model were: r _cv² = 0.691, r ² = 0.998, S _press = 0.178, s = 0.014 and F = 1080.765, while CoMSIA combined steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields: r _cv² = 0.600, r ² = 0.988, S _press = 0.206, s = 0.034 and F = 284.433. The 3D-QSAR models calculated satisfactory test set activities. The 3D-QSAR contour plots correlated strongly with the experimental data for the binding topology. For this reason, these results would be beneficial for predicting affinities with the compounds of interest, and they are advantageous for guiding the design and synthesis of new and more effective anticancer agents.     

QSPR Modeling of Gibbs Free Energy of Organic Compounds by Weighting of Nearest Neighboring Codes

We examine the encoding of chemical structure of organic compounds by Labeled Hydrogen-Filled Graphs (LHFGs). Quantitative Structure-Property Relationships (QSPR) for a representative set of 150 organic molecules have been derived by means of the optimization of correlation weights of local invariants of the LHFGs. We have tested as local invariants Morgan extended connectivity of zero- and first order, numbers of path of length 2 (P2) and valence shells of distance of 2 (S2) associated with each atom in the molecular structure, and the Nearest Neighboring Codes (NNC). The best statistical characteristics for the Gibbs free energy has been obtained for the NNC weighting. Statistical parameters corresponding to this model are the following n = 100, r² = 0.9974, s = 5.136 kJ/mol, F = 38319 (training set); n = 50, r² = 0.9990, s = 3.405 kJ/mol, F = 48717 (test set). Some possible further developments are pointed out.     

CoMFA and CoMSIA studies on a new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line

Abstract  A new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line is examined to determine the relationship between the structural properties and the biological activity of these compounds—the 3-D quantitative structure–activity relationship (3D-QSAR)—using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMFA and CoMSIA models were obtained using the atom-based alignment of 33 compounds, 22 training compounds and 11 tested compounds, and these give desirable statistics; those for the CoMFA standard model were: r _cv² = 0.691, r ² = 0.998, S _press = 0.178, s = 0.014 and F = 1080.765, while CoMSIA combined steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields: r _cv² = 0.600, r ² = 0.988, S _press = 0.206, s = 0.034 and F = 284.433. The 3D-QSAR models calculated satisfactory test set activities. The 3D-QSAR contour plots correlated strongly with the experimental data for the binding topology. For this reason, these results would be beneficial for predicting affinities with the compounds of interest, and they are advantageous for guiding the design and synthesis of new and more effective anticancer agents. Graphical abstract   A new and more effective anticancer agent of xanthone derivatives against the oral human epidermoid carcinoma (KB) cell line, as investigated by CoMFA and CoMSIA analysis     

The pressure effect study on the burning rate of ammonium nitrate-HTPB-based propellant with the influence catalysts

The burning rate of AN–HTPB-based propellant catalysed with chromium salt has been studied using conventional strand burner under the various pressure range, i.e. from atmospheric pressure to 6.897 MPa and verified with Piobert law, i.e. r = aP ⁿ . At atmospheric pressure, the burning rate AN–HTPB propellant was being accelerated with the chromium-based catalysts used. In case of lead chromate-catalysed system, burning rate was observed 2.655 times higher than burning rate (r = 0.200 mm s⁻¹) of virgin AN–HTPB propellant sample. However, the Copper chromate-catalysed propellant burned with slower rate (r = 0.160 mm s⁻¹) than the virgin AN–HTPB propellant sample. The burning rate of all catalysed propellant samples are found to be the pressure sensitive and accelerated higher with rise of pressure. The highest burning rate (r = 2.422 mm s⁻¹) was recorded with ammonium dichromate and lowest (r = 1.40 mm s⁻¹) with lead chromate-catalysed propellant sample with the rise of pressure up to 6.897 MPa at different pressures. A linear relationship was observed between the burning rate and pressure rise which followed the Piobert law, i.e. r = aP ⁿ . The pressure index (n) values of AN–HTPB-based samples were calculated higher when catalysed with ammonium dichromate, Copper Chromate, Cr₂O₃, Potassium dichromate (n = 0.525, 0.555, 0.429, and 0.408 respectively) and lower (n = 0.226) with lead chromate compared to virgin sample (n = 0.405). Higher value indicates the positive effect on accelerating the burning rate with catalyst at higher pressure ranges.     

QM/MM based 3D QSAR models for potent B-Raf inhibitors

Three dimensional (3D) quantitative structure-activity relationship studies of 37 B-Raf inhibitors, pyrazole-based derivatives, were performed. Based on the co-crystallized compound (PDB ID: 3D4Q), several alignment methods were utilized to derive reliable comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models. Receptor-guided alignment with quantum mechanics/molecular mechanics (QM/MM) minimization led to the best CoMFA model (q ² = 0.624, r ² = 0.959). With the same alignment, a statistically reliable CoMSIA model with steric, H-bond acceptor, and hydrophobic fields was also derived (q ² = 0.590, r ² = 0.922). Both models were validated with an external test set, which gave satisfactory predictive r ² values of 0.926 and 0.878, respectively. Contour maps from CoMFA and CoMSIA models revealed important structural features responsible for increasing biological activity within the active site and explained the correlation between biological activity and receptor-ligand interactions. New fragments were identified as building blocks which can replace R1-3 groups through combinatorial screening methods. By combining these fragments a compound with a high bioactivity level prediction was found. These results can offer useful information for the design of new B-Raf inhibitors.     

QSPR Modeling of Complex Stability by Optimization of Correlation Weights of the Hydrogen Bond Index and the Local Graph Invariants

The hydrogen bond index (HBI) is the global invariant of a molecular graph and equals the number of vertices representing hydrogen and nitrogen atoms. This index was considered a measure of the capability of a complex to form hydrogen bonds. Optimization of the correlation weights of the HBI and local graph invariants was used for the QSPR modeling of the stability of 110 biometal M²⁺ complexes with -amino acids and phosphate derivatives of adenosine. The statistical parameters of the best model are n = 55, r = 0.9921, s = 0.279, and F = 3328 (learning sample) and n = 55, r = 99.35, s = 0.248, and F = 4027 (control sample).     

Pharmacophore modeling, 3D-QSAR,docking, and molecular dynamics simulation on topoisomerase IV inhibitors of wild type Staphylococcus aureus

Staphylococcus aureus is a gram-positive bacterium. It is a foremost cause of skin and respiratory infections, endocarditis, osteomyelitis, Ritter’s disease, and bacteraemia. Topoisomerase enzyme is involved in preventing or correcting topological problems of overwinding or underwinding occurring in DNA before replication process. An exhaustive molecular modeling studies that includes pharmacophore modeling, ligand-based three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, molecular dynamics simulation, and ADME calculations were performed on isothiazoloquinolones derivatives which are reported as effective inhibitors against topoisomerase IV of wild type S. aureus. In pharmacophore modeling by using pharmacophore alignment and scoring engine (PHASE) a five-point model (AHHRR.3) was generated with existing compounds having statistical significant as correlation coefficient (R ² = 0.954), cross-validation coefficient (Q ² = 0.650), and F value of 130.5. Ligand-based 3D-QSAR study was applied using comparative molecular field analysis (CoMFA) with Q ² = 0.616, R ² = 0.989, and comparative molecular similarity indices analysis (CoMSIA) with Q ² = 0.510, R ² = 0.995. The predictive ability of this model was determined using a test set of molecules that gave acceptable predictive correlation (R ² Pred) values 0.55 and 0.56 for CoMFA and CoMSIA, respectively. Docking and molecular dynamic simulations were employed to position the inhibitors into protein active site to find out the most probable binding mode and most reliable conformations. Developed pharmacophore models and docking methods provide guidance to design enhanced activity molecules.

     

QSPR Modeling of Complex Stability by Correlation Weighing of the Topological and Chemical Invariants of Molecular Graphs

The concept of modeling of the complex stability based on optimization of the correlation weights of the nearest neighborhood codes (NNC), the hydrogen bond index (HBI), and the cyclicity code (CC) is described. The NNC is a local topochemical invariant of a vertex of the molecular graph whose numerical value is a function of the total number and the composition of vertices adjacent to the given vertex. The HBI is a global chemical invariant of the molecular graph calculated from the number of oxygen and nitrogen atoms. The CC is a global topological invariant of the graph equal to the number of rings present in the ligand structure. The statistical characteristics of the best model of the stability constants of the complexes are as follows: n = 75, r = 0.9738, s = 0.457, F = 1337 (training sample); n = 75, r = 0.9795, s = 0.461, F = 1724 (test sample).