Analysis of affinity and specificity in an EF-hand site using double mutant cycles

The effects of three mutations on the EF-hand Ca(2+)/Mg(2+) binding site of smooth muscle myosin regulatory light chain (RLC) were studied: D5S, in which an aspartate is replaced by a serine in position 5 of the loop; D9E, in which an aspartate is replaced by a glutamate in position 9; and D12E, in which the aspartate in position 12 is replaced by a glutamate. All possible combinations of the three mutations were produced. The single mutants D5S and D9E and the double mutant D5S/D9E have low affinity for Ca(2+). All the mutants containing mutation D12E are Ca(2+)-specific and have higher affinities than wild type, even when containing mutations D5S or D9E. All of the mutants studied have lower affinity for Mg(2+) than the wild-type protein. As expected, the changes in binding free energy that each mutant produces depend on the residues present at the other positions of the site, since the mutated positions are very close in the protein structure. Coupling energies are about the same for all pairs of mutants when binding Ca(2+), but can have different values when binding Mg(2+). D5S and D9E have a large negative coupling energy for Mg(2+) binding which suggests an interaction between these two positions. When mutation D12E is present, the coupling energy for Mg(2+) binding between D5S and D9E is much lower, suggesting that this interaction occurs only if an aspartate is in position 12. Glutamate in position 9 may be able to coordinate Mg(2+) directly in the double mutant D5S/D9E.     

The geometry of domain combination in proteins.

Most proteins in genomes are the result of the recombination of two or more domains. It has been found that if proteins are formed by a combination of domains from superfamilies A and B, then the domains may occur in the sequential order AB or BA but only in about 2% of cases do they occur in both sequential orders. The classical Rossmann domains of known structure are combined with catalytic domains from seven different superfamilies. In addition, there are eight cases where structures with both AB and BA domain combinations are known. For these two sets of structures, we analysed: (i) the relative orientation of the domains; (ii) the type of domain connection; (iii) the structure of the interdomain links; and (iv) domain function. The results of this analysis indicate that in most cases domain order is conserved because recombination of the domains has only occurred once during the course of evolution. Functional reasons become important when the domain connections are short. In seven out of the eight known cases where domains are combined in the AB and BA sequential orders they have different geometrical relationships that give them different functional properties.     

Evolution of EF-hand calcium-modulated proteins. IV. Exon shuffling did not determine the domain compositions of EF-hand Proteins

Summary In the previous three reports in this series we demonstrated that the EF-hand family of proteins evolved by a complex pattern of gene duplication, transposition, and splicing. The dendrograms based on exon sequences are nearly identical to those based on protein sequences for troponin C, the essential light chain myosin, the regulatory light chain, and calpain. This validates both the computational methods and the dendrograms for these subfamilies. The proposal of congruence for calmodulin, troponin, C, essential light chain, and regulatory light chain was confirmed. There are, however, significant differences in the calmodulin dendrograms computed from DNA and from protein sequences. In this study we find that introns are distributed throughout the EF-hand domain and the interdomain regions. Further, dendrograms based on intron type and distribution bear little resemblance to those based on protein or on DNA sequences. We conclude that introns are inserted, and probably deleted, with relatively high frequency. Further, in the EF-hand family exons do not correspond to structural domains and exon shuffling played little if any role in the evolution of this widely distributed homolog family. Calmodulin has had a turbulent evolution. Its dendrograms based on protein sequence, exon sequence, 3′-tail sequence, intron sequences, and intron positions all show significant differences.     

Hydration-coupled dynamics in proteins studied by neutron scattering and NMR: the case of the typical EF-hand calcium-binding parvalbumin

The influence of hydration on the internal dynamics of a typical EF-hand calciprotein, parvalbumin, was investigated by incoherent quasi-elastic neutron scattering (IQNS) and solid-state 13C-NMR spectroscopy using the powdered protein at different hydration levels. Both approaches establish an increase in protein dynamics upon progressive hydration above a threshold that only corresponds to partial coverage of the protein surface by the water molecules. Selective motions are apparent by NMR in the 10-ns time scale at the level of the polar lysyl side chains (externally located), as well as of more internally located side chains (from Ala and Ile), whereas IQNS monitors diffusive motions of hydrogen atoms in the protein at time scales up to 20 ps. Hydration-induced dynamics at the level of the abundant lysyl residues mainly involve the ammonium extremity of the side chain, as shown by NMR. The combined results suggest that peripheral water-protein interactions influence the protein dynamics in a global manner. There is a progressive induction of mobility at increasing hydration from the periphery toward the protein interior. This study gives a microscopic view of the structural and dynamic events following the hydration of a globular protein.     

Structural plasticity and the evolution of antibody affinity and specificity

The germline precursor to the ferrochelatase antibody 7G12 was found to bind the polyether jeffamine in addition to its cognate hapten N-methylmesoporphyrin. A comparison of the X-ray crystal structures of the ligand-free germline Fab and its complex with either hapten or jeffamine reveals that the germline antibody undergoes significant conformational changes upon the binding of these two structurally distinct ligands, which lead to increased antibody-ligand complementarity. The five somatic mutations introduced during affinity maturation lead to enhanced binding affinity for hapten and a loss in affinity for jeffamine. Moreover, a comparison of the crystal structures of the germline and affinity-matured antibodies reveals that somatic mutations not only fix the optimal binding site conformation for the hapten, but also introduce interactions that interfere with the binding of non-hapten molecules. The structural plasticity of this germline antibody and the structural effects of the somatic mutations that result in enhanced affinity and specificity for hapten likely represent general mechanisms used by the immune response, and perhaps primitive proteins, to evolve high affinity, selective receptors for so many distinct chemical structures.     

Evolution of EF-hand calcium-modulated proteins. V. The genes encoding EF-hand proteins are not clustered in mammalian genomes

Summary The chromosomal assignments of genes belonging to the EF-hand family which have a common origin are compiled in this article. So far data are available from 27 human gene loci belonging to 6 subfamilies and 8 murine loci belonging to 4 subfamilies. Chromosomal localization has been obtained by somatic-cell hybrid analysis using the Southern blot technique or PCR amplification, metaphase spread in situ hybridization, or isolation of the particular genes from chromosome-specific libraries. Except for genes of the S-100 alpha proteins which are grouped on human chromosome 1q12-25 and mouse chromosome 3, no linkage has been found for genes encoding EF-hand proteins, indicating absence of selective pressure for maintaining chromosomal clustering. Six of these genes map to known syntenic groups conserved in the human and mouse genomes. This suggests that chromosomal translocations occurred before divergence of these species. The possible significance of chromosomal positioning with respect to nearby located known genes and genetic disease loci is discussed.     

Determination of residue specificity in the EF-hand of troponin C for Ca2+ coordination, by genetic engineering.

Utilizing protein engineering of troponin C (TnC), combined with the physiology of skinned fibers, the present study sought to delineate the mechanisms for metal ion coordination and sensitivity in the sites (EF-hands) that execute the Ca2+ switch for contraction. A total TnC-encoding gene comprising multiple target sequences for restriction enzymes was synthesized, furnishing a pliant molecular handle to manipulate sites I and II in the NH2 terminus of the protein. Of the six positions (X, Y, Z, -Y, -X, and -Z) essential for metal ion chelation in a typical EF-hand, invariably the X position has aspartate, and -Z position has glutamate. In the X position of site II, mutation of aspartate for either glutamate (gamma-carboxylate) or asparagine (same side chain length as aspartate) yielded functionally inactive proteins with concomitantly diminished Ca2+ binding capacity. Similarly, in -Z position (site I), neither aspartate nor glutamine were compatible in exchange for the conserved glutamate. In contrast, for the Y coordinate of site II, a preference for asparagine comparable to that for wild-type aspartate was detected, but glutamate was impermissible. Evidently, physicochemical and steric factors both are critical in governing the mechanism for metal ion chelation in TnC in a physiological milieu. Furthermore, the findings manifest that the quaternary structure of hydrated TnC restrains the EF-hands during on-off operation of the Ca2+ switch.     

Remodeling of the AB site of rat parvalbumin and oncomodulin into a canonical EF-hand.

Parvalbumin (PV) and the homologous protein oncomodulin (OM) contain three EF-hand motifs, but the first site (AB) cannot bind Ca2+. Here we aimed to recreate the putative ancestral proteins [D19-28E]PV and [D19-28E]OM by replacing the 10-residue-long nonfunctional loop in the AB site by a 12-residue canonical loop. To create an optical conformational probe we also expressed the homologs with a F102W replacement. Unexpectedly, in none of the proteins did the mutation reactivate the AB site. The AB-remodeled parvalbumins bind two Ca2+ ions with strong positive cooperativity (nH = 2) and moderate affinity ([Ca2+]0.5 = 2 microM), compared with [Ca2+]0.5 = 37 nM and nH = 1 for the wild-type protein. Increasing Mg2+ concentrations changed nH from 2 to 0.65, but without modification of the [Ca2+]0. 5-value. CD revealed that the Ca2+ and Mg2+ forms of the remodeled parvalbumins lost one-third of their alpha helix content compared with the Ca2+ form of wild-type parvalbumin. However, the microenvironment of single Trp residues in the hydrophobic cores, monitored using intrinsic fluorescence and difference optical density, is the same. The metal-free remodeled parvalbumins possess unfolded conformations. The AB-remodeled oncomodulins also bind two Ca2+ with [Ca2+]0.5 = 43 microM and nH = 1.45. Mg2+ does not affect Ca2+ binding. Again the Ca2+ forms display two-thirds of the alpha-helical content in the wild-type, while their core is still strongly hydrophobic as monitored by Trp and Tyr fluorescence. The metal-free oncomodulins are partially unfolded and seem not to possess a hydrophobic core. Our data indicate that AB-remodeled parvalbumin has the potential to regulate cell functions, whereas it is unlikely that [D19-28E]OM can play a regulatory role in vivo. The predicted evolution of the AB site from a canonical to an abortive EF-hand may have been dictated by the need for stronger interaction with Mg2+ and Ca2+, and a high conformational stability of the metal-free forms.     

Optimizing the affinity and specificity of proteins with molecular display

Affinity maturation of receptor-ligand interactions represents an important area of academic and pharmaceutical research. Improving affinity and specificity of proteins can tailor potency for both in vivo and in vitro applications. A number of different display platforms including phage display, bacterial and yeast display, ribosome display, and mRNA display can optimize protein affinity and specificity. Here, we will review the advantages and disadvantages of these molecular display methods with a focus on their suitability for protein affinity maturation.     

Lid domain plasticity and lipid flexibility modulate enzyme specificity in human monoacylglycerol lipase

Human monoacylglycerol lipase (MAGL) is a membrane-interacting enzyme that generates pro-inflammatory signaling molecules. For this reason, MAGL inhibition is a promising strategy to treat pain, cancer, and neuroinflammatory diseases. MAGL can hydrolyze monoacylglycerols bearing an acyl chain of different lengths and degrees of unsaturation, cleaving primarily the endocannabinoid 2-arachidonoylglycerol. Importantly, the enzymatic binding site of MAGL is confined by a 75-amino-acid-long, flexible cap domain, named ‘lid domain’, which is structurally similar to that found in several other lipases. However, it is unclear how lid domain plasticity affects catalysis in MAGL. By integrating extensive molecular dynamics simulations and free-energy calculations with mutagenesis and kinetic experiments, we here define a lid-domain-mediated mechanism for substrate selection and binding in MAGL catalysis. In particular, we clarify the key role of Phe159 and Ile179, two conserved residues within the lid domain, in regulating substrate specificity in MAGL. We conclude by proposing that other structurally related lipases may share this lid-domain-mediated mechanism for substrate specificity.     

Identification of the critical structural determinants of the EF-hand domain arrangements in calcium binding proteins

EF-hand calcium binding proteins (CaBPs) share strong sequence homology, but exhibit great diversity in structure and function. Thus although calmodulin (CaM) and calcineurin B (CNB) both consist of four EF hands, their domain arrangements are quite distinct. CaM and the CaM-like proteins are characterized by an extended architecture, whereas CNB and the CNB-like proteins have a more compact form. In this study, we performed structural alignments and molecular dynamics (MD) simulations on 3 CaM-like proteins and 6 CNB-like proteins, and quantified their distinct structural and dynamical features in an effort to establish how their sequences specify their structures and dynamics. Alignments of the EF2-EF3 region of these proteins revealed that several residues (not restricted to the linker between the EF2 and EF3 motifs) differed between the two groups of proteins. A customized inverse folding approach followed by structural assessments and MD simulations established the critical role of these residues in determining the structure of the proteins. Identification of the critical determinants of the two different EF-hand domain arrangements and the distinct dynamical features relevant to their respective functions provides insight into the relationships between sequence, structure, dynamics and function among these EF-hand CaBPs.     

Identification of the critical structural determinants of the EF-hand domain arrangements in calcium binding proteins

EF-hand calcium binding proteins (CaBPs) share strong sequence homology, but exhibit great diversity in structure and function. Thus although calmodulin (CaM) and calcineurin B (CNB) both consist of four EF hands, their domain arrangements are quite distinct. CaM and the CaM-like proteins are characterized by an extended architecture, whereas CNB and the CNB-like proteins have a more compact form. In this study, we performed structural alignments and molecular dynamics (MD) simulations on 3 CaM-like proteins and 6 CNB-like proteins, and quantified their distinct structural and dynamical features in an effort to establish how their sequences specify their structures and dynamics. Alignments of the EF2-EF3 region of these proteins revealed that several residues (not restricted to the linker between the EF2 and EF3 motifs) differed between the two groups of proteins. A customized inverse folding approach followed by structural assessments and MD simulations established the critical role of these residues in determining the structure of the proteins. Identification of the critical determinants of the two different EF-hand domain arrangements and the distinct dynamical features relevant to their respective functions provides insight into the relationships between sequence, structure, dynamics and function among these EF-hand CaBPs.     

Solvation energetics and conformational change in EF-hand proteins

Calmodulin and other members of the EF-hand protein family are known to undergo major changes in conformation upon binding Ca(2+). However, some EF-hand proteins, such as calbindin D9k, bind Ca(2+) without a significant change in conformation. Here, we show the importance of a precise balance of solvation energetics to conformational change, using mutational analysis of partially buried polar groups in the N-terminal domain of calmodulin (N-cam). Several variants were characterized using fluorescence, circular dichroism, and NMR spectroscopy. Strikingly, the replacement of polar side chains glutamine and lysine at positions 41 and 75 with nonpolar side chains leads to dramatic enhancement of the stability of the Ca(2+)-free state, a corresponding decrease in Ca(2+)-binding affinity, and an apparent loss of ability to change conformation to the open form. The results suggest a paradigm for conformational change in which energetic strain is accumulated in one state in order to modulate the energetics of change to the alternative state.     

Potential influence of Asp in the Ca2+ coordination position 5 of parvalbumin on the calcium-binding affinity: a computational study

Parvalbumins (PV) are calcium-binding proteins, all sharing the common helix-loop-helix (EF-hand) motif. This motif contains a central twelve-residue Ca(2+)-binding loop with the flanking helices positioned roughly perpendicular to each other. The precise role of these coordination residues has been the subject of intense studies. In this work, we focus on the coordination position 5 in the CD Ca(2+)-binding site of silver hake parvalbumin isoform B (SHPV-B). The most common residue at site 5 of calcium-binding loop in canonical EF-hands is Asp [B.J. Marsden, G.S. Shaw, B.D. Sykes, Biochem. Cell Biol. 68 (1990) 587-601], but in the CD site of PV, this position is almost always serine (Ser). The substitution of Ser with Asp will add the 5th carboxylate residue in the CD coordination sphere. However, as predicted by the acid pair hypothesis, the Ca(2+)-binding affinity would be maximized in an EF-hand motif that has four carboxylate ligands paired along the +/-x, and +/-z-axes [R.E. Reid, R.S. Hodges, J. Theor. Biol. 84 (1980) 401-444]. Molecular dynamics simulations and free energy calculations were employed to investigate the influence of Ser to Asp mutation at position 5 on calcium-binding affinity. We found that the Asp variant exhibited remarkable stability during the entire molecular dynamics simulation, with not only the retention of the Ca(2+)-binding site, but also increased compactness in the coordination sphere. The S55D fragment also accommodated Ca(2+) well. We conclude that the reason why Asp which is the most common residue at site 5 of calcium-binding loop in canonical EF-hands has never been identified at this position experimentally for PVs might be related to its physiological functions.     

The N-terminal domain determines the affinity and specificity of H1 binding to chromatin

Linker histone H1, one of the most abundant nuclear proteins in multicellular eukaryotes, is a key component of the chromatin structure mainly due to its role in the formation and maintenance of the 30nm chromatin fiber. It has a three-domain structure; a central globular domain flanked by a short N-terminal domain and a long, highly basic C-terminal domain. Previous studies have shown that the binding abilities of H1 are at large determined by the properties of the C-terminal domain; much less attention has been paid to role of the N-terminal domain. We have previously shown that H1 can be reconstituted via cytoplasmic mRNA injection in Xenopus oocytes, cells that lack somatic H1. The heterologously expressed H1 proteins are incorporated into in vivo assembled chromatin at specific sites and the binding event is monitored as an increase in nucleosomal repeat length (NRL). Using this setup we have here compared the binding properties of wt-H1.4 and hH1.4 devoid of its N-terminal domain (ΔN-hH1.4). The ΔN-hH1.4 displays a drastically lower affinity for chromatin binding as compared to the wild type hH1.4. Our data also indicates that ΔN-hH1.4 is more prone to unspecific chromatin binding than the wild type. We conclude that the N-terminal domain of H1 is an important determinant of affinity and specificity of H1-chromatin interactions.     

Helical vs. zigzag coordination polymer: Influence of structural preference of metal-ion coordination geometry

Two new coordination polymers [(L⁶)CoCl₂]_∞ (1) and [(L⁶)CdCl₂]_∞ (2) (L⁶ = α,α′-bis(pyrazolyl)-m-xylene) have been synthesized and characterized by X-ray crystallography. Structural analysis of 1 (space group P2₁) revealed that each Co(II) center is coordinated by two chloride ions and two pyrazole nitrogens from two different L⁶ ligands, giving rise to pseudo-tetrahedral coordination environment around Co(II). For this structure 1D helical polymeric chain running along crystallographic 2₁ axis is observed. While for 2 (space group ), each Cd(II) center is coordinated by two chloride ions, two pyrazole nitrogens from two different L⁶ ligands, and a chloride ion from the adjacent layer. It provides distorted trigonal bipyramidal geometry around Cd(II), which leads to the formation of a 2D zigzag polymer. In continuation of our recent activity on inorganic crystal engineering, from the standpoint of metal-ligand coordination chemistry and C-H?Cl₂M^II hydrogen bonding, in this work we have investigated the potential of coordination tectons1 and 2 in bringing about generality in helix/zigzag network formation, strengthened primarily by C-H?Cl hydrogen-bonding interactions.     

Free-energy linkage between folding and calcium binding in EF-hand proteins

Troponin is the singular Ca²⁺-sensitive protein in the contraction of vertebrate striated muscles. Troponin C (TnC), the Ca²⁺-binding subunit of the troponin complex, has two distinct domains, C and N, which have different properties despite their extensive structural homology. In this work, we analyzed the thermodynamic stability of the isolated N-domain of TnC using a fluorescent mutant with Phe 29 replaced by Trp (F29W/N-domain, residues 1-90). The complete unfolding of the N-domain of TnC in the absence or presence of Ca²⁺ was achieved by combining high hydrostatic pressure and urea, a maneuver that allowed us to calculate the thermodynamic parameters (ΔV and ΔG_atm). In this study, we propose that part of the affinity for Ca²⁺ is contributed by the free-energy change of folding of the N- and C-domains that takes place when Ca²⁺ binds. The importance of the free-energy change for the structural and regulatory functions of the TnC isolated domains was evaluated. Our results shed light on how the coupling between folding and ion binding contributes to the fine adjustment of the affinity for Ca²⁺ in EF-hand proteins, which is crucial to function.     

Hypoxia-inducible genes encoding small EF-hand proteins in rice and tomato

Rice has evolved metabolic and morphological adaptations to low-oxygen stress to grow in submerged paddy fields. To characterize the molecular components that mediate the response to hypoxia in rice, we identified low-oxygen stress early response genes by microarray analysis. Among the highly responsive genes, five genes, OsHREF1 to OsHREF5, shared strong homology. They encoded small proteins harboring two EF-hands, typical Ca(2+)-binding motifs. Homologous genes were found in many land plants, including SlHREF in tomato, which is also strongly induced by hypoxia. SlHREF induction was detected in both roots and shoots of tomato plants under hypoxia. With the exception of OsHREF5, OsHREF expression was unaffected by drought, salinity, cold, or osmotic stress. Fluorescent signals of green fluorescent protein-fused OsHREFs were detected in the cytosol and nucleus. Ruthenium red, an inhibitor of intracellular Ca(2+) release, repressed induction of OsHREF1-4 under hypoxia. The HREFs may be related to the Ca(2+) response to hypoxia.