Ab Initio and Empirical Model MD Simulation Studies of Solvent Effects on the Properties of N-Methylacetamide along a cis−trans Isomerization Pathway

Density Functional Studies on Secondary Amides: Role of Steric Factors in Cis/Trans Isomerization

Cis/trans isomerization of amide bonds is a key step in a wide range of biological and synthetic processes. Occurring through C-N amide bond rotation, it also coincides with the activation of amides in enzymatic hydrolysis. In recently described QM studies of cis/trans isomerization in secondary amides using density functional methods, we highlighted that a peptidic prototype, such as glycylglycine methyl ester, can suitably represent the isomerization and complexities arising out of a larger molecular backbone, and can serve as the primary scaffold for model structures with different substitution patterns in order to assess and compare the steric effect of the substitution patterns. Here, we describe our theoretical assessment of such steric effects using tert-butyl as a representative bulky substitution. We analyze the geometries and relative stabilities of both trans and cis isomers, and effects on the cis/trans isomerization barrier. We also use the additivity principle to calculate absolute steric effects with a gradual increase in bulk. The study establishes that bulky substitutions significantly destabilize cis isomers and also increases the isomerization barrier, thereby synergistically hindering the cis/trans isomerization of secondary amides. These results provide a basis for the rationalization of kinetic and thermodynamic properties of peptides with potential applications in synthetic and medicinal chemistry.

Exploring Peptide⁻Solvent Interactions: A Computational Study

The dilemma of reconciling the contradictory evidence regarding the conformation of long solvated peptide chains is the so-called "reconciliation problem". Clues regarding the stability of certain conformations likely lie in the electronic structure at the peptide⁻solvent interface, but the peptide⁻solvent interaction is not fully understood. Here, we study the influence of aqueous solvent on peptide conformations by using classical molecular dynamics (MD) and quantum mechanical/molecular mechanical (QM/MM) energy calculations. The model systems include an 11-residue peptide, X 2 A 7 O 2 (XAO), where X, A, and O denote diaminobutyric acid, alanine, and ornithine, respectively, and a 9-mer (Arg-Pro-Pro-Gly-Phe-Ser-Ala-Phe-Lys). Spectroscopic and MD data present conflicting evidence regarding the structure of XAO in water; some results indicate that XAO adopts a polyproline II (P II ) conformation, whereas other findings suggest that XAO explores a range of conformations. To investigate this contradiction, we present here the results of MD simulations of XAO and the 9-mer in aqueous solution, combined with QM/MM energy calculations.

Solvation of Biomolecules by the Soft Sticky Dipole-Quadrupole-Octupole Water Model

The soft sticky dipole-quadrupole-octupole (SSDQO) potential energy function represents a water molecule by a single site with a van der Waals sphere and point multipoles. Previously, SSDQO was shown to give good properties for liquid water and solvation of simple ions and is faster than three point models. Here, SSDQO is assessed for solvating biologically relevant molecules having a multi-site, partial charge description. Monte Carlo simulations of ethanol, benzene, and N-methylacetamide in SSDQO with SPC/E moments showed the water structure was as good as in SPC/E. Thus, SSDQO is potentially useful for simulations of biological macromolecules in aqueous solution.

Estimating the "steric clash" at cis peptide bonds

To account for the scarcity of cis peptide bonds in proteins, especially in nonproline (or secondary amide) cases, a steric-clash argument is often put forward, in a scheme where the R lateral chains are facing parallel one another, and the backbone is kept in an "all- trans"-like arrangement. Although such a steric conflict can be partly relieved through proper adjustment of the backbone dihedral angles, one can try to estimate its associated energy cost. To this end, quantum-chemistry approaches using a differential-torsion protocol and bond-separation-energy analyses are applied to N-ethyl propionamide CH3-CH2-CO-NH-CH2-CH3, regarded as a model capable of exhibiting C beta...C beta interaction as in alanine succession. The calculations provide an increment of 9 kcal/mol, quite close to that obtained in the nearly isostere (gsg) rotamer of n-hexane (10 kcal/mol), suggesting the local effects induced by methyl-methyl contact are similar in both cases. Analogous treatments on larger radicals as encountered in leucine or phenylalanine dimers do not change this increment much, which therefore defines the basic reference per-plaque quota to be overcome along all- cis chains. Explicit modeling indicated it can be reduced by up to a factor of 4.