3.1: Potential Energy Surface and Bonding Interactions
- Page ID
- 398272
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)In this chapter, we return to the idea of a potential energy surface and consider models for a covalent bond. The potential energy changes as a function of the relative distance between bonded atoms. As atoms come together to form molecules, the potential energy will have a minimum at the equilibrium bond distance. For diatomic molecules, the potential energy surface is one-dimensional and depends only on the distance between the two atoms in the diatomic molecule. For polyatomic molecules with more complex structure, the potential energy surface will be a high-dimensional function of the positions all the atoms. To model the secondary structure of proteins and other macromolecules, we need potential energy functions to describe the equilibrium bond stretching of covalent bonds, the equilibrium angle bending (involving an angle between three atoms), and dihedral torsional rotations (involving four atoms). In this chapter we will look at empirical potential energy functions for these three type of interactions that govern the conformational energy of a polypeptide chain.
- Understand the difference between a harmonic potential energy and the Morse potential, and what is meant by the terms harmonic or anharmonic vibrations.
- Be able to explain and justify the different parameters used to describe single, double, and triple bonds.
- Be able to sketch the dihedral potential energy for the butane molecule.
- Understand the physical justification for the different types of empirical potential energy functions that describe bonding interactions.
The chemical bond
The starting place for understanding chemical bonding and molecular structure is the Born Oppenheimer approximation. Since atomic nuclei are much heavier than the electrons, the electrons move much faster than the nuclei. Thus, to a good approximation, the motion of the atomic nuclei and electrons can be treated separately. The energy due to the electrons in the ground state can be calculated assuming a fixed position of the nuclei. This calculation can be repeated at different values of the internuclear distance between the atoms. This gives the potential energy for the nuclear positions as a function of the internuclear distance. As a consequence of the Born Oppenheimer approximation, we can think of the atomic nuclei as moving along a potential energy surface that is a function only of the nuclei position.
We begin by considering the simplest case of the potential energy curve for a diatomic molecule, such as H2 that was introduced in Chapter II.3 in the context of transition state theory. Figure III.1.A again shows the potential energy curve for the H2 molecule. Notice that the equilibrium bond length, labeled \(r_e\), is located at the minimum of the bond potential energy curve.
![Figure3.1A.png](https://chem.libretexts.org/@api/deki/files/415004/Figure3.1A.png?revision=2&size=bestfit&width=505&height=340)
If the displacement away from the equilibrium bond length is small, the potential energy surface in the vicinity of the equilibrium bond distance can be described by a harmonic oscillator (see Appendix for the derivation)
\[V=\frac{1}{2}k \left(r-r_e\right)^2\label{EQ:harmonic_bond}\]
In Equation \ref{EQ:harmonic_bond}, \(r\) is the distance between the two atoms and \(r_e\) is the equilibrium bond distance. Figure III.1.B compares the full potential energy curve for H2 with Equation \ref{EQ:harmonic_bond}. Notice that at sufficiently small displacements from \(r_e\), vibrations in the bond are accurately described by Equation \ref{EQ:harmonic_bond}. Oscillations that are described by Equation \ref{EQ:harmonic_bond} are called harmonic. In a harmonic bond, the force that results from stretching or compressing the bond distance is proportional to the displacement from the equilibrium position (Hook’s law).
\[\begin{eqnarray}\text{Force} &=& – \frac{dV}{dr} \\[4pt] &=& – k \left(r-r_e\right)\label{EQ:harmonic_force}\end{eqnarray}\]
![Figure3.1B.png](https://chem.libretexts.org/@api/deki/files/419611/Figure3.1B.png?revision=1&size=bestfit&width=451&height=303)
Equations III.1.\ref{EQ:harmonic_bond} and III.1.\ref{EQ:harmonic_force} describe small perturbative fluctuations about the equilibrium bond distance. In this case, the bonds are said to be harmonic bonds. Large displacements of the bond away from \(r_e\) lead to anharmonic vibrations. Anharmonic vibrations describe bond stretching that is far from the equilibrium bond distance. To model anharmonic vibrations, we could consider including higher-order terms in the Taylor expansion in Appendix. Alternatively, a common potential that can model bond breaking is the Morse potential:
\[V = D_e \left(1 – e^{-a (r-r_e)}\right)^2 \label{EQ:harmonic_bond2}\]
where \(D_e\) is the bond dissociation energy, and \(a\) controls the potential width. Figure III.1.C shows a plot of the Morse potential.
![Figure3.1C.png](https://chem.libretexts.org/@api/deki/files/421563/Figure3.1C.png?revision=1&size=bestfit&width=612&height=411)
The potential energy surface for macromolecules is high-dimensional
In the previous section we described the potential energy for a diatomic molecule, such as H2, O2, N2, NO, etc …, in terms of the equilibrium bond distance \(\bf{r_e}\) and the internuclear distance \(\bf{r}\). As we saw in Chapter II.3 (see Figure II.3.C), for polyatomic molecules involving more than two bonded atoms, the potential energy surface is a function of all the internuclear distances and quickly becomes high-dimensional. For a proteins, nucleic acids, lipids, polysaccharides, or other macromolecules, if we let \(\textbf{R}\) be the set of positions of all the atoms ( \( \textbf{R}=\{\textbf{R}_1,\textbf{R}_2, …,\textbf{R}_N \} \) ), then the potential energy surface \(U(\textbf{R})\) will be a function of all atomic coordinates. Although we cannot draw in more than three dimensions, we can consider an abstract high-dimensional space where a point is given by the set of all the atomic positions: \((\textbf{R}_1,\textbf{R}_2, …,\textbf{R}_N)\). Local minima in this high-dimensional space will represent metastable conformations and local maxima represent kinetic barriers.
At constant number of particles, temperature, and volume, the probability of a given conformation is given by the normalize Boltzmann weight of the energy:
\[P(\textbf{R}) = \frac{e^{-U(\textbf{R})/k_BT}}{\int ... \int e^{-U(\textbf{R})/k_BT} \ d\textbf{R}_1, ..., d\textbf{R}_N}\label{EQ:Boltzman}\]
where the integral in the denominator is a multi-dimensional integral over all the atomic coordinates. It is common practice in computational biophysics to reduce the high-dimensional potential energy landscape to a limited number of variable that can describe the conformational state of the system, called order parameters or collective variables. An example of this was already present in Chapter II.3 where we reduce the two-dimensional potential energy surface for the SN2 substitution to a one-dimensional picture by constructing a path along the lowest energy contour from the reactant to product state as shown in Figure III.1.D
![Figure3.1D.png](https://chem.libretexts.org/@api/deki/files/421564/Figure3.1D.png?revision=1&size=bestfit&width=713&height=479)
As another example, we can consider the two dihedral angels, \(\phi\) and \(\psi\) for the Nme-Ala-Ace peptide shown here:
![Figure3.1E.png](https://chem.libretexts.org/@api/deki/files/421565/Figure3.1E.png?revision=1&size=bestfit&width=621&height=417)
Imagine we could sample randomly the equilibrium configurations drawn from a Boltzmann distribution of states given by Equation \ref{EQ:Boltzman}. In this case, we could construct a histogram of all \(\phi\), \(\psi\) angle combinations to obtain the probability \(P(\phi, \psi)\). From statistical biophysics, we can define a free energy surface from the probabilities as
\[F(\phi,\psi) = -k_BT \ln P(\phi,\psi) \label{EQ:FES}\]
The free energy surface is a low-dimensional representation of system that provides information about metastable states, their relative stability, and the free energy barriers separating them. An example two-dimensional free energy surface for the Nme-Ala-Ace peptide is shown in Figure III.1.F. From this surface, we can define the transition state as the separatrix where there is equal probability for the system to flow into either metastable basin.
![Fes_Ala.png](https://chem.libretexts.org/@api/deki/files/426807/Fes_Ala.png?revision=1&size=bestfit&width=402&height=270)
The covalent bond revisited
In molecular mechanics, we are often interested in equilibrium fluctuations of chemical bonds. For small fluctuations, it is customary to use a harmonic potential for the form of \ref{EQ:X}:
\[V=\frac{1}{2}k \left(r-r_e\right)^2 \label{EQ:X}\]
The spring constant \(k\) determines the stiffness of the bond and the value of \(r_e\) is the equilibrium bond distance. These values are parameters of the model and are usually empirically determined to fit to experimentally observed values (for example from IR spectroscopy). Each type of chemical bond will need to be parameterized to agree with experiment. For example, a carbon-carbon single bond, between two sp3 carbon atoms, such as in the ethane molecule shown in Figure III.1.G (a) will have a slightly larger value of \(r_e\) and a slightly smaller value of \(k\) as compared to a carbon-carbon double bond, between two sp2 carbon atoms, such as in the ethene molecule shown in Figure III.1.G (b).
![Figure3.1G.png](https://chem.libretexts.org/@api/deki/files/421566/Figure3.1G.png?revision=1&size=bestfit&width=507&height=341)
Table III.1.a shows the values for the spring constant and equilibrium bond distance for the different bond orders of carbon bonds.
hybridization | bond order | spring constant, k | equilibrium bond distance, re |
sp3 | 1 | 100 kcal/mol/Å2 | 1.5 Å |
sp2 | 2 | 200 kcal/mol/Å2 | 1.3 Å |
sp | 3 | 400 kcal/mol/Å2 | 1.2 Å |
Figure III.1.H shows a plot of the harmonic bond for each type of bond (single, double, and triple). We can see that as the bond order increases, the equilibrium bond distance, \(r_e\), decreases, and the stiffness of the bond, \(k\) increases.
![Figure3.1H.png](https://chem.libretexts.org/@api/deki/files/421568/Figure3.1H.png?revision=1&size=bestfit&width=416&height=280)
Angular and dihedral potentials for polyatomic molecules
In molecules with more than two bonded atoms it is possible to define a bond angle shown in Figure III.1.I. A bond angle is formed by three consecutive atoms, here labeled \(i\), \(j\), and \(k\). The corresponding bond angle is \(\theta\). If we let \(\textbf{r}_{ij}\) be the bond distance between atoms \(i\) and \(j\) and \(\textbf{r}_{jk}\) the bond distance between atoms \(j\) and \(k\), the bond angle can be expressed in terms of the atomic positions as
\[\theta = \arccos \left(\frac{\textbf{r}_{ij}\cdot \textbf{r}_{jk}}{|r_{ij}||r_{jk}|} \right)\label{EQ:Bondangle}\]
where the dot indicates the dot product between bond vectors \(\textbf{r}_{ij}\) and \(\textbf{r}_{jk}\) and \(|r_{ij}|\) and \(|r_{jk}|\) are the lengths of the bond vectors.
It is common to describe fluctuations about the equilibrium bond angle with a harmonic potential of the form:
\[V_{angle} =\frac{1}{2}k \left(\theta-\theta_e\right)^2\label{EQ:angle}\]
where \(\theta_e\) is the equilibrium bond angle and \(k\) is a force constant that parameterizes the stiffness of the bond angle.
![Figure3.1I.png](https://chem.libretexts.org/@api/deki/files/421569/Figure3.1I.png?revision=1&size=bestfit&width=552&height=371)
If we are interested in the force acting on atom \(i\), as we would need for solving Newton’s equation of motion, we can use the chain rule:
\[\frac{dV}{d r_i} = \frac{dV}{d\theta}\cdot \frac{d\theta}{d r_i}\label{EQ:chainrule}\]
where the derivate \(\frac{dV}{d\theta}\) can be computed by taking the derivative of Equation \ref{EQ:angle} with respect to \(\theta\), and the derivative \(\frac{d\theta}{d r_i}\) can be computed by taking the derivative of Equation \ref{EQ:Bondangle} with respect to \(r_i\).
For atoms with four or more bonded atoms, there is a potential energy barrier to internal rotation about the dihedral angles. As a simple example we will consider the rotation about the methyl dihedral in butane as shown below:
The lowest energy configuration is the staggered, anti-, geometry at a bond angle of 180º. There are two equivalent staggered configurations with bond angles of 300º and 60º which are 3.8 kJ mol-1 higher in energy. The large energy barrier of 19 kJ mol-1 for the eclipsed, sys-, geometry means that rotations from on gauche configuration (300º) into another gauche configuration (60º) will be a rare event. To model a periodic potential energy function as seen in Figure III.1.J, the torsional energy is given by:
\[V_{torsion} = A[1 + \cos (n \phi – \phi_0)]\label{EQ:torsion}\]
where \(\phi\) is the dihedral angle of interest and the parameter \(A\) is a constant with units of energy, \(n\) is an integer, and \(\phi_0\) is a reference dihedral angle.
![Figure_3.1J.png](https://chem.libretexts.org/@api/deki/files/426809/Figure_3.1J.png?revision=1&size=bestfit&width=552&height=371)
As an example, consider the the structure of the retinal chromophore shown in Figure III.1.K which has a conjugated \(\pi\) bond system. Absorption of a photon induces an excited state of the chromophore which is followed by the isomerization of the retinal protonated Schiff base around the C13=C14 bond (adjacent to the Schiff base group). Rotation of the dihedral bonds can be modeled using Equation \ref{EQ:torsion}. Table III.1.b reports the parameter set used in molecular dynamics simulations of the retinal Schiff base.
![Figure_retinal.png](https://chem.libretexts.org/@api/deki/files/421570/Figure_retinal.png?revision=1&size=bestfit&width=552&height=371)
Table III.1.b. The parameter set for the torsion potentials of the main polyene chain of the retinal Schiff base. Values taken from Tajkhorshid et al. 1999. The form of the torsion potential is given by Equation \ref{EQ:torsion}.
\(\phi \) | \(A \) kcal/mol | \(n\) | \( \phi_0 \) (degrees) |
C5-C6-C7-C8 | 5.62 | 2.0 | 180.0 |
C6-C7-C8-C9 | 19.99 | 2.0 | 180.0 |
C7-C8-C9-C10 | 8.515 | 2.0 | 180.0 |
C8-C9-C10-C11 | 18.64 | 2.0 | 180.0 |
C9-C10-C11-C12 | 11.25 | 2.0 | 180.0 |
C10-C11-C12-C13 | 17.54 | 2.0 | 180.0 |
C11-C12-C13-C14 | 14.15 | 2.0 | 180.0 |
C12-C13-C14-C15 | 14.73 | 2.0 | 180.0 |
C13-C14-C15-N16 | 15.215 | 2.0 | 180.0 |
C14-C15-N16-C\(_\epsilon\) | 14.38 | 2.0 | 180.0 |
Practice Problems
Coming soon ...