The Molecular Mechanism of Eukaryotic Elongation Factor 2 Kinase Activation

Calmodulin (CaM)-dependent eukaryotic elongation factor 2 kinase (eEF-2K) impedes protein synthesis through phosphorylation of eukaryotic elongation factor 2 (eEF-2). It is subject to complex regulation by multiple upstream signaling pathways, through poorly described mechanisms. Precise integration of these signals is critical for eEF-2K to appropriately regulate protein translation rates. Here, an allosteric mechanism comprising two sequential conformations is described for eEF-2K activation. First, Ca²⁺/CaM binds eEF-2K with high affinity (K_d(CaM)^app = 24 ± 5 nm) to enhance its ability to autophosphorylate Thr-348 in the regulatory loop (R-loop) by > 10⁴-fold (k_auto = 2.6 ± 0.3 s⁻¹). Subsequent binding of phospho-Thr-348 to a conserved basic pocket in the kinase domain potentially drives a conformational transition of the R-loop, which is essential for efficient substrate phosphorylation. Ca²⁺/CaM binding activates autophosphorylated eEF-2K by allosterically enhancing k_cat^app for peptide substrate phosphorylation by 10³-fold. Thr-348 autophosphorylation results in a 25-fold increase in the specificity constant (k_cat^app/K_m(Pep-S)^app), with equal contributions from k_cat^app and K_m(Pep-S)^app, suggesting that peptide substrate binding is partly impeded in the unphosphorylated enzyme. In cells, Thr-348 autophosphorylation appears to control the catalytic output of active eEF-2K, contributing more than 5-fold to its ability to promote eEF-2 phosphorylation. Fundamentally, eEF-2K activation appears to be analogous to an amplifier, where output volume may be controlled by either toggling the power switch (switching on the kinase) or altering the volume control (modulating stability of the active R-loop conformation). Because upstream signaling events have the potential to modulate either allosteric step, this mechanism allows for exquisite control of eEF-2K output.     

Stopped-Flow Fluorescence Kinetic Study of Protein Sliding and Intersegment Transfer in the Target DNA Search Process

Kinetic characterizations of protein translocation on DNA are nontrivial because the simultaneous presence of multiple different mechanisms makes it difficult to extract the information specific to a particular translocation mechanism. In this study, we have developed new approaches for the kinetic investigations of proteins' sliding and intersegment transfer (also known as “direct transfer”) in the target DNA search process. Based on the analytical expression of the mean search time for the discrete-state stochastic model, we derived analytical forms of the apparent rate constant k_app for protein-target association in systems involving competitor DNA and the intersegment transfer mechanism. Our analytical forms of k_app facilitate the experimental determination of the kinetic rate constants for intersegment transfer and sliding in the target association process. Using stopped-flow fluorescence data for the target association kinetics along with the analytical forms of k_app, we have studied the translocation of the Egr-1 zinc-finger protein in the target DNA association process. Sliding was analyzed using the DNA-length-dependent k_app data. Using the dependence of k_app on the concentration of competitor DNA, we determined the second-order rate constant for intersegment transfer. Our results indicate that a major pathway in the target association process for the Egr-1 zinc-finger protein is the one involving intersegment transfer to a nonspecific site and the subsequent sliding to the target.     

A23187—Channel behaviour: Fluorescence study

Pyranine entrapped soylipid liposomes have been used as a model system to study the proton transport across membrane in the presence of A23187, a carboxylic ionophore specific for electroneutral exchange of divalent cations. An apparent rate constant (k_app) for transport of protons has been determined from the rate of change of fluorescence intensity of pyranine by stopped flow rapid kinetics in the presence of proton gradient The variation of k_app has been studied as a function of ionophore concentration and the results have been compared with gramicidin—a well known channel former under the similar experimental conditions. The rates thus obtained showed that A23187 is not only a simple carrier but also shows channel behaviour at high concentration of ionophore.     

Kinetics of protein aggregation. Quantitative estimation of the chaperone-like activity in test-systems based on suppression of protein aggregation

The experimental data on the kinetics of irreversible aggregation of proteins caused by exposure to elevated temperatures or the action of denaturing agents (guanidine hydrochloride, urea) have been analyzed. It was shown that the terminal phase of aggregation followed, as a rule, first order kinetics. For the kinetic curves registered by an increase in the apparent absorbance (A) in time (t) the methods of estimation of the corresponding kinetic parameters A _lim and k _I (A _lim is the limiting value of A at t and k _I is the rate constant of the first order) have been proposed. Cases are revealed when the reaction rate constant k _I calculated from the kinetic curve of aggregation of the enzymes coincides with the rate constant for enzyme inactivation. Such a situation is interpreted as a case when the rate of aggregation is limited by the stage of denaturation of the enzyme. A conclusion has been made that, in order to establish the mechanism of protein aggregation, the kinetic investigations of aggregation should be carried out over a wide range of protein concentrations. The refolding experiments after denaturation of proteins by guanidine hydrochloride or urea have been also analyzed. It was shown that aggregation accompanying refolding follows first order kinetics at the final phase of the process. The model of protein refolding explaining such a kinetic regularity has been proposed. When aggregation of protein substrate follows first order kinetics, parameters A _lim and k _I may be used for the quantitative characterization of the chaperone-like activity in the test-systems based on suppression of protein aggregation.     

α-Synuclein aggregation variable temperature and variable pH kinetic data: A re-analysis using the Finke–Watzky 2-step model of nucleation and autocatalytic growth

The aggregation of proteins is believed to be intimately connected to many neurodegenerative disorders. We recently reported an “Ockham's razor”/minimalistic approach to analyze the kinetic data of protein aggregation using the Finke–Watzky (F–W) 2-step model of nucleation (A → B, rate constant k₁) and autocatalytic growth (A + B → 2B, rate constant k₂). With that kinetic model we have analyzed 41 representative protein aggregation data sets in two recent publications, including amyloid β, α-synuclein, polyglutamine, and prion proteins (Morris, A. M., et al. (2008) Biochemistry 47, 2413-2427; Watzky, M. A., et al. (2008) Biochemistry 47, 10790–10800). Herein we use the F–W model to reanalyze protein aggregation kinetic data obtained under the experimental conditions of variable temperature or pH 2.0 to 8.5. We provide the average nucleation (k₁) and growth (k₂) rate constants and correlations with variable temperature or varying pH for the protein α-synuclein. From the variable temperature data, activation parameters ΔG^‡, ΔH^‡, and ΔS^‡ are provided for nucleation and growth, and those values are compared to the available parameters reported in the previous literature determined using an empirical method. Our activation parameters suggest that nucleation and growth are energetically similar for α-synuclein aggregation (ΔG^‡_nucleation = 23(3) kcal/mol; ΔG^‡_growth = 22(1) kcal/mol at 37 °C). From the variable pH data, the F–W analyses show a maximal k₁ value at pH ~ 3, as well as minimal k₁ near the isoelectric point (pI) of α-synuclein. Since solubility and net charge are minimized at the pI, either or both of these factors may be important in determining the kinetics of the nucleation step. On the other hand, the k₂ values increase with decreasing pH (i.e., do not appear to have a minimum or maximum near the pI) which, when combined with the k₁ vs. pH (and pI) data, suggest that solubility and charge are less important factors for growth, and that charge is important in the k₁, nucleation step of α-synuclein. The chemically well-defined nucleation (k₁) rate constants obtained from the F–W analysis are, as expected, different than the 1/lag-time empirical constants previously obtained. However, k₂ × [A]₀ (where k₂ is the rate constant for autocatalytic growth and [A]₀ is the initial protein concentration) is related to the empirical constant, k_app obtained previously. Overall, the average nucleation and average growth rate constants for α-synuclein aggregation as a function of pH and variable temperature have been quantitated. Those values support the previously suggested formation of a partially folded intermediate that promotes aggregation under high temperature or acidic conditions.     

Elucidating slow binding kinetics of a protein without observable bound resonances by longitudinal relaxation NMR spectroscopy

We developed a new method to elucidate the binding kinetics k_on and k_off, and the dissociation constant K_D (=k_off/k_on), of protein-protein interactions without observable bound resonances of the protein of interest due to high molecular weight in a complex with a large target protein. In our method, k_on and k_off rates are calculated from the analysis of longitudinal relaxation rates of free resonances measured for multiple samples containing different concentration ratios of ¹⁵N-labeled protein and substoichiometric amounts of the target protein. The method is applicable to interactions that cannot be analyzed by relaxation dispersion spectroscopy due to slow interactions on millisecond to second timescale and/or minimal conformational (chemical shift) change upon binding. We applied the method to binding of the B1 domain of protein G (GB1) to immunoglobulin G, and derived the binding kinetics despite the absence of observable bound GB1 resonances.     

High affinity DNA-microtubule associated protein interaction

We have isolated the MAP/tau proteins from twice-cycled chick brain microtubule preparations and demonstrated that they are responsible for the nitrocellulose DNA binding activity we and others have measured. Using the isolated MAP/tau proteins we then measured the apparent affinity constant K_app for the homologous chick DNA interaction and found evidence for two equilibrium affinity classes-a K_app = 6 × 10⁷ M^–1, responsible for the bulk of the DNA binding activity and a small (< 10%) higher affinity K_app = 10⁸ – 10⁹ M^–1, likely due to sequence specific binding protein species. Using the same chick brain MAP-tau protein, a heterologous interaction with D. melanogaster DNA, was found to possess just the lower affinity class-K_app = 2 × 10⁷ M^–1. Under stringent binding conditions we carried out equilibrium nitrocellulose filter binding experiments in a ternary reaction mixture at constant MAP/tau protein and ³⁵S radiolabelled chick DNA concentration using increasing and excess concentrations of competitor DNAs of different sources. The order of competitor strengths found was-chick DNA > mouse DNA > D. melanogaster = E. coli. DNA. These data and specifically the homologous DNA: protein case being the strongest competitor corroborate our previous studies using total microtubule protein and provide new evidence for a conserved interaction of a small DNA sequence class with MAP/tau protein species. Moreover, these data allow us to conclude that the conserved DNA sequence: MAP/tau protein interactions do not critically depend upon any energetic feature co-involving tubulin for their properties since tubulin is absent from these preparations.     

Fitting protein-folding free energy landscape for a certain conformation to an NK fitness landscape

The NK fitness landscape is a mathematical landscape model with a parameter k that governs the degree of ruggedness of the landscape. We presented a procedure to fit a given landscape to the NK fitness landscape by introducing the “apparent k-value” k_app. In this paper, we defined the protein free energy (ΔG) landscape in amino acid sequence space, where ΔG is the folding free energy from a random coil to a “certain conformation”. Applying this landscape to our fitting procedure, we examined the statistical properties of the landscape. For calculation of a conformation energy, amino acid residues are represented by points, and interaction energies among amino acid residues are given as (1+K)-body interactions, that is, an unit of interacting (1+K) residues cooperatively contribute a single energy value to the conformational energy. Our results suggest that the apparent k-value of the free energy landscape is k_app≈K, and that the number of possible interactions among residues is unrelated to the k_app value. This leads to the inference that k_app takes values about 1-3 in real landscapes, if nature adopts two-body ∼four-body interaction energies.     

Influences of nonuniform activity distributions on the apparent maximum reaction rate and apparent Michaelis constant of immobilized enzyme reactions

The influences of nonuniform activity distribution within a porous solid support on the apparent kinetic parameters, V_m^app and K_m^app, of immobilized enzyme reactions following the Michaelis-Menten kinetics were theoretically investigated. As the enzyme is distributed to the neighborhood of the external surface of the support, V_m^app and K_m^app approach their respective intrinsic values over a wide range of substrate concentration. There is a close relationship between the nonuniform distribution and internal diffusional resistance. Changes in these two factors provide similar effects on V_m^app and K_m^app. As long as the immobilized enzyme reaction follows Michaelis-Menten kinetics, the nonuniform activity distribution never makes K_m^app less than its intrinsic value.     

An Analysis of Biomolecular Force Fields for Simulations of Polyglutamine in Solution

Polyglutamine (polyQ) peptides are a useful model system for biophysical studies of protein folding and aggregation, both for their intriguing aggregation properties and their own relevance to human disease. The genetic expansion of a polyQ tract triggers the formation of amyloid aggregates associated with nine neurodegenerative diseases. Several clearly identifiable and separable factors, notably the length of the polyQ tract, influence the mechanism of aggregation, its associated kinetics, and the ensemble of structures formed. Atomistic simulations are well positioned to answer open questions regarding the thermodynamics and kinetics of polyQ folding and aggregation. The additional, explicit representation of water permits deeper investigation of the role of solvent dynamics, and it permits a direct comparison of simulation results with infrared spectroscopy experiments. The generation of meaningful simulation results hinges on satisfying two essential criteria: achieving sufficient conformational sampling to draw statistically valid conclusions, and accurately reproducing the intermolecular forces that govern system structure and dynamics. In this work, we examine the ability of 12 biomolecular force fields to reproduce the properties of a simple, 30-residue polyQ peptide (Q₃₀) in explicit water. In addition to secondary and tertiary structure, we consider generic structural properties of polymers that provide additional dimensions for analysis of the highly degenerate disordered states of the molecule. We find that the 12 force fields produce a wide range of predictions. We identify AMBER ff99SB, AMBER ff99SB^∗, and OPLS-AA/L to be most suitable for studies of polyQ folding and aggregation.     

Purification and Biochemical Characterization of Digestive Lipase in Whiteleg Shrimp

Penaeus vannamei lipase was purified from midgut gland of whiteleg shrimp. Pure lipase (E.C. 3.1.1.3) was obtained after Superdex 200 gel filtration and Resource Q anionic exchange. The pure lipase, which is a glycosylated molecule, is a monomer having a molecular mass of about 44.8 kDa, as determined by SDS-PAGE analysis. The lipase hydrolyses short and long-chain triacylglycerols and naphthol derivates at comparable rates. A specific activity of 1787 U mg⁻¹ and 475 U mg⁻¹ was measured with triolein and tributyrin as substrates, respectively, at pH 8.0 and 30°C in the absence of colipase. The lipase showed a K _{m, app} of 3.22 mM and k _{cat, app}/K _{m, app} of 0.303 × 10³ mM⁻¹ s⁻¹ using triolein as substrate. Natural detergents, such as sodium deoxycholate, act as potent inhibitors of the lipase. This inhibition can be reversed by adding fresh oil emulsion. Result with tetrahydrolipstatin, an irreversible inhibitor, suggests that the lipase is a serine enzyme. Peptide sequences of the lipase were determined and compared with the full-length sequence of lipase which was obtained by the rapid amplification of cDNA ends method. The full cDNA of the pvl was 1,186 bp, with a deduced protein of 362 amino acids that includes a consensus sequence (GXSXG) of the lipase superfamily of α/β-hydrolase. The gene exhibits features of conserved catalytic residues and high homology with various mammalian and insect lipase genes. A potential lid sequence is suggested for pvl.     

Study of the thermal stability of D-amino acid oxidase from Trigonopsis variabilis reveals enzyme inactivation via multiple steps

The thermal stability of a highly purified preparation of D-amino acid oxidase from Trigonopsis variabilis (TvDAO), which does not show microheterogeneity due to the partial oxidation of Cys-108, was studied based on dependence of temperature (20–60°C) and protein concentration (5–100 µmol L^?1). The time courses of loss of enzyme activity in 100 mmol L^?1 potassium phosphate buffer, pH 8.0, are well described by a formal kinetic mechanism in which two parallel denaturation processes, partial thermal unfolding and dissociation of the FAD cofactor, combine to yield the overall inactivation rate. Estimates from global fitting of the data revealed that the first-order rate constant of the unfolding reaction (k_a) increased 10⁴-fold in response to an increase in temperature from 20 to 60°C. The rate constants of FAD release (k_b) and binding (k_?b) as well as the irreversible aggregation of the apo-enzyme (k_agg) were less sensitive to changes in temperature, their activation energy (E_a) being about 52 kJ mol^?1 in comparison with an E_a value of 185 kJ mol^?1 for k_a. The rate-determining step of TvDAO inactivation switched from FAD dissociation to unfolding at high temperatures. The model adequately described the effect of protein concentration on inactivation kinetics. Its predictions regarding the extent of FAD release and aggregation during thermal denaturation were confirmed by experiments. TvDAO is shown to contain two highly reactive cysteines per protein subunit whose modification with 5,5′-dithio-bis (2-nitrobenzoic acid) was accompanied by inactivation. Dithiothreitol (1 mmol L^?1) enhanced up to 10-fold the recovery of enzyme activity during ion exchange chromatography of technical-grade TvDAO. However, it did not stabilize TvDAO at all temperatures and protein concentrations, suggesting that deactivation of cysteines was not responsible for thermal denaturation.     

Exploring the active site of tripeptidyl-peptidase II through studies of pH dependence of reaction kinetics

Tripeptidyl-peptidase II (TPP II) is a subtilisin-like serine protease which forms a large enzyme complex (> 4 MDa). It is considered a potential drug target due to its involvement in specific physiological processes. However, information is scarce concerning the kinetic characteristics of TPP II and its active site features, which are important for design of efficient inhibitors. To amend this, we probed the active site by determining the pH dependence of TPP II catalysis. Access to pure enzyme is a prerequisite for kinetic investigations and herein we introduce the first efficient purification system for heterologously expressed mammalian TPP II. The pH dependence of kinetic parameters for hydrolysis of two different chromogenic substrates, Ala-Ala-Phe-pNA and Ala-Ala-Ala-pNA, was determined for murine, human and Drosophila melanogaster TPP II as well as mutant variants thereof. The investigation demonstrated that TPP II, in contrast to subtilisin, has a bell-shaped pH dependence of k_cat^app/K_M probably due to deprotonation of the N-terminal amino group of the substrate at higher pH. Since both the K_M and k_cat^app are lower for cleavage of AAA-pNA than for AAF-pNA we propose that the former can bind non-productively to the active site of the enzyme, a phenomenon previously observed with some substrates for subtilisin. Two mutant variants, H267A and D387G, showed bell-shaped pH-dependence of k_cat^app, possibly due to an impaired protonation of the leaving group. This work reveals previously unknown differences between TPP II orthologues and subtilisin as well as features that might be conserved within the entire family of subtilisin-like serine peptidases.     

The Multistep Process of Homotypic Neutrophil Aggregation: A Review of the Molecules and Effects of Hydrodynamics

Homotypic adhesion of neutrophils stimulated with chemoattractant is analogous to capture on vascular endothelium in that both processes are supported by L-selectin and β₂-integrin adhesion receptors. Under hydrodynamic shear, cell adhesion requires that receptors bind sufficient ligand over the duration of intercellular contact to withstand the hydrodynamic stresses. Using cone and plate viscometry to apply a uniform linear shear field to suspensions of neutrophils and flow cytometry to quantitate the size distribution of aggregates formed over the time course of formyl peptide stimulation, we conducted a detailed examination of the affect of shear rate and shear stress on the kinetics of cell aggregation. The efficiency of aggregate formation was fit from a mathematical model based on Smoluchowski's two-body collision theory. Over a range of venular shear rates (400–800 s^-1), β90% of the single cells are recruited into aggregates ranging from doublets to groupings larger than sextuplets. Adhesion efficiency fit to the kinetics of aggregation increased with shear rate from β20% at 100s^-1 to a maximum level of β80% at 400 s^-1. This increase to peak adhesion efficiency was dependent on L-selectin and β₂-integrin, and was resistant to shear stress up to β7 dyn/cm². When L-selectin was blocked with antibody, β₂-integrin (CD11a, b) supported adhesion at low shear rates (< 400 s^-1). Aggregates formed over the rapid phase of aggregation remain intact and resistant to shear up to 120 s. At the end of this plateau phase of stability, aggregates spontaneously dissociate back to singlets. The rate of cell disaggregation is linearly proportional to the applied shear rate. The binding kinetics of selectin and integrin appear to be optimized to function within discrete ranges of shear rate and stress, providing an intrinsic mechanism for the transition from neutrophil tethering to firm but reversible adhesion.     

Identifying protein aggregation mechanisms and quantifying aggregation rates from combined monomer depletion and continuous scattering

Parallel temperature initial rates (PTIR) from chromatographic separation of aggregating protein solutions are combined with continuous simultaneous multiple sample light scattering (SMSLS) to make quantitative deductions about protein aggregation kinetics and mechanisms. PTIR determines the rates at which initially monomeric proteins are converted to aggregates over a range of temperatures, under initial-rate conditions. Using SMSLS for the same set of conditions provides time courses of the absolute Rayleigh scattering ratio, I_R(t), from which a potentially different measure of aggregation rates can be quantified. The present report compares these measures of aggregation rates across a range of solution conditions that result in different aggregation mechanisms for anti-streptavidin (AS) immunoglobulin gamma-1 (IgG1). The results illustrate how the two methods provide complementary information when deducing aggregation mechanisms, as well as cases where they provide new mechanistic details that were not possible to deduce in previous work. Criteria are presented for when the two techniques are expected to give equivalent results for quantitative rates, the potential limitations when solution non-idealities are large, as well as a comparison of the temperature dependence of AS-IgG1 aggregation rates with published data for other antibodies.     

Biopolymer - Surfactant Interaction: 4 Kinetics of Binding of Cetyltrimethyl Ammonium Bromide with Gelatin,Hemoglobin, β-lactoglobulin and Lysozyme

Abstract

The binding of CTAB with the proteins, gelatin, hemoglobin, β-lactoglobulin and lysozyme follow first order kinetics and occurs either in two or three distinct stages. The number of stages depends on the overall configuration of the biopolymers. The denatured protein, gelatin has shown three-stage kinetics under all conditions, whereas the native proteins, hemoglobinn, β-lactoglobulin and lysozyme have exhibited two stage kinetics. Heat treated lysozyme in 8 mol dm-³ urea medium has also shown a two-stage kinetics. On the basis of non interacting binding sites on the proteins and independent sequential binding, the rates of reaction have been observed to increase with temperature and follow the trend k₁ >> k₂ > k₃. The interaction of CTA⁺ with the proteins is both electrostatic and hydrophobic. Hemoglobin has shown maximum reaction rate whereas, β-lactoglobulin has shown a minimum. The activation parameters for the kinetic process have exhibited almost non-variant Δ G? and Δ H? < T Δ S? The formation of activation complex in the Eyring model is entropy controlled so also the overall kinetics. An isokinetic entropy-enthalpy compensation phenomenon has been observed for the respective kinetic stages.     

Osmotic water permeability in glycoprotein containing liposomes

The kinetics of osmotic water permeability in proteoliposomes containing ₁-acid glycoprotein was investigated by means of stopped-flow spectrophotometry. A biphasic time-course of scattered light with time was registered. The rate constants calculated from fits to an exponential function in the first phase were proportional to the final medium osmolarity. The apparent second order rate constants K_app (Osm^-1 sec^-1) were determined at different glycoprotein concentrations in the original mixture for preparation of proteoliposomes. The value of K_app at lipid:glycoprotein weight ratio = 1 was plotted in Arrhenius coordinates. The calculated activation energy for water permeation through the lipid bilayer suggests that eventual channel mechanism may be involved due to the presence of glycoprotein molecule in the liposomes.     

Nuclear magnetic resonance studies of exchangeable protons. II. The solvent exchange rate of the indole nitrogen proton of tryptophan derivatives.

The exchange rate of the indole nitrogen proton with solvent water protons was measured as a function of pH and temperature for tryptophan, N-methyl tryptophan, glycyl tryptophan, tryptophanamide, tryptophylglycine, and tryptophylglycyl glycine. The nmr observation was by long-pulse Fourier transform methods, and kinetics were inferred from saturation recovery, H₂O transfer of saturation, and linewidth. There are observable differences between the rates of these compounds, but all are describable within a factor of two by specific acid and base-catalyzed rates k_H = 100 and k_OH = 10⁸ 1/mol-sec at 27°C. It is concluded that this behaviour is representative of this proton on the indole side chain in a random-chain peptide exposed to water.     

Diffusional and electrostatic effects on apparent kinetic parameters of reactions catalyzed by enzyme immobilized on the external surface of a support

Diffusional and electrostatic effects on the apparent maximum reaction rate V_m^app and the apparent Michaelis constant K_m^app were investigated theoretically for a system in which an enzyme immobilized on the external surface of a solid support catalyzes a reaction according to Michaelis-Menten kinetics. In such a system, the dependence of V_m^app and K_m^app on the substrate concentration can be expressed analytically. When the support and substrate carry charges of the same sign, resulting in a repulsive force between them, both V_m^app and K_m^app decrease with increasing substrate concentration, but they never decrease below the respective intrinsic values. On the other hand, when the support and substrate carry charges of opposite sign and therefore an attractive force occurs, V_m^app decreases towards its intrinsic value, while K_m^app decreases to values below its intrinsic value in the region of high substrate concentration.     

An NMR method for studying the kinetics of metal exchange in biomolecular systems

The kinetics of lanthanide (III) exchange for calcium(II) in the C-terminal EF-hand of the protein calbindin D_9khave been studied by one-dimensional (1D) stopped-flow NMR. By choosing a paramagnetic lanthanide (Ce³⁺), kinetics in the sub-second range can be easily measured. This is made possible by the fact that (i) the kinetic behaviour of hyperfine shifted signals can be monitored in 1D NMR and (ii) fast repetition rates can be employed because these hyperfine shifted signals relax fast. It is found that the Ce³⁺-Ca²⁺exchange process indeed takes place on a sub-second timescale and can be easily monitored with this technique. As the rate of calcium-cerium substitution was found not to depend on the presence of excess calcium in solution, the kinetics of the process were interpreted in terms of a bimolecular associative mechanism, and the rate constants extracted. Interestingly, the dissociative mechanism involving the apoform of the protein, which is generally assumed for metal ion exchange at protein binding sites, was not in agreement with our data.