Sergey Ovchinnikov

Robust and accurate prediction of residue-residue interactions across protein interfaces using evolutionary information.

Do the amino acid sequence identities of residues that make contact across protein interfaces covary during evolution? If so, such covariance could be used to predict contacts across interfaces and assemble models of biological complexes. We find that residue pairs identified using a pseudo-likelihood-based method to covary across protein–protein interfaces in the 50S ribosomal unit and 28 additional bacterial protein complexes with known structure are almost always in contact in the complex, provided that the number of aligned sequences is greater than the average length of the two proteins. We use this method to make subunit contact predictions for an additional 36 protein complexes with unknown structures, and present models based on these predictions for the tripartite ATP-independent periplasmic (TRAP) transporter, the tripartite efflux system, the pyruvate formate lyase-activating enzyme complex, and the methionine ABC transporter. DOI: http://dx.doi.org/10.7554/eLife.02030.001

Assessing the utility of coevolution-based residue–residue contact predictions in a sequence- and structure-rich era

Significance We develop an improved method for predicting residue–residue contacts in protein structures that achieves higher accuracy than previous methods by integrating structural context and sequence coevolution information. We then determine the conditions under which these predicted contacts are likely to be useful for structure modeling and identify more than 400 protein families where these conditions are currently met. Recently developed methods have shown considerable promise in predicting residue–residue contacts in protein 3D structures using evolutionary covariance information. However, these methods require large numbers of evolutionarily related sequences to robustly assess the extent of residue covariation, and the larger the protein family, the more likely that contact information is unnecessary because a reasonable model can be built based on the structure of a homolog. Here we describe a method that integrates sequence coevolution and structural context information using a pseudolikelihood approach, allowing more accurate contact predictions from fewer homologous sequences. We rigorously assess the utility of predicted contacts for protein structure prediction using large and representative sequence and structure databases from recent structure prediction experiments. We find that contact predictions are likely to be accurate when the number of aligned sequences (with sequence redundancy reduced to 90%) is greater than five times the length of the protein, and that accurate predictions are likely to be useful for structure modeling if the aligned sequences are more similar to the protein of interest than to the closest homolog of known structure. These conditions are currently met by 422 of the protein families collected in the Pfam database.

Protein structure determination using metagenome sequence data

Filling in the protein fold picture Fewer than a third of the 14,849 known protein families have at least one member with an experimentally determined structure. This leaves more than 5000 protein families with no structural information. Protein modeling using residue-residue contacts inferred from evolutionary data has been successful in modeling unknown structures, but it requires large numbers of aligned sequences. Ovchinnikov et al. augmented such sequence alignments with metagenome sequence data (see the Perspective by Söding). They determined the number of sequences required to allow modeling, developed criteria for model quality, and, where possible, improved modeling by matching predicted contacts to known structures. Their method predicted quality structural models for 614 protein families, of which about 140 represent newly discovered protein folds. Science, this issue p. 294; see also p. 248 Combining metagenome data with protein structure prediction generates models for 614 families with unknown structures. Despite decades of work by structural biologists, there are still ~5200 protein families with unknown structure outside the range of comparative modeling. We show that Rosetta structure prediction guided by residue-residue contacts inferred from evolutionary information can accurately model proteins that belong to large families and that metagenome sequence data more than triple the number of protein families with sufficient sequences for accurate modeling. We then integrate metagenome data, contact-based structure matching, and Rosetta structure calculations to generate models for 614 protein families with currently unknown structures; 206 are membrane proteins and 137 have folds not represented in the Protein Data Bank. This approach provides the representative models for large protein families originally envisioned as the goal of the Protein Structure Initiative at a fraction of the cost.

Large-scale determination of previously unsolved protein structures using evolutionary information

The prediction of the structures of proteins without detectable sequence similarity to any protein of known structure remains an outstanding scientific challenge. Here we report significant progress in this area. We first describe de novo blind structure predictions of unprecendented accuracy we made for two proteins in large families in the recent CASP11 blind test of protein structure prediction methods by incorporating residue–residue co-evolution information in the Rosetta structure prediction program. We then describe the use of this method to generate structure models for 58 of the 121 large protein families in prokaryotes for which three-dimensional structures are not available. These models, which are posted online for public access, provide structural information for the over 400,000 proteins belonging to the 58 families and suggest hypotheses about mechanism for the subset for which the function is known, and hypotheses about function for the remainder. DOI: http://dx.doi.org/10.7554/eLife.09248.001

Improved de novo structure prediction in CASP11 by incorporating coevolution information into Rosetta.

We describe CASP11 de novo blind structure predictions made using the Rosetta structure prediction methodology with both automatic and human assisted protocols. Model accuracy was generally improved using coevolution derived residue–residue contact information as restraints during Rosetta conformational sampling and refinement, particularly when the number of sequences in the family was more than three times the length of the protein. The highlight was the human assisted prediction of T0806, a large and topologically complex target with no homologs of known structure, which had unprecedented accuracy—<3.0 Å root‐mean‐square deviation (RMSD) from the crystal structure over 223 residues. For this target, we increased the amount of conformational sampling over our fully automated method by employing an iterative hybridization protocol. Our results clearly demonstrate, in a blind prediction scenario, that coevolution derived contacts can considerably increase the accuracy of template‐free structure modeling. Proteins 2016; 84(Suppl 1):67–75.

Architectures of Lipid Transport Systems for the Bacterial Outer Membrane.

How phospholipids are trafficked between the bacterial inner and outer membranes through the hydrophilic space of the periplasm is not known. We report that members of the mammalian cell entry (MCE) protein family form hexameric assemblies with a central channel capable of mediating lipid transport. The E.xa0coli MCE protein, MlaD, forms a ring associated with an ABC transporter complex in the inner membrane. A soluble lipid-binding protein, MlaC, ferries lipids between MlaD and an outer membrane proteinxa0complex. In contrast, EM structures of two other E.xa0coli MCE proteins show that YebT forms an elongated tube consisting of seven stacked MCE rings, and PqiB adopts a syringe-like architecture. Both YebT and PqiB create channels of sufficient length to span the periplasmic space. This work reveals diverse architectures of highly conserved protein-based channels implicated in the transport of lipids between the membranes of bacteria and some eukaryotic organelles.

Cryo-EM structure of the protein-conducting ERAD channel Hrd1 in complex with Hrd3

Misfolded endoplasmic reticulum proteins are retro-translocated through the membrane into the cytosol, where they are poly-ubiquitinated, extracted from the membrane, and degraded by the proteasome—a pathway termed endoplasmic reticulum-associated protein degradation (ERAD). Proteins with misfolded domains in the endoplasmic reticulum lumen or membrane are discarded through the ERAD-L and ERAD-M pathways, respectively. In Saccharomyces cerevisiae, both pathways require the ubiquitin ligase Hrd1, a multi-spanning membrane protein with a cytosolic RING finger domain. Hrd1 is the crucial membrane component for retro-translocation, but it is unclear whether it forms a protein-conducting channel. Here we present a cryo-electron microscopy structure of S. cerevisiae Hrd1 in complex with its endoplasmic reticulum luminal binding partner, Hrd3. Hrd1 forms a dimer within the membrane with one or two Hrd3 molecules associated at its luminal side. Each Hrd1 molecule has eight transmembrane segments, five of which form an aqueous cavity extending from the cytosol almost to the endoplasmic reticulum lumen, while a segment of the neighbouring Hrd1 molecule forms a lateral seal. The aqueous cavity and lateral gate are reminiscent of features of protein-conducting conduits that facilitate polypeptide movement in the opposite direction—from the cytosol into or across membranes. Our results suggest that Hrd1 forms a retro-translocation channel for the movement of misfolded polypeptides through the endoplasmic reticulum membrane.

Computation and Functional Studies Provide a Model for the Structure of the Zinc Transporter hZIP4

Background: ZIP transporters increase the cytosolic concentration of first row transition metals. Results: We have developed a structural model of hZIP4 by combining protein prediction methods with in situ experiments. Conclusion: Analysis of our experiments provides insight into the permeation pathway of hZIP4. Significance: Comparison of this model to membrane transporter crystal structures provides a structural linkage to MFS proteins. Members of the Zrt and Irt protein (ZIP) family are a central participant in transition metal homeostasis as they function to increase the cytosolic concentration of zinc and/or iron. However, the lack of a crystal structure hinders elucidation of the molecular mechanism of ZIP proteins. Here, we employed GREMLIN, a co-evolution-based contact prediction approach in conjunction with the Rosetta structure prediction program to construct a structural model of the human (h) ZIP4 transporter. The predicted contact data are best fit by modeling hZIP4 as a dimer. Mutagenesis of residues that comprise a central putative hZIP4 transmembrane transition metal coordination site in the structural model alter the kinetics and specificity of hZIP4. Comparison of the hZIP4 dimer model to all known membrane protein structures identifies the 12-transmembrane monomeric Piriformospora indica phosphate transporter (PiPT), a member of the major facilitator superfamily (MFS), as a likely structural homolog.

Structure of a bd oxidase indicates similar mechanisms for membrane-integrated oxygen reductases.

Peering into a membrance oxidase Microorganisms have evolved a number of enzymes to reduce oxygen and prevent oxidative stress. Cytochrome bd oxidases serve this role and also protect pathogenic bacteria from nitric acid; however, this class of enzymes so far has eluded high-resolution crystallography. Safarian et al. were able to resolve the three-dimensional structure of cytochrome bd oxidase from a thermophilic bacterium (see the Perspective by Cook and Poole). The overall structure and triangular arrangement of its heme cofactors bear little structural resemblance to those of other membrane-spanning oxidases, despite serving a similar function. Science, this issue p. 583; see also p. 518 Structural details reveal how pathogens protect against oxidative stress and nitric oxide. The cytochrome bd oxidases are terminal oxidases that are present in bacteria and archaea. They reduce molecular oxygen (dioxygen) to water, avoiding the production of reactive oxygen species. In addition to their contribution to the proton motive force, they mediate viability under oxygen-related stress conditions and confer tolerance to nitric oxide, thus contributing to the virulence of pathogenic bacteria. Here we present the atomic structure of the bd oxidase from Geobacillus thermodenitrificans, revealing a pseudosymmetrical subunit fold. The arrangement and order of the heme cofactors support the conclusions from spectroscopic measurements that the cleavage of the dioxygen bond may be mechanistically similar to that in the heme-copper–containing oxidases, even though the structures are completely different.

Origins of coevolution between residues distant in protein 3D structures

Significance Coevolution-derived contact predictions are enabling accurate protein structure modeling. However, coevolving residues are not always in contact, and this is a potential source of error in such modeling efforts. To investigate the sources of such errors and, more generally, the origins of coevolution in protein structures, we provide a global overview of the contributions to the “exceptions” to the general rule that coevolving residues are close in protein three-dimensional structures. Residue pairs that directly coevolve in protein families are generally close in protein 3D structures. Here we study the exceptions to this general trend—directly coevolving residue pairs that are distant in protein structures—to determine the origins of evolutionary pressure on spatially distant residues and to understand the sources of error in contact-based structure prediction. Over a set of 4,000 protein families, we find that 25% of directly coevolving residue pairs are separated by more than 5 Å in protein structures and 3% by more than 15 Å. The majority (91%) of directly coevolving residue pairs in the 5–15 Å range are found to be in contact in at least one homologous structure—these exceptions arise from structural variation in the family in the region containing the residues. Thirty-five percent of the exceptions greater than 15 Å are at homo-oligomeric interfaces, 19% arise from family structural variation, and 27% are in repeat proteins likely reflecting alignment errors. Of the remaining long-range exceptions (<1% of the total number of coupled pairs), many can be attributed to close interactions in an oligomeric state. Overall, the results suggest that directly coevolving residue pairs not in repeat proteins are spatially proximal in at least one biologically relevant protein conformation within the family; we find little evidence for direct coupling between residues at spatially separated allosteric and functional sites or for increased direct coupling between residue pairs on putative allosteric pathways connecting them.