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.

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.

Corticosterone targets distinct steps of synaptic transmission via concentration specific activation of mineralocorticoid and glucocorticoid receptors

Hippocampal neurons are affected by chronic stress and have a high density of cytoplasmic mineralocorticoid and glucocorticoid receptors (MR and GR). Detailed studies on the genomic effects of the stress hormone corticosterone at physiologically relevant concentrations on different steps in synaptic transmission are limited. In this study, we tried to delineate how activation of MR and GR by different concentrations of corticosterone affects synaptic transmission at various levels. The effect of 3-h corticosterone (25, 50, and 100nM) treatment on depolarization-mediated calcium influx, vesicular release and properties of miniature excitatory post-synaptic currents (mEPSCs) were studied in cultured hippocampal neurons. Activation of MR with 25nM corticosterone treatment resulted in enhanced depolarization-mediated calcium influx via a transcription-dependent process and increased frequency of mEPSCs with larger amplitude. On the other hand, activation of GR upon 100nM corticosterone treatment resulted in increase in the rate of vesicular release via the genomic actions of GR. Furthermore, GR activation led to significant increase in the frequency of mEPSCs with larger amplitude and faster decay. Our studies indicate that differential activation of the dual receptor system of MR and GR by corticosterone targets the steps in synaptic transmission differently.

Structural and kinetic studies on adenylosuccinate lyase from Mycobacterium smegmatis and Mycobacterium tuberculosis provide new insights on the catalytic residues of the enzyme

Adenylosuccinate lyase (ASL), an enzyme involved in purine biosynthesis, has been recognized as a drug target against microbial infections. In the present study, ASL from Mycobacteriumsmegmatis (MsASL) and Mycobacteriumtuberculosis (MtbASL) were cloned, purified and crystallized. The X-ray crystal structure of MsASL was determined at a resolution of 2.16 angstrom. It is the first report of an apo-ASL structure with a partially ordered active site C3 loop. Diffracting crystals of MtbASL could not be obtained and a model for its structure was derived using MsASL as a template. These structures suggest that His149 and either Lys285 or Ser279 of MsASL are the residues most likely to function as the catalytic acid and base, respectively. Most of the active site residues were found to be conserved, with the exception of Ser148 and Gly319 of MsASL. Ser148 is structurally equivalent to a threonine in most other ASLs. Gly319 is replaced by an arginine residue in most ASLs. The two enzymes were catalytically much less active compared to ASLs from other organisms. Arg319Gly substitution and reduced flexibility of the C3 loop might account for the low catalytic activity of mycobacterial ASLs. The low activity is consistent with the slow growth rate of Mycobacteria and their high GC containing genomes, as well as their dependence on other salvage pathways for the supply of purine nucleotides. Structured digital abstract andby()

C-12 Helices in Long Hybrid (alpha gamma)(n) Peptides Composed Entirely of Unconstrained Residues with Proteinogenic Side Chains

Unconstrained gamma(4) amino acid residues derived by homologation of proteinogenic amino acids facilitate helical folding in hybrid (alpha gamma)(n) sequences. The C-12 helical conformation for the decapeptide, Boc-Leu-gamma(4)(R)Val](5)-OMe, is established in crystals by X-ray diffraction. A regular C-12 helix is demonstrated by NMR studies of the 18 residue peptide, Boc-Leu-gamma(4)(AR)Val](9)-OMe, and a designed 16 residue (alpha gamma)(n) peptide, incorporating variable side chains. Unconstrained (alpha gamma)(n) peptides show an unexpectedly high propensity for helical folding in long polypeptide sequences.

Self-assembled dipeptide nanotubes constituted by flexible beta-phenylalanine and conformationally constrained alpha,beta-dehydrophenylalanine residues as drug delivery system

Peptide based self assembled nanostructures have attracted growing interest in recent years due to their numerous potential applications particularly in biomedical sciences. Di-peptide Phe-Phe was shown previously to self-assemble into nanotube like structures. In this work, we studied the affect of peptide backbone length and conformational flexibility on the self assembly process by using two dipeptides based on the Phe-Phe backbone (beta Phe-Phe and beta Phe-Delta Phe): one containing a flexible beta Phe amino acid, and the other containing both a flexible bPhe as well as a backbone constraining Alpha Phe (alpha,beta-dehydrophenylalanine) amino acid. Electron microscopy and X-ray diffraction experiments revealed that these new di-peptides can self-assemble into nanotubes having different properties than the native Phe-Phe nanotubes. These nanotubes were stable over a broad range of temperatures and the introduction of non-natural amino acids provided them with stability against the action of nonspecific proteases. Moreover, these dipeptides showed no cytotoxicity towards HeLa and L929 cells, and were able to encapsulate small drug molecules. We further showed that anticancerous drug mitoxantrone was more efficient in killing HeLa and B6F10 cells when entrapped in nanotubes as compared to free mitoxantrone. Therefore, these beta-phenylalanine and alpha, beta-dehydrophenylalanine containing dipeptide nanotubes may be useful in the development of biocompatible and proteolytically stable drug delivery vehicles.

Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety-six of these alignment positions correspond to conserved Gly residue with (phi, ) values allowed for non-glycyl residues. Reasons for this observation were investigated by in-silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (phi, ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (phi, ) values which are disallowed for Ala. In-silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.

NON-GLYCINE RESIDUES IN THE DISALLOWED REGIONS OF THE RAMACHANDRAN MAP AND THEIR NEIGHBORHOOD PREFERENCES

The allowed and the ``disallowed'' regions in the celebrated Ramachandran map (phi-psi] map) was elegantly deduced by Ramachandran, Ramakrishnan and Sasisekharan even before the protein crystal structures became available. This powerful map was derived based on rigid geometry of the peptide group and later several investigations on protein crystal structures reported the occurrence of a small fraction of the phi-psi] torsion angles in the disallowed region. The question is what factors make these residues adopt disallowed conformations? Is it driven by the necessity to maintain the overall topology or is it associated with function or is it just that the disallowed conformations are extreme limits of the allowed conformations? Today, with the availability of a large number of high resolution crystal structures, we have revisited this problem. Apart from validating some of the earlier findings such as residue propensities, preferred location in the secondary structure, we have explored their spatial neighborhood preferences using the protein structure network PSN] approach developed in our lab. Finally, the structural and functional implications of the disallowed conformations are examined.

PROPENSITIES OF AMINO ACID RESIDUES IN PROTEINS FOR DIFFERENT REGIONS OF THE RAMACHANDRAN MAP

Using a dataset of 1164 crystal structures of largely non-homologous proteins defined at a resolution of 1.5 angstrom or better, we have investigated the (phi,psi) preferences of 20 residue types by considering the residues which occur in loops. Propensities of residue types to occur in the loops with (phi,psi) values in the aa region of the Ramachandran map has a poor correlation coefficient of 0.48 to the Chou-Fasman propensities of the residue types to occur in the a-helical segments. However the correlation coefficient between propensities of residues in loops to adopt beta conformations and those in beta-sheet is much higher (0.95). These observations suggest that a-helix formation is well influenced by the local amino acid sequence while intrinsic preference of residue types for beta-sheet plays a major role in the formation of beta-sheet. The main chain polar groups of residues in loops, that can affect the (phi,psi) values, can be involved in intra-molecular hydrogen bonding. Therefore we investigated further by considering subset of residues in loops with low (0 to 2) number of intra-molecular hydrogen bonds per residue involving main chain polar atoms. For this subset, the correlation coefficients between propensities for alpha-helix and alpha(R) region and between beta-sheet and beta-region are 0.26 and 0.64 respectively. This reiterates higher intrinsic tendency of beta-region favouring residues to adopt beta-sheet than alpha(R) region favouring residues to adopt alpha-helical structure.

INDUCTION OF FOLDED STRUCTURES IN DESIGNED PEPTIDES USING CONFORMATIONALLY CONSTRAINED RESIDUES

Folding into compact globular structures, with well-defined modules of secondary structure, appears to be a characteristic of long polypeptide chains, with a specific patterning of coded amino acid residues along the length of sequence. Cooperative hydrogen bond driven secondary structure formation and solvent forces, which contribute favorably to the entropy of folding, by promoting compaction of the polymeric chain, have long been discussed as major determinants of the folding process. First principles design approaches, which use non-coded amino acids, employ an alternative structure directing strategy, by using amino acid residues which exhibit a strong conformational bias for specific regions of the Ramachandran map. This overview of ongoing studies in the authors' laboratory, attempts to explore the use of conformationally restricted amino acid residues in the design of peptides with well-defined secondary structures. Short peptides composed of 20 genetically coded amino acids usually exist in solution as an ensemble of equilibrating conformations. Apolar peptide sequences, which are readily soluble in organic solvents like chloroform and methanol, facilitate formation of structures which are predominately driven by intramolecular hydrogen bond formation. The choice of sequences containing residues with a limited range of conformational choices strongly favors formation of local turn structures, stabilized by short range intramolecular hydrogen bonds. Two residue beta-turns can nucleate either helical or hairpin folding, depending on the precise conformation of the turn segment Restriction of the conformational space available to amino acid residues is easily achieved by introduction of an additional alkyl group at the C alpha carbon atom or by side chain backbone cyclization, as in proline. Studies of synthetic sequences incorporating two prototype residues alpha-aminoisobutyric acid (Aib) and D-proline (DPro) illustrate the utility of the strategy in construction of helices and hairpins. Extensions to the design of conformationally switchable sequences and structurally defined hybrid peptides containing backbone homologated residues are also surveyed.

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.

Impact of Mutating the Key Residues of a Bifunctional 5,10-Methylenetetrahydrofolate Dehydrogenase-Cyclohydrolase from Escherichia coli on Its Activities

Methylenetetrahydrofolate dehydrogenase-cyclohydrolase (FolD) catalyzes interconversion of 5,10-methylene-tetrahydrofolate and 10-formyl-tetrahydrofolate in the one-carbon metabolic pathway. In some organisms, the essential requirement of 10-formyl-tetrahydrofolate may also be fulfilled by formyltetrahydrofolate synthetase (Fhs). Recently, we developed an Escherichia coli strain in which the folD gene was deleted in the presence of Clostridium perfringens fhs (E. coli Delta folD/p-fhs) and used it to purify FolD mutants (free from the host-encoded FolD) and determine their biological activities. Mutations in the key residues of E. coli FolD, as identified from three-dimensional structures (D121A, Q98K, K54S, Y50S, and R191E), and a genetic screen (G122D and C58Y) were generated, and the mutant proteins were purified to determine their kinetic constants. Except for the R191E and K54S mutants, others were highly compromised in terms of both dehydrogenase and cyclohydrolase activities. While the R191E mutant showed high cyclohydrolase activity, it retained only a residual dehydrogenase activity. On the other hand, the K54S mutant lacked the cyclohydrolase activity but possessed high dehydrogenase activity. The D121A and G122D (in a loop between two helices) mutants were highly compromised in terms of both dehydrogenase and cyclohydrolase activities. In vivo and in vitro characterization of wild-type and mutant (R191E, G122D, D121A, Q98K, C58Y, K54S, and Y50S) FolD together with three-dimensional modeling has allowed us to develop a better understanding of the mechanism for substrate binding and catalysis by E. coli FolD.

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.

Bioengineering of AAV2 Capsid at Specific Serine, Threonine, or Lysine Residues Improves Its Transduction Efficiency in Vitro and in Vivo

We hypothesized that the AAV2 vector is targeted for destruction in the cytoplasm by the host cellular kinase/ubiquitination/proteasomal machinery and that modification of their targets on AAV2 capsid may improve its transduction efficiency. In vitro analysis with pharmacological inhibitors of cellular serine/threonine kinases (protein kinase A, protein kinase C, casein kinase II) showed an increase (20-90%) on AAV2-mediated gene expression. The three-dimensional structure of AAV2 capsid was then analyzed to predict the sites of ubiquitination and phosphorylation. Three phosphodegrons, which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, were identified. Mutation targets comprising eight serine (S) or seven threonine (T) or nine lysine (K) residues were selected in and around phosphodegrons on the basis of their solvent accessibility, overlap with the receptor binding regions, overlap with interaction interfaces of capsid proteins, and their evolutionary conservation across AAV serotypes. AAV2-EGFP vectors with the wild-type (WT) capsid or mutant capsids (15 S/T -> alanine A] or 9 K -> arginine R] single mutant or 2 double K -> R mutants) were then evaluated in vitro. The transduction efficiencies of 11 S/T -> A and 7 K -> R vectors were significantly higher (similar to 63-90%) than the AAV2-WT vectors (similar to 30-40%). Further, hepatic gene transfer of these mutant vectors in vivo resulted in higher vector copy numbers (up to 4.9-fold) and transgene expression (up to 14-fold) than observed from the AAV2-WT vector. One of the mutant vectors, S489A, generated similar to 8-fold fewer antibodies that could be cross-neutralized by AAV2-WT. This study thus demonstrates the feasibility of the use of these novel AAV2 capsid mutant vectors in hepatic gene therapy.

Cultivation of Pleurotus spp. on a combination of anaerobically digested plant material and various agro-residues

Two species of Pleurotus, Pleurotus florida and Pleurotus flabellatus were cultivated on two agro-residues (paddy straw; PS and coir pith; CP) singly as well as in combination with biogas digester residue (BDR, main feed leaf biomass). The biological efficiency, nutritional value, composition and nutrient balance (C, N and P) achieved with these substrates were studied. The most suitable substrate that produced higher yields and biological efficiency was PS mixed with BDR followed by coir pith with BDR. Addition of BDR with agro-residues could increase mushroom yield by 20-30%. The biological efficiency achieved was high for PS + BDR (231.93% for P. florida and 209.92% for P. flabellatus) and for CP + BDR (14831% for P. florida and 188.46% for P. flabellatus). The OC (organic carbon), TKN (nitrogen) and TP (phosphate) removal of the Pleurotus spp. under investigation suggests that PS with BDR is the best substrate for growing mushroom. (C) 2015 Published by Elsevier Inc. on behalf of International Energy Initiative.