1
|
Villar JJS, Valdez ARL, Setiadi DH, Csizmadia IG, Viskolcz B, Rágyanszki A. An improved two-rotor function for conformational potential energy surfaces of 20 amino acid diamides. CAN J CHEM 2018. [DOI: 10.1139/cjc-2017-0571] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
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.
Collapse
Affiliation(s)
- John Justine S. Villar
- Scientific Computing Laboratory, Department of Computer Science, University of the Philippines, Diliman, 1101 Quezon City, Philippines
- Institute of Chemistry, University of Miskolc, 3515 Miskolc-Egyetemváros, Hungary
| | - Adrian Roy L. Valdez
- Scientific Computing Laboratory, Department of Computer Science, University of the Philippines, Diliman, 1101 Quezon City, Philippines
| | - David H. Setiadi
- Department of Chemistry, University of Toronto, Toronto, ON M5S 3H6, Canada
| | - Imre G. Csizmadia
- Institute of Chemistry, University of Miskolc, 3515 Miskolc-Egyetemváros, Hungary
- Department of Chemistry, University of Toronto, Toronto, ON M5S 3H6, Canada
| | - Béla Viskolcz
- Institute of Chemistry, University of Miskolc, 3515 Miskolc-Egyetemváros, Hungary
| | - Anita Rágyanszki
- Institute of Chemistry, University of Miskolc, 3515 Miskolc-Egyetemváros, Hungary
- Department of Chemistry, University of Toronto, Toronto, ON M5S 3H6, Canada
| |
Collapse
|
2
|
Yuan Y, Mills MJL, Popelier PLA, Jensen F. Comprehensive analysis of energy minima of the 20 natural amino acids. J Phys Chem A 2014; 118:7876-91. [PMID: 25084473 DOI: 10.1021/jp503460m] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
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.
Collapse
Affiliation(s)
- Yongna Yuan
- Manchester Institute of Biotechnology (MIB) , 131 Princess Street, Manchester M1 7DN, Great Britain and
| | | | | | | |
Collapse
|
3
|
Pal TK, Sankararamakrishnan R. Quantum chemical investigations on intraresidue carbonyl-carbonyl contacts in aspartates of high-resolution protein structures. J Phys Chem B 2010; 114:1038-49. [PMID: 20039723 DOI: 10.1021/jp909339r] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The folding and stability of a polypeptide chain are due to many different and simultaneous noncovalent interactions. Recent studies have observed several novel and counterintuitive contacts in protein structures, and the nature of interactions due to such contacts is yet to be fully elucidated. We have identified carbonyl-carbonyl intraresidue contacts in 102 Asp residues from a data set of high-resolution protein structures. At the outset, it appears that such close approach of two carbonyl oxygen atoms is energetically not favorable. We have carried out ab initio quantum chemical calculations on 10 representative examples of self-contacting Asp residues from different regions of the Ramachandran map. Potential energy scan using three levels of theory (HF, B3LYP, and MP2) and two basis sets (6-31+G* and 6-31++G**) was performed by varying the side-chain dihedral angle chi(1) while keeping all other parameters corresponding to that observed in the protein structures. We also calculated interaction energies by considering the surrounding interacting residues and water molecules. Our results show that the energy difference between a self-contacting Asp residue from the crystal structures and the minimum energy conformations is about 10-15 kcal/mol. This small energy difference is compensated by its interactions with the surrounding residues and water molecules as observed in the interaction energy analysis. The results are independent of the levels of theory used. The contacting carbonyl-carbonyl groups adopt a sheared parallel motif orientation which helps to expose both the backbone and side-chain carbonyl oxygen atoms and enable them to participate in tertiary interactions. Natural bond orbital calculations indicate that carbonyl-carbonyl groups in self-contacting Asp residues interact through n --> pi* electron delocalization. The geometry analysis and nature of chemical interactions together explain the rationale for the existence of such Asp residues in protein structures and their importance in the protein stability.
Collapse
Affiliation(s)
- Tuhin Kumar Pal
- Department of Biological Sciences and Bioengineering, Indian Institute of Technology, Kanpur, Kanpur 208016, India
| | | |
Collapse
|
4
|
Li X, Wang H, Bowen KH, Grégoire G, Lecomte F, Schermann JP, Desfrançois C. The parent anion of the RGD tripeptide: Photoelectron spectroscopy and quantum chemistry calculations. J Chem Phys 2009; 130:214301. [DOI: 10.1063/1.3137095] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
|
5
|
Echenique P, Alonso JL. Efficient model chemistries for peptides. I. General framework and a study of the heterolevel approximation in RHF and MP2 with Pople split-valence basis sets. J Comput Chem 2008; 29:1408-22. [DOI: 10.1002/jcc.20900] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
6
|
Grégoire G, Gaigeot MP, Marinica DC, Lemaire J, Schermann JP, Desfrançois C. Resonant infrared multiphoton dissociation spectroscopy of gas-phase protonated peptides. Experiments and Car–Parrinello dynamics at 300 K. Phys Chem Chem Phys 2007; 9:3082-97. [PMID: 17612732 DOI: 10.1039/b618094a] [Citation(s) in RCA: 89] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The gas-phase structures of protonated peptides are studied by means of resonant infrared multiphoton dissociation spectroscopy (R-IRMPD) performed with a free electron laser. The peptide structures and protonation sites are obtained through comparison between experimental IR spectra and their prediction from quantum chemistry calculations. Two different analyses are conducted. It is first supposed that only well-defined conformations, sufficiently populated according to a Boltzmann distribution, contribute to the observed spectra. On the contrary, DFT-based Car-Parrinello molecular dynamics simulations show that at 300 K protonated peptides no longer possess well-defined structures, but rather dynamically explore the set of conformations considered in the first conventional approach.
Collapse
Affiliation(s)
- G Grégoire
- Laboratoire de Physique des Lasers UMR 7538 Institut Galilée, Université Paris 13, Villetaneuse, France.
| | | | | | | | | | | |
Collapse
|
7
|
Sang-aroon W, Ruangpornvisuti V. Conformational study of cationic, zwitterionic, anionic species of aspartic acid, water-added forms and their protonation. A DFT method. ACTA ACUST UNITED AC 2006. [DOI: 10.1016/j.theochem.2005.10.036] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
8
|
Ramachandran backbone potential energy surfaces of aspartic acid and aspartate residues: implications on allosteric sites in receptor–ligand complexations. ACTA ACUST UNITED AC 2003. [DOI: 10.1016/j.theochem.2003.08.055] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
9
|
Koo JC, Lam JS, Chass GA, Torday LL, Varro A, Papp JG. Conformational dependence of the intrinsic acidity of the aspartic acid residue sidechain in N-acetyl-l-aspartic acid-N′-methylamide. ACTA ACUST UNITED AC 2003. [DOI: 10.1016/s0166-1280(02)00639-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
10
|
How reliable could economic Hartree–Fock computations be in studying large, folded peptides? A comparative HF and DFT case study on N- and C-protected aspartic acid. ACTA ACUST UNITED AC 2002. [DOI: 10.1016/s0166-1280(02)00579-1] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|