Computational methodology
The ChemAxon Calculator Plugins for Conformers Searching available in Marvin View 17.15.0 were used to generate 3D structures from SMILES strings, and to propose low energy conformers for structure property prediction and calculation. For the molecule considered in the current study, the lowest energy conformers were used as a starting point for the geometry optimization. The geometries of all the selected conformers were optimized with the DFTBA program. The molecular structures of the five lowest energy conformers were reoptimized by resorting to the MN12SX/Def2TZVP/H2O model chemistry. The optimized structures were confirmed to be real minima by vibrational frequency analysis (no imaginary frequency).
Results and discussion
The molecular structures of the optimized conformers of Taltobulin obtained as mentioned in the previous section, and whose graphical sketch is shown in Fig. 1, has been submitted to optimization in absence of solvent by resorting the DFTBA model available in Gaussian 09 [15] and then reoptimized using the MN12SX/Def2TZVP/H2O model chemistry mentioned in that section. Having verified that each of the structures corresponded to the minimum energy conformations by running a frequency calculation analysis, the electronic properties were determined by using the same model chemistry.
Becke has recently mentioned that the adiabatic connection and the ideas of Hohenberg, Kohn, and Sham applying only to electronic ground states is a common misconception [16]. In this regard, the HOMO–LUMO gap within the KS model represents a nice estimation of the lowest excitation energy [17]. Under this assumption, the determination of the maximum wavelength absorption of the Taltobulin peptide was performed by conducting a ground-state calculation with the aforementioned density functional at the same level of model chemistry and theory to obtain the HOMO–LUMO gap and subsequently, the \(\lambda _{max}\). The electronic energy of the neutral molecular system of Taltobulin, the HOMO and LUMO orbital energies, and the maximum absorption wavelengths \(\lambda _{max}\) calculated with the MN12SX/Def2TZVP/H2O model chemistry are − 1517.422 au, − 6.240 eV, − 1.733 eV, and 275 nm, respectively.
According to the results obtained when studying melanoidins [7,8,9,10,11,12,13] as well as peptides from marine sources [14], it can be said that the calculations performed with the MN12SX density functional render HOMO and LUMO energies that satisfy the approximate Koopmans’ theorem. Thus, the application of the KID procedure will be justified. The global reactivity descriptors Electronegativity \(\chi \) [18, 19], Global Hardness \(\eta \) [18, 19], Electrophilicity \(\omega \) [20], Electrodonating Power \(\omega ^{-}\) [21], Electroaccepting Power \(\omega ^{+}\) [21] and Net Electrophilicity \(\Delta \omega ^{\pm }\) [22] were calculated by resorting to the HOMO and LUMO energies determined with the MN12SX density functional with results being \(\chi \) = 3.986 eV, \(\eta \) = 4.507 eV, \(\omega \) = 1.763 eV, \(\omega ^{-}\) = 5.800 eV, \(\omega ^{+}\) = 1.814, and \(\Delta \omega ^{\pm }\) = 7.614 eV. The interested reader in the mathematical formulations of these reactivity descriptors is referred to the original works and to our previous research on the field [7,8,9,10,11,12,13,14]. As expected from the molecular structure of this species, its electrodonating ability is more important that its electroaccepting character.
We now turn our attention to the local descriptors of chemical reactivity, namely the Electrophilic Fukui function \(f^{-}(\mathbf {r})\) [18, 19, 23], the Nucleophilic Fukui function \(f^{+}(\mathbf {r})\) [18, 19, 23] and the Dual Descriptor (DD) \({\Delta }f(\mathbf {r})\) [24,25,26,27,28]. As for the case of the global reactivity descriptors, the interested reader in the mathematical formulations of these reactivity descriptors is referred to the original works and to our previous research on the field [7,8,9,10,11,12,13,14]. The Electrophilic Fukui functions \(f^{-}(\mathbf {r})\) and Nucleophilic Fukui functions \(f^{+}(\mathbf {r})\) for the Taltobulin peptide are shown in Fig. 2.
Within the field of Chemoinformatics applied to the discovery of new pharmaceutical drugs, it is usual to verify the drug-likeness of the involved molecules resorting to some empirical rules, as for example, the Lipinski Rule of Five (Ro5) [29]. This can be achieved by using the readily available Molinspiration software and the results for the case of Taltobulin are presented next as miLogP (the octanol/water partition coefficient) = 4.43, TPSA (the molecular polar surface area) = 98.73, nAtoms (the number of atoms of the molecule) = 34, nON (the number of hydrogen bond acceptors) = 7, nOHNH (the number of hydrogen bond donors) = 3, nviol (the number of violations of the Ro5) = 0, nrotb (the number of rotatable bonds) = 11, volume (the molecular volume) = 479.94 and MW (the molecular weight) = 473.66. Although this criteria cannot always be applied in general to peptides, it can be seen from the previous results that for Taltobulin the number of violations of the Ro5 is 0, which means that Taltobulin can be considered as a druggable molecule.
As the next step, a new task was accomplished by resorting to the same software for the determination of the bioactivity scores for different drug targets whose values for the Taltobulin peptide are GPCR Ligand = 0.43, Ion Channel Modulator = 0.15, Kinase Inhibitor = - 0.12, Nuclear Receptor Ligand = 0.19, Protease Inhibitor = 0.68 and Enzyme Inhibitor = 0.42.
These bioactivity scores for organic molecules can be interpreted as active (when the bioactivity score is > 0), moderately active (when the bioactivity score lies between − 5.0 and 0.0) and inactive (when the bioactivity score \(<-\)5.0). That means that the Taltobulin peptide can be considered a potentially bioactive as a Protease inhibitor, besides being able to act a ligand for GPCR and as an Enzyme inhibitor.