Spectral-envelope-group-delay models for transients

Transient signals such as plosives in speech or Castanets in audio do not have a specific modulation or periodic structure in time domain. However, in the spectral domain they exhibit a prominent modulation structure, which is a direct consequence of their narrow time localization. Based on this observation, a spectral-domain AM-FM model for transients is proposed. The spectral AM-FM model is built starting from real spectral zero-crossings. The AM and FM correspond to the spectral envelope (SE) and group delay (GD), respectively. Taking into account the modulation structure and spectral continuity, a local polynomial regression technique is proposed to estimate the GD function from the real spectral zeros. The SE is estimated based on the phase function computed from the estimated GD. Since the GD estimation is parametric, the degree of smoothness can be controlled directly. Simulation results based on synthetic transient signals generated using a beta density function are presented to analyze the noise-robustness of the SEGD model. Three specific applications are considered: (1) SEGD based modeling of Castanet sounds; (2) appropriateness of the model for transient compression; and (3) determining glottal closure instants in speech using a short-time SEGD model of the linear prediction residue.

Multiple molecular alterations in phosphodegrons 1-3 within AAV2 capsid demonstrates higher hepatic gene transfer efficiency

The success of AAV2 mediated hepatic gene transfer in human trials for diseases such as hemophilia has been hampered by a combination of low transduction efficiency and a robust immune response directed against these vectors. We have previously shown that AAV2 is targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal degradation machinery and modification of the serine(S)/threonine(T) kinase and lysine(K) targets on AAV capsid is beneficial. Thus targeted single mutations of S/T>A(S489A, S498A, T251A) and K>R (K532R) improved the efficiency of gene transfer in vivo as compared to wild type (WT)-AAV2 vectors (∼6-14 fold). In the present study, we evaluated if combined alteration of the phosphodegrons (PD), which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, improves further the gene transfer efficiency. Thus, we generated four multiple mutant vectors (PD: 1+3, S489A+K532R, PD: 1+3, S489A+K532R together with T251 residue which did not lie in any of the phosphodegrons but had shown increased transduction efficiency compared to the WT-AAV2 vector (∼6 fold) and was also conserved in 9 out of 10 AAV serotypes (AAV 1 to 10), PD: 1+3, S489A+K532R+S498A and a fourth combination of PD: 3, K532R+T251. We then evaluated them in vitro and in vivo and compared their gene transfer efficiency with either the WT-AAV2 or the best single mutant S489A-AAV2 vector. The novel multiple mutations on the AAV2 capsid did not affect the overall vector packaging efficiency. All the multiple AAV2 mutants showed superior transduction efficiency in HeLa cells in vitro when compared to either the WT (62-72% Vs 21%) or the single mutant S489A (62-72% Vs 50%) AAV2 vectors as demonstrated by FACS analysis (Fig. 1A). On hepatic gene transfer with 5x10^10 vgs per animal in C57BL/6 mice, all the multiple mutants showed increased transgene expression compared to either the WT-AAV2 (∼15-23 fold) or the S489A single mutant vector (∼2-3 fold) (Fig.1B and C). These novel multiple mutant AAV2 vectors also showed higher vector copy number in murine hepatocytes 4 weeks post transduction, as compared to the WT-AAV2 (∼5-6 Vs 1.4 vector copies/diploid genome) and further higher when compared to the single mutant S489A(∼5-6 fold Vs 3.8 fold) (Fig.1D). Further ongoing studies will demonstrate the therapeutic benefit of one or more of the multiple mutants vectors in preclinical models of hemophilia.

Cyclic AMP-dependent protein lysine acylation in mycobacteria regulates fatty acid and propionate metabolism

Acetylation of lysine residues is a posttranslational modification that is used by both eukaryotes and prokaryotes to regulate a variety of biological processes. Here we identify multiple substrates for the cAMP-dependent protein lysine acetyltransferase from Mycobacterium tuberculosis (KATmt). We demonstrate that a catalytically important lysine residue in a number of FadD (fatty acyl CoA synthetase) enzymes is acetylated by KATmt in a cAMP-dependent manner and that acetylation inhibits the activity of FadD enzymes. A sirtuin-like enzyme can deacetylate multiple FadDs, thus completing the regulatory cycle. Using a strain deleted for the KATmt ortholog in Mycobacterium bovis Bacillus Calmette-Guerin (BCG), we show for the first time that acetylation is dependent on intracellular cAMP levels. KATmt can utilize propionyl CoA as a substrate and, therefore, plays a critical role in alleviating propionyl CoA toxicity in mycobacteria by inactivating acyl CoA synthetase (ACS). The precision by which mycobacteria can regulate the metabolism of fatty acids in a cAMP-dependent manner appears to be unparalleled in other biological organisms and is ideally suited to adapt to the complex environment that pathogenic mycobacteria experience in the host.

C12-Helix development in ()n sequences - spectroscopic characterization of Boc-Aib-4(R)Val]-OMe oligomers

The solution conformations of the -hybrid oligopeptides Boc-Aib-4(R)Val]n-OMe (n = 1-8) in organic solvents have been probed by NMR, IR, and CD spectroscopic methods. In the solid state, this peptide series favors C12-helical conformations, which are backbone-expanded analogues of 310 helices in -peptide sequences. NMR studies of the six- (n = 3) and 16-residue (n = 8) peptides reveal that only two NH protons attached the N-terminus residues Aib(1) and 4(R)Val(2) are solvent-exposed. Sequential NiH-Ni+1H NOEs characteristic of local helical conformations are also observed at the residues. IR studies establish that chain extension leads to a large enhancement in the intensities of the hydrogen-bonded NH stretching bands (3343-3280 cm-1), which suggest elongation of intramolecularly hydrogen-bonded structures. The development of C12-helical structures upon lengthening of the sequence is supported by the NMR and IR observations. The CD spectra of the ()n peptides reveal a negative maximum at ca. 206 nm and a positive maximum at ca. 192 nm, spectral feature that are distinct from those of 310 helices in -peptides.

The Escherichia coli Phosphotyrosine Proteome Relates to Core Pathways and Virulence

While phosphotyrosine modification is an established regulatory mechanism in eukaryotes, it is less well characterized in bacteria due to low prevalence. To gain insight into the extent and biological importance of tyrosine phosphorylation in Escherichia coli, we used immunoaffinity-based phosphotyrosine peptide enrichment combined with high resolution mass spectrometry analysis to comprehensively identify tyrosine phosphorylated proteins and accurately map phosphotyrosine sites. We identified a total of 512 unique phosphotyrosine sites on 342 proteins in E. coli K12 and the human pathogen enterohemorrhagic E. coli (EHEC) O157:H7, representing the largest phosphotyrosine proteome reported to date in bacteria. This large number of tyrosine phosphorylation sites allowed us to define five phosphotyrosine site motifs. Tyrosine phosphorylated proteins belong to various functional classes such as metabolism, gene expression and virulence. We demonstrate for the first time that proteins of a type III secretion system (T3SS), required for the attaching and effacing (A/E) lesion phenotype characteristic for intestinal colonization by certain EHEC strains, are tyrosine phosphorylated by bacterial kinases. Yet, A/E lesion and metabolic phenotypes were unaffected by the mutation of the two currently known tyrosine kinases, Etk and Wzc. Substantial residual tyrosine phosphorylation present in an etk wzc double mutant strongly indicated the presence of hitherto unknown tyrosine kinases in E. coli. We assess the functional importance of tyrosine phosphorylation and demonstrate that the phosphorylated tyrosine residue of the regulator SspA positively affects expression and secretion of T3SS proteins and formation of A/E lesions. Altogether, our study reveals that tyrosine phosphorylation in bacteria is more prevalent than previously recognized, and suggests the involvement of phosphotyrosine-mediated signaling in a broad range of cellular functions and virulence.

Solution nuclear magnetic resonance structure of the GATase subunit and structural Basis of the interaction between GATase and ATPPase subunits in a two-subunit-type GMPS from methanocaldococcus jannaschii

The solution structure of the monomeric glutamine amidotransferase (GATase) subunit of the Methanocaldococcus janaschii (Mj) guanosine monophosphate synthetase (GMPS) has been determined using high-resolution nuclear magnetic resonance methods. Gel filtration chromatography and N-15 backbone relaxation studies have shown that the Mj GATase subunit is present in solution as a 21 kDa (188-residue) monomer. The ensemble of 20 lowest-energy structures showed root-mean-square deviations of 0.35 +/- 0.06 angstrom for backbone atoms and 0.8 +/- 0.06 angstrom for all heavy atoms. Furthermore, 99.4% of the backbone dihedral angles are present in the allowed region of the Ramachandran map, indicating the stereochemical quality of the structure. The core of the tertiary structure of the GATase is composed of a seven-stranded mixed beta-sheet that is fenced by five alpha-helices. The Mj GATase is similar in structure to the Pyrococcus horikoshi (Ph) GATase subunit. Nuclear magnetic resonance (NMR) chemical shift perturbations and changes in line width were monitored to identify residues on GATase that were responsible for interaction with magnesium and the ATPPase subunit, respectively. These interaction studies showed that a common surface exists for the metal ion binding as well as for the protein-protein interaction. The dissociation constant for the GATase-Mg2+ interaction has been found to be similar to 1 mM, which implies that interaction is very weak and falls in the fast chemical exchange regime. The GATase-ATPPase interaction, on the other hand, falls in the intermediate chemical exchange regime on the NMR time scale. The implication of this interaction in terms of the regulation of the GATase activity of holo GMPS is discussed.

Probing p300/CBP associated Factor (PCAF)-dependent pathways with a small molecule inhibitor

PCAF (KAT2B) belongs to the GNAT family of lysine acetyltransferases (KAT) and specifically acetylates the histone H3K9 residue and several nonhistone proteins. PCAF is also a transcriptional coactivator. Due to the lack of a PCAF KAT-specific small molecule inhibitor, the exclusive role of the acetyltransferase activity of PCAF is not well understood. Here, we report that a natural compound of the hydroxybenzoquinone class, embelin, specifically inhibits H3Lys9 acetylation in mice and inhibits recombinant PCAF-mediated acetylation with near complete specificity in vitro. Furthermore, using embelin, we have identified the gene networks that are regulated by PCAF during muscle differentiation, further highlighting the broader regulatory functions of PCAF in muscle differentiation in addition to the regulation via MyoD acetylation.

Depth: a web server to compute depth, cavity sizes, detect potential small-molecule ligand-binding cavities and predict the pK(a) of ionizable residues in proteins

Residue depth accurately measures burial and parameterizes local protein environment. Depth is the distance of any atom/residue to the closest bulk water. We consider the non-bulk waters to occupy cavities, whose volumes are determined using a Voronoi procedure. Our estimation of cavity sizes is statistically superior to estimates made by CASTp and VOIDOO, and on par with McVol over a data set of 40 cavities. Our calculated cavity volumes correlated best with the experimentally determined destabilization of 34 mutants from five proteins. Some of the cavities identified are capable of binding small molecule ligands. In this study, we have enhanced our depth-based predictions of binding sites by including evolutionary information. We have demonstrated that on a database (LigASite) of similar to 200 proteins, we perform on par with ConCavity and better than MetaPocket 2.0. Our predictions, while less sensitive, are more specific and precise. Finally, we use depth (and other features) to predict pK(a)s of GLU, ASP, LYS and HIS residues. Our results produce an average error of just <1 pH unit over 60 predictions. Our simple empirical method is statistically on par with two and superior to three other methods while inferior to only one. The DEPTH server (http://mspc.bii.a-star.edu.sg/depth/) is an ideal tool for rapid yet accurate structural analyses of protein structures.

Mechanistic features of Salmonella typhimurium propionate kinase (TdcD): Insights from kinetic and crystallographic studies

Short-chain fatty acids (SCFAs) play a major role in carbon cycle and can be utilized as a source of carbon and energy by bacteria. Salmonella typhimurium propionate kinase (StTdcD) catalyzes reversible transfer of the gamma-phosphate of ATP to propionate during L-threonine degradation to propionate. Kinetic analysis revealed that StTdcD possesses broad ligand specificity and could be activated by various SCFAs (propionate > acetate approximate to butyrate), nucleotides (ATP approximate to GTP > CTP approximate to TTP; dATP > dGTP > dCTP) and metal ions (Mg2+ approximate to Mn2+ > Co2+). Inhibition of StTdcD by tricarboxylic acid (TCA) cycle intermediates such as citrate, succinate, alpha-ketoglutarate and malate suggests that the enzyme could be under plausible feedback regulation. Crystal structures of StTdcD bound to PO4 (phosphate), AMP, ATP, Ap4 (adenosine tetraphosphate), GMP, GDP, GTP, CMP and CTP revealed that binding of nucleotide mainly involves hydrophobic interactions with the base moiety and could account for the broad biochemical specificity observed between the enzyme and nucleotides. Modeling and site-directed mutagenesis studies suggest Ala88 to be an important residue involved in determining the rate of catalysis with SCFA substrates. Molecular dynamics simulations on monomeric and dimeric forms of StTdcD revealed plausible open and closed states, and also suggested role for dimerization in stabilizing segment 235-290 involved in interfacial interactions and ligand binding. Observation of an ethylene glycol molecule bound sufficiently close to the gamma-phosphate in StTdcD complexes with triphosphate nucleotides supports direct in-line phosphoryl transfer. (C) 2013 Elsevier B.V. All rights reserved.

An automated approach to network features of protein structure ensembles

Network theory applied to protein structures provides insights into numerous problems of biological relevance. The explosion in structural data available from PDB and simulations establishes a need to introduce a standalone-efficient program that assembles network concepts/parameters under one hood in an automated manner. Herein, we discuss the development/application of an exhaustive, user-friendly, standalone program package named PSN-Ensemble, which can handle structural ensembles generated through molecular dynamics (MD) simulation/NMR studies or from multiple X-ray structures. The novelty in network construction lies in the explicit consideration of side-chain interactions among amino acids. The program evaluates network parameters dealing with topological organization and long-range allosteric communication. The introduction of a flexible weighing scheme in terms of residue pairwise cross-correlation/interaction energy in PSN-Ensemble brings in dynamical/chemical knowledge into the network representation. Also, the results are mapped on a graphical display of the structure, allowing an easy access of network analysis to a general biological community. The potential of PSN-Ensemble toward examining structural ensemble is exemplified using MD trajectories of an ubiquitin-conjugating enzyme (UbcH5b). Furthermore, insights derived from network parameters evaluated using PSN-Ensemble for single-static structures of active/inactive states of 2-adrenergic receptor and the ternary tRNA complexes of tyrosyl tRNA synthetases (from organisms across kingdoms) are discussed. PSN-Ensemble is freely available from http://vishgraph.mbu.iisc.ernet.in/PSN-Ensemble/psn_index.html.

Identification of key amino acid residues in the catalytic mechanism of diaminopropionate ammonialyase from Salmonella typhimurium

Diaminopropionate ammonialyase (DAPAL), a fold-typeII pyridoxal 5-phosphate-dependent enzyme, catalyzes the ,-elimination of diaminopropionate (DAP) to pyruvate and ammonia. DAPAL was able to utilize both d- and l-DAP as substrates with almost equal efficiency. Mutational analysis of functionally important residues such as Thr385, Asp125 and Asp194 was carried out to understand the mechanism by which the isomers are hydrolyzed. Further, the putative residues involved in the formation of disulfide bond Cys271 and Cys299 were also mutated. T385S, T385D sDAPAL were as active with dl-DAP as substrate as sDAPAL, whereas the later exhibited a threefold increase in catalytic efficiency with d-Ser as substrate. Further analysis of these mutants suggested that DAPAL might follow an anti-E-2 mechanism of catalysis that does not involve the formation of a quinonoid intermediate. Of the two mutants of Asp125, D125E showed complete loss of activity with d-DAP as substrate, whereas the reaction with l-DAP was not affected significantly, demonstrating that Asp125 was essential for abstraction of protons from the d-isomer. By contrast, mutational analysis of Asp194 showed that the residue may not be directly involved in proton abstraction from l-DAP. sDAPAL does not form a disulfide bond in solution, although the position of Cys299 and Cys271 in the modeled structure of sDAPAL favored the formation of a disulfide bond. Further, unlike eDAPAL, sDAPAL could be activated by monovalent cations. Mutation of the cysteine residues showed that Cys271 may be involved in coordinating the monovalent cation, as observed in the case of other fold-typeII enzymes.

Increasing cleavage specificity and activity of restriction endonuclease KpnI

Restriction enzyme KpnI is a HNH superfamily endonuclease requiring divalent metal ions for DNA cleavage but not for binding. The active site of KpnI can accommodate metal ions of different atomic radii for DNA cleavage. Although Mg2+ ion higher than 500 mu M mediates promiscuous activity, Ca2+ suppresses the promiscuity and induces high cleavage fidelity. Here, we report that a conservative mutation of the metal-coordinating residue D148 to Glu results in the elimination of the Ca2+-mediated cleavage but imparting high cleavage fidelity with Mg2+. High cleavage fidelity of the mutant D148E is achieved through better discrimination of the target site at the binding and cleavage steps. Biochemical experiments and molecular dynamics simulations suggest that the mutation inhibits Ca2+-mediated cleavage activity by altering the geometry of the Ca2+-bound HNH active site. Although the D148E mutant reduces the specific activity of the enzyme, we identified a suppressor mutation that increases the turnover rate to restore the specific activity of the high fidelity mutant to the wild-type level. Our results show that active site plasticity in coordinating different metal ions is related to KpnI promiscuous activity, and tinkering the metal ion coordination is a plausible way to reduce promiscuous activity of metalloenzymes.

A multi-level multi-scale approach to study essential genes in Mycobacterium tuberculosis

Background: The set of indispensable genes that are required by an organism to grow and sustain life are termed as essential genes. There is a strong interest in identification of the set of essential genes, particularly in pathogens, not only for a better understanding of the pathogen biology, but also for identifying drug targets and the minimal gene set for the organism. Essentiality is inherently a systems property and requires consideration of the system as a whole for their identification. The available experimental approaches capture some aspects but each method comes with its own limitations. Moreover, they do not explain the basis for essentiality in most cases. A powerful prediction method to recognize this gene pool including rationalization of the known essential genes in a given organism would be very useful. Here we describe a multi-level multi-scale approach to identify the essential gene pool in a deadly pathogen, Mycobacterium tuberculosis. Results: The multi-level workflow analyses the bacterial cell by studying (a) genome-wide gene expression profiles to identify the set of genes which show consistent and significant levels of expression in multiple samples of the same condition, (b) indispensability for growth by using gene expression integrated flux balance analysis of a genome-scale metabolic model, (c) importance for maintaining the integrity and flow in a protein-protein interaction network and (d) evolutionary conservation in a set of genomes of the same ecological niche. In the gene pool identified, the functional basis for essentiality has been addressed by studying residue level conservation and the sub-structure at the ligand binding pockets, from which essential amino acid residues in that pocket have also been identified. 283 genes were identified as essential genes with high-confidence. An agreement of about 73.5% is observed with that obtained from the experimental transposon mutagenesis technique. A large proportion of the identified genes belong to the class of intermediary metabolism and respiration. Conclusions: The multi-scale, multi-level approach described can be generally applied to other pathogens as well. The essential gene pool identified form a basis for designing experiments to probe their finer functional roles and also serve as a ready shortlist for identifying drug targets.

Crystal structures of SCP2-thiolases of Trypanosomatidae, human pathogens causing widespread tropical diseases: the importance for catalysis of the cysteine of the unique HDCF loop

Thiolases are essential CoA-dependent enzymes in lipid metabolism. In the present study we report the crystal structures of trypanosomal and leishmanial SCP2 (sterol carrier protein, type-2)-thiolases. Trypanosomatidae cause various widespread devastating (sub)-tropical diseases, for which adequate treatment is lacking. The structures reveal the unique geometry of the active site of this poorly characterized subfamily of thiolases. The key catalytic residues of the classical thiolases are two cysteine residues, functioning as a nucleophile and an acid/base respectively. The latter cysteine residue is part of a CxG motif. Interestingly, this cysteine residue is not conserved in SCP2-thiolases. The structural comparisons now show that in SCP2-thiolases the catalytic acid/base is provided by the cysteine residue of the HDCF motif, which is unique for this thiolase subfamily. This HDCF cysteine residue is spatially equivalent to the CxG cysteine residue of classical thiolases. The HDCF cysteine residue is activated for acid/base catalysis by two main chain NH-atoms, instead of two water molecules, as present in the CxG active site. The structural results have been complemented with enzyme activity data, confirming the importance of the HDCF cysteine residue for catalysis. The data obtained suggest that these trypanosomatid SCP2-thiolases are biosynthetic thiolases. These findings provide promise for drug discovery as biosynthetic thiolases catalyse the first step of the sterol biosynthesis pathway that is essential in several of these parasites.

Investigations of Ramachandran disallowed conformations in protein domain families

In peptide and protein structures, occurrence of (phi,psi.) angles in the disallowed region of the Ramachandran map almost always suggests local regions of error or poor accuracy. However, very rarely genuine disallowed conformations occur as noted in the current study in proteins of known structure available at ultra-high resolution (<= 1.2 (A) over circle). In the current work, extent of conservation of genuine disallowed conformations in homologous proteins of known structures has been analyzed. From a dataset of 124 protein domain families, with structure of at least one constituent member in each family available at a resolution of 1.2 (A) over circle or better, we have analyzed the conservation of 221 disallowed conformations. It is observed that the disallowed conformation is only moderately conservedin protein domain families. In the gross dataset no particular residue type adopting disallowed conformation elicit high conservation of residue type though there are alignment positions in the dataset with complete conservation of both the residue type and the disallowed conformation. Conserved disallowed conformation in protein domain families play biologically significant role in roughly 50% of the cases. The residues with the disallowed conformation or its flanking residues are often located within or around the functional site of the protein. (C) 2013 Elsevier B.V. All rights reserved.