SBION: A Program for Analyses of Salt-Bridges from Multiple Structure Files

Underlying features for the enhanced electrostatic strength of the extremophilic malate dehydrogenase interface salt-bridge compared to the mesophilic one

Malate dehydrogenase (MDH) exists in multimeric form in normal and extreme solvent conditions where residues of the interface are involved in specific interactions. The interface salt-bridge (ISB) and its microenvironment (ME) residues may have a crucial role in the stability and specificity of the interface. To gain insight into this, we have analyzed 218 ISBs from 42 interfaces of 15 crystal structures along with their sequences. Comparative analyses demonstrate that the ISB strength is ∼30 times greater in extremophilic cases than that of the normal one. To this end, the interface residue propensity, ISB design and pair selection, and ME-residue's types, i.e., type-I and type-II, are seen to be intrinsically involved. Although Type-I is a common type, Type-II appears to be extremophile-specific, where the net ME-residue count is much lower with an excessive net ME-energy contribution, which seems to be a novel interface compaction strategy. Furthermore, the interface strength can be enhanced by selecting the desired mutant from the net-energy profile of all possible mutations of an unfavorable ME-residue. The study that applies to other similar systems finds applications in protein-protein interaction and protein engineering.Communicated by Ramaswamy H. Sarma.

A computational study of structural analysis of Class I human glucose-6-phosphate dehydrogenase (G6PD) variants: Elaborating the correlation to chronic non-spherocytic hemolytic anemia (CNSHA)

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common human enzyme defect that affects more than 500 million people worldwide. Individuals affected with G6PD deficiency may occasionally suffer mild-to-severe chronic hemolytic anemia. Chronic non-spherocytic hemolytic anemia (CNSHA) is a potential result of the Class I G6PD variants. This comparative computational study attempted to correct the defect in variants structure by docking the AG1 molecule to selected Class I G6PD variants [G6PD^Nashville (Arg393His), G6PD^Alhambra (Val394Leu), and G6PD^Durham (Lys238Arg)] at the dimer interface and structural NADP⁺ binding site. It was followed by an analysis of the enzyme conformations before and after binding to the AG1 molecule using the molecular dynamics simulation (MDS) approach, while the severity of CNSHA was determined via root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), hydrogen bonds, salt bridges, radius of gyration (Rg), solvent accessible surface area analysis (SASA), and principal component analysis (PCA). The results revealed that G6PD^Nashville (Arg393His) and G6PD^Durham (Lys238Arg) had lost the direct contact with structural NADP⁺ and salt bridges at Glu419 - Arg427 and Glu206 - Lys407 were disrupted in all selected variants. Furthermore, the AG1 molecule re-stabilized the enzyme structure by restoring the missing interactions. Bioinformatics approaches were also used to conduct a detailed structural analysis of the G6PD enzyme at a molecular level to understand the implications of these variants toward enzyme function. Our findings suggest that despite the lack of treatment for G6PDD to date, AG1 remains a novel molecule that promotes activation in a variety of G6PD variants.

ProteinTools: a toolkit to analyze protein structures

The experimental characterization and computational prediction of protein structures has become increasingly rapid and precise. However, the analysis of protein structures often requires researchers to use several software packages or web servers, which complicates matters. To provide long-established structural analyses in a modern, easy-to-use interface, we implemented ProteinTools, a web server toolkit for protein structure analysis. ProteinTools gathers four applications so far, namely the identification of hydrophobic clusters, hydrogen bond networks, salt bridges, and contact maps. In all cases, the input data is a PDB identifier or an uploaded structure, whereas the output is an interactive dynamic web interface. Thanks to the modular nature of ProteinTools, the addition of new applications will become an easy task. Given the current need to have these tools in a single, fast, and interpretable interface, we believe that ProteinTools will become an essential toolkit for the wider protein research community. The web server is available at https://proteintools.uni-bayreuth.de.

Intrinsic basis of thermostability of prolyl oligopeptidase from Pyrococcus furiosus

Salt-bridges play a key role in the thermostability of proteins adapted in stress environments whose intrinsic basis remains to be understood. We find that the higher hydrophilicity of PfP than that of HuP is due to the charged but not the polar residues. The primary role of these residues is to enhance the salt-bridges and their ME. Unlike HuP, PfP has made many changes in its intrinsic property to strengthen the salt-bridge. First, the desolvation energy is reduced by directing the salt-bridge towards the surface. Second, it has made bridge-energy more favorable by recruiting energetically advantageous partners with high helix-propensity among the six possible salt-bridge pairs. Third, ME-residues that perform intricate interactions have increased their energy contribution by making major changes in their binary properties. The use of salt-bridge partners as ME-residues, and ME-residues' overlapping usage, predominant in helices, and energetically favorable substitution are some of the favorable features of PfP compared to HuP. These changes in PfP reduce the unfavorable, increase the favorable ME-energy. Thus, the per salt-bridge stability of PfP is greater than that of HuP. Further, unfavorable target ME-residues can be identified whose mutation can increase the stability of salt-bridge. The study applies to other similar systems.

AUTOMINv1.0: an automation for minimization of Protein Data Bank files and its usage

Global minimal structure of protein/enzyme is energetically compromised that maintains an intricate balance between the rigidity and the flexibility. Such a state makes it interactive to its ligand molecules. Although protein data bank files (PDB) may have achieved the state, in many situations minimization has been crucial to overcome unwanted steric clashes, and other conformational strains. It is more so, when orthologous PDB structures that are intended in a given study, show variations in resolution, R-factor, shell-water contents, loop characteristics etc. Here, a fully automated procedure of minimization would be highly useful. AUTOMINv1.0 is such an automation of minimization that runs on any number of structure files with any number of chains in them along with the inclusion of selective/non-selective shell-waters interacting with polar and or non-polar atom-types of protein. Comparison of the mean binaryitems of salt-bridges of minimized and un-minimized structures (chains > 100) of nucleoside diphosphate kinase from mimi virus shows dramatic improvements in the earlier. Again, the mean steric clashes of 2AZ3.pdb are reduced upon minimization. Remarkably, the observed steric clashes between shell-waters and atom-types of protein are seen to be removed upon minimization. Taken together, AUTOMINv1.0 is an automation of minimization that finds applications in structural bioinformatics.

Insights from the salt bridge analysis of malate dehydrogenase from H. salinarum and E.coli

Halophilic proteins have greater abundance of acidic over basic residues in sequence. In structure, the surface is decorated by negative charges, with lower content of Lysine. Using sequence BLOCKs and 3D model of malate dehydrogenase from halophilic archaea (Halobacterium salinarum; hsMDH) and X-ray structure from mesophilic bacteria (E. coli; ecMDH), we show that not only acidic and basic residues have higher mean relative abundance (MRA) and thus, impart higher polarity to the sequences, but also show their presence in the surface of the structure of hsMDH relative to its mesophilic counterpart. These observations may indicate that both the acidic and the basic residues have a concerted role in the stability of hsMDH. Analysis on salt bridges from hsMDH and ecMDH show that in the former, salt bridges are highly intricate, newly engineered and global in nature. Although, these salt bridges are abundant in hsMDH, in the active site the design remains unperturbed. In high salt where hydrophobic force is weak, these salt bridges seem to play a major role in the haloadaptation of the tertiary structure of hsMDH. This is the first report of such an observation.

Stability of buried and networked salt-bridges (BNSB)in thermophilic proteins

Thermophilic proteins function at high temperature, unlike mesophilic proteins. Thermo-stability of these proteins is due to the unique buried and networked salt-bridge (BNSB). However, it is known that the desolvation cost of BNSB is too high compared to other favorable energy terms. Nonetheless, it is known that stability is provided generally by hydrophobic isosteres without the need for BNSB. We show in this analysis that the BNSB is the optimal evolutionary design of salt-bridge to offset desolvation cost efficiently. Hence, thermophilic proteins with a high level of BNSB provide thermo-stability.

Analysis of salt-bridges in prolyl oligopeptidase from Pyrococcus furiosus and Homo sapiens

Hyper thermophilic archaea not only tolerate high temperature but also operate its biochemical machineries, normally under these conditions. However, the structural signatures in proteins that answer for the hyper thermo-stability relative to its mesophilic homologue remains poorly understood. We present comparative analyses of sequences, structures and salt-bridges of prolyl-oligopeptidase from Pyrococcus furiosus (pfPOP - PDB ID: 5T88) and human (huPOP - PDB ID: 3DDU). A similar level of hydrophobic and hydrophilic residues in pfPOP and huPOP is observed. A low level of interactions between shell-waters and atom-types in pfPOP indicated hyper thermophilic features are negligible. Salt-bridge-forming-residues (sbfrs) are high in pfPOP's core and surface (pfPOP). Increased sbfrs largely indicate specific-electrostatic is important for thermo stability in pfPOP. Four classes of sbfrs are found namely positionally non-conservative (PNCS), conservative (PCS), unchanged (PU) and interchanged (PIC) type of substitutions. PNCS-sbfrs constitutes 28% and it is associated with the topology of pfPOP at high temperature. PCS helps to increase the salt-bridge population. It is also found that PU maintains similar salt-bridges at the active site and other binding sites while PIC abolishes mesophilic patterns.

POWAINDv1.0: A Program for Protein-Water Interactions Determination

Protein is the most exposed biomolecule in the aqueous environment of the cell. Its structure maintains a delicate balance between the rigidity and the flexibility that imparts binding specificity to its substrate/ligand, etc. Intramolecular interactions of polar and non-polar groups of amino acid residues and intermolecular weak interactions between these groups and shell-waters may contribute to the overall stability of the tertiary structure. However, the question as to what are the dynamics of interactions of shell-water with respect to weak forces and atom-groups of protein (AGP), requires systematic investigations. In this end, we have developed a procedure POWAINDv1.0 that analyzes interactions of crystallographic shell-waters (CSH) in residues and AGP specific manner. The shell-water and AGP specific bridge-interactions are also extracted. Further, the program analyzes favorable and unfavorable nature of each interaction based on the actual and 75% of the sum of van der Waals (vdW) radii of interacting atoms. The EXCEL-outputs are useful in understanding the profile for AGP-CSH interactions and contribution of each component in AGP. Taken together, the program provides intricate details on CSHprotein interactions and finds application in the structural Bioinformatics

Structural insights from water-ferredoxin interaction in mesophilic algae and halophilic archaea

We analyzed the water-ferredoxin interaction in mesophilic (moderate temperature) algae (PDB ID: 1AWD) and halophilic (salt-tolerant) archaea (PDB ID: 1DOI) using POWAIND version 2.0 (a protein-water interactions calculation program). It is found that the shell water (SW) is 2.5 fold greater in halophilic ferredoxin than mesophilic ferredoxin. Water-ferredoxin interactions in the core and cavity are the signature of stability. The normalized frequency of such interactions is less in halophilic relative to mesophilic ferredoxin and the halophilic signature for stability by such interactions is negligible. However, the surface dominated with such interactions seems to be important for ferredoxin and oxido-reductase recognition.

Characterization of Ionizable Groups' Environments in Proteins and Protein-Ligand Complexes through a Statistical Analysis of the Protein Data Bank

We conduct a statistical analysis of the molecular environment of common ionizable functional groups in both protein-ligand complexes and inside proteins from the Protein Data Bank (PDB). In particular, we characterize the frequency, type, and density of the interacting atoms as well as the presence of a potential counterion. We found that for ligands, most guanidinium groups, half of primary and secondary amines, and one-fourth of imidazole neighbor a carboxylate group. Tertiary amines bind more rarely near carboxylate groups, which may be explained by a crowded neighborhood and hydrophobic character. In comparison to the environment seen by the ligands, inside proteins, an environment enriched in main-chain atoms is found, and the prevalence of direct charge neutralization by carboxylate groups is different. When the ionizable character of water molecules and phenolic or hydroxyl groups is accounted, considering a high-resolution dataset (less than 1.5 Å), charge neutralization could occur for well above 80% of the ligand functional groups considered, but for tertiary amines.

Insight into SNPs and epitopes of E protein of newly emerged genotype-I isolates of JEV from Midnapur, West Bengal, India

BACKGROUND

Japanese encephalitis virus (JEV) is a mosquito-borne flavivirus that causes Japanese Encephalitis (JE) and Acute Encephalitis Syndrome (AES) in humans. Genotype-I (as co-circulating cases with Genotype-III) was isolated in 2010 (JEV28, JEV21) and then in 2011 (JEV45) from Midnapur district, West Bengal (WB) for the first time from clinical patients who were previously been vaccinated with live attenuated SA14-14-2 strain. We apply bioinformatics and immunoinformatics on sequence and structure of E protein for analysis of crucial substitutions that might cause the genotypic transition, affecting protein-function and altering specificity of epitopes.

RESULTS

Although frequency of substitutions in E glycoprotein of JEV28, JEV21 and JEV45 isolates vary, its homologous patterns remain exactly similar as earlier Japan isolate (Ishikawa). Sequence and 3D model-structure based analyses of E protein show that only four of all substitutions are critical for genotype-I specific effect of which N103K is common among all isolates indicating its role in the transition of genotype-III to genotype-I. Predicted B-cell and T-cell epitopes are seen to harbor these critical substitutions that affect overall conformational stability of the protein. These epitopes were subjected to conservation analyses using a large set of the protein from Asian continent.

CONCLUSIONS

The study identifies crucial substitutions that contribute to the emergence of genotype-I. Predicted epitopes harboring these substitutions may alter specificity which might be the reason of reported failure of vaccine. Conservation analysis of these epitopes would be useful for design of genotype-I specific vaccine.

Discrete structural features among interface residue-level classes

BACKGROUND

Protein-protein interaction (PPI) is essential for molecular functions in biological cells. Investigation on protein interfaces of known complexes is an important step towards deciphering the driving forces of PPIs. Each PPI complex is specific, sensitive and selective to binding. Therefore, we have estimated the relative difference in percentage of polar residues between surface and the interface for each complex in a non-redundant heterodimer dataset of 278 complexes to understand the predominant forces driving binding.

RESULTS

Our analysis showed ~60% of protein complexes with surface polarity greater than interface polarity (designated as class A). However, a considerable number of complexes (~40%) have interface polarity greater than surface polarity, (designated as class B), with a significantly different p-value of 1.66E-45 from class A. Comprehensive analyses of protein complexes show that interface features such as interface area, interface polarity abundance, solvation free energy gain upon interface formation, binding energy and the percentage of interface charged residue abundance distinguish among class A and class B complexes, while electrostatic visualization maps also help differentiate interface classes among complexes.

CONCLUSIONS

Class A complexes are classical with abundant non-polar interactions at the interface; however class B complexes have abundant polar interactions at the interface, similar to protein surface characteristics. Five physicochemical interface features analyzed from the protein heterodimer dataset are discriminatory among the interface residue-level classes. These novel observations find application in developing residue-level models for protein-protein binding prediction, protein-protein docking studies and interface inhibitor design as drugs.

Linking structural features of protein complexes and biological function

Protein-protein interaction (PPI) establishes the central basis for complex cellular networks in a biological cell. Association of proteins with other proteins occurs at varying affinities, yet with a high degree of specificity. PPIs lead to diverse functionality such as catalysis, regulation, signaling, immunity, and inhibition, playing a crucial role in functional genomics. The molecular principle of such interactions is often elusive in nature. Therefore, a comprehensive analysis of known protein complexes from the Protein Data Bank (PDB) is essential for the characterization of structural interface features to determine structure-function relationship. Thus, we analyzed a nonredundant dataset of 278 heterodimer protein complexes, categorized into major functional classes, for distinguishing features. Interestingly, our analysis has identified five key features (interface area, interface polar residue abundance, hydrogen bonds, solvation free energy gain from interface formation, and binding energy) that are discriminatory among the functional classes using Kruskal-Wallis rank sum test. Significant correlations between these PPI interface features amongst functional categories are also documented. Salt bridges correlate with interface area in regulator-inhibitors (r = 0.75). These representative features have implications for the prediction of potential function of novel protein complexes. The results provide molecular insights for better understanding of PPIs and their relation to biological functions.

SBION2: Analyses of Salt Bridges from Multiple Structure Files, Version 2

Specific electrostatics (i.e. salt-bridge) includes both local and non-local interactions that contribute to the overall stability of proteins. It has been shown that a salt-bridge could either be buried or exposed, networked or isolated, hydrogen-bonded or nonhydrogen bonded, in secondary-structure or in coil, formed by single or multiple bonds. Further it could also participates either in intra- or inter-dipole interactions with preference in orientation either for basic residue at N-terminal (orientation-I) or acidic residue at N-terminal (orientation-II). In this context SBION2 is unique in that it reports above mentioned binary items in excel format along with details on intra and inter-dipole interactions and orientations. These results are suitable for post run statistical analyses involving large datasets. Reports are also made on protein-protein interactions, intervening residue distances and general residue specific salt-bridge details. A ready to use compact supplementary table is also produced. The program runs in three alternative modes. Each mode works on any number of structure files with any number of chains at any given atomic distance of ion-pair. Thus SBION2 provides intricate details on salt-bridges and finds application in structural bioinformatics.

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.