Hammad Naveed

Lipid-binding surfaces of membrane proteins: evidence from evolutionary and structural analysis.

Membrane proteins function in the diverse environment of the lipid bilayer. Experimental evidence suggests that some lipid molecules bind tightly to specific sites on the membrane protein surface. These lipid molecules often act as co-factors and play important functional roles. In this study, we have assessed the evolutionary selection pressure experienced at lipid-binding sites in a set of α-helical and β-barrel membrane proteins using posterior probability analysis of the ratio of synonymous vs. nonsynonymous substitutions (ω-ratio). We have also carried out a geometric analysis of the membrane protein structures to identify residues in close contact with co-crystallized lipids. We found that residues forming cholesterol-binding sites in both β(2)-adrenergic receptor and Na(+)-K(+)-ATPase exhibit strong conservation, which can be characterized by an expanded cholesterol consensus motif for GPCRs. Our results suggest the functional importance of aromatic stacking interactions and interhelical hydrogen bonds in facilitating protein-cholesterol interactions, which is now reflected in the expanded motif. We also find that residues forming the cardiolipin-binding site in formate dehydrogenase-N γ-subunit and the phosphatidylglycerol binding site in KcsA are under strong purifying selection pressure. Although the lipopolysaccharide (LPS)-binding site in ferric hydroxamate uptake receptor (FhuA) is only weakly conserved, we show using a statistical mechanical model that LPS binds to the least stable FhuA β-strand and protects it from the bulk lipid. Our results suggest that specific lipid binding may be a general mechanism employed by β-barrel membrane proteins to stabilize weakly stable regions. Overall, we find that the residues forming specific lipid binding sites on the surfaces of membrane proteins often experience strong purifying selection pressure.

Predicting weakly stable regions, oligomerization state, and protein-protein interfaces in transmembrane domains of outer membrane proteins

Although the structures of many β-barrel membrane proteins are available, our knowledge of the principles that govern their energetics and oligomerization states is incomplete. Here we describe a computational method to study the transmembrane (TM) domains of β-barrel membrane proteins. Our method is based on a physical interaction model, a simplified conformational space for efficient enumeration, and an empirical potential function from a detailed combinatorial analysis. Using this method, we can identify weakly stable regions in the TM domain, which are found to be important structural determinants for β-barrel membrane proteins. By calculating the melting temperatures of the TM strands, our method can also assess the stability of β-barrel membrane proteins. Predictions on membrane enzyme PagP are consistent with recent experimental NMR and mutant studies. We have also discovered that out-clamps, in-plugs, and oligomerization are 3 general mechanisms for stabilizing weakly stable TM regions. In addition, we have found that extended and contiguous weakly stable regions often signal the existence of an oligomer and that strands located in the interfaces of protein–protein interactions are considerably less stable. Based on these observations, we can predict oligomerization states and can identify the interfaces of protein–protein interactions for β-barrel membrane proteins by using either structure or sequence information. In a set of 25 nonhomologous proteins with known structures, our method successfully predicted whether a protein forms a monomer or an oligomer with 91% accuracy; in addition, our method identified with 82% accuracy the protein–protein interaction interfaces by using sequence information only when correct strands are given.

Structure-based Analysis of VDAC1 Protein: DEFINING OLIGOMER CONTACT SITES*

Background: VDAC1 was shown to undergo oligomeric assembly, an event that is coupled to apoptosis induction. Results: Structure- and computation-based predications of VDAC1 oligomerization sites were confirmed. Conclusion: Upon apoptosis induction, VDAC1 undergoes conformational changes and oligomerization. Significance: Dissection of VDAC1 dimerization/oligomerization provides structural insight into the oligomeric status of cellular VDAC1 under physiological and apoptotic conditions. The outer mitochondrial membrane protein, the voltage-dependent anion channel (VDAC), is increasingly implicated in the control of apoptosis. Oligomeric assembly of VDAC1 was shown to be coupled to apoptosis induction, with oligomerization increasing substantially upon apoptosis induction and inhibited by apoptosis blockers. In this study, structure- and computation-based selection of the predicated VDAC1 dimerization site, in combination with site-directed mutagenesis, cysteine replacement, and chemical cross-linking, were employed to identify contact sites between VDAC1 molecules in dimers and higher oligomers. The predicted weakly stable β-strands were experimentally found to represent the interfaces between VDAC1 monomers composing the oligomer. Replacing hydrophobic amino acids with charged residues in β-strands 1, 2, and 19 interfered with VDAC1 oligomerization. The proximity of β-strands 1, 2, and 19 within the VDAC1 dimer and the existence of other association sites involving β-strand 16 were confirmed when a cysteine was introduced at defined positions in cysteineless VDAC1 mutants, together with the use of cysteine-specific cross-linker bis(maleimido)ethane. Moreover, the results suggest that VDAC1 also exists as a dimer that upon apoptosis induction undergoes conformational changes and that its oligomerization proceeds through a series of interactions involving two distinct interfaces. Dissection of VDAC1 dimerization/oligomerization as presented here provides structural insight into the oligomeric status of cellular VDAC1 under physiological and apoptotic conditions.

Predicting three-dimensional structures of transmembrane domains of β-barrel membrane proteins

β-Barrel membrane proteins are found in the outer membrane of gram-negative bacteria, mitochondria, and chloroplasts. They are important for pore formation, membrane anchoring, and enzyme activity. These proteins are also often responsible for bacterial virulence. Due to difficulties in experimental structure determination, they are sparsely represented in the protein structure databank. We have developed a computational method for predicting structures of the transmembrane (TM) domains of β-barrel membrane proteins. Based on physical principles, our method can predict structures of the TM domain of β-barrel membrane proteins of novel topology, including those from eukaryotic mitochondria. Our method is based on a model of physical interactions, a discrete conformational state space, an empirical potential function, as well as a model to account for interstrand loop entropy. We are able to construct three-dimensional atomic structure of the TM domains from sequences for a set of 23 nonhomologous proteins (resolution 1.8-3.0 Å). The median rmsd of TM domains containing 75-222 residues between predicted and measured structures is 3.9 Å for main chain atoms. In addition, stability determinants and protein-protein interaction sites can be predicted. Such predictions on eukaryotic mitochondria outer membrane protein Tom40 and VDAC are confirmed by independent mutagenesis and chemical cross-linking studies. These results suggest that our model captures key components of the organization principles of β-barrel membrane protein assembly.

Computational studies of membrane proteins: models and predictions for biological understanding.

We discuss recent progresses in computational studies of membrane proteins based on physical models with parameters derived from bioinformatics analysis. We describe computational identification of membrane proteins and prediction of their topology from sequence, discovery of sequence and spatial motifs, and implications of these discoveries. The detection of evolutionary signal for understanding the substitution pattern of residues in the TM segments and for sequence alignment is also discussed. We further discuss empirical potential functions for energetics of inserting residues in the TM domain, for interactions between TM helices or strands, and their applications in predicting lipid-facing surfaces of the TM domain. Recent progresses in structure predictions of membrane proteins are also reviewed, with further discussions on calculation of ensemble properties such as melting temperature based on simplified state space model. Additional topics include prediction of oligomerization state of membrane proteins, identification of the interfaces for protein-protein interactions, and design of membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

CoCiter: An Efficient Tool to Infer Gene Function by Assessing the Significance of Literature Co-Citation

A routine approach to inferring functions for a gene set is by using function enrichment analysis based on GO, KEGG or other curated terms and pathways. However, such analysis requires the existence of overlapping genes between the query gene set and those annotated by GO/KEGG. Furthermore, GO/KEGG databases only maintain a very restricted vocabulary. Here, we have developed a tool called “CoCiter” based on literature co-citations to address the limitations in conventional function enrichment analysis. Co-citation analysis is widely used in ranking articles and predicting protein-protein interactions (PPIs). Our algorithm can further assess the co-citation significance of a gene set with any other user-defined gene sets, or with free terms. We show that compared with the traditional approaches, CoCiter is a more accurate and flexible function enrichment analysis method. CoCiter is freely available at www.picb.ac.cn/hanlab/cociter/.

Engineered oligomerization state of OmpF protein through computational design decouples oligomer dissociation from unfolding

Biogenesis of β-barrel membrane proteins is a complex, multistep, and as yet incompletely characterized process. The bacterial porin family is perhaps the best-studied protein family among β-barrel membrane proteins that allows diffusion of small solutes across the bacterial outer membrane. In this study, we have identified residues that contribute significantly to the protein-protein interaction (PPI) interface between the chains of outer membrane protein F (OmpF), a trimeric porin, using an empirical energy function in conjunction with an evolutionary analysis. By replacing these residues through site-directed mutagenesis either with energetically favorable residues or substitutions that do not occur in natural bacterial outer membrane proteins, we succeeded in engineering OmpF mutants with dimeric and monomeric oligomerization states instead of a trimeric oligomerization state. Moreover, our results suggest that the oligomerization of OmpF proceeds through a series of interactions involving two distinct regions of the extensive PPI interface: two monomers interact to form a dimer through the PPI interface near G19. This dimer then interacts with another monomer through the PPI interface near G135 to form a trimer. We have found that perturbing the PPI interface near G19 results in the formation of the monomeric OmpF only. Thermal denaturation of the designed dimeric OmpF mutant suggests that oligomer dissociation can be separated from the process of protein unfolding. Furthermore, the conserved site near G57 and G59 is important for the PPI interface and might provide the essential scaffold for PPIs.

Mechanisms of regulating cell topology in proliferating epithelia: impact of division plane, mechanical forces, and cell memory.

Regulation of cell growth and cell division has a fundamental role in tissue formation, organ development, and cancer progression. Remarkable similarities in the topological distributions were found in a variety of proliferating epithelia in both animals and plants. At the same time, there are species with significantly varied frequency of hexagonal cells. Moreover, local topology has been shown to be disturbed on the boundary between proliferating and quiescent cells, where cells have fewer sides than natural proliferating epithelia. The mechanisms of regulating these topological changes remain poorly understood. In this study, we use a mechanical model to examine the effects of orientation of division plane, differential proliferation, and mechanical forces on animal epithelial cells. We find that regardless of orientation of division plane, our model can reproduce the commonly observed topological distributions of cells in natural proliferating animal epithelia with the consideration of cell rearrangements. In addition, with different schemes of division plane, we are able to generate different frequency of hexagonal cells, which is consistent with experimental observations. In proliferating cells interfacing quiescent cells, our results show that differential proliferation alone is insufficient to reproduce the local changes in cell topology. Rather, increased tension on the boundary, in conjunction with differential proliferation, can reproduce the observed topological changes. We conclude that both division plane orientation and mechanical forces play important roles in cell topology in animal proliferating epithelia. Moreover, cell memory is also essential for generating specific topological distributions.

Outer Membrane Protein Folding and Topology from a Computational Transfer Free Energy Scale

Knowledge of the transfer free energy of amino acids from aqueous solution to a lipid bilayer is essential for understanding membrane protein folding and for predicting membrane protein structure. Here we report a computational approach that can calculate the folding free energy of the transmembrane region of outer membrane β-barrel proteins (OMPs) by combining an empirical energy function with a reduced discrete state space model. We quantitatively analyzed the transfer free energies of 20 amino acid residues at the center of the lipid bilayer of OmpLA. Our results are in excellent agreement with the experimentally derived hydrophobicity scales. We further exhaustively calculated the transfer free energies of 20 amino acids at all positions in the TM region of OmpLA. We found that the asymmetry of the Gram-negative bacterial outer membrane as well as the TM residues of an OMP determine its functional fold in vivo. Our results suggest that the folding process of an OMP is driven by the lipid-facing residues in its hydrophobic core, and its NC-IN topology is determined by the differential stabilities of OMPs in the asymmetrical outer membrane. The folding free energy is further reduced by lipid A and assisted by general depth-dependent cooperativities that exist between polar and ionizable residues. Moreover, context-dependency of transfer free energies at specific positions in OmpLA predict regions important for protein function as well as structural anomalies. Our computational approach is fast, efficient and applicable to any OMP.

Geometric order in proliferating epithelia: Impact of rearrangements and cleavage plane orientation

Regulation of proliferation is required for normal development of tissues and prevention of cancer formation. Continuous control of proliferation leads to regular shaped cells forming characteristic tissue patterns. Epithelial tissues serve as a model system for studying tissue morphogenesis. Several groups have studied epithelial morphogenesis using topological or geometric models, with various assumptions. In this study, we have developed a method to simulate the dynamic process of proliferating epithelia using an off-lattice cellular model. Our method realistically models the shape, size, geometry, lineage, cleavage plane orientation as well as topological properties of individual cells. We find that cellular rearrangements and cleavage plane orientation are critical in the formation of the observed cellular pattern of epithelia, including a high percentage of hexagons in proliferating epithelial cells. It is likely that the rearrangements and orientation of the cleavage plane reduces the overall stress on the cell. We show that a high percentage of hexagons in proliferating epithelia can be obtained using uniform growth rates, which was considered unlikely in previous studies. Our off-lattice cellular model provides an improvement over existing topological for studying the dynamics of proliferating epithelium.