Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of adenosine triphosphate (ATP) and adenosine monophosphate (AMP) to adenosine diphosphate (ADP) by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of ATP binding domain (LID) and AMP binding domain (NMP) domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motion of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis, and multidimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two-dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed root mean square deviation variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder-type theoretical model to explain the observed dynamic coupling between the motion of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations. (C) 2011 American Institute of Physics. doi:10.1063/1.3516588]

Role of Water in the Enzymatic Catalysis: Study of ATP plus AMP -> 2ADP Conversion by Adenylate Kinase

The catalytic conversion ATP + AMP -> 2ADP by the enzyme adenylate kinase (ADK) involves the binding of one ATP. molecule to the LID domain and one AMP molecule to the NMP domain. The latter is followed by a. phosphate transfer and then the release of two ADP molecules. We have computed a novel two-dimensional configurational free energy surface (2DCFES), with one reaction coordinate each for the LID and the NMP domain motions, while considering explicit water interactions. Our computed 2DCFES clearly reveals the existence of a stable half-open half-closed (HOHC) intermediate stale of the enzyme. Cycling of the enzyme through the HOHC state reduces the conformational free energy barrier for. the reaction by about 20 kJ/mol. We find that the stability of the HOHC state (missed in all earlier studies with implicit solvent model) is largely because of the increase of specific interactions of the polar amino acid side chains with water, particularly with the arginine and the histidine residues. Free energy surface of the LID domain is rather rugged, which can conveniently slow down LID's conformational motion, thus facilitating a new substrate capture after the product release in the catalytic cycle.

Energetics and Structural Characterization of the large-scale Functional Motion of Adenylate Kinase

Adenylate Kinase (AK) is a signal transducing protein that regulates cellular energy homeostasis balancing between different conformations. An alteration of its activity can lead to severe pathologies such as heart failure, cancer and neurodegenerative diseases. A comprehensive elucidation of the large-scale conformational motions that rule the functional mechanism of this enzyme is of great value to guide rationally the development of new medications. Here using a metadynamics-based computational protocol we elucidate the thermodynamics and structural properties underlying the AK functional transitions. The free energy estimation of the conformational motions of the enzyme allows characterizing the sequence of events that regulate its action. We reveal the atomistic details of the most relevant enzyme states, identifying residues such as Arg119 and Lys13, which play a key role during the conformational transitions and represent druggable spots to design enzyme inhibitors. Our study offers tools that open new areas of investigation on large-scale motion in proteins.

Structural stability of adenylate kinase from the sulfate-reducing bacteria Desulfovibrio gigas

Biophysical Chemistry 110 (2004) 83–92

Adenylate kinase phosphotransfer communicates cellular energetic signals to ATP-sensitive potassium channels

Transduction of energetic signals into membrane electrical events governs vital cellular functions, ranging from hormone secretion and cytoprotection to appetite control and hair growth. Central to the regulation of such diverse cellular processes are the metabolism sensing ATP-sensitive K+ (KATP) channels. However, the mechanism that communicates metabolic signals and integrates cellular energetics with KATP channel-dependent membrane excitability remains elusive. Here, we identify that the response of KATP channels to metabolic challenge is regulated by adenylate kinase phosphotransfer. Adenylate kinase associates with the KATP channel complex, anchoring cellular phosphotransfer networks and facilitating delivery of mitochondrial signals to the membrane environment. Deletion of the adenylate kinase gene compromised nucleotide exchange at the channel site and impeded communication between mitochondria and KATP channels, rendering cellular metabolic sensing defective. Assigning a signal processing role to adenylate kinase identifies a phosphorelay mechanism essential for efficient coupling of cellular energetics with KATP channels and associated functions.

Adenylate kinase complements nucleoside diphosphate kinase deficiency in nucleotide metabolism.

Nucleoside diphosphate (NDP) kinase is a ubiquitous nonspecific enzyme that evidently is designed to catalyze in vivo ATP-dependent synthesis of ribo- and deoxyribonucleoside triphosphates from the corresponding diphosphates. Because Escherichia coli contains only one copy of ndk, the structural gene for this enzyme, we were surprised to find that ndk disruption yields bacteria that are still viable. These mutant cells contain a protein with a small amount NDP kinase activity. The protein responsible for this activity was purified and identified as adenylate kinase. This enzyme, also called myokinase, catalyzes the reversible ATP-dependent synthesis of ADP from AMP. We found that this enzyme from E. coli as well as from higher eukaryotes has a broad substrate specificity displaying dual enzymatic functions. Among the nucleoside monophosphate kinases tested, only adenylate kinase was found to have NDP kinase activity. To our knowledge, this is the first report of NDP kinase activity associated with adenylate kinase.

Electrophoretic variation in adenylate kinase of Neisseria meningitidis is due to inter- and intraspecies recombination.

In prokaryotic and eukaryotic organisms, the electrophoretic variation in housekeeping enzymes from natural populations is assumed to have arisen by the accumulation of stochastic predominantly neutral mutations. In the naturally transformable bacterium Neisseria meningitidis, we show that variation in the electrophoretic mobility of adenylate kinase is due to inter- and intraspecies recombination rather than mutation. The nucleotide sequences of the adenylate kinase gene (adk) from isolates that express the predominant slow electrophoretic variant were rather uniform, differing in sequence at an average of 1.1% of nucleotide sites. The adk sequences of rare isolates expressing the fast migrating variant were identical to each other but had a striking mosaic structure when compared to the adk genes from strains expressing the predominant variant. Thus the sequence from the fast variants was identical to those of typical slow variants in the first 158 bp of the gene but differed by 8.4% in the rest of the gene (nt 159-636). The fast electrophoretic variant appears to have arisen by the replacement of most of the meningococcal gene with the corresponding region from the adk gene of a closely related Neisseria species. The adk genes expressing the electrophoretic variant with intermediate mobility were perfect, or almost perfect, recombinants between the adk genes expressing the fast and slow variants. Recombination may, therefore, play a major role in the generation of electrophoretically detectable variation in housekeeping enzymes of some bacterial species.

GLUT4, AMP kinase, but not the insulin receptor, are required for hepatoportal glucose sensor-stimulated muscle glucose utilization.

Recent evidence suggests the existence of a hepatoportal vein glucose sensor, whose activation leads to enhanced glucose use in skeletal muscle, heart, and brown adipose tissue. The mechanism leading to this increase in whole body glucose clearance is not known, but previous data suggest that it is insulin independent. Here, we sought to further determine the portal sensor signaling pathway by selectively evaluating its dependence on muscle GLUT4, insulin receptor, and the evolutionarily conserved sensor of metabolic stress, AMP-activated protein kinase (AMPK). We demonstrate that the increase in muscle glucose use was suppressed in mice lacking the expression of GLUT4 in the organ muscle. In contrast, glucose use was stimulated normally in mice with muscle-specific inactivation of the insulin receptor gene, confirming independence from insulin-signaling pathways. Most importantly, the muscle glucose use in response to activation of the hepatoportal vein glucose sensor was completely dependent on the activity of AMPK, because enhanced hexose disposal was prevented by expression of a dominant negative AMPK in muscle. These data demonstrate that the portal sensor induces glucose use and development of hypoglycemia independently of insulin action, but by a mechanism that requires activation of the AMPK and the presence of GLUT4.

Structural and kinetic studies of Schistosoma mansoni adenylate kinases

The human parasite Schistosoma mansoni is totally dependent on the purine salvage pathway in order to supply large quantities of purine precursors for its energy and DNA biosynthetic needs. Adenylate kinase (ADK) is responsible for the conversion of AMP (produced by the adenosine kinase reaction) into ADP, which is subsequently converted into ATP by nucleoside diphosphate kinase (NDPK). ADK and NDPK are the most active enzymes of the pathway, probably reflecting an evolutionary adaptation due to the intense use of the branch in which they participate. However, notwithstanding their importance very little information has been accumulated found regarding these enzymes. In this work two adenylate kinases from S. mansoni were cloned and heterologously expressed in Escherichia coil. The purified products were utilized in activity assays, and displayed kinetic parameters similar to the corresponding human orthologous proteins. The cytosolic S. mansoni ADK was crystallized and its structure solved allowing us to detect a difference in the nucleotide binding site when compared with the human ortholog. (C) 2012 Elsevier B.V. All rights reserved.

Structural mobility and transformations in globular proteins

Although globular proteins are endowed with well defined three-dimensional structures, they exhibit substantial mobility within the framework of the given threedimensional structure. The different types of mobility found in proteins by and large correspond to the different levels of organisational hierarchy in protein architecture. They are of considerable structural and functional significance, and can be broadly classified into(a) thermal and conformational fluctuations, (b) segmental mobility, (c) interdomain mobility and (d) intersubunit mobility. Protein crystallographic studies has provided a wealth of information on all of them. The temperature factors derived from X-ray diffraction studies provide a measure of atomic displacements caused by thermal and conformational fluctuations. The variation of displacement along the polypeptide chain have provided functionally significant information on the flexibility of different regions of the molecule in proteins such as myoglobin, lysozyme and prealbumin. Segmental mobility often involves the movement of a region or a segment of a molecule with respect to the rest, as in the transition between the apo and the holo structures of lactate dehydrogenase. It may also involve rigidification of a disordered region of the molecule as in the activation of the zymogens of serine proteases. Transitions between the apo and the holo structures of alcohol dehydrogenase,and between the free and the sugar bound forms of hexokinase, are good examples of interdomain mobility caused by hinge-bending. The capability of different domains to move semi-independently contributes greatly to the versatility of immunoglobulin molecules. Interdomain mobility in citrate synthase appears to be more complex and its study has led to an alternative description of domain closure. The classical and the most thoroughly studied case of intersubunit mobility is that in haemoglobin. The stereochemical mechanism of the action of this allosteric protein clearly brings out the functional subtilities that could be achieved through intersubunit movements. In addition to ligand binding and activation,environmental changes also often cause structural transformations. The reversible transformation between 2 Zn insulin and 4 Zn insulin is caused by changes in the ionic strength of the medium. Adenylate Kinase provides a good example for functionally significant reversible conformational transitions induced by variation in pH. Available evidences indicate that reversible structural transformations in proteins could also be caused by changes in the aqueous environment, including those in the amount of water surrounding protein molecules.

Cereulide producing B. cereus and amylosin producing B. subtilis and B. mojavensis : characterization of strains and toxigenicities

B. cereus is one of the most frequent occurring bacteria in foods . It produces several heat-labile enterotoxins and one stable non-protein toxin, cereulide (emetic), which may be pre-formed in food. Cereulide is a heat stable peptide whose structure and mechanism of action were in the past decade elucidated. Until this work, the detection of cereulide was done by biological assays. With my mentors, I developed the first quantitative chemical assay for cereulide. The assay is based on liquid chromatography (HPLC) combined with ion trap mass spectrometry and the calibration is done with valinomycin and purified cereulide. To detect and quantitate valinomycin and cereulide, their [NH4+] adducts, m/z 1128.9 and m/z 1171 respectively, were used. This was a breakthrough in the cereulide research and became a very powerful tool of investigation. This tool made it possible to prove for the first time that the toxin produced by B. cereus in heat-treated food caused human illness. Until this thesis work (Paper II), cereulide producing B. cereus strains were believed to represent a homogenous group of clonal strains. The cereulide producing strains investigated in those studies originated mostly from food poisoning incidents. We used strains of many origins and analyzed them using a polyphasic approach. We found that the cereulide producing B. cereus strains are genetically and biologically more diverse than assumed in earlier studies. The strains diverge in the adenylate kinase (adk) gene (two sequence types), in ribopatterns obtained with EcoRI and PvuII (three patterns), tyrosin decomposition, haemolysis and lecithine hydrolysis (two phenotypes). Our study was the first demonstration of diversity within the cereulide producing strains of B. cereus. To manage the risk for cereulide production in food, understanding is needed on factors that may upregulate cereulide production in a given food matrix and the environmental factors affecting it. As a contribution towards this direction, we adjusted the growth environment and measured the cereulide production by strains selected for diversity. The temperature range where cereulide is produced was narrower than that for growth for most of the producer strains. Most cereulide was by most strains produced at room temperature (20 - 23ºC). Exceptions to this were two faecal isolates which produced the same amount of cereulide from 23 ºC up until 39ºC. We also found that at 37º C the choice of growth media for cereulide production differed from that at the room temperature. The food composition and temperature may thus be a key for understanding cereulide production in foods as well as in the gut. We investigated the contents of [K+], [Na+] and amino acids of six growth media. Statistical evaluation indicated a significant positive correlation between the ratio [K+]:[Na+] and the production of cereulide, but only when the concentrations of glycine and [Na+] were constant. Of the amino acids only glycine correlated positively with high cereulide production. Glycine is used worldwide as food additive (E 640), flavor modifier, humectant, acidity regulator, and is permitted in the European Union countries, with no regulatory quantitative limitation, in most types of foods. B. subtilis group members are endospore-forming bacteria ubiquitous in the environment, similar to B. cereus in this respect. Bacillus species other than B. cereus have only sporadically been identified as causative agents of food-borne illnesses. We found (Paper IV) that food-borne isolates of B. subtilis and B. mojavensis produced amylosin. It is possible that amylosin was the agent responsible for the food-borne illness, since no other toxic substance was found in the strains. This is the first report on amylosin production by strains isolated from food. We found that the temperature requirement for amylosin production was higher for the B. subtilis strain F 2564/96, a mesophilic producer, than for B. mojavensis strains eela 2293 and B 31, psychrotolerant producers. We also found that an atmosphere with low oxygen did not prevent the production of amylosin. Ready-to-eat foods packaged in micro-aerophilic atmosphere and/or stored at temperatures above 10 °C, may thus pose a risk when toxigenic strains of B. subtilis or B. mojavensis are present.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Genetic diversity of Chinese native pigs inferred from protein electrophoresis

We examined protein polymorphism of 20 native pig breeds in China and 3 introduced pig breeds. Thirty loci have been investigated, among which six loci were found to be polymorphic. Especially, the polymorphism of malate dehydrogenase (MDH), adenylate kinase (AK), and two new alleles of adenosine deaminase (ADA) had not been reported in domestic pigs and wild pigs. The percentage of polymorphic loci (P), the mean heterozygosity (H), and the mean number of alleles (A) are 0.200, 0.065, and 1.300, respectively. The degree of genetic variability of Chinese pigs as a whole was higher than that of goats, lower than that of cattle and horses, and similar to that of sheep. Using the gene frequencies of the 30 loci, Nei's genetic distance among the 20 native breeds in China and 3 introduced pig breeds was calculated by the formula of Nei. The program NEIGHBOR in PHYLIP 3.5c was chosen to construct an UPGMA tree and a NJ tree. Our results show that, of the total genetic variation found in the native pig breeds in China, 31% (0.31) is ascribable to genetic differences among breeds. About 69% of the total genetic variation is found within breeds. Most breeds are in linkage disequilibrium. The patterns of genetic similarities between the Chinese native pig breeds were not in agreement with the proposed pig type classification.

Kinetics and Statistical Distributions of Single-Molecule Conformational Dynamics

We developed a coarse-grained yet microscopic detailed model to study the statistical fluctuations of single-molecule protein conformational dynamics of adenylate kinase. We explored the underlying conformational energy landscape and found that the system has two basins of attractions, open and closed conformations connected by two separate pathways. The kinetics is found to be nonexponential, consistent with single-molecule conformational dynamics experiments. Furthermore, we found that the statistical distribution of the kinetic times for the conformational transition has a long power law tail, reflecting the exponential density of state of the underlying landscape. We also studied the joint distribution of the two pathways and found memory effects.

Efficacy of Tenofovir 1% Vaginal Gel in Reducing the Risk of HIV-1 and HSV-2 Infection

Human Immunodeficiency Virus (HIV) is a retrovirus that can result in rare opportunistic infections occurring in humans. The onset of these infections is known as Acquired Immune Deficiency Syndrome (AIDS). Sexual transmission is responsible for the majority of infections 1, resulting in transmission of HIV due to infected semen or vaginal and cervical secretions containing infected lymphocytes. HIV microbicides are formulations of chemical or biological agents that can be applied to the vagina or rectum with the intention of reducing the acquisition of HIV. Tenofovir is an NRTI that is phosphorylated by adenylate kinase to tenofovir diphosphate, which in turn competes with deoxyadeosine 5’-triphosphate for incorporation into newly synthesized HIV DNA. Once incorporated, tenofovir diphosphate results in chain termination, thus inhibiting viral replication. Tenofovir has been formulated into a range of vaginal formulations, such as rings, tablets gels and films. It has been shown to safe and effective in numerous animal models, while demonstrating safety and acceptability in numerous human trials. The most encouraging results came from the CAPRISA 004 clinical trial which demonstrated that a 1% Tenofovir vaginal gel reduced HIV infection by approximately 39%.