Update of medRxiv. 2020 Jul 19;:

Proteomic technologies now enable the rapid quantification of thousands of proteins across genetically diverse samples. Integration of these data with systems-genetics analyses is a powerful approach to identify new regulators of economically important or disease-relevant phenotypes in various populations. In this review, we summarize the latest proteomic technologies and discuss technical challenges for their use in population studies. We demonstrate how the analysis of correlation structure and loci mapping can be used to identify genetic factors regulating functional protein networks and complex traits. Finally, we provide an extensive summary of the use of proteome-wide systems genetics throughout fungi, plant, and animal kingdoms and discuss the power of this approach to identify candidate regulators and drug targets in large human consortium studies.

Machine learning can infer how protein sequence maps to function without requiring a detailed understanding of the underlying physical or biological mechanisms. It is challenging to apply existing supervised learning frameworks to large-scale experimental data generated by deep mutational scanning (DMS) and related methods. DMS data often contain high-dimensional and correlated sequence variables, experimental sampling error and bias, and the presence of missing data. Notably, most DMS data do not contain examples of negative sequences, making it challenging to directly estimate how sequence affects function. Here, we develop a positive-unlabeled (PU) learning framework to infer sequence-function relationships from large-scale DMS data. Our PU learning method displays excellent predictive performance across ten large-scale sequence-function datasets, representing proteins of different folds, functions, and library types. The estimated parameters pinpoint key residues that dictate protein structure and function. Finally, we apply our statistical sequence-function model to design highly stabilized enzymes.

Comment in Trends Genet. 2021 May;37(5):403-404.

Redox cofactor production is integral toward antioxidant generation, clearance of reactive oxygen species, and overall tumor response to ionizing radiation treatment. To identify systems-level alterations in redox metabolism that confer resistance to radiation therapy, we developed a bioinformatics pipeline for integrating multi-omics data into personalized genome-scale flux balance analysis models of 716 radiation-sensitive and 199 radiation-resistant tumors. These models collectively predicted that radiation-resistant tumors reroute metabolic flux to increase mitochondrial NADPH stores and reactive oxygen species (ROS) scavenging. Simulated genome-wide knockout screens agreed with experimental siRNA gene knockdowns in matched radiation-sensitive and radiation-resistant cancer cell lines, revealing gene targets involved in mitochondrial NADPH production, central carbon metabolism, and folate metabolism that allow for selective inhibition of glutathione production and H2O2 clearance in radiation-resistant cancers. This systems approach represents a significant advancement in developing quantitative genome-scale models of redox metabolism and identifying personalized metabolic targets for improving radiation sensitivity in individual cancer patients.

Enzymes maintain metabolism, and their concentration affects cellular fitness: high enzyme levels are costly, and low enzyme levels can limit metabolic flux. Here, we used CRISPR interference (CRISPRi) to study the consequences of decreasing E. coli enzymes below wild-type levels. A pooled CRISPRi screen with 7,177 strains demonstrates that metabolism buffers fitness defects for hours after the induction of CRISPRi. We characterized the metabolome and proteome responses in 30 CRISPRi strains and elucidated three gene-specific buffering mechanisms: ornithine buffered the knockdown of carbamoyl phosphate synthetase (CarAB) by increasing CarAB activity, S-adenosylmethionine buffered the knockdown of homocysteine transmethylase (MetE) by de-repressing expression of the methionine pathway, and 6-phosphogluconate buffered the knockdown of 6-phosphogluconate dehydrogenase (Gnd) by activating a bypass. In total, this work demonstrates that CRISPRi screens can reveal global sources of metabolic robustness and identify local regulatory mechanisms that buffer decreases of specific enzymes. A record of this paper's transparent peer review process is included in the Supplemental Information.

Subunit vaccines induce immunity to a pathogen by presenting a component of the pathogen and thus inherently limit the representation of pathogen peptides for cellular immunity-based memory. We find that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) subunit peptides may not be robustly displayed by the major histocompatibility complex (MHC) molecules in certain individuals. We introduce an augmentation strategy for subunit vaccines that adds a small number of SARS-CoV-2 peptides to a vaccine to improve the population coverage of pathogen peptide display. Our population coverage estimates integrate clinical data on peptide immunogenicity in convalescent COVID-19 patients and machine learning predictions. We evaluate the population coverage of 9 different subunits of SARS-CoV-2, including 5 functional domains and 4 full proteins, and augment each of them to fill a predicted coverage gap.

Viruses deploy genetically encoded strategies to coopt host machinery and support viral replicative cycles. Here, we use protein structure similarity to scan for molecular mimicry, manifested by structural similarity between viral and endogenous host proteins, across thousands of cataloged viruses and hosts spanning broad ecological niches and taxonomic range, including bacteria, plants and fungi, invertebrates, and vertebrates. This survey identified over 6,000,000 instances of structural mimicry; more than 70% of viral mimics cannot be discerned through protein sequence alone. We demonstrate that the manner and degree to which viruses exploit molecular mimicry varies by genome size and nucleic acid type and identify 158 human proteins that are mimicked by coronaviruses, providing clues about cellular processes driving pathogenesis. Our observations point to molecular mimicry as a pervasive strategy employed by viruses and indicate that the protein structure space used by a given virus is dictated by the host proteome. A record of this paper's transparent peer review process is included in the Supplemental Information.