Molecular Modeling Approaches for Determining Gene Function: application to a Putative Poly-A Binding Protein from Leishmania amazonensis (LaPABP)

The great expansion in the number of genome sequencing projects has revealed the importance of computational methods to speed up the characterization of unknown genes. These studies have been improved by the use of three dimensional information from the predicted proteins generated by molecular modeling techniques. In this work, we disclose the structure-function relationship of a gene product from Leishmania amazonensis by applying molecular modeling and bioinformatics techniques. The analyzed sequence encodes a 159 aminoacids polypeptide (estimated 18 kDa) and was denoted LaPABP for its high homology with poly-A binding proteins from trypanosomatids. The domain structure, clustering analysis and a three dimensional model of LaPABP, basically obtained by homology modeling on the structure of the human poly-A binding protein, are described. Based on the analysis of the electrostatic potential mapped on the model's surface and conservation of intramolecular contacts responsible for folding stabilization we hypothesize that this protein may have less avidity to RNA than it's L. major counterpart but still account for a significant functional activity in the parasite. The model obtained will help in the design of mutagenesis experiments aimed to elucidate the mechanism of gene expression in trypanosomatids and serve as a starting point for its exploration as a potential source of targets for a rational chemotherapy.

Trypanosoma cruzi ribosomal protein S4: characterization of its coding locus, analysis of transcripts, and antigenicity of the protein

Two allelic genomic fragments containing ribosomal protein S4 encoding genes (rpS4) from Trypanosoma cruzi (CL-Brener strain) were isolated and characterized. One allele comprises two complete tandem repeats of a sequence encoding an rpS4 gene. In the other, only one rpS4 gene is found. Sequence comparison to the accessed data in the genome project database reveals that our two-copy allele corresponds to a variant haplotype. However, the deduced aminoacid sequence of all the gene copies is identical. The rpS4 transcripts processing sites were determined by comparison of genomic sequences with published cDNA data. The obtained sequence data demonstrates that rpS4 genes are expressed in epimastigotes, amastigotes, and trypomastigotes. A recombinant version of rpS4 was found to be an antigenic: it was recognized by 62.5% of the individuals with positive serology for T. cruzi and by 93.3% of patients with proven chronic chagasic disease.

Trypanosoma rangeli protein tyrosine phosphatase is associated with the parasite's flagellum

Protein tyrosine phosphatases (PTPs) play an essential role in the regulation of cell differentiation in pathogenic trypanosomatids. In this study, we describe a PTP expressed by the non-pathogenic protozoan Trypanosoma rangeli (TrPTP2). The gene for this PTP is orthologous to the T. brucei TbPTP1 and Trypanosoma cruzi (TcPTP2) genes. Cloning and expression of the TrPTP2 and TcPTP2 proteins allowed anti-PTP2 monoclonal antibodies to be generated in BALB/c mice. When expressed by T. rangeli epimastigotes and trypomastigotes, native TrPTP2 is detected as a ~65 kDa protein associated with the parasite's flagellum. Given that the flagellum is an important structure for cell differentiation in trypanosomatids, the presence of a protein responsible for tyrosine dephosphorylation in the T. rangeli flagellum could represent an interesting mechanism of regulation in this structure.

Structure-based drug design studies of the interactions ofent-kaurane diterpenes derived from Wedelia paludosa with the Plasmodium falciparumsarco/endoplasmic reticulum Ca2+-ATPase PfATP6

Malaria is responsible for more deaths around the world than any other parasitic disease. Due to the emergence of strains that are resistant to the current chemotherapeutic antimalarial arsenal, the search for new antimalarial drugs remains urgent though hampered by a lack of knowledge regarding the molecular mechanisms of artemisinin resistance. Semisynthetic compounds derived from diterpenes from the medicinal plant Wedelia paludosawere tested in silico against the Plasmodium falciparumCa2+-ATPase, PfATP6. This protein was constructed by comparative modelling using the three-dimensional structure of a homologous protein, 1IWO, as a scaffold. Compound 21 showed the best docking scores, indicating a better interaction with PfATP6 than that of thapsigargin, the natural inhibitor. Inhibition of PfATP6 by diterpene compounds could promote a change in calcium homeostasis, leading to parasite death. These data suggest PfATP6 as a potential target for the antimalarial ent-kaurane diterpenes.

Control of the phosphorylation of the astrocyte marker glial fibrillary acidic protein (GFAP) in the immature rat hippocampus by glutamate and calcium ions: possible key factor in astrocytic plasticity

The present review describes recent research on the regulation by glutamate and Ca2+ of the phosphorylation state of the intermediate filament protein of the astrocytic cytoskeleton, glial fibrillary acidic protein (GFAP), in immature hippocampal slices. The results of this research are discussed against a background of modern knowledge of the functional importance of astrocytes in the brain and of the structure and dynamic properties of intermediate filament proteins. Astrocytes are now recognized as partners with neurons in many aspects of brain function with important roles in neural plasticity. Site-specific phosphorylation of intermediate filament proteins, including GFAP, has been shown to regulate the dynamic equilibrium between the polymerized and depolymerized state of the filaments and to play a fundamental role in mitosis. Glutamate was found to increase the phosphorylation state of GFAP in hippocampal slices from rats in the post-natal age range of 12-16 days in a reaction that was dependent on external Ca2+. The lack of external Ca2+ in the absence of glutamate also increased GFAP phosphorylation to the same extent. These effects of glutamate and Ca2+ were absent in adult hippocampal slices, where the phosphorylation of GFAP was completely Ca2+-dependent. Studies using specific agonists of glutamate receptors showed that the glutamate response was mediated by a G protein-linked group II metabotropic glutamate receptor (mGluR). Since group II mGluRs do not act by liberating Ca2+ from internal stores, it is proposed that activation of the receptor by glutamate inhibits Ca2+ entry into the astrocytes and consequently down-regulates a Ca2+-dependent dephosphorylation cascade regulating the phosphorylation state of GFAP. The functional significance of these results may be related to the narrow developmental window when the glutamate response is present. In the rat brain this window corresponds to the period of massive synaptogenesis during which astrocytes are known to proliferate. Possibly, glutamate liberated from developing synapses during this period may signal an increase in the phosphorylation

Molecular biology of the kallikrein-kinin system: from structure to function

The participation of the kallikrein-kinin system, comprising the serine proteases kallikreins, the protein substrates kininogens and the effective peptides kinins, in some pathological processes like hypertension and cardiovascular diseases is still a matter of controversy. The use of different experimental set-ups in concert with the development of potent and specific inhibitors and antagonists for the system has highlighted its importance but the results still lack conclusivity. Over the last few years, transgenic and gene-targeting technologies associated with molecular biology tools have provided specific information about the elusive role of the kallikrein-kinin system in the control of blood pressure and electrolyte homeostasis. cDNA and genomic sequences for kinin receptors B2 and B1 from different species were isolated and shown to encode G-protein-coupled receptors and the structure and pharmacology of the receptors were characterized. Transgenic animals expressing an overactive kallikrein-kinin system were established to study the cardiovascular effects of these alterations and the results of these investigations further corroborate the importance of this system in the maintenance of normal blood pressure. Knockout animals for B2 and B1 receptors are available and their analysis also points to the role of these receptors in cardiovascular regulation and inflammatory processes. In this paper the most recent and relevant genetic animal models developed for the study of the kallikrein-kinin system are reviewed, and the advances they brought to the understanding of the biological role of this system are discussed.

A receptor for infectious and cellular prion protein

Prions are an unconventional form of infectious agents composed only of protein and involved in transmissible spongiform encephalopathies in humans and animals. The infectious particle is composed by PrPsc which is an isoform of a normal cellular glycosyl-phosphatidylinositol (GPI) anchored protein, PrPc, of unknown function. The two proteins differ only in conformation, PrPc is composed of 40% a helix while PrPsc has 60% ß-sheet and 20% a helix structure. The infection mechanism is trigged by interaction of PrPsc with cellular prion protein causing conversion of the latter's conformation. Therefore, the infection spreads because new PrPsc molecules are generated exponentially from the normal PrPc. The accumulation of insoluble PrPsc is probably one of the events that lead to neuronal death. Conflicting data in the literature showed that PrPc internalization is mediated either by clathrin-coated pits or by caveolae-like membranous domains. However, both pathways seem to require a third protein (a receptor or a prion-binding protein) either to make the connection between the GPI-anchored molecule to clathrin or to convert PrPc into PrPsc. We have recently characterized a 66-kDa membrane receptor which binds PrPc in vitro and in vivo and mediates the neurotoxicity of a human prion peptide. Therefore, the receptor should have a role in the pathogenesis of prion-related diseases and in the normal cellular process. Further work is necessary to clarify the events triggered by the association of PrPc/PrPsc with the receptor.

P-glycoprotein-mediated resistance to chemotherapy in cancer cells: using recombinant cytosolic domains to establish structure-function relationships

Resistance to chemotherapy in cancer cells is mainly mediated by overexpression of P-glycoprotein (Pgp), a plasma membrane ATP-binding cassette (ABC) transporter which extrudes cytotoxic drugs at the expense of ATP hydrolysis. Pgp consists of two homologous halves each containing a transmembrane domain and a cytosolic nucleotide-binding domain (NBD) which contains two consensus Walker motifs, A and B, involved in ATP binding and hydrolysis. The protein also contains an S signature characteristic of ABC transporters. The molecular mechanism of Pgp-mediated drug transport is not known. Since the transporter has an extraordinarily broad substrate specificity, its cellular function has been described as a "hydrophobic vacuum cleaner". The limited knowledge about the mechanism of Pgp, partly due to the lack of a high-resolution structure, is well reflected in the failure to efficiently inhibit its activity in cancer cells and thus to reverse multidrug resistance (MDR). In contrast to the difficulties encountered when studying the full-length Pgp, the recombinant NBDs can be obtained in large amounts as soluble proteins. The biochemical and biophysical characterization of recombinant NBDs is shown here to provide a suitable alternative route to establish structure-function relationships. NBDs were shown to bind ATP and analogues as well as potent modulators of MDR, such as hydrophobic steroids, at a region close to the ATP site. Interestingly, flavonoids also bind to NBDs with high affinity. Their binding site partly overlaps both the ATP-binding site and the steroid-interacting region. Therefore flavonoids constitute a new promising class of bifunctional modulators of Pgp.

Structure and function of the cystic fibrosis transmembrane conductance regulator

Cystic fibrosis (CF) is a lethal autosomal recessive genetic disease caused by mutations in the CF transmembrane conductance regulator (CFTR). Mutations in the CFTR gene may result in a defective processing of its protein and alter the function and regulation of this channel. Mutations are associated with different symptoms, including pancreatic insufficiency, bile duct obstruction, infertility in males, high sweat Cl-, intestinal obstruction, nasal polyp formation, chronic sinusitis, mucus dehydration, and chronic Pseudomonas aeruginosa and Staphylococcus aureus lung infection, responsible for 90% of the mortality of CF patients. The gene responsible for the cellular defect in CF was cloned in 1989 and its protein product CFTR is activated by an increase of intracellular cAMP. The CFTR contains two membrane domains, each with six transmembrane domain segments, two nucleotide-binding domains (NBDs), and a cytoplasmic domain. In this review we discuss the studies that have correlated the role of each CFTR domain in the protein function as a chloride channel and as a regulator of the outwardly rectifying Cl- channels (ORCCs).

Gap junctions in cells of the immune system: structure, regulation and possible functional roles

Gap junction channels are sites of cytoplasmic communication between contacting cells. In vertebrates, they consist of protein subunits denoted connexins (Cxs) which are encoded by a gene family. According to their Cx composition, gap junction channels show different gating and permeability properties that define which ions and small molecules permeate them. Differences in Cx primary sequences suggest that channels composed of different Cxs are regulated differentially by intracellular pathways under specific physiological conditions. Functional roles of gap junction channels could be defined by the relative importance of permeant substances, resulting in coordination of electrical and/or metabolic cellular responses. Cells of the native and specific immune systems establish transient homo- and heterocellular contacts at various steps of the immune response. Morphological and functional studies reported during the last three decades have revealed that many intercellular contacts between cells in the immune response present gap junctions or "gap junction-like" structures. Partial characterization of the molecular composition of some of these plasma membrane structures and regulatory mechanisms that control them have been published recently. Studies designed to elucidate their physiological roles suggest that they might permit coordination of cellular events which favor the effective and timely response of the immune system.

Structural, functional and immunological studies on a polymeric bacterial protein

The characterization of proteins from Brucella spp, the causative agent of brucellosis, has been the subject of intensive research. We have described an 18-kDa cytoplasmic protein of Brucella abortus and shown the potential usefulness of this protein as an antigen for the serologic diagnosis of brucellosis. The amino acid sequence of the protein showed a low but significant homology with that of lumazine synthases. Lumazine is an intermediate product in bacterial riboflavin biosynthesis. The recombinant form of the 18-kDa protein (expressed in E. coli) folds like the native Brucella protein and has lumazine-synthase enzymatic activity. Three-dimensional analysis by X-ray crystallography of the homolog Bacillus subtilis lumazine synthase has revealed that the enzyme forms an icosahedral capsid. Recombinant lumazine synthase from B. abortus was crystallized, diffracted X rays to 2.7-Å resolution at room temperature, and the structure successfully solved by molecular replacement procedures. The macromolecular assembly of the enzyme differs from that of the enzyme from B. subtilis. The Brucella enzyme remains pentameric (90 kDa) in its crystallographic form. Nonetheless, the active sites of the two enzymes are virtually identical at the structural level, indicating that inhibitors of these enzymes could be viable pharmaceuticals across a broad species range. We describe the structural reasons for the differences in their quaternary arrangement and also discuss the potential use of this protein as a target for the development of acellular vaccines.

Aggrecan structure in amphibian cartilage

The structure of the large proteoglycan present in the bullfrog epiphyseal cartilage was studied by immunochemical and biochemical methods. The isolated monomer showed a polydisperse behavior on Sepharose CL2B, with a peak at Kav = 0.14. Chondroitin sulfate chains were identified by HPLC analysis of the products formed by chondroitinase digestion and mercuric acetate treatment. These chains have approximately 38 disaccharides, a Di45:Di68 ratio of 1.6 and GalNAc4S + GalNAc4,6S are the main non-reducing terminals. Keratan sulfate was identified by the use of two monoclonal antibodies in Western blots after chondroitinase ABC treatment. A keratan sulfate-rich region (~110 kDa) was isolated by sequential treatment with chondroitinase ABC and proteases. We also employed antibodies in Western blotting experiments and showed that the full length deglycosylated core protein is about 300 kDa after SDS-PAGE. Domain-specific antibodies revealed the presence of immunoreactive sites corresponding to G1/G2 and G3 globular domains and the characterization of this large proteoglycan as aggrecan. The results indicate the high conservation of the aggrecan domain structure in this lower vertebrate.

The powerful high pressure tool for protein conformational studies

The pressure behavior of proteins may be summarized as a the pressure-induced disordering of their structures. This thermodynamic parameter has effects on proteins that are similar but not identical to those induced by temperature, the other thermodynamic parameter. Of particular importance are the intermolecular interactions that follow partial protein unfolding and that give rise to the formation of fibrils. Because some proteins do not form fibrils under pressure, these observations can be related to the shape of the stability diagram. Weak interactions which are differently affected by hydrostatic pressure or temperature play a determinant role in protein stability. Pressure acts on the 2º, 3º and 4º structures of proteins which are maintained by electrostatic and hydrophobic interactions and by hydrogen bonds. We present some typical examples of how pressure affects the tertiary structure of proteins (the case of prion proteins), induces unfolding (ataxin), is a convenient tool to study enzyme dissociation (enolase), and provides arguments to understand the role of the partial volume of an enzyme (butyrylcholinesterase). This approach may have important implications for the understanding of the basic mechanism of protein diseases and for the development of preventive and therapeutic measures.

Identification of the divergent calmodulin binding motif in yeast Ssb1/Hsp75 protein and in other HSP70 family members

Yeast soluble proteins were fractionated by calmodulin-agarose affinity chromatography and the Ca2+/calmodulin-binding proteins were analyzed by SDS-PAGE. One prominent protein of 66 kDa was excised from the gel, digested with trypsin and the masses of the resultant fragments were determined by MALDI/MS. Twenty-one of 38 monoisotopic peptide masses obtained after tryptic digestion were matched to the heat shock protein Ssb1/Hsp75, covering 37% of its sequence. Computational analysis of the primary structure of Ssb1/Hsp75 identified a unique potential amphipathic alpha-helix in its N-terminal ATPase domain with features of target regions for Ca2+/calmodulin binding. This region, which shares 89% similarity to the experimentally determined calmodulin-binding domain from mouse, Hsc70, is conserved in near half of the 113 members of the HSP70 family investigated, from yeast to plant and animals. Based on the sequence of this region, phylogenetic analysis grouped the HSP70s in three distinct branches. Two of them comprise the non-calmodulin binding Hsp70s BIP/GR78, a subfamily of eukaryotic HSP70 localized in the endoplasmic reticulum, and DnaK, a subfamily of prokaryotic HSP70. A third heterogeneous group is formed by eukaryotic cytosolic HSP70s containing the new calmodulin-binding motif and other cytosolic HSP70s whose sequences do not conform to those conserved motif, indicating that not all eukaryotic cytosolic Hsp70s are target for calmodulin regulation. Furthermore, the calmodulin-binding domain found in eukaryotic HSP70s is also the target for binding of Bag-1 - an enhancer of ADP/ATP exchange activity of Hsp70s. A model in which calmodulin displaces Bag-1 and modulates Ssb1/Hsp75 chaperone activity is discussed.

In silico analysis identifies a C3HC4-RING finger domain of a putative E3 ubiquitin-protein ligase located at the C-terminus of a polyglutamine-containing protein

Almost identical polyglutamine-containing proteins with unknown structures have been found in human, mouse and rat genomes (GenBank AJ277365, AF525300, AY879229). We infer that an identical new gene (RING) finger domain of real interest is located in each C-terminal segment. A three-dimensional (3-D) model was generated by remote homology modeling and the functional implications are discussed. The model consists of 65 residues from terminal position 707 to 772 of the human protein with a total length of 796 residues. The 3-D model predicts a ubiquitin-protein ligase (E3) as a binding site for ubiquitin-conjugating enzyme (E2). Both enzymes are part of the ubiquitin pathway to label unwanted proteins for subsequent enzymatic degradation. The molecular contact specificities are suggested for both the substrate recognition and the residues at the possible E2-binding surface. The predicted structure, of a ubiquitin-protein ligase (E3, enzyme class number 6.3.2.19, CATH code 3.30.40.10.4) may contribute to explain the process of ubiquitination. The 3-D model supports the idea of a C3HC4-RING finger with a partially new pattern. The putative E2-binding site is formed by a shallow hydrophobic groove on the surface adjacent to the helix and one zinc finger (L722, C739, P740, P741, R744). Solvent-exposed hydrophobic amino acids lie around both zinc fingers (I717, L722, F738, or P765, L766, V767, V733, P734). The 3-D structure was deposited in the protein databank theoretical model repository (2B9G, RCSB Protein Data Bank, NJ).