Transitive homology-guided structural studies lead to discovery of Cro proteins with 40% sequence identity but different folds

Fluid protein fold space and its implications

Fold-switching proteins, which remodel their secondary and tertiary structures in response to cellular stimuli, suggest a new view of protein fold space. For decades, experimental evidence has indicated that protein fold space is discrete: dissimilar folds are encoded by dissimilar amino acid sequences. Challenging this assumption, fold-switching proteins interconnect discrete groups of dissimilar protein folds, making protein fold space fluid. Three recent observations support the concept of fluid fold space: (1) some amino acid sequences interconvert between folds with distinct secondary structures, (2) some naturally occurring sequences have switched folds by stepwise mutation, and (3) fold switching is evolutionarily selected and likely confers advantage. These observations indicate that minor amino acid sequence modifications can transform protein structure and function. Consequently, proteomic structural and functional diversity may be expanded by alternative splicing, small nucleotide polymorphisms, post-translational modifications, and modified translation rates.

Identification of a covert evolutionary pathway between two protein folds

Although homologous protein sequences are expected to adopt similar structures, some amino acid substitutions can interconvert α-helices and β-sheets. Such fold switching may have occurred over evolutionary history, but supporting evidence has been limited by the: (1) abundance and diversity of sequenced genes, (2) quantity of experimentally determined protein structures, and (3) assumptions underlying the statistical methods used to infer homology. Here, we overcome these barriers by applying multiple statistical methods to a family of ~600,000 bacterial response regulator proteins. We find that their homologous DNA-binding subunits assume divergent structures: helix-turn-helix versus α-helix + β-sheet (winged helix). Phylogenetic analyses, ancestral sequence reconstruction, and AlphaFold2 models indicate that amino acid substitutions facilitated a switch from helix-turn-helix into winged helix. This structural transformation likely expanded DNA-binding specificity. Our approach uncovers an evolutionary pathway between two protein folds and provides a methodology to identify secondary structure switching in other protein families.

A review of visualisations of protein fold networks and their relationship with sequence and function

Proteins form arguably the most significant link between genotype and phenotype. Understanding the relationship between protein sequence and structure, and applying this knowledge to predict function, is difficult. One way to investigate these relationships is by considering the space of protein folds and how one might move from fold to fold through similarity, or potential evolutionary relationships. The many individual characterisations of fold space presented in the literature can tell us a lot about how well the current Protein Data Bank represents protein fold space, how convergence and divergence may affect protein evolution, how proteins affect the whole of which they are part, and how proteins themselves function. A synthesis of these different approaches and viewpoints seems the most likely way to further our knowledge of protein structure evolution and thus, facilitate improved protein structure design and prediction.

A high-throughput predictive method for sequence-similar fold switchers

Although most experimentally characterized proteins with similar sequences assume the same folds and perform similar functions, an increasing number of exceptions is emerging. One class of exceptions comprises sequence-similar fold switchers, whose secondary structures shift from α-helix <-> β-sheet through a small number of mutations, a sequence insertion, or a deletion. Predictive methods for identifying sequence-similar fold switchers are desirable because some are associated with disease and/or can perform different functions in cells. Here, we use homology-based secondary structure predictions to identify sequence-similar fold switchers from their amino acid sequences alone. To do this, we predicted the secondary structures of sequence-similar fold switchers using three different homology-based secondary structure predictors: PSIPRED, JPred4, and SPIDER3. We found that α-helix <-> β-strand prediction discrepancies from JPred4 discriminated between the different conformations of sequence-similar fold switchers with high statistical significance (P < 1.8*10-19 ). Thus, we used these discrepancies as a classifier and found that they can often robustly discriminate between sequence-similar fold switchers and sequence-similar proteins that maintain the same folds (Matthews Correlation Coefficient of 0.82). We found that JPred4 is a more robust predictor of sequence-similar fold switchers because of (a) the curated sequence database it uses to produce multiple sequence alignments and (b) its use of sequence profiles based on Hidden Markov Models. Our results indicate that inconsistencies between JPred4 secondary structure predictions can be used to identify some sequence-similar fold switchers from their sequences alone. Thus, the negative information from inconsistent secondary structure predictions can potentially be leveraged to identify sequence-similar fold switchers from the broad base of genomic sequences.

Dynamic and conformational switching in proteins

Using bioinformatic methods for treating protein dynamics, developed in earlier work, we study the relationship between sequence mobility and dynamics in proteins. It is shown that sequence mobility drives a transition between two dynamic regimes in proteins, and that the specific details of this transition differ qualitatively between α-helical proteins and those in other structural classes. We examine the possibility that conformational switching is related to dynamic switching, by considering a specific system of sequences which exhibit the switching phenomenon. It is shown that a relationship between dynamic and conformational switching is entirely plausible.

Functional and Regulatory Roles of Fold-Switching Proteins

Fold-switching proteins respond to cellular stimuli by remodeling their secondary structures and changing their functions. Whereas several previous reviews have focused on various structural, physical-chemical, and evolutionary aspects of this newly emerging class of proteins, this minireview focuses on how fold switching modulates protein function and regulates biological processes. It first compares and contrasts fold switchers with other known types of proteins. Second, it presents examples of how various proteins can change their functions through fold switching. Third, it demonstrates that fold switchers can regulate biological processes by discussing two proteins, RfaH and KaiB, whose dramatic secondary structure remodeling events directly affect gene expression and a circadian clock, respectively. Finally, this minireview discusses how the field of protein fold switching might advance.

Conformational Dynamics from Ambiguous Zinc Coordination in the RanBP2-Type Zinc Finger of RBM5

The multi-domain RNA binding protein RBM5 is a molecular signature of metastasis. RBM5 regulates alternative splicing of apoptotic genes including the cell death receptor Fas and the initiator Caspase-2. The RBM5 RanBP2-type zinc finger (Zf1) is known to specifically recognize single-stranded RNAs with high affinity. Here, we study the structure and conformational dynamics of the Zf1 zinc finger of human RBM5 using NMR. We show that the presence of a non-canonical cysteine in Zf1 kinetically destabilizes the protein. Metal-exchange kinetics show that mutation of the cysteine establishes high-affinity coordination of the zinc. Our data indicate that selection of such a structurally destabilizing mutation during the course of evolution could present an opportunity for functional adaptation of the protein.

A multifunnel energy landscape encodes the competing α-helix and β-hairpin conformations for a designed peptide

Depending on the amino acid sequence, as well as the local environment, some peptides have the capability to fold into multiple secondary structures. Conformational switching between such structures is a key element of protein folding and aggregation. Specifically, understanding the molecular mechanism underlying the transition from an α-helix to a β-hairpin is critical because it is thought to be a harbinger of amyloid assembly. In this study, we explore the energy landscape for an 18-residue peptide (DP5), designed by Araki and Tamura to exhibit equal propensities for the α-helical and β-hairpin forms. We find that the degeneracy is encoded in the multifunnel nature of the underlying free energy landscape. In agreement with experiment, we also observe that mutation of tyrosine at position 12 to serine shifts the equilibrium in favor of the α-helix conformation, by altering the landscape topography. The transition from the α-helix to the β-hairpin is a complex stepwise process, and occurs via collapsed coil-like intermediates. Our findings suggest that even a single mutation can tune the emergent features of the landscape, providing an efficient route to protein design. Interestingly, the transition pathways for the conformational switch seem to be minimally perturbed upon mutation, suggesting that there could be universal microscopic features that are conserved among different switch-competent protein sequences.

Extreme divergence between one-to-one orthologs: the structure of N15 Cro bound to operator DNA and its relationship to the λ Cro complex

The gene cro promotes lytic growth of phages through binding of Cro protein dimers to regulatory DNA sites. Most Cro proteins are one-to-one orthologs, yet their sequence, structure and binding site sequences are quite divergent across lambdoid phages. We report the cocrystal structure of bacteriophage N15 Cro with a symmetric consensus site. We contrast this complex with an orthologous structure from phage λ, which has a dissimilar binding site sequence and a Cro protein that is highly divergent in sequence, dimerization interface and protein fold. The N15 Cro complex has less DNA bending and smaller DNA-induced changes in protein structure. N15 Cro makes fewer direct contacts and hydrogen bonds to bases, relying mostly on water-mediated and Van der Waals contacts to recognize the sequence. The recognition helices of N15 Cro and λ Cro make mostly nonhomologous and nonanalogous contacts. Interface alignment scores show that half-site binding geometries of N15 Cro and λ Cro are less similar to each other than to distantly related CI repressors. Despite this divergence, the Cro family shows several code-like protein–DNA sequence covariations. In some cases, orthologous genes can achieve a similar biological function using very different specific molecular interactions.

Multistep mutational transformation of a protein fold through structural intermediates

New protein folds may evolve from existing folds through metamorphic evolution involving a dramatic switch in structure. To mimic pathways by which amino acid sequence changes could induce a change in fold, we designed two folded hybrids of Xfaso 1 and Pfl 6, a pair of homologous Cro protein sequences with ~40% identity but different folds (all-α vs. α + β, respectively). Each hybrid, XPH1 or XPH2, is 85% identical in sequence to its parent, Xfaso 1 or Pfl 6, respectively; 55% identical to its noncognate parent; and ~70% identical to the other hybrid. XPH1 and XPH2 also feature a designed hybrid chameleon sequence corresponding to the C-terminal region, which switched from α-helical to β-sheet structure during Cro evolution. We report solution nuclear magnetic resonance (NMR) structures of XPH1 and XPH2 at 0.3 Å and 0.5 Å backbone root mean square deviation (RMSD), respectively. XPH1 retains a global fold generally similar to Xfaso 1, and XPH2 retains a fold similar to Pfl 6, as measured by TM-align scores (~0.7), DALI Z-scores (7-9), and backbone RMSD (2-3 Å RMSD for the most ordered regions). However, these scores also indicate significant deviations in structure. Most notably, XPH1 and XPH2 have different, and intermediate, secondary structure content relative to Xfaso 1 and Pfl 6. The multistep progression in sequence, from Xfaso 1 to XPH1 to XPH2 to Pfl 6, thus involves both abrupt and gradual changes in folding pattern. The plasticity of some protein folds may allow for "polymetamorphic" evolution through intermediate structures.

Structural metamorphism and polymorphism in proteins on the brink of thermodynamic stability

The classical view of the structure-function paradigm advanced by Anfinsen in the 1960s is that a protein's function is inextricably linked to its three-dimensional structure and is encrypted in its amino acid sequence. However, it is now known that a significant fraction of the proteome consists of intrinsically disordered proteins (IDPs). These proteins populate a polymorphic ensemble of conformations rather than a unique structure but are still capable of performing biological functions. At the boundary, between well-ordered and inherently disordered states are proteins that are on the brink of stability, either weakly stable ordered systems or disordered but on the verge of being stable. In such marginal states, even relatively minor changes can significantly alter the energy landscape, leading to large-scale conformational remodeling. Some proteins on the edge of stability are metamorphic, with the capacity to switch from one fold topology to another in response to an environmental trigger (e.g., pH, temperature/salt, redox). Many IDPs, on the other hand, are marginally unstable such that small perturbations (e.g., phosphorylation, ligands) tip the balance over to a range of ordered, partially ordered, or even more disordered states. In general, the structural transitions described by metamorphic fold switches and polymorphic IDPs possess a number of common features including low or diminished stability, large-scale conformational changes, critical disordered regions, latent or attenuated binding sites, and expansion of function. We suggest that these transitions are, therefore, conceptually and mechanistically analogous, representing adjacent regions in the continuum of order/disorder transitions.

Evolutionary convergence and divergence in archaeal chromosomal proteins and Chromo-like domains from bacteria and eukaryotes

SH3-fold-β-barrel domains of the chromo-like superfamily recognize epigenetic marks in eukaryotic proteins. Their provenance has been placed either in archaea, based on apparent structural similarity to chromatin-compacting Sul7d and Cren7 proteins, or in bacteria based on the presence of sequence homologs. Using sequence and structural evidence we establish that the archaeal Cren7/Sul7 proteins emerged from a zinc ribbon (ZnR) ancestor. Further, we show that the ancestral eukaryotic chromo-like domains evolved from bacterial versions, likely acquired from early endosymbioses, which already possessed an aromatic cage for recognition of modified amino-groups. These bacterial versions are part of a radiation of secreted SH3-fold domains, which spawned both chromo-like domains and classical SH3 domains in the context of peptide-recognition in the peptidoglycan or the extracellular matrix. This establishes that Cren7/Sul7 converged to a “SH3”-like state from a ZnR precursor via the loss of metal-chelation and acquisition of stronger hydrophobic interactions; it is unlikely to have participated in the evolution of the chromo-like domains. We show that archaea possess several Cren7/Sul7-related proteins with intact Zn-chelating ligands, which we predict to play previously unstudied roles in chromosome segregation during cell-division comparable to the PRC barrel and CdvA domain proteins.

Lessons from making the Structural Classification of Proteins (SCOP) and their implications for protein structure modelling

The Structural Classification of Proteins (SCOP) database has facilitated the development of many tools and algorithms and it has been successfully used in protein structure prediction and large-scale genome annotations. During the development of SCOP, numerous exceptions were found to topological rules, along with complex evolutionary scenarios and peculiarities in proteins including the ability to fold into alternative structures. This article reviews cases of structural variations observed for individual proteins and among groups of homologues, knowledge of which is essential for protein structure modelling.

The Role of Evolutionary Selection in the Dynamics of Protein Structure Evolution

Homology modeling is a powerful tool for predicting a protein's structure. This approach is successful because proteins whose sequences are only 30% identical still adopt the same structure, while structure similarity rapidly deteriorates beyond the 30% threshold. By studying the divergence of protein structure as sequence evolves in real proteins and in evolutionary simulations, we show that this nonlinear sequence-structure relationship emerges as a result of selection for protein folding stability in divergent evolution. Fitness constraints prevent the emergence of unstable protein evolutionary intermediates, thereby enforcing evolutionary paths that preserve protein structure despite broad sequence divergence. However, on longer timescales, evolution is punctuated by rare events where the fitness barriers obstructing structure evolution are overcome and discovery of new structures occurs. We outline biophysical and evolutionary rationale for broad variation in protein family sizes, prevalence of compact structures among ancient proteins, and more rapid structure evolution of proteins with lower packing density.

Development of a motif-based topology-independent structure comparison method to identify evolutionarily related folds

Structure conservation, functional similarities, and homologous relationships that exist across diverse protein topologies suggest that some regions of the protein fold universe are continuous. However, the current structure classification systems are based on hierarchical organizations, which cannot accommodate structural relationships that span fold definitions. Here, we describe a novel, super-secondary-structure motif-based, topology-independent structure comparison method (SmotifCOMP) that is able to quantitatively identify structural relationships between disparate topologies. The basis of SmotifCOMP is a systematically defined super-secondary-structure motif library whose representative geometries are shown to be saturated in the Protein Data Bank and exhibit a unique distribution within the known folds. SmotifCOMP offers a robust and quantitative technique to compare domains that adopt different topologies since the method does not rely on a global superposition. SmotifCOMP is used to perform an exhaustive comparison of the known folds and the identified relationships are used to produce a nonhierarchical representation of the fold space that reflects the notion of a continuous and connected fold universe. The current work offers insight into previously hypothesized evolutionary relationships between disparate folds and provides a resource for exploring novel ones. Proteins 2016; 84:1859-1874. © 2016 Wiley Periodicals, Inc.

Designing pH induced fold switch in proteins

This work investigates the computational design of a pH induced protein fold switch based on a self-consistent mean-field approach by identifying the ensemble averaged characteristics of sequences that encode a fold switch. The primary challenge to balance the alternative sets of interactions present in both target structures is overcome by simultaneously optimizing two foldability criteria corresponding to two target structures. The change in pH is modeled by altering the residual charge on the amino acids. The energy landscape of the fold switch protein is found to be double funneled. The fold switch sequences stabilize the interactions of the sites with similar relative surface accessibility in both target structures. Fold switch sequences have low sequence complexity and hence lower sequence entropy. The pH induced fold switch is mediated by attractive electrostatic interactions rather than hydrophobic-hydrophobic contacts. This study may provide valuable insights to the design of fold switch proteins.

A λ Cro-Like Repressor Is Essential for the Induction of Conjugative Transfer of SXT/R391 Elements in Response to DNA Damage

Integrative and conjugative elements (ICEs) of the SXT/R391 family are the main contributors to acquired multidrug resistance in the seventh pandemic lineage of Vibrio cholerae, the etiological agent of the diarrheal disease cholera. Conjugative transfer of SXT/R391 ICEs is triggered by antibiotics and agents promoting DNA damage through RecA-dependent autoproteolysis of SetR, an ICE-encoded λ CI-like repressor. Here, we describe the role of CroS, a distant λ Cro homolog, as a key component contributing to the regulation of expression of the activator SetCD that orchestrates the expression of the conjugative transfer genes. We show that deletion of croS abolishes the SOS response-dependent induction of SXT despite the presence of a functional setR gene. Using quantitative reverse transcription-PCR and lacZ reporter assays, we also show that CroS represses setR and setCD expression by binding to operator sites shared with SetR. Furthermore, we provide evidence of an additional operator site bound by SetR and CroS. Finally, we show that SetCD expression generates a positive feedback loop due to SXT excision and replication in a fraction of the cell population. Together, these results refine our understanding of the genetic regulation governing the propagation of major vectors of multidrug resistance.

IMPORTANCE

Healthcare systems worldwide are challenged by an alarming drug resistance crisis caused by the massive and rapid propagation of antibiotic resistance genes and the associated emergence of multidrug-resistant pathogenic bacteria. SXT/R391 ICEs contribute to this phenomenon not only in clinical and environmental vibrios but also in several members of the family Enterobacteriaceae. We have identified and characterized here the regulator CroS as a key factor in the stimulation of conjugative transfer of these ICEs in response to DNA-damaging agents. We have also untangled conflicting evidence regarding autoactivation of transfer by the master activator of SXT/R391 ICEs, SetCD. Discovery of CroS provides a clearer and more complete understanding of the regulatory network that governs the dissemination of SXT/R391 ICEs in bacterial populations.

Nonlinearities in protein space limit the utility of informatics in protein biophysics

We examine the utility of informatic-based methods in computational protein biophysics. To do so, we use newly developed metric functions to define completely independent sequence and structure spaces for a large database of proteins. By investigating the relationship between these spaces, we demonstrate quantitatively the limits of knowledge-based correlation between the sequences and structures of proteins. It is shown that there are well-defined, nonlinear regions of protein space in which dissimilar structures map onto similar sequences (the conformational switch), and dissimilar sequences map onto similar structures (remote homology). These nonlinearities are shown to be quite common-almost half the proteins in our database fall into one or the other of these two regions. They are not anomalies, but rather intrinsic properties of structural encoding in amino acid sequences. It follows that extreme care must be exercised in using bioinformatic data as a basis for computational structure prediction. The implications of these results for protein evolution are examined.

Smooth functional transition along a mutational pathway with an abrupt protein fold switch

Recent protein design experiments have demonstrated that proteins can migrate between folds through the accumulation of substitution mutations without visiting disordered or nonfunctional points in sequence space. To explore the biophysical mechanism underlying such transitions we use a three-letter continuous protein model with seven atoms per amino acid to provide realistic sequence-structure and sequence-function mappings through explicit simulation of the folding and interaction of model sequences. We start from two 16-amino-acid sequences folding into an α-helix and a β-hairpin, respectively, each of which has a preferred binding partner with 35 amino acids. We identify a mutational pathway between the two folds, which features a sharp fold switch. By contrast, we find that the transition in function is smooth. Moreover, the switch in preferred binding partner does not coincide with the fold switch. Discovery of new folds in evolution might therefore be facilitated by following fitness slopes in sequence space underpinned by binding-induced conformational switching.

Independent Co-Option of a Tailed Bacteriophage into a Killing Complex in Pseudomonas

Competition between microbes is widespread in nature, especially among those that are closely related. To combat competitors, bacteria have evolved numerous protein-based systems (bacteriocins) that kill strains closely related to the producer. In characterizing the bacteriocin complement and killing spectra for the model strain Pseudomonas syringae B728a, we discovered that its activity was not linked to any predicted bacteriocin but is derived from a prophage. Instead of encoding an active prophage, this region encodes a bacteriophage-derived bacteriocin, termed an R-type syringacin. This R-type syringacin is striking in its convergence with the well-studied R-type pyocin of P. aeruginosa in both genomic location and molecular function. Genomic alignment, amino acid percent sequence identity, and phylogenetic inference all support a scenario where the R-type syringacin has been co-opted independently of the R-type pyocin. Moreover, the presence of this region is conserved among several other Pseudomonas species and thus is likely important for intermicrobial interactions throughout this important genus.

Evolutionary innovation is often achieved through modification of complexes or processes for alternate purposes, termed co-option. Notable examples include the co-option of a structure functioning in locomotion (bacterial flagellum) to one functioning in protein secretion (type three secretion system). Similar co-options can occur independently in distinct lineages. We discovered a genomic region in the plant pathogen Pseudomonas syringae that consists of a fragment of a bacteriophage genome. The fragment encodes only the tail of the bacteriophage, which is lethal toward strains of this species. This structure is similar to a previously described structure produced by the related species Pseudomonas aeruginosa. The two structures, however, are not derived from the same evolutionary event. Thus, they represent independent bacteriophage co-options. The co-opted bacteriophage from P. syringae is found in the genomes of many other Pseudomonas species, suggesting ecological importance across this genus.

Interplay of Electrostatics and Hydrophobic Effects in the Metamorphic Protein Human Lymphotactin

The human lymphotactin (hLtn) is a protein that features two native states both of which are physiologically relevant: it is a monomer (hLtn10) at 10 °C with 200 mM salt and a dimer (hLtn40) at 40 °C and without salt. Here we focus on the networks of electrostatic and hydrophobic interactions that display substantial changes upon the conversion from hLtn10 to hLtn40 since they are expected to modulate the relative stability of the two folds. In addition to the Arg 23-Arg 43 interaction discussed in previous work, we find several other like-charge pairs that are likely important to the stability of hLtn10. Free energy perturbation calculations are carried out to explicitly evaluate the contribution of the Arg 23-Arg 43 interaction to the hLtn10 stability. hLtn40 features a larger number of salt bridges, and a set of hydrophobic residues undergo major changes in the solvent accessible surface area between hLtn10 and hLtn40, pointing to their importance to the relative stability of the two folds. We also discuss the use of explicit and implicit solvent simulations for characterizing the conformational ensembles under different solution conditions.

Studying protein fold evolution with hybrids of differently folded homologs

To study the sequence determinants governing protein fold evolution, we generated hybrid sequences from two homologous proteins with 40% identity but different folds: Pfl 6 Cro, which has a mixed α + β structure, and Xfaso 1 Cro, which has an all α-helical structure. First, we first examined eight chimeric hybrids in which the more structurally conserved N-terminal half of one protein was fused to the more structurally divergent C-terminal half of the other. None of these chimeras folded, as judged by circular dichroism spectra and thermal melts, suggesting that both halves have strong intrinsic preferences for the native global fold pattern, and/or that the interfaces between the halves are not readily interchangeable. Second, we examined 10 hybrids in which blocks of the structurally divergent C-terminal region were exchanged. These hybrids showed varying levels of thermal stability and suggested that the key residues in the Xfaso 1 C terminus specifying the all-α fold were concentrated near the end of helix 4 in Xfaso 1, which aligns to the end of strand 2 in Pfl 6. Finally, we generated hybrid substitutions for each individual residue in this critical region and measured thermal stabilities. The results suggested that R47 and V48 were the strongest factors that excluded formation of the α + β fold in the C-terminal region of Xfaso 1. In support of this idea, we found that the folding stability of one of the original eight chimeras could be rescued by back-substituting these two residues. Overall, the results show not only that the key factors for Cro fold specificity and evolution are global and multifarious, but also that some all-α Cro proteins have a C-terminal subdomain sequence within a few substitutions of switching to the α + β fold.

Effect of single-site mutations on hydrophobic-polar lattice proteins

We developed a heuristic method for determining the ground-state degeneracy of hydrophobic-polar (HP) lattice proteins, based on Wang-Landau and multicanonical sampling. It is applied during comprehensive studies of single-site mutations in specific HP proteins with different sequences. The effects in which we are interested include structural changes in ground states, changes of ground-state energy, degeneracy, and thermodynamic properties of the system. With respect to mutations, both extremely sensitive and insensitive positions in the HP sequence have been found. That is, ground-state energies and degeneracies, as well as other thermodynamic and structural quantities, may be either largely unaffected or may change significantly due to mutation.

Evolutionary relationship of two ancient protein superfolds

Proteins are the molecular machines of the cell that fold into specific three-dimensional structures to fulfill their functions. To improve our understanding of how the structure and function of proteins arises, it is crucial to understand how evolution has generated the structural diversity we observe today. Classically, proteins that adopt different folds are considered to be nonhomologous. However, using state-of-the-art tools for homology detection, we found evidence of homology between proteins of two ancient and highly populated protein folds, the (βα)8-barrel and the flavodoxin-like fold. We detected a family of sequences that show intermediate features between both folds and determined what is to our knowledge the first representative crystal structure of one of its members, giving new insights into the evolutionary link of two of the earliest folds. Our findings contribute to an emergent vision where protein superfolds share common ancestry and encourage further approaches to complete the mapping of structure space onto sequence space.

The energy landscape of a protein switch

Protein switches are made of highly similar sequences that fold to dramatically different structures. A structural switching system with 31 sequence variants for α and α+β folds has been illustrated experimentally by He et al., Structure, 2012, 20, 283 and is investigated computationally in the present study. Methods to assign a sequence to one of the two folds are reported and analyzed. A fast and accurate protocol to identify the correct fold of the 31 sequences is based on enriching modeled structures using short molecular dynamics (MD) trajectories and scoring these structures with coarse-grained energy functions. We examine five coarse-grained energy functions and illustrate that the Hinds-Levitt potential works the best for this task. We show that enrichment by MD significantly enhances prediction accuracy. Finally, we find that melting temperature correlates well with the energy difference between the two folds (correlation coefficient ∼-0.7). The correlation reduces dramatically (∼0.4) if the absolute energy of the correct fold is considered. Moreover, prediction of melting temperature is sensitive to the structural templates. We emphasize in our analyses the use of native structures as templates since these folds are more readily available from structural biology experiments.

SCOP2 prototype: a new approach to protein structure mining

We present a prototype of a new structural classification of proteins, SCOP2 (http://scop2.mrc-lmb.cam.ac.uk/), that we have developed recently. SCOP2 is a successor to the Structural Classification of Proteins (SCOP, http://scop.mrc-lmb.cam.ac.uk/scop/) database. Similarly to SCOP, the main focus of SCOP2 is to organize structurally characterized proteins according to their structural and evolutionary relationships. SCOP2 was designed to provide a more advanced framework for protein structure annotation and classification. It defines a new approach to the classification of proteins that is essentially different from SCOP, but retains its best features. The SCOP2 classification is described in terms of a directed acyclic graph in which nodes form a complex network of many-to-many relationships and are represented by a region of protein structure and sequence. The new classification project is expected to ensure new advances in the field and open new areas of research.

Packing interface energetics in different crystal forms of the λ Cro dimer

Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site-directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson-Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ∼5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures.

A role for indels in the evolution of Cro protein folds

Insertions and deletions in protein sequences, or indels, can disrupt structure and may result in changes in protein folds during evolution or in association with alternative splicing. Pfl 6 and Xfaso 1 are two proteins in the Cro family that share a common ancestor but have different folds. Sequence alignments of the two proteins show two gaps, one at the N terminus, where the sequence of Xfaso 1 is two residues shorter, and one near the center of the sequence, where the sequence of Pfl 6 is five residues shorter. To test the potential importance of indels in Cro protein evolution, we generated hybrid variants of Pfl 6 and Xfaso 1 with indels in one or both regions, chosen according to several plausible sequence alignments. All but one deletion variant completely unfolded both proteins, showing that a longer N-terminal sequence was critical for Pfl 6 folding and a longer central region sequence was critical for Xfaso 1 folding. By contrast, Xfaso 1 tolerated a longer N-terminal sequence with little destabilization, and Pfl 6 tolerated central region insertions, albeit with substantial effects on thermal stability and some perturbation of the surrounding structure. None of the mutations appeared to convert one stable fold into the other. On the basis of this two-protein comparison, short insertion and deletion mutations probably played a role in evolutionary fold change in the Cro family, but were also not the only factors.

A polymetamorphic protein

Arc repressor is a homodimeric protein with a ribbon-helix-helix fold. A single polar-to-hydrophobic substitution (N11L) at a solvent-exposed position leads to population of an alternate dimeric fold in which 3₁₀ helices replace a β-sheet. Here we find that the variant Q9V/N11L/R13V (S-VLV), with two additional polar-to-hydrophobic surface mutations in the same β-sheet, forms a highly stable, reversibly folded octamer with approximately half the α-helical content of wild-type Arc. At low protein concentration and low ionic strength, S-VLV also populates both dimeric topologies previously observed for N11L, as judged by NMR chemical shift comparisons. Thus, accumulation of simple hydrophobic mutations in Arc progressively reduces fold specificity, leading first to a sequence with two folds and then to a manifold bridge sequence with at least three different topologies. Residues 9-14 of S-VLV form a highly hydrophobic stretch that is predicted to be amyloidogenic, but we do not observe aggregates of higher order than octamer. Increases in sequence hydrophobicity can promote amyloid aggregation but also exert broader and more complex effects on fold specificity. Altered native folds, changes in fold coupled to oligomerization, toxic pre-amyloid oligomers, and amyloid fibrils may represent a near continuum of accessible alternatives in protein structure space.

Piecing together the family portrait of TCR-CD3 complexes

The pre-T-cell receptor (TCR)-, αβTCR-, and γδTCR-CD3 complexes are members of a family of modular biosensors that are responsible for driving T-cell development, activation, and effector functions. They inform essential checkpoint decisions by relaying key information from their ligand-binding modules (TCRs) to their signaling modules (CD3γε + CD3δε and CD3ζζ) and on to the intracellular signaling apparatus. Their actions shape the T-cell repertoire, as well as T-cell-mediated immunity; yet, the mechanisms that underlie their activity remain an enigma. As with any molecular machine, understanding how they function depends upon understanding how their parts fit and work together. In the 30 years since the initial biochemical and genetic characterizations of the αβTCR, the structure and function of the individual components of these family members have been extensively characterized. Cumulatively, this information has allowed us to piece together a portrait of the αβTCR-CD3 complex and outline the form of the remaining family members. Here we review the known structural and functional characteristics of the components of these TCR-CD3 complex family members. We then discuss how these data have informed our understanding of the architecture of the αβTCR-CD3 complex as well as their implications for the other family members. The intent is to provide a framework for considering: (i) how these thematically similar complexes diverge to execute their specific functions and (ii) how our knowledge of the form and function of these distinct family members can cross-inform our understanding of the other family members.

Protein conformational switches: from nature to design

Protein conformational switches alter their shape upon receiving an input signal, such as ligand binding, chemical modification, or change in environment. The apparent simplicity of this transformation--which can be carried out by a molecule as small as a thousand atoms or so--belies its critical importance to the life of the cell as well as its capacity for engineering by humans. In the realm of molecular switches, proteins are unique because they are capable of performing a variety of biological functions. Switchable proteins are therefore of high interest to the fields of biology, biotechnology, and medicine. These molecules are beginning to be exploited as the core machinery behind a new generation of biosensors, functionally regulated enzymes, and "smart" biomaterials that react to their surroundings. As inspirations for these designs, researchers continue to analyze existing examples of allosteric proteins. Recent years have also witnessed the development of new methodologies for introducing conformational change into proteins that previously had none. Herein we review examples of both natural and engineered protein switches in the context of four basic modes of conformational change: rigid-body domain movement, limited structural rearrangement, global fold switching, and folding-unfolding. Our purpose is to highlight examples that can potentially serve as platforms for the design of custom switches. Accordingly, we focus on inducible conformational changes that are substantial enough to produce a functional response (e.g., in a second protein to which it is fused), yet are relatively simple, structurally well-characterized, and amenable to protein engineering efforts.

Mutational tipping points for switching protein folds and functions

While disordered to ordered rearrangements are relatively common, the ability of proteins to switch from one ordered fold to a completely different fold is generally regarded as rare, and few fold switches have been characterized. Here, in a designed system, we examine the mutational requirements for transitioning between folds and functions. We show that switching between monomeric 3α and 4β+α folds can occur in multiple ways with successive single amino acid changes at diverse residue positions, raising the likelihood that such transitions occur in the evolution of new folds. Even mutations on the periphery of the core can tip the balance between alternatively folded states. Ligand-binding studies illustrate that a new immunoglobulin G-binding function can be gained well before the relevant 4β+α fold is appreciably populated in the unbound protein. The results provide new insights into the evolution of fold and function.

The porous borders of the protein world

Fold switching may play a role in the evolution of new protein folds and functions. He et al., in this issue of Structure, use protein design to illustrate that the same drastic change in a protein fold can occur via multiple different mutational pathways.

Effective techniques for protein structure mining

Retrieval and characterization of protein structure relationships are instrumental in a wide range of tasks in structural biology. The classification of protein structures (COPS) is a web service that provides efficient access to structure and sequence similarities for all currently available protein structures. Here, we focus on the application of COPS to the problem of template selection in homology modeling.

Molecular simulation uncovers the conformational space of the λ Cro dimer in solution

The significant variation among solved structures of the λ Cro dimer suggests its flexibility. However, contacts in the crystal lattice could have stabilized a conformation which is unrepresentative of its dominant solution form. Here we report on the conformational space of the Cro dimer in solution using replica exchange molecular dynamics in explicit solvent. The simulated ensemble shows remarkable correlation with available x-ray structures. Network analysis and a free energy surface reveal the predominance of closed and semi-open dimers, with a modest barrier separating these two states. The fully open conformation lies higher in free energy, indicating that it requires stabilization by DNA or crystal contacts. Most NMR models are found to be unstable conformations in solution. Intersubunit salt bridging between Arg(4) and Glu(53) during simulation stabilizes closed conformations. Because a semi-open state is among the low-energy conformations sampled in simulation, we propose that Cro-DNA binding may not entail a large conformational change relative to the dominant dimer forms in solution.

Reengineering Cro protein functional specificity with an evolutionary code

Cro proteins from different lambdoid bacteriophages are extremely variable in their target consensus DNA sequences and constitute an excellent model for evolution of transcription factor specificity. We experimentally tested a bioinformatically derived evolutionary code relating switches between pairs of amino acids at three recognition helix sites in Cro proteins to switches between pairs of nucleotide bases in the cognate consensus DNA half-sites. We generated all eight possible code variants of bacteriophage λ Cro and used electrophoretic mobility shift assays to compare binding of each variant to its own putative cognate site and to the wild-type cognate site; we also tested the wild-type protein against all eight DNA sites. Each code variant showed stronger binding to its putative cognate site than to the wild-type site, except some variants containing proline at position 27; each also bound its cognate site better than wild-type Cro bound the same site. Most code variants, however, displayed poorer affinity and specificity than wild-type λ Cro. Fluorescence anisotropy assays on λ Cro and the triple code variant (PSQ) against the two cognate sites confirmed the switch in specificity and showed larger apparent effects on binding affinity and specificity. Bacterial one-hybrid assays of λ Cro and PSQ against libraries of sequences with a single randomized half-site showed the expected switches in specificity at two of three coded positions and no clear switches in specificity at noncoded positions. With a few caveats, these results confirm that the proposed Cro evolutionary code can be used to reengineer Cro specificity.

Structural gymnastics of multifunctional metamorphic proteins

The classic structure-function paradigm holds that a protein exhibits a single well-defined native state that gives rise to its biological function. Nonetheless, over the past few decades, numerous examples of proteins exhibiting biological function arising from multiple structural states of varying disorder have been identified. Most recently, several examples of 'metamorphic proteins', able to interconvert between vastly different native-like topologies under physiological conditions, have been characterised with multiple functions. In this review, we look at the concept of protein metamorphosis in relation to the current understanding of the protein structure-function landscape. Although structural dynamism observed for metamorphic proteins provides a novel source of functional versatility, the dynamic nature of the metamorphic proteins generally makes them difficult to identify and probe using conventional protein structure determination methods. However, as the existence of metamorphic proteins has now been established and techniques enabling the analysis of multiple protein conformers are improving, it is likely that this class will continue to grow in number.

Spectral analysis of a protein conformational switch

The existence of conformational switching in proteins, induced by single amino acid mutations, presents an important challenge to our understanding of the physics of protein folding. Sequence-local methods, commonly used to detect structural homology, are incapable of accounting for this phenomenon. We examine a set of proteins, derived from the G(A) and G(B) domains of Streptococcus protein G, which are known to show a dramatic conformational change as a result of single-residue replacement. It is shown that these sequences, which are almost identical locally, can have very different global patterns of physical properties. These differences are consistent with the observed complete change in conformation. These results suggest that sequence-local methods for identifying structural homology can be misleading. They point to the importance of global sequence analysis in understanding sequence-structure relationships.

Evolutionary bridges to new protein folds: design of C-terminal Cro protein chameleon sequences

Regions of amino-acid sequence that are compatible with multiple folds may facilitate evolutionary transitions in protein structure. In a previous study, we described a heuristically designed chameleon sequence (SASF1, structurally ambivalent sequence fragment 1) that could adopt either of two naturally occurring conformations (α-helical or β-sheet) when incorporated as part of the C-terminal dimerization subdomain of two structurally divergent transcription factors, P22 Cro and λ Cro. Here we describe longer chameleon designs (SASF2 and SASF3) that in the case of SASF3 correspond to the full C-terminal half of the ordered region of a P22 Cro/λ Cro sequence alignment (residues 34-57). P22-SASF2 and λ(WDD)-SASF2 show moderate thermal stability in denaturation curves monitored by circular dichroism (T(m) values of 46 and 55°C, respectively), while P22-SASF3 and λ(WDD)-SASF3 have somewhat reduced stability (T(m) values of 33 and 49°C, respectively). (13)C and (1)H NMR secondary chemical shift analysis confirms two C-terminal α-helices for P22-SASF2 (residues 36-45 and 54-57) and two C-terminal β-strands for λ(WDD)-SASF2 (residues 40-45 and 50-52), corresponding to secondary structure locations in the two parent sequences. Backbone relaxation data show that both chameleon sequences have a relatively well-ordered structure. Comparisons of (15)N-(1)H correlation spectra for SASF2 and SASF3-containing proteins strongly suggest that SASF3 retains the chameleonism of SASF2. Both Cro C-terminal conformations can be encoded in a single sequence, showing the plausibility of linking different Cro folds by smooth evolutionary transitions. The N-terminal subdomain, though largely conserved in structure, also exerts an important contextual influence on the structure of the C-terminal region.

Measuring 1HN temperature coefficients in invisible protein states by relaxation dispersion NMR spectroscopy

A method based on the Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for measuring the temperature coefficients of amide proton chemical shifts of low populated 'invisible' protein states that exchange with a 'visible' ground state on the millisecond time-scale. The utility of the approach is demonstrated with an application to an I58D mutant of the Pfl6 Cro protein that undergoes exchange between the native, folded state and a cold denatured, unfolded conformational ensemble that is populated at a level of 6% at 2.5°C. A wide distribution of amide temperature coefficients is measured for the unfolded state. The distribution is centered about -5.6 ppb/K, consistent with an absence of intra-molecular hydrogen bonds, on average. However, the large range of values (standard deviation of 2.1 ppb/K) strongly supports the notion that the unfolded state of the protein is not a true random coil polypeptide chain.

Determinants of gas-phase disassembly behavior in homodimeric protein complexes with related yet divergent structures

The overall structure of a protein-protein complex reflects an intricate arrangement of noncovalent interactions. Whereas intramolecular interactions confer secondary and tertiary structure to individual subunits, intermolecular interactions lead to quaternary structure--the ordered aggregation of separate polypeptide chains into multisubunit assemblies. The specific ensemble of noncovalent contacts dictates the stability of subunit folds, enforces protein-protein binding specificity, and determines multimer stability. Consequently, noncovalent architecture is likely to play a role in the gas-phase dissociation of these assemblies during tandem mass spectrometry (MS/MS). To further advance the applicability of MS/MS to analytical problems in structural biology, a better understanding of the interplay between the structures and fragmentation behaviors of noncovalent protein complexes is essential. The present work constitutes a systematic study of model protein homodimers (bacteriophage N15 Cro, bacteriophage λ Cro, and bacteriophage P22 Arc) with related but divergent structures, both in terms of subunit folds and protein-protein interfaces. Because each of these dimers has a well-characterized structure (solution and/or crystal structure), specific noncovalent features could be correlated with gas-phase disassembly patterns as studied by collision-induced dissociation, surface-induced dissociation, and ion mobility. Of the several respects in which the dimers differed in structure, the presence or absence of intermolecular electrostatic contacts exerted the most significant influence on the gas-phase dissociation behavior. This is attributed to the well-known enhancement of ionic interactions in the absence of bulk solvent. Because salt bridges are general contributors to both intermolecular and intramolecular stability in protein complexes, these observations are broadly applicable to aid in the interpretation or prediction of dissociation spectra for noncovalent protein assemblies.

Converting a protein into a switch for biosensing and functional regulation

Proteins that switch conformations in response to a signaling event (e.g., ligand binding or chemical modification) present a unique solution to the design of reagent-free biosensors as well as molecules whose biological functions are regulated in useful ways. The principal roadblock in the path to develop such molecules is that the majority of natural proteins do not change conformation upon binding their cognate ligands or becoming chemically modified. Herein, we review recent protein engineering efforts to introduce switching properties into binding proteins. By co-opting natural allosteric coupling, joining proteins in creative ways and formulating altogether new switching mechanisms, researchers are learning how to coax conformational changes from proteins that previously had none. These studies are providing some answers to the challenging question: how can one convert a lock-and-key binding protein into a molecular switch?

A quality metric for homology modeling: the H-factor

BACKGROUND

The analysis of protein structures provides fundamental insight into most biochemical functions and consequently into the cause and possible treatment of diseases. As the structures of most known proteins cannot be solved experimentally for technical or sometimes simply for time constraints, in silico protein structure prediction is expected to step in and generate a more complete picture of the protein structure universe. Molecular modeling of protein structures is a fast growing field and tremendous works have been done since the publication of the very first model. The growth of modeling techniques and more specifically of those that rely on the existing experimental knowledge of protein structures is intimately linked to the developments of high resolution, experimental techniques such as NMR, X-ray crystallography and electron microscopy. This strong connection between experimental and in silico methods is however not devoid of criticisms and concerns among modelers as well as among experimentalists.

RESULTS

In this paper, we focus on homology-modeling and more specifically, we review how it is perceived by the structural biology community and what can be done to impress on the experimentalists that it can be a valuable resource to them. We review the common practices and provide a set of guidelines for building better models. For that purpose, we introduce the H-factor, a new indicator for assessing the quality of homology models, mimicking the R-factor in X-ray crystallography. The methods for computing the H-factor is fully described and validated on a series of test cases.

CONCLUSIONS

We have developed a web service for computing the H-factor for models of a protein structure. This service is freely accessible at http://koehllab.genomecenter.ucdavis.edu/toolkit/h-factor.