A new scoring function and associated statistical significance for structure alignment by CE

Time-resolved comparative molecular evolution of oxygenic photosynthesis

Oxygenic photosynthesis starts with the oxidation of water to O2, a light-driven reaction catalysed by photosystem II. Cyanobacteria are the only prokaryotes capable of water oxidation and therefore, it is assumed that the origin of oxygenic photosynthesis is a late innovation relative to the origin of life and bioenergetics. However, when exactly water oxidation originated remains an unanswered question. Here we use phylogenetic analysis to study a gene duplication event that is unique to photosystem II: the duplication that led to the evolution of the core antenna subunits CP43 and CP47. We compare the changes in the rates of evolution of this duplication with those of some of the oldest well-described events in the history of life: namely, the duplication leading to the Alpha and Beta subunits of the catalytic head of ATP synthase, and the divergence of archaeal and bacterial RNA polymerases and ribosomes. We also compare it with more recent events such as the duplication of Cyanobacteria-specific FtsH metalloprotease subunits and the radiation leading to Margulisbacteria, Sericytochromatia, Vampirovibrionia, and other clades containing anoxygenic phototrophs. We demonstrate that the ancestral core duplication of photosystem II exhibits patterns in the rates of protein evolution through geological time that are nearly identical to those of the ATP synthase, RNA polymerase, or the ribosome. Furthermore, we use ancestral sequence reconstruction in combination with comparative structural biology of photosystem subunits, to provide additional evidence supporting the premise that water oxidation had originated before the ancestral core duplications. Our work suggests that photosynthetic water oxidation originated closer to the origin of life and bioenergetics than can be documented based on phylogenetic or phylogenomic species trees alone.

Early Archean origin of Photosystem II

Photosystem II is a photochemical reaction center that catalyzes the light-driven oxidation of water to molecular oxygen. Water oxidation is the distinctive photochemical reaction that permitted the evolution of oxygenic photosynthesis and the eventual rise of eukaryotes. At what point during the history of life an ancestral photosystem evolved the capacity to oxidize water still remains unknown. Here, we study the evolution of the core reaction center proteins of Photosystem II using sequence and structural comparisons in combination with Bayesian relaxed molecular clocks. Our results indicate that a homodimeric photosystem with sufficient oxidizing power to split water had already appeared in the early Archean about a billion years before the most recent common ancestor of all described Cyanobacteria capable of oxygenic photosynthesis, and well before the diversification of some of the known groups of anoxygenic photosynthetic bacteria. Based on a structural and functional rationale, we hypothesize that this early Archean photosystem was capable of water oxidation to oxygen and had already evolved protection mechanisms against the formation of reactive oxygen species. This would place primordial forms of oxygenic photosynthesis at a very early stage in the evolutionary history of life.

Synthetic STARR-seq reveals how DNA shape and sequence modulate transcriptional output and noise

The binding of transcription factors to short recognition sequences plays a pivotal role in controlling the expression of genes. The sequence and shape characteristics of binding sites influence DNA binding specificity and have also been implicated in modulating the activity of transcription factors downstream of binding. To quantitatively assess the transcriptional activity of tens of thousands of designed synthetic sites in parallel, we developed a synthetic version of STARR-seq (synSTARR-seq). We used the approach to systematically analyze how variations in the recognition sequence of the glucocorticoid receptor (GR) affect transcriptional regulation. Our approach resulted in the identification of a novel highly active functional GR binding sequence and revealed that sequence variation both within and flanking GR’s core binding site can modulate GR activity without apparent changes in DNA binding affinity. Notably, we found that the sequence composition of variants with similar activity profiles was highly diverse. In contrast, groups of variants with similar activity profiles showed specific DNA shape characteristics indicating that DNA shape may be a better predictor of activity than DNA sequence. Finally, using single cell experiments with individual enhancer variants, we obtained clues indicating that the architecture of the response element can independently tune expression mean and cell-to cell variability in gene expression (noise). Together, our studies establish synSTARR as a powerful method to systematically study how DNA sequence and shape modulate transcriptional output and noise.

The expression level of genes is controlled by transcription factors, which are proteins that bind to genomic response elements that contain their recognition DNA sequence. Importantly, genes are not simply turned on but need to be expressed at the right level. This is, at least in part, assured by the sequence composition of genomic response elements. Here, we studied how the recognition DNA sequence influences gene regulation by a transcription factor called the glucocorticoid receptor. Specifically, we developed a method to test the activity of variants in a highly parallelized setting where everything is kept identical except for the sequence of the binding site. The systematic analysis of tens of thousands of sequence variants facilitated the identification of a previously unknown sequence variant with high activity. Moreover, we report how sequence variation of the response element influences cell-to-cell variability in expression levels. Finally, we observe similar activity profiles for distinct sequence variants that share similar three-dimensional DNA shape characteristics arguing that the three-dimensional perception of DNA by the glucocorticoid receptor, modulates its activity towards individual target genes.

Discovery of peptide inhibitors targeting human programmed death 1 (PD-1) receptor

Blocking the interaction of human programmed death 1 (hPD-1) and its ligand hPD-L1 has been a promising immunotherapy in cancer treatment. In this paper, using a computational de novo peptide design method, we designed several hPD-1 binding peptides. The most potent peptide Ar5Y_4 showed a KD value of 1.38 ± 0.39 μM, comparable to the binding affinity of the cognate hPD-L1. A Surface Plasmon Resonance (SPR) competitive binding assay result indicated that Ar5Y_4 could inhibit the interaction of hPD-1/hPD-L1. Moreover, Ar5Y_4 could restore the function of Jurkat T cells which had been suppressed by stimulated HCT116 cells. Peptides described in this paper provide promising biologic candidates for cancer immunotherapy or diagnostics.

A fresh look at the evolution and diversification of photochemical reaction centers

In this review, I reexamine the origin and diversification of photochemical reaction centers based on the known phylogenetic relations of the core subunits, and with the aid of sequence and structural alignments. I show, for example, that the protein folds at the C-terminus of the D1 and D2 subunits of Photosystem II, which are essential for the coordination of the water-oxidizing complex, were already in place in the most ancestral Type II reaction center subunit. I then evaluate the evolution of reaction centers in the context of the rise and expansion of the different groups of bacteria based on recent large-scale phylogenetic analyses. I find that the Heliobacteriaceae family of Firmicutes appears to be the earliest branching of the known groups of phototrophic bacteria; however, the origin of photochemical reaction centers and chlorophyll synthesis cannot be placed in this group. Moreover, it becomes evident that the Acidobacteria and the Proteobacteria shared a more recent common phototrophic ancestor, and this is also likely for the Chloroflexi and the Cyanobacteria. Finally, I argue that the discrepancies among the phylogenies of the reaction center proteins, chlorophyll synthesis enzymes, and the species tree of bacteria are best explained if both types of photochemical reaction centers evolved before the diversification of the known phyla of phototrophic bacteria. The primordial phototrophic ancestor must have had both Type I and Type II reaction centers.

Structure of the human MLH1 N-terminus: implications for predisposition to Lynch syndrome

Mismatch repair prevents the accumulation of erroneous insertions/deletions and non-Watson-Crick base pairs in the genome. Pathogenic mutations in the MLH1 gene are associated with a predisposition to Lynch and Turcot's syndromes. Although genetic testing for these mutations is available, robust classification of variants requires strong clinical and functional support. Here, the first structure of the N-terminus of human MLH1, determined by X-ray crystallography, is described. The structure shares a high degree of similarity with previously determined prokaryotic MLH1 homologs; however, this structure affords a more accurate platform for the classification of MLH1 variants.

Solution structure and properties of AlgH from Pseudomonas aeruginosa

In Pseudomonas aeruginosa, the algH gene regulates the cellular concentrations of a number of enzymes and the production of several virulence factors, and is suggested to serve a global regulatory function. The precise mechanism by which the algH gene product, the AlgH protein, functions is unknown. The same is true for AlgH family members from other bacteria. In order to lay the groundwork for understanding the physical underpinnings of AlgH function, we examined the structure and physical properties of AlgH in solution. Under reducing conditions, results of NMR, electrophoretic mobility, and sedimentation equilibrium experiments indicate AlgH is predominantly monomeric and monodisperse in solution. Under nonreducing conditions intra and intermolecular disulfide bonds form, the latter promoting AlgH oligomerization. The high-resolution solution structure of AlgH reveals alpha/beta-sandwich architecture fashioned from ten beta strands and seven alpha helices. Comparison with available structures of orthologues indicates conservation of overall structural topology. The region of the protein most strongly conserved structurally also shows the highest amino acid sequence conservation and, as revealed by hydrogen-deuterium exchange studies, is also the most stable. In this region, evolutionary trace analysis identifies two clusters of amino acid residues with the highest evolutionary importance relative to all other AlgH residues. These frame a partially solvent exposed shallow hydrophobic cleft, perhaps identifying a site for intermolecular interactions. The results establish a physical foundation for understanding the structure and function of AlgH and AlgH family proteins and should be of general importance for further investigations of these and related proteins.

Origin and Evolution of Water Oxidation before the Last Common Ancestor of the Cyanobacteria

Photosystem II, the water oxidizing enzyme, altered the course of evolution by filling the atmosphere with oxygen. Here, we reconstruct the origin and evolution of water oxidation at an unprecedented level of detail by studying the phylogeny of all D1 subunits, the main protein coordinating the water oxidizing cluster (Mn₄CaO₅) of Photosystem II. We show that D1 exists in several forms making well-defined clades, some of which could have evolved before the origin of water oxidation and presenting many atypical characteristics. The most ancient form is found in the genome of Gloeobacter kilaueensis JS-1 and this has a C-terminus with a higher sequence identity to D2 than to any other D1. Two other groups of early evolving D1 correspond to those expressed under prolonged far-red illumination and in darkness. These atypical D1 forms are characterized by a dramatically different Mn₄CaO₅ binding site and a Photosystem II containing such a site may assemble an unconventional metal cluster. The first D1 forms with a full set of ligands to the Mn₄CaO₅ cluster are grouped with D1 proteins expressed only under low oxygen concentrations and the latest evolving form is the dominant type of D1 found in all cyanobacteria and plastids. In addition, we show that the plastid ancestor had a D1 more similar to those in early branching Synechococcus. We suggest each one of these forms of D1 originated from transitional forms at different stages toward the innovation and optimization of water oxidation before the last common ancestor of all known cyanobacteria.

Solution structure of the YTH domain in complex with N6-methyladenosine RNA: a reader of methylated RNA

N(6)A methylation is the most abundant RNA modification occurring within messenger RNA. Impairment of methylase or demethylase functions are associated with severe phenotypes and diseases in several organisms. Beside writer and eraser enzymes of this dynamic RNA epigenetic modification, reader proteins that recognize this modification are involved in numerous cellular processes. Although the precise characterization of these reader proteins remains unknown, preliminary data showed that most potential reader proteins contained a conserved YT521-B homology (YTH) domain. Here we define the YTH domain of rat YT521-B as a N(6)-methylated adenosine reader domain and report its solution structure in complex with a N(6)-methylated RNA. The structure reveals a binding preference for NGANNN RNA hexamer and a deep hydrophobic cleft for m(6)A recognition. These findings establish a molecular function for YTH domains as m(6)A reader domains and should guide further studies into the biological functions of YTH-containing proteins in m(6)A recognition.

Systematic detection of internal symmetry in proteins using CE-Symm

Symmetry is an important feature of protein tertiary and quaternary structures that has been associated with protein folding, function, evolution, and stability. Its emergence and ensuing prevalence has been attributed to gene duplications, fusion events, and subsequent evolutionary drift in sequence. This process maintains structural similarity and is further supported by this study. To further investigate the question of how internal symmetry evolved, how symmetry and function are related, and the overall frequency of internal symmetry, we developed an algorithm, CE-Symm, to detect pseudo-symmetry within the tertiary structure of protein chains. Using a large manually curated benchmark of 1007 protein domains, we show that CE-Symm performs significantly better than previous approaches. We use CE-Symm to build a census of symmetry among domain superfamilies in SCOP and note that 18% of all superfamilies are pseudo-symmetric. Our results indicate that more domains are pseudo-symmetric than previously estimated. We establish a number of recurring types of symmetry-function relationships and describe several characteristic cases in detail. With the use of the Enzyme Commission classification, symmetry was found to be enriched in some enzyme classes but depleted in others. CE-Symm thus provides a methodology for a more complete and detailed study of the role of symmetry in tertiary protein structure [availability: CE-Symm can be run from the Web at http://source.rcsb.org/jfatcatserver/symmetry.jsp. Source code and software binaries are also available under the GNU Lesser General Public License (version 2.1) at https://github.com/rcsb/symmetry. An interactive census of domains identified as symmetric by CE-Symm is available from http://source.rcsb.org/jfatcatserver/scopResults.jsp].

A horizontal alignment tool for numerical trend discovery in sequence data: application to protein hydropathy

An algorithm is presented that returns the optimal pairwise gapped alignment of two sets of signed numerical sequence values. One distinguishing feature of this algorithm is a flexible comparison engine (based on both relative shape and absolute similarity measures) that does not rely on explicit gap penalties. Additionally, an empirical probability model is developed to estimate the significance of the returned alignment with respect to randomized data. The algorithm's utility for biological hypothesis formulation is demonstrated with test cases including database search and pairwise alignment of protein hydropathy. However, the algorithm and probability model could possibly be extended to accommodate other diverse types of protein or nucleic acid data, including positional thermodynamic stability and mRNA translation efficiency. The algorithm requires only numerical values as input and will readily compare data other than protein hydropathy. The tool is therefore expected to complement, rather than replace, existing sequence and structure based tools and may inform medical discovery, as exemplified by proposed similarity between a chlamydial ORFan protein and bacterial colicin pore-forming domain. The source code, documentation, and a basic web-server application are available.

An all-atom model of the structure of human copper transporter 1

Human copper transporter 1 (hCTR1) is the major high affinity copper influx transporter in mammalian cells that also mediates uptake of the cancer chemotherapeutic agent cisplatin. A low resolution structure of hCTR1 determined by cryoelectron microscopy was recently published. Several protein structure simulation techniques were used to create an all-atom model of this important transporter using the low resolution structure as a starting point. The all-atom model provides new insights into the roles of specific residues of the N-terminal extracellular domain, the intracellular loop, and C-terminal region in metal ion transport. In particular, the model demonstrates that the central region of the pore contains four sets of methionine triads in the intramembranous region. The structure confirms that two triads of methionine residues delineate the intramembranous region of the transporter, and further identifies two additional methionine triads that are located in the extracellular N-terminal part of the transporter. Together, the four triads create a structure that promotes stepwise transport of metal ions into and then through the intramembranous channel of the transporter via transient thioether bonds to methionine residues. Putative copper-binding sites in the hCTR1 trimer were identified by a program developed by us for prediction of metal-binding sites. These sites correspond well with the known effects of mutations on the ability of the protein to transport copper and cisplatin.

Molecular architecture of the Nup84-Nup145C-Sec13 edge element in the nuclear pore complex lattice

Nuclear pore complexes (NPCs) facilitate all nucleocytoplasmic transport. These massive protein assemblies are modular, with a stable structural scaffold supporting more dynamically attached components. The scaffold is made from multiple copies of the heptameric Y complex and the heteromeric Nic96 complex. We previously showed that members of these core subcomplexes specifically share an ACE1 fold with Sec31 of the COPII vesicle coat, and we proposed a lattice model for the NPC based on this commonality. Here we present the crystal structure of the heterotrimeric 134-kDa complex of Nup84-Nup145C-Sec13 of the Y complex. The heterotypic ACE1 interaction of Nup84 and Nup145C is analogous to the homotypic ACE1 interaction of Sec31 that forms COPII lattice edge elements and is inconsistent with the alternative 'fence-like' NPC model. We construct a molecular model of the Y complex and compare the architectural principles of COPII and NPC lattices.

Chapter 4. Predicting and characterizing protein functions through matching geometric and evolutionary patterns of binding surfaces

Predicting protein functions from structures is an important and challenging task. Although proteins are often thought to be packed as tightly as solids, closer examination based on geometric computation reveals that they contain numerous voids and pockets. Most of them are of random nature, but some are binding sites providing surfaces to interact with other molecules. A promising approach for function prediction is to infer functions through discovery of similarity in local binding pockets, as proteins binding to similar substrates/ligands and carrying out similar functions have similar physical constraints for binding and reactions. In this chapter, we describe computational methods to distinguish those surface pockets that are likely to be involved in important biological functions, and methods to identify key residues in these pockets. We further describe how to predict protein functions at large scale from structures by detecting binding surfaces similar in residue make-ups, shape, and orientation. We also describe a Bayesian Monte Carlo method that can separate selection pressure due to biological function from pressure due to protein folding. We show how this method can be used to reconstruct the evolutionary history of binding surfaces for detecting similar binding surfaces. In addition, we briefly discuss how the negative image of a binding pocket can be casted, and how such information can be used to facilitate drug discovery.

Energetic profiling of protein folds

Current protein classification methods treat high-resolution structures as static entities. However, experiments have well documented the dynamic nature of proteins. With knowledge that thermodynamic fluctuations around the high-resolution structure contribute to a more physically accurate and biologically meaningful picture of a protein, the concept of a protein's energetic profile is introduced. It is demonstrated on a large scale that energetic profiles are both diagnostic of a protein fold and evolutionarily relevant. Development of Structural Thermodynamic Ensemble-based Protein Homology (STEPH), an algorithm that searches for local similarities between energetic profiles, constitutes a first step towards a long-term goal of our laboratory to integrate thermodynamic information into protein-fold classification approaches.

Multiobjective optimization in bioinformatics and computational biology

This paper reviews the application of multiobjective optimization in the fields of bioinformatics and computational biology. A survey of existing work, organized by application area, forms the main body of the review, following an introduction to the key concepts in multiobjective optimization. An original contribution of the review is the identification of five distinct "contexts," giving rise to multiple objectives: These are used to explain the reasons behind the use of multiobjective optimization in each application area and also to point the way to potential future uses of the technique.