Ring Flip of Proline Residue via the Transition State with an Envelope Conformation

Puckering transition of the proline residue along the pseudorotational path: revisited

Puckering transitions of the proline residue for Ac-Pro-X with trans and cis prolyl peptide bonds were explored along the pseudorotation phase angle using DFT methods in the gas phase and in water.

Measuring the Energy Barrier of the Structural Change That Initiates Amyloid Formation

Obtaining kinetic and thermodynamic information for protein amyloid formation can yield new insight into the mechanistic details of this biomedically important process. The kinetics of the structural change that initiates the amyloid pathway, however, has been challenging to access for any amyloid protein system. Here, using the protein β-2-microglobulin (β2m) as a model, we measure the kinetics and energy barrier associated with an initial amyloidogenic structural change. Using covalent labeling and mass spectrometry, we measure the decrease in solvent accessibility of one of β2m's Trp residues, which is buried during the initial structural change, as a way to probe the kinetics of this structural change at different temperatures and under different amyloid forming conditions. Our results provide the first-ever measure of the activation barrier for a structural change that initiates the amyloid formation pathway. The results also yield new mechanistic insight into β2m's amyloidogenic structural change, especially the role of Pro32 isomerization in this reaction.

Which DFT levels of theory are appropriate in predicting the prolyl cis–trans isomerization in solution?

DFTs were assessed for the conformational preferences of the peptides containing Pro and its derivatives in chloroform and water.

Comparative effects of trifluoromethyl- and methyl-group substitutions in proline

What is the outcome of trifluoromethyl-/methyl-substitution in each position of the proline ring? Look inside to find out.

Exploring the mechanism of isomerisation and water-migration in the water-complexes of amino-acid l-proline: electrostatic potential and vibrational analysis

This work reveals interesting pathways for water-migration and neutral ↔ zwitterionic isomerisation in the water complexes of l-proline.

Mechanistic insights into Pin1 peptidyl-prolyl cis-trans isomerization from umbrella sampling simulations

The peptidyl-proyl isomerase Pin1 plays a key role in the regulation of phospho(p)-Ser/Thr-Pro proteins, acting as a molecular timer of the cell cycle. After recognition of these motifs, Pin1 catalyzes the rapid cis-trans isomerization of proline amide bonds of substrates, contributing to maintain the equilibrium between the two conformations. Although a great interest has arisen on this enzyme, its catalytic mechanism has long been debated. Here, the cis-trans isomerization of a model peptide system was investigated by means of umbrella sampling simulations in the Pin1-bound and unbound states. We obtained free energy barriers consistent with experimental data, and identified several enzymatic features directly linked to the acceleration of the prolyl bond isomerization. In particular, an enhanced autocatalysis, the stabilization of perturbed ground state conformations, and the substrate binding in a procatalytic conformation were found as main contributions to explain the lowering of the isomerization free energy barrier.

Conformational preferences and cis-trans isomerization of L-3,4-dehydroproline residue

The conformational study of N-acetyl-N'-methylamide of L-3,4-dehydroproline (Ac-Dhp-NHMe, the Dhp dipeptide) is carried out using hybrid density functional methods with the self-consistent reaction field method in the gas phase and in solution (chloroform and water). The incorporation of a double bond between C(beta) and C(gamma) into the prolyl ring results in the puckering, backbone population, and barriers to prolyl cis-trans isomerization different from those of the Pro dipeptide. For local minima of the Dhp dipeptide in the gas phase and in water, the C(beta)-C(gamma) bonds become shorter by approximately 0.2 A and the bond angles C(alpha)-C(beta)-C(gamma) and C(beta)-C(gamma)-C(delta) are widened by approximately 8 degrees than those of the Pro dipeptide, and the puckering amplitude is computed to be 0.01-0.07 A, indicating that the 3,4-dehydroprolyl ring is quite less puckered. However, polyproline-like conformations become preferred and the relative stability of the conformation tC with a C(7) intramolecular hydrogen bond decreases as the solvent polarity increases, as found for the Pro dipeptide. The barriers to cis-trans isomerization of the Ac-Dhp peptide bond increase with the increase of solvent polarity and the isomerization is likely to proceed through the clockwise rotation in water, as found for the prolyl peptide bond. The hydrogen bond between the prolyl nitrogen and the following amide N-H group seems to contribute in stabilizing the transition state structures.

Conformational analysis of L-prolines in water

The results of the ring conformational analysis of L-proline, N-acetyl-L-proline, and trans-4-hydroxy-L-proline by NMR combined with calculations using density functional theory (DFT) and molecular dynamics (MD) are reported. Accurate values of 1H-1H J-couplings in water and other solvents have been determined. Using a two-site equilibrium model, the Cgamma-endo conformer of L-proline in water has been identified as intermediate between gammaTdelta [twist(Cgamma-endo, Cdelta-exo)] and gammaE [envelope(Cgamma-endo)] and the Cgamma-exo conformer as betaTgamma. Both conformers were equally populated at room temperature. The N-acetyl [cis-rotamer gammaTbeta(80%)/gammaE(20%) and trans-rotamer gammaTbeta(61%)/gammaE(39%)] and 4-hydroxy (gammaEpsilon) derivatives showed significant changes in both the population and the geometries of the preferred ring conformers. The combination of NMR predicted populations with the DFT B3LYP/6-311+G(2d,p)/IEFPCM calculations proved successful, resulting in fairly accurate predictions of J-couplings. Simulations using MD were mostly in favor of the two-site equilibrium model between Cgamma-endo and Cgamma-exo conformers, similar to that used for the analysis of NMR J-couplings. Various force fields examined for MD simulations failed to reproduce the ring conformational geometries and populations of L-proline in water accurately, while significantly better agreement with NMR was found for trans-N-acetyl-L-proline using GROMOS and AMBER force fields.

Backbone conformational preferences and pseudorotational ring puckering of 1-aminocyclopentane-1-carboxylic acid

We have used quantum mechanical calculations at the B3LYP/6-311G(d,p) level to determine the conformational preferences of the N-acetyl-N'-methylamide derivative of 1-aminocyclopentane-1-carboxylic acid in the gas phase, chloroform solution, and water solution. The backbone conformation of this dipeptide has been described through the dihedral angles varphi and psi, while the pseudorotational phase angle was used to define the conformation of the cyclopentane ring. Results indicate that the backbone flexibility of this amino acid is restricted by the cyclic nature of the side chain, the relative stability of the different conformations depending on the polarity of the environment. The potential energy of the pseudorotation was also studied as a function of the backbone conformation. Interestingly, the conformation of the cyclic side chain depends on the backbone arrangement. Furthermore, the number of pseudorotational states accessible at room temperature is high in all the investigated environments, especially in aqueous solution. Finally, a set of force-field parameters for classical molecular mechanics calculations was developed for the investigated amino acid. Molecular dynamics simulations in both chloroform and aqueous solutions were performed to demonstrate the reliability of such parameters.

Puckering Transition of Proline Residue in Water

The puckering transition of the proline residue with trans and cis prolyl peptide bonds was explored by optimizations along the torsion angle chi1 of the prolyl ring using quantum-chemical methods in water. By analyzing the potential energy surfaces and local minima in water, it is observed that the puckering transition of the proline residue proceeds from a down-puckered conformation to an up-puckered one and vice versa through the transition state with an envelope form having the N atom at the top of the envelope and not a planar one, as seen in the gas phase, although the backbone conformations are different in the gas phase and in water. The barriers to the puckering transition DeltaGup-->down are estimated to be 3.12 and 3.00 kcal/mol for trans and cis conformers at the B3LYP/6-311++G(d,p) level of theory in water, respectively, which are about 1.7 kcal/mol higher than those in the gas phase. Out of 2197 prolines from the 241 high-resolution PDB chains, four transition-state-like structures with the envelope ring puckering are identified. Three of them have the trans prolyl peptide bonds and one has the cis one. The favorable or steric interactions by neighboring residues may be responsible for the stabilization of these transition-state-like ring structures in the proteins.

Conformational Preferences and cis−trans Isomerization of Azaproline Residue

The conformational study of N-acetyl-N'-methylamide of azaproline (Ac-azPro-NHMe, the azPro dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the effects of the replacement of the backbone CHalpha group by the nitrogen atom on the conformational preferences and prolyl cis-trans isomerization in the gas phase and in solution (chloroform and water). The incorporation of the Nalpha atom into the prolyl ring results in the different puckering, backbone population, and barriers to prolyl cis-trans isomerization from those of Ac-Pro-NHMe (the Pro dipeptide). In particular, the azPro dipeptide has a dominant backbone conformation D (beta2) with the cis peptide bond preceding the azPro residue in both the gas phase and solution. This may be ascribed to the favorable electrostatic interaction or intramolecular hydrogen bond between the prolyl nitrogen and the amide hydrogen following the azPro residue and to the absence of the unfavorable interactions between electron lone pairs of the acetyl carbonyl oxygen and the prolyl Nalpha. This calculated higher population of the cis peptide bond is consistent with the results from X-ray and NMR experiments. As the solvent polarity increases, the conformations B and B* with the trans peptide bond become more populated and the cis population decreases more, which is opposite to the results for the Pro dipeptide. The conformation B lies between conformations D and A (alpha) and conformation B* is a mirror image of the conformation B on the phi-psi map. The barriers to prolyl cis-trans isomerization for the azPro dipeptide increase with the increase of solvent polarity, and the cis-trans isomerization proceeds through only the clockwise rotation with omega' approximately +120 degrees about the prolyl peptide bond for the azPro dipeptide in the gas phase and in solution, as seen for the Pro dipeptide. The pertinent distance d(N...H-NNHMe) and the pyramidality of imide nitrogen can describe the role of this hydrogen bond in stabilizing the transition state structure and the lower rotational barriers for the azPro dipeptide than those for the Pro dipeptide in the gas phase and in solution.

Conformational Preferences of Proline Analogues with Different Ring Size

The conformational study on L-azetidine-2-carboxylic acid (Ac-Aze-NHMe, the Aze dipeptide) and (S)-piperidine-2-carboxylic acid (Ac-Pip-NHMe, the Pip dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the differences in conformational preferences and cis-trans isomerization for proline residue and its analogues with different ring size in the gas phase and in solution (chloroform and water). The change of ring size by deleting a CH2 group from or adding a CH2 group to the prolyl ring results the remarkable changes in backbone and ring structures compared with those of the Pro dipeptide, especially in the C'-N imide bond length and the bond angles around the N-C(alpha) bond. The four-membered azetidine ring can have either puckered structure depending on the backbone structure because of the less puckered structure. The six-membered piperidine ring can adopt chair and boat conformations, but the chair conformation is more preferred than the boat conformation. These calculated preferences for puckering are consistent with experimental results from analysis of X-ray structures of Aze- and Pip-containing peptides. On going from Pro to Aze to Pip, the axiality (i.e., a tendency to adopt the axial orientation) of the NHMe group becomes stronger, which can be ascribed to reduce the steric hindrances between 1,2-substituted Ac and NHMe groups. As the solvent polarity increases, the polyproline II-like conformation becomes more populated and the relative stability of conformation tC with a C7 hydrogen bond between C'=O of the amino group and N-H of the carboxyl group decreases for both the Aze and Pip dipeptides, as seen for the Pro dipeptide. The cis population and rotational barriers for the imide bond increase with the increase of solvent polarity for both the Aze and Pip dipeptides, as seen for the Pro dipeptide. In particular, the cis-trans isomerization proceeds in common through only the clockwise rotation with omega' approximately +120 degrees about azetyl and piperidyl peptide bonds in the gas phase and in solution, as seen for alanyl and prolyl peptide bonds. The pertinent distance d(N...H-N(NHMe)) and the pyramidality of imide nitrogen can describe the role of this hydrogen bond in stabilizing the transition state structure, but the lower rotational barriers for the Aze and Pip dipeptides than those for the Pro dipeptide, which is observed from experiments, cannot be rationalized.

Conformational study of palindromic tripeptides (GPG, IPI and KPK) in HIV-1 protease--a density functional theory study

A comparative study has been carried out on three palindromic tripeptides Gly-Pro-Gly, Ile-Pro-Ile and Lys-Pro-Lys which were present in HIV protein along with their analogues applying density functional computation at B3LYP/6-31G* level of theory. Discrepancy from the structural analysis has been noted for all the systems and it was found to be more for amide capped structure at the C terminal of proline. The puckering amplitude A and Phase angle P of the pyrrolidine ring of proline in the chosen palindromic tripeptides and their analogues were calculated from the endocyclic torsion angles. The minimum energy conformers lying well within the prescribed region of proline were obtained for the derived compounds from potential energy surface scan mentioning that no role has been played by its terminal residues. This is further supported by the simulated amide bands identifying the helical structure for all three palindromic tripeptides signifying the importance of proline. The molecular properties such as stabilization energy, chemical hardness along with dipole moment were calculated and interpreted. The values of Calpha-H(s) and the peptide backbone N-Calpha-CO for all the selected conformers specify the three palindromic tripeptides to have a symmetrical achiral structure.

Conformational Preferences of Non-Prolyl and Prolyl Residues

The conformational study on Ac-Ala-NHMe (the alanine dipeptide) and Ac-Pro-NHMe (the proline dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the differences in the backbone conformational preference and the cis-trans isomerization for the non-prolyl and prolyl residues in the gas phase and in the solutions (chloroform and water). For the alanine and proline dipeptides, with the increase of solvent polarity, the populations of the conformation tC with an intramolecular C(7) hydrogen bond significantly decrease, and those of the polyproline II-like conformation tF and the alpha-helical conformation tA increase, which is in good agreement with the results from circular dichroism and NMR experiments. For both the dipeptides, as the solvent polarity increases, the relative free energy of the cis conformer to the trans conformer decreases and the rotational barrier to the cis-trans isomerization increases. It is found that the cis-trans isomerization proceeds in common through only the clockwise rotation with omega' approximately +120 degrees about the non-prolyl and prolyl peptide bonds in both the gas phase and the solutions. The pertinent distance d(N...H-N(NHMe)) can successfully describe the increase in the rotational barriers for the non-prolyl and prolyl trans-cis isomerization as the solvent polarity increases and the higher barriers for the non-prolyl residue than for the prolyl residue, as seen in experimental and calculated results. By analysis of the contributions to rotational barriers, the cis-trans isomerization for the non-prolyl and prolyl peptide bonds is proven to be entirely enthalpy driven in the gas phase and in the solutions. The calculated cis populations and rotational barriers to the cis-trans isomerization for both the dipeptides in chloroform and/or water accord with the experimental values.

Conformational Preferences of Proline Oligopeptides

A conformational study on the terminally blocked proline oligopeptides, Ac-(Pro)(n)()-NMe(2) (n = 2-5), is carried out using the ab initio Hartree-Fock level of theory with the self-consistent reaction field method in the gas phase and in solutions (chloroform, 1-propanol, and water) to explore the preference and transition between polyproline II (PPII) and polyproline I (PPI) conformations depending on the chain length, the puckering, and the solvent. The mean differences in the free energy per proline of the up-puckered conformations relative to the down-puckered conformations for both diproline and triproline increases for the PPII-like conformations and decreases for the PPI-like conformations as the solvent polarity increases. These calculated results indicate that the PPII-like structures have preferentially all-down puckerings in solutions, whereas the PPI-like structures have partially mixed puckerings. The free energy difference per proline residue between the PPII- and PPI-like structures decreases as the proline chain becomes longer in the gas phase but increases as the proline chain becomes longer in solutions and the solvent polarity increases. In particular, our calculated results indicate that each of the proline oligopeptides can exist as an ensemble of conformations with the trans and cis peptide bonds in solutions, although the PPII-like structure with all-trans peptide bonds is dominantly preferred, which is reasonably consistent with the previously observed results. In diproline Ac-(Pro)(2)-NMe(2), the rotational barrier to the cis-to-trans isomerization for the first prolyl peptide bond increases as the solvent polarity increases, whereas the rotational barrier for the second prolyl peptide bond does not show the monotonic increase as the solvent polarity increases. When the rotational barriers for these two prolyl peptide bonds were compared, it could be deduced that the conformational transition from PPI with the cis peptide bond to PPII with the trans peptide bond is initiated at the C-terminus and proceeds to the N-terminus in water. This is consistent with the results from NMR experiments on polyproline in D(2)O but opposite to the results from enzymatic hydrolysis kinetics experiments on polyproline.

Conformational Preference and Cis−Trans Isomerization of 4(R)-Substituted Proline Residues

We report here the conformational preference and prolyl cis-trans isomerization of 4(R)-substituted proline dipeptides, N-acetyl-N'-methylamides of 4(R)-hydroxy-L-proline and 4(R)-fluoro-L-proline (Ac-Hyp-NHMe and Ac-Flp-NHMe, respectively), studied at the HF/6-31+G(d), B3LYP/6-31+G(d), and B3LYP/6-311++G(d,p) levels of theory. The 4(R)-substitution by electron-withdrawing groups did not result in significant changes in backbone torsion angles as well as endocyclic torsion angles of the prolyl ring. However, the small changes in backbone torsion angles phi and psi and the decrease of bond lengths r(Cbeta-Cgamma) or r(Cgamma-Cdelta) appear to induce the increase of the relative stability of the trans up-puckered conformation and to alter the relative stabilities of transition states for prolyl cis-trans isomerization. Solvation free energies of local minima and transition states in chloroform and water were calculated using the conductor-like polarizable continuum model at the HF/6-31+G(d) level of theory. The population of trans up-puckered conformations increases in the order Ac-Pro-NHMe < Ac-Hyp-NHMe < Ac-Flp-NHMe in chloroform and water. The increase in population for trans up-puckered conformations in solution is attributed to the increase in population for the polyproline-II-like conformations with up puckering. The barriers DeltaGct++ to prolyl cis-to-trans isomerization for Ac-Hyp-NHMe and Ac-Flp-NHMe increase as the solvent polarity increases, as seen for Ac-Pro-NHMe. In particular, it was identified that the cis-trans isomerization proceeds through the clockwise rotation about the prolyl peptide bond for Ac-Hyp-NHMe and Ac-Flp-NHMe in chloroform and water, as seen for Ac-Pro-NHMe.

Antibiotic binding to monozinc CphA beta-lactamase from Aeromonas hydropila: quantum mechanical/molecular mechanical and density functional theory studies

The active-site dynamics of apo CphA beta-lactamase from Aeromonas hydropila and its complex with a beta-lactam antibiotic molecule (biapenem) are simulated using a quantum mechanical/molecular mechanical (QM/MM) method and density functional theory (DFT). The quantum region in the QM/MM simulations, which includes the Zn(II) ion and its ligands, the antibiotic molecule, the catalytic water, and an active-site histidine residue, was treated using the self-consistent charge density functional tight binding (SCC-DFTB) model. Biapenem is docked at the active site unambiguously, based on a recent X-ray structure of an enzyme-intermediate complex. The substrate is found to form the fourth ligand of the zinc ion with its 3-carboxylate oxygen and to hydrogen bond with several active-site residues. The stability of the metal-ligand bonds and the hydrogen-bond network is confirmed by 500 ps molecular dynamics simulations of both the apo enzyme and the substrate-enzyme complex. The structure and dynamics of the substrate-enzyme complex provide valuable insights into the mode of catalysis in such enzymes that is central to the bacterial resistance to beta-lactam antibiotics.

Puckering Transition of 4-Substituted Proline Residues

The puckering transition of 4-substituted proline residues by electron-withdrawing groups, i.e., 4(R)-hydroxy-L-proline (Hyp) and 4(R)-fluoro-L-proline (Flp) residues, with trans and cis prolyl peptide bonds was studied by adiabatic optimizations along the torsion angle chi1 of the prolyl ring at the HF/6-31+G(d) level. By analyzing the potential energy surface and local minima, it is observed that the puckering transition of the prolyl ring for Hyp and Flp residues proceeds from a down-puckered conformation to an up-puckered one through the transition state with an envelope form having the N atom at the top of envelope and not a planar one for both trans and cis conformers, which is the same as found for the unsubstituted proline residue. At HF/6-31+G(d) and B3LYP/6-311++G(d,p) levels, the structures of the backbone and prolyl ring for local minima of Ac-Hyp-NHMe and Ac-Flp-NHMe are quite similar to those of Ac-Pro-NHMe. However, the relative stability of the up-puckered conformation to the down-puckered one is increased for Ac-Hyp-NHMe with the cis imide bond and for Ac-Flp-NHMe with the trans and cis imide bonds. In particular, the 4(R)-substitution by hydroxy and fluorine groups has brought some structural changes in the prolyl ring of the transition states and the changes in barriers for the puckering transition. The puckering transitions for Ac-Hyp-NHMe and Ac-Flp-NHMe are proven to be predominantly electronically driven by analyzing the electronic and enthalpic contributions to the barriers, as seen for Ac-Pro-NHMe.

First Principle Computational Study on the Full Conformational Space of l-Proline Diamides

Ab initio molecular orbital computations were carried out at three levels of theory: RHF/3-21G, RHF/6-31G(d), and B3LYP/6-31G(d), on four model systems of the amino acid proline, HCO-Pro-NH2 [I], HCO-Pro-NH-Me [II], MeCO-Pro-NH2 [III], and MeCO-Pro-NH-Me [IV], representing a systematic variation in the protecting N- and C-terminal groups. Three previously located backbone conformations, gammaL, epsilonL, and alphaL, were characterized together with two ring-puckered forms syn (gauche+ = g+) or "DOWN" and anti (gauche- = g-) or "UP", as well as trans-trans, trans-cis, cis-trans, and cis-cis peptide bond isomers. The topologies of the conformational potential energy cross-sections (PECS) of the potential energy hypersurfaces (PEHS) for compounds [I]-[IV] were explored and analyzed in terms of potential energy curves (PEC), and HCO-Pro-NH2 [I] was also analyzed in terms of potential energy surfaces (PESs). Thermodynamic functions were also calculated for HCO-Pro-NH2 [I] at the CBS-4M and G3MP2 levels of theory. The study confirms that the use of the simplest model, compound [I] with P(N) = P(C) = H, along with the RHF/3-21G level of theory, is an acceptable practice for the analysis of peptide models because only minor differences in geometry and stability are observed.

Cis–trans isomerization and puckering of proline residue

We report here the results on N-acetyl-L-proline-N'-methylamide (Ac-Pro-NHMe) calculated at the HF/6-31+G(d) level with the conductor-like polarizable continuum model (CPCM) of self-consistent reaction field methods to investigate the changes of backbone and prolyl ring along the cis-trans isomerization of the prolyl peptide bond. From the potential energy surface, the barrier to ring flip from the down-puckered conformation to the up-puckered one is estimated to be 2.5 and 3.2 kcal/mol for trans and cis conformers of Ac-Pro-NHMe, respectively. In particular, the ring flip seems to be inaccessible in the intermediate regions between trans and cis conformations, because of higher barriers (approximately 13-19 kcal/mol) to rotation of the prolyl peptide bond. The torsion angles for backbone and prolyl ring vary largely around the transition states at omega' approximately 120 degrees and -70 degrees for the prolyl peptide bond. Three kinds of puckering amplitudes show the same trend of puckering along the cis-trans isomerization although their absolute values are different. In particular, trans and cis conformations have the almost same degree of puckering. The cis populations and barriers to rotation of the prolyl peptide bond for Ac-Pro-NHMe are increased with the increase of solvent polarity, which is mainly ascribed to the decreases of relative free energies for cis conformations and the increase of relative free energies for transition states.