Effects of single-tryptophan mutations on R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) is an R-plasmid-encoded enzyme that confers clinical resistance to the antibacterial drug trimethoprim. This enzyme shows no sequence or structural homology to the chromosomal DHFRs. The active form of the protein is a homotetramer possessing D(2) symmetry and a single active-site pore. Two tryptophans occur per monomer: W38 and its symmetry-related residues (W138, W238, and W338) occur at the dimer-dimer interfaces, while W45 and its symmetry-related partners (W145, W245, and W345) occur at the monomer-monomer interfaces. Two single-tryptophan mutant genes were constructed to determine the structural and functional consequences of four mutations per tetramer. The W45F mutant retains full enzyme activity and the fluorescence environment of the unmutated W38 residues clearly monitors ligand binding and a pH dependent tetramer right harpoon over left harpoon 2 dimers equilibrium. In contrast, four simultaneous W38F mutations at the dimer-dimer interfaces result in tetramer destabilization. The ensuing dimer is relatively inactive, as is dimeric wild-type R67 DHFR. A comparison of emission spectra indicates the fluorescent signal of wild-type R67 DHFR is dominated by the contribution from W38. Equilibrium unfolding/folding curves at pH 5.0, where all protein variants are dimeric, indicate the environment monitored by the W38 residue is slightly less stable than the environment monitored by the W45 residue.     

Interaction of streptavidin-based peptide-MHC oligomers (tetramers) with cell-surface TCRs

The binding of oligomeric peptide-MHC (pMHC) complexes to cell surface TCR can be considered to approximate TCR-pMHC interactions at cell-cell interfaces. In this study, we analyzed the equilibrium binding of streptavidin-based pMHC oligomers (tetramers) and their dissociation kinetics from CD8(pos) T cells from 2C-TCR transgenic mice and from T cell hybridomas that expressed the 2C TCR or a high-affinity mutant (m33) of this TCR. Our results show that the tetramers did not come close to saturating cell-surface TCR (binding only 10-30% of cell-surface receptors), as is generally assumed in deriving affinity values (K(D)), in part because of dissociative losses from tetramer-stained cells. Guided by a kinetic model, the oligomer dissociation rate and equilibrium constants were seen to depend not only on monovalent association and dissociation rates (k(off) and k(on)), but also on a multivalent association rate (μ) and TCR cell-surface density. Our results suggest that dissociation rates could account for the recently described surprisingly high frequency of tetramer-negative, functionally competent T cells in some T cell responses.     

Requirements for reliable determination of binding affinity constants by kinetic approach

The nonspecific binding (equilibrium coefficient kn) of ligand (L) and/or the incomplete recovery (alpha < 1) of the receptor-ligand (RL) complex in binding measurements, could hamper accurate determination of the association and dissociation rate constants of the R/L system. For the simplest model of R/L interaction, characterized by a bimolecular association process (rate constant k1) and a monomolecular dissociation process (rate constant k2), the consequences of kn and/or alpha neglect on k1 and k2 determination were investigated. Various situations that are especially relevant for k1 determination, were examined in which nonspecific binding was: (i) negligible relative to specific binding, or (ii) developed progressively or very rapidly in association kinetics. When only the initial kinetic phase was used, according to the situation (i.e. the nonspecific binding characteristics, and the fact that kn and/or alpha were or were not taken into account to correct the binding measurements), k1 could be accurately determined or generally slightly overestimated or slightly underestimated (in the two latter cases by factors involving mainly kn and/or alpha but not the R concentration or the R/L equilibrium association constant, K), whereas k2 should always be fairly well estimated. Consequently, for the simplest R/L systems, the k1/k2 ratio derived from such kinetic experiments should be much less susceptible to substantial underestimation than K derived from R saturation experiments [Borgna, J. Steroid Biochem. Mol. Biol. (2004)]. Kinetic experiments could also be more appropriate than R saturation experiments to detect cooperative--positive or negative--binding of L to R.     

Kinetics of gonadotropin binding by receptors of the rat testis. Analysis by a nonlinear curve-fitting method.

The kinetics of the reaction between human chorionic gonadotropin (hCG) and specific gonadotropin receptors in the rat testis were determined at 24 and 37 degrees, over a wide range of hormone concentrations. Hormone concentrations were corrected for the binding activity of the (-125I)hCG tracer preparations. Analysis of the experimental data was performed with an interactive nonlinear curve fitting program, based upon the second-order chemical kinetic differential equation. The mean values for the association rate constant (k1) were 4.7 x 10-7 M-1 min-1 at 24 degrees, and 11.0 x 10-7 M-1 min-1 at 37 degrees. At both temperatures, the values of kl were independent of hormone concentration. Initial dissociation rates were consistent with first order kinetics, with dissociation rate constant (k2) of 1.7 x 10 minus -3 and 4.6 x 10 minus -3 min minus -1 at 24 and 37 degrees, respectively. When studied over longer periods at 24 degrees, the dissociation process appeared to be multiexponential. The kinetics of degradation of (-125I)hCG and receptors were determined at both temperatures, and a mathematical model was developed by modification of the second-order chemical kinetic differential equation to take these factors into account. The application of such a model to hCG kinetic binding data demonstrated that reactant degradation had little significant effect on the derivation of the association rate constant (k1), but caused significant overestimation of the dissociation rate constant (k2) values derived from association experiments. The model was also applied by computer simulation to a theoretical analysis of the effects of degradation of free hormone and receptor sites upon kinetic and steadystate binding data. By this method, the initial velocities of hormone binding were shown to be less affected by degradation than the steady-state levels of hormone-receptor complex. Also, reactant degradation in simulated steady-state experiments caused an underestimate of the apparent equilibrium association constant, but had relatively less effect on the determination of binding site concentration.     

Towards understanding the mechanisms of molecular recognition by computer simulations of ligand-protein interactions

The thermodynamic and kinetic aspects of molecular recognition for the methotrexate (MTX)-dihydrofolate reductase (DHFR) ligand-protein system are investigated by the binding energy landscape approach. The impact of 'hot' and 'cold' errors in ligand mutations on the thermodynamic stability of the native MTX-DHFR complex is analyzed, and relationships between the molecular recognition mechanism and the degree of ligand optimization are discussed. The nature and relative stability of intermediates and thermodynamic phases on the ligand-protein association pathway are studied, providing new insights into connections between protein folding and molecular recognition mechanisms, and cooperativity of ligand-protein binding. The results of kinetic docking simulations are rationalized based on the thermodynamic properties determined from equilibrium simulations and the shape of the underlying binding energy landscape. We show how evolutionary ligand selection for a receptor active site can produce well-optimized ligand-protein systems such as MTX-DHFR complex with the thermodynamically stable native structure and a direct transition mechanism of binding from unbound conformations to the unique native structure.     

Thermodynamic studies of the reversible association of Escherichia coli ribosomal subunits.

The kinetics of association of Escherichia coli 30S and 50S ribosomal subunits have been carried out as a function of temperature after a magnesium jump from 1.5 to 3 mM. Turbidimetric recordings combined with a stopped-flow apparatus were used to follow the kinetics. The data show that the rates of formation and dissociation of the 70S particles at 3 mM Mg2+ and +25 degrees C were, respectively: k2 = 10(5) M-1 s-1, k1 = 4,5 X 10(-3) s-1; lowering the temperature decreases the rate constants with activation energies equal to E2 = 7.5 kcal/mol, E1 = 26.5 kcal/mol and enhances the association equilibrium towards the 70S species with an enthalpy change (delta H degrees assoc = -19.9 kcal/mol) dominant over the entropy change (delta S degrees assoc = -33 cal/(deg mol)). These thermodynamic parameters were compared to those obtained from studies on the interactions of codon-anticodon in yeast phenylalanine transfer RNA as well as of ribooligonucleotides. The kinetic and thermodynamic data are shown to be consistent with 16S-23S RNA interaction.     

Effect of multiple symmetries on the association of R67 DHFR subunits bearing interfacial complementing mutations

It was shown previously that complementation could be a powerful mean to probe protein-protein interactions in the normally tetrameric R67 DHFR. Indeed, mixing complementing inactive dimeric mutants produced active heterotetramers. This approach turned a homo-oligomer into a hetero-oligomer and thus allowed the use of combinatorial assays, a subtle analysis of the association forces, and a precise determination of the equilibrium dissociation constants (K(D)) by titrimetry. However, for some of the complementing pairs, the experimental data implied multiple equilibria involving heterodimers, although no monomers could be detected. Thus, the reactions involved had to be identified to elaborate a suitable model to determine the K(D) of those pairs correctly. That model suggested that homodimers associated rapidly before the protomers could be redistributed in a multiple equilibrium system. Kinetic data confirmed that view. The association data at equilibrium were analyzed by multiple curve fitting with all plausible combinations of parameters. This gave a confidence interval for K(D) that is safer than the usual 67% or 90% confidence interval. Finally, the K(D) of one specific reaction, the dissociation of a heterotetramer with the relevant symmetry into two homodimers could be determined with the relevant model for each complementing pair, although multiple equilibria were present. These K(D) can thus be used as a set of references data to test and improve theoretical methods such as association free energy calculations.     

Novel mechanism-based substrates of dihydrofolate reductase and the thermodynamics of ligand binding: A comparison of theory and experiment for 8-methylpterin and 6,8-dimethylpterin

Molecular dynamics simulation and free energy perturbation techniques have been used to study the relative binding free energies of the designed mechanism-based pterins, 8-methylpterin and 6,8-dimethylpterin, to dihydrofolate reductase (DHFR), with co-factor nicotinamide adenine dinucleotide phosphate (NADPH). The calculated free energy differences suggest that DHFR.NADPH.6,8-dimethylpterin is thermodynamically more stable than DHFR.NADPH.8-methylpterin by 2.4 kcal/mol when the substrates are protonated and by 1.3 kcal/mol when neutral. The greater binding strength of 6,8-dimethylpterin may be attributed largely to hydration effects. In terms of an appropriate model for the pH-dependent kinetic mechanism, these differences can be interpreted consistently with experimental data obtained from previous kinetic studies, i.e., 6,8-dimethylpterin is a more efficient substrate of vertebrate DHFRs than 8-methylpterin. The kinetic data suggest a value of 6.6 ± 0.2 for the pK_a of the active site Glu-30 in DHFR.NADPH. We have also used experimental data to estimate absolute values for thermodynamic dissociation constants of the active (i.e., protonated) forms of the substrates: these are of the same order as for the binding of folate (0.1–10 μM). The relative binding free energy calculated from the empirically derived dissociation constants for the protonated forms of 8-methylpterin and 6,8-dimethylpterin is 1.4 kcal/mol, a value which compares reasonably well with the theoretical value of 2.4 kcal/mol. © 1993 Wiley-Liss, Inc.     

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.     

Direct monitoring of allosteric recognition of type IIE restriction endonuclease EcoRII

EcoRII is a homodimer with two domains consisting of a DNA-binding N terminus and a catalytic C terminus and recognizes two specific sequences on DNA. It shows a relatively complicated cleavage reaction in bulk solution. After binding to either recognition site, EcoRII cleaves the other recognition site of the same DNA (cis-binding) strand and/or the recognition site of the other DNA (trans-binding) strand. Although it is difficult to separate these two reactions in bulk solution, we could simply obtain the binding and cleavage kinetics of only the cis-binding by following the frequency (mass) changes of a DNA-immobilized quartz-crystal microbalance (QCM) responding to the addition of EcoRII in aqueous solution. We obtained the maximum binding amounts (Deltam(max)), the dissociation constants (K(d)), the binding and dissociation rate constants (k(on) and k(off)), and the catalytic cleavage reaction rate constants (k(cat)) for wild-type EcoRII, the N-terminal-truncated form (EcoRII N-domain), and the mutant derivatives in its C-terminal domain (K263A and R330A). It was determined from the kinetic analyses that the N-domain, which covers the catalytic C-domain in the absence of DNA, preferentially binds to the one DNA recognition site while transforming EcoRII into an active form allosterically, and then the secondary C-domain binds to and cleaves the other recognition site of the DNA strand.     

Evolutional design of a hyperactive cysteine- and methionine-free mutant of Escherichia coli dihydrofolate reductase

We developed a strategy for finding out the adapted variants of enzymes, and we applied it to an enzyme, dihydrofolate reductase (DHFR), in terms of its catalytic activity so that we successfully obtained several hyperactive cysteine- and methionine-free variants of DHFR in which all five methionyl and two cysteinyl residues were replaced by other amino acid residues. Among them, a variant (M1A/M16N/M20L/M42Y/C85A/M92F/C152S), named as ANLYF, has an approximately seven times higher k(cat) value than wild type DHFR. Enzyme kinetics and crystal structures of the variant were investigated for elucidating the mechanism of the hyperactivity. Steady-state and transient binding kinetics of the variant indicated that the kinetic scheme of the catalytic cycle of ANLYF was essentially the same as that of wild type, showing that the hyperactivity was brought about by an increase of the dissociation rate constants of tetrahydrofolate from the enzyme-NADPH-tetrahydrofolate ternary complex. The crystal structure of the variant, solved and refined to an R factor of 0.205 at 1.9-angstroms resolution, indicated that an increased structural flexibility of the variant and an increased size of the N-(p-aminobenzoyl)-L-glutamate binding cleft induced the increase of the dissociation constant. This was consistent with a large compressibility (volume fluctuation) of the variant. A comparison of folding kinetics between wild type and the variant showed that the folding of these two enzymes was similar to each other, suggesting that the activity enhancement of the enzyme can be attained without drastic changes of the folding mechanism.     

Regulation of cAMP-dissociation kinetics in lactating rat mammary gland

The cAMP-dissociation kinetics of rat mammary gland cytosols are dependent upon the temperature of cAMP association. Dissociation rates (measured at pH 6.5, 24 degrees C) were biphasic (k = 0.08-0.23 min-1 and k = 0.02 min-1) and monophasic (k-1 = 0.02 min-1) after 0 degrees C and 24 degrees C association, respectively. The temperature-dependent change from an initial fast rate to an initial slow rate was observed at all concentrations of cAMP tested from 1 to 1000 nM. When the slow-dissociating site was associated with non-radioactive 8-bromo-cAMP, the dissociation rates of [3H]-cAMP from the remaining dissociating site was slow (k = 0.02 min-1) and fast (k = 0.05 min-1) at 24 degrees C and 0 degrees C associating rate can be converted to the slow-dissociating rate by warming. When 0.2 M sodium thiocyanate was added to the association mixture at 24 degrees C, biphasic dissociation rates of k = 0.23 min-1 and k = 0.02 min-1 were observed, suggesting that the chaotropic salt blocks the interconversion of rates. The data are consistent with the model for cAMP-dependent protein kinase which exhibits two binding sites with different affinities. The type II enzyme from mammary gland cytosol exhibits in addition the phenomenon of temperature-dependent interconversion of the two binding affinities.     

The hybridization of glyceraldehyde 3-phosphate dehydrogenases from rabbit muscle and yeast. Kinetics and thermodynamics of the reaction and isolation of the hybrid

A kinetic and thermodynamic study was made of the formation of the hybrid (R(2)Y(2)) glyceraldehyde 3-phosphate dehydrogenase from the yeast (Y(4)) and rabbit (R(4)) enzymes. The values of the thermodynamic parameters for the equilibrium between R(4), Y(4) and R(2)Y(2) suggest that the R(2)-R(2) and Y(2)-Y(2) interactions are similar. However, the failure to observe the RY(3) and R(3)Y hybrids is interpreted in terms of differences at the interfaces of the R-R and Y-Y interactions (the glyceraldehyde 3-phosphate dehydrogenase molecule being regarded as a dimer of dimers). The kinetics of formation of the R(2)Y(2) hybrid were studied and a model was proposed to account for the results. Best-fit values for the rate constants of the individual steps were evaluated by computer simulation, and the rate-limiting steps were identified as the dissociation of tetramers to dimers. It is proposed that the cleavage plane for dissociation of the tetramers corresponds to the region of low electron density through the centre of the molecule in the X-ray-crystallographic structure for human glyceraldehyde 3-phosphate dehydrogenase (Watson et al., 1972), which is probably the plane containing the Q and R axes in the lobster enzyme (Buehner et al., 1974). The R(2)Y(2) hybrid was isolated in milligram amounts by ion-exchange chromatography and its rate of reversion to the native enzyme was shown to be consistent with the kinetic model proposed from the hybrid-formation experiments.     

Binding of radioiodinated SPECT ligands to transfected cell membranes expressing single muscarinic receptor subtypes

The equilibrium dissociation constant and the kinetic rate constants were determined for the binding of (R)-[3H]3-quinuclidinyl benzilate ([3H]QNB) and [125I]3-quinuclidinyl-4-iodobenzilate ((R,R)- and (R,S)-[125I]IQNB) to transfected cell membranes expressing one single muscarinic acetylcholine receptor (mAChR) subtype. The association and dissociation kinetics for the m2 subtype were more rapid than for the m1 and m3 subtypes. The differential kinetic properties may be useful for the single photon emission computed tomographic (SPECT) evaluation of regional mAChR subtype alterations in disease states.     

Stopped-flow kinetic analysis of replication protein A-binding DNA: damage recognition and affinity for single-stranded DNA reveal differential contributions of k(on) and k(off) rate constants

Replication protein A (RPA) is a heterotrimeric protein required for many DNA metabolic functions, including replication, recombination, and nucleotide excision repair (NER). We report the pre-steady-state kinetic analysis of RPA-binding DNA substrates using a stopped-flow assay to elucidate the kinetics of DNA damage recognition. The bimolecular association rate, k(on), for RPA binding to duplex DNA substrates is greatest for a 1,3d(GXG), intermediate for a 1,2d(GpG) cisplatin-DNA adduct, and least for an undamaged duplex DNA substrate. RPA displays a decreased k(on) and an increased k(off) for a single-stranded DNA substrate containing a single 1,2d(GpG) cisplatin-DNA adduct compared with an undamaged DNA substrate. The k(on) for RPA-binding single-stranded polypyrimidine sequences appears to be diffusion-limited. There is minimal difference in k(on) for varying length DNA substrates; therefore, the difference in equilibrium binding affinity is mainly attributed to the k(off). The k(on) for a purine-rich 30-base DNA is reduced by a factor of 10 compared with a pyrimidine-rich DNA of identical length. These results provide insight into the mechanism of RPA-DNA binding and are consistent with RPA recognition of DNA-damage playing a critical role in NER.     

Stopped-flow kinetic analysis of eIF4E and phosphorylated eIF4E binding to cap analogs and capped oligoribonucleotides: evidence for a one-step binding mechanism

Recruitment of eukaryotic mRNA to the 48 S initiation complex is rate-limiting for protein synthesis under normal conditions. Binding of the 5' -terminal cap structure of mRNA to eIF4E is a critical event during this process. Mammalian eIF4E is phosphorylated at Ser-209 by Mnk1 and Mnk2 kinases. We investigated the interaction of both eIF4E and phosphorylated eIF4E (eIF4E(P)) with cap analogs and capped oligoribonucleotides by stopped-flow kinetics. For m(7)GpppG, the rate constant of association, k(on), was dependent on ionic strength, decreasing progressively up to 350 mm KCl, but the rate constant of dissociation, k(off), was independent of ionic strength. Phosphorylation of eIF4E decreased k(on) by 2.1-2.3-fold at 50-100 mm KCl but had progressively less effect at higher ionic strengths, being negligible at 350 mm. Contrary to published evidence, eIF4E phosphorylation had no effect on k(off). Several observations supported a simple one-step binding mechanism, in contrast to published reports of a two-step mechanism. The kinetic function that best fit the data changed from single- to double-exponential as the eIF4E concentration was increased. However, measuring k(off) for dissociation of a pre-formed eIF4E.m(7)GpppG complex suggested that the double-exponential kinetics were caused by dissociation of eIF4E dimers, not a two-step mechanism. Addition of a 12-nucleotide chain to the cap structure increased affinity at high ionic strength for both eIF4E (24-fold) and eIF4E(P) (7-fold), primarily due to a decrease in k(off). This suggests that additional stabilizing interactions between capped oligoribonucleotides and eIF4E, which do not occur with cap analogs alone, act to slow dissociation.     

Kinetic studies on the interaction of eglin c with human leukocyte elastase and cathepsin G

We have investigated the inhibition of human leukocyte elastase and cathepsin G by recombinant Eglin c under near physiological conditions. The association rate constants k on of Eglin c for elastase and cathepsin G were 1.3 X 10(7) M-1 s-1 and 2 X 10(6) M-1 s-1, respectively. Under identical conditions, the k on for the association of human plasma alpha 1-proteinase inhibitor with the two leukocproteinases were 2.4 X 10(7) M-1 s-1 and 10(6) M-1 s-1, respectively. The consistency of these data could be verified using a set of competition experiments. The elastase-Eglin c interaction was studied in greater detail. The dissociation rate constant k off was determined by trapping of free elastase from an equilibrium mixture of elastase and Eglin c with alpha 1-proteinase inhibitor or alpha 2-macroglobulin. The rate of dissociation was very low (k off = 3.5 X 10(-5) s-1). The calculated equilibrium dissociation constant of the complex, Ki(calc) = k off/k on, was found to be 2.7 X 10(-12) M. Ki was also measured by adding elastase to mixtures of Eglin c and substrate and determining the steady-state rates of substrate hydrolysis. The Ki determined from these experiments (7.5 X 10(-11) M) was significantly higher than Ki(calc). This discrepancy might be explained by assuming that the interaction of Eglin c with elastase involves two steps: a fast binding reaction followed by a slow isomerization step. From the above kinetic constants it may be inferred that at a therapeutic concentration of 5 X 10(-7) M, Eglin c will inhibit leukocyte elastase in one second and will bind this enzyme in a "pseudo-irreversible" manner.     

Construction of a synthetic gene for an R-plasmid-encoded dihydrofolate reductase and studies on the role of the N-terminus in the protein

R67 dihydrofolate reductase (DHFR) is a novel protein that provides clinical resistance to the antibacterial drug trimethoprim. The crystal structure of a dimeric form of R67 DHFR indicates the first 16 amino acids are disordered [Matthews et al. (1986) Biochemistry 25, 4194-4204]. To investigate whether these amino acids are necessary for protein function, the first 16 N-terminal residues have been cleaved off by chymotrypsin. The truncated protein is fully active with kcat = 1.3 s-1, Km(NADPH) = 3.0 microM, and Km(dihydrofolate) = 5.8 microM. This result suggests the functional core of the protein resides in the beta-barrel structure defined by residues 27-78. To study this protein further, synthetic genes coding for full-length and truncated R67 DHFRs were constructed. Surprisingly, the gene coding for truncated R67 DHFR does not produce protein in vivo or confer trimethoprim resistance upon Escherichia coli. Therefore, the relative stabilities of native and truncated R67 DHFR were investigated by equilibrium unfolding studies. Unfolding of dimeric native R67 DHFR is protein concentration dependent and can be described by a two-state model involving native dimer and unfolded monomer. Using absorbance, fluorescence, and circular dichroism techniques, an average delta GH2O of 13.9 kcal mol-1 is found for native R67 DHFR. In contrast, an average delta GH2O of 11.3 kcal mol-1 is observed for truncated R67 DHFR. These results indicate native R67 DHFR is 2.6 kcal mol-1 more stable than truncated protein. This stability difference may be part of the reason why protein from the truncated gene is not found in vivo in E. coli.     

Multi-method global analysis of thermodynamics and kinetics in reconstitution of monellin

Accurate and precise determinations of thermodynamic parameters of binding are important steps toward understanding many biological mechanisms. Here, a multi-method approach to binding analysis is applied and a detailed error analysis is introduced. Using this approach, the binding thermodynamics and kinetics of the reconstitution of the protein monellin have been quantitatively determined in detail by simultaneous analysis of data collected with fluorescence spectroscopy, surface plasmon resonance and isothermal titration calorimetry at 25 degrees C, pH 7.0 and 150 mM NaCl. Monellin is an intensely sweet protein composed of two peptide chains that form a single globular domain. The kinetics of the reconstitution reaction are slow, with an association rate constant, k(on) of 8.8 x 10(3) M(-1) s(-1) and a dissociation rate constant, k(off) of 3.1 x 10(-4) s(-1). The equilibrium constant K(A) is 2.8 x 10(7) M(-1) corresponding to a standard free energy of association, DeltaG degrees , of -42.5 kJ/mol. The enthalpic component, DeltaH degrees , is -18.7 kJ/mol and the entropic contribution, DeltaS degrees , is 79.8 J mol(-1) K(-1) (-TDeltaS degrees = -23.8 kJ/mol). The association of monellin is therefore a bimolecular intra-protein association whose energetics are slightly dominated by entropic factors.     

Kinetic and thermodynamic characterization of HIV-1 protease inhibitors

Interaction kinetic and thermodynamic analyses provide information beyond that obtained in general inhibition studies, and may contribute to the design of improved inhibitors and increased understanding of molecular interactions. Thus, a biosensor-based method was used to characterize the interactions between HIV-1 protease and seven inhibitors, revealing distinguishing kinetic and thermodynamic characteristics for the inhibitors. Lopinavir had fast association and the highest affinity of the tested compounds, and the interaction kinetics were less temperature-dependent as compared with the other inhibitors. Amprenavir, indinavir and ritonavir showed non-linear temperature dependencies of the kinetics. The free energy, enthalpy and entropy (DeltaG, DeltaH, DeltaS) were determined, and the energetics of complex association (DeltaG(on), DeltaH(on), DeltaS(on)) and dissociation (DeltaG(off), DeltaH(off), DeltaS(off)) were resolved. In general, the energetics for the studied inhibitors was in the same range, with the negative free energy change (DeltaG < 0) due primarily to increased entropy (DeltaS > 0). Thus, the driving force of the interaction was increased degrees of freedom in the system (entropy) rather than the formation of bonds between the enzyme and inhibitor (enthalpy). Although the DeltaG(on) and DeltaG(off) were in the same range for all inhibitors, the enthalpy and entropy terms contributed differently to association and dissociation, distinguishing these phases energetically. Dissociation was accompanied by positive enthalpy (DeltaH(off) > 0) and negative entropy (DeltaS(off) < 0) changes, whereas association for all inhibitors except lopinavir had positive entropy changes (DeltaS(on) > 0), demonstrating unique energetic characteristics for lopinavir. This study indicates that this type of data will be useful for the characterization of target-ligand interactions and the development of new inhibitors of HIV-1 protease.