Binomlabs

Analysis of neutral and oncological mutations by calculating thermodynamic characteristics.

The correlations were performed on real solutions in a biochemical laboratory.
Study of P53 protein.
Work is underway, the page is being filled as it progresses.
GO TO Main Page
CONTACTS
For the negative set, protein variants that were never found (693 variants) or found only once( 323 variants) in human cancer have been selected (no_cancer p53 variants)
[TP53_PROF:a machine learning model to predict impact of missense mutations inTP53]
A diagram of the calculated values for each set of p53 protein mutations, which we will further analyze.

Cancer-related Mutations
TP53 mutation data was aggregated from multiple studies

Effects of Common Cancer Mutations on Stability

Silent and oncological mutations in p53

Heat map of reliable electrostatic interaction energy between pairwise amino acid residues of the p53 monomer
The difference in the calculations obtained for the Positive and Negative set of p53 mutations

Interplay between physical quantities
on the affinity of the complex

The best way to improve your professional skills and increase your value
  • 45%
    Free energy
    1 year
  • 67%
    Entropy change
    3 years
  • 15%
    Entalpy change
    2 years
  • 20%
    Hydrophobicity
  • 30%
    Solvatation
  • 55%
    Stability
Characteristics of physical quantities
Relationships and influences of physical quantities among themselves and correlation with affinity.
  • Entropy change
    The entropy change accompanying a chemical reaction is defined as the difference between the sum of the entropies of all the products and the sum of the entropies of all reactants the entropy change is given by:
    dS = (cSc + dSD + • • •) - (aSA + bSB+ • • ')
  • Gibbs energy
    the thermodynamic potential that is minimized when a system reaches chemical equilibrium at constant pressure and temperature when not driven by an applied electrolytic voltage
  • Entalpy change
    The lattice energy of an ionic crystal is the enthalpy change,
    required to decompose 1 mole of the crystal into its constituent gaseous ions at any temperature T.
  • Hydrophobicity
    It is commonly understood to be the tendency of non-polar molecules to form aggregates in order to reduce their surface of contact with polar molecules such as water

Lock-and-Key:
An Entropy-Dominated Binding Process

Therefore, for the lock-and-key binding to proceed, the solvent entropy gain should be large enough to overcompensate for not only the positive enthalpy change arising from the desolvation process, but also the negative entropy change caused by the loss of rotational and translational motions of the ligand.
Indeed, the negative enthalpy change arising
from the favorable interactions (such as van der Waals forces, hydrogen bonding, electrostatic, and dipole–dipole interactions) can also contribute to the lowering of the system’s free energy, but the solvent entropy gain arising from the displacement of the water molecules plays a dominant role in lowering the free energy. Therefore, it is reasonable to conclude that the lock-and-key binding is a entropy-dominated process.
For the conformational selection binding scenario,
it is difficult to distinguish which factor (the entropy or the enthalpy) contributes more to the lowering of the system’s free energy because the large solvent entropy gain in the first step could be offset by the loss of the rotational and translational entropy and the decrease of the conformational entropy in the subsequent step, and the negative enthalpy change in the second step could be offset by the positive enthalpy changes due to the desolvation energy penalty and the disruption of the original noncovalent interactions surrounding the binding sites.

Nevertheless, the selective binding and the following conformational adjustments are dominated by the solvent entropy gain and the system enthalpy decrease, respectively, suggesting that they play a role, in a sequential manner, in lowering the system’s free energy.
In the last decades, the quest thereof has been dominated by structural and thermodynamical considerations, where candidate ligands are ranked by their equilibrium constants Kd. This constant, also known as the dissociation constant, describes the fraction of target-proteins that would be bound to the ligand at equilibrium conditions, and in particular, at constant ligand concentration. A schematic picture of protein-ligand binding is given in Figure 1, where it is seen that the equilibrium constant is given by the ratio of the kinetic rate constants koff and kon. From physical chemistry, it is known that Kd is proportional to the Boltzmann constant associated with the free energy difference of the bound and unbound states (Figure 1, left).
Figure 1: Left: schematic picture of the Protein-Ligand (P-L) binding free energy profile. The equilibrium constant Kd is fully described by the free energy difference deltaG of the unbound and bound states, describing the strength of the P-L bound complex. For this reason, Kd has served as the dominant parameter in targeted drug design. However, the binding reaction is almost always accompanied by an energy barrier which is often not negligible, and can highly impact the effectiveness of the candidate ligand in vivo (see text). Right: differential equation accompanying the bimolecular binding reaction, where [P] and [L] denote protein and ligand concentrations, respectively.
While it is true that Kd is a good indicator of binding strength, it has major shortcomings that prove to be detrimental: more than 90% of promising drug molecules fail human trials by off-target toxicity and/or bad overall target selectivity. Indeed, thermodynamic selection based on Kd neglects the existence of a free energy barrier that needs to be overcome during the binding reaction (Figure 1, left). While the final bound state at the target protein may be energetically most favourable, the energy barrier may be sufficiently high that most drug molecules bind to a – potentially toxic – off-target protein with less favourable bound state, but lower energy barrier. To this end, the notion of kinetic selectivity has gained a lot of attention in the latest years, where one tries to obtain the rate constants kon and koff. These parameters hold information about kinetic properties of the protein-ligand interaction, such as the drug residence-time, allowing better profiling of candidate drug molecules and saving costly and time-consuming in vivo experiments.
https://plato.ea.ugent.be/mp/export1.jsp?nr=25223&view=ntpbrPD
Models of drug–target binding. a Schematic diagram of a one-barrier drug–target binding free energy profile. A one-step model with one free energy barrier is used to derive the experimental rate constants. The figure and equations show how the steady-state rate constants relate to the free energy differences shown. The residence time of a drug bound to its target, τ (which is the reciprocal of the rate constant for dissociation of the drug–target complex, k off), results from the “difference” in free energy between the transition state (TS) and the bound ground state (GS), ΔG off.
The red arrows indicate that prolongation of the τ can be achieved by stabilizing the GS (increasing the magnitude of ΔG D), destabilizing the TS (increasing ΔG on) or a combination of both (i.e.,
b Diagram schematically illustrating different mechanisms of drug binding involving protein conformational changes.
R and RC denote two different conformations of the protein, the latter requires conformational changes for ligand binding.

These may occur by conformational selection (blue path) or by induced fit upon formation of an encounter complex [RL]# (red path), or by a combination of the two mechanisms. Binding proceeds through an energetically unfavorable intermediate state (TS in panel A or a local minimum in a 2 (or more)-step binding free energy profile) that, in the conformational selection and induced fit mechanisms, corresponds, respectively, to the R+L or [RL]# state of the system); the final complex is denoted by [RL].
The gray path and third equation describe the pseudo-one-step binding process shown in (a) is used to derive the experimental rate constants
[Protein conformational flexibility modulates kinetics and thermodynamics of drug binding]

Thermodynamic profiles of N-HSP90 inhibitors measured by ITC. The enthalpic and entropic components of the binding free energy are shown in a and b for for WT N-HSP90, and the L107A mutant, respectively. The dashed diagonal line (ΔH=−TΔS) divides the plot into two main areas where the enthalpy (gray) or the entropy (red) dominate the binding free energy (ΔG).
Thermodynamic profiles of N-HSP90 inhibitors
The thermodynamic profiles of the binding of the 20 resorcinol ligands to N-HSP90 obtained by isothermal titration calorimetry (ITC) are depicted in Fig. 3 and quantify the energetic differences between two states in equilibrium (free state and bound state).
Free Energy Calculations

As a reward for the highly intensive computation, the results of free energy calculations ought to be reliable and almost quantitative. The main advantages over faster scoring functions are that the free energy calculations include both the energetic (i.e., potential energy and solvation energy) and entropic (i.e., dynamics/flexibility of both protein and ligand, and solvent effects) contributions, and require no case-by-case parameter fitting
Anyway, the accurate prediction of binding free energy using the calculation methods, despite being an ambitious goal, would revolutionize their applications in basic research and in drug design and discovery.
Free energy calculations rely on the fundamental relationship between Helmholtz free energy F and the partition function Z:
F = −kBTlnZ
where kB and T is Boltzmann’s constant and temperature, respectively. When the system is treated in terms of the classic approximation of statistical thermodynamics, the partition function can be expressed as configurational integral:
where h is Planck’s constant, N is the number of atoms or particles in system, and N! is only present for indistinguishable particles.
The latter describes the interactions between the various atoms in system (i.e., potential function). The difference in free energy between two states A and B can be expressed as a ratio of their partition functions:
If the conformational sampling is carried out under constant temperature and pressure conditions (isothermal-isobaric ensemble), the Gibbs free energy can be obtained.
[Insights into Protein–Ligand Interactions: Mechanisms, Models, and Methods]

Free energy calculations of the protein-ligand binding try to compute the binding free energies based on the principles of statistical thermodynamics.
[Insights into Protein–Ligand Interactions: Mechanisms, Models, and Methods]
[Some Binding-Related Drug Properties are Dependent on Thermodynamic Signature]
Thermodynamic profiles for three pairs of HIV-1 proteinase inhibitors that vary by only a single group: (a) KNI-10033-KNI-10075 pair within which an apolar group thioether on KNI-10033 is replaced by a polar group sulfonyl to form KNI-10075; (b) KNI-10052-KNI-10054 pair within which an apolar methyl group is replaced by a polar hydroxyl group; and (c) KNI-10046-KNI-10030 pair within which a hydrogen atom on the former is replaced by an apolar methyl group to form the latter. The binding free energy (∆G), enthalpy (∆H), and entropy (T∆S) are shown
[Freire, E. The binding thermodynamics of drug candidate. In Thermodynamics and Kinetics of Drug Binding; Keserü, G.M., Swinney, D.C., Eds.; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2015; pp. 1–13.]
[Lafont, V.; Armstrong, A.A.; Ohtaka, H.; Kiso, Y.; Mario Amzel, L.; Freire, E. Compensating enthalpic and entropic changes hinder binding affinity optimization. Chem. Biol. Drug Des. 2007, 69, 413–422.]
[Kawasaki, Y.; Freire, E. Finding a better path to drug selectivity. Drug Discov. Today 2011, 16, 985–990.]
NBD-556 is a competitive inhibitor of CD4 characterized by a binding affinity of 3.7 μM. Despite the small size, NBD-556 binds with a thermodynamic signature that resembles that of CD4 (ΔH = -24.5 kcal/mol, -TΔS = 17.1 kcal/mol at 25 °C
dH, kcal/mol
dH, kcal/mol
-TΔS, kcal/mol
-TΔS, kcal/mol
The thermodynamic signatures of sCD4 and NBD-556 at 25°C. The large conformational structuring of gp120 triggered by CD4 binding is reflected in a thermodynamic signature characterized by an unusually large favorable change in enthalpy and a very large unfavorable entropy change. Except for a lower affinity, the binding of NBD-556 to gp120 is also associated with enthalpy and entropy changes similar to those observed for CD4. ΔG is represented by blue bars, ΔH by green bars and -TΔS by red bars.
Free energy calculations of protein-ligand complexes
The entropy change accompanying a chemical reaction is defined as
the difference between the sum of the entropies of all the products and
the sum of the entropies of all reactants.
the entropy change is given by:
dS = (cSc + dSD + • • •) - (aSA + bSB+ • • )
where SA, SB, etc., are the entropies per mole of the various species.

When the entropies of the individual substances correspond to a state of
unit activity they are called standard entropies and they are designated by
the symbol S°.

Again, in a reaction where all the substances involved are
at unit activity dS is written dS0, and the latter is the standard entropy
change of the reaction.
One way of obtaining these is through the third law of thermodynamics.

ENTROPY CHANGE IN CHEMICAL REACTIONS
Li and Deber [1] used circular dichroism (CD) data to rank order helical propensity of proteins within membranes. Residues such as Ile, Val and Thr, which usually exist as β-sheets in an aqueous environment, prefer an α-helical conformation in lipid membranes. Thus, the helical propensity of amino acid residues correlates with the hydrophobic nature of the side chain
Felitsky et al. introduced the use of a new parameter called “average area buried upon folding” (AABUF) [2] that explains both local contacts and long-range interactions. AABUF was used to study folding of apomyoglobin and provided additional insight into hydrophobic collapse and early folding events.
[1] A measure of helical propensity for amino acids in membrane environments
[2] Modeling transient collapsed states of an unfolded protein to provide insights into early folding events
This principle, called the third law of thermodynamics, states that the entropy of all pure crystalline solids may be taken as zero at the absolute zero of temperature.

The importance of the third law lies in the fact that it permits the calculation
of absolute values of the entropy of pure substances from thermal
data alone.

THE THIRD LAW OF THERMODYNAMICS

Further, since the third law states that for
any pure crystalline substance S = 0 at T = 0
then for any such substance the equation may be integrated
between this lower limit and any temperature T to yield:
St, known as the absolute entropy of the solid, is always a positive quantity.
All that is necessary for integration of this equation is a knowledge of the heat capacities of the solid from T = 0 to any desired temperature T
The integral in Eq. above is evaluated graphically by plotting either Cp
vs. In T or Cp/T vs. T and determining the area under the curve between
T = 0 and any temperature T. This area is then the value of the integral
and hence of ST. In practice heat.capacities are usually measured from
approximately 20°K to temperature T, and extrapolation is resorted to
from ca. 20°K to T = 0.


Such extrapolations are generally made by use of the Debye third power law for the heat capacity of solids at low temperatures, namely,

EVALUATION OF ABSOLUTE ENTROPIES

This expression substituted into Eq. above yields
and hence on integration between T=0 and Ti, the lower limit of the heat capacity data, we get
Since Cp = a(T)^3, then the value of ST1 is equal to one-third of the value
of Cp at temperature T1
A plot of Cp vs. log10 T for anhydrous sodium sulfate is shown in Figure. The area under the curve between T = 14° and
T = 298.15°K is 15.488, while the area from T = 0 to T = 14° is 0.026. Hence the absolute entropy per mole of sodium sulfate at 298.15°K is
Absolute entropies of substances that are liquid or gaseous at room temperatures can also be obtained with the third law.
Heat capacity of Na2S04 at various temperatures.
Absolute entropies are usually calculated and tabulated for 25°C and
unit activity. Values of these thermodynamic quantities in the standard
state are available at present for most elements and for many compounds.
Some of these are listed in Table:


USE OF ABSOLUTE ENTROPIES
STANDARD ABSOLUTE ENTROPIES OF ELEMENTS
AND COMPOUNDS AT 25°C
(Entropy units per gram atom or mole)
Suppose the entropy change is required for the reaction:
C(s, graphite) + 2 H2(g) + l- 0,(g) = CH30H(1)
In a similar manner may be calculated AS0 values for other reactions at 25°C provided the necessary absolute entropies are known. To obtain the entropy changes at temperatures other than 25°C:
Although entropy can be employed to measure the tendency of systems
to undergo change, under conditions most frequently encountered it is
not as convenient a quantity to use as the so-called free energy functions
The Helmholtz free energy of any system, A, is defined as A = E - TS

where E and S are, respectively, the internal energy and entropy of a
system. Since E, T, and S depend only on the state of the system, A
must also be a state function. Consequently, when the system passes
from one state to another, the change in A must be given by
THE HELMHOLTZ FREE ENERGY
Equation (2) gives the most general definition of dA. Under isothermal
conditions, when T2 = T1 = T, this equation reduces to
dA = A 2 — A i= dE - (T2S2 - T1S1) (2)
dA = dE - TdS (3)
Equation (3) allows a physical interpretation of dA. Since under isothermal conditions TdS = qr ,we have
dA = dE - qr
Hence at constant temperature the maximum work done by a system is
accomplished at the expense of a decrease in the Helmholtz free energy
of the system. This is why A is sometimes called the maximum work
content of a system.
Следовательно, при постоянной температуре максимальная работа, совершаемая системой, достигается за счет уменьшения свободной энергии Гельмгольца.
Вот почему А иногда называют максимальной работой системы.
The maximum work a process may yield is not necessarily the amount
of energy available for doing useful work, even though the process is
conducted reversibly. Of the total amount of work available, a certain
amount has to be utilized for the performance of pressure-volume work
against the atmosphere due to contraction or expansion of the system
during the process.
THE GIBBS FREE ENERGY
For any pure substance A is most conveniently expressed in terms of T and V as the independent variables:
Equation shows the dependence of A for a pure substance dn both the temperature and volume, with the individual effects of these variables
Максимальная работа, которую может произвести процесс, не обязательно равна сумме энергии, доступной для совершения полезной работы, даже если процесс проводится обратимо. Из общего объема имеющихся работ определенное количество должно быть использовано для выполнения работы давление-объем против атмосферы из-за сжатия или расширения системы во время процесса.
For a process taking place reversibly at constant temperature and pressure, and involving a volume change from V1 to V2,
the work done against the atmosphere is P(V2 — Vi) = PdV.

Since this work is accomplished at the expense of the maximum work yielded by the process, the net amount of energy available for work other than pressurevolume against the confining atmosphere must be
Net available energy at T and P
To bring out more precisely the nature of the maximum net energy
available from a process, let us define another state function F, called the Gibbs free energy, by the relation F = H - TS
The change in F between two states of a system will be, therefore,
and when the temperature is constant
An alternate, but equivalent, way of defining F is through the relation
F = A + PV
where A is the Helmholtz free energy, and P and V are the pressure and volume of the system.
which at constant pressure becomes
dF=dA+PdV
The equivalence of the two definitions may be shown by inserting Eq's. We obtain thus
Hence the total maximum work wm:
where dqr is the infinitesimal quantity of heat absorbed in a process taking place under reversible conditions at temperature T. In the case of a finite reversible change at constant temperature, dS becomes AS, dqT becomes qr
Therefore, for any isothermal reversible process in which an amount of heat qr is absorbed at temperature T, the entropy change involved is simply the absorbed heat divided by the absolute temperature. When qr is positive, i.e., heat is absorbed, AS is also positive, indicating an increase in the entropy of the system. On the other hand, when heat is evolved qT is negative and so is dS, and the system experiences a decrease in entropy.

Следовательно, для любого изотермического обратимого процесса, в котором количество тепло qr поглощается при температуре T, связанное с этим изменение энтропии представляет собой просто поглощенное тепло, деленное на абсолютную температуру. Когда qr положителен, т. е. происходит поглощение тепла, AS также положителен, что указывает на увеличение энтропии системы. С другой стороны, когда выделяется тепло, qT отрицательно, как и AS, и в системе происходит уменьшение энтропии.
The physical significance of dF at constants T and P may be obtained as follows: At constant temperature TdS = qr. Again, when the pressure is also constant, dH = dE + PdV.
dF = dE + PdV - qr
= -(qr - dE-PdV)
But, by the first law qr — dE = wm, and, therefore,
Good design is obvious. Great design is transparent.
dF = -(wm -PdV)
dF represents the maximum net energy at constant T and P available for doing useful work, the net available energy under the specified conditions results from a decrease in the free energy content of the system on passing from the initial to the final state.

It is customary to refer to the Gibbs free energy simply as the free
energy. We shall follow this procedure, and append the designation Gibbs
or Helmholtz only when necessary for clarity.
Like H, F is most conveniently expressed in terms of T and P as
independent variables.
dF представляет собой максимальную чистую энергию при постоянных T и P, доступную для совершения полезной работы, чистая доступная энергия в указанных условиях является результатом уменьшения содержания свободной энергии в системе при переходе от начального к конечному состоянию.

Свободную энергию Гиббса принято называть просто свободной энергией. Мы будем следовать этой процедуре и добавим обозначение Гиббса или Гельмгольца только тогда, когда это необходимо для ясности, подобно H, F удобнее всего выражать через T и P как независимые переменные.
A minus sign denotes, therefore, that the reaction tends to proceed spontaneously. When the tendency is from right to left, however, net work equivalent to dF has to be absorbed in order for the reaction to proceed in the direction indicated, and dF is positive. A positive sign for dF signifies, therefore, that the reaction in the given direction is not spontaneous. Finally, when the system is in equilibrium, there is no tendency to proceed in either direction, no work can be done by the system, and hence dF = 0.
These three possible conditions for the free energy change of a process at constant temperature and pressure may be summarized as follows:
When a catalyst like platinized asbestos is introduced, however, the reaction
proceeds with explosive violence. Even with a catalyst the reaction
would have been impossible had not the potentiality to react been present.
Однако при введении такого катализатора, как платинированный асбест, реакция протекает сo взрывной силой. Даже при наличии катализатора реакция была бы невозможна, если бы не существовала возможность реакции.
The arrows indicate the directions the reaction tends to follow spontaneously for the given sign of the free energy change.
A negative free energy change for a process does not necessarily mean that the process will take place. It is merely an indication that the process can occur provided the conditions are right.
Platinized asbestos is a catalyst material obtained by permeating asbestos fibers with finely divided platinum, usually via the chemical method. It is a white fibrous solid containing around 5% platinum metal by weight.
Platinized asbestos used to be sold in the past by chemical suppliers, but nowadays it's getting harder to find since most asbestos-containing products have been phased out. It is very toxic and irritant and should be handled with care. Inhalation of fine particles will lead to severe lung problems, while skin or tissue contact will cause irritation or dermatitis.
Uses
  1. Platinised asbestos is used as a catalyst in the manufacture of H2SO4.
  2. Platinised asbestos helps in the formation of SO3, from S02 and O2.
It is the sign of the free energy change which determines whether the
potentiality to react exists, and it is the magnitude of the free energy
change which tells us how large that potentiality is.
The term represents the heat interchange between the system and its surroundings when the process is conducted isothermally end reversib
When dH is greater than dF, qr is positive and energy is absorbed as heat from the surroundings. On the other hand, when dF is greater than dH, qr is negative and heat is evolved to the surroundings. Finally, in the special case when
dH = dF, heat is neither absorbed nor evolved by the system, and hence ther is no change in entropy.
— TdS or
The difference between the work so obtained and the heat which would have been liberated had the reaction been carried out completely irreversibly, as in an open beaker,
is given by the Gibbs-Helmholtz equation, and is equal to either
amount of heat qr
dF equations are useful not only for calculating the free energy change at any temperature T, but also for evaluating dS and dH.
dS is readily obtained by differentiation with respect to temperature at constant pressure of the dF expression in accordance with Eq.
Cp is the heat capacity per mole
An approach to the apolar molecule-water interaction was adopted by
W. Kauzmann. Rather than focus on the formation of crystalline hydrates, he considered the thermo dynamics of transferring an apolar residue from an apolar solvent to water. For example, he asked what the free energy change would be in transferring methane from mole fraction x in benzene to mole fraction x in water. This kind of approach basically aims at determining whether the apolar groups would prefer to reside on the protein surface—and thus interact with water—or bring themselves together in the interior of the protein—and thus create an apolar environment for themselves.
Unitary and cratic contributions to free energy of transfer
unitary and cratic contributions to free energy
Lets concider a standard free energy change

dG°=-RTlnK
The dependence of dG° on the standard state is a purely statistical effect that arises from the expression for the partial molal entropy; consider the partial molal entropy Sa of A:
where Xa is the mole fraction of A. (We are considering dilute solutions, and therefore we are ignoring activity-coefficient effects.)
The term — RInXa is called the cratic contribution to the entropy; it is a purely statistical term arising from the mixing of A with solvent molecules.
This contribution is independent of the chemical nature of A.
The term S'a is called the unitary contribution; it reflects the characteristics of A itself and its interaction with solvent.
The standard entropy change for the reaction becomes:
Equation above is the usual expression for the partial molar entropy. It can easily be derived by starting with the total differential for the Gibbs free energy:
DERIVATION OF EXPRESSION FOR PARTIAL MOLAR ENTROPY
Introductory physical chemistry textbooks derive the following relationship:

is the standard chemical potential at Xi = 1.
Using this relationship, and differentiating the preceding expression, we obtain
where Xm is the mole fraction of a molecule when its concentration is
1 mole/ liter. Because the water concentration is 55.6 molar in dilute aqueous solution, the cratic contribution to the entropy change is about — Rln56= — 8 cal /(mol deg) at 300 K, this corresponds to a cratic contribution to the free energy of about — 2.4 kcal/ mole.
If a different standard state and different concentration units were chosen, a different contribution would result. The unitary change in entropy, dSu, thus is obtained from:
As introduced above, two thermodynamic quantities, the enthalpy change and entropy change, determine the sign and magnitude of the binding free energy.

We therefore consider ΔH and ΔS as the driving factors for protein–ligand binding. The contributions of ΔH and ΔS to ΔG are closely related.

For instance, the tight binding resulting from multiple favorable noncovalent interactions between association partners will lead to a large negative enthalpy change, but this is usually accompanied by a negative entropy change due to the restriction of the mobility of the interacting partners, ultimately resulting in a medium-magnitude change in binding free energy.

Similarly, a large entropy gain is usually accompanied by an enthalpic penalty (positive enthalpy change) due to the energy required for disrupting noncovalent interactions. This phenomenon—the medium-magnitude free energy change caused by the complementary changes between enthalpy and entropy—is called the enthalpy–entropy compensation.
The main criticisms are that the compensation could be
  • a misleading interpretation of the data obtained from a relatively narrow temperature range or from a limited range for the free energies
  • the result of random experimental and systematic errors
  • and (iii) the result of data selection bias

Например, прочное связывание, возникающее в результате множественных благоприятных нековалентных взаимодействий между партнерами ассоциации, приведет к значительному отрицательному изменению энтальпии, но это обычно сопровождается отрицательным изменением энтропии из-за ограничения подвижности взаимодействующих партнеров, что в конечном итоге приводит к изменение средней величины свободной энергии связи.

Точно так же большой прирост энтропии обычно сопровождается энтальпийным штрафом (положительным изменением энтальпии) из-за энергии, необходимой для разрушения нековалентных взаимодействий. Это явление — изменение свободной энергии средней величины, вызванное дополнительными изменениями между энтальпией и энтропией, — называется энтальпийно-энтропийной компенсацией.

Lock-and-Key:
An Entropy-Dominated Binding Process

Therefore, for the lock-and-key binding to proceed, the solvent entropy gain should be large enough to overcompensate for not only the positive enthalpy change arising from the desolvation process, but also the negative entropy change caused by the loss of rotational and translational motions of the ligand.
Indeed, the negative enthalpy change arising
from the favorable interactions (such as van der Waals forces, hydrogen bonding, electrostatic, and dipole–dipole interactions) can also contribute to the lowering of the system’s free energy, but the solvent entropy gain arising from the displacement of the water molecules plays a dominant role in lowering the free energy. Therefore, it is reasonable to conclude that the lock-and-key binding is a entropy-dominated process.

Conformational Selection:
A Process in Which Entropy and Enthalpy Play Roles in a Sequential Manner

the selective binding may be dominated by the solvent entropy gain.
the presence of the conformational flexibility in the protein allows for the conformational adjustments of the residue side chaine
For the conformational selection binding scenario,
it is difficult to distinguish which factor (the entropy or the enthalpy) contributes more to the lowering of the system’s free energy because the large solvent entropy gain in the first step could be offset by the loss of the rotational and translational entropy and the decrease of the conformational entropy in the subsequent step, and the negative enthalpy change in the second step could be offset by the positive enthalpy changes due to the desolvation energy penalty and the disruption of the original noncovalent interactions surrounding the binding sites.

Nevertheless, the selective binding and the following conformational adjustments are dominated by the solvent entropy gain and the system enthalpy decrease, respectively, suggesting that they play a role, in a sequential manner, in lowering the system’s free energy.
[Insights into Protein–Ligand Interactions: Mechanisms, Models, and Methods]
[Some Binding-Related Drug Properties are Dependent on Thermodynamic Signature]
Thermodynamic profiles for three pairs of HIV-1 proteinase inhibitors that vary by only a single group: (a) KNI-10033-KNI-10075 pair within which an apolar group thioether on KNI-10033 is replaced by a polar group sulfonyl to form KNI-10075; (b) KNI-10052-KNI-10054 pair within which an apolar methyl group is replaced by a polar hydroxyl group; and (c) KNI-10046-KNI-10030 pair within which a hydrogen atom on the former is replaced by an apolar methyl group to form the latter. The binding free energy (∆G), enthalpy (∆H), and entropy (T∆S) are shown
[Freire, E. The binding thermodynamics of drug candidate. In Thermodynamics and Kinetics of Drug Binding; Keserü, G.M., Swinney, D.C., Eds.; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2015; pp. 1–13.]
[Lafont, V.; Armstrong, A.A.; Ohtaka, H.; Kiso, Y.; Mario Amzel, L.; Freire, E. Compensating enthalpic and entropic changes hinder binding affinity optimization. Chem. Biol. Drug Des. 2007, 69, 413–422.]
[Kawasaki, Y.; Freire, E. Finding a better path to drug selectivity. Drug Discov. Today 2011, 16, 985–990.]
NBD-556 is a competitive inhibitor of CD4 characterized by a binding affinity of 3.7 μM. Despite the small size, NBD-556 binds with a thermodynamic signature that resembles that of CD4 (ΔH = -24.5 kcal/mol, -TΔS = 17.1 kcal/mol at 25 °C
dH, kcal/mol
dH, kcal/mol
-TΔS, kcal/mol
-TΔS, kcal/mol
The thermodynamic signatures of sCD4 and NBD-556 at 25°C. The large conformational structuring of gp120 triggered by CD4 binding is reflected in a thermodynamic signature characterized by an unusually large favorable change in enthalpy and a very large unfavorable entropy change. Except for a lower affinity, the binding of NBD-556 to gp120 is also associated with enthalpy and entropy changes similar to those observed for CD4. ΔG is represented by blue bars, ΔH by green bars and -TΔS by red bars.
Free energy calculations of protein-ligand complexes
When the dissociation is into normal atoms, the spectroscopic dissociation energy Ds can be calculated directly from the observed frequency. However, when tfie dissociation results in one or more excited atoms, a correction has to be applied for the energy of excitation.
Typical dissociation constants for simple two-subunit systems range from:
These correspond to free energies of dimerization of —11 to —22
kcal /mole at 25°C.

When two protein subunits are brought together to form a complex, there is a loss of three translational and three rotational degrees of freedom, because the two subunits no longer can move independently. Many side chains free to move on the surfaces of the separate proteins will become immobilized at the subunit interface. These effects will lead to an intrinsic entropy loss that must be overcome by interactions between the two subunits if the complex is to be stable.
to
Much of this energy appears to come from the hydrophobic interactions
Correctly identifying the true driver mutations in a patient’s tumor is a major challenge in precision oncology. Most efforts address frequent mutations, leaving medium-and low-frequency variants mostly unaddressed. For TP53, this identification is crucial for both somatic and germline mutations, a multi organ cancer predisposition.

Here we combine calculated physical quantities and the application of machine learning methods to their processing to separate oncogenic mutations from neutral mutations of the P53 protein.

Most of p53 mutations in cancers are missense mutations, which produce the full-length mutant p53 (mutp53) protein
Fechner Correlation

The FECHNER CORRELATION command calculates the Fechner signs correlation coefficient between all the pairs of variables. Fechner correlation coefficient is used to check relationship for small samples.

How To
Run: STATISTICS->NONPARAMETRIC STATISTICS-> FECHNER CORRELATION...
Select the variables you want to correlate.
 Pairwise deletion is default for missing values removal (use the MISSING VALUES option in the PREFERENCES window to force casewise deletion).
Results:
Matrix with Fechner correlation coefficients between each pair of variables is calculated.
Fechner correlation coefficient is defined by
Negative Set of p53 mutations
p.A119P,p.A119S, p.A129G, p.A129P, p.A129S, p.A138G, p.A189S, p.C124F, p.C124W, p.C182F, p.C229G, p.C229W, p.D148A, p.D148V, p.D184A, p.D184E, p.D184V, p.D186A, p.D186E, p.D186Y p.D207A, p.D207V p.E171A, p.E171V, p.E180A, p.E180V, p.E198A, p.E198V, p.E224Q, p.E287A, ,p.F109Y, p.F113Y, p.F134Y, p.F212C, p.G108A, p.G108C, p.G112A, p.G112C, p.G112R, p.G112V, p.G117A, p.G117V, p.G117W, p.G154A, p.G154R, p.G187A, p.G226C, p.G262A, p.G262C p.G262R, p.G279A, p.H115D, p.H115L, p.H115N, p.H115P, p.H115Q, p.H115R p.H168Q, p.H178L, p.H214N, p.H233N, p.I162L, p.I195L, p.I195V, p.I255L, p.K101I, p.K101N, p.K101Q, p.K101T, p.K120T, p.L111V, p.L114F, p.L114M, p.L114S, p.L114V, p.L114W, p.L137R, p.L188Q, p.L188R, p.L201M, p.L201W, p.L206F, p.L206M, p.L206V, p.L206W, p.L252R, p.L257M, p.L264P, p.L264Q p.L264V, p.L265V, p.M169L, p.S106T, p.S116A, p.S116T, p.S116Y, p.S121A, p.S121C, p.S121T, p.S121Y, p.S127A, p.S149A, p.S149C, p.S183A, p.S183T, p.S185C, p.S185T, p.S227A, p.S227Y, p.S261I, p.T102N, p.T102P, p.T118K, p.T118P, p.T118R, p.T118S, p.T123N p.T123P, p.T123S, p.T125S, p.T140N, p.T150S, p.T170K, p.T256R, p.T284R, p.V122E, p.V122G, p.V122M, p.V147L, p.V203G, p.V225L, p.Y103C, p.Y103D, p.Y103F, p.Y103H, p.Y103N, p.Y103S, p.Y107F, p.Y107N, p.Y107S
Positive Set of p53 mutations
p.A138V, p.A159P, p.A159V, p.A161S, p.A161T, p.A276D, p.A276G, p.A276P, p.C124G, p.C135F, p.C135R, pC135S, p.C135W, p.C135Y, p.C141G, p.C141R, p.C141W, p.C141Y, p.C176F, p.C176G, p.C176R, p.C176W, p.C176Y, p.C238F, p.C238R, p.C238Y, p.C242F, p.C242G, p.C242S, p.C242Y, p.C275F, p.C275G, p.C275R, p.C275W, p.C275Y, p.C277F, p.C277Y, p.D259V, p.D259Y, p.D281E, p.D281H, p.D281N, p.D281V, p.D281Y, p.E180K, p.E224D, p.E258A, p.E258K, p.E258Q, p.E271K, p.E271V, p.E285K,p.E285V, p.E286A, p.E286G, p.E286K, p.E286Q, p.E286V, p.E287D, p.F109C, p.F109V, p.F113C, p.F113V, p.F134C, p.F134L, p.F134V, p.F270C, p.F270L, p.F270S, p.G105C, p.G105D, p.G105V, p.G154V, p.G199V, p.G244C, p.G244D, p.G244S, p.G244V, p.G245C, p.G245D, p.G245R, p.G245S, p.G245V, p.G262V, p.G266E, .G266R, p.G266V, p.G279E
The fact is that the yi values for Negative Set practically coincide with the obtained average value; as a result of the difference between the average and yi values, we get a very small value tending to zero.
Neutral Mutations
Oncological mutations
WOW! affect
Difference in calculated characteristics for Positive and Negative sets of P53 protein mutations.
To determine the correlation coefficient, we use Boolean variables and the Fechner coefficient.
rb(integr) is a Fechner coefficient
enthalpy change
N1-N60
N1-N60
Negative set
Negative set
Positive set
Positive set
Positive set
Negative set
N53-N140
N2-N89
p.A119P
p.A119S
p.A129G
p.A129P
p.A129S
p.A138G
p.A189S
p.C124F
p.C124W
p.C182F
p.C229G
p.C229W
p.D148A
p.D148V
p.D184A
p.D184E
p.D184V
p.D186A
p.D186E
p.D186Y
p.D207A
p.D207V
p.E171A
p.E171V
p.E180A
p.E180V
p.E198A
p.E198V
p.E224Q
p.E287A
p.F109Y
p.F113Y
p.F134Y
p.F212C
p.G108A
p.G108C
p.G112A
p.G112C
p.G112R
p.G112V
p.G117A
p.G117V
p.G117W
p.G154A
p.G154R
p.G187A
p.G226C
p.G262A
p.G262C
p.G262R
p.G279A
p.H115D
p.H115L
p.H115N
p.H115P
p.H115Q
p.H115R
p.H168Q
p.H178L
p.H214N
0.73993787
0.73993397
0.7399646
0.73995984
0.73995537
0.73985441
0.7399354
0.7375981
0.73758668
0.73957824
0.73991268
0.73990698
0.73998011
0.73997562
0.74010023
0.74000325
0.73998995
0.74002283
0.74000103
0.73997485
0.74000402
0.73997575
0.74002357
0.73997244
0.74015927
0.73999621
0.74005368
0.73997331
0.74003208
0.73998698
0.73999584
0.73997822
0.7399458
0.74006543
0.73999537
0.74001586
0.74005163
0.74015792
0.74014309
0.73999474
0.74011601
0.74000749
0.74000731
0.74000035
0.74007487
0.74001805
0.74003
0.74
0.74002
0.74008
0.74007
0.74
0.73998
0.73999
0.73998
0.74007
0.74013
0.74006
0.74008
0.73999
p.A138V
p.A159P
p.A159V
p.A161S
p.A161T
p.A276D
p.A276G
p.A276P
p.C124G
p.C135F
p.C135R
p.C135S
p.C135W
p.C135Y
p.C141G
p.C141R
p.C141W
p.C141Y
p.C176F
p.C176G
p.C176R
p.C176W
p.C176Y
p.C238F
p.C238R
p.C238Y
p.C242F
p.C242G
p.C242S
p.C242Y
p.C275F
p.C275G
p.C275R
p.C275W
p.C275Y
p.C277F
p.C277Y
p.D259V
p.D259Y
p.D281E
p.D281H
p.D281N
p.D281V
p.D281Y
p.E180K
p.E224D
p.E258A
p.E258K
p.E258Q
p.E271K
p.E271V
p.E285K
p.E285V
p.E286A
p.E286G
p.E286K
p.E286Q
p.E286V
p.E287D
p.F109C
0.73979859
0.73992935
0.73992999
0.73990309
0.73991134
0.739866
0.73989952
0.7398477
0.73790811
0.73676158
0.73993565
0.73758584
0.73644711
0.73674552
0.73650341
0.73801977
0.73835483
0.73771861
0.73755947
0.73683127
0.73725951
0.73756511
0.73681525
0.73657199
0.73682142
0.73757876
0.73656484
0.73762472
0.73793314
0.73740206
0.73744179
0.73944582
0.73950486
0.73972086
0.73944381
0.73941679
0.73973285
0.73972298
0.73997548
0.73997854
0.73998977
0.73998112
0.73998263
0.73997477
0.73997529
0.73999009
0.7399824
0.74008994
0.74003253
0.7402447
0.73997668
0.74017278
0.73997055
0.74000908
0.73997871
0.74014231
0.74005099
0.73996981
0.73998965
0.74006329
Negative set
Positive set
Histogram of values distribution ​​for the magnitude of enthalpy change
Сalculated data obtained for two sets of mutations. Plot of differential entropy changes for two data sets. The graphs show Boolean variables and the final Fechner coefficient.
Made on
Tilda