Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 °C in 0.1 M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge–dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds. A primary objective in this study was to investigate the influence of charge–charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the α-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge–dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge–charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pKSa−pK5K) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson–Boltzmann equation, which considers desolvation and charge–dipole interactions in addition to charge–charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge–charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.This work was supported by grants GM-37039 and GM-52483 from the National Institutes of Health (USA), grants BE-1060 and BE-1281 from the Robert A. Welch Foundation, and a grant PB-93-06777 to M.R. from the Dirección General de Investigación Cientı́fica y Técnica (Spain

Theoretically nanoscale study on ionization of muscimol nano drug in aqueous solution

In the present work, acid dissociation constant (pKa) values of muscimol derivatives were calculated using the Density Functional Theory (DFT) method. In this regard, free energy values of neutral, protonated and deprotonated species of muscimol were calculated in water at the B3LYP/6-31G(d) basis sets. The hydrogen bond formation of all species had been analyzed using the Tomasi's method. It was revealed that the theoretically calculated pKa values were in a good agreement with the existing experimental pKa values, which were determined from capillary electrophoresis, potentiometric titration and UV-visible spectrophotometric measurements

Common and Distant Structural Characteristics of Feruloyl Esterase Families from Aspergillus oryzae

Background: Feruloyl esterases (FAEs) are important biomass degrading accessory enzymes due to their capability of cleaving the ester links between hemicellulose and pectin to aromatic compounds of lignin, thus enhancing the accessibility of plant tissues to cellulolytic and hemicellulolytic enzymes. FAEs have gained increased attention in the area of biocatalytic transformations for the synthesis of value added compounds with medicinal and nutritional applications. Following the increasing attention on these enzymes, a novel descriptor based classification system has been proposed for FAEs resulting into 12 distinct families and pharmacophore models for three FAE sub-families have been developed. Methodology/Principal Findings: The feruloylome of Aspergillus oryzae contains 13 predicted FAEs belonging to six sub-families based on our recently developed descriptor-based classification system. The three-dimensional structures of the 13 FAEs were modeled for structural analysis of the feruloylome. The three genes coding for three enzymes, viz., A.O.2, A.O.8 and A.O.10 from the feruloylome of A. oryzae, representing sub-families with unknown functional features, were heterologously expressed in Pichia pastoris, characterized for substrate specificity and structural characterization through CD spectroscopy. Common feature-based pharamacophore models were developed according to substrate specificity characteristics of the three enzymes. The active site residues were identified for the three expressed FAEs by determining the titration curves of amino acid residues as a function of the pH by applying molecular simulations. Conclusions/Significance: Our findings on the structure-function relationships and substrate specificity of the FAEs of A. oryzae will be instrumental for further understanding of the FAE families in the novel classification system. The developed pharmacophore models could be applied for virtual screening of compound databases for short listing the putative substrates prior to docking studies or for post-processing docking results to remove false positives. Our study exemplifies how computational predictions can complement to the information obtained through experimental methods. © 2012 Udatha et al.published_or_final_versio

Determination of the pK(a) of the N-terminal amino group of ubiquitin by NMR

This work was supported by the Medical Research Council (grants U117533887 and grant U117565398 until March 2015) and by the Francis Crick Institute which receives its core funding from Cancer Research UK (FC001142, FC10029), the UK Medical Research Council (FC001142, FC10029), and the Wellcome Trust (FC001142, FC10029)

pK values of the ionizable groups of proteins

We have used potentiometric titrations to measure the pK values of the ionizable groups of proteins in alanine pentapeptides with appropriately blocked termini. These pentapeptides provide an improved model for the pK values of the ionizable groups in proteins. Our pK values determined in 0.1 M KCl at 25°C are: 3.67±0.03 (α-carboxyl), 3.67±0.04 (Asp), 4.25±0.05 (Glu), 6.54±0.04 (His), 8.00±0.03 (α-amino), 8.55±0.03 (Cys), 9.84±0.11 (Tyr), and 10.40±0.08 (Lys). The pK values of some groups differ from the Nozaki and Tanford (N&T) pK values often used in the literature: Asp (3.67 this work vs. 4.0 N&T); His (6.54 this work vs. 6.3 N&T); α-amino (8.00 this work vs. 7.5 N&T); Cys (8.55 this work vs. 9.5 N&T); and Tyr (9.84 this work vs. 9.6 N&T). Our pK values will be useful to those who study pK perturbations in folded and unfolded proteins, and to those who use theory to gain a better understanding of the factors that determine the pK values of the ionizable groups of proteins

The basic helix-loop-helix region of human neurogenin 1 is a monomeric natively unfolded protein which forms a

Neuronal specification is regulated by the activity of transcription factors containing the basic helix-loop-helix motif (bHLH); these regulating proteins include, among others, the neurogenin (Ngn) family, related to the atonal family of genes. Neurogenin 1 (NGN1) is a 237-residue protein that contains a bHLH domain and is involved in neuronal differentiation. In this work, we synthesized the bHLH region of NGN1 (bHLHN) comprising residues 90-150 of the full-length NGN1. The domain is a monomeric natively unfolded protein with a pH-dependent premolten globule conformation, as shown by several spectroscopic techniques (namely, NMR, fluorescence, FTIR, and circular dichroism). The unfolded character of the domain also explains, first, the impossibility of its overexpression in several Escherichia coli strains and, second, its insolubility in aqueous buffers. To the best of our knowledge, this is the first extensive study of the conformational preferences of a bHLH domain under different solution conditions. Upon binding to two DNA E-boxes, the protein forms "fuzzy" complexes (that is, the complexes were not fully folded). The affinities of bHLHN for both DNA boxes were smaller than those of other bHLH domains, which might explain why the protein-DNA complexes were not fully folded.Journal ArticleResearch Support, Non-U.S. Gov'tinfo:eu-repo/semantics/publishe