HLA-DRB1*07:01 is associated with a higher risk of asparaginase allergies.

Asparaginase is a therapeutic enzyme used to treat leukemia and lymphoma, with immune responses resulting in suboptimal drug exposure and a greater risk of relapse. To elucidate whether there is a genetic component to the mechanism of asparaginase-induced immune responses, we imputed human leukocyte antigen (HLA) alleles in patients of European ancestry enrolled on leukemia trials at St. Jude Children's Research Hospital (n = 541) and the Children's Oncology Group (n = 1329). We identified a higher incidence of hypersensitivity and anti-asparaginase antibodies in patients with HLA-DRB1*07:01 alleles (P = 7.5 × 10(-5), odds ratio [OR] = 1.64; P = 1.4 × 10(-5), OR = 2.92, respectively). Structural analysis revealed that high-risk amino acids were located within the binding pocket of the HLA protein, possibly affecting the interaction between asparaginase epitopes and the HLA-DRB1 protein. Using a sequence-based consensus approach, we predicted the binding affinity of HLA-DRB1 alleles for asparaginase epitopes, and patients whose HLA genetics predicted high-affinity binding had more allergy (P = 3.3 × 10(-4), OR = 1.38). Our results suggest a mechanism of allergy whereby HLA-DRB1 alleles that confer high-affinity binding to asparaginase epitopes lead to a higher frequency of reactions. These trials were registered at www.clinicaltrials.gov as NCT00137111, NCT00549848, NCT00005603, and NCT00075725.

Broken-Symmetry DFT Computations for the Reaction Pathway of IspH, an Iron-Sulfur Enzyme in Pathogenic Bacteria.

The recently discovered methylerythritol phosphate (MEP) pathway provides new targets for the development of antibacterial and antimalarial drugs. In the final step of the MEP pathway, the [4Fe-4S] IspH protein catalyzes the 2e(-)/2H(+) reductive dehydroxylation of (E)-4-hydroxy-3-methyl-but-2-enyl diphosphate (HMBPP) to afford the isoprenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP). Recent experiments have attempted to elucidate the IspH catalytic mechanism to drive inhibitor development. Two competing mechanisms have recently emerged, differentiated by their proposed HMBPP binding modes upon 1e(-) reduction of the [4Fe-4S] cluster: (1) a Birch reduction mechanism, in which HMBPP remains bound to the [4Fe-4S] cluster through its terminal C4-OH group (ROH-bound) until the -OH is cleaved as water; and (2) an organometallic mechanism, in which the C4-OH group rotates away from the [4Fe-4S] cluster, allowing the HMBPP olefin group to form a metallacycle complex with the apical iron (η(2)-bound). We perform broken-symmetry density functional theory computations to assess the energies and reduction potentials associated with the ROH- and η(2)-bound states implicated by these competing mechanisms. Reduction potentials obtained for ROH-bound states are more negative (-1.4 to -1.0 V) than what is typically expected of [4Fe-4S] ferredoxin proteins. Instead, we find that η(2)-bound states are lower in energy than ROH-bound states when the [4Fe-4S] cluster is 1e(-) reduced. Furthermore, η(2)-bound states can already be generated in the oxidized state, yielding reduction potentials of ca. -700 mV when electron addition occurs after rotation of the HMBPP C4-OH group. We demonstrate that such η(2)-bound states are kinetically accessible both when the IspH [4Fe-4S] cluster is oxidized and 1e(-) reduced. The energetically preferred pathway gives 1e(-) reduction of the cluster after substrate conformational change, generating the 1e(-) reduced intermediate proposed in the organometallic mechanism.

Ligand binding mode prediction by docking: mdm2/mdmx inhibitors as a case study.

The p53-binding domains of Mdm2 and Mdmx, two negative regulators of the tumor suppressor p53, are validated targets for cancer therapeutics, but correct binding poses of some proven inhibitors, particularly the nutlins, have been difficult to obtain with standard docking procedures. Virtual screening pipelines typically draw from a database of compounds represented with 1D or 2D structural information from which one or more 3D conformations must be generated. These conformations are then passed to a docking algorithm that searches for optimal binding poses on the target protein. This work tests alternative pipelines using several commonly used conformation generation programs (LigPrep, ConfGen, MacroModel, and Corina/Rotate) and docking programs (GOLD, Glide, MOE-dock, and AutoDock Vina) for their ability to reproduce known poses for a series of Mdmx and/or Mdm2 inhibitors, including several nutlins. Most combinations of these programs using default settings fail to find correct poses for the nutlins but succeed for all other compounds. Docking success for the nutlin class requires either computationally intensive conformational exploration or an "anchoring" procedure that incorporates knowledge of the orientation of the central imidazoline ring.

Incomplete folding upon binding mediates Cdk4/cyclin D complex activation by tyrosine phosphorylation of inhibitor p27 protein.

p27(Kip1) (p27), an intrinsically disordered protein, regulates the various Cdk/cyclin complexes that control cell cycle progression. The kinase inhibitory domain of p27 contains a cyclin-binding subdomain (D1), a Cdk-binding subdomain (D2), and a linker helix subdomain that connects D1 and D2. Here, we report that, despite extensive sequence conservation between Cdk4/cyclin D1 (hereafter Cdk4/cyclin D) and Cdk2/cyclin A, the thermodynamic details describing how the individual p27 subdomains contribute to equally high affinity binding to these two Cdk/cyclin complexes are strikingly different. Differences in enthalpy/entropy compensation revealed that the D2 subdomain of p27 folds incompletely when binding Cdk4/cyclin D versus Cdk2/cyclin A. Incomplete binding-induced folding exposes tyrosine 88 of p27 for phosphorylation by the nonreceptor tyrosine kinase Abl. Importantly, tyrosine phosphorylation (of p27) relieves Cdk inhibition by p27, enabling cell cycle entry. Furthermore, the interaction between a conserved hydrophobic patch on cyclin D and subdomain D1 is much weaker than that with cyclin A; consequently, a construct containing subdomains D1 and LH (p27-D1LH) does not inhibit substrate binding to Cdk4/cyclin D as it does to Cdk2/cyclin A. Our results provide a mechanism by which Cdk4 (within the p27/Cdk4/cyclin D complex) is poised to be activated by extrinsic mitogenic signals that impinge upon p27 at the earliest stage of cell division. More broadly, our results further illustrate the regulatory versatility of intrinsically disordered proteins.

Identification and characterization of the first small molecule inhibitor of MDMX.

The p53 pathway is disrupted in virtually every human tumor. In approximately 50% of human cancers, the p53 gene is mutated, and in the remaining cancers, the pathway is dysregulated by genetic lesions in other genes that modulate the p53 pathway. One common mechanism for inactivation of the p53 pathway in tumors that express wild-type p53 is increased expression of MDM2 or MDMX. MDM2 and MDMX bind p53 and inhibit its function by distinct nonredundant mechanisms. Small molecule inhibitors and small peptides have been developed that bind MDM2 in the p53-binding pocket and displace the p53 protein, leading to p53-mediated cell cycle exit and apoptosis. To date, peptide inhibitors of MDMX have been developed, but no small molecule inhibitors have been reported. We have developed biochemical and cell-based assays for high throughput screening of chemical libraries to identify MDMX inhibitors and identified the first MDMX inhibitor SJ-172550. This compound binds reversibly to MDMX and effectively kills retinoblastoma cells in which the expression of MDMX is amplified. The effect of SJ-172550 is additive when combined with an MDM2 inhibitor. Results from a series of biochemical and structural modeling studies suggest that SJ-172550 binds the p53-binding pocket of MDMX, thereby displacing p53. This lead compound is a useful chemical scaffold for further optimization of MDMX inhibitors that may eventually be used to treat pediatric cancers and various adult tumors that overexpress MDMX or have similar genetic lesions. When combined with selective MDM2 inhibitors, SJ-172550 may also be useful for treating tumors that express wild-type p53.

Density functional theory analysis of structure, energetics, and spectroscopy for the Mn-Fe active site of Chlamydia trachomatis ribonucleotide reductase in four oxidation states.

Models for the Mn-Fe active site structure of ribonucleotide reductase (RNR) from pathogenic bacteria Chlamydia trachomatis (Ct) in different oxidation states have been studied in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor like screening (COSMO) solvation model and also with finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) calculations. The detailed structures for the reduced Mn(II)-Fe(II), the met Mn(III)-Fe(III), the oxidized Mn(IV)-Fe(III) and the superoxidized Mn(IV)-Fe(IV) states are predicted. The calculated properties, including geometries, (57)Fe Mossbauer isomer shifts and quadrupole splittings, and (57)Fe and (55)Mn electron nuclear double resonance (ENDOR) hyperfine coupling constants, are compared with the available experimental data. The Mössbauer and energetic calculations show that the (mu-oxo, mu-hydroxo) models better represent the structure of the Mn(IV)-Fe(III) state than the di-mu-oxo models. The predicted Mn(IV)-Fe(III) distances (2.95 and 2.98 A) in the (mu-oxo, mu-hydroxo) models are in agreement with the extended X-ray absorption fine structure (EXAFS) experimental value of 2.92 A (Younker et al. J. Am. Chem. Soc. 2008, 130, 15022-15027). The effect of the protein and solvent environment on the assignment of the Mn metal position is examined by comparing the relative energies of alternative mono-Mn(II) active site structures. It is proposed that if the Mn(II)-Fe(II) protein is prepared with prior addition of Mn(II) or with Mn(II) richer than Fe(II), Mn is likely positioned at metal site 2, which is further from Phe127.

Quantitative structure-activity relationship (QSAR) for a series of novel cannabinoid derivatives using descriptors derived from semi-empirical quantum-chemical calculations.

Recent work implicating the cannabinoid receptors in a wide range of human pathologies has intensified the need for reliable QSAR models for drug discovery and lead optimization. Predicting the ligand selectivity of the cannabinoid CB(1) and CB(2) receptors in the absence of generally accepted models for their structures requires a ligand-based approach, which makes such studies ideally suited for quantum-chemical treatments. We present a QSAR model for ligand-receptor interactions based on quantum-chemical descriptors (an eQSAR) obtained from PM3 semi-empirical calculations for a series of phenyl-substituted cannabinoids based on a ligand with known in vivo activity against glioma [Duntsch, C.; Divi, M. K.; Jones, T.; Zhou, Q.; Krishnamurthy, M.; Boehm, P.; Wood, G.; Sills, A.; Moore. B. M., II. J. Neuro-Oncol., 2006, 77, 143] and a set of structurally similar adamantyl-substituted cannabinoids. A good model for CB(2) inhibition (R(2)=0.78) has been developed requiring only four explanatory variables derived from semi-empirical results. The role of the ligand dipole moment is discussed and we propose that the CB(2) binding pocket likely possesses a significant electric field. Describing the affinities with respect to the CB(1) receptor was not possible with the current set of ligands and descriptors, although the attempt highlighted some important points regarding the development of QSAR models.

A fluid salt-bridging cluster and the stabilization of p53.

p53 is a homotetrameric tumor suppressor protein that is found to be mutated in most human cancers. Some of these mutations, particularly mutations to R337, fall in the tetramerization domain and cause defects in tetramer formation leading to loss of function. Mutation to His at this site has been found to destabilize the tetramer in a pH-dependent fashion. In structures of the tetramerization domain determined by crystallography, R337 from one monomer makes a salt bridge with D352 from another monomer, apparently helping to stabilize the tetramer. Here we present molecular dynamics simulations of wild-type p53 and the R337His mutant at several different pH and salt conditions. We find that the 337-352 salt bridge is joined by two other charged side chains, R333 and E349. These four residues do not settle into a fixed pattern of salt bridging, but continue to exchange salt-bridging partners on the nanosecond time scale throughout the simulation. This unusual system of fluid salt bridging may explain the previous finding from alanine scanning experiments that R333 contributes significantly to protein stability, even though in the crystal structure it is extended outward into solvent. This extended conformation of R333 appears to be the result of a specific crystal contact and, this contact being absent in the simulation, R333 turns inward to join its interaction partners. When R337 is mutated to His but remains positively charged, it maintains the original interaction with D352, but the newly observed interaction with E349 is weakened, accounting for the reduced stability of R337H even under mildly acidic conditions. When this His is deprotonated, the interaction with D352 is also lost, accounting for the further destabilization observed under mildly alkaline conditions. Simulations were carried out using both explicit and implicit solvent models, and both displayed similar behavior of the fluid salt-bridging cluster, suggesting that implicit solvent models can capture at least the qualitative features of this phenomenon as well as explicit solvent. Simulations under strongly acidic conditions in implicit solvent displayed the beginnings of the unfolding process, a destabilization of the hydrophobic dimer-dimer interface. Computational alanine scanning using the molecular mechanics Poisson-Boltzmann surface area method showed significant correlation to experimental unfolding data for charged and polar residues, but much weaker correlation for hydrophobic residues.

Performance of a docking/molecular dynamics protocol for virtual screening of nutlin-class inhibitors of Mdmx.

A virtual screening protocol involving docking and molecular dynamics has been tested against the results of fluorescence polarization assays testing the potency of a series of compounds of the nutlin class for inhibition of the interaction between p53 and Mdmx, an interaction identified as a driver of certain cancers. The protocol uses a standard docking method (AutoDock) with a cutoff based on the AutoDock score (ADscore), followed by molecular dynamics simulation with a cutoff based on root-mean-square-deviation (RMSD) from the docked pose. An analysis of the experimental and computational results shows modest performance of ADscore alone, but dramatically improved performance when RMSD is also used.

Disordered p27Kip1 exhibits intrinsic structure resembling the Cdk2/cyclin A-bound conformation.

p27Kip1 (p27) influences cell division by regulating nuclear cyclin-dependent kinases. Before binding, p27 is at least partially disordered and folds upon binding its Cdk/cyclin targets. 30-40% of human proteins, including p27, are predicted to contain disordered segments, and have been termed intrinsically unstructured proteins (IUPs). Unfortunately, the inherent dynamics of IUPs hamper detailed analysis of their structure/function relationships. Here, we describe the use of molecular dynamics (MD) computations and solution NMR spectroscopy to reveal that several segments of the p27 kinase inhibitory domain (p27-KID), in addition to the previously characterized helical segment, exist as highly populated, intrinsically folded structural units (IFSUs). Several IFSUs resemble structural features of bound p27-KID, while another exhibits alternative conformations. Interestingly, the highly conserved, specificity determining segment of p27 is shown to be highly disordered. Elucidation of IFSUs within p27-KID allows consideration of their influences on the thermodynamics and kinetics of Cdk/cyclin binding. The degree to which IFSUs are populated within p27-KID is surprising and suggests that other putative IUPs contain IFSUs that may be studied using similar techniques.

Monitoring Ligand-Induced Protein Ordering in Drug Discovery.

While the gene for p53 is mutated in many human cancers causing loss of function, many others maintain a wild-type gene but exhibit reduced p53 tumor suppressor activity through overexpression of the negative regulators, Mdm2 and/or MdmX. For the latter mechanism of loss of function, the activity of endogenous p53 can be restored through inhibition of Mdm2 or MdmX with small molecules. We previously reported a series of compounds based upon the Nutlin-3 chemical scaffold that bind to both MdmX and Mdm2 [Vara, B. A. et al. (2014) Organocatalytic, diastereo- and enantioselective synthesis of nonsymmetric cis-stilbene diamines: A platform for the preparation of single-enantiomer cis-imidazolines for protein-protein inhibition. J. Org. Chem. 79, 6913-6938]. Here we present the first solution structures based on data from NMR spectroscopy for MdmX in complex with four of these compounds and compare them with the MdmX:p53 complex. A p53-derived peptide binds with high affinity (Kd value of 150nM) and causes the formation of an extensive network of hydrogen bonds within MdmX; this constitutes the induction of order within MdmX through ligand binding. In contrast, the compounds bind more weakly (Kd values from 600nM to 12µM) and induce an incomplete hydrogen bond network within MdmX. Despite relatively weak binding, the four compounds activated p53 and induced p21(Cip1) expression in retinoblastoma cell lines that overexpress MdmX, suggesting that they specifically target MdmX and/or Mdm2. Our results document structure-activity relationships for lead-like small molecules targeting MdmX and suggest a strategy for their further optimization in the future by using NMR spectroscopy to monitor small-molecule-induced protein order as manifested through hydrogen bond formation.

Proton affinity changes driving unidirectional proton transport in the bacteriorhodopsin photocycle.

Bacteriorhodopsin is the smallest autonomous light-driven proton pump. Proposals as to how it achieves the directionality of its trans-membrane proton transport fall into two categories: accessibility-switch models in which proton transfer pathways in different parts of the molecule are opened and closed during the photocycle, and affinity-switch models, which focus on changes in proton affinity of groups along the transport chain during the photocycle. Using newly available structural data, and adapting current methods of protein protonation-state prediction to the non-equilibrium case, we have calculated the relative free energies of protonation microstates of groups on the transport chain during key conformational states of the photocycle. Proton flow is modeled using accessibility limitations that do not change during the photocycle. The results show that changes in affinity (microstate energy) calculable from the structural models are sufficient to drive unidirectional proton transport without invoking an accessibility switch. Modeling studies for the N state relative to late M suggest that small structural re-arrangements in the cytoplasmic side may be enough to produce the crucial affinity change of Asp96 during N that allows it to participate in the reprotonation of the Schiff base from the cytoplasmic side. Methodologically, the work represents a conceptual advance compared to the usual calculations of pK(a) using macroscopic electrostatic models. We operate with collective states of protonation involving all key groups, rather than the individual-group pK(a) values traditionally used. When combined with state-to-state transition rules based on accessibility considerations, a model for non-equilibrium proton flow is obtained. Such methods should also be applicable to other active proton-transport systems.

Structural details, pathways, and energetics of unfolding apomyoglobin.

Protein folding is often difficult to characterize experimentally because of the transience of intermediate states, and the complexity of the protein-solvent system. Atomistic simulations, which could provide more detailed information, have had to employ highly simplified models or high temperatures, to cope with the long time scales of unfolding; direct simulation of folding is even more problematic. We report a fully atomistic simulation of the acid-induced unfolding of apomyoglobin in which the protonation of acidic side-chains to simulate low pH is sufficient to induce unfolding at room temperature with no added biasing forces or other unusual conditions; and the trajectory is validated by comparison to experimental characterization of intermediate states. Novel insights provided by their analysis include: characterization of a dry swollen globule state forming a barrier to initial unfolding or final folding; observation of cooperativity in secondary and tertiary structure formation and its explanation in terms of dielectric environments; and structural details of the intermediate and the completely unfolded states. These insights involve time scales and levels of structural detail that are presently beyond the range of experiment, but come within reach through the simulation methods described here. An implicit solvation model is used to analyze the energetics of protein folding at various pH and ionic strength values, and a reasonable estimate of folding free energy is obtained. Electrostatic interactions are found to disfavor folding.

Exploring protein native states and large-scale conformational changes with a modified generalized born model.

Implicit solvation models provide, for many applications, a reasonably accurate and computationally effective way to describe the electrostatics of aqueous solvation. Here, a popular analytical Generalized Born (GB) solvation model is modified to improve its accuracy in calculating the solvent polarization part of free energy changes in large-scale conformational transitions, such as protein folding. In contrast to an earlier GB model (implemented in the AMBER-6 program), the improved version does not overstabilize the native structures relative to the finite-difference Poisson-Boltzmann continuum treatment. In addition to improving the energy balance between folded and unfolded conformers, the algorithm (available in the AMBER-7 and NAB molecular modeling packages) is shown to perform well in more than 50 ns of native-state molecular dynamics (MD) simulations of thioredoxin, protein-A, and ubiquitin, as well as in a simulation of Barnase/Barstar complex formation. For thioredoxin, various combinations of input parameters have been explored, such as the underlying gas-phase force fields and the atomic radii. The best performance is achieved with a previously proposed modification to the torsional potential in the Amber ff99 force field, which yields stable native trajectories for all of the tested proteins, with backbone root-mean-square deviations from the native structures being approximately 1.5 A after 6 ns of simulation time. The structure of Barnase/Barstar complex is regenerated, starting from an unbound state, to within 1.9 A relative to the crystal structure of the complex.

Macroscopic electrostatic models for protonation states in proteins.

The use of macroscopic electrostatic models to calculate the relative energetics of protonation states and the pH-titration properties of ionizable groups in proteins is described. These methods treat the protein as an irregularly-shaped low-dielectric object containing embedded atomic charges immersed in a high-dielectric (solvent) medium. The energetics of altering protonation states then involves the electrostatic work of altering the embedded atomic charges. The governing electrostatic equation is either the Poisson or linearized Poisson-Boltzmann equation, which generally requires numerical solution. A tutorial approach is taken, the main aim of which is a thorough understanding of the method.

Density Functional and Electrostatic Calculations of Manganese Superoxide Dismutase Active Site Complexes in Protein Environments.

Density functional and electrostatic methods have been applied to calculate active site geometries and the redox potential of manganese superoxide dismutase (MnSOD). The initial active site clusters were built up by including only first-shell side chain ligands and then augmented by second-shell ligands. The density functional optimized Mn-ligand bond lengths for the reduced complexes in general compared fairly well with protein crystallography data; however, large deviations for calculated Mn-OH distances were found for the oxidized active site clusters. Our calculations suggest that this deviation can be attributed to the redox heterogeneity of the oxidized protein in X-ray crystallography studies. The redox potential was calculated by treating the protein environment and the solvent bulk by a semimacroscopic electrostatic model. The protein structures were taken from the Thermus thermophilus enzyme. The calculated coupled redox potentials converge toward experimental values with increasing size of the active site cluster models, and the final calculated value was +0.06 V, compared to experimental values of +0.26 V determined for Bacillus stearothermophilus and +0.31 V in Escherichia coli enzymes. Using an energy decomposition scheme, the effects of the second-shell ligands and the protein and reaction fields have been analyzed.

CuZn Superoxide Dismutase Geometry Optimization, Energetics, and Redox Potential Calculations by Density Functional and Electrostatic Methods.

The structures, energetics, and orbital- and charge-dependent properties of copper zinc superoxide dismutase (CuZnSOD) have been studied using density functional and electrostatic methods. The CuZnSOD was represented with a model consisting of copper and zinc sites connected by a bridging histidine ligand. In addition to the bridge, three histidine ligands and one water molecule were bonded to the Cu ion in the copper site as first-shell ligands. Two histidine ligands and an aspartate were coordinated to the zinc ion in the zinc site. Full optimization of the model was performed using different functionals, both local and nonlocal. Geometrical parameters calculated with the nonlocal functionals agree well with the experimental X-ray data. In our calculated results, the His61 Nepsilon-Cu bond in the active site breaks during the reduction and protonation, consistent with a number of X-ray structures and with EXAFS and NMR evidence. The reduction potential and pK(a) of the coupled electron/proton reaction catalyzed by CuZnSOD were determined using different models for the extended environment-from an electrostatic representation of continuum solvent, to the full protein/solvent environment using a Poisson-Boltzmann method. The predicted redox potential and pK(a) values determined using the model with the full protein/solvent environment are in excellent agreement with experiment. Inclusion of the full protein environment is essential for an accurate description of the redox process. Although the zinc ion does not play a direct redox role in the dismutation, its electronic contribution is very important for the catalytic mechanism.

Use of Broken-Symmetry Density Functional Theory To Characterize the IspH Oxidized State: Implications for IspH Mechanism and Inhibition.

With current therapies becoming less efficacious due to increased drug resistance, new inhibitors of both bacterial and malarial targets are desperately needed. The recently discovered methylerythritol phosphate (MEP) pathway for isoprenoid synthesis provides novel targets for the development of such drugs. Particular attention has focused on the IspH protein, the final enzyme in the MEP pathway, which uses its [4Fe-4S] cluster to catalyze the formation of the isoprenoid precursors IPP and DMAPP from HMBPP. IspH catalysis is achieved via a 2e-/2H+ reductive dehydroxylation of HMBPP; the mechanism by which catalysis is achieved, however, is highly controversial. The work presented herein provides the first step in assessing different routes to catalysis by using computational methods. By performing broken-symmetry density functional theory (BS-DFT) calculations that employ both the conductor-like screening solvation model (DFT/COSMO) and a finite-difference Poisson-Boltzmann self-consistent reaction field methodology (DFT/SCRF), we evaluate geometries, energies, and Mössbauer signatures of the different protonation states that may exist in the oxidized state of the IspH catalytic cycle. From DFT/SCRF computations performed on the oxidized state, we find a state where the substrate, HMBPP, coordinates the apical iron in the [4Fe-4S] cluster as an alcohol group (ROH) to be one of two, isoenergetic, lowest-energy states. In this state, the HMBPP pyrophosphate moiety and an adjacent glutamate residue (E126) are both fully deprotonated, making the active site highly anionic. Our findings that this low-energy state also matches the experimental geometry of the active site and that its computed isomer shifts agree with experiment validate the use of the DFT/SCRF method to assess relative energies along the IspH reaction pathway. Additional studies of IspH catalytic intermediates are currently being pursued.

Identification and characterization of an allosteric inhibitory site on dihydropteroate synthase.

The declining effectiveness of current antibiotics due to the emergence of resistant bacterial strains dictates a pressing need for novel classes of antimicrobial therapies, preferably against molecular sites other than those in which resistance mutations have developed. Dihydropteroate synthase (DHPS) catalyzes a crucial step in the bacterial pathway of folic acid synthesis, a pathway that is absent in higher vertebrates. As the target of the sulfonamide class of drugs that were highly effective until resistance mutations arose, DHPS is known to be a valuable bacterial Achilles heel that is being further exploited for antibiotic development. Here, we report the discovery of the first known allosteric inhibitor of DHPS. NMR and crystallographic studies reveal that it engages a previously unknown binding site at the dimer interface. Kinetic data show that this inhibitor does not prevent substrate binding but rather exerts its effect at a later step in the catalytic cycle. Molecular dynamics simulations and quasi-harmonic analyses suggest that the effect of inhibitor binding is transmitted from the dimer interface to the active-site loops that are known to assume an obligatory ordered substructure during catalysis. Together with the kinetics results, these structural and dynamics data suggest an inhibitory mechanism in which binding at the dimer interface impacts loop movements that are required for product release. Our results potentially provide a novel target site for the development of new antibiotics.

Analysis of the active-site mechanism of tyrosyl-DNA phosphodiesterase I: a member of the phospholipase D superfamily.

Tyrosyl-DNA phosphodiesterase I (Tdp1) is a member of the phospholipase D superfamily that hydrolyzes 3'-phospho-DNA adducts via two conserved catalytic histidines-one acting as the lead nucleophile and the second acting as a general acid/base. Substitution of the second histidine specifically to arginine contributes to the neurodegenerative disease spinocerebellar ataxia with axonal neuropathy (SCAN1). We investigated the catalytic role of this histidine in the yeast protein (His432) using a combination of X-ray crystallography, biochemistry, yeast genetics, and theoretical chemistry. The structures of wild-type Tdp1 and His432Arg both show a phosphorylated form of the nucleophilic histidine that is not observed in the structure of His432Asn. The phosphohistidine is stabilized in the His432Arg structure by the guanidinium group that also restricts the access of nucleophilic water molecule to the Tdp1-DNA intermediate. Biochemical analyses confirm that His432Arg forms an observable and unique Tdp1-DNA adduct during catalysis. Substitution of His432 by Lys does not affect catalytic activity or yeast phenotype, but substitutions with Asn, Gln, Leu, Ala, Ser, and Thr all result in severely compromised enzymes and DNA topoisomerase I-camptothecin dependent lethality. Surprisingly, His432Asn did not show a stable covalent Tdp1-DNA intermediate that suggests another catalytic defect. Theoretical calculations revealed that the defect resides in the nucleophilic histidine and that the pK(a) of this histidine is crucially dependent on the second histidine and on the incoming phosphate of the substrate. This represents a unique example of substrate-activated catalysis that applies to the entire phospholipase D superfamily.