Predicting 3D structures of transient protein-protein complexes by homology

The paper reports a homology based approach for predicting the 3D structures of full length hetero protein complexes. We have created a database of templates that includes structures of hetero protein-protein complexes as well as domain-domain structures (), which allowed us to expand the template pool up to 418 two-chain entries (at 40% sequence identity). Two protocols were tested-a protocol based on position specific Blast search (Protocol-I) and a protocol based on structural similarity of monomers (Protocol-II). All possible combinations of two monomers (350,284 pairs) in the ProtCom database were subjected to both protocols to predict if they form complexes. The predictions were benchmarked against the ProtCom database resulting to false-true positives ratios of approximately 5:1 and approximately 7:1 and recovery of 19% and 86%, respectively for protocols I and II. From 350,284 trials Protocol-I made only approximately 500 wrong predictions resulting to 0.5% error. In addition, though it was shown that artificially created domain-domain structures can in principle be good templates for modeling full length protein complexes, more sensitive methods are needed to detect homology relations. The quality of the models was assessed using two different criteria such as interfacial residues and overall RMSD. It was found that there is no correlation between these two measures. In many cases the interface residues were predicted correctly, but the overall RMSD was over 6 A and vice versa.     

Protein-protein complex structure predictions by multimeric threading and template recombination

The total number of protein-protein complex structures currently available in the Protein Data Bank (PDB) is six times smaller than the total number of tertiary structures in the PDB, which limits the power of homology-based approaches to complex structure modeling. We present a threading-recombination approach, COTH, to boost the protein complex structure library by combining tertiary structure templates with complex alignments. The query sequences are first aligned to complex templates using a modified dynamic programming algorithm, guided by ab initio binding-site predictions. The monomer alignments are then shifted to the multimeric template framework by structural alignments. COTH was tested on 500 nonhomologous dimeric proteins, which can successfully detect correct templates for 50% of the cases after homologous templates are excluded, which significantly outperforms conventional homology modeling algorithms. It also shows a higher accuracy in interface modeling than rigid-body docking of unbound structures from ZDOCK although with lower coverage. These data demonstrate new avenues to model complex structures from nonhomologous templates.     

VITAL NMR: using chemical shift derived secondary structure information for a limited set of amino acids to assess homology model accuracy

Homology modeling is a powerful tool for predicting protein structures, whose success depends on obtaining a reasonable alignment between a given structural template and the protein sequence being analyzed. In order to leverage greater predictive power for proteins with few structural templates, we have developed a method to rank homology models based upon their compliance to secondary structure derived from experimental solid-state NMR (SSNMR) data. Such data is obtainable in a rapid manner by simple SSNMR experiments (e.g., ¹³C–¹³C 2D correlation spectra). To test our homology model scoring procedure for various amino acid labeling schemes, we generated a library of 7,474 homology models for 22 protein targets culled from the TALOS+/SPARTA+ training set of protein structures. Using subsets of amino acids that are plausibly assigned by SSNMR, we discovered that pairs of the residues Val, Ile, Thr, Ala and Leu (VITAL) emulate an ideal dataset where all residues are site specifically assigned. Scoring the models with a predicted VITAL site-specific dataset and calculating secondary structure with the Chemical Shift Index resulted in a Pearson correlation coefficient (−0.75) commensurate to the control (−0.77), where secondary structure was scored site specifically for all amino acids (ALL 20) using STRIDE. This method promises to accelerate structure procurement by SSNMR for proteins with unknown folds through guiding the selection of remotely homologous protein templates and assessing model quality.     

Nitrogen Oxyanion-dependent Dissociation of a Two-component Complex That Regulates Bacterial Nitrate Assimilation

Nitrogen is an essential nutrient for growth and is readily available to microbes in many environments in the form of ammonium and nitrate. Both ions are of environmental significance due to sustained use of inorganic fertilizers on agricultural soils. Diverse species of bacteria that have an assimilatory nitrate/nitrite reductase system (NAS) can use nitrate or nitrite as the sole nitrogen source for growth when ammonium is limited. In Paracoccus denitrificans, the pathway-specific two-component regulator for NAS expression is encoded by the nasT and nasS genes. Here, we show that the putative RNA-binding protein NasT is a positive regulator essential for expression of the nas gene cluster (i.e. nasABGHC). By contrast, a nitrogen oxyanion-binding sensor (NasS) is required for nitrate/nitrite-responsive control of nas gene expression. The NasS and NasT proteins co-purify as a stable heterotetrameric regulatory complex, NasS-NasT. This protein-protein interaction is sensitive to nitrate and nitrite, which cause dissociation of the NasS-NasT complex into monomeric NasS and an oligomeric form of NasT. NasT has been shown to bind the leader RNA for nasA. Thus, upon liberation from the complex, the positive regulator NasT is free to up-regulate nas gene expression.     

Autoinduction Specificities of the Lantibiotics Subtilin and Nisin

The biosynthesis of the lantibiotics subtilin and nisin is regulated by autoinduction via two-component systems. Although subtilin is structurally closely related to nisin and contains the same lanthionine ring structure, both lantibiotics specifically autoinduce their biosynthesis. Subtilin and also the subtilin-like lantibiotics entianin and ericin autoinduce the two-component system SpaRK of Bacillus subtilis, whereas the biosynthesis of nisin is autoinduced via the two-component system NisRK of Lactococcus lactis. Autoinduction is highly specific for the respective lantibiotic and therefore of major importance for the functional expression of genetically engineered subtilin-like lantibiotics. To identify the structural features required for subtilin autoinduction, subtilin-nisin hybrids and specific point mutations of amino acid position 1 were generated. For subtilin autoinduction, the N-terminal tryptophan is the most important for full SpaK activation. The failure of subtilin to autoinduce the histidine kinase NisK mainly depends on the N-terminal tryptophan, as its single exchange to the aliphatic amino acid residues isoleucine, leucine, and valine provided NisK autoinduction. In addition, the production of subtilin variants which did not autoinduce their own biosynthesis could be rescued upon heterologous coexpression in B. subtilis DSM15029 by the autoinducing subtilin-like lantibiotic entianin.     

Universally applicable methods for monitoring response regulator aspartate phosphorylation both in vitro and in vivo using Phos-tag-based reagents

Recent development of the phosphate chelator, Phos-tag, together with Phos-tag pendant reagents, has provided new methods for detection of phosphorylated serine, threonine, tyrosine, and histidine residues in phosphoproteins. We have investigated the use of Phos-tag for detection and quantification of phospho-aspartate in response regulator proteins that function within two-component signaling systems. Alternative methods are especially important, because the labile nature of the acylphosphate bond in response regulator proteins has restricted the application of many traditional methods of phosphoprotein analysis. We demonstrate that Phos-tag gel stain can be used to detect phospho-Asp in response regulators and that Phos-tag acrylamide gel electrophoresis can be used to separate phosphorylated and unphosphorylated forms of response regulator proteins. The latter method, coupled to Western blot analysis, enables detection of specific phosphorylated proteins in complex mixtures such as cell lysates. Standards of phosphorylated proteins can be used to correct for hydrolysis of the labile phospho-Asp bond that invariably occurs during analysis. We have employed Phos-tag methods to characterize the phosphorylation state of the Escherichia coli response regulator PhoB both in vitro, using purified protein, and in vivo, by analyzing lysates of cells grown under different conditions of induction of the PhoR/PhoB phosphate assimilation pathway.     

Protein packing: dependence on protein size, secondary structure and amino acid composition

We have used the occluded surface algorithm to estimate the packing of both buried and exposed amino acid residues in protein structures. This method works equally well for buried residues and solvent-exposed residues in contrast to the commonly used Voronoi method that works directly only on buried residues. The atomic packing of individual globular proteins may vary significantly from the average packing of a large data set of globular proteins. Here, we demonstrate that these variations in protein packing are due to a complex combination of protein size, secondary structure composition and amino acid composition. Differences in protein packing are conserved in protein families of similar structure despite significant sequence differences. This conclusion indicates that quality assessments of packing in protein structures should include a consideration of various parameters including the packing of known homologous proteins. Also, modeling of protein structures based on homologous templates should take into account the packing of the template protein structure.     

How to choose templates for modeling of protein complexes: Insights from benchmarking template-based docking

Comparative docking is based on experimentally determined structures of protein-protein complexes (templates), following the paradigm that proteins with similar sequences and/or structures form similar complexes. Modeling utilizing structure similarity of target monomers to template complexes significantly expands structural coverage of the interactome. Template-based docking by structure alignment can be performed for the entire structures or by aligning targets to the bound interfaces of the experimentally determined complexes. Systematic benchmarking of docking protocols based on full and interface structure alignment showed that both protocols perform similarly, with top 1 docking success rate 26%. However, in terms of the models' quality, the interface-based docking performed marginally better. The interface-based docking is preferable when one would suspect a significant conformational change in the full protein structure upon binding, for example, a rearrangement of the domains in multidomain proteins. Importantly, if the same structure is selected as the top template by both full and interface alignment, the docking success rate increases 2-fold for both top 1 and top 10 predictions. Matching structural annotations of the target and template proteins for template detection, as a computationally less expensive alternative to structural alignment, did not improve the docking performance. Sophisticated remote sequence homology detection added templates to the pool of those identified by structure-based alignment, suggesting that for practical docking, the combination of the structure alignment protocols and the remote sequence homology detection may be useful in order to avoid potential flaws in generation of the structural templates library.     

Molecular Modeling Studies of the Interaction Between Plasmodium falciparum HslU and HslV Subunits

Abstract

The PfHslUV, a Plasmodium falciparum homolog of prokaryotic HslUV systems, is a newly identified drug target. The HslUV complex is an assembly of Heat Shock Locus gene products U and V. The formation of complete complex is essential for the proteasome to carry out its biochemical and physiological role in the parasite, namely to degrade specific target proteins in an ATP-dependent chaperone assisted manner. PfHslV subunit, a protease, exhibits increased proteolytic activity in the presence of PfHslU, the subunit believed to be responsible for allosteric activation of PfHslV. In the present work, we have employed computational methods to simulate the interaction of PfHslU and PfHslV subunits. We have used three methods—namely homology modeling, molecular docking and computational alanine scanning to model the complex, to predict the binding mode of PfHslU-V interaction and to predict the binding-energy hot-spots in protein-protein interface, respectively. The three dimensional models of PfHslV and PfHslU have been generated using MODELLER, based on the crystal structures of prokaryotic HslUV complex as templates. The modeled structures were docked using PatchDock, a geometry-based molecular docking algorithm. Finally, a three-dimensional PfHslUV complex model was generated that helped in comparing protein-protein interface characteristics with that of crystal structures of prokaryotic HslUV. Further, computational alanine scanning analysis of the generated complex was performed to calculate the binding free energy changes (ΔΔG_bind), which helped in identifying residues crucial for PfHslU and PfHslV interactions.     

Regulation of nisin biosynthesis and immunity in Lactococcus lactis 6F3.

The biosynthetic genes of the nisin-producing strain Lactococcus lactis 6F3 are organized in an operon-like structure starting with the structural gene nisA followed by the genes nisB, nisT, and nisC, which are probably involved in chemical modification and secretion of the prepeptide (G. Engelke, Z. Gutowski-Eckel, M. Hammelmann, and K.-D. Entian, Appl. Environ. Microbiol. 58:3730-3743, 1992). Subcloning of an adjacent 5-kb downstream region revealed additional genes involved in nisin biosynthesis. The gene nisI, which encodes a lipoprotein, causes increased immunity after its transformation into nisin-sensitive L. lactis MG1614. It is followed by the gene nisP, coding for a subtilisin-like serine protease possibly involved in processing of the secreted leader peptide. Adjacent to the 3' end of nisP the genes nisR and nisK were identified, coding for a regulatory protein and a histidine kinase, showing marked similarities to members of the OmpR/EnvZ-like subgroup of two-component regulatory systems. The deduced amino acid sequences of nisR and nisK exhibit marked similarities to SpaR and SpaK, which were recently identified as the response regulator and the corresponding histidine kinase of subtilin biosynthesis. By using antibodies directed against the nisin prepeptide and the NisB protein, respectively, we could show that nisin biosynthesis is regulated by the expression of its structural and biosynthetic genes. Prenisin expression starts in the exponential growth phase and precedes that of the NisB protein by approximately 30 min. Both proteins are expressed to a maximum in the stationary growth phase.     

Use of oligonucleotide probes to identify members of two-component regulatory systems in Xanthomonas campestris pathovar campestris

Summary Two-component regulatory systems comprising a sensor and a regulator protein, both with highly conserved amino acid domains, and commonly genetically linked, have been described in a range of bacterial species and are involved in sensing environmental stimuli. We used two oligonucleotide probes matching the postulated coding regions for domains of sensor and regulator proteins respectively in Xanthomonas campestris pathovar campestris (Xcc) to identify possible two-component regulatory systems in Xcc. Two different fragments of Xcc DNA with homology to both of these probes were cloned. The DNA sequence of part of one of these fragments encompassed a potential open reading frame (ORF), the predicted amino acid sequence of which had extensive homology with regulator proteins of two-component regulatory systems. Analysis of the predicted amino acid sequence for the 3 end of an adjacent ORF revealed a very high level of homology with the C-terminal end of sensor proteins. Strains of Xcc with Tn5-induced mutations in the regulator gene were affected in extracellular polysaccharide production, and also in resistance to salt and chloramphenicol. No effects of mutation in the second clone were observed.     

The arcB gene of Escherichia coli encodes a sensor-regulator protein for anaerobic repression of the arc modulon

The arcA (dye) and arcB genes of Escherichia coli are responsible for anaerobic repression of target operons and regulons of aerobic function (the arc modulon). The amino acid sequence of ArcA (Dye) indicated that it is the regulator protein of a two-component control system. Here we show that ArcB is a membrane sensor protein on the basis of its deduced amino acid sequence (778 residues), hydropathicity profile, and cellular distribution. On the carboxyl end of the ArcB sequence there is an additional domain showing homology with conserved regions of regulator proteins. Deletion into this domain destroyed ArcB function. ArcB conserved a histidine residue for autophosphorylation of the sensor proteins, and aspartic residues important for the regulator proteins.     

NMR-derived topology of a GFP-photoprotein energy transfer complex

Förster resonance energy transfer within a protein-protein complex has previously been invoked to explain emission spectral modulation observed in several bioluminescence systems. Here we present a spatial structure of a complex of the Ca²⁺-regulated photoprotein clytin with its green-fluorescent protein (cgGFP) from the jellyfish Clytia gregaria, and show that it accounts for the bioluminescence properties of this system in vitro. We adopted an indirect approach of combining x-ray crystallography determined structures of the separate proteins, NMR spectroscopy, computational docking, and mutagenesis. Heteronuclear NMR spectroscopy using variously ¹⁵N,¹³C,²H-enriched proteins enabled assignment of backbone resonances of more than 94% of the residues of both proteins. In a mixture of the two proteins at millimolar concentrations, complexation was inferred from perturbations of certain ¹H-¹⁵N HSQC-resonances, which could be mapped to those residues involved at the interaction site. A docking computation using HADDOCK was employed constrained by the sites of interaction, to deduce an overall spatial structure of the complex. Contacts within the clytin-cgGFP complex and electrostatic complementarity of interaction surfaces argued for a weak protein-protein complex. A weak affinity was also observed by isothermal titration calorimetry (K_D = 0.9 mm). Mutation of clytin residues located at the interaction site reduced the degree of protein-protein association concomitant with a loss of effectiveness of cgGFP in color-shifting the bioluminescence. It is suggested that this clytin-cgGFP structure corresponds to the transient complex previously postulated to account for the energy transfer effect of GFP in the bioluminescence of aequorin or Renilla luciferase.     

NMR and homology modeling studies of copper(II)-halocyanin from Natronobacterium pharaonis bacteria

Halocyanin from the haloalkaliphilic archaean Natronobacterium pharaonis is a peripheral membrane type 1 blue copper protein with a single polypeptide chain of 163 amino acid residues. Halocyanin participates as putative electron carrier protein associated to an electron acceptor role for a terminal oxidase and has the lowest redox potential value reported to date for a BCP. NMR studies and homology modeling calculations were performed to evaluate the electronic properties of Cu(II)-halocyanin from Natronobacterium pharaonis. The copper coordination site properties of Cu(II)-halocyanin are discussed. The ¹H NMR spectra, isotropic chemical shifts and relaxation times for halocyanin are compared with those of other BCPs such as azurin, amicyanin, plastocyanin and stellacyanin. The wild-type Cu(II)-halocyanin presents almost the same ¹H NMR spectra in comparison with Cu(II)-plastocyanin as expected from a similar coordination symmetry. However, minor differences were found. In order to get some insight on these differences, a computational model for Cu(II)-halocyanin from N. pharaonis was built. Model is based on sequential homology of halocyanin with two different families of proteins: plastocyanins and pseudoazurins. Homology modeling was performed using two different structural templates and copper ion was added for further refinement of the coordination site. Proposed structure was in good agreement with NMR experimental information and is the first three-dimensional model reported to date of an halocyanin. Small differences were found in the copper coordination site with respect to other BCP with known structure. This work is also an interesting example of expertise-driven homology modeling across different protein families.     

Accuracy of Protein-Protein Binding Sites in High-Throughput Template-Based Modeling

The accuracy of protein structures, particularly their binding sites, is essential for the success of modeling protein complexes. Computationally inexpensive methodology is required for genome-wide modeling of such structures. For systematic evaluation of potential accuracy in high-throughput modeling of binding sites, a statistical analysis of target-template sequence alignments was performed for a representative set of protein complexes. For most of the complexes, alignments containing all residues of the interface were found. The full interface alignments were obtained even in the case of poor alignments where a relatively small part of the target sequence (as low as 40%) aligned to the template sequence, with a low overall alignment identity (<30%). Although such poor overall alignments might be considered inadequate for modeling of whole proteins, the alignment of the interfaces was strong enough for docking. In the set of homology models built on these alignments, one third of those ranked 1 by a simple sequence identity criteria had RMSD<5 Å, the accuracy suitable for low-resolution template free docking. Such models corresponded to multi-domain target proteins, whereas for single-domain proteins the best models had 5 Å<RMSD<10 Å, the accuracy suitable for less sensitive structure-alignment methods. Overall, ∼50% of complexes with the interfaces modeled by high-throughput techniques had accuracy suitable for meaningful docking experiments. This percentage will grow with the increasing availability of co-crystallized protein-protein complexes.     

TASSER_low‐zsc: An approach to improve structure prediction using low z‐score–ranked templates

In a variety of threading methods, often poorly ranked (low z‐score) templates have good alignments. Here, a new method, TASSER_low‐zsc that identifies these low z‐score–ranked templates to improve protein structure prediction accuracy, is described. The approach consists of clustering of threading templates by affinity propagation on the basis of structural similarity (thread_cluster) followed by TASSER modeling, with final models selected by using a TASSER_QA variant. To establish the generality of the approach, templates provided by two threading methods, SP³ and SPARKS², are examined. The SP³ and SPARKS² benchmark datasets consist of 351 and 357 medium/hard proteins (those with moderate to poor quality templates and/or alignments) of length ≤250 residues, respectively. For SP³ medium and hard targets, using thread_cluster, the TM‐scores of the best template improve by ～4 and 9% over the original set (without low z‐score templates) respectively; after TASSER modeling/refinement and ranking, the best model improves by ～7 and 9% over the best model generated with the original template set. Moreover, TASSER_low‐zsc generates 22% (43%) more foldable medium (hard) targets. Similar improvements are observed with low‐ranked templates from SPARKS². The template clustering approach could be applied to other modeling methods that utilize multiple templates to improve structure prediction. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

The Inner Membrane Complex Sub-compartment Proteins Critical for Replication of the Apicomplexan Parasite Toxoplasma gondii Adopt a Pleckstrin Homology Fold

Toxoplasma gondii, an apicomplexan parasite prevalent in developed nations, infects up to one-third of the human population. The success of this parasite depends on several unique structures including an inner membrane complex (IMC) that lines the interior of the plasma membrane and contains proteins important for gliding motility and replication. Of these proteins, the IMC sub-compartment proteins (ISPs) have recently been shown to play a role in asexual T. gondii daughter cell formation, yet the mechanism is unknown. Complicating mechanistic characterization of the ISPs is a lack of sequence identity with proteins of known structure or function. In support of elucidating the function of ISPs, we first determined the crystal structures of representative members TgISP1 and TgISP3 to a resolution of 2.10 and 2.32 Å, respectively. Structural analysis revealed that both ISPs adopt a pleckstrin homology fold often associated with phospholipid binding or protein-protein interactions. Substitution of basic for hydrophobic residues in the region that overlays with phospholipid binding in related pleckstrin homology domains, however, suggests that ISPs do not retain phospholipid binding activity. Consistent with this observation, biochemical assays revealed no phospholipid binding activity. Interestingly, mapping of conserved surface residues combined with crystal packing analysis indicates that TgISPs have functionally repurposed the phospholipid-binding site likely to coordinate protein partners. Recruitment of larger protein complexes may also be aided through avidity-enhanced interactions resulting from multimerization of the ISPs. Overall, we propose a model where TgISPs recruit protein partners to the IMC to ensure correct progression of daughter cell formation.     

A leucine-rich repeat assembly approach for homology modeling of the human TLR5-10 and mouse TLR11-13 ectodomains

So far, 13 groups of mammalian Toll-like receptors (TLRs) have been identified. Most TLRs have been shown to recognize pathogen-associated molecular patterns from a wide range of invading agents and initiate both innate and adaptive immune responses. The TLR ectodomains are composed of varying numbers and types of leucine-rich repeats (LRRs). As the crystal structures are currently missing for most TLR ligand-binding ectodomains, homology modeling enables first predictions of their three-dimensional structures on the basis of the determined crystal structures of TLR ectodomains. However, the quality of the predicted models that are generated from full-length templates can be limited due to low sequence identity between the target and templates. To obtain better templates for modeling, we have developed an LRR template assembly approach. Individual LRR templates that are locally optimal for the target sequence are assembled into multiple templates. This method was validated through the comparison of a predicted model with the crystal structure of mouse TLR3. With this method, we also constructed ectodomain models of human TLR5, TLR6, TLR7, TLR8, TLR9, and TLR10 and mouse TLR11, TLR12, and TLR13 that can be used as first passes for a computational simulation of ligand docking or to design mutation experiments. This template assembly approach can be extended to other repetitive proteins.     

High-quality homology models derived from NMR and X-ray structures of E. coli proteins YgdK and Suf E suggest that all members of the YgdK/Suf E protein family are enhancers of cysteine desulfurases

The structural biology of proteins mediating iron-sulfur (Fe-S) cluster assembly is central for understanding several important biological processes. Here we present the NMR structure of the 16-kDa protein YgdK from Escherichia coli, which shares 35% sequence identity with the E. coli protein SufE. The SufE X-ray crystal structure was solved in parallel with the YdgK NMR structure in the Northeast Structural Genomics (NESG) consortium. Both proteins are (1) key components for Fe-S metabolism, (2) exhibit the same distinct fold, and (3) belong to a family of at least 70 prokaryotic and eukaryotic sequence homologs. Accurate homology models were calculated for the YgdK/SufE family based on YgdK NMR and SufE crystal structure. Both structural templates contributed equally, exemplifying synergy of NMR and X-ray crystallography. SufE acts as an enhancer of the cysteine desulfurase activity of SufS by SufE-SufS complex formation. A homology model of CsdA, a desulfurase encoded in the same operon as YgdK, was modeled using the X-ray structure of SufS as a template. Protein surface and electrostatic complementarities strongly suggest that YgdK and CsdA likewise form a functional two-component desulfurase complex. Moreover, structural features of YgdK and SufS, which can be linked to their interaction with desulfurases, are conserved in all homology models. It thus appears very likely that all members of the YgdK/SufE family act as enhancers of Suf-S-like desulfurases. The present study exemplifies that "refined" selection of two (or more) targets enables high-quality homology modeling of large protein families.