Cited by 114Protein Structure and DynamicsMachine Learning in BioinformaticsGenomics and Phylogenetic Studies

The DOI matches its Crossref record; 57 of the 57 references it cites were found in a registry or an index.

Large-scale discovery of protein interactions at residue resolution using co-evolution calculated from genomic sequences

Open access

Abstract

Increasing numbers of protein interactions have been identified in high-throughput experiments, but only a small proportion have solved structures. Recently, sequence coevolution-based approaches have led to a breakthrough in predicting monomer protein structures and protein interaction interfaces. Here, we address the challenges of large-scale interaction prediction at residue resolution with a fast alignment concatenation method and a probabilistic score for the interaction of residues. Importantly, this method (EVcomplex2) is able to assess the likelihood of a protein interaction, as we show here applied to large-scale experimental datasets where the pairwise interactions are unknown. We predict 504 interactions de novo in the E. coli membrane proteome, including 243 that are newly discovered. While EVcomplex2 does not require available structures, coevolving residue pairs can be used to produce structural models of protein interactions, as done here for membrane complexes including the Flagellar Hook-Filament Junction and the Tol/Pal complex.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41467-021-21636-z

Cites

57 of 57 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Scaling EVcouplings methods to full bacterial genomes. A The search problem for binary protein-protein interactions in Escherichia coli involves finding all of the estimated 104 true interactions out of 107 possible pairs. Only approximately 9% of these true direct interactions have a crystal structure solved in E. coli or a homologous structure in another organism. B Evolutionary couplings learned directly from protein sequences can resolve interfaces. Sequence alignments of both monomeric proteins are created and concatenated by the reciprocal highest identity procedure before inference of evolutionary couplings. Raw evolutionary coupling scores can be combined with features of their distribution, biochemical properties, and sequence entropy to improve inference. C A benchmark dataset of all non-redundant protein interactions with known interface structure.
Fig. 1 Scaling EVcouplings methods to full bacterial genomes. A The search problem for binary protein-protein interactions in Escherichia coli involves finding all of the estimated 104 true interactions out of 107 possible pairs. Only approximately 9% of these true direct interactions have a crystal structure solved in E. coli or a homologous structure in another organism. B Evolutionary couplings learned directly from protein sequences can resolve interfaces. Sequence alignments of both monomeric proteins are created and concatenated by the reciprocal highest identity procedure before inference of evolutionary couplings. Raw evolutionary coupling scores can be combined with features of their distribution, biochemical properties, and sequence entropy to improve inference. C A benchmark dataset of all non-redundant protein interactions with known interface structure.CC BY · PubMed Central
Fig. 2 Model predicts interacting residues and interacting proteins with high precision. A Recall on the held-out fraction of the positive benchmark set (x-axis) and false positive rate on the held-out dataset of non-interacting complexes (y-axis) at a score threshold that gives the corresponding recall. Our logistic regression model (purple) reduces the false positive rate compared to the previous EV Complex score (blue). B Prediction of protein interaction is based on the prediction of interacting residues. Number of predicted interacting residue pairs for complexes inferred to interact (purple) or not inferred to interact (gray) based on our stringent protein complex prediction threshold. C Example performance on known interaction between ABC transporter permease and ATP binding subunit (UniProt IDs: Y1470_HAEIN and Y1471_HAEIN, PDB ID: 2NQ2 chains C and A [https://www.rcsb.org/structure/2NQ2]. D Example performance on known interaction between DNA primase PriS and PriL (UniProt IDs: PRIS_SULSO and PRIL_SULSO, PDB ID: 1ZT2 chain A and B [https://www.rcsb.org/structure/1ZT2]).
Fig. 2 Model predicts interacting residues and interacting proteins with high precision. A Recall on the held-out fraction of the positive benchmark set (x-axis) and false positive rate on the held-out dataset of non-interacting complexes (y-axis) at a score threshold that gives the corresponding recall. Our logistic regression model (purple) reduces the false positive rate compared to the previous EV Complex score (blue). B Prediction of protein interaction is based on the prediction of interacting residues. Number of predicted interacting residue pairs for complexes inferred to interact (purple) or not inferred to interact (gray) based on our stringent protein complex prediction threshold. C Example performance on known interaction between ABC transporter permease and ATP binding subunit (UniProt IDs: Y1470_HAEIN and Y1471_HAEIN, PDB ID: 2NQ2 chains C and A [https://www.rcsb.org/structure/2NQ2]. D Example performance on known interaction between DNA primase PriS and PriL (UniProt IDs: PRIS_SULSO and PRIL_SULSO, PDB ID: 1ZT2 chain A and B [https://www.rcsb.org/structure/1ZT2]).CC BY · PubMed Central
Fig. 3 Discovery of hundreds of new interactions in the E. coli membrane proteome. A We searched a high-value subset of the 10∧7 possible interactions in the E. coli proteome by searching membrane compartments with themselves and with adjacent membrane compartments. B We found 504 high-scoring protein interactions in the cell envelope, including 75 with structural characterization and 186 with previous experimental evidence (and no structural characterization). C–E 3D configurations of previously structurally characterized interactions are accurately predicted by molecular docking with inferred restraints. RMSDs calculated by comparison to known structures with PDB IDs 3RKO [https://www.rcsb.org/structure/3RKO], 2WU2 [https://www.rcsb.org/structure/2WU2], and 2HQS [https://www.rcsb.org/structure/2HQS], respectively. F–H Example of three docked models of newly resolved protein complexes: BamE/MltB, YajC/FtsI, and Lnt/MurJ. Evolutionarily coupled residues used as restraints in docking are shown in magenta and connected with solid lines.
Fig. 3 Discovery of hundreds of new interactions in the E. coli membrane proteome. A We searched a high-value subset of the 10∧7 possible interactions in the E. coli proteome by searching membrane compartments with themselves and with adjacent membrane compartments. B We found 504 high-scoring protein interactions in the cell envelope, including 75 with structural characterization and 186 with previous experimental evidence (and no structural characterization). C–E 3D configurations of previously structurally characterized interactions are accurately predicted by molecular docking with inferred restraints. RMSDs calculated by comparison to known structures with PDB IDs 3RKO [https://www.rcsb.org/structure/3RKO], 2WU2 [https://www.rcsb.org/structure/2WU2], and 2HQS [https://www.rcsb.org/structure/2HQS], respectively. F–H Example of three docked models of newly resolved protein complexes: BamE/MltB, YajC/FtsI, and Lnt/MurJ. Evolutionarily coupled residues used as restraints in docking are shown in magenta and connected with solid lines.CC BY · PubMed Central
Fig. 4 Model of the flagellar hook-filament junction. A Schematic of the orientation of the bacterial flagellum. The proteins FlgL (green) and FlgK (blue) form two rings which create the junction between the hook and filament of the flagellum24. B Docked model of FlgL and FlgK using evolutionary couplings. PDB structures of homologous proteins from Salmonella typhimurium were used in docking (PDB IDs 2D4Y [https://www.rcsb.org/structure/2D4Y], and 2D4X [https://www.rcsb.org/structure/2D4X], respectively). Predicted interface residues are highlighted in purple. C A previously inferred model of the FlgK ring from Campylobacter jejuni24 was used to infer the structure of the entire hook-filament junction. Evolutionary coupled residues (purple) show the interface for FlgL ring insertion into the FlgK ring. By aligning our docked model to the C. jejuni ring, we show that an 11-mer ring of FlgL fits inside the FlgK structure.
Fig. 4 Model of the flagellar hook-filament junction. A Schematic of the orientation of the bacterial flagellum. The proteins FlgL (green) and FlgK (blue) form two rings which create the junction between the hook and filament of the flagellum24. B Docked model of FlgL and FlgK using evolutionary couplings. PDB structures of homologous proteins from Salmonella typhimurium were used in docking (PDB IDs 2D4Y [https://www.rcsb.org/structure/2D4Y], and 2D4X [https://www.rcsb.org/structure/2D4X], respectively). Predicted interface residues are highlighted in purple. C A previously inferred model of the FlgK ring from Campylobacter jejuni24 was used to infer the structure of the entire hook-filament junction. Evolutionary coupled residues (purple) show the interface for FlgL ring insertion into the FlgK ring. By aligning our docked model to the C. jejuni ring, we show that an 11-mer ring of FlgL fits inside the FlgK structure.CC BY · PubMed Central
Fig. 5 Atomic resolution model of the Tol-Pal complex. A Schematic of the proteins involved in the Tol-Pal complex2 (TolABR, CpoB, and Pal). Interactions with previously solved interfaces are shown in orange and interactions inferred by our method are shown in purple. B Complete model of the Tol-Pal complex inferred by aligning results of docked pairwise models. Note that CpoB is inferred to be a trimer in vivo but was docked as a monomer for modeling purposes C–F Residue resolution of TolB-Pal, TolB-TolA, TolR-TolB, and CpoB-TolB interfaces. The top 5 inferred interface contacts are shown in purple. Dashed lines indicate inferred contacts where one or more residues are missing from the solved structure. Structures used are 1TOL_A (TolA) [https://www.rcsb.org/structure/1TOL], 2HQS_A (TolB) [https://www.rcsb.org/structure/2HQS], 5BY4_A (TolR) [https://www.rcsb.org/structure/5BY4], 2HQS_H (Pal) [https://www.rcsb.org/structure/2HQS], and 2XDJ_A (CpoB) [https://www.rcsb.org/structure/2XDJ].
Fig. 5 Atomic resolution model of the Tol-Pal complex. A Schematic of the proteins involved in the Tol-Pal complex2 (TolABR, CpoB, and Pal). Interactions with previously solved interfaces are shown in orange and interactions inferred by our method are shown in purple. B Complete model of the Tol-Pal complex inferred by aligning results of docked pairwise models. Note that CpoB is inferred to be a trimer in vivo but was docked as a monomer for modeling purposes C–F Residue resolution of TolB-Pal, TolB-TolA, TolR-TolB, and CpoB-TolB interfaces. The top 5 inferred interface contacts are shown in purple. Dashed lines indicate inferred contacts where one or more residues are missing from the solved structure. Structures used are 1TOL_A (TolA) [https://www.rcsb.org/structure/1TOL], 2HQS_A (TolB) [https://www.rcsb.org/structure/2HQS], 5BY4_A (TolR) [https://www.rcsb.org/structure/5BY4], 2HQS_H (Pal) [https://www.rcsb.org/structure/2HQS], and 2XDJ_A (CpoB) [https://www.rcsb.org/structure/2XDJ].CC BY · PubMed Central
Fig. 6 Predicted protein interactions in eukaryotic proteomes. A Distribution of eukaryotic sequences in concatenated sequence alignments. Shown in orange are alignments that passed our sequence diversity threshold, and in gray are those that did not. B Number of correctly predicted inter-protein ECs for eukaryotic-exclusive complexes above the sequence diversity threshold. Eukaryotic-exclusive complexes are defined as complexes whose concatenated sequence alignment is at least 90% eukaryotic sequences. Inter protein ECs are defined as correct if their minimum atom distance is <8 Å. C The human spliceosome proteins Prp38 and MFAP1 have a known interface correctly predicted by our method, and have a protein interaction detected. The inter-protein ECs above the 80% precision threshold used throughout the paper are shown in purple. Known protein structure (PDB ID: 5F5U, [https://www.rcsb.org/structure/5F5U]) is used to visualize the subunits. D Schematic of predicted interaction between Lsm5 and Prp38 with top 5 inter-protein ECs.
Fig. 6 Predicted protein interactions in eukaryotic proteomes. A Distribution of eukaryotic sequences in concatenated sequence alignments. Shown in orange are alignments that passed our sequence diversity threshold, and in gray are those that did not. B Number of correctly predicted inter-protein ECs for eukaryotic-exclusive complexes above the sequence diversity threshold. Eukaryotic-exclusive complexes are defined as complexes whose concatenated sequence alignment is at least 90% eukaryotic sequences. Inter protein ECs are defined as correct if their minimum atom distance is <8 Å. C The human spliceosome proteins Prp38 and MFAP1 have a known interface correctly predicted by our method, and have a protein interaction detected. The inter-protein ECs above the 80% precision threshold used throughout the paper are shown in purple. Known protein structure (PDB ID: 5F5U, [https://www.rcsb.org/structure/5F5U]) is used to visualize the subunits. D Schematic of predicted interaction between Lsm5 and Prp38 with top 5 inter-protein ECs.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Research works printed but not linked: 1

  1. Marcotte, C. & Marcotte, E. Predicting functional linkages from gene fusion. Appl. Bioinf. 1, 1–8 (2002).

Cited by

Sorted by citations

See all 114

Metadata

Authors
Green, A. G., Elhabashy, H., Brock, K. P., Maddamsetti, R., Kohlbacher, O., & Marks, D. S.
Affiliations
  • Harvard University (US)
  • Max Planck Institute for Developmental Biology (DE)
  • Universitätsklinikum Tübingen (DE)
  • Broad Institute (US)
Corresponding author
  • Kohlbacher, Oliver
  • Marks, Debora S.
Venue
Nature Communications
Published
2021-03-02
Volume
12
Issue
1
Pages
1396-1396
ISSN
2041-1723 (ISSN-L)
OpenAlex
W3133793061
Cited by
114
References
57
Citations by year
2026: 14 · 2025: 17 · 2024: 24 · 2023: 25 · 2022: 21 · 2021: 13
FWCI
6.51
Type
journal article
Language
en
Open access
Locations
  • Nature Communications (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • RePEc: Research Papers in Economics (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
Funders
  • National Science Foundation (DGE1144152)
  • U.S. Department of Health & Human Services | National Institutes of Health (T15LM007092)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Protein Structure and Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Machine Learning in Bioinformatics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
Keywords
protein-protein interaction, membrane protein complexes
MeSH terms
Amino Acids, Bacterial Proteins, Base Sequence, Escherichia coli, Eukaryotic Cells, Membrane Proteins, Protein Binding, Genome, Bacterial, Evolution, Molecular, Proteome, Protein Interaction Mapping, Molecular Docking Simulation

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review