Hydration energy landscape of the active site cavity in cytochrome P450cam

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein–ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol⁻¹) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable. Proteins 32:381–396, 1998. © 1998 Wiley-Liss, Inc.     

WaterScore: a novel method for distinguishing between bound and displaceable water molecules in the crystal structure of the binding site of protein-ligand complexes

We have performed a multivariate logistic regression analysis to establish a statistical correlation between the structural properties of water molecules in the binding site of a free protein crystal structure, with the probability of observing the water molecules in the same location in the crystal structure of the ligand-complexed form. The temperature B-factor, the solvent-contact surface area, the total hydrogen bond energy and the number of protein-water contacts were found to discriminate between bound and displaceable water molecules in the best regression functions obtained. These functions may be used to identify those bound water molecules that should be included in structure-based drug design and ligand docking algorithms.FIGURE The binding site ( thin sticks) of penicillopepsin (3app) with its crystallographically determined water molecules ( spheres) and superimposed ligand (in thick sticks, from complexed structure 1ppk). Water molecules sterically displaced by the ligand upon complexation are shown in cyan. Bound water molecules are shown in blue. Displaced water molecules are shown in yellow. Water molecules removed from the analysis due to a lack of hydrogen bonds to the protein are shown in white. WaterScore correctly predicted waters in blue as Probability=1 to remain bound and waters in yellow as Probability<1x10(-20) to remain bound.     

Crystal structure of a catalytic site mutant of beta-amylase from Bacillus cereus var. mycoides cocrystallized with maltopentaose

The X-ray crystal structure of a catalytic site mutant of beta-amylase, E172A (Glu172 --> Ala), from Bacillus cereus var. mycoides complexed with a substrate, maltopentaose (G5), and the wild-type enzyme complexed with maltose were determined at 2.1 and 2.0 A resolution, respectively. Clear and continuous density corresponding to G5 was observed in the active site of E172A, and thus, the substrate, G5, was not hydrolyzed. All glucose residues adopted a relaxed (4)C(1) conformation, and the conformation of the maltose unit for Glc2 and Glc3 was much different from those of other maltose units, where each glucose residue of G5 is named Glc1-Glc5 (Glc1 is at the nonreducing end). A water molecule was observed 3.3 A from the C1 atom of Glc2, and 3.0 A apart from the OE1 atom of Glu367 which acts as a general base. In the wild-type enzyme-maltose complex, two maltose molecules bind at subsites -2 and -1 and at subsites +1 and +2 in tandem. The conformation of the maltose molecules was similar to that of the condensation product of soybean beta-amylase, but differed from that of G5 in E172A. When the substrate flips between Glc2 and Glc3, the conformational energy of the maltose unit was calculated to be 20 kcal/mol higher than that of the cis conformation by MM3. We suggest that beta-amylase destabilizes the bond that is to be broken in the ES complex, decreasing the activation energy, DeltaG(++), which is the difference in free energy between this state and the transition state.     

Water molecules in the nucleotide binding cleft of actin: effects on subunit conformation and implications for ATP hydrolysis

In the monomeric actin crystal structure, the positions of a highly organized network of waters are clearly visible within the active site. However, the recently proposed models of filamentous actin (F-actin) did not extend to including these waters. Since the water network is important for ATP hydrolysis, information about water position is critical to understanding the increased rate of catalysis upon filament formation. Here, we show that waters in the active site are essential for intersubdomain rotational flexibility and that they organize the active-site structure. Including the crystal structure waters during simulation setup allows us to observe distinct changes in the active-site structure upon the flattening of the actin subunit, as proposed in the Oda model for F-actin. We identify changes in both protein position and water position relative to the phosphate tail that suggest a mechanism for accelerating the rate of nucleotide hydrolysis in F-actin by stabilizing charge on the β-phosphate and by facilitating deprotonation of catalytic water.     

The free energy landscape of small molecule unbinding

The spontaneous dissociation of six small ligands from the active site of FKBP (the FK506 binding protein) is investigated by explicit water molecular dynamics simulations and network analysis. The ligands have between four (dimethylsulphoxide) and eleven (5-diethylamino-2-pentanone) non-hydrogen atoms, and an affinity for FKBP ranging from 20 to 0.2 mM. The conformations of the FKBP/ligand complex saved along multiple trajectories (50 runs at 310 K for each ligand) are grouped according to a set of intermolecular distances into nodes of a network, and the direct transitions between them are the links. The network analysis reveals that the bound state consists of several subbasins, i.e., binding modes characterized by distinct intermolecular hydrogen bonds and hydrophobic contacts. The dissociation kinetics show a simple (i.e., single-exponential) time dependence because the unbinding barrier is much higher than the barriers between subbasins in the bound state. The unbinding transition state is made up of heterogeneous positions and orientations of the ligand in the FKBP active site, which correspond to multiple pathways of dissociation. For the six small ligands of FKBP, the weaker the binding affinity the closer to the bound state (along the intermolecular distance) are the transition state structures, which is a new manifestation of Hammond behavior. Experimental approaches to the study of fragment binding to proteins have limitations in temporal and spatial resolution. Our network analysis of the unbinding simulations of small inhibitors from an enzyme paints a clear picture of the free energy landscape (both thermodynamics and kinetics) of ligand unbinding.     

Insights from atomic-resolution X-ray structures of chemically synthesized HIV-1 protease in complex with inhibitors

The human immunodeficiency virus 1 (HIV-1) protease (PR) is an aspartyl protease essential for HIV-1 viral infectivity. HIV-1 PR has one catalytic site formed by the homodimeric enzyme. We chemically synthesized fully active HIV-1 PR using modern ligation methods. When complexed with the classic substrate-derived inhibitors JG-365 and MVT-101, the synthetic HIV-1 PR formed crystals that diffracted to 1.04- and 1.2-A resolution, respectively. These atomic-resolution structures revealed additional structural details of the HIV-1 PR's interactions with its active site ligands. Heptapeptide inhibitor JG-365, which has a hydroxyethylamine moiety in place of the scissile bond, binds in two equivalent antiparallel orientations within the catalytic groove, whereas the reduced isostere hexapeptide MVT-101 binds in a single orientation. When JG-365 was converted into the natural peptide substrate for molecular dynamic simulations, we found putative catalytically competent reactant states for both lytic water and direct nucleophilic attack mechanisms. Moreover, free energy perturbation calculations indicated that the insertion of catalytic water into the catalytic site is an energetically favorable process.     

Conservation of water molecules in an antibody–antigen interaction

The solvation of the antibody–antigen Fv D1.3–lysozyme complex is investigated through a study of the conservation of water molecules in crystal structures of the wild-type Fv fragment of antibody D1.3, 5 free lysozyme, the wild-type Fv D1.3–lysozyme complex, 5 Fv D1.3 mutants complexed with lysozyme and the crystal structure of an idiotope (Fv D1.3)-abti-idiotope (Fv E5.2) complex. In all, there are 99 water molecules common to the wild-type and mutant antibody–lysozyme complexes. The antibody–lysozyme interface includes 25 well-ordered solvent molecules, conserved among the wild-type and mutant Fv D1.3–lysozyme complexes, which are bound directly or through other water molecules to both antibody and antigen. In addition to contributing hydrogen bonds to the antibody–antigen interaction the solvent molecules fill many interface cavities. Comparison with x-ray crystal structures of free Fv D1.3 and free lysozyme shows that 20 of these conserved interface waters in the complex were bound to one of the free proteins. Uo to 23 additional water molecules are also found in the antibody–antigen interface, however these waters do no bridge antibody and antigen and their temperature factors are much higher than those of the 25 well-ordered waters. Fifteen water molecules are displaced to form the complex, some of which are substituted by hydrophilic protein atoms, and 5 water molecules are added at the antibody–antigen interface with the formation of the complex. While the current crystal models of the D1.3–lysozyme complex do not demonstrate the increase in bound waters found in a physico-chemical study of the interaction at decreased water activities, the 25 well-ordered interface water contribute a net gain of 10 hydrogen bonds to complex stability.     

Crystal structure of substrate-free Pseudomonas putida cytochrome P-450

The crystal structure of Pseudomonas putida cytochrome P-450cam in the substrate-free form has been refined at 2.20-A resolution and compared to the substrate-bound form of the enzyme. In the absence of the substrate camphor, the P-450cam heme iron atom is hexacoordinate with the sulfur atom of Cys-357 providing one axial heme ligand and a water molecule or hydroxide ion providing the other axial ligand. A network of hydrogen-bonded solvent molecules occupies the substrate pocket in addition to the iron-linked aqua ligand. When a camphor molecule binds, the active site waters including the aqua ligand are displaced, resulting in a pentacoordinate high-spin heme iron atom. Analysis of the Fno camphor - F camphor difference Fourier and a quantitative comparison of the two refined structures reveal that no detectable conformational change results from camphor binding other than a small repositioning of a phenylalanine side chain that contacts the camphor molecule. However, large decreases in the mean temperature factors of three separate segments of the protein centered on Tyr-96, Thr-185, and Asp-251 result from camphor binding. This indicates that camphor binding decreases the flexibility in these three regions of the P-450cam molecule without altering the mean position of the atoms involved.     

Incorporating replacement free energy of binding‐site waters in molecular docking

Binding‐site water molecules play a crucial role in protein‐ligand recognition, either being displaced upon ligand binding or forming water bridges to stabilize the complex. However, rigorously treating explicit binding‐site waters is challenging in molecular docking, which requires to fully sample ensembles of waters and to consider the free energy cost of replacing waters. Here, we describe a method to incorporate structural and energetic properties of binding‐site waters into molecular docking. We first developed a solvent property analysis (SPA) program to compute the replacement free energies of binding‐site water molecules by post‐processing molecular dynamics trajectories obtained from ligand‐free protein structure simulation in explicit water. Next, we implemented a distance‐dependent scoring term into DOCK scoring function to take account of the water replacement free energy cost upon ligand binding. We assessed this approach in protein targets containing important binding‐site waters, and we demonstrated that our approach is reliable in reproducing the crystal binding geometries of protein‐ligand‐water complexes, as well as moderately improving the ligand docking enrichment performance. In addition, SPA program (free available to academic users upon request) may be applied in identifying hot‐spot binding‐site residues and structure‐based lead optimization. Proteins 2014; 82:1765–1776. © 2014 Wiley Periodicals, Inc.     

Active-site gorge and buried water molecules in crystal structures of acetylcholinesterase from Torpedo californica

Buried water molecules and the water molecules in the active-site gorge are analyzed for five crystal structures of acetylcholinesterase from Torpedo californica in the resolution range 2.2-2.5 A (native enzyme, and four inhibitor complexes). A total of 45 buried hydration sites are identified, which are populated with between 36 and 41 water molecules. About half of the buried water is located in a distinct region neighboring the active-site gorge. Most of the buried water molecules are very well conserved among the five structures, and have low displacement parameters, B, of magnitudes similar to those of the main-chain atoms of the central beta-sheet structure. The active-site gorge of the native enzyme is filled with over 20 water molecules, which have poor hydrogen-bond coordination with an average of 2.9 polar contacts per water molecule. Upon ligand binding, distinct groups of these water molecules are displaced, whereas the others remain in positions similar to those that they occupy in the native enzyme. Possible roles of the buried water molecules are discussed, including their possible action as a lubricant to allow large-amplitude fluctuations of the loop structures forming the gorge wall. Such fluctuations are required to facilitate traffic of substrate, products and water molecules to and from the active-site. Because of their poor coordination, the gorge water molecules can be considered as "activated" as compared to bulk water. This should allow their easy displacement by incoming substrate. The relatively loose packing of the gorge water molecules leaves numerous small voids, and more efficient space-filling by substrates and inhibitors may be a major driving force of ligand binding.     

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.     

Mutational and structural studies of the diisopropylfluorophosphatase from Loligo vulgaris shed new light on the catalytic mechanism of the enzyme

The active site, the substrate binding site, and the metal binding sites of the diisopropylfluorophosphatase (DFPase) from Loligo vulgaris have been modified by means of site-directed mutagenesis to improve our understanding of the reaction mechanism. Enzymatic characterization of mutants located in the major groove of the substrate binding pocket indicates that large hydrophobic side chains at these positions are favorable for substrate turnover. Moreover, the active site residue His287 proved to be beneficial, but not essential, for DFP hydrolysis. In most cases, hydrophobic side chains at position 287 led to significant catalytic activities although reduced relative to the wild-type enzyme. With respect to the Ca-1 binding site, where catalysis occurs, various mutants indicated that the net charge at this calcium-binding site as well as the relative positions of the charged calcium ligands is crucial for catalytic activity. The importance of the electrostatic potential at the active site was furthermore revealed by various mutations of residues lining the interior of the central water-filled tunnel, which traverses the entire protein structure. In this respect, the structural features of residue His181, which is located at the opposite end of the DFPase tunnel relative to the active site, were characterized extensively. It was concluded that a tunnel-spanning hydrogen bond network, which includes a large number of apparently slow exchanging water molecules, relays any modifications in the electrostatics of the system to the active site, thus affecting the catalytic reactivity of the enzyme.     

Enhancement of catalytic efficiency of enzymes through exposure to anhydrous organic solvent at 70 degrees C. Three-dimensional structure of a treated serine proteinase at 2.2 A resolution

The enzyme behavior in anhydrous media has important applications in biotechnology. So far chemical modifications and protein engineering have been used to alter the catalytic power of the enzymes. For the first time, it is demonstrated that an exposure of enzyme to anhydrous organic solvents at optimized high temperature enhances its catalytic power through local changes at the binding region. Six enzymes: proteinase K, wheat germ acid phosphatase, alpha-amylase, beta-glucosidase, chymotrypsin and trypsin have been exposed to acetonitrile at 70 degrees C for three hours. The activities of these enzymes were found to be considerably enhanced. In order to understand the basis of this change in the activity of these enzymes, the structure of one of these treated enzymes, proteinase K has been analyzed in detail using X-ray diffraction method. The overall structure of the enzyme is similar to the native structure in aqueous environment. The hydrogen bonding system of the catalytic triad is intact after the treatment. However, the water structure in the substrate binding site undergoes some rearrangement as some of the water molecules are either displaced or completely absent. The most striking observation concerning the water structure pertains to the complete deletion of the water molecule which occupied the position at the so-called oxyanion hole in the active site of the native enzyme. Three acetonitrile molecules were found in the present structure. All the acetonitrile molecules are located in the recognition site. The sites occupied by acetonitrile molecules are independent of water molecules. The acetonitrile molecules are involved in extensive interactions with the protein atoms. All of them are interlinked through water molecules. The methyl group of one of the acetonitrile molecules (CCN1) interacts simultaneously with the hydrophobic side chains of Leu-96, Ile-107, and Leu-133. The development of such a hydrophobic environment at the recognition site introduces a striking conformation change in Ile-107 by rotating its side chain about C(alpha)--C(beta) bond by 180 degrees to bring about the delta-methyl group within the range of attractive van der Waals interactions with the methyl group of CCN1. A similar change has earlier been observed in proteinase K when it is complexed to a substrate analog lactoferrin fragment.     

Enhancement of catalytic activity of enzymes by heating in anhydrous organic solvents: 3D structure of a modified serine proteinase at high resolution

For the first time, it is demonstrated that exposure of an enzyme to anhydrous organic solvents at optimized high temperature enhances its catalytic power through local changes at the binding region. Six enzymes, namely, proteinase K, wheat germ acid phosphatase, alpha-amylase, beta-glucosidase, chymotrypsin and trypsin were exposed to acetonitrile at 70 degrees C for three hr. The activities of these enzymes were found to be considerably enhanced. In order to understand the basis of this change in the activity of these enzymes, proteinase K was analyzed in detail using X-ray diffraction method. The overall structure of the enzyme was found to be similar to the native structure in aqueous environment. The hydrogen bonding system of the catalytic triad remained intact after the treatment. However, the water structure in the substrate binding site underwent some rearrangement as some of the water molecules were either displaced or completely absent. The most striking observation concerning the water structure was the complete deletion of the water molecule which occupied the position at the so-called oxyanion hole in the active site of the native enzyme. Three acetonitrile molecules were found in the present structure. All the acetonitrile molecules were located in the recognition site. Interlinked through water molecules, the sites occupied by acetonitrile molecules were independent of water molecules. The acetonitrile molecules are involved in extensive interactions with the protein atoms. The methyl group of one of the acetonitrile molecules (CCN1) interacts simultaneously with the hydrophobic side chains of Leu 96, Ile 107 and Leu 133. The development of such a hydrophobic environment at the recognition site introduced a striking conformation change in Ile 107 by rotating its side chain about C alpha-C beta bond by 180 degrees to bring about the delta-methyl group within the range of attractive van der Waals interactions with the methyl group of CCN1. A similar change had earlier been observed in proteinase K when it was complexed to a substrate analogue, lactoferrin fragment.     

Explaining an unusually fast parasitic enzyme: folate tail-binding residues dictate substrate positioning and catalysis in Cryptosporidium hominis thymidylate synthase

The essential enzyme TS-DHFR from Cryptosporidium hominis undergoes an unusually rapid rate of catalysis at the conserved TS domain, facilitated by two nonconserved residues, Ala287 and Ser290, in the folate tail-binding region. Mutation of these two residues to their conserved counterparts drastically affects multiple steps of the TS catalytic cycle. We have determined the crystal structures of all three mutants (A287F, S290G, and A287F/S290G) in complex with active site ligands dUMP and CB3717. The structural data show two effects of the mutations: an increased distance between the ligands in the active site and increased flexibility of the folate ligand in the partially open enzyme state that precedes conformational change to the active catalytic state. The latter effect is able to be rescued by the mutants containing the A287F mutation. In addition, the conserved water network of TS is altered in each of the mutants. The structural results point to a role of the folate tail-binding residues in closely positioning ChTS ligands and restricting ligand flexibility in the partially open state to allow for a rapid transition to the active closed state and enhanced rate of catalysis. These results provide an explanation on how folate tail-binding residues at one end of the active site affect long-range interactions throughout the TS active site and validate these residues as targets for species-specific drug design.     

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.     

Active site geometry of glucose-1-phosphate uridylyltransferase

Glucose-1-phosphate uridylyltransferase, or UGPase, catalyzes the production of UDP-glucose from glucose-1-phosphate and UTP. Because of the biological role of UDP-glucose in glycogen synthesis and in the formation of glycolipids, glycoproteins, and proteoglycans, the enzyme is widespread in nature. Recently this laboratory reported the three-dimensional structure of UGPase from Escherichia coli. While the initial X-ray analysis revealed the overall fold of the enzyme, details concerning its active site geometry were limited because crystals of the protein complexed with either substrates or products could never be obtained. In an effort to more fully investigate the active site geometry of the enzyme, UGPase from Corynebacterium glutamicum was subsequently cloned and purified. Here we report the X-ray structure of UGPase crystallized in the presence of both magnesium and UDP-glucose. Residues involved in anchoring the ligand to the active site include the polypeptide chain backbone atoms of Ala 20, Gly 21, Gly 117, Gly 180, and Ala 214, and the side chains of Glu 36, Gln 112, Asp 143, Glu 201, and Lys 202. Two magnesium ions are observed coordinated to the UDP-glucose. An alpha- and a beta-phosphoryl oxygen, three waters, and the side chain of Asp 142 ligate the first magnesium, whereas the second ion is coordinated by an alpha-phosphoryl oxygen and five waters. The position of the first magnesium is conserved in both the glucose-1-phosphate thymidylyltransferases and the cytidylyltransferases. The structure presented here provides further support for the role of the conserved magnesium ion in the catalytic mechanisms of the sugar-1-phosphate nucleotidylyltransferases.     

The structure of a glycogen phosphorylase glucopyranose spirohydantoin complex at 1.8 A resolution and 100 K: the role of the water structure and its contribution to binding.

A glucopyranose spirohydantoin (a pyranose analogue of the potent herbicide, hydantocidin) has been identified as the highest affinity glucose analogue inhibitor of glycogen phosphorylase b (GPb). In order to elucidate the structural features that contribute to the binding, the structures of GPb in the native T state conformation and in complex with glucopyranose spirohydantoin have been determined at 100 K to 2.0 A and 1.8 A resolution, respectively, and refined to crystallographic R values of 0.197 (R[free] 0.248) and 0.182 (R[free] 0.229), respectively. The low temperature structure of GPb is almost identical to that of the previously determined room temperature structure, apart from a decrease in overall atomic temperature factors ((B) room temperature GPb = 34.9 A2; (B) 100 K GPb = 23.4 A2). The glucopyranose spirohydantoin inhibitor (Ki = 3.0 microM) binds at the catalytic site and induces small changes in two key regions of the protein: the 280s loop (residues 281-286) that results in a decrease in mobility of this region, and the 380s loop (residues 377-385) that undergoes more significant shifts in order to optimize contact to the ligand. The hydantoin group, that is responsible for increasing the affinity of the glucose compound by a factor of 10(3), makes only one hydrogen bond to the protein, from one of its NH groups to the main chain oxygen of His377. The other polar groups of the hydantoin group form hydrogen bonds to five water molecules. These waters are involved in extensive networks of hydrogen bonds and appear to be an integral part of the protein structure. Analysis of the water structure at the catalytic site of the native enzyme, shows that five waters are displaced by ligand binding and that there is a significant decrease in mobility of the remaining waters on formation of the GPb-hydantoin complex. The ability of the inhibitor to exploit existing waters, to displace waters and to recruit new waters appears to be important for the high affinity of the inhibitor.     

Crystal Structures of Wild-type and Mutant Methicillin-resistant Staphylococcus aureus Dihydrofolate Reductase Reveal an Alternate Conformation of NADPH That May Be Linked to Trimethoprim Resistance

Both hospital- and community-acquired Staphylococcus aureus infections have become major health concerns in terms of morbidity, suffering and cost. Trimethoprim-sulfamethoxazole (TMP-SMZ) is an alternative treatment for methicillin-resistant S. aureus (MRSA) infections. However, TMP-resistant strains have arisen with point mutations in dihydrofolate reductase (DHFR), the target for TMP. A single point mutation, F98Y, has been shown biochemically to confer the majority of this resistance to TMP. Using a structure-based approach, we have designed a series of novel propargyl-linked DHFR inhibitors that are active against several trimethoprim-resistant enzymes. We screened this series against wild-type and mutant (F98Y) S. aureus DHFR and found that several are active against both enzymes and specifically that the meta-biphenyl class of these inhibitors is the most potent. In order to understand the structural basis of this potency, we determined eight high-resolution crystal structures: four each of the wild-type and mutant DHFR enzymes bound to various propargyl-linked DHFR inhibitors. In addition to explaining the structure-activity relationships, several of the structures reveal a novel conformation for the cofactor, NADPH. In this new conformation that is predominantly associated with the mutant enzyme, the nicotinamide ring is displaced from its conserved location and three water molecules complete a network of hydrogen bonds between the nicotinamide ring and the protein. In this new position, NADPH has reduced interactions with the inhibitor. An equilibrium between the two conformations of NADPH, implied by their occupancies in the eight crystal structures, is influenced both by the ligand and the F98Y mutation. The mutation induced equilibrium between two NADPH-binding conformations may contribute to decrease TMP binding and thus may be responsible for TMP resistance.     

Structure of native and apo carbonic anhydrase II and structure of some of its anion-ligand complexes.

In order to obtain a better structural framework for understanding the catalytic mechanism of carbonic anhydrase, a number of inhibitor complexes of the enzyme were investigated crystallographically. The three-dimensional structure of free human carbonic anhydrase II was refined at pH 7.8 (1.54 A resolution) and at pH 6.0 (1.67 A resolution). The structure around the zinc ion was identical at both pH values. The structure of the zinc-free enzyme was virtually identical with that of the native enzyme, apart from a water molecule that had moved 0.9 A to fill the space that would be occupied by the zinc ion. The complexes with the anionic inhibitors bisulfite and formate were also studied at neutral pH. Bisulfite binds with one of its oxygen atoms, presumably protonized, to the zinc ion and replaces the zinc water. Formate, lacking a hydroxyl group, is bound with its oxygen atoms not far away from the position of the non-protonized oxygen atoms of the bisulfite complex, i.e. at hydrogen bond distance from Thr199 N and at a position between the zinc ion and the hydrophobic part of the active site. The result of these and other studies have implications for our view of the catalytic function of the enzyme, since virtually all inhibitors share some features with substrate, product or expected transition states. A reaction scheme where electrophilic activation of carbon dioxide plays an important role in the hydration reaction is presented. In the reverse direction, the protonized oxygen of the bicarbonate is forced upon the zinc ion, thereby facilitating cleavage of the carbon-oxygen bond. This is achieved by the combined action of the anionic binding site, which binds carboxyl groups, the side-chain of threonine 199, which discriminates between hydrogen bond donors and acceptors, and hydrophobic interaction between substrate and the active site cavity. The required proton transfer between the zinc water and His64 can take place through water molecules 292 and 318.