Accurate prediction of interfacial residues in two-domain proteins using evolutionary information: Implications for three-dimensional modeling

With the preponderance of multidomain proteins in eukaryotic genomes, it is essential to recognize the constituent domains and their functions. Often function involves communications across the domain interfaces, and the knowledge of the interacting sites is essential to our understanding of the structure-function relationship. Using evolutionary information extracted from homologous domains in at least two diverse domain architectures (single and multidomain), we predict the interface residues corresponding to domains from the two-domain proteins. We also use information from the three-dimensional structures of individual domains of two-domain proteins to train naive Bayes classifier model to predict the interfacial residues. Our predictions are highly accurate (approximate to 85%) and specific (approximate to 95%) to the domain-domain interfaces. This method is specific to multidomain proteins which contain domains in at least more than one protein architectural context. Using predicted residues to constrain domain-domain interaction, rigid-body docking was able to provide us with accurate full-length protein structures with correct orientation of domains. We believe that these results can be of considerable interest toward rational protein and interaction design, apart from providing us with valuable information on the nature of interactions. Proteins 2014; 82:1219-1234. (c) 2013 Wiley Periodicals, Inc.

The HORMA domain: an evolutionarily conserved domain discovered in chromatin-associated proteins, has unanticipated diverse functions

The HORMA domain (for Hop1p, Rev7p and MAD2) was discovered in three chromatin-associated proteins in the budding yeast Saccharomyces cerevisiae. This domain has also been found in proteins with similar functions in organisms including plants, animals and nematodes. The HORMA domain containing proteins are thought to function as adaptors for meiotic checkpoint protein signaling and in the regulation of meiotic recombination. Surprisingly, new work has disclosed completely unanticipated and diverse functions for the HORMA domain containing proteins. A. M. Villeneuve and colleagues (Schvarzstein et al., 2013) show that meiosis-specific HORMA domain containing proteins plays a vital role in preventing centriole disengagement during Caenorhabditis elegans spermatocyte meiosis. Another recent study reveals that S. cerevisiae Atg13 HORMA domain acts as a phosphorylation-dependent conformational switch in the cellular autophagic process. (C) 2014 Elsevier B.V. All rights reserved.

Phylogenetic relationships and classification of thiolases and thiolase-like proteins of Mycobacterium tuberculosis and Mycobacterium smegmatis

Thiolases are enzymes involved in lipid metabolism. Thiolases remove the acetyl-CoA moiety from 3-ketoacyl-CoAs in the degradative reaction. They can also catalyze the reverse Claisen condensation reaction, which is the first step of biosynthetic processes such as the biosynthesis of sterols and ketone bodies. In human, six distinct thiolases have been identified. Each of these thiolases is different from the other with respect to sequence, oligomeric state, substrate specificity and subcellular localization. Four sequence fingerprints, identifying catalytic loops of thiolases, have been described. In this study genome searches of two mycobacterial species (Mycobacterium tuberculosis and Mycobacterium smegmatis), were carried out, using the six human thiolase sequences as queries. Eight and thirteen different thiolase sequences were identified in M. tuberculosis and M. smegmatis, respectively. In addition, thiolase-like proteins (one encoded in the Mtb and two in the Msm genome) were found. The purpose of this study is to classify these mostly uncharacterized thiolases and thiolase-like proteins. Several other sequences obtained by searches of genome databases of bacteria, mammals and the parasitic protist family of the Trypanosomatidae were included in the analysis. Thiolase-like proteins were also found in the trypanosomatid genomes, but not in those of mammals. In order to study the phylogenetic relationships at a high confidence level, additional thiolase sequences were included such that a total of 130 thiolases and thiolase-like protein sequences were used for the multiple sequence alignment. The resulting phylogenetic tree identifies 12 classes of sequences, each possessing a characteristic set of sequence fingerprints for the catalytic loops. From this analysis it is now possible to assign the mycobacterial thiolases to corresponding homologues in other kingdoms of life. The results of this bioinformatics analysis also show interesting differences between the distributions of M. tuberculosis and M. smegmatis thiolases over the 12 different classes. (C) 2014 Elsevier Ltd. All rights reserved.

Propensity to Form Amyloid Fibrils Is Encoded as Excitations in the Free Energy Landscape of Monomeric Proteins

Protein aggregation, linked to many of diseases, is initiated when monomers access rogue conformations that are poised to form amyloid fibrils. We show, using simulations of src SH3 domain, that mechanical force enhances the population of the aggregation-prone (N*) states, which are rarely populated under force free native conditions but are encoded in the spectrum of native fluctuations. The folding phase diagrams of SH3 as a function of denaturant concentration (C]), mechanical force (f), and temperature exhibit an apparent two-state behavior, without revealing the presence of the elusive N* states. Interestingly, the phase boundaries separating the folded and unfolded states at all C] and f fall on a master curve, which can be quantitatively described using an analogy to superconductors in a magnetic field. The free energy profiles as a function of the molecular extension (R), which are accessible in pulling experiments, (R), reveal the presence of a native-like N* with a disordered solvent-exposed amino-terminal beta-strand. The structure of the N* state is identical with that found in Fyn SH3 by NMR dispersion experiments. We show that the timescale for fibril formation can be estimated from the population of the N* state, determined by the free energy gap separating the native structure and the N* state, a finding that can be used to assess fibril forming tendencies of proteins. The structures of the N* state are used to show that oligomer formation and likely route to fibrils occur by a domain-swap mechanism in SH3 domain. (C) 2014 Elsevier Ltd. All rights reserved.

Re-analysis of cryoEM data on HCV IRES bound to 40S subunit of human ribosome integrated with recent structural information suggests new contact regions between ribosomal proteins and HCV RNA

In this study, we combine available high resolution structural information on eukaryotic ribosomes with low resolution cryo-EM data on the Hepatitis C Viral RNA (IRES) human ribosome complex. Aided further by the prediction of RNA-protein interactions and restrained docking studies, we gain insights on their interaction at the residue level. We identified the components involved at the major and minor contact regions, and propose that there are energetically favorable local interactions between 40S ribosomal proteins and IRES domains. Domain II of the IRES interacts with ribosomal proteins S5 and S25 while the pseudoknot and the downstream domain IV region bind to ribosomal proteins S26, S28 and S5. We also provide support using UV cross-linking studies to validate our proposition of interaction between the S5 and IRES domains II and IV. We found that domain IIIe makes contact with the ribosomal protein S3a (S1e). Our model also suggests that the ribosomal protein S27 interacts with domain IIIc while S7 has a weak contact with a single base RNA bulge between junction IIIabc and IIId. The interacting residues are highly conserved among mammalian homologs while IRES RNA bases involved in contact do not show strict conservation. IRES RNA binding sites for S25 and S3a show the best conservation among related viral IRESs. The new contacts identified between ribosomal proteins and RNA are consistent with previous independent studies on RNA-binding properties of ribosomal proteins reported in literature, though information at the residue level is not available in previous studies.

CLAP: A web-server for automatic classification of proteins with special reference to multi-domain proteins

Background: The function of a protein can be deciphered with higher accuracy from its structure than from its amino acid sequence. Due to the huge gap in the available protein sequence and structural space, tools that can generate functionally homogeneous clusters using only the sequence information, hold great importance. For this, traditional alignment-based tools work well in most cases and clustering is performed on the basis of sequence similarity. But, in the case of multi-domain proteins, the alignment quality might be poor due to varied lengths of the proteins, domain shuffling or circular permutations. Multi-domain proteins are ubiquitous in nature, hence alignment-free tools, which overcome the shortcomings of alignment-based protein comparison methods, are required. Further, existing tools classify proteins using only domain-level information and hence miss out on the information encoded in the tethered regions or accessory domains. Our method, on the other hand, takes into account the full-length sequence of a protein, consolidating the complete sequence information to understand a given protein better. Results: Our web-server, CLAP (Classification of Proteins), is one such alignment-free software for automatic classification of protein sequences. It utilizes a pattern-matching algorithm that assigns local matching scores (LMS) to residues that are a part of the matched patterns between two sequences being compared. CLAP works on full-length sequences and does not require prior domain definitions. Pilot studies undertaken previously on protein kinases and immunoglobulins have shown that CLAP yields clusters, which have high functional and domain architectural similarity. Moreover, parsing at a statistically determined cut-off resulted in clusters that corroborated with the sub-family level classification of that particular domain family. Conclusions: CLAP is a useful protein-clustering tool, independent of domain assignment, domain order, sequence length and domain diversity. Our method can be used for any set of protein sequences, yielding functionally relevant clusters with high domain architectural homogeneity. The CLAP web server is freely available for academic use at http://nslab.mbu.iisc.ernet.in/clap/.

Sequence and conformational preferences at termini of alpha-helices in membrane proteins: Role of the helix environment

-helices are amongst the most common secondary structural elements seen in membrane proteins and are packed in the form of helix bundles. These -helices encounter varying external environments (hydrophobic, hydrophilic) that may influence the sequence preferences at their N and C-termini. The role of the external environment in stabilization of the helix termini in membrane proteins is still unknown. Here we analyze -helices in a high-resolution dataset of integral -helical membrane proteins and establish that their sequence and conformational preferences differ from those in globular proteins. We specifically examine these preferences at the N and C-termini in helices initiating/terminating inside the membrane core as well as in linkers connecting these transmembrane helices. We find that the sequence preferences and structural motifs at capping (Ncap and Ccap) and near-helical (N' and C') positions are influenced by a combination of features including the membrane environment and the innate helix initiation and termination property of residues forming structural motifs. We also find that a large number of helix termini which do not form any particular capping motif are stabilized by formation of hydrogen bonds and hydrophobic interactions contributed from the neighboring helices in the membrane protein. We further validate the sequence preferences obtained from our analysis with data from an ultradeep sequencing study that identifies evolutionarily conserved amino acids in the rat neurotensin receptor. The results from our analysis provide insights for the secondary structure prediction, modeling and design of membrane proteins. Proteins 2014; 82:3420-3436. (c) 2014 Wiley Periodicals, Inc.

Identifying functionally important cis-peptide containing segments in proteins and their utility in molecular function annotation

Cis-peptide embedded segments are rare in proteins but often highlight their important role in molecular function when they do occur. The high evolutionary conservation of these segments illustrates this observation almost universally, although no attempt has been made to systematically use this information for the purpose of function annotation. In the present study, we demonstrate how geometric clustering and level-specific Gene Ontology molecular-function terms (also known as annotations) can be used in a statistically significant manner to identify cis-embedded segments in a protein linked to its molecular function. The present study identifies novel cis-peptide fragments, which are subsequently used for fragment-based function annotation. Annotation recall benchmarks interpreted using the receiver-operator characteristic plot returned an area-under-curve >0.9, corroborating the utility of the annotation method. In addition, we identified cis-peptide fragments occurring in conjunction with functionally important trans-peptide fragments, providing additional insights into molecular function. We further illustrate the applicability of our method in function annotation where homology-based annotation transfer is not possible. The findings of the present study add to the repertoire of function annotation approaches and also facilitate engineering, design and allied studies around the cis-peptide neighborhood of proteins.

Streptococcus pneumoniae Genome Database (SPGDB): A database for strain specific comparative analysis of Streptococcus pneumoniae genes and proteins

Streptococcus pneumoniae causes pneumonia, septicemia and meningitis. S. pneumoniae is responsible for significant mortality both in children and in the elderly. In recent years, the whole genome sequencing of various S. pneumoniae strains have increased manifold and there is an urgent need to provide organism specific annotations to the scientific community. This prompted us to develop the Streptococcus pneumoniae Genome Database (SPGDB) to integrate and analyze the completely sequenced and available S. pneumoniae genome sequences. Further, links to several tools are provided to compare the pool of gene and protein sequences, and proteins structure across different strains of S. pneumoniae. SPGDB aids in the analysis of phenotypic variations as well as to perform extensive genomics and evolutionary studies with reference to S. pneumoniae. (C) 2014 Elsevier Inc. All rights reserved.

Paralogous cAMP Receptor Proteins in Mycobacterium smegmatis Show Biochemical and Functional Divergence

The cyclic AMP receptor protein (CRP) family of transcription factors consists of global regulators of bacterial gene expression. Here, we identify two paralogous CRPs in the genome of Mycobacterium smegmatis that have 78% identical sequences and characterize them biochemically and functionally. The two proteins (MSMEG_0539 and MSMEG_6189) show differences in cAMP binding affinity, trypsin sensitivity, and binding to a CRP site that we have identified upstream of the msmeg_3781 gene. MSMEG_6189 binds to the CRP site readily in the absence of cAMP, while MSMEG_0539 binds in the presence of cAMP, albeit weakly. msmeg_6189 appears to be an essential gene, while the ?msmeg_0539 strain was readily obtained. Using promoter-reporter constructs, we show that msmeg_3781 is regulated by CRP binding, and its transcription is repressed by MSMEG_6189. Our results are the first to characterize two paralogous and functional CRPs in a single bacterial genome. This gene duplication event has subsequently led to the evolution of two proteins whose biochemical differences translate to differential gene regulation, thus catering to the specific needs of the organism.

Efficacious Gene Silencing in Serum and Significant Apoptotic Activity Induction by Survivin Downregulation Mediated by New Cationic Gemini Tocopheryl Lipids

Nonviral gene delivery offers cationic liposomes as promising instruments for the delivery of double-stranded RNA (ds RNA) molecules for successful sequence-specific gene silencing (RNA interference). The efficient delivery of siRNA (small interfering RNA) to cells while avoiding unexpected side effects is an important prerequisite for the exploitation of the power of this excellent tool. We present here six new tocopherol based cationic gemini lipids, which induce substantial gene knockdown without any obvious cytotoxicity. All the efficient coliposomal formulations derived from each of these geminis and a helper lipid, dioleoylphosphatidylethanolamine (DOPE), were well characterized using physical methods such as atomic force microscopy (AFM) and dynamic light scattering (DLS). Zeta potential measurements were conducted to estimate the surface charge of these formulations. Flow cytometric analysis showed that the optimized coliposomal formulations could transfect anti-GFP siRNA efficiently in three different GFP expressing cell lines, viz., HEK 293T, HeLa, and Caco-2, significantly better than a potent commercial standard Lipofectamine 2000 (L2K) both in the absence and in the presence of serum (FBS). Notably, the knockdown activity of coliposomes of gemini lipids was not affected even in the presence of serum (10% and 50% FBS) while it dropped down for L2K significantly. Observations under a fluorescence microscope, RT-PCR, and Western blot analysis substantiated the flow cytometry results. The efficient cellular entry of labeled siRNA in GFP expressing cells as evidenced from confocal microscopy put forward these gemini lipids among the potent lipidic carriers for siRNA. The efficient transfection capabilities were also profiled in a more relevant fashion while performing siRNA transfections against survivin (an anti-apoptotic protein) which induced substantial apoptosis. Furthermore, the survivin downregulation improved the therapeutic efficacy levels of an anticancer drug, doxorubicin, significantly. In short, the new tocopherol based gemini lipids appear to be highly promising for achieving siRNA mediated gene knockdown in various cell lines.

Remarkable visible light-triggered cytotoxicity of mitochondria targeting mixed-ligand cobalt(III) complexes of curcumin and phenanthroline bases binding to human serum albumin

Six new mixed-ligand cobalt(III) complexes of formulation Co(N-N)(2)(O-O)](ClO4)(2) (1-6), where N-N is a N,N-donor phenanthroline base, namely, 1,10-phenanthroline (phen in 1, 2), dipyrido3,2-d:2',3'-f] quinoxaline (dpq in 3, 4), and dipyrido3,2-a:2',3'-c]phenazine (dppz in 5, 6), O-O is acetylacetonate (acac in 1, 3, 5) or curcumin (bis(4-hydroxy-3-methoxyphenyl)-1,6-diene-3,5-dione, cur in 2, 4, 6), have been synthesized and characterized. The X-ray crystal structures of complex 1 (as PF6- salt, 1a) and 3 show distorted octahedral geometries formed by the CoN4O2 core. The complexes 1, 3 and 5 having the simple acac ligand are prepared as control species to understand the role of curcumin. The optimized geometries and the frontier orbitals of the curcumin complexes 2, 4, and 6 are obtained from the DFT calculations. The complexes 2, 4, and 6 having the photoactive curcumin moiety display an absorption band in the visible region near 420 nm and show remarkable photocytotoxicity in HeLa cancer cells with respective IC50 values of 7.4 mu M, 5.1 mu M and 1.6 mu M while being much less toxic in dark. MTT assay using complex 6 shows that it is not significantly photocytotoxic to MCF-10A normal cells. The control complexes having the acac ligand are non-toxic both in the presence and absence of light. The cell death is apoptotic in nature and triggered by the photogeneration of reactive oxygen species. Fluorescence imaging experiments on HeLa cells reveals that complex 6 accumulated primarily inside the mitochondria. Human serum albumin (HSA) binding experiments show that the complexes bind HSA with good affinity, but 6 binds with the highest affinity, with a K-b value of 9.8 x 10(5) M-1. Thus, complex 6 with its negligible toxicity in the dark and in normal cells but remarkable toxicity in visible light holds significant photochemotherapeutic potential.

Co-liposomes of redox-active alkyl-ferrocene modified low MW branched PEI and DOPE for efficacious gene delivery in serum

Herein, we present six new lipopolymers based on low molecular weight, branched polyethylenimine (BPEI 800 Da) which are hydrophobically modified using ferrocene terminated alkyl tails of variable lengths. The effects of degree of grafting, spacer length and the redox state of ferrocene in the lipopolymers on the self assembly properties were investigated in detail by TEM, AFM, DLS and zeta potential measurements. The assemblies displayed an oxidation induced increase in the size of the aggregates. The co-liposomes comprising the lipopolymer and a helper lipid, 1,2-dioleoyl phosphatidyl ethanolamine (DOPE), showed excellent gene (pDNA) delivery capability in a serum containing environment in two cancer cell lines (HeLa and U251 cells). Optimized formulations showed remarkably higher transfection activity than BPEI (25 kDa) and were also significantly better than a commercial transfection reagent, Lipofectamine 2000 as evidenced from both the luciferase activity and GFP expression analysis. Oxidation of ferrocene in the lipopolymers led to drastically reduced levels of gene transfection which was substantiated by reduced cellular internalization of fluorescently labelled pDNA as detected using confocal microscopy and flow cytometry. Moreover, the transfection inactive oxidized lipopolyplexes could be turned transfection active by exposure to ascorbic acid (AA) in cell culture medium during transfection. Endocytosis inhibition experiments showed that gene expression mediated by reduced formulations involved both clathrin and caveolae mediated pathways while the oxidized formulations were routed via the caveolae. Cytotoxicity assays revealed no obvious toxicity for the lipopolyplexes in the range of optimized transfection levels in both the cell lines studied. Overall, we have exploited the redox activity of ferrocene in branched PEI-based efficient polymeric gene carriers whose differential transfection activities could be harnessed for spatial or temporal cellular transfections.

Structural and functional analysis of two universal stress proteins YdaA and YnaF from Salmonella typhimurium: possible roles in microbial stress tolerance

In many organisms ``Universal Stress Proteins'' CUSPS) are induced in response to a variety of environmental stresses. Here we report the structures of two USPs, YnaF and YdaA from Salmonella typhimurium determined at 1.8 angstrom and 2.4 angstrom resolutions, respectively. YnaF consists of a single USP domain and forms a tetrameric organization stabilized by interactions mediated through chloride ions. YdaA is a larger protein consisting of two tandem USP domains. Two protomers of YdaA associate to form a structure similar to the YnaF tetramer. YdaA showed ATPase activity and an ATP binding motif G-2X-G-9X-G(S/T/N) was found in its C-terminal domain. The residues corresponding to this motif were not conserved in YnaF although YnaF could bind ATP. However, unlike YdaA, YnaF did not hydrolyse ATP in vitro. Disruption of interactions mediated through chloride ions by selected mutations converted YnaF into an ATPase. Residues that might be important for ATP hydrolysis could be identified by comparing the active sites of native and mutant structures. Only the C-terminal domain of YdaA appears to be involved in ATP hydrolysis. The structurally similar N-terminal domain was found to bind a zinc ion near the segment equivalent to the phosphate binding loop of the C-terminal domain. Mass spectrometric analysis showed that YdaA might bind a ligand of approximate molecular weight 800 daltons. Structural comparisons suggest that the ligand, probably related to an intermediate in lipid A biosynthesis, might bind at a site close to the zinc ion. Therefore, the N-terminal domain of YdaA binds zinc and might play a role in lipid metabolism. Thus, USPs appear to perform several distinct functions such as ATP hydrolysis, altering membrane properties and chloride sensing. (C) 2015 Elsevier Inc. All rights reserved.

Efficacious redox-responsive gene delivery in serum by ferrocenylated monomeric and dimeric cationic cholesterols

Herein, we present the design and synthesis of new redox-active monomeric and dimeric (gemini) cationic lipids based on ferrocenylated cholesterol derivatives for gene delivery. The cationic cholesterols are shown to be transfection efficient after being formulated with the neutral helper lipid DOPE in the presence of serum (FBS). The redox activity of the resulting co-liposomes and their lipoplexes could be regulated using the alkanyl ferrocene moiety attached to the ammonium head groups of the cationic cholesterols. Atomic force microscopy (AFM), dynamic light scattering (DLS) and zeta potential measurements were performed to characterize the co-liposomal aggregates and their complexes with pDNA. The transfection efficiency of lipoplexes could be tuned by changing the oxidation state of the ferrocene moiety. The gene transfection capability was assayed in terms of green fluorescence protein (GFP) expression using pEGFP-C3 plasmid DNA in three cell lines of different origins, namely Caco-2, HEK293T and HeLa, in the presence of serum. The vesicles possessing ferrocene in the reduced state induced an efficient transfection, even better than a commercial reagent Lipofectamine 2000 (Lipo 2000) as evidenced by flow cytometry and fluorescence microscopy. All the co-liposomes containing the oxidized ferrocene displayed diminished levels of gene expression. Gene transfection events from the oxidized co-liposomes were further potentiated by introducing ascorbic acid (AA) as a reducing agent during lipoplex incubation with cells, leading to the resumption of transfection activity. Assessment of transfection capability of both reduced and oxidized co-liposomes was also undertaken following cellular internalization of labelled pDNA using confocal microscopy and flow cytometry. Overall, we demonstrate here controlled gene transfection activities using redox-driven, transfection efficient cationic monomeric and dimeric cholesterol lipids. Such systems could be used in gene delivery applications where transfection needs to be performed spatially or temporally.