Physicochemical Interpretation of Molecular Connectivity Indices

A novel procedure for selection of molecular descriptors: QSAR model for mutagenicity of nitroaromatic compounds

Nitroaromatic compounds (NACs) stand out as pervasive organic pollutants, prompting an imperative need to investigate their hazardous effects. Computational chemistry methods play a crucial role in this exploration, offering a safer and more time-efficient approach, mandated by various legislations. In this study, our focus lay on the development of transparent, interpretable, reproducible, and publicly available methodologies aimed at deriving quantitative structure-activity relationship models and testing them by modelling the mutagenicity of NACs against the Salmonella typhimurium TA100 strain. Descriptors were selected from Mordred and RDKit molecular descriptors, along with several quantum chemistry descriptors. For that purpose, the genetic algorithm (GA), as the most widely used method in the literature, and three alternative algorithms (Boruta, Featurewiz, and ForwardSelector) combined with the forward stepwise selection technique were used. The construction of models utilized the multiple linear regression method, with subsequent scrutiny of fitting and predictive performance, reliability, and robustness through various statistical validation criteria. The models were ranked by the Multi-Criteria Decision Making procedure. Findings have revealed that the proposed methodology for descriptor selection outperforms GA, with Featurewiz showing a slight advantage over Boruta and ForwardSelector. These constructed models can serve as valuable tools for the quick and reliable prediction of NACs mutagenicity.

Correlation of the Molecular Cross-Sectional Area of Organic Monofunctional Compounds with Topological Descriptors

The molecular cross-sectional area (σ) has proved to be an interesting molecular measure not only in the field of adsorption phenomena on solids but also in biochemistry, physiology, or surfactant chemistry. The existing methods to estimate the cross-sectional areas are either not readily applicable or can only be applied to a limited number of compounds. The aim of this work was to describe a method, as general as possible, quick and easy to perform. To that end, the molecular cross-sectional areas were correlated with topological indices. The Emmett-Brunauer formula was used to calculate the reference cross-sectional areas (σEB) of 431 compounds. The correlations of the Wiener (W), hyper-Wiener (WW), Balaban (J), and Randić (χ) indices with σEB were compared for n-alkanes as well as branched and cyclic alkanes. Only the Wiener and hyper-Wiener indices correlated well with σEB, with the data being best fitted by power law regression curves. The lower degeneracy of the hyper-Wiener index did not translate into any significant gain of performance when correlated with σEB. Following the parsimony principle, the less complex Wiener index was thus selected to correlate with the σEB of compounds representing 31 other monofunctional and structural families. The integration of all the compound families into a single curve allowed a quick rough estimation of the cross-sectional areas. The specific reference equations σEB = qWp were determined for the 34 selected families, allowing the fast and reliable calculation of the cross-sectional area of any monofunctional compound.

QSAR assessment of aquatic toxicity potential of diverse agrochemicals

The fast-increasing number of commercially produced chemicals challenges the experimental ecotoxicity assessment methods, which are costly, time-consuming, and dependent on the sacrifice of animals. In this regard, Quantitative Structure-Property/Activity Relationships (QSPR/QSAR) have led the way in developing ecotoxicity assessment models. In this study, QSAR models have been developed using the pEC50 values of 82 diverse agrochemicals or agro-molecules against a planktonic crustacean Daphnia magna with easily interpretable 2D descriptors. Moreover, a link among octanol-water partition coefficient (KOW), bio-concentration factor (BCF), and critical body residue (CBR) has been addressed, and their imputation for the prediction of the toxicity endpoint (EC50) has been done with an objective of the advanced exploration of several ecotoxicological parameters for toxic chemicals. The developed partial least squares (PLS) models were validated rigorously and proved to be robust, sound, and immensely well-predictive. The final Daphnia toxicity model derived from experimental derived properties along with computed descriptors emerged better in statistical quality and predictivity than those obtained solely from computed descriptors. Additionally, the pEC50 and other important properties (log KOW, log BCF, and log CBR) for a set of external agro-molecules, not employed in model development, were predicted to show the predictive ability of the models.

Multi-Condition QSAR Model for the Virtual Design of Chemicals with Dual Pan-Antiviral and Anti-Cytokine Storm Profiles

Respiratory viruses are infectious agents, which can cause pandemics. Although nowadays the danger associated with respiratory viruses continues to be evidenced by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as the virus responsible for the current COVID-19 pandemic, other viruses such as SARS-CoV-1, the influenza A and B viruses (IAV and IBV, respectively), and the respiratory syncytial virus (RSV) can lead to globally spread viral diseases. Also, from a biological point of view, most of these viruses can cause an organ-damaging hyperinflammatory response known as the cytokine storm (CS). Computational approaches constitute an essential component of modern drug development campaigns, and therefore, they have the potential to accelerate the discovery of chemicals able to simultaneously inhibit multiple molecular and nonmolecular targets. We report here the first multicondition model based on quantitative structure-activity relationships and an artificial neural network (mtc-QSAR-ANN) for the virtual design and prediction of molecules with dual pan-antiviral and anti-CS profiles. Our mtc-QSAR-ANN model exhibited an accuracy higher than 80%. By interpreting the different descriptors present in the mtc-QSAR-ANN model, we could retrieve several molecular fragments whose assembly led to new molecules with drug-like properties and predicted pan-antiviral and anti-CS activities.

Incorporating Machine Learning in Computer-Aided Molecular Design for Fragrance Molecules

The demand for new novel flavour and fragrance (F&F) molecules has boosted the need for a systematic approach to designing fragrance molecules. However, the F&F-related industry still relies heavily on experimental approaches or on existing databases without considering the consequences resulting from changes in concentration, which could omit potential fragrances. Computer-aided molecular design (CAMD) has great potential to identify novel molecular structures to be used as fragrances. Using CAMD for this purpose requires models to predict the olfaction properties of molecules. A rough set-based machine learning (RSML) approach is used to develop an interpretable predictive model for odour characteristics in this work. New rule-based models are generated from RSML based on the dilution and a number of different topological indices which identify the structure-odour relationship of fragrance molecules. The most prominent rules are selected and formulated as constraints in a CAMD optimisation model. The combination of several rules was able to increase the coverage of different classes of molecules. To model the performance indicators that vary over a range of properties, a disjunctive programming model is also incorporated into the CAMD framework. A case study demonstrates the utilisation of this methodology to design fragrance additives in dishwashing liquid. The results illustrate the capability of the novel RSML and CAMD framework to identify potential fragrance molecules that can be used in consumer products.

PTML Modeling for Pancreatic Cancer Research: In Silico Design of Simultaneous Multi-Protein and Multi-Cell Inhibitors

Pancreatic cancer (PANC) is a dangerous type of cancer that is a major cause of mortality worldwide and exhibits a remarkably poor prognosis. To date, discovering anti-PANC agents remains a very complex and expensive process. Computational approaches can accelerate the search for anti-PANC agents. We report for the first time two models that combined perturbation theory with machine learning via a multilayer perceptron network (PTML-MLP) to perform the virtual design and prediction of molecules that can simultaneously inhibit multiple PANC cell lines and PANC-related proteins, such as caspase-1, tumor necrosis factor-alpha (TNF-alpha), and the insulin-like growth factor 1 receptor (IGF1R). Both PTML-MLP models exhibited accuracies higher than 78%. Using the interpretation from one of the PTML-MLP models as a guideline, we extracted different molecular fragments desirable for the inhibition of the PANC cell lines and the aforementioned PANC-related proteins and then assembled some of those fragments to form three new molecules. The two PTML-MLP models predicted the designed molecules as potentially versatile anti-PANC agents through inhibition of the three PANC-related proteins and multiple PANC cell lines. Conclusions: This work opens new horizons for the application of the PTML modeling methodology to anticancer research.

In Silico Drug Repurposing for Anti-Inflammatory Therapy: Virtual Search for Dual Inhibitors of Caspase-1 and TNF-Alpha

Inflammation involves a complex biological response of the body tissues to damaging stimuli. When dysregulated, inflammation led by biomolecular mediators such as caspase-1 and tumor necrosis factor-alpha (TNF-alpha) can play a detrimental role in the progression of different medical conditions such as cancer, neurological disorders, autoimmune diseases, and cytokine storms caused by viral infections such as COVID-19. Computational approaches can accelerate the search for dual-target drugs able to simultaneously inhibit the aforementioned proteins, enabling the discovery of wide-spectrum anti-inflammatory agents. This work reports the first multicondition model based on quantitative structure–activity relationships and a multilayer perceptron neural network (mtc-QSAR-MLP) for the virtual screening of agency-regulated chemicals as versatile anti-inflammatory therapeutics. The mtc-QSAR-MLP model displayed accuracy higher than 88%, and was interpreted from a physicochemical and structural point of view. When using the mtc-QSAR-MLP model as a virtual screening tool, we could identify several agency-regulated chemicals as dual inhibitors of caspase-1 and TNF-alpha, and the experimental information later retrieved from the scientific literature converged with our computational results. This study supports the capabilities of our mtc-QSAR-MLP model in anti-inflammatory therapy with direct applications to current health issues such as the COVID-19 pandemic.

Indirect-Acting Pan-Antivirals vs. Respiratory Viruses: A Fresh Perspective on Computational Multi-Target Drug Discovery

Respiratory viruses continue to afflict mankind. Among them, pathogens such as coronaviruses [including the current pandemic agent known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)] and the one causing influenza A (IAV) are highly contagious and deadly. These can evade the immune system defenses while causing a hyperinflammatory response that can damage different tissues/organs. Simultaneously targeting immunomodulatory proteins is a plausible antiviral strategy since it could lead to the discovery of indirect-acting pan-antiviral (IAPA) agents for the treatment of diseases caused by respiratory viruses. In this context, computational approaches, which are an essential part of the modern drug discovery campaigns, could accelerate the identification of multi-target immunomodulators. This perspective discusses the usefulness of computational multi-target drug discovery for the virtual screening (drug repurposing) of IAPA agents capable of boosting the immune system through the activation of the toll-like receptor 7 (TLR7) and/or the stimulator of interferon genes (STING) while inhibiting key pro-inflammatory proteins, such as caspase-1 and tumor necrosis factor-alpha (TNF-α).

QSAR models for the acute toxicity of 1,2,4-triazole fungicides to zebrafish (Danio rerio) embryos

In recent decades, the 1,2,4-triazole fungicides are widely used for crop diseases control, and their toxicity to wild lives and pollution to ecosystem have attracted more and more attention. However, how to quickly and efficiently evaluate the toxicity of these compounds to environmental organisms is still a challenge. In silico method, such like Quantitative Structure-Activity Relationship (QSAR), provides a good alternative to evaluate the environmental toxicity of a large number of chemicals. At the present study, the acute toxicity of 23 1,2,4-triazole fungicides to zebrafish (Danio rerio) embryos was firstly tested, and the LC₅₀ (median lethal concentration) values were used as the bio-activity endpoint to conduct QSAR modelling for these triazoles. After the comparative study of several QSAR models, the 2D-QSAR model was finally constructed using the stepwise multiple linear regression algorithm combining with two physicochemical parameters (logD and μ), an electronic parameter (Q_N1) and a topological parameter (X^vPC₄). The optimal model could be mathematically described as following: pLC₅₀ = -7.24-0.30X^vPC₄ + 0.76logD - 26.15Q_N1 - 0.08μ. The internal validation by leave-one-out (LOO) cross-validation showed that the R²_adj (adjusted noncross-validation squared correlation coefficient), Q² (cross-validation correlation coefficient) and RMSD (root-mean-square error) was 0.88, 0.84 and 0.17, respectively. The external validation indicated the model had a robust predictability with the q² (predictive squared correlation coefficient) of 0.90 when eliminated tricyclazole. The present study provided a potential tool for predicting the acute toxicity of new 1,2,4-triazole fungicides which contained an independent triazole ring group in their molecules to zebrafish embryos, and also provided a reference for the development of more environmentally-friendly 1,2,4-triazole pesticides in the future.

Topological Indices of Hyaluronic Acid-Paclitaxel Conjugates’ Molecular Structure in Cancer Treatment

A large number of medical experiments have confirmed that the features of drugs have a close correlation with their molecular structure. Drug properties can be obtained by studying the molecular structure of corresponding drugs. The calculation of the topological index of a drug structure enables scientists to have a better understanding of the physical chemistry and biological characteristics of drugs. In this paper, we focus on Hyaluronic Acid-Paclitaxel conjugates which are widely used in the manufacture of anticancer drugs. Several topological indices are determined by virtue of the edge-partition method, and our results remedy the lack of medicine experiments, thus providing a theoretical basis for pharmaceutical engineering.

Exploring the QSAR's predictive truthfulness of the novel N-tuple discrete derivative indices on benchmark datasets

Graph derivative indices (GDIs) have recently been defined over N-atoms (N = 2, 3 and 4) simultaneously, which are based on the concept of derivatives in discrete mathematics (finite difference), metaphorical to the derivative concept in classical mathematical analysis. These molecular descriptors (MDs) codify topo-chemical and topo-structural information based on the concept of the derivative of a molecular graph with respect to a given event (S) over duplex, triplex and quadruplex relations of atoms (vertices). These GDIs have been successfully applied in the description of physicochemical properties like reactivity, solubility and chemical shift, among others, and in several comparative quantitative structure activity/property relationship (QSAR/QSPR) studies. Although satisfactory results have been obtained in previous modelling studies with the aforementioned indices, it is necessary to develop new, more rigorous analysis to assess the true predictive performance of the novel structure codification. So, in the present paper, an assessment and statistical validation of the performance of these novel approaches in QSAR studies are executed, as well as a comparison with those of other QSAR procedures reported in the literature. To achieve the main aim of this research, QSARs were developed on eight chemical datasets widely used as benchmarks in the evaluation/validation of several QSAR methods and/or many different MDs (fundamentally 3D MDs). Three to seven variable QSAR models were built for each chemical dataset, according to the original dissection into training/test sets. The models were developed by using multiple linear regression (MLR) coupled with a genetic algorithm as the feature wrapper selection technique in the MobyDigs software. Each family of GDIs (for duplex, triplex and quadruplex) behaves similarly in all modelling, although there were some exceptions. However, when all families were used in combination, the results achieved were quantitatively higher than those reported by other authors in similar experiments. Comparisons with respect to external correlation coefficients (q²_ext) revealed that the models based on GDIs possess superior predictive ability in seven of the eight datasets analysed, outperforming methodologies based on similar or more complex techniques and confirming the good predictive power of the obtained models. For the q²_ext values, the non-parametric comparison revealed significantly different results to those reported so far, which demonstrated that the models based on DIVATI's indices presented the best global performance and yielded significantly better predictions than the 12 0-3D QSAR procedures used in the comparison. Therefore, GDIs are suitable for structure codification of the molecules and constitute a good alternative to build QSARs for the prediction of physicochemical, biological and environmental endpoints.

Improved Scaling of Molecular Network Calculations: The Emergence of Molecular Domains

The design of materials needed for the storage, delivery, and conversion of (re)useable energy is still hindered by the lack of new, hierarchical molecular screening methodologies that encode information on more than one length scale. Using a molecular network theory as a foundation, we show that to describe charge transport in disordered materials the network methodology must be scaled-up. We detail the scale-up through the use of adjacency lists and depth first search algorithms for during operations on the adjacency matrix. We consider two types of electronic acceptors, perylenediimide (PDI) and the fullerene derivative phenyl-C61-butyric acid methyl ester (PCBM), and we demonstrate that the method is scalable to length scales relevant to grain boundary and trap formations. Such boundaries lead to a decrease in the percolation ratio of PDI with system size, while the ratio for PCBM remains constant, further quantifying the stable, diverse transport pathways of PCBM and its success as a charge-accepting material.

Physico-Chemical and Structural Interpretation of Discrete Derivative Indices on N-Tuples Atoms

This report examines the interpretation of the Graph Derivative Indices (GDIs) from three different perspectives (i.e., in structural, steric and electronic terms). It is found that the individual vertex frequencies may be expressed in terms of the geometrical and electronic reactivity of the atoms and bonds, respectively. On the other hand, it is demonstrated that the GDIs are sensitive to progressive structural modifications in terms of: size, ramifications, electronic richness, conjugation effects and molecular symmetry. Moreover, it is observed that the GDIs quantify the interaction capacity among molecules and codify information on the activation entropy. A structure property relationship study reveals that there exists a direct correspondence between the individual frequencies of atoms and Hückel’s Free Valence, as well as between the atomic GDIs and the chemical shift in NMR, which collectively validates the theory that these indices codify steric and electronic information of the atoms in a molecule. Taking in consideration the regularity and coherence found in experiments performed with the GDIs, it is possible to say that GDIs possess plausible interpretation in structural and physicochemical terms.

First Multitarget Chemo-Bioinformatic Model To Enable the Discovery of Antibacterial Peptides against Multiple Gram-Positive Pathogens

Antimicrobial peptides (AMPs) have emerged as promising therapeutic alternatives to fight against the diverse infections caused by different pathogenic microorganisms. In this context, theoretical approaches in bioinformatics have paved the way toward the creation of several in silico models capable of predicting antimicrobial activities of peptides. All current models have several significant handicaps, which prevent the efficient search for highly active AMPs. Here, we introduce the first multitarget (mt) chemo-bioinformatic model devoted to performing alignment-free prediction of antibacterial activity of peptides against multiple Gram-positive bacterial strains. The model was constructed from a data set containing 2488 cases of AMPs sequences assayed against at least 1 out of 50 Gram-positive bacterial strains. This mt-chemo-bioinformatic model displayed percentages of correct classification higher than 90.00% in both training and prediction (test) sets. For the first time, two computational approaches derived from basic concepts in genetics and molecular biology were applied, allowing the calculations of the relative contributions of any amino acid (in a defined position) to the antibacterial activity of an AMP and depending on the bacterial strain used in the biological assay. The present mt-chemo-bioinformatic model constitutes a powerful tool to enable the discovery of potent and versatile AMPs.

A general ANN-based multitasking model for the discovery of potent and safer antibacterial agents

Bacteria have been one of the world's most dangerous and deadliest pathogens for mankind, nowadays giving rise to significant public health concerns. Given the prevalence of these microbial pathogens and their increasing resistance to existing antibiotics, there is a pressing need for new antibacterial drugs. However, development of a successful drug is a complex, costly, and time-consuming process. Quantitative Structure-Activity Relationships (QSAR)-based approaches are valuable tools for shortening the time of lead compound identification but also for focusing and limiting time-costly synthetic activities and in vitro/vivo evaluations. QSAR-based approaches, supported by powerful statistical techniques such as artificial neural networks (ANNs), have evolved to the point of integrating dissimilar types of chemical and biological data. This chapter reports an overview of the current research and potential applications of QSAR modeling tools toward the rational design of more efficient antibacterial agents. Particular emphasis is given to the setup of multitasking models along with ANNs aimed at jointly predicting different antibacterial activities and safety profiles of drugs/chemicals under diverse experimental conditions.

Stability prediction of Cu2+, Ni2+ and Zn2+ N-salicylidene-aminoacidato complexes by models based on connectivity index 3 χ v

This paper presents models for the estimation of stability constants (K 1 and β 2) of nickel(II), copper(II) and zinc(II) mono- and bis-complexes with 5 Schiff bases (salicylideneglycine, salicylidenealanine, salicylideneserine, salicylidenephenylalanine, and salicylidenetyrosine). The models were based on the molecular-graph theory and valence molecular connectivity index of the 3rd order, 3χv , derived from it. Univariate linear models were developed for each metal separately, while in the common models for two and three metals, the indicator variable, In, was introduced. The standard error of models for the log K 1 constant was less than 0.12, while for log β 2 models, the S.E. did not exceed 0.14.

Chemoinformatics for rational discovery of safe antibacterial drugs: simultaneous predictions of biological activity against streptococci and toxicological profiles in laboratory animals

Streptococci are a group of Gram-positive bacteria which are responsible for causing many diverse diseases in humans and other animals worldwide. The high prevalence of resistance of these bacteria to current antibacterial drugs is an alarming problem for the scientific community. The battle against streptococci by using antimicrobial chemotherapies will depend on the design of new chemicals with high inhibitory activity, having also as low toxicity as possible. Multi-target approaches based on quantitative-structure activity relationships (mt-QSAR) have played a very important role, providing a better knowledge about the molecular patterns related with the appearance of different pharmacological profiles including antimicrobial activity. Until now, almost all mt-QSAR models have considered the study of biological activity or toxicity separately. In the present study, we develop by the first time, a unified multitasking (mtk) QSAR model for the simultaneous prediction of anti-streptococci activity and toxic effects against biological models like Mus musculus and Rattus norvegicus. The mtk-QSAR model was created by using artificial neural networks (ANN) analysis for the classification of compounds as positive (high biological activity and/or low toxicity) or negative (otherwise) under diverse sets of experimental conditions. Our mtk-QSAR model, correctly classified more than 97% of the cases in the whole database (more than 11,500 cases), serving as a promising tool for the virtual screening of potent and safe anti-streptococci drugs.

Model of twelve properties of a set of organic solvents with graph-theoretical and/or experimental parameters

Twelve properties of a highly heterogeneous class of organic solvents have been modeled with a graph-theoretical molecular connectivity modified (MC) method, which allows to encode the core electrons and the hydrogen atoms. The graph-theoretical method uses the concepts of simple, general, and complete graphs, where these last types of graphs are used to encode the core electrons. The hydrogen atoms have been encoded by the aid of a graph-theoretical perturbation parameter, which contributes to the definition of the valence delta, delta(v), a key parameter in molecular connectivity studies. The model of the twelve properties done with a stepwise search algorithm is always satisfactory, and it allows to check the influence of the hydrogen content of the solvent molecules on the choice of the type of descriptor. A similar argument holds for the influence of the halogen atoms on the type of core electron representation. In some cases the molar mass, and in a minor way, special "ad hoc" parameters have been used to improve the model. A very good model of the surface tension could be obtained by the aid of five experimental parameters. A mixed model method based on experimental parameters plus molecular connectivity indices achieved, instead, to consistently improve the model quality of five properties. To underline is the importance of the boiling point temperatures as descriptors in these last two model methodologies.

Design of novel antituberculosis compounds using graph-theoretical and substructural approaches

The increasing resistance of Mycobacterium tuberculosis to the existing drugs has alarmed the worldwide scientific community. In an attempt to overcome this problem, two models for the design and prediction of new antituberculosis agents were obtained. The first used a mixed approach, containing descriptors based on fragments and the topological substructural molecular design approach (TOPS-MODE) descriptors. The other model used a combination of two-dimensional (2D) and three-dimensional (3D) descriptors. A data set of 167 compounds with great structural variability, 72 of them antituberculosis agents and 95 compounds belonging to other pharmaceutical categories, was analyzed. The first model showed sensitivity, specificity, and accuracy values above 80% and the second one showed values higher than 75% for these statistical indices. Subsequently, 12 structures of imidazoles not included in this study were designed, taking into account the two models. In both cases accuracy was 100%, showing that the methodology in silico developed by us is promising for the rational design of antituberculosis drugs.

Self-assembled monolayer of organic iodine on a Au surface for attachment of redox-active metal clusters

The attachment of a bifunctional iodo-organo-phosphinate compound to gold (Au) surfaces via chemisorption of the iodine atom is described and used to chelate a redox-active metal cluster via the phosphinate group. XPS, AFM, and electrochemical measurements show that (4-iodo-phenyl)phenyl phosphinic acid (IPPA) forms a tightly bound self-assembled monolayer (SAM) on Au surfaces. The surface coverage of an IPPA monolayer on Au was quantified by an electrochemical method and found to be 0.40 +/- 0.03 nmol/cm2, roughly corresponding to 0.4 monolayers. We show that the Au/IPPA SAM, but not the underivatized Au, adsorbs Mn4O4(Ph2PO2)6 from solution by a phosphinate exchange reaction to yield Au/IPPA/Mn4O4(Ph2PO2)5 SAM. The resulting SAM is firmly bound and not removed by sonication, as confirmed by manganese XPS (Mn 2p1/2) and by AFM. Electrochemistry confirms that Mn4O4(Ph2PO2)6 is anchored on the Au/IPPA surface and that redox chemistry can be mediated between the electrode and the surface-attached complex. Mn4O4(Ph2PO2)6 contains the reactive Mn4O46+ cubane core, a redox-active bioinspired catalyst.

Three-dimensional topographic index applied to the prediction of acyclic C5–C8 alkenes Kováts retention indices on polydimethylsiloxane and squalane columns

A novel approach is described for the prediction of gas chromatographic Kováts retention indices of 150 acyclic C5-C8 alkenes on two stationary phases (polydimethylsiloxane, PDMS, and squalane, SQ). The heuristic method was used to build multiple linear regression models using descriptors calculated by MODLESLAB software and CODESSA program. The resulting quantitative structure-retention relationship (QSRR) models were well-correlated, with predictive R2 values of 0.970 and 0.958 for retention indices on PDMS and SQ columns, respectively. 1Omegap, a three-dimensional (3D) topographic index, was found to play the most important role in the description of the chromatographic retention behavior of the alkenes in these two stationary phases. Moreover, this index could completely distinguish different isomers of alkene. Therefore, it can also be extended to distinguish different isomers of other compounds so that can well describe their quantitative structure-retention relationships.

Quantitative study of the structure-retention index relationship in the imine family

The Kováts retention index is one of the most popular descriptors of the performance of organic compounds in gas chromatography (GC). The mathematical modeling of this index is an interesting and open problem in analytical chemistry. In this paper, two models for the prediction of the Kováts retention index are presented. Topologic, topographic and quantum-chemical descriptors were used as structural descriptors. Multiple linear regression (MLR) analysis provides the first model using the forward stepwise procedure for the variable selection. For the second one, an ensemble of artificial neural network (ANN) was constructed using the pruning algorithm. Both methods were validated by an external set of compounds, by the Golbraikh and Tropsha method and by the leave-one-out (LOO) and the leave many out (LMO) procedures. The R2, RMScv and Q2, values for the training sets were 0.884, 0.589 and 0.830 for NN and 0.974, 0.417 and 0.970 for MLR models, respectively. The robustness of both models was demonstrated. Both portrait the chromatographic performance of the sample but in this case, the results of MLR equation are better than the NN ones. The MLR model is recommended because of its simplicity.

Partial least squares modelling of the acute toxicity of aliphatic compounds to Tetrahymena pyriformis

The aim of this study was to evaluate a multivariate statistical model, utilising Partial Least Squares (PLS) analysis, for the prediction of the acute toxicity of aliphatic chemicals to the ciliate Tetrahymena pyriformis. A model was developed that was capable of making a prediction regardless the mechanism of toxic action. The toxicity of 476 compounds, possessing different mechanisms of toxic action was considered. A set of 74 descriptors, including the octanol-water partition coefficient, molecular-orbital descriptors, geometrical, topological and connectivity indices, was generated. A three-component, eight-descriptor PLS model was developed. It was validated by a Y-permutation test and by simulation of external prediction for complementary subsets. A comparison with existing class or mechanism-based models, derived on the same data set, was made.

Complete graph conjecture for inner-core electrons: homogeneous index case

The complete graph conjecture that encodes the inner-core electrons of atoms with principal quantum number n >or= 2 with complete graphs, and especially with odd complete graphs, is discussed. This conjecture is used to derive new values for the molecular connectivity and pseudoconnectivity basis indices of hydrogen-suppressed chemical pseudographs. For atoms with n = 2 the new values derived with this conjecture are coincident with the old ones. The modeling ability of the new homogeneous basis indices, and of the higher-order terms, is tested and compared with previous modeling studies, which are centered on basis indices that are either based on quantum concepts or partially based on this new conjecture for the inner-core electrons. Two similar algorithms have been proposed with this conjecture, and they parallel the two "quantum" algorithms put forward by molecular connectivity for atoms with n > 2. Nine properties of five classes of compounds have been tested: the molecular polarizabilities of a class of organic compounds, the dipole moment, molar refraction, boiling points, ionization energies, and parachor of a series of halomethanes, the lattice enthalpy of metal halides, the rates of hydrogen abstraction of chlorofluorocarbons, and the pED(50) of phenylalkylamines. The two tested algorithms based on the odd complete graph conjecture give rise to a highly interesting model of the nine properties, and three of them can even be modeled by the same set of basis indices. Interesting is the role of some basis indices all along the model.

Prediction of molecular volume and surface of alkanes by molecular topology

Molecular volume and molecular surface are expressed as a function of topological degree in alkane graphs. This allows not only a straightforward approach to calculate such physicochemical magnitudes but also an interpretation of the role of the local vertex invariant (LOVI) or valence degree, delta, as well as the connectivity indices in the prediction of physicochemical properties. The interpretation is based on the concept of molecular accessibility (as introduced by Estrada, J. Phys. Chem. A 2002, 106, 9085) for which precise mathematical definitions are provided.

What are the limits of applicability for graph theoretic descriptors in QSPR/QSAR? Modeling dipole moments of aromatic compounds with TOPS-MODE descriptors

The numerous possibilities of using graph theoretic descriptors in QSPR/QSAR are analyzed, and some misunderstandings on the role of this theoretical approach in chemistry are clarified. Principal component analysis is used to obtain a property space for several physicochemical properties of aromatic compounds. It is proved that most of the QSPR applications of the graph-theoretic structure descriptors are concentrated to the description of properties in a very limited region of this property space. Here, we show that graph-theoretic approaches are also applicable to the modeling of physicochemical properties that are far away from this region traditionally studied. The molecular dipole moments of benzene derivatives, mono-, ortho-, meta-, and para-susbtituted, are modeled by using the Topological Sub-Structural Molecular Design (TOPS-MODE) approach. The TOPS-MODE approach used permits to calculate group dipole moments that are given for several substituents. The differences between these group dipoles and those obtained by simple difference between experimental values are analyzed. Some difficulties arising from this traditional way of deriving substituent constants are identified and analyzed.