Exact Solution to the Extended Zwanzig Model for Quasi-Sigmoidal Chemically Induced Denaturation Profiles: Specific Heat and Configurational Entropy

Temperature and chemically induced denaturation comprise two of themost characteristicmechanisms to achieve the passage from the native stateN to any of the unstructured statesDj in the denatured ensemble in proteins and peptides. In this work we present a full analytical solution for the configurational partition function Zqs of a homopolymer chain poly-X in the extended Zwanzig model (EZM) for a quasisigmoidal denaturation profile. This solution is built up from an EZM exact solution in the case where the fraction α of native contacts follows exact linear dependence on denaturant’s concentration ζ; thus an analytical solution for ZL in the case of an exact linear denaturation profile is also provided. A recently established connection between the number ] of potential nonnative conformations per residue and temperature-independent helical propensityω complements themodel in order to identify specific proteinogenic poly-X chains, where X represents any of the twenty naturally occurring aminoacid residues. From Zqs, equilibrium thermodynamic potentials like entropy S and average internal energy ⟨E⟩ and thermodynamic susceptibilities like specific heat Cv are calculated for poly-valine (poly-V) and poly-alanine (poly-A) chains. The influence of the rate at which native contacts denature as function of ζ on thermodynamic stability is also discussed.


Introduction
The early recognition that ordered three-dimensional macromolecular structures, later known as the native state, play a fundamental role on biological activity of proteins and peptides led to a large interest in the fundamental mechanisms that drive formation and stabilization of such structures [1][2][3][4].According to the classical view of protein folding, the native state  is a unique three-dimensional structure in which the protein displays biological activity [1,5].More recent studies, however, have shown that the structure-function paradigm may have remarkable exceptions.A large class of proteins known as intrinsically disordered proteins (IDPs) is one of the most outstanding examples of this category.IDPs can show a high level of biological activity even if the protein is in one of the   states associated with the disordered ensemble [6][7][8][9].The transition from  →   has been, accordingly, a subject of intense research.Temperature and chemically induced denaturation are two of the most common factors to induce such transitions being, the latter, probably the second most studied denaturation mechanism [10].In order to chemically induce the  →   transition, urea and guanidine hydrochloride (GdnHCl) are the most usual organic compounds [11][12][13].
Despite the large amount of theoretical and experimental work, there is still an open debate on the underlying physical mechanism driving urea and GdnHCl mediated denaturation [14][15][16].It is widely known that solvation of proteins and peptides in urea or GdnHCl decreases the free energy Δ  →  2 International Journal of Statistical Mechanics compared to free energy of solvation in water Δ ∘  →  and that such decrease is linearly related to denaturant's concentration  according to Δ  →  = Δ ∘  →  − , where  is the so-called  value [17][18][19].From theoretical or computational evaluation of  values it is expected to assess if the decrease in the free energy of solvation Δ  →  in the protein-denaturant system is due to electrostatic, hydrogen bonding, or van der Waals interactions [10,11].A second central debate refers to whether denaturant molecules act directly on protein side-chains and backbone or instead act on bulk water altering hydrogen bond network and thus lessening or screening water-protein interactions [20,21].MD simulations suggest that whereas GdnHCl seems to decrease solvation free energy by electrostatic forces, urea seems to act through van der Waals forces in the first solvation shell of protein [13,22].
From a theoretical point of view, statistical mechanics coarse-grained (SMCG) models have gained increasing popularity in order to carry out a physical-based analysis of both protein folding and protein denaturation [23][24][25][26].The central feature of SMCG methods is the use of a small set of coarse-grained variables to describe average properties of the system instead of a fully detailed atom description.Common examples of such sets of variables include, but are not restricted to, pairs of dihedral angles (  ,   ) for each residue, number of nonnative conformations ] per residue, or intrinsic helical propensity  among others [27][28][29].Such set of parameters are then used to perform statistical averages over specific ensembles following equilibrium statistical mechanics prescriptions [25,30].A pioneer example of this category of methods is the Lifson-Roig treatment of the helixcoil transition [29] in peptides in which the configurational partition function  is written in terms of a transfer matrix T which depends on the conditional probabilities (, V, ) of having a particular residue either in a coil or helical state.
Two additional recent examples of SMCG models are the elastic network model (EN) [31] and the Zwanzig model (ZW) [28].The EN model considers a polypeptide as a chain of beads representing   carbons interacting to closest neighbors via harmonic potentials of varying strengths [26,31,32].Whereas in the original version a cut-off distance  0 is considered beyond of which residues do not interact at all [31], other versions of the EN model assign a decreasing value to the elastic constant   as the distance from a reference point is increased.Albeit EN model is suitable to analyze vibrational properties of equilibrium states, the ZW model on the other hand is more appropriate to discuss equilibrium thermodynamic properties of homopolymer chains [33].In the ZW model, an ansatz is proposed a priori for the energy levels   accessible to the system thus avoiding a quantum mechanical calculation of the energy spectra of the polypeptide chain through Schrodinger's equation ĤΨ  =   Ψ  [28,34].Such energy spectra {  } are minimally frustrated with a single ground state  0 and multiple metastable states   with  ̸ = 0 [35].In addition, a coarse-grained variable ] is used to describe the number of potential nonnative conformations per residue [28].
The original ZW model considers temperature as the only driving force to induce denaturation.However, it is possible to extend it to consider the concentration  of a chemical denaturant as a contributing factor to unfolding [34].This extended Zwanzig model (EZM) uses the functional dependence of the fraction   of native contacts with respect to  to derive an analytical formula for the configurational partition function Z of the polypeptide chain in terms of an integral over microstates whose energy   depends now on , ], , and .The fact that the experimental denaturation curve   () usually follows a sigmoidal pattern hinders the possibility of a closed analytical solution for such integral [34].However, the EZM admits a full analytical solution for the case in which  depends linearly on  [36].Although proteins and peptides do not follow in most situations this type of denaturation curve under chemical denaturants, it is possible to use this zerothorder solution to build up a general solution for the sigmoidal case.
In this work we use the EZM framework to analyze thermodynamic equilibrium properties of ideal poly-X chains, where X represents any of the twenty naturally occurring aminoacid residues, under the influence of a chemical denaturation profile.In particular we focus our attention on thermodynamic susceptibilities like specific heat  V (, , ) and thermodynamic equilibrium potentials as configurational entropy S(, , ) and their connection to the mechanism of chemical denaturation imposed on the chain.In order to achieve this goal, we provide first an analytical solution for the configurational partition function Z  in the EZM model for the case in which the fraction   of native contacts is a linear function of concentration of chemical denaturant.From this result, we make an approximation for the full sigmoidal denaturation curve as a superposition of linear denaturation profiles.This quasisigmoidal representation provides us with a full analytical solution for the configurational partition function Z in EZM from which thermodynamic functions can be obtained straightforwardly.In addition, two types of denaturation mechanisms called A and B can be identified for this quasisigmoidal denaturation profile depending on the rate at which the fraction  of native structure is broken down by the chemical compound.Using a recently derived connection between helical propensity  and the number ] of potential nonnative conformations per residue [37], an association with specific poly-X chains can be accomplished.Thermodynamics for poly-alanine (poly-A) and poly-valine (poly-V) in terms of the specifics of quasisigmoidal denaturation profile are discussed.
We proceed as follows.Section 2 contains a brief and selfcontained summary of the physical foundations of the EZM.In Section 3 an analytical solution for the configurational partition function Z  in the EZM for the linear   =   () is presented.Definition and properties of a quasisigmoidal denaturation profile are defined in Section 4 together with the implementation of the linear solution to obtain a full analytical solution for the quasisigmoidal   () case.For this case, Z  is obtained in terms of standard parameters of the EZM and also in terms of the slopes   of the individual linear profiles.As we shall see, such slopes weight the different contributions of each linear denaturation profile to the full partition function.Finally, Section 5 presents the behavior of  V and configurational entropy S for two poly-X homopolymer chains with different helical propensity.A further discussion is included on how thermodynamic equilibrium properties depend on denaturation curve.

Chemically Induced Protein and
Peptide Denaturation: The Extended Zwanzig Model (EZM) Statistical mechanics based coarse-grained models (SMCG models) provide a powerful tool to unveil intrinsic chain contributions to macroscopic thermodynamic observables of protein and peptide in solution [23][24][25]38].Depending on the details included, analytically solvable SMCG models may provide the fine tuning necessary to dissect the relative contributions of each major driving force for folding or unfolding.Even they are not so copious as expected, the ZW model [28,34,37] for ideal homopolymer chains and the EN model [26] together with on-lattice models with Ising-like interactions provide some outstanding examples of this category [39,40] as mentioned before.In particular, ZW model proposes an ansatz for the energy spectrum {  } of the chain which depends essentially on the number  of residues participating in a native contact.With this ansatz, ZW model overcomes the central obstacle encountered usually in applying the methodology of equilibrium statistical mechanics to calculate thermodynamic properties, which is the finding of the energy levels { 0 ,  1 , . . ., } accessible to the system consistent with macroscopic restrictions and physical interactions encoded in the Hamiltonian operator Ĥ in quantum mechanics.Given this fact, the proposal of a phenomenological energy spectrum {  } for a chain endowed with a single ground energy state and multiple metastable states satisfies the basic expectations for minimally frustrated energy landscape in proteins and peptides [41].For a chain of  residues, the ZW model proposes that the energy spectrum  can be written as where  =   is the number of residues in a native contact,  > 0 is the energy penalty for a nonnative contact,  > 0 is the overall energy of the ground state, and  , is the Kronecker delta function.It is worth underlining that here there is only one ground state which is associated with a native state and also that folding is a cooperative effect.Only when all residues have a native contact, namely, when  = , the second term of (1) contributes.Such functional form for  implies that there is a single ground state with energy − and a set of unstable states with energies , 2, . . ., .If we assume that all residue's conformations are statistically independent and there is a number ]+1 for each one, the partition function  can be calculated readily to give where   is the degeneracy associated with energy level   and  = () −1 is the Bolztmann factor as usual.From (2) equilibrium thermodynamic functions like free energy , entropy  and specific heat  V can be calculated [21].
In its original formulation, ZW model considers temperature as the central factor driving thermodynamic phase transitions.Although this is mostly true for a large variety of protein and peptide conformational transitions, there are other physicochemical factors that may trigger similar conformational phase transition.This is particularly the case of chemical denaturation which consists of the addition of chemical compounds, usually guanidine chloride (GdnCl) or urea, to an equilibrium water-protein solution to destabilize native structure.As discussed previously, details of the physical mechanisms by which GdnHCl and urea destabilize native structure are not fully understood; however the basics have been already devised [10,14,42].For instance, for the case of urea, two central mechanisms have been proposed, a direct one in which a direct substitution of protein-protein and protein-water contacts is achieved by urea and an indirect one in which urea affects primarily the water hydrogen bond network of the first water shell around protein, thus making the exposure of nonpolar residues to protein's surface more favorable [15].Recently it has been calculated that the passage of a urea molecule from the bulk to the first shell water that immediately surrounds protein is more favorable energetically Δ = −12 kJ cal −1 than the transfer of a water molecule from the bulk to the immediate vicinity Δ = −0.25 kcal mol −1 .Thus the hypothesis of an indirect action of urea in protein denaturation is supported.The most probable picture for the action of urea is a two-step process.In a first stage formation of urea-water (U-W) contacts is achieved in the first water shell.Subsequently, when urea concentration has reached a value of approximately  = 8 Molar, formation of protein-urea (P-U) is favored over protein-protein (P-P) contacts [14].
It is possible to propose a SMCG model to include the presence of a chemical denaturating agent as a major factor in protein denaturation.To achieve this an extension of the original ZW model was proposed in the recent past [34].The corner stone of this extended Zwanzig model (EZM) considers that the action of a denaturating agent can be directly measured by its action on the fraction   of native contacts along the chain, which has to be an explicit function of the concentration  of the chemical compound.Then, the EZM assumes that   =   () (or its equivalent  = ()) which is a functional form that can be deduced experimentally.For instance, for a pure  helix-coil or  sheet-coil transition, the fraction   of residues in a native conformation at a fixed temperature can be measured by circular dichroism (CD).For most chemical and temperature-induced conformational transitions,   follows a sigmoidal relationship with either temperature or chemical concentration.Hence, for a sigmoidal denaturation profile,   () can be written as where , , and  are parameters that have to be adjusted in order to fit a specific profile.Given the fact that the energy  International Journal of Statistical Mechanics of a particular conformation is determined by the number  of native contacts and that  = (), the energy spectrum in the EZM is a function of the continuous variable .Thus the EZM partition function Z can be written as [34] where   is the degeneracy associated with the energy level () and  * is the concentration with a threshold value.Equation (4) represents the continuous version of ( Passage from a discrete model for  as in the ZW model to a continuous model for Z in the EZM is not just a simple matter of substitution of sums by integrals [34].Special care has to be taken since all variables that depend on the discrete variable  have to be correctly translated in terms of the continuous variable .For instance, the degeneracy factor   in the ZW model can be written as with (   ) being the standard binomial coefficient.However, if  is a continuous function of  as it is stated in (3) the correct translation for the degeneracy factor   in (4) implies constructing a continuous approximation to the binomial coefficient.In the same spirit, the discrete Kronecker delta function  , that represents the ground energy state in (1) and that mathematically is a well of energy can be best represented in the continuous case by -sequences   ().A particular useful -sequence   in one-dimension that describes a finite and continuous energy well is where  is any positive integer.By definition, -sequences satisfy that lim Thus, for sufficiently large ,   () is a strong peaked function [43].
The energy spectrum  as a function of denaturant's concentration is then obtained in two steps.Using remarks made above together with (3) one obtains where parameters , , and  of (3) have to be determined to fit the mathematical properties of a particular sigmoidal profile (see [34]); here  1/2 is the concentration where denaturation has reached its midpoint, that is, where ( 1/2 ) = /2 and  ⋆ is defined as  ⋆ = /√ in order to fulfill integration requirements of -sequences with  the ground state energy of the native state in the discrete version.As stated before,  is an arbitrary positive integer that can be adjusted accordingly.The partition function Z in EZM is thus obtained by substitution of ( 8) into (4).As it can be easily seen Z's integrand is a highly nonlinear function on .Then, a full analytical solution is not attainable.

Full Analytical Solution to EZM for Linear Denaturation Profiles
It is possible to develop an analytical approximation to Z taking into account that, for some values of , ,  1/2 , and ], the integrand of ( 4) has a well-defined peak at some value of 0 ⩽  ⩽  ⋆ where  ⋆ is the critical value at which secondary structure is fully unfolded.For these cases, integral in (4) can be approximated by a highly peaked single value and a zerothorder solution for Z can be obtained.Nevertheless, this approach has a strong limitation.This zeroth-order solution Z 0 retains a complex and rich dependence on ],  1/2 ,  and some other quantities of interest.However, temperature dependence essentially has been dropped out leaving only a first-order temperature on .This yields that thermodynamic functions like free energy  or specific heat  V shows no particularly interesting behavior on temperature.In fact since  V is a second derivative of Z 0 with respect to , this zerothorder approximation yields that  V = 0 giving then no further detail on the existence of thermodynamic phase transitions [34].EZM admits a full analytical solution for Z in the case of a linear denaturation profile.For such case, the number  of residues with native contacts follows a linear relationship with .This is where  and  are adjustable constants to consider specific denaturation profiles.The general tendency described in (9) dictates that the rate at which urea or GdnHCl denatures protein is constant.Fully denatured systems are reached at a particular concentration  ⋆ .Using ( 9) and ( 1) for the functional form of  as function of concentration, it is possible to obtain  as function of  as The combinatorial part of the degeneracy   = ] − (   ) can be written in terms of continuous variable  by means of a gaussian function centered at /2 which is the value where the binomial term exhibits its maximum value.This is Then, the partition function Z  in the case of a linear denaturation profile can be written as where the first term corresponds to  = 0 and the integral is performed over an infinitesimal quantity of denaturant  and afterwards a limit when  → 0 is performed.Second term in (12) can be written as If we introduce the change of variable we can write the following equivalence: Hence, Z  in the EZM for a linear denaturation profile reads as with the integral limits given by Since the integral ∫    − 2  admits a closed solution given by where Erf[] is the error function at , the function Z  can be finally written as where  = () −1 as usual is Boltzmann's factor.

Exact Solution of the EZM Model for Quasisigmoidal Denaturation Profiles
Albeit standard denaturation mechanisms usually follow a sigmoidal-like profile, the existence of an exact solution for the linear denaturation profile in EZM provides us with a powerful tool to proceed a step further towards a formulation of an analytically solvable SMCG model.While the sigmoidal denaturation profile is not analytically solvable within the EZM, it is possible to make a quite good approximation by means of a superposition of linear profiles.The idea is rather simple and it is illustrated in Figure 1.Sigmoidal denaturation profile is decomposed into the sum of three linear denaturation profiles with different slopes in such a way that sigmoidal denaturation is described as the superposition of three linear processes I, II, and III through which protein unfolds.Each process is characterized by a constant rate at which native contacts break.Initial and final processes (I and III, resp.)approximate quite well the sigmoidal profile since the derivative of hyperbolic tangent function is almost constant.However, intermediate process (II) has also a constant rate of native contact's breaking while hyperbolic tangent substantially differs from this.Nevertheless, linear and sigmoidal processes have a negative slope indicating a common decreasing mechanism.Except at the points   at which the derivative does not exist, the superposition of the processes I, II, and III approximates the rate at which native contacts are broken by denaturant.Let us consider three linear denaturation processes labeled I, II, and III according to the following equations: with { 1 ,  2 ,  3 } being their corresponding slopes and { 1 ,  2 ,  3 } being constants that ensure mathematical continuity to the superposition.At zero denaturant concentration it is expected that  = , whereas at critical concentration  ⋆ it is assumed that  = 0.At intermediate concentrations  Average rate at which native contacts break has a good correspondence between initial and final stages (processes I and III, resp.).The central discrepancy arises in the abrupt change artificially induced by the linear nature of process II.However, as we shall see, this is not a fundamental issue for the results obtained.Values used for quasisigmoidal pattern are given by  1 ( 1 = 3.5,   1 = 39/40) and  2 ( 2 = 5.2,   2 = 1/40), where process I intersects II and process II intersects process III, respectively.A chain of  = 40 residues is considered which is completely denatured at 8 Molar.
1 and  2 ,  1 and  2 residues are expected.This sets up the following equations: This defines a 6 × 6 system which has a unique solution.From this, we obtain The energy spectrum  can be written consequently as Hence, partition function Z  for the quasisigmoidal denaturation profile is defined as Each of the integrals involved in (24) has the same functional aspect as the one defined and solved previously in Section 3. Hence, an analytical full solution of each of the integrals in ( 24) can be accomplished to give an analytical expression for Z  .That is, where the quantities   are defined as follows: and quantities   are given in (22).Configurational partition function Z  defined in (25) provides us with the main frame to analyze equilibrium thermodynamic properties of polypeptide chains under the action of quasisigmoidal denaturation profiles.In order to investigate thermodynamic equilibrium properties we shall use Z  to evaluate configurational entropy S, average energy E, and thermodynamic susceptibilities like specific heat C V which contain substantial information on the thermodynamic response of the system.

Thermodynamic Equilibrium Response to Quasisigmoidal Denaturation Profiles
Given a particular quasisigmoidal denaturation profile two possible mechanisms can be identified.Figure 2 shows two different possibilities for chemical denaturation between concentrations  = 0 and  =  ⋆ = 6 M. Central difference is associated with the rate at which native contacts are broken by denaturant.As discussed previously, for a strict sigmoidal process the rate [44] at which native contacts are broken is essentially constant in the initial and final stages of denaturation, whereas the middle process has a highly nonconstant rate.A quasisigmoidal denaturation profile can approximate this property and in addition includes two variants.Figure 2(a) shows a denaturation mechanism composed by a slow and fast denaturation modes, whereas Figure 2(b) depicts a mechanism with fast and slow denaturation modes.
We shall call these two mechanisms type A and type B denaturation profiles, respectively.For type A, native contacts are broken in a slow mode initially and afterwards a fast denaturation mode is induced.In a type B process, native contacts are broken swiftly in initial stages and then a slow mode is induced.In both cases, initial and final concentrations are the same.A slow denaturation mode means that a larger amount of denaturant is needed to break down a certain fraction of native contacts in comparison to a fast denaturation mode.
If we restrict the model to helix-coil transitions, it is possible to establish a connection between the number ] of potential nonnative conformations which is a phenomenological parameter in the original ZW model and the helical propensity  of a residue to participate in a helix-coil transition [37].This relation proposes that in absence of an enthalpic energy barrier between helical and coil states, ] =  −1 .This relation is a reasonable first-order approximation since as a particular residue is more prone to participate in a helical state (native conformation), the probabilities that it participates in a coil state (nonnative conformation) clearly diminish.Such quantitative relationship enables us to particularize the solution obtained for Z  for specific homopolymer chains with the replacement ] =  −1 .Tables with measurements of   for each of the twenty or so aminoacid residues are widely known.Here, we shall use data obtained from Chakrabartty et al. [45].Figure 3 shows specific heat  V as function of temperature  for poly-alanine (poly-A) and poly-valine (poly-V) chains composed of  = 40 residues in a slow-fast denaturation (type A) and fast-slow denaturation (type B) mechanisms.According to [45] for alanine,  = 1.61 whereas valine has  = 0.23.This scale suggests that alanine is one of the strongest helix formers whereas valine is basically a helix breaker.Since ] ≃  −1 [37], poly-A chains would have a smaller number of average nonnative conformations in comparison to poly-V chains.Since helical propensity for valine,   , is smaller than helical propensity for alanine   , then ]  > ]  holds and larger peak in poly-V chain is expected.
The specific role of denaturation mode is showcased in Figure 3 as well.Figure 3(a) considers a poly-A chain subject both to a slow-fast denaturation mechanism (thick line) and a fast-slow denaturation mechanism (dotted line).It is observed that this factor plays a significative role in the thermodynamic stability of the chain.The transition temperature   varies appreciably for each case.The phase transition signaled by the peak in the specific heat  V occurs earlier for the slow-fast denaturation mode (type A) than for the fast-slow denaturation mode (type B).This phenomenon is, however, inverted in the case of the poly-V chain.In this case, the transition temperature   associated with the type B denaturation is smaller than its counterpart in the type A mechanism.The intensity in the  V maximum is also inverted.Type A mode is associated with a larger  V peak for poly-A chain whereas for poly-V chain type A denaturation mode is associated with the smaller  V peak.Within our model the difference between poly-A and poly-V chains lies exclusively in the helical propensity  and its related variable ].Hence, the number of potential nonnative conformations ] and the type of denaturation mechanism both play a fundamental role in thermodynamic stability of configurations, in particular in the transition temperature   .It is worth observing that ]  < 1 and ]  > 1; thus, ] = 1 signals a threshold on thermodynamic stability.Homopolymer chains poly-X with a larger amount of conformational degrees of freedom are expected to have higher denaturation temperatures since more thermal energy is needed to overcome nonnative conformations.Larger values of ] tend to shift   to higher values.Such observation is consistent with results shown in Figure (3).The value ] ̸ = 0 is not allowed since in that case each residue has only a single state and conformational transitions would be blocked.
For a further characterization of the equilibrium properties of the chains, we have calculated the configurational equilibrium entropy S as function of temperature for previous poly-A and poly-V chains.Figure 4 shows S() for poly-A (a) and poly-V (b) in type A and B denaturation mechanisms.Given the fact that S is a first-order derivative of the free energy F with respect to temperature, that is, a discontinuous change of  at at critical temperature   can be directed being associated with a first-order phase transition [30,46,47].The formal definition of a first-order phase transition requires that entropy undergoes a discontinuous change at the transition temperature   .However, for a finite poly-X chain, it is expected that S changes continuously.As the number of residues  increases, a true first-order phase transition is expected since the peak observed in  V becomes Figure 2: Two possible denaturation mechanisms depending on the rate at which native contacts break can be devised for a given quasisigmoidal denaturation profile.In a slow-fast mode denaturation (a), the rate at which native contacts are broken is essentially constant at the beginning and then a higher increase of the rate is observed at the end.For the fast-slow mode denaturation (b), a small amount of denaturant induces a higher rate of denaturation in the beginning.Then, the process is essentially constant up to the critical concentration.narrower and higher.Recently it has been shown that a firstorder phase transition in the limit of large  is embedded in this type of coarse-grained models [37].Figure 4 shows that although both poly-A and poly-V chains exhibit a first-order phase transition, such thermodynamic transition becomes stepper in the case of poly-V chains for the same values of parameters used.Smaller helical propensities imply then that thermodynamic transitions occur in a smaller window of temperature.For poly-V chain,  = 0.23 and the temperature range in which the transition occurs is Δ = 40 K approximately.For poly-A chains,  = 1.61 and Δ = 100 K, approximately.Entropy differences Δ between unfolded and folded structures ΔS ≃ 100 J K mol −1 for poly-A and ΔS ≃ 500 J K mol −1 suggest as well that poly-V denaturation is much more favorable in entropic terms than poly-A.This is consistent with the fact that within this Figure 4: Configurational entropy S for poly-alanine (poly-A) and poly-valine (poly-V) chains under a slow-fast denaturation mode (type A) and a fast-slow denaturation mode (type B).In both cases, the number of residues was set to  = 40,  = 0.75 × 10 3 J mol −1 , and  = 35×10 3 J mol −1 .Profile A is characterized by  1 = 0.32,  2 = 16.19, 3 = 6.25,  1 = 39, and  2 = 5 whereas profile B is defined through parameters  1 = 64,  2 = 0.44,  3 = 6,  1 = 8, and  2 = 6.Helical propensity   = 1.61 was considered for alanine while   = 0.23 was considered for valine.model poly-V residues have more nonnative conformations; thus more microstates are available to a poly-V chain in comparison with poly-A.

Conclusions and Final Remarks
A full analytical solution for the EZM model in the quasisigmoidal approximation model is presented in order to calculate equilibrium thermodynamic properties of proteinogenic homopolymer poly-X chains under the action of a chemically induced denaturation mechanism.As discussed previously, the original EZM model considers the exact sigmoidal denaturation curve common to most chemical denaturation processes.The EZM partition function is then expressed in terms of a highly nonlinear integrand on the number of microstates as function of .Given this fact, a full solution is not achievable.In this work, a quasisigmoidal representation for denaturation curve is used to obtain a full analytical solution to the EZM model.Results obtained here show that Z  can be obtained in terms of error functions Erf() and the slopes   of the individual linear denaturation mechanisms that compose the quasisigmoidal approximation.With this approximation, the EZM model can be included in the class of analytically solvable SMCG models.
For a given quasisigmoidal denaturation mechanism two submechanisms can be identified depending on the rate at which native contacts are broken by denaturant.The first process which we called type A denaturation mode is defined by a slow rate of native contacts breaking followed by a fast denaturation rate.The second process which we called type B mode is characterized by a fast followed by a slow rate of native contact breaking.Both type A and type B denaturation modes have the same initial and final stages.At  = 0,   = 1 and at  = 8 M,   = 0 in both cases.A critical value  * = 8 Molar to achieve complete denaturation was set according to reports on the literature for urea denaturation experiments, although a different value can be used.Results obtained here do not depend qualitatively on the exact value of  * .Using the fact that the number of potential nonnative conformations ] can be related thermodynamically to the helical propensity , thermodynamic properties like the specific heat C V and configurational entropy S as function of temperature were calculated for poly-alanine (poly-A) and poly-valine (poly-V) homopolymer chains.For these homopolymer chains, the influence of the type of denaturation mode was also discussed.
Specific heat  V calculations as function of temperature  and helical propensity  show that the mode of denaturation (either type A or B) does play a significant role in thermodynamic stability, specifically in the value of the transition temperature   .For poly-A chains, type B denaturation mode shows a higher degree of thermodynamic stability, that is, a larger value of the transition temperature   .However, for poly-V chains, the opposite situation is observed.This implies that although both denaturation modes have the same initial and final stages, the rate at which native contacts are broken does play a role in the location of a thermodynamic transition.Such is consistent with the fact that Z  indeed depends on the slopes  1 ,  2 , and  3 of the linear processes I, II, and III that form the quasisigmoidal denaturation profile. V curves also show that the intrinsic nature of the X residue in the poly-X chain has a role in thermodynamic properties.Since   >   and  ≃ ] −1 , poly-A chains exhibit a small International Journal of Statistical Mechanics number of nonnative microstates whereas for poly-V chains this number is larger.This is in fact what is observed in Figure 3.The intensity of the  V peak for poly-V is larger than its counterpart in poly-A.Main result obtained here is that the location of transition temperature   is both a function of the helical propensity  and denaturation mode.For instance, poly-X chains with higher propensities to undergo a random coil helix transition shall denature at lower temperatures on type A denaturation mode.
Results obtained for configurational entropy S() for both poly-A and poly-V chains confirm the fact that poly-A tends to have a smaller number of microstates than poly-V since entropy differences ΔS in poly-A and poly-V are different.For poly-A, ΔS ≃ 100 J K mol −1 whereas for poly-V, ΔS ≃ 500 J K mol −1 .Since S =  log Ω with Ω the number of microstates available, Figure 4 confirms that Ω  < Ω  as suggested by the specific heat plots.S() also contains information on the nature of the thermodynamic transition involved.Figure 4 suggests that  is the main driving factor for a steeper transition (left panel).Poly-X chains with values  < 1 display a steeper thermodynamic transition.For such chains, transition temperatures   are slightly affected by the type of transition (A or B) in contraposition to poly-X chains with  > 1 like poly-A where sensitive dependence on denaturation mode is observed.
Comparison of the results obtained in this work with MD simulations of poly-alanine chains in vacuo shows that conformation transitions do occur even in absence of solvent interactions [48].Conformational transitions predicted in this work signaled by a peak in the specific heat V (Figure 3) and configurational entropy (Figure 4) are in qualitative agreement with MD simulations regarding the fact that the steepness of the transition depends on the size of the chain [48].Smaller chains have "softer" transitions in comparison with larger ones.This fact is indeed what is expected from the general theory of equilibrium statistical mechanics [49,50] where true thermodynamic first-order phase transitions appear strictly as  → ∞ and the volume of the system remains constant.Recently sodiated poly-alanine chains Ala  Na + with  = 8, 9, 10, 11, which are poly-A chains with a Na + cation attached to the C terminus have been shown to exhibit helical stable conformations in vacuo as well suggesting that conformational transitions are robust phenomena in polypeptides with larger helical propensities even in the absence of solvent interactions [51].Free energies reported for  helical conformations of these peptides range from −30 kJ mol −1 for Ala 8 Na + up to −44.2 kJ mol −1 for Ala 10 Na + .This value can be compared with our results if we consider that the free energy change Δ can be estimated as where the internal energy change Δ = −( + ) ≃ −65 kJ mol −1 for the values considered in Figure 3 and Δ can be estimated from the results obtained in Figure 4. Since Δ ≃ 65 kJ mol −1 and Δ = 20 kJ mol −1 we can see that our results predict a value Δ ≃ −45 kJ mol −1 which is in quite good agreement with these MD results for short Alanine peptides.

Figure 1 :
Figure 1: Sigmoidal and quasisigmoidal denaturation profiles for a chemically induced protein denaturation.Quasisigmoidal denaturation profile is composed by the superposition of three linear processes called I, II, and III with different slopes (gray solid line).Average rate at which native contacts break has a good correspondence between initial and final stages (processes I and III, resp.).The central discrepancy arises in the abrupt change artificially induced by the linear nature of process II.However, as we shall see, this is not a fundamental issue for the results obtained.Values used for quasisigmoidal pattern are given by  1 ( 1 = 3.5,   1 = 39/40) and  2 ( 2 = 5.2,   2 = 1/40), where process I intersects II and process II intersects process III, respectively.A chain of  = 40 residues is considered which is completely denatured at 8 Molar.
2) written in terms of the continuous variable  instead of the discrete one .It describes a statistical average over energy conformations defined by variable  when energy spectrum  is a function also of .Standard denaturation experiments show that  usually ranges from zero to a maximum concentration  * = 8 Molar in a  →  denaturation reaction.