Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

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.     

Isolation of a single activating allosteric interaction in phosphofructokinase from Escherichia coli

Escherichia coli phosphofructokinase 1 (EcPFK) is a homotetramer with four active and four allosteric sites. Understanding of the structural basis of allosteric activation of EcPFK by MgADP is complicated by the multiplicity of binding sites. To isolate a single heterotropic allosteric interaction, hybrid tetramers were formed between wild-type and mutant EcPFK subunits in which the binding sites of the mutant subunits have decreased affinity for their respective ligands. The 1:3 (wild-type:mutant) hybrid that contained only one native active site and one native allosteric site was isolated. The affinity for the substrate fructose-6-phosphate (Fru-6-P) of a single wild-type active site is greatly decreased over that displayed by the wild-type tetramer due to the lack of homotropic activation. The free energy of activation by MgADP for this heterotropic interaction is -0.58 kcal/mol at 8.5 degrees C. This compares to -2.87 kcal/mol for a hybrid with no homotropic coupling but all four unique heterotropic interactions. Therefore, the isolated interaction contributes 20% of the total heterotropic coupling. By comparison, wild-type EcPFK exhibits a coupling free energy between Fru-6-P and MgADP of -1.56 kcal/mol under these conditions, indicating that the effects of MgADP are diminished by a homotropic activation equal to -1.3 kcal/mol. These data are not consistent with a concerted allosteric mechanism.     

High free energy of lipid/protein interaction in biological membranes

The free energy of lipid/protein interaction in biological membranes is still unknown although extensive partitioning and modelling studies have revealed many partial energetic increments. Multiple site binding kinetics are now applied to four well-studied functional membrane proteins, and mean free energy values (+/-S.D.) of -4.23+/-0.49 kcal/mol for single lipid binding sites and of -89.7+/-35.4 kcal/mol for complete lipid substitution are obtained. These high free energy values point to an important bioenergetic role of lipid/protein interaction in membrane functions.     

Mean-field analysis of protein-protein interactions

Calculations were performed on the D1.3-E5.2 antibody-antibody complex estimating the binding affinities of the wild-type and 16 alanine substitutions. Analyzed were structural models of the interfacial region containing a zinc ion and crystallographic waters. A continuum approach was used to evaluate the electrostatic free energies and the hydrophobic effect was calculated by employing a buried molecular surface area relationship. Estimates of the absolute binding affinity reproduced the experimental value within the uncertainty of assessing entropic and strain energy contributions. The best correlation for mutants with experimental data was achieved when the hydrophilicity of created cavities were considered, and yielded a correlation coefficient of 0.7 and an average error of +/-1.4 kcal/mol. Empirically fitting the free energy function produced a smaller error of +/-1.0 kcal/mol. Depending on the electrical potential and electrostatic reorganization, scaling the 'protein dielectric constant' to approximately 10 may improve the accuracy of continuum models for evaluating amino acid substitutions.     

Elsamicin A binding to DNA. A comparative thermodynamic characterization

The antitumor drug elsamicin A contains a coumarin-related chartarin chromophore that intercalates into DNA. It differs from other related molecules in its disaccharide moiety, which bears an amino sugar. Its binding to DNA was analyzed using isothermal titration calorimetry and UV thermal denaturation, and characterized thermodynamically. For the association of elsamicin A with DNA we found DeltaG degrees = -8.6 kcal mol(-1), DeltaH = -10.4 kcal mol(-1), DeltaS = -6.1 cal mol(-1) K(-1), and Kobs = 2.8(+/- 0.2) x 10(6) M(-1) at 20 degrees C in 18 mM Na+. The contributions to the free energy of binding that lead to the DNA-elsamicin complex are compared with the binding to DNA of chartreusin, another chartarin-containing drug. The results are discussed in terms of the contributions of the disaccharide moieties into the strength of binding.     

A transition state in pieces: major contributions of entropic effects to ligand binding by adenosine deaminase.

Nebularine undergoes hydration at the active site of adenosine deaminase, in a reaction analogous to a partial reaction in the displacement of ammonia from adenosine by water, to generate an inhibitory complex that captures much of the binding affinity expected of an ideal transition-state analogue. Enzyme affinities of several compounds related to nebularine 1,6-hydrate, and to its stable analog 2'-deoxycoformycin, were compared in an effort to identify the structural origins of strong binding. Binding of the stable transition-state analog inhibitor 2'-deoxycoformycin was rendered 9.8 kcal/mol less favorable by removal of substituent ribose, 9.7 kcal/mol less favorable by inversion of the 8-hydroxyl substituent of the diazepine ring, and 10.0 kcal/mol less favorable by removal of atoms 4-6 of the diazepine ring. Binding of the unstable transition-state analog nebularine hydrate was rendered at least 9.9 kcal/mol less favorable by removal of the 6-hydroxyl group and 10.2 kcal/mol less favorable by removal of atoms 1-3 of the pyrimidine ring. In each case, the enzyme exhibited only modest affinity (Kd greater than or equal to 10(-2) M) for the "missing piece", indicating that incorporation of 2 binding determinants within a single molecule permits an additional 7-12 kcal/mol of intrinsic binding energy to be manifested as observed binding energy. These results are consistent with earlier indications that adenosine deaminase may use 10.5 kcal/mol of the intrinsic free energy of binding of the two substrates to place them in positions appropriate for reaction at the active site, overcoming the unfavorable entropy change of -35 eu for the equilibrium of 1,6-hydration of purine ribonucleoside and reducing the equilibrium constant for attainment of the transition state in deamination of adenosine. Thus, adenosine deaminase may achieve up to 8 orders of magnitude of its catalytic power by converting the nonenzymatic, bimolecular, hydration reaction to a monomolecular reaction at its active site. Several new 6-substituted 1,6-dihydropurine ribonucleosides, prepared by photoaddition of formate and by low-temperature addition of organolithium reagents to a derivative of purine ribonucleoside, exhibited Ki values of 9-1400 microM against adenosine deaminase, in accord with the active site's considerable tolerance of bulky leaving groups in substrates. Inhibition by one diastereomer of 6-carboxy-1,6-dihydropurine ribonucleoside was found to be time-dependent, progressing from a weakly bound to a more strongly bound complex.     

Reciprocal cooperative effects of multiple ligand binding to pyruvate kinase

The formation of multiple ligand complexes with muscle pyruvate kinase was measured in terms of dissociation constants and the standard free energies of formation were calculated. The binding of Mn2+ to the enzyme (KA = 55 +/- 5 X 10(-6) M; deltaF degrees = -5.75 +/- 0.05 kcal/mol) and to the enzyme saturated with phosphoenolpyruvate (conditional free energy) KA' = 0.8 +/- 0.4 X 10(-6) M; deltaF degrees = -8.22 +/- 0.34 kcal/mol) has been measured under identical conditions giving a free energy of coupling, delta(deltaF degrees) = -2.47 +/- 0.34 kcal/mol. Such a large negative free energy of coupling is diagnostic of a strong positively cooperative effect in ligand binding. The binding of the substrate phosphoenolpyruvate to free enzyme and the enzyme-Mn2+ complex was, by necessity, measured by different methods. The free energy of phosphoenolpyruvate binding to free enzyme (KS = 1.58 +/- 0.10 X 10(-4)M; deltaF degrees = -5.13 +/- 0.04 kcal/mol) and to the enzyme-Mn2+ complex (K3 = 0.75 +/- 0.10 X 10(-6)M; deltaF degrees = -8.26 +/- 0.07 kcal/mol) also gives a large negative free energy of coupling, delta(deltaF degrees) = -3.16 +/- 0.08 kcal/mol. Such a large negative value confirms reciprocal binding effects between the divalent cation and the substrate phosphoenolpyruvate. The binding of Mn2+ to the enzyme-ADP complex was also investigated and a free energy of coupling, delta(deltaF degrees) = -0.08 +/- 0.08 kcal/mol, was measured, indicative of little or no cooperativity in binding. The free energy of coupling with Mn2+ and pyruvate was measured as -1.52 +/- 0.14 kcal/mol, showing a significant amount of cooperativity in ligand binding but a substantially smaller effect than that observed for phosphoenolpyruvate binding. The magnitude of the coupling free energy may be related to the role of the divalent cation in the formation of the enzyme-substrate complexes. In the absence of the activating monovalent cation, the coupling free energies for phosphoenolpyruvate and pyruvate binding decrease by 40-60% and 25%, respectively, substantiating a role for the monovalent cation in the formation of enzyme-substrate complexes with phosphoenolpyruvate and with pyruvate.     

Crystal structure and ligand-binding studies of a screened peptide complexed with streptavidin.

The thermodynamic binding parameters and crystal structure for streptavidin-peptide complexes where the peptide sequences were obtained by random screening methods are reported. The affinities between streptavidin and two heptapeptides were determined by titrating calorimetric methods [Phe-Ser-His-Pro-Gln-Asn-Thr, Ka = 7944 (+/- 224) M-1, delta G degrees = -5.32 (+/- 0.01) kcal/mol, and delta H degrees = -19.34 (+/- 0.48) kcal/mol; His-Asp-His-Pro-Gln-Asn-Leu, Ka = 3542 (+/- 146) M-1, delta G degrees = -4.84 (+/- 0.03) kcal/mol, and delta H degrees = -19.00 (+/- 0.64) kcal/mol]. The crystal structure of streptavidin complexed with one of these peptides has been determined at 2.0-A resolution. The peptide (Phe-Ser-His-Pro-Gln-Asn-Thr) binds in a turn conformation with the histidine, proline, and glutamine side chains oriented inward at the biotin-binding site. A water molecule is immobilized between the histidine and glutamine side chains of the peptide and an aspartic acid side chain of the protein. Although some of the residues that participate in binding biotin also interact with the screened peptide, the peptide adopts an alternate method of utilizing binding determinants in the biotin-binding site of streptavidin.     

Linkage between cooperative oxygenation and subunit assembly of cobaltous human hemoglobin.

The thermodynamic linkage between cooperative oxygenation and dimer-tetramer subunit assembly has been determined for cobaltous human hemoglobin in which iron(II) protoporphyrin IX is replaced by cobalt(II) protoporphyrin IX. The equilibrium parameters of the linkage system were determined by global nonlinear least-squares regression of oxygenation isotherms measured over a range of hemoglobin concentrations together with the deoxygenated dimer-tetramer assembly free energy determined independently from forward and reverse reaction rates. The total cooperative free energy of tetrameric cobalt hemoglobin (over all four binding steps) is found to be 1.84 (+/- 0.13) kcal, compared with the native ferrous hemoglobin value of 6.30 (+/- 0.14) kcal. Detailed investigation of stepwise cooperativity effects shows the following: (1) The largest change occurs at the first ligation step and is determined on model-independent grounds by knowledge of the intermediate subunit assembly free energies. (2) Cooperativity in the shape of the tetrameric isotherm occurs mainly during the middle two steps and is concomitant with the release of quaternary constraints. (3) Although evaluation of the pure tetrameric isotherm portrays identical binding affinity between the last two steps, this apparent noncooperativity is the result of a "hidden" oxygen affinity enhancement at the last step of 0.48 (+/- 0.12) kcal. This quaternary enhancement energy is revealed by the difference in subunit assembly free energies of the triply and fully ligated species and is manifested visually by the oxygenation isotherms at high versus low hemoglobin concentration. (4) Cobaltous hemoglobin dimers exhibit apparent anticooperativity of 0.49 (+/- 0.16) kcal (presumed to arise from heterogeneity of subunit affinities).(ABSTRACT TRUNCATED AT 250 WORDS)     

Structural analysis of the carbohydrate backbone of Vibrio parahaemolyticus O2 lipopolysaccharides

A structural investigation has been carried out on the carbohydrate backbone of Vibrio parahaemolyticus O2 lipopolysaccharides (LPS) isolated by dephosphorylation, O-deacylation and N-deacylation. The carbohydrate backbone is a short-chain saccharide consisting of nine monosaccharide units i.e., 1 mol each of D-galactose (Gal), D-glucose (Glc), D-glucuronic acid (GlcA), L-glycero-D-manno-heptose (L,D-Hep), D-glycero-D-manno-heptose (D,D-Hep), 3-deoxy-D-manno-oct-2-ulosonic acid (Kdo), 5,7-diacetamido-3,5,7,9-tetradeoxy-D-glycero-D-galacto-non-2-ulosonic acid (NonlA), and 2 mol of 2-amino-2-deoxy-D-glucose (D-glucosamine, GlcN). Based on the data obtained by NMR spectroscopy, fast-atom bombardment mass spectrometry (FABMS) and methylation analysis, a structure was elucidated for the carbohydrate backbone of O2 LPS. In the native O2 LPS, the 2-amino-2-deoxy-D-glucitol (GlcN-ol) at the reducing end of the nonasaccharide is present as GlcN. The lipid A backbone is a beta-D-GlcN-(1-->6)-D-GlcN disaccharide as is the case for many Gram-negative bacterial LPS. The lipid A proximal Kdo is substituted by the distal part of the carbohydrate chain at position-5. In the native O2 LPS, D-galacturonic acid, which is liberated from LPS by mild acid treatment or by dephosphorylation in hydrofluoric acid, is present although its binding position is unknown at present.     

Mutation in the glmS gene responsible for the cell wall synthesis increases resistance to the herbicide amitrole in the cyanobacterium Synechocystis sp. PCC6803

A DNA fragment transforming the cells of the cyanobacterium Synechocystis sp. PCC 6803 to amitrole (3-amino-1,2,4-triazole) resistance was cloned from the Atr2 mutant resistant to this herbicide. The transforming activity of the cloned fragment was shown to be related to the missence-mutation "Val250-->Leu250" in the glmS gene encoding glucosamine-6-phosphate synthase, a key enzyme of cell wall synthesis. The amino acid substitution is localized in the central nonconservative part of the GlmS protein, far from two reaction centers positioned at the ends of a polypeptide. It is suggested that the mutant protein has lost sensitivity to amitrole. In the wild type strain, this herbicide causes conditional glucosamine auxotrophy (exogenous glucosamine restores ability of the cells to row in the presence of the lethal amitrole concentrations). Val250 is proposed to be allosteric binding site of AM in the GlmS protein of cyanobacterium.     

Alternative route for biosynthesis of amino sugars in Escherichia coli K-12 mutants by means of a catabolic isomerase.

By inserting a lambda placMu bacteriophage into gene glmS encoding glucosamine 6-phosphate synthetase (GlmS), the key enzyme of amino sugar biosynthesis, a nonreverting mutant of Escherichia coli K-12 that was strictly dependent on exogenous N-acetyl-D-glucosamine or D-glucosamine was generated. Analysis of suppressor mutations rendering the mutant independent of amino sugar supply revealed that the catabolic enzyme D-glucosamine-6-phosphate isomerase (deaminase), encoded by gene nagB of the nag operon, was able to fulfill anabolic functions in amino sugar biosynthesis. The suppressor mutants invariably expressed the isomerase constitutively as a result of mutations in nagR, the locus for the repressor of the nag regulon. Suppression was also possible by transformation of glmS mutants with high-copy-number plasmids expressing the gene nagB. Efficient suppression of the glmS lesion, however, required mutations in a second locus, termed glmX, which has been localized to 26.8 min on the standard E. coli K-12 map. Its possible function in nitrogen or cell wall metabolism is discussed.     

Analysis of the binding energies of testosterone, 5alpha-dihydrotestosterone,androstenedione and dehydroepiandrosterone sulfate with an antitestosterone antibody

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the binding of testosterone (TES), 5alpha-dihydrotestosterone (5ADHT), androstenedione (AND), and dehydroepiandrosterone sulfate (DHEAS) to the monoclonal antitestosterone antibody 3-C(4)F(5). The relative binding free energy of TES and AND was also calculated with free energy perturbation (FEP) simulations. The antibody 3-C(4)F(5) has a relatively high affinity (3 x 10(8) M(-1)) and on overall good binding profile for testosterone but its cross-reactivity with DHEAS has been the main reason for the failure to use this antibody in clinical immunoassays. The relative binding free energies obtained with the MM-PBSA method were 1.5 kcal/mol for 5ADHT, 3.8 kcal/mol for AND, and 4.3 kcal/mol for DHEAS, as compared to TES. When a water molecule of the ligand binding site, observed in the antibody-TES crystal structure, was explicitly included in MM-PBSA calculations, the relative binding energies were 3.4, 4.9, and 5.4 kcal/mol for 5ADHT, AND, and DHEAS, respectively. The calculated numbers are in correct order but larger than the corresponding experimental energies of 1.3, 1.5, and 2.6 kcal/mol, respectively. The fact that the MM-PBSA method reproduced the relative binding free energies of DHEAS, a steroid having a negatively charged sulfate group, and the neutrally charged TES, 5ADHT, and AND in satisfactory agreement with experiment shows the robustness of the method in predicting relative binding affinities. The 800-ps FEP simulations predicted that the antibody 3-C(4)F(5) binds TES 1.3 kcal/mol tighter than AND. Computational mutagenesis of selected amino acid residues of the ligand binding site revealed that the lower affinities of AND and DHEAS as compared to TES are due to a combined effect of several residues, each contributing a small fraction to the tighter binding of TES. An exception to this is Tyr99H, whose mutation to Ala lowered the binding of DHEAS 0.7 kcal/mol more than the binding of TES. This is probably due to the hydrogen bonding interaction formed between the OH group of Tyr99H and the sulfate group of DHEAS. Computational mutagensis data also showed that the affinity of the steroids to the antitestosterone antibody 3-C(4)F(5) would be enhanced if Trp47H were repositioned so that it would make more extensive contacts with the bound ligands. In addition, the binding of steroids to antitestosterone, antiprogesterone, and antiestradiol antibodies is discussed.     

Flanking DNA-sequences contribute to the specific binding of cI-repressor and OR1.

The binding of cI-repressor to a series of mutant operators containing OR1 of the right operator of bacteriophage lambda was investigated. Sites OR2 and/or OR3 were inactivated by either point or deletion mutations. The free energy of binding repressor to OR1 in the wildtype operator, delta G1, is -13.7 +/- 0.3 kcal/mol. delta G1 determined for an OR2- operator created by a single point mutation in OR2 is -13.6 +/- 0.2 kcal/mol. In contrast, delta G1 for the binding of repressor to a cloned synthetic OR1 operator containing only 24 bp of lambda sequence is -12.2 +/- 0.1 kcal/mol. When sequence 5' to OR1 is present, the binding affinity increases to -13.0 +/- 0.1 kcal/mol. In addition, the proximity of OR1 to a fragment-end decreases delta G1 from -13.7 to -12.3 +/- 0.1 kcal/mol. These results suggest that the DNA sequence outside the 17 bp OR1 binding-site contributes to the specific binding of cI-repressor.     

The crystal and solution studies of glucosamine-6-phosphate synthase from Candida albicans

Glucosamine 6-phosphate (GlcN-6-P) synthase is an ubiquitous enzyme that catalyses the first committed step in the reaction pathway that leads to formation of uridine 5'-diphospho-N-acetyl-D-glucosamine (UDP-GlcNAc), a precursor of macromolecules that contain amino sugars. Despite sequence similarities, the enzyme in eukaryotes is tetrameric, whereas in prokaryotes it is a dimer. The activity of eukaryotic GlcN-6-P synthase (known as Gfa1p) is regulated by feedback inhibition by UDP-GlcNAc, the end product of the reaction pathway, whereas in prokaryotes the GlcN-6-P synthase (known as GlmS) is not regulated at the post-translational level. In bacteria and fungi the enzyme is essential for cell wall synthesis. In human the enzyme is a mediator of insulin resistance. For these reasons, Gfa1p is a target in anti-fungal chemotherapy and in therapeutics for type-2 diabetes. The crystal structure of the Gfa1p isomerase domain from Candida albicans has been analysed in complex with the allosteric inhibitor UDP-GlcNAc and in the presence of glucose 6-phosphate, fructose 6-phosphate and an analogue of the reaction intermediate, 2-amino-2-deoxy-d-mannitol 6-phosphate (ADMP). A solution structure of the native Gfa1p has been deduced using small-angle X-ray scattering (SAXS). The tetrameric Gfa1p can be described as a dimer of dimers, with each half similar to the related enzyme from Escherichia coli. The core of the protein consists of the isomerase domains. UDP-GlcNAc binds, together with a metal cation, in a well-defined pocket on the surface of the isomerase domain. The residues responsible for tetramerisation and for binding UDP-GlcNAc are conserved only among eukaryotic sequences. Comparison with the previously studied GlmS from E. coli reveals differences as well as similarities in the isomerase active site. This study of Gfa1p focuses on the features that distinguish it from the prokaryotic homologue in terms of quaternary structure, control of the enzymatic activity and details of the isomerase active site.     

A novel empirical free energy function that explains and predicts protein-protein binding affinities

A free energy function can be defined as a mathematical expression that relates macroscopic free energy changes to microscopic or molecular properties. Free energy functions can be used to explain and predict the affinity of a ligand for a protein and to score and discriminate between native and non-native binding modes. However, there is a natural tension between developing a function fast enough to solve the scoring problem but rigorous enough to explain and predict binding affinities. Here, we present a novel, physics-based free energy function that is computationally inexpensive, yet explanatory and predictive. The function results from a derivation that assumes the cost of polar desolvation can be ignored and that includes a unique and implicit treatment of interfacial water-bridged interactions. The function was parameterized on an internally consistent, high quality training set giving R2=0.97 and Q2=0.91. We used the function to blindly and successfully predict binding affinities for a diverse test set of 31 wild-type protein-protein and protein-peptide complexes (R2=0.79, rmsd=1.2 kcal mol(-1)). The function performed very well in direct comparison with a recently described knowledge-based potential and the function appears to be transferable. Our results indicate that our function is well suited for solving a wide range of protein/peptide design and discovery problems.     

On the multiple functional roles of the active site histidine in catalysis and allosteric regulation of Escherichia coli glucosamine 6-phosphate deaminase

The active site of glucosamine-6-phosphate deaminase (EC 3.5.99.6, formerly 5.3.1.10) from Escherichia coli was first characterized on the basis of the crystallographic structure of the enzyme bound to the competitive inhibitor 2-amino-2-deoxy-glucitol 6-phosphate. The structure corresponds to the R allosteric state of the enzyme; it shows the side-chain of His143 in close proximity to the O5 atom of the inhibitor. This arrangement suggests that His143 could have a role in the catalysis of the ring-opening step of glucosamine 6-phosphate whose alpha-anomer is the true substrate. The imidazole group of this active-site histidine contacts the carboxy groups from Glu148 and Asp141, via its Ndelta1 atom [Oliva et al. (1995) Structure 3, 1323-1332]. These interactions change in the T state because the side chain of Glu148 moves toward the allosteric site, leaving at the active site the dyad Asp141-His143 [Horjales et al. (1999) Structure 7, 527-536]. In this research, a dual approach using site-directed mutagenesis and controlled chemical modification of histidine residues has been used to investigate the role of the active-site histidine. Our results support a multifunctional role of His143; in the forward reaction, it is involved in the catalysis of the ring-opening step of the substrate, glucosamine 6-P. In the reverse reaction, the substrate fructose 6-P binds in its open chain, carbonylic form. The role of His143 in the binding of both glucosamine 6-P and reaction intermediates in their extended-chain forms was demonstrated by binding experiments using the reaction intermediate analogue, 2-amino-2-deoxy-D-glucitol 6-phosphate. His143 was also shown to be a critical residue for the conformational coupling between active and allosteric sites. From the pH dependence of the reactivity of the active site histidine to diethyl dicarbonate, we observed a pK(a) change of 1.2 units to the acid side when the enzyme undergoes the allosteric T to R transition during which the side chain of Glu148 moves toward the active site. The kinetic study of the Glu148-Gln mutant deaminase shows that the loss of the carboxy group and its replacement with the corresponding amide modifies the k(cat) versus pH profile of the enzyme, suggesting that the catalytic step requiring the participation of His143 has become rate-limiting. This, in turn, indicates that the interaction Glu148-His143 in the wild-type enzyme in the R state contributes to make the enzyme functional over a wide pH range.     

Ligand interactions of the cation-dependent mannose 6-phosphate receptor. Comparison with the cation-independent mannose 6-phosphate receptor

The interactions of the bovine cation-dependent mannose 6-phosphate receptor with monovalent and divalent ligands have been studied by equilibrium dialysis. This receptor appears to be a homodimer or a tetramer. Each mole of receptor monomer bound 1.2 mol of the monovalent ligands, mannose 6-phosphate and pentamannose phosphate with Kd values of 8 X 10(-6) M and 6 X 10(-6) M, respectively and 0.5 mol of the divalent ligand, a high mannose oligosaccharide with two phosphomonoesters, with a Kd of 2 X 10(-7) M. When Mn2+ was replaced by EDTA in the dialysis buffer, the Kd for pentamannose phosphate was 2.5 X 10(-5) M. By measuring the affinity of the cation-dependent and cation-independent mannose 6-phosphate receptors for a variety of mannose 6-phosphate analogs, we conclude that the 6-phosphate and the 2-hydroxyl of mannose 6-phosphate each contribute approximately 4-5 kcal/mol of Gibb's free energy to the binding reaction. Neither receptor appears to interact substantially with the anomeric oxygen of mannose 6-phosphate. The receptors differ in that the cation-dependent receptor displays no detectable affinity for N-acetylglucosamine 1'-(alpha-D-methylmannopyranose 6-monophosphate) whereas this ligand binds to the cation-independent receptor with a poor, but readily measurable Kd of about 0.1 mM. The spacing of the mannose 6-phosphate-binding sites relative to each other may also differ for the two receptors.