Next Article in Journal
Computational Studies of a Mechanism for Binding and Drug Resistance in the Wild Type and Four Mutations of HIV-1 Protease with a GRL-0519 Inhibitor
Next Article in Special Issue
Development of an in Silico Model of DPPH• Free Radical Scavenging Capacity: Prediction of Antioxidant Activity of Coumarin Type Compounds
Previous Article in Journal
Phenolic Compounds in the Potato and Its Byproducts: An Overview
Previous Article in Special Issue
Fast Modeling of Binding Affinities by Means of Superposing Significant Interaction Rules (SSIR) Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Oscar Martínez-Santiago
1,2,
Yovani Marrero-Ponce
1,3,4,5,*,
Stephen J. Barigye
1,6,
Huong Le Thi Thu
7,
F. Javier Torres
4,8,
Cesar H. Zambrano
4,8,
Jorge L. Muñiz Olite
9,
Maykel Cruz-Monteagudo
10,
Ricardo Vivas-Reyes
11,12,
Liliana Vázquez Infante
1,2 and
Luis M. Artiles Martínez
1
1
Computer-Aided Molecular “Biosilico” Discovery and Bioinformatic Research International Network (CAMD-BIR IN), Cumbayá-Tumbaco, Quito 170184, Ecuador
2
Department of Chemical Science, Faculty of Chemistry-Pharmacy, Universidad Central “Martha Abreu” de Las Villas, Santa Clara 54830, Villa Clara, Cuba
3
Universidad San Francisco de Quito (USFQ), Grupo de Medicina Molecular y Traslacional (MeM&T), Colegio de Ciencias de la Salud (COCSA), Escuela de Medicina, Edificio de Especialidades Médicas, Hospital de los Valles, Av. Interoceánica Km 12 ½—Cumbayá, Quito 170157, Ecuador
4
Universidad San Francisco de Quito (USFQ), Instituto de Simulación Computacional (ISC-USFQ), Diego de Robles y vía Interoceánica, Quito 170157, Ecuador
5
Grupo de Investigación Microbiología y Ambiente (GIMA), Programa de Bacteriología, Facultad Ciencias de la Salud, Universidad de San Buenaventura, Calle Real de Ternera, Cartagena de Indias, Bolívar 130010, Colombia
6
Departamento de Química, Universidade Federal de Lavras (UFLA), Caixa Postal 3037, Lavras 37200-000, MG, Brazil
7
School of Medicine and Pharmacy, Vietnam National University, Hanoi (VNU) 144 Xuan Thuy, Cau Giay, Hanoi 100000, Vietnam
8
Universidad San Francisco de Quito (USFQ), Grupo de Química Computacional y Teórica (QCT-USFQ), Departamento de Ingeniería Química, Diego de Robles y Vía Interoceánica, Quito 170157, Ecuador
9
Grupo de Investigación en Estudios Químicos y Biológicos, Facultad de Ciencias Básicas, Universidad Tecnológica de Bolívar (UTB), Parque Industrial y Tecnológico Carlos Vélez Pombo Km 1 vía Turbaco, Cartagena de Indias, Bolívar 130010, Colombia
10
Instituto de Investigaciones Biomédicas (IIB), Universidad de Las Américas (UDLA), Quito 170513, Ecuador
11
Grupo de Química Cuántica y Teórica, Facultad de Ciencias, Universidad de Cartagena, Cartagena de Indias, Bolívar 130001, Colombia
12
Grupo CipTec, Fundación Universitaria Tecnológico de Comfenalco, Facultad de Ingenierías, Programa de Ingeniería de Procesos, Cartagena de Indias, Bolívar 130001, Colombia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(6), 812; https://doi.org/10.3390/ijms17060812
Submission received: 31 January 2016 / Revised: 27 April 2016 / Accepted: 4 May 2016 / Published: 27 May 2016

Abstract

:
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.

Graphical Abstract

1. Introduction

Many mathematical invariants used in the codification of chemical information of molecular structures have in the recent years gained important utility in several research fields [1,2,3]. These invariants are more advantageous relative to physicochemical parameters customarily employed in describing, for instance, the hydrophobic, steric and/or electronic effects due modifications of substituents in a molecule (e.g., the Hammett’s sigma constant); because they are able to quantify greater chemical information on molecules and usually yield better performance in studies on structure-property/activity relationships, similarity/diversity, virtual screening, among others. These aforementioned mathematical invariants are known as Molecular Descriptors (MDs).
Formally, MDs may be defined as the final result of a logical and mathematical procedure in which the chemical information codified in a symbolic representation of the molecule is transformed into a number [4].
When MDs are based on the topology of the molecular structure, these are denominated as Topological Indices (TIs), which are derived from a graph-theoretical invariant and are able to codify information on molecular connectivity [5]. In other words, the TIs are numeric representations of the molecular structure and should encode useful structural characteristics such as: size, symmetry, ramifications, cycles, type of atoms, as well as the multiplicity of bonds in a chemical structure.
One of the 13 properties proposed by Randic as desirable for any MD is the direct structural interpretation [6]. Additionally, the fifth aspect taken into account to validate a Quantitative Structure-Active Relationship QSAR model [7] is related with its interpretation, which is directly related to the understanding of MDs in structural and/or physicochemical terms. However, the majority of the research efforts have been focused on applying these MDs (TIs) in a significant number of in silico molecular modeling tasks and the virtual screening of chemical compounds of interest. Efforts destined to the direct interpretation of the majority of the graph-theoretical invariants in structural and/or physicochemical terms are generally rare [8]. This fact is clearly portrayed in statement by Hoffmann that [9], “In many interesting areas of chemistry we are approaching predictability, but I would claim, not understanding”.
MDs will probably play an increasing role in computational, theoretical and medicinal chemistry. In fact, the availability of a large number of theoretical MDs containing diverse sources of chemical information would be useful to comprehend better the relationship between molecular structures and experimental evidence, taking advantage of the increasingly more powerful methods, computational algorithms and fast computers [4].
The previous affirmation demonstrates the increasing need to find some nexus among modern TIs with familiar concepts in established sciences, in order to optimize known procedures, group similar methods and create new MDs that codify orthogonal chemical information, and ultimately enhance their utility and scope of application.
Recently, some of the authors of the present work defined a new family of TIs based on the concept of the discrete derivative (specifically, derivative of a molecular graph), which was denominated Graph Derivative Indices (GDI) [10,11]. These GDIs have been applied in several Quantitative Structure/Property Relationship (QSAR/QSPR) studies showing satisfactory results [10,11,12]. These indices have been defined as global and local invariants for atoms and/or group of atoms and are based on the concept of discrete derivative of a Graph “G” with respect to an event “S” over duplex, triplex and quadruplex relations of atoms (vertexes) [10,12].
The concept of the graph-based derivative with respect to a given event “S” was proposed initially by V. A. Gorbátov in 1988 [13]. To evaluate the derivative of a group of elements belonging to a graph it is necessary to know the individual and the reciprocal participation frequencies of the graph elements (vertexes), in the set of conditions (sub-graphs), for which the event is true [10,11,13]. The derivative values characterize the non-homogeneous participation degree of the groups of elements from a graph in a given event [10,13].
The derivative value over n-elements of a G can be obtained as:
G S ( m 1 , m 2 , ... m n ) = 1 f m 1 , m 2 , ... m n ( i f i 2 i 1 , i 2 i 1 i 2 f i 1 i 2 + ... + ( 1 ) α + 1 α i 1 , i 2 ... i α i 1 i 2 , ... , i α 1 i α f i 1 i 2 ... i α + ... + ( 1 ) n + 1 n i 1 , i 2 , ... , i n i 1 i 2 , ... , i n 1 i n f i 1 i 2 ... i n )
If we wish to know the derivative over a pair (i,j) of elements from a graph (duplex), the Equation (1) may be simplified to:
G S ( i   j ) = f i 2 f i j + f j f i j
The Equation (1) (and Equation (2)) is the main expression for the underlying concept of the GDIs, thus a reliable interpretation of its elements is necessary to get an adequate understanding of the indices in terms of physical observables associated to the molecular structure.
In this report, a deduction of the graph discrete derivative, analogous with the classical derivative from mathematical analysis will be presented. This analysis will be used as basis for getting physicochemical and structural interpretation of the derivative values over atom-pairs. Besides, a group of carefully designed experiments to corroborate the veracity of the propositions made in the GDIs interpretation will be reported.

2. Graph Derivative Indices

2.1. Discrete Derivative. Deduction and Analogy with the Classical Derivative Concept

Isaac Newton is credited with the pioneer development of modern-day differential calculus, built on earlier work by Fermat, Barrow and Descartes. He defined the derivative of the function y = f ( x ) at the point x, as the limit of the differences relation as Δ x 0 :
dy dx = lim Δ x o Δy Δx = lim Δ x 0 f ( x + Δ x ) f ( x ) Δ x
Note that this notation was proposed by Leibnitz [14]. In the mathematical analysis the derivative characterizes a function’s variation degree when a small variation in its argument is made, that is, the classical definition of the derivative is based on the concept of the limit [15]. In discrete mathematics the definition of the limit does not hold due to the non-continuous nature of the functions, therefore it is impossible to establish an exact application of the derivative concept from continuous to discrete mathematics [13]. However to solve optimization problems in discrete mathematics, the discrete derivatives is introduced based on the use of the frequency of letters in words of a ψ model. In order to illustrate this concept, a graph G composed of four vertexes and five edges, as shown in Figure 1, is taken as an example. A new event “S” is defined (Figure 1B) so that it is possible to obtain fragments of G in “5” sub-graphs (words from model-collection of conditions for which the event is true).
It is feasible to characterize the participation frequency of conditions (letters, atoms) in the collections of conditions (words, sub-structures) in which the event is true using the corresponding set of conditions. For vertexes “a” and “b” the corresponding frequencies would be fa = 2, fb = 3. With this scheme is also possible to determine the simultaneous inclusion’s frequency for vertexes “a” and “b”, i.e., the Reciprocal Frequency (fab). For the example from Figure 1, the fab = 1. Based on latter analysis, we may conclude that for any pair of vertexes (letters) i & j from a G, the individual frequencies are greater or equal to the reciprocal frequencies between them, that is f i f i j and f j f i j [10,11,13].
Up to now two sub-graph sets are defined, namely Ma and Mb, obtained by graph fragmentation, taking as base the application of a given event. Both sets may or may not have coincident elements (sub-graphs that contain “a” and “b” vertexes at the same time). The number of coincident sub-graphs is quantified by the reciprocal frequency of vertexes “a” and “b” (fab) and expresses the magnitude of the separation of these sets. Therefore, fab is similar to the increase or perturbation ∆x, defined in classical mathematics.
Maintaining the analogy with mathematical analysis, it is possible to evaluate de variation’s degree as the symmetrical difference between sets, as shown in Equation 4:
M a Δ M b = ( M a \ M b ) ( M b \ M a )
If this expression is written taking as base the individual and reciprocal frequencies, it would be:
M a Δ M b = ( f a f a b ) + ( f b f a b )
M a Δ M b = f a 2 f a b + f b
Dividing Equation (6) by fab, which is representative of the perturbation degree from a set to another, and expressed in terms of the differentiation with respect to the event, Equation (7) is obtained:
d G d S ( a   b ) = f a 2 f a b + f b f a b
The mathematical Expression (7) is denominated as the graph derivative with respect to an event over a pair of vertexes (duplex relations), and it defines the non-uniform participation degree in the “S” event of the pair (a,b) belonging to G. In other words, Equation (7) expresses the heterogeneity degree of pairs of components constituting G with respect to any previous event [10,11,13] and it can be interpreted as a non-directed weighed graph < V , ( U , P ) > whose bearer corresponds with that of a model determined by this event and vertexes (vi, vj), weighted by the ratio of the incompatible frequency [(fi − fij) + (fj − fij)] with the compatible frequency fij in the event, and with a particularity that [13]:
( v i , v j ) U ,  if  dG dS ( v i , v j ) =  
( v i , v j ) U ,  if  dG dS ( v i , v j ) = finite magnitude different of zero
( v i = v j ) ,  if  dG dS ( v i , v j ) = 0
A scheme that summarizes the analogy between discrete and continuous derivatives is shown in Figure 2. Note that various events may be applied to a given model; thus allowing for varied information to be retrieved from the chemical structure, and ultimately yielding different MDs [16]. Recently 12 theoretical events were introduced and these include: connected sub-graphs (S), walks of length k (K), Sach’s sub-graph (H), Quantum (Q), Terminal Path (T), path-vertex incidence (V), Multiplicity (M), MACC fingerprint (C), Sub-structures fingerprint (B), E-State fingerprint (E), Alog P (A), and Refractivity (R). The application of these events has allowed for the obtaining of varied information and generating a high number of MDs that characterize the chemical structure from dissimilar perspectives. These events have been successfully applied in several applications and are grouped in three clusters: Topological events (S, K, H, Q, T, V and M), Fingerprint-based events (C, B and E) and Physico-chemical events (A and R) [3,16].
It is also important to analyze the characteristics of the Relations Frequency Matrix (F), given that its elements make possible the evaluation of the discrete derivatives. The matrix F is an n × n symmetrical matrix (where, n is the number of atoms in the molecule), whose diagonal elements are denominated as Individual Frequencies (fi) of each element in G, and the off-diagonal elements correspond to the Reciprocal Frequencies.
To determine the derivatives over n-elements as shown in the Equation (1), F would be an n-dimensional matrix (or so-called hypermatrix). In a previous report, the GDIs have been generalized for calculating derivatives over n-tuples of atoms [12], and a similar theoretical scaffold was also applied in generalizing the GT-STAF (acronym for Graph Theoretical Thermodynamic STAte Functions) indices based on information theory [17], From this point onwards the attention will be focused on the graph derivative over duplex relations, which will be calculated when the molecule is fragmented according to an event criteria S and using order 1 substructures. This easy description of the molecular structure allows arriving at conclusions which will serve as a base for the interpretation of the values of the discrete derivative over a pair of atoms and posteriorly generalized to n-tuples of atoms.

2.2. Graph Derivative Indices (GDI). Application to Chemical Codification

The main goal of this manuscript is to find the structural and/or physicochemical interpretation of GDIs. Although GDIs have been successfully defined and applied in several studies [10], it is necessary to remember some medullar aspects of their theory and as well as the corresponding mathematical algorithms used in the description of the organic structure.
Let’s apply the previously discussed aspects in an example using the molecular structure of 2-amino-5-vinylfurane (Figure 3).
The event connected sub-graphs (S) considers the formation of the molecular structure (G) taking as base molecular fragments (sub-graphs with different orders and types according to Kier-Hall nomenclature, which is Path, Cluster, Path-Cluster and Chain) [10,11,16], By using the scheme of the Figure 2 and applying the event (S) for the order 1 (pairs of connected atoms or individual bonds), the following relations frequency matrix is obtained:
F = 1 1 0 0 0 0 0 0 1 3 1 0 0 1 0 0 0 1 2 1 0 0 0 0 0 0 1 2 1 0 0 0 0 0 0 1 3 1 1 0 0 1 0 0 1 2 0 0 0 0 0 0 1 0 2 1 0 0 0 0 0 0 1 1
For evaluating which bonds contribute in a greater measure, from a topological point of view, to the formation of the structure, the derivative over the pairs of atoms is calculated. For our previous example, the derivative values are:
G S ( 1 , 2 ) = 1 2 ( 1 ) + 3 1 = 2.00   G S ( 2 , 3 ) = 3 2 ( 1 ) + 2 1 = 3.00
Analogically, for the rest of the connected pairs their derivative values are:
G S ( 2 , 6 ) = 3.00 ;   G S ( 3 , 4 ) = 2.00 ;   G S ( 4 , 5 ) = 3.00 ;
G S ( 5 , 6 ) = 3.00 ; G S ( 5 , 7 ) = 3.00 ;   G S ( 7 , 8 ) = 1.00  
Analyzing these derivative values, it is observed that the most influent pairs according to the chosen event (i.e., the most contributing bonds in the formation of the molecular structure) are 2-3, 2-6, 4-5, 5-6 and 5-7, respectively. This is a logical result because the atoms 2 and 5 present the highest number of connections in the molecular graph G (Figure 3B). The GDI values obtained for all atom pairs can be organized as a matrix D = [ G S ( i , j ) ] n × n , where it is possible to obtain individual atomic local indices by adding all the values of the derivative for the atom i, which is equivalent to adding all the elements from each row or column from the D matrix. This atomic index Δ i = i = 1 n G S ( i , j ) is a Local Vertexes Invariant (LOVI) [4,11,12,18,19]. The atomic indices Δi for a given molecular structure with n atoms may be expressed as a LOVIs vector, (VL). In this sense, the LOVIs vector for the molecule of Figure 3 is: V L = [ 2 8 5 5 9 6 4 1 ] 1 × 8 .
However, VL only takes into account information on the connectivity’s degree of the different atoms in the structure without considering the type of bonded atom. It is important to apply a weighting scheme to the atoms in a compound to yield a description much closer to the molecular structure reality. In the GDI this information is introduced through three different schemes, as it will be explained below.

2.2.1. Atomic Differentiation

The matrix treatment of graphs does not guarantee the chemical differentiation of the elements in the G; thus, methods for differentiating the atoms of different nature in the molecule have been introduced. As a preliminary step, an atomic weight is assigned to each atom, through the expression:
ϑ i = P i δ i
where, δ i is the vertex degree for atom i in the molecular structure and Pi is a property that characterizes each atom. With this definition a weight vector V w = [ ϑ i ] 1 × n may be constructed, which contains information on each type of atom present in the molecule and positioned within a particular electronic environment. For practical purposes, the resulting values may be written as a row vector, column vector or as a W diagonal matrix, which is applied in three different schemes [10].
Weighting in the incidence matrix: A weighted incidence matrix may be obtained as a result of the multiplication of the Q incidence matrix with the weighted matrix W ( Q   ×   W = Q w ). Subsequently, the graph’s derivation process is performed as previously described. A same result is obtained if we introduce the direct weighting of each atom (and/or atom-pairs) in the frequency matrix and ultimately yielding a weighted LOVIs vector (wVL-f). For the molecule in the Figure 3, the weighted LOVIs vector, using Pauling’s electronegativity as atoms weighting scheme, is V L f w = [ 4.92 ,   11.29 ,   9.55 ,   10.61 ,   4.92 ,   3.60 ,   3.75 ,   0.83 ] .
Weighting in the derivative matrix: Once we have the D matrix, it may be multiplied with V w vector and a new weighted LOVIs vector (wVL-d) is obtained. For molecular structure of 2-amino-5-vinylfurane, this vector would be V L d w = [ 3.61 ,   14.40 ,   3.82 ,   10.57 ,   3.61 ,   1.27 ,   3.19 ,   0.85 ]   .
Weighting in the LOVI vector: The product of the multiplication of VL × W yields a new vector and its component also represent weighted LOVIs   ( V L w = V L   ×   W ). For the molecule from Figure 3 V L w = [ 4.25 ,   5.10 ,   10.95 ,   5.74 ,   4.25 ,   6.38 ,   3.40 ,   1.27 ] .
The definition of local atomic descriptors, following the previous prescription, is in harmony with one of the 13 properties proposed by Randic that new TIs should possess [4,6].

2.2.2. Total and/or Local (Group or Atom-Type) Description

In analogy with the Molecular Orbitals Theory, which states that molecular orbitals can be obtained as a linear combination of atomic orbitals [20,21,22], MDs associated with the whole molecular structure (or to a group of atoms in the molecule) can be calculated through the mathematical invariants shown in Table 1.
In Table 1, the invariants (aggregation operators) are classified in four groups: (1) norms; (2) means; (3) statistical parameters and (4) “classic” invariants. Let’s compute some total and local invariants (unsaturated bonds and heteroatoms) for the 2-amino-5-vinylfurane molecule using the LOVIs of the unweighted atoms as an example:
Total Invariants: N1 = 40 (N1 means Manhattan norm (sum of atom-level descriptors (LOVIs), A (Harmonic Mean) = 5 and R (Range) = 8 (see Table 1 for more details).
Local Invariants over unsaturated bonds and heteroatoms: N1 (Unsaturated bonds) = 32 and N2 (heteroatoms) = 6.3245.
With this procedure, a great number of possibilities in the description of important chemical information on the molecular structure may be explored (see Table 1 for more details). In previous papers, all these invariants were introduced and applied to obtain total and local indices over groups and/or atom types, generalizing the traditional way for obtaining the global indices [10,11,12]. Other families of indices recently defined by our research group (Computer-Aided Molecular “Biosilico” Discovery and Bioinformatic Research (CAMD-BIR) International Network) also make use of these invariants and have been applied to many structure-activity/property relationship experiments with relevant results [3,17,23,24].

3. Structural and Physicochemical Interpretation of GDIs

The GDIs have been previously used in several theoretical applications, providing relevant results fundamentally in QSAR/QSPR studies. [10,11,12] However, little effort has been destined to the interpretation of these MDs in structural and/or physicochemical terms [10,11,12,16]. In this section, the results of different experiments designed to demonstrate the reliability of GDIs in structural and physicochemical hypotheses are presented. The GDIs interpretation will be developed using three different approaches, which will be detailed for every case study.
Regarding the computational details, the calculation of GDIs for all the databases was performed using the open-source, java-based software denominated as DIVATI (acronym for DIscrete DeriVAtive Type Indices), a new module of ToMoCoMD (Topological Molecular Computational Design) software [25].
QSAR/QSPR models were developed by using Multiple Linear Regression (MLR) with the MobyDigs software [26]. This program allows obtaining MLR equations using the genetic algorithm [27] as MD selection method. Additionally, it allows for several model validation procedures such as: internal cross validation (Q2loo), external validation (Q2ext), bootstrapping (Q2boot), and Y-randomization, as well as prediction analysis.

3.1. Structural Interpretation

Some of the 13 characteristics, proposed by M. Randic, that new TIs should ideally possess include: possible structural interpretation, isomers recognition, possibility for local definition, as well as present a correct dependence with size and gradual change with structural variations, among others [28]. These properties are closely related and all of them indicate the direct correspondence that should exist among topological indices calculated for a molecule and its structure. If this condition is achieved the MDs can, at least in principle, describe any chemical information extracted from the connectivity of the molecular structure.
All calculations for the structural interpretation experiments performed in this section were made for derivatives over pairs of atoms with respect to the connected sub-graphs event (each event yields different LOVI values, but the interrelations among them are similar), using the generalized matrix and Pauling’s electronegativity over the incidence matrix as the weighting scheme The LOVIs for each atom in the molecular structures employed in the present section, as well as the corresponding Minkowsky norm for p = 1(equivalent to the summation operator), the arithmetic mean, and the range (see Table 1) are depicted in Table 2, Table 3, Table 4 and Table 5.

3.1.1. Differentiation among Homologs, Multiple Bond Codification and Positions

From Table 2, it can be inferred that LOVIs for all the different atoms decrease as one get closer to the center of the carbon chain, belonging to the considered linear alkanes and their homologs (for a sake of simplicity, the hydrogen atoms are no taken into account in the representation of the G).
The value of every atom of the chain also increases as one moves from an inferior homolog to a superior one, maintaining the regularity among the LOVIs of atoms in the same molecule. Besides, the total invariants shown in this Table reflect a regular increase in their respective values when the number of methylene groups –CH2– in the structure increases from one homolog to another. The increase of the electronic density and its position in a specific region can be also quantified using this mathematical description. It is observed that atoms connected by multiple bonds also increase their LOVI values when the number of bonds among them increases. Consequently, the total (global) invariants from this molecule increase when the number of atoms increases. As illustrated these total invariants can be used for differentiating position isomers because their values decrease in the measure that the multiple bond is more embedded in the molecular structure.

3.1.2. Differentiation among Chain Isomers

In Table 3, are the results of an experiment analogous to the previous one but now using chain isomers of saturated hydrocarbons with six carbon atoms. Table 2 shows how the variations in the carbon chain from the molecule cause variations in the LOVIs values from each atom, according to the modifications experimented in the molecular structure, from one isomer to another.
It is interesting to note the peculiar variation of the different invariants when the quantities of ramifications increase, and consequently the length of the main chain decreases. Note that the invariants decrease their values, which is a logical result since LOVIs also experiment a decrease in their values in the measure that molecules are more ramified.

3.1.3. Codification of the Presence of Different Functions

Until this point, only molecules comprised exclusively of carbon and hydrogen have been described. Here, the capability of GDIs for describing structures containing different organic functions is evaluated.
The organic functions play a fundamental role in the properties of bio-active molecules; therefore their description cannot be ignored. In Table 4 the results of the calculation of the atomic and total indices for a series of propane derivatives is summarized. As in the previous cases, the electronegativity according to Pauling’s scale is employed as label for differentiating the atoms in the molecule.
It is observed that when the electronegativity of atoms located in “X” position increases, the LOVIs values also increase and the total invariants increase in almost a regular way. It is interesting to detail that the increase observed in LOVI values, from a molecule to another, is proportional to the increase in the property used for characterizing the atoms. Besides, it should be pointed out that this increase in the LOVI values due to the presence of different X groups is identical to the summation of the increase of the rest of the chain.

3.1.4. Cycle, Conjugation and Aromaticity

Table 5 shows a set of cyclic structures and some few acyclic chains. The purpose of choosing these systems is to determine if GDIs are able to retrieve information about the conjugation, aromaticity and the presence of cycles, during the numerical description of the molecular structure.
From a detailed analysis of data in Table 5, it can be observed that LOVIs from a cyclic structure are smaller in comparison to the values obtained for the corresponding acyclic structures with an isolated unsaturated bond. In the measure that the number of unsaturated bonds increase, the LOVIs values decrease, but these decrements are more significant when electronic conjugation is present in the structure, meaning that the aromatic molecules are characterized by the lowest values. The presence of heteroatoms in the systems increases the LOVIs values for carbon atoms compared to counterparts without heteroatoms. Moreover, atoms bound to heteroatoms present a smaller index value than the other ones, which is an indication that these atoms have a lower electronic density provoked by the direct interaction with a more electronegative atom. This result is consistent with the chemical reality of heteroatomic systems.

3.1.5. Molecular Symmetry Described in Terms of GDIs

A more complete structural description is obtained when the symmetry of the molecules is incorporated in the description achieved for the GDIs as previously shown.
In this sense, Shannon’s Entropy is used as a symmetry index (see Equation (9)) [4]:
S I = i = 1 G ρ ( g i ) N l o g 2 ρ ( g i ) N l o g 2 N
where ρ ( g i ) is the cardinality of the equivalent classes, which in our particular case is the number of equivalent LOVIs and N is the total number of atoms in the molecule. The symmetry index may be interpreted as a measure of the bidimensional symmetry of an organic molecule, and it takes values between 0 and 1 ( 0 S I 1 ). A value of 0 corresponds to molecules with high 2D symmetry, whereas a value of 1 is associated with asymmetric molecules.
Table 6 shows the SI values computed for ten organic molecules, taking as base the cardinality of the LOVIs obtained with GDIs for derivatives over pairs of atoms with respect to the connected sub-graph event and by weighting the incidence matrix with Pauling’s electronegativity.
From the analysis of Table 6, it can be inferred that the resulting SI are logical and acceptable if it is considered that this index contains information on the symmetry of the molecules based on the local vertex invariants calculated for these 2D structures.

3.2. Reactivity Based on Geometrical Effects

Consider a molecule fragmented with respect to the connected sub-graph (S) for order 1. Note that the individual frequencies are the number of total connections for each atom in the model generated by the considered event. For this particular case (S event for order 1), the individual frequencies correspond to the number of sigma bonds for each atom. In electronic terms this frequency may be written as:
f i 1 = V i t H i π i n i
f i 1 = V i H S π i n i = σ i
where, V i t , H i ,   V i H S , π i ,    n i , and σ i are the total valence, the number of hydrogen atoms bond, the valence number with suppressed hydrogen atoms, the quantity of π bonds, the number of lone pairs associated to atom i and the number of sigma bonds from the atom i, respectively.
According to the latter expression, it can be observed that the individual frequency of order 1 can be simply defined as the number of sigma bonds of atom i. In the case where the generalized frequency matrix (including kth-order sub-graphs, k > 1) is employed, the interpretation of the individual frequency holds but should be identified as the number of connections between atoms separated by distance k, depending on the model generated by the chosen event and the maximum fragment order.
In this sense, the derivative evaluated for frequencies of order 1, over a pair of bonding atoms may be written as:
d G d S ( i   j ) = ( V i H S + V j H S 2 ) 2 π n i n j
d G d S ( i   j ) = A 2 π n i n j
For carbon atoms in an organic molecule Equation (13) can be simplified to:
d G d S ( i   j ) = A 2 π
By analyzing Equations (12)–(14), it is evident that the derivative codifies information on the total number of sigma type connections associated to a pair of atoms (i j) connected at a specific distance. In this sense, the derivative can be defined as measure of the “accessibility” of a given region of a molecule.
The LOVI for each atom i can be defined as the multiplicity of its sigma bonds affected or influenced by the respective surrounding multiplicities. This definition can be mathematically written as:
Δ i = n σ i + j = 1 n σ j 2 n
The term 2n appears from the expression of the derivative and it has the effect of cancelling out the elements repeated in both groups of bonds.
Based on these concepts, the LOVI value of an atom can be regarded as the degree of participation of this atom in the formation of the molecular structure; therefore, the accessibility ε i to an atom can be evaluated as the inverse of its LOVI because it is proportional to the total area in which an atom is likely to be attacked by any external entity. In this sense, the interaction capacity of two atoms located in different molecules can be evaluated as the product of their accessibilities ( ρ i j = ε i ε j ) , defining the interaction capacity of an atom i with any other atom j as the addition of accessibility contributions of the atom i with the rest ( I i = j = 1 n ρ i j ) . Generalizing this idea, the interaction capacity of two molecules ( Ω ), understood as the total contribution of all the individual capacities of the constituent atoms determined from their accessibilities, should be proportional to the interaction probability of two molecules, which can be mathematically expressed as:
P ~ Ω = M ^ ( j = 1 n ρ i j )
where, M ^ is an operator, which involves the operation of the individual interaction capacity from each atom (M would be the sum, product, etc.), P is the real probability of interaction between two molecules and Ω is the theoretical probability evaluated by taking as base the LOVIs from each atom. The real value of the intermolecular interaction probability P depends on the structures (symmetry, form, size and distribution). These structural parameters are carefully quantified by GDIs as was showed in the previous epigraph.
In Collisions Kinetic Theory, the magnitude known as the Steric Factor (F), is considered as a consequence of the analysis of incongruences between theory and the experimentation results and it expresses the probability of the corresponding geometrical configuration during interactions [29]. Posterior deductions showed that the Steric Factor from Collisions Kinetic Theory is analogous to the e Δ S / R term, which is related with the probability of existence of a viable configuration able to interact with another specie [29] and Δ S is the variation of entropy of the active complex and R a constant ( R = 8.314  J · mol 1 · K 1 ) .
The LOVIs, their inverses, and the parameters introduced in both kinetic theories, have a direct relationship with the probability of intermolecular interaction based on structural configurations during the interaction. In this sense, it is possible find analogies between these kinetic parameters and the local and total indices derived from the corresponding mathematical procedures implicit in GDIs calculation.
To evaluate the veracity of the aforementioned ideas, second order dimerization reactions, in gaseous states were studied as an example. Table 7 summarizes the ΔS# data for each studied reaction [30].
In each case, the total interaction capacity Ω was evaluated and it was related with the activation entropy for the process of dimerization of diverse unsaturated molecules, yielding a correlation of 97.58%. The operator taken to evaluate the totality of the atomic contributions was the standard deviation, which is a measure of the dispersion of the values.
The activation entropy is a parameter associated with the probability of an adequate interaction between two spaces, taking as reference the structural organization of the created transition state, which is product of a suitable interaction. In this sense, it is logical to find a strong relation between parameters obtained from GDIs calculation and this magnitude. This example can be taken as base for understanding the molecular graph derivative from a kinetic point of view, which is closely related with the possibility of interaction of the molecules based on their corresponding geometrical structures.

3.2.1. Accessibility as Measure of the Interaction’s Capacity

Accessibility is defined as the possibility of entrance or access capacity, in this case, to a specific region in a molecule. The Kier and Hall [31,32] molecular connectivity indices have been recently studied by Estrada in terms of the Relative Bond Accessibility Area (RBA), Cij, which is expressed in square Randics (R2) [33,34]. The reactivity of an atom is related with its accessibility, and the accessibility can by quantified from topological features of the molecular structure that contains this atom [35,36].
As was mentioned in the previous section, the inverse of a LOVI may be considered as a measure of the availability of a particular zone of a molecule (fully characterized by a pair of atoms) to interact with an external agent. From this definition, it can be stated that this MD is closely related to the steric effects controlling molecular reactivity or simply steric reactivity. The steric reactivity associated to a pair of bound atoms can be evaluated as ( 1 Δ i Δ Δ j ) , where Δ i is the LOVI of atom i. In Table 8, steric reactivity indices computed for a collection of molecules are summarized together with RBA in R2 calculated by Estrada for these molecules [33,34].
Both sets of data from Table 8 are plotted in Figure 4, where a qualitative agreement between the two curves may be observed. Note that for the pair of atoms (a,b) in the molecules 4 and 5 (2,2-dimethylpropane and butane, respectively) the same RBA values were obtained using the connectivity index, while the evaluation of the steric reactivity with GDIs calculations yielded better differentiation more consistent with the chemical reality because it is anticipated that the a–b bond in the 2,2-dimetylpropane molecule would be less accessible due to greater steric hindrance than the a–b bond in the butane molecule. This challenge of index degeneration in situations like these is not encountered with GDIs.
Another important aspect is the infinite value obtained for the trivial case of the ethane molecule. Although this value is useless for further statistic or algebraic developments, it is logical if we take into account that it is evaluating the accessibility of an entity external to a particular bond in each molecule. With this idea it plausible that the ethane molecule has infinite accessibility possibilities given that it is constituted by only one bond (G with suppressed H-atoms were considered in all cases). To visualize this effect in Figure 4, an R-GDI value equal to 1.5 was arbitrarily assigned to this molecule to give the idea that it is superior and thus allowing for the visualization of the regularities of the rest without affecting the scale. Note, however, that this does not mean that GDIs would not have the capacity of codifying the structure of one molecule as simple as ethane, because this infinite value is only given for the order 1, which is the configuration used for simplifying the physical interpretation of the GDIs that may later be generalized to more complex systems (i.e., using sub-graphs of superior orders in the event S, derivatives over n-elements, using other events, and so on).

3.2.2. Specific Reaction Rate and Its Relation with GDI

The specific reaction rate constant is the speed of chemical transformation when the concentration of all the reagents is equal to 1 mol/L [29,30]. Thus, the rate constant is a magnitude proportional to the reactivity of the molecules that participate in a reaction, with the temperature of the system being constant [29].
In the previous sections the derivative over a pair of bonded atoms was defined as the quantity of sigma connections of both atoms and was expressed as a difference of electronic contributions that express the part of the total electronic density destined to connections or sigma bonds. In this sense, the derivative over a pair of atoms i and j can also be understood as a measure of the interaction capacity of a bond with a neighboring molecule.
To evaluate the previous affirmation a data comprised of 34 derivatives of 2-vinylfuranes was employed, for which the specific rate constant for the nucleophylic addition to the double exocyclic bond from this molecules with mercaptoacetic acid has been reported [2]. The best one variable regression shows that there exists a moderate relation between the property and the GDIs. However, taking into account that the best models reported by other authors for describing this property use seven variables [2,37], it is noteworthy that with only one descriptor (derivative over the double exocyclic bond) it is possible to explain approximately 80% of the variance of the experimental property. Figure 5 shows the correspondence between the derivative values and the experimental values from logK.

3.3. Interpretation of GDIs in Electronics Terms

The nature of the atoms and the molecules is determined by their electronic structure [21], thus, by describing the dynamics, distribution and the electronic energy of these systems, it is feasible to establish a useful nexus for the better understanding of molecular structures and/or the methods used to codify them.
The detailed structural analysis in the previous epigraph showed that there is a relationship between GDIs and the influence of the electronic richness of atoms bound to a specific center in a molecule. Moreover, the obtained result demonstrated the possibility of a physicochemical interpretation in kinetic terms, on the basis of the quantification of the interaction potential of a molecule described by means of GDI and LOVIs obtained in geometrical framework. All these aspects are easily assimilable if we take into account the individual frequencies and the derivatives over atom pairs expressed in electronic terms as explained in the previous section. Figure 6 illustrates the electronic decomposition of the individual frequencies of order 1 and the derivatives of order 1 for a pair of bonded atoms.
Up to this point, only the quantity and distribution of bonds in the molecular structure has been taken into account for the GDIs interpretation. However, the knowledge about these bonds and their distribution around each vertex determines the electronic density in the environment of each atomic nucleus in the molecular structure, and thus motivating a deeper analysis of the GDIs from a Quantum Mechanics perspective.
From Figure 6, the orbital description of the frequencies allows for two observations: (1) if the frequency is equal to the number of sigma bonds of an atom, it then also quantifies the hybridization of the atom and the symmetry of the sigma electronic distribution and (2) the parameters to the left suggest a separation of the electronic terms corresponding to the quantification of σ, π and non-shared electrons, which allows for the evaluation of the separated influence of these types of electrons, as it will be corroborated in the next part of the present study regarding the relation of the frequencies with Hückel’s free valence. It can be pointed out that only the frequency of order 1 is similar to the number of sigma bonds belonging to an atom. Nonetheless, it is anticipated that a similar concept holds in an abstract model based higher order frequencies obtained according to different fragmentation approaches.
For the case of the derivative, it can by described as a quantity related to the part of the total valence of the bonding atoms, which is destined to the establishment of connections in a molecular network. Equally, the separation of the π electrons offers a differentiated treatment to each type of electrons (it makes reference to the electrons that form part of the covalent bonds, π- or σ-type and non-shared electrons), considering the basic differences among their characteristics and expressing one in function of others, which shows the relation among them.
The fact that both the frequency and the derivative can be can be expressed in terms of the distribution of electrons around one atom or bond approximates to the idea of chemical reactivity, now understood from the perspective of the electronic interaction capacity. In next sections, carefully designed experiments will show the strength of these basic local descriptors from the discrete derivative algorithm in describing the electronic characteristics of atoms and molecules.

3.3.1. GDIs, Chemical Reactivity and Relation with Hückel’s Free Valence

A possible approximation for studying the chemical reactivity is to determine the degree in which the atoms in a molecule are united to the adjacent atoms, which is relative to their theoretical bonding capacity [38], The degree in which one atom is united to its neighbors can be calculated summing all the values of the bond orders of that atom. If the sum of all the bond orders is subtracted from the value of the highest bonding capacity, we would obtain the free valence:
F r = Highest bonding capacity j = 1 n ρ i j
For conjugated hydrocarbons, the free valence index can be evaluated as:
F r = 4.732 j = 1 n ρ i j
where, ρ i j is the bond order for the atoms (i j), determined by Hückel molecular orbitals method [38,39].
The individual frequency of each atom evaluates the number of connections of an atom in a specific model generated by an adequate event for describing the molecular structure. In the particular case of the event S, if only the order 1 matrix is taken into account, it is possible to relate this frequency with the quantity of sigma bonds of an atom.
Reorganizing the terms for the individual frequency shown in the Figure 6, we obtain:
V i T f i 1 = H i + π i
The term to the left of the Equation (19) describes the part of the total valence dedicated to forming π type bonds. The right term shows the number of π reactive electrons and the number of H-atoms bonded to each carbon. Both members from that equation show huge similarities with the Equations (17) and (18), although from different optics but explaining the same characteristic: chemical reactivity of each carbon atom in the molecular system, therefore, it should be hoped that:
F r ~ ( 4 f i ) = ( H i + π i )
Taking into account the equation 19, a study was performed on the behavior of the free valence calculated by the Hückel molecular orbitals method for specific atoms with different environments in a dataset composed of 19 conjugated hydrocarbons, with some of them being aromatic. This chemical dataset was primarily used by Kier and Hall [39] with the goal of finding a relationship between the free valence and the Electropologic State index, which is a descriptor regularly employed in correlations with electronic properties of different molecular systems and known to yield good results [39,40].
Figure 7 shows the free valence calculated by the Hückel’s molecular orbitals method (blue line), the Electropologic State index (red line) and the derivative indices taking as base the individual frequencies of the atoms. A uniform variation between the values obtained for Hückel’s free valence and the GDIs is observed and this evidently shows that the GDIs quantify electronic environments of atoms in the molecules and reinforcing the understanding of the GDIs in reactivity terms. In previous studies, Kier and Hall grouped the atom types in these molecules in three categories: carbon atoms with one hydrogen, carbon atoms with two hydrogens and carbon atoms without hydrogen [39]. Likewise, a linear relationship is found in the present experiment:
F r 1 f i 1 + C t
where, C t is the proportionality constant which adopts different values for each type of atom from the system and it allows finding reactivity values closer to the values evaluated by Huckel’s method. Table 9 shows the six types of atoms that can appear in conjugated systems conformed by only one carbon and hydrogen.
The Equation (21) shows a simple relation to evaluate Hückel’s free valence taking as base the order 1 frequencies with the S event. It is obvious that in the Expression (21), the goal is not the substitution of the real calculation of Hückel’s free valence, but the use of this expression allows establishing in a faster and simple way an estimate of the chemical reactivity of any carbon atom in a conjugated molecule without evaluating the bond order, for which it is necessary to know the coefficients of each atomic wave function that participates in the linear combination to form the molecular wave function.
Another advantage of the Equations (20) and (21) is that they demonstrate a clear existence of a nexus between the electronic environment of the atoms and their corresponding atomic indices evaluated with derivative indices, reaffirming the notion that graph derivative indices describe electronic properties of atoms in a molecule and consequently their electro-chemical reactivity.

3.3.2. Electronic Interpretation

In the beginnings of 1950s, it was discovered that the resonance frequency of a nuclide depends not only of its magnetogyric ratio and the intensity of the magnetic field B0, but that it also depends of the electronic environment where the nuclide is located [41]. For one nuclide in a specific substance, there will be as many resonance frequencies as the electronic environments. This phenomenon known as chemical displacement (shift), is the base of the NMR (Nuclear Magnetic Resonance) [22,41] chemical applications. The chemical shift is a descriptor of the electronic characteristics of each atom from a molecule. In the following experiments we intend to find linear relationships between the chemical shift of some active nuclides in NMR, with the goal of discovering in what measure the electronic information of atoms and molecules is codified in the structural descriptions based on the GDI concept.

QSPRs of Chemical Shift of 17O-NMR for Aldehydes and Ketones

For this analysis a data of aldehydes and ketones was used, which have been previously studied by Kier and Hall [40] with the Electropologic State index. All molecules are aliphatic and their chemical shifts of the 17O have been reported in the literature (see Table 10).
The variables used for linear regression were obtained by calculating the duplex derivative with respect to 10 different events (using chemical, physical and graph-theoretical atom-labels) and several norms, means, statistical and classic invariants as total and local MDs. For this experiment the MobyDigs software was used [26]. Figure 8A shows the performance of one variable regressions built for each of the events in the present study, according to the cross validation parameter Q2Loo. As can be observed, the events with the best correlations for the studied property are Sub-Structure (B) and Multiplicity (M), respectively. This observation is logical considering the fact that first one is a fingerprint-based event [3,16], and it conforms the incidence matrix only with substructures with functional groups and/or atom types of chemical interest, while the second one is a topological description of the connections at one step topological distances and their multiplicities (simple, double and triple bonds between pairs of atoms in a molecule) [3,16]. These two events yield matrices that reflect the electronic richness of the molecule, fragmented in individual bonds and their multiplicities. The best regression (based on the sub-structure event) obtained in the present study is shown below (see Equation (22)):
δ = 592.94 ( ± 1.58 ) 8.62 ( ± 0.38 ) [ A N / I n ( I S ) ] B D f
R2 = 98.65%, Q2Loo = 98.1%, Q2Boot = 98.21%, s = 2.304, scv = 2.41, ysc = −0.051, F = 512.25.
As it can be observed from the statistics of this equation, a good correlation (R2 = 98.65% and Q2Loo = 98.10%) is obtained between the calculated MDs and the experimental chemical shift for 17O. An analysis of the descriptor contained in this model shows that it expresses arithmetic mean of the LOVIs of the carbon and oxygen from carbonyl group (IS: unsaturated bond), weighted by the valence degree, which is a topological expression of each atom [4]. This is a logical result taking into consideration that the value of the chemical shift of oxygen mainly depends on the electronic environment of this atom as well as of the influence of the electronic density of its unique adjacent atom (carbonyl carbon); therefore, the average is a direct quantitative measure of the electronic richness in the model, which is the main factor that influences the chemical shift of the nuclide of 17O.

QSPRs of Chemical Shift of 17O-NMR for Eithers

A similar study was performed using a dataset composed of 10 aliphatic ethers (Table 11), for which the chemical shift of 17O-NMR was reported [40]. The best one-variable model obtained in this study is shown below (Equation (23)):
δ = 105.54 ( ± 2.82 ) 3.64 ( ± 0.109 ) [ A T / I n ( H T ) ] M D f
R2 = 99.28%, Q2Loo = 98.94%, Q2Boot = 99.05%, s = 3.588, scv = 3.891, ysc = −0.046, F = 1102.54.
As can be seen, this model is statistically robust (R2 = 99.28% and Q2 = 98.94%). It is especially interesting that the variable that best correlates with the modeled property is the one that uses the topological polar surface area (T) to weight the oxygen LOVI. This property is an expression of the electronic and steric environment of the oxygen nucleus and the LOVI explains the influence of the groups adjacent to oxygen. Therefore, it is logical that this local atomic index has greater influence in modeling the studied property.
Figure 8B shows the performance of one variable models built for the events in the present study, according to the cross validation parameter Q2Loo. As can be observed, the best events for describing the electronic properties of atoms in these molecules are: multiplicity (M) and sub-structure events (B), respectively. Good performance is also achieved with connected sub-graphs (S) and Sach’s sub-graphs (H).
Table 11 shows the values of the variable that best correlates with the mentioned property. In the measure that chemical shift values increase, the LOVIs values decrease in a regular way, having a negative contribution in the model.
This linear relationship between the LOVIs from the oxygen atom and its chemical shift in NMR implicates that the LOVIs from the oxygen atom codify electronic information varying in an almost uniform way with the chemical shift values.

Comparison between Real Spectrums and Atomic GDIs for Alkanes

To comprehend better the relationship of the GDIs with the electronic properties of the atoms of a molecule, some simple alkane molecules have been selected and their protonic Nuclear Magnetic Resonance spectrum (1H-NMR) predicted. Posteriorly, these are superimposed with the LOVIs values computed for the atoms in each molecule. The spectra were predicted using the ChemDraw program [42] and in both examples a good estimation was achieved. The selected molecules were methylbutane and 2,2,3-trimethylpentane, both molecules have all their sp3 carbon with different environments. Figure 9 shows the correspondence obtained between the chemical shift in ppm from each proton and the inverse of the LOVI value for each carbon corresponding to that proton (or group of protons). It can be observed that there is an unequivocal numerical proximity between both groups of values.
If the chemical shift is a numerical expression for electronic density of a nuclei as well as the surrounding electronic environment, influenced by the electronic richness of the adjacent atoms, then a linear relationship between the chemical shift and the GDIs may be found, which means that our MDs codify steric and electronic information of the atoms in the molecules. Figure 9B and Figure 10B in each previous case are an approximation of the real spectrum (without taking into account the signal’s multiplicities), where the similarity of both spectra may be observed: the spectrum obtained from the chemical shift and the one obtained from the calculated GDIs. However, when a more realistic spectrum is needed, additional considerations are required: Calculation of the integrated intensity and multiplicity.
The Integrated Intensity (quantities of hydrogen bonded to the atom (i), which provokes the signal) is computed according to the Equation (24):
NH = 4 – δi
On the other hand, the signal’s Multiplicity is determined by Equation (25):
Mi = Σ(4 − δj) + 1
where δ is the vertex degree of an atom and the atoms in the sum are those that possess derivative values different from zero with the atom i, for the order 1, according to the Connected Sub-graph (CS) or Multiplicity (M) event.

Description of Global Electronic Properties. Energy of Resonance

The resonance energy or stabilization energy by resonance is the difference between the energy corresponding to the structure with double (or triple) rigid bonds established in positions located in molecules with alternate unsaturated bonds (the most probable Kekule’s structure) and the real energy of the substance. The latter (i.e., real energy) is less than the former and this decrement is associated with electronic delocalization [43,44]. In this study the correlation between the resonance energy from a dataset composed of 17 aromatic molecules and values calculated by GDIs was determined using several atomic labels (chemical, physical and graph-theoretical atomic properties) and 10 events. The existence of a linear correlation between the previous mentioned property and GDIs implies that they are able to characterize electronic densities and their delocalization capacity (which would corroborate with the study in Epigraph 3.1.4).
The best one and two variable models together with their corresponding statistical parameters are shown in Equations (26) and (27), respectively:
E R = 4.16 ( ± 2.77 ) + 0.48 ( ± 0.02 ) [ N 1 L / I n ] M D f
R2 = 98.25, s = 5.286, Q2Loo = 97.41, sCV = 6.049, Q2Boot = 97.60, ysc = −0.038, F = 844.03.
E R = 4.24 ( ± 2.22 ) + 8.69 ( ± 1.58 ) [ C C / P d N 1 ( K ) ] H D f 0.94 ( ± 0.02 ) [ A G / P l C 5 ( S ) ( I S ) ] H D f
R2 = 99.24, s = 3.613, Q2Loo = 98.99, sCV = 3.764, Q2Boot = 99.01, ysc = 0.017, F = 910.02.
The satisfactory statistical parameters from these models demonstrate the close linear relationship between GDIs and the resonance energy for this group of molecules. The Equation (26) with only one variable is able to explain more than the 98% of the variance of the property. Table 12 shows the experimental and calculated values with the Equations (26) and (27), and their corresponding residuals. It is interesting to point out that the invariant that best correlates with this electronic property is the norm 1, which is the linear combination of the individual LOVI values.
In the previous experiment it was demonstrated that the atomic index values have a close relationship with the electronic properties of atoms in their molecular environment; therefore, as expected one total descriptor such as the norm 1 (N1, see Table 1 for more details) adequately codifies all the contributions and it expresses the electronic characteristics of the molecule, showing also a good correlation with the resonance energy, which is a property that express the electronic behavior product of conjugation. The two variable model contains a local descriptor for atoms forming part of unsaturated bonds. This variable quantifies the electronic effects of all the atoms with these characteristics in the molecule. Both descriptors from Equation (27) were derived with respect to the Sach’s sub-graphs (H) event, because its fragments only take into account the sub-graphs of order 1 and ring sub-graphs [3,16,45,46,47].
The present study adds value to the previous experiments, which show the capacity of quantifying the electronic environment in atoms and molecules with this mathematical approach. The Figure 11 shows the regression graphs and the comparative behavior between experimental and theoretical values from the Equation (27).

4. Conclusions

The capacity of the GDIs to extract relevant structural information of organic molecules and consequently express it in terms of atomic local and total values has been demonstrated. It was observed that GDIs codify information on molecular symmetry, allow for the characterization of molecular structures with different sizes, are sensitive to structural ramifications, adequately characterize differences in the electronic densities of atoms at different positions in the structures, including cases of conjugated systems. Additionally, these TIs take into account the presence of heteroatoms and how they affect the electronic environment of the molecule.
Taking in consideration the regularity and coherence found among GDIs and each one of the structures described by this method, it may be affirmed that the GDIs possess direct structural interpretation allowing for greater comprehension of the chemical information codified [28]. Additionally, it was demonstrated that there exists a relationship between GDIs and the geometric reactivity, seen as a combination of the accessibility to the molecular structure and the activation entropy in the interaction process.
The transformation of the frequency and derivative parameters in electronic terms revealed that GDIs locally codify the characteristics of the electronic distribution of atoms and bonds and can be expressed as the electronic reactivity of atoms and molecules. The application of several mathematic operators to obtain global and local indices over a group of atoms, as well as a combination of these in linear models are an expression of more complex molecular properties such as the energy, analogous to the function of the operators employed in quantum mechanics.

Author Contributions

Oscar Martínez-Santiago, Yovani Marrero-Ponce and Stephen J. Barigye proposed the theory of the GDI indices and designed the theoretical experiment; Oscar Martínez-Santiago and Yovani Marrero-Ponce supervised the QSPR modeling on the all chemical datasets; Luis M. Artiles Martínez, Yovani Marrero-Ponce and Oscar Martínez-Santiago proposed the analogy between the classical and discrete derivative; Liliana Vázquez Infante, Oscar Martínez-Santiago, Cesar H. Zambrano, Huong Le Thi Thu, Jorge L. Muñiz Olite and Maykel Cruz-Monteagudo worked in the QSPR modeling on all chemical datasets; Oscar Martínez-Santiago, Yovani Marrero-Ponce, F. Javier Torres and Ricardo Vivas-Reyes analyzed the results and proposed the interpretation of the GDIs; Oscar Martínez-Santiago wrote the paper; Finally, F. Javier Torres, Ricardo Vivas-Reyes and Maykel Cruz-Monteagudo contributed to the revision of the manuscript too. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Putz, M.V.; Lacrămă, A.-M. Introducing spectral structure activity relationship (S-SAR) analysis. Application to ecotoxicology. Int. J. Mol. Sci. 2007, 8, 363–391. [Google Scholar] [CrossRef]
  2. Estrada, E.; Molina, E. Novel local (fragment-based) topological molecular descriptors for QSPR/QSAR and molecular design. J. Mol. Graph. Model 2001, 20, 54–64. [Google Scholar] [CrossRef]
  3. Barigye, S.J.; Marrero-Ponce, Y.; Martínez López, Y.; Martínez-Santiago, O.; Torrens, F.; Domenech, R.G.; Galvez, J. Event-based criteria in GT-STAF information indices: Theory, exploratory diversity analysis and QSPR applications. SAR QSAR Environ. Res. 2013, 24, 3–34. [Google Scholar] [CrossRef] [PubMed]
  4. Todeschini, R.; Consonni, V. Molecular Descriptors for Chemoinformatics; Wiley-VCH: Hoboken, NJ, USA, 2009; Volume 1–2. [Google Scholar]
  5. Devillers, J.; Balaban, A.T. Topological Indices and Related Descriptors in QSAR and QSPR; Gordon and Breach: Amsterdam, The Netherlands, 1999. [Google Scholar]
  6. Gutman, T.; Ruscic, B.; Trinajstic, N.; Wilcox, C.F. Graph Theory and Molecular Orbitals. XII. Acyclic Polyenes. J. Chem. Phys. 1975, 62, 3399. [Google Scholar] [CrossRef]
  7. Organisation for Economic Co-operation and Development. Guidance Document on the Validation of (Quantitative) Structure-Activity Relationship [(Q)SAR] Models; OECD: Paris, France, 2007; p. 154. [Google Scholar]
  8. Barigye, S.J.; Marrero-Ponce, Y.; Zupan, J.; Pérez-Giménez, F.; Freitas, M.P. Structural and physicochemical interpretation of GT-STAF information theory-based indices. Bull. Chem. Soc. Jpn. 2015, 88, 97–109. [Google Scholar] [CrossRef]
  9. Hoffmann, R. Qualitative thinking in the age of modern computational chemistry—Or what Lionel Salem knows. J. Mol. Struct. THEOCHEM 1998, 424, 3. [Google Scholar] [CrossRef]
  10. Marrero-Ponce, Y.; Martínez-Santiago, O.; Martínez López, Y.; Barigye, S.J.; Torrens, F. Derivatives in discrete mathematics: A novel graph-theoretical invariant for generating new 2/3D molecular descriptors. I. Theory and QSPR application. J. Comput. Aided Mol. Des. 2012, 26, 1229–1246. [Google Scholar] [CrossRef] [PubMed]
  11. Martínez-Santiago, O.; Millán Cabrera, R.; Marrero-Ponce, Y.; Barigye, S.J.; Martínez-López, Y.; Torrens, F.; Pérez-Giménez, F. Discrete derivatives for atom-pairs as a novel graph-theoretical invariant for generating new molecular descriptors: Orthogonality, interpretation and QSARs/QSPRs on benchmark databases. J. Mol. Inform. 2014, 33, 343–368. [Google Scholar] [CrossRef]
  12. Martínez Santiago, O.; Marrero-Ponce, Y.; Millán Cabrera, R.; Barigye, S.; Martínez-López, Y.; Artiles Martínez, L.M.; Guerra de León, J.O.; Perez-Giménez, F.; Torrens, F. Extending graph (discrete) derivative descriptors to n-tuple atom-relations. MATCH Commun. Math. Comput. Chem. 2015, 73, 397–420. [Google Scholar]
  13. Gorbátov, V.A. Fundamentos de la Matemática Discreta; Editorial Mir: Moscú, Russia, 1988. [Google Scholar]
  14. Avery, J. The Quantum Theory of Atoms, Molecules and Photons; McGraw-Hill Book Company (UK) Ltd.: London, UK, 1972. [Google Scholar]
  15. Cockett, M.; Doggett, G. Maths for Chemists; Royal Society of Chemistry: Cambridge, UK, 2003; Volume 1. [Google Scholar]
  16. Martínez Santiago, O.; Millán Cabrera, R.; Marrero-Ponce, Y.; Barigye, S.J.; Martínez López, Y.; Torrens, F. Topo-Chemical Extended Indices Derived from Event-Based Graph Derivative. Theory, Analysis and QSPR Validation. Curr. Pharm. Des. in revision.
  17. Barigye, S.J.; Marrero-Ponce, Y.; Martínez-López, Y.; Artiles Martínez, L.M.; Pino-Urias, R.W.; Martínez-Santiago, O.; Torrens, F. Relations frequency hypermatrices in mutual, conditional and joint entropy-based information indices. J. Comput. Chem. 2013, 34, 259–274. [Google Scholar] [CrossRef] [PubMed]
  18. Balaban, A.T. Numerical Modelling of Chemical Structures: Local Graph Invariants and Topological Indices. In Graph Theory and Topology in Chemistry; King, R.B., Rouvray, D.H., Eds.; Elsevier: Amsterdam, The Netherlands, 1987; pp. 159–176. [Google Scholar]
  19. Ivanciuc, O.; Balaban, T.S.; Balaban, A.T. Design of topological indices. Part 4. Reciprocal distance matrix, related local vertex invariants and topological indices. J. Math. Chem. 1993, 12, 309–318. [Google Scholar] [CrossRef]
  20. Magnasco, V. Models for Bonding in Chemistry; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2010. [Google Scholar]
  21. Levine, I.N. Química Cuántica; Pearson Educación: Madrid, Spain, 2001. [Google Scholar]
  22. Atkins, P.; Friedman, R. Molecular Quantum Mechanics; Oxford University Press: New York, NY, USA, 2005. [Google Scholar]
  23. García-Jacas, C.R.; Marrero-Ponce, Y.; Barigye, S.J.; Valdés-Martiní, J.R.; Rivera-Borroto, O.M.; Verbel, O. N-linear algebraic maps to codify chemical structures: Is a suitable generalization to the atom-pairs approaches? Curr. Drug Metab. 2014, 15, 441–469. [Google Scholar] [CrossRef] [PubMed]
  24. García-Jacas, C.R.; Marrero-Ponce, Y.; Acevedo-Martínez, L.; Barigye, S.J.; Valdés-Martiní, J.R.; Contreras-Torres, E. QuBiLS-MIDAS: A parallel free-software for molecular descriptors computation based on multi-linear algebraic maps. J. Comput. Chem. 2014, 35, 1395–1409. [Google Scholar] [CrossRef] [PubMed]
  25. Marrero-Ponce, Y.; Martínez Santiago, O.; Barigye, S.J. Divati 1.0; Unit of Computer-Aided Molecular “Biosilico” Discovery and Bioinformatic Research (CAMD-BIR Unit): Santa Clara, Cuba, 2013. [Google Scholar]
  26. Todeschini, R.; Ballabio, D.; Consonni, V.; Mauri, A.; Pavan, M. Mobydigs. Version 1.0; Talete SRL Ed.: Milano, Italy, 2004. [Google Scholar]
  27. Leardi, R. Nature-Inspired Methods in Chemometrics: Genetic Algorithms and Artificial Neural Networks; Department of Pharmaceutical and Food Chemistry and Technology, University of Genova: Genova, Italy, 2003. [Google Scholar]
  28. Randic, M. Generalized molecular descriptors. J. Math. Chem. 1991, 7, 155. [Google Scholar] [CrossRef]
  29. Guerasimov, Y.; Dreving, V.; Eriomin, E.; Kiseliov, A.; Lebedev, V.; Panchenkov, G.; Shliguin, A. Curso de Química Física Tomo II; Editorial Mir: Moscú, Russia, 1971. [Google Scholar]
  30. Frost, A.A.; Pearson, R.G. Kinetics and Mechanism, 2nd ed.; Wiley: New York, NY, USA, 1961. [Google Scholar]
  31. Kier, L.B.; Hall, L.H. Molecular connectivity: Intermolecular accessibility and encounter simulation. J. Mol. Graph. Model. 2011, 20, 76–83. [Google Scholar] [CrossRef]
  32. Kier, L.B.; Hall, L.H. Intermolecular accessibility: The meaning of molecular connectivity. J. Chem. Inf. Comput. Sci. 2000, 40, 792–795. [Google Scholar] [CrossRef] [PubMed]
  33. Estrada, E. The Structural Interpretation of the Randic Index. Internet Electron. J. Mol. Des. 2002, 1, 360–366. [Google Scholar]
  34. Estrada, E. Physicochemical interpretation of molecular connectivity indices. J. Phys. Chem. A 2002, 106, 9085–9091. [Google Scholar] [CrossRef]
  35. Tudoran, M.A.; Putz, M.V. Molecular graph theory: From adjacency information to colored topology by chemical reactivity. Curr. Org. Chem. 2015, 19, 359–386. [Google Scholar] [CrossRef]
  36. Putz, M.V.; Tudoran, M.A.; Ori, O. Topological organic chemistry: From distance matrix to Timisoara eccentricity. Curr. Org. Chem. 2015, 19, 249–273. [Google Scholar] [CrossRef]
  37. Estrada, E.; Molina, E. 3D conectivity indices in QSPR/QSAR studies. J. Chem. Inf. Comput. Sci. 2001, 41, 791. [Google Scholar] [CrossRef] [PubMed]
  38. Ortiz del Toro, P.J.; Pérez Martínez, C.S. Química Cuántica. Elementos de Estructura Molecular; MES (acronym of Ministerio de Eduación Superior): Ciudad de la Habana, Cuba, 1984. [Google Scholar]
  39. Kier, L.B.; Hall, L.H. The e-state as an extended free valence. J. Chem. Inf. Comput. Sci. 1997, 37, 548–552. [Google Scholar] [CrossRef]
  40. Kier, L.B.; Hall, L.H. Molecular Structure Description. The Electrotopological State; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  41. Pérez Martínez, C.S.; Ortiz del Toro, P.J.; Alonso Becerra, E. Resonancia Magnética Nuclear; MES (acronym of Ministerio de Eduación Superior): Ciudad de la Habana, Cuba, 1983. [Google Scholar]
  42. Cambridgesoft-Corporation. Chemdraw; Cambridgesoft-Corporation: Cambridge, MA, USA, 2003. [Google Scholar]
  43. Morrison, R.T.; Boyd, R.N. Organic Chemistry, 7th ed.; Prentice-Hall Inc.: New Jersey, USA, 1992. [Google Scholar]
  44. Putz, M.V. Compactness aromaticity of atoms in molecules. Int. J. Mol. Sci. 2010, 11, 1269–1310. [Google Scholar] [CrossRef] [PubMed]
  45. Sachs, H. Beziehungen Zwischen den in Einem Graphen Enthaltenen Kreisen und Seinem Characteristischen Polynom; Debrecen, Hungary, 1964; Volume 11. [Google Scholar]
  46. Skvortsova, M.I.; Stankevich, I.V. Eigenvectors of weighted graphs: A supplement to Sachs’ theorem. J. Mol. Struct. 2005, 719, 213–223. [Google Scholar] [CrossRef]
  47. Gutman, I. Impact of the Sachs theorem on theoretical chemistry: A participant’s testimony. MATCH Commun. Math. Comput. Chem. 2003, 48, 17–34. [Google Scholar]
Figure 1. Example of graph G (A), fragmentation according to an event “S” (B) and the sub-graphs sets for vertexes a and b (C).
Figure 1. Example of graph G (A), fragmentation according to an event “S” (B) and the sub-graphs sets for vertexes a and b (C).
Ijms 17 00812 g001
Figure 2. Scheme of the analogy between derivatives and their mathematical development. (A) Obtaining of the discrete derivative over a pair of elements i j from a graph G; (B) Algebraic development of the process for obtaining the discrete derivative over pairs of vertexes; (C) Obtaining of the classical derivative from the mathematical analysis.
Figure 2. Scheme of the analogy between derivatives and their mathematical development. (A) Obtaining of the discrete derivative over a pair of elements i j from a graph G; (B) Algebraic development of the process for obtaining the discrete derivative over pairs of vertexes; (C) Obtaining of the classical derivative from the mathematical analysis.
Ijms 17 00812 g002
Figure 3. (A) Molecular structure of 2-amino-5-vinylfurane; (B) Corresponding graph with arbitrary numeration.
Figure 3. (A) Molecular structure of 2-amino-5-vinylfurane; (B) Corresponding graph with arbitrary numeration.
Ijms 17 00812 g003
Figure 4. Graphical behavior of the steric reactivity based in LOVIs calculated by GDIs and the Relative Bond Accessibility Area (RBA) proposed by Estrada for evaluating the accessibility.
Figure 4. Graphical behavior of the steric reactivity based in LOVIs calculated by GDIs and the Relative Bond Accessibility Area (RBA) proposed by Estrada for evaluating the accessibility.
Ijms 17 00812 g004
Figure 5. Regularity in the Derivative (over exocyclic double bond) variation and the logK values.
Figure 5. Regularity in the Derivative (over exocyclic double bond) variation and the logK values.
Ijms 17 00812 g005
Figure 6. Own Frequency and duplex derivative in electronic terms.
Figure 6. Own Frequency and duplex derivative in electronic terms.
Ijms 17 00812 g006
Figure 7. Behavior of the free valence, the E-States and GDIs for atoms from 19 conjugated molecules.
Figure 7. Behavior of the free valence, the E-States and GDIs for atoms from 19 conjugated molecules.
Ijms 17 00812 g007
Figure 8. Development of the one-variable linear regression models obtained for each event: (A) Ethers; (B) Aldehydes and Ketones.
Figure 8. Development of the one-variable linear regression models obtained for each event: (A) Ethers; (B) Aldehydes and Ketones.
Ijms 17 00812 g008aIjms 17 00812 g008b
Figure 9. Methylbutane. (A) Equivalence between LOVIs values and the chemical shift in ppm; (B) Quantity of protons that provoke the signal vs. LOVIs and ppm.
Figure 9. Methylbutane. (A) Equivalence between LOVIs values and the chemical shift in ppm; (B) Quantity of protons that provoke the signal vs. LOVIs and ppm.
Ijms 17 00812 g009
Figure 10. 2,2,3-Trimethylpentane. (A) Equivalence between LOVIs values and the chemical shift in ppm; (B) Number of protons that provoke the signal vs. LOVIs and ppm.
Figure 10. 2,2,3-Trimethylpentane. (A) Equivalence between LOVIs values and the chemical shift in ppm; (B) Number of protons that provoke the signal vs. LOVIs and ppm.
Ijms 17 00812 g010
Figure 11. Regression and prediction graphs for the Equation (27).
Figure 11. Regression and prediction graphs for the Equation (27).
Ijms 17 00812 g011
Table 1. Invariants functions employed to derive molecular descriptors (total and local) from local vertex invariants (LOVIs).
Table 1. Invariants functions employed to derive molecular descriptors (total and local) from local vertex invariants (LOVIs).
No.GroupNameIDFormula a
1Norms (Metrics)Minkowski norms ( p = 1) Manhattan normN1 N 1 = i = 1 n | L i |
2Minkowski norms ( p = 2) Euclidean normN2 N 2 = i = 1 n | L i | 2
3Minkowski norms ( p = 3)N3 N 3 = i = 1 n | L i | 3 3
4Chebyshev distanceNI N I = lim n [ i = 1 n L i n ] 1 n
5Mean (first statistical moment)Geometric MeanG G = i = 1 n L i n
6Arithmetic Mean (Power mean of degree α = 1)M (or M1) M α = ( L 1 α + L 2 α + + L n α n ) 1 α
7Quadratic Mean (Power mean of degree α = 2)P2 (or M2)
8Power mean of degree α = 3P3 (or M3)
9Harmonic Mean (Power mean of degree α = −1)A (or M−1)
10Statistical (highest statistical moments)VarianceV V = i = 1 n ( L i M ) 2 n 1 ,   M: arithmetic mean
11SkewnessS S = n ( X 3 ) / [ ( n 1 ) ( n 2 ) ( D E ) 3 ] , n, number of vertices. X 3 = i = 1 n ( L i M ) 3 , M: arithmetic mean, DE, standard deviation
12KurtosisK K = [ n ( n + 1 ) ( X 4 ) 3 ( X 2 ) ( X 2 ) ( n 1 ) ] [ ( n 1 ) ( n 2 ) ( n 3 ) ( DE ) 4 , n, number of vertices; X j = i = 1 n ( L i M ) j , M: arithmetic mean, DE, standard deviation
13Standard DeviationDE D E = i = 1 n ( L i M ) 2 n 1
14Variation CoefficientCV C V = D E M
15RangeR R = L m a x L m i n
16Percentile 25Q1 Q 1 = ( N 4 + 1 2 ) , N: Li number
17Percentile 50Q2 Q 2 = ( N 2 + 1 2 ) , N: Li number
18Percentile 75Q3 Q 3 = ( 3 N 4 + 1 2 ) , N: Li number
19Inter-quartile RangeI50 I 50 = Q 3 Q 1
20Maximum valueMXMX = Li max
21Minimum valueMNMN = Li min
22Classical (Invariants)AutocorrelationACk A C k = i = 1 n j 1 n L i × L j · [ δ ( d i j , k ) ] , k = 1 , 2 , 7
23GravitationalGIk G I k = 1 n i = 1 n j = 1 n L i L j d i j k · [ δ ( d i j , k ) ] , k = 1 , 2 , 7
24Total sum at lag kTSk T S k = i = 1 n j = 1 n L i j · [ δ ( d i j , k ) ] , k = 1 , 2 , 7
25Kier-Hall connectivityCNm K H t m = i = 1 n ( i = 1 n k L i , w ) k λ
where, k is the number of sub-graphs, nk is the number of atoms in a fragment, λ is equal to ½, m and t are the sub-graph order and type, respectively
26Mean Information ContentMI M I = i = 1 A N g N 0 l o g 2 N g N 0
where, Ng is the number of atoms with the same LOVI value; No, is the number of atoms in a molecule
27Total Information ContentTI T I = N 0 l o g 2 N 0 g = 1 G N g l o g 2 N g
28Standardized Information ContentSI S I = T I N 0 l o g 2 N 0
29Electrotopological state (E-state index)ES S I = I i + Δ I i = I i + j = 1 n I i I j ( d i j + 1 ) 2  
where, Ii is the intrinsic state of the ith atom and ΔIi is the field effect on the ith atom calculated as perturbation of the Ii of ith atom by all other atoms in the molecule, dij is the topological distance between the ith and the jth atoms, and n is the number of atoms. The exponent k is 2.
30Ivanciuc-Balaban Type-IndicesIB J k = n 2 Δ B n + C + 1 i = 1 n 1 j = i + 1 n α i j ( L i × L j ) 1 2
where, the summation goes over all pairs of atoms but only pairs of adjacent atoms are accounted for by means of the elements αij of the adjacency matrix. The n, B, and C are the number of atoms, bonds, and rings (cyclomatic number), respectively.
a The formulae used in these invariants, are simplified forms of general equations given that the vector   y ¯ is constituted of the coordinates of the origin. For example, in the case of the Euclidean norm (N2), the general formula is: || x ¯ || 2 = i = 1 n ( x i y i ) 2 + ( x j y j ) 2 + ( x z y z ) 2 . However given that   y ¯ = (0, 0, 0), this formula reduces to || x ¯ || 2 = i = 1 n | x i | 2 .
Table 2. Codification of the chain size, multiple bonds and their positions in the molecule.
Table 2. Codification of the chain size, multiple bonds and their positions in the molecule.
MoleculeLOVIs ValuesTotal Invariants
C1C2C3C4C5C6N1ARA
Ijms 17 00812 i00122420
Ijms 17 00812 i002646165.332
Ijms 17 00812 i00311.176.176.1711.1734.678.675
Ijms 17 00812 i00417.339.337.339.3317.3360.6712.3310
Ijms 17 00812 i00524.413.329.589.5813.3224.494.615.7714.82
Ijms 17 00812 i00626.6316.419.9210.0313.2329.57105.7917.6319.65
Ijms 17 00812 i00735.4321.0611.4612.0815.5737.15132.7522.1325.69
Ijms 17 00812 i00827.913.910.2512.3117.2126.27107.8317.9717.65
Ijms 17 00812 i00928.5522.3816.1711.4215.2332.15125.920.9820.73
Ijms 17 00812 i01026.7313.6112.2112.2113.6126.73105.117.5214.53
Ijms 17 00812 i01129.6514.48161614.4829.65120.2720.0415.17
Ijms 17 00812 i01226.0515.6310.8310.3813.5331.9108.3218.0521.53
Ijms 17 00812 i01322.8812.389912.3822.8888.5314.7613.88
Table 3. Differentiation among chain isomers.
Table 3. Differentiation among chain isomers.
MoleculeLOVIs ValuesTotal Invariants
C1C2C3C4C5C6N1ARA
Ijms 17 00812 i01424.413.329.589.5813.3224.494.615.7714.82
Ijms 17 00812 i01516.7911.018.2611.6220.2216.7984.7014.1211.96
Ijms 17 00812 i01618.049.799.59.7918.041580.1413.368.54
Ijms 17 00812 i01714.3310.3510.3514.3514.3314.3378133.975
Ijms 17 00812 i01813.62513.259.34716.2213.6313.6379.6913.286.88
Table 4. Evaluation of some organic molecules with heteroatoms.
Table 4. Evaluation of some organic molecules with heteroatoms.
MoleculeType of Atom in XLOVIs ValuesTotal Invariants
C1C2C3XN1ARA
Ijms 17 00812 i019CH311.176.176.1711.1734.678.675.00
NH211.376.876.6412.5437.419.355.90
OH11.697.447.0213.8139.959.996.79
F11.997.907.3214.8842.0910.527.56
Table 5. Codification of cyclization, conjugation and aromaticity.
Table 5. Codification of cyclization, conjugation and aromaticity.
MoleculeLOVIs ValuesTotal Invariants
C1C2N1ARA
Ijms 17 00812 i02016.298.7657.7211.548.68
Ijms 17 00812 i0218.978.9754.949.160.58
Ijms 17 00812 i02210.335.3331.337.835
Ijms 17 00812 i0239.108.9755.59.250.71
Ijms 17 00812 i0248.398.3950.348.390
Ijms 17 00812 i02517.2014.93146.8614.698.02
Ijms 17 00812 i0267.5812.3340.510.138.75
Ijms 17 00812 i0275.675.6722.675.670
Ijms 17 00812 i0289.689.6855.509.250.71
Ijms 17 00812 i0298.468.4951.268.540.39
Ijms 17 00812 i0307.737.9941.608.322.44
Table 6. Symmetry Index for some organic molecules.
Table 6. Symmetry Index for some organic molecules.
MoleculesSIMoleculesSI
Ijms 17 00812 i0311.00 Ijms 17 00812 i0321.00
Ijms 17 00812 i0330.75 Ijms 17 00812 i0340.6554
Ijms 17 00812 i0350.60 Ijms 17 00812 i0360.5870
Ijms 17 00812 i0370.5793 Ijms 17 00812 i0380.50
Ijms 17 00812 i0390.4911 Ijms 17 00812 i0400.00
Table 7. Activation entropy and molecular interaction capacity for the dimerization of second order for unsaturated molecules in gaseous state.
Table 7. Activation entropy and molecular interaction capacity for the dimerization of second order for unsaturated molecules in gaseous state.
MoleculeΩ (GDI)ΔS#
Isobutylene1.19988−6
2-Isoprene0.40035−8
Ethylene0−12
1,3-Butadiene0−13
1,3-Pentadiene−0.1948−14
Propylene−0.5501−15
Cyclopentadiene−2.1654−26
Table 8. Accessibility to the chemical bond (Connectivity Indices and Graph Derivative Indices (GDIs)).
Table 8. Accessibility to the chemical bond (Connectivity Indices and Graph Derivative Indices (GDIs)).
NameFormulaR-GDIRBA (R2)
1: ethane Ijms 17 00812 i041Infinite1
2: propane Ijms 17 00812 i0420.50.7071
3: methyl propane Ijms 17 00812 i0430.08330.5774
4: 2,2-dimetylpropane Ijms 17 00812 i0440.02780.5
5: n-butane Ijms 17 00812 i0450.11110.5
6: methyl butane Ijms 17 00812 i0460.03570.4082
7: 2,2-dimethyl butane Ijms 17 00812 i0470.01540.3536
8: 2,3-dimethyl butane Ijms 17 00812 i0480.01560.3333
9: 2,2,3-trimethyl butane Ijms 17 00812 i0490.00790.2887
10: 2,2,3,3-tetramethyl butane Ijms 17 00812 i0500.00440.25
Table 9. C t values for each type of carbon atoms.
Table 9. C t values for each type of carbon atoms.
Type of Atom H 2 C = CH = CH > C = Ijms 17 00812 i051 Ijms 17 00812 i052
C t 0.1364−0.0371−0.0778−0.1211−0.1401.380
Table 10. Chemical Shift ( δ ) in 17O-NMR and GDI values.
Table 10. Chemical Shift ( δ ) in 17O-NMR and GDI values.
No.CompoundsE-State a [ A N / I n ( I S ) ] B D f bδ (ppm) cδp (ppm) d
1CH3CHO8.8060.167592.0591.12
2C2H5CHO9.1741.319579.5582.29
3i-C3H7CHO9.5052.069574.557 5.24
4(CH3)2CO9.4443.250569.0564.41
5CH3COC2H59.8134.000557.0558.66
6i-C3H7COCH310.1444.417557.0554.57
7(C2H5)2CO10.1815.000547.0550.41
8i-C3H7COC2H510.5125.833543.5542.41
9(i-C3H7)CO10.8436.667535.0535.77
a E-State Index; b Arithmetic mean of LOVIs from oxygen and carbon atoms in double bond; c Chemical shifts of the 17O-MNR; d Chemical shifts of the 17O calculated by Equation (22).
Table 11. Chemical Shift ( δ ) in 17O-NMR and GDI values.
Table 11. Chemical Shift ( δ ) in 17O-NMR and GDI values.
No.CompoundE-State a
[ A T / I n ( H T ) ] M D / f
b
δ (ppm) cδp (ppm) d
1Methoxymethane4.2042.58−52.2−53.12
2Methoxyethane4.5435.74−22.5−22.64
32-Methoxypropane4.7528.80−2−1.56
42-Methoxy-2-methylpropane4.9425.158.59.36
5Ethoxyethane4.8328.906.57.72
62-Ethoxypropane5.0421.962828.75
72-Ethoxy-2-methylpropane5.2318.3140.539.50
82-Isopropoxypropane5.2515.0252.550.84
92-Isopropoxy-2-methylpropane5.4411.3762.562.59
102-t-Butoxy-2-methylpropane5.637.727676.37
a E-State Index; b LOVI of oxygen atoms; c Chemical shifts of the 17O-MNR; d Chemical shifts of the 17O calculated by Equation (23).
Table 12. Experimental resonance energy calculated by the GDI-based models.
Table 12. Experimental resonance energy calculated by the GDI-based models.
MoleculesERexpERcal (Ec. 24)Res. aERcal (Ec. 25)Res. a
Benzene3628.61−7.3935.170.83
Naphthalene6157.79−3.2157.53.5
Anthracene8386.973.9780.772.23
Phenanthrene9186.25−4.7588.072.93
Styrene3839.471.4729.548.46
Stilbene7474.810.8170.953.05
Biphenyl7165.94−5.0669.541.46
Butadiene3.512.298.794.57−1.07
Fluorene7678.492.4983.19−7.19
3,5-Triphenylbenzene149140.6−8.4147.611.39
Toluene3533.14−1.8638.01−3.01
O-Xylene3536.951.9541.5−6.5
Diphenylmethane6765.97−1.0371.6−4.6
Naphthacene110116.156.15111.58−1.58
Chrysene116.5115.26−1.24115.321.18
Pyrene108.9107.11−1.79108.05−0.05
Perylene126.3135.399.09126.75−0.45
a Residual (ERexp–Ercal).

Share and Cite

MDPI and ACS Style

Martínez-Santiago, O.; Marrero-Ponce, Y.; Barigye, S.J.; Le Thi Thu, H.; Torres, F.J.; Zambrano, C.H.; Muñiz Olite, J.L.; Cruz-Monteagudo, M.; Vivas-Reyes, R.; Vázquez Infante, L.; et al. Physico-Chemical and Structural Interpretation of Discrete Derivative Indices on N-Tuples Atoms. Int. J. Mol. Sci. 2016, 17, 812. https://doi.org/10.3390/ijms17060812

AMA Style

Martínez-Santiago O, Marrero-Ponce Y, Barigye SJ, Le Thi Thu H, Torres FJ, Zambrano CH, Muñiz Olite JL, Cruz-Monteagudo M, Vivas-Reyes R, Vázquez Infante L, et al. Physico-Chemical and Structural Interpretation of Discrete Derivative Indices on N-Tuples Atoms. International Journal of Molecular Sciences. 2016; 17(6):812. https://doi.org/10.3390/ijms17060812

Chicago/Turabian Style

Martínez-Santiago, Oscar, Yovani Marrero-Ponce, Stephen J. Barigye, Huong Le Thi Thu, F. Javier Torres, Cesar H. Zambrano, Jorge L. Muñiz Olite, Maykel Cruz-Monteagudo, Ricardo Vivas-Reyes, Liliana Vázquez Infante, and et al. 2016. "Physico-Chemical and Structural Interpretation of Discrete Derivative Indices on N-Tuples Atoms" International Journal of Molecular Sciences 17, no. 6: 812. https://doi.org/10.3390/ijms17060812

APA Style

Martínez-Santiago, O., Marrero-Ponce, Y., Barigye, S. J., Le Thi Thu, H., Torres, F. J., Zambrano, C. H., Muñiz Olite, J. L., Cruz-Monteagudo, M., Vivas-Reyes, R., Vázquez Infante, L., & Artiles Martínez, L. M. (2016). Physico-Chemical and Structural Interpretation of Discrete Derivative Indices on N-Tuples Atoms. International Journal of Molecular Sciences, 17(6), 812. https://doi.org/10.3390/ijms17060812

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop