Computational prediction of the pKas of small peptides through Conceptual DFT descriptors

Quantitative In Silico Prediction of the Rate of Protodeboronation by a Mechanistic Density Functional Theory-Aided Algorithm

Computational reaction prediction has become a ubiquitous task in chemistry due to the potential value accurate predictions can bring to chemists. Boronic acids are widely used in industry; however, understanding how to avoid the protodeboronation side reaction remains a challenge. We have developed an algorithm for in silico prediction of the rate of protodeboronation of boronic acids. A general mechanistic model devised through kinetic studies of protodeboronation was found in the literature and forms the foundation on which the algorithm presented in this work is built. Protodeboronation proceeds through 7 distinct pathways, though for any particular boronic acid, only a subset of mechanistic pathways are active. The rate of each active mechanistic pathway is linearly correlated with its characteristic energy difference, which in turn can be determined using Density Functional Theory. We validated the algorithm using leave-one-out cross-validation on a data set of 50 boronic acids and made a further 50 rate predictions on academically and industrially important boronic acids out of sample. We believe this work will provide great assistance to chemists performing reactions that feature boronic acids, such as Suzuki-Miyaura and Chan-Evans-Lam couplings.

An insight into reactivity and bioactivity properties of quorum sensing peptides against PDE10A: a computational peptidology approach

Peptides are currently the most promising lead molecules. Quorum sensing peptides have a variety of structural features and are regularly exposed to post-translational modifications. Antiparkinsonian drugs lose their efficacy after a long period of use, and patients develop motor problems such as drug-induced dyskinesia (DIDs). The interaction between PDE10A and cAMP is necessary for dopamine neurotransmission and may play a role in Parkinson's disease pathogenesis. cAMP and cGMP are cyclic nucleotides that act as secondary messengers in the signal transduction pathway, influencing a range of CNS activities. PDE enzymes hydrolyze phosphodiester bonds to break down cAMP and cGMP, allowing them to control intracellular levels of these second messengers effectively. PDE expression, and hence cyclic nucleotide levels and their downstream targets, may change with age and in numerous age-related illnesses, including Parkinson's disease, according to mounting evidence. At the peak of dyskinesias, cyclic nucleotide levels were lower, and using phosphodiesterase inhibitors before antiparkinsonian medicines reduced the severity of dyskinesias. In a recent study, PapRIV was found to have the ability to activate BV-2 microglia cells, indicating that this quorum sensing peptide may play a role in gut-brain contact. As a result of the current in silico work, mainly focused on QSPs as a lead molecule for inhibiting PDE10A, the SRNAT QSP sequence has been a potent molecule in molecular docking and molecular dynamics simulations. Furthermore, we can test the efficiency of therapeutic components in vitro and in vivo utilizing this computational approach against PDE10A.

Virtual Prospection of Marine Cyclopeptides as Therapeutics by Means of Conceptual DFT and Computational ADMET

Bioactive peptides are chemical compounds created through the covalent bonding of amino acids, known as amide or peptide bonds. Due to their unusual chemistry and various biological effects, marine bioactive peptides have garnered considerable research. The effectiveness of a bioactive marine peptide is attributed to its structural features, such as amino acid content and sequence, which vary depending on the degree of action. Cyclic peptides combine several favorable properties such as good binding affinity, target selectivity and low toxicity that render them an attractive modality for the development of therapeutics. The apratoxins are a class of molecules formed by a series of cyclic depsipeptides with potent cytotoxic activities. The objective of this research is to pursue a computational prospection of the molecular structures and properties of several cylopeptides of marine origin with potential therapeutic applications. The methodology will be based on the determination of the chemical reactivity descriptors of the studied molecules through the consideration of the Conceptual DFT model and validation of a particular model chemistry, MN12SX/Def2TZVP/H2O. These studies will be complemented by a determination of the pharmacokinetics and ADMET parameters by resorting to certain cheminformatics tools.

Conceptual DFT-Based Computational Peptidology, Pharmacokinetics Study and ADMET Report of the Veraguamides A–G Family of Marine Natural Drugs

As a continuation of our research on the chemical reactivity, pharmacokinetics and ADMET properties of cyclopeptides of marine origin with potential therapeutic abilities, in this work our already presented integrated molecular modeling protocol has been used for the study of the chemical reactivity and bioactivity properties of the Veraguamides A–G family of marine natural drugs. This protocol results from the estimation of the conceptual density functional theory (CDFT) chemical reactivity descriptors together with several chemoinformatics tools commonly considered within the process of development of new therapeutic drugs. CP-CDFT is a branch of computational chemistry and molecular modeling dedicated to the study of peptides, and it is a protocol that allows the estimation with great accuracy of the CDFT-based reactivity descriptors and the associated physical and chemical properties, which can aid in determining the ability of the studied peptides to behave as potential useful drugs. Moreover, the superiority of the MN12SX density functional over other long-range corrected density functionals for the prediction of chemical and physical properties in the presence of water as the solvent is clearly demonstrated. The research was supplemented with an investigation of the bioactivity of the molecular systems and their ADMET (absorption, distribution, metabolism, excretion, and toxicity) parameters, as is customary in medicinal chemistry. Some instances of the CDFT-based chemical reactivity descriptors’ capacity to predict the pKas of peptides as well as their potential as AGE inhibitors are also shown.

An integrated molecular modeling protocol for drug screening based on conceptual density functional theory and chemoinformatics for the study of marine cyclopeptides

An integrated molecular modeling protocol resulting from the combination of conceptual density functional theory (CDFT) chemical reactivity descriptors with several chemoinformatics tools has been used for the study of the chemical reactivity and bioactivity properties of a group of marine cyclic peptides. CP-CDFT is a branch of computational chemistry and molecular modeling dedicated to the study of peptides. The protocol allowed the estimation of the CDFT-based reactivity indices together with the associated physicochemical parameters that can help to identify the ability of the studied peptides to behave as potential useful drugs. This was complemented with an analysis of the bioactivity and pharmacokinetics parameters related to the ADMET (absorption, distribution, metabolism, excretion, and toxicity) features. Some examples related to the ability of the CDFT-based chemical reactivity descriptors for the prediction of the pKas of the peptides as well as their potential as AGE inhibitors are also presented.

Effect of the aqueous matrix on the inactivation of E. coli by permaleic acid

In this work permaleic acid (PMA) was investigated as possible disinfecting agent and compared to peracetic acid (PAA) in real tap water and wastewater. Preliminary tests in lysogeny broth (LB) were also performed. PMA was synthesized from maleic anhydride and hydrogen peroxide and, for the first time, its antimicrobial activity was evaluated with respect to the growth inhibition of E. coli. The effect of the pH and bivalent ions, typically occurring in real water matrices (namely, Mg²⁺ and Ca²⁺), was also investigated. pKa values for PMA were calculated for the first time by DFT calculations. The concentration of bivalent ions strongly affected disinfection efficiency with PMA (Ca²⁺=0.33 mgL^-1 and Mg²⁺=0.35 mgL^-1: 100% E. coli reduction > log 5; Ca²⁺=13.3 mg L^-1 and Mg²⁺=25.6 mg L^-1: E. coli reduction < log 1, after 60 min), and such results were supported by DFT modelling outcomes (pKa₂ of PMA 7.3) and disinfection tests in presence of EDTA chelating agent. More alkaline pH conditions drastically decreased PMA disinfection (pH = 5: > log 5 E.coli reduction; pH = 9: < log 1 E.coli reduction, after 60 min). PMA disinfection efficiency is strongly affected by the target water quality, the concentration of metal bivalent ions and the initial pH.

Virtual Screening of Marine Natural Compounds by Means of Chemoinformatics and CDFT-Based Computational Peptidology

This work presents the results of a computational study of the chemical reactivity and bioactivity properties of the members of the theopapuamides A-D family of marine peptides by making use of our proposed methodology named Computational Peptidology (CP) that has been successfully considered in previous studies of this kind of molecular system. CP allows for the determination of the global and local descriptors that come from Conceptual Density Functional Theory (CDFT) that can give an idea about the chemical reactivity properties of the marine natural products under study, which are expected to be related to their bioactivity. At the same time, the validity of the procedure based on the adoption of the KID (Koopmans In DFT) technique, as well as the MN12SX/Def2TZVP/H2O model chemistry is successfully verified. Together with several chemoinformatic tools that can be used to improve the process of virtual screening, some additional properties of these marine peptides are identified related to their ability to behave as useful drugs. With the further objective of analyzing their bioactivity, some useful parameters for future QSAR studies, their predicted biological targets, and the ADMET (Absorption, Distribution, Metabolism, Excretion and Toxicity) parameters related to the theopapuamides A-D pharmacokinetics are also reported.

Conceptual DFT-Based Computational Peptidology of Marine Natural Compounds: Discodermins A-H

A methodology based on the concepts that arise from Density Functional Theory named Conceptual Density Functional Theory (CDFT) was chosen for the calculation of some global and local reactivity descriptors of the Discodermins A-H family of marine peptides through the consideration of the KID (Koopmans in DFT) technique that was successfully used in previous studies of this kind of molecular systems. The determination of active sites of the studied molecules for different kinds of reactivities was achieved by resorting to some CDFT-based descriptors like the Fukui functions as well as the Parr functions derived from Molecular Electron Density Theory (MEDT). A few properties identified with their ability to behave as a drug and the bioactivity of the peptides considered in this examination were acquired by depending on a homology model by studying the correlation with the known bioactivity of related molecules in their interaction with various biological receptors. With the further object of analyzing their bioactivity, some parameters of usefulness for future QSAR studies, their predicted biological targets, and the ADME (Absorption, Distribution, Metabolism, and Excretion) parameters related to the Discodermins A-H pharmacokinetics are also reported.

Chemical reactivity and bioactivity properties of pyrazinamide analogs of acetylsalicylic acid and salicylic acid using conceptual density functional theory

Conventional drugs used to treat Tuberculosis (TB) are becoming ineffective due to the occurrence of multiple drug resistant strains of tuberculosis (TB). This has made the TB disease a a serious global health dilemma. Hence, there is desperate necessity for the advancement of new drugs. In this work, the chemical reactivity and bioactivity of several analogs ofpyrazinamide (PZA) were investigated. PZA is one of the first-line of drugs used to treat tuberculosis and is a key contributor to shortening the treatment time for the disease. Chemical reactivity descriptors of pyrazinamide (PZA) and its analogs of acetylsalicyclic acid and salicyclic acid were investigated using conceptual density functional theory in water as a solvent at the MN12SX/Def2TZVP level of theory. Results have shown that all PZA analogs have improved their global and local reactivity indeces as compared to pyrazinamide based on its electronegativity, electrodonating power, electroaccepting power, eletrophilicity, global hardness and dual descriptor condensed fukui indexes. Moreover, their pKa values are slightly higher than PZA. In terms of its drug-likeness, all PZA analogs passed the Lipinski's Rule of Five criteria. Furthermore, their bioactivity scores are significantly better than pyrazinamide indicating good reaction to G-Protein Coupled Receptor (GPCR) ligands, kinase inhibitors, ion channel modulators, nuclear receptor ligands, protease inhibitors and other enzyme targets. Overall, the PZA analogs are found to be promising anti-tuberculosis drugs. Based on global and local reactivity descriptors, pKa and bioactivity scores, PZA analog of 5-n-Octanoylsalicylic acid is the most reactive among the PZA analogs tested.

Conceptual DFT as a Novel Chemoinformatics Tool for Studying the Chemical Reactivity Properties of the Amatoxin Family of Fungal Peptides

The chemical structures and molecular reactivities of the Amatoxin group of fungi-derived peptides have been determined by means of the consideration of a model chemistry that has been previously validated as well-behaved for our purposes. The reactivity descriptors were calculated on the basis of a methodological framework built around the concepts that are the outcome of the so called Conceptual Density Functional Theory (CDFT). This procedure in connection with the different Fukui functions allowed to identify the chemically active regions within the molecules. By considering a simple protocol designed by our research group for the estimation of the pKa of peptides through the information coming from the chemical hardness, these property has been established for the different molecular systems explored in this research. The information reported through this work could be of interest for medicinal chemistry researchers in using this knowledge for the design of new medicines based on the studied peptides or as a help for the understanding of the toxicity mechanisms exerted by them.

Chemical-Reactivity Properties, Drug Likeness, and Bioactivity Scores of Seragamides A–F Anticancer Marine Peptides: Conceptual Density Functional Theory Viewpoint

A methodology based on concepts that arose from Density Functional Theory (CDFT) was chosen for the calculation of global and local reactivity descriptors of the Seragamide family of marine anticancer peptides. Determination of active sites for the molecules was achieved by resorting to some descriptors within Molecular Electron Density Theory (MEDT) such as Fukui functions. The pKas of the six studied peptides were established using a proposed relationship between this property and calculated chemical hardness. The drug likenesses and bioactivity properties of the peptides considered in this study were obtained by resorting to a homology model by comparison with the bioactivity of related molecules in their interaction with different receptors. With the object of analyzing the concept of drug repurposing, a study of potential AGE-inhibition abilities of Seragamides peptides was pursued by comparison with well-known drugs that are already available as pharmaceuticals.

Calculation of the Global and Local Conceptual DFT Indices for the Prediction of the Chemical Reactivity Properties of Papuamides A-F Marine Drugs

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 structures of the Papuamide family of marine peptides. A methodology based on Conceptual Density Functional Theory (CDFT) was chosen for the determination of the reactivity descriptors. The molecular active sites were associated with the active regions of the molecules related to the nucleophilic and electrophilic Parr functions. Finally, the drug-likenesses and the bioactivity scores for the Papuamide peptides were predicted through a homology methodology relating them with the calculated reactivity descriptors, while other properties such as the pKas were determined following a methodology developed by our group.

Chemical reactivity and bioactivity properties of the Phallotoxin family of fungal peptides based on Conceptual Peptidology and DFT study

A methodology based on the concepts that arise from Density Functional Theory (CDFT) was chosen for the calculation of the global and local reactivity descriptors of the Phallotoxin family of fungal peptides. The determination of the active sites for the molecules has been achieved by resorting some descriptors within Molecular Electron Density Theory (MEDT) like the Dual Descriptor and the Parr functions. Phallosacin has been found as the most reactive of the peptides on the basis of the calculated Global Reactivity Descriptors. The pKas of the seven studied peptides were established using a proposed relationship between this property and the calculated Global Hardness. The bioactivity properties of the peptides considered in this study were obtained by resorting to a homology model by comparison with the bioactivity of related molecules in their interaction with different receptors.

Computational Peptidology Assisted by Conceptual Density Functional Theory for the Study of Five New Antifungal Tripeptides

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 structures of a group of five new antifungal tripeptides, namely (2R)-2-[(2S)-2-[(2S)-2-amino-3-phenylpropanamido]propanamido]-5-[(diaminomethylidene)amino]pentanoic acid, (2S)-2-[(2S)-2-[(2S)-2-amino-3-phenyl propanamido]propanamido]-3-(4-hydroxyphenyl)propanoic acid, (2S)-2-[(2S)-2-[(2S)-2-amino-3-phenylpropanamido]-3-methylbutanamido]-3-(4-hydroxyphenyl)propanoic acid, (2R)-2-[(2S)-2-[(2S)-2-amino-3-phenylpropanamido]-3-(1H-indol-3-yl)propanamido]-3-sulfanylpropanoic acid, and (2S)-2-[(2S)-2-[(2S)-2-amino-3-phenylpropanamido]-3-(1H-indol-3-yl)propanamido]-3-(4-hydroxyphenyl)propanoic acid, according to their amino acid sequences. A methodology based on conceptual density functional theory was chosen for the determination of the reactivity descriptors. The molecular active sites were associated with the active regions of the molecules that were associated with the nucleophilic, electrophilic, and radical Fukui functions. Additionally, the pK a values for the different peptides are predicted with great accuracy, which constitutes a useful knowledge for the process of drug design. Finally, the bioactivity scores for the new antifungal peptides are predicted through a homology methodology relating them with the calculated reactivity descriptors.

CDFT-Based Reactivity Descriptors as a Useful MEDT Chemoinformatics Tool for the Study of the Virotoxin Family of Fungal Peptides

Virotoxins are monocyclic peptides formed by at least five different compounds: alaviroidin, viroisin, deoxoviroisin, viroidin and deoxovirodin. These are toxic peptides singularly found in Amanita virosa mushrooms. Here we perform computational studies on the structural and electronic conformations of these peptides using the MN12SX/Def2TZVP/H2O chemistry model to investigate their chemical reactivity. CDFT-based descriptors (for Conceptual Density Functional Theory) (e.g., Parr functions and Nucleophilicity) are also considered. At the same time, other properties (e.g., pKas) will be determined and used to study virotoxins solubility and to inform decisions about repurposing these agents in medicinal chemistry.

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

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.

A computational scheme of pKa values based on the three-dimensional reference interaction site model self-consistent field theory coupled with the linear fitting correction scheme

A scheme for quantitatively computing the acid dissociation constant, pKa, of hydrated molecules is proposed. It is based on the three-dimensional reference interaction site model self-consistent field (3D-RISM-SCF) theory coupled with the linear fitting correction (LFC) scheme. In LFC/3D-RISM-SCF, pKa values of target molecules are evaluated using the Gibbs energy difference between the protonated and unprotonated states calculated by 3D-RISM-SCF and the parameters fitted by the LFC scheme to the experimental values of training set systems. The pKa values computed by LFC/3D-RISM-SCF show quantitative agreement with the experimental data.

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

The MN12SX density functional, in connection with the Def2TZVP basis set, was assessed, together with the SMD solvation model (Solvation Model based on the Density), for calculation of the molecular properties and structure of a group of peptides of marine origin named Mirabamides A⁻H. All the chemical reactivity descriptors for the systems were calculated via Conceptual Density Functional Theory (CDFT). The active sites suitable for nucleophilic, electrophilic, and radical attacks were chosen by linking them with the Fukui function indices, nucleophilic and electrophilic Parr functions, and condensed Dual Descriptor Δ f ( r ) , respectively. Additionally, the p K a values for the different peptides are predicted with great accuracy as well as the ability of the studied molecule in acting as an efficient inhibitor of the formation of Advanced Glycation Endproducts (AGEs), which constitutes a useful knowledge for the development of drugs for fighting Diabetes, Alzheimer and Parkinson diseases. Finally, the bioactivity scores for the Mirabamides A⁻H are predicted through different methodologies.

Blue M2: an intermediate melanoidin studied via conceptual DFT

In this computational study, ten density functionals, viz. CAM-B3LYP, LC-ω PBE, M11, M11L, MN12L, MN12SX, N12, N12SX, ω B97X, and ω B97XD, related to the Def2TZVP basis sets, are assessed together with the SMD solvation model for calculation of the molecular properties and structure of blue-M2 intermediate melanoidin pigment. All the chemical reactivity descriptors for the system are calculated via conceptual density functional theory (DFT). The active sites suitable for nucleophilic, electrophilic, and radical attacks are selected by linking them with the Fukui function indices, electrophilic Parr functions, and condensed dual descriptors Δf(r), respectively. The prediction of the maximum absorption wavelength is considerably accurate relative to its experimental value. The study reveals that the MN12SX and N12SX density functionals are the most appropriate density functionals for predicting the chemical reactivity of the molecule under study.

Conceptual DFT Study of the Local Chemical Reactivity of the Colored BISARG Melanoidin and Its Protonated Derivative

This computational study assessed eight fixed RSH (range-separated hybrid) density functionals that include CAM-B3LYP, LC-ωPBE, M11, MN12SX, N12SX, ωB97, ωB97X, and ωB97XD related to the Def2TZVP basis sets together with the SMD solvation model in the calculation the molecular structure and reactivity properties of the BISARG intermediate melanoidin pigment (5-(2-(E)-(Z)-5-[(2-furyl)methylidene]-3-(4-acetylamino-4-carboxybutyl)-2-imino-1,3-dihydroimidazol-4-ylideneamino(E)-4-[(2-furyl)methylidene]-5-oxo-1H-imidazol-1-yl)-2-acetylaminovaleric acid) and its protonated derivative, BISARG(p). The chemical reactivity descriptors for the systems were calculated via the Conceptual Density Functional Theory. The choice of active sites applicable to nucleophilic, electrophilic as well as radical attacks were made by linking them with Fukui functions indices, electrophilic and nucleophilic Parr functions, and the condensed Dual Descriptor Δf(r). The study found the MN12SX and N12SX density functionals to be the most appropriate in predicting the chemical reactivity of the molecular systems under study starting from the knowledge of the HOMO, LUMO, and HOMO-LUMO gap energies.

Molecular Reactivity and Absorption Properties of Melanoidin Blue-G1 through Conceptual DFT

This computational study presents the assessment of eleven density functionals that include CAM-B3LYP, LC-wPBE, M11, M11L, MN12L, MN12SX, N12, N12SX, wB97, wB97X and wB97XD related to the Def2TZVP basis sets together with the Solvation Model Density (SMD) solvation model in calculating the molecular properties and structure of the Blue-G1 intermediate melanoidin pigment. The chemical reactivity descriptors for the system are calculated via the conceptual Density Functional Theory (DFT). The choice of the active sites related to the nucleophilic, electrophilic, as well as radical attacks is made by linking them with the Fukui function indices, the electrophilic Parr functions and the condensed dual descriptor Δ f ( r ) . The prediction of the maximum absorption wavelength tends to be considerably accurate relative to its experimental value. The study found the MN12SX and N12SX density functionals to be the most appropriate density functionals in predicting the chemical reactivity of the studied molecule.