Molecular self-assembled helix peptide nanotubes based on some amino acid molecules and their dipeptides: molecular modeling studies.

Abstract

Context: The paper considers the features of the structure and dipole moments of several amino acids and their dipeptides which play an important role in the formation of the peptide nanotubes based on them. The influence of the features of their chirality (left L and right D) and the alpha-helix conformations of amino acids are taken into account. In particular, amino acids with aromatic rings, such as phenylalanine (Phe/F), and branched-chain amino acids (BCAAs)-leucine (Leu/L) and isoleucine (Ile/I)-as well as corresponding dipeptides (diphenylalanine (FF), dileucine (LL), and diisoleucine (II)) are considered. The main features and properties of these dipeptide structures and peptide nanotubes (PNTs), based on them, are investigated using computational molecular modeling and quantum-chemical semi-empirical calculations. Their polar, piezoelectric, and photoelectronic properties and features are studied in detail. The results of calculations of dipole moments and polarization, as well as piezoelectric coefficients and band gap width, for different types of helical peptide nanotubes are presented. The calculated values of the chirality indices of various nanotubes are given, depending on the chirality of the initial dipeptides-the results obtained are consistent with the law of changes in the type of chirality as the hierarchy of molecular structures becomes more complex. The influence of water molecules in the internal cavity of nanotubes on their physical properties is estimated. A comparison of the results of these calculations by various computational methods with the available experimental data is presented and discussed.

Method: The main tool for molecular modeling of all studied nanostructures in this work was the HyperChem 8.01 software package. The main approach used here is the Hartree-Fock (HF) self-consistent field (SCF) with various quantum-chemical semi-empirical methods (AM1, PM3, RM1) in the restricted Hartree-Fock (RHF) and in the unrestricted Hartree-Fock (UHF) approximations. Optimization of molecular systems and the search for their optimal geometry is carried out in this work using the Polak-Ribeire algorithm (conjugate gradient method), which determines the optimized geometry at the point of their minimum total energy. For such optimized structures, dipole moments D and electronic energy levels (such as E_HOMO and E_LUMO), as well as the band gap E_g = E_LUMO - E_HOMO, were then calculated. For each optimized molecular structure, the volume was calculated using the QSAR program implemented also in the HyperChem software package.