The evolutionary origins and catalytic importance of conserved electrostatic networks within TIM-barrel proteins

Conservation of function is the basic tenet of protein evolution. Conservation of key electrostatic properties is a frequently employed mechanism that leads to conserved function. In a previous report, we identified several conserved electrostatic properties in four protein families and one functionally diverse enzyme superfamily. In this report, we demonstrate the evolutionary and catalytic importance of electrostatic networks in three ubiquitous metabolic enzymes: triosephosphate isomerase, enolase, and transaldolase. Evolutionary importance is demonstrated using phylogenetic motifs (sequence fragments that parallel the overall familial phylogeny). Phylogenetic motifs frequently correspond to both catalytic residues and conserved interactions that fine-tune catalytic residue pKa values. Further, in the case of triosephosphate isomerase, quantitative differences in the catalytic Glu169 pKa values parallel subfamily differentiation. Finally, phylogenetic motifs are shown to structurally cluster around the active sites of eight different TIM-barrel families. Depending upon the mechanistic requisites of each reaction catalyzed, interruptions to the canonical fold may or may not be identified as phylogenetic motifs.     

The role of surface electrostatics on the stability,function and regulation of human cystathionine β-synthase,a complex multidomain and oligomeric protein

Human cystathionine β-synthase (hCBS) is a key enzyme of sulfur amino acid metabolism, controlling the commitment of homocysteine to the transsulfuration pathway and antioxidant defense. Mutations in hCBS cause inherited homocystinuria (HCU), a rare inborn error of metabolism characterized by accumulation of toxic homocysteine in blood and urine. hCBS is a complex multidomain and oligomeric protein whose activity and stability are independently regulated by the binding of S-adenosyl-methionine (SAM) to two different types of sites at its C-terminal regulatory domain. Here we study the role of surface electrostatics on the complex regulation and stability of hCBS using biophysical and biochemical procedures. We show that the kinetic stability of the catalytic and regulatory domains is significantly affected by the modulation of surface electrostatics through noticeable structural and energetic changes along their denaturation pathways. We also show that surface electrostatics strongly affect SAM binding properties to those sites responsible for either enzyme activation or kinetic stabilization. Our results provide new insight into the regulation of hCBS activity and stability in vivo with implications for understanding HCU as a conformational disease. We also lend experimental support to the role of electrostatic interactions in the recently proposed binding modes of SAM leading to hCBS activation and kinetic stabilization.     

Classification of protein sequences by homology modeling and quantitative analysis of electrostatic similarity

Protein electrostatics plays a key role in ligand binding and protein-protein interactions. Therefore, similarities or dissimilarities in electrostatic potentials can be used as indicators of similarities or dissimilarities in protein function. We here describe a method to compare the electrostatic properties within protein families objectively and quantitatively. Three-dimensional structures are built from database sequences by comparative modeling. Molecular potentials are then computed for these with a continuum solvation model by finite difference solution of the Poisson-Boltzmann equation or analytically as a multipole expansion that permits rapid comparison of very large datasets. This approach is applied to 104 members of the Pleckstrin homology (PH) domain family. The deviation of the potentials of the homology models from those of the corresponding experimental structures is comparable to the variation of the potential in an ensemble of structures from nuclear magnetic resonance data or between snapshots from a molecular dynamics simulation. For this dataset, the results for analysis of the full electrostatic potential and the analysis using only monopole and dipole terms are very similar. The electrostatic properties of the PH domains are generally conserved despite the extreme sequence divergence in this family. Notable exceptions from this conservation are seen for PH domains linked to a Db1 homology (DH) domain and in proteins with internal PH domain repeats.     

The Two Sides of Complement C3d: Evolution of Electrostatics in a Link between Innate and Adaptive Immunity

The interaction between complement fragment C3d and complement receptor 2 (CR2) is a key aspect of complement immune system activation, and is a component in a link between innate and adaptive immunities. The complement immune system is an ancient mechanism for defense, and can be found in species that have been on Earth for the last 600 million years. However, the link between the complement system and adaptive immunity, which is formed through the association of the B-cell co-receptor complex, including the C3d-CR2 interaction, is a much more recent adaptation. Human C3d and CR2 have net charges of −1 and +7 respectively, and are believed to have evolved favoring the role of electrostatics in their functions. To investigate the role of electrostatics in the function and evolution of human C3d and CR2, we have applied electrostatic similarity methods to identify regions of evolutionarily conserved electrostatic potential based on 24 homologues of complement C3d and 4 homologues of CR2. We also examine the effects of structural perturbation, as introduced through molecular dynamics and mutations, on spatial distributions of electrostatic potential to identify perturbation resistant regions, generated by so-called electrostatic “hot-spots”. Distributions of electrostatic similarity based on families of perturbed structures illustrate the presence of electrostatic “hot-spots” at the two functional sites of C3d, while the surface of CR2 lacks electrostatic “hot-spots” despite its excessively positive nature. We propose that the electrostatic “hot-spots” of C3d have evolved to optimize its dual-functionality (covalently attaching to pathogen surfaces and interaction with CR2), which are both necessary for the formation B-cell co-receptor complexes. Comparison of the perturbation resistance of the electrostatic character of the homologues of C3d suggests that there was an emergence of a new role of electrostatics, and a transition in the function of C3d, after the divergence of jawless fish.     

Additive Modulation of DNA-DNA Interactions by Interstitial Ions

Quantitative understanding of biomolecular electrostatics, particularly involving multivalent ions and highly charged surfaces, remains lacking. Ion-modulated interactions between nucleic acids provide a model system in which electrostatics plays a dominant role. Using ordered DNA arrays neutralized by spherical cobalt³⁺ hexammine and Mg²⁺ ions, we investigate how the interstitial ions modulate DNA-DNA interactions. Using methods of ion counting, osmotic stress, and x-ray diffraction, we systematically determine thermodynamic quantities, including ion chemical potentials, ion partition, DNA osmotic pressure and force, and DNA-DNA spacing. Analyses of the multidimensional data provide quantitative insights into their interdependencies. The key finding of this study is that DNA-DNA forces are observed to linearly depend on the partition of interstitial ions, suggesting the dominant role of ion-DNA coupling. Further implications are discussed in light of physical theories of electrostatic interactions and like-charge attraction.     

Leveraging enzyme structure-function relationships for functional inference and experimental design: the structure-function linkage database

The study of mechanistically diverse enzyme superfamilies-collections of enzymes that perform different overall reactions but share both a common fold and a distinct mechanistic step performed by key conserved residues-helps elucidate the structure-function relationships of enzymes. We have developed a resource, the structure-function linkage database (SFLD), to analyze these structure-function relationships. Unique to the SFLD is its hierarchical classification scheme based on linking the specific partial reactions (or other chemical capabilities) that are conserved at the superfamily, subgroup, and family levels with the conserved structural elements that mediate them. We present the results of analyses using the SFLD in correcting misannotations, guiding protein engineering experiments, and elucidating the function of recently solved enzyme structures from the structural genomics initiative. The SFLD is freely accessible at http://sfld.rbvi.ucsf.edu.     

Electrostatics in computational protein design

Catalytic activity and protein-protein recognition have proven to be significant challenges for computational protein design. Electrostatic interactions are crucial for these and other protein functions, and therefore accurate modeling of electrostatics is necessary for successfully advancing protein design into the realm of protein function. This review focuses on recent progress in modeling electrostatic interactions in computational protein design, with particular emphasis on continuum models.     

Conserved signature proposed for folding in the lipocalin superfamily

We systematically identify a group of evolutionarily conserved residues proposed for folding in a model beta-barrel superfamily, the lipocalins. The nature of conservation at the structural level is defined and we show that the conserved residues are involved in a network of interactions that form the core of the fold. Exploratory kinetic studies are conducted with a model superfamily member, human serum retinol-binding protein, to examine their role. The present results, coupled with key experimental studies conducted with another lipocalin beta-lactoglobulin, suggest that the evolutionarily conserved regions fold on a faster folding time-scale than the non-conserved regions.     

Electrostatic Effects in Filamentous Protein Aggregation

Electrostatic forces play a key role in mediating interactions between proteins. However, gaining quantitative insights into the complex effects of electrostatics on protein behavior has proved challenging, due to the wide palette of scenarios through which both cations and anions can interact with polypeptide molecules in a specific manner or can result in screening in solution. In this article, we have used a variety of biophysical methods to probe the steady-state kinetics of fibrillar protein self-assembly in a highly quantitative manner to detect how it is modulated by changes in solution ionic strength. Due to the exponential modulation of the reaction rate by electrostatic forces, this reaction represents an exquisitely sensitive probe of these effects in protein-protein interactions. Our approach, which involves a combination of experimental kinetic measurements and theoretical analysis, reveals a hierarchy of electrostatic effects that control protein aggregation. Furthermore, our results provide a highly sensitive method for the estimation of the magnitude of binding of a variety of ions to protein molecules.     

Mutational analysis of conserved outer sphere arginine residues of chalcone synthase

Chalcone synthase (CHS), a key enzyme in flavonoid biosynthesis, catalyses sequential decarboxylative condensations of p-coumaroyl-CoA with three malonyl-CoA molecules and cyclizes the resulting tetraketide intermediate to produce chalcone. Phenylglyoxal, an Arg selective reagent, was found to inactivate the enzyme, although no Arg is found at the active site. Conserved, non-active site Arg residues of CHS were individually mutated and the results were discussed in the context of the 3D structure of CHS. Arg199 and Arg350 were shown to provide important interactions to maintain the structural integrity and foldability of the enzyme. Arg68, Arg172 and Arg328 interact with highly conserved Gln33/Phe215, Glu380 and Asp311/Glu314, respectively, thus helping position the catalytic Cys-His-Asn triad and the (372)GFGPG loop in correct topology at the active site. In particular, a mutation of Arg172 resulted in selective impairment in the cyclization activities of CHS and stilbene synthase, a related enzyme that catalyses a different cyclization of the same tetraketide intermediate. These Arg residues and their interactions are well conserved in other enzymes of the CHS superfamily, suggesting that they may serve similar functions in other enzymes. Mutations of Arg68 and Arg328 had been found in mutant plants that showed impaired CHS activity.     

Structural basis for the activation of phenylalanine in the non-ribosomal biosynthesis of gramicidin S.

The non-ribosomal synthesis of the cyclic peptide antibiotic gramicidin S is accomplished by two large multifunctional enzymes, the peptide synthetases 1 and 2. The enzyme complex contains five conserved subunits of approximately 60 kDa which carry out ATP-dependent activation of specific amino acids and share extensive regions of sequence similarity with adenylating enzymes such as firefly luciferases and acyl-CoA ligases. We have determined the crystal structure of the N-terminal adenylation subunit in a complex with AMP and L-phenylalanine to 1.9 A resolution. The 556 amino acid residue fragment is folded into two domains with the active site situated at their interface. Each domain of the enzyme has a similar topology to the corresponding domain of unliganded firefly luciferase, but a remarkable relative domain rotation of 94 degrees occurs. This conformation places the absolutely conserved Lys517 in a position to form electrostatic interactions with both ligands. The AMP is bound with the phosphate moiety interacting with Lys517 and the hydroxyl groups of the ribose forming hydrogen bonds with Asp413. The phenylalanine substrate binds in a hydrophobic pocket with the carboxylate group interacting with Lys517 and the alpha-amino group with Asp235. The structure reveals the role of the invariant residues within the superfamily of adenylate-forming enzymes and indicates a conserved mechanism of nucleotide binding and substrate activation.     

A unified view of the role of electrostatic interactions in modulating the gating of Cys loop receptors

In the Cys loop superfamily of ligand-gated ion channels, a global conformational change, initiated by agonist binding, results in channel opening and the passage of ions across the cell membrane. The detailed mechanism of channel gating is a subject that has lent itself to both structural and electrophysiological studies. Here we defined a gating interface that incorporates elements from the ligand binding domain and transmembrane domain previously reported as integral to proper channel gating. An overall analysis of charged residues within the gating interface across the entire superfamily showed a conserved charging pattern, although no specific interacting ion pairs were conserved. We utilized a combination of conventional mutagenesis and the high precision methodology of unnatural amino acid incorporation to study extensively the gating interface of the mouse muscle nicotinic acetylcholine receptor. We found that charge reversal, charge neutralization, and charge introduction at the gating interface are often well tolerated. Furthermore, based on our data and a reexamination of previously reported data on gamma-aminobutyric acid, type A, and glycine receptors, we concluded that the overall charging pattern of the gating interface, and not any specific pairwise electrostatic interactions, controls the gating process in the Cys loop superfamily.     

Electrostatic Interactions in the Binding Pathway of a Transient Protein Complex Studied by NMR and Isothermal Titration Calorimetry

Much of our knowledge of protein binding pathways is derived from extremely stable complexes that interact very tightly, with lifetimes of hours to days. Much less is known about weaker interactions and transient complexes because these are challenging to characterize experimentally. Nevertheless, these types of interactions are ubiquitous in living systems. The combination of NMR relaxation dispersion Carr–Purcell–Meiboom–Gill (CPMG) experiments and isothermal titration calorimetry allows the quantification of rapid binding kinetics for complexes with submillisecond lifetimes that are difficult to study using conventional techniques. We have used this approach to investigate the binding pathway of the Src homology 3 (SH3) domain from the Fyn tyrosine kinase, which forms complexes with peptide targets whose lifetimes are on the order of about a millisecond. Long range electrostatic interactions have been shown to play a critical role in the binding pathways of tightly binding complexes. The role of electrostatics in the binding pathways of transient complexes is less well understood. Similarly to previously studied tight complexes, we find that SH3 domain association rates are enhanced by long range electrostatics, whereas short range interactions are formed late in the docking process. However, the extent of electrostatic association rate enhancement is several orders of magnitudes less, whereas the electrostatic-free basal association rate is significantly greater. Thus, the SH3 domain is far less reliant on electrostatic enhancement to achieve rapid association kinetics than are previously studied systems. This suggests that there may be overall differences in the role played by electrostatics in the binding pathways of extremely stable versus transient complexes.     

Automated computational framework for the analysis of electrostatic similarities of proteins

Charge plays an important role in protein-protein interactions. In the case of excessively charged proteins, their electrostatic potentials contribute to the processes of recognition and binding with other proteins or ligands. We present an automated computational framework for determining the contribution of each charged amino acid to the electrostatic properties of proteins, at atomic resolution level. This framework involves computational alanine scans, calculation of Poisson-Boltzmann electrostatic potentials, calculation of electrostatic similarity distances (ESDs), hierarchical clustering analysis of ESDs, calculation of solvation free energies of association, and visualization of the spatial distributions of electrostatic potentials. The framework is useful to classify families of mutants with similar electrostatic properties and to compare them with the parent proteins in the complex. The alanine scan mutants introduce perturbations in the local electrostatic properties of the proteins and aim in delineating the contribution of each mutated amino acid in the spatial distribution of electrostatic potential, and in biological function when electrostatics is a dominant contributing factor in protein-protein interactions. The framework can be used to design new proteins with tailored electrostatic properties, such as immune system regulators, inhibitors, and vaccines, and in guiding experimental studies. We present an example for the interaction of the immune system protein C3d (the d-fragment of complement protein C3) with its receptor CR2, and we discuss our data in view of a binding site controversy.     

Glutathione-Binding Site of a Bombyx mori Theta-Class Glutathione Transferase

The glutathione transferase (GST) superfamily plays key roles in the detoxification of various xenobiotics. Here, we report the isolation and characterization of a silkworm protein belonging to a previously reported theta-class GST family. The enzyme (bmGSTT) catalyzes the reaction of glutathione with 1-chloro-2,4-dinitrobenzene, 1,2-epoxy-3-(4-nitrophenoxy)-propane, and 4-nitrophenethyl bromide. Mutagenesis of highly conserved residues in the catalytic site revealed that Glu66 and Ser67 are important for enzymatic function. These results provide insights into the catalysis of glutathione conjugation in silkworm by bmGSTT and into the metabolism of exogenous chemical agents.     

The role of electrostatics in the interaction of the inhibitory region of troponin I with troponin C

We have addressed the electrostatic interactions occurring between the inhibitory region of cardiac troponin I with the C-lobe of troponin C using scanning glycine mutagenesis of the inhibitory region. We report variations in the electric potentials due to mutation of charged residues within this complex based upon the solved NMR structure (1OZS). These results demonstrate the importance of electrostatics within this complex, and it is proposed that electrostatic interactions are integral to the formation and function of larger ternary troponin complexes. To address this hypothesis, we report (15)N NMR relaxation measurements, which suggest that, within a ternary complex involving the C-lobe and the N-terminal region of troponin I (residues 34-71), the inhibitory region maintains the electrostatic interactions with the E-helix of the C-lobe as observed within the binary complex. These results imply that, in solution, the cardiac troponin complex behaves in a manner consistent with that of the crystal structure of the skeletal isoform (1YTZ). A cardiac troponin complex possessing domain orientations similar to that of the skeletal isoform provides structural insights into altered troponin I activities as observed for the familial hypertrophic cardiomyopathy mutation R144G and phosphorylation of Thr142.     

Role of electrostatic interactions in PDZ domain ligand recognition

PDZ domains are protein-protein interaction modules that normally recognize short C-terminal peptides. The apparent requirement for a ligand with a free terminal carboxylate group has led to the proposal that electrostatic interactions with the terminus play a significant role in recognition. However, this model has been called into question by the more recent finding that PDZ domains can recognize some internal peptide motifs that occur within a specific secondary structure context. Although these motifs bind at the same interface, they lack a terminal charge. Here we have investigated the role of electrostatics in PDZ-mediated recognition in the mouse alpha1-syntrophin PDZ domain by examining the salt dependence of binding to both terminal and internal ligands and the effects of mutating a conserved basic residue previously proposed to play a role in electrostatic recognition. These studies indicate that direct electrostatic interactions with the peptide terminus do not play a significant energetic role in binding. Additional chemical modification studies of the peptide terminus support a model in which steric and hydrogen bonding complementarity play a primary role in recognition specificity. Peptides with a free carboxy terminus, or presented within a specific structural context, can satisfy these requirements.