Single Molecule Thermodynamics Hypothesis of Protein Folding and Drug Design

Share and Cite:

Yi, F. (2019) Single Molecule Thermodynamics Hypothesis of Protein Folding and Drug Design. Journal of Biosciences and Medicines, 7, 164-172. doi: 10.4236/jbm.2019.711015.

1. Introduction: A Single Molecule View of Protein Folding

To resolve the protein folding problem, that is: predicting the native structure and describing the folding dynamics, we must work with the fundamental physical law that directly governs protein folding process. That law is the Thermodynamic Principle of Protein Folding [1], it is just the Second Law of Thermodynamics since Anfinsen and others already shown that the folding process is spontaneous. In the protein folding case, the second law is that the Gibbs free energy achieves a minimum at the native structure.

Therefore, we have to figure out what is the Gibbs free energy. The question is, is there a Gibbs free energy function whose variables are all possible conformations of a given protein molecule? Or is it only a Gibbs free energy difference between the folded ensemble of protein molecules and its counterpart, the unfolded ensemble? The former is a single molecule view coming from contemplating a protein molecule comes out of ribosome and changes its conformations until it achieves its native structure; the latter is an ensemble view coming from staring at a tube of purified protein solution and trying to figure out the collective behaviours of the protein molecules in the solution while the solution is going towards equilibrium.

Via quantum statistics applied to a tiny thermodynamic system S X which is tailor made for the conformation X and its immediate physiological environment E U , we have derived the conformational Gibbs free energy (CGFE) function G ( X ; E U , U ) for globular proteins U [2] [3] [4].

Applying CGFE function we translate the thermodynamic principle of protein folding into the Single Molecule Thermodynamic Hypothesis (SMTH) of Protein Folding: putting a protein molecule U in an environment E , a stable conformation X E (may not be unique) of U must be a minimizer (local or global) of the CGFE function G ( X ; E , U ) . In particular, the gradient vanishes at X E , G ( X E , E , U ) = 0 .

Another big question in protein folding is that is there a folding force? Leventhal in 1969 [5] has shown by contradiction that there must be a folding force, otherwise if the folding process were only random, it would have taken a time span longer than the Earth’s age.

The CGFE function also gives us the deterministic part of the folding force F i acting on an atom a i of U . It is F i ( X ) = x i G ( X ; E U , U ) .

A scientific hypothesis has to give verifiable predictions to let people confirm or refute it. We suggest two verifiable predictions of the single molecule thermodynamic hypothesis.

1) ab initio predictions of native structures of globular proteins: the native structure X U of the protein U is a (local or global) minimizer of the CGFE function:

G ( X U ; E U , U ) = min X U X U G ( X ; E U , U ) , forsomeneighbourhood U of X U . (1)

Especially, since G ( X ; E U , U ) is smooth, G ( X U ; E U , U ) = 0 . In case X U is a global minimizer, then

G ( X U ; E U , U ) = min X X U G ( X ; E U , U ) . (2)

Here X denotes a conformation and X U is the set of all possible conformations of U .

Therefore, for a globular protein U (we know the formula of G ( X ; E U , U ) ), the prediction of native structure is reduced to a pure mathematical problem, the minimization problem of a known smooth function. A mathematical theorem guarantees that this problem has a solution, i.e., mimimizers always exist, the real task is to find them and determine which is the native structure, a hard programming problem.

2) For a globular protein U , starting from any initial conformation X 0 , there is a folding path X ( t ) = X ( t ; X 0 ) satisfying X ( t 0 ; X 0 ) = X 0 , and for t t 0 , the following Langevin equation:

m i d 2 x i ( t ) d t 2 = F total = x i G ( X ( t ) ; E U , U ) η i d x i ( t ) d t + F i ( t ) , i = 1, , n . (3)

Here η i is the solvent friction. The random force F i ( t ) is caused by occasionally bumping into another non-solvent molecule. Because of it, we do not have a completely deterministic folding path. Again, mathematical theorems guarantee that such folding path X ( t ; X 0 ) exists. Moreover, mathematics also tells us that it is highly depending on its initial conformation X 0 , so an important issue is to know the protein’s initial conformation as it is out of ribosome.

If the two predictions are positively verified, then we can say that theoretically the protein folding problem is resolved, at least for globular proteins.

2. The CGFE Function for Globular Proteins

To explain our CGFE function, we start with a conformation X . A protein U consists of n atoms ( a 1 , , a n ) , a conformation of U can be expressed as a point of 3n-dimensional space 3 n , X = ( x 1 , , x n ) , where x i 3 is the nuclear position of a i in the 3-dimensional space 3 . Not all points in 3 n are conformations of U , bond lengths, bond angles, and van der Waals distances in general, are natural constraints. Denote X U as the set of all possible conformations of U . A conformational function of U is a function f : X U . For example, all force fields used in molecular dynamics simulations are conformational functions.

The 3-dimensional conformation P X of U is P X = i = 1 n C i ( x i ) 3 , where C i 3 is the shape of the atom a i in U , C i ( x i ) indicates it has been congruently moved to the nuclear position x i . C i exists, will change with X , see [6]. A very good approximation to C i ( x i ) is a solid ball B ( x i , r i ) , centred at x i 3 with radius r i that is the a i ’s van der Waals radius. For simplicity, we adopt that C i ( x i ) = B ( x i , r i ) .

In natural, and even in most of artificial environments of protein folding, the immediate environment of a globular protein U is just one layer of water molecules surrounding the conformation P X = i = 1 n B ( x i , r i ) . This is true even the protein molecule is inside a crystal [7].

P X plus the one layer water molecules consists of a tiny thermodynamic system S X , tailor made for P X . As an open thermodynamics system, S X has a Gibbs free energy G ( S X ) . The CGFE function G ( ; E U , U ) : X U then is defined as G ( X ; E U , U ) = G ( S X ) .

Between P X and the layer of water molecules, is an interface M X , for example, the solvent accessible surface P X . The expression of G ( S X ) is via global geometric features of M X and its surface chemical potentials. A protein molecule has many moieties or atom groups, some are charged, some are polar, others are non-polar. They can be classified into hydrophobicity classes H i , 1 i H , H > 1 , from most hydrophobic (non-polar) to the most hydrophilic (polar or charged). An atom a i H j if it belongs to a moiety of class H j . Define P X , i = a j H i B ( x j , r j ) , then P X = i = 1 n P X , i . The space containing water molecules in S X , the ring X = S X \ P X ¯ , is decomposed into H parts (not necessarily connected) via the distance function d i s t ( x , P X , i ) = min y P X , i | x y | , see Figure 1.

X , i = { x X : d i s t ( x , P X , i ) d i s t ( x , P X , j ) , forany j i } , i = 1, , H . (4)

The interface M X then is decomposed accordingly

M X = i = 1 H M X , i , M X , i = M X X , i . (5)

A water molecule in X , i will touch M X , i , so it will be attracted ( H i is charged or polar) or repulsed ( H i is non-polar) by P X . So the same water molecule in different X , i has different chemical potentials μ i . For non-polar H i , μ i > 0 , for charged or polar H i , μ i < 0 . Thus there is a 1 < k < H , such that

μ 1 > μ 2 > > μ k > 0 > μ k + 1 > > μ H . (6)

is called hydrophobic surface, is called hydrophilic surface.

Since water molecules and electrons can enter or leave, is an open thermodynamic system. Hence, the ensemble used here is the grand canonical ensemble. Applying quantum statistics to the open system, the number of water molecules in turns to be a Hermitian operator with mean value. The same is true to the number of electrons. The Gibbs free energy of then is


Since is fixed, here the kinetic energy operator vanishes.

Figure 1. Water molecules are contained inside, note it is not necessarily connected. is an interface.

Since every water molecule in has contact with the surface, is proportional to the area. Therefore, there are, such that


Let be the coordinates of one electron in, the coordinates of water molecules in. Then the electronic density distribution function ((1.3) of ([6], page 6)) gives


The is obtained by mean value theorem of integrals.

Let be the domain enclosed by (), we have roughly the volume, where is the diameter of a water molecule and is the area of. Then taking the mean in (9) we have


Let and. Substitute (8) and (10) into (7), the conformational Gibbs free energy function of a globular protein in its physiological environment is


The function is smooth, i.e., the first and second derivatives exist in except at points such that for some. But such case cannot happen for a conformation, because of the van der Waals distances must be positive. Therefore, on, is smooth.

3. Explanations of Protein Folding

A scientific hypothesis has to be able to explain natural and artificial phenomena. We will explain several phenomena of protein folding, unfolding, and docking, and suggest an application to drug design, according to the SMTH.

3.1. What the CGFE Function in (11) Reveal?

It is well known ([8] and [9]) that native structures of globular proteins have three important global geometric features. Comparing to unfolded conformations, they have: 1) Smaller volume; 2) Smaller surface area; and 3) Compactly packed hydrophobic cores.

Hence folding towards native structure, the volume and area are going to shrink. Then the first two terms in (11) must have positive coefficients, i.e., , otherwise shrinking volume and area would have enlarged the Gibbs free energy.

As the hydrophobic core, look at (11), the hydrophobic surfaces have positive chemical potentials, thus shrinking will reduce the Gibbs free energy. On the other hand, for hydrophilic surface, , therefore, enlarging will reduce the Gibbs free energy. Enlarging of hydrophilic surface area is equivalent to shrinking the hydrophobic surface area, since is also shrinking towards the native structure.

According to the SMTH, predicting native structure of a globular protein is to minimize the CGFE function (11) as in (1) and (2). As analysed above, it is essentially making an ever better hydrophobic core by shrinking volume, area, and hydrophobic area simultaneously and cohesively. This is the “cooperativity” searched in [10], of “the concurrent participation of different regions of the biomolecule to promote and sustain intramolecular or intermolecular interactions”.

One may ask that where are hydrogen bonds in (11)? The answer is that secondary structures and hydrogen bonds are products of minimizing the CGFE function to find the native structures. A judicious examination of an exhaustive PDB sample of small soluble globular proteins of moderate size (residues) showed that the hydrophobic collapsing is coupled with backbone hydrogen-bond formation [11]. In [12], we neglected the volume and area, only shrank the hydrophobic surface area (equivalent to make better hydrophobic core), hydrogen bonds, secondary structures such as helices, strands, and turns, duly appeared with statistical significance.

The explanation is: proteins in their physiological environment are special among polymers. Polymers do not have specified structures, proteins have native structures in their physiological environment. Why? Globular proteins’ peptide chains are special, folding in their physiological environment, while collapsing to hydrophobic cores the residues are putting just in places to be able to form secondary structures and hydrogen bonds simultaneously. Evolution selects the very few peptide chains to be foldable globular proteins. In fact, randomly picking a 400 residue peptide chain, the probability that it is a protein’s peptide chain is at most 10−460 [2]. In any computer, 10−460 is zero.

3.2. Explanation of Denaturation

According to the pioneer research of denaturation [13], all known denature phenomena are caused by environment change from to some. Denaturation, or unfolding, is the same as folding, only in a different environment and with a different CGFE function. According to the SMTH, the unfolding will end at minimizers’s (may not be unique, local or global), of.

Experiments show that the difference between folding and unfolding is that folding leads to a unique native structure, unfolding leads to many different stable conformations [1].

Either in folding or unfolding, a conformation moves along a folding (unfolding) path which satisfies an equation of motion, the Langevin Equation (3), with deterministic forces and. This may explain why folding leads to a unique native structure and unfolding leads to many stable conformations. The random forces in the Langevin equation for folding and unfolding may be also different, in a denatured environment randomly bumping upon some other molecules will happen more often.

Moreover, the initial conformation of the folding or unfolding path also determines where the path ends. Any local minimizer has a domain of attractive basin (U in (1)) such that any initial conformation in will fold to. Because of evolutional selection for a protein’s peptide chain made it fit in the protein’s physiological environment, the attractive basin of is large enough to contain all the initial conformations freshly come out of ribosome, thus even has more than one minimizer, the native structure is the unique folding result.

But is not in denaturation, the initial conformation is the unique native structure? Why it unfolds to many stable denatured conformations? It is because catastrophe, a phenomenon often happens in nature, a description of it is given in [14]: “Catastrophe theory is concerned with the mathematical modelling of sudden changes―so called ‘catastrophes’―in the behaviour of natural systems, which can appear as a consequence of continuous changes of the system parameters”.

Actually, from physiological environment to the denaturation environment, there must be a family of environments connecting them. Except and, these environments are not in equilibrium or quasi-static so the function is not well defined. Hence, although the parameter t varies continuously, catastrophe does happen so that various copies of the same native structure suddenly changed to different structures and when the environment finally changes to the denaturation environment, these changed structures became different initial conformations of denaturation paths under.

3.3. Explanation of Docking

Docking is trying to bound two molecules to form a stable complex.

Let and be two molecules (proteins or others), and and. The 3-dimensional conformations and are contained in their tailor made thermodynamic systems and respectively in their common physiological environment. In Figure 2, in the beginning,. Now suppose that there is a congruence brings to, i.e.,. Van der Waals repulse tells us that it is always 1) and 2). Thus, we define and

Figure 2. (a) Two independent molecules; (b) Invading each other; (c) Forming new minimiser.

binding if and only if 1) and 2).

If and binding, there is a binding energy depending on the net effect of the way of bring to close to. Neglecting, consider the Gibbs free energy for the conformation of the “molecule”), then


If is just water, then is given by (11) for the conformation of the “molecule”. If and are not binding, then, thus (12) is true for any.

In general, even and are minimizers of and respectively, will not be a minimizer of .

3.4. Drug Design

We address a question asked in [10], “the rational drug designer faces a many-body problem: the interactions between the protein target and the drug/ligand involve more than groups matched up in a pairwise fashion at the target-ligand interface... what sort of many-body problem is the drug designer facing and how can this knowledge play advantageously to address the major therapeutic imperatives of today and tomorrow”?

In fact, a drug is much smaller than its target, a globular protein. Thus, the binding energy will be easier to figure out than that between two large proteins. Then the stable structures of the drug/target complex will be (local or global) minimizers of in (12). Thus the drug’s efficiency and safety depend on the properties of minimizers of . This might be an answer to what is the “many-body problem” in [10] and also suggests a particular solution.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.


[1] Anfinsen, C.B. (1973) Principles That Govern the Folding of Protein Chains. Science, 181, 223-230.
[2] Fang, Y. Gibbs Free Energy Formula for Protein Folding. In: Morales-Rodriguez, R., Ed., Thermodynamics—Fundamentals and Its Application in Science, 47-82.
[3] Fang, Y. (2014) A Gibbs Free Energy Formula for Protein Folding Derived from Quantum Statistics. Science China, Physics, Mechanics & Astronomy, 57, 1547-1551.
[4] Fang, Y. (2015) Thermodynamic Principle Revisited: Theory of Protein Folding. Advances in Bioscience and Biotechnology, 6, 37-48.
[5] Levinthal, C. (1969) How to Fold Graciously. Mössbauer Spectroscopy in Biological Systems Proceedings. Proceedings of a Meeting Held at Allerton House, Monticello, Illinois. University of Illinois Bulletin, 67, 22-26.
[6] Bader, R.F.W. (1990) Atoms in Molecules: A Quantum Theory. Clarendon Press, Oxford.
[7] Lattman, E.E. and Loll, P.J. (2008) Protein Crystallography: A Concise Guide. The Johns Hopkins University Press, Baltimore.
[8] Novotny, J., Bruccoleri, R. and Karplus, M. (1984) An Analysis of Incorrectly Folded Protein Models. Implications for Structure Predictions. J. Mol. Biol., 177, 787-818.
[9] Novotny, J., Rashin, A.A. and Bruccoleri, R. (1988) Criteria That Discriminate between Native Proteins and Incorrectly Folded Models. Proteins, 4, 19-30.
[10] Fernádez Stigliano, A. (2015) Biomolecular Interfaces: Interactions, Functions and Drug Design. Springer International Publishing.
[11] Fernandez, A., Kardos, J. and Goto, Y. (2003) Protein Folding: Could Hydrophobic Collapse Be Coupled with Hydrogen-Bond Formation? FEBS Letters, 536, 187-192.
[12] Fang, Y. and Jing, J. (2010) Geometry, Thermodynamics, and Protein. Journal of Theoretical Biology, 262, 382-390.
[13] Wu, H. (1931) Studies of Denaturation of Proteins XIII. A Theory of Denaturation. Chinese J. Physiol., 5, 321-344.
[14] Sanns, W. (2009) Catastrophe Theory. In: Meyers, R.A., Ed., Encyclopedia of Complexity and System Science, Springer, Vol. 4, 703-719.

Copyright © 2023 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.