Quaternary structure constraints on evolutionary sequence divergence

The structurally constrained protein evolution (SCPE) model simulates protein divergence considering protein structure explicitly. The model is based on the observation that protein structure is more conserved during evolution than the sequences encoding for that structure. In the previous work, the SCPE model considered only the tertiary structure. Here we show that the performance of the model is enhanced when the oligomeric structure is taken into account. Our results agree with recent evolutionary studies of oligomeric proteins, which show that conservation of the quaternary structure imposes additional constraints on sequence divergence. The incorporation of protein-protein interactions into protein evolution models may be important in the study of quaternary protein structures and complex protein assemblies.     

Mechanism of chaperone function in small heat shock proteins: dissociation of the HSP27 oligomer is required for recognition and binding of destabilized T4 lysozyme

Mammalian small heat shock proteins (sHSP) form polydisperse and dynamic oligomers that undergo equilibrium subunit exchange. Current models of their chaperone activity hypothesize that recognition and binding of protein non-native states involve changes in the oligomeric state. The equivalent thermodynamic representation is a set of three coupled equilibria that includes the sHSP oligomeric equilibrium, the substrate folding equilibrium, and the equilibrium binding between the sHSP and the substrate non-native states. To test this hypothesis and define the binding-competent oligomeric state of human Hsp27, we have perturbed the two former equilibria and quantitatively determined the consequences on binding. The substrate is a set of T4 lysozyme (T4L) mutants that bind under conditions that favor the folded state over the unfolded state by 10(2)-10(4)-fold. The concentration-dependent oligomer equilibrium of Hsp27 was perturbed by mutations that alter the relative stability of two major oligomeric states including phosphorylation-mimicking mutations that result in the dissociation to a small multimer over a wide range of concentrations. Correlation of binding isotherms with size exclusion chromatography analysis of the Hsp27 oligomer equilibrium demonstrates that the multimer is the binding-competent state. Binding occurs through two modes, each characterized by different affinity and number of binding sites, and results in T4L.Hsp27 complexes of different hydrodynamic properties. Mutants of the Hsp27 phosphorylation mimic that reverse the reduction in oligomer size also reduce the extent of T4L binding. Taken together, these results suggest a central role for the oligomeric equilibrium in regulating the chaperone activity of sHSP. The mutants identify sequence features important for modulating this equilibrium.     

Role of the conserved SRLFDQFFG region of alpha-crystallin, a small heat shock protein. Effect on oligomeric size, subunit exchange, and chaperone-like activity

Small heat shock proteins (sHsps) are necessary for several cellular functions and in stress tolerance. Most sHsps are oligomers; intersubunit interactions leading to changes in oligomeric structure and exposure of specific regions may modulate their functioning. Many sHsps, including alpha A- and alpha B-crystallin, contain a well conserved SRLFDQFFG sequence motif in the N-terminal region. Sequence-based prediction shows that it exhibits helical propensity with amphipathic character, suggesting that it plays a critical role in the structure and function of alpha-crystallins. In order to investigate the role of this motif in the structure and function of sHsps, we have made constructs deleting this sequence from alpha A- and alpha B-crystallin, overexpressed, purified, and studied these engineered proteins. Circular dichroism spectroscopic studies show changes in tertiary and secondary structure on deletion of the sequence. Glycerol density gradient centrifugation and dynamic light scattering studies show that the multimeric size of the mutant proteins is significantly reduced, indicating a role for this motif in higher order organization of the subunits. Both deletion mutants exhibit similar oligomeric size and increased chaperone-like activity. Urea-induced denaturation study shows that the SRLFDQFFG sequence contributes significantly to the structural stability. Fluorescence resonance energy transfer studies show that the rate of exchange of the subunits in the alpha Adel-crystallin oligomer is higher compared with that in the alpha A-crystallin oligomer, suggesting that this region contributes to the oligomer dynamics in addition to the higher order assembly and structural stability. Thus, our study shows that the SRLFDQFFG sequence is one of the critical motifs in structure-function regulation of alpha A- and alpha B-crystallin.     

Dissociative mechanism for irreversible thermal denaturation of oligomeric proteins

Protein stability is a fundamental characteristic essential for understanding conformational transformations of the proteins in the cell. When using protein preparations in biotechnology and biomedicine, the problem of protein stability is of great importance. The kinetics of denaturation of oligomeric proteins may have characteristic properties determined by the quaternary structure. The kinetic schemes of denaturation can include the multiple stages of conformational transitions in the protein oligomer and stages of reversible dissociation of the oligomer. In this case, the shape of the kinetic curve of denaturation or the shape of the melting curve registered by differential scanning calorimetry can vary with varying the protein concentration. The experimental data illustrating dissociative mechanism for irreversible thermal denaturation of oligomeric proteins have been summarized in the present review. The use of test systems based on thermal aggregation of oligomeric proteins for screening of agents possessing anti-aggregation activity is discussed.     

The theoretical analysis of kinetic behaviour of “hysteretic” allosteric enzymes V. Relaxation kinetics of dissociating enzyme systems

The expressions for relaxation time as a function of enzyme and specific ligand concentration are deduced for dissociating enzyme system 2p ? P (P is enzyme oligomer which is able to dissociate reversibly forming two identical halves p). It is assumed that ligand binding sites are equivalent and independent in each oligomeric enzyme form and the equilibrium between oligomeric forms develops rather slowly in comparison with the rate of the binding of the ligand. The kinetics of relaxation of the dissociating enzyme system 2p ? P with progressive change of the rate constants for association of oligomeric form p has been analysed in graphic form. The situations when one of the oligomeric enzyme forms is not able to bind the ligand are also considered. The principles of the analysis of relaxation kinetics of dissociating enzyme systems 2p ? P are discussed.     

Oligomerization triggers binding of a Bacillus thuringiensis Cry1Ab pore-forming toxin to aminopeptidase N receptor leading to insertion into membrane microdomains

Bacillus thuringiensis Cry1A toxins, in contrast to other pore-forming toxins, bind two putative receptor molecules, aminopeptidase N (APN) and cadherin-like proteins. Here we show that Cry1Ab toxin binding to these two receptors depends on the toxins' oligomeric structure. Toxin monomeric structure binds to Bt-R1, a cadherin-like protein, that induces proteolytic processing and oligomerization of the toxin (Gomez, I., Sanchez, J., Miranda, R., Bravo A., Soberon, M., FEBS Lett. (2002) 513, 242-246), while the oligomeric structure binds APN, which drives the toxin into the detergent-resistant membrane (DRM) microdomains causing pore formation. Cleavage of APN by phospholipase C prevented the location of Cry1Ab oligomer and Bt-R1 in the DRM microdomains and also attenuates toxin insertion into membranes despite the presence of Bt-R1. Immunoprecipitation experiments demonstrated that initial Cry1Ab toxin binding to Bt-R1 is followed by binding to APN. Also, immunoprecipitation of Cry1Ab toxin-binding proteins using pure oligomeric or monomeric structures showed that APN was more efficiently detected in samples immunoprecipitated with the oligomeric structure, while Bt-R1 was preferentially detected in samples immunoprecipitated with the monomeric Cry1Ab. These data agrees with the 200-fold higher apparent affinity of the oligomer than that of the monomer to an APN enriched protein extract. Our data suggest that the two receptors interact sequentially with different structural species of the toxin leading to its efficient membrane insertion.     

Protein-protein interactions in microsomal cytochrome P-450 isozyme LM2 and their effect on substrate binding

The effects of protein-protein interactions and substrate binding on the structure of the active site of rabbit liver microsomal cytochrome P-450 LM2 have been analyzed by resonance Raman spectroscopy of the monomeric and oligomeric protein in solution. Also H2O2-dependent catalytic activities of the two states have been compared. The two vinyl substituents of the heme exhibit different orientations, as indicated by the frequencies and intensities of their stretching vibrations. One group lies in the plane of the heme and remains unchanged in the two states of cytochrome P-450 LM2, the other is tilted out of the plane. The tilting angle in oligomers was smaller than in monomers. These vinyl stretching modes together with some porphyrin modes, were found to be sensitive indicators of the quaternary structure and of substrate binding. In both the oligomer and the monomer, substrate binding causes changes of the relative intensities of some porphyrin modes and the vinyl stretching vibrations which may reflect modifications of the electronic transitions due to hydrophobic interactions between the bound substrate and the heme. In contrast to the monomeric cytochrome P-450 LM2, benzphetamine binding to the oligomers of this isozyme additionally produces a shift of the spin-state equilibrium. This indicates that in the oligomer the substrate-binding pocket is converted by protein-protein interaction to a structure that forces substrates to interfere with the sixth ligands, inducing an increase of the five-coordinated high-spin configuration. In the monomer the substrate-binding pocket can accommodate benzphetamine without affecting the spin state. Binding of imidazole to the monomeric and oligomeric cytochrome P-450 LM2 produces essentially the same resonance Raman spectra. Apparently the replacement of the native sixth ligand by imidazole disturbs the structure of the active site in such a way that it becomes insensitive to protein-protein interactions. H2O2-dependent N-demethylation of benzphetamine and aniline p-hydroxylation by cytochrome P-450 LM2 did not depend on its state of aggregation.     

Compensatory Evolution of a WW Domain Variant Lacking the Strictly Conserved Trp Residue

Replacement of conserved amino acid residues during evolution of proteins can lead to divergence and the formation of new families with novel functions, but is often deleterious to both protein structure and function. Using the WW domain, we experimentally examined whether and to what degree second-site mutations can compensate for the reduction of function and loss of structure that accompany substitution of a strictly conserved amino acid residue. The W17F mutant of the WW domain, with substitution of the most strictly conserved Trp residue, is known to lack a specific three-dimensional structure and shows reduced binding affinity in comparison to the wild type. To obtain second-site revertants, we performed a selection experiment based on the proline-rich peptide (PY ligand) binding affinity using the W17F mutant as the initial sequence. After selection by ribosome display, we were able to select revertants that exhibited a maximum ninefold higher affinity to the PY ligand than the W17F mutant and showed an even better affinity than the wild type. In addition, we found that the functional restoration resulted in increased binding specificity in selected revertants, and the structures were more compact, with increased amounts of secondary structure, in comparison to the W17F mutant. Our results suggest that the defective structure and function of the proteins caused by mutations in highly conserved residues occurring through divergent evolution not only can be restored but can be further improved by compensatory mutations.     

Structural changes of the Cry1Ac oligomeric pre-pore from bacillus thuringiensis induced by N-acetylgalactosamine facilitates toxin membrane insertion

The primary action of Cry toxins produced by Bacillus thuringiensis is to lyse midgut epithelial cells in their target insect by forming lytic pores. The toxin-receptor interaction is a complex process, involving multiple interactions with different receptor and carbohydrate molecules. It has been proposed that Cry1A toxins sequentially interact with a cadherin receptor, leading to the formation of a pre-pore oligomer structure, and that the oligomeric structure binds to glycosylphosphatidyl-inositol-anchored aminopeptidase-N (APN) receptor. The Cry1Ac toxin specifically recognizes the N-acetylgalactosamine (GalNAc) carbohydrate present in the APN receptor from Manduca sexta larvae. In this work, we show that the Cry1Ac pre-pore oligomer has a higher binding affinity with APN than the monomeric toxin. The effects of GalNAc binding on the toxin structure were studied in the monomeric Cry1Ac, in the soluble pre-pore oligomeric structure, and in its membrane inserted state by recording the fluorescence status of the tryptophan (W) residues. Our results indicate that the W residues of Cry1Ac have a different exposure to the solvent when compared with that of the closely related Cry1Ab toxin. GalNAc binding specifically affects the exposure of W545 in the pre-pore oligomer in contrast to the monomer where GalNAc binding did not affect the fluorescence of the toxin. These results indicate a subtle conformational change in the GalNAc binding pocket in the pre-pore oligomer that could explain the increased binding affinity of the Cry1Ac pre-pore to APN. Although our analysis did not reveal major structural changes in the pore-forming domain I upon GalNAc binding, it showed that sugar interaction enhanced membrane insertion of soluble pre-pore oligomeric structure. Therefore, the data presented here permits to propose a model in which the interaction of Cry1Ac pre-pore oligomer with APN receptor facilitates membrane insertion and pore formation.     

Structure-based ligand binding sites of protein p14.5, a member of protein family YER057c/YIL051c/YjgF

Seventeen mutants with one, two or three amino acids substitutions of human protein p14.5, homologue to well-known tumor antigen from goat liver UK114 and a member of proteins YER057c/YIL051c/YjgF family, have been used for structure-functional relation studies and ligand binding analysis using cross-linking by triacryloyl-hexahydro-s-triazine (TAT), size exclusion chromatography, free fatty acid and 8-anilino-1-naphthalenesulfonic acid (ANS) binding assays. Amino acids having the most significant impact on the ligand binding activity have been determined: R107, N93, Y21 and F89. Arginine 107 has been identified as the most accessible amino acid in the cleft. Trimeric structure of protein p14.5 has been confirmed as being essential for stoichiometric small ligand binding activity and oligomeric structure of p14. Ligand binding activity may be related with the biological functions of these proteins, which still are not understood well.     

Collective dynamics differentiates functional divergence in protein evolution

Protein evolution is most commonly studied by analyzing related protein sequences and generating ancestral sequences through Bayesian and Maximum Likelihood methods, and/or by resurrecting ancestral proteins in the lab and performing ligand binding studies to determine function. Structural and dynamic evolution have largely been left out of molecular evolution studies. Here we incorporate both structure and dynamics to elucidate the molecular principles behind the divergence in the evolutionary path of the steroid receptor proteins. We determine the likely structure of three evolutionarily diverged ancestral steroid receptor proteins using the Zipping and Assembly Method with FRODA (ZAMF). Our predictions are within ∼2.7 Å all-atom RMSD of the respective crystal structures of the ancestral steroid receptors. Beyond static structure prediction, a particular feature of ZAMF is that it generates protein dynamics information. We investigate the differences in conformational dynamics of diverged proteins by obtaining the most collective motion through essential dynamics. Strikingly, our analysis shows that evolutionarily diverged proteins of the same family do not share the same dynamic subspace, while those sharing the same function are simultaneously clustered together and distant from those, that have functionally diverged. Dynamic analysis also enables those mutations that most affect dynamics to be identified. It correctly predicts all mutations (functional and permissive) necessary to evolve new function and ∼60% of permissive mutations necessary to recover ancestral function.     

Fold independent structural comparisons of protein-ligand binding sites for exploring functional relationships

The rapid growth in protein structural data and the emergence of structural genomics projects have increased the need for automatic structure analysis and tools for function prediction. Small molecule recognition is critical to the function of many proteins; therefore, determination of ligand binding site similarity is important for understanding ligand interactions and may allow their functional classification. Here, we present a binding sites database (SitesBase) that given a known protein-ligand binding site allows rapid retrieval of other binding sites with similar structure independent of overall sequence or fold similarity. However, each match is also annotated with sequence similarity and fold information to aid interpretation of structure and functional similarity. Similarity in ligand binding sites can indicate common binding modes and recognition of similar molecules, allowing potential inference of function for an uncharacterised protein or providing additional evidence of common function where sequence or fold similarity is already known. Alternatively, the resource can provide valuable information for detailed studies of molecular recognition including structure-based ligand design and in understanding ligand cross-reactivity. Here, we show examples of atomic similarity between superfamily or more distant fold relatives as well as between seemingly unrelated proteins. Assignment of unclassified proteins to structural superfamiles is also undertaken and in most cases substantiates assignments made using sequence similarity. Correct assignment is also possible where sequence similarity fails to find significant matches, illustrating the potential use of binding site comparisons for newly determined proteins.     

The biological functions of low-frequency vibrations (phonons) 5. A phenomenological theory

Low-frequency internal motions of a biomacromolecule are thought to possess significant biological function from the dynamic point of view. In this paper, a general phenomenological theory is established by which it is clearly verified that low-frequency resonance plays a central role in the energy transmission required during the cooperative interaction between subunits in a protein oligomer. According to the present theory, it is found that the energy transmission between a pair of diagonal subunits in a protein oligomer with a polygon arrangement is the most efficient, so as to in a sense further predict that after a ligand is bound to a subunit by random collision, its diagonal subunit in the same protein oligomer will possess the greatest probability of binding with the next ligand. Furthermore, based on the concept of the 'resonance-controlled trigger' derived from the phenomenological theory, it is feasible to estimate the lower time limit of allosteric transition from one subunit to the other. Such a time limit depends on the dominant low-frequency mode of each subunit, the ratio of the coupling force constant to the corresponding inherent force constant, as well as the geometrical arrangement of subunits in a protein oligomer. So far none of the allosteric transitions observed in proteins has exceeded the time limit as defined here, indicating a logical consistency between our theory and the experiments.     

Oligomerization triggers binding of a Bacillus thuringiensis Cry1Ab pore-forming toxin to aminopeptidase N receptor leading to insertion into membrane microdomains

Bacillus thuringiensis Cry1A toxins, in contrast to other pore-forming toxins, bind two putative receptor molecules, aminopeptidase N (APN) and cadherin-like proteins. Here we show that Cry1Ab toxin binding to these two receptors depends on the toxins' oligomeric structure. Toxin monomeric structure binds to Bt-R₁, a cadherin-like protein, that induces proteolytic processing and oligomerization of the toxin (Gómez, I., Sánchez, J., Miranda, R., Bravo A., Soberón, M., FEBS Lett. (2002) 513, 242-246), while the oligomeric structure binds APN, which drives the toxin into the detergent-resistant membrane (DRM) microdomains causing pore formation. Cleavage of APN by phospholipase C prevented the location of Cry1Ab oligomer and Bt-R₁ in the DRM microdomains and also attenuates toxin insertion into membranes despite the presence of Bt-R₁. Immunoprecipitation experiments demonstrated that initial Cry1Ab toxin binding to Bt-R₁ is followed by binding to APN. Also, immunoprecipitation of Cry1Ab toxin-binding proteins using pure oligomeric or monomeric structures showed that APN was more efficiently detected in samples immunoprecipitated with the oligomeric structure, while Bt-R₁ was preferentially detected in samples immunoprecipitated with the monomeric Cry1Ab. These data agrees with the 200-fold higher apparent affinity of the oligomer than that of the monomer to an APN enriched protein extract. Our data suggest that the two receptors interact sequentially with different structural species of the toxin leading to its efficient membrane insertion.     

Searching distant homologs of the regulatory ACT domain in phenylalanine hydroxylase

High sequence divergence, evolutionary mobility, and superfold topology characterize the ACT domain. Frequently found in multidomain proteins, these domains induce allosteric effects by binding a regulatory ligand usually to an ACT domain dimer interface. In mammalian phenylalanine hydroxylase (PAH), no contacts are formed between ACT domains, and the domain promotes an allosteric effect despite the apparent lack of ligand binding. The increased functional scenario of this abundant domain encouraged us to search for distant homologs, aiming to enhance the understanding of the ACT domain in general and the ACT domain of PAH in particular. The PDB was searched using the FATCAT server with the ACT domain of PAH as a query. The hits that were confirmed by the SSAP algorithm were divided into known ACT domains (KADs) and potential ACT domains (PADs). The FATCAT/SSAP procedure recognized most of the established KADs, as well 18 so far unrecognized non-redundant PADs with extremely low sequence identities and high divergence in functionality and oligomerization. However, analysis of the structural similarity provides remarkable clustering of the proteins according to similarities in ligand binding. Despite enormous sequence divergence and high functional variability, there is a common regulatory theme among these domains. The results reveal the close relationships of the ACT domain of PAH with amino acid binding and metallobinding ACT domains and with acylphosphatase.     

The highly dynamic oligomeric structure of bradavidin II is unique among avidin proteins

Bradavidin II is a biotin‐binding protein from Bradyrhizobium japonicum that resembles chicken avidin and bacterial streptavidin. A biophysical characterization was carried out using dynamic light scattering, native mass spectrometry, differential scanning calorimetry, and isothermal titration calorimetry combined with structural characterization using X‐ray crystallography. These observations revealed that bradavidin II differs from canonical homotetrameric avidin protein family members in its quaternary structure. In contrast with the other avidins, bradavidin II appears to have a dynamic (transient) oligomeric state in solution. It is monomeric at low protein concentrations but forms higher oligomeric assemblies at higher concentrations. The crystal structure of bradavidin II revealed an important role for Phe42 in shielding the bound ligand from surrounding water molecules, thus functionally replacing the L7,8 loop essential for tight ligand binding in avidin and streptavidin. This bradavidin II characterization opens new avenues for oligomerization‐independent biotin‐binding protein development.     

Genome analysis of the unicellular green alga Chlamydomonas reinhardtii Indicates an ancient evolutionary origin for key pattern recognition and cell-signaling protein families

The evolution of specific cell signaling and adhesion domains may have played an important role in the transition to a multicellular existence in the metazoans. Genomic analysis indicates that several signaling domains predominately found in animals are also present in the unicellular green alga, Chlamydomonas reinhardtii. A large group of proteins is present, containing scavenger receptor cysteine-rich (SRCR) and C-type lectin domains, which function in ligand binding and play key roles in the innate immune system of animals. Chlamydomonas also contains a large family of putative tyrosine kinases, suggesting an important role for phosphotyrosine signaling in the green algae. These important signaling domains may therefore be widespread among eukaryotes and most probably evolved in ancestral eukaryotes before the divergence of the Opisthokonts (the animal and fungal lineage).     

Sequence–function relationships in folding upon binding

Folding coupled to binding is ubiquitous in biology. Nevertheless, the relationship of sequence to function for protein segments that undergo coupled binding and folding remains to be determined. Specifically, it is not known if the well-established rules that govern protein folding and stability are relevant to ligand-linked folding transitions. Upon small ligand biotinoyl-5′-AMP (bio-5′-AMP) binding the Escherichia coli protein BirA undergoes a disorder-to-order transition that results in formation of a network of packed hydrophobic side chains. Ligand binding is also allosterically coupled to protein association, with bio-5′-AMP binding enhancing the dimerization free energy by −4.0 kcal/mol. Previous studies indicated that single alanine replacements in a three residue hydrophobic cluster that contributes to the larger network disrupt cluster formation, ligand binding, and allosteric activation of protein association. In this work, combined equilibrium and kinetic measurements of BirA variants with alanine substitutions in the entire hydrophobic network reveal large functional perturbations resulting from any single substitution and highly non-additive effects of multiple substitutions. These substitutions also disrupt ligand-linked folding. The combined results suggest that, analogous to protein folding, functional disorder-to-order linked to binding requires optimal packing of the relevant hydrophobic side chains that contribute to the transition. The potential for many combinations of residues to satisfy this requirement implies that, although functionally important, segments of homologous proteins that undergo folding linked to binding can exhibit sequence divergence.