PocketAlign A Novel Algorithm for Aligning Binding Sites in Protein Structures

A fundamental task in bioinformatics involves a transfer of knowledge from one protein molecule onto another by way of recognizing similarities. Such similarities are obtained at different levels, that of sequence, whole fold, or important substructures. Comparison of binding sites is important to understand functional similarities among the proteins and also to understand drug cross-reactivities. Current methods in literature have their own merits and demerits, warranting exploration of newer concepts and algorithms, especially for large-scale comparisons and for obtaining accurate residue-wise mappings. Here, we report the development of a new algorithm, PocketAlign, for obtaining structural superpositions of binding sites. The software is available as a web-service at http://proline.physicslisc.emetin/pocketalign/. The algorithm encodes shape descriptors in the form of geometric perspectives, supplemented by chemical group classification. The shape descriptor considers several perspectives with each residue as the focus and captures relative distribution of residues around it in a given site. Residue-wise pairings are computed by comparing the set of perspectives of the first site with that of the second, followed by a greedy approach that incrementally combines residue pairings into a mapping. The mappings in different frames are then evaluated by different metrics encoding the extent of alignment of individual geometric perspectives. Different initial seed alignments are computed, each subsequently extended by detecting consequential atomic alignments in a three-dimensional grid, and the best 500 stored in a database. Alignments are then ranked, and the top scoring alignments reported, which are then streamed into Pymol for visualization and analyses. The method is validated for accuracy and sensitivity and benchmarked against existing methods. An advantage of PocketAlign, as compared to some of the existing tools available for binding site comparison in literature, is that it explores different schemes for identifying an alignment thus has a better potential to capture similarities in ligand recognition abilities. PocketAlign, by finding a detailed alignment of a pair of sites, provides insights as to why two sites are similar and which set of residues and atoms contribute to the similarity.

Dynamic and thermodynamic anomalies of water at low temperatures: from bulk water to reverse micelles and DNA hydration layer

Liquid water is known to exhibit remarkable thermodynamic and dynamic anomalies, ranging from solvation properties in supercritical state to an apparent divergence of the linear response functions at a low temperature. Anomalies in various dynamic properties of water have also been observed in the hydration layer of proteins, DNA grooves and inside the nanocavity, such as reverse micelles and nanotubes. Here we report studies on the molecular origin of these anomalies in supercooled water, in the grooves of DNA double helix and reverse micelles. The anomalies have been discussed in terms of growing correlation length and intermittent population fluctuation of 4- and 5-coordinated species. We establish correlation between thermodynamic response functions and mean squared species number fluctuation. Lifetime analysis of 4- and 5-coordinated species reveals interesting differences between the role of the two species in supercooled and constrained water. The nature and manifestations of the apparent and much discussed liquid-liquid transition under confinement are found to be markedly different from that in the bulk. We find an interesting `faster than bulk' relaxation in reverse micelles which we attribute to frustration effects created by competition between the correlations imposed by surface interactions and that imposed by hydrogen bond network of water.

Prediction of activities of oxygen in dilute quaternary solutions using binary data

Equations are developed for predicting the activity coefficients of oxygen dissolved in ternary liquid alloys. These are extensions of earlier treatments, and are based on a model in which each oxygen atom is assumed to make four bonds with neighboring metal atoms. It is also postulated that the strong oxygen-metal bonds distort the electronic configuration around the metal atoms bonded to oxygen, and that the quantitative reduction of the strength of bonds made by these atoms with all of the adjacent metal atoms is equivalent to a factor of approximately two. The predictions of the quasichemical equation which is derived agree satisfactorily with the partial molar free energies of oxygen in Ag-Cu-Sn solutions at 1200°C reported in literature. An extension of this treatment to multicomponent solutions is also indicated.

Evidence of residual entropy in the cubic spinel Zn2TiO4

The standard free energies of formation of Zn2Ti04 and ZnTi03 have been determined in the temperature range 930° to i ioo'x from electromotive force measurements on reversible solid oxide galvanic cells;Ag-5at%znll I Pt, + CaO-Zr02 ZnO I II Ag-5at%Zn Y20r Th02 CaO-Zr02 + ,Pt Zn2Ti04+ ZnTi03 and II Ag-5at%Zn CaO-Zr02 + ,Pt ZnTi03+ Ti02 The values may be expressed by the equations,2ZnO (wurtz) + Ti02(rut) -> Zn2Ti04(sp), f:!:.Go = -750-2-46T (±75)cal;ZnO(wurtz) +Ti02(rut) -> ZnTi03(ilmen) ,f:!:.Co = -]600-0·]99T(±50)cal.Combination of the free energy values with the calorimetric heat of formation, and low-temperature and high-temperature heat capacity of Zn2Ti04 reported in literature, suggests a residual entropy of ],9 (±0·6) cal K-1 mol ? for the cubic spinel. Ideal mixing of Zn2+ and Ti4+ ions on the octahedral sites would result in a configurational contribution to the entropy of 2· 75 cal K-1 rnol ".The difference is indicative of short-range ordering of cations on octahedral sites.

Activities and their relation to cation distribution in NiAl2O4 MgAl2O4 spinel solid solutions

The activity of NiAl2O4 in NiAl2O4MgAl2O4 solid solutions has been measured by using a solid oxide galvanic cell of the type, Pt, Ni + NiAl2O4 + Al2O3(α)/CaOZrO2/Ni + NixMg1−xAl2O4 + Al2O3(α). Pt, in the temperature range 750–1150°C. The activities in the spinel solid solutions show negative deviations from Raoult's law. The cation distribution in the solid solutions has been calculated using site preference energies independent of composition for Ni2+, Mg2+, and Al3+ ions obtained from crystal field theory and measured cation disorder in pure NiAl2O4 and MgAl2O4, and assumi g ideal mixing of cations on the tetrahedral and octahedral positions. The calculated values correctly predict the decrease in the fraction, α, of Ni2+ ions on tetrahedral sites for 1>x>0.25, observed by Porta et al. [J. Solid State Chem.11, 135 (1974)] but do not support their tentative evidence for an increase in α for x < 0.25. The measured excess free energy of mixing can be completely accounted for by using either the calculated or the measured cation distributions. This suggests that the Madelung energy is approximately a linear function of composition in the solid solutions. The composition of NiOMgO solid solutions in equilibrium with NiAl2O4MgAl2O4 solid solutions has been calculated from the results and information available in literature.

Stability numbers for an unsupported vertical circular excavation in c-phi soil

By using an axisymmetric lower bound finite element limit analysis formulation, the stability numbers (gamma H/C) for an unsupported vertical circular excavation in a cohesive-frictional soil have been generated. The numerical results are obtained for values of normalized excavation height (H/b) and friction angle (phi) greater than those considered previously in the literature. The results compare well with those available in literature. The stability numbers presented in this note would be beneficial from a design point of view. (C) 2011 Elsevier Ltd. All rights reserved.

Tar characterisation in new generation agro-residue gasifiers - Cyclone and downdraft open top twin air entry systems

This paper deals with the characterisation of tar from two configurations of bioresidue thermochemical conversion reactors designed for producer gas based power generation systems. The pulverised fuel reactor is a cyclone system (R1) and the solid bioresidue reactor (denoted R2) is an open top twin air entry system both at 75-90 kg/h capacity (to generate electricity similar to 100 kVA). The reactor, R2, has undergone rigorous test in a major Indo-Swiss programme for the tar quantity at various conditions. The former is a recent technology development. Tars collected from these systems by a standard tar collection apparatus at the laboratory at Indian Institute of Science have been analysed at the Royal Institute of Technology (KTH), Sweden. The results of these analyses show that these thermochemical conversion reactors behave differently from the earlier reactors reported in literature in so far as tar generation is concerned. The extent of tar in hot gas is about 700-800 ppm for R1 and 70-100 ppm for R2. The amounts of the major compounds - naphthalene and phenol-are much lower that what is generally understood to happen in the gasifiers in Europe. It is suggested that the longer residence times at high temperatures allowed for in these reactors is responsible for this behavior. It is concluded the new generation reactor concepts extensively tried out at lower power levels hold promise for high power atmospheric gasification systems for woody as well as pulverisable bioresidues.

From anomalies in neat liquid to structure, dynamics and function in the biological world

Water brings its remarkable thermodynamic and dynamic anomalies in the pure liquid state to biological world where water molecules face a multitude of additional interactions that frustrate its hydrogen bond network. Yet the water molecules participate and control enormous number of biological processes in manners which are yet to be understood at a molecular level. We discuss thermodynamics, structure, dynamics and properties of water around proteins and DNA, along with those in reverse micelles. We discuss the roles of water in enzyme kinetics, in drug-DNA intercalation and in kinetic-proof reading ( the theory of lack of errors in biosynthesis). We also discuss how water may play an important role in the natural selection of biomolecules. (C) 2011 Elsevier B. V. All rights reserved.

Modeling of ternary solubilities of solids in supercritical carbon dioxide in the presence of cosolvents or cosolutes

The investigation of ternary solubilities of solids is essential for the efficient design of extraction processes. The ternary solubilities of solids for cosolvent and cosolute systems are complex functions of temperature, pressure and cosolvent/cosolute composition. The intermolecular interactions between the molecules have a significant role in the solubilities of mixed solids in SCCO2 and cosolvent ternary systems. Two model equations were developed for ternary SCCO2 + cosolvent/cosolute systems by using association and activity coefficient models. Both the model equations consist of five adjustable parameters and correlate the ternary solubilities of solids in terms of temperature, pressure, density and cosolvent/cosolute composition. The model equation for cosolvent systems correlated 43 solid pollutants-cosolvent-SCCO2, while the model equation for cosolute systems correlated 19 solute-cosolute-SCCO2 systems available in literature. The average AARD of the model equations are 4.73% and 4.87% for cosolvent ternary systems and mixed solids in SCCO2, respectively. (C) 2011 Elsevier B.V. All rights reserved.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K ion are comparable to that observed in bulk salt solutions. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4717710]

Phase inversion in agitated liquid-liquid dispersions: Anomalous effect of electrolyte

An imbalance between breakup and coalescence of drops in turbulent liquid-liquid dispersions leads to inversion of phases the dispersed phase becomes continuous and vice versa. An increase in the rate of coalescence of drops is expected to decrease the dispersed phase fraction at which inversion occurs. In the present work, we increased the rate of coalescence of drops by adding electrolyte to pure liquid-liquid dispersions. The experiments carried out for three representative liquid-liquid systems show that contrary to the expectation the addition of an electrolyte increases the dispersed phase fraction at which inversion occurs for both, oil-in-water and water-in-oil dispersions. The step-down experiments confirm that the addition of the electrolyte increases the rate of coalescence of drops in lean dispersions under the same conditions, thereby confirming an anomalous effect of the presence of an electrolyte on the stability of dispersions. (C) 2012 Elsevier Ltd. All rights reserved.

Comparison of tertiary structures of proteins in protein-protein complexes with unbound forms suggests prevalence of allostery in signalling proteins

Abstract: Background: Most signalling and regulatory proteins participate in transient protein-protein interactions during biological processes. They usually serve as key regulators of various cellular processes and are often stable in both protein-bound and unbound forms. Availability of high-resolution structures of their unbound and bound forms provides an opportunity to understand the molecular mechanisms involved. In this work, we have addressed the question "What is the nature, extent, location and functional significance of structural changes which are associated with formation of protein-protein complexes?" Results: A database of 76 non-redundant sets of high resolution 3-D structures of protein-protein complexes, representing diverse functions, and corresponding unbound forms, has been used in this analysis. Structural changes associated with protein-protein complexation have been investigated using structural measures and Protein Blocks description. Our study highlights that significant structural rearrangement occurs on binding at the interface as well as at regions away from the interface to form a highly specific, stable and functional complex. Notably, predominantly unaltered interfaces interact mainly with interfaces undergoing substantial structural alterations, revealing the presence of at least one structural regulatory component in every complex. Interestingly, about one-half of the number of complexes, comprising largely of signalling proteins, show substantial localized structural change at surfaces away from the interface. Normal mode analysis and available information on functions on some of these complexes suggests that many of these changes are allosteric. This change is largely manifest in the proteins whose interfaces are altered upon binding, implicating structural change as the possible trigger of allosteric effect. Although large-scale studies of allostery induced by small-molecule effectors are available in literature, this is, to our knowledge, the first study indicating the prevalence of allostery induced by protein effectors. Conclusions: The enrichment of allosteric sites in signalling proteins, whose mutations commonly lead to diseases such as cancer, provides support for the usage of allosteric modulators in combating these diseases.

Interdiffusion in Nb-Mo, Nb-Ti and Nb-Zr systems

Diffusion couple technique is used to study interdiffusion in Nb-Mo, Nb-Ti and Nb-Zr systems. Interdiffusion coefficients at different temperatures and compositions are determined using the relation developed by Wagner. The change in activation energy for interdiffusion with composition is determined. Further, impurity diffusion coefficient of the species are determined and compared with the available data in literature.