Comparative kinomics of human and chimpanzee reveal unique kinship and functional diversity generated by new domain combinations

Background: Phosphorylation by protein kinases is a common event in many cellular processes. Further, many kinases perform specialized roles and are regulated by non-kinase domains tethered to kinase domain. Perturbation in the regulation of kinases leads to malignancy. We have identified and analysed putative protein kinases encoded in the genome of chimpanzee which is a close evolutionary relative of human. Result: The shared core biology between chimpanzee and human is characterized by many orthologous protein kinases which are involved in conserved pathways. Domain architectures specific to chimp/human kinases have been observed. Chimp kinases with unique domain architectures are characterized by deletion of one or more non-kinase domains in the human kinases. Interestingly, counterparts of some of the multi-domain human kinases in chimp are characterized by identical domain architectures but with kinase-like non-kinase domain. Remarkably, out of 587 chimpanzee kinases no human orthologue with greater than 95% sequence identity could be identified for 160 kinases. Variations in chimpanzee kinases compared to human kinases are brought about also by differences in functions of domains tethered to the catalytic kinase domain. For example, the heterodimer forming PB1 domain related to the fold of ubiquitin/Ras-binding domain is seen uniquely tethered to PKC-like chimpanzee kinase. Conclusion: Though the chimpanzee and human are evolutionary very close, there are chimpanzee kinases with no close counterpart in the human suggesting differences in their functions. This analysis provides a direction for experimental analysis of human and chimpanzee protein kinases in order to enhance our understanding on their specific biological roles.

Evolution of domain combinations in protein kinases and its implications for functional diversity

Protein kinases phosphorylating Ser/Thr/Tyr residues in several cellular proteins exert tight control over their biological functions. They constitute the largest protein family in most eukaryotic species. Protein kinases classified based on sequence similarity in their catalytic domains, cluster into subfamilies, which share gross functional properties. Many protein kinases are associated or tethered covalently to domains that serve as adapter or regulatory modules,naiding substrate recruitment, specificity, and also serve as scaffolds. Hence the modular organisation of the protein kinases serves as guidelines to their functional and molecular properties. Analysis of genomic repertoires of protein kinases in eukaryotes have revealed wide spectrum of domain organisation across various subfamilies of kinases. Occurrence of organism-specific novel domain combinations suggests functional diversity achieved by protein kinases in order to regulate variety of biological processes. In addition, domain architecture of protein kinases revealed existence of hybrid protein kinase subfamilies and their emerging roles in the signaling of eukaryotic organisms. In this review we discuss the repertoire of non-kinase domains tethered to multi-domain kinases in the metazoans. Similarities and differences in the domain architectures of protein kinases in these organisms indicate conserved and unique features that are critical to functional specialization. (C) 2009 Elsevier Ltd. All rights reserved.

The repertoire of protein kinases encoded in the draft version of the human genome: atypical variations and uncommon domain combinations

Background: Phosphorylation by protein kinases is central to cellular signal transduction. Abnormal functioning of kinases has been implicated in developmental disorders and malignancies. Their activity is regulated by second messengers and by the binding of associated domains, which are also influential in translocating the catalytic component to their substrate sites, in mediating interaction with other proteins and carrying out their biological roles. Results: Using sensitive profile-search methods and manual analysis, the human genome has been surveyed for protein kinases. A set of 448 sequences, which show significant similarity to protein kinases and contain the critical residues essential for kinase function, have been selected for an analysis of domain combinations after classifying the kinase domains into subfamilies. The unusual domain combinations in particular kinases suggest their involvement in ubiquitination pathways and alternative modes of regulation for mitogen-activated protein kinase kinases (MAPKKs) and cyclin-dependent kinase (CDK)-like kinases. Previously unexplored kinases have been implicated in osteoblast differentiation and embryonic development on the basis of homology with kinases of known functions from other organisms. Kinases potentially unique to vertebrates are involved in highly evolved processes such as apoptosis, protein translation and tyrosine kinase signaling. In addition to coevolution with the kinase domain, duplication and recruitment of non-catalytic domains is apparent in signaling domains such as the PH, DAG-PE, SH2 and SH3 domains. Conclusions: Expansion of the functional repertoire and possible existence of alternative modes of regulation of certain kinases is suggested by their uncommon domain combinations. Experimental verification of the predicted implications of these kinases could enhance our understanding of their biological roles.

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.

Analogous Synaptic Plasticity Profiles Emerge from Disparate Channel Combinations

An open question within the Bienenstock-Cooper-Munro theory for synaptic modification concerns the specific mechanism that is responsible for regulating the sliding modification threshold (SMT). In this conductance-based modeling study on hippocampal pyramidal neurons, we quantitatively assessed the impact of seven ion channels (R- and T-type calcium, fast sodium, delayed rectifier, A-type, and small-conductance calcium-activated (SK) potassium and HCN) and two receptors (AMPAR and NMDAR) on a calcium-dependent Bienenstock-Cooper-Munro-like plasticity rule. Our analysis with R- and T-type calcium channels revealed that differences in their activation-inactivation profiles resulted in differential impacts on how they altered the SMT. Further, we found that the impact of SK channels on the SMT critically depended on the voltage dependence and kinetics of the calcium sources with which they interacted. Next, we considered interactions among all the seven channels and the two receptors through global sensitivity analysis on 11 model parameters. We constructed 20,000 models through uniform randomization of these parameters and found 360 valid models based on experimental constraints on their plasticity profiles. Analyzing these 360 models, we found that similar plasticity profiles could emerge with several nonunique parametric combinations and that parameters exhibited weak pairwise correlations. Finally, we used seven sets of virtual knock-outs on these 360 models and found that the impact of different channels on the SMT was variable and differential. These results suggest that there are several nonunique routes to regulate the SMT, and call for a systematic analysis of the variability and state dependence of the mechanisms underlying metaplasticity during behavior and pathology.

Do carboximide-carboxylic acid combinations form co-crystals? The role of hydroxyl substitution on the formation of co-crystals and eutectics

Carboxylic acids, amides and imides are key organic systems which provide understanding of molecular recognition and binding phenomena important in biological and pharmaceutical settings. In this context, studies of their mutual interactions and compatibility through co-crystallization may pave the way for greater understanding and new applications of their combinations. Extensive co-crystallization studies are available for carboxylic acid/amide combinations, but only a few examples of carboxylic acid/imide co-crystals are currently observed in the literature. The non-formation of co-crystals for carboxylic acid/imide combinations has previously been rationalized, based on steric and computed stability factors. In the light of the growing awareness of eutectic mixtures as an alternative outcome in co-crystallization experiments, the nature of various benzoic acid/cyclic imide combinations is established in this paper. Since an additional functional group can provide sites for new intermolecular interactions and, potentially, promote supramolecular growth into a co-crystal, benzoic acids decorated with one or more hydroxyl groups have been systematically screened for co-crystallization with one unsaturated and two saturated cyclic imides. The facile formation of an abundant number of hydroxybenzoic acid/cyclic carboximide co-crystals is reported, including polymorphic and variable stoichiometry co-crystals. In the cases where co-crystals did not form, the combinations are shown invariably to result in eutectics. The presence or absence and geometric disposition of hydroxyl functionality on benzoic acid is thus found to drive the formation of co- crystals or eutectics for the studied carboxylic acid/imide combinations.

Free vibration of a kinked cantilever with attached masses

In this paper free vibration characteristics of a centrally kinked cantilever beam of unit mass carrying masses at the kink (m(k)) and at the tip (m(t)) are analyzed. Frequency factors are presented for the first two modes for different combinations of m(k),m(t) and the kink angle delta. A relationship of the form f(m(k),m(t), delta) = m(k) + m(t)(4 + 10/3 cos delta+ 2/3 cos(2) delta)=const appears to give the same fundamental frequency for a given delta and different combinations of [m(k), m(t)]. Mode shapes as well as bending moments at the support and at the kink are also discussed. The utility of a discrete beam model in understanding the free vibration characteristics is also highlighted.

Numerical prediction of load-displacement behaviors of adhesively bonded joints at different extension rates and temperatures

The present work focuses on simulation of nonlinear mechanical behaviors of adhesively bonded DLS (double lap shear) joints for variable extension rates and temperatures using the implicit ABAQUS solver. Load-displacement curves of DLS joints at nine combinations of extension rates and environmental temperatures are initially obtained by conducting tensile tests in a UTM. The joint specimens are made from dual phase (DP) steel coupons bonded with a rubber-toughened adhesive. It is shown that the shell-solid model of a DLS joint, in which substrates are modeled with shell elements and adhesive with solid elements, can effectively predict the mechanical behavior of the joint. Exponent Drucker-Prager or Von Mises yield criterion together with nonlinear isotropic hardening is used for the simulation of DLS joint tests. It has been found that at a low temperature (-20 degrees C), both Von Mises and exponent Drucker-Prager criteria give close prediction of experimental load-extension curves. However. at a high temperature (82 degrees C), Von Mises condition tends to yield a perceptibly softer joint behavior, while the corresponding response obtained using exponent Drucker-Prager criterion is much closer to the experimental load-displacement curve.

Implicit and multigrid procedures for steady-state computations with upwind algorithms

A computational study for the convergence acceleration of Euler and Navier-Stokes computations with upwind schemes has been conducted in a unified framework. It involves the flux-vector splitting algorithms due to Steger-Warming and Van Leer, the flux-difference splitting algorithms due to Roe and Osher and the hybrid algorithms, AUSM (Advection Upstream Splitting Method) and HUS (Hybrid Upwind Splitting). Implicit time integration with line Gauss-Seidel relaxation and multigrid are among the procedures which have been systematically investigated on an individual as well as cumulative basis. The upwind schemes have been tested in various implicit-explicit operator combinations such that the optimal among them can be determined based on extensive computations for two-dimensional flows in subsonic, transonic, supersonic and hypersonic flow regimes. In this study, the performance of these implicit time-integration procedures has been systematically compared with those corresponding to a multigrid accelerated explicit Runge-Kutta method. It has been demonstrated that a multigrid method employed in conjunction with an implicit time-integration scheme yields distinctly superior convergence as compared to those associated with either of the acceleration procedures provided that effective smoothers, which have been identified in this investigation, are prescribed in the implicit operator.

Simple Coding Scheme For Modular Arithmetic

A dual representation scheme for performing arithmetic modulo an arbitrary integer M is presented. The coding scheme maps each integer N in the range 0 <= N < M into one of two representations, each being identified by its most significant bit. The encoding of numbers is straightforward and the problem of checking for unused combinations is eliminated.

Base-base interactions in nucleic acids containing A-T base pairs: Structure of poly[d(A-T)]

The Watson-Crick type of base pairing is considered to be mandatory for the formation of duplex DNA. However, conformational calculations carried out in our laboratory, have shown that some combinations of backbone torsion angles and sugar pucker lead to duplexes with Hoogsteen type of base pairing also. Here we present the results of energy calculations performed on A-T containing doublet sequences in the D-form with both Hoogsteen and Watson-Crick type of base pairing and the 3 viable models for the A-T containing polynucleotide duplex poly[d(A-T)].

The credit risk(+) model with general sector correlations

We consider an enhancement of the credit risk+ model to incorporate correlations between sectors. We model the sector default rates as linear combinations of a common set of independent variables that represent macro-economic variables or risk factors. We also derive the formula for exact VaR contributions at the obligor level.

Stereochemical studies on cyclic peptides. XIII. Energy minimization studies on cyclic hexapeptides having hydrogen bonds

The basic cyclic hexapeptide conformations which accommodate hydrogen bonded β and γ turns in the backbone have been worked out using stereochemical criteria and energy minimization procedures. It was found that cyclic hexapeptides can be made up of all possible combinations of 4 ± 1 hydrogen bonded types I, I', II and II' β turns, giving rise to symmetric conformations having twofold and inversion symmetries as well as nonsymmetric structures. Conformations having exclusive features of 3 ± 1 hydrogen bonded γ turns were found to be possible in threefold and S6 symmetric cyclic hexapeptides. The results show that the cyclic hexapeptides formed by the linking of two β turn tripeptide fragments differ mainly in (a) the hydrogen bonding scheme present in the β turn tripeptides and (b) the conformation at the α-carbon atoms where the two tripeptide fragments link. The different hydrogen bonding schemes found in the component β turns are: 1) a β turn with only a 4 ± 1 hydrogen bond, 2) a type I or I' β turn with 4 ± 1 and 3 ± 1 hydrogen bonds occurring in a bifurcated form and 3) a type II or II' β turn having both the 4 ± 1 and the 3 ± 1 hydrogen bonds with the same acceptor oxygen atom. The conformation at the linking α-carbon atoms was found to lie either in the extended region or in the 3 ± 1 hydrogen bonded γ turn or inverse γ turn regions. Further, the threefold and the S6 symmetric conformations have three γ turns interleaved by three extended regions or three inverse γ turns, respectively. The feasibility of accommodating alanyl residues of both isomeric forms in the CHP minima has been explored. Finally, the available experimental data are reviewed in the light of the present results.

Dry sliding wear behaviour of magnesium alloy based hybrid composites in the longitudinal direction

In the present investigation, the wear behaviour of a creep-resistant AE42 magnesium alloy and its composites reinforced with Saffil short fibres and SiC particles in various combinations is examined in the longitudinal direction i.e., the plane containing random fibre orientation is perpendicular to the steel counter-face. Wear tests are conducted on a pin-on-disc set-up under dry sliding condition having a constant sliding velocity of 0.837 m/s for a constant sliding distance of 2.5 km in the load range of 10-40 N. It is observed that the wear rate increases with increase in load for the alloy and the composites, as expected. Wear rate of the composites is lower than the alloy and the hybrid composites exhibit a lower wear rate than the Saffil short fibres reinforced composite at all the loads. Therefore, the partial replacement of Saffil short fibres by an equal volume fraction of SiC particles not only reduces the cost but also improves the wear resistance of the composite. Microstructural investigation of the surface and subsurface of the worn pin and wear debris is carried out to explain the observed results and to understand the wear mechanisms. It is concluded that the presence of SiC particles in the hybrid composites improves the wear resistance because these particles remain intact and retain their load bearing capacity even at the highest load employed, they promote the formation of iron-rich transfer layer and they also delay the fracture of Saffil short fibres to higher loads. Under the experimental conditions used in the present investigation, the dominant wear mechanism is found to be abrasion for the AE42 alloy and its composites. It is accompanied by severe plastic deformation of surface layers in case of alloy and by the fracture of Saffil short fibres as well as the formation of iron-rich transfer layer in case of composites.

Impression creep behaviour of magnesium alloy-based hybrid composites in the transverse direction

The creep behaviour of a creep-resistant AE42 magnesium alloy reinforced with Saffil short fibres and SiC particulates in various combinations has been investigated in the transverse direction, i.e., the plane containing random fibre orientation was perpendicular to the loading direction, in the temperature range of 175-300 degrees C at the stress levels ranging from 60 to 140 MPa using impression creep test technique. Normal creep behaviour, i.e., strain rate decreasing with strain and then reaching a steady state, is observed at 175 degrees C at all the stresses employed, and up to 80 MPa stress at 240 degrees C. A reverse creep behaviour, i.e., strain rate increasing with strain, then reaching a steady state and then decreasing, is observed above 80 MPa stress at 240 degrees C and at all the stress levels at 300 degrees C. This pattern remains the same for all the composites employed. The reverse creep behaviour is found to be associated with fibre breakage. The apparent stress exponent is found to be very high for all the composites. However, after taking the threshold stress into account, the true stress exponent is found to range between 4 and 7, which suggests viscous glide and dislocation climb being the dominant creep mechanisms. The apparent activation energy Q(C) was not calculated due to insufficient data at any stress level either for normal or reverse creep behaviour. The creep resistance of the hybrid composites is found to be comparable to that of the composite reinforced with 20% Saffil short fibres alone at all the temperatures and stress levels investigated. The creep rate of the composites in the transverse direction is found to be higher than the creep rate in the longitudinal direction reported in a previous paper.