Essential dynamics and sidechain hydrogen bond cluster studies on eosinophil cationic protein

Eosinophil Cationic Protein (ECP) is a member of RNase A superfamily which carries out the obligatory catalytic role of cleaving RNA. It is involved in a variety of biological functions. Molecular dynamics simulations followed by essential dynamics analysis on this protein are carried out with the goal of gaining insights into the dynamical properties at atomic level. The top essential modes contribute to subspaces and to the transition phase. Further, the sidechain-sidechain/sidechain-mainchain hydrogen bond clusters are analyzed in the top modes, and compared with those of crystal structure. The role of residues identified by these methods is discussed in the context of concerted motion, structure and stability of the protein.

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Endonuclease Active Site Plasticity Allows DNA Cleavage with Diverse Alkaline Earth and Transition Metal Ions

A majority of enzymes show a high degree of specificity toward a particular metal ion in their catalytic reaction. However, Type II restriction endonuclease (REase) R.KpnI, which is the first member of the HNH superfamily of REases, exhibits extraordinary diversity in metal ion dependent DNA cleavage. Several alkaline earth and transition group metal ions induce high fidelity and promiscuous cleavage or inhibition depending upon their concentration. The metal ions having different ionic radii and co-ordination geometries readily replace each other from the enzyme's active site, revealing its plasticity. Ability of R KpnI to cleave DNA with both alkaline earth and transition group metal ions having varied ionic radii could imply utilization of different catalytic site(s). However, mutation of the invariant His residue of the HNH motif caused abolition of the enzyme activity with all of the cofactors, indicating that the enzyme follows a single metal ion catalytic mechanism for DNA cleavage. Indispensability of His in nucleophile activation together with broad cofactor tolerance of the enzyme indicates electrostatic stabilization function of metal ions during catalysis. Nevertheless, a second metal ion is recruited at higher concentrations to either induce promiscuity or inhibit the DNA cleavage. Regulation of the endonuclease activity and fidelity by a second metal ion binding is a unique feature of R.KpnI among REases and HNH nucleases. The active site plasticity of R.KpnI opens up avenues for redesigning cofactor specificities and generation of mutants specific to a particular metal ion.

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background: Bacteria such as Escherichia coli and Salmonella typhimurium can utilize acetate as the sole source of carbon and energy. Acetate kinase (AckA) and phosphotransacetylase (Pta), key enzymes of acetate utilization pathway, regulate flux of metabolites in glycolysis, gluconeogenesis, TCA cycle, glyoxylate bypass and fatty acid metabolism. Results: Here we report kinetic characterization of S. typhimurium AckA (StAckA) and structures of its unliganded (Form-I, 2.70 angstrom resolution) and citrate-bound (Form-II, 1.90 angstrom resolution) forms. The enzyme showed broad substrate specificity with k(cat)/K-m in the order of acetate > propionate > formate. Further, the K-m for acetyl-phosphate was significantly lower than for acetate and the enzyme could catalyze the reverse reaction (i.e. ATP synthesis) more efficiently. ATP and Mg2+ could be substituted by other nucleoside 5'-triphosphates (GTP, UTP and CTP) and divalent cations (Mn2+ and Co2+), respectively. Form-I StAckA represents the first structural report of an unliganded AckA. StAckA protomer consists of two domains with characteristic beta beta beta alpha beta alpha beta alpha topology of ASKHA superfamily of proteins. These domains adopt an intermediate conformation compared to that of open and closed forms of ligand-bound Methanosarcina thermophila AckA (MtAckA). Spectroscopic and structural analyses of StAckA further suggested occurrence of inter-domain motion upon ligand-binding. Unexpectedly, Form-II StAckA structure showed a drastic change in the conformation of residues 230-300 compared to that of Form-I. Further investigation revealed electron density corresponding to a citrate molecule in a pocket located at the dimeric interface of Form-II StAckA. Interestingly, a similar dimeric interface pocket lined with largely conserved residues could be identified in Form-I StAckA as well as in other enzymes homologous to AckA suggesting that ligand binding at this pocket may influence the function of these enzymes. Conclusions: The biochemical and structural characterization of StAckA reported here provides insights into the biochemical specificity, overall fold, thermal stability, molecular basis of ligand binding and inter-domain motion in AckA family of enzymes. Dramatic conformational differences observed between unliganded and citrate-bound forms of StAckA led to identification of a putative ligand-binding pocket at the dimeric interface of StAckA with implications for enzymatic function.

Improved Detection of Remote Homologues Using Cascade PSI-BLAST: Influence of Neighbouring Protein Families on Sequence Coverage

Background: Development of sensitive sequence search procedures for the detection of distant relationships between proteins at superfamily/fold level is still a big challenge. The intermediate sequence search approach is the most frequently employed manner of identifying remote homologues effectively. In this study, examination of serine proteases of prolyl oligopeptidase, rhomboid and subtilisin protein families were carried out using plant serine proteases as queries from two genomes including A. thaliana and O. sativa and 13 other families of unrelated folds to identify the distant homologues which could not be obtained using PSI-BLAST. Methodology/Principal Findings: We have proposed to start with multiple queries of classical serine protease members to identify remote homologues in families, using a rigorous approach like Cascade PSI-BLAST. We found that classical sequence based approaches, like PSI-BLAST, showed very low sequence coverage in identifying plant serine proteases. The algorithm was applied on enriched sequence database of homologous domains and we obtained overall average coverage of 88% at family, 77% at superfamily or fold level along with specificity of similar to 100% and Mathew's correlation coefficient of 0.91. Similar approach was also implemented on 13 other protein families representing every structural class in SCOP database. Further investigation with statistical tests, like jackknifing, helped us to better understand the influence of neighbouring protein families. Conclusions/Significance: Our study suggests that employment of multiple queries of a family for the Cascade PSI-BLAST searches is useful for predicting distant relationships effectively even at superfamily level. We have proposed a generalized strategy to cover all the distant members of a particular family using multiple query sequences. Our findings reveal that prior selection of sequences as query and the presence of neighbouring families can be important for covering the search space effectively in minimal computational time. This study also provides an understanding of the `bridging' role of related families.

Mapping the apoptosis inducing domain of an immunomodulatory protein: glycodelin A

Glycodelin A (GdA) is a dimeric glycoprotein synthesized by the human endometrium under progesterone regulation. Based on the high sequence similarity with beta-lactoglobulin, it is placed under the lipocalin superfamily. The protein is one of the local immunomodulators present at the feto-maternal interface which affects both the innate as well as the acquired arms of the immune system, thereby bringing about successful establishment and progression of pregnancy. Our previous studies revealed that the domain responsible for the immunosuppressive activity of glycodelin lies on its protein backbone and the glycans modulate the same. This study attempts to further delineate the apoptosis inducing region of GdA. Our results demonstrate that the stretch of amino acid sequence between Met24 to Leu105 is necessary and sufficient to inhibit proliferation of T cells and induce apoptosis in them. Further, within this region the key residues involved in harboring the activity were shown to be present between Asp52 and Ser65.

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.

Conformational Analysis of a 20-Membered Cyclic Peptide Disulfide from Conus virgo with a WPW Segment: Evidence for an Aromatic-Proline Sandwich

A novel peptide containing a single disulfide bond, CIWPWC (Vi804), has been isolated and characterised from the venom of the marine cone snail, Conus virgo. A precursor polypeptide sequence derived from complementary DNA, corresponding to the M-superfamily conotoxins, has been identified. The identity of the synthetic and natural peptide sequence has been established. A detailed analysis of the conformation in solution is reported for Vi804 and a synthetic analogue, (CIWPWC)-W-D ((D)W3-Vi804), in order to establish the structure of the novel WPW motif, which occurs in the context of a 20-membered macrocyclic disulfide. Vi804 exists exclusively in the cis W3P4 conformer in water and methanol, whereas (D)W3-Vi804 occurs exclusively as the trans conformer. NMR spectra revealed a W3P4 typeVI turn in Vi804 and a typeII turn in the analogue peptide, (D)W3-Vi804. The extremely high-field chemical shifts of the proline ring protons, together with specific nuclear Overhauser effects, are used to establish a conformation in which the proline ring is sandwiched between the flanking Trp residues, which emphasises a stabilising role for the aromatic-proline interactions, mediated predominantly by dispersion forces.

The Non-catalytic ``Cap Domain'' of a Mycobacterial Metallophosphoesterase Regulates Its Expression and Localization in the Cell

Despite highly conserved core catalytic domains, members of the metallophosphoesterase (MPE) superfamily perform diverse and crucial functions ranging from nucleotide and nucleic acid metabolism to phospholipid hydrolysis. Unique structural elements outside of the catalytic core called ``cap domains'' are thought to provide specialization to these enzymes; however, no directed study has been performed to substantiate this. The cap domain of Rv0805, an MPE from Mycobacterium tuberculosis, is located C-terminal to its catalytic domain and is dispensable for the catalytic activity of this enzyme in vitro. We show here that this C-terminal extension (CTE) mediates in vivo localization of the protein to the cell membrane and cell wall as well as modulates expression levels of Rv0805 in mycobacteria. We also demonstrate that Rv0805 interacts with the cell wall of mycobacteria, possibly with the mycolyl-arabinogalactan-peptidoglycan complex, by virtue of its C terminus, a hitherto unknown property of this MPE. Using a panel of mutant proteins, we identify interactions between active site residues of Rv0805 and the CTE that determine its association with the cell wall. Finally, we show that Rv0805 and a truncated mutant devoid of the CTE produce different phenotypic effects when expressed in mycobacteria. Our study thus provides a detailed dissection of the functions of the cap domain of an MPE and suggests that the repertoire of cellular functions of MPEs cannot be understood without exploring the modulatory effects of these subdomains.

Enriching the annotation of Mycobacterium tuberculosis H37Rv proteome using remote homology detection approaches: Insights into structure and function

The availability of the genome sequence of Mycobacterium tuberculosis H37Rv has encouraged determination of large numbers of protein structures and detailed definition of the biological information encoded therein; yet, the functions of many proteins in M. tuberculosis remain unknown. The emergence of multidrug resistant strains makes it a priority to exploit recent advances in homology recognition and structure prediction to re-analyse its gene products. Here we report the structural and functional characterization of gene products encoded in the M. tuberculosis genome, with the help of sensitive profile-based remote homology search and fold recognition algorithms resulting in an enhanced annotation of the proteome where 95% of the M. tuberculosis proteins were identified wholly or partly with information on structure or function. New information includes association of 244 proteins with 205 domain families and a separate set of new association of folds to 64 proteins. Extending structural information across uncharacterized protein families represented in the M. tuberculosis proteome, by determining superfamily relationships between families of known and unknown structures, has contributed to an enhancement in the knowledge of structural content. In retrospect, such superfamily relationships have facilitated recognition of probable structure and/or function for several uncharacterized protein families, eventually aiding recognition of probable functions for homologous proteins corresponding to such families. Gene products unique to mycobacteria for which no functions could be identified are 183. Of these 18 were determined to be M. tuberculosis specific. Such pathogen-specific proteins are speculated to harbour virulence factors required for pathogenesis. A re-annotated proteome of M. tuberculosis, with greater completeness of annotated proteins and domain assigned regions, provides a valuable basis for experimental endeavours designed to obtain a better understanding of pathogenesis and to accelerate the process of drug target discovery. (C) 2014 Elsevier Ltd. All rights reserved.

Metallophosphoesterases: structural fidelity with functional promiscuity

Calcineurin-like metallophosphoesterases (MPEs) form a large superfamily of binuclear metal-ion-centre-containing enzymes that hydrolyse phosphomono-, phosphodi-or phosphotri-esters in a metal-dependent manner. The MPE domain is found in Mre11/SbcD DNA-repair enzymes, mammalian phosphoprotein phosphatases, acid sphingomyelinases, purple acid phosphatases, nucleotidases and bacterial cyclic nucleotide phosphodiesterases. Despite this functional diversity, MPEs show a remarkably similar structural fold and active-site architecture. In the present review, we summarize the available structural, biochemical and functional information on these proteins. We also describe how diversification and specialization of the core MPE fold in various MPEs is achieved by amino acid substitution in their active sites, metal ions and regulatory effects of accessory domains. Finally, we discuss emerging roles of these proteins as non-catalytic protein-interaction scaffolds. Thus we view the MPE superfamily as a set of proteins with a highly conserved structural core that allows embellishment to result in dramatic and niche-specific diversification of function.

A sleep-inducing peptide from the venom of the Indian cone snail Conus araneosus

The marine snail Conus araneosus has unusual significance due to its confined distribution to coastal regions of southeast India and Sri Lanka. Due to its relative scarceness, this species has been poorly studied. In this work, we characterized the venom of C. araneosus to identify new venom peptides. We identified 14 novel compounds. We determined amino acid sequences from chemically-modified and unmodified crude venom using liquid chromatography-electrospray ionization mass spectrometry and matrix assisted laser desorption ionization time-of-flight mass spectrometry. Ten sequences showed six Cys residues arranged in a pattern that is most commonly associated with the M-superfamily of conotoxins. Four other sequences had four Cys residues in a pattern that is most commonly associated with the T-superfamily of conotoxins. The post-translationally modified residue (pyroglutamate) was determined at the N-terminus of two sequences, ar3h and ar3i respectively. In addition, two sequences, ar3g and ar3h were C-terminally amidated. At a dose of 2 nmol, peptide ar3j elicited sleep when injected intraperitoneally into mice. To our knowledge, this is the first report of a peptide from a molluscivorous cone snail with sleep-inducing effects in mice. The novel peptides characterized herein extend the repertoire of unique peptides derived from cone snails and may add value to the therapeutic promise of conotoxins. (C) 2015 Elsevier Ltd. All rights reserved.

Structures of substrate- and nucleotide-bound propionate kinase from Salmonella typhimurium: substrate specificity and phosphate-transfer mechanism

Kinases are ubiquitous enzymes that are pivotal to many biochemical processes. There are contrasting views on the phosphoryl-transfer mechanism in propionate kinase, an enzyme that reversibly transfers a phosphoryl group from propionyl phosphate to ADP in the final step of non-oxidative catabolism of L-threonine to propionate. Here, X-ray crystal structures of propionate- and nucleotide-bound Salmonella typhimurium propionate kinase are reported at 1.8-2.0 angstrom resolution. Although the mode of nucleotide binding is comparable to those of other members of the ASKHA superfamily, propionate is bound at a distinct site deeper in the hydrophobic pocket defining the active site. The propionate carboxyl is at a distance of approximate to 5 angstrom from the -phosphate of the nucleotide, supporting a direct in-line transfer mechanism. The phosphoryl-transfer reaction is likely to occur via an associative S(N)2-like transition state that involves a pentagonal bipyramidal structure with the axial positions occupied by the nucleophile of the substrate and the O atom between the - and the -phosphates, respectively. The proximity of the strictly conserved His175 and Arg236 to the carboxyl group of the propionate and the -phosphate of ATP suggests their involvement in catalysis. Moreover, ligand binding does not induce global domain movement as reported in some other members of the ASKHA superfamily. Instead, residues Arg86, Asp143 and Pro116-Leu117-His118 that define the active-site pocket move towards the substrate and expel water molecules from the active site. The role of Ala88, previously proposed to be the residue determining substrate specificity, was examined by determining the crystal structures of the propionate-bound Ala88 mutants A88V and A88G. Kinetic analysis and structural data are consistent with a significant role of Ala88 in substrate-specificity determination. The active-site pocket-defining residues Arg86, Asp143 and the Pro116-Leu117-His118 segment are also likely to contribute to substrate specificity.

Draft genome of a commonly misdiagnosed multidrug resistant pathogen Candida auris

Background: Candida auris is a multidrug resistant, emerging agent of fungemia in humans. Its actual global distribution remains obscure as the current commercial methods of clinical diagnosis misidentify it as C. haemulonii. Here we report the first draft genome of C. auris to explore the genomic basis of virulence and unique differences that could be employed for differential diagnosis. Results: More than 99.5 % of the C. auris genomic reads did not align to the current whole (or draft) genome sequences of Candida albicans, Candida lusitaniae, Candida glabrata and Saccharomyces cerevisiae; thereby indicating its divergence from the active Candida clade. The genome spans around 12.49 Mb with 8527 predicted genes. Functional annotation revealed that among the sequenced Candida species, it is closest to the hemiascomycete species Clavispora lusitaniae. Comparison with the well-studied species Candida albicans showed that it shares significant virulence attributes with other pathogenic Candida species such as oligopeptide transporters, mannosyl transfersases, secreted proteases and genes involved in biofilm formation. We also identified a plethora of transporters belonging to the ABC and major facilitator superfamily along with known MDR transcription factors which explained its high tolerance to antifungal drugs. Conclusions: Our study emphasizes an urgent need for accurate fungal screening methods such as PCR and electrophoretic karyotyping to ensure proper management of fungemia. Our work highlights the potential genetic mechanisms involved in virulence and pathogenicity of an important emerging human pathogen namely C. auris. Owing to its diversity at the genomic scale; we expect the genome sequence to be a useful resource to map species specific differences that will help develop accurate diagnostic markers and better drug targets.

Differential dynamics of the serotonin(1A) receptor in membrane bilayers of varying cholesterol content revealed by all atom molecular dynamics simulation

The serotonin(1A) receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin(1A) receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3s, to analyze the effect of cholesterol on the structure and dynamics of the serotonin(1A) receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin(1A) receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.