The Three-dimensional Solution Structure of NaD1, a New Floral Defensin from Nicotiana alata and its Application to a Homology Model of the Crop Defense Protein alfAFP

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of NaD1, a novel antifungal and insecticidal protein isolated from the flowers of Nicotiana alata. NaD1 is a basic, cysteine-rich protein of 47 residues and is the first example of a plant defensin from flowers to be characterized structurally. Its three-dimensional structure consists of an a-helix and a triple-stranded anti-parallel beta-sheet that are stabilized by four intramolecular disulfide bonds. NaD1 features all the characteristics of the cysteine-stabilized up motif that has been described for a variety of proteins of differing functions ranging from antibacterial insect defensins and ion channel-perturbing scorpion toxins to an elicitor of the sweet taste response. The protein is biologically active against insect pests, which makes it a potential candidate for use in crop protection. NaD1 shares 31% sequence identity with alfAFP, an antifungal protein from alfalfa that confers resistance to a fungal pathogen in transgenic potatoes. The structure of NaD1 was used to obtain a homology model of alfAFP, since NaD1 has the highest level of sequence identity with alfAFP of any structurally characterized antifungal defensin. The structures of NaD1 and alfAFP were used in conjunction with structure - activity data for the radish defensin Rs-AFP2 to provide an insight into structure-function relationships. In particular, a putative effector site was identified in the structure of NaD1 and in the corresponding homology model of alfAFP. (C) 2002 Elsevier Science Ltd. All rights reserved.

Structural characterization of respiratory syncytial virus fusion inhibitor escape mutants: homology model of the F protein and a syncytium formation assay

Respiratory syncytial virus (RSV) is a ubiquitous human pathogen and the leading cause of lower respiratory tract infections in infants. Infection of cells and subsequent formation of syncytia occur through membrane fusion mediated by the RSV fusion protein (RSV-F). A novel in vitro assay of recombinant RSV-F function has been devised and used to characterize a number of escape mutants for three known inhibitors of RSV-F that have been isolated. Homology modeling of the RSV-F structure has been carried out on the basis of a chimera derived from the crystal structures of the RSV-F core and a fragment from the orthologous fusion protein from Newcastle disease virus (NDV). The structure correlates well with the appearance of RSV-F in electron micrographs, and the residues identified as contributing to specific binding sites for several monoclonal antibodies are arranged in appropriate solvent-accessible clusters. The positions of the characterized resistance mutants in the model structure identify two promising regions for the design of fusion inhibitors. (C) 2003 Elsevier Science (USA). All rights reserved.

Enzymatic characterization and homology model of a catalytically active recombinant West Nile virus NS3 protease

West Nile Virus (WNV) is a mosquito-borne flavivirus with a rapidly expanding global distribution. Infection causes severe neurological disease and fatalities in both human and animal hosts. The West Nile viral protease (NS2B-NS3) is essential for post-translational processing in host-infected cells of a viral polypeptide precursor into structural and functional viral proteins, and its inhibition could represent a potential treatment for viral infections. This article describes the design, expression, and enzymatic characterization of a catalytically active recombinant WNV protease, consisting of a 40-residue component of cofactor NS2B tethered via a noncleavable nonapeptide (G(4)SG(4)) to the N-terminal 184 residues of NS3. A chromogenic assay using synthetic para-nitroanilide (pNA) hexapeptide substrates was used to identify optimal enzyme-processing conditions (pH 9.5, I < 0.1 M, 30% glycerol, 1 mM CHAPS), preferred substrate cleavage sites, and the first competitive inhibitor (Ac-FASGKR- H, IC50 &SIM; 1 μM). A putative three-dimensional structure of WNV protease, created through homology modeling based on the crystal structures of Dengue-2 and Hepatitis C NS3 viral proteases, provides some valuable insights for structure-based design of potent and selective inhibitors of WNV protease.

Isolation and characterization of a cone snail protease with homology to CRISP proteins of the pathogenesis-related protein superfamily

The pathogenesis-related (PR) protein superfamily is widely distributed in the animal, plant, and fungal kingdoms and is implicated in human brain tumor growth and plant pathogenesis. The precise biological activity of PR proteins, however, has remained elusive. Here we report the characterization, cloning and structural homology modeling of Tex31 from the venom duct of Conus textile. Tex31 was isolated to >95% purity by activity-guided fractionation using a para-nitroanilide substrate based on the putative cleavage site residues found in the propeptide precursor of conotoxin TxVIA. Tex31 requires four residues including a leucine N-terminal of the cleavage site for efficient substrate processing. The sequence of Tex31 was determined using two degenerate PCR primers designed from N-terminal and tryptic digest Edman sequences. A BLAST search revealed that Tex31 was a member of the PR protein superfamily and most closely related to the CRISP family of mammalian proteins that have a cysteine-rich C-terminal tail. A homology model constructed from two PR proteins revealed that the likely catalytic residues in Tex31 fall within a structurally conserved domain found in PR proteins. Thus, it is possible that other PR proteins may also be substrate-specific proteases.

Molecular cloning and homology modelling of a subtilisin-like serine protease from the marine fungus, Engyodontium album BTMFS10

An alkaline protease gene (Eap) was isolated for the first time from a marine fungus, Engyodontium album. Eap consists of an open reading frame of 1,161 bp encoding a prepropeptide consisting of 387 amino acids with a calculated molecular mass of 40.923 kDa. Homology comparison of the deduced amino acid sequence of Eap with other known proteins indicated that Eap encode an extracellular protease that belongs to the subtilase family of serine protease (Family S8). A comparative homology model of the Engyodontium album protease (EAP) was developed using the crystal structure of proteinase K. The model revealed that EAP has broad substrate specificity similar to Proteinase K with preference for bulky hydrophobic residues at P1 and P4. Also, EAP is suggested to have two disulfide bonds and more than two Ca2? binding sites in its 3D structure; both of which are assumed to contribute to the thermostable nature of the protein.

Molecular cloning and homology modelling of a subtilisin-like serine protease from the marine fungus, Engyodontium album BTMFS10

An alkaline protease gene (Eap) was isolated for the first time from a marine fungus, Engyodontium album. Eap consists of an open reading frame of 1,161 bp encoding a prepropeptide consisting of 387 amino acids with a calculated molecular mass of 40.923 kDa. Homology comparison of the deduced amino acid sequence of Eap with other known proteins indicated that Eap encode an extracellular protease that belongs to the subtilase family of serine protease (Family S8). A comparative homology model of the Engyodontium album protease (EAP) was developed using the crystal structure of proteinase K. The model revealed that EAP has broad substrate specificity similar to Proteinase K with preference for bulky hydrophobic residues at P1 and P4. Also, EAP is suggested to have two disulfide bonds and more than two Ca2? binding sites in its 3D structure; both of which are assumed to contribute to the thermostable nature of the protein.

Solution structure of the human signaling protein RACK1

Background: The adaptor protein RACK1 (receptor of activated kinase 1) was originally identified as an anchoring protein for protein kinase C. RACK1 is a 36 kDa protein, and is composed of seven WD repeats which mediate its protein-protein interactions. RACK1 is ubiquitously expressed and has been implicated in diverse cellular processes involving: protein translation regulation, neuropathological processes, cellular stress, and tissue development. Results: In this study we performed a biophysical analysis of human RACK1 with the aim of obtaining low resolution structural information. Small angle X-ray scattering (SAXS) experiments demonstrated that human RACK1 is globular and monomeric in solution and its low resolution structure is strikingly similar to that of an homology model previously calculated by us and to the crystallographic structure of RACK1 isoform A from Arabidopsis thaliana. Both sedimentation velocity and sedimentation equilibrium analytical ultracentrifugation techniques showed that RACK1 is predominantly a monomer of around 37 kDa in solution, but also presents small amounts of oligomeric species. Moreover, hydrodynamic data suggested that RACK1 has a slightly asymmetric shape. The interaction of RACK1 and Ki1/57 was tested by sedimentation equilibrium. The results suggested that the association between RACK1 and Ki-1/57(122-413) follows a stoichiometry of 1:1. The binding constant (KB) observed for RACK1-Ki-1/57(122-413) interaction was of around (1.5 +/- 0.2) x 10(6) M(-1) and resulted in a dissociation constant (KD) of (0.7 +/- 0.1) x 10(-6) M. Moreover, the fluorescence data also suggests that the interaction may occur in a cooperative fashion. Conclusion: Our SAXS and analytical ultracentrifugation experiments indicated that RACK1 is predominantly a monomer in solution. RACK1 and Ki-1/57(122-413) interact strongly under the tested conditions.

Characterization of the retinoid orphan-related receptor-alpha coactivator binding interface: A structural basis for ligand-independent transcription

The retinoid orphan-related receptor-alpha (RORalpha) is a member of the ROR subfamily of orphan receptors and acts as a constitutive activator of transcription in the absence of exogenous ligands. To understand the basis of this activity, we constructed a homology model of Rill using the closely related TRbeta as a template. Molecular modeling suggested that bulky hydrophobic side chains occupy the RORa ligand cavity leaving a small but distinct cavity that may be involved in receptor stabilization. This model was subject to docking simulation with a receptor-interacting peptide from the steroid receptor coactivator, GR-interacting protein-1, which delineated a coactivator binding surface consisting of the signature motif spanning helices 3-5 and helix 12 [activation function 2 (AF2)]. Probing this surface with scanning alanine mutagenesis showed structural and functional equivalence between homologous residues of RORalpha and TRbeta. This was surprising (given that Rill is a ligand-independent activator, whereas TRbeta has an absolute requirement for ligand) and prompted us to use molecular modeling to identify differences between Rill and TRbeta in the way that the All helix interacts with the rest of the receptor. Modeling highlighted a nonconserved amino acid in helix 11 of RORa (Phe491) and a short-length of 3.10 helix at the N terminus of AF2 which we suggest i) ensures that AF2 is locked permanently in the holoconformation described for other liganded receptors and thus 2) enables ligand-independent recruitment of coactivators. Consistent with this, mutation of RORa Phe491 to either methionine or alanine (methionine is the homologous residue in TRbeta), reduced and ablated transcriptional activation and recruitment of coactivators, respectively. Furthermore, we were able to reconstitute transcriptional activity for both a deletion mutant of Ill lacking All and Phe491 Met, by overexpression of a GAL-AF2 fusion protein, demonstrating ligand-independent recruitment of AF2 and a role for Phe491 in recruiting AF2.

Targeting aquaporin function: potent inhibition of aquaglyceroporin-3 by a gold-based compound.

Aquaporins (AQPs) are membrane channels that conduct water and small solutes such as glycerol and are involved in many physiological functions. Aquaporin-based modulator drugs are predicted to be of broad potential utility in the treatment of several diseases. Until today few AQP inhibitors have been described as suitable candidates for clinical development. Here we report on the potent inhibition of AQP3 channels by gold(III) complexes screened on human red blood cells (hRBC) and AQP3-transfected PC12 cells by a stopped-flow method. Among the various metal compounds tested, Auphen is the most active on AQP3 (IC(50) = 0.8±0.08 µM in hRBC). Interestingly, the compound poorly affects the water permeability of AQP1. The mechanism of gold inhibition is related to the ability of Au(III) to interact with sulphydryls groups of proteins such as the thiolates of cysteine residues. Additional DFT and modeling studies on possible gold compound/AQP adducts provide a tentative description of the system at a molecular level. The mapping of the periplasmic surface of an homology model of human AQP3 evidenced the thiol group of Cys40 as a likely candidate for binding to gold(III) complexes. Moreover, the investigation of non-covalent binding of Au complexes by docking approaches revealed their preferential binding to AQP3 with respect to AQP1. The high selectivity and low concentration dependent inhibitory effect of Auphen (in the nanomolar range) together with its high water solubility makes the compound a suitable drug lead for future in vivo studies. These results may present novel metal-based scaffolds for AQP drug development.

Megazol and its bioisostere 4H-1,2,4-triazole: comparing the trypanocidal, cytotoxic and genotoxic activities and their in vitro and in silico interactions with the Trypanosoma brucei nitroreductase enzyme

Megazol (7) is a 5-nitroimidazole that is highly active against Trypanosoma cruzi and Trypanosoma brucei, as well as drug-resistant forms of trypanosomiasis. Compound 7 is not used clinically due to its mutagenic and genotoxic properties, but has been largely used as a lead compound. Here, we compared the activity of 7 with its 4H-1,2,4-triazole bioisostere (8) in bloodstream forms of T. brucei and T. cruzi and evaluated their activation by T. brucei type I nitroreductase (TbNTR) enzyme. We also analysed the cytotoxic and genotoxic effects of these compounds in whole human blood using Comet and fluorescein diacetate/ethidium bromide assays. Although the only difference between 7 and 8 is the substitution of sulphur (in the thiadiazole in 7) for nitrogen (in the triazole in 8), the results indicated that 8 had poorer antiparasitic activity than 7 and was not genotoxic, whereas 7 presented this effect. The determination of Vmax indicated that although 8 was metabolised more rapidly than 7, it bounds to the TbNTR with better affinity, resulting in equivalent kcat/KM values. Docking assays of 7 and 8 performed within the active site of a homology model of the TbNTR indicating that 8 had greater affinity than 7.

Modeling of the TCR-MHC-peptide complex.

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling.

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

Computational toolbox towards evolutionary domain mapping of membrane proteins

Membrane proteins account for about 20% to 30% of all proteins encoded in a typical genome. They play central roles in multiple cellular processes mediating the interaction of the cell with its surrounding. Over 60% of all drug targets contain a membrane domain. The experimental difficulties of obtaining a crystal structural severely limits our ability or understanding of membrane protein function. Computational evolutionary studies of proteins are crucial for the prediction of 3D structures. In this project, we construct a tool able to quantify the evolutionary positive selective pressure on each residue of membrane proteins through maximum likelihood phylogeny reconstruction. The conservation plot combined with a structural homology model is also a potent tool to predict those residues that have essentials roles in the structure and function of a membrane protein and can be very useful in the design of validation experiments.

Prediction of 3D structure and ligand-binding properties : ABC transporters and flavodiiron proteins from Synechocystis sp. strain PCC 6803

Cyanobacteria are unicellular, non-nitrogen-fixing prokaryotes, which perform photosynthesis similarly as higher plants. The cyanobacterium Synechocystis sp. strain PCC 6803 is used as a model organism in photosynthesis research. My research described herein aims at understanding the function of the photosynthetic machinery and how it responds to changes in the environment. Detailed knowledge of the regulation of photosynthesis in cyanobacteria can be utilized for biotechnological purposes, for example in the harnessing of solar energy for biofuel production. In photosynthesis, iron participates in electron transfer. Here, we focused on iron transport in Synechocystis sp. strain PCC 6803 and particularly on the environmental regulation of the genes encoding the FutA2BC ferric iron transporter, which belongs to the ABC transporter family. A homology model built for the ATP-binding subunit FutC indicates that it has a functional ATPbinding site as well as conserved interactions with the channel-forming subunit FutB in the transporter complex. Polyamines are important for the cell proliferation, differentiation and apoptosis in prokaryotic and eukaryotic cells. In plants, polyamines have special roles in stress response and in plant survival. The polyamine metabolism in cyanobacteria in response to environmental stress is of interest in research on stress tolerance of higher plants. In this thesis, the potd gene encoding an polyamine transporter subunit from Synechocystis sp. strain PCC 6803 was characterized for the first time. A homology model built for PotD protein indicated that it has capability of binding polyamines, with the preference for spermidine. Furthermore, in order to investigate the structural features of the substrate specificity, polyamines were docked into the binding site. Spermidine was positioned very similarly in Synechocystis PotD as in the template structure and had most favorable interactions of the docked polyamines. Based on the homology model, experimental work was conducted, which confirmed the binding preference. Flavodiiron proteins (Flv) are enzymes, which protect the cell against toxicity of oxygen and/or nitric oxide by reduction. In this thesis, we present a novel type of photoprotection mechanism in cyanobacteria by the heterodimer of Flv2/Flv4. The constructed homology model of Flv2/Flv4 suggests a functional heterodimer capable of rapid electron transfer. The unknown protein sll0218, encoded by the flv2-flv4 operon, is assumed to facilitate the interaction of the Flv2/Flv4 heterodimer and energy transfer between the phycobilisome and PSII. Flv2/Flv4 provides an alternative electron transfer pathway and functions as an electron sink in PSII electron transfer.

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).