Mapping and molecular modeling of S-adenosyl-L-methionine binding sites in N-methyltransferase domains of the multifunctional polypeptide cyclosporin synthetase

We employed a highly specific photoaffinity labeling procedure, using 14C-labeled S-adenosyl-L-methionine (AdoMet) to define the chemical structure of the AdoMet binding centers on cyclosporin synthetase (CySyn). Tryptic digestion of CySyn photolabeled with either [methyl-14C]AdoMet or [carboxyl-14C]AdoMet yielded the sequence H2N-Asn-Asp-Gly-Leu-Glu-Ser-Tyr-Val-Gly-Ile-Glu-Pro-Ser-Arg-COOH (residues 10644-10657), situated within the N-methyltransferase domain of module 8 of CySyn. Radiosequencing detected Glu10654 and Pro10655 as the major sites of derivatization. [carboxyl-14C]AdoMet in addition labeled Tyr10650. Chymotryptic digestion generated the radiolabeled peptide H2N-Ile-Gly-Leu-Glu-Pro-Ser-Gln-Ser-Ala-Val-Gln-Phe-COOH, corresponding to amino acids 2125-2136 of the N-methyltransferase domain of module 2. The radiolabeled amino acids were identified as Glu2128 and Pro2129, which are equivalent in position and function to the modified residues identified with tryptic digestions in module 8. Homology modeling of the N-methyltransferase domains indicates that these regions conserve the consensus topology of the AdoMet binding fold and consensus cofactor interactions seen in structurally characterized AdoMet-dependent methyltransferases. The modified sequence regions correspond to the motif II consensus sequence element, which is involved in directly complexing the adenine and ribose components of AdoMet. We conclude that the AdoMet binding to nonribosomal peptide synthetase N-methyltransferase domains obeys the consensus cofactor interactions seen among most structurally characterized low molecular weight AdoMet-dependent methyltransferases. <br />

Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the ��2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the ��2 adrenergic crystal structure by optimizing the binding site using ��2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

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.

Adaptive evolution of color vision as seen through the eyes of butterflies

Butterflies and primates are interesting for comparative color vision studies, because both have evolved middle- (M) and long-wavelength- (L) sensitive photopigments with overlapping absorbance spectrum maxima (lambda(max) values). Although positive selection is important for the maintenance of spectral variation within the primate pigments, it remains an open question whether it contributes similarly to the diversification of butterfly pigments. To examine this issue, we performed epimicrospectrophotometry on the eyes of five Limenitis butterfly species and found a 31-nm range of variation in the lambda(max) values of the L-sensitive photopigments (514-545 nm). We cloned partial Limenitis L opsin gene sequences and found a significant excess of replacement substitutions relative to polymorphisms among species. Mapping of these L photopigment lambda(max) values onto a phylogeny revealed two instances within Lepidoptera of convergently evolved L photopigment lineages whose lambda(max) values were blue-shifted. A codon-based maximum-likelihood analysis indicated that, associated with the two blue spectral shifts, four amino acid sites (Ile17Met, Ala64Ser, Asn70Ser, and Ser137Ala) have evolved substitutions in parallel and exhibit significant d(N)/d(S) >1. Homology modeling of the full-length Limenitis arthemis astyanax L opsin placed all four substitutions within the chromophore-binding pocket. Strikingly, the Ser137Ala substitution is in the same position as a site that in primates is responsible for a 5- to 7-nm blue spectral shift. Our data show that some of the same amino acid sites are under positive selection in the photopigments of both butterflies and primates, spanning an evolutionary distance >500 million years.

Proteolytic activation of Chlamydia trachomatis HTRA is mediated by PDZ1 domain interactions with protease domain loops L3 and LC and beta strand ��5

Chlamydia trachomatis is a bacterial pathogen responsible for one of the most prevalent sexually transmitted infections worldwide. Its unique development cycle has limited our understanding of its pathogenic mechanisms. However, CtHtrA has recently been identified as a potential C. trachomatis virulence factor. CtHtrA is a tightly regulated quality control protein with a monomeric structural unit comprised of a chymotrypsin-like protease domain and two PDZ domains. Activation of proteolytic activity relies on the C-terminus of the substrate allosterically binding to the PDZ1 domain, which triggers subsequent conformational change and oligomerization of the protein into 24-mers enabling proteolysis. This activation is mediated by a cascade of precise structural arrangements, but the specific CtHtrA residues and structural elements required to facilitate activation are unknown. Using in vitro analysis guided by homology modeling, we show that the mutation of residues Arg362 and Arg224, predicted to disrupt the interaction between the CtHtrA PDZ1 domain and loop L3, and between loop L3 and loop LD, respectively, are critical for the activation of proteolytic activity. We also demonstrate that mutation to residues Arg299 and Lys160, predicted to disrupt PDZ1 domain interactions with protease loop LC and strand ��5, are also able to influence proteolysis, implying their involvement in the CtHtrA mechanism of activation. This is the first investigation of protease loop LC and strand ��5 with respect to their potential interactions with the PDZ1 domain. Given their high level of conservation in bacterial HtrA, these structural elements may be equally significant in the activation mechanism of DegP and other HtrA family members.

Structural studies on members of the picornavirus superfamily

The pathogenic members of the picornavirus superfamily have adverse effects on humans, their crops and their livestock. As structure is related to function, detailed structural studies on these viruses are important not only for fundamental understanding of the viral life cycle, but also for the rational design of vaccines and inhibitors for disease control. These viruses have positive sense, single-stranded RNA genomes enclosed in a protein capsid. X-ray crystallography and cryo-electron microscopy studies have revealed that the isometric members of this group have icosahedrally-symmetric capsids made up of 60 copies of each of the structural proteins. The members that infect animal cells often employ one or more cellular receptors to facilitate cell entry which in some cases is known to initiate the uncoating sequence of the genome. The nature of the interactions between individual viruses and alternative cellular receptors has rarely been probed. The capsid assembly of the members of the picornavirus superfamily is considered to be cooperative and the interactions of RNA and capsid proteins are thought to play an important role in orchestrating virus assembly. The major aims of this thesis were to solve the structures of blackcurrant reversion virus (BRV), human parechovirus 1 (HPEV1) and coxsackievirus A7 (CAV7), as well as the structure of HPEV1 complexed with two of its cellular receptors using cryo-electron microscopy, three-dimensional image reconstruction and homology modeling. Each of the selected viruses represents a taxonomic group where little or no structural data was previously available. The results enabled the detailed comparison of the new structures to those of known picornaviruses, the identification of surface-exposed epitopes potentially important for host interaction, the mapping of RNA-capsid protein interactions and the elucidation of the basis for the specificity of two different receptor molecules for the same capsid. This work will form the basis for further studies on the influence of RNA on parechovirus assembly as a potential target for drug design.

Identification of Critical Residues of the Mycobacterial Dephosphocoenzyme A Kinase by Site-Directed Mutagenesis

Dephosphocoenzyme A kinase performs the transfer of the c-phosphate of ATP to dephosphocoenzyme A, catalyzing the last step of coenzyme A biosynthesis. This enzyme belongs to the P-loop-containing NTP hydrolase superfamily, all members of which posses a three domain topology consisting of a CoA domain that binds the acceptor substrate, the nucleotide binding domain and the lid domain. Differences in the enzymatic organization and regulation between the human and mycobacterial counterparts, have pointed out the tubercular CoaE as a high confidence drug target (HAMAP database). Unfortunately the absence of a three-dimensional crystal structure of the enzyme, either alone or complexed with either of its substrates/regulators, leaves both the reaction mechanism unidentified and the chief players involved in substrate binding, stabilization and catalysis unknown. Based on homology modeling and sequence analysis, we chose residues in the three functional domains of the enzyme to assess their contributions to ligand binding and catalysis using site-directed mutagenesis. Systematically mutating the residues from the P-loop and the nucleotide-binding site identified Lys14 and Arg140 in ATP binding and the stabilization of the phosphoryl intermediate during the phosphotransfer reaction. Mutagenesis of Asp32 and Arg140 showed catalytic efficiencies less than 5-10% of the wild type, indicating the pivotal roles played by these residues in catalysis. Non-conservative substitution of the Leu114 residue identifies this leucine as the critical residue from the hydrophobic cleft involved in leading substrate, DCoA binding. We show that the mycobacterial enzyme requires the Mg2+ for its catalytic activity. The binding energetics of the interactions of the mutant enzymes with the substrates were characterized in terms of their enthalpic and entropic contributions by ITC, providing a complete picture of the effects of the mutations on activity. The properties of mutants defective in substrate recognition were consistent with the ordered sequential mechanism of substrate addition for CoaE.

AlignHUSH: Alignment of HMMs using structure and hydrophobicity information

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

A Bifunctional Enzyme That Has Both Monoacylglycerol Acyltransferase and Acyl Hydrolase Activities

Monoacylglycerol acyltransferase (MGAT) catalyzes the synthesis of diacylglycerol, the precursor of triacylglycerol biosynthesis and an important signaling molecule. Here, we describe the isolation and characterization of the peanut (Arachis hypogaea) MGAT gene. The soluble enzyme utilizes invariant histidine-62 and aspartate-67 residues of the acyltransferase motif for its MGAT activity. A sequence analysis revealed the presence of a hydrolase (GXSXG) motif, and enzyme assays revealed the presence of monoacylglycerol (MAG) and lysophosphatidylcholine (LPC) hydrolytic activities, indicating the bifunctional nature of the enzyme. The overexpression of the MGAT gene in yeast (Saccharomyces cerevisiae) caused an increase in triacylglycerol accumulation. Similar to the peanut MGAT, the Arabidopsis (Arabidopsis thaliana) homolog (At1g52760) also exhibited both acyltransferase and hydrolase activities. Interestingly, the yeast homolog lacks the conserved HX4D motif, and it is deficient in the acyltransferase function but exhibits MAG and LPC hydrolase activities. This study demonstrates the presence of a soluble MGAT/hydrolase in plants. The predicted three-dimensional homology modeling and substrate docking suggested the presence of two separate substrate (MAG and LPC)-binding sites in a single polypeptide. Our study describes a soluble bifunctional enzyme that has both MGAT and hydrolase functions.

Cloning, expression, purification, and biochemical characterisation of the FIC motif containing protein of Mycobacterium tuberculosis

The role of FIC (Filamentation induced by cAMP)(2) domain containing proteins in the regulation of many vital pathways, mostly through the transfer of NMPs from NTPs to specific target proteins (NMPylation), in microorganisms, higher eukaryotes, and plants is emerging. The identity and function of FIC domain containing protein of the human pathogen, Mycobacterium tuberculosis, remains unknown. In this regard, M. tuberculosis fic gene (Mtfic) was cloned, overexpressed, and purified to homogeneity for its biochemical characterisation. It has the characteristic FIC motif, HPFREGNGRSTR (HPFxxGNGRxxR), spanning 144th to 155th residue. Neither the His-tagged nor the GST-tagged MtFic protein, overexpressed in Escherichia coil, nor expression of Mtfic in Mycobacterium smegmatis, yielded the protein in the soluble fraction. However, the maltose binding protein (MBP) tagged MtFic (MBP-MtFic) could be obtained partly in the soluble fraction. The cloned, overexpressed, and purified recombinant MBP-MtFic showed conversion of ATP, GTP, CTP, and UTP into AMP. GMP, CMP, and UMP, respectively. Sequence alignment with several FIC motif containing proteins, complemented with homology modeling on the FIC motif containing protein, VbhT of Bartonella schoenbuchensis as the template, showed conservation and interaction of residues constituting the FIC domain. Site-specific mutagenesis of the His144, or Glu148, or Asn150 of the FIC motif, or of Arg87 residue that constitutes the FIC domain, or complete deletion of the FIC motif, abolished the NTP to NMP conversion activity. The design of NMP formation assay using the recombinant, soluble MtFic would enable identification of its target substrate for NMPylation. (C) 2012 Elsevier Inc. All rights reserved.

Heat and SDS insensitive NDK dimers are largely stabilised by hydrophobic interaction to form functional hexamer in Mycobacterium smegmatis

The primary structure and function of nucleoside diphosphate kinase (NDK), a substrate non-specific enzyme involved in the maintenance of nucleotide pools is also implicated to play pivotal roles in many other cellular processes. NDK is conserved from bacteria to human and forms a homotetramer or hexamer to exhibit its biological activity. However, the nature of the functional oligomeric form of the enzyme differs among different organisms. The functional form of NDKs from many bacterial systems, including that of the human pathogen, Mycobacterium tuberculosis (MtuNDK), is a hexamer, although some bacterial NDKs are tetrameric in nature. The present study addresses the oligomeric property of MsmNDK and how a dimer, the basic subunit of a functional hexamer, is stabilized by hydrogen bonds and hydrophobic interactions. Homology modeling was generated using the three-dimensional structure of MtuNDK as a template; the residues interacting at the monomer-monomer interface of MsmNDK were mapped. Using recombinant enzymes of wild type, catalytically inactive mutant, and monomer-monomer interactive mutants of MsmNDK, the stability of the dimer was verified under heat, SDS, low pH, and methanol. The predicted residues (Gln17, Ser24 and Glu27) were engaged in dimer formation, however the mutated proteins retained the ATPase and GTPase activity even after introducing single (MsmNDK- Q17A, MsmNDK-E27A, and MsmNDK-E27Q) and double (MsmNDK-E27A/Q17A) mutation. However, the monomer monomer interaction could be abolished using methanol, indicating the stabilization of the monomer-monomer interaction by hydrophobic interaction.

A Mycobacterial Phosphoribosyltransferase Promotes Bacillary Survival by Inhibiting Oxidative Stress and Autophagy Pathways in Macrophages and Zebrafish

Mycobacterium tuberculosis employs various strategies to modulate host immune responses to facilitate its persistence in macrophages. The M. tuberculosis cell wall contains numerous glycoproteins with unknown roles in pathogenesis. Here, by using Concanavalin A and LC-MS analysis, we identified a novel mannosylated glycoprotein phosphoribosyltransferase, encoded by Rv3242c from M. tuberculosis cell walls. Homology modeling, bioinformatic analyses, and an assay of phosphoribosyltransferase activity in Mycobacterium smegmatis expressing recombinant Rv3242c (MsmRv3242c) confirmed the mass spectrometry data. Using Mycobacterium marinum-zebrafish and the surrogate MsmRv3242c infection models, we proved that phosphoribosyltransferase is involved in mycobacterial virulence. Histological and infection assays showed that the M. marinum mimG mutant, an Rv3242c orthologue in a pathogenic M. marinum strain, was strongly attenuated in adult zebrafish and also survived less in macrophages. In contrast, infection with wild type and the complemented Delta mimG: Rv3242c M. marinum strains showed prominent pathological features, such as severe emaciation, skin lesions, hemorrhaging, and more zebrafish death. Similarly, recombinant Msm Rv3242c bacteria showed increased invasion in non-phagocytic epithelial cells and longer intracellular survival in macrophages as compared with wild type and vector control M. smegmatis strains. Further mechanistic studies revealed that the Rv3242c- and mimG-mediated enhancement of intramacrophagic survival was due to inhibition of autophagy, reactive oxygen species, and reduced activities of superoxide dismutase and catalase enzymes. Infection with MsmRv3242c also activated the MAPK pathway, NF-kappa B, and inflammatory cytokines. In summary, we show that a novel mycobacterial mannosylated phosphoribosyltransferase acts as a virulence and immunomodulatory factor, suggesting that it may constitute a novel target for antimycobacterial drugs.

Structure Prediction of G-Protein Coupled Receptors

<p>G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.</p> <p>In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.</p> <p>To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.</p> <p>Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.</p>

Cloning,Expression and Sequence Analysis of A Luciferase Gene from the Chinese Firefly Pyrocoelia pygidialis

���������������������������������������(������������)������������mRNA������������������������������������������������������.���������������������������cDNA���������1647���������,������548������������������.������������������������������������������������������:������������������������������������Lampyris noctiluca,L.turkestanicus���Nyctophila cf.caucasica���������������������������������97.8%������������������.���������������������������������������������������������,���������������:���������������Lampyris+Nyctophila������������,������������������������������������������������������������.���������������������������������������������������������������70kDa,���������������������������������������������������.���������������������������������������������,������������������������������������������������������������������������������������������5������������,���������������������������������������������������������������������������������������.���������������������������������3������������������������������������,���13������������������������,���������������������������������.���������������������������������������������������������������������������������������������������������������������������������������������.

Visible Volume: a Robust Measure for Protein Structure Characterization

We propose a new characterization of protein structure based on the natural tetrahedral geometry of the �� carbon and a new geometric measure of structural similarity, called visible volume. In our model, the side-chains are replaced by an ideal tetrahedron, the orientation of which is fixed with respect to the backbone and corresponds to the preferred rotamer directions. Visible volume is a measure of the non-occluded empty space surrounding each residue position after the side-chains have been removed. It is a robust, parameter-free, locally-computed quantity that accounts for many of the spatial constraints that are of relevance to the corresponding position in the native structure. When computing visible volume, we ignore the nature of both the residue observed at each site and the ones surrounding it. We focus instead on the space that, together, these residues could occupy. By doing so, we are able to quantify a new kind of invariance beyond the apparent variations in protein families, namely, the conservation of the physical space available at structurally equivalent positions for side-chain packing. Corresponding positions in native structures are likely to be of interest in protein structure prediction, protein design, and homology modeling. Visible volume is related to the degree of exposure of a residue position and to the actual rotamers in native proteins. In this article, we discuss the properties of this new measure, namely, its robustness with respect to both crystallographic uncertainties and naturally occurring variations in atomic coordinates, and the remarkable fact that it is essentially independent of the choice of the parameters used in calculating it. We also show how visible volume can be used to align protein structures, to identify structurally equivalent positions that are conserved in a family of proteins, and to single out positions in a protein that are likely to be of biological interest. These properties qualify visible volume as a powerful tool in a variety of applications, from the detailed analysis of protein structure to homology modeling, protein structural alignment, and the definition of better scoring functions for threading purposes.