Titration_DB: storage and analysis of NMR-monitored protein pH titration curves

The influence of random-coil chemical shifts on the assessment of structural propensities in folded proteins and IDPs

In studying secondary structural propensities of proteins by nuclear magnetic resonance (NMR) spectroscopy, secondary chemical shifts (SCSs) serve as the primary atomic scale observables. For SCS calculation, the selection of an appropriate random coil chemical shift (RCCS) dataset is a crucial step, especially when investigating intrinsically disordered proteins (IDPs). The scientific literature is abundant in such datasets, however, the effect of choosing one over all the others in a concrete application has not yet been studied thoroughly and systematically. Hereby, we review the available RCCS prediction methods and to compare them, we conduct statistical inference by means of the nonparametric sum of ranking differences and comparison of ranks to random numbers (SRD-CRRN) method. We try to find the RCCS predictors best representing the general consensus regarding secondary structural propensities. The existence and the magnitude of resulting differences on secondary structure determination under varying sample conditions (temperature, pH) are demonstrated and discussed for globular proteins and especially IDPs.

Molecular Dynamics and Structural Studies of Zinc Chloroquine Complexes

Chloroquine (CQ) is a first-choice drug against malaria and autoimmune diseases. It has been co-administered with zinc against SARS-CoV-2 and soon dismissed because of safety issues. The structural features of Zn-CQ complexes and the effect of CQ on zinc distribution in cells are poorly known. In this study, state-of-the-art computations combined with experiments were leveraged to solve the structural determinants of zinc-CQ interactions in solution and the solid state. NMR, ESI-MS, and X-ray absorption and diffraction methods were combined with ab initio molecular dynamics calculations to address the kinetic lability of this complex. Within the physiological pH range, CQ binds Zn²⁺ through the quinoline ring nitrogen, forming [Zn(CQH)Cl_x(H₂O)_3-x]^(3+)-x (x = 0, 1, 2, and 3) tetrahedral complexes. The Zn(CQH)Cl₃ species is stable at neutral pH and at high chloride concentrations typical of the extracellular medium, but metal coordination is lost at a moderately low pH as in the lysosomal lumen. The pentacoordinate complex [Zn(CQH)(H₂O)₄]³⁺ may exist in the absence of chloride. This in vitro/in silico approach can be extended to other metal-targeting drugs and bioinorganic systems.

Fabrication of Polysaccharide-Based Halochromic Nanofibers via Needle-Less Electrospinning and Their Characterization: A Study of the Leaching Effect

Responsive materials, i.e., smart materials, have the ability to change their physical or chemical properties upon certain external signals. The development of nanofibrous halochromic materials, specifically combining the pH-sensitive functionality and unique nanofiber properties, could yield interesting new applications, especially when the common problem of dye leaching is successfully tackled. Therefore, in this article, we studied the fabrication process of polysaccharide-based halochromic nanofibrous materials by using a combination of various halochromic dyes (bromothymol blue, bromocresol green, and thymol blue) and cellulose acetate in a spinning solution using a one-pot strategy. The inhibition of leaching was addressed by using a complexing agent: poly-diallyl-dimethylammonium chloride (PDADMAC). The preparation of hybrid spinning solutions, their characterization, and ability to form continuous nanofibers were studied using a high production needle-less electrospinning system. The produced hybrid solutions and nanofibers were characterized, in terms of their rheological properties, chemical structure, morphology, and functionality. Fabricated nanofibrous halochromic structures show a clear color change upon exposure to different pH values, as well as the reduced leaching of dyes, upon the addition of a complexing agent. The leaching decreased by 61% in the case of bromocresol green, while, in the case of bromothymol blue and thymol blue, the leaching was reduced by 95 and 99%, respectively.

Detecting aspartate isomerization and backbone cleavage after aspartate in intact proteins by NMR spectroscopy

The monitoring of non-enzymatic post-translational modifications (PTMs) in therapeutic proteins is important to ensure drug safety and efficacy. Together with methionine and asparagine, aspartic acid (Asp) is very sensitive to spontaneous alterations. In particular, Asp residues can undergo isomerization and peptide-bond hydrolysis, especially when embedded in sequence motifs that are prone to succinimide formation or when followed by proline (Pro). As Asp and isoAsp have the same mass, and the Asp-Pro peptide-bond cleavage may lead to an unspecific mass difference of + 18 Da under native conditions or in the case of disulfide-bridged cleavage products, it is challenging to directly detect and characterize such modifications by mass spectrometry (MS). Here we propose a 2D NMR-based approach for the unambiguous identification of isoAsp and the products of Asp-Pro peptide-bond cleavage, namely N-terminal Pro and C-terminal Asp, and demonstrate its applicability to proteins including a therapeutic monoclonal antibody (mAb). To choose the ideal pH conditions under which the NMR signals of isoAsp and C-terminal Asp are distinct from other random coil signals, we determined the pK_a values of isoAsp and C-terminal Asp in short peptides. The characteristic ¹H-¹³C chemical shift correlations of isoAsp, N-terminal Pro and C-terminal Asp under standardized conditions were used to identify these PTMs in lysozyme and in the therapeutic mAb rituximab (MabThera) upon prolonged storage under acidic conditions (pH 4-5) and 40 °C. The results show that the application of our 2D NMR-based protocol is straightforward and allows detecting chemical changes of proteins that may be otherwise unnoticed with other analytical methods.

In vivo detection of teriflunomide-derived fluorine signal during neuroinflammation using fluorine MR spectroscopy

Background: Magnetic resonance imaging (MRI) is indispensable for diagnosing neurological conditions such as multiple sclerosis (MS). MRI also supports decisions regarding the choice of disease-modifying drugs (DMDs). Determining in vivo tissue concentrations of DMDs has the potential to become an essential clinical tool for therapeutic drug monitoring (TDM). The aim here was to examine the feasibility of fluorine-19 (¹⁹F) MR methods to detect the fluorinated DMD teriflunomide (TF) during normal and pathological conditions.

Methods: We used ¹⁹F MR spectroscopy to detect TF in the experimental autoimmune encephalomyelitis (EAE) mouse model of multiple sclerosis (MS) in vivo. Prior to the in vivo investigations we characterized the MR properties of TF in vitro. We studied the impact of pH and protein binding as well as MR contrast agents.

Results: We could detect TF in vivo and could follow the ¹⁹F MR signal over different time points of disease. We quantified TF concentrations in different tissues using HPLC/MS and showed a significant correlation between ex vivo TF levels in serum and the ex vivo¹⁹F MR signal.

Conclusion: This study demonstrates the feasibility of ¹⁹F MR methods to detect TF during neuroinflammation in vivo. It also highlights the need for further technological developments in this field. The ultimate goal is to add ¹⁹F MR protocols to conventional ¹H MRI protocols in clinical practice to guide therapy decisions.

In vitro singlet state and zero-quantum encoded magnetic resonance spectroscopy: Illustration with N-acetyl-aspartate

Magnetic resonance spectroscopy (MRS) allows the analysis of biochemical processes non-invasively and in vivo. Still, its application in clinical diagnostics is rare. Routine MRS is limited to spatial, chemical and temporal resolutions of cubic centimetres, mM and minutes. In fact, the signal of many metabolites is strong enough for detection, but the resonances significantly overlap, exacerbating identification and quantification. Besides, the signals of water and lipids are much stronger and dominate the entire spectrum. To suppress the background and isolate selected signals, usually, relaxation times, J-coupling and chemical shifts are used. Here, we propose methods to isolate the signals of selected molecular groups within endogenous metabolites by using long-lived spin states (LLS). We exemplify the method by preparing the LLSs of coupled protons in the endogenous molecules N-acetyl-L-aspartic acid (NAA). First, we store polarization in long-lived, double spin states, followed by saturation pulses before the spin order is converted back to observable magnetization or double quantum filters to suppress background signals. We show that LLS and zero-quantum coherences can be used to selectively prepare and measure the signals of chosen metabolites or drugs in the presence of water, inhomogeneous field and highly concentrated fatty solutions. The strong suppression of unwanted signals achieved allowed us to measure pH as a function of chemical shift difference.

PKAD: a database of experimentally measured pKa values of ionizable groups in proteins

Ionizable residues play key roles in many biological phenomena including protein folding, enzyme catalysis and binding. We present PKAD, a database of experimentally measured pKas of protein residues reported in the literature or taken from existing databases. The database contains pKa data for 1350 residues in 157 wild-type proteins and for 232 residues in 45 mutant proteins. Most of these values are for Asp, Glu, His and Lys amino acids. The database is available as downloadable file as well as a web server (http://compbio.clemson.edu/pkad). The PKAD database can be used as a benchmarking source for development and improvement of pKa's prediction methods. The web server provides additional information taken from the corresponding structures and amino acid sequences, which allows for easy search and grouping of the experimental pKas according to various biophysical characteristics, amino acid type and others.

Modeling pH-Dependent NMR Chemical Shift Perturbations in Peptides

Modeling the pH dependence of protein and peptide chemical shifts outside the range of physiological values (6.5-7) is key to understanding structure-function relationships of these systems. These capabilities are largely not available in current chemical shift prediction software. In this study, we utilize a combination of molecular dynamics and quantum mechanics to investigate the through-space and through-bond contributions of protonation-dependent chemical shift perturbations (CSPs) in model tripeptides. By altering the protonation state of the titratable group in the tripeptides, we observe a notable difference in the conformational ensembles and attendantly compute significant CSPs for all nuclei near the site of protonation. We thus demonstrate the ability to recapitulate experimental pH-dependent CSPs with good agreement (R = 0.85, 0.99, and 0.98 for ¹³C, ¹⁵N, and ¹H, respectively). Broadly, we provide the groundwork for incorporating pH effects into empirical and semiempirical chemical shift predictors.

Helix-loop-helix peptide foldamers and their use in the construction of hydrolase mimetics

De novo designed helix-loop-helix peptide foldamers containing cis-2-aminocyclopentanecarboxylic acid residues were evaluated for their conformational stability and possible use in enzyme mimetic development. The correlation between hydrogen bond network size and conformational stability was demonstrated through CD and NMR spectroscopies. Molecules incorporating a Cys/His/Glu triad exhibited enzyme-like hydrolytic activity.

A Targeted Mutation Identified through pKa Measurements Indicates a Postrecruitment Role for Fis1 in Yeast Mitochondrial Fission

The tail-anchored protein Fis1 is implicated as a passive tether in yeast mitochondrial fission. We probed the functional role of Fis1 Glu-78, whose elevated side chain pKa suggests participation in protein interactions. Fis1 binds partners Mdv1 or Dnm1 tightly, but mutation E78A weakens Fis1 interaction with Mdv1, alters mitochondrial morphology, and abolishes fission in a growth assay. In fis1Δ rescue experiments, Fis1-E78A causes a novel localization pattern in which Dnm1 uniformly coats the mitochondria. By contrast, Fis1-E78A at lower expression levels recruits Dnm1 into mitochondrial punctate structures but fails to support normal fission. Thus, Fis1 makes multiple interactions that support Dnm1 puncta formation and may be essential after this step, supporting a revised model for assembly of the mitochondrial fission machinery. The insights gained by mutating a residue with a perturbed pKa suggest that side chain pKa values inferred from routine NMR sample pH optimization could provide useful leads for functional investigations.

A Soluble, Folded Protein without Charged Amino Acid Residues

Charges are considered an integral part of protein structure and function, enhancing solubility and providing specificity in molecular interactions. We wished to investigate whether charged amino acids are indeed required for protein biogenesis and whether a protein completely free of titratable side chains can maintain solubility, stability, and function. As a model, we used a cellulose-binding domain from Cellulomonas fimi, which, among proteins of more than 100 amino acids, presently is the least charged in the Protein Data Bank, with a total of only four titratable residues. We find that the protein shows a surprising resilience toward extremes of pH, demonstrating stability and function (cellulose binding) in the pH range from 2 to 11. To ask whether the four charged residues present were required for these properties of this protein, we altered them to nontitratable ones. Remarkably, this chargeless protein is produced reasonably well in Escherichia coli, retains its stable three-dimensional structure, and is still capable of strong cellulose binding. To further deprive this protein of charges, we removed the N-terminal charge by acetylation and studied the protein at pH 2, where the C-terminus is effectively protonated. Under these conditions, the protein retains its function and proved to be both soluble and have a reversible folding-unfolding transition. To the best of our knowledge, this is the first time a soluble, functional protein with no titratable side chains has been produced.

Role of Deprotonation Free Energies in pKa Prediction and Molecule Ranking

A computationally efficient classical molecular simulation technique is derived for ranking the pKa values of a set of chemically similar congeneric molecules in an implicit solvent model of water. This uses the deprotonation free energy of the titratable group in the gas and aqueous phases obtained by thermodynamic integration (TI). For a series of alcohols and acids a strong linear correlation is demonstrated between the experimental pKa and the deprotonation free energy difference in the gas and liquid phases. These calculations also show that classical TI is more efficient than slow-growth TI in calculating deprotonation free energies for the series of molecules considered herein.

Asymmetric protonation of EmrE

The small multidrug resistance transporter EmrE is a homodimer that uses energy provided by the proton motive force to drive the efflux of drug substrates. The pKa values of its "active-site" residues--glutamate 14 (Glu14) from each subunit--must be poised around physiological pH values to efficiently couple proton import to drug export in vivo. To assess the protonation of EmrE, pH titrations were conducted with (1)H-(15)N TROSY-HSQC nuclear magnetic resonance (NMR) spectra. Analysis of these spectra indicates that the Glu14 residues have asymmetric pKa values of 7.0 ± 0.1 and 8.2 ± 0.3 at 45°C and 6.8 ± 0.1 and 8.5 ± 0.2 at 25°C. These pKa values are substantially increased compared with typical pKa values for solvent-exposed glutamates but are within the range of published Glu14 pKa values inferred from the pH dependence of substrate binding and transport assays. The active-site mutant, E14D-EmrE, has pKa values below the physiological pH range, consistent with its impaired transport activity. The NMR spectra demonstrate that the protonation states of the active-site Glu14 residues determine both the global structure and the rate of conformational exchange between inward- and outward-facing EmrE. Thus, the pKa values of the asymmetric active-site Glu14 residues are key for proper coupling of proton import to multidrug efflux. However, the results raise new questions regarding the coupling mechanism because they show that EmrE exists in a mixture of protonation states near neutral pH and can interconvert between inward- and outward-facing forms in multiple different protonation states.

The pKa value and accessibility of cysteine residues are key determinants for protein substrate discrimination by glutaredoxin

The enzyme glutaredoxin catalyzes glutathione exchange, but little is known about its interaction with protein substrates. Very different proteins are substrates in vitro, and the enzyme seems to have low requirements for specific protein interactions. Here we present a systematic investigation of the interaction between human glutaredoxin 1 and glutathionylated variants of a single model protein. Thus, single cysteine variants of acyl-coenzyme A binding protein were produced creating a set of substrates in the same protein background. The rate constants for deglutathionylation differ by more than 2 orders of magnitude between the best (k1 = 1.75 × 10(5) M(-1) s(-1)) and the worst substrate (k1 = 4 × 10(2) M(-1) s(-1)). The pKa values of the substrate cysteine residues were determined by NMR spectroscopy and found to vary from 8.2 to 9.9. Rates of glutaredoxin 1-catalyzed deglutathionylation were assessed with respect to substrate cysteine pKa values, cysteine residue accessibility, local stability, and backbone dynamics. Good substrates are characterized by a combination of high accessibility of the glutathionylated site and low pKa of the cysteine residue.

Protein dielectric constants determined from NMR chemical shift perturbations

Understanding the connection between protein structure and function requires a quantitative understanding of electrostatic effects. Structure-based electrostatic calculations are essential for this purpose, but their use has been limited by a long-standing discussion on which value to use for the dielectric constants (ε(eff) and ε(p)) required in Coulombic and Poisson-Boltzmann models. The currently used values for ε(eff) and ε(p) are essentially empirical parameters calibrated against thermodynamic properties that are indirect measurements of protein electric fields. We determine optimal values for ε(eff) and ε(p) by measuring protein electric fields in solution using direct detection of NMR chemical shift perturbations (CSPs). We measured CSPs in 14 proteins to get a broad and general characterization of electric fields. Coulomb's law reproduces the measured CSPs optimally with a protein dielectric constant (ε(eff)) from 3 to 13, with an optimal value across all proteins of 6.5. However, when the water-protein interface is treated with finite difference Poisson-Boltzmann calculations, the optimal protein dielectric constant (ε(p)) ranged from 2 to 5 with an optimum of 3. It is striking how similar this value is to the dielectric constant of 2-4 measured for protein powders and how different it is from the ε(p) of 6-20 used in models based on the Poisson-Boltzmann equation when calculating thermodynamic parameters. Because the value of ε(p) = 3 is obtained by analysis of NMR chemical shift perturbations instead of thermodynamic parameters such as pK(a) values, it is likely to describe only the electric field and thus represent a more general, intrinsic, and transferable ε(p) common to most folded proteins.

Mass spectrometry and NMR analysis of ligand binding by human liver fatty acid binding protein

Human liver fatty acid binding protein (hL-FABP) is the most abundant cytosolic protein in the liver. This protein plays important roles associated to partitioning of fatty acids (FAs) to specific metabolic pathways, nuclear signaling and protection against oxidative damage. The protein displays promiscuous binding properties and can bind two internal ligands, unlike FABPs from other tissues. Different topologies for the ligand located in the more accessible site have been reported, with either a 'head-in' or 'head-out' orientation of the carboxylate end. Electrospray-ionization mass spectrometry and nuclear magnetic resonance titrations are employed here in order to investigate in further detail the binding properties of this system, the equilibria established in solution and the pH dependence of the complexes. The results are consistent with two binding sites with different affinity and a unique head-out topology for the second molecule of either ligand. Competition experiments indicate a higher affinity for oleic acid relative to palmitic acid at each binding site.

Highly perturbed pKa values in the unfolded state of hen egg white lysozyme

The majority of pK(a) values in protein unfolded states are close to the amino acid model pK(a) values, thus reflecting the weak intramolecular interactions present in the unfolded ensemble of most proteins. We have carried out thermal denaturation measurements on the WT and eight mutants of HEWL from pH 1.5 to pH 11.0 to examine the unfolded state pK(a) values and the pH dependence of protein stability for this enzyme. The availability of accurate pK(a) values for the folded state of HEWL and separate measurements of mutant-induced effects on the folded state pK(a) values, allows us to estimate the pK(a) values of seven acidic residues in the unfolded state of HEWL. Asp-48 and Asp-66 display pK(a) values of 2.9 and 3.1 in our analysis, thus representing the most depressed unfolded state pK(a) values observed to date. We observe a strong correlation between the folded state pK(a) values and the unfolded state pK(a) values of HEWL, thus suggesting that the unfolded state of HEWL possesses a large degree of native state characteristics.

Elucidation of isoform-dependent pH sensitivity of troponin i by NMR spectroscopy

Myocardial ischemia is characterized by reduced blood flow to cardiomyocytes, which can lead to acidosis. Acidosis decreases the calcium sensitivity and contractile efficiency of cardiac muscle. By contrast, skeletal and neonatal muscles are much less sensitive to changes in pH. The pH sensitivity of cardiac muscle can be reduced by replacing cardiac troponin I with its skeletal or neonatal counterparts. The isoform-specific response of troponin I is dictated by a single histidine, which is replaced by an alanine in cardiac troponin I. The decreased pH sensitivity may stem from the protonation of this histidine at low pH, which would promote the formation of electrostatic interactions with negatively charged residues on troponin C. In this study, we measured acid dissociation constants of glutamate residues on troponin C and of histidine on skeletal troponin I (His-130). The results indicate that Glu-19 comes in close contact with an ionizable group that has a pK(a) of ∼6.7 when it is in complex with skeletal troponin I but not when it is bound to cardiac troponin I. The pK(a) of Glu-19 is decreased when troponin C is bound to skeletal troponin I and the pK(a) of His-130 is shifted upward. These results strongly suggest that these residues form an electrostatic interaction. Furthermore, we found that skeletal troponin I bound to troponin C tighter at pH 6.1 than at pH 7.5. The data presented here provide insights into the molecular mechanism for the pH sensitivity of different muscle types.

Selective observation of Asp and Glu resonances in (13)CO detected experiments

We present a filter element to observe exclusively the resonances of Asp and Glu residues in several (13)C detected experiments, CACO, CBCACO, and CCCO. As in these experiments the carboxylate resonances appear in the directly detected dimension, it is possible to measure their chemical shifts with a high precision. Therefore, the experiments can be of great utility in the accurate determination of the pKa 's of these residues. Furthermore, the experiments can be applied in the study of deuterated proteins, where the usual experiments for pKa determination cannot be used. Finally, the good resolution with which the carboxylate spectrum is obtained allows observing the coupling between the carboxyl carbon and the backbone CO in Asp residues, what provides information on their side chain conformation.

Dissecting electrostatic interactions in Bacillus circulans xylanase through NMR-monitored pH titrations

NMR-monitored pH titration curves of proteins provide a rich source of structural and electrostatic information. Although relatively straightforward to measure, interpreting pH-dependent chemical shift changes to obtain site-specific acid dissociation constants (pK (A) values) is challenging. In order to analyze the biphasic titrations exhibited by the side chain (13)C(γ) nuclei of the nucleophilic Glu78 and general acid/base Glu172 in Bacillus circulans xylanase, we have revisited the formalism for the ionization equilibria of two coupled acidic residues. In general, fitting NMR-monitored pH titration curves for such a system will only yield the two macroscopic pK (A) values that reflect the combined effects of both deprotonation reactions. However, through the use of mutations complemented with ionic strength-dependent measurements, we are able to extract the four microscopic pK (Ai) values governing the branched acid/base equilibria of Glu78 and Glu172 in BcX. These data, confirmed through theoretical calculations, help explain the pH-dependent mechanism of this model GH11 xylanase by demonstrating that the kinetically determined pK (A) values and hence catalytic roles of these two residues result from their electrostatic coupling.

On the development of protein pKa calculation algorithms

Protein pK(a) calculation methods are developed partly to provide fast non-experimental estimates of the ionization constants of protein side chains. However, the most significant reason for developing such methods is that a good pK(a) calculation method is presumed to provide an accurate physical model of protein electrostatics, which can be applied in methods for drug design, protein design, and other structure-based energy calculation methods. We explore the validity of this presumption by simulating the development of a pK(a) calculation method using artificial experimental data derived from a human-defined physical reality. We examine the ability of an RMSD-guided development protocol to retrieve the correct (artificial) physical reality and find that a rugged optimization landscape and a huge parameter space prevent the identification of the correct physical reality. We examine the importance of the training set in developing pK(a) calculation methods and investigate the effect of experimental noise on our ability to identify the correct physical reality, and find that both effects have a significant and detrimental impact on the physical reality of the optimal model identified. Our findings are of relevance to all structure-based methods for protein energy calculations and simulation, and have large implications for all types of current pK(a) calculation methods. Our analysis furthermore suggests that careful and extensive validation on many types of experimental data can go some way in making current models more realistic.

Electrostatics in proteins and protein-ligand complexes

Accurate computational methods for predicting electrostatic energies are of major importance for our understanding of protein energetics in general for computer-aided drug design as well as for the design of novel biocatalysts and protein therapeutics. Electrostatic energies are of particular importance in such applications as virtual screening, drug design and protein-protein docking due to the high charge density of protein ligands and small-molecule drugs, and the frequent protonation state changes observed when drugs bind to their protein targets. Therefore, the development of a reliable and fast algorithm for the evaluation of electrostatic free energies, as an important contributor to the overall protein energy function, has been the focus for many scientists over the past three decades. In this review we describe the current state-of-the-art in modeling electrostatic effects in proteins and protein-ligand complexes. We focus mainly on the merits and drawbacks of the continuum methodology, and speculate on future directions in refining algorithms for calculating electrostatic energies in proteins using experimental data.

Calculating pKa values in the cAMP-dependent protein kinase: the effect of conformational change and ligand binding

The conformational change observed upon ligand binding and phosphorylation for the cAMP-dependent protein kinase (protein kinase A-PKA) is of high importance for the regulation of its activity. We calculate pKa values and net charges for 18 3D structures of PKA in various conformations and liganded states to examine the role of electrostatics in ligand binding and activation. We find that the conformational change of PKA takes place without any significant net proton uptake/release at all pH values, thus indicating that PKA has evolved to reduce any pH-dependent barriers to the conformational motion. We furthermore find that the binding of ligands induces large changes in the net charge of PKA at most pH values, but significantly, we find that the net charge difference at physiological pH is close to zero, thus indicating that the active-site pKa values have been preorganized for substrate binding. We are unable to unequivocally resolve the identity of the groups responsible for determining the pH-activity profile of PKA but speculate that the titration of Lys 168 or the titration of ATP itself could be responsible for the loss of activity at high pH values. Finally, we examine the effect of point mutations on the pKa values of the PKA catalytic residues and find these to be relatively insensitive to both noncharge-altering and charge-altering mutations.

The molecular basis of distinct aggregation pathways of islet amyloid polypeptide

Abnormal aggregation of islet amyloid polypeptide (IAPP) into amyloid fibrils is a hallmark of type 2 diabetes. In this study, we investigated the initial oligomerization and subsequent addition of monomers to growing aggregates of human IAPP at the residue-specific level using NMR, atomic force microscopy, mass spectroscopy, and computational simulations. We found that in solution IAPPs rapidly associate into transient low-order oligomers such as dimers and trimers via interactions between histidine 18 and tyrosine 37. This initial event is proceeded by slow aggregation into higher-order spherical oligomers and elongated fibrils. In these two morphologically distinct types of aggregates IAPPs adopt structures with markedly different residual flexibility. Here we show that the anti-amyloidogenic compound resveratrol inhibits oligomerization and amyloid formation via binding to histidine 18, supporting the finding that this residue is crucial for on-pathway oligomer formation.

Remeasuring HEWL pK(a) values by NMR spectroscopy: methods, analysis, accuracy, and implications for theoretical pK(a) calculations

Site-specific pK(a) values measured by NMR spectroscopy provide essential information on protein electrostatics, the pH-dependence of protein structure, dynamics and function, and constitute an important benchmark for protein pK(a) calculation algorithms. Titration curves can be measured by tracking the NMR chemical shifts of several reporter nuclei versus sample pH. However, careful analysis of these curves is needed to extract residue-specific pK(a) values since pH-dependent chemical shift changes can arise from many sources, including through-bond inductive effects, through-space electric field effects, and conformational changes. We have re-measured titration curves for all carboxylates and His 15 in Hen Egg White Lysozyme (HEWL) by recording the pH-dependent chemical shifts of all backbone amide nitrogens and protons, Asp/Glu side chain protons and carboxyl carbons, and imidazole protonated carbons and protons in this protein. We extracted pK(a) values from the resulting titration curves using standard fitting methods, and compared these values to each other, and with those measured previously by ¹H NMR (Bartik et al., Biophys J 1994;66:1180–1184). This analysis gives insights into the true accuracy associated with experimentally measured pK(a) values. We find that apparent pK(a) values frequently differ by 0.5–1.0 units depending upon the nuclei monitored, and that larger differences occasionally can be observed. The variation in measured pK(a) values, which reflects the difficulty in fitting and assigning pH-dependent chemical shifts to specific ionization equilibria, has significant implications for the experimental procedures used for measuring protein pK(a) values, for the benchmarking of protein pK(a) calculation algorithms, and for the understanding of protein electrostatics in general.

Variation and analysis of second-sphere interactions and axial histidinate character in c-type cytochromes

The electron-donating properties of the axial His ligand to heme iron in cytochromes c (cyts c) are found to be correlated with the midpoint reduction potential (E(m)) in variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht cyt c(552)) in which mutations have been made in and near the Cys-X-X-Cys-His (CXXCH) heme-binding motif. To probe the strength of the His-Fe(III) interaction, we have measured (13)C nuclear magnetic resonance (NMR) chemical shifts for (13)CN(-) bound to heme iron trans to the axial His in Ht Fe(III) cyt c(552) variants. We observe a linear relationship between these (13)C chemical shifts and E(m), indicating that the His-Fe(III) bond strength correlates with E(m). To probe a conserved hydrogen bonding interaction between the axial His Hdelta1 and the backbone carbonyl of a Pro residue, we measured the pK(a) of the axial His Hdelta1 proton (pK(a(2))), which we propose to relate to the His-Fe(III) interaction, reduction potential, and local electrostatic effects. The observed linear relationship between the axial His (13)Cbeta chemical shift and E(m) is proposed to reflect histidinate (anionic) character of the ligand. A linear relationship also is seen between the average heme methyl (1)H chemical shift and E(m) which may reflect variation in axial His electron-donating properties or in the ruffling distortion of the heme plane. In summary, chemical shifts of the axial His and exogenous CN(-) bound trans to His are shown to be sensitive probes of the His-Fe(III) interaction in variants of Ht cyt c(552) and display trends that correlate with E(m).

Capturing, sharing and analysing biophysical data from protein engineering and protein characterization studies

Large amounts of data are being generated annually on the connection between the sequence, structure and function of proteins using site-directed mutagenesis, protein design and directed evolution techniques. These data provide the fundamental building blocks for our understanding of protein function, molecular biology and living organisms in general. However, much experimental data are never deposited in databases and is thus 'lost' in journal publications or in PhD theses. At the same time theoretical scientists are in need of large amounts of experimental data for benchmarking and calibrating novel predictive algorithms, and theoretical progress is therefore often hampered by the lack of suitable data to validate or disprove a theoretical assumption. We present PEAT (Protein Engineering Analysis Tool), an application that integrates data deposition, storage and analysis for researchers carrying out protein engineering projects or biophysical characterization of proteins. PEAT contains modules for DNA sequence manipulation, primer design, fitting of biophysical characterization data (enzyme kinetics, circular dichroism spectroscopy, NMR titration data, etc.), and facilitates sharing of experimental data and analyses for a typical university-based research group. PEAT is freely available to academic researchers at http://enzyme.ucd.ie/PEAT.