Correction to the Peptide Linkage

The residues in peptides are joined together by peptide linkages, -HNCO-. These linkages are almost flat, and normally adopt a trans configuration, the hydrogen and oxygen atoms being on opposite sides of the C-N bond. Experimentally, the barrier to interconversion in N-methyl acetamide is about 14 kcal/mol, but all four methods within MOPAC predict a significantly lower barrier, PM3 giving the lowest value.

The low barrier can be traced to the tendency of semiempirical methods to give pyramidal nitrogens. The degree to which pyramidalization of the nitrogen atom is preferred can be seen in the series of compounds given in the Table.

 

Table: Comparison of Observed and Calculated Pyramidalization of Nitrogen

Compound

MNDO

AM1

PM3

PM6

Exp

Ammonia

Py

Py

Py

 

Py

Aniline

Py

Py

Py

 

Py

Formamide

Py

Flat

Py

 

Flat

Acetamide

Py

Flat

Py

 

Flat

N-methyl formamide

Py

Flat

Py

 

Flat

N-methyl acetamide

Flat

Flat

Py

 

Flat

k

6.17

3.32

7.19

   

To correct this, a molecular-mechanics correction has been applied. This consists of identifying the -R-HNCO- unit, and adding a torsion potential of form:

begin{displaymath}k sin ! ^{2} thetaend{displaymath}

where θ is the X-N-C-O angle, X=R or H, and k varies from method to method. Values of k in kcal/mol are also given in the Table. This has two effects: there is a force constraining the nitrogen to be planar, and the HNCO barrier in N-methyl acetamide is raised to 14.00 kcal/mol. When the MM correction is in place, the nitrogen atom for all methods for the last three compounds shown above is planar. The correction should be user-transparent.