Posts Tagged ‘energy’

Conformational analysis and enzyme activity: models for amide hydrolysis.

Sunday, April 12th, 2009

The diagram below summarizes an interesting result recently reported by Hanson and co-workers (DOI: 10.1021/jo800706y. At ~neutral pH, compound 13 hydrolyses with a half life of 21 minutes, whereas 14 takes 840 minutes. Understanding this difference in reactivity may allow us to understand why some enzymes can catalyze the hydrolysis of peptides with an acceleration of up to twelve orders of magnitude.

Models for peptide cleavage.

Models for peptide cleavage.

The secret to understanding this behaviour lies in a technique known as conformational analysis, for which Derek Barton was awarded a Nobel prize. Indeed, the very molecules for which he first developed his technique were the decalins, of which molecule  13 is an example of a cis-decalin and 14 a trans-decalin. Barton’s insight was to recognize that both types of ring prefer to exist in chair conformations rather than the alternative boat shape.

The technique pioneered by Barton for estimating the energies of these various conformations is called Molecular Mechanics, and can be used to explain the difference in reactivity. Considering first molecule 13, one can calculate its molecular mechanics energy for two conformations, differing in whether the N-alkyl sidechain is equatorial (left) or axial (right).

Cis amide

Cis amide. Click for Equatorial 3D.

The equatorial form (green box) comes out about 5 kcal/mol lower in energy than the axial (red box). One can also calculate the energy of the product, which arises from the OH attacking the carbon of the amide (dashed lines), evicting ammonia, and forming a cyclic lactone. Here, the most stable product (by ~10 kcal/mol) is again that resulting from the green bond forming. From the simple relationship ΔG = -RT Ln K (where K describes the position of the equatorial/axial equilibrium), one can conclude that the ratio equatorial/axialis ~4000, i.e. the favoured reaction arises from the most abundant reactant.

Trans amide

Trans amide. Click for 3D.

With the trans amide, the equatorial conformation (green box) is around 3 kcal/mol lower than the axial (red box), but now the most stable lactone product (by ~ 3 kcal/mol) arises (green bond) from the less stableaxial reactant. For reaction to occur, the equatorial reactant has to first isomerise to the axial, which imposes a ~3 kcal/mol penalty on the reaction. This is enough to slow the rate of the reaction significantly compared to the un-penalised cis-decalin reaction.

 

How do molecules interact with each other?

Sunday, April 12th, 2009

Understanding how molecules interact (bind) with each other when in close proximity is essential in all areas of chemistry. One specific example of this need is for the molecule shown below.

The Pirkle reagent

The Pirkle reagent

This is the so-called Pirkle Reagent and is much used to help resolve the two enantiomers of a racemic mixture, particularly drug molecules. The reagent binds to each enantiomer of a racemic drug differently, and this difference can be exploited by using e.g. column chromatography to separate the two forms. The conventional wisdom is that such chiral recognition occurs via a three-point binding model. In other words, at least three different interactions must occur between the Pirkle reagent and the drug to allow such chiral recognition.

So how do we identify where these bindings might occur? A good place to start is to look at the self-binding of the Pirkle reagent, in other words, how does it interact with itself in the crystal state? An X-ray structure of the pure enantiomer of the Pirkle reagent shows that it binds with itself to form a loose dimer. We are now in a position to analyze exactly how this binding occurs. To do this, we are going to invoke a technique known as Atoms-in-molecules or AIM. This effectively looks at the curvature of the electron density in the dimer, and from the characteristics of this curvature, identifies a series of so called critical points, or regions where the first derivative of the electron density (referred to as ρ(r) ) with respect to the geometry is zero. These critical points come in four varieties only;

  1. A nuclear critical point, which almost exactly corresponds to where the nuclei are
  2. A bond critical point, which is the key to understanding not only where actual bonds are in the molecule, but also a range of weaker interactions which are conventionally not graced with the term bond, but which nevertheless can be essential in understanding how to molecules interact weakly with each other.
  3. The remaining two types of critical point relate to rings and cages, and we will not be concerned further with them here.

The electron density required for this analysis could in principle come from the X-ray measurements themselves, but it is not easy to acquire this to the required accuracy (although it can be done). In this case, it is easier (and probably no less accurate) to calculate the density from a DFT-based quantum mechanical calculation. The result of this is shown below.

Pirkle dimer. Click on image to obtain model

Pirkle dimer. Click for 3D.


The light blue spheres show the position of selected bond critical points or BCPs in the AIM analysis. So what do they tell us about how two molecules of Pirkle molecule interact with each other? Three different points labelled 1-3 are highlighted for discussion.

  1. Points 1 connect the hydrogen of the OH group with the carbons of the π-face of the anthracene ring (the left ring of the molecule as shown above). This is an unusual type of interaction known as a π-facial hydrogen bond, and it has only been recognized as such in the last 30 years. Note that this interaction is not restricted to occur just between a pair of atoms, but can involve more (in this case almost a whole benzene ring). By finding the value of the electron density ρ(r) at this BCP, one can estimate the energy of interaction resulting from its formation. In this case, ρ(r) ~ 0.014 au, and comparison with other types of hydrogen bond suggests that this value corresponds to an interaction energy of around 2.5 kcal/mol. This is a little weaker than a conventional OH…O hydrogen bond, but is still quite significant. Two of these interactions occur in this Pirkle dimer.
  2. Points 2 are equally unexpected. They connect the oxygen of the same OH group involved in the previous interaction, and one of the ring C-H groups. Again, that C-H…O groups can interact has only been recognized relatively recently. The value of ρ(r) of ~ 0.018 indicates a hydrogen bond strength of ~3 kcal/mol, again hardly insignificant.
  3. There are four specific interactions of the final type 3. These occur in the region of overlap of the two anthracene rings, and these are referred to as π-π stacking interactions. Again, the ρ(r) of ~ 0.005, calibrated against known systems, suggests that each is individually worth around 1 kcal/mol.

So adding up all eight interactions indicates that the two molecules of the Pirkle reagent have an interaction energy of around 15 kcal/mol resulting just from these weak bonds (there are other types of interactions between two molecules known as dispersion forces, which also contribute), and which together provide more than enough free energy to overcome the entropy required to bring the two molecules together.

Armed with tools such as AIM, one can now be more confident in analyzing the various terms that contribute to two molecules interacting with each other, and in the case of chiral molecules, how these interactions may result in chiral recognitions.

The SN-1 Reaction live!

Friday, April 3rd, 2009