The relationship between classification of multi-domain proteins using an alignment-free approach and their functions: a case study with immunoglobulins

Establishing functional relationships between multi-domain protein sequences is a non-trivial task. Traditionally, delineating functional assignment and relationships of proteins requires domain assignments as a prerequisite. This process is sensitive to alignment quality and domain definitions. In multi-domain proteins due to multiple reasons, the quality of alignments is poor. We report the correspondence between the classification of proteins represented as full-length gene products and their functions. Our approach differs fundamentally from traditional methods in not performing the classification at the level of domains. Our method is based on an alignment free local matching scores (LMS) computation at the amino-acid sequence level followed by hierarchical clustering. As there are no gold standards for full-length protein sequence classification, we resorted to Gene Ontology and domain-architecture based similarity measures to assess our classification. The final clusters obtained using LMS show high functional and domain architectural similarities. Comparison of the current method with alignment based approaches at both domain and full-length protein showed superiority of the LMS scores. Using this method we have recreated objective relationships among different protein kinase sub-families and also classified immunoglobulin containing proteins where sub-family definitions do not exist currently. This method can be applied to any set of protein sequences and hence will be instrumental in analysis of large numbers of full-length protein sequences.

Alkyl chain substituted 1,9-pyrazoloanthrones exhibit prominent inhibitory effect on c-Jun N-terminal kinase (JNK)

N-Alkyl substituted pyrazoloanthrone derivatives were synthesized, characterized and tested for their in vitro inhibitory activity over c-Jun N-terminal kinase (JNK). Among the tested molecules, a few derivatives showed significant inhibitory activity against JNK with minimal off-target effect on other mitogen-activated protein kinase (MAP kinase) family members such as MEK1/2 and MKK3,6. These results suggested that N-alkyl (propyl and butyl) bearing pyrazoloanthrone scaffolds provide promising therapeutic inhibitors for JNK in regulating inflammation associated disorders.

The HORMA domain: an evolutionarily conserved domain discovered in chromatin-associated proteins, has unanticipated diverse functions

The HORMA domain (for Hop1p, Rev7p and MAD2) was discovered in three chromatin-associated proteins in the budding yeast Saccharomyces cerevisiae. This domain has also been found in proteins with similar functions in organisms including plants, animals and nematodes. The HORMA domain containing proteins are thought to function as adaptors for meiotic checkpoint protein signaling and in the regulation of meiotic recombination. Surprisingly, new work has disclosed completely unanticipated and diverse functions for the HORMA domain containing proteins. A. M. Villeneuve and colleagues (Schvarzstein et al., 2013) show that meiosis-specific HORMA domain containing proteins plays a vital role in preventing centriole disengagement during Caenorhabditis elegans spermatocyte meiosis. Another recent study reveals that S. cerevisiae Atg13 HORMA domain acts as a phosphorylation-dependent conformational switch in the cellular autophagic process. (C) 2014 Elsevier B.V. All rights reserved.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

Trm1p, a Zn(II)(2)Cys(6)-type transcription factor, is essential for the transcriptional activation of genes of methanol utilization pathway, in Pichia pastoris

The zinc finger transcription factors Mxr1p and Rop are key regulators of methanol metabolism in the methylotrophic yeast, Pichia pastoris, while Trm1p and Trm2p regulate methanol metabolism in Candida boidinii. Here, we demonstrate that Trm1p is essential for the expression of genes of methanol utilization (mut) pathway in P. pastoris as well. Expression of AOXI and other genes of mut pathway is severely compromised in P. pastoris Delta Trm1 strain resulting in impaired growth on media containing methanol as the sole source of carbon. Trm1p localizes to the nucleus of cells cultured on glucose or methanol. The zinc finger domain of Mxr1p but not Trm1p binds to AOXI promoter sequences in vitro, indicating that these two positive regulators act by different mechanisms. We conclude that both Trm1p and Mxr1p are essential for the expression of genes of mut pathway in P. pastoris and the mechanism of transcriptional regulation of mut pathway may be similar in P. pastoris and C. boidinii. (C) 2014 Elsevier Inc. All rights reserved.

Structure, interactions and evolutionary implications of a domain-swapped lectin dimer from Mycobacterium smegmatis

Crystal structure determination of the lectin domain of MSMEG_3662 from Mycobacterium smegmatis and its complexes with mannose and methyl-alpha-mannose, the first effort of its kind on a mycobacterial lectin, reveals a structure very similar to beta-prism II fold lectins from plant sources, but with extensive unprecedented domain swapping in dimer formation. The two subunits in a dimer often show small differences in structure, but the two domains, not always related by 2-fold symmetry, have the same structure. Each domain carries three sugar-binding sites, similar to those in plant lectins, one on each Greek key motif. The occurrence of beta-prism II fold lectins in bacteria, with characteristics similar to those from plants, indicates that this family of lectins is of ancient origin and had evolved into a mature system before bacteria and plants diverged. In plants, the number of binding sites per domain varies between one and three, whereas the number is two in the recently reported lectin domains from Pseudomonas putida and Pseudomonas aeruginosa. An analysis of the sequences of the lectins and the lectin domains shows that the level of sequence similarity among the three Greek keys in each domain has a correlation with the number of binding sites in it. Furthermore, sequence conservation among the lectins from different species is the highest for that Greek key which carries a binding site in all of them. Thus, it would appear that carbohydrate binding influences the course of the evolution of the lectin.

Grain and the concomitant ferroelectric domain size dependent physical properties of Ba0.85Ca0.15Zr0.1Ti0.9O3 ceramics fabricated using powders derived from oxalate precursor route

Fine powders comprising nanocrystallites of Ba0.85Ca0.15Zr0.1Ti0.9O3 (BCZT) were synthesized via oxalate precursor method, which facilitated to obtain homogenous and large grain sized ceramics at a lower sintering temperature. The compacted powders were sintered at various temperatures in the range of 1200 degrees C-1500 degrees C for an optimized duration of 10 h. Interestingly the one that was sintered at 1450 degrees C/10 h exhibited well resolved Morphotrophic Phase Boundary. The average grain size associated with this sample was 30 mu m accompanied by higher domain density mostly with 90 degrees twinning. These were believed to have significant contribution towards obtaining large strain of about 0.2% and piezoelectric coefficient as high as 563 pC/N. The maximum force that was generated by BCZT ceramic (having 30 mu m grain size) was found to be 161 MPa, which is much higher than that of known actuator materials such as PZT (40MPa) and NKN-5-LT (7 MPa). (C) 2014 AIP Publishing LLC.

Exact Phase Retrieval for a Class of 2-D Parametric Signals

We address the problem of two-dimensional (2-D) phase retrieval from magnitude of the Fourier spectrum. We consider 2-D signals that are characterized by first-order difference equations, which have a parametric representation in the Fourier domain. We show that, under appropriate stability conditions, such signals can be reconstructed uniquely from the Fourier transform magnitude. We formulate the phase retrieval problem as one of computing the parameters that uniquely determine the signal. We show that the problem can be solved by employing the annihilating filter method, particularly for the case when the parameters are distinct. For the more general case of the repeating parameters, the annihilating filter method is not applicable. We circumvent the problem by employing the algebraically coupled matrix pencil (ACMP) method. In the noiseless measurement setup, exact phase retrieval is possible. We also establish a link between the proposed analysis and 2-D cepstrum. In the noisy case, we derive Cramer-Rao lower bounds (CRLBs) on the estimates of the parameters and present Monte Carlo performance analysis as a function of the noise level. Comparisons with state-of-the-art techniques in terms of signal reconstruction accuracy show that the proposed technique outperforms the Fienup and relaxed averaged alternating reflections (RAAR) algorithms in the presence of noise.

Conservation Properties of the Trapezoidal Rule in Linear Time Domain Analysis of Acoustics and Structures

The trapezoidal rule, which is a special case of the Newmark family of algorithms, is one of the most widely used methods for transient hyperbolic problems. In this work, we show that this rule conserves linear and angular momenta and energy in the case of undamped linear elastodynamics problems, and an ``energy-like measure'' in the case of undamped acoustic problems. These conservation properties, thus, provide a rational basis for using this algorithm. In linear elastodynamics problems, variants of the trapezoidal rule that incorporate ``high-frequency'' dissipation are often used, since the higher frequencies, which are not approximated properly by the standard displacement-based approach, often result in unphysical behavior. Instead of modifying the trapezoidal algorithm, we propose using a hybrid finite element framework for constructing the stiffness matrix. Hybrid finite elements, which are based on a two-field variational formulation involving displacement and stresses, are known to approximate the eigenvalues much more accurately than the standard displacement-based approach, thereby either bypassing or reducing the need for high-frequency dissipation. We show this by means of several examples, where we compare the numerical solutions obtained using the displacement-based and hybrid approaches against analytical solutions.

Typical and atypical domain combinations in human protein kinases: functions, disease causing mutations and conservation in other primates

Ser/Thr and Tyr protein kinases orchestrate many signalling pathways and hence loss in this balance leads to many disease phenotypes. Due to their high abundance, diversity and importance, efforts have been made in the past to classify kinases and annotate their functions at both gross and fine levels. These kinases are conventionally classified into subfamilies based on the sequences of catalytic domains. Usually the domain architecture of a full-length kinase is consistent with the subfamily classification made based on the sequence of kinase domain. Important contributions of modular domains to the overall function of the kinase are well known. Recently occurrence of two kinds of outlier kinases-''Hybrid'' and ``Rogue'' has been reported. These show considerable deviations in their domain architectures from the typical domain architecture known for the classical kinase subfamilies. This article provides an overview of the different subfamilies of human kinases and the role of non-kinase domains in functions and diseases. Importantly this article provides analysis of hybrid and rogue kinases encoded in the human genome and highlights their conservation in closely related primate species. These kinases are examples of elegant rewiring to bring about subtle functional differences compared to canonical variants.

WNT-Inflammasome Signaling Mediates NOD2-Induced Development of Acute Arthritis in Mice

In addition to its role in innate immunity, the intracellular pathogen sensor nucleotide-binding oligomerization domain 2 (NOD2) has been implicated in various inflammatory disorders, including the development of acute arthritis. However, the molecular mechanisms involved in the development of NOD2-responsive acute arthritis are not clear. In this study, we demonstrate that NOD2 signals to a cellular protein, Ly6/PLAUR domain-containing protein 6, in a receptor-interacting protein kinase 2-TGF-beta-activated kinase 1-independent manner to activate the WNT signaling cascade. Gain- or loss-of-function of the WNT signaling pathway in an in vivo experimental mouse arthritis model or in vitro systems established the role for WNT-responsive X-linked inhibitor of apoptosis during the development of acute arthritis. Importantly, WNT-stimulated X-linked inhibitor of apoptosis mediates the activation of inflammasomes. The subsequent caspase-1 activation and IL-1 beta secretion together contribute to the phenotypic character of the inflammatory condition of acute arthritis. Thus, identification of a role for WNT-mediated inflammasome activation during NOD2 stimulation serves as a paradigm to understand NOD2-associated inflammatory disorders and develop novel therapeutics.

Novel anti IGFBP2 single chain variable fragment inhibits glioma cell migration and invasion

Insulin like growth factor binding protein 2 (IGFBP2) is highly up regulated in glioblastoma (GBM) tissues and has been one of the prognostic indicators. There are compelling evidences suggesting important roles for IGFBP2 in glioma cell proliferation, migration and invasion. Extracellular IGFBP2 through its carboxy terminal arginine glycine aspartate (RGD) motif can bind to cell surface alpha 5 beta 1 integrins and activate pathways downstream to integrin signaling. This IGFBP2 activated integrin signaling is known to play a crucial role in IGFBP2 mediated invasion of glioma cells. Hence a molecular inhibitor of carboxy terminal domain of IGFBP2 which can inhibit IGFBP2-cell surface interaction is of great therapeutic importance. In an attempt to develop molecular inhibitors of IGFBP2, we screened single chain variable fragment (scFv) phage display libraries, Tomlinson I (Library size 1.47 x 10(8)) and Tomlinson J (Library size 1.37 x 10(8)) using human recombinant IGFBP2. After screening we obtained three IGFBP2 specific binders out of which one scFv B7J showed better binding to IGFBP2 at its carboxy terminal domain, blocked IGFBP2-cell surface association, reduced activity of matrix metalloprotease 2 in the conditioned medium of glioma cells and inhibited IGFBP2 induced migration and invasion of glioma cells. We demonstrate for the first time that in vitro inhibition of extracellular IGFBP2 activity by using human scFv results in significant reduction of glioma cell migration and invasion. Therefore, the inhibition of IGFBP2 can serve as a potential therapeutic strategy in the management of GBM.

Interferon-Gamma and Nitric Oxide Synthase 2 Mediate the Aggregation of Resident Adherent Peritoneal Exudate Cells: Implications for the Host Response to Pathogens

Interferon-gamma (Ifn gamma), a key macrophage activating cytokine, plays pleiotropic roles in host immunity. In this study, the ability of Ifn gamma to induce the aggregation of resident mouse adherent peritoneal exudate cells (APECs), consisting primarily of macrophages, was investigated. Cell-cell interactions involve adhesion molecules and, upon addition of Ifn gamma, CD11b re-localizes preferentially to the sites of interaction on APECs. A functional role of CD11b in enhancing aggregation is demonstrated using Reopro, a blocking reagent, and siRNA to Cd11b. Studies with NG-methyl-L-arginine (LNMA), an inhibitor of Nitric oxide synthase (Nos), NO donors, e.g., S-nitroso-N-acetyl-DL-penicillamine (SNAP) or Diethylenetriamine/ nitric oxide adduct (DETA/NO), and Nos2(-/-) mice identified Nitric oxide (NO) induced by Ifn gamma as a key regulator of aggregation of APECs. Further studies with Nos2(-/-) APECs revealed that some Ifn. responses are independent of NO: induction of MHC class II and CD80. On the other hand, Nos2 derived NO is important for other functions: motility, phagocytosis, morphology and aggregation. Studies with cytoskeleton depolymerizing agents revealed that Ifn gamma and NO mediate the cortical stabilization of Actin and Tubulin which contribute to aggregation of APECs. The biological relevance of aggregation of APECs was delineated using infection experiments with Salmonella Typhimurium (S. Typhimurium). APECs from orally infected, but not uninfected, mice produce high amounts of NO and aggregate upon ex vivo culture in a Nos2-dependent manner. Importantly, aggregated APECs induced by Ifn gamma contain fewer intracellular S. Typhimurium compared to their single counterparts post infection. Further experiments with LNMA or Reopro revealed that both NO and CD11b are important for aggregation; in addition, NO is bactericidal. Overall, this study elucidates novel roles for Ifn gamma and Nos2 in regulating Actin, Tubulin, CD11b, motility and morphology during the aggregation response of APECs. The implications of aggregation or ``group behavior'' of APECs are discussed in the context of host resistance to infectious organisms.

Bioengineering of AAV2 Capsid at Specific Serine, Threonine, or Lysine Residues Improves Its Transduction Efficiency in Vitro and in Vivo

We hypothesized that the AAV2 vector is targeted for destruction in the cytoplasm by the host cellular kinase/ubiquitination/proteasomal machinery and that modification of their targets on AAV2 capsid may improve its transduction efficiency. In vitro analysis with pharmacological inhibitors of cellular serine/threonine kinases (protein kinase A, protein kinase C, casein kinase II) showed an increase (20-90%) on AAV2-mediated gene expression. The three-dimensional structure of AAV2 capsid was then analyzed to predict the sites of ubiquitination and phosphorylation. Three phosphodegrons, which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, were identified. Mutation targets comprising eight serine (S) or seven threonine (T) or nine lysine (K) residues were selected in and around phosphodegrons on the basis of their solvent accessibility, overlap with the receptor binding regions, overlap with interaction interfaces of capsid proteins, and their evolutionary conservation across AAV serotypes. AAV2-EGFP vectors with the wild-type (WT) capsid or mutant capsids (15 S/T -> alanine A] or 9 K -> arginine R] single mutant or 2 double K -> R mutants) were then evaluated in vitro. The transduction efficiencies of 11 S/T -> A and 7 K -> R vectors were significantly higher (similar to 63-90%) than the AAV2-WT vectors (similar to 30-40%). Further, hepatic gene transfer of these mutant vectors in vivo resulted in higher vector copy numbers (up to 4.9-fold) and transgene expression (up to 14-fold) than observed from the AAV2-WT vector. One of the mutant vectors, S489A, generated similar to 8-fold fewer antibodies that could be cross-neutralized by AAV2-WT. This study thus demonstrates the feasibility of the use of these novel AAV2 capsid mutant vectors in hepatic gene therapy.

Targeted Modifications in Adeno-Associated Virus Serotype 8 Capsid Improves Its Hepatic Gene Transfer Efficiency In Vivo

Recombinant adeno-associated virus vectors based on serotype 8 (AAV8) have shown significant promise for liver-directed gene therapy. However, to overcome the vector dose dependent immunotoxicity seen with AAV8 vectors, it is important to develop better AAV8 vectors that provide enhanced gene expression at significantly low vector doses. Since it is known that AAV vectors during intracellular trafficking are targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal machinery, we modified specific serine/threonine kinase or ubiquitination targets on the AAV8 capsid to augment its transduction efficiency. Point mutations at specific serine (S)/threonine (T)/lysine (K) residues were introduced in the AAV8 capsid at the positions equivalent to that of the effective AAV2 mutants, generated successfully earlier. Extensive structure analysis was carried out subsequently to evaluate the structural equivalence between the two serotypes. scAAV8 vectors with the wild-type (WT) and each one of the S/T -> Alanine (A) or K-Arginine (R) mutant capsids were evaluated for their liver transduction efficiency in C57BL/6 mice in vivo. Two of the AAV8-S -> A mutants (S279A and S671A), and a K137R mutant vector, demonstrated significantly higher enhanced green fluorescent protein (EGFP) transcript levels (similar to 9- to 46-fold) in the liver compared to animals that received WT-AAV8 vectors alone. The best performing AAV8 mutant (K137R) vector also had significantly reduced ubiquitination of the viral capsid, reduced activation of markers of innate immune response, and a concomitant two-fold reduction in the levels of neutralizing antibody formation in comparison to WT-AAV8 vectors. Vector bio-distribution studies revealed that the K137R mutant had a significantly higher and preferential transduction of the liver (106 vs. 7.7 vector copies/mouse diploid genome) when compared to WT-AAV8 vectors. To further study the utility of the K137R-AAV8 mutant in therapeutic gene transfer, we delivered human coagulation factor IX (h. FIX) under the control of liver-specific promoters (LP1 or hAAT) into C57BL/6 mice. The circulating levels of h. FIX: Ag were higher in all the K137R-AAV8 treated groups up to 8 weeks post-hepatic gene transfer. These studies demonstrate the feasibility of the use of this novel AAV8 vectors for potential gene therapy of hemophilia B.