Study of pK Values and Effective Dielectric Constants of Ionizable Residues in Pentapeptides and in Staphylococcal Nuclease (SNase) Using a Mean-Field Approach

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Towards the validation of "in silico" models and physicochemical filters to identify and characterize new chemical entities

Summary: Lipophilicity plays an important role in the determination and the comprehension of the pharmacokinetic behavior of drugs. It is usually expressed by the partition coefficient (log P) in the n-octanol/water system. The use of an additional solvent system (1,2-dichlorethane/water) is necessary to obtain complementary information, as the log Poct values alone are not sufficient to explain ail biological properties. The aim of this thesis is to develop tools allowing to predict lipophilicity of new drugs and to analyze the information yielded by those log P values. Part I presents the development of theoretical models used to predict lipophilicity. Chapter 2 shows the necessity to extend the existing solvatochromic analyses in order to predict correctly the lipophilicity of new and complex neutral compounds. In Chapter 3, solvatochromic analyses are used to develop a model for the prediction of the lipophilicity of ions. A global model was obtained allowing to estimate the lipophilicity of neutral, anionic and cationic solutes. Part II presents the detailed study of two physicochemical filters. Chapter 4 shows that the Discovery RP Amide C16 stationary phase allows to estimate lipophilicity of the neutral form of basic and acidic solutes, except of lipophilic acidic solutes. Those solutes present additional interactions with this particular stationary phase. In Chapter 5, 4 different IANI stationary phases are investigated. For neutral solutes, linear data are obtained whatever the IANI column used. For the ionized solutes, their retention is due to a balance of electrostatic and hydrophobie interactions. Thus no discrimination is observed between different series of solutes bearing the same charge, from one column to an other. Part III presents two examples illustrating the information obtained thanks to Structure-Properties Relationships (SPR). Comparing graphically lipophilicity values obtained in two different solvent systems allows to reveal the presence of intramolecular effects .such as internai H-bond (Chapter 6). SPR is used to study the partitioning of ionizable groups encountered in Medicinal Chemistry (Chapter7). Résumé La lipophilie joue un .rôle important dans la détermination et la compréhension du comportement pharmacocinétique des médicaments. Elle est généralement exprimée par le coefficient de partage (log P) d'un composé dans le système de solvants n-octanol/eau. L'utilisation d'un deuxième système de solvants (1,2-dichloroéthane/eau) s'est avérée nécessaire afin d'obtenir des informations complémentaires, les valeurs de log Poct seules n'étant pas suffisantes pour expliquer toutes les propriétés biologiques. Le but de cette thèse est de développer des outils permettant de prédire la lipophilie de nouveaux candidats médicaments et d'analyser l'information fournie par les valeurs de log P. La Partie I présente le développement de modèles théoriques utilisés pour prédire la lipophilie. Le chapitre 2 montre la nécessité de mettre à jour les analyses solvatochromiques existantes mais inadaptées à la prédiction de la lipophilie de nouveaux composés neutres. Dans le chapitre 3, la même méthodologie des analyses solvatochromiques est utilisée pour développer un modèle permettant de prédire la lipophilie des ions. Le modèle global obtenu permet la prédiction de la lipophilie de composés neutres, anioniques et cationiques. La Partie II présente l'étude approfondie de deux filtres physicochimiques. Le Chapitre 4 montre que la phase stationnaire Discovery RP Amide C16 permet la détermination de la lipophilie de la forme neutre de composés basiques et acides, à l'exception des acides très lipophiles. Ces derniers présentent des interactions supplémentaires avec cette phase stationnaire. Dans le Chapitre 5, 4 phases stationnaires IAM sont étudiées. Pour les composés neutres étudiés, des valeurs de rétention linéaires sont obtenues, quelque que soit la colonne IAM utilisée. Pour les composés ionisables, leur rétention est due à une balance entre des interactions électrostatiques et hydrophobes. Donc aucune discrimination n'est observée entre les différentes séries de composés portant la même charge d'une colonne à l'autre. La Partie III présente deux exemples illustrant les informations obtenues par l'utilisation des relations structures-propriétés. Comparer graphiquement la lipophilie mesurée dans deux différents systèmes de solvants permet de mettre en évidence la présence d'effets intramoléculaires tels que les liaisons hydrogène intramoléculaires (Chapitre 6). Cette approche des relations structures-propriétés est aussi appliquée à l'étude du partage de fonctions ionisables rencontrées en Chimie Thérapeutique (Chapitre 7) Résumé large public Pour exercer son effet thérapeutique, un médicament doit atteindre son site d'action en quantité suffisante. La quantité effective de médicament atteignant le site d'action dépend du nombre d'interactions entre le médicament et de nombreux constituants de l'organisme comme, par exemple, les enzymes du métabolisme ou les membranes biologiques. Le passage du médicament à travers ces membranes, appelé perméation, est un paramètre important à optimiser pour développer des médicaments plus puissants. La lipophilie joue un rôle clé dans la compréhension de la perméation passive des médicaments. La lipophilie est généralement exprimée par le coefficient de partage (log P) dans le système de solvants (non miscibles) n-octanol/eau. Les valeurs de log Poct seules se sont avérées insuffisantes pour expliquer la perméation à travers toutes les différentes membranes biologiques du corps humain. L'utilisation d'un système de solvants additionnel (le système 1,2-dichloroéthane/eau) a permis d'obtenir les informations complémentaires indispensables à une bonne compréhension du processus de perméation. Un grand nombre d'outils expérimentaux et théoriques sont à disposition pour étudier la lipophilie. Ce travail de thèse se focalise principalement sur le développement ou l'amélioration de certains de ces outils pour permettre leur application à un champ plus large de composés. Voici une brève description de deux de ces outils: 1)La factorisation de la lipophilie en fonction de certaines propriétés structurelles (telle que le volume) propres aux composés permet de développer des modèles théoriques utilisables pour la prédiction de la lipophilie de nouveaux composés ou médicaments. Cette approche est appliquée à l'analyse de la lipophilie de composés neutres ainsi qu'à la lipophilie de composés chargés. 2)La chromatographie liquide à haute pression sur phase inverse (RP-HPLC) est une méthode couramment utilisée pour la détermination expérimentale des valeurs de log Poct.

Caracterização ácido-base da superfície de espécies mistas da alga Spirulina através de titulação potenciométrica e modelo de distribuição de sítios discretos

Acid base properties of mixed species of the microalgae Spirulina were studied by potentiometric titration in medium of 0.01 and 0.10 mols L-1 NaNO3 at 25.0±0.10 C using modified Gran functions or nonlinear regression techniques for data fitting. The discrete site distribution model was used, permitting the characterization of five classes of ionizable sites in both ionic media. This fact suggests that the chemical heterogeneity of the ionizable sites on the cell surface plays a major role on the acid-base properties of the suspension in comparison to electrostatic effects due to charge-charge interactions. The total of ionizable sites were 1.75±0.10 and 1.86±0.20 mmolsg-1 in ionic media of 0.01 and 0.10 mols L-1 NaNO3, respectively. A major contribution of carboxylic groups was observed with an average 34 and 22% of ionizable sites being titrated with conditional pcKa of 4.0 and 5.4, respectively. The remaining 44% of ionizable sites were divided in three classes with averaged conditional pcKa of 6.9, 8.7 and 10.12, which may be assigned respectively to imidazolic, aminic, and phenolic functionalities.

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.

Purification, characterization and sequencing of the major beta-1,3-glucanase from the midgut of Tenebrio molitor larvae

The major beta-1,3-glucanase from Tenebrio molitor (TLam) was purified to homogeneity (yield, 6%; enrichment, 113 fold; specific activity, 4.4 U/mg). TLam has a molecular weight of 50 kDa and a pH optimum of 6. It is an encloglucanase that hydrolyzes beta-1,3-glucans as laminarin and yeast beta-1,3-1,6-glucan, but is inactive toward other polysaccharides (as unbranched beta-1,3-glucans or mixed beta-1,3-1,4-glucan from cereals) or disaccharides. The enzyme is not inhibited by high substrate concentrations and has low processivity (0.6). TLam has two ionizable groups involved in catalysis, and His, Tyr and Arg residues plus a divalent ion at the active site. A Cys residue important for TLam activity is exposed after laminarin binding. The cDNA coding for this enzyme was cloned and sequenced. It belongs to glycoside hydrolase family 16, and is related to other insect glucanases and glucan-binding proteins. Sequence analysis and homology modeling allowed the identification of some residues (E174, E179, H204, Y304, R127 and R181) at the active site of the enzyme, which may be important for TLam activity. TLam efficiently lyses fungal cells, suggesting a role in making available walls and cell contents to digestion and in protecting the midgut from pathogen infections. (C) 2009 Elsevier Ltd. All rights reserved.

Determinação do grau de ionização de aminoácidos polares carregados

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Determinação do gradiente de pH na região interna a duas membranas em equilíbrio com eletrólito

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

pH-dependent distribution of chlorin e6 derivatives across phospholipid bilayers probed by NMR spectroscopy

The pH-dependent membrane adsorption and distribution of three chlorin derivatives, chlorin e6 (CE), rhodin G7 (RG), and monoaspartyl-chlorin e6 (MACE), in the physiological pH range (pH 6-8) were probed by NMR spectroscopy. Unilamellar vesicles consisting of dioleoyl-phosphatidyl-choline (DOPC) were used as membrane models. The chlorin derivatives were characterized with respect to their aggregation behavior, the pK(a) values of individual carboxylate groups, the extent of membrane adsorption, and their flip-flop rates across the bilayer membrane for pH 6-8. External membrane adsorption was found to be lower for RG than for CE and MACE. Both electrostatic interactions and the extent of aggregation seemed to be the main determinants of membrane adsorption. Rate constants for chlorin transfer across the membrane were found to correlate strongly with the pH of the surrounding medium, in particular, for CE and RG. In acidic solution, CE and RG transfer across the membrane was strongly accelerated, and in basic solution, all compounds were retained, mostly in the outer monolayer. In contrast, MACE flip-flop across the membrane remained very low even at pH 6. The protonation of ionizable groups is suggested to be a major determinant of chlorin transfer rates across the bilayer. pK(a) values of CE and RG were found to be between 6 and 8, and two of the carboxylate groups in MACE had pK(a) values below 6. For CE and RG, the kinetic profiles at acidic pH indicated that the initial fast membrane distribution was followed by secondary steps that are discussed in this article.

Role of hydrophobic interactions and desolvation in determining the structural properties of a model alpha beta peptide.

Model AB, a 20-amino acid peptide that was designed to adopt an alpha beta tertiary structure stabilized by hydrophobic interactions between residues in adjacent helical and extended segments, exhibited large pKa shifts of several ionizable groups and slow hydrogen/deuterium exchange rates of nearly all the peptide amide groups [Butcher, D. J., Bruch, M. D. & Moe, G. T. (1995) Biopolymers 36, 109-120]. These properties, which depend on structure and hydration, are commonly observed in larger proteins but are quite unusual for small peptides. To identify which of several possible features of the peptide design are most important in determining these properties, several closely related analogs of Model AB were characterized by CD and NMR spectroscopy. The results show that hydrophobic interactions between adjacent helical and extended segments are structure-determining and have the additional effect of altering water-peptide interactions over much of the peptide surface. These results may have important implications for understanding mechanisms of protein folding and for the design of independently folding peptides.

Electrostatic Contributions to the Thermodynamics of Ribonuclease P Protein Folding

Electrostatic interactions are of fundamental importance in determining the structure and stability of macromolecules. For example, charge-charge interactions modulate the folding and binding of proteins and influence protein solubility. Electrostatic interactions are highly variable and can be both favorable and unfavorable. The ability to quantify these interactions is challenging but vital to understanding the detailed balance and major roles that they have in different proteins and biological processes. Measuring pKa values of ionizable groups provides a sensitive method for experimentally probing the electrostatic properties of a protein.

pKa values report the free energy of site-specific proton binding and provide a direct means of studying protein folding and pH-dependent stability. Using a combination of NMR, circular dichroism, and fluorescence spectroscopy along with singular value decomposition, we investigated the contributions of electrostatic interactions to the thermodynamic stability and folding of the protein subunit of Bacillus subtilis ribonuclease P, P protein. Taken together, the results suggest that unfavorable electrostatics alone do not account for the fact that P protein is intrinsically unfolded in the absence of ligand because the pKa differences observed between the folded and unfolded state are small. Presumably, multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of p H and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups.

A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: Proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of pH and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups. © 2011 American Institute of Physics.

A direct potentiometric titration study of the dissociation of humic acid with selectively blocked functional groups

A direct potentiometric titration method was applied to commercial and soil humic acids in order to determine their carboxyl and phenol group concentrations and apparent and intrinsic pK. In that context, acid-base properties of humic acids are interpreted by selective blocking of carboxylic and phenolic groups by esterification and acetylation. Differences in underivatized and derivatized HA's acid-base properties are ascribed to carboxyl and phenol groups influence on total humic acidity. Potentiometric data were treated with the modified Henderson-Hasselbalch equation. Infra red results, the acidic group contents and the average values of apparent and intrinsic pK for underivatized and derivatized HAs confirmed the selectivity of esterification derivatization method. After blocking of the functional groups, the values of acidic group contents decreased, while the value of apparent pK increased after derivatization. Phenol groups cannot be specifically identified by the acetylation method, due to low selectivity of the acetylation method.

Design, Synthesis And Biological Evaluation Of Hybrid Bioisoster Derivatives Of N-acylhydrazone And Furoxan Groups With Potential And Selective Anti-trypanosoma Cruzi Activity.

Hybrid bioisoster derivatives from N-acylhydrazones and furoxan groups were designed with the objective of obtaining at least a dual mechanism of action: cruzain inhibition and nitric oxide (NO) releasing activity. Fifteen designed compounds were synthesized varying the substitution in N-acylhydrazone and in furoxan group as well. They had its anti-Trypanosoma cruzi activity in amastigotes forms, NO releasing potential and inhibitory cruzain activity evaluated. The two most active compounds (6, 14) both in the parasite amastigotes and in the enzyme contain the nitro group in para position of the aromatic ring. The permeability screening in Caco-2 cell and cytotoxicity assay in human cells were performed for those most active compounds and both showed to be less cytotoxic than the reference drug, benznidazole. Compound 6 was the most promising, since besides activity it showed good permeability and selectivity index, higher than the reference drug. Thereby the compound 6 was considered as a possible candidate for additional studies.

Comparative Cytogenetics Of Physalaemus albifrons And Physalaemus Cuvieri Species Groups (anura, leptodactylidae).

Recently, Physalaemus albifrons (Spix, 1824) was relocated from the Physalaemus cuvieri group to the same group as Physalaemus biligonigerus (Cope, 1861), Physalaemus marmoratus (Reinhardt & Lütken, 1862) and Physalaemus santafecinus Barrio, 1965. To contribute to the analysis of this proposition, we studied the karyotypes of Physalaemus albifrons, Physalaemus santafecinus and three species of the Physalaemus cuvieri group. The karyotype of Physalaemus santafecinus was found to be very similar to those of Physalaemus biligonigerus and Physalaemus marmoratus, which were previously described. A remarkable characteristic that these three species share is a conspicuous C-band that extends from the pericentromeric region almost to the telomere in the short arm of chromosome 3. This characteristic is not present in the Physalaemus albifrons karyotype and could be a synapomorphy of Physalaemus biligonigerus, Physalaemus marmoratus and Physalaemus santafecinus. The karyotype of Physalaemus santafecinus is also similar to those of Physalaemus marmoratus and Physalaemus biligonigerus owing to the presence of several terminal C-bands and the distal localization of the NOR in a small metacentric chromosome. In contrast, the Physalaemus albifrons karyotype has no terminal C-bands and its NOR is located interstitially in the long arm of submetacentric chromosome 8. The NOR-bearing chromosome of Physalaemus albifrons very closely resembles those found in Physalaemus albonotatus (Steindachner, 1864), Physalaemus cuqui Lobo, 1993 and some populations of Physalaemus cuvieri Fitzinger, 1826. Additionally, the Physalaemus albifrons karyotype has an interstitial C-band in chromosome 5 that has been exclusively observed in species of the Physalaemus cuvieri group. Therefore, we were not able to identify any chromosomal feature that supports the reallocation of Physalaemus albifrons.