PHYSICO: An UNIX based Standalone Procedure for Computation of Individual and Group Properties of Protein Sequences

In the genomic and proteomic era, efficient and automated analyses of sequence properties of protein have become an important task in bioinformatics. There are general public licensed (GPL) software tools to perform a part of the job. However, computations of mean properties of large number of orthologous sequences are not possible from the above mentioned GPL sets. Further, there is no GPL software or server which can calculate window dependent sequence properties for a large number of sequences in a single run. With a view to overcome above limitations, we have developed a standalone procedure i.e. PHYSICO, which performs various stages of computation in a single run based on the type of input provided either in RAW-FASTA or BLOCK-FASTA format and makes excel output for: a) Composition, Class composition, Mean molecular weight, Isoelectic point, Aliphatic index and GRAVY, b) column based compositions, variability and difference matrix, c) 25 kinds of window dependent sequence properties. The program is fast, efficient, error free and user friendly. Calculation of mean and standard deviation of homologous sequences sets, for comparison purpose when relevant, is another attribute of the program; a property seldom seen in existing GPL softwares.

Availability

PHYSICO is freely available for non-commercial/academic user in formal request to the corresponding author ni.ca.vinurub.hcetoib@eejrenabka     

Salt-Bridge Energetics in Halophilic Proteins

Halophilic proteins have greater abundance of acidic over basic and very low bulky hydrophobic residues. Classical electrostatic stabilization was suggested as the key determinant for halophilic adaptation of protein. However, contribution of specific electrostatic interactions (i.e. salt-bridges) to overall stability of halophilic proteins is yet to be understood. To understand this, we use Adaptive-Poison-Boltzmann-Solver Methods along with our home-built automation to workout net as well as associated component energy terms such as desolvation energy, bridge energy and background energy for 275 salt-bridges from 20 extremely halophilic proteins. We then perform extensive statistical analysis on general and energetic attributes on these salt-bridges. On average, 8 salt-bridges per 150 residues protein were observed which is almost twice than earlier report. Overall contributions of salt-bridges are −3.0 kcal mol⁻¹. Majority (78%) of salt-bridges in our dataset are stable and conserved in nature. Although, average contributions of component energy terms are equal, their individual details vary greatly from one another indicating their sensitivity to local micro-environment. Notably, 35% of salt-bridges in our database are buried and stable. Greater desolvation penalty of these buried salt-bridges are counteracted by stable network salt-bridges apart from favorable equal contributions of bridge and background terms. Recruitment of extensive network salt-bridges (46%) with a net contribution of −5.0 kcal mol⁻¹ per salt-bridge, seems to be a halophilic design wherein favorable average contribution of background term (−10 kcal mol⁻¹) exceeds than that of bridge term (−7 kcal mol⁻¹). Interiors of proteins from halophiles are seen to possess relatively higher abundance of charge and polar side chains than that of mesophiles which seems to be satisfied by cooperative network salt-bridges. Overall, our theoretical analyses provide insight into halophilic signature in its specific electrostatic interactions which we hope would help in protein engineering and bioinformatics studies.     

Understanding antibody–antigen associations by molecular dynamics simulations: Detection of important intra- and inter-molecular salt bridges

1 NSec molecular dynamics (MD) simulation of anti-hen egg white antibody, HyHEL63 (HH63), complexed with HEL reveals important molecular interactions, not revealed in its X-ray crystal structure. These molecular interactions were predicted to be critical for the complex formation, based on structure–function studies of this complex and 3-other anti-HEL antibodies, HH8, HH10 and HH26, HEL complexes. All four antibodies belong to the same structural family, referred to here as HH10 family. Ala scanning results show that they recognize ‘coincident epitopes’. 1 NSec explicit, with periodic boundary condition, MD simulation of HH63-HEL reveals the presence of functionally important salt-bridges. Around 200 ps in vacuo and an additional 20 ps explicit simulation agree with the observations from 1 Nsec simulation. Intra-molecular salt-bridges predicted to play significant roles in the complex formation, were revealed during MD simulation. A very stabilizing salt-bridge network, and another intra-molecular salt-bridge, at the binding site of HEL, revealed during the MD simulation, is proposed to predipose binding site geometry for specific binding. All the revealed salt-bridges are present in one or more of the other three complexes and/or involve “hot-spot” epitope and paratope residues. Most of these charged epitope residues make large contribution to the binding free energy. The “hot spot” epitope residue Lys97_Y, which significantly contributes to the free energy of binding in all the complexes, forms an intermolecular salt-bridge in several MD conformers. Our earlier computations have shown that this inter-molecular salt-bridge plays a significant role in determining specificity and flexibility of binding in the HH8-HEL and HH26-HEL complexes. Using a robust criterion of salt-bridge detection, this inter-molecular salt-bridge was detected in the native structures of the HH8-HEL and HH26-HEL complexes, but was not revealed in the crystal structure of HH63-HEL complex. The electrostatic strength of this revealed salt-bridge was very strong. During 1 Nsec MD simulation this salt-bridge networks with another inter-molecular salt-bridge to form an inter-molecular salt-bridge triad. Participation of Lys97_Y in the formation of inter-molecular triad further validates the functional importance of Lys97_Y in HH63-HEL associations. These results demonstrate that many important structural details of biomolecular interactions can be better understood when studied in a dynamic environment, and that MD simulations can complement and expand information obtained from static X-ray structure. This study also highlights “hot-spot” molecular interactions in HyHEL63-HEL complex. The publisher or recipient acknowledges right of the U.S. Government to retain a non exclusive, royalty-free license in and to any copyright covering the article.     

Salt-bridges in the microenvironment of stable protein structures

Salt-bridges (sb) play an important role in the folding and stability of proteins. This is deduced from the evaluation of net energy in the microenvironments (ME, residues that are 4 Å away from positive and negative partners of salt-bridge and interact with them). MEs act as a determinant of net-energy due to the intrinsic features in the sequence. The stability of extremophilic proteins is due to the presence of favorable residues at the ME without any unfavorable residues. We studied a dataset of four structures from the protein data bank (PDB) and a homology model (1HM5) to gain insights on this issue. Data shows that the presence of isolated charges and polar residues in the core of extremophilic proteins helps in the formation of stable salt-bridges with reduced desolvation. Thus, site-specific mutations with favorable residues at the ME will help to develop thermo stable proteins with strong salt bridges.     

Structural and mechanistic exploration of acid resistance: kinetic stability facilitates evolution of extremophilic behavior

Kinetically stable proteins are unique in that their stability is determined solely by kinetic barriers rather than by thermodynamic equilibria. To better understand how kinetic stability promotes protein survival under extreme environmental conditions, we analyzed the unfolding behavior and determined the structure of Nocardiopsis alba Protease A (NAPase), an acid-resistant, kinetically stable protease, and compared these results with a neutrophilic homolog, α-lytic protease (αLP). Although NAPase and αLP have the same number of acid-titratable residues, kinetic studies revealed that the height of the unfolding free energy barrier for NAPase is less sensitive to acid than that of αLP, thereby accounting for NAPase's improved tolerance of low pH. A comparison of the αLP and NAPase structures identified multiple salt-bridges in the domain interface of αLP that were relocated to outer regions of NAPase, suggesting a novel mechanism of acid stability in which acid-sensitive electrostatic interactions are rearranged to similarly affect the energetics of both the native state and the unfolding transition state. An acid-stable variant of αLP in which a single interdomain salt-bridge is replaced with a corresponding intradomain NAPase salt-bridge shows a dramatic > 15-fold increase in acid resistance, providing further evidence for this hypothesis. These observations also led to a general model of the unfolding transition state structure for αLP protease family members in which the two domains separate from each other while remaining relatively intact themselves. These results illustrate the remarkable utility of kinetic stability as an evolutionary tool for developing longevity over a broad range of harsh conditions.     

Comparative studies on ion-pair energetic,distribution among three domains of life: Archaea,eubacteria, and eukarya

Salt-bridges play a unique role in the structural and functional stability of proteins, especially under harsh environments. How these salt-bridges contribute to the overall thermodynamic stability of protein structure and function across different domains of life is elusive still date. To address the issue, statistical analyses on the energies of salt-bridges, involved in proteins' structure and function, are performed across three domains of life, that is, archaea, eubacteria, and eukarya. Results show that although the majority of salt-bridges are stable and conserved, yet the stability of archaeal proteins (∆∆G_net = −5.06 ± 3.8) is much more than that of eubacteria (∆∆G_net = −3.7 ± 2.9) and eukarya (∆∆G_net = −3.54 ± 3.1). Unlike earlier study with archaea, in eukarya and eubacteria, not all buried salt-bridge in our dataset are stable. Buried salt-bridges play surprising role in protein stability, whose variations are clearly observed among these domains. Greater desolvation penalty of buried salt-bridges is compensated by stable network of salt-bridges apart from equal contribution of bridge and background energy terms. On the basis proteins' secondary structure, topology, and evolution, our observation shows that salt-bridges when present closer to each other in sequence tend to form a greater number. Overall, our comparative study provides insight into the role of specific electrostatic interactions in proteins from different domains of life, which we hope, would be useful for protein engineering and bioinformatics study.     

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.     

Use of models in understanding the mechanism of action of haemoglobin

Models with three, four and eight salt-bridges have been used to study the mechanism of action of haemoglobin. Both side chains forming a salt-bridge, i.e. the proton acceptor and the proton donor, are postulated to change pK on ligation of oxygen. The eight salt-bridge model is able to predict, as a unified theory, both the degree of oxygenation and the Bohr effect at any PH and p_O2 value; this has not been done by any other published model. The predicted pK values for the Borh groups corresponde well with those measured experiemntally. This model predicts the pK values of those side chains responsible for the acid Bohr effect, suggesting that these correspond to the proton acceptors of the salt-bridges. The model also fulfils the condition of linearity between the fractional degree of oxygenation and fractional number of protons released. It is postulated that there is a gradual change in structure on going from deoxy to oxyhaemoglobin, due to the rupture of salt-bridges. The path folowed during this process will be both pH and p_O2 dependent. A formula describing the number of intact or broken salt-bridges as a function of pH and p_O2 was developed. This formula shows that the fractional number of broken salt-bridges reaches a minimum value of 0.2 at around pH 6.3 in the absence of oxygen. However, if oxygen is added, this fractional number approaches 1.0 soon after the partial pressure of oxygen goes above 40 mm Hg.     

Structure and Binding Interface of the Cytosolic Tails of αXβ2 Integrin

Background

Integrins are signal transducer proteins involved in a number of vital physiological processes including cell adhesion, proliferation and migration. Integrin molecules are hetero-dimers composed of two distinct subunits, α and β. In humans, 18 α and 8 β subunits are combined into 24 different integrin molecules. Each of the subunit comprises a large extracellular domain, a single pass transmembrane segment and a cytosolic tail (CT). The CTs of integrins are vital for bidirectional signal transduction and in maintaining the resting state of the receptors. A large number of intracellular proteins have been found to interact with the CTs of integrins linking integrins to the cytoskeleton.

Methodology/Principal Findings

In this work, we have investigated structure and interactions of CTs of the leukocyte specific integrin αXβ2. We determined the atomic resolution structure of a myristoylated CT of αX in perdeuterated dodecylphosphocholine (DPC) by NMR spectroscopy. Our results reveal that the 35-residue long CT of αX adopts an α-helical conformation for residues F4-N17 at the N-terminal region. The remaining residues located at the C-terminal segment of αX delineate a long loop of irregular conformations. A segment of the loop maintains packing interactions with the helical structure by an extended non-polar surface of the αX CT. Interactions between αX and β2 CTs are demonstrated by ¹⁵N-¹H HSQC NMR experiments. We find that residues constituting the polar face of the helical conformation of αX are involved in interactions with the N-terminal residues of β2 CT. A docked structure of the CT complex indicates that a network of polar and/or salt-bridge interactions may sustain the heteromeric interactions.

Conclusions/Significance

The current study provides important insights into the conservation of interactions and structures among different CTs of integrins.     

The regulatory alpha subunit of phosphorylase kinase may directly participate in the binding of glycogen phosphorylase

The yeast two-hybrid screen has been used to identify potential regions of interaction of the largest regulatory subunit, , of phosphorylase kinase (PhK) with two fragments of its protein substrate, glycogen phosphorylase b (Phb). One fragment, corresponding to residues 17-484 (PhbN"), contained the regulatory domain of the protein, but in missing the first 16 residues was devoid of the sole phosphorylation site of Phb, Ser14; the second fragment corresponded to residues 485-843 (PhbC) and contained the catalytic domain of Phb. Truncation fragments of the subunit were screened for interactions against these two substrate fragments. PhbC was not found to interact with any constructs; however, PhbN" interacted with a region of (residues 864-1014) that is near the phosphorylatable region of that subunit. PhbN" was also screened for interactions against a variety of fragments of the catalytic subunit of PhK; however, no interactions were detected, even with fulllength . Our results support the idea that amino acid residues proximal to the convertible serine of Phb are important for its specific interaction with the catalytic subunit of PhK, but that regions distinct from the convertible serine residue of Phb and from the catalytic domain of PhK may also be involved in the interaction of these two proteins.     

Identification of the amino acid residues involved in the hemolytic activity of the Cucumaria echinata lectin CEL-III

Background

CEL-III is a hemolytic lectin isolated from the sea cucumber Cucumaria echinata that shows Ca^2 +-dependent Gal/GalNAc-binding specificity. This lectin is composed of two carbohydrate-recognition domains (domains 1 and 2) and an oligomerization domain (domain 3) that facilitates CEL-III assembly in the target cell membrane to form ion-permeable pores.

Methods

Several amino acid residues in domain 3 were replaced by alanine, and hemolytic activity of the mutants was examined.

Results

K344A, K351A, K405A, K420A and K425A showed marked increases in activity. In particular, K405A had activity that was 360-fold higher than the wild-type recombinant CEL-III and 3.6-fold higher than the native protein purified from sea cucumber. Since these residues appear to play roles in the stabilization of domain 3 through ionic and hydrogen bonding interactions with other residues, the mutations of these residues presumably lead to destabilization of domain 3, which consequently induces the oligomerization of the protein through association of domain 3 in the membrane. In contrast, K338A, R378A and R408A mutants exhibited a marked decrease in hemolytic activity. Since these residues are located on the surface of domain 3 without significant interactions with other residue, they may be involved in the interaction with components of the target cell membrane.

Conclusions

Several amino acid residues, especially basic residues, are found to be involved in the hemolytic activity as well as the oligomerization ability of CEL-III.

General significance

The results provide important clues to the membrane pore-forming mechanism of CEL-III, which is also related to that of bacterial pore-forming toxins.     

A homologous mapping method for three-dimensional reconstruction of protein networks reveals disease-associated mutations

Background

One of the crucial steps toward understanding the associations among molecular interactions, pathways, and diseases in a cell is to investigate detailed atomic protein-protein interactions (PPIs) in the structural interactome. Despite the availability of large-scale methods for analyzing PPI networks, these methods often focused on PPI networks using genome-scale data and/or known experimental PPIs. However, these methods are unable to provide structurally resolved interaction residues and their conservations in PPI networks.

Results

Here, we reconstructed a human three-dimensional (3D) structural PPI network (hDiSNet) with the detailed atomic binding models and disease-associated mutations by enhancing our PPI families and 3D–domain interologs from 60,618 structural complexes and complete genome database with 6,352,363 protein sequences across 2274 species. hDiSNet is a scale-free network (γ?=?2.05), which consists of 5177 proteins and 19,239 PPIs with 5843 mutations. These 19,239 structurally resolved PPIs not only expanded the number of PPIs compared to present structural PPI network, but also achieved higher agreement with gene ontology similarities and higher co-expression correlation than the ones of 181,868 experimental PPIs recorded in public databases. Among 5843 mutations, 1653 and 790 mutations involved in interacting domains and contacting residues, respectively, are highly related to diseases. Our hDiSNet can provide detailed atomic interactions of human disease and their associated proteins with mutations. Our results show that the disease-related mutations are often located at the contacting residues forming the hydrogen bonds or conserved in the PPI family. In addition, hDiSNet provides the insights of the FGFR (EGFR)-MAPK pathway for interpreting the mechanisms of breast cancer and ErbB signaling pathway in brain cancer.

Conclusions

Our results demonstrate that hDiSNet can explore structural-based interactions insights for understanding the mechanisms of disease-associated proteins and their mutations. We believe that our method is useful to reconstruct structurally resolved PPI networks for interpreting structural genomics and disease associations.

     

Stabilizing salt-bridge enhances protein thermostability by reducing the heat capacity change of unfolding

Most thermophilic proteins tend to have more salt bridges, and achieve higher thermostability by up-shifting and broadening their protein stability curves. While the stabilizing effect of salt-bridge has been extensively studied, experimental data on how salt-bridge influences protein stability curves are scarce. Here, we used double mutant cycles to determine the temperature-dependency of the pair-wise interaction energy and the contribution of salt-bridges to ΔC_p in a thermophilic ribosomal protein L30e. Our results showed that the pair-wise interaction energies for the salt-bridges E6/R92 and E62/K46 were stabilizing and insensitive to temperature changes from 298 to 348 K. On the other hand, the pair-wise interaction energies between the control long-range ion-pair of E90/R92 were negligible. The ΔC_p of all single and double mutants were determined by Gibbs-Helmholtz and Kirchhoff analyses. We showed that the two stabilizing salt-bridges contributed to a reduction of ΔC_p by 0.8–1.0 kJ mol⁻¹ K⁻¹. Taken together, our results suggest that the extra salt-bridges found in thermophilic proteins enhance the thermostability of proteins by reducing ΔC_p, leading to the up-shifting and broadening of the protein stability curves.     

NMR Structure of Integrin α4 Cytosolic Tail and Its Interactions with Paxillin

Background

Integrins are a group of transmembrane signaling proteins that are important in biological processes such as cell adhesion, proliferation and migration. Integrins are α/β hetero-dimers and there are 24 different integrins formed by specific combinations of 18 α and 8 β subunits in humans. Generally, each of these subunits has a large extracellular domain, a single pass transmembrane segment and a cytosolic tail (CT). CTs of integrins are important in bidirectional signal transduction and they associate with a large number of intracellular proteins.

Principal Findings

Using NMR spectroscopy, we determined the 3-D structure of the full-length α4 CT (Lys968-Asp999) and characterize its interactions with the adaptor protein paxillin. The α4 CT assumes an overall helical structure with a kink in its membrane proximal region. Residues Gln981-Asn997 formed a continuous helical conformation that may be sustained by potential ionic and/or hydrogen bond interactions and packing of aromatic-aliphatic side-chains. ¹⁵N-¹H HSQC NMR experiments reveal interactions of the α4 CT C-terminal region with a fragment of paxillin (residues G139-K277) that encompassed LD2-LD4 repeats. Residues of these LD repeats including their adjoining linkers showed α4 CT binding-induced chemical shift changes. Furthermore, NMR studies using LD-containing peptides showed predominant interactions between LD3 and LD4 of paxillin and α4 CT. Docked structures of the α4 CT with these LD repeats suggest possible polar and/or salt-bridge and non-polar packing interactions.

Significance

The current study provides molecular insights into the structural diversity of α CTs of integrins and interactions of integrin α4 CT with the adaptor protein paxillin.     

Fast photochemical oxidation of proteins for comparing solvent-accessibility changes accompanying protein folding: Data processing and application to barstar

Mass spectrometry-based protein footprinting reveals regional and even amino-acid structural changes and fills the gap for many proteins and protein interactions that cannot be studied by X-ray crystallography or NMR spectroscopy. Hydroxyl radical-mediated labeling has proven to be particularly informative in this pursuit because many solvent-accessible residues can be labeled by OH in a protein or protein complex, thus providing more coverage than does specific amino-acid modifications. Finding all the OH-labeling sites requires LC/MS/MS analysis of a proteolyzed sample, but data processing is daunting without the help of automated software. We describe here a systematic means for achieving a comprehensive residue-resolved analysis of footprinting data in an efficient manner, utilizing software common to proteomics core laboratories. To demonstrate the method and the utility of OH-mediated labeling, we show that FPOP easily distinguishes the buried and exposed residues of barstar in its folded and unfolded states. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.     

The PRC-barrel: a widespread,conserved domain shared by photosynthetic reaction center subunits and proteins of RNA metabolism

Background

The H subunit of the purple bacterial photosynthetic reaction center (PRC-H) is important for the assembly of the photosynthetic reaction center and appears to regulate electron transfer during the reduction of the secondary quinone. It contains a distinct cytoplasmic β-barrel domain whose fold has no close structural relationship to any other well known β-barrel domain.

Results

We show that the PRC-H β-barrel domain is the prototype of a novel superfamily of protein domains, the PRC-barrels, approximately 80 residues long, which is widely represented in bacteria, archaea and plants. This domain is also present at the carboxyl terminus of the pan-bacterial protein RimM, which is involved in ribosomal maturation and processing of 16S rRNA. A family of small proteins conserved in all known euryarchaea are composed entirely of a single stand-alone copy of the domain. Versions of this domain from photosynthetic proteobacteria contain a conserved acidic residue that is thought to regulate the reduction of quinones in the light-induced electron-transfer reaction. Closely related forms containing this acidic residue are also found in several non-photosynthetic bacteria, as well as in cyanobacteria, which have reaction centers with a different organization. We also show that the domain contains several determinants that could mediate specific protein-protein interactions.

Conclusions

The PRC-barrel is a widespread, ancient domain that appears to have been recruited to a variety of biological systems, ranging from RNA processing to photosynthesis. Identification of this versatile domain in numerous proteins could aid investigation of unexplored aspects of their biology.

     

Energetic contribution of solvent-exposed ion pairs to alpha-helix structure.

Understanding the role of amino acid side-chain interactions in forming secondary structure in proteins is useful for deciphering how proteins fold and for predicting folded structures of proteins from their sequence. Analysis of the secondary structure as a function of pH in two designed synthetic peptides with identical composition but different sequences, affords a quantitative estimate of the free energy contribution of a single ion pair to the stability of an isolated alpha-helix. One peptide contains repeated blocks of Glu4Lys4. The second has repeated blocks of Glu2Lys2. The former contains significant helical structure at neutral pH while the latter has none, based on ultraviolet light circular dichroism measurements and 1H nuclear magnetic resonance spectroscopy. The difference is attributed to formation of helix-stabilizing salt-bridges between Glu- and Lys+ spaced at i, i + 4 intervals in the former peptide. The free energy of formation of a single Glu(-)-Lys+ salt-bridge can be evaluated by using a statistical model of the helix-coil transition that explicitly includes salt-bridges: the result is -0.50(+/- 0.05) kcal/mol at 4 degrees C and neutral pH in 10 mM salt, in agreement with a value derived for a single salt-bridge in a helix on the surface of a globular protein.