Autocatalytic cooperativity and self-regulation of ATPase pumps in membrane active transport

We investigate the effect of autocatalysis on the conformational changes of membrane pumps during active transport driven by ATP. The translocation process is described by means of an alternating access model. The usual kinetic scheme is extended by introducing autocatalytic steps and allowing for dynamic formation of enzyme complexes. The usual features of cooperative models are recovered, i.e., sigmoid shapes of flux versus concentration curves. We show also that two autocatalytic steps lead to a mechanism of inhibition by the substrate as experimentally observed for some ATPase pumps. In addition, when the formation of enzyme complexes is allowed, the model exhibits a multiple stationary states regime, which can be related to a self-regulation mechanism of the active transport in biological systems. Correspondence to: G. Weissmüller     

Evidence for two distinct conformations of the Escherichia coli mannitol permease that are important for its transport and phosphorylation functions

Column chromatography of the Escherichia coli mannitol permease (mannitol-specific enzyme II of the phosphotransferase system) in the presence of deoxycholate has revealed that the active permease can exist in at least two association states with apparent molecular weights consistent with a monomer and a dimer. The monomeric conformation is favored by the presence of mannitol and by the phosphoenolpyruvate (PEP)-dependent phosphorylation of the protein. The dimer is stabilized by inorganic phosphate (Pi), which also stimulates phospho-exchange between mannitol and mannitol 1-phosphate (a partial reaction in the overall PEP-dependent phosphorylation of mannitol). Kinetic analysis of the phospho-exchange reaction revealed that Pi stimulates phospho-exchange by increasing the Vmax of the reaction. A kinetic model for mannitol permease function is presented involving both conformations of the permease. The monomer (or a less-stable conformation of the dimer) is hypothesized to be involved in the initial mannitol-binding and PEP-dependent phosphorylation steps, while the stably associated dimer is suggested to participate in later steps involving direct phosphotransfer between the permease, mannitol and mannitol 1-phosphate.     

The dimer state of GyrB is an active form: implications for the initial complex assembly and processive strand passage

In a previous study, we presented the dimer structure of DNA gyrase B' domain (GyrB C-terminal domain) from Mycobacterium tuberculosis and proposed a 'sluice-like' model for T-segment transport. However, the role of the dimer structure is still not well understood. Cross-linking and analytical ultracentrifugation experiments showed that the dimer structure exists both in the B' protein and in the full-length GyrB in solution. The cross-linked dimer of GyrB bound GyrA very weakly, but bound dsDNA with a much higher affinity than that of the monomer state. Using cross-linking and far-western analyses, the dimer state of GyrB was found to be involved in the ternary GyrA-GyrB-DNA complex. The results of mutational studies reveal that the dimer structure represents a state before DNA cleavage. Additionally, these results suggest that the dimer might also be present between the cleavage and reunion steps during processive transport.     

Thermodynamics of the ATPase cycle of GlcV, the nucleotide-binding domain of the glucose ABC transporter of sulfolobus solfataricus

ATP-binding cassette transporters drive the transport of substrates across the membrane by the hydrolysis of ATP. They typically have a conserved domain structure with two membrane-spanning domains that form the transport channel and two cytosolic nucleotide-binding domains (NBDs) that energize the transport reaction. Binding of ATP to the NBD monomer results in formation of a NBD dimer. Hydrolysis of the ATP drives the dissociation of the dimer. The thermodynamics of distinct steps in the ATPase cycle of GlcV, the NBD of the glucose ABC transporter of the extreme thermoacidophile Sulfolobus solfataricus, were studied by isothermal titration calorimetry using the wild-type protein and two mutants, which are arrested at different steps in the ATP hydrolytic cycle. The G144A mutant is unable to dimerize, while the E166A mutant is defective in dimer dissociation. The ATP, ADP, and AMP-PNP binding affinities, stoichiometries, and enthalpies of binding were determined at different temperatures. From these data, the thermodynamic parameters of nucleotide binding, NBD dimerization, and ATP hydrolysis were calculated. The data demonstrate that the ATP hydrolysis cycle of isolated NBDs consists of consecutive steps where only the final step of ADP release is energetically unfavorable.     

The role of competition in creating multiple steady states and favouring evolution of prebiotic polymers

In a prebiotic context the assumption is made that a mutation leads to the formation of a new polymer displaying autocatalytic and cross-catalytic properties. By the use of a simple chemical model it is shown that, despite these favourable attributes, the new species either remains at low concentration or accumulates, unless a competition is imposed between the new polymer and its precursor. In the latter case, the well-known situation of multiple steady states and hysteresis is established and the new “improved” species, under specific conditions, may become the dominant one. The need of competition suggests an early compartmentation of the system.     

Functional organization of single and paired V(D)J cleavage complexes

RAG-1 and RAG-2 initiate V(D)J recombination by binding to specific recognition sequences (RSS) and then cleave the DNA in two steps: nicking and hairpin formation. Recent work has established that a dimer of RAG-1 and either one or two monomers of RAG-2 bind to a single RSS, but the enzymatic contributions of the RAG molecules within this nucleoprotein complex and its functional organization have not been elucidated. Using heterodimeric protein preparations containing both wild-type and catalytically deficient RAG-1 molecules, we found that one active monomer is sufficient for both nicking and hairpin formation at a single RSS, demonstrating that a single active site can carry out both cleavage steps. Furthermore, the mutant heterodimers efficiently cleaved both RSS in a synaptic complex. These results strongly suggest that two RAG-1 dimers are responsible for RSS cleavage in a synaptic complex, with one monomer of each dimer catalyzing both nicking and hairpin formation at each RSS.     

ATP-driven MalK dimer closure and reopening and conformational changes of the "EAA" motifs are crucial for function of the maltose ATP-binding cassette transporter (MalFGK2)

We have investigated conformational changes of the purified maltose ATP-binding cassette transporter (MalFGK(2)) upon binding of ATP. The transport complex is composed of a heterodimer of the hydrophobic subunits MalF and MalG constituting the translocation pore and of a homodimer of MalK, representing the ATP-hydrolyzing subunit. Substrate is delivered to the transporter in complex with periplasmic maltose-binding protein (MalE). Cross-linking experiments with a variant containing an A85C mutation within the Q-loop of each MalK monomer indicated an ATP-induced shortening of the distance between both monomers. Cross-linking caused a substantial inhibition of MalE-maltose-stimulated ATPase activity. We further demonstrated that a mutation affecting the "catalytic carboxylate" (E159Q) locks the MalK dimer in the closed state, whereas a transporter containing the "ABC signature" mutation Q140K permanently resides in the resting state. Cross-linking experiments with variants containing the A85C mutation combined with cysteine substitutions in the conserved EAA motifs of MalF and MalG, respectively, revealed close proximity of these residues in the resting state. The formation of a MalK-MalG heterodimer remained unchanged upon the addition of ATP, indicating that MalG-EAA moves along with MalK during dimer closure. In contrast, the yield of MalK-MalF dimers was substantially reduced. This might be taken as further evidence for asymmetric functions of both EAA motifs. Cross-linking also caused inhibition of ATPase activity, suggesting that transporter function requires conformational changes of both EAA motifs. Together, our data support ATP-driven MalK dimer closure and reopening as crucial steps in the translocation cycle of the intact maltose transporter and are discussed with respect to a current model.     

A nonuniform stepping mechanism for E. coli UvrD monomer translocation along single-stranded DNA

E. coli UvrD is an SF1 helicase involved in several DNA metabolic processes. Although a UvrD dimer is needed for helicase activity, a monomer can translocate with 3' to 5' directionality along single-stranded DNA, and this ATP-dependent translocation is likely involved in RecA displacement. In order to understand how the monomeric translocase functions, we have combined fluorescence stopped-flow kinetic methods with recently developed analysis methods to determine the kinetic mechanism, including ATP coupling stoichiometry, for UvrD monomer translocation along ssDNA. Our results suggest that the macroscopic rate of UvrD monomer translocation is not limited by each ATPase cycle but rather by a slow step (pause) in each translocation cycle that occurs after four to five rapid 1 nt translocation steps, with each rapid step coupled to hydrolysis of one ATP. These results suggest a nonuniform stepping mechanism that differs from either a Brownian motor or previous structure-based inchworm mechanisms.     

Picosecond laser fluorometry of FAD of D-amino acid oxidase-benzoate complex

Formation of a complex of D-amino acid oxidase (D-amino acid:O2 oxidoreductase (deaminating), EC 1.4.3.3) and benzoate, an enzyme-substrate complex model, was studied by measuring the fluorescence life-time of the coenzyme FAD of the complex by using a mode-locked Nd:YAG laser and a streak camera. The value of lifetime was 60 +/- 10 ps in the monomer of the complex and it was extremely short (much less than 5 ps) in the dimer of the complex. Since the values of fluorescence lifetime of the coenzyme are 130 ps in the monomeric form of free enzyme and 40 ps in the dimeric form of free enzyme, the decrease in the lifetime upon complex formation with benzoate is slight in the monomer (reduced to one-half) whereas marked in the dimer (reduced to less than 1/10). By analyzing the fluorescence decay curve, a dissociation constant of the monomer-dimer equilibrium of the complex was evaluated to be 0.4 +/- 0.3 microM, which is much smaller than that in free enzyme. Fluorescence analysis under steady state excitation revealed that the apparent dissociation constant (K) of FAD from the enzyme was decreased by 1:1000 upon the complex formation. Relative quantum yield of the fluorescence of FAD in the complex to that of free FAD exhibited appreciable dependence on the complex concentration: greater in the monomer and less in the dimer. These results suggest that a molecular interaction between FAD and amino acid residue(s) is strengthened by the complex formation, which contributes to a remarkable conformational change in the protein moiety of the complex.     

Kinetic mechanism for formation of the active,dimeric UvrD helicase-DNA complex

Escherichia coli UvrD protein is a 3' to 5' SF1 helicase required for DNA repair as well as DNA replication of certain plasmids. We have shown previously that UvrD can self-associate to form dimers and tetramers in the absence of DNA, but that a UvrD dimer is required to form an active helicase-DNA complex in vitro. Here we have used pre-steady state, chemical quenched flow methods to examine the kinetic mechanism for formation of the active, dimeric helicase-DNA complex. Experiments were designed to examine the steps leading to formation of the active complex, separate from the subsequent DNA unwinding steps. The results show that the active dimeric complex can form via two pathways. The first, faster path involves direct binding to the DNA substrate of a pre-assembled UvrD dimer (dimer path), whereas the second, slower path proceeds via sequential binding to the DNA substrate of two UvrD monomers (monomer path), which then assemble on the DNA to form the dimeric helicase. The rate-limiting step within the monomer pathway involves dimer assembly on the DNA. These results show that UvrD dimers that pre-assemble in the absence of DNA are intermediates along the pathway to formation of the functional dimeric UvrD helicase.     

Impossibility of multiple steady states in a family of active membrane transport models

To study the dynamical behavior of active membrane transport models, Vieira and Bisch proposed a complex chemical network (model 3) with two cycles. One cycle involves monomers as pump units while the other cycle uses dimers. In their work, the stoichiometric network analysis was used to study the stability of steady states and the bifurcation analysis was done through numerical methods. They concluded that the possibility of multiple steady states in the model 3 could not be discarded. Here, a zero eigenvalue analysis is applied to prove the impossibility of multiple positive steady states in the model 3. (A positive steady state is one for which all species have positive concentrations.) Moreover, the result is generalized to its family networks. Received: 6 April 1998 / Revised version: 16 October 1998 / Accepted: 28 October 1998     

Disruption of the bacteriophage T4 Mre11 dimer interface reveals a two-state mechanism for exonuclease activity

The Mre11-Rad50 (MR) complex is a central player in DNA repair and is implicated in the processing of DNA ends caused by double strand breaks. Recent crystal structures of the MR complex suggest that several conformational rearrangements occur during its ATP hydrolysis cycle. A comparison of the Mre11 dimer interface from these structures suggests that the interface is dynamic in nature and may adopt several different arrangements. To probe the functional significance of the Mre11 dimer interface, we have generated and characterized a dimer disruption Mre11 mutant (L101D-Mre11). Although L101D-Mre11 binds to Rad50 and dsDNA with affinity comparable with the wild-type enzyme, it does not activate the ATP hydrolysis activity of Rad50, suggesting that the allosteric communication between Mre11 and Rad50 has been interrupted. Additionally, the dsDNA exonuclease activity of the L101D-MR complex has been reduced by 10-fold under conditions where processive exonuclease activity is required. However, we unexpectedly found that under steady state conditions, the nuclease activity of the L101D-MR complex is significantly greater than that of the wild-type complex. Based on steady state and single-turnover nuclease assays, we have assigned the rate-determining step of the steady state nuclease reaction to be the productive assembly of the complex at the dsDNA end. Together, our data suggest that the Mre11 dimer interface adopts at least two different states during the exonuclease reaction.     

Intramolecular dynamics of low molecular weight protein tyrosine phosphatase in monomer-dimer equilibrium studied by NMR: a model for changes in dynamics upon target binding

Low molecular weight protein tyrosine phosphatase (LMW-PTP) dimerizes in the phosphate-bound state in solution with a dissociation constant of K(d)=1.5(+/-0.1)mM and an off-rate on the order of 10(4)s(-1). 1H and 15N NMR chemical shifts identify the dimer interface, which is in excellent agreement with that observed in the crystal structure of the dimeric S19A mutant. Two tyrosine residues of each molecule interact with the active site of the other molecule, implying that the dimer may be taken as a model for a complex between LMW-PTP and a target protein. 15N relaxation rates for the monomeric and dimeric states were extrapolated from relaxation data acquired at four different protein concentrations. Relaxation data of satisfactory precision were extracted for the monomer, enabling model-free analyses of backbone fluctuations on pico- to nanosecond time scales. The dimer relaxation data are of lower quality due to extrapolation errors and the possible presence of higher-order oligomers at higher concentrations. A qualitative comparison of order parameters in the monomeric and apparent dimeric states shows that loops forming the dimer interface become rigidified upon dimerization. Qualitative information on monomer-dimer exchange and intramolecular conformational exchange was obtained from the concentration dependence of auto- and cross-correlated relaxation rates. The loop containing the catalytically important Asp129 fluctuates between different conformations in both the monomeric and dimeric (target bound) states. The exchange rate compares rather well with that of the catalyzed reaction step, supporting existing hypotheses that catalysis and enzyme dynamics may be coupled. The side-chain of Trp49, which is important for substrate specificity, exhibits conformational dynamics in the monomer that are largely quenched upon formation of the dimer, suggesting that binding is associated with the selection of a single side-chain conformer.     

Cargo-binding makes a wild-type single-headed myosin-VI move processively

Class VI myosin is an intracellular vesicle and organelle transporter that moves along actin filaments in a direction opposite to most other known myosin classes. The myosin-VI was expected to form a dimer to move processively along actin filaments with a hand-over-hand mechanism like other myosin organelle transporters. Recently, however, wild-type myosin-VI was demonstrated to be monomer and single-headed, casting a doubt on its processivity. By using single molecule techniques, we show that green-fluorescent-protein-tagged single-headed, wild-type myosin-VI does not move processively. However, when coupled to 200-nm polystyrene beads (comparable to intracellular vesicles in size) at a ratio of one head per bead, single-headed myosin-VI moves processively with large (40-nm) steps. The characteristics of this monomer-driven movement were different to that of artificial dimer-driven movement: Compared to the artificial dimer, the monomer-bead complex had a reduced stall force (1 pN compared to 2 pN), an average run length 2.5-fold shorter (91 nm compared to 220 nm) and load-dependent step size. Furthermore, we found that a monomer-bead complex moved more processively in a high viscous solution (40-fold higher than water) similar to cellular environment. Because the diffusion constant of the bead is 60-fold lower than myosin-VI heads alone in water, we propose a model in which the bead acts as a diffusional anchor for the myosin-VI, enhancing its rebinding following detachment and supporting processive movement of the bead-monomer complexes. Although a single-headed myosin-VI was able to move processively with a large cargo, the travel distance was rather short. Multiple molecules may be involved in the cargo transport for a long travel distance in cells.     

Opening of the ADP-bound active site in the ABC transporter ATPase dimer: evidence for a constant contact, alternating sites model for the catalytic cycle

ABC transporters are ubiquitous, ATP-dependent transmembrane pumps. The mechanism by which ATP hydrolysis in the nucleotide-binding domain (NBD) effects conformational changes in the transmembrane domain that lead to allocrite translocation remains largely unknown. A possible aspect of this mechanism was suggested by previous molecular dynamics simulations of the MJ0796 NBD dimer, which revealed a novel, nucleotide-dependent intrasubunit conformational change involving the relative rotation of the helical and catalytic subdomains. Here, we find that in four of five simulations of the ADP/ATP-bound dimer, the relative rotation of the helical and catalytic subdomains in the ADP-bound monomer results in opening of the ADP-bound active site, probably sufficient or close to sufficient to allow nucleotide exchange. We also observe that in all five simulations of the ADP/ATP-bound dimer, the intimate contact of the LSGGQ signature sequence with the ATP gamma-phosphate is weakened by the intrasubunit conformational change within the ADP-bound monomer. We discuss how these results support a constant contact model for the function of the NBD dimer in contrast to switch models, in which the NBDs are proposed to fully disassociate during the catalytic cycle.     

Oscillatory isozymes as the simplest model for coupled biochemical oscillators

We analyze a simple model for two autocatalytic reactions catalyzed by two distinct isozymes transforming, with different kinetic properties, a given substrate into the same product. This two-variable system can be viewed as the simplest model of chemically coupled biochemical oscillators. Phase-plane analysis indicates how the kinetic differences between the two enzymes give rise to complex oscillatory phenomena such as the coexistence of a stable steady state and a stable limit cycle, or the co-existence of two simultaneously stable oscillatory regimes (birhythmicity). The model allows one to verify a previously proposed conjecture for the origin of birhythmicity. In other conditions, the system admits multiple oscillatory domains as a function of a control parameter whose variation gives rise to markedly different types of oscillations. The latter behavior provides an explanation for the occurrence of multiple modes of oscillations in thalamic neurons.     

Microtubules as mechanical force sensors

Microtubules are polymers of tubulin subunits (dimers) arranged on a hexagonal lattice. Each tubulin dimer comprises two monomers, the alpha-tubulin and beta-tubulin, and can be found in two states. In the first state a mobile negative charge is located into the alpha-tubulin monomer and in the second into the beta-tubulin monomer. Each tubulin dimer is modeled as an electrical dipole coupled to its neighbors by electrostatic forces. The location of the mobile charge in each dimer depends on the location of the charges in the dimer's neighborhood. Mechanical forces that act on the microtubule affect the distances between the dimers and alter the electrostatic potential. Changes in this potential affect the mobile negative charge location in each dimer and the charge distribution in the microtubule. The net effect is that mechanical forces affect the charge distribution in microtubules. We propose to exploit this effect and use microtubules as mechanical force sensors. We model each dimer as a two-state quantum system and, following the quantum computation paradigm, we use discrete quantum random walk on the hexagonal microtubule lattice to determine the charge distribution. Different forces applied on the microtubule are modeled as different coin biases leading to different probability distributions of the quantum walker location, which are directly connected to different charge distributions. Simulation results show that there is a strong indication that microtubules can be used as mechanical force sensors and that they can also detect the force directions and magnitudes.     

Calcium-induced strain in the monomer promotes dimerization in neural cadherin

Cadherins are cell adhesion proteins that are important for tissue formation and integrity. Cell-cell adhesion occurs through the formation of the strand-crossover dimer between identical cadherins on the surface of neighboring cells. The strand-crossover dimer forms exclusively between their EC1 domains via swapping of the βA sheet by undocking the conserved tryptophan 2, W2, from its own hydrophobic pocket and docking it into the hydrophobic pocket of its adhesive partner. An interesting aspect of the system is the fact that critical noncovalent interactions in the monomer re-form in the dimer. Thus, as these noncovalent interactions are conserved, what drives the formation of dimer? Moreover, why is dimer formation calcium-dependent? Thus, to probe the structural and energetic effects of calcium on the noncovalent interactions that are necessary for dimer formation, we performed spectroscopic, stability, and assembly studies of wild-type and two mutants, W2A and E89A, of neural (N-) cadherin. We find that while the ionic interaction involving E89 has a minimal effect on the general stability of the closed conformation of the βA sheet, the hydrophobic interaction involving W2 is the source of the calcium requirement for adhesive dimer formation. The binding of calcium creates strain in the W2-hydrophbic pocket interaction through direct connection of E11 at the C-terminus of the βA sheet to calcium. To overcome this unfavorable condition in the monomer, N-cadherin forms a dimer. Taken together, our data provide a thermodynamic basis for the calcium dependence of strand-crossover dimer formation in N-cadherin.