Chemical synthesis and evaluation of a backbone‐cyclized minimized 2‐helix Z‐domain

The Z‐molecule is a small, engineered IgG‐binding affinity protein derived from the immunoglobulin‐binding domain B of Staphylococcus aureus protein A. The Z‐domain consists of 58 amino acids forming a well‐defined antiparallel three‐helix structure. Two of the three helices are involved in ligand binding, whereas the third helix provides structural support to the three‐helix bundle. The small size and the stable three‐helix structure are two attractive properties comprised in the Z‐domain, but a further reduction in size of the protein is valuable for several reasons. Reduction in size facilitates synthetic production of any protein‐based molecule, which is beneficial from an economical viewpoint. In addition, a smaller protein is easier to manipulate through chemical modifications. By omitting the third stabilizing helix from the Z‐domain and joining the N‐ and C‐termini by a native peptide bond, the affinity protein obtains the advantageous properties of a smaller scaffold and in addition becomes resistant to exoproteases. We here demonstrate the synthesis and evaluation of a novel cyclic two‐helix Z‐domain. The molecule has retained affinity for its target protein, is resistant to heat treatment, and lacks both N‐ and C‐termini. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Functional properties of calbindin D9K mutants with exchanged Ca2+ binding sites

We have constructed three different engineered proteins based on calbindin D9K by either exchanging the two calcium binding sites within the protein or making the amino acid sequence of the two calcium binding sites identical. The individual calcium binding sites of the engineered proteins retain most of their ion binding characteristics as well as the basal structure of their Ca2+ ligand sphere in the new environment. Even the protein with its sites interchanged, a mutation involving 30 amino acids out of a total of 75, still binds calcium with an affinity as high as that of many natural EF-hand proteins.     

Crystal structure of the HNF4 alpha ligand binding domain in complex with endogenous fatty acid ligand

HNF4 alpha is an orphan member of the nuclear receptor family with prominent functions in liver, gut, kidney and pancreatic beta cells. We have solved the x-ray crystal structure of the HNF4 alpha ligand binding domain, which adopts a canonical fold. Two conformational states are present within each homodimer: an open form with alpha helix 12 (alpha 12) extended and collinear with alpha 10 and a closed form with alpha 12 folded against the body of the domain. Although the protein was crystallized without added ligands, the ligand binding pockets of both closed and open forms contain fatty acids. The carboxylic acid headgroup of the fatty acid ion pairs with the guanidinium group of Arg(226) at one end of the ligand binding pocket, while the aliphatic chain fills a long, narrow channel that is lined with hydrophobic residues. These findings suggest that fatty acids are endogenous ligands for HNF4 alpha and establish a framework for understanding how HNF4 alpha activity is enhanced by ligand binding and diminished by MODY1 mutations.     

Beyond the Binding Site: The Role of the β2 – β3 Loop and Extra-Domain Structures in PDZ Domains

A general paradigm to understand protein function is to look at properties of isolated well conserved domains, such as SH3 or PDZ domains. While common features of domain families are well understood, the role of subtle differences among members of these families is less clear. Here, molecular dynamics simulations indicate that the binding mechanism in PSD95-PDZ3 is critically regulated via interactions outside the canonical binding site, involving both the poorly conserved loop and an extra-domain helix. Using the CRIPT peptide as a prototypical ligand, our simulations suggest that a network of salt-bridges between the ligand and this loop is necessary for binding. These contacts interconvert between each other on a time scale of a few tens of nanoseconds, making them elusive to X-ray crystallography. The loop is stabilized by an extra-domain helix. The latter influences the global dynamics of the domain, considerably increasing binding affinity. We found that two key contacts between the helix and the domain, one involving the loop, provide an atomistic interpretation of the increased affinity. Our analysis indicates that both extra-domain segments and loosely conserved regions play critical roles in PDZ binding affinity and specificity.     

Comparative normal mode analysis of LFA-1 integrin I-domains

The conformational dynamics of the Inserted domain (I-domain) from the lymphocyte function-associated antigen-1 (LFA-1) was investigated by normal mode analysis of multiple structures of the low, intermediate, and high affinity states. LFA-1 is an integrin expressed on leukocytes and is of critical importance in adhesion reactions, like antigen-specific responses, homing, and diapedesis. The main ligand binding site of LFA-1 is the I-domain, which recognizes intercellular adhesion molecules (ICAMs), members of the immunoglobulin superfamily. From experimental crystal structures, a large-scale conformational change of, among others, the α7 helix of the I-domain has been observed leading to the proposal that these structural changes are linked to the conformational regulation of LFA-1. The results from the present calculations show that structural changes of the α7 helix consistent with those observed in the crystal structures are significantly sampled by the low frequency modes. This was found to be particularly true for the low affinity state of the I-domain, indicating that low frequency motions favor the conformational transition implicated in activation. However, beyond the simple downward shift of the helix implied by the crystal structures, the calculations further show that there is a noticeable swinging-out motion of the helix. The consequences of this motion are discussed in the context of integrin activation and inhibition. Moreover, significant changes in the atomic-level dynamics and in long-range correlated motions of the I-domain were found to occur upon binding of the natural ligand ICAM. These changes were more local upon binding of an allosteric inhibitor. The present study opens the question of how changes in dynamics may contribute to the long-range transmission of signal upon ICAM binding by the LFA-1 I-domain.     

Design and folding of a multidomain protein

To test whether the folding process of a large protein can be understood on the basis of the folding behavior of the domains that constitute it, we coupled two well-studied small -helical proteins, the B-domain of protein A (60 amino acids) and Rd-apocytochrome b562 (Rd-apocyt b562, 106 amino acids), by fusing the C-terminal helix of the B-domain of protein A with the N-terminal helix of Rd-apocyt b562 without changing their hydrophobic core residues. The success of the design was confirmed by determining the structure of the engineered protein with multidimensional NMR methods. Kinetic studies showed that the logarithms of the folding/unfolding rate constants of the engineered protein are linearly dependent on concentrations of guanidinium chloride in the measurable range from 1.7 to 4 M. Their slopes (m-values) are close to those of Rd-apocyt b562. In addition, the 1H-15N HSQC spectrum taken at 1.5 M guanidinium chloride reveals that only the Rd-apocyt b562 domain in the designed protein remained folded. These results suggest that the two domains have weak energetic coupling. Interestingly, the redesigned protein folds faster than Rd-apocyt b562, suggesting that the fused helix stabilizes the rate-limiting transition state.     

Design of improved protein inhibitors of HIV-1 cell entry: Optimization of electrostatic interactions at the binding interface

Continuum electrostatic methods are a powerful tool for the analysis and design of biomolecular complexes, with methodologies that allow for the detailed analysis of the electrostatic contributions to binding affinities and procedures for computing the properties of electrostatically optimal ligands. We have applied these methods to the design of improved inhibitors of HIV-1 cell entry. HIV infection of a cell requires viral-cell membrane fusion, an event partially driven by a large-scale conformational change in the viral membrane glycoprotein gp41. This transformation involves the docking of a helix from the C-terminal region of three gp41 chains against a pre-formed trimeric-coiled coil; several protein constructs that inhibit membrane fusion act by binding to an isolated C-terminal helix and blocking the formation of the fusogenic structure. A detailed analysis of the electrostatic contributions to the binding of one such inhibitor (5-Helix) to a C-terminal helix was performed using the X-ray crystal structure of the core of the HIV-1 gp41 ectodomain as a structural model, and several residues on 5-Helix that make substantial contributions to binding, both favorable and unfavorable, were identified. An electrostatic affinity optimization methodology was applied to the side chains of 5-Helix, with the results showing that significant improvements in binding affinity are possible if the electrostatic contributions to the binding free energy are optimized. Several mutations accessible by experimental methods are suggested, with calculated improvements in binding affinity of as much as 500-fold and greater.     

Insights into Lou Gehrig's disease from the structure and instability of the A4V mutant of human Cu,Zn superoxide dismutase

We have engineered enhanced DNA-binding function into the a1 homeodomain by making changes in a loop distant from the DNA-binding surface. Comparison of the free and bound a1 structures suggested a mechanism linking van der Waals stacking changes in this loop to the ordering of a final turn in the DNA-binding helix of a1. Inspection of the protein sequence revealed striking differences in amino acid identity at positions 24 and 25 compared to related homeodomain proteins. These positions lie in the loop connecting helix-1 and helix-2, which is involved in heterodimerization with the alpha 2 protein. A series of single and double amino acid substitutions (a1-Q24R, a1-S25Y, a1-S25F and a1-Q24R/S25Y) were engineered, expressed and purified for biochemical and biophysical study. Calorimetric measurements and HSQC NMR spectra confirm that the engineered variants are folded and are equally or more stable than the wild-type a1 homeodomain. NMR analysis of a1-Q24R/S25Y demonstrates that the DNA recognition helix (helix-3) is extended by at least one turn as a result of the changes in the loop connecting helix-1 and helix-2. As shown by EMSA, the engineered variants bind DNA with enhanced affinity (16-fold) in the absence of the alpha 2 cofactor and the variant alpha 2/a1 heterodimers bind cognate DNA with specificity and affinity reflective of the enhanced a1 binding affinity. Importantly, in vivo assays demonstrate that the a1-Q24R/S25Y protein binds with fivefold greater affinity than wild-type a1 and is able to partially suppress defects in repression by alpha 2 mutants. As a result of these studies, we show how subtle differences in residues at a surface distant from the functional site code for a conformational switch that allows the a1 homeodomain to become active in DNA binding in association with its cofactor alpha 2.     

Stability of Ligand-induced Protein Conformation Influences Affinity in Maltose-binding Protein

Our understanding of what determines ligand affinity of proteins is poor, even with high-resolution structures available. Both the non-covalent ligand–protein interactions and the relative free energies of available conformations contribute to the affinity of a protein for a ligand. Distant, non-binding site residues can influence the ligand affinity by altering the free energy difference between a ligand-free and ligand-bound conformation. Our hypothesis is that when different ligands induce distinct ligand-bound conformations, it should be possible to tweak their affinities by changing the free energies of the available conformations. We tested this idea for the maltose-binding protein (MBP) from Escherichia coli. We used single-molecule Förster resonance energy transfer (smFRET) to distinguish several unique ligand-bound conformations of MBP. We engineered mutations, distant from the binding site, to affect the stabilities of different ligand-bound conformations. We show that ligand affinity can indeed be altered in a conformation-dependent manner. Our studies provide a framework for the tuning of ligand affinity, apart from modifying binding site residues.     

Beyond the binding site: the role of the β₂-β₃ loop and extra-domain structures in PDZ domains

A general paradigm to understand protein function is to look at properties of isolated well conserved domains, such as SH3 or PDZ domains. While common features of domain families are well understood, the role of subtle differences among members of these families is less clear. Here, molecular dynamics simulations indicate that the binding mechanism in PSD95-PDZ3 is critically regulated via interactions outside the canonical binding site, involving both the poorly conserved β?-β? loop and an extra-domain helix. Using the CRIPT peptide as a prototypical ligand, our simulations suggest that a network of salt-bridges between the ligand and this loop is necessary for binding. These contacts interconvert between each other on a time scale of a few tens of nanoseconds, making them elusive to X-ray crystallography. The loop is stabilized by an extra-domain helix. The latter influences the global dynamics of the domain, considerably increasing binding affinity. We found that two key contacts between the helix and the domain, one involving the β?-β? loop, provide an atomistic interpretation of the increased affinity. Our analysis indicates that both extra-domain segments and loosely conserved regions play critical roles in PDZ binding affinity and specificity.     

The effect of denaturants on protein thermal stability analyzed through a theoretical model considering multiple binding sites

A novel mathematical development applied to protein ligand binding thermodynamics is proposed, which allows the simulation, and therefore the analysis of the effects of multiple and independent binding sites to the Native and/or Unfolded protein conformations, with different binding constant values. Protein stability is affected when it binds to a small number of high affinity ligands or to a high number of low affinity ligands. Differential scanning calorimetry (DSC) measures released or absorbed energy of thermally induced structural transitions of biomolecules. This paper presents the general theoretical development for the analysis of thermograms of proteins obtained for n-ligands bound to the native protein and m-ligands bound to their unfolded form. In particular, the effect of ligands with low affinity and with a high number of binding sites (n and/or m > 50) is analyzed. If the interaction with the native form of the protein is the one that predominates, they are considered stabilizers and if the binding with the unfolded species predominates, it is expected a destabilizing effect. The formalism presented here can be adapted to fitting routines in order to simultaneously obtain the unfolding energy and ligand binding energy of the protein. The effect of guanidinium chloride on bovine serum albumin thermal stability, was successfully analyzed with the model considering low number of middle affinity binding sites to the native state and a high number of weak binding sites to the unfolded state.     

Characterization of an imidazoline/guanidinium receptive site distinct from the alpha 2-adrenergic receptor

alpha 2-Adrenergic receptors recognize a number of molecules with diverse chemical structures, including the yohimban diastereoisomers yohimbine and rauwolscine, catecholamines, guanidinium analogs, and imidazolines, such as clonidine. The affinity of the receptor protein for some of these ligands can vary by 10-100-fold among various tissues and species, suggesting a heterogeneous class of binding sites. Certain cellular effects elicited by the compounds possessing an imidazoline or guanidinium moiety may actually be mediated by a membrane receptor distinct from the alpha 2-adrenergic receptor. To determine whether this imidazoline/guanidinium receptive site (IGRS) and the alpha 2-adrenergic receptor represent distinct proteins, we solubilized and partially characterized the two binding sites in rabbit kidney. This tissue expresses both alpha 2-adrenergic receptors and high affinity imidazoline/guanidinium binding sites, the latter which are rauwolscine-insensitive but can be identified with the benzodioxan [3H]idazoxan. The IGRS and alpha 2-adrenergic receptor in rabbit kidney exhibit distinct ligand recognition properties, which are maintained after solubilization and partial purification. In addition, the two receptors can be physically separated by heparin-agarose or lectin affinity chromatography indicating that the two binding sites are distinct entities. [3H]Idazoxan binding is trypsin-sensitive, indicating that the IGRS is a protein rather than a lipid component of the plasma membrane. [3H]Idazoxan binding is not inhibited by endogenous agonists for known neurotransmitter receptors. However, the IGRS does recognize clonidine-displacing substance, a small non-catechol compound isolated from calf brain, suggesting the existence of a previously uncharacterized hormonal/neurotransmitter receptor system.     

Ligand binding properties of the human erythropoietin receptor extracellular domain expressed in Escherichia coli.

We developed an assay to directly measure the ligand binding properties of the cloned human erythropoietin receptor (EpoR). The cDNA encoding the extracellular domain of the human EpoR was amplified by polymerase chain reaction and ligated into the prokaryotic expression vector pGEX3X. Synthesis in Escherichia coli was induced and a soluble glutathione S-transferase fusion protein, EREx, was purified by erythropoietin affinity chromatography. Purified EREx was bound to GSH agarose beads and used in a solid phase ligand binding assay. Specific binding of 125I-erythropoietin to EREx beads was demonstrated. A single affinity class (Kd = 1.5 nM) of the binding site was evident on Scatchard analysis. The Kd of this site is quantitatively equivalent to that of the "low" affinity cellular binding site. Kinetic analysis of ligand binding to EREx revealed both the on and off rates to be rapid, with t1/2 of 60 and 40 s, respectively. EREx ligand binding exhibits no obvious metal ion dependence or cross-competition by other hemopoietins. Antibodies to EREx block the binding of erythropoietin to the cellular EpoR. We conclude that the 66-kDa EpoR protein is capable of specific ligand binding and that no covalent modifications or associated molecules are required for this interaction. We speculate that the "high" affinity cellular binding site (Kd less than 0.2 nM) results from the interaction of the EpoR with another molecule, either additional EpoR or associated subunits, that decreases the ligand off rate.     

Zinc coordination geometry and ligand binding affinity: the structural and kinetic analysis of the second-shell serine 228 residue and the methionine 180 residue of the aminopeptidase from Vibrio proteolyticus

The chemical properties of zinc make it an ideal metal to study the role of coordination strain in enzymatic rate enhancement. The zinc ion and the protein residues that are bound directly to the zinc ion represent a functional charge/dipole complex, and polarization of this complex, which translates to coordination distortion, may tune electrophilicity, and hence, reactivity. Conserved protein residues outside of the charge/dipole complex, such as second-shell residues, may play a role in supporting the electronic strain produced as a consequence of functional polarization. To test the correlation between charge/dipole polarity and ligand binding affinity, structure-function studies were carried out on the dizinc aminopeptidase from Vibrio proteolyticus. Alanine substitutions of S228 and M180 resulted in catalytically diminished enzymes whose crystal structures show very little change in the positions of the metal ions and the protein residues. However, more detailed inspections of the crystal structures show small positional changes that account for differences in the zinc ion coordination geometry. Measurements of the binding affinity of leucine phosphonic acid, a transition state analogue, and leucine, a product, show a correlation between coordination geometry and ligand binding affinity. These results suggest that the coordination number and polarity may tune the electrophilicity of zinc. This may have provided the evolving enzyme with the ability to discriminate between reaction coordinate species.     

NMR structure and dynamics of the chimeric protein SH3-F2

In order to further elucidate structural and dynamic principles of protein self-organization and protein-ligand interactions, a new chimeric protein was designed and a genetically engineered construct was created. SH3-F2 amino acid sequence consists of polyproline ligand mgAPPLPPYSA, GG linker, and the sequence of spectrin SH3 domain circular permutant S19-P20s. Structural and dynamic properties of the protein were studied with high-resolution NMR. According to NMR data, the tertiary structure of the chimeric protein SH3-F2 has a topology that is typical for SH3 domains in the complex with the ligand forming polyproline type II helix located in the conservative region of binding in the orientation II. The polyproline ligand closely adjoins with the protein globule and is stabilized by hydrophobic interactions. However, the interactions of the ligand and the part of globule related to SH3 domain is not too large, because the analysis of protein dynamical characteristics points to the low amplitude, high-frequency ligand tumbling relative to the slow intramolecular motions of the main globule. The constructed chimera allows carrying out further structural and thermodynamic investigations of polyproline helix properties and its interaction with regulatory domains.     

Protein stabilization and destabilization by guanidinium salts

Preferential interactions of bovine serum albumin were measured with guanidine sulfate, guanidine acetate, and guanidine hydrochloride. The results showed an increasing preferential hydration with increasing salt concentration for the sulfate, positive preferential salt binding for the hydrochloride, and an intermediate situation for the acetate. These results correlate well with the known effects of the three salts on protein stability, namely, the stabilizing effect of guanidine sulfate and the denaturing effect of guanidine hydrochloride. Comparison of guanidinium and magnesium salts indicated that the substitution of guanidinium ion for Mg2+ decreases the preferential hydration and increases the preferential salt binding, suggesting that the perturbation by guanidinium ion binding of the surface free energy is greater than that by Mg2+ ion. It was concluded that guanidine salts are not a special class, but their activity toward proteins is modulated by the same fine balance between hydration and salt binding to protein as in the case of other salts, with the second factor being stronger in guanidine salts.     

Manipulating cofactor binding thermodynamics in an artificial oxygen transport protein

We report the mutational analysis of an artificial oxygen transport protein, HP7, which operates via a mechanism akin to that of human neuroglobin and cytoglobin. This protein destabilizes one of two heme-ligating histidine residues by coupling histidine side chain ligation with the burial of three charged glutamate residues on the same helix. Replacement of these glutamate residues with alanine, which is uncharged, increases the affinity of the distal histidine ligand by a factor of 13. Paradoxically, it also decreases heme binding affinity by a factor of 5 in the reduced state and 60 in the oxidized state. Application of a three-state binding model, in which an initial pentacoordinate binding event is followed by a protein conformational change to hexacoordinate, provides insight into the mechanism of this seemingly counterintuitive result: the initial pentacoordinate encounter complex is significantly destabilized by the loss of the glutamate side chains, and the increased affinity for the distal histidine only partially compensates for that. These results point to the importance of considering each oxidation and conformational state in the design of functional artificial proteins.     

Extrusion of the C-terminal helix in navel orangeworm moth pheromone-binding protein (AtraPBP1) controls pheromone binding

The navel orangeworm, Amyelois transitella (Walker), is an agricultural insect pest that can be controlled by disrupting male–female communication with sex pheromones, a technique known as mating disruption. Insect pheromone-binding proteins (PBPs) provide fast transport of hydrophobic pheromones through aqueous sensillar lymph and promote sensitive delivery of pheromones to receptors. Here we present a mutational analysis on a PBP from A. transitella (AtraPBP1) to evaluate how the C-terminal helix in this protein controls pheromone binding as a function of pH. Pheromone binds tightly to AtraPBP1 at neutral pH, but the binding is much weaker at pH below 5. Deletion of the entire C-terminal helix (residues 129–142) causes more than 100-fold increase in pheromone-binding affinity at pH 5 and only a 1.5-fold increase at pH 7. A similar pH-dependent increase in pheromone binding is also seen for the H80A/H95A double mutant that promotes extrusion of the C-terminal helix by disabling salt bridges at each end of the helix. The single mutants (H80A and H95A) also exhibit pheromone binding at pH below 5, but with ∼2-fold weaker affinity. NMR and circular dichroism data demonstrate a large overall structural change in each of these mutants at pH 4.5, indicating an extrusion of the C-terminal helix that profoundly affects the overall structure of the low pH form. Our results confirm that sequestration of the C-terminal helix at low pH as seen in the recent NMR structure may serve to block pheromone binding. We propose that extrusion of these C-terminal residues at neutral pH (or by the mutations in this study) exposes a hydrophobic cleft that promotes high affinity pheromone binding.     

Localized influence of 2'-hydroxyl groups and helix geometry on protein recognition in the RNA major groove

The local geometry of a DNA helix can influence protein recognition, but the sequence-specific features that contribute to helix structure are not fully understood, and even less is known about how RNA helix geometry may affect protein recognition. To begin to understand how local or global helix structure may influence binding in an RNA model system, we generated a series of DNA analogues of HIV and BIV TAR RNAs in which ribose sugars were systematically substituted in and around the known binding sites for argininamide and a BIV Tat arginine-rich peptide, respectively, and measured their corresponding binding affinities. For each TAR interaction, binding occurs in the RNA major groove with high specificity, whereas binding to the all-DNA analogue is weak and nonspecific. Relatively few substitutions are needed to convert either DNA analogue of TAR into a high-affinity binder, with the ribose requirements being restricted largely to regions that directly contact the ligand. Substitutions at individual positions show up to 70-fold differences in binding affinity, even at adjacent base pairs, while two base pairs at the core of the BIV Tat peptide-RNA interface are largely unaffected by deoxyribose substitution. These results suggest that the helix geometries and unique conformational features required for binding are established locally and are relatively insulated from effects more than one base pair away. It seems plausible that arginine-rich peptides are able to adapt to a mosaic helical architecture in which segments as small as single base steps may be considered as modular recognition units.