Assignment of the 1H and 13C NMR signals of some hydroxyphenylcoumarins.
The complete 1H and 13C NMR assignments of six hydroxyphenylcoumarins have been made using 1D and 2D NMR techniques, including COSY, HMQC and HMBC experiments.
Carbon Isotopes↗
Biomedical subjects
Eugenio Uriarte
Publications and source records attributed to Eugenio Uriarte.
ANN-QSAR model for selection of anticancer leads from structurally heterogeneous series of compounds.
Developing a model for predicting anticancer activity of any classes of organic compounds based on molecular structure is very important goal for medicinal chemist. Different molecular descriptors can be used to solve this problem. Stochastic molecular descriptors so-called the MARCH-INSIDE approach, shown to be very successful in drug design. Nevertheless, the structural diversity of compounds is so vast that we may need non-linear models such as artificial neural networks (ANN) instead of linear ones. SmartMLP-ANN analysis used to model the anticancer activity of organic compounds has shown high average accuracy of 93.79% (train performance) and predictability of 90.88% (validation performance) for the 8:3-MLP topology with different training and predicting series. This ANN model favourably compares with respect to a previous linear discriminant analysis (LDA) model [H. González-Díaz et al., J. Mol. Model 9 (2003) 395] that showed only 80.49% of accuracy and 79.34% of predictability. The present SmartMLP approach employed shorter training times of only 10h while previous models give accuracies of 70-89% only after 25-46 h of training. In order to illustrate the practical use of the model in bioorganic medicinal chemistry, we report the in silico prediction, and in vitro evaluation of six new synthetic tegafur analogues having IC(50) values in a broad range between 37.1 and 138 microgmL(-1) for leukemia (L1210/0) and human T-lymphocyte (Molt4/C8, CEM/0) cells. Theoretical predictions coincide very well with experimental results.
Animals↗
Unified QSAR approach to antimicrobials. Part 2: predicting activity against more than 90 different species in order to halt antibacterial resistance.
There are many different kinds of pathogenic bacteria species with very different susceptibility profiles to different antibacterial drugs. One limitation of QSAR models is that they consider the biological activity of drugs against only one species of bacteria. In a previous paper, we developed a unified Markov model to describe the biological activity of different drugs tested in the literature against some antimicrobial species. Consequently, predicting the probability with which a drug is active against different species of bacteria with a single unified model is a goal of major importance. The work described here develops the unified Markov model to describe the biological activity of more than 70 drugs from the literature tested against 96 species of bacteria. We applied linear discriminant analysis (LDA) to classify drugs as active or inactive against the different tested bacterial species. The model correctly classified 199 out of 237 active compounds (83.9%) and 168 out of 200 inactive compounds (84%). Overall training predictability was 84% (367 out of 437 cases). Validation of the model was carried out using an external predicting series, with the model classifying 202 out of 243 (i.e., 83.13%) of the compounds. In order to show how the model functions in practice, a virtual screening was carried out and the model recognized as active 84.5% (480 out of 568) antibacterial compounds not used in the training or predicting series. The current study is an attempt to calculate within a unified framework the probabilities of antibacterial action of drugs against many different species.
Algorithms↗
QSAR study of anticoccidial activity for diverse chemical compounds: prediction and experimental assay of trans-2-(2-nitrovinyl)furan.
In this work we report a QSAR model that discriminates between chemically heterogeneous classes of anticoccidial and non-anticoccidial compounds. For this purpose we used the Markovian Chemicals in silico Design (MARCH-INSIDE) approach J. Mol. Mod.2002, 8, 237-245; J. Mol. Mod.2003, 9, 395-407]. Linear discriminant analysis allowed us to fit the discriminant function. This function correctly classifies 86.67% of anticoccidial compounds and 96.23% of inactive compounds in the training series. Overall classification is 94.12%. We validated the model by means of an external predicting series, with 86.96% of global predictability. Remarkably, the present model is based on topological as well as configuration-dependent molecular descriptors. Therefore, the model performs timely calculations and allows discrimination between Z/E and chiral isomers. Finally, to exemplify the use of the model in practice we report the prediction and experimental assay of trans-2-(2-nitrovinyl)furan. It is notable that lesion control was 72.86% at mg/kg of body weight with respect to 60% at 125 mg/kg for amprolium (control drug). The back-projection map for this compound predicts a high level of importance for the double bond and for the nitro group in the trans position. We conclude that the MARCH-INSIDE approach enables the accurate fast track identification of anticoccidial hits. Moreover, trans-2-(2-nitrovinyl)furan seems to be a promising drug for the treatment of coccidiosis.
Animals↗
Simple stochastic fingerprints towards mathematical modeling in biology and medicine 2. Unifying Markov model for drugs side effects.
Most of present mathematical models for biological activity consider just the molecular structure. In the present article we pretend extending the use of Markov chain models to define novel molecular descriptors, which consider in addition other parameters like target site or biological effect. Specifically, this mathematical model takes into consideration not only the molecular structure but the specific biological system the drug affects too. Herein, a general Markov model is developed that describes 19 different drugs side effects grouped in eight affected biological systems for 178 drugs, being 270 cases finally. The data was processed by linear discriminant analysis (LDA) classifying drugs according to their specific side effects, forward stepwise was fixed as strategy for variables selection. The average percentage of good classification and number of compounds used in the training/predicting sets were 100/95.8% for endocrine manifestations, (18 out of 18)/(13 out of 14); 90.5/92.3% for gastrointestinal manifestations, (38 out of 42)/(30 out of 32); 88.5/86.5% for systemic phenomena, (23 out of 26)/(17 out of 20); 81.8/77.3% for neurological manifestations, (27 out of 33)/(19 out of 25); 81.6/86.2% for dermal manifestations, (31 out of 38)/(25 out of 29); 78.4/85.1% for cardiovascular manifestation, (29 out of 37)/(24 out of 28); 77.1/75.7% for breathing manifestations, (27 out of 35)/(20 out of 26) and 75.6/75% for psychiatric manifestations, (31 out of 41)/(23 out of 31). Additionally a back-projection analysis (BPA) was carried out for two ulcerogenic drugs to prove in structural terms the physical interpretation of the models obtained. This article develops a mathematical model that encompasses a large number of drugs side effects grouped in specifics biological systems using stochastic absolute probabilities of interaction ((A)pi(k)(j)) by the first time.
Drug-Related Side Effects and Adverse Reactions↗
New furan side tetracyclic allopsoralen derivatives: synthesis and photobiological evaluation.
Novel tetracyclic allopsoralen derivatives characterized by the condensation of a fourth cyclohexenylic (5-7) or benzenic (8-10) ring at the furan side and a methoxy (5 and 8), a hydroxy (6 and 9), or a dimethylaminopropoxy (7 and 10) side chain in the 10 position of the chromophore were prepared. Compounds 7 and 10 showed a strong photoantiproliferative activity, up to 3 orders of magnitude higher than that of the photochemotherapeutic drug 8-methoxypsoralen (8-MOP). The investigation into the mechanism of action demonstrated for 10 the capacity to intercalate between DNA base pairs in the ground state, to give rise to a covalent photoaddition upon UVA irradiation, and to inhibit polymerase chain reaction (PCR) in a sequence-specific manner. Conversely, compound 7 showed a limited capacity to form an intercalative complex and the lack of ability to photoadd to the macromolecule, thus revealing a novel and unusual behavior for an allopsoralen derivative.
Animals↗
Unify QSAR approach to antimicrobials. Part 1: predicting antifungal activity against different species.
Most of up-to-date reported molecular descriptors encode only information about the molecular structure. In previous papers, we have extended stochastic descriptors to encode additional information such as target site, partition system, or biological species [Bioorg. Med. Chem. Lett.2005, 15, 551; Bioorg. Med. Chem. 2005, 13, 1119]. This work develops an unify Markov model to describe with a single linear equation the biological activity of 74 drugs tested in the literature against some of the fungi species selected from a list of 87 species (491 cases in total). The data were processed by linear discriminant analysis (LDA) classifying drugs as active or non-active against the different tested fungi species. The model correctly classifies 338 out of 368 active compounds (91.85%) and 89 out of 123 non-active compounds (72.36%). Overall training predictability was 86.97% (427 out of 491 compounds). Validation of the model was carried out by means of leave-species-out (LSO) procedure. After elimination step-by-step of all drugs tested against one specific species, we record the percentage of good classification of leave-out compounds (LSO-predictability). In addition, robustness of the model to the elimination of the compounds (LSO-robustness) was considered. This aspect was considered as the variation of the percentage of good classification of the modified model (Delta) in LSO with respect to the original one. Average LSO-predictability was 86.41+/-0.95% (average+/-SD) and Delta = -0.55%, being 6 the average number of drugs tested against each fungi species. Results for some of the 87 studied species were Candida albicans: 43 tested compounds, 100% of LSO-predictability, Delta = -3.49%; Candida parapsilosis 23, 100%, Delta = -0.86%; Aspergillus fumigatus 21, 95.20%, Delta = 0.05%; Microsporum canis 12, 91.60%, Delta = -2.84%; Trichophyton mentagrophytes 11, 100%, Delta = -0.51%; Cryptococcus neoformans 10, 90%, Delta = -0.90%. The present one is the first reported unify model that allows one predicting antifungal activity of any organic compound against a very large diversity of fungi pathogens.
Antifungal Agents↗
Probabilistic neural network model for the in silico evaluation of anti-HIV activity and mechanism of action.
A theoretical model has been developed that discriminates between active and nonactive drugs against HIV-1 with four different mechanisms of action for the active drugs. The model was built up using a probabilistic neural network (PNN) algorithm and a database of 2720 compounds. The model showed an overall accuracy of 97.34% in the training series, 85.12% in the selection series, and 84.78% in an external prediction series. The model not only correctly classified a very heterogeneous series of organic compounds but also discriminated between very similar active/nonactive chemicals that belong to the same family of compounds. More specifically, the model recognized 96.02% of nonactive compounds, 94.24% of active compounds that inhibited reverse transcriptase, 97.24% of protease inhibitors, 97.14% of virus uncoating inhibitors, and 90.32% of integrase inhibitors. The results indicate that this approach may represent a powerful tool for modeling large databases in QSAR with applications in medicinal chemistry.
Algorithms↗
A QSAR model for in silico screening of MAO-A inhibitors. Prediction, synthesis, and biological assay of novel coumarins.
This work explores the potential of the MARCH-INSIDE methodology to seek a QSAR for MAO-A inhibitors from a heterogeneous series of compounds. A Markov model was used to quickly calculate the molecular electron delocalization, polarizability, refractivity, and n-octanol/water partition coefficients for a series of 1406 active/nonactive compounds. LDA was subsequently used to fit a classification function. The model showed 92.8% and 91.8% global accuracy and predictability in training and validation studies. This QSAR model was validated through a virtual screening of a series of coumarin derivatives. The 15 selected compounds were prepared and evaluated as in vitro MAO-A inhibitors. The theoretical prediction was compared with the experimental results and the model correctly predicted 13 compounds with only two mistakes on compounds with activities very close to the cutoff point established for the model. Consequently, this method represents a useful tool for the "in silico" screening of MAO-A inhibitors.
Algorithms↗
Novel 2D maps and coupling numbers for protein sequences. The first QSAR study of polygalacturonases; isolation and prediction of a novel sequence from Psidium guajava L.
The development of 2D graph-theoretic representations for DNA sequences was very important for qualitative and quantitative comparison of sequences. Calculation of numeric features for these representations is useful for DNA-QSAR studies. Most of all graph-theoretic representations identify each one of the four bases with a unitary walk in one axe direction in the 2D space. In the case of proteins, twenty amino acids instead of four bases have to be considered. This fact has limited the introduction of useful 2D Cartesian representations and the corresponding sequences descriptors to encode protein sequence information. In this study, we overcome this problem grouping amino acids into four groups: acid, basic, polar and non-polar amino acids. The identification of each group with one of the four axis directions determines a novel 2D representation and numeric descriptors for proteins sequences. Afterwards, a Markov model has been used to calculate new numeric descriptors of the protein sequence. These descriptors are called herein the sequence 2D coupling numbers (zeta(k)). In this work, we calculated the zeta(k) values for 108 sequences of different polygalacturonases (PGs) and for 100 sequences of other proteins. A Linear Discriminant Analysis model derived here (PG=5.36.zeta1-3.98.zeta3-42.21) successfully discriminates between PGs and other proteins. The model correctly classified 100% of a subset of 81 PGs and 75 non-PG proteins sequences used to train the model. The model also correctly classified 51 out of 52 (98.07%) of proteins sequences used as external validation series. The uses of different group of amino acids and/or axes orientation give different results, so it is suggested to be explored for other databases. Finally, to illustrates the use of the model we report the isolation and prediction of the PG action for a novel sequence (AY908988) isolated by our group from Psidium guajava L. This prediction coincides very well with sequence alignment results found by the BLAST methodology. These findings illustrate the possibilities of the sequence descriptors derived for this novel 2D sequence representation in proteins sequence QSAR studies.
Algorithms↗
Assignment of the 1H and 13C NMR signals of some benzofurocoumarins.
The complete 1H and 13C NMR assignment of four 6,7-benzo-fused furocoumarins (1-4) and three 3,4-benzo-fused furocoumarins (5-7) has been performed using 1D and 2D NMR techniques, including COSY, HMQC and HMBC experiments.
Carbon Isotopes↗
Effect of protein backbone folding on the stability of protein-ligand complexes.
The role played by the degree of folding of protein backbones in explaining the binding energetics of protein-ligand interactions has been studied. We analyzed the protein/peptide interactions in the RNase-S system in which amino acids at two positions of the peptide S have been mutated. The global degree of folding of the protein S correlates in a significant way with the free energy and enthalpy of the protein-peptide interactions. A much better correlation is found with the local contribution to the degree of folding of one amino acid residue: Thr36. This residue is shown to have a destabilizing interaction with Lys41, which interacts directly with peptide S. Another system, consisting of the interactions of small organic molecules with HIV-1 protease was also studied. In this case, the global change in the degree of folding of the protease backbone does not explain the binding energetics of protein-ligand interactions. However, a significant correlation is observed between the free energy of binding and the contribution of two amino acid residues in the HVI-1 protease: Gly49 and Ile66. In general, it was observed that the changes in the degree of folding are not restricted to the binding site of the protein chain but are distributed along the whole protein backbone. This study provides a basis for further consideration of the degree of folding as a parameter for empirical structural parametrizations of the binding energetics of protein folding and binding.
Amino Acids↗
Synthesis and anti-HIV activity of novel cyclopentenyl nucleoside analogues of 8-azapurine.
Novel nucleoside analogues of structure 3-5 were synthesized starting from (+/-)-cis-2-amino-3-cyclopentenylmethanol (1). The chlorine derivative 3 inhibited both HIV-1 and HIV-2 replication in MT-4 cells with IC(50) values of 10.67 microM and of 13.79 microM, respectively.
Anti-HIV Agents↗
Proteins QSAR with Markov average electrostatic potentials.
Classic physicochemical and topological indices have been largely used in small molecules QSAR but less in proteins QSAR. In this study, a Markov model is used to calculate, for the first time, average electrostatic potentials xik for an indirect interaction between aminoacids placed at topologic distances k within a given protein backbone. The short-term average stochastic potential xi1 for 53 Arc repressor mutants was used to model the effect of Alanine scanning on thermal stability. The Arc repressor is a model protein of relevance for biochemical studies on bioorganics and medicinal chemistry. A linear discriminant analysis model developed correctly classified 43 out of 53, 81.1% of proteins according to their thermal stability. More specifically, the model classified 20/28, 71.4% of proteins with near wild-type stability and 23/25, 92.0% of proteins with reduced stability. Moreover, predictability in cross-validation procedures was of 81.0%. Expansion of the electrostatic potential in the series xi0, xi1, xi2, and xi3, justified the use of the abrupt truncation approach, being the overall accuracy >70.0% for xi0 but equal for xi1, xi2, and xi3. The xi1 model compared favorably with respect to others based on D-Fire potential, surface area, volume, partition coefficient, and molar refractivity, with less than 77.0% of accuracy [Ramos de Armas, R.; González-Díaz, H.; Molina, R.; Uriarte, E. Protein Struct. Func. Bioinf.2004, 56, 715]. The xi1 model also has more tractable interpretation than others based on Markovian negentropies and stochastic moments. Finally, the model is notably simpler than the two models based on quadratic and linear indices. Both models, reported by Marrero-Ponce et al., use four-to-five time more descriptors. Introduction of average stochastic potentials may be useful for QSAR applications; having xik amenable physical interpretation and being very effective.
Markov Chains↗
QSAR study for mycobacterial promoters with low sequence homology.
The general belief is that quantitative structure-activity relationship (QSAR) techniques work only for small molecules and, protein sequences or, more recently, DNA sequences. However, with non-branched graph for proteins and DNA sequences the QSAR often have to be based on powerful non-linear techniques such as support vector machines. In our opinion, linear QSAR models based on RNA could be useful to assign biological activity when alignment techniques fail due to low sequence homology. The idea bases the high level of branching for the RNA graph. This work introduces the so-called Markov electrostatic potentials (k)xi(M) as a new class of RNA 2D-structure descriptors. Subsequently, we validate these molecular descriptors solving a QSAR classification problem for mycobacterial promoter sequences (mps), which constitute a very low sequence homology problem. The model developed (mps=-4.664.(0)xi(M)+0. 991.(1)xi(M)-2.432) was intended to predict whether a naturally occurring sequence is an mps or not on the basis of the calculated (k)xi(M) value for the corresponding RNA secondary structure. The RNA-QSAR approach recognises 115/135mps (85.2%) and 100% of control sequences. Average predictability and robustness were greater than 95%. A previous non-linear model predicts mps with a slightly higher accuracy (97%) but uses a very large parameter space for DNA sequences. Conversely, the (k)xi(M)-based RNA-QSAR encodes more structural information and needs only two variables.
Bacterial Proteins↗
Design, synthesis, and vasorelaxant and platelet antiaggregatory activities of coumarin-resveratrol hybrids.
We have synthesized the coumarin-resveratrol hybrid 4 and its dimethoxy derivative 3 by a very direct synthetic route involving a Pechmann procedure. Compound 4 has also been synthesized by an alternative route (Perkin), which also allowed the synthesis of compounds 9-13. In addition, we have evaluated the potential vasorelaxant activity of the new compounds in endothelium-containing rat aorta rings pre-contracted with noradrenaline, as well as the inhibitory effects on platelet aggregation induced by thrombin in washed human platelets. The compounds reported here relaxed vascular smooth muscle and inhibited platelet aggregation with a profile similar to that of trans-resveratrol (t-RESV) and, in some cases, showed activity higher than that of the natural compound. This is the case for compound 13, which has a vasorelaxant activity that is twice as high as that of t-resveratrol and a platelet antiaggregant activity that is six times higher. These results suggest that these novel compounds may have potential as structural templates for the design and subsequent development of new vasodilatory and platelet antiaggregatory drugs.
Animals↗
Stochastic entropy QSAR for the in silico discovery of anticancer compounds: prediction, synthesis, and in vitro assay of new purine carbanucleosides.
A Markov model based QSAR is introduced for the rational selection of anticancer compounds. The model discriminates 90.3% of 226 structurally heterogeneous anticancer/non-anticancer compounds in training series. External validation series were used to validate the model; the 91.8% containing 85 compounds, not considered to fit the model, were correctly classified. The model developed is afterwards used in a simulation of a virtual search for anticancer compounds never considered either in training or in predicting series. The 87.7% of the 213 anticancer compounds used in this simulated search were correctly classified. The model also shows high values for specificity (0.89), sensitivity (0.91), and Mathews correlation coefficient (0.79). In addition, the present model compares better-to-similar with respect to other four models elsewhere reported if one takes into consideration 26 comparison parameters. Finally, we exemplify the use of the model in practice with the design of a new series of carbanucleosides. The compounds evaluated with the model were synthesized and experimentally assayed for their antitumor effects on the proliferation of murine leukemia cells (L1210/0) and human T-lymphocyte cells (CEM/0 and Molt4/C8). The more interesting activity was detected for the compound 5a with a predicted probability of 80.2% and IC(50) = 27.0, 27.2, and 29.4 microM, respectively, against the above-mentioned cellular lines. These values are comparable to those for the control compound Ara-A.
Animals↗
Folding degrees of azurins and pseudoazurins. Implications for structure and function.
A quantitative measure of the degree of folding of azurins and pseudoazurins has been made. We have found that the reduction potential of azurins and pseudoazurins is a function of the contribution to the degree of folding of His117, a key amino acid in electron transfer which is directly bonded to copper in these proteins. The folding degree of His117 explains 95% of the variance in the experimental values of the reduction potential of azurins and pseudoazurins. The change in the folding degree of this amino acid influences several geometric parameters of the main backbones of these proteins. Among them, the angle formed between N(His117)...Cu...S(Cys112), which plays an important role in electron transport, but not the N(His117)...Cu distance, shows some non-linear correlation with the reduction potential of azurins and pseudoazurins. However, it is only able to explain less than 75% in the variance of the reduction potential of these proteins instead of the 95% explained by the folding degree of His117.
Algorithms↗