Next Article in Journal
Zinc Porphyrin-Functionalized Fullerenes for the Sensitization of Titania as a Visible-Light Active Photocatalyst: River Waters and Wastewaters Remediation
Next Article in Special Issue
In Silico Evaluation of Ibuprofen and Two Benzoylpropionic Acid Derivatives with Potential Anti-Inflammatory Activity
Previous Article in Journal
Metabolic Profiling of Pitaya (Hylocereus polyrhizus) during Fruit Development and Maturation
Previous Article in Special Issue
Novel Approach for the Search for Chemical Scaffolds with Activity at Both Acetylcholinesterase and the α7 Nicotinic Acetylcholine Receptor: A Perspective on Scaffolds with Dual Activity for the Treatment of Neurodegenerative Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chemical Reactivity Theory and Empirical Bioactivity Scores as Computational Peptidology Alternative Tools for the Study of Two Anticancer Peptides of Marine Origin

by
Juan Frau
1,†,
Norma Flores-Holguín
2,† and
Daniel Glossman-Mitnik
1,2,*,†
1
Departament de Química, Universitat de les Illes Balears, 07122 Palma de Mallorca, Spain
2
Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Miguel de Cervantes 120, Complejo Industrial Chihuahua, Chihuahua, Chih 31136, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 4 February 2019 / Revised: 5 March 2019 / Accepted: 19 March 2019 / Published: 21 March 2019
(This article belongs to the Special Issue Computational Approaches for Drug Discovery)

Abstract

:
This work presents an account of the reactivity behavior of the anticancer marine drugs, Soblidotin and Tasidotin, based on the calculation of the global and local descriptors resulting from Chemical Reactivity Theory (CRT), also known as Conceptual DFT, for their consideration as a useful complement to approximations based on Molecular Docking. The information on the global and local reactivity descriptors of the Soblidotin and Tasidotin molecules, obtained through our proposed methodology, may be used for the design of new pharmaceutical analogs by relying on the chemical interactions between these peptides and their protein-type biological receptors. It can be concluded that the CRT approximation to the global and local chemical reactivity, based on the descriptors, can provide interesting information for the consideration of both molecules as potential therapeutic drugs. This is complemented by a study on Advanced Glycation Endproduct (AGE) inhibition, by comparison with the usual molecular systems considered for the task, as a re-purposing study. Finally, the bioactivity scores for Soblidotin and Tasidotin are predicted through an empirical procedure, based on comparison with molecular structures with well-known pharmacological properties.

Graphical Abstract

1. Introduction

There is considerable interest in finding new molecular systems that could be a basis for the design of novel therapeutic drugs. For the accomplishment of this objective, researchers have carried out searches for these potential therapeutic molecules in nature, based on experience acquired in recent decades. The sea is one of those places in which searches are carried out, as it is inhabited by an enormous variety of animal, plant, and microbial life, which generates diverse molecules and remains relatively unexplored.
Among the chemical species that can be found in marine life, peptides are one of the most important chemical groups. Peptides are molecules formed by the combination of natural amino acids and some of their variants, and represent ideal chemical species for the design of therapeutic drugs, since they possess interesting properties inherent to the molecular systems, ranging from small amino acids to large proteins, with which they are often compared.
Soblidotin, also known as Auristatin PE, has been investigated due its potential anti-tumor effects [1] and is a good example of a sea-derived peptide. According to the National Cancer Institute (NCI) Thesaurus, Tasidotin is a third-generation, synthetic, water-soluble, pentapeptide analog of the marine depsipeptide Dolastatin 15, with potential anti-mitotic and anti-neoplastic activities.
When one is working with Computational Chemistry and Molecular Modeling, there is not a universal methodology that could be applied to the entire spectra of known or unknown molecular systems. As such, our own devised methodology cannot be considered useful for all systems. The workaround that researchers in this field have found relies on studying, with a particular methodology, different but related families of molecules, in order to see if the proposed methodology can be applied with the same degree of success to the different group of molecules, and that the results obtained with a particular family of molecules are good not just by chance or some kind of serendipity. Indeed, the larger the number of studies with different groups of molecules, the greater the validity of the used methodologies.
Thus, the objective of this work is to perform a comparative study of the chemical reactivity of the Soblidotin and Tasidotin peptides, of marine origin, by resorting to Chemical Reactivity Theory (CRT) [2], in addition to an investigation of their potential AGEs inhibition abilities, following a protocol devised by our group [3,4,5,6,7,8,9,10,11,12]. Moreover, the bioactivity scores for Soblidotin and Tasidotin are predicted through an empirical procedure, based on comparison with molecular structures with well-known pharmacological properties.

2. Computational Methodology

Rather than resorting to the usual methodologies in Medicinal Chemistry, related to Molecular Docking, QSAR, and QSPR, we prefer to study the potential usefulness of the peptides by considering Computational Chemistry and Molecular Modelling, based on CRT and some semiempirical relationships relying on collected databases.
As with our recent studies [3,4,5,6,7,8,9,10], the computational tasks in this work have been done by considering the popular Gaussian 09 software [13]. Following the conclusions obtained from those studies, the MN12SX density functional [14] is chosen, because it can be considered to be well-behaved; according to our proposed, above-mentioned criteria. Accordingly, the calculation of the electronic properties used a model chemistry, based on the mentioned density functional in connection with the Def2TZVP basis set, while a smaller Def2SVP was considered for the prediction of the most stable structures [15,16]. In order to obtain accurate results, all calculations were performed using water as the solvent, simulated with the SMD model [17].

3. Results and Discussion

The molecular structures of the Soblidotin and Tasidotin peptides, which are depicted in Figure 1, were optimized in absence of solvent by resorting to the semiempirical DFTBA model (available in Gaussian 09), starting from the five most stable conformers selected from a pre-optimization accomplished by means of Molecular Mechanics techniques [18,19,20,21,22]. The resulting conformers were processed, as is customary within Computational Chemistry, to obtain the desired calculated properties, according to the techniques mentioned in the previous section.
The calculation of the maximum wavelength absorption of the Soblidotin and Tasidotin peptides was performed by conducting ground-state determinations with the same model chemistry. These results, together with the total, HOMO, and LUMO electronic energies, as well as the HOMO–LUMO gap, are shown in Table 1.
To get a glimpse of the electronic distribution around the atoms of both molecules, graphical sketches of the total electronic densities for Soblidotin and Tasidotin are presented in Figure 2.

3.1. Calculation of Global Reactivity Descriptors

The definitions for the chemical reactivity descriptors, coming from CRT, are provided in terms of the energies of the HOMO and LUMO, ϵ H and ϵ L , respectively, as: The electronegativity χ 1 2 ( ϵ L + ϵ H ) [23,24], the global hardness η ( ϵ L ϵ H ) [23,24], the electrophilicity ω = ( ϵ L + ϵ H ) 2 / 4 ( ϵ L ϵ H ) [25], the electrodonating power ω ( 3 ϵ H + ϵ L ) 2 / 16 η [26], the electroaccepting power ω + ( ϵ H + 3 ϵ L ) 2 / 16 η [26], and the net electrophilicity Δ ω ± = ω + + ω [27]. The results of their calculations for the Soblidotin and Tasidotin peptides are presented in Table 2.
We can see, from Table 2, that, for both peptides, their electronating ability is more important that their electroaccepting ability, although there was not a significant difference between the values of the global reactivity descriptors of each of them.

3.2. Determination of the Local Reactivity Descriptors

If we now focus on the local reactivity descriptors coming from CRT, then the definitions will be: f + ( r ) = ρ N + 1 ( r ) ρ N ( r ) for the Nucleophilic Fukui Function [23,24], f ( r ) = ρ N ( r ) ρ N 1 ( r ) for the Electrophilic Fukui Function [23,24], and Δ f ( r ) = f ( r ) / N υ ( r ) for the Dual Descriptor [2,28,29,30,31,32], considering the electronic densities for a system with N + 1 , N, and N 1 electrons, respectively.
The Electrophilic Fukui functions f ( r ) and Nucleophilic Fukui functions f + ( r ) for the anticancer marine peptides Soblidotin and Tasidotin are shown in Figure 3.

3.3. Calculation of the pKas of the Antimicrobial Peptides of Marine Origin and Quantification of the AGEs Inhibitor Ability

For the consideration of natural products to be useful for applications in Medicinal Chemistry, the knowledge of their pKas is very important. For large molecules, such as the peptides studied in this work, the experimental determination of their pKa is sometimes difficult. However, we have already shown that the pKa of a peptide can be predicted by applying a relationship with the global hardness η [11]. Thus, the values of the pKas of Soblidotin and Tasidotin, determined through that relationship, are presented in Table 3:
It is our belief that these results could be of interest for processes, involving drug design and development, using these peptides.
The Maillard reaction is, in fact, a collection of chained chemical reactions that start from the interaction between a reducing carbonyl and the free amino group of a peptide or protein. This chain of reactions leads to a final group of molecules, known as Advanced Glycation Endproducts (or AGEs), which are considered to be the main reasons for the development of some diseases, such as Diabetes, Alzheimers, and Parkinsons [33].
Many compounds have been devised as drugs with the goal of inhibiting the formation of AGEs, which include Pyridoxamine, Aminoguanidine, Carnosine, Metformin, Pioglitazone, and Tenilsetam [34,35]. Previous results from our group on the glycation process led us to conclude that the nucleophilicity N [36] may be a good descriptor for the AGE inhibition ability, and a nice definition for it could be based on the the inverse of the net electrophilicity Δ ω ± .
If we consider that this relationship is valid, then, based on the calculated global descriptors for Soblidotin and Tasidotin, it should be possible to predict their AGE inhibition ability, by comparison with the mentioned drugs:
ALT-946 > Aminoguanidine > Metformin > Carnosine > Tasidotin >
> Soblidotin > Tenilsetam > Pyridoxamine > > Pioglitazone
It can be seen that the studied peptides possess AGE inhibitor abilities, which are lower than that of other molecular systems already on the market [34,35].

3.4. Bioactivity Scores

When looking for the potential therapeutic capacities of new molecules, an interesting approach is an empirical procedure based on comparison with the properties of molecular structures with well-known pharmacological properties. These values are known as Bioactivity Scores, and can be easily estimated by resorting to the free online Molinspiration software, as shown in Table 4.
The results from Table 4 can be interpreted by considering the factors for evaluation of the Bioactivity Scores, which are shown in Table 5.

4. Conclusions

The outcomes of this work have resulted in the determination of important information about the chemical reactivity of two peptides of marine origin with potential therapeutic properties as anticancer drugs, Soblidotin and Tasidotin. This has been done by resorting to Chemical Reactivity Theory and some empirical relationships for the prediction of the pKa values, for their AGE inhibition abilities, and for the prediction of their Bioactivity Scores (as a measure of bioactivity).
The total electronic densities for both peptides have been displayed in the form of graphical sketches, as well as the HOMO and LUMO orbitals, as an approximation to the electrophilic and nucleophilic Fukui functions.
All these computed values could be of interest in the design and development of new pharmaceutical drugs, as well as for the practice of Medicinal Chemistry.

Author Contributions

D.G.-M. conceived and designed the research and headed, wrote, and revised the manuscript, while J.F. and N.F.-H. contributed to the writing and the revision of the article.

Funding

This research was funded by Consejo Nacional de Ciencia y Tecnología (CONACYT, Mexico) through Grant 219566-2014 for Basic Science Research and by the Ministerio de Economía y Competitividad (MINECO) and the European Fund for Regional Development (FEDER) (CTQ2014-55835-R).

Acknowledgments

Daniel Glossman-Mitnik conducted this work as a Visiting Lecturer at the University of the Balearic Islands, from which support is gratefully acknowledged. Norma Flores-Holguín and Daniel Glossman-Mitnik are researchers of CIMAV and CONACYT.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mohammad, R.; Limvarapuss, C.; Wall, N.; Hamdy, N.; Beck, F.; Pettit, G.; Al-Katib, A. A new tubulin polymerization inhibitor, auristatin PE, induces tumor regression in a human waldenstrom’s macroglobulinemia xenograft model. Int. J. Oncol. 1999, 15, 367–439. [Google Scholar] [CrossRef] [PubMed]
  2. Toro-Labbé, A. (Ed.) Theoretical Aspects of Chemical Reactivity; Elsevier Science: Amsterdam, The Netherlands, 2007. [Google Scholar]
  3. Frau, J.; Glossman-Mitnik, D. Molecular reactivity and absorption properties of melanoidin blue-G1 through conceptual DFT. Molecules 2018, 23, 559. [Google Scholar] [CrossRef] [PubMed]
  4. Frau, J.; Glossman-Mitnik, D. Conceptual DFT study of the local chemical reactivity of the dilysyldipyrrolones A and B intermediate melanoidins. Theor. Chem. Acc. 2018, 137, 1210. [Google Scholar] [CrossRef]
  5. Frau, J.; Glossman-Mitnik, D. Conceptual DFT study of the local chemical reactivity of the colored BISARG melanoidin and its protonated derivative. Front. Chem. 2018, 6, 1–9. [Google Scholar] [CrossRef] [PubMed]
  6. Frau, J.; Glossman-Mitnik, D. Molecular reactivity of some maillard reaction products studied through conceptual DFT. Contemp. Chem. 2018, 1, 1–14. [Google Scholar]
  7. Frau, J.; Glossman-Mitnik, D. Computational study of the chemical reactivity of the blue-M1 intermediate melanoidin. Comput. Theor. Chem. 2018, 1134, 22–29. [Google Scholar] [CrossRef]
  8. Frau, J.; Glossman-Mitnik, D. Chemical reactivity theory applied to the calculation of the local reactivity descriptors of a colored maillard reaction product. Chem. Sci. Int. J. 2018, 22, 1–14. [Google Scholar] [CrossRef]
  9. Frau, J.; Glossman-Mitnik, D. Blue M2: An intermediate melanoidin studied via conceptual DFT. J. Mol. Model. 2018, 24, 1–13. [Google Scholar] [CrossRef]
  10. Frau, J.; Flores-Holguín, N.; Glossman-Mitnik, D. Chemical reactivity properties, pKa values, AGEs inhibitor abilities and bioactivity scores of the mirabamides A–H peptides of marine origin studied by means of conceptual DFT. Mar. Drugs 2018, 16, 302. [Google Scholar] [CrossRef]
  11. Frau, J.; Hernández-Haro, N.; Glossman-Mitnik, D. Computational prediction of the pKas of small peptides through conceptual DFT descriptors. Chem. Phys. Lett. 2017, 671, 138–141. [Google Scholar] [CrossRef]
  12. Frau, J.; Glossman-Mitnik, D. Chemical reactivity theory study of advanced glycation endproduct inhibitors. Molecules 2017, 22, 226. [Google Scholar] [CrossRef]
  13. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09 Revision E.01; Gaussian Inc.: Wallingford, UK, 2016. [Google Scholar]
  14. Peverati, R.; Truhlar, D.G. Screened-exchange density functionals with broad accuracy for chemistry and solid-state physics. Phys. Chem. Chem. Phys. 2012, 14, 16187–16191. [Google Scholar] [CrossRef]
  15. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef]
  16. Weigend, F. Accurate Coulomb-fitting basis sets for H to R. Phys. Chem. Chem. Phys. 2006, 8, 1057–1065. [Google Scholar] [CrossRef]
  17. Marenich, A.; Cramer, C.; Truhlar, D. Universal solvation model based on solute electron density and a continuum model of the solvent defined by the bulk dielectric constant and stomic surface tensions. J. Phys. Chem. B 2009, 113, 6378–6396. [Google Scholar] [CrossRef]
  18. Halgren, T.A. Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J. Comput. Chem. 1996, 17, 490–519. [Google Scholar] [CrossRef]
  19. Halgren, T.A. Merck molecular force field. II. MMFF94 van der Waals and electrostatic parameters for intermolecular interactions. J. Comput. Chem. 1996, 17, 520–552. [Google Scholar] [CrossRef]
  20. Halgren, T.A. MMFF VI. MMFF94s option for energy minimization studies. J. Comput. Chem. 1999, 20, 720–729. [Google Scholar] [CrossRef]
  21. Halgren, T.A.; Nachbar, R.B. Merck molecular force field. IV. Conformational energies and geometries for MMFF94. J. Comput. Chem. 1996, 17, 587–615. [Google Scholar] [CrossRef]
  22. Halgren, T.A. Merck molecular force field. V. Extension of MMFF94 using experimental data, additional computational data, and empirical rules. J. Comput. Chem. 1996, 17, 616–641. [Google Scholar] [CrossRef]
  23. Parr, R.; Yang, W. Density-Functional Theory of Atoms and Molecules; Oxford University Press: New York, NY, USA, 1989. [Google Scholar]
  24. Geerlings, P.; De Proft, F.; Langenaeker, W. Conceptual density functional theory. Chem. Rev. 2003, 103, 1793–1873. [Google Scholar] [CrossRef]
  25. Parr, R.; Szentpaly, L.; Liu, S. Electrophilicity index. J. Am. Chem. Soc. 1999, 121, 1922–1924. [Google Scholar] [CrossRef]
  26. Gázquez, J.; Cedillo, A.; Vela, A. Electrodonating and electroaccepting powers. J. Phys. Chem. 2007, 111, 1966–1970. [Google Scholar] [CrossRef]
  27. Chattaraj, P.; Chakraborty, A.; Giri, S. Net electrophilicity. J. Phys. Chem. 2009, 113, 10068–10074. [Google Scholar] [CrossRef]
  28. Morell, C.; Grand, A.; Toro-Labbé, A. New dual descriptor for chemical reactivity. J. Phys. Chem. 2005, 109, 205–212. [Google Scholar] [CrossRef]
  29. Morell, C.; Grand, A.; Toro-Labbé, A. Theoretical support for using the Δf(r) descriptor. Chem. Phys. Lett. 2006, 425, 342–346. [Google Scholar] [CrossRef]
  30. Martínez-Araya, J.I. Revisiting Caffeate’s capabilities as a complexation agent to silver cation in mining processes by means of the dual descriptor—A conceptual DFT approach. J. Mol. Model. 2012, 18, 4299–4307. [Google Scholar] [CrossRef]
  31. Martínez-Araya, J.I. Explaining reaction mechanisms using the dual descriptor: A complementary tool to the molecular electrostatic potential. J. Mol. Model. 2012, 19, 2715–2722. [Google Scholar] [CrossRef]
  32. Martínez-Araya, J.I. Why is the dual descriptor a more accurate local reactivity descriptor than Fukui functions? J. Math. Chem. 2015, 53, 451–465. [Google Scholar] [CrossRef]
  33. Ahmed, N. Advanced glycation endproducts—Role in pathology of diabetic complications. Diabetes Res. Clin. Pract. 2005, 67, 3–21. [Google Scholar] [CrossRef]
  34. Rahbar, S.; Figarola, J.L. Novel inhibitors of advanced glycation endproducts. Arch. Biochem. Biophys. 2003, 419, 63–79. [Google Scholar] [CrossRef]
  35. Peyroux, J.; Sternberg, M. Advanced glycation endproducts (AGEs): Pharmacological inhibition in diabetes. Pathol. Biol. 2006, 54, 405–419. [Google Scholar] [CrossRef]
  36. Domingo, L.R.; Perez, P. The nucleophilicity N index in organic chemistry. Org. Biomol. Chem. 2011, 9, 7168–7175. [Google Scholar] [CrossRef]
Figure 1. Graphical sketches of the molecular structures of the anticancer marine peptides: (a) Soblidotin and (b) Tasidotin.
Figure 1. Graphical sketches of the molecular structures of the anticancer marine peptides: (a) Soblidotin and (b) Tasidotin.
Molecules 24 01115 g001
Figure 2. Graphical sketches of the total electronic densities of the anticancer marine peptides: (a) Soblidotin and (b) Tasidotin.
Figure 2. Graphical sketches of the total electronic densities of the anticancer marine peptides: (a) Soblidotin and (b) Tasidotin.
Molecules 24 01115 g002
Figure 3. Graphical sketches of the HOMO and LUMO of the Soblidotin and Tasidotin marine peptides as approximations to the electrophilic (left column) and nucleophilic (right column) Fukui functions.
Figure 3. Graphical sketches of the HOMO and LUMO of the Soblidotin and Tasidotin marine peptides as approximations to the electrophilic (left column) and nucleophilic (right column) Fukui functions.
Molecules 24 01115 g003
Table 1. Electronic energies of the neutral molecular systems (in au) of the anticancer marine peptides, Soblidotin and Tasidotin; the HOMO and LUMO orbital energies, as well as the HOMO–LUMO gap (in eV), and the maximum absorption wavelengths λ m a x (in nm).
Table 1. Electronic energies of the neutral molecular systems (in au) of the anticancer marine peptides, Soblidotin and Tasidotin; the HOMO and LUMO orbital energies, as well as the HOMO–LUMO gap (in eV), and the maximum absorption wavelengths λ m a x (in nm).
MoleculeTotal Electronic EnergyHOMOLUMOHOMO–LUMO Gap λ max
Soblidotin−2250.5737−6.3060−0.73535.5707223
Tasidotin−1958.1339−6.2790−0.64905.6300220
Table 2. Global reactivity descriptors of the anticancer marine peptides Soblidotin and Tasidotin.
Table 2. Global reactivity descriptors of the anticancer marine peptides Soblidotin and Tasidotin.
MoleculeElectronegativityGlobal HardnessElectrophilicity
Soblidotin3.52065.57071.1125
Tasidotin3.46405.63001.0657
MoleculeElectrodonating PowerElectroaccepting PowerNet Electrophilicity
Soblidotin4.33340.81285.1463
tasidotin4.21520.75124.9964
Table 3. pKa values for the anticancer marine peptides Soblidotin and Tasidotin.
Table 3. pKa values for the anticancer marine peptides Soblidotin and Tasidotin.
MoleculepKa
Soblidotin11.70
Tasidotin11.65
Table 4. Empirical Bioactivity Scores of the anticancer marine peptides, Soblidotin and Tasidotin.
Table 4. Empirical Bioactivity Scores of the anticancer marine peptides, Soblidotin and Tasidotin.
MoleculeGPCR LigandIon Channel ModulatorKinase InhibitorNuclear Receptor LigandProtease InhibitorEnzyme Inhibitor
Soblidotin−0.30−1.30−1.07−1.340.19−0.83
Tasidotin0.19−0.43−0.230.360.51−0.13
Table 5. Factors for Evaluation of the Empirical Bioactivity Scores.
Table 5. Factors for Evaluation of the Empirical Bioactivity Scores.
ActiveModerately ActiveInactive
>0Between −5.0 and 0.0<−5.0

Share and Cite

MDPI and ACS Style

Frau, J.; Flores-Holguín, N.; Glossman-Mitnik, D. Chemical Reactivity Theory and Empirical Bioactivity Scores as Computational Peptidology Alternative Tools for the Study of Two Anticancer Peptides of Marine Origin. Molecules 2019, 24, 1115. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules24061115

AMA Style

Frau J, Flores-Holguín N, Glossman-Mitnik D. Chemical Reactivity Theory and Empirical Bioactivity Scores as Computational Peptidology Alternative Tools for the Study of Two Anticancer Peptides of Marine Origin. Molecules. 2019; 24(6):1115. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules24061115

Chicago/Turabian Style

Frau, Juan, Norma Flores-Holguín, and Daniel Glossman-Mitnik. 2019. "Chemical Reactivity Theory and Empirical Bioactivity Scores as Computational Peptidology Alternative Tools for the Study of Two Anticancer Peptides of Marine Origin" Molecules 24, no. 6: 1115. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules24061115

Article Metrics

Back to TopTop