Conceptual DFT as a chemoinformatics tool for the study of the Taltobulin anticancer peptide

Objective A well-behaved model chemistry previously validated for the study of the chemical reactivity of peptides was considered for the calculation of the molecular properties and structure of the Taltobulin anticancer peptide. A methodology based on Conceptual Density Functional Theory (CDFT) was chosen for the determination of the reactivity descriptors. Results The molecular active sites were associated with the active regions of the molecule were associated with the nucleophilic and electrophilic Fukui functions. Finally, the bioactivity scores for the Taltobulin peptide are predicted through a homology methodology relating them with the calculated reactivity descriptors.


Introduction
The biodiversity of the marine environment and the associated chemical diversity constitute a practically unlimited resource of new antitumor agents in the field of the development of marine bioactive substances [1].
Hemiasterlins comprise a small family of naturally occurring N-methylated tripeptides with highly alkylated unnatural amino acids, was originally isolated from the sponge Hemiasterella minor (class Demospongiae, order Hadromedidia, family Hemiasterllidae). Hemiasterlins act as potent tumor growth inhibitors [2,3].
A synthetic analogue of hemiasterlin, taltobulin (HTI-286) has been synthesized in which the 3-substituted indole ring was replaced by phenyl group. Taltobulin inhibits the polymerization of purified tubulin and disrupts microtubule organization in cells. Then, it is considered as a potent inhibitor of proliferation and has substantially less interaction with multidrug resistance protein than currently used antimicrotubule agents [4].
Assuming that an understanding of the chemical interactions is essential for the development of new pharmaceutical drugs, in this work we will be studying the chemical reactivity properties of Talbodulin by resorting to the Conceptual DFT methodology, which will allow the determination of the global properties as well as the local properties for the prediction of the active reaction sites, both electrophilic and nucleophilic. In a similar way, the descriptors of bioactivity (bioactivity scores) will be established through a procedure described in the literature [5,6] trying to relate them with the global and local CDFT reactivity descriptors that result from a calculation protocol based on DFT already validated by our group in previous research [7][8][9][10][11][12][13][14].

Computational methodology
The 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 max . The electronic energy of the neutral molecular system of Taltobulin, the HOMO and LUMO orbital energies, and the maximum absorption wavelengths max calculated with the MN12SX/ Def2TZVP/H2O model chemistry are − 1517.422 au, − 6.240 eV, − 1.733 eV, and 275 nm, respectively.
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)  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  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.

Limitations
In this research note, we have presented the results of a study of the chemical reactivity of a the Taltobulin anticancer peptide on the basis of the Conceptual DFT as a tool to explain the molecular interactions, the molecular properties related to bioavailability and the bioactivity scores.
However, this information is insufficient to be considered for peptidomimetics studies and additional results from other peptides will be necessary. Authors' contributions DGM conceived and designed the research and headed, wrote and revised the manuscript, while JF and NFH contributed to the writing and the revision of the article. All authors read and approved the final manuscript.