Structure-based pKa prediction provides a thermodynamic basis for the role of histidines in pH-induced conformational transitions in dengue virus

pH-induced conformational changes in dengue virus (DENV) are critical to its ability to infect host cells. The envelope protein heterodimers that make up the viral envelope shift from a dimer to a trimer conformation at low-pH during membrane fusion. Previous studies have suggested that the ionization of histidine residues at low-pH is central to this pH-induced conformational change. We sought out to use molecular modeling with structure-based pKa prediction to provide a quantitative basis for the role of histidines in pH-induced conformational changes and identify which histidine residues were primarily responsible for this transition. We combined existing crystallographic and cryo-electron microscopy data to construct templates of the dimer and trimer conformations for the mature and immature virus. We then generated homology models for the four DENV serotypes and carried out structure-based pKa prediction using Rosetta. Our results showed that the pKa values of a subset of conserved histidines in DENV successfully capture the thermodynamics necessary to drive pH-induced conformational changes during fusion. Here, we identified the structural determinants underlying these pKa values and compare our findings with previous experimental results.     

Calculation of the stability of mini-duplexes of DNA in an electrolyte solution using a molecular mechanics method in approximating a statistical model of the environment

Approximation of the statistical model of environment (SME) to estimate the energy of the macromolecule in electrolyte solution has been developed and used for calculating the conformational energy of nucleic acids by means of the molecular mechanics method. Calculation of base pairs opening delta Hcalop enthalpies and enthalpies of DNA miniduplexes dissociation delta Hcaldis were performed for 10 types of diduplexes. The approximation SME enable to perform calculations of the absolute base-dependent values of delta Hcalop which coincide in the range of 1 kcal/mol with the experimental base-dependent values of the helix-coil transition enthalpies. Values of dissociation enthalpies delta Hcalop greater than delta Hcaldis for all miniduplexes, the difference of delta Hcalop--delta Hcaldis determine the base-dependent energy of the helix-coil boundary. The values of activation barriers for the strands dissociation delta H d not equal to congruent to 8 kcal/mol and association delta H not equal to as congruent to 4 kcal/mol were obtained for GG/CC and AA/TT duplexes. It is concluded that the approximation SME enables to increase substantially the accuracy of the calculations of the macromolecule conformational rearrangement enthalpies in the electrolytes solution.     

Interaction of Trimeresurus flavoviridis phospholipase A2 and its fragment with calcium ion

Dimeric T. flavoviridis phospholipase A2 has been studied in terms of the interaction with essential Ca2+ by equilibrium gel filtration, ultraviolet difference spectroscopy, fluorescence measurements, and chemical modifications with p-bromophenacyl bromide. The subunit bound to Ca2+ with a 1:1 molar ratio and no cooperative binding was observed. The hypochromic effect produced upon the binding of Ca2+ is due to perturbation of (a) specific tryptophan residue(s) located in the vicinity of the active site and appears to be characteristic of this enzyme. On the basis of the pH dependence of the dissociation constants, it has been found that the alpha-amino group (pKa 8.7) controls the binding of Ca2+. Deprotonation of the alpha-amino group is possibly accompanied by conformational transition to the active form which is able to bind Ca2+. This is in contrast to the case of bovine pancreatic phospholipase A2 in which Asp-49 (pKa 5.2) is responsible for the metal ion binding (Fleer et al. (1981) Eur. J. Biochem. 113, 283-288). Des-octapeptide(1-8)-phospholipase A2 (L-fragment) was found to be capable of binding Ca2+ under the control of a group with a pKa of 7.6. This pKa value was similar to an apparent pKa of 7.5 determined for the histidine residue in the active site of the native enzyme by way of p-bromophenacyl bromide modification. It appears that the N-terminal (octapeptide) sequence affects the binding mode of Ca2+, possibly because of conformational transition arising from its removal. The reinvestigation showed that the N-terminal octapeptide sequence is Gly-Leu-Trp-Gln-Phe-Glu-Asn-Met.     

The anomalous pKa of Tyr-9 in glutathione S-transferase A1-1 catalyzes product release

The pKa of the catalytic Tyr-9 in glutathione S-transferase (GST) A1-1 is lowered from 10.3 to approximately 8.1 in the apoenzyme and approximately 9.0 with a GSH conjugate bound at the active site. However, a clear functional role for the unusual Tyr-9 pKa has not been elucidated. GSTA1-1 also includes a dynamic C terminus that undergoes a ligand-dependent disorder-to-order transition. Previous studies suggest a functional link between Tyr-9 ionization and C-terminal dynamics. Here we directly probe the role of Tyr-9 ionization in ligand binding and C-terminal conformation. An engineered mutant of rGSTA1-1, W21F/F222W, which contains a single Trp at the C terminus, was used as a fluorescent reporter of pH-dependent C-terminal dynamics. This mutant exhibited a pH-dependent change in Trp-222 emission properties consistent with changes in C-terminal solvation or conformation. The apparent pKa values for the conformational transition were 7.9 +/- 0.1 and 9.3 +/- 0.1 for the apoenzyme and ligand-bound enzyme, respectively, in excellent agreement with the pKa for Tyr-9 in these states. The Y9F/W21F/F222W mutant, however, exhibited no such pH-dependent changes. Time-resolved fluorescence anisotropy studies revealed a ligand-dependent, Tyr-9-dependent, change in the order parameter of Trp-222. However, no pH dependence was observed. In equilibrium and pre-steady-state ligand binding studies, product conjugate had a decreased equilibrium binding affinity (KD), concomitant with increased binding and dissociation rates, at higher pH values. Furthermore, the recovered pKa values for the pH-dependent microscopic rate constants ranged from 7.7 to 8.4, also in agreement with the pKa of Tyr-9. In contrast, the Y9F/W21F/F222W mutant had no pH-dependent transition in KD or rate constants for ligand binding or dissociation. The combined results indicate that the macroscopic populations of "open" and "closed" states of the C terminus are not determined solely by the ionization state of Tyr-9. However, the rates of transition between these states are faster for the ionized Tyr-9. The ionized Tyr-9 states provide a parallel pathway for product dissociation, which is kinetically and thermodynamically favored. In silico kinetic models further support the functional role for the parallel dissociation pathway provided by ionized Tyr-9.     

Conformational stability of dimeric proteins: quantitative studies by equilibrium denaturation.

The conformational stability of dimeric globular proteins can be measured by equilibrium denaturation studies in solvents such as guanidine hydrochloride or urea. Many dimeric proteins denature with a 2-state equilibrium transition, whereas others have stable intermediates in the process. For those proteins showing a single transition of native dimer to denatured monomer, the conformational stabilities, delta Gu (H2O), range from 10 to 27 kcal/mol, which is significantly greater than the conformational stability found for monomeric proteins. The relative contribution of quaternary interactions to the overall stability of the dimer can be estimated by comparing delta Gu (H2O) from equilibrium denaturation studies to the free energy associated with simple dissociation in the absence of denaturant. In many cases the large stabilization energy of dimers is primarily due to the intersubunit interactions and thus gives a rationale for the formation of oligomers. The magnitude of the conformational stability is related to the size of the polypeptide in the subunit and depends upon the type of structure in the subunit interface. The practical use, interpretation, and utility of estimation of conformational stability of dimers by equilibrium denaturation methods are discussed.     

Effects of conformational selectivity and of overlapping kinetically influential ionizations on the characteristics of pH-dependent enzyme kinetics. Implications of free-enzyme pKa variability in reactions of papain for its catalytic mechanism.

The effects of selection by a small molecule, when binding to a protein, of a particular conformation from an equilibrium stereopopulation on the characteristics of the pH-dependence of reaction with a reactivity probe or substrate were determined by analysis of an appropriate kinetic model. For reaction in one protonic state containing an equilibrium mixture of two conformational isomers, the pH-second-order rate constant (k) profile is of conventional sigmoidal form. The apparent pKa value is a composite of the pKa values of the two conformational states. The value of pKapp. for a given enzyme under given experimental conditions will always be the same (provided that the site-specificity assumed in the model is maintained) irrespective of whether only one conformation reacts or both react, with the same or with different rate constants. The experimentally determined pH-independent rate constant (kapp.) is an average of the reactivities of the two conformational states weighted in favour of the predominant form. The presence of an additional but unreactive conformational state also affects the value of kapp. The possibility that overlapping acid dissociations that affect the reactivity of the enzyme might provide pH-k profiles often indistinguishable in practice from simple sigmoidal dissociation curves and subject to variability in apparent pKa values was evaluated by a simulation study. If two reactive protonic states of the enzyme respond differently to changes in the structure of the substrate or site-specific reactivity probe, differences in apparent pKa values of up to approx. 1 unit can be exhibited without deviation from sigmoidal behaviour being reliably observed. Differences in apparent pKa values observed in some site-specific reactions of papain and their possible consequences for its catalytic mechanism are discussed.     

Calorimetric determination of linkage effects involving an acyl-enzyme intermediate

Enthalpy changes of alpha-chymotrypsin acylation by 3-(2-furyl)acryloylimidazole (FAI) were calorimetrically determined as a function of pH. By observing the functional dependence of acylation enthalpies on buffer ionization heats, a complex pH profile was obtained describing proton release accompanying formation of acyl-enzyme. A pKa of 4.0 for FAI ionization and apparent pKa values of 6.8, 7.55 and 8.8 on the enzyme were used to account for the proton release data. A model which accounts for the proton release behavior was used to fit the acylation enthalpy data and values for the apparent dissociation enthalpies of the groups involved were obtained along with a pH-independent intrinsic enthalpy of acylation. This model suggests a group with an apparent pK = 6.8 and delta Hion = 8.7 kcal/mol which is perturbed to a pK of 7.55 and delta Hion = 7.6 kcal/mol on attachment of the acyl moiety to the enzyme. The apparent ionization enthalpy change for the active-inactive transition (pK3 = 8.8; delta H = 3.0 kcal/mol) corresponds with that calculated from the data of Fersht (J. Mol. Biol. 64 (1972) 497). The pH-independent intrinsic enthalpy of acylation (delta H = -7.9 kcal/mol) is corrected for group ionizations linked to the acylation process. Consequently, it more closely reflects molecular processes of interest such as substrate binding, covalent bond rearrangement, and product release.     

Free energy of transition for the individual alkaline conformers of yeast iso-1-cytochrome c

Direct protein electrochemistry was used to obtain the thermodynamic parameters of transition from the native (state III) to the alkaline (state IV) conformer for untrimethylated Saccharomyces cerevisiae iso-1-cytochrome c expressed in E. coli and its single and multiple lysine-depleted variants. In these variants, one or more of the lysine residues involved in axial Met substitution (Lys72, Lys73, and Lys79) was mutated to alanine. The aim of this work is to determine the thermodynamic affinity of each of the substituting lysines for the heme iron and evaluate the interplay of enthalpic and entropic factors. The equilibrium constants for the deprotonation reaction of Lys72, 73, and 79 were computed for the minimized MD average structures of the wild-type and mutated proteins, applying a modified Tanford-Kirkwood calculation. Solvent accessibility calculations for the substituting lysines in all variants were also performed. The transition enthalpy and entropy values within the protein series show a compensatory behavior, typical of a process involving extensive solvent reorganization effects. The experimental and theoretical data indicate that Lys72 most readily deprotonates and replaces M80 as the axial heme iron ligand, whereas Lys73 and Lys79 show comparably higher pKa values and larger transition free energies. A good correlation is found within the series between the lowest calculated Lys pKa value and the corresponding experimental pKa value, which can be interpreted as indicative of the deprotonating lysine itself acting as the triggering group for the conformational transition. The triple Lys to Ala mutant, in which no lysine residues are available for heme iron binding, features transition thermodynamics consistent with a hydroxide ion replacing the axial methionine ligand.     

Hydroxyl-ion-induced subunit dissociation of east cytoplasmic pyruvate decarboxylase. A circular dichroism study

Cytoplasmic pyruvate decarboxylase (EC 4.1.1.1, from Saccharomyces carlsbergensis) exhibits in its circular dichroic spectrum in the 250--320-nm range a multiple two-signal band. This couplet disappears on increasing the pH up to pH 8.5. Two classes of two protons each can be quantified by these spectral changes. The first class dissociates rapidly and the apparent pK is 7.84. The thermodynamic data are delta G = 87.7 kJ mol-1, delta H = + 56.0 kJ mol-1, delta S = - 108 J mol-1 K-1, very characteristic for the deprotonation of an amino-acid side chain. The second class of the protons has the following thermodynamic data: delta G = 88.3 kJ mol-1, delta H = - 64.3 kJ mol-1, delta S = - 520 J mol-1 K-1 which, in conjunction with kinetic reasoning and in view of enzyme stoichiometry and symmetry, suggests a conformational equilibrium exposing the second two protons. Th enzyme dissociates into two dimeric subunits. This dissociation step is considered to be rate-determining for the overall process. The data are kp = 1.4 . 10(-3), delta H not equal to = + 128.3 kJ mol-1, delta S not equal = + 136 J mol-1 K-1. If there is a conformational equilibrium, the rate constant of product formation kp will be modified by a factor beta = kc/(1 + Kc) (0 < beta less than or equal to 1) where Kc is the conformational equilibrium constant. The subunit dissociation appears to be controlled by the enthalpy of activation indicating that a number of interactions, i.e. ionic, hydrogen and hydrophobic bridges, are to be broken. Optimal conditions for the preparation of the apo-enzyme are derived from the data.     

Ionizable groups linked to the reaction of 2,2'-dithiobispyridine with hemoglobin.

The pH-dependence of the second-order rate-constant for the reaction of 2,2'-dithiobispyridine with the CysF9(93) beta sulphydryl group of hemoglobin in the R quaternary structure is analyzed in terms of a tentative model based on the observation that this sulphydryl exists as a mixture of two tertiary conformations in dynamic equilibrium. For the four aquomethemoglobins studied (human A and S, dog and rabbit), the equation derived from this model gives a better fit than a simpler equation based on the assumption of only one tertiary conformation. For the corresponding carbonmonoxyhemoglobins the simpler equation gives a better fit. The dog and rabbit oxy and azidomet data are better fitted by the model equation, whereas the data for the corresponding human A and S derivatives are better fitted by the simpler equation. From the analysis pKa values of 6.1 and 8.7 are obtained for the ionization of groups coupled to the presumed conformational transition. The pKa of 6.1 is assigned to HisHC3(146) beta; the pKa of 8.7 is assigned to the CysF9(93) beta sulphydryl group in its external conformation. It is estimated that the pKa of this sulphydryl may be as high as 12.9 in its internal conformation.     

Stability of recombinant Lys25-ribonuclease T1

The conformational stability of recombinant Lys25-ribonuclease T1 has been determined by differential scanning microcalorimetry (DSC), UV-monitored thermal denaturation measurements, and isothermal Gdn.HCl unfolding studies. Although rather different extrapolation procedures are involved in calculating the Gibbs free energy of stabilization, there is fair agreement between the delta G degrees values derived from the three different experimental techniques at pH 5, theta = 25 degrees C: DSC, 46.6 +/- 2.1 kJ/mol; UV melting curves, 48.7 +/- 5 kJ/mol; Gdn.HCl transition curves, 40.8 +/- 1.5 kJ/mol. Thermal unfolding of the enzyme is a reversible process, and the ratio of the van't Hoff and calorimetric enthalpy, delta HvH/delta Hcal, is 0.97 +/- 0.06. This result strongly suggests that the unfolding equilibrium of Lys25-ribonuclease T1 is adequately described by a simple two-state model. Upon unfolding the heat capacity increases by delta Cp degrees = 5.1 +/- 0.5 kJ/(mol.K). Similar values have been found for the unfolding of other small proteins. Surprisingly, this denaturational heat capacity change practically vanishes in the presence of moderate NaCl concentrations. The molecular origin of this effect is not clear; it is not observed to the same extent in the unfolding of bovine pancreatic ribonuclease A, which was employed in control experiments. NaCl stabilizes Lys25-ribonuclease T1. The transition temperature varies with NaCl activity in a manner that suggests two limiting binding equilibria to be operative. Below approximately 0.2 M NaCl activity unfolding is associated with dissociation of about one ion, whereas above that concentration about four ions are released in the unfolding reaction.(ABSTRACT TRUNCATED AT 250 WORDS)     

Crystal structure analysis of oxidized Pseudomonas aeruginosa azurin at pH 5.5 and pH 9.0. A pH-induced conformational transition involves a peptide bond flip

The X-ray crystal structure of recombinant wild-type azurin from Pseudomonas aeruginosa was determined by difference Fourier techniques using phases derived from the structure of the mutant His35Leu. Two data sets were collected from a single crystal of oxidized azurin soaked in mother liquor buffered at pH 5.5 and pH 9.0, respectively. Both data sets extend to 1.93 A resolution. The two pH forms were refined independently to crystallographic R-factors of 17.6% (pH 5.5) and 17.5% (pH 9.0). The conformational transition previously attributed to the protonation/deprotonation of residue His35 (pKa(red) = 7.3, pKa(ox) = 6.2), which lies in a crevice of the protein close to the copper binding site, involves a concomitant Pro36-Gly37 main-chain peptide bond flip. At the lower pH, the protonated imidazole N delta 1 of His35 forms a strong hydrogen bond with the carbonyl oxygen from Pro36, while at alkaline pH the deprotonated N delta 1 acts as an acceptor of a weak hydrogen bond from HN Gly37. The structure of the remainder of the azurin molecule, including the copper binding site, is not significantly affected by this transition.     

Electrostatic interaction between anions bound to site I and the retinal Schiff base of halorhodopsin

The influence of different anions on the deprotonation of the retinal Schiff base of halorhodopsin in the dark was investigated. We find that a large number of anions cause a significant increase of the pKa of the Schiff base, an effect attributed to binding to "site I" on the protein. The concentration dependencies of the spectroscopic shifts associated with the changes of the pKa yielded dissociation constants (and thus binding energies) for the anions, which were related to the Stokes radii. The data fit the predictions of electrostatic interaction between the anions and the positive charge associated with site I, if the latter is located within a few angstroms from the surface of the protein. The specificity of site I toward various anions is quantitatively explained by the differences in the change of Born energy upon transfer of the anions from water to the binding site. The changes in the deprotonation energy of the Schiff base upon the binding of anions, delta delta Gdeprot, could be calculated from the delta pKa at infinite anion concentration. Unexpectedly, the delta delta Gdeprot values were remarkably close to the energies of binding to site I. Thus, site I and the Schiff base are strongly electrostatically coupled, either because of close proximity or because of the possibility of allosteric energy transfer between them.     

Kinesin undergoes a 9 S to 6 S conformational transition.

Addition of NaCl or KCl in the presence of 50 nM ATP induces a shift in the sedimentation coefficient (apparent S20,w) of kinesin from 9.4 S at low ionic strength to 6.5 S at high ionic strength. The midpoint for the transition occurs at ionic strength values of 0.39, 0.25, and 0.18 for pH values of 6.3, 6.9, and 8.3, respectively. Gel filtration experiments indicate that the transition to the 6.5 S species is accompanied by a decrease in the diffusion coefficient. Under all conditions which were tested, the 64-kDa beta subunits comigrate with the 120-kDa alpha subunits without any evidence for dissociation of the alpha 2 beta 2 complex. These results are consistent with the change in sedimentation coefficient being due to a conformational transition between a folded form at low ionic strength and an extended form at high ionic strength. This conformational transition is not significantly affected by the nature of the nucleotide bound at the active site since similar results are obtained both in the presence of excess EDTA, which removes the bound ADP, and after replacement of the bound ADP with adenosine 5'-(beta,gamma-imino)triphosphate. The alpha 2 form of kinesin, which lacks the beta subunits, undergoes a similar transition between a 6.7 S form at low ionic strength and a 5.1 S form at high ionic strength with a midpoint for the transition at an ionic strength of 0.5 at pH 6.9. Electron microscopic observation also indicates a transition between a folded conformation at low ionic strength and an extended conformation at high ionic strength for both the alpha 2 beta 2 and alpha 2 species.     

Mechanism and dynamics of conformational ordering in xanthan polysaccharide

The thermally induced order-disorder transition of xanthan (extracellular bacterial polysaccharide from Xanthomonas campestris) has been investigated by optical rotation, differential scanning calorimetry, stopped-flow reaction kinetics and low-angle laser light scattering, and the results have been analysed in terms of Zimm -Bragg helix-coil transition theory. The reciprocal of the transition midpoint temperature (Tm) varies linearly with the logarithm of cation (K+) the salt dependence of Tm, is in agreement with Manning polyelectrolyte theory the ordered structure. The associated increase in cation binding, calculated from the salt dependence of tm, is in agreement with the Manning polyelectrolyte theory for one of the candidate structures from X-ray diffraction, a 5(1) single helix stabilized by packing of side-chains along the polymer backbone, but not for the alternative double-helix structure that has also been proposed. At each salt concentration, the two fundamental parameters of the Zimm -Bragg theory, s and sigma, were calculated. The equilibrium constant for growth of the ordered structure (s) is derived directly from calorimetric measurement of transition enthalpy (delta Hcal ), and sigma, which quantifies the relative instability of the helix nucleus, is derived from the ratio of delta Hcal to the apparent transition enthalpy (delta Happ ) obtained by van't Hoff analysis of the optical rotation data. The temperature course of conformational ordering calculated theoretically is in good quantitative agreement with experimental results from both optical rotation and scanning calorimetry. The calculated average length of stable, ordered chain-sequences increases with decreasing temperature, but equals or exceeds the total chain length from light scattering only at temperatures more than approximately equal to 70 K below Tm, suggesting that ordered and disordered regions may co-exist within the same xanthan molecule. Consistent with this interpretation, the observed rate of conformational ordering increases sharply under conditions where the starting solution for dynamic measurements is partially ordered, suggesting that ordered sequences within each chain may act as helix nuclei for adjacent disordered regions, so that helix growth, rather than the slower nucleation process, becomes rate limiting.     

Conformational studies of nucleic acids: IV. The conformational energetics of oligonucleotides: d(ApApApA) and ApApApA

Utilizing a new method for modeling furanose pseudorotation (D. A. Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A. Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA. The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: delta (C5'-C4'-C3'-O3')-chi (O4'-C1'-N9-C4), where the energetic basis for an approximately linear delta-chi correlation can be seen: zeta (C3'-O3'-P-O5')-delta, in which the experimentally observed linear correlation between zeta and delta in DNA(220 degrees less than zeta less than 280 degrees) is clearly predicted; zeta-epsilon (C4'-C3'-O3'-P), which shows that epsilon increases with decreasing zeta less than 260 degrees; alpha (O3'-P-O5'-C5')-gamma (O5'-C5'-C4'-C3') where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion delta is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in alpha and gamma. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.     

Energetics of superhelicity and of B-Z transitions in superhelical DNA

The linking difference, alpha, imposed upon a superhelically constrained DNA molecule must be partitioned between twisting and bending deformations. Transitions to alternative secondary structures can occur at susceptible sites, altering the local molecular twist by an amount delta Twtrans. That part of the linking difference not accommodated in this way, the residual linking difference alpha res, must be manifested as smooth torsional and flexural deformations of secondary structure. The competition among the alternative ways of accommodating the imposed linking difference alpha determines a stressed equilibrium state. The superhelical free energy, G(alpha), is the excess free energy of the equilibrium state at linking difference alpha above that of the relaxed state under identical conditions. In this paper a method is described by which the free energies associated both to linking, G(alpha), and to residual linking differences can be determined from data on superhelical conformational transitions. The application of this approach to previously published experimental data on the B-Z transition suggests that the free energy associated with alpha res is about 30% larger at substantial superhelicities than it is near the relaxed state. At the onset of transition the functional form of G(alpha) is shown to change in a manner dependent upon the length of the Z-susceptible site.     

Effect of methylamine and plasmin on the conformation of human alpha 2-macroglobulin as revealed by differential scanning calorimetric analysis.

Differential scanning calorimetric analysis was used as a probe of the conformational alteration in human alpha 2-macroglobulin (AM) upon its complex formation with methylamine and with the protease, human plasmin. The slow electrophoretic form of AM displayed a single thermal transition, characterized by a temperature midpoint (Tm) of 65.8 +/- 0.3 degrees, a calorimetric enthalpy (delta Hc) of 2,550 +/- 150 kcal/mol and a van't Hoff enthalpy (delta Hvh) of 140 kcal/mol. In the presence of sufficient methylamine to irreversibly disrupt the four thiol ester bonds in AM, a single thermal transition was obtained, characterized by a Tm of 62.8 +/- 0.3 degrees, a delta Hc of 1,700 +/- 100 kcal/mol, and a delta Hvh of 169 kcal/mol. These data suggest that a major conformational alteration is produced in AM upon complex formation with methylamine. When plasmin interacts with AM, the resulting thermogram displays Tm values for AM of 68-69 degrees and 77 degrees, also suggestive of a large conformational alteration in AM. However, this latter alteration appears dissimilar to the change induced by methylamine.     

Two steps in the transition between the native and acid states of bovine α-lactalbumin detected by circular polarization of luminescence: Evidence for a premolten globule state?

A few studies indirectly support the existence of an intermediate in the transition of Ca(2+)-saturated bovine alpha-lactalbumin (alpha-LA) from the native (N) to the acidic (A) state, known as the molten globule state. However, direct experimental evidence for the appearance of this intermediate has not been obtained. The signal of circular polarization of luminescence (CPL) is sensitive to fine conformational transitions because of its susceptibility to changes in the environmental asymmetry of fluorescent chromophores in their excited electronic states. In the present study, CPL measurements were applied using the intrinsic tryptophan fluorescence of alpha-LA as well as the fluorescence of 8-anilino-1-naphthalenesulfonic acid (ANS) bound to alpha-LA. CPL of tryptophan and ANS was measured in the pH range of 2.5-6 in order to find direct experimental evidence for the proposed intermediate. CPL (characterized by the emission anisotropy factor, g(em)) depends on the asymmetry of the protein molecular structure in the environment of the tryptophan and the ANS chromophores in the excited electronic state. The pH dependence of both the gab, absorption anisotropy factor determined by CD, and the ANS steady state fluorescence, showed a single transition at pH 3-3.7 as already reported elsewhere. This transition was interpreted as being a result of a change of the alpha-LA tertiary structure, which resulted in a loss of asymmetry of the environment of both the tryptophan residues and the ANS hydrophobic binding sites. The pH dependence of the tryptophan and ANS g(em) showed an additional conformational transition at pH 4-5, which coincided with the pKa of Ca2+ dissociation (pKa 5), as predicted by Permyakov et al. (1981, Biochem Biophys Res Commun 100:191-197). The titration curve showed that there is a pH range between 3.7 and 4.1 in which alpha-LA exists in an intermediate state between the N- and A-state. We suggest that the intermediate is the premolten globule state characterized by a reduced Ca2+ binding to the alpha-LA, native-like tertiary structure, and reduced asymmetric fluctuation of the tertiary structure on the nanosecond time scale. This intermediate resembles the "critical activated state" theoretically deduced by Kuwajima et al. (1989, J Mol Biol 206:547-561). The present study demonstrates the power of CPL measurements for the investigation of folding/unfolding transitions in proteins.     

Effect of monomer composition on proton dissociation of weak polyacids

The study of the proton dissociation process of weak polyacids (eg α carboxylic poly(monoprotic)acid) is based on the knowledge of the change in electrostatic free energy, G(el), as a function of the variation of the number of charges on the polymer chain. The original treatment proposed by Manning can be used to describe the proton dissociation process of weak poly(monoprotic)acids, in the absence of pH-induced conformational transitions. In order to describe the alpha dependence of pKa of weak co-poly(monoprotic)acids containing two different acidic groups in different amounts along the polymer chain, a simple modification of the model is proposed. Abbreviations used: Due to the difficulty of using non-ASCI characters in the main text, we have used the following abbreviations: G(el) Gel 10-pH 10-pH pKa pKa alpha α alphan αn dG(el)/dalpha δGel/δα KA and KB KA and KB CA and CB CA and CB βA and βB βA and βB δpK δpK §max ξmax This revised version was published online in November 2006 with corrections to the Cover Date.