Ranking models of transmembrane β-barrel proteins using Z-coordinate predictions

Outer membrane β-barrel structure prediction through the lens of AlphaFold2

Most proteins found in the outer membrane of gram-negative bacteria share a common domain: the transmembrane β-barrel. These outer membrane β-barrels (OMBBs) occur in multiple sizes and different families with a wide range of functions evolved independently by amplification from a pool of homologous ancestral ββ-hairpins. This is part of the reason why predicting their three-dimensional (3D) structure, especially by homology modeling, is a major challenge. Recently, DeepMind's AlphaFold v2 (AF2) became the first structure prediction method to reach close-to-experimental atomic accuracy in CASP even for difficult targets. However, membrane proteins, especially OMBBs, were not abundant during their training, raising the question of how accurate the predictions are for these families. In this study, we assessed the performance of AF2 in the prediction of OMBBs and OMBB-like folds of various topologies using an in-house-developed tool for the analysis of OMBB 3D structures, and barrOs. In agreement with previous studies on other membrane protein classes, our results indicate that AF2 predicts transmembrane β-barrel structures at high accuracy independently of the use of templates, even for novel topologies absent from the training set. These results provide confidence on the models generated by AF2 and open the door to the structural elucidation of novel transmembrane β-barrel topologies identified in high-throughput OMBB annotation studies or designed de novo.

General features of transmembrane beta barrels from a large database

Large datasets contribute new insights to subjects formerly investigated by exemplars. We used coevolution data to create a large, high-quality database of transmembrane β-barrels (TMBB). By applying simple feature detection on generated evolutionary contact maps, our method (IsItABarrel) achieves 95.88% balanced accuracy when discriminating among protein classes. Moreover, comparison with IsItABarrel revealed a high rate of false positives in previous TMBB algorithms. In addition to being more accurate than previous datasets, our database (available online) contains 1,938,936 bacterial TMBB proteins from 38 phyla, respectively, 17 and 2.2 times larger than the previous sets TMBB-DB and OMPdb. We anticipate that due to its quality and size, the database will serve as a useful resource where high-quality TMBB sequence data are required. We found that TMBBs can be divided into 11 types, three of which have not been previously reported. We find tremendous variance in proteome percentage among TMBB-containing organisms with some using 6.79% of their proteome for TMBBs and others using as little as 0.27% of their proteome. The distribution of the lengths of the TMBBs is suggestive of previously hypothesized duplication events. In addition, we find that the C-terminal β-signal varies among different classes of bacteria though its consensus sequence is LGLGYRF. However, this β-signal is only characteristic of prototypical TMBBs. The ten non-prototypical barrel types have other C-terminal motifs, and it remains to be determined if these alternative motifs facilitate TMBB insertion or perform any other signaling function.

Computational structure prediction provides a plausible mechanism for electron transfer by the outer membrane protein Cyc2 from Acidithiobacillus ferrooxidans

Cyc2 is the key protein in the outer membrane of Acidithiobacillus ferrooxidans that mediates electron transfer between extracellular inorganic iron and the intracellular central metabolism. This cytochrome c is specific for iron and interacts with periplasmic proteins to complete a reversible electron transport chain. A structure of Cyc2 has not yet been characterized experimentally. Here we describe a structural model of Cyc2, and associated proteins, to highlight a plausible mechanism for the ferrous iron electron transfer chain. A comparative modeling protocol specific for trans membrane beta barrel (TMBB) proteins in acidophilic conditions (pH ~ 2) was applied to the primary sequence of Cyc2. The proposed structure has three main regimes: Extracellular loops exposed to low-pH conditions, a TMBB, and an N-terminal cytochrome-like region within the periplasmic space. The Cyc2 model was further refined by identifying likely iron and heme docking sites. This represents the first computational model of Cyc2 that accounts for the membrane microenvironment and the acidity in the extracellular matrix. This approach can be used to model other TMBBs which can be critical for chemolithotrophic microbial growth.

Enhancing the Prediction of Transmembrane β-Barrel Segments with Chain Learning and Feature Sparse Representation

Transmembrane β-barrels (TMBs) are one important class of membrane proteins that play crucial functions in the cell. Membrane proteins are difficult wet-lab targets of structural biology, which call for accurate computational prediction approaches. Here, we developed a novel method named MemBrain-TMB to predict the spanning segments of transmembrane β-barrel from amino acid sequence. MemBrain-TMB is a statistical machine learning-based model, which is constructed using a new chain learning algorithm with input features encoded by the image sparse representation approach. We considered the relative status information between neighboring residues for enhancing the performance, and the matrix of features was translated into feature image by sparse coding algorithm for noise and dimension reduction. To deal with the diverse loop length problem, we applied a dynamic threshold method, which is particularly useful for enhancing the recognition of short loops and tight turns. Our experiments demonstrate that the new protocol designed in MemBrain-TMB effectively helps improve prediction performance.

All-atom 3D structure prediction of transmembrane β-barrel proteins from sequences

Transmembrane β-barrels (TMBs) carry out major functions in substrate transport and protein biogenesis but experimental determination of their 3D structure is challenging. Encouraged by successful de novo 3D structure prediction of globular and α-helical membrane proteins from sequence alignments alone, we developed an approach to predict the 3D structure of TMBs. The approach combines the maximum-entropy evolutionary coupling method for predicting residue contacts (EVfold) with a machine-learning approach (boctopus2) for predicting β-strands in the barrel. In a blinded test for 19 TMB proteins of known structure that have a sufficient number of diverse homologous sequences available, this combined method (EVfold_bb) predicts hydrogen-bonded residue pairs between adjacent β-strands at an accuracy of ∼70%. This accuracy is sufficient for the generation of all-atom 3D models. In the transmembrane barrel region, the average 3D structure accuracy [template-modeling (TM) score] of top-ranked models is 0.54 (ranging from 0.36 to 0.85), with a higher (44%) number of residue pairs in correct strand-strand registration than in earlier methods (18%). Although the nonbarrel regions are predicted less accurately overall, the evolutionary couplings identify some highly constrained loop residues and, for FecA protein, the barrel including the structure of a plug domain can be accurately modeled (TM score = 0.68). Lower prediction accuracy tends to be associated with insufficient sequence information and we therefore expect increasing numbers of β-barrel families to become accessible to accurate 3D structure prediction as the number of available sequences increases.

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1-2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug-specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans-membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α-helical MPs as well as β-barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge-based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade.

Charge asymmetry in the proteins of the outer membrane

MOTIVATION

Outer membrane beta-barrels (OMBBs) are the proteins found in the outer membrane of bacteria, mitochondria and chloroplasts. There are thousands of beta-barrels reported in genomic databases with ∼2-3% of the genes in gram-negative bacteria encoding these proteins. These proteins have a wide variety of biological functions including active and passive transport, cell adhesion, catalysis and structural anchoring. Of the non-redundant OMBB structures in the Protein Data Bank, half have been solved during the past 5 years. This influx of information provides new opportunities for understanding the chemistry of these proteins. The distribution of charges in proteins in the outer membrane has implications for how the mechanism of outer membrane protein insertion is understood. Understanding the distribution of charges might also assist in organism selection for the heterologous expression of mitochondrial OMBBs.

RESULTS

We find a strong asymmetry in the charge distribution of these proteins. For the outward-facing residues of the beta-barrel within regions of similar amino acid density for both membrane leaflets, the external side of the outer membrane contains almost three times the number of charged residues as the internal side of the outer membrane. Moreover, the lipid bilayer of the outer membrane is asymmetric, and the overall preference for amino acid types to be in the external leaflet of the membrane correlates roughly with the hydrophobicity of the membrane lipids. This preference is demonstrably related to the difference in lipid composition of the external and internal leaflets of the membrane.