Acetylcholinesterase (AChE) is an enzyme broadly distributed in many species, including parasites. It occurs in multiple molecular forms that differ in their quaternary structure and mode of anchoring to the cell surface. This review summarizes biochemical and immunological investigations carried out in our laboratories on AChE of the helmint, Schistosoma mansoni. AChE appears in S. mansoni in two principal molecular forms, both globular, with sedimentation coefficients of ∼6.5 and 8 S. On the basis of their substrate specificity and sensitivity to inhibitors, both are “true” acetylcholinesterases. Approximately half of the AChE activity of S. mansoni is located on the outer surface of the parasite, attached to the tegumental membrane via a covalently attached glycosylphosphatidylinositol anchor. The remainder is located within the parasite, mainly associated with muscle tissue. Whereas the internal enzyme is most likely involved in termination of neurotransmission at cholinergic synapses, the role of the surface enzyme remains to be established; there are, however, indications that it is involved in signal transduction. The two forms of AChE differ in their heparin-binding properties, only the internal 8 S form of the AChE being retained on a heparin column. The two forms differ also in their immunological specificity, since they are selectively recognized by different monoclonal antibodies. Polyclonal antibodies raised against S. mansoni AChE purified by affinity chromatography are specific for the parasite AChE, reacting with both molecular forms, but do not recognize AChE from other species. They interact with the surface-localized enzyme on the intact organism, and produce almost total complement-dependent killing of the parasite. S. mansoni AChE is thus demonstrated to be a functional protein, involved in multifaceted activities, which can serve as a suitable candidate for diagnostic purposes, vaccine development, and drug design.

Recent developments on virus-like particles have demonstrated their potential in transfecting eucaryotic cells. In the case of particles based on the major coat protein VP1 of polyoma virus, transfection occurs via binding of VP1 to sialic acids. Since sialic acid is present on almost every eucaryotic cell line, this results in an unspecific cell targeting. Generation of a cell-type specificity of this system would imply the presentation of a new function on the surface of VP1. To analyze whether a new functional protein can be placed on VP1, we inserted dihydrofolate reductase from Escherichia coli as a model protein. The effect of such an insertion on both VP1 and the inserted protein was investigated, respectively. The function of VP1, like the formation of pentameric capsomers and its ability to assemble into capsids, was not influenced by the insertion. The inserted dihydrofolate reductase showed major changes when compared to the wild-type form. The thermal stability of the enzyme was dramatically reduced in the fusion protein; nevertheless, the dihydrofolate reductase proved to be a fully active enzyme with only slightly increased KM values for its substrates. This model system provides the basis for further modifications of the VP1 protein to achieve an altered surface of VP1 with new properties.

Protein phosphatase 1 (PP1) is widely distributed among tissues and species and acts as a regulator of many important cellular processes. By targeting the catalytic part of PP1 (PP1C) toward particular loci and substrates, regulatory subunits constitute key elements conferring specificity to the holoenzyme. Here, we report the identification of an (α/β)8-barrel-like structure within the N-ter stretch of the human PP1 regulatory subunit hGM, which is part of the family of diverse proteins associated with glycogen metabolism. Protein homology modeling gave rise to a three-dimensional (3D) model for the 381 N-ter residue stretch of hGM, based on sequence similarity with Streptomyces olivochromogenes xylose isomerase, identified by using FASTA. The alignment was subsequently extended by using hydrophobic cluster analysis. The homology-derived model includes the putative glycogen binding area located within the 142–230 domain of hGM as well as a structural characterization of the PP1C interacting domain (segment 51–67). Refinement of the latter by molecular dynamics afforded a topology that is in agreement with previous X-ray studies (Egloff et al., 1997). Finite difference Poisson–Boltzmann calculations performed on the interacting domains of PP1C and hGM confirm the complementarity of the local electrostatic potentials of the two partners. This work highlights the presence of a conserved fold among distant species (mammalian, Caenorhabditis elegans, yeast) and, thus, emphasizes the involvement of PP1 in crucial basic cellular functions.

Calerythrin is a 20 kDa calcium-binding protein isolated from gram-positive bacterium Saccharopolyspora erythraea. Based on amino acid sequence homology, it has been suggested that calerythrin belongs to the family of invertebrate sarcoplasmic EF-hand calcium-binding proteins (SCPs), and therefore it is expected to function as a calcium buffer. NMR spectroscopy was used to obtain structural information on the protein in solution. Backbone and side chain 1H, 13C, and 15N assignments were obtained from triple resonance experiments HNCACB, HN(CO)CACB, HNCO, CC(CO)NH, and [15N]-edited TOCSY, and HCCH-TOCSY. Secondary structure was determined by using secondary chemical shifts and characteristic NOEs. In addition, backbone N-H residual dipolar couplings were measured from a spin-state selective [1H, 15N] correlation spectrum acquired from a sample dissolved in a dilute liquid crystal. Four EF-hand motifs with characteristic helix-loop-helix patterns were observed. Three of these are typical calcium-binding EF-hands, whereas site 2 is an atypical nonbinding site. The global fold of calerythrin was assessed by dipolar couplings. Measured dipolar couplings were compared with values calculated from four crystal structures of proteins with sequence homology to calerythrin. These data allowed us to recognize an overall similarity between the folds of calerythrin and sarcoplasmic calcium-binding proteins from the sandworm Nereis diversicolor and the amphioxus Branchiostoma lanceolatum.

The neurotrophins are growth factors that are involved in the development and survival of neurons. Neurotrophin release by a target tissue results in neuron growth along the neurotrophin concentration gradient, culminating in the eventual innervation of the target tissue. These activities are mediated through trk cell surface receptors. We have determined the structures of the heterodimer formed between brain-derived neurotrophic factor (BDNF) and neurotrophin 4 (NT4), as well as the structure of homodimer of NT4. We also present the structure of the Neurotrophin 3 homodimer, which is refined to higher resolution than previously published. These structures provide the first views of the architecture of the NT4 protomer. Comparison of the surface of a model of the BDNF homodimer with the structures of the neurotrophin homodimers reveals common features that may be important in the binding between the neurotrophins and their receptors. In particular, there exists an analogous region on the surface of each neurotrophin that is likely to be involved in trk receptor binding. Variations in sequence on the periphery of this common region serve to confer trk receptor specificity.

We present direct evidence for a change in protein structural specificity due to hydrophobic core packing. High resolution structural analysis of a designed core variant of ubiquitin reveals that the protein is in slow exchange between two conformations. Examination of side-chain rotamers indicates that this dynamic response and the lower stability of the protein are coupled to greater strain and mobility in the core. The results suggest that manipulating the level of side-chain strain may be one way of fine tuning the stability and specificity of proteins.

α1-Microglobulin (α1m) is an electrophoretically heterogeneous plasma protein. It belongs to the lipocalin superfamily, a group of proteins with a three-dimensional (3D) structure that forms an internal hydrophobic ligand-binding pocket. α1m carries a covalently linked unidentified chromophore that gives the protein a characteristic brown color and extremely heterogeneous optical properties. Twenty-one different colored tryptic peptides corresponding to residues 88–94, 118–121, and 122–134 of human α1m were purified. In these peptides, the side chains of Lys92, Lys118, and Lys130 carried size heterogeneous, covalently attached, unidentified chromophores with molecular masses between 122 and 282 atomic mass units (amu). In addition, a previously unknown uncolored lipophilic 282 amu compound was found strongly, but noncovalently associated with the colored peptides. Uncolored tryptic peptides containing the same Lys residues were also purified. These peptides did not carry any additional mass (i.e., chromophore) suggesting that only a fraction of the Lys92, Lys118, and Lys130 are modified. The results can explain the size, charge, and optical heterogeneity of α1m. A 3D model of α1m, based on the structure of rat epididymal retinoic acid-binding protein (ERABP), suggests that Lys92, Lys118, and Lys130 are semiburied near the entrance of the lipocalin pocket. This was supported by the fluorescence spectra of α1m under native and denatured conditions, which indicated that the chromophores are buried, or semiburied, in the interior of the protein. In human plasma, approximately 50% of α1m is complex bound to IgA. Only the free α1m carried colored groups, whereas α1m linked to IgA was uncolored.

The binding free energies of four inhibitors to bovine β-trypsin are calculated. The inhibitors use either ornithine, lysine, or arginine to bind to the S1 specificity site. The electrostatic contribution to binding free energy is calculated by solving the finite difference Poisson–Boltzmann equation, the contribution of nonpolar interactions is calculated using a free energy-surface area relationship and the loss of conformational entropy is estimated both for trypsin and ligand side chains. Binding free energy values are of a reasonable magnitude and the relative affinity of the four inhibitors for trypsin is correctly predicted. Electrostatic interactions are found to oppose binding in all cases. However, in the case of ornithine- and lysine-based inhibitors, the salt bridge formed between their charged group and the partially buried carboxylate of Asp189 is found to stabilize the complex. Our analysis reveals how the molecular architecture of the trypsin binding site results in highly specific recognition of substrates and inhibitors. Specifically, partially burying Asp189 in the inhibitor-free enzyme decreases the penalty for desolvation of this group upon complexation. Water molecules trapped in the binding interface further stabilize the buried ion pair, resulting in a favorable electrostatic contribution of the ion pair formed with ornithine and lysine side chains. Moreover, all side chains that form the trypsin specificity site are partially buried, and hence, relatively immobile in the inhibitor-free state, thus reducing the entropic cost of complexation. The implications of the results for the general problem of recognition and binding are considered. A novel finding in this regard is that like charged molecules can have electrostatic contributions to binding that are more favorable than oppositely charged molecules due to enhanced interactions with the solvent in the highly charged complex that is formed.

Sequential 1H-NMR assignments of mouse [Cd7]-metallothionein-1 (MT1) have been carried out by standard homonuclear NMR methods and the use of an accordion–heteronuclear multiple quantum correlation (HMQC) experiment for establishing the metal, 113Cd2+, to cysteine connectivities. The three-dimensional structure was then calculated using the distance constraints from two-dimensional nuclear Overhauser effect (NOE) spectroscopy spectra and the Cys–Cd connectivities as input for a distance geometry-dynamical simulated annealing protocol in X-PLOR 3.851. Similar to the mammalian MT2 isoforms, the homologous primary structure of MT1 suggested two separate domains, each containing one metal cluster. Because there were no interdomain constraints, the structure calculation for the N-terminal β- and the C-terminal α-domain were carried out separately. The structures are based on 409 NMR constraints, consisting of 381 NOEs and 28 cysteine–metal connectivities. The only elements of regular secondary structure found were two short stretches of 310 helices along with some half-turns in the α-domain. Structural comparison with rat liver MT2 showed high similarity, with the β-domain structure in mouse MT1 showing evidence of increased flexibility compared to the same domain in MT2. The latter was reflected by the presence of fewer interresidue NOEs, no slowly exchanging backbone amide protons, and enhanced cadmium–cadmium exchange rates found in the β-domain of MT1.

Pyrazole and its 4-alkyl substituted derivatives are potent inhibitors for many alcohol dehydrogenases. However, the human σσ isoenzyme exhibits a 580-fold lower affinity for 4-methylpyrazole than does the human β1β1 isoenzyme, with which it shares 69% sequence identity. In this study, structural and kinetic studies were utilized in an effort to identify key structural features that affect the binding of 4-methylpyrazole in human alcohol dehydrogenase isoenzymes. We have extended the resolution of the human σσ alcohol dehydrogenase (ADH) isoenzyme to 2.5 Å resolution. Comparison of this structure to the human β1β1 isoenzyme structure indicated that the side-chain position for Met141 in σσ ADH might interfere with 4-methylpyrazole binding. Mutation of Met141 in σσ ADH to Leu (σ141L) lowers the Ki for 4-methylpyrazole from 350 to 10 μM, while having a much smaller effect on the Ki for pyrazole. Thus, the mutagenesis results show that the residue at position 141, which lines the substrate-binding pocket at a position close to the methyl group of 4-methylpyrazole, directly affects the binding of the inhibitor. To rule out nonspecific structural changes due to the mutation, the X-ray structure of the σ141L mutant enzyme was determined to 2.4 Å resolution. The three-dimensional structure of the mutant enzyme is identical to the wild-type enzyme, with the exception of the residue at position 141. Thus, the differences in 4-methylpyrazole binding between the mutant and wild-type σσ ADH isoenzymes can be completely ascribed to the local changes in the topology of the substrate binding site, and provides an explanation for the class-specific differences in 4-methylpyrazole binding to the human ADH isoenzymes.

Microcalorimetry has been used to measure the stabilities of mutational variants of yeast iso-1 cytochrome c in which F82 and L85 have been replaced by other hydrophobic amino acids. Specifically, F82 has been replaced by Y and L85 by A. The double mutant F82Y, L85A iso-1 has also been studied, and the mutational perturbations are compared to those for the two single mutants, F82Y iso-1 and L85A iso-1. Results are interpreted in terms of known crystallographic structures. The data show that (1) the destabilization of the mutant proteins is similar in magnitude to that which is theoretically predicted by the more obvious mutation-induced structural effects; (2) the free energy of destabilization of the double mutant, F82Y, L85A iso-1, is less than the sum of those of the two single mutants, almost certainly because, in the double mutant, the -OH group of Y82 is able to protrude into the cavity formed by the L85A substitution. The more favorable structural accommodation of the new -OH group in the double mutant leads to additional stability through (1) further decreases in the volumes of internal cavities and (2) formation of an extra protein–protein hydrogen bond.

The crystal structure of the S642A mutant of mitochondrial aconitase (mAc) with citrate bound has been determined at 1.8 Å resolution and 100 K to capture this binding mode of substrates to the native enzyme. The 2.0 Å resolution, 100 K crystal structure of the S642A mutant with isocitrate binding provides a control, showing that the Ser → Ala replacement does not alter the binding of substrates in the active site. The aconitase mechanism requires that the intermediate product, cis-aconitate, flip over by 180° about the Cα–Cβ double bond. Only one of these two alternative modes of binding, that of the isocitrate mode, has been previously visualized. Now, however, the structure revealing the citrate mode of binding provides direct support for the proposed enzyme mechanism.

The human fibrinogen γ-chain C-terminal segment functions as the platelet integrin binding site as well as the Factor XIIIa cross-linking substrate and thus plays an important role in blood clot formation and stabilization. The three-dimensional structure of this segment has been determined using carrier protein driven crystallization. The C-terminal segment, γ-(398–411), was attached to a linker sequence at the C-terminus of glutathione S-transferase and the structure of this fusion protein determined at 1.8 Å resolution. Functional studies of the chimeric protein demonstrate that the fibrinogen sequence in the presence of the carrier protein retains its specific functions as ligand for platelet integrin αIIb β3 (gpIIb/IIIa) and as a cross-linking substrate for Factor XIIIa. The structure obtained for the fibrinogen γ-chain segment is not affected by crystal packing and can provide the missing links to the recently reported model of cross-linked fibrin.

We have determined the three-dimensional structure of the potassium channel inhibitor HsTX1, using nuclear magnetic resonance and molecular modeling. This protein belongs to the scorpion short toxin family, which essentially contains potassium channel blockers of 29 to 39 amino acids and three disulfide bridges. It is highly active on voltage-gated Kv1.3 potassium channels. Furthermore, it has the particularity to possess a fourth disulfide bridge. We show that HsTX1 has a fold similar to that of the three-disulfide–bridged toxins and conserves the hydrophobic core found in the scorpion short toxins. Thus, the fourth bridge has no influence on the global conformation of HsTX1. Most residues spatially analogous to those interacting with voltage-gated potassium channels in the three-disulfide–bridged toxins are conserved in HsTX1. Thus, we propose that Tyr21, Lys23, Met25, and Asn26 are involved in the biological activity of HsTX1. As an additional positively charged residue is always spatially close to the aromatic residue in toxins blocking the voltage-gated potassium channels, and as previous mutagenesis experiments have shown the critical role played by the C-terminus in HsTX1, we suggest that Arg33 is also important for the activity of the four disulfide-bridged toxin. Docking calculations confirm that, if Lys23 and Met25 interact with the GYGDMH motif of Kv1.3, Arg33 can contact Asp386 and, thus, play the role of the additional positively charged residue of the toxin functional site. This original configuration of the binding site of HsTX1 for Kv1.3, if confirmed experimentally, offers new structural possibilities for the construction of a molecule blocking the voltage-gated potassium channels.

The monoclonal antibody 1696, directed against the HIV-1 protease, displays strong inhibitory effects toward the catalytic activity of the enzyme of both the HIV-1 and HIV-2 isolates. This antibody cross-reacts with peptides that include the N-terminus of the enzyme, a region that is well conserved in sequence among different viral strains and which, furthermore, is crucial for homodimerization to the active enzymatic form. This observation, as well as antigen-binding studies in the presence of an active site inhibitor, suggest that 1696 inhibits the HIV protease by destabilizing its active homodimeric form. To characterize further how the antibody 1696 inhibits the HIV-1 and HIV-2 proteases, we have solved the crystal structure of its Fab fragment by molecular replacement and refined it at 3.0 Å resolution. The antigen binding site has a deep cavity at its center, which is lined mainly by acidic and hydrophobic residues, and is large enough to accommodate several antigen residues. The structure of the Fab 1696 could form a starting basis for the design of alternative HIV protease-inhibiting molecules of broad specificity.

Direct thermodynamic and kinetic investigations of the binding of nucleotides to the nucleoside monophosphate (NMP) site of NMP kinases have not been possible so far because a spectroscopic probe was not available. By coupling a fluorescent N-methylanthraniloyl- (mant) group to the β-phosphate of CDP via a butyl linker, a CDP analogue [(Pβ)MABA-CDP] was obtained that still binds specifically to the NMP site of UmpKdicty, because the base and the ribose moieties, which are involved in specific interactions, are not modified. This allows the direct determination of binding constants for its substrates in competition experiments.

The motor domain regions of three novel members of the kinesin superfamily TLKIF1, TLKIFC, and TLBIMC were identified in a thermophilic fungus Thermomyces lanuginosus. Based on sequence similarity, they were classified as members of the known kinesin families Unc104/KIF1, KAR3, and BIMC. TLKIF1 was subsequently expressed in Escherichia coli. The expression level was high, and the protein was mostly soluble, easy to purify, and enzymatically active. TLKIF1 is a monomeric kinesin motor, which in a gliding motility assay displays a robust plus-directed microtubule movement up to 2 μm/s. The discovery of TLKIF1 also demonstrates that a family of kinesin motors not previously found in fungi may in fact be used in this group of organisms.

In the fold recognition approach to structure prediction, a sequence is tested for compatibility with an already known fold. For membrane proteins, however, few folds have been determined experimentally. Here the feasibility of computing the vast majority of likely membrane protein folds is tested. The results indicate that conformation space can be effectively sampled for small numbers of helices. The vast majority of potential monomeric membrane protein structures can be represented by about 30-folds for three helices, but increases exponentially to about 1,500,000 folds for seven helices. The generated folds could serve as templates for fold recognition or as starting points for conformational searches that are well distributed throughout conformation space.

The accelerated pace of genomic sequencing has increased the demand for structural models of gene products. Improved quantitative methods are needed to study the many systems (e.g., macromolecular assemblies) for which data are scarce. Here, we describe a new molecular dynamics method for protein structure determination and molecular modeling. An energy function, or database potential, is derived from distributions of interatomic distances obtained from a database of known structures. X-ray crystal structures are refined by molecular dynamics with the new energy function replacing the Van der Waals potential. Compared to standard methods, this method improved the atomic positions, interatomic distances, and side-chain dihedral angles of structures randomized to mimic the early stages of refinement. The greatest enhancement in side-chain placement was observed for groups that are characteristically buried. More accurate calculated model phases will follow from improved interatomic distances. Details usually seen only in high-resolution refinements were improved, as is shown by an R-factor analysis. The improvements were greatest when refinements were carried out using X-ray data truncated at 3.5 Å. The database potential should therefore be a valuable tool for determining X-ray structures, especially when only low-resolution data are available.

The reaction mechanism of the catalytic phosphoryl transfer of cAMP-dependent protein kinase (cAPK) was investigated by semi-empirical AM1 molecular orbital computations of an active site model system derived from the crystal structure of the catalytic subunit of the enzyme. The activation barrier is calculated as 20.7 kcal mol−1 and the reaction itself to be exothermic by 12.2 kcal mol−1. The active site residue Asp166, which was often proposed to act as a catalytic base, does not accept a proton in any of the reaction steps. Instead, the hydroxyl hydrogen of serine is shifted to the simultaneously transferred phosphate group of ATP. Although the calculated transition state geometry indicates an associative phosphoryl transfer, no concentration of negative charge is found. To study the influence of protein mutations on the reaction mechanism, we compared two-dimensional energy hypersurfaces of the protein kinase wild-type model and a corresponding mutant in which Asp166 was replaced by alanine. Surprisingly, they show similar energy profiles despite the experimentally known decrease of catalytic activity for corresponding mutants. Furthermore, a model structure was examined, where the charged NH3 group of Lys168 was replaced by a neutral methyl group. The energetic hypersurface of this hypothetical mutant shows two possible pathways for phosphoryl transfer, which both require significantly higher activation energies than the other systems investigated, while the energetic stabilization of the reaction product is similar in all systems. As the position of the amino acid side chains and the substrate peptide is virtually unchanged in all model systems, our results suggest that the exchange of Asp166 by other amino acid is less important to the phosphoryl transfer itself, but crucial to maintain the configuration of the active site in vivo. The positively charged side chain of Lys168, however, is necessary to stabilize the intermediate reaction states, particularly the side chain of the substrate peptide.