Thermal coefficients of the methyl groups within ubiquitin

Physiological processes such as protein folding and molecular recognition are intricately linked to their dynamic signature, which is reflected in their thermal coefficient. In addition, the local conformational entropy is directly related to the degrees of freedom, which each residue possesses within its conformational space. Therefore, the temperature dependence of the local conformational entropy may provide insight into understanding how local dynamics may affect the stability of proteins. Here, we analyze the temperature dependence of internal methyl group dynamics derived from the cross-correlated relaxation between dipolar couplings of two CH bonds within ubiquitin. Spanning a temperature range from 275 to 308 K, internal methyl group dynamics tend to increase with increasing temperature, which translates to a general increase in local conformational entropy. With this data measured over multiple temperatures, the thermal coefficient of the methyl group order parameter, the characteristic thermal coefficient, and the local heat capacity were obtained. By analyzing the distribution of methyl group thermal coefficients within ubiquitin, we found that the N-terminal region has relatively high thermostability. These results indicate that methyl groups contribute quite appreciably to the total heat capacity of ubiquitin through the regulation of local conformational entropy.     

A thermodynamical study of the clockwork hypothesis proposed by E. Schrödinger

In the present study, the 'clockwork' hypothesis proposed by Schr?dinger was examined from the viewpoint of thermodynamics. Firstly, noticing a unidirectional transfer of entropy in a heat engine, the logic was briefly explained about a close relation between this entropy transfer and an irreversible cycle performed by a working body. Next, paying attention to two fundamental differences between a heat engine and a biological system, we considered an isolated system Asigma consisting of three one-component systems (Ai, A, Ao) and noted a case that the same molecules as the component ones flowed quasistatically into Ai from the outside. Then, the unidirectional flows of the molecules, energy and entropy, which were induced by the above inflow in Asigma, were formulated on the basis of the equilibrium thermodynamics for an open system. Furthermore, it was clarified that the fundamental equation for these flows is the Schr?dinger inequality and that the necessary-sufficient condition for this inequality is the existence of an irreversible cycle performed by A. Here A corresponds to a working body in a heat engine. It was, thus, concluded that the 'clockwork' hypothesis by Schr?dinger is considered to be reasonable for a biological system composed of various irreversible subsystems.     

The effect of alumina nanoparticles on the thermal properties of PMMA: a molecular dynamics simulation

Metal oxides, as one of the most promising flame retardant additives, improve the fire retardant and the thermal stability properties of polymers. In the present study, molecular dynamics (MD) simulations based on the united atom model were applied to study the effect of alumina nanoparticles on the density, thermal conductivity, heat capacity, and thermal diffusivity of isotactic poly(methyl methacrylate) (is-PMMA). Thermal diffusivity of PMMA and PMMA/alumina nanocomposite were investigated through calculating thermal conductivity, density and heat capacity in the range of 300–700?K. Heat capacity can be calculated using fluctuations properties of energy. Thermal conductivity was calculated through the nonequilibrium molecular dynamics (NEMD) simulation by Fourier’s law approach. Our results show that the addition of alumina nanoparticles decreases the heat capacity and increases the glass transition temperature (T_g), thermal conductivity and thermal diffusivity of the PMMA. Therefore, the addition of alumina nanoparticles to PMMA improves the fire retardancy of the polymer. In addition, we illustrate the links between the intermolecular and bulk properties of PMMA in the presence of the alumina nanoparticles.     

Studies on the interaction of the food colorant tartrazine with double stranded deoxyribonucleic acid

Interaction of the food additive tartrazine with double-stranded DNA was studied by spectroscopic and calorimetric techniques. Absorbance studies revealed that tartrazine exhibited hypochromism in the presence of DNA without any bathochromic effects. Minor groove displacement assay of DAPI and Hoechst 33258 suggested that tartrazine binds in the minor groove of DNA. The complexation was predominantly entropy driven with a smaller but favorable enthalpic contribution to the standard molar Gibbs energy. The equilibrium constant was evaluated to be (3.68?±?.08)?×?10⁴?M^?1 at 298.15?K. The negative standard molar heat capacity value along with an enthalpy–entropy compensation phenomenon proposed the involvement of dominant hydrophobic forces in the binding process. Tartrazine enhanced the thermal stability of DNA by 7.53?K under saturation conditions.     

Stability of yeast iso-1-ferricytochrome c as a function of pH and temperature.

Absorbance-detected thermal denaturation studies of the C102T variant of Saccharomyces cerevisiae iso-1-ferricytochrome c were performed between pH 3 and 5. Thermal denaturation in this pH range is reversible, shows no concentration dependence, and is consistent with a 2-state model. Values for free energy (delta GD), enthalpy (delta HD), and entropy (delta SD) of denaturation were determined as functions of pH and temperature. The value of delta GD at 300 K, pH 4.6, is 5.1 +/- 0.3 kcal mol-1. The change in molar heat capacity upon denaturation (delta Cp), determined by the temperature dependence of delta HD as a function of pH (1.37 +/- 0.06 kcal mol-1 K-1), agrees with the value determined by differential scanning calorimetry. pH-dependent changes in the Soret region indicate that a group or groups in the heme environment of the denatured protein, probably 1 or both heme propionates, ionize with a pK near 4. The C102T variant exhibits both enthalpy and entropy convergence with a delta HD of 1.30 kcal mol-1 residue-1 at 373.6 K and a delta SD of 4.24 cal mol-1 K-1 residue-1 at 385.2 K. These values agree with those for other single-domain, globular proteins.     

Common mechanism of thermostability in small α- and β-proteins studied by molecular dynamics

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (T_m) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to T_m. Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH − TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.     

Forming and maintaining a heat engine for quantum biology

Chemical reactions upholding biological functions and structures are the process of measurement taking place among the participating chemical reactants. Chemical reactions occurring in thermal environments are either endothermic or exothermic. In particular, exothermic reactions that can live with temperature gradients of exogenous origin could potentially be competent enough to synthesize a robust quantum as a heat engine. Molecular organizations leading to the origin of the phenomenon of life might have been associated with the emergence of a quantum coherence embodied in a robust heat engine feeding on quantum decoherence. Evolutionary maintenance of a robust quantum heat engine, once appeared, can further be empowered by the build-up of temperature gradients of endogenous origin. Biology enriches the repertoire of quantum mechanics so as to include a robust heat engine as a legitimate member of a quantum in addition to the already established member of a quantum including an atom, molecule, and macromolecule.     

The thermal properties of binary structure sI clathrate hydrate from molecular dynamics simulation

In this work, we present temperature dependence of lattice parameter and normalised lattice parameter in the atmospheric pressure and 120 bar and also pressure dependence of unit cell volume and normalised unit cell volume at 150 and 250?K for variety guests with different size, polarity and guest–host hydrogen bonding capability such as trimethylene oxide (TMO), ethylene oxide (EO), formaldehyde (FA), cyclobutane (CB), cyclopropane (CP) and ethane (Et) in the large cages with CH₄ in small cages of sI clathrate hydrates by molecular dynamics simulations. The obtained values of lattice parameters for the guest species are compatible with the experimental values. These clathrate hydrates are simulated with TIP4P/ice four-site water potential. Herein, isobaric thermal expansivity and isothermal compressibility are calculated at a temperature range of 50–250?K and a wide pressure range. These structural properties have been compared for guests which they are isoelectronic and have similar masses but with different size and polarity. We use molecular dynamics simulations to relate microscopic guest properties, like guest–host hydrogen bonding to macroscopic sI clathrate hydrate properties. The temperature dependence of thermodynamic properties such as constant-volume and constant-pressure heat capacity is presented in the atmospheric pressure for these guest species.     

Molecular dynamics simulations of the thermal stability of tteRBP and ecRBP

Molecular dynamics simulations were performed for investigating the thermal stability of the extremely thermophilic Thermoanaerobacter tengcongensis ribose binding protein (tteRBP) and the mesophilic homologous Escherichia coli ribose binding protein (ecRBP). The simulations for the two proteins were carried out under the room temperature (300?K) and the optimal activity temperature (tteRBP 375?K and ecRBP 329?K), respectively. The comparative analyses of the trajectories show that the two proteins have stable overall structures at the two temperatures; further analyses indicate that they both have strong side-chain interactions and different backbone flexibilities at the different temperatures. The tteRBP 375?K and ecRBP 329?K have stronger internal motion and higher flexibility than tteRBP 300?K and ecRBP 300?K, respectively, it is noted that the flexibility of tteRBP is much higher than that of ecRBP at the two temperatures. Therefore, tteRBP 375?K can adapt to high temperature due to its higher flexibility of backbone. Combining with the researches by Cuneo et al., it is concluded that the side-chain interactions and flexibility of backbone are both the key factors to maintain thermal stability of the two proteins.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:22     

Comparative mechanical unfolding studies of spectrin domains R15, R16 and R17

Spectrins belong to repetitive three-helix bundle proteins that have vital functions in multicellular organisms and are of potential value in nanotechnology. To reveal the unique physical features of repeat proteins we have studied the structural and mechanical properties of three repeats of chicken brain α-spectrin (R15, R16 and R17) at the atomic level under stretching at constant velocities (0.01, 0.05 and 0.1?Å·ps^?1) and constant forces (700 and 900?pN) using molecular dynamics (MD) simulations at T?=?300?K. 114 independent MD simulations were performed and their analysis has been done. Despite structural similarity of these domains we have found that R15 is less mechanically stable than R16, which is less stable than R17. This result is in agreement with the thermal unfolding rates. Moreover, we have observed the relationship between mechanical stability, flexibility of the domains and the number of aromatic residues involved in aromatic clusters.     

Rare Dissipative Transitions Punctuate the Initiation of Chemical Denaturation in Proteins

Protein unfolding dynamics are bound by their degree of entropy production, a quantity that relates the amount of heat dissipated by a nonequilibrium process to a system’s forward and time-reversed trajectories. We here explore the statistics of heat dissipation that emerge in protein molecules subjected to a chemical denaturant. Coupling large molecular dynamics datasets and Markov state models with the theory of entropy production, we demonstrate that dissipative processes can be rigorously characterized over the course of the urea-induced unfolding of the protein chymotrypsin inhibitor 2. By enumerating full entropy production probability distributions as a function of time, we first illustrate that distinct passive and dissipative regimes are present in the denaturation dynamics. Within the dissipative dynamical region, we next find that chymotrypsin inhibitor 2 is strongly driven into unfolded states in which the protein’s hydrophobic core has been penetrated by urea molecules and disintegrated. Detailed analyses reveal that urea’s interruption of key hydrophobic contacts between core residues causes many of the protein’s native structural features to dissolve.     

Comparative simulation studies of native and single-site mutant human beta-defensin-1 peptides

Human defensins play important roles in a broad range of biological functions, such as microbial defense and immunity. Yet, little is known about their molecular properties, i.e. secondary structure stability, structural variability, important side chain interactions, surface charge distribution, and resistance to thermal fluctuations, and how these properties are related to their functions. To assess these factors, we studied the native human β-defensin-1 monomer and dimer as well as several single-site mutants using molecular dynamics simulations. The results showed that disulfide bonds are important determinants in maintaining the defensins’ structural integrity, as no structural transitions were observed at 300?K and only minor structural unfolding was detected upon heating to 500?K. The α-helix was less thermally stable than the core β-sheet structure held together by hydrogen bonds and hydrophobic interactions. The monomer α-helix stability was directly correlated, whereas the end-to-end distance was inversely correlated to the experimentally measured β-defensin-1 chemotactic activity, in the order: mutant 2 (Gln24Glu)?>?mutant 3 (Lys31Ala)?=?wild type?>?mutant 1 (Asn4Ala). The structural stability of the β-defensin-1 dimer species exhibited an inverse correlation to their chemotactic activity. In dimers formed by mutants 2 and 3, we observed sliding of one monomer upon the surface of the other in the absence of unbinding. This dynamic sliding feature may enhance the molecular oligomerization of β-defensin-1 peptides contributing to their antibacterial activity. It could also help these peptides orient correctly in the CC chemokine receptor 6 binding site, thereby initiating their chemotactic activity. In agreement with this notion, the remarkable sliding behavior was observed only for the mutants with the highest chemotactic activity.     

Water penetration and escape in proteins

The kinetics of water penetration and escape in cytochrome c (cyt c) is studied by molecular dynamics (MD) simulations at various temperatures. Water molecules that penetrate the protein interior during the course of an MD simulation are identified by monitoring the number of water molecules in the first coordination shell (within 3.5 A) of each water molecule in the system. Water molecules in the interior of cyt c have 0-3 water molecules in their first hydration shell and this coordination number persists for extended periods of time. At T = 300 K we identify over 200 events in which water molecules penetrate the protein and reside inside for at least 5 picoseconds (ps) within a 1.5 nanoseconds (ns) time period. Twenty-seven (27) water molecules reside for at least 300 ps, 17 water molecules reside in the protein interior for times longer than 500 ps, and two interior water molecules do not escape; at T = 360 K one water molecule does not escape; at 430 K all water molecules exchange. Some of the internal water molecules show mean square displacements (MSD) of 1 A2 characteristic of structural waters. Others show MSD as large as 12 A2, suggesting that some of these water molecules occupy transient cavities and diffuse extensively within the protein. Motions of protein-bound water molecules are rotationally hindred, but show large librations. Analysis of the kinetics of water escape in terms of a survival time correlation function shows a power law behavior in time that can be interpreted in terms of a broad distribution of energy barriers, relative to kappa BT, for water exchange. At T = 300 K estimates of the roughness of the activation energy distribution is 4-10 kJ/mol (2-4 kappa BT). Activation enthalpies for water escape are 6-23 kJ/mol. The difference in activation entropies between fast exchanging (0.01 ns) and slow exchanging (0.1-1 ns) water molecules is -27 J/K/mol. Dunitz (Science 1997;264:670.) has estimated the maximum entropy loss of a water molecule due to binding to be 28 J/K/mol. Therefore, our results suggest that the entropy of interior water molecules is similar to entropy of bulk water.     

Vacuum Ultraviolet Studies on the Nature of the Radiation Inactivation of Trypsin

Trypsin, in powder form and in frozen D(2)O-glucose solutions, at temperatures from 100 degrees to 300 degrees K, was excited with vacuum ultraviolet and near ultraviolet radiation to determine how absorbed photon energy is partitioned into radiative, nonradiative and/or inactivating processes; at 300 degrees K most of the absorbed energy is not reemitted, so that it (0.98-0.99 for excitation at 120 nm and 0.75-0.90 at 280 nm) is potentially available for inactivation. Since the effects of excitation wavelength and temperature on the emission quenching yields are generally different from those on the inactivation yields of dry trypsin, the mere retention of quenched energy by an enzyme does not necessarily lead to its inactivation. Thus, as predicted previously, the radiation inactivation of trypsin must proceed by rather specific mechanisms which undoubtedly depend upon environment-sensitive processes, since the nature of the molecular environment can modify the partitioning of energy so significantly; for example, there are differences in the phosphorescence-to-fluorescence ratio, in the activation energy for quenching, and in the lifetimes and kinetics of the decay of phosphorescence when trypsin in frozen glasses and dry trypsin are excited by various wavelengths of ultraviolet radiation.     

A new theoretical approach to biological self-assembly

Upon biological self-assembly, the number of accessible translational configurations of water in the system increases considerably, leading to a large gain in water entropy. It is important to calculate the solvation entropy of a biomolecule with a prescribed structure by accounting for the change in water–water correlations caused by solute insertion. Modeling water as a dielectric continuum is not capable of capturing the physical essence of the water entropy effect. As a reliable tool, we propose a hybrid of the angle-dependent integral equation theory combined with a multipolar water model and a morphometric approach. Using our methods wherein the water entropy effect is treated as the key factor, we can elucidate a variety of processes such as protein folding, cold, pressure, and heat denaturating of a protein, molecular recognition, ordered association of proteins such as amyloid fibril formation, and functioning of ATP-driven proteins.     

The nature of the free energy barriers to two-state folding

We introduce a simple procedure to analyze the temperature dependence of the folding and unfolding rates of two-state proteins. We start from the simple transition-state-like rate expression: k = D(eff)exp(-DeltaG(TS)/RT), in which upper and lower bounds for the intra-chain effective diffusion coefficient (D(eff)) are obtained empirically using the timescales of elementary processes in protein folding. From the changes in DeltaG(TS) as a function of temperature, we calculate enthalpies and heat capacities of activation, together with the more elusive entropies of activation. We then estimate the conformational entropy of the transition state by extrapolation to the temperature at which the solvation entropy vanishes by cancellation between polar and apolar terms. This approach is based on the convergence temperatures for the entropy of solvating apolar (approximately 385 K) and polar groups (approximately 335 K), the assumption that the structural properties of the transition state are somewhere in between the unfolded and folded states, and the established relationship between observed heat capacity and solvent accessibility.1 To circumvent the lack of structural information about transition states, we use the empirically determined heat capacities of activation as constraints to identify the extreme values of the transition state conformational entropy that are consistent with experiment. The application of this simple approach to six two-state folding proteins for which there is temperature-dependent data available in the literature provides important clues about protein folding. For these six proteins, we obtain an average equilibrium cost in conformational entropy of -4.3 cal x mol(-1)K(-1)per residue, which is in close agreement to previous empirical and computational estimates of the same quantity. Furthermore, we find that all these proteins have a conformationally diverse transition state, with more than half of the conformational entropy of the unfolded state. In agreement with predictions from theory and computer simulations, the transition state signals the change from a regime dominated by loss in conformational entropy to one driven by the gain in stabilization free energy (i.e., including protein interactions and solvation effects). Moreover, the height of the barrier is determined by how much stabilization free energy is realized at that point, which is related to the relative contribution of local versus non-local interactions. A remarkable observation is that the fraction of conformational entropy per residue that is present in the transition state is very similar for the six proteins in this study. Based on this commonality, we propose that the observed change in thermodynamic regime is connected to a change in the pattern of structure formation: from one driven by formation of pairwise interactions to one dominated by coupling of the networks of interactions involved in forming the protein core. In this framework, the barrier to two-state folding is crossed when the folding protein reaches a "critical native density" that allows expulsion of remaining interstitial water and consolidation of the core. The principle of critical native density should be general for all two-state proteins, but can accommodate different folding mechanisms depending on the particularities of the structure and sequence.     

Experimental and computational studies on the binding of diazinon to human serum albumin

In the present research, the binding properties of diazinon (DZN), as an organophosphorus herbicide, to human serum albumin (HSA) were investigated using combination of spectroscopic, electrochemistry, and molecular modeling techniques. Changes in the UV–Vis and FT-IR spectra were observed upon ligand binding along with a significant degree of tryptophan fluorescence quenching on complex formation. The obtained results from spectroscopic and electrochemistry experiments along with the computational studies suggest that DZN binds to residues located in subdomains IIA of HSA with binding constant about 1410.9 M^?1 at 300 K. From the thermodynamic parameters calculated according to the van’t Hoff equation, the enthalpy change ΔH° and entropy change ΔS° were found to be ?16.695 and 0.116 KJ/mol K, respectively. The primary binding pattern is determined by hydrophobic interaction and hydrogen binding occurring in so-called site I of HSA. DZN could slightly alter the secondary structure of HSA. All of experimental results are supported by computational techniques such as docking and molecular dynamics simulation using a HSA crystal model.     

Base Pairing at the Stem-loop Junction in the SL1 Kissing Complex of HIV-1 RNA: A Thermodynamic Study Probed by Molecular Dynamics Simulation

Abstract

The thermodynamics of the opening/closure process of a GC base pair located at the stem-loop junction of the SL1 sequence from HIV-1_Lai genomic RNA was investigated in the context of a loop-loop homodimer (or kissing complex) using molecular dynamics simulation. The free energy, enthalpy and entropy changes for the closing reaction are 0 kcal·mol^?1, ?11 kcal·mol^?1and ?0.037 kcal·mol^?1-K^?1 at 300° K respectively. Furthermore it is found that the free energy change is the same for the formation of a 11 nucleotide loop closed with UG and for the formation of a 9 nucleotide loop closed with GC. Our study evidences the high flexibility of the nucleotides at the stem-loop junction explaining the weak stability of this structure.     

Comparative molecular dynamics simulation studies for determining factors contributing to the thermostability of chemotaxis protein “CheY”

Comparative molecular dynamics simulations of chemotaxis protein “CheY” from thermophilic origin Thermotoga maritima and its mesophilic counterpart Salmonella enterica have been performed for 10?ns each at 300 and 350?K, and 20?ns each at 400 and 450?K. The trajectories were analyzed in terms of different factors like root-mean-square deviation, root-mean-square fluctuation, radius of gyration, solvent accessible surface area, H-bonds, salt bridge content, and protein–solvent interactions which indicate distinct differences between the two of them. The two proteins also follow dissimilar unfolding pathways. The overall flexibility calculated by the trace of the diagonalized covariance matrix displays similar flexibility of both the proteins near their optimum growth temperatures. However, at higher temperatures mesophilic protein shows increased overall flexibility than its thermophilic counterpart. Principal component analysis also indicates that the essential subspaces explored by the simulations of two proteins at different temperatures are nonoverlapping and they show significantly different directions of motion. However, there are significant overlaps within the trajectories and similar direction of motions are observed for both proteins at 300?K. Overall, the mesophilic protein leads to increased conformational sampling of the phase space than its thermophilic counterpart. This is the first ever study of thermostability of CheY protein homologs by using protein dynamism as a main impact. Our study might be used as a model for studying the molecular basis of thermostability of two homologous proteins from two organisms living at different temperatures with less visible differences.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.