Thermodynamics of site-specific small molecular ion interactions with DNA duplex: a molecular dynamics study

The stability and dynamics of a double-stranded DNA (dsDNA) is affected by the preferential occupancy of small monovalent molecular ions. Small metal and molecular ions such as sodium and alkyl ammonium have crucial biological functions in human body, affect the thermodynamic stability of the duplex DNA and exhibit preferential binding. Here, using atomistic molecular dynamics simulations, we investigate the preferential binding of metal ion such as Na⁺ and molecular ions such as tetramethyl ammonium (TMA⁺) and 2-hydroxy-N,N,N-trimethylethanaminium (CHO⁺) to double-stranded DNA. The thermodynamic driving force for a particular molecular ion-DNA interaction is determined by decomposing the free energy of binding into its entropic and enthalpic contributions. Our simulations show that each of these molecular ions preferentially binds to the minor groove of the DNA and the extent of binding is highest for CHO⁺. The ion binding processes are found to be entropically favourable. In addition, the contribution of hydrophobic effects towards the entropic stabilisation (in case of TMA⁺) and the effect of hydrogen bonding contributing to enthalpic stabilisation (in case of CHO⁺) have also been investigated.     

Enthalpie and Entropie Contributions in the Transesterification of Sucrose: Computational Study of Lipases and Subtilisin

Abstract

Transesterification of sucrose with fatty acids catalyzed by subtilisin Carlsberg occurs with regioselectivity that is different from that in lipases. Thermomyces lanuginosus lipase (TlL) and Candida antarctica lipase B (CALB) catalyze synthesis at positions 6 and 6′, with differing abilities, while subtilisin catalysis leads to the l′-acylated sucrose. The catalytic machinery in lipases is approximately mirrored in subtilisins but different pocket morphologies including size, shape, and rearrangement of the catalytic elements underlies the differing regioselectivities. The thermodynamic consequences of these differences on the above reactions have been explored systematically using computational methods, determining the free energies of interaction of the putative transition-state adducts. Analysis of the conformers with the lowest transition state energies (protein-ligand interactions and vibrational entropy contributions) indicates that enthalpic factors control specificities in lipases while entropic factors are more important in subtilisin.     

The thermodynamics of protein interactions with essential first row transition metals

BackgroundThe binding of metal ions to proteins is a crucial process required for their catalytic activity, structural stability and/or functional regulation. Isothermal titration calorimetry provides a wealth of fundamental information which when combined with structural data allow for a much deeper understanding of the underlying molecular mechanism.Scope of reviewA rigorous understanding of any molecular interaction requires in part an in-depth quantification of its thermodynamic properties. Here, we provide an overview of recent studies that have used ITC to quantify the interaction of essential first row transition metals with relevant proteins and highlight major findings from these thermodynamic studies.General significanceThe thermodynamic characterization of metal ion–protein interactions is one important step to understanding the role that metal ions play in living systems. Such characterization has important implications not only to elucidating proteins' structure-function relationships and biological properties but also in the biotechnology sector, medicine and drug design particularly since a number of metal ions are involved in several neurodegenerative diseases.Major conclusionsIsothermal titration calorimetry measurements can provide complete thermodynamic profiles of any molecular interaction through the simultaneous determination of the reaction binding stoichiometry, binding affinity as well as the enthalpic and entropic contributions to the free energy change thus enabling a more in-depth understanding of the nature of these interactions. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

Structure,dynamics, and molecular inhibition of the Staphylococcus aureus m1A22-tRNA methyltransferase TrmK

The enzyme m¹A22-tRNA methyltransferase (TrmK) catalyzes the transfer of a methyl group to the N1 of adenine 22 in bacterial tRNAs. TrmK is essential for Staphylococcus aureus survival during infection but has no homolog in mammals, making it a promising target for antibiotic development. Here, we characterize the structure and function of S. aureus TrmK (SaTrmK) using X-ray crystallography, binding assays, and molecular dynamics simulations. We report crystal structures for the SaTrmK apoenzyme as well as in complexes with methyl donor SAM and co-product product SAH. Isothermal titration calorimetry showed that SAM binds to the enzyme with favorable but modest enthalpic and entropic contributions, whereas SAH binding leads to an entropic penalty compensated for by a large favorable enthalpic contribution. Molecular dynamics simulations point to specific motions of the C-terminal domain being altered by SAM binding, which might have implications for tRNA recruitment. In addition, activity assays for SaTrmK-catalyzed methylation of A22 mutants of tRNA^Leu demonstrate that the adenine at position 22 is absolutely essential. In silico screening of compounds suggested the multifunctional organic toxin plumbagin as a potential inhibitor of TrmK, which was confirmed by activity measurements. Furthermore, LC-MS data indicated the protein was covalently modified by one equivalent of the inhibitor, and proteolytic digestion coupled with LC-MS identified Cys92 in the vicinity of the SAM-binding site as the sole residue modified. These results identify a cryptic binding pocket of SaTrmK, laying a foundation for future structure-based drug discovery.     

Kinetic and thermodynamic analysis of dimerization inhibitors binding to HIV protease monomers by surface plasmon resonance

The analysis of kinetic and thermodynamic parameters of binding of peptide and nonpeptide dimerization inhibitors of HIV protease (HIVp) to the enzyme monomers immobilized on an optical chip has been studied by surface plasmon resonance. The molecular interactions were investigated at different inhibitor concentrations (0–80 μM) and temperatures (15–35°C). Determination of kinetic (k _on, k _off), equilibrium (K _d), and thermodynamic (ΔG, ΔH, and -TΔS) has shown that both inhibitors are characterized by similar interaction parameters and the entropic term (-TΔS) of about −20 kcal/mol is the main driving force for the HIVp complex formation with the inhibitors, while the positive value (14 kcal/mol) of the enthalpic term (ΔH) counteracted the complex formation.     

Molecular mechanisms of temperature adaptation in fish myofibrillar adenosine triphosphatases

Summary Studies have been carried out on the Mg²⁺ Ca²⁺-myofibrillar ATPase from the muscles of fish adapted to different environmental temperatures. The thermal stability of the ATPase is strongly correlated with mean habitat temperature. Activities of Antarctic fish ATPases are significantly higher at low temperatures than those of temperate and tropical water species. The effects of ionic strength on ATPase activity have also been studied. The Gibbs free energy of activation (G ^#) was found to increase and enzyme activity decrease with increasing ionic strength within the physiological temperature range of each species. Significantly lower values of G ^#, of around 1 Kcal/mole, are obtained for the ATPase of cold-adapted compared to tropical fish. Enthalpic and entropic activation energies were also reduced in the cold adapted ATPases. It is postulated that the reduction of the enthalpic activation term in the cold adapted enzyme confers the advantage of reducing the temperature sensitivity of the rate limiting step thus partly compensating for the low heat content of the cellular environment. Possible molecular mechanisms of temperature compensation in fish myofibrillar ATPase are discussed.     

Drug‐binding energetics of human α‐1‐acid glycoprotein assessed by isothermal titration calorimetry and molecular docking simulations

This study utilizes sensitive, modern isothermal titration calorimetric methods to characterize the microscopic thermodynamic parameters that drive the binding of basic drugs to α‐1‐acid glycoprotein (AGP) and thereby rationalize the thermodynamic data in relation to docking models and crystallographic structures of the drug–AGP complexes. The binding of basic compounds from the tricyclic antidepressant series, together with miaserine, chlorpromazine, disopyramide and cimetidine, all displayed an exothermically driven binding interaction with AGP. The impact of protonation/deprotonation events, ionic strength, temperature and the individual selectivity of the A and F1*S AGP variants on drug‐binding thermodynamics was characterized. A correlation plot of the thermodynamic parameters for all of the test compounds revealed that an enthalpy–entropy compensation is in effect. The exothermic binding energetics of the test compounds were driven by a combination of favorable (negative) enthalpic (?Hº) and favorable (positive) entropic (?Sº) contributions to the Gibbs free energy (?Gº). Collectively, the data imply that the free energies that drive drug binding to AGP and its relationship to drug serum residency evolve from the complex interplay of enthalpic and entropic forces from interactions with explicit combinations of hydrophobic and polar side‐chain sub‐domains within the multi‐lobed AGP ligand binding cavity.Copyright © 2012 John Wiley & Sons, Ltd.     

Purification and Some Properties of Carboxylesterase from Arthrobacter globiformis; Stereoselective Hydrolysis of Ethyl Chrysanthemate

An esterase with excellent stereoselectivity for (+)-trans-ethyl chrysanthemate was purified to homogeneity from Arthrobacter globiformis SC-6-98-28. The purified enzyme hydrolyzed a mixture of ethyl chrysanthemate isomers stereoselectively to produce (+)-trans-acid with 100% stereoisomeric purity. The apparent molecular weight of the purified enzyme was 43,000 on SDS–polyacrylamide gel electrophoresis, and 94,000 on gel filtration chromatography. The optimum conditions for the ester hydrolysis were pH 10.0 at 45°C. The purified esterase hydrolyzed short-chain fatty acid esters, but did not have detectable activity on long-chain water-insoluble fatty acid esters. The enzyme activity was inbibited by diisopropyl fluorophosphate and phenylmethylsulfonyl fluoride.     

Long-Timescale Molecular-Dynamics Simulations of the Major Urinary Protein Provide Atomistic Interpretations of the Unusual Thermodynamics of Ligand Binding

The mouse major urinary protein (MUP) has proved to be an intriguing test bed for detailed studies on protein-ligand recognition. NMR, calorimetric, and modeling investigations have revealed that the thermodynamics of ligand binding involve a complex interplay between competing enthalpic and entropic terms. We performed six independent, 1.2 μs molecular-dynamics simulations on MUP—three replicates on the apo-protein, and three on the complex with the pheromone isobutylmethoxypyrazine. Our findings provide the most comprehensive picture to date of the structure and dynamics of MUP, and how they are modulated by ligand binding. The mechanical pathways by which amino acid side chains can transmit information regarding ligand binding to surface loops and either increase or decrease their flexibility (entropy-entropy compensation) are identified. Dewetting of the highly hydrophobic binding cavity is confirmed, and the results reveal an aspect of ligand binding that was not observed in earlier, shorter simulations: bound ligand retains extensive rotational freedom. Both of these features have significant implications for interpretations of the entropic component of binding. More generally, these simulations test the ability of current molecular simulation methods to produce a reliable and reproducible picture of protein dynamics on the microsecond timescale.     

Collapse Dynamics of Single Proteins Extended by Force

Single-molecule force spectroscopy has opened up new approaches to the study of protein dynamics. For example, an extended protein folding after an abrupt quench in the pulling force was shown to follow variable collapse trajectories marked by well-defined stages that departed from the expected two-state folding behavior that is commonly observed in bulk. Here, we explain these observations by developing a simple approach that models the free energy of a mechanically extended protein as a combination of an entropic elasticity term and a short-range potential representing enthalpic hydrophobic interactions. The resulting free energy of the molecule shows a force-dependent energy barrier of magnitude, ΔE = ɛ(F − F_c)^3/2, separating the enthalpic and entropic minima that vanishes at a critical force F_c. By solving the Langevin equation under conditions of a force quench, we generate folding trajectories corresponding to the diffusional collapse of an extended polypeptide. The predicted trajectories reproduce the different stages of collapse, as well as the magnitude and time course of the collapse trajectories observed experimentally in ubiquitin and I27 protein monomers. Our observations validate the force-clamp technique as a powerful approach to determining the free-energy landscape of proteins collapsing and folding from extended states.     

Effects of nonspecific ion-protein interactions on the redox chemistry of cytochrome c

 The effects of the ionic atmosphere on the enthalpic and entropic contributions to the reduction potential of native (state III) beef heart cytochrome c have been determined through variable-temperature direct electrochemistry experiments. At neutral or slightly alkaline pH values, from 5 to 50  °C, the reduction enthalpy and entropy become less negative with decreasing ionic strength. The reduction entropy extrapolated at null ionic strength is approximately zero, indicating that, in the absence of the screening effects of the salt ions on the network of the electrostatic interactions at the protein-solvent interface, the solvation properties and the conformational flexibility of the two redox states are comparable. The moderate decrease in E°′ observed with increasing ionic strength [ΔE°′_IS =(E°′) _{I =0.1 M}–(E°′) _{I =0 M}=–0.035 V at 25  °C], once the compensating enthalpic and entropic effects of the salt-induced changes in the hydrogen bonding within the hydration sphere of the molecule in the two redox states are factorized out, results in being ultimately determined by the stabilizing enthalpic effect of the negatively charged ionic atmosphere on the ferri form. At pH 9, the ionic strength dependence of the reduction termodynamics of cytochrome c follows distinctive patterns, possibly as a result of specific binding of the hydroxide ion to the protein. A decrease in ionic strength at constant pH, as well as a pH increase at constant ionic strength, induces a depression of the temperature of the transition from the low-T to high-T conformer of cytochrome c, which suggests that a temperature-induced decrease in the pK _a for a residue deprotonation is the key event of this conformational change. Received: 7 April 1999 / Accepted: 19 July 1999     

Substrate specificity and nucleotides binding properties of NM23H2/nucleoside diphosphate kinase homolog from <Emphasis Type="Italic">Plasmodium falciparum</Emphasis>

Nucleoside diphosphate kinases (NDKs) play a key role in maintaining the intracellular energy resources as well as the balance of nucleotide pools. Recently, attention has been directed to NDKs owing to its role in activating various chemotherapeutic agents. The binding affinity of different nucleotides with P. falciparum NDK was varied according to the following order ADP ~ GDP > dGDP > dADP > dTDP > CDP > dCDP > UDP. The binding of purines nucleotides was stronger than pyrimidines. Furthermore, PfNDK showed more preferences to ribonucleotides over deoxyribonucleotides. Pyrimidines showed lower negative free energy compared with that of purines. The interaction of all nucleotides showed favorable enthalpic and entropic terms. However, the enthalpic terms were the main deriving forces for purine nucleotides, while the entropic contributions were the predominant forces for pyrimidines. Interestingly, TDP showed marked affinity and more favorable enthalpic and less entropic contributions. In conclusion, the size of nucleotide was the critical factor in PfNDK ligand affinity.     

Kinetic and thermodynamic studies of peptidyltransferase in ribosomes from the extreme thermophile Thermus thermophilus

Throughout evolution, emerging organisms survived by adapting existing biochemical processes to new reaction conditions. Simple protein enzymes balanced changes in structural stability with changes that permitted optimal catalysis by adjustments in both entropic and enthalpic contributions to the free energy of activation for the reaction. Study of adaptive mechanisms by large multicomponent enzymes such as the ribosome has been largely unexplored. Here we have determined the kinetic and thermodynamic parameters of peptidyltransferase in ribosomes from the extreme thermophile Thermus thermophilus. Activity of thermophilic enzymes can be assayed over a wide range of temperatures, enabling one to measure accurate catalytic rates and determine enthalpic and entropic contributions to the free energy of activation of the reaction. Differences in the reaction conditions used here and in published studies on mesophilic ribosomes prevent direct comparison, but our data on Thermus ribosomes suggest that these ribosomes have adapted to changing environments using the same strategies as simple protein enzymes, balancing stability and flexibility without loss of catalytic rate. This strategy must be a very ancient process, perhaps first used by primitive ribosomes in the RNA World.     

Constituents of the Essential Oil of Chrysanthemum japonense var. debile

One hundred fifty strains of actinomycetes were isolated from soils on plate cultures containing beet arabinan as the sole carbon source. About one-third of the culture fluids were found to have arabinosidase activity. A wild-type strain, Streptomyces sp. No. 17-1, was selected as the best producer of arabinosidase. The highest enzymatic activity was obtained in the culture fluid when the initial pH was adjusted to 9.0. An α-l-arabinofuranosidase was highly purified from the culture filtrate of No. 17-1 by combining column chromatography on DEAE-cellulose, gel filtration on Sephadex G-100, and isoelectric focusing. The molecular weight of the purified enzyme was estimated to be about 92, 000, and its isoelectric point was pH 4.4. The enzymatic activity was maximum at pH 6.0 and was completely inhibited by Hg²⁺. The apparent Km value of the enzyme for p-nitrophenyl-α-l-arabinofuranoside was determined to be 3.6 mM.     

Spatially resolved free-energy contributions of native fold and molten-globule-like Crambin

The folding stability of a protein is governed by the free-energy difference between its folded and unfolded states, which results from a delicate balance of much larger but almost compensating enthalpic and entropic contributions. The balance can therefore easily be shifted by an external disturbance, such as a mutation of a single amino acid or a change of temperature, in which case the protein unfolds. Effects such as cold denaturation, in which a protein unfolds because of cooling, provide evidence that proteins are strongly stabilized by the solvent entropy contribution to the free-energy balance. However, the molecular mechanisms behind this solvent-driven stability, their quantitative contribution in relation to other free-energy contributions, and how the involved solvent thermodynamics is affected by individual amino acids are largely unclear. Therefore, we addressed these questions using atomistic molecular dynamics simulations of the small protein Crambin in its native fold and a molten-globule-like conformation, which here served as a model for the unfolded state. The free-energy difference between these conformations was decomposed into enthalpic and entropic contributions from the protein and spatially resolved solvent contributions using the nonparametric method Per|Mut. From the spatial resolution, we quantified the local effects on the solvent free-energy difference at each amino acid and identified dependencies of the local enthalpy and entropy on the protein curvature. We identified a strong stabilization of the native fold by almost 500 kJ mol⁻¹ due to the solvent entropy, revealing it as an essential contribution to the total free-energy difference of (53 ± 84) kJ mol⁻¹. Remarkably, more than half of the solvent entropy contribution arose from induced water correlations.     

Enzymatic characterization of an extracellular glucoamylase from Tricholoma matsutake and its cloning and secretory expression in Pichia pastoris

ABSTRACT

A glucoamylase from the ectomycorrhizal fungus Tricholoma matsutake (TmGLA) was purified 33.2-fold to homogeneity as a single monomeric glycoprotein with a molecular mass of 63.9 kDa. Maximum activity was observed at 60°C and pH 5.0. The enzyme is active down to 50°C and in the pH range of 4.0–6.0, and its activity is strongly inhibited by Ag⁺. It degrades α-1,4- and α-1,6-glycosidic linkages in various polysaccharides. Its gene (TmGlu1) was cloned using information from the enzyme’s internal amino acid sequences and the whole genome sequence of T. matsutake NBRC 30605. The deduced amino acid sequence showed clear homology with those of GH family 15 proteins. Pichia pastoris transformed with TmGlu1 secreted the active enzyme in a glycosylated form, and its characteristics were the same as the native enzyme.     

Binding Thermodynamics of Ferredoxin:NADP Reductase: Two Different Protein Substrates and One Energetics

The thermodynamics of the formation of binary and ternary complexes between Anabaena PCC 7119 FNR and its substrates, NADP⁺ and Fd, or Fld, has been studied by ITC. Despite structural dissimilarities, the main difference between Fd and Fld binding to FNR relates to hydrophobicity, reflected in different binding heat capacity and number of water molecules released from the interface. At pH 8, the formation of the binary complexes is both enthalpically and entropically driven, accompanied by the protonation of at least one ionizable group. His²⁹⁹ FNR has been identified as the main responsible for the proton exchange observed. However, at pH 10, where no protonation occurs and intrinsic binding parameters can be obtained, the formation of the binary complexes is entropically driven, with negligible enthalpic contribution. Absence of the FMN cofactor in Fld does not alter significantly the strength of the interaction, but considerably modifies the enthalpic and entropic contributions, suggesting a different binding mode. Ternary complexes show negative cooperativity (6-fold and 11-fold reduction in binding affinity, respectively), and an increase in the enthalpic contribution (more favorable) and a decrease in the entropic contribution (less favorable), with regard to the binary complexes energetics.     

Spatially resolved free-energy contributions of native fold and molten-globule-like Crambin

The folding stability of a protein is governed by the free-energy difference between its folded and unfolded states, which results from a delicate balance of much larger but almost compensating enthalpic and entropic contributions. The balance can therefore easily be shifted by an external disturbance, such as a mutation of a single amino acid or a change of temperature, in which case the protein unfolds. Effects such as cold denaturation, in which a protein unfolds because of cooling, provide evidence that proteins are strongly stabilized by the solvent entropy contribution to the free-energy balance. However, the molecular mechanisms behind this solvent-driven stability, their quantitative contribution in relation to other free-energy contributions, and how the involved solvent thermodynamics is affected by individual amino acids are largely unclear. Therefore, we addressed these questions using atomistic molecular dynamics simulations of the small protein Crambin in its native fold and a molten-globule-like conformation, which here served as a model for the unfolded state. The free-energy difference between these conformations was decomposed into enthalpic and entropic contributions from the protein and spatially resolved solvent contributions using the nonparametric method Per|Mut. From the spatial resolution, we quantified the local effects on the solvent free-energy difference at each amino acid and identified dependencies of the local enthalpy and entropy on the protein curvature. We identified a strong stabilization of the native fold by almost 500 kJ mol⁻¹ due to the solvent entropy, revealing it as an essential contribution to the total free-energy difference of (53 ± 84) kJ mol⁻¹. Remarkably, more than half of the solvent entropy contribution arose from induced water correlations.