Melting Regime of the Anionic Phospholipid DMPG: New Lamellar Phase and Porous Bilayer Model

Aqueous dispersions of the anionic phospholipid dimyristoyl phosphatidylglycerol (DMPG) at pH above the apparent pK of DMPG and concentrations in the interval 70-300 mM have been investigated by small (SAXS) and wide-angle X-ray scattering, differential scanning calorimetry, and polarized optical microscopy. The order. disorder transition of the hydrocarbon chains occurs along an interval of about 10 degrees C (between T(m)(on) similar to 20 degrees C and T(m)(off) similar to 30 degrees C). Such melting regime was previously characterized at lower concentrations, up to 70 mM DMPG, when sample transparency was correlated with the presence of pores across the bilayer. At higher concentrations considered here, the melting regime persists but is not transparent. Defined SAXS peaks appear and a new lamellar phase L(p) with pores is proposed to exist above 70 mM DMPG, starting at similar to 23 degrees C (similar to 3 degrees C above T(m)(on)) and losing correlation after T(m)(off). A new model for describing the X-ray scattering of bilayers with pores, presented here, is able to explain the broad band attributed to in-plane correlation between pores. The majority of cell membranes have a net negative charge, and the opening of pores across the membrane tuned by ionic strength, temperature, and lipid composition is likely to have biological relevance.

First observation of the Cabibbo-suppressed decays Xi(+)(c)->Sigma(+)pi(-)pi(+) and Xi(+)(c)->Sigma(-)pi(+)pi(+) and measurement of their branching ratios

We report the first observation of two Cabibbo-suppressed ecay modes, Xi(+)(c) -> Sigma(+)pi(-)pi(+) and Xi c+ -> Sigma(-)pi(+)pi(+). We observe 59 +/- 14 over a background of 87, and 22 +/- 8 over a background of 13 events, respectively, for the signals. The data were accumulated using the SELEX spectrometer during the 1996-1997 fixed target run at Fermilab, chiefly from a 600 Gev/c Sigma(-) beam. The branching ratios of the decays relative to the Cabibbo-favored Xi c+ -> Xi(-)pi(+)pi(+) are measured to be B(Xi(+)(c) -> Sigma(+)pi(-)pi(+))/B(Xi(+)(c) -> Xi(-)pi(+)pi(+)) = 0.48 +/- 0.20, and B(Xi(+)(c) -> Sigma(-)pi(+)pi(+))/B(Xi(+)(c) -> Sigma(-)pi(+)pi(+)) = 0.18 +/- 0.09, respectively. We also report branching ratios for the same decay modes of the Delta(+)(c) relative to Delta(+)(c) -> pK(-)pi(+.) (C) 2008 Elsevier B.V. All rights reserved.

Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Catalytic mechanism of inulinase from Arthrobacter sp S37

Detailed catalytic roles of the conserved Glu323, Asp460, and Glu519 of Arthrobacter sp. S37 inulinase (EnIA), a member of the glycoside hydrolase family 32, were investigated by site-directed mutagenesis and pH-dependence studies of the enzyme efficiency and homology modeling were carried out for EnIA and for D460E mutant. The enzyme efficiency (k(cat)/K-m) of the E323A and E519A mutants was significantly lower than that of the wild-type due to a substantial decrease in k(cat), but not due to variations in K-m, consistent with their putative roles as nucleophile and acid/base catalyst, respectively. The D460A mutant was totally inactive, whereas the D460E and D460N mutants were active to some extent, revealing Asp460 as a catalytic residue and demonstrating that the presence of a carboxylate group in this position is a prerequisite for catalysis. The pH-dependence studies indicated that the pK(a) of the acid/base catalyst decreased from 9.2 for the wild-type enzyme to 7.0 for the D460E mutant, implicating Asp460 as the residue that interacts with the acid/base catalyst Glu519 and elevates its pK(a). Homology modeling and molecular dynamics simulation of the wild-type enzyme and the D460E mutant shed light on the structural roles of Glu323, Asp460, and Glu519 in the catalytic activity of the enzyme. (C) 2008 Elsevier Inc. All rights reserved.

pH-Sensitive Binding of Cytochrome c to the Inner Mitochondrial Membrane. Implications for the Participation of the Protein in Cell Respiration and Apoptosis

Cytochrome c exhibits two positively charged sites: site A containing lysine residues with high pK(a) values and site L containing ionizable groups with pK(aobs),values around 7.0. This protein feature implies that cytochrome c can participate in the fusion of mitochondria and have its detachment from the inner membrane regulated by cell acidosis and alkalosis. In this study, We demonstrated that both horse and tuna cytochrome c exhibited two types of binding to inner mitochondrial membranes that contributed to respiration: a high-affinity and low-efficiency pi-I-independent binding (microscopic dissociation constant K(sapp2), similar to 10 nM) and a low-affinity and high-efficiency pH-dependent binding that for horse cytochrome c had a pK(a) of similar to 6.7. For tuna cytochrome c (Lys22 and His33 replaced with Asn and Trp, respectively), the effect of pH on K(sapp1), was less striking than for the horse heme protein, and both tuna and horse cytochrome c had closed K(sapp1) values at pH 7.2 and 6.2, respectively. Recombinant mutated cytochrome c H26N and H33N also restored the respiration of the cytochrome c-depleted mitoplast in a pH-dependent manner. Consistently, the detachment of cytochrome c from nondepleted mitoplasts was favored by alkalinization, suggesting that site Lionization influences the participation of cytochrome c in the respiratory chain and apoptosis.

Effects of Micelles and Vesicles on the Oximolysis of p-Nitrophenyl Diphenyl Phosphate: A Model System for Surfactant-Based Skin-defensive Formulations against Organophosphates

The rates of oximolysis of p-nitrophenyl diphenyl phosphate (PNPDPP) by Acetophenoxime; 10-phenyl-10-hydi-oxyiminodecanoic acid; 4-(9-carboxynonanyl)-1-(9-carboxy-1-hydroyiminononanyl) benzene; 1-dodecyl-2-[(hydroxyimino)methyl]-pyridinium chloride (IV) and N-methylpyridinium-2-aldoxime chloride were determined in micelles of N-hexadecyl-N,N,N-trimethylammonium chloride (CTAC), N-hexadecyl-N,N-dimethylammonium propanesulfonate and dioctadecyldimethylammonium chloride (DODAC) vesicles. The effects of CTAC micelles and DODAC vesicles on the rates of oxymolysis of O,O-Diethyl O-(4-nitrophenyl) phosphate (paraoxon) by oxime IV were also determined. Analysis of micellar and vesicular effects on oximolysis of PNPDPP, using pseudophase or pseudophase with explicit consideration of ion exchange models, required the determination of the aggregate`s effects on the pK(a), of oximes and on the rates of PNPDPP hydrolysis. All aggregates increased the rate of oximolysis of PNPDPP and the results were analyzed quantitatively. In particular, DODAC vesicles catalyzed the reaction and increased the rate of oximolysis of PNPDPP by IV several million fold at pH`s compatible with pharmaceutical formulations. The rate increase produced by DODAC vesicles on the rate of oximolysis paraoxon by IV demonstrates the pharmaceutical potential of this system, since the substrate is used as an agricultural defensive agent and the surfactant is extensively employed in cosmetic formulations. (C) 2008 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 98:1040-1052, 2009

1,4-Diamino-2-butanone, a wide-spectrum microbicide, yields reactive species by metal-catalyzed oxidation

The alpha-aminoketone 1,4-diamino-2-butanone (DAB), a putrescine analogue, is highly toxic to various microorganisms, including Trypanosoma cruzi. However, little is known about the molecular mechanisms underlying DAB`s cytotoxic properties. We report here that DAB (pK(a) 7.5 and 9.5) undergoes aerobic oxidation in phosphate buffer, pH 7.4, at 37 degrees C, catalyzed by Fe(II) and Cu(II) ions yielding NH(4)(+) ion, H(2)O(2), and 4-amino-2-oxobutanal (oxoDAB). OxoDAB, like methylglyoxal and other alpha-oxoaldehydes, is expected to cause protein aggregation and nucleobase lesions. Propagation of DAB oxidation by superoxide radical was confirmed by the inhibitory effect of added SOD (50 U ml(-1)) and stimulatory effect of xanthine/xanthine oxidase, a source of superoxide radical. EPR spin trapping studies with 5,5-dimethyl-1-pyrroline-1-oxide (DMPO) revealed an adduct attributable to DMPO-HO(center dot), and those with alpha-(4-pyridyl-1-oxide)-N-tert-butylnitrone or 3,5-dibromo-4-nitrosobenzenesulfonic acid, a six-line adduct assignable to a DAB(center dot) resonant enoyl radical adduct. Added horse spleen ferritin (HoSF) and bovine apo-transferrin underwent oxidative changes in tryptophan residues in the presence of 1.0-10 mM DAB. Iron release from HoSF was observed as well. Assays performed with fluorescein-encapsulated liposomes of cardiolipin and phosphatidylcholine (20:80) incubated with DAB resulted in extensive lipid peroxidation and consequent vesicle permeabilization. DAB (0-10 mM) administration to cultured LLC-MK2 epithelial cells caused a decline in cell viability, which was inhibited by preaddition of either catalase (4.5 mu M) or aminoguanidine (25 mM). Our findings support the hypothesis that DAB toxicity to several pathogenic microorganisms previously described may involve not only reported inhibition of polyamine metabolism but also DAB pro-oxidant activity. (C) 2011 Elsevier Inc. All rights reserved.

Horseradish peroxidase compound I as a tool to investigate reactive protein-cysteine residues: from quantification to kinetics

Proteins containing reactive cysteine residues (protein-Cys) are receiving increased attention as mediators of hydrogen peroxide signaling. These proteins are mainly identified by mining the thiol proteomes of oxidized protein-Cys in cells and tissues. However, it is difficult to determine if oxidation occurs through a direct reaction with hydrogen peroxide or by thiol-disulfide exchange reactions. Kinetic studies with purified proteins provide invaluable information about the reactivity of protein-Cys residues with hydrogen peroxide. Previously, we showed that the characteristic UV-Vis spectrum of horseradish peroxidase compound I, produced from the oxidation of horseradish peroxidase by hydrogen peroxide, is a simple, reliable, and useful tool to determine the second-order rate constant of the reaction of reactive protein-Cys with hydrogen peroxide and peroxynitrite. Here, the method is fully described and extended to quantify reactive protein-Cys residues and micromolar concentrations of hydrogen peroxide. Members of the peroxiredoxin family were selected for the demonstration and validation of this methodology. In particular, we determined the pK(a) of the peroxidatic thiol of rPrx6 (5.2) and the second-order rate constant of its reactions with hydrogen peroxide ((3.4 +/- 0.2) x 10(7) M(-1) s(-1)) and peroxynitrite ((3.7 +/- 0.4) x 10(5) M(-1) s(-1)) at pH 7.4 and 25 degrees C. (C) 2011 Elsevier Inc. All rights reserved.

Acetyl Radical Production by the Methylglyoxal-Peroxynitrite System: A Possible Route for L-Lysine Acetylation

Methylglyoxal is an a-oxoaldehyde putatively produced in excess from triose phosphates, aminoacetone, and acetone in some disorders, particularly in diabetes. Here, we investigate the nucleophilic addition of ONOO(-), known as a potent oxidant and nucleophile, to methylglyoxal, yielding an acetyl radical intermediate and ultimately formate and acetate ions. The rate of ONOO(-) decay in the presence of methylglyoxal [k(2,app) = (1.0 +/- 0.1) x 10(3) M(-1) s(-1); k(2) approximate to 1.0 x 10(5) M(-1) s(-1)] at pH 7.2 and 25 degrees C was found to be faster than that reported with monocarbonyl substrates (k(2) < 10(3) M(-1) diacetyl (k(2) = 1.0 x 10(4) M(-1) s(-1)), or CO(2) (k(2) = 3-6 x 10(4) M(-1) s(-1)). The pH profile of the methylglyoxal peroxynitrite reaction describes an ascendant curve with an inflection around pH 7.2, which roughly coincides with the pK(a) values of both ONOOH and H(2)PO(4)(-) ion. Electron paramagnetic resonance spin trapping experiments with 2-methyl-2-nitrosopropane revealed concentration-dependent formation of an adduct that can be attributed to 2-methyl-2-nitrosopropane-CH(3)CO(center dot) (a(N) = 0.83 mT). Spin trapping with 3,5-dibromo-4-nitrosobenzene sulfonate gave a signal that could be assigned to a methyl radical adduct [a(N) = 1.41 mT; a(H) = 1.35 mT; a(H(m)) = 0.08 mT]. The 2-methyl-2-nitrosopropane-CH(3)CO(center dot) adduct could also be observed by replacement of ONOO(-) with H(2)O(2), although at much lower yields. Acetyl radicals could be also trapped by added L-lysine as indicated by the presence of W-acetyl-L-lysine in the spent reaction mixture. This raises the hypothesis that ONOO(-)/H(2)O(2) in the presence of methylglyoxal is endowed with the potential to acetylate proteins in post-translational processes.

Investigating the formation and the properties of monoalkyl carbonates in aqueous medium using capillary electrophoresis with capacitively coupled contactless conductivity detection

Although alkyl carbonic acids (ACAs) and their salts are referred to as instable species in aqueous medium, we demonstrate that a monoalkyl carbonate (MAC) can in fact be easily formed from bicarbonate and an alcohol even in the presence of a high amount of water. A CE system with two capacitively coupled contactless conductivity detectors (C(4)Ds) was used to obtain different parameters about these species and their reactions. Based on the mobilities obtained for a series of alcohols ranging from 1 to 5 carbons, the coefficients of diffusion and the hydrodynamic radii were calculated. When compared with the equivalent carboxylates, MACs have radii systematically smaller. Although the precise pK(a) values of the ACAs could not be obtained, because of the fast decomposition in acid medium, it was possible, for the first time, to show that they are below 4.0. This result suggests that the acidity of an ACA is quite similar to the first hydrogen of H(2)CO(3). Using a new approach to indirectly calibrate the C(4)D, the kinetic constants and the equilibrium constants of formation were also obtained. The results suggest that the increase in the chain length makes the MACs less stable and more inert.

Solvatochromism in Binary Mixtures: First Report on a Solvation Free Energy Relationship between Solvent Exchange Equilibrium Constants and the Properties of the Medium

We have employed UV-vis spectroscopy in order to investigate details of the solvation of six solvatochromic indicators, hereafter designated as ""probes"", namely, 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr,, respectively. These can be divided into three pairs, each includes two probes of similar pK(a) in water and different lipophilicity. Solvation has been studied in binary mixtures, BMs, of water, W, with 12 protic organic solvents, S, including mono- and bifunctional alcohols (2-alkoxyethanoles, unsaturated and chlorinated alcohols). Each medium was treated as a mixture of S, W, and a complex solvent, S-W, formed by hydrogen bonding. Values of lambda(max) (of the probe intramolecular charge transfer) were converted into empirical polarity scales, E(T)(probe) in kcal/mol, whose values were correlated with the effective mole fraction of water in the medium, chi w(effective). This correlation furnished three equilibrium constants for the exchange of solvents in the probe solvation shell; phi(W/S) (W substitutes S): phi(S-W/W) (S-W substitutes W), and phi(S-W/S) (S-W substitutes S), respectively. The values of these constants depend on the physicochemical properties of the probe and the medium. We tested, for the first time, the applicability of a new solvation free energy relationship: phi = constant + a alpha(BM) + b beta(BM) + s(pi*(BM) + d delta) + p log P(BM), where a, b, s, and p are regression coefficients alpha(BM), beta(BM), and pi*(BM) are solvatochromic parameters of the BM, delta is a correction term for pi*, and log P is an empirical scale of lipophilicity. Correlations were carried out with two-, three-, and four-medium descriptors. In all cases, three descriptors gave satisfactory correlations; use of four parameters gave only a marginal increase of the goodness of fit. For phi(W/S), the most important descriptor was found to be the lipophilicity of the medium; for phi(S-W/W) and phi(S-W/S), solvent basicity is either statistically relevant or is the most important descriptor. These responses are different from those of E(T)(probe) of many solvatochromic indicators in pure solvents, where the importance of solvent basicity is usually marginal, and can be neglected.

Solvation in Pure Liquids: What Can Be Learned from the Use of Pairs of Indicators?

The solvation of six solvatochromic probes in a large number of solvents (33-68) was examined at 25 degrees C. The probes employed were the following: 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenyl pyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2), respectively. Of these, MePMBr is a novel compound. They can be grouped in three pairs, each with similar pK(a) in water but with different molecular properties, for example, lipophilicity and dipole moment. These pairs are formed by RB and MePM; QB and MePMBr; WB and MePMBr(2), respectively. Theoretical calculations were carried out in order to calculate their physicochemical properties including bond lengths, dihedral angles, dipole moments, and wavelength of absorption of the intramolecular charge-transfer band in four solvents, water, methanol, acetone, and DMSO, respectively. The data calculated were in excellent agreement with available experimental data, for example, bond length and dihedral angles. This gives credence to the use of the calculated properties in explaining the solvatochromic behaviors observed. The dependence of an empirical solvent polarity scale E(T)(probe) in kcal/mol on the physicochemical properties of the solvent (acidity, basicity, and dipolarity/polarizability) and those of the probes (pK(a), and dipole moment) was analyzed by using known multiparameter solvation equations. For each pair of probes, values of E(T)(probe) (for example, E(T)(MePM) versus E(T)(RB)) were found to be linearly correlated with correlation coefficients, r, between 0.9548 and 0.9860. For the mercyanine series, the values of E(T)(probe) also correlated linearly, with (r) of 0.9772 (MePMBr versus MePM) and 0.9919 (MePMBr(2) versus MePM). The response of each pair of probes (of similar pK(a)) to solvent acidity is the same, provided that solute-solvent hydrogen-bonding is not seriously affected by steric crowding (as in case of RB). We show, for the first time, that the response to solvent dipolarity/polarizability is linearly correlated to the dipole moment of the probes. The successive introduction of bromine atoms in MePM (to give MePMBr, then MePMBr(2)) leads to the following linear decrease: pK(a) in water, length of the phenolate oxygen-carbon bond, length of the central ethylenic bond, susceptibility to solvent acidity, and susceptibility to solvent dipolarity/polarizability. Thus studying the solvation of probes whose molecular structures are varied systematically produces a wealth of information on the effect of solute structure on its solvation. The results of solvation of the present probes were employed in order to test the goodness of fit of two independent sets of solvent solvatochromic parameters.

Ion pair stabilization effects on a series of procaine structural analogs

In this work, a series of 10 structural procaine analogs have been synthesized in order to investigate the structural features affecting the stability of ion pair formation and its influence on the lipophilicity of ionizable compounds. The structural variation within this series was focused on the terminal nitrogen substituents and on the intermediate chain linkage nature. The hydrophobic parameters log P(n) and log P(i) (partition coefficient of the neutral and ionic species, respectively), as well as the ionization constants pK(a) and pK(a)(oct), were obtained from log D-pH profiles measured at pH values ranging from 2 to 12. The difference between log P(i) and log P(n) values (i.e. difflog P) of each prepared compound was considered a measure of the stability of ion pair formation. In this set, the difflog P values varied nearly over one log unit, ranging from -2.40 to -3.37. It has been observed that the presence of hydrogen bonding groups (especially donor) and low steric hindrance around the terminal amine ionizable group increases the relative lipophilicity of the ionic species as compared to the corresponding neutral species. These results were interpreted as due to the increased stability of ion pairs of the compounds bearing these structural features. (C) 2010 Elsevier B.V. All rights reserved.

Light-Induced Oxidation of Unsaturated Lipids as Sensitized by Flavins

Triplet-excited riboflavin ((3)RF*) was found by laser flash photolysis to be quenched by polyunsaturated fatty acid methyl esters in tert-butanol/water (7:3, v/v) in a second-order reaction with k similar to 3.0 x 10(5) L mol(-1) s(-1) at 25 degrees C for methyl linoleate and 3.1 x 10(6) L mol(-1) s(-1), with Delta H double dagger = 22.6 kJ mol(-1) and Delta S double dagger = -62.3 J K(-1) mol(-1), for methyl linolenate in acetonitrile/water (8:2, v/v). For methyl oleate, k was <10(4) L mol(-1) s(-1). For comparison, beta-casein was found to have a rate constant k similar to 4.9 x 10(8) L mol(-1) s(-1). Singlet-excited flavin was not quenched by the esters as evidenced by insensitivity of steady-state fluorescence to their presence. Density functional theory (DFT) calculations showed that electron transfer from unsaturated fatty acid esters to triplet-excited flavins is endergonic, while a formal hydrogen atom transfer is exergonic (Delta G(HAT)degrees = -114.3, -151.2, and -151.2 kJ mol(-1) for oleate, linoleate, and linolenate, respectively, in acetonitrile). The reaction is driven by acidity of the lipid cation radical for which a pK(a) similar to -0.12 was estimated by DFT calculations. Absence of electrochemical activity in acetonitrile during cyclic voltammetry up to 2.0 V versus NHE confirmed that Delta G(ET)degrees > 0 for electron transfer. Interaction of methyl esters with (3)RF* is considered as initiation of the radical chain, which is subsequently propagated by combination reactions with residual oxygen. In this respect, carbon-centered and alkoxyl radicals were detected using the spin trapping technique in combination with electron paramagnetic resonance spectroscopy. Moreover, quenching of 3RF* yields, directly or indirectly, radical species which are capable of initiating oxidation in unsaturated fatty acid methyl esters. Still, deactivation of triplet-excited flavins by lipid derivatives was slower than by proteins (factor up to 10(4)), which react preferentially by electron transfer. Depending on the reaction environment in biological systems (including food), protein radicals are expected to interfere in the mechanism of light-induced lipid oxidation.

A quantum chemical study on a set of non-imidazole H(3) antihistamine molecules

Molecular orbital calculations were carried out on a set of 28 non-imidazole H(3) antihistamine compounds using the Hartree-Fock method in order to investigate the possible relationships between electronic structural properties and binding affinity for H3 receptors (pK(i)). It was observed that the frontier effective-for-reaction molecular orbital (FERMO) energies were better correlated with pK(i) values than highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values. Exploratory data analysis through hierarchical cluster (HCA) and principal component analysis (PCA) showed a separation of the compounds in two sets, one grouping the molecules with high pK(i) values, the other gathering low pK(i) value compounds. This separation was obtained with the use of the following descriptors: FERMO energies (epsilon(FERMO)), charges derived from the electrostatic potential on the nitrogen atom (N(1)), electronic density indexes for FERMO on the N(1) atom (Sigma((FERMO))c(i)(2)). and electrophilicity (omega`). These electronic descriptors were used to construct a quantitative structure-activity relationship (QSAR) model through the partial least-squares (PLS) method with three principal components. This model generated Q(2) = 0.88 and R(2) = 0.927 values obtained from a training set and external validation of 23 and 5 molecules, respectively. After the analysis of the PLS regression equation and the values for the selected electronic descriptors, it is suggested that high values of FERMO energies and of Sigma((FERMO))c(i)(2), together with low values of electrophilicity and pronounced negative charges on N(1) appear as desirable properties for the conception of new molecules which might have high binding affinity. 2010 Elsevier Inc. All rights reserved.