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.

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.

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.

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.

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.

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.

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.

Catalysis and sulfa drug resistance in dihydropteroate synthase.

The sulfonamide antibiotics inhibit dihydropteroate synthase (DHPS), a key enzyme in the folate pathway of bacteria and primitive eukaryotes. However, resistance mutations have severely compromised the usefulness of these drugs. We report structural, computational, and mutagenesis studies on the catalytic and resistance mechanisms of DHPS. By performing the enzyme-catalyzed reaction in crystalline DHPS, we have structurally characterized key intermediates along the reaction pathway. Results support an S(N)1 reaction mechanism via formation of a novel cationic pterin intermediate. We also show that two conserved loops generate a substructure during catalysis that creates a specific binding pocket for p-aminobenzoic acid, one of the two DHPS substrates. This substructure, together with the pterin-binding pocket, explains the roles of the conserved active-site residues and reveals how sulfonamide resistance arises.

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.

Mössbauer properties of the diferric cluster and the differential iron(II)-binding affinity of the iron sites in protein R2 of class Ia Escherichia coli ribonucleotide reductase: a DFT/electrostatics study.

The R2 subunit of class-Ia ribonucleotide reductase (RNR) from Escherichia coli (E. coli) contains a diiron active site. Starting from the apo-protein and Fe(II) in solution at low Fe(II)/apoR2 ratios, mononuclear Fe(II) binding is observed indicating possible different Fe(II) binding affinities for the two alternative sites. Further, based on their Mössbauer spectroscopy and two-iron-isotope reaction experiments, Bollinger et al. (J. Am. Chem. Soc., 1997, 119, 5976-5977) proposed that the site Fe1, which bonds to Asp84, should be associated with the higher observed (57)Fe Mössbauer quadrupole splitting (2.41 mm s(-1)) and lower isomer shift (0.45 mm s(-1)) in the Fe(III)Fe(III) state, site Fe2, which is further from Tyr122, should have a greater affinity for Fe(II) binding than site Fe1, and Fe(IV) in the intermediate X state should reside at site Fe2. In this paper, using density functional theory (DFT) incorporated with the conductor-like screening (COSMO) solvation model and with the finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) methodologies, we have demonstrated that the observed large quadrupole splitting for the diferric state R2 does come from site Fe1(III) and it is mainly caused by the binding position of the carboxylate group of the Asp84 sidechain. Further, a series of active site clusters with mononuclear Fe(II) binding at either site Fe1 or Fe2 have been studied, which show that with a single dielectric medium outside the active site quantum region, there is no energetic preference for Fe(II) binding at one site over another. However, when including the explicit extended protein environment in the PB-SCRF model, the reaction field favors the Fe(II) binding at site Fe2 rather than at site Fe1 by ~9 kcal mol(-1). Therefore our calculations support the proposal of the previous Mössbauer spectroscopy and two-iron-isotope reaction experiments by Bollinger et al.

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.

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.

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.

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.

Experimental and DFT studies: novel structural modifications greatly enhance the solvent sensitivity of live cell imaging dyes.

Structural modifications of previously reported merocyanine dyes (Toutchkine, A.; Kraynov, V.; Hahn, K. J. Am. Chem. Soc. 2003, 125, 4132-4145) were found to greatly enhance the solvent dependence of their absorbance and fluorescence emission maxima. Density functional theory (DFT) calculations have been performed to understand the differences in optical properties between the new and previously synthesized dyes. Absorption and emission energies were calculated for several new dyes using DFT vertical self-consistent reaction field (VSCRF) methods. Geometries of ground and excited states were optimized with a conductor-like screening model (COSMO) and self-consistent field (SCF) methods. The new dyes have enhanced zwitterionic character in the ground state and much lower polarity in the excited state, as shown by the DFT-VSCRF calculations. Consistently, the position of the absorption bands are strongly blue-shifted in more polar solvent (methanol compared to benzene), as predicted by the DFT spectral calculations. Inclusion of explicit H-bonding solvent molecules within the quantum model further enhances the predicted shifts and is consistent with the observed spectral broadening. Smaller but significant spectral shifts in polar versus nonpolar solvent are predicted and observed for emission bands. The new dyes show large fluorescence quantum yields in polar hydrogen-bonding solvents; qualitatively, the longest bonds along the conjugated chain at the excited S1 state minimum are shorter in the more polar solvent, inhibiting photoisomerization. The loss of photostability of the dyes is a consequence of the reaction with and electron transfer to singlet oxygen, starting oxidative dye cleavage. The calculated vertical ionization potentials of three dyes I-SO, AI-SO(4), and AI-BA(4) in benzene and methanol are consistent with their relative photobleaching rates; the charge distributions along the conjugated chains for the three dyes are similarly predictive of higher reaction rates for AI-SO(4) and AI-BA(4) than for I-SO. Time-dependent DFT calculations were also performed on AI-BA(4); these were less accurate than the VSCRF method in predicting the absorption energy shift from benzene to methanol.

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.

Model for proton transport coupled to protein conformational change: application to proton pumping in the bacteriorhodopsin photocycle.

A modeling method is presented for protein systems in which proton transport is coupled to conformational change, as in proton pumps and in motors driven by the proton-motive force. Previously developed methods for calculating pKa values in proteins using a macroscopic dielectric model are extended beyond the equilibrium case to a master-equation model for the time evolution of the system through states defined by ionization microstate and a discrete set of conformers. The macroscopic dielectric model supplies free energy changes for changes of protonation microstate, while the method for obtaining the energetics of conformational change and the relaxation rates, the other ingredients needed for the master equation, are system dependent. The method is applied to the photoactivated proton pump, bacteriorhodopsin, using conformational free energy differences from experiment and treating relaxation rates through three adjustable parameters. The model is found to pump protons with an efficiency relatively insensitive to parameter choice over a wide range of parameter values, and most of the main features of the known photocycle from very early M to the return to the resting state are reproduced. The boundaries of these parameter ranges are such that short-range proton transfers are faster than longer-range ones, which in turn are faster than conformational changes. No relaxation rates depend on conformation. The results suggest that an "accessibility switch", while not ruled out, is not required and that vectorial proton transport can be achieved through the coupling of the energetics of ionization and conformational states.

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.

On the role of the conserved aspartate in the hydrolysis of the phosphocysteine intermediate of the low molecular weight tyrosine phosphatase.

The usual rate-determining step in the catalytic mechanism of the low molecular weight tyrosine phosphatases involves the hydrolysis of a phosphocysteine intermediate. To explain this hydrolysis, general base-catalyzed attack of water by the anion of a conserved aspartic acid has sometimes been invoked. However, experimental measurements of solvent deuterium kinetic isotope effects for this enzyme do not reveal a rate-limiting proton transfer accompanying dephosphorylation. Moreover, base activation of water is difficult to reconcile with the known gas-phase proton affinities and solution phase pK(a)'s of aspartic acid and water. Alternatively, hydrolysis could proceed by a direct nucleophilic attack by a water molecule. To understand the hydrolysis mechanism, we have used high-level density functional methods of quantum chemistry combined with continuum electrostatics models of the protein and the solvent. Our calculations do not support a catalytic activation of water by the aspartate. Instead, they indicate that the water oxygen directly attacks the phosphorus, with the aspartate residue acting as a H-bond acceptor. In the transition state, the water protons are still bound to the oxygen. Beyond the transition state, the barrier to proton transfer to the base is greatly diminished; the aspartate can abstract a proton only after the transition state, a result consistent with experimental solvent isotope effects for this enzyme and with established precedents for phosphomonoester hydrolysis.

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.