Magnesium spin state determines the energetics of an adenosinetriphosphate molecule

The molecular dynamics method (density functional theory) DFT:B3LYP (6-3IG** basis set, t = 310 K) was used to study interactions between a molecule of adenosinetriphosphate (ATP) (ATP subsystem) and the [Mg(H₂O)₆]²⁺ magnesium cofactor (Mg subsystem) in an aqueous medium simulated by 78 water molecules in the singlet (S) and triplet (T) states. Potential energy surfaces (PESs) for the S (lowest in energy) and T states (highest in energy) are significantly separated in space. Motion along them directs the Mg complex either to oxygen atoms of the γ-β-phosphate groups (O1–O2) (S state of PES) or to oxygen atoms of the β-α-phosphate groups (O2–O3) (T state of PES). Chelation of the γ-β- and β-α-phosphates leads to formation of a stable low-energy ([Mg(H₂O)₄-(OI-O2)ATP]^2?) complex or a metastable high-energy ([Mg(H₂O)₂-(O2–O3)ATP]^2?) complex, respectively, which differ in number of water molecules surrounding the Mg atom. Intersection of two T PESs is accompanied by formation of an unstable state characterized by redistribution of spins between the Mg and ATP subsystems. This state, being sensitive to interaction with the Mg nuclear spin (²⁵Mg), induces an unpaired electron spin, which initiates the ATP cleavage by the ion-radical mechanism, yielding a reactive ion radical of adenosinemonophosphate (·AMP^?), which was earlier found experimentally by the method of chemically induced dynamic nuclear polarization (CIDNP). Biological aspects of the results obtained are discussed.     

Mg spin determines adenosinetriphosphate energetics

The molecular dynamics method DFT:B3LYP (6-31G** basis set, T = 310 K) was used to study interactions between adenosinetriphosphate (ATP), ATP subsystem, and magnesium cofactor [Mg(H2O)6]2+, Mg subsystem, in water environment modeled with 78 water molecules in singlet (S) and triplet (T) states. The lowest in energy singlet (S) and triplet (T) potential energy surfaces, PESs, are remarkably separated in space and direct the Mg cofactor towards the gamma-beta-phosphate oxygens (O1-O2), S path, or towards the beta-alpha-phosphate oxygens (O2-O3), T path. Chelation of the gamma-beta-phosphates and beta2-alpha-phosphates ends, respectively, in the formation of stable, low-energy, ([Mg(H2O)4-(O1-O2)ATP]2-) and metastable, high-energy, ([Mg(H2O)2-(O2-O3)ATP]2-) chelates, differing in the number of water molecules around the Mg. Intersection between the two T PESs produces an unstable state, a result of spin redistribution between the Mg and ATP subsystems. This state, which is sensitive to a hyperfine interaction with the Mg nuclear spin, 25Mg, reveals an unpaired electron spin and initiates the ATP cleavage along the ion-radical path, yielding a highly reactive adenosinemonophosphate ion-radical, *AMP-, earlier observed in the CIDNP (Chemically Induced Dynamic Nuclear Polarization) experiment (A.A. Tulub, 2006). Biological consequences of the findings are discussed.     

New horizons of adenosinetriphosphate energetics arising from interaction with magnesium cofactor

MD DFT:B3LYP (6-31G** basis set, T = 310 K) method is used to study interactions [singlet (S) and triplet (T) reaction paths] between adenosinetriphosphate, ATP⁴⁻, and [Mg(H₂O)₆]²⁺ in water environment, modeled with 78 water molecules. Computations reveal the appearance of low and high-energy states (stable, quasi-stable, and unstable), assigned to different spin symmetries. At the initial stage of interaction, ATP donates a part of its negative charge to the Mg complex making the Mg slightly charged. As a result, the original octahedral Mg complex loses two (S state) or four (T state) water molecules. Moving along S or T potential energy surfaces (PESs), Mg(H₂O)₄ or Mg(H₂O)₂ display different ways of complexation with ATP. S path favors the formation of a stable chelate with the O1–O2 fragment of ATP triphosphate tail, whereas T path favors producing a single-bonded complex with the O2. The latter, being unstable, undergoes a further conversion into a spin-separated complex, also unstable, and two metastable S complexes, which finally arise in two stable, low-energy and high-energy, chelates. The spin-separated complex experiences rapid decomposition resulting in the production of a highly reactive adenosinemonophosphate ion-radical •AMP, early observed in the CIDNP experiment (Tulub 2006). Biological consequences of the findings are discussed.     

The role of magnesium for geometry and charge in GTP hydrolysis, revealed by quantum mechanics/molecular mechanics simulations

The coordination of the magnesium ion in proteins by triphosphates plays an important role in catalytic hydrolysis of GTP or ATP, either in signal transduction or energy conversion. For example, in Ras the magnesium ion contributes to the catalysis of GTP hydrolysis. The cleavage of GTP to GDP and P(i) in Ras switches off cellular signaling. We analyzed GTP hydrolysis in water, Ras, and Ras·Ras-GTPase-activating protein using quantum mechanics/molecular mechanics simulations. By comparison of the theoretical IR-difference spectra for magnesium ion coordinated triphosphate to experimental ones, the simulations are validated. We elucidated thereby how the magnesium ion contributes to catalysis. It provides a temporary storage for the electrons taken from the triphosphate and it returns them after bond cleavage and P(i) release back to the diphosphate. Furthermore, the Ras·Mg(2+) complex forces the triphosphate into a stretched conformation in which the β- and γ-phosphates are coordinated in a bidentate manner. In this conformation, the triphosphate elongates the bond, which has to be cleaved during hydrolysis. Furthermore, the γ-phosphate adopts a more planar structure, driving the conformation of the molecule closer to the hydrolysis transition state. GTPase-activating protein enhances these changes in GTP conformation and charge distribution via the intruding arginine finger.     

The effects of magnesium ions on the interactions of ox brain and liver glutamate dehydrogenase with ATP and GTP.

The effects of ATP and GTP on the activities of ox liver and brain glutamate dehydrogenase were determined in the absence and presence of added Mg2+ ions. Although GTP was an inhibitor of the enzyme reaction assayed in the direction of NAD+ reduction, the magnesium complex of this nucleotide had no effect on the activity. Similarly the magnesium complex of ATP was without effect on the activity of the enzyme although the free nucleotide was an activator. These results suggest that it is important to take account of magnesium complex formation when considering the regulatory actions of these nucleotides.     

Coherent singlet-triplet transitions at the initial step of tubulin assembly into microtubules

Quantum chemistry calculations [DFT-B3LYP QM/MM method, 6-31G** basis set, + ab initio molecular dynamics] were used to study the action of Mg2+ on tubulin properties. It was shown that the hydration of the guanosine triphosphate-tubulin forms a protein zone structure, which includes a electron-occupied zone and a conductivity zone. The binding of Mg2+ to guanosine triphosphate-tubulin results in the unpairing of electrons in the occupied zone (triplet state formation) followed by their transition to the conductivity zone in which the inversion of spin occurs (singlet state formation). The formation of triplet state is the initial step in the subsequent protein dynamics in the picosecond range of time. The dynamics shows up as a coherent oscillating transition of tubulin between the triplet and singlet states, which is evidence of a simultaneous adjustment between nuclear and electron configurations of the protein (ab initio molecular dynamics calculations). The barrier between the triplet and singlet states does not exceed 0.60 kcal x mol(-1). The barrier overcome is considered as electron tunneling through the Fermi surface, which separates the occupied and conductivity zones. Zone formation occurs in the presence of the shell of biological water surrounding the protein.     

The cytochrome c oxidase-azide-nitric oxide complex as a model for the oxygen-binding site

The complex of cytochrome c oxidase with NO and azide has been studied by EPR at 9.2 and 35 GHz. This complex which shows delta ms = 2 EPR triplet and strong anisotropic signals, due to the interaction of cytochrome a2+3 X NO (S = 1/2) and Cu2+B (S = 1/2), is photodissociable . Its action spectrum is similar to that of cytochrome a2+3 X NO with bands at 430, 560 and 595 nm, but shows an additional band in the near ultraviolet region. The quantum yield of the photodissociation process of cytochrome a2+3 X NO in the metal pair appears to depend on the redox state of CuB. When the photolysed sample was warmed to 77 K, a complex was observed with the EPR parameters of cytochrome a3+3 - N-3 - Cu1 +B (S = 1/2). This process of electron and ligand transfer can be reversed by heating the sample to 220 K. It is suggested that in the triplet species azide is bound to Cu2+B whereas NO is bridged between Cu2+B and the haem iron of the cytochrome a2+3. The complex has a triplet ground state and a singlet excited state with an exchange interaction J = -7.1 cm-1 between both spins. The anisotropy in the EPR spectra is mainly due to a magnetic dipole-dipole interaction between cytochrome a2+3 X NO and Cu2+B. From simulations of the triplet EPR spectra obtained at 9 and 35 GHz, a value for the distance between the nitroxide radical and Cu2+B of 0.33 nm was found. A model of the NO binding in the cytochrome a3-Cu pair shows a distance between the haem iron of cytochrome a3 and CuB of 0.45 nm. It is concluded that the cytochrome a3-CuB pair forms a cage in which the dioxygen molecule is bidentate coordinated to the two metals during the catalytic reaction.     

The single chlorophyll a molecule in the cytochrome b6f complex: unusual optical properties protect the complex against singlet oxygen

The cytochrome b(6)f complex of oxygenic photosynthesis mediates electron transfer between the reaction centers of photosystems I and II and facilitates coupled proton translocation across the membrane. High-resolution x-ray crystallographic structures (Kurisu et al., 2003; Stroebel et al., 2003) of the cytochrome b(6)f complex unambiguously show that a Chl a molecule is an intrinsic component of the cytochrome b(6)f complex. Although the functional role of this Chl a is presently unclear (Kuhlbrandt, 2003), an excited Chl a molecule is known to produce toxic singlet oxygen as the result of energy transfer from the excited triplet state of the Chl a to oxygen molecules. To prevent singlet oxygen formation in light-harvesting complexes, a carotenoid is typically positioned within approximately 4 A of the Chl a molecule, effectively quenching the triplet excited state of the Chl a. However, in the cytochrome b(6)f complex, the beta-carotene is too far (> or =14 Angstroms) from the Chl a for effective quenching of the Chl a triplet excited state. In this study, we propose that in this complex, the protection is at least partly realized through special arrangement of the local protein structure, which shortens the singlet excited state lifetime of the Chl a by a factor of 20-25 and thus significantly reduces the formation of the Chl a triplet state. Based on optical ultrafast absorption difference experiments and structure-based calculations, it is proposed that the Chl a singlet excited state lifetime is shortened due to electron exchange transfer with the nearby tyrosine residue. To our knowledge, this kind of protection mechanism against singlet oxygen has not yet been reported for any other chlorophyll-containing protein complex. It is also reported that the Chl a molecule in the cytochrome b(6)f complex does not change orientation in its excited state.     

Molecular dynamics simulations show that bound Mg2+ contributes to amino acid and aminoacyl adenylate binding specificity in aspartyl-tRNA synthetase through long range electrostatic interactions

Molecular recognition between the aminoacyl-tRNA synthetase enzymes and their cognate amino acid ligands is essential for the faithful translation of the genetic code. In aspartyl-tRNA synthetase (AspRS), the co-substrate ATP binds preferentially with three associated Mg2+ cations in an unusual, bent geometry. The Mg2+ cations play a structural role and are thought to also participate catalytically in the enzyme reaction. Co-binding of the ATP x Mg3(2+) complex was shown recently to increase the Asp/Asn binding free energy difference, indicating that amino acid discrimination is substrate-assisted. Here, we used molecular dynamics free energy simulations and continuum electrostatic calculations to resolve two related questions. First, we showed that if one of the Mg2+ cations is removed, the Asp/Asn binding specificity is strongly reduced. Second, we computed the relative stabilities of the three-cation complex and the 2-cation complexes. We found that the 3-cation complex is overwhelmingly favored at ordinary magnesium concentrations, so that the protein is protected against the 2-cation state. In the homologous LysRS, the 3-cation complex was also strongly favored, but the third cation did not affect Lys binding. In tRNA-bound AspRS, the single remaining Mg2+ cation strongly favored the Asp-adenylate substrate relative to Asn-adenylate. Thus, in addition to their structural and catalytic roles, the Mg2+ cations contribute to specificity in AspRS through long range electrostatic interactions with the Asp side chain in both the pre- and post-adenylation states.     

On the roles of magnesium and spermidine in the isoleucyl-tRNA synthetase reaction. Analysis of the reaction mechanism by total rate equations

The reaction of isoleucyl-tRNA synthetase from Escherichia coli B was analysed by deriving total steady-state rate equations for the ATP/PPi exchange reaction and for the aminoacylation of tRNA, and by fitting these rate equations to series of experimental results. The analysis suggests that (a) a Mg2+ inhibits the aminoacylation of tRNA but not the activation of the amino acid. In the chosen mechanism, this enzyme-bound Mg2+ is required at the activation step. (b) Another Mg2+ is required at ATP, but the MgATP apparently can be replaced by the spermidine.ATP complex. Spermidine.ATP is a weaker substrate. The role of spermidine.ATP is especially suggested by the relative rates of the aminoacylation of tRNA when the spermidine and magnesium concentrations are varied. The aminoacylation measurements still suggest that (c) two (or more) Mg2+ are bound to the tRNA molecule and are required for enzyme activity at the transfer step, and that these Mg2+ can be replaced by spermidines.     

Spin polymerization of DNA/RNA nucleotides

The Car-Parrinello Molecular Dynamics (CPMD) has been used to study the ion-radical (IR) polymerization (triplet (T) and singlet (S/T₀) states) of adenine mononucleotides upon interaction with Mg²⁺(H₂O)₂-ATP⁴⁻. It has been found that the IR polymerization occurs only upon Mg²⁺(H₂O)₂-ATP⁴⁻ excitation into a T state (the Franck-Condon or femtosecond laser excitation); it naturally occurs in the dark with DNA polymerase or another Mg-holoenzyme upon interaction of Mg with two Asp residues. The IR path affects only the HO-C_3′ group of ribose, leaving the HO-C_2′ group inactive. The IR polymerization starts with the homolytic removal of the hydrogen atom from the HO-C_3′ group and its transfer onto the hydroxyl radical ·OH, a product of the ATP cleavage, which yields a water molecule. A further progress of the reaction involves interaction between two ion-radicals ·ATP. The reaction is sensitive to the recombination of ·OH and ·ATP. It is mostly suppressed by the appearance of identically directed electron spins on both radicals (the radical pair in the T₀ state) in the vicinity of the HO-C_3′ group and not suppressed in the vicinity of the HO-C_2′ group (the spins in the radical pair are oppositely directed, the radical pair in the T₀ state), making the latter inert on the IR polymerization, but allowing it to be active in the ionic (hydrolytic) polymerization.     

Theoretical studies of the ATP hydrolysis mechanism of myosin

The ATP hydrolysis mechanism of myosin was studied using quantum chemical (QM) and molecular dynamics calculations. The initial model compound for QM calculations was constructed on the basis of the energy-minimized structure of the myosin(S1dc)-ATP complex, which was determined by molecular mechanics calculations. The result of QM calculations suggested that the ATP hydrolysis mechanism of myosin consists of a single elementary reaction in which a water molecule nucleophilically attacked gamma-phosphorus of ATP. In addition, we performed molecular dynamics simulations of the initial and final states of the ATP hydrolysis reaction, that is, the myosin-ATP and myosin-ADP.Pi complexes. These calculations revealed roles of several amino acid residues (Lys185, Thr186, Ser237, Arg238, and Glu459) in the ATPase pocket. Lys185 maintains the conformation of beta- and gamma-phosphate groups of ATP by forming the hydrogen bonds. Thr186 and Ser237 are coordinated to a Mg(2+) ion, which interacts with the phosphates of ATP and therefore contributes to the stabilization of the ATP structure. Arg238 and Glu459, which consisted of the gate of the ATPase pocket, retain the water molecule acting on the hydrolysis at the appropriate position for initiating the hydrolysis.     

Binding of magnesium to myosin subfragment-1 ATPase

Tyr 180 of chicken breast muscle alkali light chain A1 was nitrated with tetranitromethane. The nitroA1 was incorporated into chicken breast muscle subfragment-1 (S-1) by exchange with the intrinsic alkali light chain. In the presence of adenylylimidodiphosphate (AMPPNP) or ADP, the S-1 containing nitroA1 showed a difference visible absorption spectrum by Mg2+ or Ca2+. The difference spectrum has a trough around 435 nm, indicating a blue shift of the absorption spectrum due to the nitrophenol chromophore of the modified A1. The plot of delta A at 435 nm versus concentration of free Mg2+ fitted a single binding curve, independent of the total concentration of AMPPNP. These results reveal that free Mg2+ binds to the active site of S-1 ATPase, but not as Mg-AMPPNP complex. The dissociation constants of magnesium from S-1 complex were different with the two nucleotides and were 1.25 X 10(-8) M and 1.24 X 10(-7) with AMPPNP and ADP, respectively. The difference spectrum was also obtained in the presence of ATP. The delta epsilon value after adding ATP changed with the ATPase reaction. The steady state rate of S-1 ATPase was measured at various concentrations of free Mg2+. The dissociation constant of magnesium from the steady state complex, EPADP(a), was estimated as 6 X 10(-8) M. These results suggest that the affinity of magnesium at the active site of ATPase changes with the intermediate states of ATPase reaction. The affinity of calcium was lower than that of magnesium.     

A quantum chemical investigation of the electronic structure of thionine

We have examined the electronic and molecular structure of 3,7-diaminophenothiazin-5-ium dye (thionine) in the electronic ground state and in the lowest excited states. The electronic structure was calculated using a combination of density functional theory and multi-reference configuration interaction (DFT/MRCI). Equilibrium geometries were optimized employing (time-dependent) density functional theory (B3LYP functional) combined with the TZVP basis set. Solvent effects were estimated using the COSMO model and micro-hydration with up to five explicit water molecules. Our calculated electronic energies are in good agreement with experimental data. We find the lowest excited singlet and triplet states at the ground state geometry to be of π→π* (S(1), S(2), T(1), T(2)) and n→π* (S(3), T(3)) character. This order changes when the molecular structure in the electronically excited states is relaxed. Geometry relaxation has almost no effect on the energy of the S(1) and T(1) states (~0.02 eV). The relaxation effects on the energies of S(2) and T(2) are moderate (0.14-0.20 eV). The very small emission energy results in a very low fluorescence rate. While we were not able to locate the energetic minimum of the S(3) state, we found a non-planar minimum for the T(3) state with an energy which is very close to the energy of the S(1) minimum in the gas phase (0.04 eV above). When hydration effects are taken into account, the n→π* states S(3) and T(3) are strongly blueshifted (0.33 and 0.46 eV), while the π→π* states are only slightly affected (<0.06 eV).     

The use of nucleotide phosphorothioate diastereomers to define the structure of metal-nucleotide bound to GTP-AMP and ATP-AMP phosphotransferases from beef-heart mitochondria

The diastereomers of adenosine 5'-O-[1-thio]triphosphate (ATP[alpha S]) and adenosine 5'-O-[2-thio]triphosphate (ATP[beta S]) were utilized to seek unambiguous assignment of Mg2+ coordination to ATP when bound to ATP-AMP phosphotransferase from beef heart mitochondria (AK2). Similarly, the diastereomers of guanosine 5'-O-[thio]triphosphate (GTP[alpha S]) and guanosine 5'-O-[2-thio]triphosphate (GTP[beta S]) were utilized to seek unambiguous assignment of Mg2+ coordination to GTP when bound to GTP-AMP phosphotransferase from beef heart mitochondria (AK3). Furthermore the diastereomers of guanosine 5'-O-[1-thio]diphosphate (GDP-[alpha S]) have been used to assign Mg2+ coordination to GDP when bound to AK3. The ratios (V for isomer Sp)/(V for isomer Rp) obtained in the presence of Mg2+ and Cd2+ are compared to those already published for ATP-AMP phosphotransferases from pig muscle (AK1) [Kalbitzer et al. (1983) Eur. J. Biochem. 133, 221-227] and from baker's yeast (AKy) [Tomasselli and Noda (1983) Eur. J. Biochem. 132, 109-115]. In all cases, coordination of Mg2+ to the beta-phosphate via the pro-R oxygen is present, as shown by reversal of specificity for the diastereomers of ATP [beta S] or GTP [beta S] respectively on changing the metal ion. In contrast, there is no reversal of specificity for the diastereomers of ATP [alpha S] or GTP[alpha S], or for GDP[alpha S] in the case of AK3 for the reverse reaction, indicating that there is no interaction of the metal with the alpha-phosphate group. The observed stereospecificity for the alpha-thiophosphate is consistent with the assumption of an interaction of the pro-R oxygen of the alpha-phosphate group with the enzyme.     

Identification by time-resolved EPR of the peridinins directly involved in chlorophyll triplet quenching in the peridinin-chlorophyll a-protein from Amphidinium carterae

The mechanism of triplet-triplet energy transfer in the peridinin-chlorophyll-protein (PCP) from Amphidinium carterae was investigated by time-resolved EPR (TR-EPR). The approach exploits the concept of spin conservation during triplet-triplet energy transfer, which leads to spin polarization conservation in the observed TR-EPR spectra. The acceptor (peridinin) inherits the polarization of the donor (chlorophyll) in a way which depends on the relative geometrical arrangement of the donor-acceptor couple. Starting from the initially populated chlorophyll triplet state and taking the relative positions among Chls and peridinins from the X-ray structure of PCP, we calculated the expected triplet state polarization of any peridinin in the complex. Comparison with the experimental data allowed us to propose a path for triplet quenching in the protein. The peridinin-chlorophyll pair directly involved in the triplet-triplet energy transfer coincides with the one having the shortest center to center distance. A water molecule, which is coordinated to the central Mg atom of the Chl, is also placed in close contact with the peridinin. The results support the concept of localization of the triplet state mainly in one specific peridinin in each of the two pigment subclusters related by a pseudo C2 symmetry.     

Investigation of the substrate structure and metal cofactor requirements of the rat liver mitochondrial ATP synthase/ATPase complex

The F1 moiety of the rat liver mitochondrial ATP synthase/ATPase complex contains as isolated 2 mol Mg2+/mol F1, 1 mol of which is nonexchangeable and the other which is exchangeable (N. Williams, J. Hullihen, and P.L. Pedersen, (1987) Biochemistry 26, 162-169). In addition, the enzyme binds 1 mol ADP/mol F1 and 3 mol AMP.PNP, the latter of which can bind in complex formation with divalent cation and displace the Mg2+ at the exchangeable site. Thus, in terms of ligand binding sites the fully loaded rat liver F1 complex contains 3 mol MgAMP.PNP, 1 mol ADP, and 1 mol Mg2+. In this study we have used several metal ATP complexes or analogs thereof to gain further insight into the ligand binding domains of rat liver F1 and the mechanism by which it catalyzes ATP hydrolysis in soluble and membrane bound form. Studies with LaATP confirmed that MgATP is the most likely substrate for rat liver F1, and provided evidence that the enzyme may contain additional Mg2+ binding sites, undetected in previous studies of F1-ATPases, that are required for catalytic activity. Thus, F1 containing the thermodynamically stable LaATP complex in place of MgATP requires added Mg2+ to induce ATP hydrolysis. As Mg2+ cannot readily displace La2+ under these conditions there appears to be a catalytically important class of Mg2+ binding sites on rat liver F1, distinct from the nonexchangeable Mg2+ site and the sites involved in binding MgATP. Additional studies carried out with exchange inert metal-nucleotide complexes involving rhodium and the Mg2+ and Cd2+ complexes of ATP beta S and ATP alpha S imply that the rate-limiting step in the ATPase reaction pathway occurs subsequent to the P gamma-O-P beta bond cleavage steps, perhaps at the level of Mg(ADP)(Pi) hydrolysis or MgADP release. Evidence is presented that Mg2+ remains coordinated to the leaving group of the reaction, i.e., the beta phosphoryl group. Finally, in contrast to soluble F1, F1 bound to F0 in the inner mitochondrial membrane failed to discriminate between the Mg2+ complexes of the ATP beta S isomers. This indicates that a fundamental difference may exist between the catalytic or kinetic mechanism of F1 and the more physiologically intact F0F1 complex.     

Global conformation of the Escherichia coli replication factor DnaC protein in absence and presence of nucleotide cofactors

Global conformational and oligomeric states of the Escherichia coli replicative factor DnaC protein in the absence and presence of magnesium and nucleotide cofactors, ATP and ADP, and their fluorescent analogues, MANT-ATP and MANT-ADP, have been examined using analytical sedimentation velocity and time-dependent fluorescence anisotropy techniques. In solution, the DnaC protein exists exclusively as a monomer over a large protein concentration range. The value of s(degrees) (20, w)= 2.45 +/- 0.07 S indicates that the protein molecule has an elongated shape. When modeled as a prolate ellipsoid of revolution, the hydrated DnaC protein has an axial ratio of 4.0 +/- 0.6 with long axis a = 112 A and the short axis b = 28 A, respectively. The presence of magnesium or nucleotide cofactors, ATP or ADP, does not affect the global conformation of the protein and its monomeric state. These data indicate that recently found cooperative interactions between the DnaC molecules, in the complex with the DnaB helicase, are induced by the binding to the helicase, i.e., they are not the intrinsic property of the DnaC protein. Fluorescence anisotropy decays of the DnaC-MANT-ATP and DnaC-MANT-ADP complexes indicate that the protein has a rigid global structure on the nanosecond time scale, little affected by the nucleotide cofactors. Nevertheless, the complex with ATP has a more flexible structure, while the complex with ADP is more rigid, with the protein molecule assuming a more elongated shape. Magnesium exerts control only on the complex with the ATP analogue. In the absence of magnesium, the ATP analogue is firmly held in the binding site. In the presence of Mg(2+), this fixed location is released and the analogue is allowed to assume a flexible conformational state. The significance of the results for the functioning of the DnaC protein is discussed.     

QM/MM modeling the Ras-GAP catalyzed hydrolysis of guanosine triphosphate

The mechanism of the hydrolysis reaction of guanosine triphosphate (GTP) by the protein complex Ras-GAP (p21(ras) - p120(GAP)) has been modeled by the quantum mechanical-molecular mechanical (QM/MM) and ab initio quantum calculations. Initial geometry configurations have been prompted by atomic coordinates of a structural analog (PDBID:1WQ1). It is shown that the minimum energy reaction path is consistent with an assumption of two-step chemical transformations. At the first stage, a unified motion of Arg789 of GAP, Gln61, Thr35 of Ras, and the lytic water molecule results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low-barrier transition state TS1. At the second stage, Gln61 abstracts and releases protons within the subsystem including Gln61, the lytic water molecule and the gamma-phosphate group of GTP through the corresponding transition state TS2. Direct quantum calculations show that, in this particular environment, the reaction GTP + H(2)O --> GDP + H(2)PO(4) (-) can proceed with reasonable activation barriers of less than 15 kcal/mol at every stage. This conclusion leads to a better understanding of the anticatalytic effect of cancer-causing mutations of Ras, which has been debated in recent years.     

Inhibition of (Na+,K+)-ATPase by magnesium ions and inorganic phosphate and release of these ligands in the cycles of ATP hydrolysis

The kinetic data of magnesium and inorganic phosphate inhibition of the (Na+,K+)-dependent ATP hydrolysis are consistent with a model where both ligands act independently and their release in the ATPase cycle is an ordered process where inorganic phosphate is released first. The effects of magnesium on the stimulation of the ATPase activity by Na+, K+ and ATP, and the inhibition of that activity by inorganic phosphate, are consistent with Mg2+ acting both as a 'product' and as a dead-end inhibitor. The dead-end Mg-enzyme complex would be produced with an enzyme form located downstream in the reaction sequence from the point where Mg2+ acts as a 'product' inhibitor. In the absence of K+, Mg2+ inhibition was reduced when either Na+ or ATP concentrations were increased well beyond those concentrations needed to saturate their high-affinity sites. This ATP effect suggests that the dead-end Mg-enzyme complex formation is affected by the speed of the E2-E1 conformational change. The present model is consistent with the formation of an Mg-phosphoenzyme complex insensitive to K+ which could become K+-sensitive in the presence of high Na+ concentrations. These Mg-enzyme complexes appear as intermediaries in the Na+-ATPase activity found in the absence of external Na+ and K+. These results can be interpreted on the basis of Mg2+ binding to a single site in the enzyme molecule. In addition, these experiments provide kinetic evidence indicating that the stimulation by external Na+ of the ATPase activity in the absence of K+ is due to a K+-like action of Na+ on the external K+ sites.