Effect of the phosphate substrate on drug-inhibitor binding to human purine nucleoside phosphorylase

The thermodynamics of the drug-inhibitors acyclovir, ganciclovir, and 9-benzylguanine binding to human purine nucleoside phosphorylase (hsPNP) were determined from isothermal titration calorimetry as a function of the substrate phosphate ion (Pi) concentration from 0 to 0.125 M and temperature from 15 °C to 35 °C. At 25 °C and with an increase in the Pi concentration from 0 to 50 mM, acyclovir binding becomes more entropically-driven and ganciclovir binding becomes more enthalpically-driven. At 25 °C, the tighter 9-benzylguanine binding reaction goes from an enthalpically-driven reaction in the absence of Pi to an entropically-driven reaction at 10 mM Pi, and the enthalpically-driven nature of the binding reaction is restored at 75 mM Pi. Since the dependencies of the driving-nature of the binding reactions on Pi concentration can be simulated by Pi binding to its catalytic site, it is believed that bound Pi affects the interactions of the side-chains with the ribose catalytic site. However, the binding constants are unaffected by change in the bound Pi concentration because of enthalpy-entropy compensation. The enzymatic activity of hsPNP was determined by an ITC-based assay employing 7-methylguanosine and Pi as the substrates. The heat of reaction determined from the assay increased by 7.5 kJ mol⁻¹ with increase in Pi concentration from 50 to 100 mM and is attributed to weak binding of the Pi to a secondary regulatory site. Although the binding constants of acyclovir and ganciclovir at 20 μM hsPNP were in agreement with the inverse inhibition constants determined from the ITC enzyme inhibition assays at 60 nM, the binding constant of 9-benzylguanine, which interacts with Phe159 from an adjacent subunit, decreased from 5.62 × 10⁵ M⁻¹ to 1.14 × 10⁵ M⁻¹. This reduction in the 9-benzylguanine binding affinity along with a 7-fold increase in the specific activity of hsPNP at 14.5 nM results from partial dissociation of the hsPNP trimer into monomers below the 60 nM level.     

Structure of a mutant human purine nucleoside phosphorylase with the prodrug, 2‐fluoro‐2′‐deoxyadenosine and the cytotoxic drug, 2‐fluoroadenine

A double mutant of human purine nucleoside phosphorylase (hDM) with the amino acid mutations Glu201Gln:Asn243Asp cleaves adenosine‐based prodrugs to their corresponding cytotoxic drugs. When fused to an anti‐tumor targeting component, hDM is targeted to tumor cells, where it effectively catalyzes phosphorolysis of the prodrug, 2‐fluoro‐2′‐deoxyadenosine (F‐dAdo) to the cytotoxic drug, 2‐fluoroadenine (F‐Ade). This cytotoxicity should be restricted only to the tumor microenvironment, because the endogenously expressed wild type enzyme cannot use adenosine‐based prodrugs as substrates. To gain insight into the interaction of hDM with F‐dAdo, we have determined the crystal structures of hDM with F‐dAdo and F‐Ade. The structures reveal that despite the two mutations, the overall fold of hDM is nearly identical to the wild type enzyme. Importantly, the residues Gln201 and Asp243 introduced by the mutation form hydrogen bond contacts with F‐dAdo that result in its binding and catalysis. Comparison of substrate and product complexes suggest that the side chains of Gln201 and Asp243 as well as the purine base rotate during catalysis possibly facilitating cleavage of the glycosidic bond. The two structures suggest why hDM, unlike the wild‐type enzyme, can utilize F‐dAdo as substrate. More importantly, they provide a critical foundation for further optimization of cleavage of adenosine‐based prodrugs, such as F‐dAdo by mutants of human purine nucleoside phosphorylase.     

The structure of the large regulatory α subunit of phosphorylase kinase examined by modeling and hydrogen‐deuterium exchange

Phosphorylase kinase (PhK), a 1.3 MDa regulatory enzyme complex in the glycogenolysis cascade, has four copies each of four subunits, (αβγδ)₄, and 325 kDa of unique sequence (the mass of an αβγδ protomer). The α, β and δ subunits are regulatory, and contain allosteric activation sites that stimulate the activity of the catalytic γ subunit in response to diverse signaling molecules. Due to its size and complexity, no high resolution structures have been solved for the intact complex or its regulatory α and β subunits. Of PhK's four subunits, the least is known about the structure and function of its largest subunit, α. Here, we have modeled the full‐length α subunit, compared that structure against previously predicted domains within this subunit, and performed hydrogen‐deuterium exchange on the intact subunit within the PhK complex. Our modeling results show α to comprise two major domains: an N‐terminal glycoside hydrolase domain and a large C‐terminal importin α/β‐like domain. This structure is similar to our previously published model for the homologous β subunit, although clear structural differences are present. The overall highly helical structure with several intervening hinge regions is consistent with our hydrogen‐deuterium exchange results obtained for this subunit as part of the (αβγδ)₄ PhK complex. Several low exchanging regions predicted to lack ordered secondary structure are consistent with inter‐subunit contact sites for α in the quaternary structure of PhK; of particular interest is a low‐exchanging region in the C‐terminus of α that is known to bind the regulatory domain of the catalytic γ subunit.     

Direct evidence by H/D exchange and ESI-MS for transient unproductive domain interaction in the refolding of an antibody scFv fragment

The refolding kinetics of a single-chain Fv (scFv) fragment, derived from a stabilized mutant of the phosphorylcholine binding antibody McPC603, was investigated by H/D exchange and ESI-MS and compared with the folding kinetics of its constituting domains V(H) and V(L). Both V(H) and V(L) adopt essentially native-like exchange protection within the dead time of the manual-mixing H/D exchange experiment (10 s) and in the case of V(L), which contains two cis-prolines in the native conformation, this fast protection is independent of proline cis/trans isomerization. At the earliest time point resolvable by manual mixing, fewer deuterons are protected in the scFv fragment than in the two isolated domains together, despite the fact that the scFv fragment is significantly more stable than V(L) and V(H). Full H/D exchange protection in the scFv fragment is gained on a time scale of minutes. This means that the domains in the scFv fragment do not refold independently. Rather, they associate prematurely and in nonnative form, a kinetic trap. Unproductive domain association is observed both after equilibrium- and short-term denaturation. For the equilibrium-denatured scFv fragment, whose native structure formation is dependent on a cis conformation of an interface proline in V(L), this cis/trans isomerization reaction proceeds about one order in magnitude more slowly than the escape from the trap to a conformation where full H/D exchange protection is already achieved. We interpret these data in terms of a general kinetic scheme involving intermediates with and without domain association.     

Cardiolipin interacts with beta‐2‐glycoprotein I and forms an open conformation–Mechanisms analyzed using hydrogen/deuterium exchange

Beta‐2‐glycoprotein I (β₂GPI) is the major antigen for the antiphospholipid antibodies in the antiphospholipid syndrome. The exposed epitope in domain I of β₂GPI can be recognized by the anti‐β₂GPI antibody. Here, we prepared the anionic di‐oleoyl‐phosphatidylserine (DOPS) and cardiolipin (CL) liposomes to interact with the β₂GPI. The conformational changes of β₂GPI upon binding with the liposomes were analyzed using hydrogen/deuterium exchange mass spectrometry. The exchange level of sequences 21–27 significantly increased after β₂GPI had interacted with DOPS. This change indicated a reduced interaction between domain I and domain V, inferring to a protrusion of the sequences 21–27 from the ring conformation. After β₂GPI had interacted with CL for 30 min, the exchange levels in 4 of the 5 domains increased significantly. The deuteration levels of sequences 1–20, 21–27, 196–205, 273–279 and 297–306 increased, suggesting that these regions had become more exposed, and the domain I was no longer in contact with domain V. The increasing deuteration levels in sequences 70–86, 153–162, 191–198, 196–205 and 273–279 indicated β₂GPI undergoing conformational changes to expose these inner regions, suggesting a structural transition. Overall, DOPS and CL induced minor conformational changes of β₂GPI at sequences 21–27 and forms an intermediate conformation after 10 min of interaction. After a complete protein–lipid interaction, high negatively charged CL membrane induced a major conformation transition of β₂GPI.     

Detection of early unfolding events in a dimeric protein by amide proton exchange and native electrospray mass spectrometry

Oligomeric proteins generally undergo unfolding through a dissociation/denaturation mechanism wherein the subunits first dissociate and then unfold. This mechanism can be detected by the fact that the proteins exhibit a concentration dependence of the denaturation curve. However, the concentration dependence does not answer the question of whether there are thermally induced conformational changes that facilitate subunit dissociation. To fully probe these mechanisms it is desirable to have an analytical approach that is capable of measuring both subunit dissociation and protein denaturation in a highly sensitive manner. In this article, we demonstrate that the combined use of native mass spectrometry to detect subunit mixing, and amide hydrogen/deuterium exchange to detect transient unfolding events can provide a very unique insight into the pre‐melting transitions in a protein oligomer. Both methods keep an isotopic record of each transformation event, without the dependence on equilibrium of the unfolding reaction. Here, we use a combined form of H/D exchange/mass spectrometry and isotopic labeling/native electrospray mass spectrometry to study the pre‐unfolding events of Bacillus subtilis NAD⁺ synthetase, a symmetrical dimer protein, which plays a vital role in the lifecycle of the bacteria. In the experimental outcome provided, we were able to clearly illustrate that at elevated temperatures, the NAD synthetase dimer undergoes reversible dissociation without monomer unfolding, while at temperatures where monomer unfolding is observed to take place, the rate of dimer dissociation still yet exceeds the rate of unfolding. Information provided by combining these two mass spectrometric methods was found to be very robust, and allowed us to establish an NAD synthetase unfolding model, where primary dissociation occurs prior to the complete unfolding of the NAD⁺ synthetase.     

Deuterium exchange and mass spectrometry reveal the interaction differences of two synthetic modulators of RXRalpha LBD

Protein amide hydrogen/deuterium (H/D) exchange was used to compare the interactions of two antagonists, UVI 2112 and UVI 3003, with that of the agonist, 9-cis-retinoic acid, upon binding to the human retinoid X receptor alpha ligand-binding domain (hRXRalpha LBD) homodimer. Analysis of the H/D content by mass spectrometry showed that in comparison to 9-cis-retinoic acid, the antagonists provide much greater protection toward deuterium exchange-in throughout the protein, suggesting that the protein-antagonist complex adopts a more restricted conformation or ensemble of conformations in which solvent accesses to amide protons are reduced. A comparison between the two antagonists shows that UVI 3003 is more protective in the C-terminal region due to the extra hydrophobic interactions derived from the atoms in the benzene ring of the carboxylic acid chain. It was less protective within regions comprising peptides 271-278 and 326-330 due to differences in conformational orientation, and/or shorter carboxylic acid chain length. Decreased deuterium exchange-in in the segment 234-239 where the residues do not involve interactions with the ligand was observed with the two antagonists, but not with 9-cis-RA. The amide protons of helix 12 of the agonist- or antagonist-occupied protein in solution have the same deuterium exchange rates as the unliganded protein, supporting a suggestion made previously that helix 12 can cover the occupied binding cavity only with the cofactor present to adjust its location.     

Two different proteins that compete for binding to thrombin have opposite kinetic and thermodynamic profiles

Thrombin binds thrombomodulin (TM) at anion binding exosite 1, an allosteric site far from the thrombin active site. A monoclonal antibody (mAb) has been isolated that competes with TM for binding to thrombin. Complete binding kinetic and thermodynamic profiles for these two protein-protein interactions have been generated. Binding kinetics were measured by Biacore. Although both interactions have similar K(D)s, TM binding is rapid and reversible while binding of the mAb is slow and nearly irreversible. The enthalpic contribution to the DeltaG(bind) was measured by isothermal titration calorimetry and van't Hoff analysis. The contribution to the DeltaG(bind) from electrostatic steering was assessed from the dependence of the k(a) on ionic strength. Release of solvent H(2)O molecules from the interface was assessed by monitoring the decrease in amide solvent accessibility at the interface upon protein-protein binding. The mAb binding is enthalpy driven and has a slow k(d). TM binding appears to be entropy driven and has a fast k(a). The favorable entropy of the thrombin-TM interaction seems to be derived from electrostatic steering and a contribution from solvent release. The two interactions have remarkably different thermodynamic driving forces for competing reactions. The possibility that optimization of binding kinetics for a particular function may be reflected in different thermodynamic driving forces is discussed.     

PAS domain allostery and light-induced conformational changes in photoactive yellow protein upon I2 intermediate formation, probed with enhanced hydrogen/deuterium exchange mass spectrometry

Photoactive yellow protein (PYP) is a small bacterial photoreceptor that undergoes a light-activated reaction cycle. PYP is also the prototypical Per-Arnt-Sim (PAS) domain. PAS domains, found in diverse multi-domain proteins from bacteria to humans, mediate protein-protein interactions and function as sensors and signal transducers. Here, we investigate conformational and dynamic changes in solution in wild-type PYP upon formation of the long-lived putative signaling intermediate I2 with enhanced hydrogen/deuterium exchange mass spectrometry (DXMS). The DXMS results showed that the central beta-sheet remains stable but specific external protein segments become strongly deprotected. Light-induced disruption of the dark-state hydrogen bonding network in I2 produces increased flexibility and opening of PAS core helices alpha3 and alpha4, releases the beta4-beta5 hairpin, and propagates conformational changes to the central beta-sheet. Surprisingly, the first approximately 10 N-terminal residues, which are essential for fast dark-state recovery from I2, become more protected. By combining the DXMS results with our crystallographic structures, which reveal detailed changes near the chromophore but limited protein conformational change, we propose a mechanism for I2 state formation. This mechanism integrates the results from diverse biophysical studies of PYP, and links an allosteric T to R-state conformational transition to three pathways for signal propagation within the PYP fold. On the basis of the observed changes in PYP plus commonalities shared among PAS domain proteins, we further propose that PAS domains share this conformational mechanism, which explains the versatile signal transduction properties of the structurally conserved PYP/PAS module by framework-encoded allostery.     

cGMP-binding prepares PKG for substrate binding by disclosing the C-terminal domain

Type I cyclic guanosine 3′,5′-monophosphate (cGMP)-dependent protein kinase (PKG) is involved in the nitric oxide/cGMP signaling pathway. PKG has been identified in many different species, ranging from unicelõlular organisms to mammals. The enzyme serves as one of the major receptor proteins for intracellular cGMP and controls a variety of cellular responses, ranging from smooth-muscle relaxation to neuronal synaptic plasticity. In the absence of a crystal structure, the three-dimensional structure of the homodimeric 152-kDa kinase PKG is unknown; however, there is evidence that the kinase adopts a distinct cGMP-dependent active conformation when compared to the inactive conformation. We performed mass-spectrometry-based hydrogen/deuterium exchange experiments to obtain detailed information on the structural changes in PKG Iα induced by cGMP activation. Site-specific exchange measurements confirmed that the autoinhibitory domain and the hinge region become more solvent exposed, whereas the cGMP-binding domains become more protected in holo-PKG (dimeric PKG saturated with four cGMP molecules bound). More surprisingly, our data revealed a specific disclosure of the substrate-binding region of holo-PKG, shedding new light into the kinase-activation process of PKG.     

Conserved and nonconserved features of the folding pathway of hisactophilin, a beta-trefoil protein

Based on previous studies of interleukin-1beta (IL-1beta) and both acidic and basic fibroblast growth factors (FGFs), it has been suggested that the folding of beta-trefoil proteins is intrinsically slow and may occur via the formation of essential intermediates. Using optical and NMR-detected quenched-flow hydrogen/deuterium exchange methods, we have measured the folding kinetics of hisactophilin, another beta-trefoil protein that has < 10% sequence identity and unrelated function to IL-1beta and FGFs. We find that hisactophilin can fold rapidly and with apparently two-state kinetics, except under the most stabilizing conditions investigated where there is evidence for formation of a folding intermediate. The hisactophilin intermediate has significant structural similarities to the IL-1beta intermediate that has been observed experimentally and predicted theoretically using a simple, topology-based folding model; however, it appears to be different from the folding intermediate observed experimentally for acidic FGF. For hisactophilin and acidic FGF, intermediates are much less prominent during folding than for IL-1beta. Considering the structures of the different beta-trefoil proteins, it appears that differences in nonconserved loops and hydrophobic interactions may play an important role in differential stabilization of the intermediates for these proteins.     

H/D exchange mass spectrometry and statistical coupling analysis reveal a role for allostery in a ferredoxin-dependent bifurcating transhydrogenase catalytic cycle

Recent investigations into ferredoxin-dependent transhydrogenases, a class of enzymes responsible for electron transport, have highlighted the biological importance of flavin-based electron bifurcation (FBEB). FBEB generates biomolecules with very low reduction potential by coupling the oxidation of an electron donor with intermediate potential to the reduction of high and low potential molecules. Bifurcating systems can generate biomolecules with very low reduction potentials, such as reduced ferredoxin (Fd), from species such as NADPH. Metabolic systems that use bifurcation are more efficient and confer a competitive advantage for the organisms that harbor them. Structural models are now available for two NADH-dependent ferredoxin-NADP⁺ oxidoreductase (Nfn) complexes. These models, together with spectroscopic studies, have provided considerable insight into the catalytic process of FBEB. However, much about the mechanism and regulation of these multi-subunit proteins remains unclear. Using hydrogen/deuterium exchange mass spectrometry (HDX-MS) and statistical coupling analysis (SCA), we identified specific pathways of communication within the model FBEB system, Nfn from Pyrococus furiosus, under conditions at each step of the catalytic cycle. HDX-MS revealed evidence for allosteric coupling across protein subunits upon nucleotide and ferredoxin binding. SCA uncovered a network of co-evolving residues that can provide connectivity across the complex. Together, the HDX-MS and SCA data show that protein allostery occurs across the ensemble of iron?sulfur cofactors and ligand binding sites using specific pathways that connect domains allowing them to function as dynamically coordinated units.     

Molecular dissection of the inhibitor binding pocket of mitotic kinesin Eg5 reveals mutants that confer resistance to antimitotic agents

The mitotic kinesin Eg5 plays an essential role in establishing the bipolar spindle. Recently, several antimitotic inhibitors have been shown to share a common binding region on Eg5. Considering the importance of Eg5 as a potential drug target for cancer chemotherapy it is essential to understand the molecular mechanism, by which these agents block Eg5 activity, and to determine the "key residues" crucial for inhibition. Eleven residues in the inhibitor binding pocket were mutated and the effects were monitored by kinetic analysis and mass spectrometry. Mutants R119A, D130A, P131A, I136A, V210A, Y211A and L214A abolish the inhibitory effect of monastrol. Results for W127A and R221A are less striking, but inhibitor constants are still considerably modified compared to wild-type Eg5. Only one residue, Leu214, was found to be essential for inhibition by STLC. W127A, D130A, V210A lead to increased K(i)(app) values, but binding of STLC is still tight. R119A, P131A, Y211A and R221A convert STLC into a classical rather than a tight-binding inhibitor with increased inhibitor constants. These results demonstrate that monastrol and STLC interact with different amino acids within the same binding region, suggesting that this site is highly flexible to accommodate different types of inhibitors. The drug specificity is due to multiple interactions not only with loop L5, but also with residues located in helices alpha2 and alpha3. These results suggest that tumour cells might develop resistance to Eg5 inhibitors, by expressing Eg5 point mutants that retain the enzyme activity, but prevent inhibition, a feature that is observed for certain tubulin inhibitors.     

Cooperative unfolding of a metastable serpin to a molten globule suggests a link between functional and folding energy landscapes

Alpha-1 antitrypsin (alpha(1)-AT) is a member of the serpin class of protease inhibitors, and folds to a metastable state rather than its thermodynamically most stable native state. Upon cleavage by a target protease, alpha(1)-AT undergoes a dramatic conformational change to a stable form, translocating the bound protease more than 70 A to form an inhibitory protease-serpin complex. Numerous mutagenesis studies on serpins have demonstrated the trade-off between the stability of the metastable state on the one hand and the inhibitory efficiency on the other. Studies of the equilibrium unfolding of serpins provide insight into this connection between structural plasticity and metastability. We studied equilibrium unfolding of wild-type alpha(1)-AT using hydrogen-deuterium/exchange mass spectrometry to characterize the structure and the stability of an equilibrium intermediate that was observed in low concentrations of denaturant in earlier studies. Our results show that the intermediate observed at low concentrations of denaturant has no protection from hydrogen-deuterium exchange, indicating a lack of stable structure. Further, differential scanning calorimetry of alpha(1)-AT at low concentrations of denaturant shows no heat capacity peak during thermal denaturation, indicating that the transition from the intermediate to the unfolded state is not a cooperative first-order-like phase transition.. Our results show that the unfolding of alpha(1)-AT involves a cooperative transition to a molten globule form, followed by a non-cooperative transition to a random-coil form as more guanidine is added. Thus, the entire alpha(1)-AT molecule consists of one cooperative structural unit rather than multiple structural domains with different stabilities. Furthermore, our results together with previous mutagenesis studies suggest a possible link between an equilibrium molten globule and a functional intermediate that may be populated during the protease inhibition.     

Metabolite induction via microorganism co-culture: A potential way to enhance chemical diversity for drug discovery

Microorganisms have a long track record as important sources of novel bioactive natural products, particularly in the field of drug discovery. While microbes have been shown to biosynthesize a wide array of molecules, recent advances in genome sequencing have revealed that such organisms have the potential to yield even more structurally diverse secondary metabolites. Thus, many microbial gene clusters may be silent under standard laboratory growth conditions. In the last ten years, several methods have been developed to aid in the activation of these cryptic biosynthetic pathways. In addition to the techniques that demand prior knowledge of the genome sequences of the studied microorganisms, several genome sequence-independent tools have been developed. One of these approaches is microorganism co-culture, involving the cultivation of two or more microorganisms in the same confined environment. Microorganism co-culture is inspired by the natural microbe communities that are omnipresent in nature. Within these communities, microbes interact through signaling or defense molecules. Such compounds, produced dynamically, are of potential interest as new leads for drug discovery. Microorganism co-culture can be achieved in either solid or liquid media and has recently been used increasingly extensively to study natural interactions and discover new bioactive metabolites. Because of the complexity of microbial extracts, advanced analytical methods (e.g., mass spectrometry methods and metabolomics) are key for the successful detection and identification of co-culture-induced metabolites.