Dimension reduction in conformational analysis: a two-rotor mathematical model of amino acid diamide conformational potential energy surface

Abstract

The conformational potential energy surface (PES) of a molecule provides insights into the relative stability of the possible foldamers. However, the time and space complexity of electronic structure calculations, commonly used to generate PES, increases exponentially with an increasing number of atoms. The use of mathematical functions to model the topology of conformational PES is an alternative to more computer-intensive quantum chemical calculations, but the choice and complexity of functions used are crucial in achieving more accurate results. This paper presents a method to illustrate the topology of amino acid diamide PESs through a linear combination of a Fourier series and a mixture of Gaussian functions. Results yield a significantly small error, with an average RMSE of 4.9946 kJ mol−1 for all fits, which suggest that these functions may be used to represent the topology of the PESs, with around twofold order of magnitude decrease in computational time, with respect to DFT electronic structure calculations. This study ultimately aims to provide a foundation for a framework on building polypeptide PES from individual amino acid PESs.