Decorsin is an antagonist of integrin αIIb β3 and a potent platelet aggregation inhibitor. A synthetic gene encoding decorsin, originally isolated from the leech Macrobdella decora, was designed, constructed, and expressed in Escherichia coli. The synthetic gene was fused to the stII signal sequence and expressed under the transcriptional control of the E. coli alkaline phosphatase promoter. The protein was purified by size-exclusion filtration of the periplasmic contents followed by reversed-phase high-performance liquid chromatography. Purified recombinant decorsin was found to be indistinguishable from leech-derived decorsin based on amino acid composition, mass spectral analysis, and biological activity assays. Complete sequential assignments of 1H and proton bound 13C resonances were established. Stereospecific assignments of 21 of 25 nondegenerate β-methylene groups were determined. The RGD adhesion site recognized by integrin receptors was found at the apex of a most exposed hairpin loop. The dynamic behavior of decorsin was analyzed using several independent NMR parameters. Although the loop containing the RGD sequence is the most flexible one in decorsin, the conformation of the RGD site itself is more restricted than in other proteins with similar activities.

The lipase from Pseudomonas cepacia represents a widely applied catalyst for highly enantioselective resolution of chiral secondary alcohols. While its stereopreference is determined predominantly by the substrate structure, stereoselectivity depends on atomic details of interactions between substrate and lipase. Thirty secondary alcohols with published E values using P. cepacia lipase in hydrolysis or esterification reactions were selected, and models of their octanoic acid esters were docked to the open conformation of P. cepacia lipase. The two enantiomers of 27 substrates bound preferentially in either of two binding modes: the fast-reacting enantiomer in a productive mode and the slow-reacting enantiomer in a nonproductive mode. Nonproductive mode of fast-reacting enantiomers was prohibited by repulsive interactions. For the slow-reacting enantiomers in the productive binding mode, the substrate pushes the active site histidine away from its proper orientation, and the distance d(HNε − Oalc) between the histidine side chain and the alcohol oxygen increases. d(HNε − Oalc) was correlated to experimentally observed enantioselectivity: in substrates for which P. cepacia lipase has high enantioselectivity (E > 100), d(HNε − Oalc) is >2.2 Å for slow-reacting enantiomers, thus preventing efficient catalysis of this enantiomer. In substrates of low enantioselectivity (E < 20), the distance d(HNε − Oalc) is less than 2.0 Å, and slow- and fast-reacting enantiomers are catalyzed at similar rates. For substrates of medium enantioselectivity (20 < E < 100), d(HNε − Oalc) is around 2.1 Å. This simple model can be applied to predict enantioselectivity of P. cepacia lipase toward a broad range of secondary alcohols.

The crystal structure of ligand-free tryptophanyl-tRNA synthetase (TrpRS) was solved at 2.9 Å using a combination of molecular replacement and maximum-entropy map/phase improvement. The dimeric structure (R = 23.7, Rfree = 26.2) is asymmetric, unlike that of the TrpRS tryptophanyl–5′AMP complex (TAM; Doublié S, Bricogne G, Gilmore CJ, Carter CW Jr, 1995, Structure 3:17–31). In agreement with small-angle solution X-ray scattering experiments, unliganded TrpRS has a conformation in which both monomers open, leaving only the tryptophan-binding regions of their active sites intact. The amino terminal αA-helix, TIGN, and KMSKS signature sequences, and the distal helical domain rotate as a single rigid body away from the dinucleotide-binding fold domain, opening the AMP binding site, seen in the TAM complex, into two halves. Comparison of side-chain packing in ligand-free TrpRS and the TAM complex, using identification of nonpolar nuclei (Ilyin VA, 1994, Protein Eng 7:1189–1195), shows that significant repacking occurs between three relatively stable core regions, one of which acts as a bearing between the other two. These domain rearrangements provide a new structural paradigm that is consistent in detail with the “induced-fit” mechanism proposed for TyrRS by Fersht et al. (Fersht AR, Knill-Jones JW, Beduelle H, Winter G, 1988, Biochemistry 27:1581–1587). Coupling of ATP binding determinants associated with the two catalytic signature sequences to the helical domain containing the presumptive anticodon-binding site provides a mechanism to coordinate active-site chemistry with relocation of the major tRNA binding determinants.

By using a novel consensus approach, we have previously managed to generate a fully synthetic phytase, consensus phytase-1, that was 15–26 °C more thermostable than the parent fungal phytases used in its design (Lehmann et al., 2000). We now sought to use the backbone of consensus phytase-1 and to modify its catalytic properties. This was done by replacing a considerable part of the active site (i.e., all the divergent residues) with the corresponding residues of Aspergillus niger NRRL 3135 phytase, which displays pronounced differences in specific activity, substrate specificity, and pH-activity profile. For the new protein termed consensus phytase-7, a major—although not complete—shift in catalytic properties was observed, demonstrating that rational transfer of favorable catalytic properties from one phytase to another is possible by using this approach. Although the exchange of the active site was associated with a 7.6 °C decrease in unfolding temperature (Tm) as measured by differential scanning calorimetry, consensus phytase-7 still was >7 °C more thermostable than all wild-type ascomycete phytases known to date. Thus, combination of the consensus approach with the selection of a “preferred” active site allows the design of a thermostabilized variant of an enzyme family of interest that (most closely) matches the most favorable catalytic properties found among its family members.

Aspartic proteinases share a conserved network of hydrogen bonds (termed “fireman's grip”), which involves the hydroxyl groups of two threonine residues in the active site Asp-Thr-Gly triplets (Thr26 in the case of human immunodeficiency virus type 1 (HIV-1) PR). In the case of retroviral proteinases (PRs), which are active as symmetrical homodimers, these interactions occur at the dimer interface. For a systematic analysis of the “fireman's grip,” Thr26 of HIV-1 PR was changed to either Ser, Cys, or Ala. The variant enzymes were tested for cleavage of HIV-1 derived peptide and polyprotein substrates. PR(T26S) and PR(T26C) showed similar or slightly reduced activity compared to wild-type HIV-1 PR, indicating that the sulfhydryl group of cysteine can substitute for the hydroxyl of the conserved threonine in this position. PR(T26A), which lacks the “fireman's grip” interaction, was virtually inactive and was monomeric in solution at conditions where wild-type PR exhibited a monomer–dimer equilibrium. All three mutations had little effect when introduced into only one chain of a linked dimer of HIV-1 PR. In this case, even changing both Thr residues to Ala yielded residual activity suggesting that the “fireman's grip” is not essential for activity but contributes significantly to dimer formation. Taken together, these results indicate that the “fireman's grip” is crucial for stabilization of the retroviral PR dimer and for overall stability of the enzyme.

The ligand-binding domain of the human low-density lipoprotein receptor consists of seven modules, each of 40–45 residues. In the presence of calcium, these modules adopt a common polypeptide fold with three conserved disulfide bonds. A concatemer of the first and second modules (LB1–2) folds efficiently in the presence of calcium ions, forming the same disulfide connectivities as in the isolated modules. The three-dimensional structure of LB1–2 has now been solved using two-dimensional 1H NMR spectroscopy and restrained molecular dynamics calculations. No intermodule nuclear Overhauser effects were observed, indicating the absence of persistent interaction between them. The near random-coil NH and Hα chemical shifts and the low φ and ψ angle order parameters of the four-residue linker suggest that it has considerable flexibility. The family of LB1–2 structures superimposed well over LB1 or LB2, but not over both modules simultaneously. LB1 and LB2 have a similar pattern of calcium ligands, but the orientations of the indole rings of the tryptophan residues W23 and W66 differ, with the latter limiting solvent access to the calcium ion. From these studies, it appears that although most of the modules in the ligand-binding region of the receptor are joined by short segments, these linkers may impart considerable flexibility on this region.

A direct conflict between the stabilization free energy parameters of cytochrome c determined by optical methods and by hydrogen exchange (HX) is quantitatively explained when the partially folded intermediates seen by HX are taken into account. The results support the previous HX measurements of intermediate populations, show how intermediates can elude the standard melting analysis, and illustrate how they confuse the analysis when they are significantly populated within the melting transition region.

Polymerization of metallothioneins is one of the usually encountered puzzles during the research process of metallothioneins' structure and function. Our work focuses on the cysteine independently occurred polymerization from metallothioneins monomers in different milieus, while it leaves out the aggregation caused by the oxidation of cysteine, because the latter circumstance is the result of purification lapsus. After the purification of metallothioneins monomers, a dynamic light-scattering technique is used to detect the polymerized states of rabbit liver metallothionein I and II in different buffers, which is the first systematical detection of polymerized states of metallothioneins in solutions. The effects of different compositions of each buffer are discussed in details. Steric complementarity, hydrophobic, and electrostatic interaction characteristics are studied, following the modeling of monomers and relevant polymers of rat metallothionein II, rabbit liver metallothionein I and II. These theoretical calculations are the first complete computer simulations on different factors affecting metallothioneins' polymerization. A molecular recognition mechanism of metallothioneins' polymerization in solutions is proposed on the bases of experimental results and theoretical calculations. Preliminary X-ray studies of two crystal forms of rabbit liver metallothionein II are compared with the crystal structure of rat metallothionein II, and the polymerized states in crystal packing are discussed with the knowledge of polymerization of metallothioneins in solutions. The hypothesis, which is consistent with theoretical calculations and experimental results, is expected to construct a connection between the biochemical characteristics and physiological functions of metallothioneins, and this research may give some enlightenment to the topics of protein polymerizations.

Blue copper proteins are type-I copper-containing redox proteins whose role is to shuttle electrons from an electron donor to an electron acceptor in bacteria and plants. A large amount of experimental data is available on blue copper proteins; however, their functional characterization is hindered by the complexity of redox processes in biological systems. We describe here the application of a semiquantitative method based on a comparative analysis of molecular interaction fields to gain insights into the recognition properties of blue copper proteins. Molecular electrostatic and hydrophobic potentials were computed and compared for a set of 33 experimentally-determined structures of proteins from seven blue copper subfamilies, and the results were quantified by means of similarity indices. The analysis provides a classification of the blue copper proteins and shows that (1) comparison of the molecular electrostatic potentials provides useful information complementary to that highlighted by sequence analysis; (2) similarities in recognition properties can be detected for proteins belonging to different subfamilies, such as amicyanins and pseudoazurins, that may be isofunctional proteins; (3) dissimilarities in interaction properties, consistent with experimentally different binding specificities, may be observed between proteins belonging to the same subfamily, such as cyanobacterial and eukaryotic plastocyanins; (4) proteins with low sequence identity, such as azurins and pseudoazurins, can have sufficient similarity to bind to similar electron donors and acceptors while having different binding specificity profiles.

We have reported variants of the M13 bacteriophage major coat protein (P8) that enable high copy display of monomeric and oligomeric proteins, such as human growth hormone and steptavidin, on the surface of phage particles (Sidhu SS, Weiss GA, Wells JA. 2000. High copy display of large proteins on phage for functional selections. J Mol Biol 296:487–495). Here, we explore how an optimized P8 variant (opti-P8) could evolve the ability to efficiently display a protein fused to its N-terminus. Reversion of individual opti-P8 residues back to the wild-type P8 residue identifies a limited set of hydrophobic residues responsible for the high copy protein display. These hydrophobic amino acids bracket a conserved hydrophobic face on the P8 alpha helix thought to be in contact with the phage coat. Mutations additively combine to promote high copy protein display, which was further enhanced by optimization of the linker between the phage coat and the fusion protein. These data are consistent with a model in which protein display-enhancing mutations allow for better packing of the fusion protein into the phage coat. The high tolerance for phage coat protein mutations observed here suggests that filamentous phage coat proteins could readily evolve new capabilities.

Inspection of high resolution three-dimensional (3D) structures from the protein database reveals an increasing number of cis-Xaa-Pro and cis-Xaa-Yaa peptide bonds. However, we are still far from being able to predict whether these bonds will remain cis upon single-site substitution of Pro or Yaa and/or cleavage of a peptide bond close to it in the sequence. We have chosen oxidized Escherichia coli thioredoxin (Trx), a member of the Trx superfamily with a single α/β domain and cis P76 to determine the effect of single-site substitution and/or cleavage on this isomer. Standard two-dimensional (2D) NMR analysis were performed on cleaved Trx (1–73/74–108) and its P76A variant. Analysis of the NOE connectivities indicates remarkable similarity between the secondary and supersecondary structure of the noncovalent complexes and Trx. Analysis of the 2D version of the HCCH-TOCSY and HMQC-NOESY-HMQC and 13C-filtered HMQC-NOESY spectra of cleaved Trx with uniformly 13C-labeled I75 and P76 shows surprising conservation of both cis P76 and packing of I75 against W31. A similar NMR analysis of its P76A variant provides no evidence for cis A76 and shows only subtle local changes in both the packing of I75 and the interstrand connectivities between its most protected hydrophobic strands (β2 and β4). Indeed, a molecular simulation model for the trans P76A variant of Trx shows only subtle local changes around the substitution site. In conclusion, cleavage of R73 is insufficient to provoke cis/trans isomerization of P76, but cleavage and single-site substitution (P76A) favors the trans isomer.

The protease domain of coagulation factor VIIa (FVIIa) is homologous to trypsin with a similar active site architecture. The catalytic function of FVIIa is regulated by allosteric modulations induced by binding of divalent metal ions and the cofactor tissue factor (TF). To further elucidate the mechanisms behind these transformations, the effects of Zn2+ binding to FVIIa in the free form and in complex with TF were investigated. Equilibrium dialysis suggested that two Zn2+ bind with high affinity to FVIIa outside the N-terminal γ-carboxyglutamic acid (Gla) domain. Binding of Zn2+ to FVIIa, which was influenced by the presence of Ca2+, resulted in decreased amidolytic activity and slightly reduced affinity for TF. After binding to TF, FVIIa was less susceptible to zinc inhibition. Alanine substitutions for either of two histidine residues unique for FVIIa, His216, and His257, produced FVIIa variants with decreased sensitivity to Zn2+ inhibition. A search for putative Zn2+ binding sites in the crystal structure of the FVIIa protease domain was performed by Grid calculations. We identified a pair of Zn2+ binding sites in the Glu210–Glu220 Ca2+ binding loop adjacent to the so-called activation domain canonical to serine proteases. Based on our results, we propose a model that describes the conformational changes underlying the Zn2+-mediated allosteric down-regulation of FVIIa's activity.

This paper presents an analytically tractable model that captures the most elementary aspect of the protein folding problem, namely that both the energy and the entropy decrease as a protein folds. In this model, the system diffuses within a sphere in the presence of an attractive spherically symmetric potential. The native state is represented by a small sphere in the center, and the remaining space is identified with unfolded states. The folding temperature, the time-dependence of the populations, and the relaxation rate are calculated, and the folding dynamics is analyzed for both golf-course and funnel-like energy landscapes. This simple model allows us to illustrate a surprising number of concepts including entropic barriers, transition states, funnels, and the origin of single exponential relaxation kinetics.

Distant homologies between proteins are often discovered only after three-dimensional structures of both proteins are solved. The sequence divergence for such proteins can be so large that simple comparison of their sequences fails to identify any similarity. New generation of sensitive alignment tools use averaged sequences of entire homologous families (profiles) to detect such homologies. Several algorithms, including the newest generation of BLAST algorithms and BASIC, an algorithm used in our group to assign fold predictions for proteins from several genomes, are compared to each other on the large set of structurally similar proteins with little sequence similarity. Proteins in the benchmark are classified according to the level of their similarity, which allows us to demonstrate that most of the improvement of the new algorithms is achieved for proteins with strong functional similarities, with almost no progress in recognizing distant fold similarities.

It is also shown that details of profile calculation strongly influence its sensitivity in recognizing distant homologies. The most important choice is how to include information from diverging members of the family, avoiding generating false predictions, while accounting for entire sequence divergence within a family. PSI-BLAST takes a conservative approach, deriving a profile from core members of the family, providing a solid improvement without almost any false predictions. BASIC strives for better sensitivity by increasing the weight of divergent family members and paying the price in lower reliability. A new FFAS algorithm introduced here uses a new procedure for profile generation that takes into account all the relations within the family and matches BASIC sensitivity with PSI-BLAST like reliability.

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.

We have crystallized Drosophila melanogaster acetylcholinesterase and solved the structure of the native enzyme and of its complexes with two potent reversible inhibitors, 1,2,3,4-tetrahydro-N-(phenylmethyl)-9-acridinamine and 1,2,3,4-tetrahydro-N-(3-iodophenyl-methyl)-9-acridinamine—all three at 2.7 Å resolution. The refined structure of D. melanogaster acetylcholinesterase is similar to that of vertebrate acetylcholinesterases, for example, human, mouse, and fish, in its overall fold, charge distribution, and deep active-site gorge, but some of the surface loops deviate by up to 8 Å from their position in the vertebrate structures, and the C-terminal helix is shifted substantially. The active-site gorge of the insect enzyme is significantly narrower than that of Torpedo californica AChE, and its trajectory is shifted several angstroms. The volume of the lower part of the gorge of the insect enzyme is ∼50% of that of the vertebrate enzyme. Upon binding of either of the two inhibitors, nine aromatic side chains within the active-site gorge change their conformation so as to interact with the inhibitors. Some differences in activity and specificity between the insect and vertebrate enzymes can be explained by comparison of their three-dimensional structures.

The crystal structures of four active site-directed thrombin inhibitors, 1–4, in a complex with human α-thrombin have been determined and refined at up to 2.0 Å resolution using X-ray crystallography. These compounds belong to a structurally novel family of inhibitors based on a 2,3-disubstituted benzo[b]thiophene structure. Compared to traditional active-site directed inhibitors, the X-ray crystal structures of these complexes reveal a novel binding mode. Unexpectedly, the lipophilic benzo[b]thiophene nucleus of the inhibitor appears to bind in the S1 specificity pocket. At the same time, the basic amine of the C-3 side chain of the inhibitor interacts with the mostly hydrophobic proximal, S2, and distal, S3, binding sites. The second, basic amine side chain at C-2 was found to point away from the active site, occupying a location between the S1 and S′1 sites. Together, the aromatic rings of the C-2 and C-3 side chains sandwich the indole ring of Trp60D contained in the thrombin S2 insertion loop defined by the sequence “Tyr-Pro-Pro-Trp.” [The thrombin residue numbering used in this study is equivalent to that reported for chymotrypsinogen (Hartley BS, Shotton DM, 1971, The enzymes, vol. 3. New York: Academic Press. pp 323–373).] In contrast to the binding mode of more classical thrombin inhibitors (d-Phe-Pro-Arg-H, NAPAP, Argatroban), this novel class of benzo[b]thiophene derivatives does not engage in hydrogen bond formation with Gly216 of the thrombin active site. A detailed analysis of the three-dimensional structures not only provides a clearer understanding of the interaction of these agents with thrombin, but forms a foundation for rational structure-based drug design. The use of the data from this study has led to the design of derivatives that are up to 2,900-fold more potent than the screening hit 1.

The active sites and substrate bindings of Rhizobium trifolii molonyl-CoA synthetase (MCS) catalyzing the malonyl-CoA formation from malonate and CoA have been determined based on NMR spectroscopy, site-directed mutagenesis, and comparative modeling methods. The MCS-bound conformation of malonyl-CoA was determined from two-dimensional–transferred nuclear Overhauser effect spectroscopy data. MCS protein folds into two structural domains and consists of 16 α-helices, 24 β-strands, and several long loops. The core active site was determined as a wide cleft close to the end of the small C-terminal domain. The catalytic substrate malonate is placed between ATP and His206 in the MCS enzyme, supporting His206 in its catalytic role as it generates reaction intermediate, malonyl-AMP. These findings are strongly supported by previous biochemical data, as well as by the site-directed mutagenesis data reported here. This structure reveals the biochemical role as well as the substrate specificity that conservative residues of adenylate-forming enzymes have.

Human serum albumin (HSA) interacts with a vast array of chemically diverse ligands at specific binding sites. To pinpoint the essential structural elements for the formation of the warfarin binding site on human serum albumin, a defined set of five recombinant proteins comprising combinations of domains and/or subdomains of the N-terminal part were prepared and characterized by biochemical standard procedures, tryptophanyl fluorescence, and circular dichroic measurements, indicating well-preserved secondary and tertiary structures. Affinity constants for binding to warfarin were estimated by fluorescence titration experiments and found to be highest for HSA-DOM I-II and HSA, followed by HSA-DOM IB-II, HSA-DOM II, and HSA-DOM I-IIA. In addition, ultraviolet difference spectroscopy and induced circular dichroism experiments were carried out to get an in depth understanding of the binding mechanism of warfarin to the fragments as stand-alone proteins. This systematic study indicates that the primary warfarin binding site is centered in subdomain IIA with indispensable structural contributions of subdomain IIB and domain I, while domain III is not involved in this binding site, underlining the great potential that lies in the use of combinations of recombinant fragments for the study and accurate localization of ligand binding sites on HSA.

Structural genomic projects envision almost routine protein structure determinations, which are currently imaginable only for small proteins with molecular weights below 25,000 Da. For larger proteins, structural insight can be obtained by breaking them into small segments of amino acid sequences that can fold into native structures, even when isolated from the rest of the protein. Such segments are autonomously folding units (AFU) and have sizes suitable for fast structural analyses. Here, we propose to expand an intuitive procedure often employed for identifying biologically important domains to an automatic method for detecting putative folded protein fragments. The procedure is based on the recognition that large proteins can be regarded as a combination of independent domains conserved among diverse organisms. We thus have developed a program that reorganizes the output of BLAST searches and detects regions with a large number of similar sequences. To automate the detection process, it is reduced to a simple geometrical problem of recognizing rectangular shaped elevations in a graph that plots the number of similar sequences at each residue of a query sequence. We used our program to quantitatively corroborate the premise that segments with conserved sequences correspond to domains that fold into native structures. We applied our program to a test data set composed of 99 amino acid sequences containing 150 segments with structures listed in the Protein Data Bank, and thus known to fold into native structures. Overall, the fragments identified by our program have an almost 50% probability of forming a native structure, and comparable results are observed with sequences containing domain linkers classified in SCOP. Furthermore, we verified that our program identifies AFU in libraries from various organisms, and we found a significant number of AFU candidates for structural analysis, covering an estimated 5 to 20% of the genomic databases. Altogether, these results argue that methods based on sequence similarity can be useful for dissecting large proteins into small autonomously folding domains, and such methods may provide an efficient support to structural genomics projects.