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.

Amino acid interaction preferences in helical membrane proteins

Membrane proteins are involved in a number of important biological functions. Yet, they are poorly understood from the structure and folding point of view. The external environment being drastically different from that of globular proteins, the intra-protein interactions in membrane proteins are also expected to be different. Hence, statistical potentials representing the features of inter-residue interactions based exclusively on the structures of membrane proteins are much needed. Currently, a reasonable number of structures are available, making it possible to undertake such an analysis on membrane proteins. In this study we have examined the inter-residue interaction propensities of amino acids in the membrane spanning regions of the alpha-helical membrane (HM) proteins. Recently we have shown that valuable information can be obtained on globular proteins by the evaluation of the pair-wise interactions of amino acids by classifying them into different structural environments, based on factors such as the secondary structure or the number of contacts that a residue can make. Here we have explored the possible ways of classifying the intra-protein environment of HM proteins and have developed scoring functions based on different classification schemes. On evaluation of different schemes, we find that the scheme which classifies amino acids to different intra-contact environment is the most promising one. Based on this classification scheme, we also redefine the hydrophobicity scale of amino acids in HM proteins.

Probing the Allosteric Mechanism in Pyrrolysyl-tRNA Synthetase Using Energy-Weighted Network Formalism

Pyrrolysyl-tRNA synthetase (PyIRS) is an atypical enzyme responsible for charging tRNA(Pyl) with pyrrolysine, despite lacking precise tRNA anticodon recognition. This dimeric protein exhibits allosteric regulation of function, like any other tRNA synthetases. In this study we examine the paths of allosteric communication at the atomic level, through energy-weighted networks of Desulfitobacterium hafniense PyIRS (DhPyIRS) and its complexes with tRNA(Pyl) and activated pyrrolysine. We performed molecular dynamics simulations of the structures of these complexes to obtain an ensemble conformation-population perspective. Weighted graph parameters relevant to identifying key players and ties in the context of social networks such as edge/node betweenness, closeness index, and the concept of funneling are explored in identifying key residues and interactions leading to shortest paths of communication in the structure networks of DhPylRS. Further, the changes in the status of important residues and connections and the costs of communication due to ligand induced perturbations are evaluated. The optimal, suboptimal, and preexisting paths are also investigated. Many of these parameters have exhibited an enhanced asymmetry between the two subunits of the dimeric protein, especially in the pretransfer complex, leading us to conclude that encoding of function goes beyond the sequence/structure of proteins. The local and global perturbations mediated by appropriate ligands and their influence on the equilibrium ensemble of conformations also have a significant role to play in the functioning of proteins. Taking a comprehensive view of these observations, we propose that the origin of many functional aspects (allostery rand half-sites reactivity in the case of DhPyIRS) lies in subtle rearrangements of interactions and dynamics at a global level.

The destabilization of human GCAP1 by a proline to leucine mutation might cause cone-rod dystrophy

Guanylate cyclase activating protein-1 (GCAP1) is required for activation of retinal guanylate cyclase-1 (RetGC1), which is essential for recovery of photoreceptor cells to the dark state. In this paper, experimentally derived observations are reported that help in explaining why a proline→leucine mutation at position 50 of human GCAP1 results in cone–rod dystrophy in a family carrying this mutation. The primary amino acid sequence of wild-type GCAP1 was mutated using site-directed mutagenesis to give a leucine at position 50. In addition, serine replaced a glutamic acid residue at position 6 to promote N‐terminal myristoylation, yielding the construct GCAP1 E6S/P50L. The enzyme was over-expressed in Escherichia coli cells, isolated and purified before being used in assays with RetGC1, characterized by circular dichroism (CD) spectroscopy, and investigated for protease resistance and thermal stability. Assays of cyclic guanosine monophosphate (cGMP) synthesis from guanosine triphosphate by RetGC1 in the presence of E6S/P50L showed that E6S/P50L could activate RetGC1 and displayed similar calcium sensitivity to wild-type GCAP1. In addition, E6S/P50L and wild-type GCAP1 possess similar CD spectra. However, there was a marked increase in the susceptibility to protease degradation and also a reduction in the thermal stability of E6S/P50L as observed by both the cGMP assay and CD spectroscopy. It is therefore suggested that although GCAP1 E6S/P50L has a similar activity and calcium dependency profile to the wild-type GCAP1, its lower stability could reduce its cellular concentration, which would in turn alter [Ca2+] and result in death of cells.

A gradient PCR-based screen for use in site-directed mutagenesis

Site-directed mutagenesis is widely used to study protein and nucleic acid structure and function. Despite recent advancements in the efficiency of procedures for site-directed mutagenesis, the fraction of site-directed mutants by most procedures rarely exceeds 50% on a routine basis and is never 100%. Hence it is typically necessary to sequence two or three clones each time a site-directed mutant is constructed. We describe a simple and robust gradient-PCR-based screen for distinguishing site-directed mutants from the starting, unmutated plasmid. The procedure can use either purified plasmid DNA or colony PCR, starting from a single colony. The screen utilizes the primer used for mutagenesis and a common outside primer that can be used for all other mutants constructed with the same template. Over 30 site-specific mutants in a variety of templates were successfully screened and all of the mutations detected were subsequently confirmed by DNA sequencing. A single base pair mismatch could be detected in an oligonucleotide of 36 bases. Detection efficiency was relatively independent of starting template concentration and the nature of the outside primer used. (C) 2003 Elsevier Science (USA). All rights reserved.

Mutagenicity associated with O6-methylguanine-DNA damage and mechanism of nucleotide flipping by AGT during repair

Methylated guanine damage at O6 position (i.e. O6MG) is dangerous due to its mutagenic and carcinogenic character that often gives rise to G:C-A:T mutation. However, the reason for this mutagenicity is not known precisely and has been a matter of controversy. Further, although it is known that O6-alkylguanine-DNA alkyltransferase (AGT) repairs O6MG paired with cytosine in DNA, the complete mechanism of target recognition and repair is not known completely. All these aspects of DNA damage and repair have been addressed here by employing high level density functional theory in gas phase and aqueous medium. It is found that the actual cause of O6MG mediated mutation may arise due to the fact that DNA polymerases incorporate thymine opposite to O6MG, misreading the resulting O6MG:T complex as an A:T base pair due to their analogous binding energies and structural alignments. It is further revealed that AGT mediated nucleotide flipping occurs in two successive steps. The intercalation of the finger residue Arg 128 into the DNA double helix and its interaction with the O6MG: C base pair followed by rotation of the O6MG nucleotide are found to be crucial for the damage recognition and nucleotide flipping.

Entrapment of a Water Wire in a Hydrophobic Peptide Channel with an Aromatic Lining

A one-dimensional water wire has been characterized by X-ray diffraction in single crystals of the tripeptide Ac-Phe-Pro-Trp-OMe. Crystals in the hexagonal space group P6(5) reveal a central hydrophobic channel lined by aromatic residues which entraps an approximately linear array of hydrogen bonded water molecules. The absence of any significant van der Waals contact with the channel walls suggests that the dominant interaction between the ``water wire'' and ``peptide nanotube'' is electrostatic in origin. An energy difference of 16 KJmol(-1) is estimated for the distinct orientations of the water wire dipole with respect to the macrodipole of the peptide nanotube. The structural model suggests that Grotthuss type proton conduction may, through constricted hydrophobic channels, be facilitated by concerted, rotational reorientation of water molecules.

AlignHUSH: Alignment of HMMs using structure and hydrophobicity information

Background: Sensitive remote homology detection and accurate alignments especially in the midnight zone of sequence similarity are needed for better function annotation and structural modeling of proteins. An algorithm, AlignHUSH for HMM-HMM alignment has been developed which is capable of recognizing distantly related domain families The method uses structural information, in the form of predicted secondary structure probabilities, and hydrophobicity of amino acids to align HMMs of two sets of aligned sequences. The effect of using adjoining column(s) information has also been investigated and is found to increase the sensitivity of HMM-HMM alignments and remote homology detection. Results: We have assessed the performance of AlignHUSH using known evolutionary relationships available in SCOP. AlignHUSH performs better than the best HMM-HMM alignment methods and is observed to be even more sensitive at higher error rates. Accuracy of the alignments obtained using AlignHUSH has been assessed using the structure-based alignments available in BaliBASE. The alignment length and the alignment quality are found to be appropriate for homology modeling and function annotation. The alignment accuracy is found to be comparable to existing methods for profile-profile alignments. Conclusions: A new method to align HMMs has been developed and is shown to have better sensitivity at error rates of 10% and above when compared to other available programs. The proposed method could effectively aid obtaining clues to functions of proteins of yet unknown function. A web-server incorporating the AlignHUSH method is available at http://crick.mbu.iisc.ernet.in/similar to alignhush/

Denaturant mediated unfolding of both native and molten globule states of maltose binding protein are accompanied by large [Delta]Cp's

Maltose binding protein (MBP) is a large, monomeric two domain protein containing 370 amino acids. In the absence of denaturant at neutral pH, the protein is in the native state, while at pH 3.0 it forms a molten globule. The molten globule lacks a tertiary circular dichroism signal but has secondary structure similar to that of the native state. The molten globule binds 8-anilino-1-naphthalene sulfonate (ANS). The unfolding thermodynamics of MBP at both pHs were measured by carrying out a series of isothermal urea melts at temperatures ranging from 274–329 K. At 298 K, values of [Delta]G°, [Delta]Cp, and Cm were 3.1 ± 0.2 kcal mol−1, 5.9 ± 0.8 kcal mol−1 K−1 (15.9 cal (mol-residue)−1 K−1), and 0.8 M, respectively, at pH 3.0 and 14.5 ± 0.4 kcal mol−1, 8.3 ± 0.7 kcal mol−1 K−1 (22.4 kcal (mol-residue)−1 K−1), and 3.3 M, respectively, at pH 7.1. Guanidine hydrochloride denaturation at pH 7.1 gave values of [Delta]G° and [Delta]Cp similar to those obtained with urea. The m values for denaturation are strongly temperature dependent, in contrast to what has been previously observed for small globular proteins. The value of [Delta]Cp per mol-residue for the molten globule is comparable to corresponding values of [Delta]Cp for the unfolding of typical globular proteins and suggests that it is a highly ordered structure, unlike molten globules of many small proteins. The value of [Delta]Cp per mol-residue for the unfolding of the native state is among the highest currently known for any protein.

Crystal structures of open and closed forms of D-serine deaminase from Salmonella typhimurium - implications on substrate specificity and catalysis

Metabolism of D-amino acids is of considerable interest due to their key importance in cell structure and function. Salmonella typhimurium D-serine deaminase (StDSD) is a pyridoxal 5' phosphate (PLP) dependent enzyme that catalyses degradation of D-Ser to pyruvate and ammonia. The first crystal structure of D-serine deaminase described here reveals a typical Foldtype II or tryptophan synthase beta subunit fold of PLP-dependent enzymes. Although holoenzyme was used for crystallization of both wild-type StDSD (WtDSD) and selenomethionine labelled StDSD (SeMetDSD), significant electron density was not observed for the cofactor, indicating that the enzyme has a low affinity for the cofactor under crystallization conditions. Interestingly, unexpected conformational differences were observed between the two structures. The WtDSD was in an open conformation while SeMetDSD, crystallized in the presence of isoserine, was in a closed conformation suggesting that the enzyme is likely to undergo conformational changes upon binding of substrate as observed in other Foldtype II PLP-dependent enzymes. Electron density corresponding to a plausible sodium ion was found near the active site of the closed but not in the open state of the enzyme. Examination of the active site and substrate modelling suggests that Thr166 may be involved in abstraction of proton from the C alpha atom of the substrate. Apart from the physiological reaction, StDSD catalyses a, b elimination of D-Thr, D-Allothr and L-Ser to the corresponding alpha-keto acids and ammonia. The structure of StDSD provides a molecular framework necessary for understanding differences in the rate of reaction with these substrates.

Cloning, expression, purification and preliminary X-ray diffraction studies of a putative Mycobacterium smegmatis thiolase

Thiolases are important in fatty-acid degradation and biosynthetic pathways. Analysis of the genomic sequence of Mycobacterium smegmatis suggests the presence of several putative thiolase genes. One of these genes appears to code for an SCP-x protein. Human SCP-x consists of an N-terminal domain (referred to as SCP2 thiolase) and a C-terminal domain (referred as sterol carrier protein 2). Here, the cloning, expression, purification and crystallization of this putative SCP-x protein from M. smegmatis are reported. The crystals diffracted X-rays to 2.5 angstrom resolution and belonged to the triclinic space group P1. Calculation of rotation functions using X-ray diffraction data suggests that the protein is likely to possess a hexameric oligomerization with 32 symmetry which has not been observed in the other six known classes of this enzyme.

Synthesis, biological studies of linear and branched arabinofuranoside-containing glycolipids and their interaction with surfactant protein A

Oligoarabinofuranoside-containing glycolipids relevant to mycobacterial cell wall components were synthesized in order to understand the functional roles of such glycolipids. A series of linear tetra-, hexa-, octa-and a branched heptasaccharide oligoarabinofuranosides, with 1 -> 2 and 1 -> 5 a-linkages between the furanoside residues, were synthesized by chemical methods from readily available monomer building blocks. Upon the synthesis of glycolipids, constituted with a double alkyl chain-substituted sn-glycerol core and oligosaccharide fragments, biological studies were performed to identify the effect of synthetic glycolipids on the biofilm formation and sliding motilities of Mycobacterium smegmatis. Synthetic glycolipids and arabinofuranosides displayed an inhibitory effect on the growth profile, but mostly on the biofilm formation and maturation. Similarly, synthetic compounds also influenced the sliding motility of the bacteria. Further, biophysical studies were undertaken, so as to identify the interactions of the glycolipids with a pulmonary surfactant protein, namely surfactant protein A (SP-A), with the aid of the surface plasmon resonance technique. Specificities of each glycolipid interacting with SP-A were thus evaluated. From this study, glycolipids were found to exhibit higher apparent association constants than the corresponding oligosaccharide portion alone, without the double alkyl group-substituted glycerol core.

Ligand-modulated Parallel Mechanical Unfolding Pathways of Maltose-binding Proteins

Protein folding and unfolding are complex phenomena, and it is accepted that multidomain proteins generally follow multiple pathways. Maltose-binding protein (MBP) is a large (a two-domain, 370-amino acid residue) bacterial periplasmic protein involved in maltose uptake. Despite the large size, it has been shown to exhibit an apparent two-state equilibrium unfolding in bulk experiments. Single-molecule studies can uncover rare events that are masked by averaging in bulk studies. Here, we use single-molecule force spectroscopy to study the mechanical unfolding pathways of MBP and its precursor protein (preMBP) in the presence and absence of ligands. Our results show that MBP exhibits kinetic partitioning on mechanical stretching and unfolds via two parallel pathways: one of them involves a mechanically stable intermediate (path I) whereas the other is devoid of it (path II). The apoMBP unfolds via path I in 62% of the mechanical unfolding events, and the remaining 38% follow path II. In the case of maltose-bound MBP, the protein unfolds via the intermediate in 79% of the cases, the remaining 21% via path II. Similarly, on binding to maltotriose, a ligand whose binding strength with the polyprotein is similar to that of maltose, the occurrence of the intermediate is comparable (82% via path I) with that of maltose. The precursor protein preMBP also shows a similar behavior upon mechanical unfolding. The percentages of molecules unfolding via path I are 53% in the apo form and 68% and 72% upon binding to maltose and maltotriose, respectively, for preMBP. These observations demonstrate that ligand binding can modulate the mechanical unfolding pathways of proteins by a kinetic partitioning mechanism. This could be a general mechanism in the unfolding of other large two-domain ligand-binding proteins of the bacterial periplasmic space.