Temperature dependent free energy surface of polymer folding from equilibrium and quench studies

Langevin dynamics simulation studies have been employed to calculate the temperature dependent free energy surface and folding characteristics of a 500 monomer long linear alkane (polyethylene) chain with a realistic interaction potential. Both equilibrium and temperature quench simulation studies have been carried out. Using the shape anisotropy parameter (S) of the folded molecule as the order parameter, we find a weakly first order phase transition between the high-temperature molten globule and low-temperature rodlike crystalline states separated by a small barrier of the order of k(B)T. Near the melting temperature (580 K), we observe an intriguing intermittent fluctuation with pronounced ``1/f noise characteristics'' between these two states with large difference in shape and structure. We have also studied the possibilities of different pathways of folding to states much below the melting point. At 300 K starting from the all-trans linear configuration, the chain folds stepwise into a very regular fourfold crystallite with very high shape anisotropy. Whereas, when quenched from a high temperature (900 K) random coil regime, we identify a two step transition from the random coiled state to a molten globulelike state and, further, to a anisotropic rodlike state. The trajectory reveals an interesting coupling between the two order parameters, namely, radius of gyration (R-g) and the shape anisotropy parameter (S). The rodlike final state of the quench trajectory is characterized by lower shape anisotropy parameter and significantly larger number of gauche defects as compared to the final state obtained through equilibrium simulation starting from all-trans linear chain. The quench study shows indication of a nucleationlike pathway from the molten globule to the rodlike state involving an underlying rugged energy landscape. (C) 2010 American Institute of Physics. doi:10.1063/1.3509398]

Critical phenomena in polymer solutions: Scaling of the free energy

The thermodynamics of monodisperse solutions of polymers in the neighborhood of the phase separation temperature is studied by means of Wilson’s recursion relation approach, starting from an effective ϕ4 Hamiltonian derived from a continuum model of a many‐chain system in poor solvents. Details of the chain statistics are contained in the coefficients of the field variables ϕ, so that the parameter space of the Hamiltonian includes the temperature, coupling constant, molecular weight, and excluded volume interaction. The recursion relations are solved under a series of simplifying assumptions, providing the scaling forms of the relevant parameters, which are then used to determine the scaling form of the free energy. The free energy, in turn, is used to calculate the other singular thermodynamic properties of the solution. These are characteristically power laws in the reduced temperature and molecular weight, with the temperature exponents being the same as those of the 3d Ising model. The molecular weight exponents are unique to polymer solutions, and the calculated values compare well with the available experimental data.

A critical assessment of the standard molar gibbs free energy of formation of NiWO4

Three independent studies have been reported on the free energy of formation of NiWO4. Results of these measurements are analyzed by the �third-law� method, using thermal functions for NiWO4 derived from both low and high temperature heat capacity measurements. Values for the standard molar enthalpy of formation of NiWO4 at 298·15 K obtained from �third-law� analysis are compared with direct calorimetric determinations. Only one set of free energy measurements is found to be compatible with calorimetric enthalpies of formation. The selected value for ?f H m 0 (NiWO4, cr, 298·15 K) is the average of the three calorimetric measurements, using both high temperature solution and combustion techniques, and the compatible free energy determination. A new set of evaluated data for NiWO4 is presented.

Solubilization of n-decane into n-decyl-polyoxyethylene surfactants: A cavity forming free energy based model

The present paper reports the results of a theoretical study of the forces and factors driving the solubilization of n-alkane solubilizates into the micellar core of some non-ionic surfactants, based on a micellar model which includes the cavity forming free energy as a component of micellization. The solubilizate is n-decane and the non-ionic surfactants considered are n-decyl-polyoxyethylene surfactants. The extent of solubilization, i.e. the mole fraction of the solubilizate within the core has been calculated. The results indicate that the incorporated solubilizate has more translational and rotational degrees of freedom as compared to those of the tail parts of the surfactants present in the core. This drives the total free energy of aggregation after solubilization into a more favourable direction. The results are in fair agreement with the experimental results.

Time scales for transitions between free-energy minima of a dense hard-sphere system

Time scales associated with activated transitions between glassy metastable states of a free-energy functional appropriate for a dense hard-sphere system are calculated by using a new Monte Carlo method for the local density variables. In particular, we calculate the time the system, initially placed in a shallow glassy minimum of the free-energy, spends in the neighborhood of this minimum before making a transition to the basin of attraction of another free-energy minimum. This time scale is found to increase as the average density is increased. We find a crossover density near which this time scale increases very sharply and becomes longer than the longest times accessible in our simulation. This time scale does not show any evidence of increasing with sample size

DNA Free Energy-Based Promoter Prediction and Comparative Analysis of Arabidopsis and Rice Genomes

The cis-regulatory regions on DNA serve as binding sites for proteins such as transcription factors and RNA polymerase. The combinatorial interaction of these proteins plays a crucial role in transcription initiation, which is an important point of control in the regulation of gene expression. We present here an analysis of the performance of an in silico method for predicting cis-regulatory regions in the plant genomes of Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) on the basis of free energy of DNA melting. For protein-coding genes, we achieve recall and precision of 96% and 42% for Arabidopsis and 97% and 31% for rice, respectively. For noncoding RNA genes, the program gives recall and precision of 94% and 75% for Arabidopsis and 95% and 90% for rice, respectively. Moreover, 96% of the false-positive predictions were located in noncoding regions of primary transcripts, out of which 20% were found in the first intron alone, indicating possible regulatory roles. The predictions for orthologous genes from the two genomes showed a good correlation with respect to prediction scores and promoter organization. Comparison of our results with an existing program for promoter prediction in plant genomes indicates that our method shows improved prediction capability.

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Estimation of surface free energy of pyrites by contact angle measurements

The surface properties of coal-pyrite play a major role in determining its separation from coal in processes such as flotation. The solution pH is an important parameter in determining the surface properties of both coal and coal-pyrite such as surface free energy and zeta-potential. In the present investigation, the effect of pH on the surface free energy of pyrites from different sources was studied. The surface free energy of solids is made up of two components, i.e. the dispersive surface free energy and the acid-base interaction energy. Various methods have been used by previous researchers to evaluate these two components for different solids. In the present study, a new approach was developed and used to study the surface free energy of pyrite surfaces. Results indicate that the dispersion surface free energy of various pyrites is independent of pH while the acid-base interaction energy is strongly dependent on the pH. The acid-base interaction energy is different for each pyrite sample and also the change with pH varies with the type of pyrite. Coal-pyrite was found to be more hydrophobic than ore-pyrite which may be attributed to the presence of carbon in coal-pyrites. The acid-base interaction energy varied little with pH for coal pyrites than ore-pyrite. Comparison of acid-base interaction energy with zeta-potential measurements shows a good correlation between the minimum in acid-base interaction energy and the pHpzc.

Understanding glassy dynamics from the free-energy landscape

Many interesting features of the dynamics of simple liquids near the glass transition may be understood in terms of properties of the free-energy landscape obtained from numerical studies of a model free-energy functional. Main results obtained from this approach are summarized and a list of references to relevant publications is provided. (C) 2002 Elsevier Science B.V. All rights reserved.

Docking and free energy simulations to predict conformational domains involved in hCG-LH receptor interactions using recombinant antibodies

Single chain fragment variables (ScFvs) have been extensively employed in studying the protein-protein interactions. ScFvs derived from phage display libraries have an additional advantage of being generated against a native antigen, circumventing loss of information on conformational epitopes. In the present study, an attempt has been made to elucidate human chorionic gonadotropin (hCG)-luteinizing hormone (LH) receptor interactions by using a neutral and two inhibitory ScFvs against hCG. The objective was to dock a computationally derived model of these ScFvs onto the crystal structure of hCG and understand the differential roles of the mapped epitopes in hCG-LH receptor interactions. An anti-hCG ScFv, whose epitope was mapped previously using biochemical tools, served as the positive control for assessing the quality of docking analysis. To evaluate the role of specific side chains at the hCG-ScFv interface, binding free energy as well as residue interaction energies of complexes in solution were calculated using molecular mechanics Poisson-Boltzmann/surface area method after performing the molecular dynamic simulations on the selected hCG-ScFv models and validated using biochemical and SPR analysis. The robustness of these calculations was demonstrated by comparing the theoretically determined binding energies with the experimentally obtained kinetic parameters for hCG-ScFv complexes. Superimposition of hCG-ScFv model onto a model of hCG complexed with the 51-266 residues of LH receptor revealed importance of the residues previously thought to be unimportant for hormone binding and response. This analysis provides an alternate tool for understanding the structure-function analysis of ligand-receptor interactions. Proteins 2011;79:3108-3122. (C) 2011 Wiley-Liss, Inc.

Potentiometric determination of the Gibbs free energy of formation of cadmium and magnesium chromites

The standard Gibbs free energy of formation of magnesium and cadmiumchromites have been determined by potentiometric measurements on reversiblesolid-state electrochemical cells [dformula (Au-5%Cd, , Au-5%Cd; Pt, + , CaO-ZrO[sub 2], + ,Pt; CdO, , CdCr[sub 2]O[sub 4] + Cr[sub 2]O[sub 3])] in the temperature range 500°–730°C, and [dformula Pt, Cr + Cr[sub 2]O[sub 3]/Y[sub 2]O[sub 3]-ThO[sub 2]/Cr + MgCr[sub 2]O[sub 4] + MgO, Pt] in the temperature range 800°–1200°C. The temperature dependence of the freeenergies of formation of the ternary compounds can be represented by theequations [dformula CdO(r.s.) + Cr[sub 2]O[sub 3](cor) --> CdCr[sub 2]O[sub 4](sp)] [dformula Delta G[sup 0] = - 42,260 + 7.53T ([plus-minus]400) J] and [dformula MgO(r.s.) + Cr[sub 2]O[sub 3](cor) --> MgCr[sub 2]O[sub 4](sp)] [dformula Delta G[sup 0] = - 45,200 + 5.36T ([plus-minus]400) J] The entropies of formation of these spinels are discussed in terms of cationdisorder and extent of reduction of Cr3+ ions to Cr2+ ions. Thermodynamicdata on the chromates of cadmium and magnesium are derived by combiningthe results obtained in this study with information available in the literatureon high temperature, high pressure phase equilibria in the systems CdO-Cr2O3-O2 and MgO-Cr2O3-O2.

Gibbs free energy of formation of magnesium stannate

he thermodynamic properties of the spinel Mg2SnO4 have been determined by emf measurements on the solid oxide galvanic cell,View the MathML source in the temperature range 600 to 1000°C. The Gibbs' free energy of formation of Mg2SnO4 from the component oxides can be expressed as View the MathML source,View the MathML source These values are in good agreement with the information obtained by Jackson et al. [Earth Planet. Sci. Lett.24, 203 (1974)] on the high pressure decomposition of magnesium stannate into component oxides at different temperatures. The thermodynamic data suggest that the spinel phase is entropy stabilized, and would be unstable below 207 (±25)°C at atmospheric pressure. Based on the information obtained in this study and trends in the stability of aluminate and chromite spinels, it can be deduced that the stannates of nickel and copper(II) are unstable.

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.