Membrane protein-regulated networks across human cancers.

Alterations in membrane proteins (MPs) and their regulated pathways have been established as cancer hallmarks and extensively targeted in clinical applications. However, the analysis of MP-interacting proteins and downstream pathways across human malignancies remains challenging. Here, we present a systematically integrated method to generate a resource of cancer membrane protein-regulated networks (CaMPNets), containing 63,746 high-confidence protein-protein interactions (PPIs) for 1962 MPs, using expression profiles from 5922 tumors with overall survival outcomes across 15 human cancers. Comprehensive analysis of CaMPNets links MP partner communities and regulated pathways to provide MP-based gene sets for identifying prognostic biomarkers and druggable targets. For example, we identify CHRNA9 with 12 PPIs (e.g., ERBB2) can be a therapeutic target and find its anti-metastasis agent, bupropion, for treatment in nicotine-induced breast cancer. This resource is a study to systematically integrate MP interactions, genomics, and clinical outcomes for helping illuminate cancer-wide atlas and prognostic landscapes in tumor homo/heterogeneity.

Relationships between residue Voronoi volume and sequence conservation in proteins.

BACKGROUND: Functional and biophysical constraints can cause different levels of sequence conservation in proteins. Previously, structural properties, e.g., relative solvent accessibility (RSA) and packing density of the weighted contact number (WCN), have been found to be related to protein sequence conservation (CS). The Voronoi volume has recently been recognized as a new structural property of the local protein structural environment reflecting CS. However, for surface residues, it is sensitive to water molecules surrounding the protein structure. Herein, we present a simple structural determinant termed the relative space of Voronoi volume (RSV); it uses the Voronoi volume and the van der Waals volume of particular residues to quantify the local structural environment. METHODS: RSV (range, 0-1) is defined as (Voronoi volume-van der Waals volume)/Voronoi volume of the target residue. The concept of RSV describes the extent of available space for every protein residue. RESULTS: RSV and Voronoi profiles with and without water molecules (RSVw, RSV, VOw, and VO) were compared for 554 non-homologous proteins. RSV (without water) showed better Pearson's correlations with CS than did RSVw, VO, or VOw values. The mean correlation coefficient between RSV and CS was 0.51, which is comparable to the correlation between RSA and CS (0.49) and that between WCN and CS (0.56). CONCLUSIONS: RSV is a robust structural descriptor with and without water molecules and can quantitatively reflect evolutionary information in a single protein structure. Therefore, it may represent a practical structural determinant to study protein sequence, structure, and function relationships.

On the relationship between residue structural environment and sequence conservation in proteins.

Residues that are crucial to protein function or structure are usually evolutionarily conserved. To identify the important residues in protein, sequence conservation is estimated, and current methods rely upon the unbiased collection of homologous sequences. Surprisingly, our previous studies have shown that the sequence conservation is closely correlated with the weighted contact number (WCN), a measure of packing density for residue's structural environment, calculated only based on the Cα positions of a protein structure. Moreover, studies have shown that sequence conservation is correlated with environment-related structural properties calculated based on different protein substructures, such as a protein's all atoms, backbone atoms, side-chain atoms, or side-chain centroid. To know whether the Cα atomic positions are adequate to show the relationship between residue environment and sequence conservation or not, here we compared Cα atoms with other substructures in their contributions to the sequence conservation. Our results show that Cα positions are substantially equivalent to the other substructures in calculations of various measures of residue environment. As a result, the overlapping contributions between Cα atoms and the other substructures are high, yielding similar structure-conservation relationship. Take the WCN as an example, the average overlapping contribution to sequence conservation is 87% between Cα and all-atom substructures. These results indicate that only Cα atoms of a protein structure could reflect sequence conservation at the residue level.

Ubiquitin-coated nanodiamonds bind to autophagy receptors for entry into the selective autophagy pathway.

Selective macroautophagy/autophagy plays a pivotal role in the processing of foreign pathogens and cellular components to maintain homeostasis in human cells. To date, numerous studies have demonstrated the uptake of nanoparticles by cells, but their intracellular processing through selective autophagy remains unclear. Here we show that carbon-based nanodiamonds (NDs) coated with ubiquitin (Ub) bind to autophagy receptors (SQSTM1 [sequestosome 1], OPTN [optineurin], and CALCOCO2/NDP52 [calcium binding and coiled-coil domain 2]) and are then linked to MAP1LC3/LC3 (microtubule-associated protein 1 light chain 3) for entry into the selective autophagy pathway. NDs are ultimately delivered to lysosomes. Ectopically expressed SQSTM1-green fluorescence protein (GFP) could bind to the Ub-coated NDs. By contrast, the Ub-associated domain mutant of SQSTM1 (ΔUBA)-GFP did not bind to the Ub-coated NDs. Chloroquine, an autophagy inhibitor, prevented the ND-containing autophagosomes from fusing with lysosomes. Furthermore, autophagy receptors OPTN and CALCOCO2/NDP52, involved in the processing of bacteria, were found to be involved in the selective autophagy of NDs. However, ND particles located in the lysosomes of cells did not induce mitotic blockage, senescence, or cell death. Single ND clusters in the lysosomes of cells were observed in the xenografted human lung tumors of nude mice. This study demonstrated for the first time that Ub-coated nanoparticles bind to autophagy receptors for entry into the selective autophagy pathway, facilitating their delivery to lysosomes.

MIB: Metal Ion-Binding Site Prediction and Docking Server.

The structure of a protein determines its biological function(s) and its interactions with other factors; the binding regions tend to be conserved in sequence and structure, and the interacting residues involved are usually in close 3D space. The Protein Data Bank currently contains more than 110 000 protein structures, approximately one-third of which contain metal ions. Identifying and characterizing metal ion-binding sites is thus essential for investigating a protein's function(s) and interactions. However, experimental approaches are time-consuming and costly. The web server reported here was built to predict metal ion-binding residues and to generate the predicted metal ion-bound 3D structure. Binding templates have been constructed for regions that bind 12 types of metal ion-binding residues have been used to construct binding templates. The templates include residues within 3.5 Šof the metal ion, and the fragment transformation method was used for structural comparison between query proteins and templates without any data training. Through the adjustment of scoring functions, which are based on the similarity of structure and binding residues. Twelve kinds of metal ions (Ca2+, Cu2+, Fe3+, Mg2+, Mn2+, Zn2+, Cd2+, Fe2+, Ni2+, Hg2+, Co2+, and Cu+) binding residues prediction are supported. MIB also provides the metal ions docking after prediction. The MIB server is available at http://bioinfo.cmu.edu.tw/MIB/ .

Sequence Conservation, Radial Distance and Packing Density in Spherical Viral Capsids.

The conservation level of a residue is a useful measure about the importance of that residue in protein structure and function. Much information about sequence conservation comes from aligning homologous sequences. Profiles showing the variation of the conservation level along the sequence are usually interpreted in evolutionary terms and dictated by site similarities of a proper set of homologous sequences. Here, we report that, of the viral icosahedral capsids, the sequence conservation profile can be determined by variations in the distances between residues and the centroid of the capsid - with a direct inverse proportionality between the conservation level and the centroid distance - as well as by the spatial variations in local packing density. Examining both the centroid and the packing density models against a dataset of 51 crystal structures of nonhomologous icosahedral capsids, we found that many global patterns and minor features derived from the viral structures are consistent with those present in the sequence conservation profiles. The quantitative link between the level of conservation and structural features like centroid-distance or packing density allows us to look at residue conservation from a structural viewpoint as well as from an evolutionary viewpoint.

(PS)2: protein structure prediction server version 3.0.

Protein complexes are involved in many biological processes. Examining coupling between subunits of a complex would be useful to understand the molecular basis of protein function. Here, our updated (PS)(2) web server predicts the three-dimensional structures of protein complexes based on comparative modeling; furthermore, this server examines the coupling between subunits of the predicted complex by combining structural and evolutionary considerations. The predicted complex structure could be indicated and visualized by Java-based 3D graphics viewers and the structural and evolutionary profiles are shown and compared chain-by-chain. For each subunit, considerations with or without the packing contribution of other subunits cause the differences in similarities between structural and evolutionary profiles, and these differences imply which form, complex or monomeric, is preferred in the biological condition for the subunit. We believe that the (PS)(2) server would be a useful tool for biologists who are interested not only in the structures of protein complexes but also in the coupling between subunits of the complexes. The (PS)(2) is freely available at http://ps2v3.life.nctu.edu.tw/.

The conservation profile of a protein bears the imprint of the molecule that is evolutionarily coupled to the protein.

The conservation profile of a protein is a curve of the conservation levels of amino acids along the sequence. Biologists are usually more interested in individual points on the curve (namely, the conserved amino acids) than the overall shape of the curve. Here, we show that the conservation curves of proteins bear the imprints of molecules that are evolutionarily coupled to the proteins. Our method is based on recent studies that a sequence conservation profile is quantitatively linked to its structural packing profile. We find that the conservation profiles of nucleic acid (NA) binding proteins are better correlated with the packing profiles of the protein-NA complexes than those of the proteins alone. This indicates that a nucleic acid binding protein evolves to accommodate the nucleic acid in such a way that the residues involved in binding have their conservation levels closely coupled with the specific nucleotides.

On the packing density of the unbound protein-protein interaction interface and its implications in dynamics.

BACKGROUND: Characterizing the interface residues will help shed light on protein-protein interactions, which are involved in many important biological processes. Many studies focus on characterizing sequence or structure features of protein interfaces, but there are few studies characterizing the dynamics of interfaces. Therefore, we would like to know whether there is any specific dynamics pattern in the protein-protein interaction interfaces. Thermal fluctuation is an important dynamical property for a residue, and could be quickly estimated by local packing density without large computation since studies have showen closely relationship between these two properties. Therefore, we divided surface of an unbound subunit (free protein subunits before they are involved in forming the protein complexes) into several separate regions, and compared their average thermal fluctuations of different regions in order to characterize the dynamics pattern in unbound protein-protein interaction interfaces. RESULTS: We used weighted contact numbers (WCN), a parameter-free method to quantify packing density, to estimate the thermal fluctuations of residues in the interfaces. By analyzing the WCN distributions of interfaces in unbound subunits from 1394 non-homologous protein complexes, we show that the residues in the central regions of interfaces have higher packing density (i.e. more rigid); on the other hand, residues surrounding the central regions have smaller packing density (i.e. more flexible). The distinct distributions of packing density, suggesting distinct thermal fluctuation, reveals specific dynamics pattern in the interface of unbound protein subunits. CONCLUSIONS: We found general trend that the unbound protein-protein interaction interfaces consist of rigid residues in the central regions, which are surrounded by flexible residues. This finding suggests that the dynamics might be one of the important features for the formation of protein complexes.

Protein thermal stability enhancement by designing salt bridges: a combined computational and experimental study.

Protein thermal stability is an important factor considered in medical and industrial applications. Many structural characteristics related to protein thermal stability have been elucidated, and increasing salt bridges is considered as one of the most efficient strategies to increase protein thermal stability. However, the accurate simulation of salt bridges remains difficult. In this study, a novel method for salt-bridge design was proposed based on the statistical analysis of 10,556 surface salt bridges on 6,493 X-ray protein structures. These salt bridges were first categorized based on pairing residues, secondary structure locations, and Cα-Cα distances. Pairing preferences generalized from statistical analysis were used to construct a salt-bridge pair index and utilized in a weighted electrostatic attraction model to find the effective pairings for designing salt bridges. The model was also coupled with B-factor, weighted contact number, relative solvent accessibility, and conservation prescreening to determine the residues appropriate for the thermal adaptive design of salt bridges. According to our method, eight putative salt-bridges were designed on a mesophilic ß-glucosidase and 24 variants were constructed to verify the predictions. Six putative salt-bridges leaded to the increase of the enzyme thermal stability. A significant increase in melting temperature of 8.8, 4.8, 3.7, 1.3, 1.2, and 0.7°C of the putative salt-bridges N437K-D49, E96R-D28, E96K-D28, S440K-E70, T231K-D388, and Q277E-D282 was detected, respectively. Reversing the polarity of T231K-D388 to T231D-D388K resulted in a further increase in melting temperatures by 3.6°C, which may be caused by the transformation of an intra-subunit electrostatic interaction into an inter-subunit one depending on the local environment. The combination of the thermostable variants (N437K, E96R, T231D and D388K) generated a melting temperature increase of 15.7°C. Thus, this study demonstrated a novel method for the thermal adaptive design of salt bridges through inference of suitable positions and substitutions.

Local packing density is the main structural determinant of the rate of protein sequence evolution at site level.

Functional and biophysical constraints result in site-dependent patterns of protein sequence variability. It is commonly assumed that the key structural determinant of site-specific rates of evolution is the Relative Solvent Accessibility (RSA). However, a recent study found that amino acid substitution rates correlate better with two Local Packing Density (LPD) measures, the Weighted Contact Number (WCN) and the Contact Number (CN), than with RSA. This work aims at a more thorough assessment. To this end, in addition to substitution rates, we considered four other sequence variability scores, four measures of solvent accessibility (SA), and other CN measures. We compared all properties for each protein of a structurally and functionally diverse representative dataset of monomeric enzymes. We show that the best sequence variability measures take into account phylogenetic tree topology. More importantly, we show that both LPD measures (WCN and CN) correlate better than all of the SA measures, regardless of the sequence variability score used. Moreover, the independent contribution of the best LPD measure is approximately four times larger than that of the best SA measure. This study strongly supports the conclusion that a site's packing density rather than its solvent accessibility is the main structural determinant of its rate of evolution.

CELLO2GO: a web server for protein subCELlular LOcalization prediction with functional gene ontology annotation.

CELLO2GO (http://cello.life.nctu.edu.tw/cello2go/) is a publicly available, web-based system for screening various properties of a targeted protein and its subcellular localization. Herein, we describe how this platform is used to obtain a brief or detailed gene ontology (GO)-type categories, including subcellular localization(s), for the queried proteins by combining the CELLO localization-predicting and BLAST homology-searching approaches. Given a query protein sequence, CELLO2GO uses BLAST to search for homologous sequences that are GO annotated in an in-house database derived from the UniProt KnowledgeBase database. At the same time, CELLO attempts predict at least one subcellular localization on the basis of the species in which the protein is found. When homologs for the query sequence have been identified, the number of terms found for each of their GO categories, i.e., cellular compartment, molecular function, and biological process, are summed and presented as pie charts representing possible functional annotations for the queried protein. Although the experimental subcellular localization of a protein may not be known, and thus not annotated, CELLO can confidentially suggest a subcellular localization. CELLO2GO should be a useful tool for research involving complex subcellular systems because it combines CELLO and BLAST into one platform and its output is easily manipulated such that the user-specific questions may be readily addressed.

A mechanistic stress model of protein evolution accounts for site-specific evolutionary rates and their relationship with packing density and flexibility.

BACKGROUND: Protein sites evolve at different rates due to functional and biophysical constraints. It is usually considered that the main structural determinant of a site's rate of evolution is its Relative Solvent Accessibility (RSA). However, a recent comparative study has shown that the main structural determinant is the site's Local Packing Density (LPD). LPD is related with dynamical flexibility, which has also been shown to correlate with sequence variability. Our purpose is to investigate the mechanism that connects a site's LPD with its rate of evolution. RESULTS: We consider two models: an empirical Flexibility Model and a mechanistic Stress Model. The Flexibility Model postulates a linear increase of site-specific rate of evolution with dynamical flexibility. The Stress Model, introduced here, models mutations as random perturbations of the protein's potential energy landscape, for which we use simple Elastic Network Models (ENMs). To account for natural selection we assume a single active conformation and use basic statistical physics to derive a linear relationship between site-specific evolutionary rates and the local stress of the mutant's active conformation.We compare both models on a large and diverse dataset of enzymes. In a protein-by-protein study we found that the Stress Model outperforms the Flexibility Model for most proteins. Pooling all proteins together we show that the Stress Model is strongly supported by the total weight of evidence. Moreover, it accounts for the observed nonlinear dependence of sequence variability on flexibility. Finally, when mutational stress is controlled for, there is very little remaining correlation between sequence variability and dynamical flexibility. CONCLUSIONS: We developed a mechanistic Stress Model of evolution according to which the rate of evolution of a site is predicted to depend linearly on the local mutational stress of the active conformation. Such local stress is proportional to LPD, so that this model explains the relationship between LPD and evolutionary rate. Moreover, the model also accounts for the nonlinear dependence between evolutionary rate and dynamical flexibility.

Site-specific structural constraints on protein sequence evolutionary divergence: local packing density versus solvent exposure.

Protein sequences evolve under selection pressures imposed by functional and biophysical requirements, resulting in site-dependent rates of amino acid substitution. Relative solvent accessibility (RSA) and local packing density (LPD) have emerged as the best candidates to quantify structural constraint. Recent research assumes that RSA is the main determinant of sequence divergence. However, it is not yet clear which is the best predictor of substitution rates. To address this issue, we compared RSA and LPD with site-specific rates of evolution for a diverse data set of enzymes. In contrast with recent studies, we found that LPD measures correlate better than RSA with evolutionary rate. Moreover, the independent contribution of RSA is minor. Taking into account that LPD is related to backbone flexibility, we put forward the possibility that the rate of evolution of a site is determined by the ease with which the backbone deforms to accommodate mutations.

Functional investigation of transmembrane helix 3 in H⁺-translocating pyrophosphatase.

H⁺-translocating pyrophosphatase (H⁺-PPase, EC 3.6.1.1) plays an important role in acidifying vacuoles by transporting protons across membranes at the expense of pyrophosphate (PP(i)) hydrolysis. Vigna radiata H⁺-PPase (VrH⁺-PPase) contains 16 transmembrane helices (TMs). The hydrophobicity of TM3 is relatively lower than that of most other TMs, and the amino acids in this TM are highly conserved in plants. Furthermore, TM5 and -6, which are the core TMs involving in H⁺-PPase functions, are near TM3. It is thus proposed that TM3 is associated with H⁺-PPase activity. To address this possibility, site-directed mutagenesis was applied in this investigation to determine the role of TM3 in VrH⁺-PPase. Upon alanine/serine substitution, T138 and S142, whose side chains face toward the center TMs, were found to be involved in efficient proton transport. G149/S153 and G160/A164 pairs at the crucial termini of the two GxxxG-like motifs are indispensable in maintaining enzymatic activities and conformational stability. Moreover, stability in the vicinity surrounding G149 is pivotal for efficient expression. S153, M161 and A164 are critical for the K⁺-mediated stimulation of H⁺-PPase. Taken together, our results demonstrate that TM3 plays essential roles in PP(i) hydrolysis, proton transport, expression, and K⁺ stimulation of H⁺-PPase.

Computational analysis of KRAS mutations: implications for different effects on the KRAS p.G12D and p.G13D mutations.

BACKGROUND: The issue of whether patients diagnosed with metastatic colorectal cancer who harbor KRAS codon 13 mutations could benefit from the addition of anti-epidermal growth factor receptor therapy remains under debate. The aim of the current study was to perform computational analysis to investigate the structural implications of the underlying mutations caused by c.38G>A (p.G13D) on protein conformation. METHODS: Molecular dynamics (MD) simulations were performed to understand the plausible structural and dynamical implications caused by c.35G>A (p.G12D) and c.38G>A (p.G13D). The potential of mean force (PMF) simulations were carried out to determine the free energy profiles of the binding processes of GTP interacting with wild-type (WT) KRAS and its mutants (MT). RESULTS: Using MD simulations, we observed that the root mean square deviation (RMSD) increased as a function of time for the MT c.35G>A (p.G12D) and MT c.38G>A (p.G13D) when compared with the WT. We also observed that the GTP-binding pocket in the c.35G>A (p.G12D) mutant is more open than that of the WT and the c.38G>A (p.G13D) proteins. Intriguingly, the analysis of atomic fluctuations and free energy profiles revealed that the mutation of c.35G>A (p.G12D) may induce additional fluctuations in the sensitive sites (P-loop, switch I and II regions). Such fluctuations may promote instability in these protein regions and hamper GTP binding. CONCLUSIONS: Taken together with the results obtained from MD and PMF simulations, the present findings implicate fluctuations at the sensitive sites (P-loop, switch I and II regions). Our findings revealed that KRAS mutations in codon 13 have similar behavior as KRAS WT. To gain a better insight into why patients with metastatic colorectal cancer (mCRC) and the KRAS c.38G>A (p.G13D) mutation appear to benefit from anti-EGFR therapy, the role of the KRAS c.38G>A (p.G13D) mutation in mCRC needs to be further investigated.

On the relationship between the sequence conservation and the packing density profiles of the protein complexes.

We have recently showed that the weighted contact number profiles (or the packing density profiles) of proteins are well correlated with those of the corresponding sequence conservation profiles. The results suggest that a protein structure may contain sufficient information about sequence conservation comparable to that derived from multiple homologous sequences. However, there are ambiguities concerning how to compute the packing density of the subunit of a protein complex. For the subunits of a complex, there are different ways to compute its packing density--one including the packing contributions of the other subunits and the other one excluding their contributions. Here we selected two sets of enzyme complexes. Set A contains complexes with the active sites comprising residues from multiple subunits, while set B contains those with the active sites residing on single subunits. In Set A, if the packing density profile of a subunit is computed considering the contributions of the other subunits of the complex, it will agree better with the sequence conservation profile. But in Set B the situations are reversed. The results may be due to the stronger functional and structural constraints on the evolution processes on the complexes of Set A than those of Set B to maintain the enzymatic functions of the complexes. The comparison of the packing density and the sequence conservation profiles may provide a simple yet potentially useful way to understanding the structural and evolutionary couplings between the subunits of protein complexes.

Deriving correlated motions in proteins from X-ray structure refinement by using TLS parameters.

Dynamic information in proteins may provide valuable information for understanding allosteric regulation of protein complexes or long-range effects of the mutations on enzyme activity. Experimental data such as X-ray B-factors or NMR order parameters provide a convenient estimate of atomic fluctuations (or atomic auto-correlated motions) in proteins. However, it is not as straightforward to obtain atomic cross-correlated motions in proteins - one usually resorts to more sophisticated computational methods such as Molecular Dynamics, normal mode analysis or atomic network models. In this report, we show that atomic cross-correlations can be reliably obtained directly from protein structure using X-ray refinement data. We have derived an analytic form of atomic correlated motions in terms of the original TLS parameters used to refine the B-factors of X-ray structures. The correlated maps computed using this equation are well correlated with those of the method based on a mechanical model (the correlation coefficient is 0.75) for a non-homologous dataset comprising 100 structures. We have developed an approach to compute atomic cross-correlations directly from X-ray protein structure. Being in analytic form, it is fast and provides a feasible way to compute correlated motions in proteins in a high throughput way. In addition, avoiding sophisticated computational operations; it provides a quick, reliable way, especially for non-computational biologists, to obtain dynamics information directly from protein structure relevant to its function.

CPred: a web server for predicting viable circular permutations in proteins.

Circular permutation (CP) is a protein structural rearrangement phenomenon, through which nature allows structural homologs to have different locations of termini and thus varied activities, stabilities and functional properties. It can be applied in many fields of protein research and bioengineering. The limitation of applying CP lies in its technical complexity, high cost and uncertainty of the viability of the resulting protein variants. Not every position in a protein can be used to create a viable circular permutant, but there is still a lack of practical computational tools for evaluating the positional feasibility of CP before costly experiments are carried out. We have previously designed a comprehensive method for predicting viable CP cleavage sites in proteins. In this work, we implement that method into an efficient and user-friendly web server named CPred (CP site predictor), which is supposed to be helpful to promote fundamental researches and biotechnological applications of CP. The CPred is accessible at http://sarst.life.nthu.edu.tw/CPred.

pKNOT v.2: the protein KNOT web server.

Knotted proteins have recently received lots of attention due to their interesting topological novelty as well as its puzzling folding mechanisms. We previously published a pKNOT server, which provides a structural database of knotted proteins, analysis tools for detecting and analyzing knotted regions from structures as well as a Java-based 3D graphics viewer for visualizing knotted structures. However, there lacks a convenient platform performing similar tasks directly from 'protein sequences'. In the current version of the web server, referred to as pKNOT v.2, we implement a homology modeling tool such that the server can now accept protein sequences in addition to 3D structures or Protein Data Bank (PDB) IDs and return knot analysis. In addition, we have updated the database of knotted proteins from the current PDB with a combination of automatic and manual procedure. We believe that the updated pKNOT server with its extended functionalities will provide better service to biologists interested in the research of knotted proteins. The pKNOT v.2 is available from http://pknot.life.nctu.edu.tw/.