Free energy approximations in simple lattice proteins

This work addresses the question of whether it is possible to define simple pairwise interaction terms to approximate free energies of proteins or polymers. Rather than ask how reliable a potential of mean force is, one can ask how reliable it could possibly be. In a two-dimensional, infinite lattice model system one can calculate exact free energies by exhaustive enumeration. A series of approximations were fitted to exact results to assess the feasibility and utility of pairwise free energy terms. Approximating the true free energy with pairwise interactions gives a poor fit with little transferability between systems of different size. Adding extra artificial terms to the approximation yields better fits, but does not improve the ability to generalize from one system size to another. Furthermore, one cannot distinguish folding from nonfolding sequences via the approximated free energies. Most usefully, the methodology shows how one can assess the utility of various terms in lattice protein/polymer models. (C) 2001 American Institute of Physics.

Mechanical alloying of MoSi2 with ternary alloying elements. Part 2 - computer simulation

A computer model of the mechanical alloying process has been developed to simulate phase formation during the mechanical alloying of Mo and Si elemental powders with a ternary addition of Al, Mg, Ti or Zr. Using the Arhennius equation, the model balances the formation rates of the competing reactions that are observed during milling. These reactions include the formation of tetragonal C11(b) MOSi2 (t-MoSi2) by combustion, the formation of the hexagonal C40 MoSi2 polymorph (h-MoSi2), the transformation of the tetragonal to the hexagonal form, and the recovery of t-MoSi2 from h-MoSi2 and deformed t-MoSi2. The addition of the ternary additions changes the free energy of formation of the associated MoSi2 alloys, i.e. Mo(Si, Al)(2), Mo(Mg, Al)(2), (Mo, Ti)Si-2 (Mo, Zr)Si-2 and (Mo, Fe)Si-2, respectively. Variation of the energy of formation alone is sufficient for the simulation to accurately model the observed phase formation. (C) 2003 Elsevier B.V. All rights reserved.

A methodology for building energy modelling and calibration in warm climates

In the last 7 years, a method has been developed to analyse building energy performance using computer simulation, in Brazil. The method combines analysis of building design plans and documentation, walk-through visits, electric and thermal measurements and the use of an energy simulation tool (DOE-2.1E code), The method was used to model more than 15 office buildings (more than 200 000 m(2)), located between 12.5degrees and 27.5degrees South latitude. The paper describes the basic methodology, with data for one building and presents additional results for other six cases. (C) 2002 Elsevier Science Ltd. All rights reserved.

Adsorption in slit-like pores of activated carbons: Improvement of the Horvath and Kawazoe method

Anew thermodynamic approach has been developed in this paper to analyze adsorption in slitlike pores. The equilibrium is described by two thermodynamic conditions: the Helmholtz free energy must be minimal, and the grand potential functional at that minimum must be negative. This approach has led to local isotherms that describe adsorption in the form of a single layer or two layers near the pore walls. In narrow pores local isotherms have one step that could be either very sharp but continuous or discontinuous benchlike for a definite range of pore width. The latter reflects a so-called 0 --> 1 monolayer transition. In relatively wide pores, local isotherms have two steps, of which the first step corresponds to the appearance of two layers near the pore walls, while the second step corresponds to the filling of the space between these layers. All features of local isotherms are in agreement with the results obtained from the density functional theory and Monte Carlo simulations. The approach is used for determining pore size distributions of carbon materials. We illustrate this with the benzene adsorption data on activated carbon at 20, 50, and 80 degreesC, argon adsorption on activated carbon Norit ROX at 87.3 K, and nitrogen adsorption on activated carbon Norit R1 at 77.3 K.

Tuning the metal-to-metal charge transfer energy of cyano-bridged dinuclear complexes

The metal-to-metal charge transfer (MMCT) transitions of a series of Class II mixed valence dinuclear complexes bearing cyano bridging ligands may be varied systematically by variations to either the hexacyanometallate(II) donor or Co-III acceptor moieties. Specifically, the new dinuclear species trans-[(LCoNCFe)-Co-14S(CN)(5)](-) (L-14S = 6-methyl-1,11-diaza-4,8-dithia- cyclotetradecane-6-amine) and trans-[(LCoNCRu)-Co-14(CN)(5)]-(L-14 = 6-methyl-1,4,8,11-tetraazacyclotetradecane-6-amine) have been prepared and their spectroscopic and electrochemical properties are compared with the relative trans-[(LCoNCFe)-Co-14(CN)(5)](-). The crystal structures of Na{trans-[(LCoNCFe)-Co-14S(CN)(5)]}.51/2H(2)O.1/2EtOH, Na{trans-[(LCoNCRu)-Co-14(CN)(5)]}.3H(2)O and Na{trans-[(LCoNCRu)-Co-14(CN)(5)]}.8H(2)O are also reported. The ensuing changes to the MMCT energy have been examined within the framework of Hush theory, and it was found that the free energy change between the redox isomers was the dominant effect in altering the energy of the MMCT transition.

Application of density functional theory to equilibrium adsorption of argon and nitrogen on amorphous silica surface

We present a new version of non-local density functional theory (NL-DFT) adapted to description of vapor adsorption isotherms on amorphous materials like non-porous silica. The novel feature of this approach is that it accounts for the roughness of adsorbent surface. The solid–fluid interaction is described in the same framework as in the case of fluid–fluid interactions, using the Weeks–Chandler–Andersen (WCA) scheme and the Carnahan–Starling (CS) equation for attractive and repulsive parts of the Helmholtz free energy, respectively. Application to nitrogen and argon adsorption isotherms on non-porous silica LiChrospher Si-1000 at their boiling points, recently published by Jaroniec and co-workers, has shown an excellent correlative ability of our approach over the complete range of pressures, which suggests that the surface roughness is mostly the reason for the observed behavior of adsorption isotherms. From the analysis of these data, we found that in the case of nitrogen adsorption short-range interactions between oxygen atoms on the silica surface and quadrupole of nitrogen molecules play an important role. The approach presented in this paper may be further used in quantitative analysis of adsorption and desorption isotherms in cylindrical pores such as MCM-41 and carbon nanotubes.

Adsorption of N-2, Ar, acetone, chloroform and acetone-chloroform mixture on carbon black in the framework of molecular layer structure theory (MLST)

Adsorption of pure nitrogen, argon, acetone, chloroform and acetone-chloroform mixture on graphitized thermal carbon black is considered at sub-critical conditions by means of molecular layer structure theory (MLST). In the present version of the MLST an adsorbed fluid is considered as a sequence of 2D molecular layers, whose Helmholtz free energies are obtained directly from the analysis of experimental adsorption isotherm of pure components. The interaction of the nearest layers is accounted for in the framework of mean field approximation. This approach allows quantitative correlating of experimental nitrogen and argon adsorption isotherm both in the monolayer region and in the range of multi-layer coverage up to 10 molecular layers. In the case of acetone and chloroform the approach also leads to excellent quantitative correlation of adsorption isotherms, while molecular approaches such as the non-local density functional theory (NLDFT) fail to describe those isotherms. We extend our new method to calculate the Helmholtz free energy of an adsorbed mixture using a simple mixing rule, and this allows us to predict mixture adsorption isotherms from pure component adsorption isotherms. The approach, which accounts for the difference in composition in different molecular layers, is tested against the experimental data of acetone-chloroform mixture (non-ideal mixture) adsorption on graphitized thermal carbon black at 50 degrees C. (C) 2005 Elsevier Ltd. All rights reserved.

Modeling of adsorption and nucleation in infinite cylindrical pores by two-dimensional density functional theory

In this paper, we present an analysis of argon adsorption in cylindrical pores having amorphous silica structure by means of a nonlocal density functional theory (NLDFT). In the modeling, we account for the radial and longitudinal density distributions, which allow us to consider the interface between the liquidlike and vaporlike fluids separated by a hemispherical meniscus in the canonical ensemble. The Helmholtz free energy of the meniscus was determined as a function of pore diameter. The canonical NLDFT simulations show the details of density rearrangement at the vaporlike and liquidlike spinodal points. The limits of stability of the smallest bridge and the smallest bubble were also determined with the canonical NLDFT. The energy of nucleation as a function of the bulk pressure and the pore diameter was determined with the grand canonical NLDFT using an additional external potential field. It was shown that the experimentally observed reversibility of argon adsorption isotherms at its boiling point up to the pore diameter of 4 nm is possible if the potential barrier of 22kT is overcome due to density fluctuations.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Effect of adsorption deformation on thermodynamic characteristics of a fluid in slit pores at sub-critical conditions

Solvation. pressure due to adsorption of fluids in porous materials is the cause of elastic deformation of an adsorbent, which is accessible to direct experimental measurements. Such a deformation contributes to the Helmholtz free energy of the whole adsorbent-adsorbate system due to accumulation of compression or tension energy by the solid. It means that in the general case the solid has to be considered as not solely a source of the external potential field for the fluid confined in the pore volume, but also as thermodynamically nonmert component of the solid-fluid system. We present analysis of nitrogen adsorption isotherms and heat of adsorption in slit graphitic pores accounting for the adsorption deformation by means of nonlocal density functional theory. (c) 2006 Elsevier Ltd. All rights reserved.

Equation of state of a superfluid Fermi gas in the BCS-BEC crossover

We present a theory for a superfluid Fermi gas near the BCS-BEC crossover, including pairing fluctuation contributions to the free energy similar to that considered by Nozieres and Schmitt-Rink for the normal phase. In the strong coupling limit, our theory is able to recover the Bogoliubov theory of a weakly interacting Bose gas with a molecular scattering length very close to the known exact result. We compare our results with recent Quantum Monte Carlo simulations both for the ground state and at finite temperature. Excellent agreement is found for all interaction strengths where simulation results are available.

In vivo protein cyclization promoted by a circularly permuted Synechocystis sp PCC6803 DnaB mini-intein

A synthetic Synechocystis sp. PCC6803 DnaB split mini-intein gene was constructed for the in vivo cyclization of recombinant proteins expressed in Escherichia coli. The system was used to cyclize the NH2-terminal domain of E. coli DnaB, the structure of which had been determined previously by NMR spectroscopy. Cyclization was found to proceed efficiently, with little accumulation of precursor, and the product was purified in high yield. The solution structure of cyclic DnaB-N is not significantly different from that of linear DnaB-N and it unfolds reversibly at temperatures similar to14 degreesC higher. Improved hydrogen bonding was observed in the first and last helices, and the length of the last helix was increased, while the 9-amino acid linker used to join the NH2 and COOH termini was found to be highly mobile. The measured thermodynamic stabilization of the structure (DeltaDeltaG approximate to 2 kcal/mol) agrees well with the value estimated from the reduced conformational entropy in the unfolded form. Simple polymer theory can be used to predict likely free energy changes resulting from protein cyclization and how the stabilization depends on the size of the protein and the length of the linker used to connect the termini.

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10degreesC by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG = 2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding I measured by stopped-flow circular dichroism. rate of less than 4 min(-1) The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E. coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding. (C) 2005 Elsevier Ltd. All rights reserved.

Effects of molecular crowding by saccharides on alpha-chymotrypsin dimerization

Given the importance of protein complexes as therapeutic targets, it is necessary to understand the physical chemistry of these interactions under the crowded conditions that exist in cells. We have used sedimentation equilibrium to quantify the enhancement of the reversible homodimerization of alpha-chymotrypsin by high concentrations of the osmolytes glucose, sucrose, and raffinose. In an attempt to rationalize the ostuolyte-mediated stabilization of the a-chymotrypsin homodimer, we have used models based on binding interactions (transfer-free energy analysis) and steric interactions (excluded volume theory) to predict the stabilization. Although transfer-free energy analysis predicts reasonably well the relatively small stabilization observed for complex formation between cytochrome c and cytochrome c peroxidase, as well as that between bobtail quail lysozyme and a monoclonal Fab fragment, it underestimates the sugar-mediated stabilization of the alpha-chymotrypsin dimer. Although predictions based on excluded volume theory overestimate the stabilization, it would seem that a major determinant in the observed stabilization of the a-chymotrypsin homodimer is the thermodynamic nonideality arising from molecular crowding by the three small sugars.

Factors affecting the rate of phosphocreatine resynthesis following intense exercise

Within the skeletal muscle cell at the onset of muscular contraction, phosphocreatine (PCr) represents the most immediate reserve for the rephosphorylation of adenosine triphosphate (ATP). As a result, its concentration can be reduced to less than 30% of resting levels during intense exercise. As a fall in the level of PCr appears to adversely affect muscle contraction, and therefore power output in a subsequent bout, maximising the rate of PCr resynthesis during a brief recovery period will be of benefit to an athlete involved in activities which demand intermittent exercise. Although this resynthesis process simply involves the rephosphorylation of creatine by aerobically produced ATP (with the release of protons), it has both a fast and slow component, each proceeding at a rate that is controlled by different components of the creatine kinase equilibrium. The initial fast phase appears to proceed at a rate independent of muscle pH. Instead, its rate appears to be controlled by adenosine diphosphate (ADP) levels; either directly through its free cytosolic concentration, or indirectly, through its effect on the free energy of ATP hydrolysis. Once this fast phase of recovery is complete, there is a secondary slower phase that appears almost certainly rate-dependant on the return of the muscle cell to homeostatic intracellular pH. Given the importance of oxidative phosphorylation in this resynthesis process, those individuals with an elevated aerobic power should be able to resynthesise PCr at a more rapid rate than their sedentary counterparts. However, results from studies that have used phosphorus nuclear magnetic resonance (P-31-NMR) spectroscopy, have been somewhat inconsistent with respect to the relationship between aerobic power and PCr recovery following intense exercise. Because of the methodological constraints that appear to have limited a number of these studies, further research in this area is warranted.