Mutation of conserved histidines alters tertiary structure and nanomechanics of consensus ankyrin repeats

The conserved TPLH tetrapeptide motif of ankyrin repeats (ARs) plays an important role in stabilizing AR proteins, and histidine (TPLH)-to-arginine (TPLR) mutations in this motif have been associated with a hereditary human anemia, spherocytosis. Here, we used a combination of atomic force microscopy-based single-molecule force spectroscopy and molecular dynamics simulations to examine the mechanical effects of His → Arg substitutions in TPLH motifs in a model AR protein, NI6C. Our molecular dynamics results show that the mutant protein is less mechanically stable than the WT protein. Our atomic force microscopy results indicate that the mechanical energy input necessary to fully unfold the mutant protein is only half of that necessary to unfold the WT protein (53 versus 106 kcal/mol). In addition, the ability of the mutant to generate refolding forces is also reduced. Moreover, the mutant protein subjected to cyclic stretch-relax measurements displays mechanical fatigue, which is absent in the WT protein. Taken together, these results indicate that the His → Arg substitutions in TPLH motifs compromise mechanical properties of ARs and suggest that the origin of hereditary spherocytosis may be related to mechanical failure of ARs.     

Folding kinetics of the cooperatively folded subdomain of the IκBα ankyrin repeat domain

The ankyrin repeat (AR) domain of IκBα consists of a cooperative folding unit of roughly four ARs (AR1-AR4) and of two weakly folded repeats (AR5 and AR6). The kinetic folding mechanism of the cooperative subdomain, IκBα_67-206, was analyzed using rapid mixing techniques. Despite its apparent architectural simplicity, IκBα_67-206 displays complex folding kinetics, with two sequential on-pathway high-energy intermediates. The effect of mutations to or away from the consensus sequences of ARs on folding behavior was analyzed, particularly the GXTPLHLA motif, which have not been examined in detail previously. Mutations toward the consensus generally resulted in an increase in folding stability, whereas mutations away from the consensus resulted in decreased overall stability. We determined the free energy change upon mutation for three sequential transition state ensembles along the folding route for 16 mutants. We show that folding initiates with the formation of the interface of the outer helices of AR3 and AR4, and then proceeds to consolidate structure in these repeats. Subsequently, AR1 and AR2 fold in a concerted way in a single kinetic step. We show that this mechanism is robust to the presence of AR5 and AR6 as they do not strongly affect the folding kinetics. Overall, the protein appears to fold on a rather smooth energy landscape, where the folding mechanism conforms a one-dimensional approximation. However, we note that the AR does not necessarily act as a single folding element.     

Dissection of protein-protein interaction and CDK4 inhibition in the oncogenic versus tumor suppressing functions of gankyrin and P16

Protein-protein interactions usually involve a large number of residues; thus it is difficult to elucidate functional and structural roles of specific residues located in the interface. This problem is particularly challenging for ankyrin repeat proteins (ARs), which consist of linear arrays of small repeating units and play critical roles in almost every life process via protein-protein interactions, because the residues involved are discontinuously dispersed in both the ARs and their partners. Our previous studies showed that while both specific CDK4 inhibitor p16INK4A (P16) and gankyrin bind to cyclin-dependent kinase 4 (CDK4) in similar fashion, only P16 inhibits the kinase activity of CDK4. While this could explain why P16 is a tumor suppressor and gankyrin is oncogenic, the structural basis of these contrasting properties was unknown. Here we show that a double mutant of gankyrin, L62H/I79D, inhibits the kinase activity of CDK4, similar to P16, and such CDK4-inhibtory activity is associated with the I79D but not L62H mutation. In addition, mutations at I79 and L62 bring about a moderate decrease in the stability of gankyrin. Further structural and biophysical analyses suggest that the substitution of Ile79 with Asp leads to local conformational changes in loops I-III of gankyrin. Taken together, our results allow the dissection of the "protein-protein binding" and "CDK4 inhibition" functions of P16, show that the difference between tumor suppressing and oncogenic functions of P16 and gankyrin, respectively, mainly resides in a single residue, and provide structural insight to the contrasting biological functions of the two AR proteins.     

Characterization and further stabilization of designed ankyrin repeat proteins by combining molecular dynamics simulations and experiments

Multiple molecular dynamics simulations with explicit solvent at room temperature and at 400 K were carried out to characterize designed ankyrin repeat (AR) proteins with full-consensus repeats. Using proteins with one to five repeats, the stability of the native structure was found to increase with the number of repeats. The C-terminal capping repeat, originating from the natural guanine-adenine-binding protein, was observed to denature first in almost all high-temperature simulations. Notably, a stable intermediate is found in experimental equilibrium unfolding studies of one of the simulated consensus proteins. On the basis of simulation results, this intermediate is interpreted to represent a conformation with a denatured C-terminal repeat. To validate this interpretation, constructs without C-terminal capping repeat were prepared and did not show this intermediate in equilibrium unfolding experiments. Conversely, the capping repeats were found to be essential for efficient folding in the cell and for avoiding aggregation, presumably because of their highly charged surface. To design a capping repeat conferring similar solubility properties yet even higher stability, eight point mutations adapting the C-cap to the consensus AR and adding a three-residue extension at the C-terminus were predicted in silico and validated experimentally. The in vitro full-consensus proteins were also compared with a previously published designed AR protein, E3_5, whose internal repeats show 80% identity in primary sequence. A detailed analysis of the simulations suggests that networks of salt bridges between β-hairpins, as well as additional interrepeat hydrogen bonds, contribute to the extraordinary stability of the full consensus.     

Contribution of hydrogen bonding to the conformational stability of ribonuclease T1.

For 30 years, the prevailing view has been that the hydrophobic effect contributes considerably more than hydrogen bonding to the conformational stability of globular proteins. The results and reasoning presented here suggest that hydrogen bonding and the hydrophobic effect make comparable contributions to the conformational stability of ribonuclease T1 (RNase T1). When RNase T1 folds, 86 intramolecular hydrogen bonds with an average length of 2.95 A are formed. Twelve mutants of RNase T1 [Tyr----Phe (5), Ser----Ala (3), and Asn----Ala (4)] have been prepared that remove 17 of the hydrogen bonds with an average length of 2.93 A. On the basis of urea and thermal unfolding studies of these mutants, the average decrease in conformational stability due to hydrogen bonding is 1.3 kcal/mol per hydrogen bond. This estimate is in good agreement with results from several related systems. Thus, we estimate that hydrogen bonding contributes about 110 kcal/mol to the conformational stability of RNase T1 and that this is comparable to the contribution of the hydrophobic effect. Accepting the idea that intramolecular hydrogen bonds contribute 1.3 +/- 0.6 kcal/mol to the stability of systems in an aqueous environment makes it easier to understand the stability of the "molten globule" states of proteins, and the alpha-helical conformations of small peptides.     

Molecular dynamics simulations reveal insights into key structural elements of adenosine receptors

The crystal structure of the human A(2A) adenosine receptor, a member of the G protein-coupled receptor (GPCR) family, is used as a starting point for the structural characterization of the conformational equilibrium around the inactive conformation of the human A(2) (A(2A) and A(2B)) adenosine receptors (ARs). A homology model of the closely related A(2B)AR is reported, and the two receptors were simulated in their apo form through all-atom molecular dynamics (MD) simulations. Different conditions were additionally explored in the A(2A)AR, including the protonation state of crucial histidines or the presence of the cocrystallized ligand. Our simulations reveal the role of several conserved residues in the ARs in the conformational equilibrium of the receptors. The "ionic lock" absent in the crystal structure of the inactive A(2A)AR is rapidly formed in the two simulated receptors, and a complex network of interacting residues is presented that further stabilizes this structural element. Notably, the observed rotameric transition of Trp6.48 ("toggle switch"), which is thought to initiate the activation process in GPCRs, is accompanied by a concerted rotation of the conserved residue of the A(2)ARs, His6.52. This new conformation is further stabilized in the two receptors under study by a novel interaction network involving residues in transmembrane (TM) helices TM5 (Asn5.42) and TM3 (Gln3.37), which resemble the conformational changes recently observed in the agonist-bound structure of β-adrenoreceptors. Finally, the interaction between Glu1.39 and His7.43, a pair of conserved residues in the family of ARs, is found to be weaker than previously thought, and the role of this interaction in the structure and dynamics of the receptor is thoroughly examined. All these findings suggest that, despite the commonalities with other GPCRs, the conformational equilibrium of ARs is also modulated by specific residues of the family.     

Molecular dynamics study of the stabilities of consensus designed ankyrin repeat proteins

Two designed ankyrin repeat (AR) proteins (E3_5 and E3_19) are high homologous (with about 87% sequence identity) and their crystal structures have a Calpha atom-positional root-mean-square difference of about 0.14 nm. However, it was found that E3_5 is considerably more stable than E3_19 in guanidinium hydrochloride and thermal denaturation experiments. With the goal of providing insights into the various factors contributing to the stabilities of the designed AR proteins and suggesting possible mutations to enhance their stabilities, homology modeling and molecular dynamics (MD) simulations with explicit solvent have been performed. Because the crystal structure of E3_19 was solved later than that of E3_5, a homology model of E3_19 based on the crystal structure of E3_5 was also used in the simulations. E3_5 shows a very stable trajectory in both crystal and solution simulations. In contrast, the C-terminal repeat of E3_19 unfolds in the simulations starting from either the modeled structure or the crystal structure, although it has a sequence identical to that of E3_5. A continuum electrostatic model was used to estimate the effect of single mutations on protein stability and to study the interaction between the internal ARs and the C-terminal capping AR. Mutations involving charged residues were found to have large effects on stability. Due to the difference in charge distribution in the internal ARs of E3_19 and E3_5, their interaction with the C-terminal capping AR is less favorable in E3_19. The simulation trajectories suggest that the stability of the designed AR proteins can be increased by optimizing the electrostatic interactions within and between the different repeats.     

A C-H triplebond O hydrogen bond stabilized polypeptide chain reversal motif at the C terminus of helices in proteins

The serendipitous observation of a C-H cdots, three dots, centered O hydrogen bond mediated polypeptide chain reversal in synthetic peptide helices has led to a search for the occurrence of a similar motif in protein structures. From a dataset of 634 proteins, 1304 helices terminating in a Schellman motif have been examined. The C-H triplebond O interaction between the T-4 C(alpha)H and T+1 Cz doublebond O group (C triplebond O< or =3.5A) becomes possible only when the T+1 residue adopts an extended beta conformation (T is defined as the helix terminating residue adopting an alpha(L) conformation). In all, 111 examples of this chain reversal motif have been identified and the compositional and conformational preferences at positions T-4, T, and T+1 determined. A marked preference for residues like Ser, Glu and Gln is observed at T-4 position with the motif being further stabilized by the formation of a side-chain-backbone O triplebond H-N hydrogen bond involving the side-chain of residue T-4 and the N-H group of residue T+3. In as many as 57 examples, the segment following the helix was extended with three to four successive residues in beta conformation. In a majority of these cases, the succeeding beta strand lies approximately antiparallel with the helix, suggesting that the backbone C-H triplebond O interactions may provide a means of registering helices and strands in an antiparallel orientation. Two examples were identified in which extended registry was detected with two sets of C-H cdots, three dots, centered O hydrogen bonds between (T-4) C(alpha)H triplebond O (T+1) and (T-8) C(alpha)H triplebondC doublebond O (T+3).     

A 1.3-A structure of zinc-bound N-terminal domain of calmodulin elucidates potential early ion-binding step

Calmodulin (CaM) is a 16.8-kDa calcium-binding protein involved in calcium-signal transduction. It is the canonical member of the EF-hand family of proteins, which are characterized by a helix-loop-helix calcium-binding motif. CaM is composed of N- and C-terminal globular domains (N-CaM and C-CaM), and within each domain there are two EF-hand motifs. Upon binding calcium, CaM undergoes a significant, global conformational change involving reorientation of the four helix bundles in each of its two domains. This conformational change upon ion binding is a key component of the signal transduction and regulatory roles of CaM, yet the precise nature of this transition is still unclear. Here, we present a 1.3-Å structure of zinc-bound N-terminal calmodulin (N-CaM) solved by single-wavelength anomalous diffraction phasing of a selenomethionyl N-CaM. Our zinc-bound N-CaM structure differs from previously reported CaM structures and resembles calcium-free apo-calmodulin (apo-CaM), despite the zinc binding to both EF-hand motifs. Structural comparison with calcium-free apo-CaM, calcium-loaded CaM, and a cross-linked calcium-loaded CaM suggests that our zinc-bound N-CaM reveals an intermediate step in the initiation of metal ion binding at the first EF-hand motif. Our data also suggest that metal ion coordination by two key residues in the first metal-binding site represents an initial step in the conformational transition induced by metal binding. This is followed by reordering of the N-terminal region of the helix exiting from this first binding loop. This conformational switch should be incorporated into models of either stepwise conformational transition or flexible, dynamic energetic state sampling-based transition.     

Conformational equilibria of valine studied by dynamics simulation.

The conformational probability distribution of a valine residue in the valine dipeptide and of the valine side chain in an alpha-helix, as well as the change in helix stability for replacing alanine with valine, has been calculated by molecular dynamics simulations of explicitly hydrated systems: dipeptide, tetrapeptide and 10-, 14- and 18-residue oligoalanine helices. All computed free-energy differences are means from at least eight separate slow-growth simulations, four in each direction and are reported with their root-mean-square deviations. Different values for the change in free energy of folding (delta delta G degrees) have been calculated with the use of forcefields having an all-atom and a central-atom representation of methyl groups, etc. The value obtained with the all-atom forcefield agrees well with new experimental values (3 kJ/mol = 0.7 kcal/mol). Furthermore, the most stable valine side-chain rotamer in the helix is different for these two representations. The most stable rotamer for the all atom conformation is the same one that predominates for valines in alpha-helices in proteins of known conformation. The lower conformational freedom of the valine side chain in the helix contributes 1 kJ/mol to the difference in stability computed with the all-atom potential; unfavorable interactions of the side chain with helix, even in the most stable conformation, further increase delta delta G degrees.     

Long time molecular dynamic simulation on the agonist binding and activation of the β2-adrenergic receptor

β₂-Adrenergic receptor (β₂AR) plays a critical role in mediating the effects of catecholamine hormones. Due to the flexibility of the structure of its active state, study of agonist–β2AR is usually performed by molecular dynamic (MD) simulation. In this study we show the representative characteristics of agonist binding and activation on β₂AR by MD simulation. The binding site and the conformational changes in the specific regions of β2AR are reasonable which confirmed the conclusion that agonist–β2AR reached its active-like state. We have studied the disruption of non-covalent intramolecular interactions, including the conserved DRY motif, the rotamer toggle switch and the ionic lock, the cytoplasmic ends of transmembranes 5 and 6, and some water-mediated hydrogen bond network regions. We conclude that agonist induced β₂AR to its active conformation, or at least the active-like conformation.     

CD-n.m.r. study of the solution conformation of bradykinin analogs containing alpha-aminoisobutyric acid

The conformation in aqueous solution of several alpha-aminoisobutyric acid (AIB)-containing analogs of bradykinin (BK) has been probed by complementary CD and 1H n.m.r. measurements. The conclusion reached is that substitution of AIB for Pro2 and/or Pro3 in BK stabilizes a degree of beta-turn conformation in the N-terminal tetrapeptide moiety of the resulting analogs. Changing the solvent from water to DMSO or TFE further enhances the contribution of particular hydrogen bonded structures to the time-averaged conformation of these peptides. Bradykinin and [AIB7]-BK adopt similar hydrogen bonded conformations in TFE, apparently with a contribution from a beta-turn involving their common Arg1-Pro2-Pro3-Gly4 moiety. The contrasting biological activities of BK and its AIB-analogs are considered in terms of the conformational analogy between the AIB-residue and cis' Pro and the propensity for a beta-turn at the N-terminus of the peptide.     

Molecular dynamics simulations to explore the active/inactive conformers of guinea pig β2 adrenoceptor for the selective design of agonists or antagonists

It is well known that guinea pig β₂ adrenoceptors (Gβ₂ARs) and human β₂ adrenoceptors (Hβ₂ARs) have structural similarity. However, only one conformational state of Gβ₂ARs has been studied – the putative inactive state. As adrenoceptors have a repertoire of conformations, and there is evidence that a certain conformation is stabilised as a ligand approaches, the aim of this study was to build four models of Gβ₂ARs by using putative active/inactive Hβ₂AR conformers as a template. We evaluated the accuracy of these models in regard to the binding mode and affinity values of a set of known β₂AR ligands through docking and molecular dynamics simulations. During docking simulations, ligands reached Gβ₂AR sites similar to those reported for Hβ₂ARs. The greatest differences between conformational states were found in the domains (TM5 and TM6) previously suggested as being key to ligand recognition. The coefficients of determination between experimental and calculated affinity values were near to but less than 0.66 in all cases. The highest values were for agonists on the active models and antagonists on the inactive model. The four Gβ₂AR models proved useful for analysing agonist/antagonist activity. The results suggest that the selection of an adequate model is dependent on the intrinsic activity of a given ligand.     

Designed armadillo repeat proteins as general peptide-binding scaffolds: consensus design and computational optimization of the hydrophobic core

Armadillo repeat proteins are abundant eukaryotic proteins involved in several cellular processes, including signaling, transport, and cytoskeletal regulation. They are characterized by an armadillo domain, composed of tandem armadillo repeats of approximately 42 amino acids, which mediates interactions with peptides or parts of proteins in extended conformation. The conserved binding mode of the peptide in extended form, observed for different targets, makes armadillo repeat proteins attractive candidates for the generation of modular peptide-binding scaffolds. Taking advantage of the large number of repeat sequences available, a consensus-based approach combined with a force field-based optimization of the hydrophobic core was used to derive soluble, highly expressed, stable, monomeric designed proteins with improved characteristics compared to natural armadillo proteins. These sequences constitute the starting point for the generation of designed armadillo repeat protein libraries for the selection of peptide binders, exploiting their modular structure and their conserved binding mode.     

Conformational energy analysis of the leucine repeat regions of C/EBP,GCN4, and the proteins of themyc,jun, andfos oncogenes

It has been recently proposed that certain DNA binding proteins (including C/EBP, GCN4 and themyc, jun, andfos oncogene proteins) share a common structural motif based on helix-promoting regions containing heptad repeat sequences of leucines. It has been suggested that this structure is critical to the biological activity of these proteins, since it facilitates the formation of functional dimers held together by interdigitating leucine side-chains along the hydrophobic interfaces between long α-helical regions of the polypeptide chains in a configuration termed the “leucine zipper.” In this paper, conformational energy analysis is used to determine the preferred three-dimensional structures of the leucine repeat regions of these proteins. The results indicate that, in all cases, the global minimum energy conformation for these regions is an amphipathic α-helix with the leucine side-chains arrayed on one side in such a way to favor “leucine zipper” dimerization. Furthermore, amino acid substitutions in these regions (such as Pro for Leu), that are known to inhibit dimer formation and prevent DNA binding, are found to produce significant conformational changes that disrupt the amphipathic helical structure. Thus, these results provide support for the proposed “leucine zipper” configuration as a critical structural feature of this class of DNA binding proteins.

     

Designing repeat proteins: well-expressed,soluble and stable proteins from combinatorial libraries of consensus ankyrin repeat proteins

We describe an efficient way to generate combinatorial libraries of stable, soluble and well-expressed ankyrin repeat (AR) proteins. Using a combination of sequence and structure consensus analyses, we designed a 33 amino acid residue AR module with seven randomized positions having a theoretical diversity of 7.2x10(7). Different numbers of this module were cloned between N and C-terminal capping repeats, i.e. ARs designed to shield the hydrophobic core of stacked AR modules. In this manner, combinatorial libraries of designed AR proteins consisting of four to six repeats were generated, thereby potentiating the theoretical diversity. All randomly chosen library members were expressed in soluble form in the cytoplasm of Escherichia coli in amounts up to 200 mg per 1 l of shake-flask culture. Virtually pure proteins were obtained in a single purification step. The designed AR proteins are monomeric and display CD spectra identical with those of natural AR proteins. At the same time, our AR proteins are highly thermostable, with T(m) values ranging from 66 degrees C to well above 85 degrees C. Thus, our combinatorial library members possess the properties required for biotechnological applications. Moreover, the favorable biophysical properties and the modularity of the AR fold may account, partly, for the abundance of natural AR proteins.     

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6 ± 0.3 kcal/mol per -CH₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH₂- group on folding contributes, on average, 1.1 ± 0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG_tr values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60 ± 4% and hydrogen bonds contribute 40 ± 4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.