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.

Rapid synthesis of room temperature ferromagnetic Ag-doped LaMnO3 perovskite phases by the solution combustion method

We report the rapid solution combustion synthesis and characterization of Ag-substituted LaMnO3 phases at relatively low temperature using oxalyl dihydrazide, as fuel. Structural parameters were refined by the Rietveld method using powder X-ray diffraction data. While the parent LaMnO3 crystallizes in the orthorhombic structure, the Ag-substituted compounds crystallize in the rhombohedral symmetry. On increasing Ag-content, unit cell volume and Mn-O-Mn bond angle decreases. The Fourier transform infra red spectrum shows two absorption bands corresponding to Mn-O stretching vibration (v(s) mode) and Mn-O-Mn deformation vibration (v(b) mode) around 600 cm(-1) and 400 cm(-1) for the compositions x = 0.0, 0.05 and 0.10, respectively. Electrical resistivity measurements reveal that composition-controlled metal to insulator transition, with the maximum metal to insulator being 280 K for the composition La0.75Ag0.25MnO3. Increase in magnetic moment was observed with increase in Ag-content. The maximum magnetic moment of 35 emu/g was observed for the composition La0.80Ag0.20MnO3. (C) 2010 Elsevier Ltd. All rights reserved.

A formula for the solution of general Abel integral equation

A new formula for the solution of the general Abel Integral equation is derived, and an important special case is checked with the known result.

On the solution of the equation ut+unux+H(x, t, u)=0

We consider the equation u(t) + u(n)u(x) + H(x, t, u) = 0 and derive a transformation relating it to u(t) + u(n)u(x) = 0. Special cases of the equation appearing in applications are discussed. Initial value problems and asymptotic behaviour of the solution are studied.

On the solution of the problem of scattering of surface water waves by a sharp discontinuity in the surface boundary conditions

Closed-form analytical expressions are derived for the reflection and transmission coefficients for the problem of scattering of surface water waves by a sharp discontinuity in the surface-boundary-conditions, for the case of deep water. The method involves the use of the Havelock-type expansion of the velocity potential along with an analysis to solve a Carleman-type singular integral equation over a semi-infinite range. This method of solution is an alternative to the Wiener-Hopf technique used previously.

Structural basis for the functional and inhibitory mechanisms of beta-hydroxyacyl-acyl carrier protein dehydratase (FabZ) of Plasmodium falciparum

The beta-hydroxyacyl-acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) catalyzes the third and important reaction of the fatty acid elongation cycle. The crystal structure of PfFabZ is available in hexameric (active) and dimeric (inactive) forms. However, PfFabZ has not been crystallized with any bound inhibitors until now. We have designed a new condition to crystallize PfFabZ with its inhibitors bound in the active site, and determined the crystal structures of four of these complexes. This is the first report on any FabZ enzyme with active site inhibitors that interact directly with the catalytic residues. Inhibitor binding not only stabilized the substrate binding loop but also revealed that the substrate binding tunnel has an overall shape of ``U''. In the crystal structures, residue Phe169 located in the middle of the tunnel was found to be in two different conformations, open and closed. Thus, Phe169, merely by changing its side chain conformation, appears to be controlling the length of the tunnel to make it suitable for accommodating longer substrates. The volume of the substrate binding tunnel is determined by the sequence as well as by the conformation of the substrate binding loop region and varies between organisms for accommodating fatty acids of different chain lengths. This report on the crystal structures of the complexes of PfFabZ provides the structural basis of the inhibitory mechanism of the enzyme that could be used to improve the potency of inhibitors against an important component of fatty acid synthesis common to many infectious organisms. (C) 2011 Elsevier Inc. All rights reserved.

Phase-plane techniques for the solution of transient-stability problems

The paper presents a graphical-numerical method for determining the transient stability limits of a two-machine system under the usual assumptions of constant input, no damping and constant voltage behind transient reactance. The method presented is based on the phase-plane criterion,1, 2 in contrast to the usual step-by-step and equal-area methods. For the transient stability limit of a two-machine system, under the assumptions stated, the sum of the kinetic energy and the potential energy, at the instant of fault clearing, should just be equal to the maximum value of the potential energy which the machines can accommodate with the fault cleared. The assumption of constant voltage behind transient reactance is then discarded in favour of the more accurate assumption of constant field flux linkages. Finally, the method is extended to include the effect of field decrement and damping. A number of examples corresponding to each case are worked out, and the results obtained by the proposed method are compared with those obtained by the usual methods.

An (O)over-bar(m(2)n) randomized algorithm to compute a minimum cycle basis of a directed graph

We consider the problem of computing a minimum cycle basis in a directed graph G. The input to this problem is a directed graph whose arcs have positive weights. In this problem a {- 1, 0, 1} incidence vector is associated with each cycle and the vector space over Q generated by these vectors is the cycle space of G. A set of cycles is called a cycle basis of G if it forms a basis for its cycle space. A cycle basis where the sum of weights of the cycles is minimum is called a minimum cycle basis of G. The current fastest algorithm for computing a minimum cycle basis in a directed graph with m arcs and n vertices runs in O(m(w+1)n) time (where w < 2.376 is the exponent of matrix multiplication). If one allows randomization, then an (O) over tilde (m(3)n) algorithm is known for this problem. In this paper we present a simple (O) over tilde (m(2)n) randomized algorithm for this problem. The problem of computing a minimum cycle basis in an undirected graph has been well-studied. In this problem a {0, 1} incidence vector is associated with each cycle and the vector space over F-2 generated by these vectors is the cycle space of the graph. The fastest known algorithm for computing a minimum cycle basis in an undirected graph runs in O(m(2)n + mn(2) logn) time and our randomized algorithm for directed graphs almost matches this running time.

On the solution of bivariate population balance equations for aggregation: Pivotwise expansion of solution domain

The solution of a bivariate population balance equation (PBE) for aggregation of particles necessitates a large 2-d domain to be covered. A correspondingly large number of discretized equations for particle populations on pivots (representative sizes for bins) are solved, although at the end only a relatively small number of pivots are found to participate in the evolution process. In the present work, we initiate solution of the governing PBE on a small set of pivots that can represent the initial size distribution. New pivots are added to expand the computational domain in directions in which the evolving size distribution advances. A self-sufficient set of rules is developed to automate the addition of pivots, taken from an underlying X-grid formed by intersection of the lines of constant composition and constant particle mass. In order to test the robustness of the rule-set, simulations carried out with pivotwise expansion of X-grid are compared with those obtained using sufficiently large fixed X-grids for a number of composition independent and composition dependent aggregation kernels and initial conditions. The two techniques lead to identical predictions, with the former requiring only a fraction of the computational effort. The rule-set automatically reduces aggregation of particles of same composition to a 1-d problem. A midway change in the direction of expansion of domain, effected by the addition of particles of different mean composition, is captured correctly by the rule-set. The evolving shape of a computational domain carries with it the signature of the aggregation process, which can be insightful in complex and time dependent aggregation conditions. (c) 2012 Elsevier Ltd. All rights reserved.

Structural basis for the redox sensitivity of the Mycobacterium tuberculosis SigK-RskA sigma-anti-sigma complex

The host-pathogen interactions in Mycobacterium tuberculosis infection are significantly influenced by redox stimuli and alterations in the levels of secreted antigens. The extracyto-plasmic function (ECF) sigma factor sigma(K) governs the transcription of the serodominant antigens MPT70 and MPT83. The cellular levels of sigma(K) are regulated by the membrane-associated anti-sigma(K) (RskA) that localizes sigma(K) in an inactive complex. The crystal structure of M. tuberculosis sigma(K) in complex with the cytosolic domain of RskA (RskAcyto) revealed a disulfide bridge in the -35 promoter-interaction region of sigma(K). Biochemical experiments reveal that the redox potential of the disulfide-forming cysteines in sigma(K) is consistent with its role as a sensor. The disulfide bond in sigma(K) influences the stability of the sigma(K)-RskA(cyto) complex but does not interfere with sigma(K)-promoter DNA interactions. It is noted that these disulfide-forming cysteines are conserved across homologues, suggesting that this could be a general mechanism for redox-sensitive transcription regulation.

Molecular Basis for the Differential Quinolone Susceptibility of Mycobacterial DNA Gyrase

DNA gyrase is a type II topoisomerase that catalyzes the introduction of negative supercoils in the genomes of eubacteria. Fluoroquinolones (FQs), successful as drugs clinically, target the enzyme to trap the gyrase-DNA complex, leading to the accumulation of double-strand breaks in the genome. Mycobacteria are less susceptible to commonly used FQs. However, an 8-methoxy-substituted FQ, moxifloxacin (MFX), is a potent antimycobacterial, and a higher susceptibility of mycobacterial gyrase to MFX has been demonstrated. Although several models explain the mechanism of FQ action and gyrase-DNA-FQ interaction, the basis for the differential susceptibility of mycobacterial gyrase to various FQs is not understood. We have addressed the basis of the differential susceptibility of the gyrase and revisited the mode of action of FQs. We demonstrate that FQs bind both Escherichia coli and Mycobacterium tuberculosis gyrases in the absence of DNA and that the addition of DNA enhances the drug binding. The FQs bind primarily to the GyrA subunit of mycobacterial gyrase, while in E. coli holoenzyme is the target. The binding of MFX to GyrA of M. tuberculosis correlates with its effectiveness as a better inhibitor of the enzyme and its efficacy in cell killing.

Structural basis for the catalytic mechanism of homoserine dehydrogenase

Homoserine dehydrogenase (HSD) is an oxidoreductase in the aspartic acid pathway. This enzyme coordinates a critical branch point of the metabolic pathway that leads to the synthesis of bacterial cell-wall components such as L-lysine and m-DAP in addition to other amino acids such as L-threonine, L-methionine and L-isoleucine. Here, a structural rationale for the hydride-transfer step in the reaction mechanism of HSD is reported. The structure of Staphylococcus aureus HSD was determined at different pH conditions to understand the basis for the enhanced enzymatic activity at basic pH. An analysis of the crystal structure revealed that Lys105, which is located at the interface of the catalytic and cofactor-binding sites, could mediate the hydride-transfer step of the reaction mechanism. The role of Lys105 was subsequently confirmed by mutational analysis. Put together, these studies reveal the role of conserved water molecules and a lysine residue in hydride transfer between the substrate and the cofactor.

Prediction of Indian rainfall during the summer monsoon season on the basis of links with equatorial Pacific and Indian Ocean climate indices

Interannual variation of Indian summer monsoon rainfall (ISMR) is linked to El Nino-Southern oscillation (ENSO) as well as the Equatorial Indian Ocean oscillation (EQUINOO) with the link with the seasonal value of the ENSO index being stronger than that with the EQUINOO index. We show that the variation of a composite index determined through bivariate analysis, explains 54% of ISMR variance, suggesting a strong dependence of the skill of monsoon prediction on the skill of prediction of ENSO and EQUINOO. We explored the possibility of prediction of the Indian rainfall during the summer monsoon season on the basis of prior values of the indices. We find that such predictions are possible for July-September rainfall on the basis of June indices and for August-September rainfall based on the July indices. This will be a useful input for second and later stage forecasts made after the commencement of the monsoon season.