Conformational species of gramicidin A in non-polar solvent. A kinetic and thermodynamic treatment in the absence and presence of phosphatidylcholine as studied by high-performance liquid chromatography

A kinetic and thermodynamic study has been carried out to characterize quantitatively the conformational equilibrium of gramicidin A (GA) in tetrahydrofuran at different peptide concentrations in the absence and presence of egg yolk phosphatidylcholine by using size-exclusion high-performance liquid chromatographic analysis. In the absence of lipid, the experimental data fit a simple dimer-monomer equilibrium, the rate and equilibrium constants for the dissociation process being (1.6 +/- 0.7) X 10(-7) s-1 and (8.5 +/- 0.3) X 10(-6) M, respectively. A higher extent of monomerization and a decrease in the time required for reaching equilibrium are detected in the presence of phospholipid, the kinetic and thermodynamic effects depending on both lipid and GA concentrations. In order to account for these observations a cyclic equilibrium mechanism is proposed which is analysed in terms of four conformational species, namely, free monomer, free dimer, lipid-bound monomer and lipid-bound dimer. The results obtained are discussed in relation to recent literature data on lipid-protein interactions.     

Enantioselective interactions of dextromethorphan and levomethorphan with the alpha 3 beta 4-nicotinic acetylcholine receptor: comparison of chromatographic and functional data

The enantioselectivity of the interaction of dextromethorphan (DM) and levomethorphan (LM) with an immobilized alpha 3 beta 4 subtype of the nicotinic acetylcholine receptor (nAChR) liquid chromatographic stationary phase has been compared to DM- and LM-induced non-competitive blockade of nicotine-stimulated 86Rb(+) efflux from cells expressing the alpha 3 beta 4-nAChR. The association rate constants (k(on)) and dissociation rate constants (k(off)) for the formation of the DM and LM complexes with the nAChR were determined using non-linear chromatographic techniques and the k(off) value for DM (1.01+/-0.01 s(-1)) was significantly lower than the k(off) for LM (1.55+/-0.002s(-1)) while the k(on) values did not significantly differ (23.66+/-0.61 and 18.61+/-0.36 microM(-1)s(-1), respectively). In thermodynamic studies using the van't Hoff approach, the enthalpy change (Delta H degrees) of the DM-nAChR complex was 330 calmol(-1) more stable than the LM-nAChR complex, while there was no significant difference in the entropy change (DeltaS degrees ). In the functional in vitro cell-based studies, there was no significant difference in the observed IC(50) values for DM (10.1+/-1.01 microM) and LM (10.9+/-1.08 microM), but the recovery from the DM-induced blockade was slower than the recovery from LM-induced blockade; after 7 min: 38.25+/-15.46% recovery from DM blockade, 63.30+/-16.08% from LM blockade; after 4h: 76.20+/-4.51% recovery from DM blockade and 93.12+/-8.76% from LM blockade. The enantioselective differences in the functional effects are consistent with the chromatographic and thermodynamic data and indicate that this difference is due to increased stability of the DM-nAChR complex. The results suggest that the chromatographic approach can be used to probe the interaction of non-competitive inhibitors (NCIs) with nAChRs and to predict relative duration of functional blockades.     

Steady-state kinetics of electron transfer through cytochrome chain of uncoupled submitochondrial particles. II. Influence of pH on kinetics of electron transfer.

pH Dependences of steady-state kinetic parameters of cytochrome chains of submitochondrial particles have been studies. It has been shown that the lifetimes of activated states (tau) of the pairs of cytochromes b leads to c1 and a leads to a3 have different pH dependences; those for the c1 leads to c and c leads to a cytochrome pairs being similar. The rate constants for the non-activated state of the respiratory chains decreased for the b leads to c1 pair and increased for the a leads to a3 pair when the pH value was increased. The values of pK calculated from these dependences for the pairs b leads to c1 and a leads to a3 were 7.2 and 8.9, respectively. It has been supposed that the ratio of activated to non-activated electron carriers may be controlled by the local pH value in the mitochondrial membrane, the latter being dependent upon the rate of electron transfer. The kinetic model based on this assumption allows one to explain the experimental dependences on pH of the rate constants for cytochromes b leads to c, and a leads to a3. The values of the diffusion rate constants for H+ and OH- ions in the mitochondrial membrane estimated from these kinetic data obtained in this study were 10(4)--10(5) s-1 and 10(2)--10(3) s-1, respectively.     

Considering the oligonucleotides secondary structures at thermodynamic and kinetic analysis of the DNA-duplexes formation

The influence of secondary structures of DNA oligonucleotides on thermodynamics and kinetics at the formation of their bimolecular complexes (duplexes) has been studied. The models considering inherent secondary structures of duplex components and their influence on quantitative thermodynamic and kinetic characteristics of the duplexes have been developed. The values of thermodynamic impacts given by individual structural elements of the double helix have been shown to depend on hairpin structuring of the duplex components. The "concentration" method to consider oligonucleotides intramolecular structure with thermodynamic parameters of bimolecular duplex formation has been proposed. According to stop-flow measurements, the observed values of association and dissociation constants are influenced by the presence of inherent structures in duplex components. The influence observed is increased with the lowering of the sample temperature. The analysis of experimental data involving the developed models provides the possibility to determine "proper" kinetic constants for the helix-to-coil transition. The difference between observed and calculated rate constants can amount up to two or more orders of magnitude.     

The mechanism of skeletal muscle myosin ATPase. The mechanism of ADP interaction

The mechanism of interaction between ADP and the myosin active center has been studied using a transient kinetic technique. The results show that the interaction of ADP with the myosin active center is a homogeneous process independent of the association state of the active centers; namely, whether ADP interacts with the monomeric myosin subfragment-1, or with the dimeric forms heavy meromyosin and myosin. The kinetics of the interaction conforms to a simple two-step reaction mechanism for ADP binding. The kinetic and thermodynamic constants for this mechanism have been determined. In addition, analysis of the binding isotherm indicates that the two active sites in heavy meromyosin and myosin function as identical and independent sites.     

Considering the oligonucleotide secondary structures in thermodynamic and kinetic analysis of DNA duplex formation

The influence of secondary structures of DNA oligonucleotides on thermodynamics and kinetics at the formation of their bimolecular complexes (duplexes) has been studied. The models considering inherent secondary structures of duplex components and their influence on quantitative thermodynamic and kinetic characteristics of the duplexes have been developed. The values of thermodynamic impacts given by individual structural elements of the double helix have been shown to depend on hairpin structuring of the duplex components. The «concentration» method to consider oligonucleotides intramolecular structure with thermodynamic parameters of bimolecular duplex formation has been proposed. According to stop-flow measurements, the observed values of association and dissociation constants are influenced by the presence of inherent structures in duplex components. The influence observed is increased with the lowering of the sample temperature. The analysis of experimental data involving the developed models provides the possibility to determine «proper» kinetic constants for the helix-to-coil transition. The difference between observed and calculated rate constants can amount up to two or more orders of magnitude.     

Global parameter optimization for cardiac potassium channel gating models.

Quantitative ion channel model evaluation requires the estimation of voltage dependent rate constants. We have tested whether a unique set of rate constants can be reliably extracted from nonstationary macroscopic voltage clamp potassium current data. For many models, the rate constants derived independently at different membrane potentials are not unique. Therefore, our approach has been to use the exponential voltage dependence predicted from reaction rate theory (Stevens, C. F. 1978. Biophys. J. 22:295-306; Eyring, H., S. H. Lin, and S. M. Lin. 1980. Basic Chemical Kinetics. Wiley and Sons, New York) to couple the rate constants derived at different membrane potentials. This constrained the solution set of rate constants to only those that also obeyed this additional set of equations, which was sufficient to obtain a unique solution. We have tested this approach with data obtained from macroscopic delayed rectifier potassium channel currents in voltage-clamped guinea pig ventricular myocyte membranes. This potassium channel has relatively simple kinetics without an inactivation process and provided a convenient system to determine a globally optimized set of voltage-dependent rate constants for a Markov kinetic model. The ability of the fitting algorithm to extract rate constants from the macroscopic current data was tested using "data" synthesized from known rate constants. The simulated data sets were analyzed with the global fitting procedure and the fitted rate constants were compared with the rate constants used to generate the data. Monte Carlo methods were used to examine the accuracy of the estimated kinetic parameters. This global fitting approach provided a useful and convenient method for reliably extracting Markov rate constants from macroscopic voltage clamp data over a broad range of membrane potentials. The limitations of the method and the dependence on initial guesses are described.     

Thermodynamic and kinetic characterisation of individual haems in multicentre cytochromes c3

The characterisation of individual centres in multihaem proteins is difficult due to the similarities in the redox and spectroscopic properties of the centres. NMR has been used successfully to distinguish redox centres and allow the determination of the microscopic thermodynamic parameters in several multihaem cytochromes c(3) isolated from different sulphate-reducing bacteria. In this article we show that it is also possible to discriminate the kinetic properties of individual centres in multihaem proteins, if the complete microscopic thermodynamic characterisation is available and the system displays fast intramolecular equilibration in the time scale of the kinetic experiment. The deconvolution of the kinetic traces using a model of thermodynamic control provides a reference rate constant for each haem that does not depend on driving force and can be related to structural factors. The thermodynamic characterisation of three tetrahaem cytochromes and their kinetics of reduction by sodium dithionite are reported in this paper. Thermodynamic and kinetic data were fitted simultaneously to a model to obtain microscopic reduction potentials, haem-haem and haem-proton interacting potentials, and reference rate constants for the haems. The kinetic information obtained for these cytochromes and recently published data for other multihaem cytochromes is discussed with respect to the structural factors that determine the reference rates. The accessibility for the reducing agent seems to play an important role in controlling the kinetic rates, although is clearly not the only factor.     

Ionization dependence of camphor binding and spin conversion of the complex between cytochrome P-450 and camphor. Kinetic and static studies at sub-zero temperatures.

The kinetic rate constants of formation and dissociation of the cytochrome-P-450 - camphor complex (Fe3+-RH) have been obtained by low-temperature (+ 5 degrees C to -20 degrees C) stopped-flow experiments. Simiarly the high-spin/low-spin equilibrium of this complex has been studied as a function of temperature and protonic activity. Both the camphor-binding mechanism and the high-spin/low-spin thermodynamic parameters of Fe3+-RH depend on the protonic activity of the medium in the physiological pH range. The binding rate constants are shown to depend on the ionization of a residue of the protein, probably a histidine. Linear enthalpy-entropy compensation is observed for the camphor binding as well as for the spin-state transition. A camphor-binding-induced change of the electrostatic potential is discussed.     

Modelling, steady state analysis and optimization of the catalytic efficiency of the triosephosphate isomerase

In the present work we have modelled and optimized the reaction mechanism of the triose phosphate isomerase (TIM) enzyme (E.C. 5.3.1.1). For this purpose we have used an approach that combines the S-system representation within the power law formalism and linear programming techniques. By this means we have explored those rate constants whose alterations are likely to improve the catalytic efficiency of the enzyme and investigated the available room for optimization in different metabolic conditions. The role and plausibility of the different types of mutations on the evolution of this enzyme have also been considered. Steady state sensitivity analysis was carried out and a new set of aggregated logarithmic gains was defined in order to quantify the responses of the system to changes in groups of rate constants that could be explained in terms of mutations affecting the catalytic properties of the enzyme. Evaluation of these logarithmic gains at different levels of saturation and disequilibrium ratios enabled us to reach conclusions about the meaning and role of the diffusion limitation terms. The catalytic efficiency of the monoenzymatic system was optimized through changes in the kinetic rate constants within different sets of restrictions ranging from thermodynamic or kinetic to evolutionary ones. Results showed that, at very different conditions, there is still room for improvement in the TIM enzyme. Thus, in a wide range of metabolically significant values of the disequilibrium ratio there is a minimal variation in the optimal profile that yields 2.1 times the velocity of the basal states. Though most of this increase is accounted for by the increase of the second order constants (that could have already reached a theoretical maximum) significant increases (10%) in catalytic efficiencies are obtained by changes of the internal steps only. Besides these new findings our optimization approach has been able to reproduce results obtained with other approaches.     

Reaction-kinetic parameters of glycidamide as determinants of mutagenic potency

Values for reaction-kinetic parameters of electrophiles can be used to predict mutagenic potency. One approach employs the Swain-Scott relationship for comparative kinetic studies of electrophilic agents reacting with nucleophiles. In this way glycidamide (GA), the putatively mutagenic/carcinogenic metabolite of acrylamide, was assessed by determining the rates of reaction with different nucleophiles. The rate constants (kNu) were determined using the "supernucleophile" cob(I)alamin [Cbl(I)] as an analytical tool. The Swain-Scott parameters for GA were compared with those of ethylene oxide (EO). The substrate constants, s values, for GA and for EO were found to be 1.0 and 0.93, respectively. The reaction rates at low values of nucleophilic strength (n=1-3), corresponding to oxygens in DNA, were determined to be 2-3.5 times higher for GA compared to EO. GA was also more reactive than EO towards other nucleophiles (n=0-6.4). The mutagenic potency of GA was determined in Chinese hamster ovary cells (hprt mutations in CHO-AA8 cells per dose unit with gamma-radiation as reference standard). The potency of GA was estimated to be about three mutations per 10(5) cells and mMh corresponding to about 40 rad-equ./mMh. A preliminary comparison of the mutagenic potency (per mMh and as rad-equivalents) of GA and EO shows an approximately seven times higher potency for GA. A higher mutagenic potency of GA compared to EO is compatible with expectation from reaction-kinetic data of the two compounds. The data confirmed that GA is not a strong mutagen, which is in line with what is expected for simple oxiranes. The present study shows the value of cob(I)alamin for the determination of reaction-kinetic parameters and their use for prediction of mutagenic potency.     

Analysis of the interaction of gallic acid and myoglobin by UV-vis absorption spectroscopy

UV-vis absorption spectroscopy has been used to analyze the interaction of myoglobin (Мb) and gallic acid (GA). The binding constants (4.38 × 10⁴ M^–1 at 298.15 K and 0.42 × 10⁴ М^–1 at 308.15K), the number of binding sites (h = 1.0), and the thermodynamic parameters of binding (ΔH, ΔS, and ΔG) have been determined. Hydrogen bonds have been shown to play a major role in the stabilization of the GA–Мb complexes. GA binding led to slight changes in the electronic state of the heme ring of the protein.     

Kinetic study of the reduction mechanism for Desulfovibrio gigas cytochrome c3.

The kinetic aspects of the reduction process in cytochrome c3 from Desulfovibrio gigas have been investigated over a wide range of pH values ranging between pH 5.8 and pH 9.8. The data have been analyzed in the framework of an I2H4 interaction network coupled to a proton-linked equilibrium between two tertiary structures (Cornish-Bowden, A. & Koshland, D.E. Jr (1970) J. Biol. Chem. 245, 6241-6250). The kinetic rate constants for the reduction of the four hemes for the two tertiary conformations have been characterized in the framework of the thermodynamic network obtained from the equilibrium analysis (Coletta, M., Catarino, T., LeGall, J.J. & Xavier, A.V. (1991) Eur. J. Biochem. 202, 1101-1106). The intrinsic reduction rate constants determined by reaction with sodium dithionite for two hemes (namely heme 4 and heme 1) are significantly faster than those for the other two heme residues. In view of the equilibrium redox properties, heme 4 (with the fastest reduction rate) may then work as the kinetic electron-capturing site for the electrons from sodium dithionite. The transfer to hemes 2 and 3 then occurs by virtue of their free-energy levels at equilibrium. At our experimental conditions, there is also transfer of electrons to hemes 2 and 3 from heme 1, which is reduced at a slower rate than heme 4, thus contributing to the biphasic kinetics observed for the overall process. The kinetic parameters obtained are discussed in terms of the mechanism proposed for the coupling between the electron and proton transfer, as induced by the heme/heme cooperativity network.     

Determination of binding constant of transcription factor AP-1 and DNA. Application of inhibitors.

The equilibrium binding and association kinetics of the fos-jun dimer (basic and leucine zipper domain) to the AP-1 DNA were studied using a quantitative assay. The basic-region and leucine zipper (bZip) domain of c-fos was expressed as a fusion protein with glutathione S-transferase, and it was bound to glutathione-agarose. The GST-fused fos bZip region was allowed to form a heterodimer with the bZip domain of c-jun, to which radiolabeled AP-1 nucleotides were added. After thorough washing, the gel-bound radioactivity was counted. The binding and dissociation rate constants (k(1) and k-(1)) of the fos-jun dimer and DNA could be obtained from a time-course experiment. The association binding constant (K(1)) was determined using both a thermodynamic equation and kinetic parameters. Nordihydroguaiaretic acid (NDGA), momordin I, natural product inhibitors of the fos-jun/DNA complex formation, was applied to this jun-GST-fused fos system and it was found to decrease the apparent equilibrium binding of dimer and DNA. The thermodynamic constant of dimer and inhibitor binding was also determined.     

Horseradish peroxidase catalyzed hydroxylation of phenol: II. Kinetic model

A numeric kinetic model of the horseradish peroxidase catalyzed hydroxylation of phenol is proposed to complete the previous thermodynamic analysis. As previously stated, the basic role of HRP is to catalyze the production of DHF* radicals. These further form hydroxyl radicals that hydroxylate phenol via noncatalyzed reactions. The transient differential equations of the model are solved numerically. Several kinetic constants are adjusted to fit basic experimental data. This set of values is then kept constant to simulate additive experiments carried out under different conditions. Predictions of the model concerning the effects of HRP concentration, temperature variation, and presence of catalase and superoxide dismutase are consistent with the experimental results. The quantitative kinetic approach consequently fully confirmed the previous thermodynamic conclusions.     

Kinetic and thermodynamic aspects of lipid translocation in biological membranes

A theoretical analysis of the lipid translocation in cellular bilayer membranes is presented. We focus on an integrative model of active and passive transport processes determining the asymmetrical distribution of the major lipid components between the monolayers. The active translocation of the aminophospholipids phosphatidylserine and phosphatidylethanolamine is mathematically described by kinetic equations resulting from a realistic ATP-dependent transport mechanism. Concerning the passive transport of the aminophospholipids as well as of phosphatidylcholine, sphingomyelin, and cholesterol, two different approaches are used. The first treatment makes use of thermodynamic flux-force relationships. Relevant forces are transversal concentration differences of the lipids as well as differences in the mechanical states of the monolayers due to lateral compressions. Both forces, originating primarily from the operation of an aminophospholipid translocase, are expressed as functions of the lipid compositions of the two monolayers. In the case of mechanical forces, lipid-specific parameters such as different molecular surface areas and compression force constants are taken into account. Using invariance principles, it is shown how the phenomenological coefficients depend on the total lipid amounts. In a second approach, passive transport is analyzed in terms of kinetic mechanisms of carrier-mediated translocation, where mechanical effects are incorporated into the translocation rate constants. The thermodynamic as well as the kinetic approach are applied to simulate the time-dependent redistribution of the lipid components in human red blood cells. In the thermodynamic model the steady-state asymmetrical lipid distribution of erythrocyte membranes is simulated well under certain parameter restrictions: 1) the time scales of uncoupled passive transbilayer movement must be different among the lipid species; 2) positive cross-couplings of the passive lipid fluxes are needed, which, however, may be chosen lipid-unspecifically. A comparison of the thermodynamic and the kinetic approaches reveals that antiport mechanisms for passive lipid movements may be excluded. Simulations with kinetic symport mechanisms are in qualitative agreement with experimental data but show discrepancies in the asymmetrical distribution for sphingomyelin.     

Affinity chromatographic purification of papain.

The reinvestigation of the affinity chromatographic method of purifying papain has been carried out. It has been reported that papain could be purified by taking advantage of the affinity of the enzyme for the insolubilized peptide inhibitor, agarose-Gly-Gly-Tyr(Bz)-Arg. Using pure tetrapeptide obtained commercially and standard coupling procedures, a significant purification of papain could not be achieved. Both active and nonactivatible enzyme bound to a column prepared in this manner were eluted together by the use of deionized water. An affinity medium with properties similar to those reported by Blumberg et al. was obtained by removal of the benzyl group on tyrosine prior to coupling with agarose. The deprotected tetrapeptide was also synthesized by an independent route and inhibition constants for the binding of the protected and deprotected tetrapeptide to papain were determined in kinetic experiments.     

Inter-flavin electron transfer in cytochrome P450 reductase - effects of solvent and pH identify hidden complexity in mechanism

This study on human cytochrome P450 reductase (CPR) presents a comprehensive analysis of the thermodynamic and kinetic effects of pH and solvent on two- and four-electron reduction in this diflavin enzyme. pH-dependent redox potentiometry revealed that the thermodynamic equilibrium between various two-electron reduced enzyme species (FMNH*,FADH*; FMN,FADH2; FMNH2,FAD) is independent of pH. No shift from the blue, neutral di-semiquinone (FMNH*,FADH*) towards the red, anionic species is observed upon increasing the pH from 6.5 to 8.5. Spectrophotometric analysis of events following the mixing of oxidized CPR and NADPH (1 to 1) in a stopped-flow instrument demonstrates that the establishment of this thermodynamic equilibrium becomes a very slow process at elevated pH, indicative of a pH-gating mechanism. The final level of blue di-semiquinone formation is found to be pH independent. Stopped-flow experiments using excess NADPH over CPR provide evidence that both pH and solvent significantly influence the kinetic exposure of the blue di-semiquinone intermediate, yet the observed rate constants are essentially pH independent. Thus, the kinetic pH-gating mechanism under stoichiometric conditions is of no significant kinetic relevance for four-electron reduction, but rather modulates the observed semiquinone absorbance at 600 nm in a pH-dependent manner. The use of proton inventory experiments and primary kinetic isotope effects are described as kinetic tools to disentangle the intricate pH-dependent kinetic mechanism in CPR. Our analysis of the pH and isotope dependence in human CPR reveals previously hidden complexity in the mechanism of electron transfer in this complex flavoprotein.