Nature of Glycine and Its α-Carbon Radical in Aqueous Solution: A Theoretical Investigation

Modeling the Initiation Phase of the Catalytic Cycle in the Glycyl-Radical Enzyme Benzylsuccinate Synthase

The reaction of benzylsuccinate synthase, the radical-based addition of toluene to a fumarate cosubstrate, is initiated by hydrogen transfer from a conserved cysteine to the nearby glycyl radical in the active center of the enzyme. In this study, we analyze this step by comprehensive computer modeling, predicting (i) the influence of bound substrates or products, (ii) the energy profiles of forward- and backward hydrogen-transfer reactions, (iii) their kinetic constants and potential mechanisms, (iv) enantiospecificity differences, and (v) kinetic isotope effects. Moreover, we support several of the computational predictions experimentally, providing evidence for the predicted H/D-exchange reactions into the product and at the glycyl radical site. Our data indicate that the hydrogen transfer reactions between the active site glycyl and cysteine are principally reversible, but their rates differ strongly depending on their stereochemical orientation, transfer of protium or deuterium, and the presence or absence of substrates or products in the active site. This is particularly evident for the isotope exchange of the remaining protium atom of the glycyl radical to deuterium, which appears dependent on substrate or product binding, explaining why the exchange is observed in some, but not all, glycyl-radical enzymes.

Radical Stabilization Energies for Enzyme Engineering: Tackling the Substrate Scope of the Radical Enzyme QueE

Experimental assessment of catalytic reaction mechanisms and profiles of radical enzymes can be severely challenging due to the reactive nature of the intermediates and sensitivity of cofactors such as iron-sulfur clusters. Here, we present an enzyme-directed computational methodology for the assessment of thermodynamic reaction profiles and screening for radical stabilization energies (RSEs) for the assessment of catalytic turnovers in radical enzymes. We have applied this new screening method to the radical S-adenosylmethione enzyme 7-carboxy-7-deazaguanine synthase (QueE), following a detailed molecular dynamics (MD) analysis that clarifies the role of both specific enzyme residues and bound Mg²⁺, Ca²⁺, or Na⁺. The MD simulations provided the basis for a statistical approach to sample different conformational outcomes. RSE calculation at the M06-2X/6-31+G* level of theory provided the most computationally cost-effective assessment of enzyme-based energies, facilitated by an initial triage using semiempirical methods. The impact of intermolecular interactions on RSE was clearly established, and application to the assessment of potential alternative substrates (focusing on radical clock type rearrangements) proposes a selection of carbon-substituted analogues that would react to afford cyclopropylcarbinyl radical intermediates as candidates for catalytic turnover by QueE.

Toward theoretical terahertz spectroscopy of glassy aqueous solutions: partially frozen solute-solvent couplings of glycine in water

The molecular-level understanding of THz spectra of aqueous solutions under ambient conditions has been greatly advanced in recent years. Here, we go beyond previous analyses by performing ab initio molecular dynamics simulations of glycine in water with artificially frozen solute or solvent molecules, respectively, while computing the total THz response as well as its decomposition into mode-specific resonances based on the "supermolecular solvation complex" technique. Clamping the water molecules and keeping glycine moving breaks the coupling of glycine to the structural dynamics of the solvent, however, the polarization and dielectric solvation effects in the static solvation cage are still at work since the full electronic structure of the quenched solvent is taken into account. The complementary approach of fixing glycine reveals both the dynamical and electronic response of the solvation cage at the level of its THz response. Moreover, to quantitatively account for the electronic contribution solely due to solvent embedding, the solute species is "vertically desolvated", thus preserving the fully coupled solute-solvent motion in terms of the solute's structural dynamics in solution, while its electronic structure is no longer subject to solute-solvent polarization and charge transfer effects. When referenced to the free simulation of Gly(aq), this three-fold approach allows us to decompose the THz spectral contributions due to the correlated solute-solvent dynamics into entirely structural and purely electronic effects. Beyond providing hitherto unknown insights, the observed systematic changes of THz spectra in terms of peak shifts and lineshape modulations due to conformational freezing and frozen solvation cages might be useful to investigate the solvation of molecules in highly viscous H-bonding solvents such as ionic liquids and even in cryogenic ices as relevant to polar stratospheric and dark interstellar clouds.

Solvation of the Glycyl Radical

The effect of adding explicit water molecules to the neutral (N) and zwitterionic (Z) forms of the glycyl radical has been examined. The results show that a minimum of three water molecules is required to stabilize the Z radical as a local minimum, with an energy gap of 123 kJ mol^-1 between the N and Z forms at this point, in favor of the N form. Increasing the number of water molecules to ∼20 leads to a converged Z-N energy difference of ∼50 kJ mol^-1 still in favor of the N form, even though the radical is not considered fully solvated from a structural point of view. Thus, energetic convergence is determined mainly by solvation of the polar functional groups, and a complete coverage of the entire molecule is not necessary. Because aqueous closed-shell glycine exists as a zwitterion while aqueous glycyl radical prefers the neutral form, the conversion between the two necessitates a change along the hydrogen-abstraction reaction pathway. In this regard, the transition structure for α-hydrogen abstraction by the ·OH radical has greater resemblance to glycine than to the glycyl radical. Overall, the barrier for hydrogen abstraction from Z glycine is larger than that from the N isomer, and this might act to provide some protection against radical damage to the free amino acid in the (aqueous) biological environment.

Radical Reaction Control in the AdoMet Radical Enzyme CDG Synthase (QueE): Consolidate, Destabilize, Accelerate

Controlling radical intermediates and thus catalysing and directing complex radical reactions is a central feature of S-adensosylmethionine (SAM)-dependent radical enzymes. We report ab initio and DFT calculations highlighting the specific influence of ion complexation, including Mg2+ , identified as a key catalytic component on radical stability and reaction control in 7-carboxy-7-deazaguanine synthase (QueE). Radical stabilisation energies (RSEs) of key intermediates and radical clock-like model systems of the enzyme-catalysed rearrangement of 6-carboxytetrahydropterin (CPH4), reveals a directing role of Mg2+ in destabilising both the substrate-derived radical and corresponding side reactions, with the effect that the experimentally-observed rearrangement becomes dominant over possible alternatives. Importantly, this is achieved with minimal disruption of the thermodynamics of the substrate itself, affording a novel mechanism for an enzyme to both maintain binding potential and accelerate the rearrangement step. Other mono and divalent ions were probed with only dicationic species achieving the necessary radical conformation to facilitate the reaction.

Site specificity on OH α-H abstraction reaction for a zwitterionic β-hairpin peptide in aqueous solution: a theoretical investigation

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.