HpTX2 is a toxin from the venom of Heteropoda venatoria spider that has been demonstrated to bind on Kv4.2 potassium channel. We have determined the solution structure of recombinant HpTX2 by use of conventional two-dimensional NMR techniques followed by distance-geometry and molecular dynamics. The calculated structure belongs to the Inhibitory Cystin Knot structural family that consists in a compact disulfide-bonded core, from which four loops emerge. A poorly defined two-stranded antiparallel β-sheet (residues 20–23 and 25–28) is detected. Analysis of the electrostatic charge anisotropy allows us to propose a functional map of HpTX2 different from the one described for κ-conotoxin PVIIA, but strongly related to the one of charybdotoxin. The orientation of the dipole moment of HpTX2 emerges through K27 which could therefore be the critical lysine residue. Close to this lysine are a second basic residue, R23, an aromatic cluster (F7, W25, W30) and an hydrophobic side chain (L24). The high density in aromatic side chains of the putative functional surface as well as the lack of an asparagine is proposed to be the structural basis of the specificity of HpTX2 toward Kv4.2 channel.

The ba3-type cytochrome c oxidase from Thermus thermophilus is known as a two subunit enzyme. Deduced from the crystal structure of this enzyme, we discovered the presence of an additional transmembrane helix “subunit IIa” spanning the membrane. The hydrophobic N-terminally blocked protein was isolated in high yield using high-performance liquid chromatography. Its complete amino acid sequence was determined by a combination of automated Edman degradation of both the deformylated and the cyanogen bromide cleaved protein and automated C-terminal sequencing of the native protein.

The molecular mass of 3,794 Da as determined by MALDI-MS and by ESI requires the N-terminal methionine to be formylated and is in good agreement with the value calculated from the formylmethionine containing sequence (3,766.5 Da + 28 Da = 3,794.5 Da). This subunit consits of 34 residues forming one helix across the membrane (Lys5–Ala34), which corresponds in space to the first transmembrane helix of subunit II of the cytochrome c oxidases from Paracoccus denitrificans and bovine heart, however, with opposite polarity. It is 35% identical to subunit IV of the ba3-cytochrome oxidase from Natronobacterium pharaonis.

The open reading frame encoding this new subunit IIa (cbaD) is located upstream of cbaB in the same operon as the genes for subunit I (cbaA) and subunit II (cbaB).

We describe the design of Escherichia coli cells that synthesize a structurally perfect, recombinant cytochrome c from the Thermus thermophilus cytochrome c552 gene. Key features are (1) construction of a plasmid-borne, chimeric cycA gene encoding an Escherichia coli-compatible, N-terminal signal sequence (MetLysIleSerIleTyrAlaThrLeuAlaAlaLeuSerLeuAlaLeuProAlaGlyAla) followed by the amino acid sequence of mature Thermus cytochrome c552; and (2) coexpression of the chimeric cycA gene with plasmid-borne, host-specific cytochrome c maturation genes (ccmABCDEFGH). Approximately 1 mg of purified protein is obtained from 1 L of culture medium. The recombinant protein, cytochrome rsC552, and native cytochrome c552 have identical redox potentials and are equally active as electron transfer substrates toward cytochrome ba3, a Thermus heme-copper oxidase. Native and recombinant cytochromes c were compared and found to be identical using circular dichroism, optical absorption, resonance Raman, and 500 MHz 1H-NMR spectroscopies. The 1.7 Å resolution X-ray crystallographic structure of the recombinant protein was determined and is indistinguishable from that reported for the native protein (Than, ME, Hof P, Huber R, Bourenkov GP, Bartunik HD, Buse G, Soulimane T, 1997, J Mol Biol 271:629–644). This approach may be generally useful for expression of alien cytochrome c genes in E. coli.

The crystal structure of the photoprotein obelin (22.2 kDa) from Obelia longissima has been determined and refined to 1.7 Å resolution. Contrary to the prediction of a peroxide, the noncovalently bound substrate, coelenterazine, has only a single oxygen atom bound at the C2-position. The protein-coelenterazine 2-oxy complex observed in the crystals is photo-active because, in the presence of calcium ion, bioluminescence emission within the crystal is observed. This structure represents only the second de novo protein structure determined using the anomalous scattering signal of the sulfur substructure in the crystal. The method used here is theoretically different from that used for crambin in 1981 (4.72 kDa) and represents a significant advancement in protein crystal structure determination.

Calbindin D28k is a highly conserved Ca2+-binding protein abundant in brain and sensory neurons. The 261-residue protein contains six EF-hands packed into one globular domain. In this study, we have reconstituted calbindin D28k from two fragments containing three EF-hands each (residues 1–132 and 133–261, respectively), and from other combinations of small and large fragments. Complex formation is studied by ion-exchange and size-exclusion chromatography, electrophoresis, surface plasmon resonance, as well as circular dichroism (CD), fluorescence, and NMR spectroscopy. Similar chromatographic behavior to the native protein is observed for reconstituted complexes formed by mixing different sets of complementary fragments, produced by introducing a cut between EF-hands 1, 2, 3, or 4. The C-terminal half (residues 133–261) appears to have a lower intrinsic stability compared to the N-terminal half (residues 1–132). In the presence of Ca2+, NMR spectroscopy reveals a high degree of structural similarity between the intact protein and the protein reconstituted from the 1–132 and 133–261 fragments. The affinity between these two fragments is 2 × 107 M−1, with association and dissociation rate constants of 2.7 × 104 M−1 s−1 and 1.4 × 10−3 s−1, respectively. The complex formed in the presence of Ca2+ is remarkably stable towards unfolding by urea and heat. Both the complex and intact protein display cold and heat denaturation, although residual α-helical structure is seen in the urea denatured state at high temperature. In the absence of Ca2+, the fragments do not recombine to yield a complex resembling the intact apo protein. Thus, calbindin D28k is an example of a protein that can only be reconstituted in the presence of bound ligand. The α-helical CD signal is increased by 26% after addition of Ca2+ to each half of the protein. This suggests that Ca2+-induced folding of the fragments is important for successful reconstitution of calbindin D28k.

Chaperonins cpn60/cpn10 (GroEL/GroES in Escherichia coli) assist folding of nonnative polypeptides. Folding of the chaperonins themselves is distinct in that it entails assembly of a sevenfold symmetrical structure. We have characterized denaturation and renaturation of the recombinant human chaperonin 10 (cpn10), which forms a heptamer. Denaturation induced by chemical denaturants urea and guanidine hydrochloride (GuHCl) as well as by heat was monitored by tyrosine fluorescence, far-ultraviolet circular dichroism, and cross-linking; all denaturation reactions were reversible. GuHCl-induced denaturation was found to be cpn10 concentration dependent, in accord with a native heptamer to denatured monomer transition. In contrast, urea-induced denaturation was not cpn10 concentration dependent, suggesting that under these conditions cpn10 heptamers denature without dissociation. There were no indications of equilibrium intermediates, such as folded monomers, in either denaturant. The different cpn10 denatured states observed in high [GuHCl] and high [urea] were supported by cross-linking experiments. Thermal denaturation revealed that monomer and heptamer reactions display the same enthalpy change (per monomer), whereas the entropy-increase is significantly larger for the heptamer. A thermodynamic cycle for oligomeric cpn10, combining chemical denaturation with the dissociation constant in absence of denaturant, shows that dissociated monomers are only marginally stable (3 kJ/mol). The thermodynamics for co-chaperonin stability appears conserved; therefore, instability of the monomer could be necessary to specify the native heptameric structure.

Peptide GFSKAELAKARAAKRGGY folds in an α-helical conformation that is stabilized by formation of a hydrophobic staple motif and an N-terminal capping box (Munoz V, Blanco FJ, Serrano L, 1995, Struct Biol 2:380–385). To investigate backbone and side-chain internal motions within the helix and hydrophobic staple, residues F2, A5, L7, A8, and A10 were selectively 13C- and 15N-enriched and NMR relaxation experiments were performed in water and in water/trifluoroethanol (TFE) solution at four Larmor frequencies (62.5, 125, 150, and 200 MHz for 13C). Relaxation data were analyzed using the model free approach and an anisotropic diffusion model. In water, angular variances of motional vectors range from 10 to 20° and backbone φ,ψ bond rotations for helix residues A5, L7, A8, and A10 are correlated indicating the presence of Cα-H, Cα-Cβ, and N-H rocking-type motions along the helix dipole axis. L7 side-chain CβH2 and CγH motions are also correlated and as motionally restricted as backbone CαH, suggesting considerable steric hindrance with neighboring groups. In TFE which stabilizes the fold, internal motional amplitudes are attenuated and rotational correlations are increased. For the side chain of hydrophobic staple residue F2, wobbling-in-a-cone type motions dominate in water, whereas in TFE, the Cβ-Cγ bond and phenyl ring fluctuate more simply about the Cα-Cβ bond. These data support the Daragan–Mayo model of correlated bond rotations (Daragan VA, Mayo KH, 1996, J Phys Chem 100:8378–8388) and contribute to a general understanding of internal motions in peptides and proteins.

The DNA-binding domain of the yeast heat shock transcription factor (HSF) contains a strictly conserved proline that is at the center of a kink. To define the role of this conserved proline-centered kink, we replaced the proline with a number of other residues. These substitutions did not diminish the ability of the full-length protein to support growth of yeast or to activate transcription, suggesting that the proline at the center of the kink is not conserved for function. The stability of the isolated mutant DNA-binding domains was unaltered from the wild-type, so the proline is not conserved to maintain the stability of the protein. The crystal structures of two of the mutant DNA-binding domains revealed that the helices in the mutant proteins were still kinked after substitution of the proline, suggesting that the proline does not cause the α-helical kink. So why are prolines conserved in this and the majority of other kinked α-helices if not for structure, function, or stability? The mutant DNA-binding domains are less soluble than wild-type when overexpressed. In addition, the folding kinetics, as measured by stopped-flow fluorescence, is faster for the mutant proteins. These two results support the premise that the presence of the proline is critical for the folding pathway of HSF's DNA-binding domain. The finding may also be more general and explain why kinked helices maintain their prolines.

The formation of the N-terminal β-hairpin of ubiquitin is thought to be an early event in the folding of this small protein. Previously, we have shown that a peptide corresponding to residues 1–17 of ubiquitin folds autonomously and is likely to have a native-like hairpin register. To investigate the causes of the stability of this fold, we have made mutations in the amino acids at the apex of the turn. We find that in a peptide where Thr9 is replaced by Asp, U(1–17)T9D, the native conformation is stabilized with respect to the wild-type sequence, so much so that we are able to characterize the structure of the mutant peptide fully by NMR spectroscopy. The data indicate that U(1–17)T9D peptide does indeed form a hairpin with a native-like register and a type I turn with a G1 β-bulge, as in the full-length protein. The reason for the greater stability of the U(1–17)T9D mutant remains uncertain, but there are nuclear Overhauser effects between the side chains of Asp9 and Lys11, which may indicate that a charge–charge interaction between these residues is responsible.

Previous conformational analysis of 10-residue linear peptides enabled us to identify some cross-strand side-chain interactions that stabilize β-hairpin conformations. The stabilizing influence of these interactions appeared to be greatly reduced when the interaction was located at the N- and C-termini of these 10-residue peptides. To investigate the effect of the position relative to the turn of favorable interactions on β-hairpin formation, we have designed two 15-residue β-hairpin forming peptides with the same residue composition and differing only in the location of two residues within the strand region. The conformational properties of these two peptides in aqueous solution were studied by 1H and 13C NMR. Differences in the conformational behavior of the two designed 15-residue peptides suggest that the influence of stabilizing factors for β-hairpin formation, in particular, cross-strand side-chain interactions, depends on their proximity to the turn. Residues adjacent to the turn are most efficient in that concern. This result agrees with the proposal that the turn region acts as the driving force in β-hairpin folding.

The Escherichia coli RecA protein triggers both DNA repair and mutagenesis in a process known as the SOS response. The 81-residue E. coli protein DinI inhibits activity of RecA in vivo. The solution structure of DinI has been determined by multidimensional triple resonance NMR spectroscopy, using restraints derived from two sets of residual dipolar couplings, obtained in bicelle and phage media, supplemented with J couplings and a moderate number of NOE restraints. DinI has an α/β fold comprised of a three-stranded β-sheet and two α-helices. The β-sheet topology is unusual: the central strand is flanked by a parallel and an antiparallel strand and the sheet is remarkably flat. The structure of DinI shows that six negatively charged Glu and Asp residues on DinI's kinked C-terminal α-helix form an extended, negatively charged ridge. We propose that this ridge mimics the electrostatic character of the DNA phospodiester backbone, thereby enabling DinI to compete with single-stranded DNA for RecA binding. Biochemical data confirm that DinI is able to displace ssDNA from RecA.

The SH3 domain, comprised of approximately 60 residues, is found within a wide variety of proteins, and is a mediator of protein–protein interactions. Due to the large number of SH3 domain sequences and structures in the databases, this domain provides one of the best available systems for the examination of sequence and structural conservation within a protein family. In this study, a large and diverse alignment of SH3 domain sequences was constructed, and the pattern of conservation within this alignment was compared to conserved structural features, as deduced from analysis of eighteen different SH3 domain structures. Seventeen SH3 domain structures solved in the presence of bound peptide were also examined to identify positions that are consistently most important in mediating the peptide-binding function of this domain. Although residues at the two most conserved positions in the alignment are directly involved in peptide binding, residues at most other conserved positions play structural roles, such as stabilizing turns or comprising the hydrophobic core. Surprisingly, several highly conserved side-chain to main-chain hydrogen bonds were observed in the functionally crucial RT-Src loop between residues with little direct involvement in peptide binding. These hydrogen bonds may be important for maintaining this region in the precise conformation necessary for specific peptide recognition. In addition, a previously unrecognized yet highly conserved β-bulge was identified in the second β-strand of the domain, which appears to provide a necessary kink in this strand, allowing it to hydrogen bond to both sheets comprising the fold.

We develop a protocol for estimating the free energy difference between different conformations of the same polypeptide chain. The conformational free energy evaluation combines the CHARMM force field with a continuum treatment of the solvent. In almost all cases studied, experimentally determined structures are predicted to be more stable than misfolded “decoys.” This is due in part to the fact that the Coulomb energy of the native protein is consistently lower than that of the decoys. The solvation free energy generally favors the decoys, although the total electrostatic free energy (sum of Coulomb and solvation terms) favors the native structure. The behavior of the solvation free energy is somewhat counterintuitive and, surprisingly, is not correlated with differences in the burial of polar area between native structures and decoys. Rather, the effect is due to a more favorable charge distribution in the native protein, which, as is discussed, will tend to decrease its interaction with the solvent. Our results thus suggest, in keeping with a number of recent studies, that electrostatic interactions may play an important role in determining the native topology of a folded protein. On this basis, a simplified scoring function is derived that combines a Coulomb term with a hydrophobic contact term. This function performs as well as the more complete free energy evaluation in distinguishing the native structure from misfolded decoys. Its computational efficiency suggests that it can be used in protein structure prediction applications, and that it provides a physically well-defined alternative to statistically derived scoring functions.

We have modeled the structure of human lymphotactin (hLpnt), by homology modeling and molecular dynamics simulations. This chemokine is unique in having a single disulfide bond and a long C-terminal tail. Because other structural classes of chemokines have two pairs of Cys residues, compared to one in Lpnt, and because it has been shown that both disulfide bonds are required for stability and function, the question arises how the Lpnt maintains its structural integrity. The initial structure of hLpnt was constructed by homology modeling. The first 63 residues in the monomer of hLpnt were modeled using the structure of the human CC chemokine, RANTES, whose sequence appeared most similar. The structure of the long C-terminal tail, missing in RANTES, was taken from the human muscle fatty-acid binding protein. In a Protein Data Bank search, this protein was found to contain a sequence that was most homologous to the long tail. Consequently, the modeled hLpnt C-terminal tail consisted of both α-helical and β-motifs. The complete model of the hLpnt monomer consisted of two α-helices located above the five-stranded β-sheet. Molecular dynamics simulations of the solvated initial model have indicated that the stability of the predicted fold is related to the geometry of Pro78. The five-stranded β-sheet appeared to be preserved only when Pro78 was modeled in the cis conformation. Simulations were also performed both for the C-terminal truncated forms of the hLpnt that contained one or two (CC chemokine-like) disulfide bonds, and for the chicken Lpnt (cLpnt). Our MD simulations indicated that the turn region (T30–G34) in hLpnt is important for the interactions with the receptor, and that the long C-terminal region stabilizes both the turn (T30–G34) and the five-stranded β-sheet. The major conclusion from our theoretical studies is that the lack of one disulfide bond and the extension of the C-terminus in hLptn are mutually complementary. It is very likely that removal of two Cys residues sufficiently destabilizes the structure of a chemokine molecule, particularly the core β-sheet, to abolish its biological function. However, this situation is rectified by the long C-terminal segment. The role of this long region is most likely to stabilize the first β-turn region and α-helix H1, explaining how this chemokine can function with a single disulfide bond.

Metabotropic glutamate receptors (mGluRs) belong to the family 3 of G-protein-coupled receptors. On these proteins, agonist binding on the extracellular domain leads to conformational changes in the 7-transmembrane domains required for G-protein activation. To elucidate the structural features that might be responsible for such an activation mechanism, we have generated models of the amino terminal domain (ATD) of type 4 mGluR (mGlu4R). The fold recognition search allowed the identification of three hits with a low sequence identity, but with high secondary structure conservation: leucine isoleucine valine-binding protein (LIVBP) and leucine-binding protein (LBP) as already known, and acetamide-binding protein (AmiC). These proteins are characterized by a bilobate structure in an open state for LIVBP/LBP and a closed state for AmiC, with ligand binding in the cleft. Models for both open and closed forms of mGlu4R ATD have been generated. ACPT-I (1-aminocyclopentane 1,3,4-tricarboxylic acid), a selective agonist, has been docked in the two models. In the open form, ACPT-I is only bound to lobe I through interactions with Lys74, Arg78, Ser159, and Thr182. In the closed form, ACPT-I is trapped between both lobes with additional binding to Tyr230, Asp312, Ser313, and Lys317 from lobe II. These results support the hypothesis that mGluR agonists bind a closed form of the ATDs, suggesting that such a conformation of the binding domain corresponds to the active conformation.

Presentation of peptides derived from endogenous proteins by class I major histocompatibility complex molecules is essential both for immunological self-tolerance and induction of cytotoxic T-cell responses against intracellular parasites. Despite frequent and diverse post-translational modification of eukaryotic cell proteins, very few class I-bound peptides with post-translationally modified residues are known. Here we describe a natural dodecamer ligand of HLA-B39 (B*3910) derived from an RNA-binding nucleoprotein that carried NG,NG-dimethyl-Arg. Although common among RNA-binding proteins, this modification was not previously known among natural class I ligands. The sequence of this peptide was determined by Edman degradation and electrospray ion trap mass spectrometry. The fragmentation pattern of the dimethyl-Arg side chain observed with this latter technique allowed us to unambiguously assign the isomeric form of the modified residue. The post-translationally modified ligand was a prominent component (1–2%) of the B*3910-bound peptide repertoire. The dimethyl-Arg residue was located in a central position of the peptide, amenable to interacting with T-cell receptors, and most other residues in the middle region of the peptide were Gly. These structural features strongly suggest that the post-translationally modified residue may have a major influence on the antigenic properties of this natural ligand.

The three-dimensional structures of cystatins, and other evidence, suggest that the flexible N-terminal region of these inhibitors may bind to target proteinases independent of the two rigid hairpin loops forming the remainder of the inhibitory surface. In an attempt to demonstrate such two-step binding, which could not be identified in previous kinetics studies, we introduced a cysteine residue before the N-terminus of cystatin A and labeled this residue with fluorescent probes. Binding of AANS- and AEDANS-labeled cystatin A to papain resulted in ∼4-fold and 1.2-fold increases of probe fluorescence, respectively, reflecting the interaction of the N-terminal region with the enzyme. Observed pseudo-first-order rate constants, measured by the loss of papain activity in the presence of a fluorogenic substrate, for the reaction of the enzyme with excess AANS-cystatin A increased linearly with the concentration of the latter. In contrast, pseudo-first-order rate constants, obtained from measurements of the change of probe fluorescence with either excess enzyme or labeled inhibitor, showed an identical hyperbolic dependence on the concentration of the reactant in excess. This dependence demonstrates that the binding occurs in two steps, and implies that the labeled N-terminal region of cystatin A interacts with the proteinase in the second step, subsequent to the hairpin loops. The comparable affinities and dissociation rate constants for the binding of labeled and unlabeled cystatin A to papain indicate that the label did not appreciably perturb the interaction, and that unlabeled cystatin therefore also binds in a similar two-step manner. Such independent binding of the N-terminal regions of cystatins to target proteinases after the hairpin loops may be characteristic of most cystatin–proteinase reactions.

The reaction mechanism of phosphoryl transfer catalyzed by UMP/CMP-kinase from Dictyostelium discoideum was investigated by semiempirical AM1 molecular orbital computations of an active site model system derived from crystal structures that contain a transition state analog or a bisubstrate inhibitor. The computational results suggest that the nucleoside monophosphate must be protonated for the forward reaction while it is unprotonated in the presence of aluminium fluoride, a popular transition state analog for phosphoryl transfer reactions. Furthermore, a compactification of the active site model system during the reaction and for the corresponding complex containing AlF3 was observed. For the active site residues that are part of the LID domain, conformational flexibility during the reaction proved to be crucial. On the basis of the calculations, a concerted phosphoryl transfer mechanism is suggested that involves the synchronous shift of a proton from the monophosphate to the transferred PO3-group. The proposed mechanism is thus analogous to the phosphoryl transfer mechanism in cAMP-dependent protein kinase that phosphorylates the hydroxyl groups of serine residues.

The Q7K/L33I/L63I HIV-1 protease mutant was expressed in Escherichia coli and the effect of binding a substrate-analog inhibitor, acetyl-pepstatin, was investigated by fluorescence spectroscopy and molecular dynamics. The dimeric enzyme has four intrinsic tryptophans, located at positions 6 and 42 in each monomer. Fluorescence spectra and acrylamide quenching experiments show two differently accessible Trp populations in the apoenzyme with kq1 = 6.85 × 109 M−1 s−1 and kq2 = 1.88 × 109 M−1 s−1, that merge into one in the complex with kq = 1.78 × 109 M−1 s−1.

500 ps trajectory analysis of Trp χ1/χ2 rotameric interconversions suggest a model to account for the observed Trp fluorescence. In the simulations, Trp6/Trp6B rotameric interconversions do not occur on this timescale for both HIV forms. In the apoenzyme simulations, however, both Trp42s and Trp42Bs are flipping between χ1/χ2 states; in the complexed form, no such interconverions occur. A detailed investigation of the local Trp environments sampled during the molecular dynamics simulation suggests that one of the apoenzyme Trp42B rotameric interconversions would allow indole-quencher contact, such as with nearby Tyr59. This could account for the short lifetime component. The model thus interprets the experimental data on the basis of the conformational fluctuations of Trp42s alone. It suggests that the rotameric interconversions of these Trps, located relatively far from the active site and at the very start of the flap region, becomes restrained when the apoenzyme binds the inhibitor. The model is thus consistent with associating components of the fluorescence decay in HIV-1 protease to ground state conformational heterogeneity.

Conformational changes are essential for the activity of many proteins. If, or how fast, internal fluctuations are related to slow conformational changes that mediate protein function is not understood. In this study, we measure internal fluctuations of the transport protein lactose permease in the presence and absence of substrate by tryptophan fluorescence spectroscopy. We demonstrate that nanosecond fluctuations of α-helices are enhanced when the enzyme transports substrate. This correlates with previously published kinetic data from transport measurements showing that millisecond conformational transitions of the substrate-loaded carrier are faster than those in the absence of substrate. These findings corroborate the hypothesis of the hierarchical model of protein dynamics that predicts that slow conformational transitions are based on fast, thermally activated internal motions.