Solving Maximum Clique Problem for Protein Structure Similarity

Computing the similarity between two protein structures is a crucial task in molecular biology, and has been extensively investigated. Many protein structure comparison methods can be modeled as maximum weighted clique problems in specific k-partite graphs, referred here as alignment graphs. In this paper we present both a new integer programming formulation for solving such clique problems and a dedicated branch and bound algorithm for solving the maximum cardinality clique problem. Both approaches have been integrated in VAST, a software for aligning protein 3D structures largely used in the National Center for Biotechnology Information, an original clique solver which uses the well known Bron and Kerbosch algorithm (BK). Our computational results on real protein alignment instances show that our branch and bound algorithm is up to 116 times faster than BK.

The accuracy of protein structure alignment servers

Background: Protein structural alignment is one of the most fundamental and crucial areas of research in the domain of computational structural biology. Comparison of a protein structure with known structures helps to classify it as a new or belonging to a known group of proteins. This, in turn, is useful to determine the function of protein, its evolutionary relationship with other protein molecules and grasping principles underlying protein architecture and folding. Results: A large number of protein structure alignment methods are available. Each protein structure alignment tool has its own strengths andweaknesses that need to be highlighted.We compared and presented results of six most popular and publically available servers for protein structure comparison. These web-based servers were compared with the respect to functionality (features provided by these servers) and accuracy (how well the structural comparison is performed). The CATH was used as a reference. The results showed that overall CE was top performer. DALI and PhyreStorm showed similar results whereas PDBeFold showed the lowest performance. In case of few secondary structural elements, CE, DALI and PhyreStorm gave 100% success rate. Conclusion: Overall none of the structural alignment servers showed 100% success rate. Studies of overall performance, effect of mainly alpha and effect of mainly beta showed consistent performance. CE, DALI, FatCat and PhyreStorm showed more than 90% success rate.

Preferential assembly of epithelial sodium channel (ENaC) subunits in Xenopus oocytes: role of furin-mediated endogenous proteolysis.

The epithelial sodium channel (ENaC) is preferentially assembled into heteromeric alphabetagamma complexes. The alpha and gamma (not beta) subunits undergo proteolytic cleavage by endogenous furin-like activity correlating with increased ENaC function. We identified full-length subunits and their fragments at the cell surface, as well as in the intracellular pool, for all homo- and heteromeric combinations (alpha, beta, gamma, alphabeta, alphagamma, betagamma, and alphabetagamma). We assayed corresponding channel function as amiloride-sensitive sodium transport (I(Na)). We varied furin-mediated proteolysis by mutating the P1 site in alpha and/or gamma subunit furin consensus cleavage sites (alpha(mut) and gamma(mut)). Our findings were as follows. (i) The beta subunit alone is not transported to the cell surface nor cleaved upon assembly with the alpha and/or gamma subunits. (ii) The alpha subunit alone (or in combination with beta and/or gamma) is efficiently transported to the cell surface; a surface-expressed 65-kDa alpha ENaC fragment is undetected in alpha(mut)betagamma, and I(Na) is decreased by 60%. (iii) The gamma subunit alone does not appear at the cell surface; gamma co-expressed with alpha reaches the surface but is not detectably cleaved; and gamma in alphabetagamma complexes appears mainly as a 76-kDa species in the surface pool. Although basal I(Na) of alphabetagamma(mut) was similar to alphabetagamma, gamma(mut) was not detectably cleaved at the cell surface. Thus, furin-mediated cleavage is not essential for participation of alpha and gamma in alphabetagamma heteromers. Basal I(Na) is reduced by preventing furin-mediated cleavage of the alpha, but not gamma, subunits. Residual current in the absence of furin-mediated proteolysis may be due to non-furin endogenous proteases.

Oligomerization of the alpha 1a- and alpha 1b-adrenergic receptor subtypes. Potential implications in receptor internalization.

We combined biophysical, biochemical, and pharmacological approaches to investigate the ability of the alpha 1a- and alpha 1b-adrenergic receptor (AR) subtypes to form homo- and hetero-oligomers. Receptors tagged with different epitopes (hemagglutinin and Myc) or fluorescent proteins (cyan and green fluorescent proteins) were transiently expressed in HEK-293 cells either individually or in different combinations. Fluorescence resonance energy transfer measurements provided evidence that both the alpha 1a- and alpha 1b-AR can form homo-oligomers with similar transfer efficiency of approximately 0.10. Hetero-oligomers could also be observed between the alpha 1b- and the alpha 1a-AR subtypes but not between the alpha 1b-AR and the beta2-AR, the NK1 tachykinin, or the CCR5 chemokine receptors. Oligomerization of the alpha 1b-AR did not require the integrity of its C-tail, of two glycophorin motifs, or of the N-linked glycosylation sites at its N terminus. In contrast, helix I and, to a lesser extent, helix VII were found to play a role in the alpha 1b-AR homo-oligomerization. Receptor oligomerization was not influenced by the agonist epinephrine or by the inverse agonist prazosin. A constitutively active (A293E) as well as a signaling-deficient (R143E) mutant displayed oligomerization features similar to those of the wild type alpha 1b-AR. Confocal imaging revealed that oligomerization of the alpha1-AR subtypes correlated with their ability to co-internalize upon exposure to the agonist. The alpha 1a-selective agonist oxymetazoline induced the co-internalization of the alpha 1a- and alpha 1b-AR, whereas the alpha 1b-AR could not co-internalize with the NK1 tachykinin or CCR5 chemokine receptors. Oligomerization might therefore represent an additional mechanism regulating the physiological responses mediated by the alpha 1a- and alpha 1b-AR subtypes.

Syntaxin1a variants lacking an N-peptide or bearing the LE mutation bind to Munc18a in a closed conformation.

In neurons, soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) proteins drive the fusion of synaptic vesicles to the plasma membrane through the formation of a four-helix SNARE complex. Members of the Sec1/Munc18 protein family regulate membrane fusion through interactions with the syntaxin family of SNARE proteins. The neuronal protein Munc18a interacts with a closed conformation of the SNARE protein syntaxin1a (Syx1a) and with an assembled SNARE complex containing Syx1a in an open conformation. The N-peptide of Syx1a (amino acids 1-24) has been implicated in the transition of Munc18a-bound Syx1a to Munc18a-bound SNARE complex, but the underlying mechanism is not understood. Here we report the X-ray crystal structures of Munc18a bound to Syx1a with and without its native N-peptide (Syx1aΔN), along with small-angle X-ray scattering (SAXS) data for Munc18a bound to Syx1a, Syx1aΔN, and Syx1a L165A/E166A (LE), a mutation thought to render Syx1a in a constitutively open conformation. We show that all three complexes adopt the same global structure, in which Munc18a binds a closed conformation of Syx1a. We also identify a possible structural connection between the Syx1a N-peptide and SNARE domain that might be important for the transition of closed-to-open Syx1a in SNARE complex assembly. Although the role of the N-peptide in Munc18a-mediated SNARE complex assembly remains unclear, our results demonstrate that the N-peptide and LE mutation have no effect on the global conformation of the Munc18a-Syx1a complex.

An external site controls closing of the epithelial Na+ channel ENaC.

Members of the ENaC/degenerin family of ion channels include the epithelial sodium channel (ENaC), acid-sensing ion channels (ASICs) and the nematode Caenorhabditis elegans degenerins. These channels are activated by a variety of stimuli such as ligands (ASICs) and mechanical forces (degenerins), or otherwise are constitutively active (ENaC). Despite their functional heterogeneity, these channels might share common basic mechanisms for gating. Mutations of a conserved residue in the extracellular loop, namely the 'degenerin site' activate all members of the ENaC/degenerin family. Chemical modification of a cysteine introduced in the degenerin site of rat ENaC (betaS518C) by the sulfhydryl reagents MTSET or MTSEA, results in a approximately 3-fold increase in the open probability. This effect is due to an 8-fold shortening of channel closed times and an increase in the number of long openings. In contrast to the intracellular gating domain in the N-terminus which is critical for channel opening, the intact extracellular degenerin site is necessary for normal channel closing, as illustrated by our observation that modification of betaS518C destabilises the channel closed state. The modification by the sulfhydryl reagents is state- and size-dependent consistent with a conformational change of the degenerin site during channel opening and closing. We propose that the intracellular and extracellular modulatory sites act on a common channel gate and control the activity of ENaC at the cell surface.

C3-symmetric peptide scaffolds are functional mimetics of trimeric CD40L.

Interaction between CD40, a member of the tumor necrosis factor receptor (TNFR) superfamily, and its ligand CD40L, a 39-kDa glycoprotein, is essential for the development of humoral and cellular immune responses. Selective blockade or activation of this pathway provides the ground for the development of new treatments against immunologically based diseases and malignancies. Like other members of the TNF superfamily, CD40L monomers self-assemble around a threefold symmetry axis to form noncovalent homotrimers that can each bind three receptor molecules. Here, we report on the structure-based design of small synthetic molecules with C3 symmetry that can mimic CD40L homotrimers. These molecules interact with CD40, compete with the binding of CD40L to CD40, and reproduce, to a certain extent, the functional properties of the much larger homotrimeric soluble CD40L. Architectures based on rigid C3-symmetric cores may thus represent a general approach to mimicking homotrimers of the TNF superfamily.

In vivo activation of PPAR target genes by RXR homodimers.

The ability of a retinoid X receptor (RXR) to heterodimerize with many nuclear receptors, including LXR, PPAR, NGF1B and RAR, underscores its pivotal role within the nuclear receptor superfamily. Among these heterodimers, PPAR:RXR is considered an important signalling mediator of both PPAR ligands, such as fatty acids, and 9-cis retinoic acid (9-cis RA), an RXR ligand. In contrast, the existence of an RXR/9-cis RA signalling pathway independent of PPAR or any other dimerization partner remains disputed. Using in vivo chromatin immunoprecipitation, we now show that RXR homodimers can selectively bind to functional PPREs and induce transactivation. At the molecular level, this pathway requires stabilization of the homodimer-DNA complexes through ligand-dependent interaction with the coactivator SRC1 or TIF2. This pathway operates both in the absence and in the presence of PPAR, as assessed in cells carrying inactivating mutations in PPAR genes and in wild-type cells. In addition, this signalling pathway via PPREs is fully functional and can rescue the severe hypothermia phenotype observed in fasted PPARalpha-/- mice. These observations have important pharmacological implications for the development of new rexinoid-based treatments.

Biological activity of ectodysplasin A is conditioned by its collagen and heparan sulfate proteoglycan-binding domains.

Mutations in the TNF family ligand EDA1 cause X-linked hypohidrotic ectodermal dysplasia (XLHED), a condition characterized by defective development of skin appendages. The EDA1 protein displays a proteolytic processing site responsible for its conversion to a soluble form, a collagen domain, and a trimeric TNF homology domain (THD) that binds the receptor EDAR. In-frame deletions in the collagen domain reduced the thermal stability of EDA1. Removal of the collagen domain decreased its activity about 100-fold, as measured with natural and engineered EDA1-responsive cell lines. The collagen domain could be functionally replaced by multimerization domains or by cross-linking antibodies, suggesting that it functions as an oligomerization unit. Surprisingly, mature soluble EDA1 containing the collagen domain was poorly active when administered in newborn, EDA-deficient (Tabby) mice. This was due to a short stretch of basic amino acids located at the N terminus of the collagen domain that confers EDA1 with proteoglycan binding ability. In contrast to wild-type EDA1, EDA1 with mutations in this basic sequence was a potent inducer of tail hair development in vivo. Thus, the collagen domain activates EDA1 by multimerization, whereas the proteoglycan-binding domain may restrict the distribution of endogeneous EDA1 in vivo.

Formation of virus-like clusters is an intrinsic property of the tumor necrosis factor family member BAFF (B cell activating factor).

The oligomeric state of BAFF (B cell activing factor), a tumor necrosis factor (TNF) family cytokine that plays a critical role in B cell development and survival, has been the subject of recent debate. Myc-tagged BAFF starting at residue Gln136 was previously reported to crystallize as trimers at pH 4.5, whereas a histidine-tagged construct of BAFF, starting at residue Ala134, formed a virus-like cluster containing 60 monomers when crystallized at pH 9.0. The formation of the BAFF 60-mer was pH dependent, requiring pH >or= 7.0. More recently, 60-mer formation was suggested to be artificially induced by the histidine tag, and it was proposed that BAFF, like all other TNF family members, is trimeric. We report here that a construct of BAFF with no amino-terminal tag (Ala134-BAFF) can form a 60-mer in solution. Using size exclusion chromatography and static light scattering to monitor trimer to 60-mer ratios in BAFF preparations, we find that 60-mer formation is pH-dependent and requires histidine 218 within the DE loop of BAFF. Biacore measurements established that the affinity of Ala134-BAFF for the BAFF receptor BAFFR/BR3 is similar to that of myc-Gln136-BAFF, which is exclusively trimeric in solution. However, Ala134-BAFF is more efficacious than myc-Gln136-BAFF in inducing B cell proliferation in vitro. We additionally show that BAFF that is processed and secreted by 293T cells transfected with full-length BAFF, or by a histiocytic lymphoma cell line (U937) that expresses BAFF endogenously, forms a pH-dependent 60-mer in solution. Our results indicate that the formation of the 60-mer in solution by the BAFF extracellular domain is an intrinsic property of the protein, and therefore that this more active form of BAFF may be physiologically relevant.

Paired activation of two components within muscarinic M3 receptor dimers is required for recruitment of beta-arrestin-1 to the plasma membrane.

beta-Arrestins regulate the functioning of G protein-coupled receptors in a variety of cellular processes including receptor-mediated endocytosis and activation of signaling molecules such as ERK. A key event in these processes is the G protein-coupled receptor-mediated recruitment of beta-arrestins to the plasma membrane. However, despite extensive knowledge in this field, it is still disputable whether activation of signaling pathways via beta-arrestin recruitment entails paired activation of receptor dimers. To address this question, we investigated the ability of different muscarinic receptor dimers to recruit beta-arrestin-1 using both co-immunoprecipitation and fluorescence microscopy in COS-7 cells. Experimentally, we first made use of a mutated muscarinic M(3) receptor, which is deleted in most of the third intracellular loop (M(3)-short). Although still capable of activating phospholipase C, this receptor loses almost completely the ability to recruit beta-arrestin-1 following carbachol stimulation in COS-7 cells. Subsequently, M(3)-short was co-expressed with the M(3) receptor. Under these conditions, the M(3)/M(3)-short heterodimer could not recruit beta-arrestin-1 to the plasma membrane, even though the control M(3)/M(3) homodimer could. We next tested the ability of chimeric adrenergic muscarinic alpha(2)/M(3) and M(3)/alpha(2) heterodimeric receptors to co-immunoprecipitate with beta-arrestin-1 following stimulation with adrenergic and muscarinic agonists. beta-Arrestin-1 co-immunoprecipitation could be induced only when carbachol or clonidine were given together and not when the two agonists were supplied separately. Finally, we tested the reciprocal influence that each receptor may exert on the M(2)/M(3) heterodimer to recruit beta-arrestin-1. Remarkably, we observed that M(2)/M(3) heterodimers recruit significantly greater amounts of beta-arrestin-1 than their respective M(3)/M(3) or M(2)/M(2) homodimers. Altogether, these findings provide strong evidence in favor of the view that binding of beta-arrestin-1 to muscarinic M(3) receptors requires paired stimulation of two receptor components within the same receptor dimer.

Thermal denaturation and aggregation of hemoglobin of Glossoscolex paulistus in acid and neutral media

The thermal denaturation and aggregation of the HbGp, in the oxy- and cyanomet-forms, was investigated by DSC, AUC, DLS, optical absorption and CD, in the pH range from 5.0 to 7.0. Oxy-HbGp has a denaturation process partially reversible and dependent on the temperature. DSC melting curve is characterized by a single peak with Tc value of 333.4±0.2K for oxy-HbGp, while two peaks with Tc values of 332.2±0.1 and 338.4±0.2K are observed for cyanomet-HbGp, at pH 7.0. In acidic pH oxy- and cyanomet-HbGp are more stable showing higher Tc values and aggregation. AUC data show that, HbGp, at pH 7.0, upon denaturation, remains undissociated at 323K, presenting oligomeric dissociation at 333 (12±3% of tetramer and 88±5% of whole HbGp) and 343K (70±5% of monomer and 30±2% of trimer). DLS data show that the lag period before aggregation is dependent on the temperature and HbGp concentration. Optical absorption and CD results show that the increase of temperature leads to the oxy-HbGp oxidation and aggregation, above 331K, in acidic pH. CD data, for HbGp, present a greater thermal stability in acid medium than at neutral pH, with similar Tc values for both oxidation forms. Our data are consistent with previous studies and represents an advance in understanding the thermal stability of oligomeric HbGp structure. © 2012 Elsevier B.V.

Activation of DegP chaperone-protease via formation of large cage-like oligomers upon binding to substrate proteins.

Cells use molecular chaperones and proteases to implement the essential quality control mechanism of proteins. The DegP (HtrA) protein, essential for the survival of Escherichia coli cells at elevated temperatures with homologues found in almost all organisms uniquely has both functions. Here we report a mechanism for DegP to activate both functions via formation of large cage-like 12- and 24-mers after binding to substrate proteins. Cryo-electron microscopic and biochemical studies revealed that both oligomers are consistently assembled by blocks of DegP trimers, via pairwise PDZ1-PDZ2 interactions between neighboring trimers. Such interactions simultaneously eliminate the inhibitory effects of the PDZ2 domain. Additionally, both DegP oligomers were also observed in extracts of E. coli cells, strongly implicating their physiological importance.

Rubisco in complex with Rubisco large subunit methyltransferase.

SET domain protein lysine methyltransferases (PKMT) are a structurally unique class of enzymes that catalyze the specific methylation of lysine residues in a number of different substrates. Especially histone-specific SET domain PKMTs have received widespread attention because of their roles in the regulation of epigenetic gene expression and the development of some cancers. Rubisco large subunit methyltransferase (RLSMT) is a chloroplast-localized SET domain PKMT responsible for the formation of trimethyl-lysine-14 in the large subunit of Rubisco, an essential photosynthetic enzyme. Here, we have used cryoelectron microscopy to produce an 11-A density map of the Rubisco-RLSMT complex. The atomic model of the complex, obtained by fitting crystal structures of Rubisco and RLSMT into the density map, shows that the extensive contact regions between the 2 proteins are mainly mediated by hydrophobic residues and leucine-rich repeats. It further provides insights into potential conformational changes that may occur during substrate binding and catalysis. This study presents the first structural analysis of a SET domain PKMT in complex with its intact polypeptide substrate.

Subnanometer-resolution structures of the grass carp reovirus core and virion.

Grass carp reovirus (GCRV) is a member of the Aquareovirus genus of the family Reoviridae, a large family of double-stranded RNA (dsRNA) viruses infecting plants, insects, fishes and mammals. We report the first subnanometer-resolution three-dimensional structures of both GCRV core and virion by cryoelectron microscopy. These structures have allowed the delineation of interactions among the over 1000 molecules in this enormous macromolecular machine and a detailed comparison with other dsRNA viruses at the secondary-structure level. The GCRV core structure shows that the inner proteins have strong structural similarities with those of orthoreoviruses even at the level of secondary-structure elements, indicating that the structures involved in viral dsRNA interaction and transcription are highly conserved. In contrast, the level of similarity in structures decreases in the proteins situated in the outer layers of the virion. The proteins involved in host recognition and attachment exhibit the least similarities to other members of Reoviridae. Furthermore, in GCRV, the RNA-translocating turrets are in an open state and lack a counterpart for the sigma1 protein situated on top of the close turrets observed in mammalian orthoreovirus. Interestingly, the distribution and the organization of GCRV core proteins resemble those of the cytoplasmic polyhedrosis virus, a cypovirus and the structurally simplest member of the Reoviridae family. Our results suggest that GCRV occupies a unique structure niche between the simpler cypoviruses and the considerably more complex mammalian orthoreovirus, thus providing an important model for understanding the structural and functional conservation and diversity of this enormous family of dsRNA viruses.