Desolvation Penalty for Burying Hydrogen-Bonded Peptide Groups in Protein Folding

Ordered-domain unfolding of thermophilic isolated β subunit ATP synthase

The flexibility of the ATP synthase's β subunit promotes its role in the ATP synthase rotational mechanism, but its domains stability remains unknown. A reversible thermal unfolding of the isolated β subunit (Tβ) of the ATP synthase from Bacillus thermophilus PS3, tracked through circular dichroism and molecular dynamics, indicated that Tβ shape transits from an ellipsoid to a molten globule through an ordered unfolding of its domains, preserving the β-sheet residual structure at high temperature. We determined that part of the stability origin of Tβ is due to a transversal hydrophobic array that crosses the β-barrel formed at the N-terminal domain and the Rossman fold of the nucleotide-binding domain (NBD), while the helix bundle of the C-terminal domain is the less stable due to the lack of hydrophobic residues, and thus the more flexible to trigger the rotational mechanism of the ATP synthase.

Physical driving force of actomyosin motility based on the hydration effect

We propose a driving force hypothesis based on previous thermodynamics, kinetics and structural data as well as additional experiments and calculations presented here on water-related phenomena in the actomyosin systems. Although Szent-Györgyi pointed out the importance of water in muscle contraction in 1951, few studies have focused on the water science of muscle because of the difficulty of analyzing hydration properties of the muscle proteins, actin, and myosin. The thermodynamics and energetics of muscle contraction are linked to the water-mediated regulation of protein-ligand and protein-protein interactions along with structural changes in protein molecules. In this study, we assume the following two points: (1) the periodic electric field distribution along an actin filament (F-actin) is unidirectionally modified upon binding of myosin subfragment 1 (M or myosin S1) with ADP and inorganic phosphate Pi (M.ADP.Pi complex) and (2) the solvation free energy of myosin S1 depends on the external electric field strength and the solvation free energy of myosin S1 in close proximity to F-actin can become the potential force to drive myosin S1 along F-actin. The first assumption is supported by integration of experimental reports. The second assumption is supported by model calculations utilizing molecular dynamics (MD) simulation to determine solvation free energies of a small organic molecule and two small proteins. MD simulations utilize the energy representation method (ER) and the roughly proportional relationship between the solvation free energy and the solvent-accessible surface area (SASA) of the protein. The estimated driving force acting on myosin S1 is as high as several piconewtons (pN), which is consistent with the experimentally observed force.

Thermodynamic characterization of the unfolding of the prion protein

The prion protein appears to be unusually susceptible to conformational change, and unlike nearly all other proteins, it can easily be made to convert to alternative misfolded conformations. To understand the basis of this structural plasticity, a detailed thermodynamic characterization of two variants of the mouse prion protein (moPrP), the full-length moPrP (23-231) and the structured C-terminal domain, moPrP (121-231), has been carried out. All thermodynamic parameters governing unfolding, including the changes in enthalpy, entropy, free energy, and heat capacity, were found to be identical for the two protein variants. The N-terminal domain remains unstructured and does not interact with the C-terminal domain in the full-length protein at pH 4. Moreover, the enthalpy and entropy of unfolding of moPrP (121-231) are similar in magnitude to values reported for other proteins of similar size. However, the protein has an unusually high native-state heat capacity, and consequently, the change in heat capacity upon unfolding is much lower than that expected for a protein of similar size. It appears, therefore, that the native state of the prion protein undergoes substantial fluctuations in enthalpy and hence, in structure.

Contribution of hydrogen bonds to protein stability

Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein.

Protein dynamics are influenced by the order of ligand binding to an antibiotic resistance enzyme

The aminoglycoside N3 acetyltransferase-IIIb (AAC) is responsible for conferring bacterial resistance to a variety of aminoglycoside antibiotics. Nuclear magnetic resonance spectroscopy and dynamic light scattering analyses revealed a surprising result; the dynamics of the ternary complex between AAC and its two ligands, an antibiotic and coenzyme A, are dependent upon the order in which the ligands are bound. Additionally, two structurally similar aminoglycosides, neomycin and paromomycin, induce strikingly different dynamic properties when they are in their ternary complexes. To the best of our knowledge, this is the first example of a system in which two identically productive pathways of forming a simple ternary complex yield significant differences in dynamic properties. These observations emphasize the importance of the sequence of events in achieving optimal protein-ligand interactions and demonstrate that even a minor difference in molecular structure can have a profound effect on biochemical processes.

A large solvent isotope effect on protein association thermodynamics

Solvent reorganization can contribute significantly to the energetics of protein-protein interactions. However, our knowledge of the magnitude of the energetic contribution is limited, in part, by a dearth of quantitative experimental measurements. The biotin repressor forms a homodimer as a prerequisite to DNA binding to repress transcription initiation. At 20 °C, the dimerization reaction, which is thermodynamically coupled to binding of a small ligand, bio-5'-AMP, is characterized by a Gibbs free energy of -7 kcal/mol. This modest net dimerization free energy reflects underlying, very large opposing enthalpic and entropic driving forces of 41 ± 3 and -48 ± 3 kcal/mol, respectively. The thermodynamics have been interpreted as indicating coupling of solvent release to dimerization. In this work, this interpretation has been investigated by measuring the effect of replacing H2O with D2O on the dimerization thermodynamics. Sedimentation equilibrium measurements performed at 20 °C reveal a solvent isotope effect of -1.5 kcal/mol on the Gibbs free energy of dimerization. Analysis of the temperature dependence of the reaction in D2O indicates enthalpic and entropic contributions of 28 and -37 kcal/mol, respectively, considerably smaller than the values measured in H2O. These large solvent isotope perturbations to the thermodynamics are consistent with a significant contribution of solvent release to the dimerization reaction.

Extended solvent-contact model for protein solvation: test cases for dipeptides

Solvation effects are critically important in the structural stabilization and functional optimization of proteins. Here, we propose a new solvation free energy function for proteins, and test its applicability in predicting the solvation free energies of dipeptides. The present solvation model involves the improvement of the previous solvent-contact model assuming that the molecular solvation free energy could be given by the sum over the individual atomic contributions. In addition to the existing solvent-contact term, the modified solvation free energy function includes the self-solvation term that reflects the effects of intramolecular interactions in the solute molecule on solute-solvent interactions. Four kinds of atomic parameters should be determined in this solvation model: atomic fragmental volume, maximum atomic occupancy, atomic solvation, and atomic self-solvation parameters. All of these parameters for 16 atom types are optimized with a standard genetic algorithm in such a way to minimize the difference between the solvation free energies of dipeptides obtained from high-level quantum chemical calculations and those predicted by the solvation free energy function. The solvation free energies of dipeptides estimated from the new solvation model are in good agreement with the quantum chemical results. Therefore, the optimized solvation free energy function is expected to be useful for examining the structural and energetic features of proteins in aqueous solution.

New solvation free energy function comprising intermolecular solvation and intramolecular self-solvation terms

Solvation free energy is a fundamental thermodynamic quantity that should be determined to estimate various physicochemical properties of a molecule and the desolvation cost for its binding to macromolecular receptors. Here, we propose a new solvation free energy function through the improvement of the solvent-contact model, and test its applicability in estimating the solvation free energies of organic molecules with varying sizes and shapes. This new solvation free energy function is constructed by combining the existing solute-solvent interaction term with the self-solvation term that reflects the effects of intramolecular interactions on solvation. Four kinds of atomic parameters should be determined in this solvation model: atomic fragmental volume, maximum atomic occupancy, atomic solvation, and atomic self-solvation parameters. All of these parameters for total 37 atom types are optimized by the operation of a standard genetic algorithm in such a way to minimize the difference between the experimental solvation free energies and those calculated by the solvation free energy function for 362 organic molecules. The solvation free energies estimated from the new solvation model compare well with the experimental results with the associated squared correlation coefficients of 0.88 and 0.85 for training and test sets, respectively. The present solvation model is thus expected to be useful for estimating the solvation free energies of organic molecules.

Properties of hydrophobic free energy found by gas-liquid transfer

The hydrophobic free energy in current use is based on transfer of alkane solutes from liquid alkanes to water, and it has been argued recently that these values are incorrect and should be based instead on gas-liquid transfer data. Hydrophobic free energy is measured here by gas-liquid transfer of hydrocarbon gases from vapor to water. The new definition reduces more than twofold the values of the apparent hydrophobic free energy. Nevertheless, the newly defined hydrophobic free energy is still the dominant factor that drives protein folding as judged by ΔCp, the change in heat capacity, found from the free energy change for heat-induced protein unfolding. The ΔCp for protein unfolding agrees with ΔCp values for solvating hydrocarbon gases and disagrees with ΔCp for breaking peptide hydrogen bonds, which has the opposite sign. The ΔCp values for the enthalpy of liquid-liquid and gas-liquid transfer are similar. The plot of free energy against the apparent solvent-exposed surface area is given for linear alkanes, but only for a single conformation, the extended conformation, of these flexible-chain molecules. The ability of the gas-liquid hydrophobic factor to predict protein stability is tested and reasonable agreement is found, using published data for the dependences on temperature of the unfolding enthalpy of ribonuclease T1 and the solvation enthalpies of the nonpolar and polar groups.

Gas-liquid transfer data used to analyze hydrophobic hydration and find the nature of the Kauzmann-Tanford hydrophobic factor

Hydrophobic free energy for protein folding is currently measured by liquid-liquid transfer, based on an analogy between the folding process and the transfer of a nonpolar solute from water into a reference solvent. The second part of the analogy (transfer into a nonaqueous solvent) is dubious and has been justified by arguing that transfer out of water probably contributes the major part of the free energy change. This assumption is wrong: transfer out of water contributes no more than half the total, often less. Liquid-liquid transfer of the solute from water to liquid alkane is written here as the sum of 2 gas-liquid transfers: (i) out of water into vapor, and (ii) from vapor into liquid alkane. Both gas-liquid transfers have known free energy values for several alkane solutes. The comparable values of the two different transfer reactions are explained by the values, determined in 1991 for three alkane solutes, of the cavity work and the solute-solvent interaction energy. The transfer free energy is the difference between the positive cavity work and the negative solute-solvent interaction energy. The interaction energy has similar values in water and liquid alkane that are intermediate in magnitude between the cavity work in water and in liquid alkane. These properties explain why the transfer free energy has comparable values (with opposite signs) in the two transfers. The current hydrophobic free energy is puzzling and poorly defined and needs a new definition and method of measurement.

Antibiotic selection by the promiscuous aminoglycoside acetyltransferase-(3)-IIIb is thermodynamically achieved through the control of solvent rearrangement

The results presented here show the first known observation of opposite signs of change in heat capacity (ΔC(p)) of two structurally similar ligands binding to the same protein site. Neomycin and paromomycin are aminoglycoside antibiotics that are substrates for the resistance-conferring enzyme, the aminoglycoside acetyltransferase-(3)-IIIb (AAC). These antibiotics are identical to one another except at the 6' position where neomycin has an amine and paromomycin has a hydroxyl. The opposite trends in ΔC(p) of binding of these two drugs to AAC suggest a differential exposure of nonpolar amino acid side chains. Nuclear magnetic resonance experiments further demonstrate significantly different changes in AAC upon interaction with neomycin and paromomycin. Experiments in H(2)O and D(2)O reveal the first observed temperature dependence of solvent and vibrational contributions to ΔC(p). Coenzyme A significantly influences these effects. Together, the data suggest that AAC exploits solvent properties to facilitate favorable thermodynamic selection of antibiotics.