Different amino acid substitutions at the same position in rhodopsin lead to distinct phenotypes

PURPOSE: Identification of a novel rhodopsin mutation in a family with retinitis pigmentosa and comparison of the clinical phenotype to a known mutation at the same amino acid position. METHODS: Screening for mutations in rhodopsin was performed in 78 patients with retinitis pigmentosa. All exons and flanking intronic regions were amplified by PCR, sequenced, and compared to the reference sequence derived from the National Center for Biotechnology Information (NCBI, Bethesda, MD) database. Patients were characterized clinically according to the results of best corrected visual acuity testing (BCVA), slit lamp examination (SLE), funduscopy, Goldmann perimetry (GP), dark adaptometry (DA), and electroretinography (ERG). Structural analyses of the rhodopsin protein were performed with the Swiss-Pdb Viewer program available on-line (http://www.expasy.org.spdvbv/ provided in the public domain by Swiss Institute of Bioinformatics, Geneva, Switzerland). RESULTS: A novel rhodopsin mutation (Gly90Val) was identified in a Swiss family of three generations. The pedigree indicated autosomal dominant inheritance. No additional mutation was found in this family in other autosomal dominant genes. The BCVA of affected family members ranged from 20/25 to 20/20. Fundus examination showed fine pigment mottling in patients of the third generation and well-defined bone spicules in patients of the second generation. GP showed concentric constriction. DA demonstrated monophasic cone adaptation only. ERG revealed severely reduced rod and cone signals. The clinical picture is compatible with retinitis pigmentosa. A previously reported amino acid substitution at the same position in rhodopsin leads to a phenotype resembling night blindness in mutation carriers, whereas patients reported in the current study showed the classic retinitis pigmentosa phenotype. The effect of different amino acid substitutions on the three-dimensional structure of rhodopsin was analyzed by homology modeling. Distinct distortions of position 90 (shifts in amino acids 112 and 113) and additional hydrogen bonds were found. CONCLUSIONS: Different amino acid substitutions at position 90 of rhodopsin can lead to night blindness or retinitis pigmentosa. The data suggest that the property of the substituted amino acid distinguishes between the phenotypes.

Discovery of Novel Adenosine Receptor Agonists That Exhibit Subtype Selectivity

A series of N6-bicyclic and N6-(2-hydroxy)cyclopentyl derivatives of adenosine were synthesized as novel A1R agonists and their A1R/A2R selectivity assessed using a simple yeast screening platform. We observed that the most selective, high potency ligands were achieved through N6-adamantyl substitution in combination with 5′-N-ethylcarboxamido or 5′-hydroxymethyl groups. In addition, we determined that 5′-(2-fluoro)thiophenyl derivatives all failed to generate a signaling response despite showing an interaction with the A1R. Some selected compounds were also tested on A1R and A3R in mammalian cells revealing that four of them are entirely A1R-selective agonists. By using in silico homology modeling and ligand docking, we provide insight into their mechanisms of recognition and activation of the A1R. We believe that given the broad tissue distribution, but contrasting signaling profiles, of adenosine receptor subtypes, these compounds might have therapeutic potential.

New cytochromes P450: Tale of two variants

Cytochrome P450s, a superfamily of heme enzymes found in most living organisms. They are responsible for metabolism of many therapeutic drugs, industrial pollutants, carcinogens, and additives to foodstuffs, as well as some endogenous compounds including fatty acids and steroids. First pass drug metabolism studies represent mainly liver and small intestine elimination, and are viewed as the standard to predict therapeutic outcome. However, drug plasma levels determined after administration do not always correlate with therapeutic efficacy of the drug. Therefore, a possible explanation may come by understanding drug metabolism in extrahepatic tissues and/or at the site of drug action. Identification and characterization of novel tissue specific isoforms of P450 generated by alternative splicing of known P450 genes or as yet unidentified genes is essential to predict pharmacological outcome of drugs or the fate of a carcinogen that act at sites remote from liver. ^ Using RT-PCR, brain-specific cytochrome P450s were detected in samples of human autopsy brain. So far, we have identified two human brain variants including P450 2D7 and P450 1A1. We have shown the presence of the P450 1A1 brain specific splice variant in African Americans, Caucasians and Indians albeit different patterns of liver to brain variant ratio were seen distributed throughout each population. Interestingly, the splice variant was detected only in the brain but not in any other tissues from the same individual. Homology modeling was used to compare the variant 3D structure to the liver form structure and differences in the substrate access channels and substrate binding sites were noticed. Automated computational docking was used to predict the metabolic fate of the potent carcinogenic substrate, benzo[a]pyrene. P450 1A1 brain variant showed no binding orientations that could produce the active metabolite, whereas P450 1A1 liver form did reveal orientations capable of generating active carcinogenic product. In vitro P32 labeling studies verified the docking predictions. Therefore, the data support the hypothesis that P450 brain splice variants mediate the metabolism of xenobiotics by mechanisms distinct from the well-studied liver counterparts. ^

An integrated molecular biology and computational approach to CYP1A1 expression in human brain

The cytochromes P450 comprise a superfamily of heme-containing mono-oxygenases. These enzymes metabolize numerous xenobiotics, but also play a role in metabolism of endogenous compounds. The P450 1A1 enzyme generally metabolizes polycyclic aromatic hydrocarbons, and its expression can be induced by aryl hydrocarbon receptor (AhR) activation. CYP1A1 is an exception to the generality that the majority of CYPs demonstrate highest expression in liver; CYP1Al is present in numerous extrahepatic tissues, including brain. This P450 has been observed in two forms, wildtype (WT) and brain variant (BV), arising from alternatively spliced mRNA transcripts. The CYP1A1 BV mRNA presented an exon deletion and was detected in human brain but not liver tissue of the same individuals. ^ Quantitative PCR analyses were performed to determine CYP1A1 WT and BV transcript expression levels in normal, bipolar disorder or schizophrenic groups. In our samples, we show that CYP1A1 BV mRNA, when present, is found alongside the full-length form. Furthermore, we demonstrate a significant decrease in expression of CYP1A1 in patients with bipolar disorder or schizophrenia. The expression level was not influenced by post-mortem interval, tissue pH, age, tobacco use, or lifetime antipsychotic medication load. ^ There is no indication of increased brain CYP1A1 expression in normal smokers versus non-smokers in these samples. We observed slightly increased CYP1A1 expression only in bipolar and schizophrenic smokers versus non-smokers. This may be indicative of complex interactions between neuronal chemical environments and AhR-mediated CYP1A1 induction in brain. ^ Structural homology modeling demonstrated that P450 1A1 BV has several alterations to positions/orientations of substrate recognition site residues compared to the WT isoform. Automated substrate docking was employed to investigate the potential binding of neurological signaling molecules and neurotropic drugs, as well as to differentiate specificities of the two P450 1A1 isoforms. We consistently observed that the BV isoform produced energetically favorable substrate dockings in orientations not observed for the same substrate in the WT isoform. These results demonstrated that structural differences, namely an expanded substrate access channel and active site, confer greater capacity for unique compound docking positions suggesting a metabolic profile distinct from the wildtype form for these test compounds. ^

The chitinase PfCHT1 from the human malaria parasite Plasmodium falciparum lacks proenzyme and chitin-binding domains and displays unique substrate preferences

Within hours after the ingestion of a blood meal, the mosquito midgut epithelium synthesizes a chitinous sac, the peritrophic matrix. Plasmodium ookinetes traverse the peritrophic matrix while escaping the mosquito midgut. Chitinases (EC 3.2.1.14) are critical for parasite invasion of the midgut: the presence of the chitinase inhibitor, allosamidin, in an infectious blood meal prevents oocyst development. A chitinase gene, PgCHT1, recently has been identified in the avian malaria parasite P. gallinaceum. We used the sequence of PgCHT1 to identify a P. falciparum chitinase gene, PfCHT1, in the P. falciparum genome database. PfCHT1 differs from PgCHT1 in that the P. falciparum gene lacks proenzyme and chitin-binding domains. PfCHT1 was expressed as an active recombinant enzyme in Escherichia coli. PfCHT1 shares with PgCHT1 a substrate preference unique to Plasmodium chitinases: the enzymes cleave tri- and tetramers of GlcNAc from penta- and hexameric oligomers and are unable to cleave smaller native chitin oligosaccharides. The pH activity profile of PfCHT1 and its IC50 (40 nM) to allosamidin are distinct from endochitinase activities secreted by P. gallinaceum ookinetes. Homology modeling predicts that PgCHT1 has a novel pocket in the catalytic active site that PfCHT1 lacks, which may explain the differential sensitivity of PfCHT1 and PgCHT1 to allosamidin. PfCHT1 may be the ortholog of a second, as yet unidentified, chitinase gene of P. gallinaceum. These results may allow us to develop novel strategies of blocking human malaria transmission based on interfering with P. falciparum chitinase.

RegA proteins from phage T4 and RB69 have conserved helix–loop groove RNA binding motifs but different RNA binding specificities

The RegA proteins from the bacteriophage T4 and RB69 are translational repressors that control the expression of multiple phage mRNAs. RegA proteins from the two phages share 78% sequence identity; however, in vivo expression studies have suggested that the RB69 RegA protein binds target RNAs with a higher affinity than T4 RegA protein. To study the RNA binding properties of T4 and RB69 RegA proteins more directly, the binding sites of RB69 RegA protein on synthetic RNAs corresponding to the translation initiation region of two RB69 target genes were mapped by RNase protection assays. These assays revealed that RB69 RegA protein protects nucleotides –9 to –3 (relative to the start codon) on RB69 gene 44, which contains the sequence GAAAAUU. On RB69 gene 45, the protected site (nucleotides –8 to –3) contains a similar purine-rich sequence: GAAAUA. Interestingly, T4 RegA protein protected the same nucleotides on these RNAs. To examine the specificity of RNA binding, quantitative RNA gel shift assays were performed with synthetic RNAs corresponding to recognition elements (REs) in three T4 and three RB69 mRNAs. Comparative gel shift assays demonstrated that RB69 RegA protein has an ∼7-fold higher affinity for T4 gene 44 RE RNA than T4 RegA protein. RB69 RegA protein also binds RB69 gene 44 RE RNA with a 4-fold higher affinity than T4 RegA protein. On the other hand, T4 RegA exhibited a higher affinity than RB69 RegA protein for RB69 gene 45 RE RNA. With respect to their affinities for cognate RNAs, both RegA proteins exhibited the following hierarchy of affinities: gene 44 > gene 45 > regA. Interestingly, T4 RegA exhibited the highest affinity towards RB69 gene 45 RE RNA, whereas RB69 RegA protein had the highest affinity for T4 gene 44 RE RNA. The helix–loop groove RNA binding motif of T4 RegA protein is fully conserved in RB69 RegA protein. However, homology modeling of the structure of RB69 RegA protein reveals that the divergent residues are clustered in two areas of the surface, and that there are two large areas of high conservation near the helix–loop groove, which may also play a role in RNA binding.

Demonstration of the extrinsic coagulation pathway in teleostei: Identification of zebrafish coagulation factor VII

It is not known whether the mammalian mechanism of coagulation initiation is conserved in fish. Identification of factor VII is critical in providing evidence for such a mechanism. A cDNA was cloned from a zebrafish (teleost) library that predicted a protein with sequence similarity to human factor VII. Factor VII was shown to be present in zebrafish blood and liver by Western blot analysis and immunohistochemistry. Immunodepletion of factor VII from zebrafish plasma selectively inhibited thromboplastin-triggered thrombin generation. Heterologous expression of zebrafish factor VII demonstrated a secreted protein (50 kDa) that reconstituted thromboplastin-triggered thrombin generation in immunodepleted zebrafish plasma. These results suggest conservation of the extrinsic coagulation pathway between zebrafish and humans and add credence to the zebrafish as a model for mammalian hemostasis. The structure of zebrafish factor VIIa predicted by homology modeling was consistent with the overall three-dimensional structure of human factor VIIa. However, amino acid disparities were found in the epidermal growth factor-2/serine protease regions that are present in the human tissue factor–factor VIIa contact surface, suggesting a structural basis for the species specificity of this interaction. In addition, zebrafish factor VII demonstrates that the Gla-EGF-EGF-SP domain structure, which is common to coagulation factors VII, IX, X, and protein C, was present before the radiation of the teleosts from the tetrapods. Identification of zebrafish factor VII significantly narrows the evolutionary window for development of the vertebrate coagulation cascade and provides insight into the structural basis for species specificity in the tissue factor–factor VIIa interaction.

Shared fine specificity between T-cell receptors and an antibody recognizing a peptide/major histocompatibility class I complex.

Cytotoxic T cells recognize mosaic structures consisting of target peptides embedded within self-major histocompatibility complex (MHC) class I molecules. This structure has been described in great detail for several peptide-MHC complexes. In contrast, how T-cell receptors recognize peptide-MHC complexes have been less well characterized. We have used a complete set of singly substituted analogs of a mouse MHC class I, Kk-restricted peptide, influenza hemagglutinin (Ha)255-262, to address the binding specificity of this MHC molecule. Using the same peptide-MHC complexes we determined the fine specificity of two Ha255-262-specific, Kk-restricted T cells, and of a unique antibody, pSAN, specific for the same peptide-MHC complex. Independently, a model of the Ha255-262-Kk complex was generated through homology modeling and molecular mechanics refinement. The functional data and the model corroborated each other showing that peptide residues 1, 3, 4, 6, and 7 were exposed on the MHC surface and recognized by the T cells. Thus, the majority, and perhaps all, of the side chains of the non-primary anchor residues may be available for T-cell recognition, and contribute to the stringent specificity of T cells. A striking similarity between the specificity of the T cells and that of the pSAN antibody was found and most of the peptide residues, which could be recognized by the T cells, could also be recognized by the antibody.

Repacking protein cores with backbone freedom: structure prediction for coiled coils.

Progress in homology modeling and protein design has generated considerable interest in methods for predicting side-chain packing in the hydrophobic cores of proteins. Present techniques are not practically useful, however, because they are unable to model protein main-chain flexibility. Parameterization of backbone motions may represent a general and efficient method to incorporate backbone relaxation into such fixed main-chain models. To test this notion, we introduce a method for treating explicitly the backbone motions of alpha-helical bundles based on an algebraic parameterization proposed by Francis Crick in 1953 [Crick, F. H. C. (1953) Acta Crystallogr. 6, 685-689]. Given only the core amino acid sequence, a simple calculation can rapidly reproduce the crystallographic main-chain and core side-chain structures of three coiled coils (one dimer, one trimer, and one tetramer) to within 0.6-A root-mean-square deviations. The speed of the predictive method [approximately 3 min per rotamer choice on a Silicon Graphics (Mountain View, CA) 4D/35 computer] permits it to be used as a design tool.

Comparison of the defensive elicitation activity of cerato-populin and cerato-platanin: a structural model

Plants can defend themselves from potential pathogenic microorganisms relying on a complex interplay of signaling pathways: activation of the MAPK cascade, transcription of defense related genes, production of reactive oxygen species, nitric oxide and synthesis of other defensive compounds such as phytoalexins. These events are triggered by the recognition of pathogen’s effectors (effector-triggered immunity) or PAMPs (PAMP-triggered immunity). The Cerato Platanin Family (CPF) members are Cys-rich proteins secreted and localized on fungal cell walls, involved in several aspects of fungal development and pathogen-host interactions. Although more than hundred genes of the CPF have been identified and analyzed, the structural and functional characterization of the expressed proteins has been restricted only to few members of the family. Interestingly, those proteins have been shown to bind chitin with diverse affinity and after foliar treatment they elicit defensive mechanisms in host and non-host plants. This property turns cerato platanins into interesting candidates, worth to be studied to develop new fungal elicitors with applications in sustainable agriculture. This study focus on cerato-platanin (CP), core member of the family and on the orthologous cerato-populin (Pop1). The latter shows an identity of 62% and an overall homology of 73% with respect to CP. Both proteins are able to induce MAPKs phosphorylation, production of reactive oxygen species and nitric oxide, overexpression of defense’s related genes, programmed cell death and synthesis of phytoalexins. CP, however, when compared to Pop1, induces a faster response and, in some cases, a stronger activity on plane leaves. Aim of the present research is to verify if the dissimilarities observed in the defense elicitation activity of these proteins can be associated to their structural and dynamic features. Taking advantage of the available CP NMR structure, Pop1’s 3D one was obtained by homology modeling. Experimental residual dipolar couplings and 1H, 15N, 13C resonance assignments were used to validate the model. Previous works on CPF members, addressed the highly conserved random coil regions (loops b1-b2 and b2-b3) as sufficient and necessary to induce necrosis in plants’ leaves: that region was investigated in both Pop1 and CP. In the two proteins the loops differ, in their primary sequence, for few mutations and an insertion with a consequent diversification of the proteins’ electrostatic surface. A set of 2D and 3D NMR experiments was performed to characterize both the spatial arrangement and the dynamic features of the loops. NOE data revealed a more extended network of interactions between the loops in Pop1 than in CP. In addition, in Pop1 we identified a salt bridge Lys25/Asp52 and a strong hydrophobic interaction between Phe26/Trp53. These structural features were expected not only to affect the loops’ spatial arrangement, but also to reduce the degree of their conformational freedom. Relaxation data and the order parameter S2 indeed highlighted reduced flexibility, in particular for loop b1-b2 of Pop1. In vitro NMR experiments, where Pop1 and CP were titrated with oligosaccharides, supported the hypothesis that the loops structural and dynamic differences may be responsible for the different chitin-binding properties of the two proteins: CP selectively binds tetramers of chitin in a shallow groove on one side of the barrel defined by loops b1-b2, b2-b3 and b4-b5, Pop1, instead, interacts in a non-specific fashion with oligosaccharides. Because the region involved in chitin-binding is also responsible for the defense elicitation activity, possibly being recognized by plant's receptors, it is reasonable to expect that those structural and dynamic modifications may also justify the different extent of defense elicitation. To test that hypothesis, the initial steps of a protocol aimed to the identify a receptor for CP, in silico, are presented.

Determination of alpha-conotoxin binding modes on neuronal nicotinic acetylcholine receptors

alpha-Conotoxins, from cone snails, and alpha-neurotoxins, from snakes, are competitive inhibitors of nicotinic acetylcholine receptors (nAChRs) that have overlapping binding sites in the ACh binding pocket. These disulphide-rich peptides are used extensively as tools to localize and pharmacologically characterize specific nAChRs subtypes. Recently, a homology model based on the high-resolution structure of an ACh binding protein (AChBP) allowed the three-fingered alpha-neurotoxins to be docked onto the alpha7 nAChR. To investigate if alpha-conotoxins interact with the nAChR in a similar manner, we built homology models of human alpha7 and alpha3beta2 nAChRs, and performed docking simulations of alpha-conotoxins ImI, PnIB, PnIA and MII using the program GOLD. Docking revealed that alpha-conotoxins have a different mode of interaction compared with alpha-neurotoxins, with surprisingly few nAChR residues in common between their overlapping binding sites. These docking experiments show that Imi and PnIB bind to the ACh binding pocket via a small cavity located above the beta9/beta10 hairpin of the (+)alpha7 nAChR subunit. Interestingly, PnIB, PnIA and MII were found to bind in a similar location on alpha7 or alpha3beta2 receptors mostly through hydrophobic interactions, while ImI bound further from the ACh binding pocket, mostly through electrostatic interactions. These findings, which distinguish alpha-conotoxin and alpha-neurotoxin binding modes, have implications for the rational design of selective nAChR antagonists. Copyright (C) 2004 John Wiley Sons, Ltd.

Site-directed mutagenesis and kinetic studies of the West Nile virus NS3 protease identify key enzyme-substrate interactions

The flavivirus West Nile virus (WNV) has spread rapidly throughout the world in recent years causing fever, meningitis, encephalitis, and fatalities. Because the viral protease NS2B/NS3 is essential for replication, it is attracting attention as a potential therapeutic target, although there are currently no antiviral inhibitors for any flavivirus. This paper focuses on elucidating interactions between a hexapeptide substrate (Ae-KPGLKR-p-nitroanilide) and residues at S1 and S2 in the active site of WNV protease by comparing the catalytic activities of selected mutant recombinant proteases in vitro. Homology modeling enabled the predictions of key mutations in VWNV NS3 protease at S1 (V115A/F, D129A/ E/N, S135A, Y150A/F, S160A, and S163A) and S2 (N152A) that might influence substrate recognition and catalytic efficiency. Key conclusions are that the substrate P1 Arg strongly interacts with S1 residues Asp-129, Tyr-150, and Ser-163 and, to a lesser extent, Ser-160, and P2 Lys makes an essential interaction with Asn-152 at S2. The inferred substrate-enzyme interactions provide a basis for rational protease inhibitor design and optimization. High sequence conservation within flavivirus proteases means that this study may also be relevant to design of protease inhibitors for other flavivirus proteases.

Molecular recognition of an RNA trafficking element by heterogeneous nuclear ribonucleoprotein A2

Heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multitasking protein involved in RNA packaging, alternative splicing of pre-mRNA. telomere maintenance, cytoplasmic RNA trafficking, and translation. It binds short segments of single-stranded nucleic acids, including the A2RE11 RNA element that is necessary and sufficient for cytoplasmic transport of a subset of rnRNAs in oligodendrocytes and neurons. We have explored the structures of hnRNP A2, its RNA recognition motifs (RRMs) and Gly-rich module, and the RRM complexes with A2RE11. Circular dichroism spectroscopy showed that the secondary structure of the first 189 residues of hnRNP A2 parallels that of the tandem beta alpha beta beta alpha beta RRMs of its paralogue, hnRNP A1, previously deduced from X-ray diffraction studies. The unusual GRD was shown to have substantial beta-sheet and beta-turn structure. Sedimentation equilibrium and circular dichroism results were consistent with the tandem RRM region being monomeric and supported earlier evidence for the binding of two A2RE11 oligoribonucleotides to this domain, in contrast to the protein dimer formed by the complex of hnRNP A1 with the telomeric ssDNA repeat. A three-dimensional structure for the N-terminal, two-RRM-containing segment of hnRNP A2 was derived by homology modeling. This structure was used to derive a model for the complex with A2RE11 using the previously described interaction of pairs of stacked nucleotides with aromatic residues on the RRM beta-sheet platforms, conserved in other RRM-RNA complexes, together with biochemical data and molecular dynamics-based observations of inter-RRM mobility.

In silico identification of novel G protein coupled receptors

G-protein coupled receptors (GPCRs) are a superfamily of membrane integral proteins responsible for a large number of physiological functions. Approximately 50% of marketed drugs are targeted toward a GPCR. Despite showing a high degree of structural homology, there is a large variance in sequence within the GPCR superfamily which has lead to difficulties in identifying and classifying potential new GPCR proteins. Here the various computational techniques that can be used to characterize a novel GPCR protein are discussed, including both alignment-based and alignment-free approaches. In addition, the application of homology modeling to building the three-dimensional structures of GPCRs is described.

Filogenia molecular das enzimas isocitrato liase e malato sintase e sua evolução em Viridiplanta

The plant metabolism consists of a complex network of physical and chemical events resulting in photosynthesis, respiration, synthesis and degradation of organic compounds. This is only possible due to the different kinds of responses to many environmental variations that a plant could be subject through evolution, leading also to conquering new surroundings. The glyoxylate cycle is a metabolic pathway found in glyoxysomes plant, which has unique role in the seedling establishment. Considered as a variation of the citric acid cycle, it uses an acetyl coenzyme A molecule, derived from lipids beta-oxidation to synthesize compounds which are used in carbohydrate synthesis. The Malate synthase (MLS) and Isocitrate lyase (ICL) enzyme of this cycle are unique and essential in regulating the biosynthesis of carbohydrates. Because of the absence of decarboxylation steps as rate-limiting steps, detailed studies of molecular phylogeny and evolution of these proteins enables the elucidation of the effects of this route presence in the evolutionary processes involved in their distribution across the genome from different plant species. Therefore, the aim of this study was to establish a relationship between the molecular evolution of the characteristics of enzymes from the glyoxylate cycle (isocitrate lyase and malate synthase) and their molecular phylogeny, among green plants (Viridiplantae). For this, amino acid and nucleotide sequences were used, from online repositories as UniProt and Genbank. Sequences were aligned and then subjected to an analysis of the best-fit substitution models. The phylogeny was rebuilt by distance methods (neighbor-joining) and discrete methods (maximum likelihood, maximum parsimony and Bayesian analysis). The identification of structural patterns in the evolution of the enzymes was made through homology modeling and structure prediction from protein sequences. Based on comparative analyzes of in silico models and from the results of phylogenetic inferences, both enzymes show significant structure conservation and their topologies in agreement with two processes of selection and specialization of the genes. Thus, confirming the relevance of new studies to elucidate the plant metabolism from an evolutionary perspective