Functional characterisation of urease accessory protein G (ureG) from potato

The activation of the nickel metalloenzyme urease is a complex process. In bacteria, several urease accessory proteins are essential for incorporation of nickel into the active centre of urease. Comparatively little is known about the activation process and the proteins involved in plants. We cloned five different cDNAs encoding isoforms of urease accessory protein G (ureG) in potato. The 5-coding region of these cDNAs is highly polymorphic within Solanum tuberosum ssp. tuberosum, containing mainly a simple sequence repeat encoding histidine and aspartate. Mapping on an ultrahigh-density map of the potato genome and Southern blot analysis showed that the isoforms arise from allelic differences of a single-copy gene which was located on chromosome 2. Expression analysis at the mRNA and protein levels indicated the presence of ureG in almost all tissues examined, consistent with the ubiquitous expression of urease. An attempt to correlate urease activity with ureG expression levels in different tissues was made. Allelic copies of ureG were expressed in a tissue-specific manner. UreG from potato and the Klebsiella aerogenes urease operon defective in bacterial ureG were co-expressed in Escherichia coli. The plant gene complements the K. aerogenes ureG mutation, demonstrating that it encodes a urease accessory protein and indicating a structural conservation between the plant and the bacterial urease activation complexes.     

Nickel binding properties of Helicobacter pylori UreF,an accessory protein in the nickel-based activation of urease

Helicobacter pylori UreF (HpUreF) is involved in the insertion of Ni²⁺ in the urease active site. The recombinant protein in solution is a dimer characterized by an extensive α-helical structure and a well-folded tertiary structure. HpUreF binds two Ni²⁺ ions per dimer, with a micromolar dissociation constant, as shown by calorimetry. X-ray absorption spectroscopy indicated that the Ni²⁺ ions reside in a five-coordinate pyramidal geometry comprising exclusively N/O-donor ligands derived from the protein, including one or two histidine imidazole and carboxylate ligands. Binding of Ni²⁺ does not affect the solution properties of the protein. Mutation to alanine of His229 and/or Cys231, a pair of residues located on the protein surface that interact with H. pylori UreD, altered the affinity of the protein for Ni²⁺. This result, complemented by the findings from X-ray absorption spectroscopy, indicates that the Ni²⁺ binding site involves His229, and that Cys231 has an indirect structural role in metal binding. An in vivo assay of urease activation demonstrated that H229A HpUreF, C231A HpUreF, and H229/C231 HpUreF are significantly less competent in this process, suggesting a role for a Ni²⁺ complex with UreF in urease maturation. This hypothesis was supported by calculations revealing the presence of a tunnel that joins the Cys-Pro-His metal binding site on UreG and an opening on the UreD surface, passing through UreF close to His229 and Cys231, in the structure of the H. pylori UreDFG complex. This tunnel could be used to transfer nickel into the urease active site during apoenzyme-to-holoenzyme activation.     

Crystal structure of a ferredoxin reductase for the CYP199A2 system from Rhodopseudomonas palustris

Cytochrome P450‐199A2 from Rhodopseudomonas palustris oxidizes para‐substituted benzoic acids and may play a role in lignin and aromatic acid degradation pathways in the bacterium. CYP199A2 has an associated [2Fe‐2S] ferredoxin, palustrisredoxin (Pux) but not a ferredoxin reductase. A genome search identified the palustrisredoxin reductase (PuR) gene. PuR was produced in Escherichia coli and shown to be a flavin‐dependent protein that supports efficient electron transfer from NADH to Pux, thus reconstituting CYP199A2 monooxygenase activity (k_cat = 37.9 s^–1 with 4‐methoxybenzoic acid). The reduction of Pux by PuR shows K_m = 4.2 μM and k_cat = 262 s^–1 in 50 mM Tris, pH 7.4. K_m is increased to 154 μM in the presence of 200 mM KCl, indicating the importance of ionic interactions in PuR/Pux binding. The crystal structure of PuR has been determined at 2.2 Å resolution and found to be closely related to that of other oxygenase‐coupled NADH‐dependent ferredoxin reductases. Residues on the surface that had been proposed to be involved in ferredoxin reductase‐ferredoxin binding are conserved in PuR. However, Lys328 in PuR lies over the FAD isoalloxazine ring and, together with His11 and Gln41, render the electrostatic potential of the surface more positive and may account for the greater involvement of electrostatic interactions in ferredoxin binding by PuR. Consistent with these observations the K328G mutation weakened Pux binding and virtually eliminated the dependence of PuR/Pux binding on salt concentration, thus confirming that the FAD si side surface in the vicinity of Lys328 is the ferredoxin binding site. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural characterization of the tunnels of Mycobacterium tuberculosis truncated hemoglobin N from molecular dynamics simulations

The structure of oxygenated trHbN from Mycobacterium tuberculosis shows an extended heme distal hydrogen‐bond network that includes Tyr33(B10), Gln58(E11), and the bound O₂. In addition, trHbN structure shows a network of hydrophobic cavities organized in two orthogonal branches. In the present work, the structure and the dynamics of oxygenated and deoxygenated trHbN in explicit water was investigated from 100 ns molecular dynamics (MD) simulations. Results show that, depending on the presence or the absence of a coordinated O₂, the Tyr33(B10) and Gln58(E11) side chains adopt two different configurations in concert with hydrogen bond network rearrangement. In addition, our data indicate that Tyr33(B10) and Gln58(E11) control the dynamics of Phe62(E15). In deoxy‐trHbN, Phe62(E15) is restricted to one conformation. Upon O₂ binding, the conformation of Gln58(E11) changes and residue Phe62(E15) fluctuates between two conformations. We also conducted a systematic study of trHbN tunnels by analyzing thousands of MD snapshots with CAVER. The results show that tunnel formation is the result of the dynamic reshaping of short‐lived hydrophobic cavities. The analyses indicate that the presence of these cavities is likely linked to the rigid structure of trHbN and also reveal two tunnels, EH and BE, that link the protein surface to the buried distal heme pocket and not present in the crystallographic structure. The cavities are sufficiently large to accomodate and store ligands. Tunnel dynamics in trHbN was found to be controlled by the side‐chain conformation of the Tyr33(B10), Gln58(E11), and Phe62(E15) residues. Importantly, in contrast to recently published works, our extensive systematic studies show that the presence or absence of a coordinated dioxygen does not control the opening of the long tunnel but rather the opening of the EH tunnel. In addition, the data lead to new and distinctly different conclusion on the impact of the Phe62(E15) residue on trHbN tunnels. We propose that the EH and the long tunnels are used for apolar ligands storage. The trajectories bring important new structural insights related to trHbN function and to ligand diffusion in proteins. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Adherence protects the binding sites of Helicobacter pylori urease from acid-induced damage

Colonization by Helicobacter pylori partly depends on acid-dependent adherence by urease to gastric mucin. To further verify the relevance of urease adherence to colonization, the influence of acidity on the binding sites of H. pylori urease was investigated. When enzyme-based in vitro ligand capture assays were used, the effect of acidity on the binding site of H. pylori urease was determined against a backdrop medium consisting of acidic buffers simulating the luminal side of gastric mucus. A high degree of stability was exhibited by adherent urease, suggesting a pivotal role by the denatured enzyme in the persistence of the bacterium within the acidified compartment of gastric mucus.     

Subunit-specific backbone NMR assignments of a 64 kDa trp repressor/DNA complex: A role for N-terminal residues in tandem binding

Deuterium decoupled, triple resonance NMR spectroscopy was used to analyze complexes of ²H,¹⁵N,¹³C labelled intact and (des2–7) trp repressor (2–7 trpR) from E. coli bound in tandem to an idealized 22 basepair trp operator DNA fragment and the corepressor 5-methyltryptophan. The DNA sequence used here binds two trpR dimers in tandem resulting in chemically nonequivalent environments for the two subunits of each dimer. Sequence- and subunit-specific NMR resonance assignments were made for backbone ¹HN, ¹⁵N, ¹³C positions in both forms of the protein and for¹³ C in the intact repressor. The differences in backbone chemical shifts between the two subunits within each dimer of 2–7 trpR reflect dimer-dimer contacts involving the helix-turn-helix domains and N-terminal residues consistent with a previously determined crystal structure [Lawson and Carey (1993) Nature, 366, 178–182]. Comparison of the backbone chemical shifts of DNA-bound 2–7 trpR with those of DNA-bound intact trpR reveals significant changes for those residues involved in N-terminal-mediated interactions observed in the crystal structure. In addition, our solution NMR data contain three sets of resonances for residues 2–12 in intact trpR suggesting that the N-terminus has multiple conformations in the tandem complex. Analysis of C chemical shifts using a chemical shift index (CSI) modified for deuterium isotope effects has allowed a comparison of the secondary structure of intact and 2–7 tprR. Overall these data demonstrate that NMR backbone chemical shift data can be readily used to study specific structural details of large protein complexes.     

A model-based proposal for the role of UreF as a GTPase-activating protein in the urease active site biosynthesis

UreF is a protein that plays a role in the in vivo urease activation as a chaperone involved in the insertion of two Ni(2+) ions in the apo-urease active site. The molecular details of this process are unknown. In the absence of any molecular information on the UreF protein class, and as a step toward the comprehension of the relationships between UreF function and structure, we applied a structural modeling approach to infer useful biochemical knowledge on Bacillus pasteurii UreF (BpUreF). Similarity searches and multiple alignment of UreF protein sequences indicated that this class of proteins has a low homology with proteins of known structure. Fold recognition methods were therefore used to identify useful protein structural templates to model the structure of BpUreF. In particular, the templates belong to the class of GTPase-activating proteins. Modeling of BpUreF based on these templates was performed using the program MODELLER. The structure validation yielded good statistics, indicating that the model is plausible. This result suggests a role for UreF in urease active site biosynthesis as a regulator of the activity of UreG, a small G protein involved in the in vivo apo-urease activation process and established to catalyze GTP hydrolysis.     

Zn2+-linked dimerization of UreG from Helicobacter pylori, a chaperone involved in nickel trafficking and urease activation

The biosynthesis of the active metal-bound form of the nickel-dependent enzyme urease involves the formation of a lysine-carbamate functional group concomitantly with the delivery of two Ni(2+) ions into the precast active site of the apoenzyme and with GTP hydrolysis. In the urease system, this role is performed by UreG, an accessory protein belonging to the group of homologous P-loop GTPases, often required to complete the biosynthesis of nickel-enzymes. This study is focused on UreG from Helicobacter pylori (HpUreG), a bacterium responsible for gastric ulcers and cancer, infecting large part of the human population, and for which urease is a fundamental virulence factor. The soluble HpUreG was expressed in E. coli and purified to homogeneity. On-line size exclusion chromatography and light scattering indicated that apo-HpUreG exists as a monomer in solution. Circular dichroism, which demonstrated the presence of a well-defined secondary structure, and NMR spectroscopy, which revealed a large number of residues that appear structured on the basis of their backbone amide proton chemical shift dispersion, indicated that, at variance with other UreG proteins so far characterized, this protein is significantly folded in solution. The amino acid sequence of HpUreG is 29% identical to that of HypB from Methanocaldococcus jannaschii, a dimeric zinc-binding GTPase involved in the in vivo assembly of [Ni,Fe]-hydrogenase. A homology-based molecular model of HpUreG was calculated, which allowed us to identify structural and functional features of the protein. Isothermal titration microcalorimetry demonstrated that HpUreG specifically binds 0.5 equivalents of Zn(2+) per monomer (K(d) = 0.33 +/- 0.03 microM), whereas it has 20-fold lower affinity for Ni(2+) (K(d) = 10 +/- 1 microM). Zinc ion binding (but not Ni(2+) binding) causes protein dimerization, as confirmed using light scattering measurements. The structural rearrangement occurring upon Zn(2+)-binding and consequent dimerization was evaluated using circular dichroism and fluorescence spectroscopy. Fully conserved histidine and cysteine residues were identified and their role in zinc binding was verified by site-directed mutagenesis and microcalorimetry. The results are analyzed and discussed with respect to analogous examples of GTPases in nickel metabolism.     

Predicting permanent and transient protein–protein interfaces

Protein–protein interactions (PPIs) are involved in diverse functions in a cell. To optimize functional roles of interactions, proteins interact with a spectrum of binding affinities. Interactions are conventionally classified into permanent and transient, where the former denotes tight binding between proteins that result in strong complexes, whereas the latter compose of relatively weak interactions that can dissociate after binding to regulate functional activity at specific time point. Knowing the type of interactions has significant implications for understanding the nature and function of PPIs. In this study, we constructed amino acid substitution models that capture mutation patterns at permanent and transient type of protein interfaces, which were found to be different with statistical significance. Using the substitution models, we developed a novel computational method that predicts permanent and transient protein binding interfaces (PBIs) in protein surfaces. Without knowledge of the interacting partner, the method uses a single query protein structure and a multiple sequence alignment of the sequence family. Using a large dataset of permanent and transient proteins, we show that our method, BindML+, performs very well in protein interface classification. A very high area under the curve (AUC) value of 0.957 was observed when predicted protein binding sites were classified. Remarkably, near prefect accuracy was achieved with an AUC of 0.991 when actual binding sites were classified. The developed method will be also useful for protein design of permanent and transient PBIs. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Using collections of structural models to predict changes of binding affinity caused by mutations in protein–protein interactions

Protein–protein interactions (PPIs) in all the molecular aspects that take place both inside and outside cells. However, determining experimentally the structure and affinity of PPIs is expensive and time consuming. Therefore, the development of computational tools, as a complement to experimental methods, is fundamental. Here, we present a computational suite: MODPIN, to model and predict the changes of binding affinity of PPIs. In this approach we use homology modeling to derive the structures of PPIs and score them using state‐of‐the‐art scoring functions. We explore the conformational space of PPIs by generating not a single structural model but a collection of structural models with different conformations based on several templates. We apply the approach to predict the changes in free energy upon mutations and splicing variants of large datasets of PPIs to statistically quantify the quality and accuracy of the predictions. As an example, we use MODPIN to study the effect of mutations in the interaction between colicin endonuclease 9 and colicin endonuclease 2 immune protein from Escherichia coli. Finally, we have compared our results with other state‐of‐art methods.     

Functional protein domains from the thermally driven motion of polypeptide chains: A proposal

It is proposed that the thermally driven motion of certain polypeptide chains, including those that are part of an otherwise stable folded protein, produces time-averaged three-dimensional domains that confer unique functions to a protein. These domains may be controlled by collapsing the polypeptide into an enthalpically favored structure, or extending it into an entropically dominated form. In the extended form, these domains occupy a relatively large space, which may be used to regulate protein–protein interactions and confer mechanical properties to proteins. This “entropic bristle” model makes several predictions about the structure and properties of these domains, and the predictions are used to reevaluate a range of biophysical studies on proteins. The outcome of the analysis suggests that the entropic bristle can be used to explain a wide range of disparate and apparently unrelated experimental observations. Proteins 32:223–228, 1998. © 1998 Wiley-Liss, Inc.     

Benchmarking of structure refinement methods for protein complex models

Protein structure docking is the process in which the quaternary structure of a protein complex is predicted from individual tertiary structures of the protein subunits. Protein docking is typically performed in two main steps. The subunits are first docked while keeping them rigid to form the complex, which is then followed by structure refinement. Structure refinement is crucial for a practical use of computational protein docking models, as it is aimed for correcting conformations of interacting residues and atoms at the interface. Here, we benchmarked the performance of eight existing protein structure refinement methods in refinement of protein complex models. We show that the fraction of native contacts between subunits is by far the most straightforward metric to improve. However, backbone dependent metrics, based on the Root Mean Square Deviation proved more difficult to improve via refinement.     

Target of rapamycin FATC domain as a general membrane anchor: The FKBP‐12 like domain of FKBP38 as a case study

Increased efforts have been undertaken to better understand the formation of signaling complexes at cellular membranes. Since the preparation of proteins containing a transmembrane domain or a prenylation motif is generally challenging an alternative membrane anchoring unit that is easy to attach, water‐soluble and binds to different membrane mimetics would find broad application. The 33‐residue long FATC domain of yeast TOR1 (y1fatc) fulfills these criteria and binds to neutral and negatively charged micelles, bicelles, and liposomes. As a case study, we fused it to the FKBP506‐binding region of the protein FKBP38 (FKBP38‐BD) and used ¹H–¹⁵N NMR spectroscopy to characterize localization of the chimeric protein to micelles, bicelles, and liposomes. Based on these and published data for y1fatc, its use as a C‐terminally attachable membrane anchor for other proteins is compatible with a wide range of buffer conditions (pH circa 6–8.5, NaCl 0 to >150 mM, presence of reducing agents, different salts such as MgCl₂ and CaCl₂). The high water‐solubility of y1fatc enables its use for titration experiments against a membrane‐localized interaction partner of the fused target protein. Results from studies with peptides corresponding to the C‐terminal 17–11 residues of the 33‐residue long domain by 1D ¹H NMR and CD spectroscopy indicate that they still can interact with membrane mimetics. Thus, they may be used as membrane anchors if the full y1fatc sequence is disturbing or if a chemically synthesized y1fatc peptide shall be attached by native chemical ligation, for example, unlabeled peptide to ¹⁵N‐labeled target protein for NMR studies.     

A simple contact mapping algorithm for identifying potential peptide mimetics in protein–protein interaction partners

A simple, static contact mapping algorithm has been developed as a first step at identifying potential peptide biomimetics from protein interaction partner structure files. This rapid and simple mapping algorithm, “OpenContact” provides screened or parsed protein interaction files based on specified criteria for interatomic separation distances and interatomic potential interactions. The algorithm, which uses all‐atom Amber03 force field models, was blindly tested on several unrelated cases from the literature where potential peptide mimetics have been experimentally developed to varying degrees of success. In all cases, the screening algorithm efficiently predicted proposed or potential peptide biomimetics, or close variations thereof, and provided complete atom‐atom interaction data necessary for further detailed analysis and drug development. In addition, we used the static parsing/mapping method to develop a peptide mimetic to the cancer protein target, epidermal growth factor receptor. In this case, secondary, loop structure for the peptide was indicated from the intra‐protein mapping, and the peptide was subsequently synthesized and shown to exhibit successful binding to the target protein. The case studies, which all involved experimental peptide drug advancement, illustrate many of the challenges associated with the development of peptide biomimetics, in general. Proteins 2014; 82:2253–2262. © 2014 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Knitting and untying the protein network: Modulation of protein ensembles as a therapeutic strategy

Proteins constitute the working machinery and structural support of all organisms. In performing a given function, they must adopt highly specific structures that can change with their level of activity, often through the direct or indirect action of other proteins. Indeed, proteins typically function within an ensemble, rather than individually. Hence, they must be sufficiently flexible to interact with each other and execute diverse tasks. The discovery that errors within these groups can ultimately cause disease has led to a paradigm shift in drug discovery, from an emphasis on single protein targets to a holistic approach whereby entire ensembles are targeted.     

Importance of the interaction protein–protein of the CaM–PDE1A and CaM–MLCK complexes in the development of new anti‐CaM drugs

Protein–protein interactions play central roles in physiological and pathological processes. The bases of the mechanisms of drug action are relevant to the discovery of new therapeutic targets. This work focuses on understanding the interactions in protein–protein–ligands complexes, using proteins calmodulin (CaM), human calcium/calmodulin‐dependent 3′,5′‐cyclic nucleotide phosphodiesterase 1A active human (PDE1A), and myosin light chain kinase (MLCK) and ligands αII–spectrin peptide (αII–spec), and two inhibitors of CaM (chlorpromazine (CPZ) and malbrancheamide (MBC)). The interaction was monitored with a fluorescent biosensor of CaM (hCaM M124C–mBBr). The results showed changes in the affinity of CPZ and MBC depending on the CaM–protein complex under analysis. For the Ca²⁺–CaM, Ca²⁺–CaM–PDE1A, and Ca²⁺–CaM–MLCK complexes, CPZ apparent dissociation constants (K_ds) were 1.11, 0.28, and 0.55 μM, respectively; and for MBC K_ds were 1.43, 1.10, and 0.61 μM, respectively. In competition experiments the addition of calmodulin binding peptide 1 (αII–spec) to Ca²⁺–hCaM M124C–mBBr quenched the fluorescence (K_d = 2.55 ± 1.75 pM) and the later addition of MBC (up to 16 μM) did not affect the fluorescent signal. Instead, the additions of αII–spec to a preformed Ca²⁺–hCaM M124C–mBBr–MBC complex modified the fluorescent signal. However, MBC was able to displace the PDE1A and MLCK from its complex with Ca²⁺–CaM. In addition, docking studies were performed for all complexes with both ligands showing an excellent correlation with experimental data. These experiments may help to explain why in vivo many CaM drugs target prefer only a subset of the Ca²⁺–CaM regulated proteins and adds to the understanding of molecular interactions between protein complexes and small ligands. Copyright © 2013 John Wiley & Sons, Ltd.     

Functional evolution within a protein superfamily

The ability to predict and characterize distributions of reactivities over families and even superfamilies of proteins opens the door to an array of analyses regarding functional evolution. In this article, insights into functional evolution in the Kazal inhibitor superfamily are gained by analyzing and comparing predicted association free energy distributions against six serine proteinases, over a number of groups of inhibitors: all possible Kazal inhibitors, natural avian ovomucoid first and third domains, and sets of Kazal inhibitors with statistically weighted combinations of residues. The results indicate that, despite the great hypervariability of residues in the 10 proteinase-binding positions, avian ovomucoid third domains evolved to inhibit enzymes similar to the six enzymes selected, whereas the orthologous first domains are not inhibitors of these enzymes on purpose. Hypervariability arises because of similarity in energetic contribution from multiple residue types; conservation is in terms of functionality, with "good" residues, which make positive or less deleterious contributions to the binding, selected more frequently, and yielding overall the same distributional characteristics. Further analysis of the distributions indicates that while nature did optimize inhibitor strength, the objective may not have been the strongest possible inhibitor against one enzyme but rather an inhibitor that is relatively strong against a number of enzymes.