De novo inference of protein function from coarse‐grained dynamics

Inference of molecular function of proteins is the fundamental task in the quest for understanding cellular processes. The task is getting increasingly difficult with thousands of new proteins discovered each day. The difficulty arises primarily due to lack of high‐throughput experimental technique for assessing protein molecular function, a lacunae that computational approaches are trying hard to fill. The latter too faces a major bottleneck in absence of clear evidence based on evolutionary information. Here we propose a de novo approach to annotate protein molecular function through structural dynamics match for a pair of segments from two dissimilar proteins, which may share even <10% sequence identity. To screen these matches, corresponding 1 µs coarse‐grained (CG) molecular dynamics trajectories were used to compute normalized root‐mean‐square‐fluctuation graphs and select mobile segments, which were, thereafter, matched for all pairs using unweighted three‐dimensional autocorrelation vectors. Our in‐house custom‐built forcefield (FF), extensively validated against dynamics information obtained from experimental nuclear magnetic resonance data, was specifically used to generate the CG dynamics trajectories. The test for correspondence of dynamics‐signature of protein segments and function revealed 87% true positive rate and 93.5% true negative rate, on a dataset of 60 experimentally validated proteins, including moonlighting proteins and those with novel functional motifs. A random test against 315 unique fold/function proteins for a negative test gave >99% true recall. A blind prediction on a novel protein appears consistent with additional evidences retrieved therein. This is the first proof‐of‐principle of generalized use of structural dynamics for inferring protein molecular function leveraging our custom‐made CG FF, useful to all. Proteins 2014; 82:2443–2454. © 2014 Wiley Periodicals, Inc.     

Mutational Analysis of the Shigella flexneri O-Antigen Polymerase Wzy: Identification of Wzz-Dependent Wzy Mutants

The O-antigen (Oag) component of lipopolysaccharide (LPS) is a major virulence determinant of Shigella flexneri and is synthesized by the O-antigen polymerase, Wzy_Sf. Oag chain length is regulated by chromosomally encoded Wzz_Sf and pHS-2 plasmid-encoded Wzz_pHS2. To identify functionally important amino acid residues in Wzy_Sf, random mutagenesis was performed on the wzy_Sf gene in a pWaldo-TEV-GFP plasmid, followed by screening with colicin E2. Analysis of the LPS conferred by mutated Wzy_Sf proteins in the wzy_Sf-deficient (Δwzy) strain identified 4 different mutant classes, with mutations found in periplasmic loop 1 (PL1), PL2, PL3, and PL6, transmembrane region 2 (TM2), TM4, TM5, TM7, TM8, and TM9, and cytoplasmic loop 1 (CL1) and CL5. The association of Wzy_Sf and Wzz_Sf was investigated by transforming these mutated wzy_Sf plasmids into a wzy_Sf- and wzz_Sf-deficient (Δwzy Δwzz) strain. Comparison of the LPS profiles in the Δwzy and Δwzy Δwzz backgrounds identified Wzy_Sf mutants whose polymerization activities were Wzz_Sf dependent. Colicin E2 and bacteriophage Sf6c sensitivities were consistent with the LPS profiles. Analysis of the expression levels of the Wzy_Sf-GFP mutants in the Δwzy and Δwzy Δwzz backgrounds identified a role for Wzz_Sf in Wzy_Sf stability. Hence, in addition to its role in regulating Oag modal chain length, Wzz_Sf also affects Wzy_Sf activity and stability.     

Stringent response in Vibrio cholerae: genetic analysis of spoT gene function and identification of a novel (p)ppGpp synthetase gene

RelA and SpoT of Gram-negative organisms critically regulate cellular levels of (p)ppGpp. Here, we have dissected the spoT gene function of the cholera pathogen Vibrio cholerae by extensive genetic analysis. Unlike Escherichia coli , V. cholerae Δ relA Δ spoT cells accumulated (p)ppGpp upon fatty acid or glucose starvation. The result strongly suggests RelA-SpoT-independent (p)ppGpp synthesis in V. cholerae . By repeated subculturing of a V. cholerae Δ relA Δ spoT mutant, a suppressor strain with (p)ppGpp⁰ phenotype was isolated. Bioinformatics analysis of V. cholerae whole genome sequence allowed identification of a hypothetical gene ( VC1224 ), which codes for a small protein (∼29 kDa) with a (p)ppGpp synthetase domain and the gene is highly conserved in vibrios; hence it has been named relV . Using E. coli Δ relA or Δ relA Δ spoT mutant we showed that relV indeed codes for a novel (p)ppGpp synthetase. Further analysis indicated that relV gene of the suppressor strain carries a point mutation at nucleotide position 676 of its coding region (Δ relA Δ spoT relV676 ), which seems to be responsible for the (p)ppGpp⁰ phenotype. Analysis of a V. cholerae Δ relA Δ spoT Δ relV triple mutant confirmed that apart from canonical relA and spoT genes, relV is a novel gene in V. cholerae responsible for (p)ppGpp synthesis.     

Mutational analysis of the (p)ppGpp synthetase activity of the Rel enzyme of Mycobacterium tuberculosis

Rel_Mtb, a GTP pyrophosphokinase encoded by the Mycobacterium tuberculosis (Mtb) genome, catalyzes synthesis of (p)ppGpp from ATP and GDP(GTP) and its hydrolysis to GDP(GTP) and pyrophosphate to mediate stringent response, which helps bacteria to survive during nutrient limitation. Like other members of Rel_Spo homologs, Rel_Mtb has four distinct domains: HD, Rel_Spo (RSD), TGS and ACT. The N-terminal HD and RSD are responsible for (p)ppGpp hydrolysis and synthesis, respectively. In this study, we have dissected the rel _Mtb gene function and determined the minimal region essential for (p)ppGpp synthetic activity. The Rel_Mtb and its truncated derivatives were expressed from an arabinose inducible promoter (P _BAD ), and in vivo functional analyses were done in a (p)ppGpp null Escherichia coli strain. Our results indicate that only 243 amino acids (188–430 residues) containing fragment are sufficient for Rel_Mtb (p)ppGpp synthetic activity. The results were further confirmed by in vitro assays using purified proteins. We further characterized the RSD of Rel_Mtb by substituting several conserved amino acids with structurally related residues and identified six such residues, which appeared to be critical for maintaining its catalytic activity. Furthermore, we have also extended our analysis to an RSD encoding gene rv1366 of Mtb, and experimental results indicated that the encoded protein Rv1366 is unable to synthesize (p)ppGpp.     

How does Sec63 affect the conformation of Sec61 in yeast?

The Sec complex catalyzes the translocation of proteins of the secretory pathway into the endoplasmic reticulum and the integration of membrane proteins into the endoplasmic reticulum membrane. Some substrate peptides require the presence and involvement of accessory proteins such as Sec63. Recently, a structure of the Sec complex from Saccharomyces cerevisiae, consisting of the Sec61 channel and the Sec62, Sec63, Sec71 and Sec72 proteins was determined by cryo-electron microscopy (cryo-EM). Here, we show by co-precipitation that the Sec61 channel subunit Sbh1 is not required for formation of stable Sec63-Sec61 contacts. Molecular dynamics simulations started from the cryo-EM conformation of Sec61 bound to Sec63 and of unbound Sec61 revealed how Sec63 affects the conformation of Sec61 lateral gate, plug, pore region and pore ring diameter via three intermolecular contact regions. Molecular docking of SRP-dependent vs. SRP-independent signal peptide chains into the Sec61 channel showed that the pore regions affected by presence/absence of Sec63 play a crucial role in positioning the signal anchors of SRP-dependent substrates nearby the lateral gate.     

Binding affinity and in vitro cytotoxicity of harmaline targeting different motifs of nucleic acids: An ultimate drug designing approach

The work focuses towards interaction of harmaline, with nucleic acids of different motifs by multispectroscopic and calorimetric techniques. Findings of this study suggest that binding constant varied in the order single‐stranded (ss) poly(A) > double‐stranded calf thymus (CT) DNA > double‐stranded poly(G)·poly(C) > clover leaf tRNA^Phe. Prominent structural changes of ss poly(A), CT DNA, and poly(G)· poly(C) with concomitant induction of optical activity in the bound achiral alkaloid molecule was observed, while with tRNA^Phe, very weak induced circular dichroism perturbation was seen. The interaction was predominantly exothermic, enthalpy driven, and entropy favored with CT DNA and poly(G)·poly(C), while it was entropy driven with poly(A) and tRNA^Phe. Intercalated state of harmaline inside poly(A), CT DNA, and poly(G)·poly(C) was shown by viscometry, ferrocyanide quenching, and molecular docking. All these findings unequivocally pointed out preference of harmaline towards ss poly(A) inducing self‐structure formation. Furthermore, harmaline administration caused a significant decrease in proliferation of HeLa and HepG₂ cells with GI₅₀ of 28μM and 11.2μM, respectively. Nucleic acid fragmentation, cellular ultramorphological changes, decreased mitochondrial membrane potential, upregulation of p53 and caspase 3, generation of reactive oxygen species, and a significant increase in the G₂/M population made HepG₂ more prone to apoptosis than are HeLa cells.