Rapid Evolution of Plethodontid Modulating Factor, a Hypervariable Salamander Courtship Pheromone, is Driven by Positive Selection

Sexual communication in plethodontid salamanders is mediated by a proteinaceous pheromone that a male delivers to a female during courtship, boosting her receptivity. The pheromone consists of three proteins from three unrelated protein families. These proteins are among a small group of pheromones known to affect female receptivity in vertebrates. Previously, we showed that the genes of two of these proteins (PRF and SPF) are prone to incessant evolution driven by positive selection, presumably as a consequence of coevolution with female receptors. In this report, we focus on the evolution of the third pheromone protein gene family, plethodontid modulating factor (PMF), to determine whether it shows the same pattern of diversification. We used RT-PCR in mental gland cDNA to survey PMF sequences from three genera of plethodontid salamanders (27 spp.) to measure rates of evolution, level of gene diversification, modes of selection, and types of amino acid substitution. Like PRF and SPF, PMF is produced by a multigene family characterized by gene duplication and high levels of polymorphism. PMF evolution is rapid, incessant, and driven by positive selection. PMF is more extreme in these dimensions than both PRF and SPF. Nestled within this extraordinary variation, however, is a signature of purifying selection, acting to preserve important structural and biochemical features of the PMF protein (i.e., secretion signal, cysteine residues, and pI). Although a pattern of persistent diversification exists at the molecular level, the morphological and behavioral aspects of the pheromone delivery system show evolutionary stasis over millions of years.     

Peptide Modulation of Class I Major Histocompatibility Complex Protein Molecular Flexibility and the Implications for Immune Recognition

T cells use the αβ T cell receptor (TCR) to recognize antigenic peptides presented by class I major histocompatibility complex proteins (pMHCs) on the surfaces of antigen-presenting cells. Flexibility in both TCRs and peptides plays an important role in antigen recognition and discrimination. Less clear is the role of flexibility in the MHC protein; although recent observations have indicated that mobility in the MHC can impact TCR recognition in a peptide-dependent fashion, the extent of this behavior is unknown. Here, using hydrogen/deuterium exchange, fluorescence anisotropy, and structural analyses, we show that the flexibility of the peptide binding groove of the class I MHC protein HLA-A*0201 varies significantly with different peptides. The variations extend throughout the binding groove, impacting regions contacted by TCRs as well as other activating and inhibitory receptors of the immune system. Our results are consistent with statistical mechanical models of protein structure and dynamics, in which the binding of different peptides alters the populations and exchange kinetics of substates in the MHC conformational ensemble. Altered MHC flexibility will influence receptor engagement, impacting conformational adaptations, entropic penalties associated with receptor recognition, and the populations of binding-competent states. Our results highlight a previously unrecognized aspect of the “altered self” mechanism of immune recognition and have implications for specificity, cross-reactivity, and antigenicity in cellular immunity.     

Tunicates Illuminate the Enigmatic Evolution of Chordate Metallothioneins by Gene Gains and Losses,Independent Modular Expansions,and Functional Convergences

To investigate novel patterns and processes of protein evolution, we have focused in the metallothioneins (MTs), a singular group of metal-binding, cysteine-rich proteins that, due to their high degree of sequence diversity, still represents a “black hole” in Evolutionary Biology. We have identified and analyzed more than 160 new MTs in nonvertebrate chordates (especially in 37 species of ascidians, 4 thaliaceans, and 3 appendicularians) showing that prototypic tunicate MTs are mono-modular proteins with a pervasive preference for cadmium ions, whereas vertebrate and cephalochordate MTs are bimodular proteins with diverse metal preferences. These structural and functional differences imply a complex evolutionary history of chordate MTs—including de novo emergence of genes and domains, processes of convergent evolution, events of gene gains and losses, and recurrent amplifications of functional domains—that would stand for an unprecedented case in the field of protein evolution.     

KRAB型锌指蛋白的进化及在物种演化中的功能

C2H2型锌指蛋白家族是目前发现的哺乳动物中最大的转录/转录调控因子家族,由一小群古老的含有真核锌指结构的转录因子经过多次基因复制和功能分化演化而来。KRAB型锌指蛋白(KRAB-containing zinc finger proteins, KRAB-ZFPs)作为C2H2型锌指蛋白家族中最大的亚家族,最早出现在四足脊椎动物,并随物种的进化数量快速增长,在人类中占据C2H2型锌指蛋白的60%左右。在物种演化中,进化压力主要改变KRAB-ZFPs的DNA结合能力,而KRAB-ZFPs介导的转录抑制能力则稳定存在。同时,多种KRAB-ZFPs能够与KRAB相关蛋白1(KRAB-associated protein 1, KAP1)协同作用沉默哺乳动物中反转录元件的活性,并与之协同进化,严格限制反转录原件的跳跃能力。本文综述了KRAB-ZFPs的数量倍增、锌指结构的灵活多变、KRAB-ZFPs/KAP1的转录抑制能力和反转录元件的跳跃性在促进哺乳动物调控网络的差异、基因组稳定性的变化和物种进化中的作用,旨在进一步揭示KRAB-ZFPs在推动物种稳定演化中的特点和功能。     

Crystal Structure of the Src Family Kinase Hck SH3-SH2 Linker Regulatory Region Supports an SH3-dominant Activation Mechanism

Most mammalian cell types depend on multiple Src family kinases (SFKs) to regulate diverse signaling pathways. Strict control of SFK activity is essential for normal cellular function, and loss of kinase regulation contributes to several forms of cancer and other diseases. Previous x-ray crystal structures of the SFKs c-Src and Hck revealed that intramolecular association of their Src homology (SH) 3 domains and SH2 kinase linker regions has a key role in down-regulation of kinase activity. However, the amino acid sequence of the Hck linker represents a suboptimal ligand for the isolated SH3 domain, suggesting that it may form the polyproline type II helical conformation required for SH3 docking only in the context of the intact structure. To test this hypothesis directly, we determined the crystal structure of a truncated Hck protein consisting of the SH2 and SH3 domains plus the linker. Despite the absence of the kinase domain, the structures and relative orientations of the SH2 and SH3 domains in this shorter protein were very similar to those observed in near full-length, down-regulated Hck. However, the SH2 kinase linker adopted a modified topology and failed to engage the SH3 domain. This new structure supports the idea that these noncatalytic regions work together as a “conformational switch” that modulates kinase activity in a manner unique to the SH3 domain and linker topologies present in the intact Hck protein. Our results also provide fresh structural insight into the facile induction of Hck activity by HIV-1 Nef and other Hck SH3 domain binding proteins and implicate the existence of innate conformational states unique to individual Src family members that “fine-tune” their sensitivities to activation by SH3-based ligands.     

ECOD: An Evolutionary Classification of Protein Domains

Understanding the evolution of a protein, including both close and distant relationships, often reveals insight into its structure and function. Fast and easy access to such up-to-date information facilitates research. We have developed a hierarchical evolutionary classification of all proteins with experimentally determined spatial structures, and presented it as an interactive and updatable online database. ECOD (Evolutionary Classification of protein Domains) is distinct from other structural classifications in that it groups domains primarily by evolutionary relationships (homology), rather than topology (or “fold”). This distinction highlights cases of homology between domains of differing topology to aid in understanding of protein structure evolution. ECOD uniquely emphasizes distantly related homologs that are difficult to detect, and thus catalogs the largest number of evolutionary links among structural domain classifications. Placing distant homologs together underscores the ancestral similarities of these proteins and draws attention to the most important regions of sequence and structure, as well as conserved functional sites. ECOD also recognizes closer sequence-based relationships between protein domains. Currently, approximately 100,000 protein structures are classified in ECOD into 9,000 sequence families clustered into close to 2,000 evolutionary groups. The classification is assisted by an automated pipeline that quickly and consistently classifies weekly releases of PDB structures and allows for continual updates. This synchronization with PDB uniquely distinguishes ECOD among all protein classifications. Finally, we present several case studies of homologous proteins not recorded in other classifications, illustrating the potential of how ECOD can be used to further biological and evolutionary studies.     

Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.     

The Crystal Structures of Yeast Get3 Suggest a Mechanism for Tail-Anchored Protein Membrane Insertion

Tail-anchored (TA) proteins represent a unique class of membrane proteins that contain a single C-terminal transmembrane helix. The post-translational insertion of the yeast TA proteins into the ER membrane requires the Golgi ER trafficking (GET) complex which contains Get1, Get2 and Get3. Get3 is an ATPase that recognizes and binds the C-terminal transmembrane domain (TMD) of the TA proteins. We have determined the crystal structures of Get3 from two yeast species, S. cerevisiae and D. hansenii, respectively. These high resolution crystal structures show that Get3 contains a nucleotide-binding domain and a “finger” domain for binding the TA protein TMD. A large hydrophobic groove on the finger domain of S. cerevisiae Get3 structure might represent the binding site for TMD of TA proteins. A hydrophobic helix from a symmetry-related Get3 molecule sits in the TMD-binding groove and mimics the TA binding scenario. Interestingly, the crystal structures of the Get3 dimers from S. cerevisiae and D. hansenii exhibit distinct conformations. The S. cerevisiae Get3 dimer structure does not contain nucleotides and maintains an “open” conformation, while the D. hansenii Get3 dimer structure binds ADP and stays in a “closed” conformation. We propose that the conformational changes to switch the Get3 between the open and closed conformations may facilitate the membrane insertions for TA proteins.     

Crystal structure and molecular dynamics of human POLDIP2, a multifaceted adaptor protein in metabolism and genome stability

Polymerase δ‐interacting protein 2 (POLDIP2, PDIP38) is a multifaceted, “moonlighting” protein, involved in binding protein partners from many different cellular processes, including mitochondrial metabolism and DNA replication and repair. How POLDIP2 interacts with many different proteins is unknown. Towards this goal, we present the crystal structure of POLDIP2 to 2.8 Å, which exhibited a compact two‐domain β‐strand‐rich globular structure, confirmed by circular dichroism and small angle X‐ray scattering approaches. POLDIP2 comprised canonical DUF525 and YccV domains, but with a conserved domain linker packed tightly, resulting in an “extended” YccV module. A central channel was observed, which we hypothesize could influence structural changes potentially mediated by redox conditions, following observation of a modified cysteine residue in the channel. Unstructured regions were rebuilt by ab initio modelling to generate a model of full‐length POLDIP2. Molecular dynamics simulations revealed a highly dynamic N‐terminal region tethered to the YccV‐domain by an extended linker, potentially facilitating interactions with distal binding partners. Models of POLDIP2 complexed with two of its partners, PrimPol and PCNA, indicated that dynamic flexibility of the POLDIP2 N‐terminus and loop regions likely mediate protein interactions.     

The Structure of the Cytomegalovirus-Encoded m04 Glycoprotein,a Prototypical Member of the m02 Family of Immunoevasins

The ability of CMVs to evade the immune system of the host is dependent on the expression of a wide array of glycoproteins, many of which interfere with natural killer cell function. In murine CMV, two large protein families mediate this immune-evasive function. Although it is established that the m145 family members mimic the structure of MHC-I molecules, the structure of the m02 family remains unknown. The most extensively studied m02 family member is m04, a glycoprotein that escorts newly assembled MHC-I molecules to the cell surface, presumably to avoid “missing self” recognition. Here we report the crystal structure of the m04 ectodomain, thereby providing insight into this large immunoevasin family. m04 adopted a β-sandwich immunoglobulin variable (Ig-V)-like fold, despite sharing very little sequence identity with the Ig-V superfamily. In addition to the Ig-V core, m04 possesses several unique structural features that included an unusual β-strand topology, a number of extended loops and a prominent α-helix. The m04 interior was packed by a myriad of hydrophobic residues that form distinct clusters around two conserved tryptophan residues. This hydrophobic core was well conserved throughout the m02 family, thereby indicating that murine CMV encodes a number of Ig-V-like molecules. We show that m04 binds a range of MHC-I molecules with low affinity in a peptide-independent manner. Accordingly, the structure of m04, which represents the first example of an murine CMV encoded Ig-V fold, provides a basis for understanding the structure and function of this enigmatic and large family of immunoevasins.     

The C-Terminal V5 Domain of Protein Kinase Cα Is Intrinsically Disordered,with Propensity to Associate with a Membrane Mimetic

The C-terminal V5 domain is one of the most variable domains in Protein Kinase C isoforms (PKCs). V5 confers isoform specificity on its parent enzyme through interactions with isoform-specific adaptor proteins and possibly through specific intra-molecular interactions with other PKC domains. The structural information about V5 domains in solution is sparse. The objective of this work was to determine the conformational preferences of the V5 domain from the α isoform of PKC (V5α) and evaluate its ability to associate with membrane mimetics. We show that V5α and its phosphorylation-mimicking variant, dmV5α, are intrinsically disordered protein domains. Phosphorylation-mimicking mutations do not alter the overall conformation of the polypeptide backbone, as evidenced by the local nature of chemical shift perturbations and the secondary structure propensity scores. However, the population of the “cis-trans” conformer of the Thr⁶³⁸-Pro⁶³⁹-Pro⁶⁴⁰ turn motif, which has been implicated in the down-regulation of PKCα via peptidyl-prolyl isomerase Pin1, increases in dmV5α, along with the conformational flexibility of the region between the turn and hydrophobic motifs. Both wild type and dmV5α associate with micelles made of a zwitterionic detergent, n-dodecylphosphocholine. Upon micelle binding, V5α acquires a higher propensity to form helical structures at the conserved “NFD” motif and the entire C-terminal third of the domain. The ability of V5α to partition into the hydrophobic micellar environment suggests that it may serve as a membrane anchor during the PKC maturation process.     

Protein evolution within a structural space

Understanding of the evolutionary origins of protein structures represents a key component of the understanding of molecular evolution as a whole. Here we seek to elucidate how the features of an underlying protein structural “space” might impact protein structural evolution. We approach this question using lattice polymers as a completely characterized model of this space. We develop a measure of structural comparison of lattice structures that is analogous to the one used to understand structural similarities between real proteins. We use this measure of structural relatedness to create a graph of lattice structures and compare this graph (in which nodes are lattice structures and edges are defined using structural similarity) to the graph obtained for real protein structures. We find that the graph obtained from all compact lattice structures exhibits a distribution of structural neighbors per node consistent with a random graph. We also find that subgraphs of 3500 nodes chosen either at random or according to physical constraints also represent random graphs. We develop a divergent evolution model based on the lattice space which produces graphs that, within certain parameter regimes, recapitulate the scale-free behavior observed in similar graphs of real protein structures.     

3D Bite Modeling and Feeding Mechanics of the Largest Living Amphibian,the Chinese Giant Salamander Andrias davidianus (Amphibia:Urodela)

Biting is an integral feature of the feeding mechanism for aquatic and terrestrial salamanders to capture, fix or immobilize elusive or struggling prey. However, little information is available on how it works and the functional implications of this biting system in amphibians although such approaches might be essential to understand feeding systems performed by early tetrapods. Herein, the skull biomechanics of the Chinese giant salamander, Andrias davidianus is investigated using 3D finite element analysis. The results reveal that the prey contact position is crucial for the structural performance of the skull, which is probably related to the lack of a bony bridge between the posterior end of the maxilla and the anterior quadrato-squamosal region. Giant salamanders perform asymmetrical strikes. These strikes are unusual and specialized behavior but might indeed be beneficial in such sit-and-wait or ambush-predators to capture laterally approaching prey. However, once captured by an asymmetrical strike, large, elusive and struggling prey have to be brought to the anterior jaw region to be subdued by a strong bite. Given their basal position within extant salamanders and their “conservative” morphology, cryptobranchids may be useful models to reconstruct the feeding ecology and biomechanics of different members of early tetrapods and amphibians, with similar osteological and myological constraints.     

BCL::Score—Knowledge Based Energy Potentials for Ranking Protein Models Represented by Idealized Secondary Structure Elements

The topology of most experimentally determined protein domains is defined by the relative arrangement of secondary structure elements, i.e. α-helices and β-strands, which make up 50–70% of the sequence. Pairing of β-strands defines the topology of β-sheets. The packing of side chains between α-helices and β-sheets defines the majority of the protein core. Often, limited experimental datasets restrain the position of secondary structure elements while lacking detail with respect to loop or side chain conformation. At the same time the regular structure and reduced flexibility of secondary structure elements make these interactions more predictable when compared to flexible loops and side chains. To determine the topology of the protein in such settings, we introduce a tailored knowledge-based energy function that evaluates arrangement of secondary structure elements only. Based on the amino acid C_β atom coordinates within secondary structure elements, potentials for amino acid pair distance, amino acid environment, secondary structure element packing, β-strand pairing, loop length, radius of gyration, contact order and secondary structure prediction agreement are defined. Separate penalty functions exclude conformations with clashes between amino acids or secondary structure elements and loops that cannot be closed. Each individual term discriminates for native-like protein structures. The composite potential significantly enriches for native-like models in three different databases of 10,000–12,000 protein models in 80–94% of the cases. The corresponding application, “BCL::ScoreProtein,” is available at www.meilerlab.org.     

Ca2+ Binding Enhanced Mechanical Stability of an Archaeal Crystallin

Structural topology plays an important role in protein mechanical stability. Proteins with β-sandwich topology consisting of Greek key structural motifs, for example, I27 of muscle titin and ¹⁰FNIII of fibronectin, are mechanically resistant as shown by single-molecule force spectroscopy (SMFS). In proteins with β-sandwich topology, if the terminal strands are directly connected by backbone H-bonding then this geometry can serve as a “mechanical clamp”. Proteins with this geometry are shown to have very high unfolding forces. Here, we set out to explore the mechanical properties of a protein, M-crystallin, which belongs to β-sandwich topology consisting of Greek key motifs but its overall structure lacks the “mechanical clamp” geometry at the termini. M-crystallin is a Ca²⁺ binding protein from Methanosarcina acetivorans that is evolutionarily related to the vertebrate eye lens β and γ-crystallins. We constructed an octamer of crystallin, (M-crystallin)₈, and using SMFS, we show that M-crystallin unfolds in a two-state manner with an unfolding force ∼90 pN (at a pulling speed of 1000 nm/sec), which is much lower than that of I27. Our study highlights that the β-sandwich topology proteins with a different strand-connectivity than that of I27 and ¹⁰FNIII, as well as lacking “mechanical clamp” geometry, can be mechanically resistant. Furthermore, Ca²⁺ binding not only stabilizes M-crystallin by 11.4 kcal/mol but also increases its unfolding force by ∼35 pN at the same pulling speed. The differences in the mechanical properties of apo and holo M-crystallins are further characterized using pulling speed dependent measurements and they show that Ca²⁺ binding reduces the unfolding potential width from 0.55 nm to 0.38 nm. These results are explained using a simple two-state unfolding energy landscape.     

Assessing the accuracy of ancestral protein reconstruction methods

The phylogenetic inference of ancestral protein sequences is a powerful technique for the study of molecular evolution, but any conclusions drawn from such studies are only as good as the accuracy of the reconstruction method. Every inference method leads to errors in the ancestral protein sequence, resulting in potentially misleading estimates of the ancestral protein's properties. To assess the accuracy of ancestral protein reconstruction methods, we performed computational population evolution simulations featuring near-neutral evolution under purifying selection, speciation, and divergence using an off-lattice protein model where fitness depends on the ability to be stable in a specified target structure. We were thus able to compare the thermodynamic properties of the true ancestral sequences with the properties of “ancestral sequences” inferred by maximum parsimony, maximum likelihood, and Bayesian methods. Surprisingly, we found that methods such as maximum parsimony and maximum likelihood that reconstruct a “best guess” amino acid at each position overestimate thermostability, while a Bayesian method that sometimes chooses less-probable residues from the posterior probability distribution does not. Maximum likelihood and maximum parsimony apparently tend to eliminate variants at a position that are slightly detrimental to structural stability simply because such detrimental variants are less frequent. Other properties of ancestral proteins might be similarly overestimated. This suggests that ancestral reconstruction studies require greater care to come to credible conclusions regarding functional evolution. Inferred functional patterns that mimic reconstruction bias should be reevaluated.