Structural and functional consequences of amino acid substitutions in the second conserved loop of Escherichia coli adenylate kinase

Molecular Docking Evaluation of (E)-5-arylidene-2-thioxothiazolidin-4-one Derivatives as Selective Bacterial Adenylate Kinase Inhibitors

Multi-drug resistant microorganism infections with emerging problems that require not only a prevention strategy, but also the development of new inhibitory compounds. Six previously synthesized 5-arylidene-2-thioxothiazolidin-4-one derivatives 1a–f, were screened for inhibitory activity on adenylate kinases of different origins by molecular docking. The compounds 1c and 1d were the most efficient inhibitors of bacterial and some archean adenylate kinases. Hydrogen bond interactions were observed with the residues belonging to the ATP binding site. Moreover human adenylate kinases are poor targets, suggesting that this selectivity offers promising prospectives for refining the structure of our compounds.

Dynamics in Thermotoga neapolitana adenylate kinase: 15N relaxation and hydrogen-deuterium exchange studies of a hyperthermophilic enzyme highly active at 30 degrees C

Backbone conformational dynamics of Thermotoga neapolitana adenylate kinase in the free form (TNAK) and inhibitor-bound form (TNAK*Ap5A) were investigated at 30 degrees C using (15)N NMR relaxation measurements and NMR monitored hydrogen-deuterium exchange. With kinetic parameters identical to those of Escherichia coli AK (ECAK) at 30 degrees C, TNAK is a unique hyperthermophilic enzyme. These catalytic properties make TNAK an interesting and novel model to study the interplay between protein rigidity, stability, and activity. Comparison of fast time scale dynamics (picosecond to nanosecond) in the open and closed states of TNAK and ECAK at 30 degrees C reveals a uniformly higher rigidity across all domains of TNAK. Within this framework of a rigid TNAK structure, several residues located in the AMP-binding domain and in the core-lid hinge regions display high picosecond to nanosecond time scale flexibility. Together with the recent comparison of ECAK dynamics with those of hyperthermophilic Aquifex aeolicus AK (AAAK), our results provide strong evidence for the role of picosecond to nanosecond time scale fluctuations in both stability and activity. In the slow time scales, TNAK's increased rigidity is not uniform but localized in the AMP-binding and lid domains. The core domain amides of ECAK and TNAK in the open and closed states show comparable protection against exchange. Significantly, the hinges framing the lid domain show similar exchange data in ECAK and TNAK open and closed forms. Our NMR relaxation and hydrogen-deuterium exchange studies therefore suggest that TNAK maintains high activity at 30 degrees C by localizing flexibility to the hinge regions that are key to facilitating conformational changes.

Crystal structure and nucleotide binding of the Thermus thermophilus RNA helicase Hera N-terminal domain

DEAD box RNA helicases use the energy of ATP hydrolysis to unwind double-stranded RNA regions or to disrupt RNA/protein complexes. A minimal RNA helicase comprises nine conserved motifs distributed over two RecA-like domains. The N-terminal domain contains all motifs involved in nucleotide binding, namely the Q-motif, the DEAD box, and the P-loop, as well as the SAT motif, which has been implicated in the coordination of ATP hydrolysis and RNA unwinding. We present here the crystal structure of the N-terminal domain of the Thermus thermophilus RNA helicase Hera in complex with adenosine monophosphate (AMP). Upon binding of AMP the P-loop adopts a partially collapsed or half-open conformation that is still connected to the DEAD box motif, and the DEAD box in turn is linked to the SAT motif via hydrogen bonds. This network of interactions communicates changes in the P-loop conformation to distant parts of the helicase. The affinity of AMP is comparable to that of ADP and ATP, substantiating that the binding energy from additional phosphate moieties is directly converted into conformational changes of the entire helicase. Importantly, the N-terminal Hera domain forms a dimer in the crystal similar to that seen in another thermophilic prokaryote. It is possible that this mode of dimerization represents the prototypic architecture in RNA helicases of thermophilic origin.

The crystal structure of Mycobacterium tuberculosis adenylate kinase in complex with two molecules of ADP and Mg2+ supports an associative mechanism for phosphoryl transfer

The crystal structure of Mycobacterium tuberculosis adenylate kinase (MtAK) in complex with two ADP molecules and Mg2+ has been determined at 1.9 A resolution. Comparison with the solution structure of the enzyme, obtained in the absence of substrates, shows significant conformational changes of the LID and NMP-binding domains upon substrate binding. The ternary complex represents the state of the enzyme at the start of the backward reaction (ATP synthesis). The structure is consistent with a direct nucleophilic attack of a terminal oxygen from the acceptor ADP molecule on the beta-phosphate from the donor substrate, and both the geometry and the distribution of positive charge in the active site support the hypothesis of an associative mechanism for phosphoryl transfer.

Associative mechanism for phosphoryl transfer: a molecular dynamics simulation of Escherichia coli adenylate kinase complexed with its substrates

The ternary complex of Escherichia coli adenylate kinase (ECAK) with its substrates adenosine monophosphate (AMP) and Mg-ATP, which catalyzes the reversible transfer of a phosphoryl group between adenosine triphosphate (ATP) and AMP, was studied using molecular dynamics. The starting structure for the simulation was assembled from the crystal structures of ECAK complexed with the bisubstrate analog diadenosine pentaphosphate (AP(5)A) and of Bacillus stearothermophilus adenylate kinase complexed with AP(5)A, Mg(2+), and 4 coordinated water molecules, and by deleting 1 phosphate group from AP(5)A. The interactions of ECAK residues with the various moieties of ATP and AMP were compared to those inferred from NMR, X-ray crystallography, site-directed mutagenesis, and enzyme kinetic studies. The simulation supports the hypothesis that hydrogen bonds between AMP's adenine and the protein are at the origin of the high nucleoside monophosphate (NMP) specificity of AK. The ATP adenine and ribose moieties are only loosely bound to the protein, while the ATP phosphates are strongly bound to surrounding residues. The coordination sphere of Mg(2+), consisting of 4 waters and oxygens of the ATP beta- and gamma-phosphates, stays approximately octahedral during the simulation. The important role of the conserved Lys13 in the P loop in stabilizing the active site by bridging the ATP and AMP phosphates is evident. The influence of Mg(2+), of its coordination waters, and of surrounding charged residues in maintaining the geometry and distances of the AMP alpha-phosphate and ATP beta- and gamma-phosphates is sufficient to support an associative reaction mechanism for phosphoryl transfer.

Effects of crowding on the thermal stability of heterogeneous protein solutions

Crowding can substantially affect the transition of a protein between its native (N) and unfolded (U) states via volume exclusion effects. Also, it influences considerably the aggregation (A) of unfolded proteins. To examine the details, we developed an approach for computing the kinetic rates of the process N <--> U --> A in which the concentration of the protein is explicitly taken into account. We then compute the relative change with temperature of the protein denaturation for various fractional volume occupancies and partition of proteins in solution. The analysis indicates that, in protein solutions in which the average distance between proteins is comparable with the radius of gyration of an unfolded protein, steric effects increase the stability of the proteins which are in compact, native states. In heterogeneous protein solutions containing various types of proteins with different thermal stabilities, the unfolding of the most thermolabile proteins will increase the stability of the other proteins. The results shed light on the way proteins change the thermal stability of a cell as they unfold and aggregate. This study may be valuable in questions related to the dynamics of thermal injuries.

Large-scale prediction of protein geometry and stability changes for arbitrary single point mutations

We have developed a method to both predict the geometry and the relative stability of point mutants that may be used for arbitrary mutations. The geometry optimization procedure was first tested on a new benchmark of 2141 ordered pairs of X-ray crystal structures of proteins that differ by a single point mutation, the largest data set to date. An empirical energy function, which includes terms representing the energy contributions of the folded and denatured proteins and uses the predicted mutant side chain conformation, was fit to a training set consisting of half of a diverse set of 1816 experimental stability values for single point mutations in 81 different proteins. The data included a substantial number of small to large residue mutations not considered by previous prediction studies. After removing 22 (approximately 2%) outliers, the stability calculation gave a standard deviation of 1.08 kcal/mol with a correlation coefficient of 0.82. The prediction method was then tested on the remaining half of the experimental data, giving a standard deviation of 1.10 kcal/mol and covariance of 0.66 for 97% of the test set. A regression fit of the energy function to a subset of 137 mutants, for which both native and mutant structures were available, gave a prediction error comparable to that for the complete training set with predicted side chain conformations. We found that about half of the variation is due to conformation-independent residue contributions. Finally, a fit to the experimental stability data using these residue parameters exclusively suggests guidelines for improving protein stability in the absence of detailed structure information.

Determination of intramolecular distance distribution during protein folding on the millisecond timescale

A method for determination of transient (on the millisecond timescale) intramolecular distance distributions (IDDs) by time-resolved dynamic non-radiative excitation energy transfer measurements was developed. The time-course of the development of the IDD between residues 73 and 203 in the CORE domain of Escherichia coli adenylate kinase throughout refolding from the GuHCl-induced denatured state was determined. The mean of the apparent IDD reduced to a value close to its magnitude in the native protein, within 2 ms (the dead-time of the instrument). At that time the width of that distribution was rather large (16+/-2 A). The large width implies that the intramolecular diffusion coefficient of the labeled segment does not exceed 10(-7) cm(2)/second. In a second slower phase of the refolding transition, the width was reduced to its native value (6+/-4 A).

Isolation and characterization of mutated mitochondrial GTP:AMP phosphotransferase

GTP:AMP phosphotransferase (adenylate kinase isozyme 3, AK3) mutants were obtained by using the ability of AK3 to complement a temperature-sensitive mutation of Escherichia coli adenylate kinase (AKe). Five mutants, P16L, G19S, G91D, G91S, and P93L, had mutation sites located at two loops that are involved in substrate binding of the enzyme. P16L and G19S bearing changes at the first loop showed reduced affinity for both GTP and AMP, the extent of reduction being slightly higher for GTP than AMP. In contrast, G91S and P93L having alterations at the second loop had lower affinities for AMP. Only the alterations at the second loop strongly influenced the Vmax value of the enzyme. Another mutant, D163N, had a substitution at the site forming a salt bridge in adenylate kinase isozyme 1 (AK1), which influenced the Vmax as well as the Km values for both substrates. The kinetic characteristics of these mutants were comparable to those of the corresponding AK1 or AKe mutants. Furthermore, from the results of mutations T201P and T201A that had alterations in all the kinetic parameters of AK3 and from a comparison with the structure and the kinetic parameters of AKe, we expect that a residue(s) around Thr201 is involved in recognition of the base of nucleoside triphosphate.

Characterization of metal and nucleotide liganded forms of adenylate kinase by electrospray ionization mass spectrometry

Complexes of adenylate kinase from Escherichia coli, Bacillus subtilis, and Bacillus stearothermophilus with the bisubstrate nucleotide analog P1,P5-di(adenosine 5')-pentaphosphate and with metal ions (Zn2+ and/or Mg2+) were analyzed by electrospray ionization mass spectrometry. P1,P5-di(adenosine 5')-pentaphosphate. adenylate kinase complex was detected in the positive mode at pH as low as 3.8. Binding of nucleotide to adenylate kinase stabilizes the overall structure of the protein and preserves the Zn2+ chelated form of the enzyme from the gram-positive organisms. In this way, it is possible in a single mass spectrometry experiment to screen metal-chelating adenylate kinases, without use of radioactively labeled compounds. Binding of Mg2+ to enzyme via P1,P5-di(adenosine 5')-pentaphosphate was also demonstrated by mass spectrometry. Although no amino acid side chain in adenylate kinase is supposed to interact with Mg2+, Asp93 in porcine muscle cytosolic enzyme, equivalent to Asp84 in the E. coli adenylate kinase, was proposed to stabilize the nucleotide.Mg2+ complex via water molecules.

Mutation at a single acidic amino acid enhances the halophilic behaviour of malate dehydrogenase from Haloarcula marismortui in physiological salts

In a statistical analysis of the amino acid compositions of 26 halophilic proteins, 24 showed an increase in acidic amino acids and a decrease in basic ones when compared to their non-halophilic homologues. The role of acidic residues in halophilic adaptation was investigated by site-directed mutagenesis of malate dehydrogenase (MalDH) from Haloarcula marismortui. In all of 40 non-halophilic homologous proteins, the position aligned with E243 in halophilic MalDH is occupied by a non-acidic amino acid, most frequently by arginine. The E243R mutant of halophilic MalDH was constructed, over-expressed in Escherichia coli, renatured and purified. Its salt-dependent catalytic activity was not affected compared to the wild-type enzyme and both proteins have the same Km values for their substrates. The resistance to denaturation of the mutant was compared to that of the wild-type protein in different physiological salt (NaCl or KCl) and temperature conditions and interpreted in terms of classical quasi-thermodynamic parameters. The mutant is more halophilic than the wild-type protein; it is more sensitive to temperature and requires significantly higher concentrations of NaCl or KCl for equivalent stability. These results highlight the role of acidic amino acids in halophilic behaviour and are in agreement with a model in which these amino acids act cooperatively to organise hydrated ion binding to the protein.

Primary structure of the hydrogenosomal adenylate kinase of Trichomonas vaginalis and its phylogenetic relationships

Hydrogenosomal adenylate kinase of the amitochondriate protist, Trichomonas vaginalis, has been purified and the sequence of its 39 amino-terminal residues established. Based on this sequence and a conserved internal region of the enzyme, a probe was obtained by DNA polymerase chain reaction and used to isolate a genomic DNA clone containing the gene of this enzyme. This gene exists probably as a single copy in T. vaginalis and is not interrupted by introns. The open reading frame obtained codes for a large type adenylate kinase with a mature molecular mass of 24.5 kDa. The T. vaginalis enzyme is homologous with adenylate kinases of other eukaryotes and eubacteria. Strongly conserved parts and residues of the molecule are conserved also in this enzyme. Phylogenetic trees obtained with various methods placed the T. vaginalis adenylate kinase close to the point where the different subfamilies of this enzyme branch from each other, indicating that the T. vaginalis enzyme has no close relationship to any of these subfamilies and that it separated early from other adenylate kinases. The conceptual translation predicts the existence of an amino-terminal nonapeptide absent from the protein purified from hydrogenosomes, similar to the processed amino-terminal extensions of other hydrogenosomal proteins. These extensions have been considered as putative targeting and import signals.

Structural and physico-chemical characteristics of Bordetella pertussis adenylate kinase, a tryptophan-containing enzyme

The adk gene from the Gram-negative pathogen Bordetella pertussis was cloned by complementing the thermosensitive Escherichia coli adk strain CR341T28. B. pertussis adenylate kinase is a 218-amino-acid protein that has high similarity with adenylate kinase from Escherichia coli and Hemophilus influenzae (57%). A distinct characteristic of enzyme from B. pertussis, not found in other bacterial adenylate kinases, is the presence of a tryptophan residue at position 185. Although distant from the catalytic site, this single tryptophan serves as a convenient probe for monitoring the binding of nucleotide substrates or analogs to the enzyme. Differential scanning calorimetry and equilibrium unfolding experiments in guanidine.HCl indicate similar stabilities for adenylate kinase from B. pertussis and E. coli. An extensive comparison between physico-chemical properties of adenylate kinase from B. pertussis and the enzyme from E. coli showed that the kinetic and structural properties of the two enzymes are very similar. However, infrared spectroscopy has allowed to identify small but significant differences in the secondary structure of the two proteins.

Mechanism of adenylate kinase. What can be learned from a mutant enzyme with minor perturbation in kinetic parameters?

The structural and functional roles of threonine-23 in the chicken muscle adenylate kinase (AK) were investigated by site-directed mutagenesis coupled with proton nuclear magnetic resonance (NMR) and phosphorus stereochemistry. The residue is potentially important because it is conserved among all types of AK and is part of the consensus P-loop sequence, 15GXPGXGKGT23. A mutant enzyme T23A (replacing threonine-23 with alanine) was constructed. Analyses of conformational stability and proton NMR indicate that the side chain of this residue contributes little to the structure of AK, which suggests that the side chain of Thr-23 does not play a structural role. The steady-state kinetic data of the mutant enzyme T23A showed no change in kcat and only 5-7-fold increases in Km and dissociation constants. Such minor changes in kinetic data are insufficient to suggest a functional role of Thr-23. However, two-dimensional NMR analyses of WT.MgAP5A and T23A.MgAP5A complexes indicated that the side chain of Thr-23 is in proximity to the adenine ring of the ATP moiety in the WT.MgAP5A complex in solution. In addition, T23A showed a significant perturbation in the stereospecificity toward the diastereomers of (Rp)- and (Sp)-adenosine 5'-(1-thiotriphosphate) (ATP alpha S), with the Rp/Sp ratio increased from < 0.02 in wild-type to 0.37 in T23A. Detailed 31P NMR analysis indicated that the stereospecificity at the AMP site was not perturbed. These results suggest that the side chain of Thr-23 is involved in catalysis, most likely via a hydrogen bonding interaction Thr-OH...O-P alpha(ATP).(ABSTRACT TRUNCATED AT 250 WORDS)