Conserved water molecules in a large family of microbial ribonucleases.

We systematically analyzed the crystallographically determined water molecules of all known structures of RNase T1 and compared them to the ordered solvent in a large number of related microbial nucleases. To assess the crystallographers' impact on the interpretation of the solvent structure, we independently refined five validation structures from diffraction data derived from five isomorphous crystals of RNase T1. We also compared the positions of water molecules found in 11 published isomorphous RNase T1 inhibitor complexes. These data suggest that the positions of most of the waters located on the surface of a protein and that are well-determined in the experimental electron density maps are determined primarily by crystal packing forces. Water molecules with less well-defined electron density are in general unique to one or a small number of crystal structures. Only a small number of the well-defined waters are found to be independent of the crystal environment. These waters have a low accessible surface area and B-factor, and tend to be conserved in the crystal structures of a number of evolutionary related ribonucleases as well. A single water molecule is found conserved in all known microbial ribonucleases.     

Novel X-ray sequences and crystal structures of Persian and Starry sturgeon methemoglobins: Highlighting the role of heme pocket waters in causing autoxidation

Although there is a high sequence similarity between mammalian and fish hemoglobin (Hb), the oxidation and heme loss rates can vary greatly between them such that fish Hbs oxidise much more rapidly than mammalian Hbs. There is to date no sequence or structural data for any sturgeon Hb to reveal the level of autoxidation in these fish. In this study, novel high resolution X-ray sequences and crystal structures of methemoglobin (Met-Hb) from two sturgeon fish including Persian sturgeon (Acipenser percisus) and Starry sturgeon (Acipenser stellatus) belonging to the Caspian sea has been determined. A comprehensive sequence and structure comparison between these sturgeon Met-Hbs and a number of non-sturgeon and normal and sickle cell anaemia human Hb in varying heme states has been carried out highlighting (i) the structural variability in the heme propionate groups; (ii) the existence of certain residues or their displacement and shift in the heme pocket allowing entry of water molecules into the heme pocket; (iii) the importance of the number of water molecules in the heme pocket; (iv) the hydrogen bonding between oxygens of A and D propionate groups and that of waters in the heme pocket; and (v) the role of heme binding waters causing oxidative stress and heme autoxidation.     

Structure-based identification of binding sites,native ligands and potential inhibitors for G-protein coupled receptors

G-protein coupled receptors (GPCRs) are the largest family of cell-surface receptors involved in signal transmission. Drugs associated with GPCRs represent more than one fourth of the 100 top-selling drugs and are the targets of more than half of the current therapeutic agents on the market. Our methodology based on the internal coordinate mechanics (ICM) program can accurately identify the ligand-binding pocket in the currently available crystal structures of seven transmembrane (7TM) proteins [bacteriorhodopsin (BR) and bovine rhodopsin (bRho)]. The binding geometry of the ligand can be accurately predicted by ICM flexible docking with and without the loop regions, a useful finding for GPCR docking because the transmembrane regions are easier to model. We also demonstrate that the native ligand can be identified by flexible docking and scoring in 1.5% and 0.2% (for bRho and BR, respectively) of the best scoring compounds from two different types of compound database. The same procedure can be applied to the database of available chemicals to identify specific GPCR binders. Finally, we demonstrate that even if the sidechain positions in the bRho binding pocket are entirely wrong, their correct conformation can be fully restored with high accuracy (0.28 A) through the ICM global optimization with and without the ligand present. These binding site adjustments are critical for flexible docking of new ligands to known structures or for docking to GPCR homology models. The ICM docking method has the potential to be used to "de-orphanize" orphan GPCRs (oGPCRs) and to identify antagonists-agonists for GPCRs if an accurate model (experimentally and computationally validated) of the structure has been constructed or when future crystal structures are determined.     

Three-dimensional solution structure of the HIV-1 protease complexed with DMP323, a novel cyclic urea-type inhibitor, determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of the HIV-1 protease homodimer, MW 22.2 kDa, complexed to a potent, cyclic urea-based inhibitor, DMP323, is reported. This is the first solution structure of an HIV protease/inhibitor complex that has been elucidated. Multidimensional heteronuclear NMR spectra were used to assemble more than 4,200 distance and angle constraints. Using the constraints, together with a hybrid distance geometry/simulated annealing protocol, an ensemble of 28 NMR structures was calculated having no distance or angle violations greater than 0.3 A or 5 degrees, respectively. Neglecting residues in disordered loops, the RMS deviation (RMSD) for backbone atoms in the family of structures was 0.60 A relative to the average structure. The individual NMR structures had excellent covalent geometry and stereochemistry, as did the restrained minimized average structure. The latter structure is similar to the 1.8-A X-ray structure of the protease/DMP323 complex (Chang CH et al., 1995, Protein Science, submitted); the pairwise backbone RMSD calculated for the two structures is 1.22 A. As expected, the mismatch between the structures is greatest in the loops that are disordered and/or flexible. The flexibility of residues 37-42 and 50-51 may be important in facilitating substrate binding and product release, because these residues make up the respective hinges and tips of the protease flaps. Flexibility of residues 4-8 may play a role in protease regulation by facilitating autolysis.     

Crystal structure of the core domain of RhoE/Rnd3: a constitutively activated small G protein

We report the 2.1 A crystal structure of the core G protein domain of the unusual Rho family member RhoE/Rnd3 in complex with endogenous GTP and magnesium. Unlike other small G proteins, RhoE, along with two other proteins Rnd1/Rho6 and Rnd2/RhoN, does not hydrolyze GTP. The main reason for this is the presence of serines in the positions equivalent to Ala59 and Gln61 in Ras. The structure shows that there are still water molecules in similar positions to the waters thought to be involved in the hydrolysis reaction in other G proteins. The structure suggests three not necessarily exclusive explanations for the lack of hydrolysis. The lack of the conserved glutamine raises the energy of the transition state inhibiting hydrolysis. The serines may restrain the waters from moving closer to the GTP, a step that is required to attain the transition state. They also stabilize the GTP-bound conformation of switch II and could prevent conformational changes required during hydrolysis. By superposition of the RhoE structure on structures of Rho family proteins in complex with binding partners, we make predictions on RhoE interactions with these partners.     

Solvent structure in crystals of trypsin determined by X-ray and neutron diffraction.

The solvent structure in orthorhombic crystals of bovine trypsin has been independently determined by X-ray diffraction to 1.35 A resolution and by neutron diffraction to 2.1 A resolution. A consensus model of the water molecule positions was obtained using oxygen positions identified in the electron density map determined by X-ray diffraction, which were verified by comparison to D2O-H2O difference neutron scattering density. Six of 184 water molecules in the X-ray structure, all with B-factors greater than 50 A2, were found to be spurious after comparison with neutron results. Roughly two-thirds of the water of hydration expected from thermodynamic data for proteins was localized by neutron diffraction; approximately one-half of the water of hydration was located by X-ray diffraction. Polar regions of the protein are well hydrated, and significant D2O-H2O difference density is seen for a small number of water molecules in a second shell of hydration. Hydrogen bond lengths and angles calculated from unconstrained refinement of water positions are distributed about values typically seen in small molecule structures. Solvent models found in seven other bovine trypsin and trypsinogen and rat trypsin structures determined by X-ray diffraction were compared. Internal water molecules are well conserved in all trypsin structures including anionic rat trypsin, which is 65% homologous to bovine trypsin. Of the 22 conserved waters in trypsin, 19 were also found in trypsinogen, suggesting that they are located in regions of the apoprotein that are structurally conserved in the transition to the mature protein. Seven waters were displaced upon activation of trypsinogen. Water structure at crystal contacts is not generally conserved in different crystal forms. Three groups of integral structural water molecules are highly conserved in all solvent structures, including a spline of water molecules inserted between two beta-strands, which may resemble an intermediate in the formation of beta sheets during the folding of a protein.     

Novel DNA binding motifs in the DNA repair enzyme endonuclease III crystal structure.

The 1.85 A crystal structure of endonuclease III, combined with mutational analysis, suggests the structural basis for the DNA binding and catalytic activity of the enzyme. Helix-hairpin-helix (HhH) and [4Fe-4S] cluster loop (FCL) motifs, which we have named for their secondary structure, bracket the cleft separating the two alpha-helical domains of the enzyme. These two novel DNA binding motifs and the solvent-filled pocket in the cleft between them all lie within a positively charged and sequence-conserved surface region. Lys120 and Asp138, both shown by mutagenesis to be catalytically important, lie at the mouth of this pocket, suggesting that this pocket is part of the active site. The positions of the HhH motif and protruding FCL motif, which contains the DNA binding residue Lys191, can accommodate B-form DNA, with a flipped-out base bound within the active site pocket. The identification of HhH and FCL sequence patterns in other DNA binding proteins suggests that these motifs may be a recurrent structural theme for DNA binding proteins.     

Crystal structure of mouse carnitine octanoyltransferase and molecular determinants of substrate selectivity

Carnitine acyltransferases have crucial functions in fatty acid metabolism. Members of this enzyme family show distinctive substrate preferences for short-, medium- or long-chain fatty acids. The molecular mechanism for this substrate selectivity is not clear as so far only the structure of carnitine acetyltransferase has been determined. To further our understanding of these important enzymes, we report here the crystal structures at up to 2.0-A resolution of mouse carnitine octanoyltransferase alone and in complex with the substrate octanoylcarnitine. The structures reveal significant differences in the acyl group binding pocket between carnitine octanoyltransferase and carnitine acetyltransferase. Amino acid substitutions and structural changes produce a larger hydrophobic pocket that binds the octanoyl group in an extended conformation. Mutation of a single residue (Gly-553) in this pocket can change the substrate preference between short- and medium-chain acyl groups. The side chains of Cys-323 and Met-335 at the bottom of this pocket assume dual conformations in the substrate complex, and mutagenesis studies suggest that the Met-335 residue is important for catalysis.     

On the contribution of water-mediated interactions to protein-complex stability

Protein-water interactions have long been recognized as a major determinant of chain folding, conformational stability, binding specificity and catalysis. However, the detailed effects of water on stabilizing protein-protein interactions remain elusive. A way to test experimentally the contribution of water-mediated interactions is by applying double mutant cycle analysis on pairs of residues that do not form direct interactions, but are bridged by water. Seven such interactions within the interface between TEM1 and BLIP proteins were evaluated. No significant interaction free energy was found between either of them. Water can bridge interactions, but also stabilize the structure of the monomer. To distinguish between these, we performed a bioinformatic analysis using AQUAPROT (http://bioinfo.weizmann.ac.il/aquaprot) to determine the degree of water conservation between the bound and unbound states. 29 structures of twelve complexes and 20 related monomers were analyzed. Of the 262 water molecules located within the interfaces, 145 were conserved between the unbound and bound structures. Strikingly, all 50 buried or partially buried waters in the monomer structures were conserved at the same location in the bound structures. Thus, buried waters have an important role in stabilizing the monomer fold rather than contributing to protein-protein binding, and are not replaced by residues from the incoming protein. Taking together the experimental and bioinformatics evidence suggests that exposed waters within the interface may be good sites for protein engineering, while buried or mostly buried waters should be left unchanged.     

Structure of a second crystal form of Bence-Jones protein Loc: strikingly different domain associations in two crystal forms of a single protein

We have determined the structure of the immunoglobulin light-chain dimer Loc in a second crystal form that was grown from distilled water. The crystal structure was determined to 2.8-A resolution; the R factor is 0.22. The two variable domains are related by local 2-fold axes and form an antigen binding "pocket". The variable domain-variable domain interaction observed in this crystal form differs from the one exhibited by the protein when crystallized from ammonium sulfate in which the two variable domains formed a protrusion (Chang et al., 1985). The structure attained in the distilled water crystals is similar to, but not identical with, the one observed for the Mcg light-chain dimer in crystals grown from ammonium sulfate. Thus, two strikingly different structures were attained by this multisubunit protein in crystals grown under two different, commonly used, crystallization techniques. The quaternary interactions exhibited by the protein in the two crystal forms are sufficiently different to suggest fundamentally different interpretations of the structural basis for the function of this protein. This observation may have general implications regarding the use of single crystallographic determinations for detailed identification of structural and functional relationships. On the other hand, proteins whose structures can be altered by manipulation of crystallization conditions may provide useful systems for study of fundamental structural chemistry.     

The three-dimensional structure of recombinant bovine chymosin at 2.3 A resolution

The crystal structure of recombinant bovine chymosin (EC 3.4.23.4; renin), which was cloned and expressed in Escherichia coli, has been determined using X-ray data extending to 2.3 A resolution. The crystals of the enzyme used in this study belong to the space group I222 with unit cell dimensions alpha = 72.7 A, b = 80.3 A, and c = 114.8 A. The structure was solved by the molecular replacement method and was refined by a restrained least-squares procedure. The crystallographic R factor is 0.165 and the deviation of bond distances from ideality is 0.020 A. The resulting model includes all 323 amino acid residues, as well as 297 water molecules. The enzyme has an irregular shape with approximate maximum dimensions of 40 x 50 x 65 A. The secondary structure consists primarily of parallel and antiparallel beta-strands with a few short alpha-helices. The enzyme can be subdivided into N- and C-terminal domains which are separated by a deep cleft containing the active aspartate residues Asp-34 and Asp-216. The amino acid residues and waters at the active site form an extensive hydrogen-bonded network which maintains the pseudo 2-fold symmetry of the entire structure. A comparison of recombinant chymosin with other acid proteinases reveals the high degree of structural similarity with other members of this family of proteins as well as the subtle differences which make chymosin unique. In particular, Tyr-77 of the flap region of chymosin does not hydrogen bond to Trp-42 but protrudes out in the P1 pocket forming hydrophobic interactions with Phe-119 and Leu-32. This may have important implications concerning the mechanism of substrate binding and substrate specificity.     

A new protein binding pocket similarity measure based on comparison of clouds of atoms in 3D: application to ligand prediction

Background  

Predicting which molecules can bind to a given binding site of a protein with known 3D structure is important to decipher the protein function, and useful in drug design. A classical assumption in structural biology is that proteins with similar 3D structures have related molecular functions, and therefore may bind similar ligands. However, proteins that do not display any overall sequence or structure similarity may also bind similar ligands if they contain similar binding sites. Quantitatively assessing the similarity between binding sites may therefore be useful to propose new ligands for a given pocket, based on those known for similar pockets.     

Interactions of different inhibitors with active-site aspartyl residues of HIV-1 protease and possible relevance to pepsin

The importance of the active site region aspartyl residues 25 and 29 of the mature HIV-1 protease (PR) for the binding of five clinical and three experimental protease inhibitors [symmetric cyclic urea inhibitor DMP323, nonhydrolyzable substrate analog (RPB) and the generic aspartic protease inhibitor acetyl-pepstatin (Ac-PEP)] was assessed by differential scanning calorimetry. DeltaT(m) values, defined as the difference in T(m) for a given protein in the presence and absence of inhibitor, for PR with DRV, ATV, SQV, RTV, APV, DMP323, RPB, and Ac-PEP are 22.4, 20.8, 19.3, 15.6, 14.3, 14.7, 8.7, and 6.5 degrees C, respectively. Binding of APV and Ac-PEP is most sensitive to the D25N mutation, as shown by DeltaT(m) ratios [DeltaT(m)(PR)/DeltaT(m)(PR(D25N))] of 35.8 and 16.3, respectively, whereas binding of DMP323 and RPB (DeltaT(m) ratios of 1-2) is least affected. Binding of the substrate-like inhibitors RPB and Ac-PEP is nearly abolished (DeltaT(m)(PR)/DeltaT(m)(PR(D29N)) > or = 44) by the D29N mutation, whereas this mutation only moderately affects binding of the smaller inhibitors (DeltaT(m) ratios of 1.4-2.2). Of the nine FDA-approved clinical HIV-1 protease inhibitors screened, APV, RTV, and DRV competitively inhibit porcine pepsin with K(i) values of 0.3, 0.6, and 2.14 microM, respectively. DSC results were consistent with this relatively weak binding of APV (DeltaT(m) 2.7 degrees C) compared with the tight binding of Ac-PEP (DeltaT(m) > or = 17 degrees C). Comparison of superimposed structures of the PR/APV complex with those of PR/Ac-PEP and pepsin/pepstatin A complexes suggests a role for Asp215, Asp32, and Ser219 in pepsin, equivalent to Asp25, Asp25', and Asp29 in PR in the binding and stabilization of the pepsin/APV complex.     

Exploration of the P6/P7 region of the peptide-binding site of the human class II major histocompatability complex protein HLA-DR1

Crystal structures of the class II major histocompatibilty complex (MHC) protein, HLA-DR1, generally show a tight fit between MHC and bound peptide except in the P6/P7 region of the peptide-binding site. In this region, there is a shallow water-filled pocket underneath the peptide and between the pockets that accommodate the P6 and P7 side chains. We investigated the properties of this pocket with the idea of engineering substitutions into the corresponding region of peptide antigens to increase their binding affinity for HLA-DR1. We investigated d-amino acids and N-alkyl modifications at both the P6 and P7 positions of the peptide and found that binding of peptides to HLA-DR1 could be increased by incorporating an N-methyl substitution at position 7 of the peptide. The crystal structure of HLA-DR1 bound to a peptide containing a P7 N-methyl alanine was determined. The N-methyl group orients in the P6/P7 pocket, displacing one of the waters usually bound in this pocket. The structure shows that the substitution does not alter the conformation of the bound peptide, which adopts the usual polyproline type II helix. An antigenic peptide carrying the N-methyl modification is taken up by antigen-presenting cells and loaded onto endogenous class II MHC molecules for presentation, and the resultant MHC-peptide complexes activate antigen-specific T-cells. These results suggest a possible strategy for increasing the affinity of weakly immunogenic peptides that might be applicable to the development of vaccines and diagnostic reagents.     

DRoP: Automated detection of conserved solvent-binding sites on proteins

Water and ligand binding play critical roles in the structure and function of proteins, yet their binding sites and significance are difficult to predict a priori. Multiple solvent crystal structures (MSCS) is a method where several X-ray crystal structures are solved, each in a unique solvent environment, with organic molecules that serve as probes of the protein surface for sites evolved to bind ligands, while the first hydration shell is essentially maintained. When superimposed, these structures contain a vast amount of information regarding hot spots of protein-protein or protein-ligand interactions, as well as conserved water-binding sites retained with the change in solvent properties. Optimized mining of this information requires reliable structural data and a consistent, objective analysis tool. Detection of related solvent positions (DRoP) was developed to automatically organize and rank the water or small organic molecule binding sites within a given set of structures. It is a flexible tool that can also be used in conserved water analysis given multiple structures of any protein independent of the MSCS method. The DRoP output is an HTML format list of the solvent sites ordered by conservation rank in its population within the set of structures, along with renumbered and recolored PDB files for visualization and facile analysis. Here, we present a previously unpublished set of MSCS structures of bovine pancreatic ribonuclease A (RNase A) and use it together with published structures to illustrate the capabilities of DRoP.     

Presence of ferric hydroxide clusters in mutants of Haemophilus influenzae ferric ion-binding protein A

The periplasmic iron binding protein plays an essential role in the iron uptake pathway of Gram-negative pathogenic bacteria from the Pasteurellaceae and Neisseriaceae families and is critical for survival of these pathogens within the host. In this study, we report the crystal structures of two mutant forms of ferric ion-binding protein A (FbpA) from Haemophilus influenzae with bound multinuclear oxo-metal clusters. Crystals of site-directed mutants in the metal or anion binding ligands contain protein in the open conformation, and two mutant FbpAs, H9A and N175L, contain different cluster arrangements in the iron-binding pocket. The iron clusters are anchored by binding to the two tyrosine ligands (Tyr195 and Tyr196) positioned at the vertex of the iron-binding pocket but are not coordinated by the other metal binding ligands. Our results suggest that the metal clusters may have formed in situ, suggesting that the mutant FbpAs may serve as a simple model for protein-mediated mineralization.     

Refined 2.3 A X-ray crystal structure of bovine thrombin complexes formed with the benzamidine and arginine-based thrombin inhibitors NAPAP, 4-TAPAP and MQPA. A starting point for improving antithrombotics.

Well-diffracting crystals of bovine epsilon-thrombin in complex with several "non-peptidic" benzamidine and arginine-based thrombin inhibitors have been obtained by co-crystallization. The 2.3 A crystal structures of three complexes formed either with NAPAP, 4-TAPAP, or MQPA, were solved by Patterson search methods and refined to crystallographic R-values of 0.167 to 0.178. The active-site environment of thrombin is only slightly affected by binding of the different inhibitors; in particular, the exposed "60-insertion loop" essentially maintains its typical projecting structure. The D-stereoisomer of NAPAP and the L-stereoisomer of MQPA bind to thrombin with very similar conformations, as previously inferred from their binding to bovine trypsin; the arginine side-chain of the latter inserts into the specificity pocket in a "non-canonical" manner. The L-stereoisomer of 4-TAPAP, whose binding geometry towards trypsin was only poorly defined, is bound to the thrombin active-site in a compact conformation. In contrast to NAPAP, the distal p-amidino/guanidino groups of 4-TAPAP and MQPA do not interact with the carboxylate group of Asp189 in the thrombin specificity pocket in a "symmetrical" twin N-twin O manner, but through "lateral" single N-twin O contacts; in contrast to the p-amidino group of 4-TAPAP, however, the guanidyl group of MQPA packs favourably in the pocket due to an elaborate hydrogen bond network, which includes two entrapped water molecules. These thrombin structures confirm previous conclusions of the important role of the intermolecular hydrogen bonds formed with Gly216, and of the good sterical fit of the terminal bulky hydrophobic inhibitor groups with the hydrophobic aryl binding site and the S2-cavity, respectively, for tight thrombin active site binding of these non-peptidic inhibitors. These accurate crystal structures are presumed to be excellent starting points for the design and the elaboration of improved antithrombotics.     

Potent cyclic urea HIV protease inhibitors with benzofused heterocycles as P2/P2′ groups

A series of benzofused heterocycles was examined to replace the metabolically unstable benzyl alcohol P2/P2′ groups of DMP 323 (1). Extremely potent inhibitors of HIV protease (Ki < 0.01 nM) and excellent antiviral activity (IC₉₀ = 8 nM) were found. An X-ray crystal structure of benzimidazolone 5a bound to HIV protease showed H-bonds to Asp30 and a bridging water molecule to Gly48.     

Structural analyses of covalent enzyme-substrate analog complexes reveal strengths and limitations of de novo enzyme design

We report the cocrystal structures of a computationally designed and experimentally optimized retro-aldol enzyme with covalently bound substrate analogs. The structure with a covalently bound mechanism-based inhibitor is similar to, but not identical with, the design model, with an RMSD of 1.4 Å over active-site residues and equivalent substrate atoms. As in the design model, the binding pocket orients the substrate through hydrophobic interactions with the naphthyl moiety such that the oxygen atoms analogous to the carbinolamine and β-hydroxyl oxygens are positioned near a network of bound waters. However, there are differences between the design model and the structure: the orientation of the naphthyl group and the conformation of the catalytic lysine are slightly different; the bound water network appears to be more extensive; and the bound substrate analog exhibits more conformational heterogeneity than typical native enzyme–inhibitor complexes. Alanine scanning of the active-site residues shows that both the catalytic lysine and the residues around the binding pocket for the substrate naphthyl group make critical contributions to catalysis. Mutating the set of water-coordinating residues also significantly reduces catalytic activity. The crystal structure of the enzyme with a smaller substrate analog that lacks naphthyl ring shows the catalytic lysine to be more flexible than in the naphthyl–substrate complex; increased preorganization of the active site would likely improve catalysis. The covalently bound complex structures and mutagenesis data highlight the strengths and weaknesses of the de novo enzyme design strategy.