Comparison of the structures and the crystal contacts of trypanosomal triosephosphate isomerase in four different crystal forms

The Potential of Secondary Metabolites from Plants as Drugs or Leads against Protozoan Neglected Diseases-Part III: In-Silico Molecular Docking Investigations

Malaria, leishmaniasis, Chagas disease, and human African trypanosomiasis continue to cause considerable suffering and death in developing countries. Current treatment options for these parasitic protozoal diseases generally have severe side effects, may be ineffective or unavailable, and resistance is emerging. There is a constant need to discover new chemotherapeutic agents for these parasitic infections, and natural products continue to serve as a potential source. This review presents molecular docking studies of potential phytochemicals that target key protein targets in Leishmania spp., Trypanosoma spp., and Plasmodium spp.

Crystallization of dienelactone hydrolase in two space groups: structural changes caused by crystal packing

Dienelactone hydrolase (DLH) is a monomeric protein with a simple α/β-hydrolase fold structure. It readily crystallizes in space group P2₁2₁2₁ from either a phosphate or ammonium sulfate precipitation buffer. Here, the structure of DLH at 1.85 Å resolution crystallized in space group C2 with two molecules in the asymmetric unit is reported. When crystallized in space group P2₁2₁2₁ DLH has either phosphates or sulfates bound to the protein in crucial locations, one of which is located in the active site, preventing substrate/inhibitor binding. Another is located on the surface of the enzyme coordinated by side chains from two different molecules. Crystallization in space group C2 from a sodium citrate buffer results in new crystallographic protein-protein interfaces. The protein backbone is highly similar, but new crystal contacts cause changes in side-chain orientations and in loop positioning. In regions not involved in crystal contacts, there is little change in backbone or side-chain configuration. The flexibility of surface loops and the adaptability of side chains are important factors enabling DLH to adapt and form different crystal lattices.

Wildtype and engineered monomeric triosephosphate isomerase from Trypanosoma brucei: partitioning of reaction intermediates in D2O and activation by phosphite dianion

Product yields for the reactions of (R)-glyceraldehyde 3-phosphate (GAP) in D2O at pD 7.9 catalyzed by wildtype triosephosphate isomerase from Trypanosoma brucei brucei (Tbb TIM) and a monomeric variant (monoTIM) of this wildtype enzyme were determined by (1)H NMR spectroscopy and were compared with the yields determined in earlier work for the reactions catalyzed by TIM from rabbit and chicken muscle [O'Donoghue, A. C., Amyes, T. L., and Richard, J. P. (2005), Biochemistry 44, 2610 - 2621]. Three products were observed from the reactions catalyzed by TIM: dihydroxyacetone phosphate (DHAP) from isomerization with intramolecular transfer of hydrogen, d-DHAP from isomerization with incorporation of deuterium from D2O into C-1 of DHAP, and d-GAP from incorporation of deuterium from D2O into C-2 of GAP. The yield of DHAP formed by intramolecular transfer of hydrogen decreases from 49% for the muscle enzymes to 40% for wildtype Tbb TIM to 34% for monoTIM. There is no significant difference in the ratio of the yields of d-DHAP and d-GAP for wildtype TIM from muscle sources and Trypanosoma brucei brucei, but partitioning of the enediolate intermediate of the monoTIM reaction to form d-DHAP is less favorable ((k(C1))(D)/(k(C2))(D) = 1.1) than for the wildtype enzyme ((k(C1))(D)/(k(C2))(D) = 1.7). Product yields for the wildtype Tbb TIM and monoTIM-catalyzed reactions of glycolaldehyde labeled with carbon-13 at the carbonyl carbon ([1-(13)C]-GA) at pD 7.0 in the presence of phosphite dianion and in its absence were determined by (1)H NMR spectroscopy [Go, M. K., Amyes, T. L., and Richard, J. P. (2009) Biochemistry 48, 5769-5778]. There is no detectable difference in the yields of the products of wildtype muscle and Tbb TIM-catalyzed reactions of [1-(13)C]-GA in D2O. The kinetic parameters for phosphite dianion activation of the reactions of [1-(13)C]-GA catalyzed by wildtype Tbb TIM are similar to those reported for the enzyme from rabbit muscle [Amyes, T. L. and Richard, J. P. (2007) Biochemistry 46, 5841-5854], but there is no detectable dianion activation of the reaction catalyzed by monoTIM. The engineered disruption of subunit contacts at monoTIM causes movement of the essential side chains of Lys-13 and His-95 away from the catalytic active positions. We suggest that this places an increased demand that the intrinsic binding energy of phosphite dianion be utilized to drive the change in the conformation of monoTIM back to the active structure for wildtype TIM.

The Limit of Accuracy of Protein Modeling: Influence of Crystal Packing on Protein Structure

The size of the protein database (PDB) makes it now feasible to arrive at statistical conclusions regarding structural effects of crystal packing. These effects are relevant for setting upper practical limits of accuracy on protein modeling. Proteins whose crystals have more than one molecule in the asymmetric unit or whose structures were determined at least twice by X-ray crystallography were paired and their differences analyzed. We demonstrate a clear influence of crystal environment on protein structure, including backbone conformations, hinge-like motions and side-chain conformations. The positions of surface water molecules tend to be variable in different crystal environments while those of ligands are not. Structures determined by independent groups vary more than structures determined by the same authors. The use of different refinement methods is a major source for this effect. Our pair-wise analysis derives a practical limit to the accuracy of protein modeling. For different crystal forms, the limit of accuracy (C(alpha), root-mean-square deviation (RMSD)) is approximately 0.8A for the entire protein, which includes approximately 0.3A due to crystal packing. For organized secondary elements, the upper limit of C(alpha) RMSD is 0.5-0.6A while for loops or protein surface it reaches 1.0A. Twenty percent of exposed side- chains exhibit different chi(1+2) conformations with approximately half of the effect also resulting from crystal packing. A web based tool for analysis and graphic presentation of surface areas of crystal contacts is available (http://ligin.weizmann.ac.il/cryco).

Closed Conformation of the Active Site Loop of Rabbit Muscle Triosephosphate Isomerase in the Absence of Substrate: Evidence of Conformational Heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.

On the role of the crystal environment in determining protein side-chain conformations

The role of crystal packing in determining the observed conformations of amino acid side-chains in protein crystals is investigated by (1) analysis of a database of proteins that have been crystallized in different unit cells (space group or unit cell dimensions) and (2) theoretical predictions of side-chain conformations with the crystal environment explicitly represented. Both of these approaches indicate that the crystal environment plays an important role in determining the conformations of polar side-chains on the surfaces of proteins. Inclusion of the crystal environment permits a more sensitive measurement of the achievable accuracy of side-chain prediction programs, when validating against structures obtained by X-ray crystallography. Our side-chain prediction program uses an all-atom force field and a Generalized Born model of solvation and is thus capable of modeling simple packing effects (i.e. van der Waals interactions), electrostatic effects, and desolvation, which are all important mechanisms by which the crystal environment impacts observed side-chain conformations. Our results are also relevant to the understanding of changes in side-chain conformation that may result from ligand docking and protein-protein association, insofar as the results reveal how side-chain conformations change in response to their local environment.

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.

Effect of pH and salt bridges on structural assembly: molecular structures of the monomer and intertwined dimer of the Eps8 SH3 domain

The SH3 domain of Eps8 was previously found to form an intertwined, domain-swapped dimer. We report here a monomeric structure of the EPS8 SH3 domain obtained from crystals grown at low pH, as well as an improved domain-swapped dimer structure at 1.8 A resolution. In the domain-swapped dimer the asymmetric unit contains two "hybrid-monomers." In the low pH form there are two independently folded SH3 molecules per asymmetric unit. The formation of intermolecular salt bridges is thought to be the reason for the formation of the dimer. On the basis of the monomer SH3 structure, it is argued that Eps8 SH3 should, in principle, bind to peptides containing a PxxP motif. Recently it was reported that Eps8 SH3 binds to a peptide with a PxxDY motif. Because the "SH3 fold" is conserved, alternate binding sites may be possible for the PxxDY motif to bind. The strand exchange or domain swap occurs at the n-src loops because the n-src loops are flexible. The thermal b-factors also indicate the flexible nature of n-src loops and a possible handle for domain swap initiation. Despite the loop swapping, the typical SH3 fold in both forms is conserved structurally. The interface of the acidic form of SH3 is stabilized by a tetragonal network of water molecules above hydrophobic residues. The intertwined dimer interface is stabilized by hydrophobic and aromatic stacking interactions in the core and by hydrophilic interactions on the surface.

Backbone dipoles generate positive potentials in all proteins: origins and implications of the effect

Asymmetry in packing the peptide amide dipole results in larger positive than negative regions in proteins of all folding motifs. The average side chain potential in 305 proteins is 109 +/- 30 mV (2. 5 +/- 0.7 kcal/mol/e). Because the backbone has zero net charge, the non-zero potential is unexpected. The larger oxygen at the negative and smaller proton at the positive end of the amide dipole yield positive potentials because: 1) at allowed phi and psi angles residues come off the backbone into the positive end of their own amide dipole, avoiding the large oxygen; and 2) amide dipoles with their carbonyl oxygen surface exposed and amine proton buried make the protein interior more positive. Twice as many amides have their oxygens exposed than their amine protons. The distribution of acidic and basic residues shows the importance of the bias toward positive backbone potentials. Thirty percent of the Asp, Glu, Lys, and Arg are buried. Sixty percent of buried residues are acids, only 40% bases. The positive backbone potential stabilizes ionization of 20% of the acids by >3 pH units (-4.1 kcal/mol). Only 6.5% of the bases are equivalently stabilized by negative regions. The backbone stabilizes bound anions such as phosphates and rarely stabilizes bound cations.

Catalytic competence of O-acetylserine sulfhydrylase in the crystal probed by polarized absorption microspectrophotometry

The reactions of the pyridoxal 5'-phosphate-dependent enzyme O-acetylserine sulfhydrylase with the substrate O-acetyl-L-serine and substrate analogs have been investigated in the crystalline state by single-crystal polarized absorption microspectrophotometry. This approach has allowed us to examine the catalytic competence of the enzyme in different crystalline states, one of which was used to determine the three-dimensional structure; experimental conditions were defined for the accumulation of catalytic intermediates in the crystal suitable for crystallographic analyses.O-Acetyl-L-serine reacts with the enzyme in one of the crystal forms leading via a beta-elimination reaction to the accumulation of the alpha-aminoacrylate Schiff base, absorbing maximally at 320 and 470 nm, as in solution. The dissociation constant for the alpha-aminoacrylate Schiff base is in the millimolar range, 500-fold higher than in solution, suggesting that crystal lattice interactions may oppose functionally relevant conformational changes. The dissociation constant exhibits a bell-shaped dependence on pH centered at pH 7. At this pH the alpha-aminoacrylate species slowly decays with time (30% decrease in 24 hours). The alpha-aminoacrylate intermediate readily reacts with sodium azide, an analog of sulfide, the natural nucleophilic agent, to give a new amino acid and the native enzyme, indicating that the crystalline enzyme catalyzes the overall beta-replacement reaction as in solution. In other crystal forms, including that used for the X-ray investigation, O-acetyl-L-serine either has an even higher dissociation constant or causes crystal damage upon binding. When the crystalline enzyme reacts with either L-cysteine or L-serine, the external aldimine intermediate is formed. The dissociation constants for both substrate analogs are closer to those observed in solution and are modulated by pH as in solution. Findings demonstrate that O-acetylserine sulfhydrylase is catalytically competent in the crystal although some regions of the molecule, likely involved in an open-closed transition induced by O-acetyl-L-serine binding, may have a limited flexibility. The accumulation in the crystal of both the external aldimine and the alpha-aminoacrylate intermediate makes feasible their structural determination and, therefore, the elucidation of the catalytic pathway at the molecular level.

Three-dimensional structures of single-chain Fv-neuraminidase complexes

The structure of the complex between a recombinant single-chain Fv construct of antibody NC10 with a five-residue peptide linker between VH and VL (termed scFv(5)), and its antigen, tetrameric neuraminidase from influenza virus (NA), has been determined and refined at 2.5 A resolution. The antibody-antigen binding interface is very similar to that of a similar NC10 scFv-NA complex in which the scFv has a 15-residue peptide linker (scFv(15)), and the NC10 Fab-NA complex. However, scFv(5) and scFv(15) have different stoichiometries in solution. While scFv(15) is predominantly monomeric in solution, scFv(5) forms dimers exclusively, because the five-residue linker is not long enough to permit VH and VL domains from the same polypeptide associating and forming an antigen-binding site. Upon forming a complex with NA, scFv(15) forms a approximately 300 kDa complex corresponding to one NA tetramer binding four scFv(15) monomers, while scFv(5) forms a approximately 590 kDa complex, corresponding to two NA tetramers crosslinked by four bivalent scFv(5) dimers. However, the dimeric scFv(5) in the scFv(5)-NA crystals does not crosslink NA tetramers, and modelling studies indicate that it is not possible to pack four dimeric and simultaneously bivalent scFvs between the NA tetramers with only a five-residue linker between VH and VL. The inability arises from the exacting requirement to orient the two antigen-binding surfaces to bind the tetrameric NA antigen while avoiding steric clashes with NC10 scFv(5) dimers bound to other sites on the NA tetramer. The utility of bivalent or bifunctional scFvs with short linkers may therefore be restricted by the steric constraints imposed by binding multivalent antigens.

Extent and nature of contacts between protein molecules in crystal lattices and between subunits of protein oligomers

A survey was compiled of several characteristics of the intersubunit contacts in 58 oligomeric proteins, and of the intermolecular contracts in the lattice for 223 protein crystal structures. The total number of atoms in contact and the secondary structure elements involved are similar in the two types of interfaces. Crystal contact patches are frequently smaller than patches involved in oligomer interfaces. Crystal contacts result from more numerous interactions by polar residues, compared with a tendency toward nonpolar amino acids at oligomer interfaces. Arginine is the only amino acid prominent in both types of interfaces. Potentials of mean force for residue-residue contacts at both crystal and oligomer interfaces were derived from comparison of the number of observed residue-residue interactions with the number expected by mass action. They show that hydrophobic interactions at oligomer interfaces favor aromatic amino acids and methionine over aliphatic amino acids; and that crystal contacts form in such a way as to avoid inclusion of hydrophobic interactions. They also suggest that complex salt bridges with certain amino acid compositions might be important in oligomer formation. For a protein that is recalcitrant to crystallization, substitution of lysine residues with arginine or glutamine is a recommended strategy.

Prediction of protein side-chain rotamers from a backbone-dependent rotamer library: a new homology modeling tool

Modeling by homology is the most accurate computational method for translating an amino acid sequence into a protein structure. Homology modeling can be divided into two sub-problems, placing the polypeptide backbone and adding side-chains. We present a method for rapidly predicting the conformations of protein side-chains, starting from main-chain coordinates alone. The method involves using fewer than ten rotamers per residue from a backbone-dependent rotamer library and a search to remove steric conflicts. The method is initially tested on 299 high resolution crystal structures by rebuilding side-chains onto the experimentally determined backbone structures. A total of 77% of chi1 and 66% of chi(1 + 2) dihedral angles are predicted within 40 degrees of their crystal structure values. We then tested the method on the entire database of known structures in the Protein Data Bank. The predictive accuracy of the algorithm was strongly correlated with the resolution of the structures. In an effort to simulate a realistic homology modeling problem, 9424 homology models were created using three different modeling strategies. For prediction purposes, pairs of structures were identified which shared between 30% and 90% sequence identity. One strategy results in 82% of chi1 and 72% chi(1 + 2) dihedral angles predicted within 40 degrees of the target crystal structure values, suggesting that movements of the backbone associated with this degree of sequence identity are not large enough to disrupt the predictive ability of our method for non-native backbones. These results compared favorably with existing methods over a comprehensive data set.

Protein crystallization

Crystallization is necessary to obtain the three-dimensional structure of proteins and nucleic acids; it often represents the bottleneck in structure determination. Our understanding of crystallization mechanisms is still incomplete. In this review, we emphasize fundamental aspects of the crystallization process. Protein-protein contacts in crystals are complex, involving a delicate balance of specific and nonspecific interactions. Depending on solution conditions, these interactions can lead to nucleation of crystals or to amorphous aggregation; this stage of crystallization has been successfully studied by light scattering. Post-nucleation crystal growth may proceed by mechanisms involving crystal defects or two-dimensional nucleation, as observed by atomic force and interference microscopy. Cessation of growth has been observed but remains incompletely understood. Impurities may play important roles during all stages of crystallization. Phase diagrams can guide optimization of conditions for nucleation and subsequent crystal growth; a theoretical understanding relating these to the intermolecular interactions is beginning to develop.

The 1.5-A resolution crystal structure of bacterial luciferase in low salt conditions

Bacterial luciferase is a flavin monooxygenase that catalyzes the oxidation of a long-chain aldehyde and releases energy in the form of visible light. A new crystal form of luciferase cloned from Vibrio harveyi has been grown under low-salt concentrations, which diffract x-rays beyond 1.5-A resolution. The x-ray structure of bacterial luciferase has been refined to a conventional R-factor of 18.2% for all recorded synchrotron data between 30.0 and 1.50-A resolution. Bacterial luciferase is an alpha-beta heterodimer, and the individual subunits fold into a single domain (beta/alpha)8 barrel. The high resolution structure reveals a non-prolyl cis peptide bond that forms between Ala74 and Ala75 in the alpha subunit near the putative active site. This cis peptide bond may have functional significance for creating a cavity at the active site. Bacterial luciferase employs reduced flavin as a substrate rather than a cofactor. The structure presented was determined in the absence of substrates. A comparison of the structural similarities between luciferase and a nonfluorescent flavoprotein, which is expressed in the lux operon of one genus of bioluminescent bacteria, suggests that the two proteins originated from a common ancestor. However, the flavin binding sites of the nonfluorescent protein are likely not representative of the flavin binding site on luciferase. The structure presented here will furnish a detailed molecular model for all bacterial luciferases.

Two crystal structures of the leupeptin-trypsin complex

Three-dimensional structures of trypsin with the reversible inhibitor leupeptin have been determined in two different crystal forms. The first structure was determined at 1.7 A resolution with R-factor = 17.7% in the trigonal crystal space group P3(1)21, with unit cell dimensions of a = b = 55.62 A, c = 110.51 A. The second structure was determined at a resolution of 1.8 A with R-factor = 17.5% in the orthorhombic space group P2(1)2(1)2(1), with unit cell dimensions of a = 63.69 A, b = 69.37 A, c = 63.01 A. The overall protein structure is very similar in both crystal forms, with RMS difference for main-chain atoms of 0.27 A. The leupeptin backbone forms four hydrogen bonds with trypsin and a fifth hydrogen bond interaction is mediated by a water molecule. The aldehyde carbonyl of leupeptin forms a covalent bond of 1.42 A length with side-chain oxygen of Ser-195 in the active site. The reaction of trypsin with leupeptin proceeds through the formation of stable tetrahedral complex in which the hemiacetal oxygen atom is pointing out of the oxyanion hole and forming a hydrogen bond with His-57.