Thermodynamic and kinetic determinants of Thermotoga maritima cold shock protein stability: a structural and dynamic analysis

The cold shock protein (CSP) from hyperthermophile Thermotoga maritima (TmCSP) is only marginally stable (DeltaG(T(opt)) = 0.3 kcal/mol) at 353 K, the optimum environmental temperature (T(opt)) for T. maritima. In comparison, homologous CSPs from E. coli (DeltaG(T(opt)) = 2.2 kcal/mol) and B. subtilis (DeltaG(T(opt)) = 1.5 kcal/mol) are at least five times more stable at 310 K, the T(opt) for the mesophiles. Yet at the room temperature, TmCSP is more stable (DeltaG(T(R)) = 4.7 kcal/mol) than its homologues (DeltaG(T(R)) = 3.0 kcal/mol for E. coli CSP and DeltaG(T(R)) = 2.1 kcal/mol for B. subtilis CSP). This unique observation suggests that kinetic, rather than thermodynamic, barriers toward unfolding might help TmCSP native structure at high temperatures. Consistently, the unfolding rate of TmCSP is considerably slower than its homologues. High temperature (600 K) complete unfolding molecular dynamics (MD) simulations of TmCSP support our hypothesis and reveal an unfolding scheme unique to TmCSP. For all the studied homologues of TmCSP, the unfolding process first starts at the C-terminal region and N-terminal region unfolds in the end. But for TmCSP, both the terminals resist unfolding for consistently longer simulation times and, in the end, unfold simultaneously. In TmCSP, the C-terminal region is better fortified and has better interactions with the N-terminal region due to the charged residues, R2, E47, E49, H61, K63, and E66, being in spatial vicinity. The electrostatic interactions among these residues are unique to TmCSP. Consistently, the room temperature MD simulations show that TmCSP is more rigid at its N- and C-termini as compared to its homologues from E. coli, B. subtilis, and B. caldolyticus.     

The glass transition temperatures of amorphous trehalose-water mixtures and the mobility of water: an experimental and in silico study

Isothermal-isobaric molecular dynamics simulations are used to calculate the specific volume of models of trehalose and three amorphous trehalose-water mixtures (2.9%, 4.5% and 5.3% (w/w) water, respectively) as a function of temperature. Plots of specific volume versus temperature exhibit a characteristic change in slope when the amorphous systems change from the glassy to the rubbery state and the intersection of the two regression lines provides an estimate of the glass transition temperature T(g). A comparison of the calculated and experimental T(g) values, as obtained from differential scanning calorimetry, shows that despite the predicted values being systematically higher (about 21-26K), the trend and the incremental differences between the T(g) values have been computed correctly: T(g)(5.3%(w/w))相似文献   

Molecular modelling of the mass density of single proteins

Using molecular dynamics (MD) simulations, the density of single proteins and its temperature dependence was modelled starting from the experimentally determined protein structure and a generic, transferable force field, without the need of prior parameterization. Although all proteins consist of the same 20 amino acids, their density in aqueous solution varies up to 10% and the thermal expansion coefficient up to twofold. To model the protein density, systematic MD simulations were carried out for 10 proteins with a broad range of densities (1.32-1.43 g/cm(3)) and molecular weights (7-97 kDa). The simulated densities deviated by less than 1.4% from their experimental values that were available for four proteins. Further analyses of protein density showed that it can be essentially described as a consequence of amino acid composition. For five proteins, the density was simulated at different temperatures. The simulated thermal expansion coefficients ranged between 4.3 and 7.1 × 10(-4) K(-1) and were similar to the experimentally determined values of ribonuclease-A and lysozyme (deviations of 2.4 and 14.6%, respectively). Further analyses indicated that the thermal expansion coefficient is linked to the temperature dependence of atomic fluctuations: proteins with a high thermal expansion coefficient show a low increase in flexibility at increasing temperature. A low increase in atomic fluctuations with temperature has been previously described as a possible mechanism of thermostability. Thus, a high thermal expansion coefficient might contribute to protein thermostability.     

Application of isotropic periodic sum method for 4-pentyl-4′-cyanobiphenyl liquid crystal

In future large-scale molecular dynamics (MD) simulations that will use parallel computing, the isotropic periodic sum (IPS) method is expected to effectively reduce the cost of interaction calculations while maintaining adequate accuracy. To assess the accuracy of this method in estimating low-charge-density polymer systems, we performed atomistic MD simulations of the bulk state of liquid crystal systems based on 4-pentyl-4′-cyanobiphenyl (5CB). In conditions of 270 K ≤ T ≤ 320 K and a normal pressure, the temperature dependence of the density, potential energy and order parameter was estimated using the IPS and Ewald sum method. The results of the IPS method and Ewald sum were consistent within the range of error. In conditions close to the phase transition point, however, the averaged values of potential energy and order parameter had a small difference. We concluded that the fundamental physical properties for the bulk state of 5CB systems are determined reasonably by using the IPS method, at least in conditions that are not close to the phase transition point.     

Rheological characterization of the galactomannan from Leucaena leucocephala seed

A water soluble galactomannan isolated from Leucaena leucocephala seeds gave an intrinsic viscosity of 3.5dl/g and viscosity average molecular mass, M(v), of 6.98×10(5)g/mol. This was in reasonably good agreement with the value of the weight average molecular mass, M(w), of 5.44±0.20×10(5)g/mol determined by GPC-MALLS coupled to RI. The onset of polymer coil overlap occurred at c*[η] of 2.1, with slope of 3.0 above and 1.3 below the point of polymer coil overlap. The shear viscosity of the polysaccharide was temperature dependent and decreased with increasing temperature. The activation energy for viscous flow of 3.0% polysaccharide concentration obtained by Arrhenius plot of zero shear viscosity as a function of temperature was 26.4kJ/mol. Both the storage modulus (G') and loss modulus (G″) showed strong dependence on frequency indicating the presence of entangled coils. The Cox-Merz plot gave close superimposition of the complex and shear viscosities.     

Physical properties, crystallization, and spherulite growth of linear and 3-arm poly(L-lactide)s

The physical properties, crystallization, and spherulite growth behavior and mechanism of linear and 3-arm poly(L-lactide) [i.e., poly(L-lactic acid) (PLLA)] have been investigated using absolute molecular weight as a molecular index. The branching reduces the chain mobility of PLLA and must be excluded from the crystalline regions. The former factor gives the higher glass transition temperature (T(g)) and starting temperature for thermal degradation (T(d,S)) of 3-arm PLLA compared with those of linear PLLA. On the other hand, both the former and the latter factors lead to the higher cold crystallization temperature (T(cc)), the longer induction period for spherulite growth (t(i)), the lower melting temperature (T(m)), crystallinitiy (X(c)), and radius growth rate of the spherulties (G) for the 3-arm PLLA compared with those for the linear PLLA. The G of 3-arm PLLA showed the vague dependence on number-average molecular weight (M(n)), probably because the branching effect was balanced with the molecular weight effect. At the M(n) exceeding critical values, the linear and 3-arm PLLA crystallize in regime II or regime III kinetics, depending on crystallization temperature (T(c)). In contrast, at the M(n) below critical values, the linear and 3-arm PLLA crystallize according solely to regime III and regime II kinetics, respectively, for all the T(c).     

A test of improved force field parameters for urea: molecular-dynamics simulations of urea crystals

Molecular-dynamics (MD) simulations of urea crystals of different shapes (cubic, rectangular prismatic, and sheet) have been performed using our previously published force field for urea. This force field has been validated by calculating values for the cohesive energy, sublimation temperature, and melting point from the MD data. The cohesive energies computed from simulations of cubic and rectangular prismatic urea crystals in vacuo at 300 K agreed very well with the experimental sublimation enthalpies reported at 298 K. We also found very good agreement between the melting points as observed experimentally and from simulations. Annealing the crystals just below the melting point leads to reconstruction to form crystal faces that are consistent with experimental observations. The simulations reveal a melting mechanism that involves surface (corner/edge) melting well below the melting point, and rotational disordering of the urea molecules in the corner/edge regions of the crystal, which then facilitates the translational motion of these molecules.     

Low-temperature energy transfer in LHC-II trimers from the Chl a/b light-harvesting antenna of photosystem II.

Temperature dependence in electronic energy transfer steps within light-harvesting antenna trimers from photosystem II was investigated by studying Chl a pump-probe anisotropy decays at several wavelengths from 675 to 682 nm. The anisotropy lifetime is markedly sensitive to temperature at the longest wavelengths (680-682 nm), increasing by factors of 5 to 6 as the trimers are cooled from room temperature to 13 K. The temperature dependence is muted at 677 and 675 nm. This behavior is modeled using simulations of temperature-broadened Chl a absorption and fluorescence spectra in spectral overlap calculations of Förster energy transfer rates. In this model, the 680 nm anisotropy decays are dominated by uphill energy transfers from 680 nm Chl a pigments at the red edge of the LHC-II spectrum; the 675 nm anisotropy decays reflect a statistical average of uphill and downhill energy transfers from 676-nm pigments. The measured temperature dependence is consistent with essentially uncorrelated inhomogeneous broadening of donor and acceptor Chl a pigments.     

Spatial glass transition temperature variations in polymer glass: application to a maltodextrin-water system

A model was developed to predict spatial glass transition temperature (T(g)) distributions in glassy maltodextrin particles during transient moisture sorption. The simulation employed a numerical mass transfer model with a concentration dependent apparent diffusion coefficient (D(app)) measured using Dynamic Vapor Sorption. The mass average moisture content increase and the associated decrease in T(g) were successfully modeled over time. Large spatial T(g) variations were predicted in the particle, resulting in a temporary broadening of the T(g) region. Temperature modulated differential scanning calorimetry confirmed that the variation in T(g) in nonequilibrated samples was larger than in equilibrated samples. This experimental broadening was characterized by an almost doubling of the T(g) breadth compared to the start of the experiment. Upon reaching equilibrium, both the experimental and predicted T(g) breadth contracted back to their initial value.     

High critical temperature above T(g) may contribute to the stability of biological systems

In this study, we characterized the molecular mobility around T(g) in sugars, poly-L-lysine and dry desiccation-tolerant biological systems, using ST-EPR, (1)H-NMR, and FTIR spectroscopy, to understand the nature and composition of biological glasses. Two distinct changes in the temperature dependence of the rotational correlation time (tau(R)) of the spin probe 3-carboxy-proxyl or the second moment (M(2)) were measured in sugars and poly-L-lysine. With heating, the first change was associated with the melting of the glassy state (T(g)). The second change (T(c)), at which tau(R) abruptly decreased over several orders of magnitude, was found to correspond with the so-called cross-over temperature, where the dynamics changed from solid-like to liquid-like. The temperature interval between T(g) and T(c) increased in the order of sucrose < trehalose < raffinose 50 degrees C, implying that the stability above T(g) improved in the same order. These differences in temperature-dependent mobilities above T(g) suggest that proteins rather than sugars play an important role in the intracellular glass formation. The exceptionally high T(c) of intracellular glasses is expected to provide excellent long-term stability to dry organisms, maintaining a slow molecular motion in the cytoplasm even at temperatures far above T(g).     

A thermodynamic molecular switch in biological systems: ribonuclease S' fragment complementation reactions

It is well known that essentially all biological systems function over a very narrow temperature range. Most typical macromolecular interactions show DeltaH degrees (T) positive (unfavorable) and a positive DeltaS degrees (T) (favorable) at low temperature, because of a positive (DeltaCp degrees /T). Because DeltaG degrees (T) for biological systems shows a complicated behavior, wherein DeltaG degrees (T) changes from positive to negative, then reaches a negative value of maximum magnitude (favorable), and finally becomes positive as temperature increases, it is clear that a deeper-lying thermodynamic explanation is required. This communication demonstrates that the critical factor is a temperature-dependent DeltaCp degrees (T) (heat capacity change) of reaction that is positive at low temperature but switches to a negative value at a temperature well below the ambient range. Thus the thermodynamic molecular switch determines the behavior patterns of the Gibbs free energy change and hence a change in the equilibrium constant, K(eq), and/or spontaneity. The subsequent, mathematically predictable changes in DeltaH degrees (T), DeltaS degrees (T), DeltaW degrees (T), and DeltaG degrees (T) give rise to the classically observed behavior patterns in biological reactivity, as may be seen in ribonuclease S' fragment complementation reactions.     

Characterization of molecular mobility in seed tissues: an electron paramagnetic resonance spin probe study.

The relationship between molecular mobility (tauR) of the polar spin probe 3-carboxy-proxyl and water content and temperature was established in pea axes by electron paramagnetic resonance (EPR) and saturation transfer EPR. At room temperature, tauR increased during drying from 10(-11) s at 2.0 g water/g dry weight to 10(-4) s in the dry state. At water contents below 0.07 g water/g dry weight, tauR remained constant upon further drying. At the glass transition temperature, tauR was constant at approximately 10(-4) s for all water contents studied. Above Tg, isomobility lines were found that were approximately parallel to the Tg curve. The temperature dependence of tauR at all water contents studied followed Arrhenius behavior, with a break at Tg. Above Tg the activation energy for rotational motion was approximately 25 kJ/mol compared to 10 kJ/mol below Tg. The temperature dependence of tauR could also be described by the WLF equation, using constants deviating considerably from the universal constants. The temperature effect on tauR above Tg was much smaller in pea axes, as found previously for sugar and polymer glasses. Thus, although glasses are present in seeds, the melting of the glass by raising the temperature will cause only a moderate increase in molecular mobility in the cytoplasm as compared to a huge increase in amorphous sugars.     

Calorimetric study of crystal growth of ice in hydrated methemoglobin and of redistribution of the water clusters formed on melting the ice.

Calorimetric studies of the melting patterns of ice in hydrated methemoglobin powders containing between 0.43 and 0.58 (g water)/(g protein), and of their dependence on annealing at subzero temperatures and on isothermal treatment at ambient temperature are reported. Cooling rates were varied between approximately 1500 and 5 K min-1 and heating rate was 30 K min-1. Recrystallization of ice during annealing is observed at T > 228 K. The melting patterns of annealed samples are characteristically different from those of unannealed samples by the shifting of the melting temperature of the recrystallized ice fraction to higher temperatures toward the value of "bulk" ice. The "large" ice crystals formed during recrystallization melt on heating into "large" clusters of water whose redistribution and apparent equilibration is followed as a function of time and/or temperature by comparison with melting endotherms. We have also studied the effect of cooling rate on the melting pattern of ice with a methemoglobin sample containing 0.50 (g water)/(g protein), and we surmise that for this hydration cooling at rates of > or = approximately 150 K min-1 preserves on the whole the distribution of water molecules present at ambient temperature.     

The effect of gamma-irradiation on the temperature dependence of D.C. electrical conductivity of dry bone.

The effect of gamma-irradiation with doses from 10 to 500 kGy on the electrical conductivity (g) of dry bone was studied. Temperature measurement of electrical conductivity were made from 393 to 533 K. The dependence obtained indicates the increase in g with temperature. An increase in irradiation dose resulted in a decreased g value for each dose up to temperature 462 K. Temperature 462 K was interpreted as the temperature of collagen melting point in dry bone. Above 462 K, g values were dose independent. A dose of 500 kGy shifted the melting point to lower temperature. In addition, the activation energy for the charge conduction process was calculated. Obtained values for electrical conductivity and activation energy were typical for dielectrics and indicated degradation of the organic component of bone.     

The reversible two-state unfolding of a monocot mannose-binding lectin from garlic bulbs reveals the dominant role of the dimeric interface in its stabilization

Allium sativum agglutinin (ASAI) is a heterodimeric mannose-specific bulb lectin possessing two polypeptide chains of molecular mass 11.5 and 12.5 kDa. The thermal unfolding of ASAI, characterized by differential scanning calorimetry and circular dichroism, shows it to be highly reversible and can be defined as a two-state process in which the folded dimer is converted directly to the unfolded monomers (A2 if 2U). Its conformational stability has been determined as a function of temperature, GdnCl concentration, and pH using a combination of thermal and isothermal GdnCl-induced unfolding monitored by DSC, far-UV CD, and fluorescence, respectively. Analyses of these data yielded the heat capacity change upon unfolding (DeltaC(p) and also the temperature dependence of the thermodynamic parameters, namely, DeltaG, DeltaH, and DeltaS. The fit of the stability curve to the modified Gibbs-Helmholtz equation provides an estimate of the thermodynamic parameters DeltaH(g), DeltaS(g), and DeltaC(p) as 174.1 kcal x mol(-1), 0.512 kcal x mol(-1) x K(-1), and 3.41 kcal x mol(-1) x K(-1), respectively, at T(g) = 339.4 K. Also, the free energy of unfolding, DeltaG(s), at its temperature of maximum stability (T(s) = 293 K) is 13.13 kcal x mol(-1). Unlike most oligomeric proteins studied so far, the lectin shows excellent agreement between the experimentally determined DeltaC(p) (3.2 +/- 0.28 kcal x mol(-1) x K(-1)) and those evaluated from a calculation of its accessible surface area. This in turn suggests that the protein attains a completely unfolded state irrespective of the method of denaturation. The absence of any folding intermediates suggests the quaternary interactions to be the major contributor to the conformational stability of the protein, which correlates well with its X-ray structure. The small DeltaC(p) for the unfolding of ASAI reflects a relatively small, buried hydrophobic core in the folded dimeric protein.     

Structure prediction and functional analyses of a thermostable lipase obtained from Shewanella putrefaciens

Previous experimental studies on thermostable lipase from Shewanella putrefaciens suggested the maximum activity at higher temperatures, but with little information on its conformational profile. In this study, the three-dimensional structure of lipase was predicted and a 60 ns molecular dynamics (MD) simulation was carried out at temperatures ranging from 300 to 400 K to better understand its thermostable nature at the molecular level. MD simulations were performed in order to predict the optimal activity of thermostable lipase. The results suggested strong conformational temperature dependence. The thermostable lipase maintained its bio-active conformation at 350 K during the 60 ns MD simulations.     

Molecular dynamics simulation of the size-dependent morphological stability of cubic shape silver nanoparticles

The morphological stability of sharp-edged silver nanoparticles is examined by the classical molecular dynamics (MD) simulations. The crystalline structure and the perfect fcc atom packing of a series of silver nanocubes (AgNC) of different sizes varying from 63 up to 1099 atoms are compared against quasi-spherical nanoparticles of the same sizes at temperature 303 K. Our MD simulations demonstrate that starting from the preformed perfect crystalline structures the cubic shape is preserved for AgNCs composed of 365–1099 atoms. Surprisingly, the rapid loss of the cubic shape morphology and transformation into the non-fcc-structure are found for smaller AgNCs composed of less than ~256 atoms. No such loss of the preformed crystalline structure is seen for quasi-spherical nanoparticles composed of 38–1007 atoms. The analysis of the temperature dependence and the binding energy of outermost Ag surface atoms suggests that the loss of the perfect cubic shape, rounding and smoothing of sharp edges and corners are driven by the tendency towards the increase in their coordination number. In addition, we revealed that AgNC₁₀₉₉ partially loses its sharp edges and corners in the aqueous environment; however, the polymer coating with poly(vinyl alcohol) (PVA) was able to preserve the well-defined cubic morphology. Finally, these results help improve the understanding of the role of surface capping agents in solution phase synthesis of Ag nanocubes.     

Effect of temperature on elastic properties of CNT-polyethylene nanocomposite and its interface using MD simulations

This paper investigates the effect of temperature on the elastic modulus of carbon nanotube-polyethylene (CNT-PE) nanocomposite and its interface using molecular dynamics (MD) simulations, by utilizing the second-generation polymer consistent force field (PCFF). Two CNTs—armchair and zigzag—were selected as reinforcing nano-fillers, and amorphous PE was used as the polymer matrix. For atomistic modelling of the nanocomposite, the commercially available code Materials Studio 8.0 was used and all other MD simulations were subsequently performed using the open source code Large-Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS). To obtain the elastic modulus of the nanocomposite, stress-strain curves were drawn at different temperatures by performing uniaxial deformation tests on the nanocomposite material, whereas the curvatures of the interfacial interaction energy vs. strain curves were utilized to obtain Young’s modulus of the interface. In addition, the glass transition temperatures of the polymer matrix and nanocomposites were also evaluated using density-temperature curves. Based on the results, it is concluded that, irrespective of temperature condition, a nanocomposite reinforced with CNT of larger chirality (i.e., armchair) yields a higher value of Young’s modulus of the nanocomposite and its interface. It was also found that, at the phase transition (from a glassy to a rubbery state) temperature (i.e., glass transition temperature), Young’s moduli of the polymer matrix, nanocomposite, and its interface drop suddenly. The results obtained from MD simulations were verified with results obtained from continuum-based rule-of-mixtures.