Gas Phase Conformations of Selenocysteine and Related Ions: A Comprehensive Theoretical Study

Identifying the protonation site and the scope of non-proline cis-peptide bond conformations: a first-principles study on protonated oligopeptides

The existence of non-proline cis-peptide bond conformations of protonated triglycine proposed by us has been verified through a recent IR-IR double resonance experiment. However, the scope of such unique structures in protonated oligopeptides and whether protonation at amide oxygen is more stable than that at traditional amino nitrogen remain unsolved. In this study, the most stable conformers of a series of protonated oligopeptides were fully searched. Our findings reveal that the special cis-peptide bond structure appears with high energies for diglycine and is energetically less favored for tetra- and pentapeptides, while it acts as the global minimum only for tripeptides. To explore the formation mechanism of the cis-peptide bond, electrostatic potential analysis, and intramolecular interactions were analyzed. Advanced theoretical calculations confirmed that amino nitrogen is still preferred as the protonated site in most cases except glycylalanylglycine(GAG). The energy difference between the two protonated isomers of GAG is only 0.03 kcal mol^-1, indicating that the tripeptide is most likely to be protonated on the amide oxygen first. We also conducted chemical (infrared (IR)) and electronic (X-ray photoelectron spectra (XPS) and near-edge X-ray absorption fine structure spectra (NEXAFS)) structure calculations of these peptides to identify their notable differences unambiguously. This study thus provides valuable information for exploring the scope of cis-peptide bond conformation and the competition between two different protonated ways.

Theoretical Spectroscopic Studies on Chemical and Electronic Structures of Selenocysteine and Pyrrolysine

The chemical and electronic structures of the 21st and 22nd proteinogenic amino acid selenocysteine (Sec), pyrrolysine (Pyl), and their derivatives (deprotonated and protonated ions) were extensively characterized for the first time. Through the fragment based step-by-step research on their potential energy surface (PES), electronic energies of the most stable conformers of Sec, Pyl and the related ions were finally determined at the advanced CBS-QB3 and DSD-PBEP86-D3(BJ)/aug-cc-pVTZ levels, respectively, with the identification of many new low-energy conformers. The infrared spectra (IR) at 298 K of the most abundant conformers in different forms were scaled by comparison with the anharmonic frequency calculations and analyzed comparing with the experimental spectra of similar molecules. The characteristic soft X-ray spectra (including X-ray photoelectron spectra (XPS) and near-edge X-ray absorption fine-structure spectra (NEXAFS)) of the most stable conformers at 498 K were also simulated. In particular, the two possible protonated configurations of Pyl can be clearly distinguished by their different spectral features. Furthermore, a small binding energy intersection appeared around 293 eV at the C 1s edge between the canonical and protonated Pyl conformers, which is different from all the previous studies. This work thus filled the gap in our knowledge by providing detailed information on the chemical and electronic structures of Sec and Pyl and will be a useful guidance for future experimental research.

An improved two-rotor function for conformational potential energy surfaces of 20 amino acid diamides

Predicting the three-dimensional structure of a protein from its amino acid sequence requires a complete understanding of the molecular forces that influences the protein folding process. Each possible conformation has its corresponding potential energy, which characterizes its thermodynamic stability. This is needed to identify the primary intra- and inter-molecular interactions, so that we can reduce the dimensionality of the problem, and create a relatively simple representation of the system. Investigating this problem using quantum chemical methods produces accurate results; however, this also entails large computational resources. In this study, an improved two-rotor potential energy function is proposed to represent the backbone interactions in amino acids through a linear combination of a Fourier series and a mixture of Gaussian functions. This function is applied to approximate the 20 amino acid diamide Ramachandran-type PESs, and results yielded an average RMSE of 2.36 kJ mol−1, which suggest that the mathematical model precisely captures the general topology of the conformational potential energy surface. Furthermore, this paper provides insights on the conformational preferences of amino acid diamides through local minima geometries and energy ranges, using the improved mathematical model. The proposed mathematical model presents a simpler representation that attempts to provide a framework on building polypeptide models from individual amino acid functions, and consequently, a novel method for rapid but accurate evaluation of potential energies for biomolecular simulations.

Identification of the protonation site of gaseous triglycine: the cis-peptide bond conformation as the global minimum

Extensive ab initio investigations have been performed to characterize stable conformers of protonated triglycine (GGGH) in the gas phase.

Comprehensive analysis of energy minima of the 20 natural amino acids

Energy minima of the 20 natural amino acids (capped by a peptide bond at both the N and the C termini, CH3-C(═O)-N(H)-(H)Cα(R)-C(═O)-N(H)-CH3), were obtained by ab initio geometry optimization. Starting with a large number of minima, quickly generated by MarvinView, geometry optimization at the HF/6-31G(d,p) level of theory reduced the number of minima, followed by further optimization at the B3LYP/apc-1 and MP2/cc-pVDZ levels, which caused some minima to disappear and some stable minima to migrate on the Ramachandran map. There is a relation between the number of minima and the size and the flexibility of the side chain. The energy minima of the 20 amino acids are mainly located in the regions of βL, γL, δL, and αL of the Ramachandran map. Multipole moments of atoms occurring in the fragment [-NH-Cα-C(═O)-] common to all 20 amino acids were calculated at the three levels of theory mentioned above. The near parallelism in behavior of these moments between levels of theory is beneficial toward estimating moments with the more expensive B3LYP and MP2 methods from data calculated with the cheaper HF method. Finally, we explored the transferability of properties between different amino acids: the bond length and angles of the common fragment [-NH-Cα(HαCβ)-C'(═O)-] in all amino acids except Gly and Pro. All bond lengths are highly transferable between different amino acids, and the standard deviations are small.