Publications

Background
A major bottleneck in our understanding of the molecular underpinnings of life is the assignment of function to proteins. While molecular experiments provide the most reliable annotation of proteins, their relatively low throughput and restricted purview have led to an increasing role for computational function prediction. However, assessing methods for protein function prediction and tracking progress in the field remain challenging.

Results
We conducted the second critical assessment of functional annotation (CAFA), a timed challenge to assess computational methods that automatically assign protein function. We evaluated 126 methods from 56 research groups for their ability to predict biological functions using Gene Ontology and gene-disease associations using Human Phenotype Ontology on a set of 3681 proteins from 18 species. CAFA2 featured expanded analysis compared with CAFA1, with regards to data set size, variety, and assessment metrics. To review progress in the field, the analysis compared the best methods from CAFA1 to those of CAFA2.

Conclusions
The top-performing methods in CAFA2 outperformed those from CAFA1. This increased accuracy can be attributed to a combination of the growing number of experimental annotations and improved methods for function prediction. The assessment also revealed that the definition of top-performing algorithms is ontology specific, that different performance metrics can be used to probe the nature of accurate predictions, and the relative diversity of predictions in the biological process and human phenotype ontologies. While there was methodological improvement between CAFA1 and CAFA2, the interpretation of results and usefulness of individual methods remain context-dependent.

Poly-proline type II (PPII) helical PXXP motifs are the recognition elements for a variety of protein–protein interactions that are critical for cellular signaling. Despite development of protocols for locking peptides into α-helical and β-strand conformations, there remains a lack of analogous methods for generating mimics of PPII helical structures. We describe herein a strategy to enforce PPII helical secondary structure in the 19-residue TrpPlexus miniature protein. Through sequence variation, we showed that a network of cation–π interactions could drive the formation of PPII helical conformations for both peptide and N-substituted glycine peptoid residues. The achievement of chemically diverse PPII helical scaffolds provides a new route towards discovering peptidomimetic inhibitors of protein–protein interactions mediated by PXXP motifs.

Protein secondary structures serve as geometrically constrained scaffolds for the display of key interacting residues at protein interfaces. Given the critical role of secondary structures in protein folding and the dependence of folding propensities on backbone dihedrals, secondary structure is expected to influence the identity of residues that are important for complex formation. Counter to this expectation, we find that a narrow set of residues dominates the binding energy in protein–protein complexes independent of backbone conformation. This finding suggests that the binding epitope may instead be substantially influenced by the side-chain conformations adopted. We analyzed side-chain conformational preferences in residues that contribute significantly to binding. This analysis suggests that preferred rotamers contribute directly to specificity in protein complex formation and provides guidelines for peptidomimetic inhibitor design.

The significant rotational energy barrier about the stereogenic carbon–carbon bond of axially chiral 3,3′-bipyrroles has been investigated by electronic circular dichroism (ECD) spectroscopy, time dependent HPLC analysis, and computational modeling. The results elucidate pathways and transition states involved in configurational inversion, thereby confirming that 3,3′-bipyrrole derivatives can exist in stable and isolable atropisomeric forms.

GIANT API provides biomedical researchers programmatic access to tissue-specific and global networks in humans and model organisms, and associated tools, which includes functional re-prioritization of existing genome-wide association study (GWAS) data. Using tissue-specific interaction networks, researchers are able to predict relationships between genes specific to a tissue or cell lineage, identify the changing roles of genes across tissues and uncover disease-gene associations. Additionally, GIANT API enables computational tools like NetWAS, which leverages tissue-specific networks for re-prioritization of GWAS results. The web services covered by the API include 144 tissue-specific functional gene networks in human, global functional networks for human and six common model organisms and the NetWAS method. GIANT API conforms to the REST architecture, which makes it stateless, cacheable and highly scalable. It can be used by a diverse range of clients including web browsers, command terminals, programming languages and standalone apps for data analysis and visualization. The API is freely available for use at http://giant-api.princeton.edu.

Metabolic network rewiring is the rerouting of metabolism through the use of alternate enzymes to adjust pathway flux and accomplish specific anabolic or catabolic objectives. Here, we report the first characterization of two parallel pathways for the breakdown of the short chain fatty acid propionate in Caenorhabditis elegans. Using genetic interaction mapping, gene co-expression analysis, pathway intermediate quantification and carbon tracing, we uncover a vitamin B12-independent propionate breakdown shunt that is transcriptionally activated on vitamin B12 deficient diets, or under genetic conditions mimicking the human diseases propionic- and methylmalonic acidemia, in which the canonical B12-dependent propionate breakdown pathway is blocked. Our study presents the first example of transcriptional vitamin-directed metabolic network rewiring to promote survival under vitamin deficiency. The ability to reroute propionate breakdown according to B12 availability may provide C. elegans with metabolic plasticity and thus a selective advantage on different diets in the wild.

During class switch recombination (CSR), B cells replace the Igh Cμ or δ exons with another downstream constant region exon (CH), altering the antibody isotype. CSR occurs through the introduction of AID-mediated double-strand breaks (DSBs) in switch regions and subsequent ligation of broken ends. Here, we developed an assay to investigate the dynamics of DSB formation in individual cells. We demonstrate that the upstream switch region Sμ is first targeted during recombination and that the mechanism underlying this control relies on 53BP1. Surprisingly, regulation of break order occurs through residual binding of 53BP1 to chromatin before the introduction of damage and independent of its established role in DNA repair. Using chromosome conformation capture, we show that 53BP1 mediates changes in chromatin architecture that affect break order. Finally, our results explain how changes in Igh architecture in the absence of 53BP1 could promote inversional rearrangements that compromise CSR.

Th17 lymphocytes protect mucosal barriers from infections, but also contribute to multiple chronic inflammatory diseases. Their differentiation is controlled by RORγt, a ligand-regulated nuclear receptor. We identified the DEAD-box RNA helicase DDX5 as a RORγt partner that coordinates transcription of selective Th17 genes and is required for Th17-mediated inflammatory pathologies. Surprisingly, the ability of DDX5 to interact with RORγt and co-activate its targets depends on its intrinsic RNA helicase activity and binding of a conserved nuclear long noncoding RNA (lncRNA), Rmrp, which is mutated in Cartilage-Hair Hypoplasia (CHH) patients. A targeted Rmrp mutation in mice, corresponding to one in CHH patients, abrogated the lncRNA’s chromatin recruitment, ability to potentiate DDX5-RORγt interaction and RORγt target gene transcription. Elucidation of the link between Rmrp and the DDX5-RORγt complex reveals a role for RNA helicases and lncRNAs in tissue-specific transcriptional regulation and promises new opportunities for therapeutic intervention in Th17-dependent diseases.

V(D)J recombination relies on the presence of proximal enhancers that activate the antigen receptor (AgR) loci in a lineage- and stage-specific manner. Unexpectedly, we find that both active and inactive AgR enhancers cooperate to disseminate their effects in a localized and long-range manner. Here, we demonstrate the importance of short-range contacts between active enhancers that constitute an Igk super-enhancer in B cells. Deletion of one element reduces the interaction frequency between other enhancers in the hub, which compromises the transcriptional output of each component. Furthermore, we establish that, in T cells, long-range contact and cooperation between the inactive Igk enhancer MiEκ and the active Tcrb enhancer Eβ alters enrichment of CBFβ binding in a manner that impacts Tcrb recombination. These findings underline the complexities of enhancer regulation and point to a role for localized and long-range enhancer-sharing between active and inactive elements in lineage- and stage-specific control.

Selecting regularization parameters in penalized high-dimensional graphical models in a principled, data-driven, and computationally efficient manner continues to be one of the key challenges in high-dimensional statistics. We present substantial computational gains and conceptual generalizations of the Stability Approach to Regularization Selection (StARS), a state-of-the-art graphical model selection scheme. Using properties of the Poisson-Binomial distribution and convex non-asymptotic distributional modeling we propose lower and upper bounds on the StARS graph regularization path which results in greatly reduced computational cost without compromising regularization selection. We also generalize the StARS criterion from single edge to induced subgraph (graphlet) stability. We show that simultaneously requiring edge and graphlet stability leads to superior graph recovery performance independent of graph topology. These novel insights render Gaussian graphical model selection a routine task on standard multi-core computers.