Simplified methods for pKa and acid pH-dependent stability estimation in proteins: removing dielectric and counterion boundaries

Much computational research aimed at understanding ionizable group interactions in proteins has focused on numerical solutions of the Poisson-Boltzmann (PB) equation, incorporating protein exclusion zones for solvent and counterions in a continuum model. Poor agreement with measured pKas and pH-dependent stabilities for a (protein, solvent) relative dielectric boundary of (4,80) has lead to the adoption of an intermediate (20,80) boundary. It is now shown that a simple Debye-Huckel (DH) calculation, removing both the low dielectric and counterion exclusion regions associated with protein, is equally effective in general pKa calculations. However, a broad-based discrepancy to measured pH-dependent stabilities is maintained in the absence of ionizable group interactions in the unfolded state. A simple model is introduced for these interactions, with a significantly improved match to experiment that suggests a potential utility in predicting and analyzing the acid pH-dependence of protein stability. The methods are applied to the relative pH-dependent stabilities of the pore-forming domains of colicins A and N. The results relate generally to the well-known preponderance of surface ionizable groups with solvent-mediated interactions. Although numerical PB solutions do not currently have a significant advantage for overall pKa estimations, development based on consideration of microscopic solvation energetics in tandem with the continuum model could combine the large deltapKas of a subset of ionizable groups with the overall robustness of the DH model.     

Calculating pH-dependent free energy of proteins by using Monte Carlo protonation probabilities of ionizable residues

Protein folding, stability, and function are usually influenced by pH. And free energy plays a fundamental role in analysis of such pH-dependent properties. Electrostatics-based theoretical framework using dielectric solvent continuum model and solving Poisson-Boltzmann equation numerically has been shown to be very successful in understanding the pH-dependent properties. However, in this approach the exact computation of pH-dependent free energy becomes impractical for proteins possessing more than several tens of ionizable sites (e.g.>30), because exact evaluation of the partition function requires a summation over a vast number of possible protonation microstates. Here we present a method which computes the free energy using the average energy and the protonation probabilities of ionizable sites obtained by the well-established Monte Carlo sampling procedure. The key feature is to calculate the entropy by using the protonation probabilities. We used this method to examine a well-studied protein (lysozyme) and produced results which agree very well with the exact calculations. Applications to the optimum pH of maximal stability of proteins and protein–DNA interactions have also resulted in good agreement with experimental data. These examples recommend our method for application to the elucidation of the pH-dependent properties of proteins.     

The ATP-dependent Clp protease in chloroplasts of higher plants

The best-known proteases in plastids are those that belong to families common to eubacteria. One of the first identified was the ATP-dependent caseinolytic protease (Clp), whose structure and function have been well characterized in Escherichia coli . Plastid Clp proteins in higher plants are surprisingly numerous and diverse, with at least 16 distinct Clp proteins in the model plant Arabidopsis thaliana . Multiple paralogues exist for several of the different types of plastid Clp protein, with the most extreme being five for the proteolytic subunit ClpP. Both biochemical and genetic studies have recently begun to reveal the intricate structural interactions between the various Clp proteins, and their importance for chloroplast function and plant development. Much of the recent data suggests that the function of many of the Clp proteins probably affects more specific processes within chloroplasts, in addition to the more general 'housekeeping' role previously assumed.     

pH-Dependent Interactions Guide the Folding and Gate the Transmembrane Pore of the β-Barrel Membrane Protein OmpG

The physical interactions that switch the functional state of membrane proteins are poorly understood. Previously, the pH-gating conformations of the β-barrel forming outer membrane protein G (OmpG) from Escherichia coli have been solved. When the pH changes from neutral to acidic the flexible extracellular loop L6 folds into and closes the OmpG pore. Here, we used single-molecule force spectroscopy to structurally localize and quantify the interactions that are associated with the pH-dependent closure. At acidic pH, we detected a pH-dependent interaction at loop L6. This interaction changed the (un)folding of loop L6 and of β-strands 11 and 12, which connect loop L6. All other interactions detected within OmpG were unaffected by changes in pH. These results provide a quantitative and mechanistic explanation of how pH-dependent interactions change the folding of a peptide loop to gate the transmembrane pore. They further demonstrate how the stability of OmpG is optimized so that pH changes modify only those interactions necessary to gate the transmembrane pore.     

Structure and mechanism of the saposin-like domain of a plant aspartic protease

Many plant aspartic proteases contain an additional sequence of ∼100 amino acids termed the plant-specific insert, which is involved in host defense and vacuolar targeting. Similar to all saposin-like proteins, the plant-specific insert functions via protein-membrane interactions; however, the structural basis for such interactions has not been studied, and the nature of plant-specific insert-mediated membrane disruption has not been characterized. In the present study, the crystal structure of the saposin-like domain of potato aspartic protease was resolved at a resolution of 1.9 Å, revealing an open V-shaped configuration similar to the open structure of human saposin C. Notably, vesicle disruption activity followed Michaelis-Menten-like kinetics, a finding not previously reported for saposin-like proteins including plant-specific inserts. Circular dichroism data suggested that secondary structure was pH-dependent in a fashion similar to influenza A hemagglutinin fusion peptide. Membrane effects characterized by atomic force microscopy and light scattering indicated bilayer solubilization as well as fusogenic activity. Taken together, the present study is the first report to elucidate the membrane interaction mechanism of plant saposin-like domains whereby pH-dependent membrane interactions resulted in bilayer fusogenic activity that probably arose from a viral type pH-dependent helix-kink-helix motif at the plant-specific insert N terminus.     

A fast and accurate computational approach to protein ionization

We report a very fast and accurate physics-based method to calculate pH-dependent electrostatic effects in protein molecules and to predict the pK values of individual sites of titration. In addition, a CHARMm-based algorithm is included to construct and refine the spatial coordinates of all hydrogen atoms at a given pH. The present method combines electrostatic energy calculations based on the Generalized Born approximation with an iterative mobile clustering approach to calculate the equilibria of proton binding to multiple titration sites in protein molecules. The use of the GBIM (Generalized Born with Implicit Membrane) CHARMm module makes it possible to model not only water-soluble proteins but membrane proteins as well. The method includes a novel algorithm for preliminary refinement of hydrogen coordinates. Another difference from existing approaches is that, instead of monopeptides, a set of relaxed pentapeptide structures are used as model compounds. Tests on a set of 24 proteins demonstrate the high accuracy of the method. On average, the RMSD between predicted and experimental pK values is close to 0.5 pK units on this data set, and the accuracy is achieved at very low computational cost. The pH-dependent assignment of hydrogen atoms also shows very good agreement with protonation states and hydrogen-bond network observed in neutron-diffraction structures. The method is implemented as a computational protocol in Accelrys Discovery Studio and provides a fast and easy way to study the effect of pH on many important mechanisms such as enzyme catalysis, ligand binding, protein-protein interactions, and protein stability.     

On the pH regulation of plant aquaporins

The majority of plants are unable to evade unfavorable conditions such as flooding, salinity, or drought. Therefore, a fine-tuned water homeostasis appears to be of crucial importance for plant survival, and it was assumed that aquaporins play a significant role in these processes. Regulation of plant aquaporin conductivity was suggested to be achieved by a gating mechanism that involves protein phosphorylation under drought stress conditions and protonation after cytosolic acidification during flooding. The effect of protein phosphorylation or protonation of aquaporins was studied on two plasma membrane intrinsic proteins, NtPIP2;1 and NtAQP1 from tobacco, which were heterologously expressed in yeast. Our results on mutated aquaporins with serine-to-alanine exchange indicate that phosphorylation of the two key serine residues did not affect the pH-dependent modification of water permeability. Protonation on a conserved histidine residue decreased water conductivity of NtPIP2;1. Although cells expressing NtPIP2;1 with a replacement of the histidine by an alanine were found to be pH-insensitive with regard to water permeability, these maintain high water transport rates, similar to those obtained under acidic conditions. The data clearly support the role of histidine at 196 as a component of pH-dependent modification of aquaporin-facilitated water transport. The predictions of combined effects from phosphorylation at conserved serines and histidine protonation were not supported by the results of functional analysis. The obtained results challenge the gating model as a general regulation mechanism for plant plasma membrane aquaporins.     

Recent development and application of constant pH molecular dynamics

Solution pH is a critical environmental factor for chemical and biological processes. Over the last decade, significant efforts have been made in the development of constant pH molecular dynamics (pHMD) techniques for gaining detailed insights into pH-coupled dynamical phenomena. In this article, we review the advancement of this field in the past five years, placing a special emphasis on the development of the all-atom continuous pHMD technique. We discuss various applications, including the prediction of pK_a shifts for proteins, nucleic acids and surfactant assemblies, elucidation of pH-dependent population shifts, protein–protein and protein–RNA binding, as well as the mechanisms of pH-dependent self-assembly and phase transitions of surfactants and peptides. We also discuss future directions for the further improvement of the pHMD techniques.     

Adenovirus complex structures

Adenovirus has, for a long time, been a model system for understanding complex virus structure, assembly and interference in host cell processes. Recent structures of adenoviral capsid proteins critical for cell entry have given new insights into both interactions with host cell receptors and inter-capsid protein interactions, which determine the capsid architecture. Such studies are of importance in engineering adenovirus for use in various gene transfer applications. Remarkable and unexpected similarities have been revealed between the cell-attachment proteins and primary receptors of adenovirus and the unrelated reovirus, and between the capsid proteins and architecture of adenovirus, the enveloped bacteriophage PRD1 and other large DNA viruses.     

An integrated model of thermodynamic-hemodynamic-pharmacokinetic system and its application on decoupling control of intracranial temperature and pressure in brain hypothermia treatment

Brain hypothermia treatment (BHT) is an intensive care characterized by simultaneous managements of various vital signs, such as intracranial temperature (ICT) and pressure (ICP), of the severe neuropatient. Medical treatments including therapeutic ambient cooling and diuresis are separately carried out based on the experience of the medical staff involved in the clinical management of various pathophysiological processes, such as thermodynamics, hemodynamics and pharmacokinetics. However, no special attention has been paid to the interactions among these subsystems in therapeutic hypothermia because of the lack of theoretical knowledge. Therefore, quantitative analyses using an integrated model of various physiological processes and their interactions are of pressing need. In the present paper, we propose a general compartmental model to describe the pathophysiological processes of the three aforementioned dynamics, on account of the dynamical analogy of temperature, pressure and concentration. The model is verified by the agreement of model-based simulation results with clinical evidence. Based on responses of the integrated model to various stimuli, a transfer function matrix is identified to linearly approximate the characteristic interrelationships between medical treatments (ambient cooling and diuresis) and the vital signs (ICT and ICP). Then a controller that decouples ambient cooling and diuresis is proposed for efficient management of ICT and ICP, enhancement of hypothermic decompression and reduction of diuretic dosage. Decoupling control simulation indicates that ICT and ICP of the integrated model, representing a patient under BHT, can be simultaneously regulated by a single PID controller for ambient cooling and another for diuresis. The proposed decoupler effectively establishes hypothermic decompression, reduces the dosage of diuretic and improves ICP management. Theoretical analyses of the integrated model and decoupling control of ICT and ICP provide insights into the intensive care of various pathophysiological processes in patients undergoing BHT.     

Ionization, partitioning, and dynamics of tryptophan octyl ester: implications for membrane-bound tryptophan residues.

The presence of tryptophan residues as intrinsic fluorophores in most proteins makes them an obvious choice for fluorescence spectroscopic analyses of such proteins. Membrane proteins have been reported to have a significantly higher tryptophan content than soluble proteins. The role of tryptophan residues in the structure and function of membrane proteins has attracted a lot of attention. Tryptophan residues in membrane proteins and peptides are believed to be distributed asymmetrically toward the interfacial region. Tryptophan octyl ester (TOE) is an important model for membrane-bound tryptophan residues. We have characterized this molecule as a fluorescent membrane probe in terms of its ionization, partitioning, and motional characteristics in unilamellar vesicles of dioleoylphosphatidylcholine. The ionization property of this molecule in model membranes has been studied by utilizing its pH-dependent fluorescence characteristics. Analysis of pH-dependent fluorescence intensity and emission maximum shows that deprotonation of the alpha-amino group of TOE occurs with an apparent pKa of approximately 7.5 in the membrane. The fluorescence lifetime of membrane-bound TOE also shows pH dependence. The fluorescence lifetimes of TOE have been interpreted by using the rotamer model for the fluorescence decay of tryptophan. Membrane/water partition coefficients of TOE were measured in both its protonated and deprotonated forms. No appreciable difference was found in its partitioning behavior with ionization. Analysis of fluorescence polarization of TOE as a function of pH showed that there is a decrease in polarization with increasing pH, implying more rotational freedom on deprotonation. This is further supported by pH-dependent red edge excitation shift and the apparent rotational correlation time of membrane-bound TOE. TOE should prove useful in monitoring the organization and dynamics of tryptophan residues incorporated into membranes.     

Interactions between PEG and type I soluble tumor necrosis factor receptor: modulation by pH and by PEGylation at the N terminus

The effects of polyethylene glycol (PEG) on protein structure and the molecular details that regulate its association to polypeptides are largely unknown. These issues were addressed using type I soluble tumor necrosis factor receptor (sTNF-RI) as a model system. Changes in solution viscosity established that a truncated form of sTNF-RI bound free PEG in a pH-dependent manner. Above pH 5.3, the viscosity escalated as the pH increased, while no effect occurred below pH 5.0. Conjugation of 2 kD, 5 kD, or 20 kD PEG to the N terminus attenuated the viscosity at the higher pH values. Tryptophan phosphorescence spectroscopy correlated changes in the protein structure about Trp-107, at the C terminus, with the pH-dependent and PEGylation-dependent attenuation of the viscosity. The results indicate that specific interactions between PEG and the truncated form of sTNF-RI are elicited by an increased flexibility of the truncated protein combined perhaps with removal of steric or charge barriers. Covalently bound PEG at the N terminus reduced the protein affinity for the free polymer and induced a more rigid and polar configuration around Trp-107. Deprotonation of His-105, which is perpendicular to Trp-107, was integral to the binding mechanism producing a pH-dependent switching mechanism. These findings stress the importance of surface charge and structural plasticity in determining macromolecular binding affinities and demonstrate the ability of conjugated PEG to modify the localized surface structure in proteins away from the site of conjugation.     

Modeling protein recognition of carbohydrates

We have a limited understanding of the details of molecular recognition of carbohydrates by proteins, which is critical to a multitude of biological processes. Furthermore, carbohydrate-modifying proteins such as glycosyl hydrolases and phosphorylases are of growing importance as potential drug targets. Interactions between proteins and carbohydrates have complex thermodynamics, and in general the specific positioning of only a few hydroxyl groups determines their binding affinities. A thorough understanding of both carbohydrate and protein structures is thus essential to predict these interactions. An atomic-level view of carbohydrate recognition through structures of carbohydrate-active enzymes complexed with transition-state inhibitors reveals some of the distinctive molecular features unique to protein-carbohydrate complexes. However, the inherent flexibility of carbohydrates and their often water-mediated hydrogen bonding to proteins makes simulation of their complexes difficult. Nonetheless, recent developments such as the parameterization of specific force fields and docking scoring functions have greatly improved our ability to predict protein-carbohydrate interactions. We review protein-carbohydrate complexes having defined molecular requirements for specific carbohydrate recognition by proteins, providing an overview of the different computational techniques available to model them.     

Membrane insertion pathway of annexin B12: thermodynamic and kinetic characterization by fluorescence correlation spectroscopy and fluorescence quenching

Experimental determination of the free energy stabilizing the structure of membrane proteins in their native lipid environment is undermined by the lack of appropriate methods and suitable model systems. Annexin B12 (ANX) is a soluble protein which reversibly inserts into lipid membranes under mildly acidic conditions, which makes it a good experimental model for thermodynamic studies of folding and stability of membrane proteins. Here we apply fluorescence correlation spectroscopy for quantitative analysis of ANX partitioning into large unilamellar vesicles containing either 25% or 75% anionic lipids. Membrane binding of ANX results in changes of autocorrelation time and amplitude, both of which are used in quantitative analysis. The thermodynamic scheme describing acid-induced membrane interactions of ANX considers two independent processes: pH-dependent formation of a membrane-competent form near the membrane interface and its insertion into the lipid bilayer. Our novel fluorescence lifetime topology method demonstrates that the insertion proceeds via an interfacial refolded intermediate state, which can be stabilized by anionic lipids. Lipid titration measurements are used to determine the free energy of both transmembrane insertion and interfacial penetration, which are found to be similar, approximately -10-12 kcal/mol. The formation of the membrane-competent form, examined in a lipid saturation experiment, was found to depend on the local proton concentration near the membrane interface, occurring with pK = 4.3 and involving the protonation of two residues. Our results demonstrate that fluorescence correlation spectroscopy is a convenient tool for the quantitative characterization of the energetics of transmembrane insertion and that pH-triggered ANX insertion is a useful model for studying the thermodynamic stability of membrane proteins.     

Protein complementation as tool for studying protein-protein interactions in living cells

The association and dissociation of protein-protein complexes play an important role in various processes in living cells. The disruption of protein-protein interactions is observed in various pathologies. The study of the nature of these interactions will contribute to a better understanding of the molecular basis of the pathogenesis of the disease and the development of new approaches to therapy. Now there is a set of methods that allow one to reveal and analyze the interaction of proteins in vitro. However, more accurate data can be obtained by studying protein-protein interactions in vivo. One of a few prospective methods is based on the effect of the complementation of fragments of reporter proteins. These reporter systems are based on the change in the fluorescent properties or enzymatic activity of the proteins that can be measured using colorimetric, fluorescent, or other substrates. The principle of the complementation is widely used to analyze protein interactions, to determine of order of interaction of protein partners in different signaling pathways, as well as in high-performance screening studies for detecting and mapping previously unknown protein-protein interactions. The possibilities of existing complementation reporter systems allow one to solve problems that are far beyond the simple registration of the interactions of two or more proteins.     

pH-Dependent enzyme deactivation models

A pH-dependent "series-type" enzyme deactivation model using rapid protonation and deprotonation equilibria and the relatively slower inactivation rates is presented. From the enzyme activity-time trajectories at different pH the models presented permit the evaluation of some of the protonation and inactivation rate constants as well as the specific activities of the different enzyme forms. pH dependence of enzyme deactivations may also exhibit deactivation disguised kinetics. Three different examples of pH-dependent enzyme deactivations available in the literature are appropriately modeled to indicate the general applicability of the model. The model presented is consistent with the data and provides mechanistic insights into the pH-dependent deactivation of different enzymes.