Studies on active site mutants of P-falciparum adenylosuccinate synthetase: Insights into enzyme catalysis and activation

Adenylosuccinate synthetase catalyzes a reversible reaction utilizing IMP, GTP and aspartate in the presence of Mg2+ to form adenylosuccinate, GDP and inorganic phosphate. Comparison of similarly liganded complexes of Plasmodium falciparum, mouse and Escherichia coil AdSS reveals H-bonding interactions involving nonconserved catalytic loop residues (Asn429, Lys62 and Thr307) that are unique to the parasite enzyme. Site-directed mutagenesis has been used to examine the role of these interactions in catalysis and structural organization of P. falciparum adenylosuccinate synthetase (PfAdSS). Mutation of Asn429 to Val, Lys62 to Leu and Thr307 to Val resulted in an increase in K-m values for IMP, GTP and aspartate, respectively along with a 5 fold drop in the k(cat) value for N429V mutant suggesting the role of these residues in ligand binding and/or catalysis. We have earlier shown that the glycolytic intermediate, fructose 1,6 bisphosphate, which is an inhibitor of mammalian AdSS is an activator of the parasite enzyme. Enzyme kinetics along with molecular docking suggests a mechanism for activation wherein F16BP seems to be binding to the Asp loop and inducing a conformation that facilitates aspartate binding to the enzyme active site. Like in other AdSS, a conserved arginine residue (Arg155) is involved in dimer crosstalk and interacts with IMP in the active site of the symmetry related subunit of PfAdSS. We also report on the iochemical characterization of the arginine mutants (R155L, R155K and R155A) which suggests that unlike in E. coil AdSS, Arg155 in PfAdSS influences both ligand binding and catalysis. (C) 2010 Elsevier B.V. All rights reserved.

Studies on Stone Age Housepits in Fennoscandia (4000 - 2000 cal BC) : Changes in ground plan, site location and degree of sedentism

Housepits have a remarkably short research history as compared to Fennoscandian archaeological research on the Stone Age in general. The current understanding of the numbers and the distribution of Stone Age housepits in the Nordic countries has, for the most part, been shaped by archaeological studies carried out over the last twenty to thirty years. The main subjects of this research are Neolithic housepits, which are archaeological remains of semi-subterranean pithouses. This dissertation consists of five peer-reviewed articles and a synthesis paper. The articles deal with the development of housepits as seen in the data gathered from Finland (the Lake Saimaa area and south-eastern Finland) and Russia (the Karelian Isthmus). This synthesis expands the discussion of the changes observed in the Papers to include Fennoscandian housepit research as a whole. Certain changes in the size, shape, environmental location, and clustering of housepits extended into various cultures and ecological zones in northern Fennoscandia. Previously, the evolution of housepits has been interpreted to have been caused by the adaptation of Neolithic societies to prevailing environmental circumstances or to re-organization following contacts with the agrarian Corded Ware/Battle Axe Cultures spreading to North. This dissertation argues for two waves of change in the pithouse building tradition. Both waves brought with them certain changes in the pithouses themselves and in the practices of locating the dwellings in the environment/landscape. The changes in housepits do not go hand in hand with other changes in material culture, nor are the changes restricted to certain ecological environments. Based on current information, it appears that the changes relate primarily to the spread of new concepts of housing and possibly to new technology, as opposed to representing merely a local response to environmental factors. This development commenced already before the birth of the Corded Ware/Battle Axe Cultures. Therefore, the changes are argued to have resulted from the spreading of new ideas through the same networks that actively distributed commodities, exotic goods, and raw materials over vast areas between the southern Baltic Sea, the north-west Russian forest zone, and Fennoscandia.

Understorey vegetation as an indicator of forest site potential in Southern Finland.

The relationship between site characteristics and understorey vegetation composition was analysed with quantitative methods, especially from the viewpoint of site quality estimation. Theoretical models were applied to an empirical data set collected from the upland forests of southern Finland comprising 104 sites dominated by Scots pine (Pinus sylvestris L.), and 165 sites dominated by Norway spruce (Picea abies (L.) Karsten). Site index H100 was used as an independent measure of site quality. A new model for the estimation of site quality at sites with a known understorey vegetation composition was introduced. It is based on the application of Bayes' theorem to the density function of site quality within the study area combined with the species-specific presence-absence response curves. The resulting posterior probability density function may be used for calculating an estimate for the site variable. Using this method, a jackknife estimate of site index H100 was calculated separately for pine- and spruce-dominated sites. The results indicated that the cross-validation root mean squared error (RMSEcv) of the estimates improved from 2.98 m down to 2.34 m relative to the "null" model (standard deviation of the sample distribution) in pine-dominated forests. In spruce-dominated forests RMSEcv decreased from 3.94 m down to 3.16 m. In order to assess these results, four other estimation methods based on understorey vegetation composition were applied to the same data set. The results showed that none of the methods was clearly superior to the others. In pine-dominated forests, RMSEcv varied between 2.34 and 2.47 m, and the corresponding range for spruce-dominated forests was from 3.13 to 3.57 m.

Analysis of forest soil chemistry and hydrology with a dynamic model ACIDIC.

In this study we analyze how the ion concentrations in forest soil solution are determined by hydrological and biogeochemical processes. A dynamic model ACIDIC was developed, including processes common to dynamic soil acidification models. The model treats up to eight interacting layers and simulates soil hydrology, transpiration, root water and nutrient uptake, cation exchange, dissolution and reactions of Al hydroxides in solution, and the formation of carbonic acid and its dissociation products. It includes also a possibility to a simultaneous use of preferential and matrix flow paths, enabling the throughfall water to enter the deeper soil layers in macropores without first reacting with the upper layers. Three different combinations of routing the throughfall water via macro- and micropores through the soil profile is presented. The large vertical gradient in the observed total charge was simulated succesfully. According to the simulations, gradient is mostly caused by differences in the intensity of water uptake, sulfate adsorption and organic anion retention at the various depths. The temporal variations in Ca and Mg concentrations were simulated fairly well in all soil layers. For H+, Al and K there were much more variation in the observed than in the simulated concentrations. Flow in macropores is a possible explanation for the apparent disequilibrium of the cation exchange for H+ and K, as the solution H+ and K concentrations have great vertical gradients in soil. The amount of exchangeable H+ increased in the O and E horizons and decreased in the Bs1 and Bs2 horizons, the net change in whole soil profile being a decrease. A large part of the decrease of the exchangeable H+ in the illuvial B horizon was caused by sulfate adsorption. The model produces soil water amounts and solution ion concentrations which are comparable to the measured values, and it can be used in both hydrological and chemical studies of soils.

Interactions of methoxyamine with pyridoxal-5'-phosphate-Schiff's base at the active site of sheep liver serine hydroxymethyltransferase

The mechanism of interaction of methoxyamine with sheep liver serine hydroxymethyltransferase (EC 2.1.2.1) (SHMT) was established by measuring changes in enzyme activity, visible absorption spectra, circular dichroism and fluorescence, and by evaluating the rate constant by stopped-flow spectrophotometry. Methoxyamine can be considered as the smallest substituted aminooxy derivative of hydroxylamine. It was a reversible noncompetitive inhibitor (Ki = 25 microM) of SHMT similar to O-amino-D-serine. Like in the interaction of O-amino-D-serine and aminooxyacetic acid, the first step in the reaction was very fast. This was evident by the rapid disappearance of the enzyme-Schiff base absorbance at 425 nm with a rate constant of 1.3 x 10(3) M-1 sec-1 and CD intensity at 430 nm. Concomitantly, there was an increase in absorbance at 388 nm (intermediate I). The next step in the reaction was the unimolecular conversion (1.1 x 10(-3) sec-1) of this intermediate to the final oxime absorbing at 325 nm. The identity of the oxime was established by its characteristic fluorescence emission at 460 nm when excited at 360 nm and by high performance liquid chromatography. These results highlight the specificity in interactions of aminooxy compounds with sheep liver serine hydroxymethyltransferase and that the carboxyl group of the inhibitors enhances the rate of the initial interaction with the enzyme.

Source and target enzyme signature in serine protease inhibitor active site sequences

Amino acid sequences of proteinaceous proteinase inhibitors have been extensively analysed for deriving information regarding the molecular evolution and functional relationship of these proteins. These sequences have been grouped into several well defined families. It was found that the phylogeny constructed with the sequences corresponding to the exposed loop responsible for inhibition has several branches that resemble those obtained from comparisons using the entire sequence. The major branches of the unrooted tree corresponded to the families to which the inhibitors belonged. Further branching is related to the enzyme specificity of the inhibitor. Examination of the active site loop sequences of trypsin inhibitors revealed that there are strong preferences for specific amino acids at different positions of the loop. These preferences are inhibitor class specific. Inhibitors active against more than one enzyme occur within a class and confirm to class specific sequence in their loops. Hence, only a few positions in the loop seem to determine the specificity. The ability to inhibit the same enzyme by inhibitors that belong to different classes appears to be a result of convergent evolution

Communication pathways from the anti-codon region to the aminoacylation site in Methionyl tRNA Synthetase: Molecular dynamics simulations and the structure network analysis

The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by ?70 in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross-correlations between the residues in MetRS have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.

Sensitivity of wetland methane emissions to model assumptions: application and model testing against site observations

Abstract. Methane emissions from natural wetlands and rice paddies constitute a large proportion of atmospheric methane, but the magnitude and year-to-year variation of these methane sources is still unpredictable. Here we describe and evaluate the integration of a methane biogeochemical model (CLM4Me; Riley et al., 2011) into the Community Land Model 4.0 (CLM4CN) in order to better explain spatial and temporal variations in methane emissions. We test new functions for soil pH and redox potential that impact microbial methane production in soils. We also constrain aerenchyma in plants in always-inundated areas in order to better represent wetland vegetation. Satellite inundated fraction is explicitly prescribed in the model because there are large differences between simulated fractional inundation and satellite observations. A rice paddy module is also incorporated into the model, where the fraction of land used for rice production is explicitly prescribed. The model is evaluated at the site level with vegetation cover and water table prescribed from measurements. Explicit site level evaluations of simulated methane emissions are quite different than evaluating the grid cell averaged emissions against available measurements. Using a baseline set of parameter values, our model-estimated average global wetland emissions for the period 1993–2004 were 256 Tg CH4 yr−1, and rice paddy emissions in the year 2000 were 42 Tg CH4 yr−1. Tropical wetlands contributed 201 Tg CH4 yr−1, or 78 % of the global wetland flux. Northern latitude (>50 N) systems contributed 12 Tg CH4 yr−1. We expect this latter number may be an underestimate due to the low high-latitude inundated area captured by satellites and unrealistically low high-latitude productivity and soil carbon predicted by CLM4. Sensitivity analysis showed a large range (150–346 Tg CH4 yr−1) in predicted global methane emissions. The large range was sensitive to: (1) the amount of methane transported through aerenchyma, (2) soil pH (± 100 Tg CH4 yr−1), and (3) redox inhibition (± 45 Tg CH4 yr−1).

Crucial contribution of the multiple copies of the initiator tRNA genes in the fidelity of tRNA(fMet) selection on the ribosomal P-site in Escherichia coli

The accuracy of the initiator tRNA (tRNA(fMet)) selection in the ribosomal P-site is central to the fidelity of protein synthesis. A highly conserved occurrence of three consecutive G-C base pairs in the anticodon stem of tRNA(fMet) contributes to its preferential selection in the P-site. In a genetic screen, using a plasmid borne copy of an inactive tRNA(fMet) mutant wherein the three G-C base pairs were changed, we isolated Escherichia coli strains that allow efficient initiation with the tRNA(fMet) mutant. Here, extensive characterization of two such strains revealed novel mutations in the metZWV promoter severely compromising tRNA(fMet) levels. Low cellular abundance of the chromosomally encoded tRNA(fMet) allows efficient initiation with the tRNA(fMet) mutant and an elongator tRNA(Gln), revealing that a high abundance of the cellular tRNA(fMet) is crucial for the fidelity of initiator tRNA selection on the ribosomal P-site in E. coli. We discuss possible implications of the changes in the cellular tRNA(fMet) abundance in proteome remodeling.

Identification of Critical Residues of the Mycobacterial Dephosphocoenzyme A Kinase by Site-Directed Mutagenesis

Dephosphocoenzyme A kinase performs the transfer of the c-phosphate of ATP to dephosphocoenzyme A, catalyzing the last step of coenzyme A biosynthesis. This enzyme belongs to the P-loop-containing NTP hydrolase superfamily, all members of which posses a three domain topology consisting of a CoA domain that binds the acceptor substrate, the nucleotide binding domain and the lid domain. Differences in the enzymatic organization and regulation between the human and mycobacterial counterparts, have pointed out the tubercular CoaE as a high confidence drug target (HAMAP database). Unfortunately the absence of a three-dimensional crystal structure of the enzyme, either alone or complexed with either of its substrates/regulators, leaves both the reaction mechanism unidentified and the chief players involved in substrate binding, stabilization and catalysis unknown. Based on homology modeling and sequence analysis, we chose residues in the three functional domains of the enzyme to assess their contributions to ligand binding and catalysis using site-directed mutagenesis. Systematically mutating the residues from the P-loop and the nucleotide-binding site identified Lys14 and Arg140 in ATP binding and the stabilization of the phosphoryl intermediate during the phosphotransfer reaction. Mutagenesis of Asp32 and Arg140 showed catalytic efficiencies less than 5-10% of the wild type, indicating the pivotal roles played by these residues in catalysis. Non-conservative substitution of the Leu114 residue identifies this leucine as the critical residue from the hydrophobic cleft involved in leading substrate, DCoA binding. We show that the mycobacterial enzyme requires the Mg2+ for its catalytic activity. The binding energetics of the interactions of the mutant enzymes with the substrates were characterized in terms of their enthalpic and entropic contributions by ITC, providing a complete picture of the effects of the mutations on activity. The properties of mutants defective in substrate recognition were consistent with the ordered sequential mechanism of substrate addition for CoaE.

Functional Analysis of an Acid Adaptive DNA Adenine Methyltransferase from Helicobacter pylori 26695

HP0593 DNA-(N-6-adenine)-methyltransferase (HP0593 MTase) is a member of a Type III restriction-modification system in Helicobacter pylori strain 26695. HP0593 MTase has been cloned, overexpressed and purified heterologously in Escherichia coli. The recognition sequence of the purified MTase was determined as 5'-GCAG-3' and the site of methylation was found to be adenine. The activity of HP0593 MTase was found to be optimal at pH 5.5. This is a unique property in context of natural adaptation of H. pylori in its acidic niche. Dot-blot assay using antibodies that react specifically with DNA containing m6A modification confirmed that HP0593 MTase is an adenine-specific MTase. HP0593 MTase occurred as both monomer and dimer in solution as determined by gel-filtration chromatography and chemical-crosslinking studies. The nonlinear dependence of methylation activity on enzyme concentration indicated that more than one molecule of enzyme was required for its activity. Analysis of initial velocity with AdoMet as a substrate showed that two molecules of AdoMet bind to HP0593 MTase, which is the first example in case of Type III MTases. Interestingly, metal ion cofactors such as Co2+, Mn2+, and also Mg2+ stimulated the HP0593 MTase activity. Preincubation and isotope partitioning analyses clearly indicated that HP0593 MTase-DNA complex is catalytically competent, and suggested that DNA binds to the MTase first followed by AdoMet. HP0593 MTase shows a distributive mechanism of methylation on DNA having more than one recognition site. Considering the occurrence of GCAG sequence in the potential promoter regions of physiologically important genes in H. pylori, our results provide impetus for exploring the role of this DNA MTase in the cellular processes of H. pylori.

Binding of active site directed ligands to phospholipase A2: Implications on the molecular constraints and catalytic mechanism

Molecular constraints for the localization of active site directed ligands (competitive inhibitors and substrates) in the active site of phospholipase A2 (PLA2) are characterized. Structure activity relationships with known inhibitors suggest that the head : group interactions dominate the selectivity as well as a substantial part of the affinity. The ab initio fitting of the amide ligands in the active site was carried out to characterize the head group interactions. Based on a systematic coordinate space search, formamide is docked with known experimental constraints such as coordination of the carbonyl group to Ca2+ and hydrogen bond between amide nitrogen and ND1 of His48. An optimal position for a bound water molecule is identified and its significance for the catalytic mechanism is postulated. Unlike the traditional ''pseudo-triad'' mechanism, the ''Ca-coordinatedoxyanion'' mechanism proposed here invokes activation of the catalytic water to form the oxyanion in the coordination sphere of calcium. As it attacks the carbonyl carbon of the ester, a near-tetrahedral intermediate is formed. As the second proton of the catalytic water is abstracted by the ester oxygen, its reorientation and simultaneous cleavage form hydrogen bond with ND1 of His48. In this mechanism of esterolysis, a catalytic role for the water co-ordinated to Ca2+ is recognised.

Interactions of L-serine at the active site of serine hydroxymethyltransferases: induction of thermal stability

Serine hydroxymethyltransferase (SHMT), EC 2.1.2.1, exhibits broad substrate and reaction specificity. In addition to cleaving many 3-hydroxyamino acids to glycine and an aldehyde, the enzyme also catalyzed the decarboxylation, transamination and racemization of several substrate analogues of amino acids. To elucidate the mechanism of interaction of substrates, especially L-serine with the enzyme, a comparative study of interaction of L-serine with the enzyme from sheep liver and Escherichia coli, was carried out. The heat stability of both the enzymes was enhanced in the presence of serine, although to different extents. Thermal denaturation monitored by spectral changes indicated an alteration in the apparent T, of sheep liver and E. coli SHMTs from 55 +/- 1 degrees C to 72 +/- 3 degrees C at 40 mM serine and from 67 +/- 1 degrees C to 72 +/- 1 degrees C at 20 mM serine, respectively. Using stopped flow spectrophotometry k values of (49 +/- 5)(.)10(-3) s(-1) and (69 +/- 7).10(-3) s(-1) for sheep liver and E. coli enzymes were determined at 50 mM serine. The binding of serine monitored by intrinsic fluorescence and sedimentation velocity measurements indicated that there was no generalized change in the structure of both proteins. However, visible CD measurements indicated a change in the asymmetric environment of pyridoxal 5'-phosphate at the active site upon binding of serine to both the enzymes. The formation of an external aldimine was accompanied by a change in the secondary structure of the enzymes monitored by far UV-CD spectra. Titration microcalorimetric studies in the presence of serine (8 mM) also demonstrated a single class of binding and the conformational changes accompanying the binding of serine to the enzyme resulted in a more compact structure leading to increased thermal stability of the enzyme.

Identification of the active-site peptide of 2,3-dihydroxybenzoic acid decarboxylase from Aspergillus oryzae

The non-oxidative decarboxylation of aromatic acids is a poorly understood reaction. The transformation of 2,3-dihydroxybenzoic acid to catechol in the fungal metabolism of indole is a prototype of such a reaction. 2,3-Dihydroxybenzoic acid decarboxylase (EC 4.1.1.46) which catalyzes this reaction was purified to homogeneity from anthranilate induced cultures of Aspergillus oryzae using affinity chromatography. The enzyme did not require cofactors like NAD(+), PLP, TPP or metal ions for its activity. There was no spectral evidence for the presence of enzyme bound cofactors. The preparation, which was adjudged homogeneous by the criteria of SDS-PAGE, sedimentation analysis and N-terminal analysis, was characterized for its physicochemical and kinetic parameters. The enzyme was inactivated by group-specific modifiers like diethyl pyrocarbonate (DEPC) and N-ethylmaleimide (NEM). The kinetics of inactivation by DEPC suggested the presence of a single class of essential histidine residues, the second order rate constant of inactivation for which was 12.5 M(-1) min(-1). A single class of cysteine residues was modified by NEM with a second order rate constant of 33 M(-1) min(-1). Substrate analogues protected the enzyme against inactivation by both DEPC and NEM, suggesting the Location of the essential histidine and cysteine to be at the active site of the enzyme. The incorporation of radiolabelled NEM in a differential labelling experiment was 0.73 mol per mol subunit confirming the presence of a single essential cysteine per active-site. Differentially labelled enzyme was enzymatically cleaved and the peptide bearing the label was purified and sequenced. The active-site peptide LLGLAETCK and the N-terminal sequence MLGKIALEEAFALPRFEEKT did not bear any similarity to sequences reported in the Swiss-Prot Protein Sequence Databank, a reflection probably of the unique primary structure of this novel enzyme. The sequences reported in this study will appear in the Swiss-Prot Protein Sequence Databank under the accession number P80402.

Probing the role of the fully conserved Cys126 in triosephosphate isomerase by site-specific mutagenesis - distal effects on dimer stability

Cys126 is a completely conserved residue in triosephosphate isomerase that is proximal to the active site but has been ascribed no specific role in catalysis. A previous study of the C126S and C126A mutants of yeast TIM reported substantial catalytic activity for the mutant enzymes, leading to the suggestion that this residue is implicated in folding and stability [Gonzalez-Mondragon E et al. (2004) Biochemistry43, 3255–3263]. We re-examined the role of Cys126 with the Plasmodium falciparum enzyme as a model. Five mutants, C126S, C126A, C126V, C126M, and C126T, were characterized. Crystal structures of the 3-phosphoglycolate-bound C126S mutant and the unliganded forms of the C126S and C126A mutants were determined at a resolution of 1.7–2.1 Å. Kinetic studies revealed an approximately five-fold drop in kcat for the C126S and C126A mutants, whereas an approximately 10-fold drop was observed for the other three mutants. At ambient temperature, the wild-type enzyme and all five mutants showed no concentration dependence of activity. At higher temperatures (> 40 °C), the mutants showed a significant concentration dependence, with a dramatic loss in activity below 15 μm. The mutants also had diminished thermal stability at low concentration, as monitored by far-UV CD. These results suggest that Cys126 contributes to the stability of the dimer interface through a network of interactions involving His95, Glu97, and Arg98, which form direct contacts across the dimer interface.