Site Specificity on OH α-H Abstraction Reaction for a Zwitterionic β-Hairpin Peptide in Aqueous Solution: A Theoretical Investigation

Abstract

Protein backbone oxidation was investigated by studying the α-H abstraction reaction in a ß-hairpin peptide, called Chignolin (PDB ID 1UAO), with density functional theory calculation at B3LYP/6-31G(d,p) without any constraint. In order to stabilize the zwitterionic form of Chignolin with the salt bridges, the effects of aqueous solution were implemented by using microsolvation combined with a conductor-like polarizable continuum model (CPCM). Comparison between three glycine residues located at three different sites in Chignolin was used to examine the possible site specificity of this backbone oxidation. To construct the reaction profile of these α-H abstraction reactions, the pre- and postreactive complexes along with their associated transition states were located and verified with the intrinsic reaction coordinate (IRC) method. The bond dissociation energy and reaction rates of these OH α-H abstraction reactions were calculated with transition state theory. The differences in this abstraction reaction between the neutral and zwitterionic forms of Chignolin were also compared. A molecular dynamics simulation was implemented to study the explicit solvation effect on the abstracted Chignolin structure. The range of the simulation time scale covers from femtoseconds to microseconds, i.e., from onset of the abstraction to the abstracted products reaching thermal equilibrium. Our results show that there are three kinds of site-specificity in this abstraction reaction. The reactivity and stability of the abstraction products and their abstraction modes are all dependent on the location where OH attacks. Furthermore, the free energy landscapes of these abstraction products are distinctively different. This may imply that the pathological disorders or diseases caused by this type of radicals are also dependent on the abstraction location