The period gene of the German cockroach and its novel linking power between vertebrate and invertebrate

Circadian clock protein PERIOD (PER) is essential for the endogenous clockworks in diverse lineages within Metazoa, but the protein sequences, the clock protein interactions, and the photic responses are variant and different between vertebrate and invertebrate PER homologs. Here we identified the German cockroach PER homologs and found it could bridge the huge phylogenetic gap and make possible a more precise protein sequence comparison between vertebrate and invertebrate PER homologs. Seven blocks of conserved regions (c1-c7) interspersed within PER proteins were defined, and two new significant homologies were found in the upstream portion of c3 region and in the c7 region, respectively. In addition, we found all dipteran insects PER homologs lack the c7 region and its following amino acid residues. Our results not only reveal the homology and divergence, but also imply the constraint and plasticity of divergent PER proteins during the course of evolution. These findings lay a solid foundation for understanding the general and divergent properties of circadian clockworks in diverse lineages within Metazoa.     

Improved scoring function for comparative modeling using the M4T method

Improvements in comparative protein structure modeling for the remote target-template sequence similarity cases are possible through the optimal combination of multiple template structures and by improving the quality of target-template alignment. Recently developed MMM and M4T methods were designed to address these problems. Here we describe new developments in both the alignment generation and the template selection parts of the modeling algorithms. We set up a new scoring function in MMM to deliver more accurate target-template alignments. This was achieved by developing and incorporating into the composite scoring function a novel statistical pairwise potential that combines local and non-local terms. The non-local term of the statistical potential utilizes a shuffled reference state definition that helped to eliminate most of the false positive signal from the background distribution of pairwise contacts. The accuracy of the scoring function was further increased by using BLOSUM mutation table scores.     

Potential for dramatic improvement in sequence alignment against structures of remote homologous proteins by extracting structural information from multiple structure alignment

A novel method has been developed for acquiring the correct alignment of a query sequence against remotely homologous proteins by extracting structural information from profiles of multiple structure alignment. A systematic search algorithm combined with a group of score functions based on sequence information and structural information has been introduced in this procedure. A limited number of top solutions (15,000) with high scores were selected as candidates for further examination. On a test-set comprising 301 proteins from 75 protein families with sequence identity less than 30%, the proportion of proteins with completely correct alignment as first candidate was improved to 39.8% by our method, whereas the typical performance of existing sequence-based alignment methods was only between 16.1% and 22.7%. Furthermore, multiple candidates for possible alignment were provided in our approach, which dramatically increased the possibility of finding correct alignment, such that completely correct alignments were found amongst the top-ranked 1000 candidates in 88.3% of the proteins. With the assistance of a sequence database, completely correct alignment solutions were achieved amongst the top 1000 candidates in 94.3% of the proteins. From such a limited number of candidates, it would become possible to identify more correct alignment using a more time-consuming but more powerful method with more detailed structural information, such as side-chain packing and energy minimization, etc. The results indicate that the novel alignment strategy could be helpful for extending the application of highly reliable methods for fold identification and homology modeling to a huge number of homologous proteins of low sequence similarity. Details of the methods, together with the results and implications for future development are presented.     

Identification and in silico analysis of a new group of double-histone fold-containing proteins

The double-histone fold is a rare protein fold in which two consecutive regions characterized by the typical structure of histones assemble together, thus giving a histone pseudodimer. Previously, this fold was found in a few prokaryotic histones and in the regulatory region of guanine–nucleotide exchange factors of the Sos family. Standard methods of sequence comparison did not allow us to find new proteins containing a histone pseudodimer, as previously reported (Sondermann et al. 2003). However, a deeper investigation of protein sequences showed that the two histone folds included in Sos proteins share significant sequence similarity with nucleosomal histones. On the basis of this observation, we applied a specific strategy of sequence-homology search, which led to the identification of a new group of histone pseudodimers in Cca3 and proteins similar to Cca3 (Cca3S). A homology model of the histone pseudodimer included in rat Cca3 was constructed. A subsequent structure–function relationship study revealed that the histone pseudodimers included in Cca3 and Cca3S proteins, but not those present in Sos proteins, could retain the ability of mediating protein–DNA interactions, and could consequently act as DNA-binding modules. Figure a and b Graphical representation of the statistical parameters (Pscore and Parea, see [20]) on which the prediction of DNA-binding site is based. Black crosses indicate the Pscore and Parea values calculated for 63 representative dsDNA-binding proteins, while the red asterisks refer to the values of the same parameters for Cca3 histone pseudodimer model (a), and for the amino-terminal domain of hSos1 (b). Only the proteins with Pscore > 0.12 and Parea > 250 (thus included in the upper right region of the graph) are considered dsDNA-binding proteins. c Localization of the predicted DNA-binding surface (in blue) on the rat Cca3 model     

Ligand specificity of odorant receptors

Odorant receptors belong to class A of the G protein-coupled receptors (GPCRs) and detect a large number of structurally diverse odorant molecules. A recent structural bioinformatic analysis suggests that structural features are conserved across class A of GPCRs in spite of their low sequence identity. Based on this work, we have aligned the sequences of 29 ORs for which ligand binding data are available. Recent site-directed mutagenesis experiments on one such receptor (MOR174-9) provide information that helped to identify nine amino-acid residues involved in ligand binding. Our modeling provides a rationale for amino acids in equivalent positions in most of the odorant receptors considered and helps to identify other amino acids that could be important for ligand binding. Our findings are consistent with most of the previous models and allow predictions for site-directed mutagenesis experiments, which could also validate our model.     

Sub-classification of response regulators using the surface characteristics of their receiver domains

The omnipresent bacterial switch known as a two-component system is comprised of a response regulator and a sensor kinase with which it interacts. Sensor kinases have been classified and further sub-classified into groups based on their sequence similarity, loop lengths and domain organization. Response regulators have been classified predominantly by the identity and function of their output domains. Here, comparative based homology modeling of the receiver domains of the OmpR sub-family of response regulators in Bacillus subtilis and Escherichia coli suggests further sub-classification is possible. A color-coded scale is used to show trends in surface hydrophobicity. For the OmpR receiver domains modeled these trends allow further sub-classification. The specific surface regions used for this sub-classification procedure correlate with clusters of residues that are important for interaction with cognate four helix bundle HisKA/Hpt domains.     

Improving the quality of twilight-zone alignments

Several recent publications illustrated advantages of using sequence profiles in recognizing distant homologies between proteins. At the same time, the practical usefulness of distant homology recognition depends not only on the sensitivity of the algorithm, but also on the quality of the alignment between a prediction target and the template from the database of known proteins. Here, we study this question for several supersensitive protein algorithms that were previously compared in their recognition sensitivity (Rychlewski et al., 2000). A database of protein pairs with similar structures, but low sequence similarity is used to rate the alignments obtained with several different methods, which included sequence-sequence, sequence-profile, and profile-profile alignment methods. We show that incorporation of evolutionary information encoded in sequence profiles into alignment calculation methods significantly increases the alignment accuracy, bringing them closer to the alignments obtained from structure comparison. In general, alignment quality is correlated with recognition and alignment score significance. For every alignment method, alignments with statistically significant scores correlate with both correct structural templates and good quality alignments. At the same time, average alignment lengths differ in various methods, making the comparison between them difficult. For instance, the alignments obtained by FFAS, the profile-profile alignment algorithm developed in our group are always longer that the alignments obtained with the PSI-BLAST algorithms. To address this problem, we develop methods to truncate or extend alignments to cover a specified percentage of protein lengths. In most cases, the elongation of the alignment by profile-profile methods is reasonable, adding fragments of similar structure. The examples of erroneous alignment are examined and it is shown that they can be identified based on the model quality.     

Conservation of Amino Acid Sequences Between Human and Porcine Choline Acetyltransferase

Abstract: The amino acid sequence of 11 peptides generated from human placental choline acetyltransferase was compared to the corresponding amino acid sequences predicted from the nucleotide sequence of a recently cloned porcine choline acetyltransferase cDNA. These peptides, which were generated by cyanogen bromide cleavage or tryptic digestion, accounted for 23% of the amino acids in the enzyme. Of the 145 amino acids sequenced eight differed between the two species, yielding an identity of 94% over the regions sampled.
Of the eight amino acids that differed six could represent single base changes in the DNA sequence. These findings demonstrate strong sequence similarity between porcine and human choline acetyltransferase and indicate that they are closely related evolutionarily.     

Computational analyses of JAK1 kinase domain: subtle changes in the catalytic cleft influence inhibitor specificity

The Janus kinases (JAKs) are a family of intracellular non-receptor tyrosine kinases which transmit signals by phosphorylation of downstream substrates. A myriad of cytokines can trigger the JAK-STAT pathway which influences immune response, embryonic development, and cellular transformation. Here, we built a comparative model for Jak1 based on the crystal structure of Jak2 (PDB code:2B7A) and Jak3 (PDB code: 1YVJ) using the InsightII package. 3D-Profile and stereochemical analysis further verified the validity of the proposed structure. Adenosine 5′-triphosphate (ATP) was then docked into its catalytic cleft. Although the shape of Jak1 kinase cleft is fairly similar to that of Jak3, we observed minute changes in the key residues of the binding interface which influenced the docking of a specific Jak3 inhibitor, WHI-P131. Superimposition of the interface residues suggested that substitution of Asp 99 (Jak3) into Glu 101 (Jak1) generated steric hindrance and a Tyr 91 to Phe 93 switch altered the shape of catalytic cleft which collectively prohibited the inhibitor binding. Furthermore, in-silico mutagenesis of these two residues back to Asp and Tyr enabled Jak1 to accommodate WHI-P131. These results may provide clues for the design and optimization of selective kinase inhibitors.     

Functional interactions between A' helices in the C-linker of open CNG channels

Cyclic nucleotide-gated (CNG) channels are nonselective cation channels that are activated by the direct binding of the cyclic nucleotides cAMP and cGMP. The region linking the last membrane-spanning region (S6) to the cyclic nucleotide binding domain in the COOH terminus, termed the C-linker, has been shown to play an important role in coupling cyclic nucleotide binding to opening of the pore. In this study, we explored the intersubunit proximity between the A' helices of the C-linker regions of CNGA1 in functional channels using site-specific cysteine substitution. We found that intersubunit disulfide bonds can be formed between the A' helices in open channels, and that inducing disulfide bonds in most of the studied constructs resulted in potentiation of channel activation. This suggests that the A' helices of the C-linker regions are in close proximity when the channel is in the open state. Our finding is not compatible with a homology model of the CNGA1 C-linker made from the recently published X-ray crystallographic structure of the hyperpolarization-activated, cyclic nucleotide-modulated (HCN) channel COOH terminus, and leads us to suggest that the C-linker region depicted in the crystal structure may represent the structure of the closed state. The opening conformational change would then involve a movement of the A' helices from a position parallel to the axis of the membrane to one perpendicular to the axis of the membrane.     

Predicted structures of two proteins involved in human diseases

Structures of 79 proteins involved in human diseases were predicted by sequence alignments with structural templates. The predicted structures for ALDP and CSA, proteins responsible for adrenoleukodystrophy and the Cockayne syndrome, respectively, were analyzed to elucidate the molecular basis of disease mutations. In particular we positioned residue P484 of ALDP in the homodimer interface. This positioning is consistent with a recent experimental finding that the mutation P484R significantly decreases the self-interaction of ALDP and suggests that the disease mechanism of this mutation lies in the impaired ALDP dimerization. We identified two new WD repeats in CSA and suggest that one of these forms part of the interaction surface with other proteins.     

Homology Models of the Tetramerization Domain of Six Eukaryotic Voltage-gated Potassium Channels Kv1.1-Kv1.6

Abstract

The homology models of the tetramerization (T1) domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were constructed based on the crystal structure of the Shaker T1 domain. The results of amino acid sequence alignment indicate that the T1 domains of these K⁺ channels are highly conserved, with the similarities varying from 77% between Shaker and Kv1.6 to 93% between Kv1.2 and Kv1.3. The homology models reveal that the T1 domains of these Kv channels exhibit similar folds as those of Shaker K⁺ channel. These models also show that each T1 monomer consists of three distinct layers, with N-terminal layer 1 and C- terminal layer 3 facing the cytoplasm and the membrane, respectively. Layer 2 exhibits the highest structural conservation because it is located around the central hydrophobic core. For each Kv channel, four identical subunits assemble into the homotetramer architecture around a four-fold axis through the hydrogen bonds and salt bridges formed by 15 highly conserved polar residues. The narrowest opening of the pore is formed by the four conserved residues corresponding to R115 of the Shaker T1 domain. The homology models of these Kv T1 domains provide particularly attractive targets for further structure-based studies.     

Block of sodium channels by divalent mercury: role of specific cysteinyl residues in the P-loop region

Divalent mercury (Hg(2+)) blocked human skeletal Na(+) channels (hSkM1) in a stable dose-dependent manner (K(d) = 0.96 microM) in the absence of reducing agent. Dithiothreitol (DTT) significantly prevented Hg(2+) block of hSkM1, and Hg(2+) block was also readily reversed by DTT. Both thimerosal and 2,2'-dithiodipyridine had little effect on hSkM1; however, pretreatment with thimerosal attenuated Hg(2+) block of hSkM1. Y401C+E758C rat skeletal muscle Na(+) channels (mu1) that form a disulfide bond spontaneously between two cysteines at the 401 and 758 positions showed a significantly lower sensitivity to Hg(2+) (K(d) = 18 microM). However, Y401C+E758C mu1 after reduction with DTT had a significantly higher sensitivity to Hg(2+) (K(d) = 0.36 microM) than wild-type hSkM1. Mutants C753Amu1 (K(d) = 8.47 microM) or C1521A mu1 (K(d) = 8.63 microM) exhibited significantly lower sensitivity to Hg(2+) than did wild-type hSkM1, suggesting that these two conserved cysteinyl residues of the P-loop region may play an important role in the Hg(2+) block of the hSkM1 isoform. The heart Na(+) channel (hH1) was significantly more sensitive to low-dose Hg(2+) (K(d) = 0.43 microM) than was hSkM1. The C373Y hH1 mutant exhibited higher resistance (K(d) = 1.12 microM) to Hg(2+) than did wild-type hH1. In summary, Hg(2+) probably inhibits the muscle Na(+) channels at more than one cysteinyl residue in the Na(+) channel P-loop region. Hg(2+) exhibits a lower K(d) value (<1. 23 microM) for inhibition by forming a sulfur-Hg-sulfur bridge, as compared to reaction at a single cysteinyl residue with a higher K(d) value (>8.47 microM) by forming sulfur-Hg(+) covalently. The heart Na(+) channel isoform with more than two cysteinyl residues in the P-loop region exhibits an extremely high sensitivity (K(d) < 0. 43 microM) to Hg(+), accounting for heart-specific high sensitivity to the divalent mercury.     

Characterization of hemoglobin from the lizard Uromastix hardwickii

Hemoglobin from the tropic lizard Uromastix hardwickii was isolated. Chain separations were studied, and the whole carboxymethylated globin was cleaved with trypsin. Peptides were pre-fractionated by exclusion chromatography and finally purified by reversed phase high-performance liquid chromatography. Amino acid sequence analysis permitted ordering of peptides in alpha- and beta-chains by homology with known structures in other hemoglobins. Results show large structural variations (about 50% homology between Uromastix and viper alpha-chains) and suggest chain heterogeneity with the presence of at least two types of both the alpha- and beta-chains in the preparations.     

The pore-lining region of shaker voltage-gated potassium channels: comparison of beta-barrel and alpha-helix bundle models.

Although there is a large body of site-directed mutagenesis data that identify the pore-lining sequence of the voltage-gated potassium channel, the structure of this region remains unknown. We have interpreted the available biochemical data as a set of topological and orientational restraints and employed these restraints to produce molecular models of the potassium channel pore region, H5. The H5 sequence has been modeled either as a tetramer of membrane-spanning beta-hairpins, thus producing an eight-stranded beta-barrel, or as a tetramer of incompletely membrane-spanning alpha-helical hairpins, thus producing an eight-staved alpha-helix bundle. In total, restraints-directed modeling has produced 40 different configurations of the beta-barrel model, each configuration comprising an ensemble of 20 structures, and 24 different configurations of the alpha-helix bundle model, each comprising an ensemble of 24 structures. Thus, over 1300 model structures for H5 have been generated. Configurations have been ranked on the basis of their predicted pore properties and on the extent of their agreement with the biochemical data. This ranking is employed to identify particular configurations of H5 that may be explored further as models of the pore-lining region of the voltage-gated potassium channel pore.     

Modeling the E. coli 4-hydroxybenzoic acid oligoprenyltransferase (ubiA transferase) and characterization of potential active sites

4-Hydroxybenzoate oligoprenyltransferase of E. coli, encoded in the gene ubiA, is an important key enzyme in the biosynthetic pathway to ubiquinone. It catalyzes the prenylation of 4-hydroxybenzoic acid in position 3 using an oligoprenyl diphosphate as a second substrate. Up to now, no X-ray structure of this oligoprenyltransferase or any structurally related enzyme is known. Knowledge of the tertiary structure and possible active sites is, however, essential for understanding the catalysis mechanism and the substrate specificity.With homology modeling techniques, secondary structure prediction tools, molecular dynamics simulations, and energy optimizations, a model with two putative active sites could be created and refined. One active site selected to be the most likely one for the docking of oligoprenyl diphosphate and 4-hydroxybenzoic acid is located near the N-terminus of the enzyme. It is widely accepted that residues forming an active site are usually evolutionary conserved within a family of enzymes. Multiple alignments of a multitude of related proteins clearly showed 100% conservation of the amino acid residues that form the first putative active site and therefore strongly support this hypothesis. However, an additional highly conserved region in the amino acid sequence of the ubiA enzyme could be detected, which also can be considered a putative (or rudimentary) active site. This site is characterized by a high sequence similarity to the aforementioned site and may give some hints regarding the evolutionary origin of the ubiA enzyme.Semiempirical quantum mechanical PM3 calculations have been performed to investigate the thermodynamics and kinetics of the catalysis mechanism. These results suggest a near S_N1 mechanism for the cleavage of the diphosphate ion from the isoprenyl unit. The 4-hydroxybenzoic acid interestingly appears not to be activated as benzoate anion but rather as phenolate anion to allow attack of the isoprenyl cation to the phenolate, which appeared to be the rate limiting step of the whole process according to our quantum chemical calculations. Our models are a basis for developing inhibitors of this enzyme, which is crucial for bacterial aerobic metabolism. Figure Structure of the model of ubiA oligoprenyltransferase derived from the photosynthetic reaction center (1PRC). Putative active amino acid residues and substrates are shown as capped sticks to describe their location and geometry in the putative active sites. The violet spheres identify Mg²⁺.This revised version was published online in April 2005 with corrections to Table 3 and the page make-up.     

Homology-based modeling of 3D structures of protein-protein complexes using alignments of modified sequence profiles

Customary practice in predicting 3D structures of protein-protein complexes is employment of various docking methods when the structures of separate monomers are known a priori. The alternative approach, i.e. the template-based prediction with pure sequence information as a starting point, is still considered as being inferior mostly due to presumption that the pool of available structures of protein-protein complexes, which can serve as putative templates, is not sufficiently large. Recently, however, several labs have developed databases containing thousands of 3D structures of protein-protein complexes, which enable statistically reliable testing of homology-based algorithms. In this paper we report the results on homology-based modeling of 3D structures of protein complexes using alignments of modified sequence profiles. The method, called HOMology-BAsed COmplex Prediction (HOMBACOP), has two distinctive features: (I) extra weight on aligning interfacial residues in the dynamical programming algorithm, and (II) increased gap penalties for the interfacial segments. The method was tested against our recently developed ProtCom database and against the Boston University protein-protein BENCHMARK. In both cases, models generated were compared to the models built on basis of customarily protein structure initiative (PSI)-BLAST sequence alignments. It was found that existence of homologous (by the means of PSI-BLAST) templates (44% of cases) enables both methods to produce models of good quality, with the profiles method outperforming the PSI-BLAST models (with respect to the percentage of correctly predicted residues on the complex interface and fraction of native interfacial contacts). The models were evaluated according to the CAPRI assessment criteria and about two thirds of the models were found to fall into acceptable and medium-quality categories. The same comparison of a larger set of 463 protein complexes showed again that profiles generate better models. We further demonstrate, using our ProtCom database, the suitability of the profile alignment algorithm in detecting remote homologues between query and template sequences, where the PSI-BLAST method fails.     

Hematopoietic cytokines: similarities and differences in the structures, with implications for receptor binding.

Crystal and NMR structures of helical cytokines--interleukin-4 (IL-4), granulocyte-macrophage colony-stimulating factor (GM-CSF), and interleukin-2 (IL-2)--have been compared. Root mean square deviations in the C alpha coordinates for the conserved regions of the helices were 1-2 A between different cytokines, about twice the differences observed for independently determined crystal and solution structures of IL-4. Considerable similarity in amino acid sequence in the areas expected to interact with the receptors was detected, and the available mutagenesis data for these cytokines were correlated with structure conservation. Models of cytokine-receptor interactions were postulated for IL-4 based on its structure as well as on the published structure of human growth hormone interacting with its receptors (de Vos, A.M., Ultsch, M., & Kossiakoff, A.A., 1992, Science 255, 306-312). Patches of positively charged residues on the surfaces of helices C and D of IL-4 may be responsible for the interactions with the negatively charged residues found in the complementary parts of the IL-4 receptors.     

Analysis of a three-dimensional structure of human acidic mammalian chitinase obtained by homology modeling and ligand binding studies

The three-dimensional (3D) model of the human acidic mammalian chitinase (hAMCase) was constructed based on the crystal structure of the human chitotriosidase (EC 3.2.1.44, PDB code 1HKK) by using InsightII/Homology module. With the aid of molecular mechanics and molecular dynamics methods, the last refined model was obtained and further assessed by Profile-3D and Procheck, which confirms that the refined model is reliable. Furthermore, the docking results of the ligands (allosamidin and NAG₂) into the active site of hAMCase indicate that allosamidin is a more preferred ligand than NAG₂, and that Glu119 forms hydrogen bond with allosamidin, which is in good agreement with the experimental results. From the docking studies, we also suggest that Trp10, Glu49, Asp192, and Glu276 in hAMCase are four important determinant residues in binding as they have strong van-der-Waals and electrostatic interactions with the ligand, respectively.