Next Article in Journal
Essential Oil Composition and Bioactivities of Waldheimia glabra (Asteraceae) from Qinghai-Tibet Plateau
Next Article in Special Issue
A Molecular Electron Density Theory Study of the Reactivity of Azomethine Imine in [3+2] Cycloaddition Reactions
Previous Article in Journal
Evaluation of Maltose-Based Cationic Liposomes with Different Hydrophobic Tails for Plasmid DNA Delivery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Prediction of the Protonation Sites of Ac-Lys-(Ala)n-Lys-NH2 Peptides through Conceptual DFT Descriptors

by
Sebastián Sastre
1,
Juan Frau
1 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, Centro de Investigación en Materiales Avanzados, Departamento de Medio Ambiente y Energía, Chihuahua, Chih 31136, Mexico
*
Author to whom correspondence should be addressed.
Molecules 2017, 22(3), 458; https://doi.org/10.3390/molecules22030458
Submission received: 6 February 2017 / Accepted: 10 March 2017 / Published: 13 March 2017

Abstract

:
Six density functionals (M11, M11L, MN12L, MN12SX, N12, and N12SX) in connection with the Def2TZVP basis set and the SMD solvation model (water as a solvent) have been assessed for the calculation of the molecular structure and properties of several peptides with the general formula Ac-Lys-(Ala) n -Lys-NH2, with n = 0 to 5. The chemical reactivity descriptors for these systems have been calculated through Conceptual density functional theory (DFT). The active sites for nucleophilic and electrophilic attacks have been chosen by relating them to the Fukui function indices, the condensed dual descriptor Δf ( r ) , and the electrophilic Parr functions. The results allowed the prediction of the protonation sites of the peptides and rendered a qualitative explanation of the difference in pK a of the two Lys groups.

Graphical Abstract

1. Introduction

Conceptual density functional theory (DFT, also known as chemical reactivity theory) is a powerful tool for the prediction, analysis, and interpretation of the outcome of chemical reactions [1,2,3,4].
Recently, Luis R. Domingo [5] proposed a new theory for the study of the reactivity in organic chemistry, which he named molecular electron density theory (MEDT). In this new theory, the capability for changes in electron density is responsible for the molecular reactivity [5], while the electron density distribution at the ground state is responsible for physical and chemical molecular properties.
Within MEDT (which encompasses conceptual DFT), several new concepts and reactivity descriptors have been defined, such as the global electron density transfer (GEDT), the nucleophilicity N index, and local condensed descriptors like the electrophilic P k - and nucleophilic P k + Parr functions [6]. The use of the Parr functions is becoming widespread for the understanding of organic chemistry processes, as evidenced by some recent publications related to Diels–Alder and cycloadditions reactions [7,8,9].
An interesting chemical reaction amenable to being studied through conceptual DFT is protein glycation. This reaction is initiated by a nucleophilic addition (nonenzymatic glycation or Maillard reaction), a reaction between a free amino group from a protein and a carbonyl group from a reducing sugar to form a freely reversible Schiff base. Glycated proteins can undergo further reactions, giving rise to some structures called advanced glycation endproducts (AGEs) [10].
Indeed, the glycation sites of an amino acid, peptide, or protein may also be the preferred sites for protonation. The reason for this assertion is because from a chemical point of view, the interaction of the carbon atom of a reducing carbonyl compound with an amino nitrogen atom is the same kind of reaction that occurs when a proton interacts with that nitrogen atom of the amino group. For a peptide or protein with many glycation or protonation sites, there is a rank of preferred reaction sites that depends on the pH, the solvent, and the conformational structure of the molecular system. This means that that if we know the ordering of the protonation of the peptide or protein, it will be possible to discern (at least qualitatively) the pK a of the Lys residues.
In some recent interesting works [11,12], Makowska et al. studied the acidic–basic properties of alanine-based peptides containing acidic and basic side chains. In particular, they experimentally estimated the influence of the length of the alanine spacer on the pK a of the Lys residues of Ac-Lys-(Ala)n-Lys-NH2 peptides (wth n = 0, 1, 2, , 5) [11].
Following the pioneering work of Parr et al. [1], a number of useful concepts have been derived from the analysis of the density of any molecular system through DFT. These concepts that allow a researcher to make qualitative predictions about the chemical reactivity of a given system can also be quantified, and are collectively known as conceptual DFT Descriptors.
Therefore, we believe that it could be of interest to apply the concepts of density functional theory to the study of those peptides in order to determine if it is possible to discern between the potential protonation sites and to obtain a qualitative estimation of the pK a differences between them.
In order to obtain quantitative values of the conceptual DFT descriptors, it is necessary to resort to the Kohn–Sham theory through calculations of the molecular density, the energy of the system, and the orbital energies—specifically those related to the frontier orbitals; that is, HOMO and LUMO [13,14,15,16,17,18].
Although the foundations of DFT have established that a universal density functional must exist and that all of the properties of the system can be obtained through calculations with this functional, in practice one must resort to some of the approximate density functionals that have been developed in the last thirty years. Because these are approximate functionals (i.e., not universal functionals), many of them are good for predicting some properties and others are good for other properties. Sometimes you can find density functionals that are excellent for describing the properties of a given molecular system with a particular functional group, but it is necessary to resort to other density functionals for a different functional group that you want to include in the molecular system under study.
When one is dealing with the study of chemical reactivity (i.e., a process that involves the transference of electrons), it is normal to perform calculations not only of the ground state, but also for open systems like the radical cation and radical anion. These systems are often difficult to converge giving trustworthy results—specially if diffuse functions must be included in the basis set [13,14,15,16,17,18]. For this reason, it is convenient to have a method that can give all information that one needs directly from the results of the calculation of the ground state of the molecular system under study. In particular, one may want to obtain the ionization potential (I) and electron affinity (A) of the system avoiding the calculation of the anion and cation radicals. Indeed, the link for this is given by the so-called Koopmans’ theorem [15,16,17,18].
Koopmans’ theorem is not valid within DFT. However, from an empirical and practical point of view, it is meaningful to follow the procedure of assigning KS HOMO as equal to and opposite of the vertical ionization potential, ϵ H = –I, and the KS LUMO as equal to and opposite of the vertical electron affinity, ϵ L = –A. We have coined the acronym KID (“Koopmans in DFT”) for this empirical procedure.
This means that the goodness of a given density functional for the purpose of predicting the Conceptual DFT descriptors directly from the properties of the neutral molecule can be estimated by checking how well it follows the KID procedure. Thus, it will be interesting to consider several recent density functionals that have shown great accuracy across a broad spectrum of databases in chemistry and physics [19] to evaluate their performance in fulfilling this practical technique.
The objective of this work is twofold: (i) to conduct a comparative study of the performance of some density functionals from the Minnesota family for reproducing chemical reactivity descriptors of the Ala-spaced Lys-containing peptides proposed by Makowska et al. [12] within the KID formalism; and (ii) to predict the preferred protonation sites that could be give a qualitative estimation of the pK a of the Lys residues using condensed dual descriptors.

2. Theoretical Background

As this work is part of an ongoing project, the theoretical background is similar to that presented in previous research [20,21,22,23,24,25,26,27] and will be shown here for the sake of completeness. Within the conceptual framework of DFT [2,28], the chemical potential μ is defined as:
μ = E N v ( r ) = - χ ,
where χ is the electronegativity, while the global hardness η is
η = 2 E N 2 v ( r ) .
Using a finite difference approximation and the KID procedure, the above expressions can be written as:
μ = - 1 2 ( I + A ) 1 2 ( ϵ L + ϵ H ) = χ K ,
η = ( I - A ) ( ϵ L - ϵ H ) = η K ,
where ϵ H and ϵ L are the energies of the highest occupied and the lowest unoccupied molecular orbitals (HOMO and LUMO), respectively. The electrophilicity index ω has been defined as [29]:
ω = μ 2 2 η = ( I + A ) 2 4 ( I - A ) ( ϵ L + ϵ H ) 2 4 ( ϵ L - ϵ H ) = ω K .
The Fukui function is defined in terms of the derivative of ρ ( r ) with respect to N [2]:
f ( r ) = ρ ( r ) N υ ( r )
The function f ( r ) reflects the ability of a molecular site to accept or donate electrons. High values of f ( r ) are related to a high reactivity at point r [2]. By applying a finite difference approximation to Equation (6), two definitions of Fukui functions depending on total electronic densities are obtained:
f + ( r ) = ρ N + 1 ( r ) - ρ N ( r ) ,
f - ( r ) = ρ N ( r ) - ρ N - 1 ( r ) ,
where ρ N + 1 ( r ) , ρ N ( r ) , and ρ N - 1 ( r ) are the electronic densities at point r for the system with N + 1 , N, and N - 1 electrons, respectively. The first one, f + ( r ) , has been associated to reactivity for a nucleophilic attack, so it measures the intramolecular reactivity at the site r towards a nucleophilic reagent. The second one, f - ( r ) , has been associated to reactivity for an electrophilic attack, so this function measures the intramolecular reactivity at the site r toward an electrophilic reagent [28].
Morell et al. [3,30,31,32,33,34,35] proposed a local reactivity descriptor (LRD) which is called the dual descriptor (DD) f ( 2 ) ( r ) f ( r ) . The definition of f ( r ) is shown as indicated by Morell et al. [30,31]:
f ( r ) = f ( r ) N υ ( r ) .
The dual descriptor can be condensed over the atomic sites: when f k > 0 , the process is driven by a nucleophilic attack on atom k, and then that atom acts as an electrophilic species; conversely, when f k < 0 , the process is driven by an electrophilic attack over atom k, and therefore atom k acts as a nucleophilic species.
In an interesting and very enlightening article [36], Jorge Ignacio Martínez-Araya has shown why the dual descriptor is a more accurate local reactivity descriptor than the Fukui functions. He has shown that although the Fukui function has the capability of revealing nucleophilic and electrophilic regions on a molecule, the dual descriptor is able to unambiguously expose truly nucleophilic and electrophilic regions, but along with the latter, dual descriptor is less affected by the lack of relaxation terms than the Fukui function when the frontier molecular orbital approximation is applied, and that this implies that the dual descriptor can be considered a more reliable descriptor to measure local reactivity than Fukui function. For this reason, we have chosen the DD to fulfil the objectives of our study instead of the Fukui functions individually.
In 2013, Domingo proposed the Parr functions P( r ) [37,38], which are given by the following equations:
P - ( r ) = ρ s r c ( r ) ( for electrophilic attacks ) ,
P + ( r ) = ρ s r a ( r ) ( for nucleophilic attacks ) ,
which are related to the atomic spin density (ASD) at the r atom of the radical cation or anion of a given molecule, respectively. The ASD over each atom of the radical cation and radical anion of the molecule gives the local nucleophilic P k - and electrophilic P k + Parr functions of the neutral molecule [6].

3. Settings and Computational Methods

Following the lines of our previous work [20,21,22,23,24,25,26,27], the computational studies were performed with the Gaussian 09 [39] series of programs with density functional methods as implemented in the computational package. The equilibrium geometries of the molecules were determined by means of the gradient technique. The force constants and vibrational frequencies were determined by computing analytical frequencies on the stationary points obtained after the optimization to check if there were true minima. The basis set used in this work was Def2SVP for geometry optimization and frequencies, while Def2TZVP was considered for the calculation of the electronic properties [40,41].
For the calculation of the molecular structure and properties of the studied systems, we chose six density functionals from the latest Minnesota density functionals family which consistently provide satisfactory results for several structural and thermodynamic properties [19]: M11, which is a is a range-separated hybrid meta-GGA (generalized gradient approximation) [42]; M11L, which is a dual-range local meta-GGA [43]; MN12L, which is a nonseparable local meta-NGA (nonseparable gradient approximation) [44]; MN12SX, which is a range-separated hybrid nonseparable meta-NGA [45]; N12, which is a nonseparable gradient approximation [46]; and N12SX, which is a range-separated hybrid nonseparable gradient approximation [45]. In these functionals, GGA stands for generalized gradient approximation (in which the density functional depends on the up and down spin densities and their reduced gradient) and NGA stands for nonseparable gradient approximation (in which the density functional depends on the up/down spin densities and their reduced gradient, and also adopts a nonseparable form). All the calculations were performed in the presence of water as a solvent by doing Integral Equation Formalism–Polarized Continuum Model (IEF-PCM) computations according to the solvation model density (SMD) solvation model [47].

4. Results and Discussion

The molecular structures of the six peptides with the general formula Ac-Lys-(Ala) n -Lys-NH2, with n = 0 to 5 were built by means of the Avogadro 1.2.0 program [48,49]. The resulting structures were pre-optimized by finding the most stable conformers through a random sampling with molecular mechanics techniques and a consideration of all the torsional angles through a random sampling with molecular mechanics techniques and a consideration of all the torsional angles through the general AMBER force field [50]. The structures of the resulting conformers were then reoptimized with the six density functionals mentioned in the previous section in conjunction with the Def2SVP basis set and the SMD solvation model, using water as a solvent. As an illustration of the optimized peptides, the geometrical conformation of Ac-Lys-(Ala) n -Lys-NH2 with n = 1 (KAK) is presented in Figure 1.
The HOMO and LUMO orbital energies (in eV), vertical ionization potentials I, and electron affinities A (in eV), global electronegativity χ, total chemical hardness η, and global electrophilicity ω of the peptides calculated with the six density functionals and the Def2TZVP basis set using water as as solvent simulated with the SMD parametrization of the IEF-PCM model are presented in Tables S1A to S6A of the Electronic Supplementary Information (ESI). The upper part of the tables shows the results derived assuming the KID procedure (hence the subscript K) and the lower part shows the results derived from the calculated ΔSCF vertical energies.
With the objective of analyzing our results in order to verify the fulfillment of the KID procedure, we previously designed several descriptors that relate the results obtained through the HOMO and LUMO calculations with those obtained by means of the vertical I and A with a ΔSCF procedure. The first three descriptors are related to the simplest fulfillment of Koopmans’ theorem by relating ϵ H with –I, ϵ L with –A, and the behavior of them in the description of the band gap: J I = | ϵ H + E g s ( N - 1 ) - E g s ( N ) | , J A = | ϵ L + E g s ( N ) - E g s ( N + 1 ) | , and J H L = J I 2 + J A 2 .
Next, we consider four other descriptors that analyze how well the studied density functionals are useful for the prediction of the electronegativity χ, the global hardness η, and the global electrophilicity ω, and for a combination of these conceptual DFT descriptors, considering only the energies of the HOMO and LUMO or the vertical I and A: J χ = | χ - χ K | , J η = | η - η K | , J ω = | ω - ω K | , and J D 1 = J χ 2 + J η 2 + J ω 2 , where D1 stands for the first group of conceptual DFT descriptors.
The results of the calculations of J I , J A , J H L , J χ , J η , J ω , and J D 1 for the six peptides considered in this work are displayed in Tables S1B to S6B of the Electronic Supplementary Information (ESI).
Based on the results for the descriptors presented on Tables S1B to S6B of the ESI, we have compiled the average values for each density functional on the whole group of peptides, and the calculated results are displayed in Table 1.
As can be seen from Table 1, the KID procedure holds with great accuracy for the MN12SX and N12SX density functionals, which are range-separated hybrid meta-NGA and range-separated hybrid NGA density functionals, respectively. Indeed, the values of J I , J A , and J H L are not exactly zero. However, the results are impressive, particularly for the N12SX density functional.
It is interesting to see that the same density functionals also fulfill the KID procedure for the other descriptors—namely J χ , J η , J ω , and J D 1 . These results are very important, because they show that it is not enough to rely only on J I , J A , and J H L . For example, if we consider only J χ , for all of the density functionals studied in this work, the values are very close to zero. As for the other descriptors, only the MN12SX and N12SX density functionals show this behavior. This means that the results for J χ could be due to a fortuitous cancellation of errors.
An important fact is that although the range-separated hybrid NGA and range-separated hybrid meta-NGA density functionals can be useful for the calculation of the conceptual DFT descriptors, it is not the same for the range-separated hybrid GGA (M11) density functional. An inspection of Table S1A of the ESI shows that this is due to the fact that this functional inadequately describes the energy of the LUMO, leading to negative values of A, which are in contradiction with the ΔSCF results.
The condensed Fukui functions and condensed dual descriptor have been calculated using the AOMix molecular analysis program [51,52] starting from single-point energy calculations.
The condensed dual descriptors f k and the electrophilic Parr functions P k - (mpa) and P k - (mpa) over the N atoms of the Lys residues of the KK, KAK, KA2K, KA3K, KA4K, and KA5K peptides of the Lys residues of the KK, KAK, KA2K, KA3K, KA4K, and KA5K peptides calculated with the MN12SX and N12SX density functionals and the Def2TZVP basis set using water as as solvent simulated with the SMD parametrization of the IEF-PCM model are shown in Table 2.
As the studied peptides are end-capped, it would be expected that the chemical reactivities of the Lys nitrogen atoms will be different. In fact, the values presented in Table 2 correspond to the condensed dual descriptors and electrophilic Parr functions over the nitrogen atoms of the Lys closer to the end-terminal Ac- group (Lys1). The values of the descriptors over the nitrogen atom of the other Lys residue (Lys2) are not shown because they are almost zero in all cases. This means that the preferred glycation site will be Lys1, but it will also be the preferred protonation site. Therefore, it can be predicted that the pK a of Lys 1 will be greater than the pK a of Lys2, although it is not possible to assign a particular value to both pK a s. Moreover, although there is not a clear trend, it can be observed from the results in Table 2 for the condensed dual descriptor and electrophilic Parr functions (either from a Mulliken population analysis (MPA) or Hirshfeld population analysis (HPA)) that the value of the descriptors are larger for KA4K and KA5K than for the other peptides. Thus, it can be concluded that the effect of the alanine spacers will be that of increasing the difference between the pK a of Lys1 and and that of Lys2. This can probably be explained by considering that the increase of the number of Ala units leads to geometrical conformations of the peptides, making the terminal Lys residues be apart from each other and avoiding the competition between them for the reacting proton H + .

5. Conclusions

The Minnesota family of density functionals (M11, M11L, MN12L, MN12SX, N12, N12SX) has been tested for the fulfillment of the empirical KID procedure by comparison of the HOMO- and LUMO-derived values with those obtained through a ΔSCF technique. It has been shown that the range-separated hybrid meta-NGA density functional (MN12SX) and the range-separated hybrid NGA density functional (N12SX) are the best for the accomplishment of this objective. As such, they are a good alternative to those density functionals whose behavior have been tuned through a gap-fitting procedure and a good prospect for their usefulness in the description of the chemical reactivity of molecular systems of large size.
From the whole of the results presented in this work, it can be seen that the sites of interaction of the studied peptides can be predicted by using DFT-based reactivity descriptors such as the condensed dual descriptor and Parr functions calculations. These descriptors were used in the characterization and successful description of the preferred reactive sites and provide a firm explanation for the reactivity of those molecules.
Moreover, the preferred glycation and protonation sites could be identified, and helped to understand the effect of the alanine spacer on such properties as the pK a of the different Lys residues. As a result, the difference between the pK a of the Lys residue closer to the Ac- end-capped site and the pK a of the other Lys residue will be always greater than zero, and will be increased as the number of alanine spacers grows.

Supplementary Materials

The following are available online. ESI.pdf : Electronic Supplementary Information.

Acknowledgments

This work has been partially supported by CIMAV, SC and Consejo Nacional de Ciencia y Tecnología (CONACYT, Mexico) through Grant 219566/2014 for Basic Science Research and Grant 265217/2016 for a Foreign Sabbatical Leave. Daniel Glossman-Mitnik conducted this work while a Sabbatical Fellow at the University of the Balearic Islands from which support is gratefully acknowledged. This work was cofunded by the Ministerio de Economía y Competitividad (MINECO) and the European Fund for Regional Development (FEDER) (CTQ2014-55835-R).

Author Contributions

Daniel Glossman-Mitnik conceived and designed the research and headed, wrote and revised the manuscript, while Sebastián Sastre and Juan Frau contributed to the calculations and to the writing and the revision of the article. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest regarding the publication of this manuscript.

Abbreviations

The following abbreviations are used in this manuscript:
DFTDensity Functional Theory
KSKohn-Sham
EElectronic Energy
IIonization Potential
AElectron Affinity
NNumber of Electrons
AGEAdvanced Glycation Endproducts
MEDTMolecular Electron Density Theory
SCFSelf Consistent Field
KIDKoopmans in DFT
HOMOHigher Ocuppied Molecular Orbital
LUMOLower Unocuppied Molecular Orbital
IEF-PCMIntegral Equation Formalism–Polarized Continuum Model
SMDSolvation Model Density

References

  1. Parr, R.; Yang, W. Density-Functional Theory of Atoms and Molecules; Oxford University Press: New York, NY, USA, 1989. [Google Scholar]
  2. Geerlings, P.; de Proft, F.; Langenaeker, W. Conceptual Density Functional Theory. Chem. Rev. 2003, 103, 1793–1873. [Google Scholar] [CrossRef] [PubMed]
  3. Toro-Labbé, A. Theoretical Aspects of Chemical Reactivity; Elsevier Science: Amsterdam, The Netherlands, 2007; Volume 19. [Google Scholar]
  4. Chattaraj, P. Chemical Reactivity Theory—A Density Functional View; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2009. [Google Scholar]
  5. Domingo, L.R. Molecular Electron Density Theory: A Modern View of Reactivity in Organic Chemistry. Molecules 2016, 21, 1–15. [Google Scholar] [CrossRef] [PubMed]
  6. Domingo, L.R.; Ríos-Gutiérrez, M.; Pérez, P. Applications of the Conceptual Density Functional Theory Indices to Organic Chemistry Reactivity. Molecules 2016, 21, 1–22. [Google Scholar] [CrossRef] [PubMed]
  7. Jasinski, R. Searching for Zwitterionic Intermediates in Hetero Diels-Alder Reactions between Methyl α,p-dinitrocinnamate and Vinyl-Alkyl Ethers. Comput. Theor. Chem. 2014, 1046, 93–98. [Google Scholar] [CrossRef]
  8. Jasiński, R.; Mróz, K.; Kacka, A. Experimental and Theoretical DFT Study on Synthesis of Sterically Crowded 2,3,3,(4)5-Tetrasubstituted-4-nitroisoxazolidines via 1,3-Dipolar Cycloaddition Reactions Between Ketonitrones and Conjugated Nitroalkenes. J. Heterocycl. Chem. 2016, 53, 1424–1429. [Google Scholar] [CrossRef]
  9. Jasiński, R.; Jasińska, E.; Dresler, E. A DFT Computational Study of the Molecular Mechanism of [3+2] Cycloaddition Reactions between Nitroethene and Benzonitrile N-oxides. J. Mol. Mod. 2017, 23, 1–9. [Google Scholar] [CrossRef] [PubMed]
  10. Ahmed, N. Advanced Glycation Endproducts—Role in Pathology of Diabetic Complications. Diabetes Res. Clin. Pract. 2005, 67, 3–21. [Google Scholar] [CrossRef] [PubMed]
  11. Makowska, J.; Bagińska, K.; Liwo, A.; Chmurzynski, L.; Scheraga, H. Acidic-Basic Properties of Three Alanine-Based Peptides Containing Acidic and Basic Side Chains: Comparison between Theory and Experiment. Biopolymers 2008, 90, 724–732. [Google Scholar] [CrossRef] [PubMed]
  12. Makowska, J.; Liwo, A.; Chmurzynski, L.; Scheraga, H.A. Influence of the Length of the Alanine Spacer on the Acidic-Basic Properties of the Ac-Lys-(Ala)n-Lys-NH2 Peptides (n = 0, 1, 2, , 5). J. Solut. Chem. 2012, 41, 1738–1746. [Google Scholar] [CrossRef] [PubMed]
  13. Huzinaga, S.; Andzelm, J.; Klobukowski, M.; Radzio-Audzelm, E.; Sakai, Y.; Tatewaki, H. Gaussian Basis Sets for Molecular Calculations; Elsevier: Amsterdam, The Netherlands, 1984. [Google Scholar]
  14. Easton, R.; Giesen, D.; Welch, A.; Cramer, C.; Truhlar, D. The MIDI! Basis Set for Quantum Mechanical Calculations of Molecular Geometries and Partial Charges. Theor. Chem. Acc. 1996, 93, 281–301. [Google Scholar] [CrossRef]
  15. Lewars, E. Computational Chemistry—Introduction to the Theory and Applications of Molecular and Quantum Mechanics; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2003. [Google Scholar]
  16. Young, D. Computational Chemistry—A Practical Guide for Applying Techniques to Real-World Problems; John Wiley & Sons: New York, NY, USA, 2001. [Google Scholar]
  17. Jensen, F. Introduction to Computational Chemistry, 2nd ed.; John Wiley & Sons: Chichester, UK, 2007. [Google Scholar]
  18. Cramer, C. Essentials of Computational Chemistry—Theories and Models, 2nd ed.; John Wiley & Sons: Chichester, UK, 2004. [Google Scholar]
  19. Peverati, R.; Truhlar, D.G. Quest for a Universal Density Functional: The Accuracy of Density Functionals Across a Broad Spectrum of Databases in Chemistry and Physics. Philos. Trans. R. Soc. A 2014, 372, 1–52. [Google Scholar] [CrossRef] [PubMed]
  20. Glossman-Mitnik, D. A Comparison of the Chemical Reactivity of Naringenin Calculated with M06 Family of Density Functionals. Chem. Cent. J. 2013, 7, 155–161. [Google Scholar] [CrossRef] [PubMed]
  21. Martínez-Araya, J.I.; Salgado-Morán, G.; Glossman-Mitnik, D. Computational Nanochemistry Report on the Oxicams—Conceptual DFT and Chemical Reactivity. J. Phys. Chem. B 2013, 117, 6639–6651. [Google Scholar] [CrossRef] [PubMed]
  22. Glossman-Mitnik, D. Computational Nanochemistry Study of the Chemical Reactivity Properties of the Rhodamine B Molecule. Procedia Comput. Sci. 2013, 18, 816–825. [Google Scholar] [CrossRef]
  23. Martínez-Araya, J.I.; Salgado-Morán, G.; Glossman-Mitnik, D. Computational Nutraceutics: Chemical Reactivity Properties of the Flavonoid Naringin by Means of Conceptual DFT. J. Chem. 2013. [Google Scholar] [CrossRef]
  24. Glossman-Mitnik, D. Chemical Reactivity Theory within DFT Applied to the Study of the Prunin Flavonoid. Eur. Int. J. Sci. Technol. 2014, 3, 195–207. [Google Scholar]
  25. Glossman-Mitnik, D. Computational Chemistry of Natural Products: A Comparison of the Chemical Reactivity of Isonaringin Calculated with the M06 Family of Density Functionals. J. Mol. Mod. 2014, 20, 1–7. [Google Scholar] [CrossRef] [PubMed]
  26. Frau, J.; Muñoz, F.; Glossman-Mitnik, D. A Molecular Electron Density Theory Study of the Chemical Reactivity of Cis- and Trans-Resveratrol. Molecules 2016, 21, 1–13. [Google Scholar] [CrossRef] [PubMed]
  27. Frau, J.; Glossman-Mitnik, D. Chemical Reactivity Theory Study of Advanced Glycation Endproduct Inhibitors. Molecules 2017, 22, 1–9. [Google Scholar] [CrossRef] [PubMed]
  28. Parr, R.; Yang, W. Density Functional Approach to the Frontier-Electron Theory of Chemical Reactivity. J. Am. Chem. Soc. 1984, 106, 4049–4050. [Google Scholar] [CrossRef]
  29. Parr, R.; Szentpaly, L.; Liu, S. Electrophilicity Index. J. Am. Chem. Soc. 1999, 121, 1922–1924. [Google Scholar] [CrossRef]
  30. Morell, C.; Grand, A.; Toro-Labbé, A. New Dual Descriptor for Chemical Reactivity. J. Phys. Chem. A 2005, 109, 205–212. [Google Scholar] [CrossRef] [PubMed]
  31. 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]
  32. Cárdenas, C.; Rabi, N.; Ayers, P.; Morell, C.; Jaramillo, P.; Fuentealba, P. Chemical Reactivity Descriptors for Ambiphilic Reagents: Dual Descriptor, Local Hypersoftness, and Electrostatic Potential. J. Phys. Chem. A 2009, 113, 8660–8667. [Google Scholar] [CrossRef] [PubMed]
  33. Ayers, P.; Morell, C.; de Proft, F.; Geerlings, P. Understanding the Woodward-Hoffmann Rules by Using Changes in Electron Density. Chem. Eur. J. 2007, 13, 8240–8247. [Google Scholar] [CrossRef] [PubMed]
  34. Morell, C.; Ayers, P.; Grand, A.; Gutiérrez-Oliva, S.; Toro-Labbé, A. Rationalization of the Diels-Alder Reactions through the Use of the Dual Reactivity Descriptor Δf(r). Phys. Chem. Chem. Phys. 2008, 10, 7239–7246. [Google Scholar] [CrossRef] [PubMed]
  35. Morell, C.; Hocquet, A.; Grand, A.; Jamart-Grégoire, B. A Conceptual DFT Study of Hydrazino Peptides: Assessment of the Nucleophilicity of the Nitrogen Atoms by Means of the Dual Descriptor Δf(r). J. Mol. Struct. THEOCHEM 2008, 849, 46–51. [Google Scholar] [CrossRef]
  36. 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]
  37. Domingo, L.R.; Pérez, P.; Sáez, J. Understanding the Local Reactivity in Polar Organic Reactions through Electrophilic and Nucleophilic Parr Functions. RSC Adv. 2013, 3, 1486–1494. [Google Scholar] [CrossRef]
  38. Chamorro, E.; Pérez, P.; Domingo, L.R. On the Nature of Parr Functions to Predict the Most Reactive Sites along Organic Polar Reactions. Chem. Phys. Lett. 2013, 582, 141–143. [Google Scholar] [CrossRef]
  39. 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 D.01; Gaussian Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  40. 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] [PubMed]
  41. Weigend, F. Accurate Coulomb-fitting Basis Sets for H to R. Phys. Chem. Chem. Phys. 2006, 8, 1057–1065. [Google Scholar] [CrossRef] [PubMed]
  42. Peverati, R.; Truhlar, D.G. Improving the Accuracy of Hybrid Meta-GGA Density Functionals by Range Separation. J. Phys. Chem. Lett. 2011, 2, 2810–2817. [Google Scholar] [CrossRef]
  43. Peverati, R.; Truhlar, D.G. M11-L: A Local Density Functional That Provides Improved Accuracy for Electronic Structure Calculations in Chemistry and Physics. J. Phys. Chem. Lett. 2012, 3, 117–124. [Google Scholar] [CrossRef]
  44. Peverati, R.; Truhlar, D.G. An Improved and Broadly Accurate Local Approximation to the Exchange-Correlation Density Functional: the MN12-L Functional for Electronic Structure Calculations in Chemistry and Physics. Phys. Chem. Chem. Phys. 2012, 14, 13171–13174. [Google Scholar] [CrossRef] [PubMed]
  45. 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] [PubMed]
  46. Peverati, R.; Truhlar, D.G. Exchange-Correlation Functional with Good Accuracy for Both Structural and Energetic Properties while Depending Only on the Density and Its Gradient. J. Chem. Theory Comput. 2012, 8, 2310–2319. [Google Scholar] [CrossRef] [PubMed]
  47. 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 Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378–6396. [Google Scholar] [CrossRef] [PubMed]
  48. Avogadro: An Open-Source Molecular Builder and Visualization Tool—Version 1.2.0. Available online: http://avogadro.openmolecules.net (accessed on 6 October 2016).
  49. Hanweel, M.; Curtis, D.E.; Lonie, D.; Vandermeersch, T.; Zurek, E.; Hutchison, G. Avogadro: An Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J. Cheminf. 2012, 4, 1–17. [Google Scholar] [CrossRef] [PubMed]
  50. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and Testing of a General AMBER Force Field. J. Comp. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef] [PubMed]
  51. Gorelsky, S. AOMix Program for Molecular Orbital Analysis—Version 6.5; University of Ottawa: Ottawa, ON, Canada, 2011. [Google Scholar]
  52. Gorelsky, S.; Lever, A. Electronic Structure and Spectra of Ruthenium Diimine Complexes by Density Functional Theory and INDO/S - Comparison of the Two Methods. J. Organomet. Chem. 2001, 635, 187–196. [Google Scholar] [CrossRef]
  • Sample Availability: Not available.
Figure 1. Geometrical conformation for the KAK peptide(Ac-Lys-(Ala) n -Lys-NH2 with n = 1).
Figure 1. Geometrical conformation for the KAK peptide(Ac-Lys-(Ala) n -Lys-NH2 with n = 1).
Molecules 22 00458 g001
Table 1. Average descriptors J I , J A , J H L , J χ , J η , J ω , J D 1 for the KK, KAK, KA2K, KA3K, KA4K and KAK5 peptides calculated with the M11, M11L, MN12L, MN12SX, N12, and N12SX density functionals and the Def2TZVP basis set using water as as solvent simulated with the solvation model density (SMD) parametrization of the Integral Equation Formalism–Polarized Continuum Model (IEF-PCM) model.
Table 1. Average descriptors J I , J A , J H L , J χ , J η , J ω , J D 1 for the KK, KAK, KA2K, KA3K, KA4K and KAK5 peptides calculated with the M11, M11L, MN12L, MN12SX, N12, and N12SX density functionals and the Def2TZVP basis set using water as as solvent simulated with the solvation model density (SMD) parametrization of the Integral Equation Formalism–Polarized Continuum Model (IEF-PCM) model.
J I J A J HL J χ J η J ω J D 1
M112.792.383.670.205.370.375.19
M11L1.660.071.670.801.730.171.91
MN12L1.870.141.880.872.010.82.20
MN12SX0.130.080.150.090.060.050.14
N121.400.221.440.591.620.101.74
N12SX0.000.070.070.030.070.030.08
Table 2. Condensed dual descriptors and electrophilic Parr functions P k - for the KK, KAK, KA2K, KA3K, KA4K, and KA5K peptides calculated with the MN12SX and N12SX density functionals and the Def2TZVP basis set using water as as solvent simulated with the SMD parametrization of the IEF-PCM model. MPA: Mulliken Population Analysis; HPA: Hirshfeld Population Analysis.
Table 2. Condensed dual descriptors and electrophilic Parr functions P k - for the KK, KAK, KA2K, KA3K, KA4K, and KA5K peptides calculated with the MN12SX and N12SX density functionals and the Def2TZVP basis set using water as as solvent simulated with the SMD parametrization of the IEF-PCM model. MPA: Mulliken Population Analysis; HPA: Hirshfeld Population Analysis.
MN12SXN12SX
f k P k - (mpa)P k - (hpa) f k P k - (mpa)P k - (hpa)
KK−0.630.760.74−0.690.760.73
KAK−0.620.760.73−0.690.760.73
KA2K−0.710.810.78−0.740.780.73
KA3K−0.630.760.74−0.690.760.73
KA4K−0.710.810.78−0.740.810.78
KA5K−0.710.810.78−0.740.810.78

Share and Cite

MDPI and ACS Style

Sastre, S.; Frau, J.; Glossman-Mitnik, D. Computational Prediction of the Protonation Sites of Ac-Lys-(Ala)n-Lys-NH2 Peptides through Conceptual DFT Descriptors. Molecules 2017, 22, 458. https://doi.org/10.3390/molecules22030458

AMA Style

Sastre S, Frau J, Glossman-Mitnik D. Computational Prediction of the Protonation Sites of Ac-Lys-(Ala)n-Lys-NH2 Peptides through Conceptual DFT Descriptors. Molecules. 2017; 22(3):458. https://doi.org/10.3390/molecules22030458

Chicago/Turabian Style

Sastre, Sebastián, Juan Frau, and Daniel Glossman-Mitnik. 2017. "Computational Prediction of the Protonation Sites of Ac-Lys-(Ala)n-Lys-NH2 Peptides through Conceptual DFT Descriptors" Molecules 22, no. 3: 458. https://doi.org/10.3390/molecules22030458

APA Style

Sastre, S., Frau, J., & Glossman-Mitnik, D. (2017). Computational Prediction of the Protonation Sites of Ac-Lys-(Ala)n-Lys-NH2 Peptides through Conceptual DFT Descriptors. Molecules, 22(3), 458. https://doi.org/10.3390/molecules22030458

Article Metrics

Back to TopTop