Amal Kumar Bandyopadhyay

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−1. 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−1 per salt-bridge, seems to be a halophilic design wherein favorable average contribution of background term (−10 kcal mol−1) exceeds than that of bridge term (−7 kcal mol−1). 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.

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

Salt-bridge and network salt-bridge are specific electrostatic interactions that contribute to the overall stability of proteins. In hierarchical protein folding model, these interactions play crucial role in nucleation process. The advent and growth of protein structure database and its availability in public domain made an urgent need for context dependent rapid analysis of salt-bridges. While these analyses on single protein is cumbersome and time-consuming, batch analyses need efficient software for rapid topological scan of a large number of protein for extracting details on (i) fraction of salt-bridge residues (acidic and basic). (ii) Chain specific intra-molecular salt-bridges, (iii) inter-molecular salt-bridges (protein-protein interactions) in all possible binary combinations (iv) network salt-bridges and (v) secondary structure distribution of salt-bridge residues. To the best of our knowledge, such efficient software is not available in public domain. At this juncture, we have developed a program i.e. SBION which can perform all the above mentioned computations for any number of protein with any number of chain at any given distance of ion-pair. It is highly efficient, fast, error-free and user friendly. Finally we would say that our SBION indeed possesses potential for applications in the field of structural and comparative bioinformatics studies. Availability SBION is freely available for non-commercial/academic institutions on formal request to the corresponding author (akbanerjee@biotech.buruniv.ac.in).

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. Availability SBION2 is freely available at http://sourceforge.net/projects/sbion2/ for academic users

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 akbanerjee@biotech.buruniv.ac.in

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

BackgroundJapanese 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.ResultsAlthough 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.ConclusionsThe 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.

PHYSICO2: an UNIX based standalone procedure for computation of physicochemical, window-dependent and substitution based evolutionary properties of protein sequences along with automated block preparation tool, version 2.

Automated genome sequencing procedure is enriching the sequence database very fast. To achieve a balance between the entry of sequences in the database and their analyses, efficient software is required. In this end PHYSICO2, compare to earlier PHYSICO and other public domain tools, is most efficient in that it i] extracts physicochemical, window-dependent and homologousposition-based-substitution (PWS) properties including positional and BLOCK-specific diversity and conservation, ii] provides users with optional-flexibility in setting relevant input-parameters, iii] helps users to prepare BLOCK-FASTA-file by the use of Automated Block Preparation Tool of the program, iv] performs fast, accurate and user-friendly analyses and v] redirects itemized outputs in excel format along with detailed methodology. The program package contains documentation describing application of methods. Overall the program acts as efficient PWS-analyzer and finds application in sequence-bioinformatics. Availability PHYSICO2: is freely available at http://sourceforge.net/projects/physico2/ along with its documentation at https://sourceforge.net/projects/physico2/files/Documentation.pdf/download for all users.

ADSBET2: Automated Determination of Salt-Bridge Energy-Terms version 2.

Component (bridge: ΔΔGbrd, background: ΔΔGprot, desolvation: ΔΔGdsolv) and net (ΔΔGnet) energy-terms of salt-bridge-structure (SBS) are auto-generated by the program ADSBET that makes use of general purpose Adaptive Poison Boltzmann Solver (APBS) method. While the procedure reports gross energy terms (Kcal Mol-1), report on bond-multiplicity corrected normalized energyterms (Kcal Mol-1 Bond-1) along with their accessibility (ASA) in monomer, isolated-SBS (ISBS) and networked-SBS (NSBS) format would be very useful for statistical comparison among SBSs and understanding their location in protein structure. In this end, ADSBET2 potentially incorporates these features along with additional model for side-chain. Gross and normalized energy-terms are redirected in monomer, ISBS and NSBS format along with their ASA informations. It works on any number of SBSs for any number of structure files present in a database. Taken together, ADSBET2 has been suitable for statistical analyses of SBSs energetics and finds applications in protein engineering and structural bioinformatics. Availability ADSBET2 is freely available at http://sourceforge.net/projects/ADSBET2/ for all users.

PROPAB: Computation of Propensities and Other Properties from Segments of 3D structure of Proteins

Residues in allelic positions, in the local segment of aligned sequences of proteins show wide variations. Here, we describe PROPAB that computes the propensity tables for helix, strand and coil types from multiple 3D structure files following ab initio statistical procedure. It also classifies them in range specific and chain specific manners. It further computes percentage composition and physicochemical properties along with residues propensities. It also prepares FASTA files for different segments (helix, strand and coil) in the exact order that they follow in the sequence. Representative analyses on orthologous (homologous across species) proteins demonstrate wide segmental variations of physicochemical properties. Such variations provide insights to relate the adaptation of these proteins in a given functional constraint under diverse environmental conditions. Thus, the program finds applications in the structural and evolutionary analysis of proteins. Availability: PROPAB is freely available at http://sourceforge.net/projects/propab/for worldwide user.

Synthesis, molecular modeling of N-acyl benzoazetinones and their docking simulation on fungal modeled target

ABSTRACT A series of stable N-acyl benzoazetinones have been synthesized in moderate to good yields (58–85%) from easily available substrates such as 2-(N-acyl) amino benzoic acids through intramolecular amidation under mild conditions. These geometry-optimized benzoazetinones were docked in the model target of P450, class CYP53A15, a benzoate 4-monooxygenase abundantly found in the genome of ascomycetes and Basidiomycetes classes of pathogenic fungi. Low per residue root-mean-square deviation (RMSD) of modeled structure of the enzyme indicated similar topology as template (4D6Z.pdb). Observed score judges site-specific docking, and the interaction of quantum mechanically optimized benzoazetinone derivatives with the target enzyme. These results suggest that 3i is the best antifungal agent. The specific hydrophobic substituent in the benzoazetinones contributed to the stability of ligand–target complex. Overall, the study provided insight into the specificity of the site-specific interactions, thereby, facilitating the possibility of development of broad-spectrum antifungal agents against opportunistic and infectious fungi. GRAPHICAL ABSTRACT

Substitutional Analysis of Orthologous Protein Families Using BLOCKS

Orthologous proteins, form due to divergence of parental sequence, perform similar function under different environmental and biological conditions. Amino acid changes at locus specific positions form hetero-pairs whose role in BLOCK evolution is yet to be understood. We involve eight protein BLOCKs of known divergence rate to gain insight into the role of hetero-pairs in evolution. Our procedure APBEST uses BLOCK-FASTA file to extract BLOCK specific evolutionary parameters such as dominantly used hetero-pair (D), usage of hetero-pairs (E), non-conservative to conservative substitution ratio (R), maximally-diverse residue (MDR), residue (RD) and class (CD) specific diversity. All these parameters show BLOCK specific variation. Conservative nature of D points towards restoration of function of BLOCK. While E sets the upper-limit of usage of hereto-pairs, strong correlation of R with divergence-rate indicates that the later is directly dependent on non-conservative substitutions. The observation that MDR, measure of positional diversity, occupy very limited positions in BLOCK indicates accommodation of diversity is positionally restricted. Overall, the study extract observed hetero-pair related quantitative and multi-parametric details of BLOCK, which finds application in evolutionary biology.