HARMONY: a server for the assessment of protein structures

Protein structure validation is an important step in computational modeling and structure determination. Stereochemical assessment of protein structures examine internal parameters such as bond lengths and Ramachandran (phi, psi) angles. Gross structure prediction methods such as inverse folding procedure and structure determination especially at low resolution can sometimes give rise to models that are incorrect due to assignment of misfolds or mistracing of electron density maps. Such errors are not reflected as strain in internal parameters. HARMONY is a procedure that examines the compatibility between the sequence and the structure of a protein by assigning scores to individual residues and their amino acid exchange patterns after considering their local environments. Local environments are described by the backbone conformation, solvent accessibility and hydrogen bonding patterns. We are now providing HARMONY through a web server such that users can submit their protein structure files and, if required, the alignment of homologous sequences. Scores are mapped on the structure for subsequent examination that is useful to also recognize regions of possible local errors in protein structures. HARMONY server is located at http://caps.ncbs.res.in/harmony/

Characterization of a mercury-reducing Bacillus cereus strain isolated from the Pulicat Lake sediments, south east coast of India

Pulicat Lake sediments are often severely polluted with the toxic heavy metal mercury. Several mercury-resistant strains of Bacillus species were isolated from the sediments and all the isolates exhibited broad spectrum resistance (resistance to both organic and inorganic mercuric compounds). Plasmid curing assay showed that all the isolated Bacillus strains carry chromosomally borne mercury resistance. Polymerase chain reaction and southern hybridization analyses using merA and merB3 gene primers/probes showed that five of the isolated Bacillus strains carry sequences similar to known merA and merB3 genes. Results of multiple sequence alignment revealed 99% similarity with merA and merB3 of TnMERI1 (class II transposons). Other mercury resistant Bacillus species lacking homology to these genes were not able to volatilize mercuric chloride, indicating the presence of other modes of resistance to mercuric compounds.

Splittings of free groups, normal forms and partitions of ends

Splittings of a free group correspond to embedded spheres in the 3-manifold M = # (k) S (2) x S (1). These can be represented in a normal form due to Hatcher. In this paper, we determine the normal form in terms of crossings of partitions of ends corresponding to normal spheres, using a graph of trees representation for normal forms. In particular, we give a constructive proof of a criterion determining when a conjugacy class in pi (2)(M) can be represented by an embedded sphere.

DNA Clasping by Mycobacterial HU: The C-Terminal Region of HupB Mediates Increased Specificity of DNA Binding

Background: HU a small, basic, histone like protein is a major component of the bacterial nucleoid. E. coli has two subunits of HU coded by hupA and hupB genes whereas Mycobacterium tuberculosis (Mtb) has only one subunit of HU coded by ORF Rv2986c (hupB gene). One noticeable feature regarding Mtb HupB, based on sequence alignment of HU orthologs from different bacteria, was that HupB(Mtb) bears at its C-terminal end, a highly basic extension and this prompted an examination of its role in Mtb HupB function. Methodology/Principal Findings: With this objective two clones of Mtb HupB were generated; one expressing full length HupB protein (HupB(Mtb)) and another which expresses only the N terminal region (first 95 amino acid) of hupB (HupB(MtbN)). Gel retardation assays revealed that HupBMtbN is almost like E. coli HU (heat stable nucleoid protein) in terms of its DNA binding, with a binding constant (K-d) for linear dsDNA greater than 1000 nM, a value comparable to that obtained for the HU alpha alpha and HU alpha beta forms. However CTR (C-terminal Region) of HupB(Mtb) imparts greater specificity in DNA binding. HupB(Mtb) protein binds more strongly to supercoiled plasmid DNA than to linear DNA, also this binding is very stable as it provides DNase I protection even up to 5 minutes. Similar results were obtained when the abilities of both proteins to mediate protection against DNA strand cleavage by hydroxyl radicals generated by the Fenton's reaction, were compared. It was also observed that both the proteins have DNA binding preference for A: T rich DNA which may occur at the regulatory regions of ORFs and the oriC region of Mtb. Conclusions/Significance: These data thus point that HupB(Mtb) may participate in chromosome organization in-vivo, it may also play a passive, possibly an architectural role.

Characterization of a mercury-reducing Bacillus cereus strain isolated from the Pulicat Lake sediments, south east coast of India

Pulicat Lake sediments are often severely polluted with the toxic heavy metal mercury. Several mercury-resistant strains of Bacillus species were isolated from the sediments and all the isolates exhibited broad spectrum resistance (resistance to both organic and inorganic mercuric compounds). Plasmid curing assay showed that all the isolated Bacillus strains carry chromosomally borne mercury resistance. Polymerase chain reaction and southern hybridization analyses using merA and merB3 gene primers/probes showed that five of the isolated Bacillus strains carry sequences similar to known merA and merB3 genes. Results of multiple sequence alignment revealed 99% similarity with merA and merB3 of TnMERI1 (class II transposons). Other mercury resistant Bacillus species lacking homology to these genes were not able to volatilize mercuric chloride, indicating the presence of other modes of resistance to mercuric compounds.

Classification of Protein Kinases on the Basis of Both Kinase and Non-Kinase Regions

Background: Protein phosphorylation is a generic way to regulate signal transduction pathways in all kingdoms of life. In many organisms, it is achieved by the large family of Ser/Thr/Tyr protein kinases which are traditionally classified into groups and subfamilies on the basis of the amino acid sequence of their catalytic domains. Many protein kinases are multidomain in nature but the diversity of the accessory domains and their organization are usually not taken into account while classifying kinases into groups or subfamilies. Methodology: Here, we present an approach which considers amino acid sequences of complete gene products, in order to suggest refinements in sets of pre-classified sequences. The strategy is based on alignment-free similarity scores and iterative Area Under the Curve (AUC) computation. Similarity scores are computed by detecting common patterns between two sequences and scoring them using a substitution matrix, with a consistent normalization scheme. This allows us to handle full-length sequences, and implicitly takes into account domain diversity and domain shuffling. We quantitatively validate our approach on a subset of 212 human protein kinases. We then employ it on the complete repertoire of human protein kinases and suggest few qualitative refinements in the subfamily assignment stored in the KinG database, which is based on catalytic domains only. Based on our new measure, we delineate 37 cases of potential hybrid kinases: sequences for which classical classification based entirely on catalytic domains is inconsistent with the full-length similarity scores computed here, which implicitly consider multi-domain nature and regions outside the catalytic kinase domain. We also provide some examples of hybrid kinases of the protozoan parasite Entamoeba histolytica. Conclusions: The implicit consideration of multi-domain architectures is a valuable inclusion to complement other classification schemes. The proposed algorithm may also be employed to classify other families of enzymes with multidomain architecture.

Directional solidification of aluminum-nickel eutectic alloys using electroslag remelting

Attempts were made to produce directionally solidified, specifically grain aligned Al-6 wt pct Ni eutectic alloy using a laboratory scale ESR unit. For this purpose sand cast alloy electrodes were electroslag remelted under different mold conditions. The grain structure of the ingots obtained from these meltings showed that insulated silica molds gave the best vertical alignment of grains along the length of the ingot. The NiAl3 fibers within the grains tended to fan out and there was only a preferred alignment of fibers along the growth direction under the conditions of our experiments. The ESR parameters most suitable for vertical alignment of eutectic grains have been identified. In some electroslag remelting trials ingots were grown on a seed ingot. This resulted in a fewer vertical grains compared to the case when no seed ingot was used. The sand cast specimen of the eutectic exhibited a maximum tensile strength of around 88.2 MN/m2 (9.0 kg/mm2) whereas conventional ESR using water cooled mold gave strength value of 98.0 MN/m2 (10 kg/mm2). The directionally solidified ESR material showed longitudinal tensile strength as high as 213.7 MN/m2 (21.8 kg/mm2) which could be further increased to 220.6 MN/m2 (22.5 kg/mm2) by using the seed ingot. The average growth rate was varied between 5 to 25 mm/min during electroslag remelting in this study. The flow stresses, tangent modulus and ultimate tensile strength of directionally solidified eutectic increased with increasing growth rates.

Field-tunable stochasticity in the magnetization reversal of a cylindrical nanomagnet

The nature of magnetization reversal in an isolated cylindrical nanomagnet has been studied employing time-resolved magnetoresistance measurement. We find that the reversal mode is highly stochastic, occurring either by multimode or single-step switching. Intriguingly, the stochasticity was found to depend on the alignment of the driving magnetic field to the long axis of the nanowires, where predominantly multimode switching gives way to single-step switching behavior as the field direction is rotated from parallel to transverse with respect to the nanowire axis.

Out-of-equilibrium microrheology using optical tweezers to probe directional viscoelastic properties under shear

Many wormlike micellar systems exhibit appreciable shear thinning due to shear-induced alignment. As the micelles get aligned introducing directionality in the system, the viscoelastic properties are no longer expected to be isotropic. An optical-tweezers-based active microrheology technique enables us to probe the out-of-equilibrium rheological properties of a wormlike micellar system simultaneously along two orthogonal directions-parallel to the applied shear, as well as perpendicular to it. While the displacements of a trapped bead in response to active drag force carry signature of conventional shear thinning, its spontaneous position fluctuations along the perpendicular direction manifest an orthogonal shear thickening, an effect hitherto unobserved. Copyright (C) EPLA, 2010

A Flexible Class of Regenerating Codes for Distributed Storage

In the distributed storage setting introduced by Dimakis et al., B units of data are stored across n nodes in the network in such a way that the data can be recovered by connecting to any k nodes. Additionally one can repair a failed node by connecting to any d nodes while downloading at most beta units of data from each node. In this paper, we introduce a flexible framework in which the data can be recovered by connecting to any number of nodes as long as the total amount of data downloaded is at least B. Similarly, regeneration of a failed node is possible if the new node connects to the network using links whose individual capacity is bounded above by beta(max) and whose sum capacity equals or exceeds a predetermined parameter gamma. In this flexible setting, we obtain the cut-set lower bound on the repair bandwidth along with a constructive proof for the existence of codes meeting this bound for all values of the parameters. An explicit code construction is provided which is optimal in certain parameter regimes.

Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to a-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the ``structurally variable'' regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of `variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.

Influence of poling methods on the orientational dynamics of 2-methyl-4-nitro-aniline in poly(methyl methacrylate) studied by second harmonic generation

The magnitude and stability of the induced dipolar orientation of 2-methyl-4-nitroaniline (MNA)/poly(methyl methacrylate) (PMMA) guest/host system is investigated. The chromophores are aligned using both the corona discharge and contact electrode poling techniques. The magnitude of order parameter (also an indicator for the second order nonlinear susceptibility) is measured by recording absorbances of the poled (by the two different techniques) and unpoled PMMA films at different concentrations of MNA. Under the same conditions the corona poling technique creates a higher alignment of molecules along the field direction. The time dependence of the second harmonic intensity of the MNA/PMMA film prepared by the two techniques can be described by a Kohlrausch-Williams-Watts stretched exponential. The temperature dependence of the decay time constant is found to generally follow a modified Williams-Landel-Ferry (WLF) or Vogel-Tamann-Fulcher (VTF) equation. The glass transition temperature seems to be the single most important parameter for determining the relaxation time tau(T).

Highly Efficient and Thermally Stable Second-Order Nonlinear Optical Chromophores and Electrooptic Polymers

Organic polymeric electro-optic (E-O) materials have attracted significant attention because of their potential use as fast and efficient components of integrated photonic devices (1,2). However, the practical application of these materials in optical devices is somewhat limited by the stringent material requirements imposed by the device design, fabrication processes and operating environments. Among the various material requirements, the most notable ones are large electro-optic coefficients (r(33)) and high thermal stability (3). The design of poled polymeric materials with high electro-optic activity (r(33)) involves the optimization of the percent incorporation of efficient (large beta mu) second order nonlinear optical (NLO) chromophores into the polymer matrices and the effective creation of poling-induced non-centrosymmetric structures. The factors that affect the material stability are a) the inherent thermal stability of the NLO chromophores, b) the chemical stability of the NLO chromophores during the polymer processing conditions, and c) the long-term dipolar alignment stability at high temperatures. Although considerable progress has been made in achieving these properties (4), organic polymeric materials suitable for practical E-O device applications are yet to be developed. This chapter highlights some of our approaches in the optimization of molecular and material nonlinear optical and thermal properties.

Spin-Selective Correlation Experiment for Measurement of Long-Range J Couplings and for Assignment of (R/S) Enantiomers from the Residual Dipolar Couplings and DFT

We report the C-HETSERF experiment for determination of long- and short-range homo- and heteronuclear scalar couplings ((n)J(HH) and (n)J(XH), n >= 1) of organic molecules with a low sensitivity dilute heteronucleus in natural abundance. The method finds significant advantage in measurement of relative signs of long-range heteronuclear total couplings in chiral organic liquid crystal. The advantage of the method is demonstrated for the measurement of residual dipolar couplings (RDCs) in enantiomers oriented in the chiral liquid crystal with a focus to unambiguously assign R/S designation in a 2D spectrum. The alignment tensor calculated from the experimental RDCs and with the computed structures of enantiomers obtained by DFT calculations provides the size of the back-calculated RDCs. Smaller root-mean-square deviations (rmsd) between experimental and calculated RDCs indicate better agreement with the input structure and its correct designation of the stereogenic center.

BSSB: BLAST Server for Structural Biologists

The Basic Local Alignment Search Tool (BLAST) is one of the most widely used sequence alignment programs with which similarity searches, for both protein and nucleic acid sequences, can be performed against large databases at high speed. A large number of tools exist for processing BLAST output, but none of them provide three-dimensional structure visualization. This shortcoming has been addressed in the proposed tool BLAST Server for Structural Biologists (BSSB), which maps a BLAST output onto the three-dimensional structure of the subject protein. The three-dimensional structure of the subject protein is represented using a three-color coding scheme (identical: red; similar: yellow; and mismatch: white) based on the pairwise alignment obtained. Thus, the user will be able to visualize a possible three-dimensional structure for the query protein sequence. This information can be used to gain a deeper insight into the sequence-structure correlation. Furthermore, the additional structure-level information enables the user to make coherent and logical decisions regarding the type of input model structure or fragment that can be used for molecular replacement calculations. This tool is freely available to all users at http://bioserver1.physics.iisc.ernet.in/bssb/.