Conserved network properties of helical membrane protein structures and its implication for improving membrane protein homology modeling at the twilight zone

GRAFENE: Graphlet-based alignment-free network approach integrates 3D structural and sequence (residue order) data to improve protein structural comparison

Initial protein structural comparisons were sequence-based. Since amino acids that are distant in the sequence can be close in the 3-dimensional (3D) structure, 3D contact approaches can complement sequence approaches. Traditional 3D contact approaches study 3D structures directly and are alignment-based. Instead, 3D structures can be modeled as protein structure networks (PSNs). Then, network approaches can compare proteins by comparing their PSNs. These can be alignment-based or alignment-free. We focus on the latter. Existing network alignment-free approaches have drawbacks: 1) They rely on naive measures of network topology. 2) They are not robust to PSN size. They cannot integrate 3) multiple PSN measures or 4) PSN data with sequence data, although this could improve comparison because the different data types capture complementary aspects of the protein structure. We address this by: 1) exploiting well-established graphlet measures via a new network alignment-free approach, 2) introducing normalized graphlet measures to remove the bias of PSN size, 3) allowing for integrating multiple PSN measures, and 4) using ordered graphlets to combine the complementary PSN data and sequence (specifically, residue order) data. We compare synthetic networks and real-world PSNs more accurately and faster than existing network (alignment-free and alignment-based), 3D contact, or sequence approaches.

Residue centrality in alpha helical polytopic transmembrane protein structures

Transmembrane proteins serve as receptors, transporters or as enzymes. They mediate a broad range of fundamental cellular activities including signal transduction, cell trafficking and photosynthesis. In this study, we analyzed the significance of central residues in the polytopic transmembrane proteins. Each protein is represented as an undirected graph, where residues represent nodes and inter-residue interactions as the edges. Residue centrality was calculated by removing the nodes and its corresponding edges from the protein contact network. Results revealed that 80% of the predicted central residues had normalized conservation values below the mean since they were slowly evolving conserved sites. We also found that 56% of amino acids were interacting with the ligand molecules and metal ions. Predicted central residues in the polytopic transmembrane proteins were found to account for 84% of binding and active site amino acids. From mutation sensitivity analysis, it was observed that 89% of central residues had deleterious mutations whose probabilities were greater than their mean value. Interestingly, we find that z-score values of each amino acid positively correlate with the conservation scores and also with the degrees of each node. Results show that 87% of central residues are hub residues.

HYDROGEN BOND NETWORK OF 1-ALKYL-3-METHYLIMIDAZOLIUM IONIC LIQUIDS: A NETWORK THEORY ANALYSIS

Hydrogen bond is a key factor in the determination of structures and properties of room-temperature ionic liquids. Connections of these hydrogen bonds form a network. In this work, we analyzed the hydrogen bond network of 1-alkyl-3-methylimidazolium ionic liquids using network theory. A two-dimensional view of the hydrogen bond network has been generated, the connection pattern shown that the average length of line shape connection is 2.44 to 2.77 for six 1-alkyl-3-methylimidazolium ionic liquids, and the connection patterns are different for short and long alkyl side chain length. The degree of each ion was calculated and analyzed. The nodes with zero degree were adopted to detect the boundary of the clusters in the ionic liquids, which have no hydrogen bond connected with neighbor ions. This work indicates that the network analysis method is useful for understanding and predicting the structure and function of RTILs.

Developing a high-quality scoring function for membrane protein structures based on specific inter-residue interactions

Membrane proteins are of particular biological and pharmaceutical importance, and computational modeling and structure prediction approaches play an important role in studies of membrane proteins. Developing an accurate model quality assessment program is of significance to the structure prediction of membrane proteins. Few such programs are proposed that can be applied to a broad range of membrane protein classes and perform with high accuracy. We developed a new model scoring function Interaction-based Quality assessment (IQ), based on the analysis of four types of inter-residue interactions within the transmembrane domains of helical membrane proteins. This function was tested using three high-quality model sets: all 206 models of GPCR Dock 2008, all 284 models of GPCR Dock 2010, and all 92 helical membrane protein models of the HOMEP set. For all three sets, the scoring function can select the native structures among all of the models with the success rates of 93, 85, and 100% respectively. For comparison, these three model sets were also adopted for a recently published model assessment program for membrane protein structures, ProQM, which gave the success rates of 85, 79, and 92% separately. These results suggested that IQ outperforms ProQM when only the transmembrane regions of the models are considered. This scoring function should be useful for the computational modeling of membrane proteins.