Quantum Mechanics/molecular Mechanics Research Articles

Study of Interactions of Nucleoside Anticancer Drugs, Capecitabine and Gemcitabine, with SWNT and BNNT using Molecular and Quantum Mechanical Calculations

Abstract: Using single-wall carbon nanotubes (SWCNTs) and boron nitride nanotubes (BNNTs), this study evaluated the interactions between Capecitabine (CAP) and Gemcitabine (GEM). Molecular mechanics and quantum mechanics were used in the analysis. : The interaction between CP and GM with SWCNTs and BNNTs under various solvents was analyzed using the self-consistent reaction field (SCRF) as a model and DFT as the analytical method. Additionally, the effect of temperature on the stability of the molecule interactions during the experiment was examined. The thermodynamic properties of the title compounds were analyzed based on theoretical calculations. It included the calculation of Frontier Molecular Orbitals (FMOs) and Total Density of States (DOS). In addition to studying the ionization potential (I), we also examined the other molecular properties of the structures, such as electrophilicity (ω), electron affinity (A), chemical hardness (η), and electronic chemical potential (μ). To investigate the interaction between CAP and GEM with SWCNTs and BNNTs, molecular mechanics methods (MM) including AMBER, OPLS, CHARMM, and MM+ force fields, were employed. Monte Carlo simulation techniques were used to calculate the results at different temperatures. : The effects of the liquid phase and mixed solvent media with varying dielectric constants (Water, DMSO, Methanol, and Ethanol) on the interaction between CAP and GEM using the force fields described above were examined in this study. The correlation between data generated by Monte Carlo, Quantum Mechanics, and Molecular Mechanics was demonstrated. There was a striking similarity in the thermodynamic properties and conformer populations of all three materials.

QM/MM Study of the Catalytic Mechanism and Substrate Specificity of the Aromatic Substrate C-Methyltransferase Fur6.

In the field of medical chemistry and other organic chemistry, introducing a methyl group into a designed position has been difficult to achieve. However, owing to the vigorous developments in the field of enzymology, methyltransferases are considered potential tools for addressing this problem. Within the methyltransferase family, Fur6 catalyzes the methylation of C3 of 1,2,4,5,7-pentahydroxynaphthalene (PHN) using S-adenosyl-l-methionine (SAM) as the methyl donor. Here, we report the catalytic mechanism and substrate specificity of Fur6 based on computational studies. Our molecular dynamics (MD) simulation studies reveal the reactive form of PHN and its interactions with the enzyme. Our hybrid quantum mechanics/molecular mechanics (QM/MM) calculations suggest the reaction pathway of the methyl transfer step in which the energy barrier is 8.6 kcal mol-1. Our free-energy calculations with a polarizable continuum model (PCM) indicate that the final deprotonation step of the methylated intermediate occurs after it is ejected into the water solvent from the active center pocket of Fur6. Additionally, our studies on the protonation states, the highest occupied molecular orbital (HOMOs), and the energy barriers of the methylation reaction for the analogs of PHN demonstrate the mechanism of the specificity to PHN. Our study provides valuable insights into Fur6 chemistry, contributing to a deeper understanding of molecular mechanisms and offering an opportunity to engineer the enzyme to achieve high yields of the desired product(s).

Deciphering the structure of deep eutectic solvents: A computational study from the solute's viewpoint

In the last few years, Deep Eutectic Solvents (DESs) have emerged as sustainable alternatives to traditional organic solvents. The high degree of freedom when tailoring DES structures represents, at the same time, their most praised and most complicated feature. Indeed, given the enormous number of possible combinations, the selection of the most suitable DES for a given application cannot be based on a trial-and-error approach; therefore, reliable computational tools are needed to fully exploit DESs potentialities. In this work, we propose the first computational protocol to investigate absorption properties of solutes dissolved in DES. The protocol combines an accurate sampling of the solute-solvent phase-space by means of Molecular Dynamics (MD) and a fully atomistic Quantum Mechanics/Molecular Mechanics (QM/MM) approach to describe intermolecular interactions and spectral properties. In this way, specific interactions such as hydrogen bonding, which characterize DES complex interaction patterns, are properly modelled. The robustness and reliability of the method are proved by comparing computed data with experimental spectra of 2-hydroxymethylfurfural and syringic, vanillic, p-coumaric, gallic, and caffeic acids dissolved in DES and water. These molecules have been chosen because of their relevance in the frame of biomass valorization. Remarkably, the protocol allows getting insights into DES/water mixtures and opens up to the prediction of aqueous DES behaviour of the water threshold which causes the switch from water-in-DES to aqueous electrolyte-like environments.

Surface-Accelerated String Method for Locating Minimum Free Energy Paths.

We present a surface-accelerated string method (SASM) to efficiently optimize low-dimensional reaction pathways from the sampling performed with expensive quantum mechanical/molecular mechanical (QM/MM) Hamiltonians. The SASM accelerates the convergence of the path using the aggregate sampling obtained from the current and previous string iterations, whereas approaches like the string method in collective variables (SMCV) or the modified string method in collective variables (MSMCV) update the path only from the sampling obtained from the current iteration. Furthermore, the SASM decouples the number of images used to perform sampling from the number of synthetic images used to represent the path. The path is optimized on the current best estimate of the free energy surface obtained from all available sampling, and the proposed set of new simulations is not restricted to being located along the optimized path. Instead, the umbrella potential placement is chosen to extend the range of the free energy surface and improve the quality of the free energy estimates near the path. In this manner, the SASM is shown to improve the exploration for a minimum free energy pathway in regions where the free energy surface is relatively flat. Furthermore, it improves the quality of the free energy profile when the string is discretized with too few images. We compare the SASM, SMCV, and MSMCV using 3 QM/MM applications: a ribozyme methyltransferase reaction using 2 reaction coordinates, the 2'-O-transphosphorylation reaction of Hammerhead ribozyme using 3 reaction coordinates, and a tautomeric reaction in B-DNA using 5 reaction coordinates. We show that SASM converges the paths using roughly 3 times less sampling than the SMCV and MSMCV methods. All three algorithms have been implemented in the FE-ToolKit package made freely available.

Reconciling TD-DFT and CASPT2 electronic structure methods for describing the photophysics of DNA.

Time-dependent density functional theory (TD-DFT) and multiconfigurational second-order perturbation theory (CASPT2) are two of the most widely used methods to investigate photoinduced dynamics in DNA-based systems. These methods sometimes give diverse dynamics in physiological environments usually modeled by quantum mechanics/molecular mechanics (QM/MM) protocol. In this work, we demonstrate for the uridine test case that the underlying topology of the potential energy surfaces of electronic states involved in photoinduced relaxation is similar in both electronic structure methods. This is verified by analyzing surface-hopping dynamics performed at the QM/MM level on aqueous solvated uridine at TD-DFT and CASPT2 levels. By constraining the dynamics to remain on state we observe similar fluctuations in energy and relaxation lifetimes in surface-hopping dynamics in both TD-DFT and experimentally validated CASPT2 methods. This finding calls for a systematic comparison of the ES potential energy surfaces of DNA and RNA nucleosides at the single- and multi-reference levels of theory. The anomalous long excited state lifetime at the TD-DFT level is explained by trapping due to the tendency of TD-DFT in QM/MM schemes with electrostatic embedding to underestimate the energy of the state leading to a wrong energetic order. A study of the FC energetics suggests that improving the description of the surrounding environment through polarizable embedding or by the expansion of QM layer with hydrogen-bonded waters helps restore the correct state order at TD-DFT level. Thus by combining TDDFT with an accurate modeling of the environment, TD-DFT is positioned as the standout protocol to model photoinduced dynamics in DNA-based aggregates and multimers.

Catalytic divergencies in the mechanism of L-arginine hydroxylating nonheme iron enzymes.

Many enzymes in nature utilize a free arginine (L-Arg) amino acid to initiate the biosynthesis of natural products. Examples include nitric oxide synthases, which generate NO from L-Arg for blood pressure control, and various arginine hydroxylases involved in antibiotic biosynthesis. Among the groups of arginine hydroxylases, several enzymes utilize a nonheme iron(II) active site and let L-Arg react with dioxygen and α-ketoglutarate to perform either C3-hydroxylation, C4-hydroxylation, C5-hydroxylation, or C4-C5-desaturation. How these seemingly similar enzymes can react with high specificity and selectivity to form different products remains unknown. Over the past few years, our groups have investigated the mechanisms of L-Arg-activating nonheme iron dioxygenases, including the viomycin biosynthesis enzyme VioC, the naphthyridinomycin biosynthesis enzyme NapI, and the streptothricin biosynthesis enzyme OrfP, using computational approaches and applied molecular dynamics, quantum mechanics on cluster models, and quantum mechanics/molecular mechanics (QM/MM) approaches. These studies not only highlight the differences in substrate and oxidant binding and positioning but also emphasize on electronic and electrostatic differences in the substrate-binding pockets of the enzymes. In particular, due to charge differences in the active site structures, there are changes in the local electric field and electric dipole moment orientations that either strengthen or weaken specific substrate C-H bonds. The local field effects, therefore, influence and guide reaction selectivity and specificity and give the enzymes their unique reactivity patterns. Computational work using either QM/MM or density functional theory (DFT) on cluster models can provide valuable insights into catalytic reaction mechanisms and produce accurate and reliable data that can be used to engineer proteins and synthetic catalysts to perform novel reaction pathways.

Probing Nonadiabaticity of Proton-Coupled Electron Transfer in Ribonucleotide Reductase.

The enzyme ribonucleotide reductase, which is essential for DNA synthesis, initiates the conversion of ribonucleotides to deoxyribonucleotides via radical transfer over a 32 Å pathway composed of proton-coupled electron transfer (PCET) reactions. Previously, the first three PCET reactions in the α subunit were investigated with hybrid quantum mechanical/molecular mechanical (QM/MM) free energy simulations. Herein, the fourth PCET reaction in this subunit between C439 and guanosine diphosphate (GDP) is simulated and found to be slightly exoergic with a relatively high free energy barrier. To further elucidate the mechanisms of all four PCET reactions, we analyzed the vibronic and electron-proton nonadiabaticities. This analysis suggests that interfacial PCET between Y356 and Y731 is vibronically and electronically nonadiabatic, whereas PCET between Y731 and Y730 and between C439 and GDP is fully adiabatic and PCET between Y730 and C439 is in the intermediate regime. These insights provide guidance for selecting suitable rate constant expressions for these PCET reactions.

Molecular mechanisms and evolutionary robustness of a color switch in proteorhodopsins.

Proteorhodopsins are widely distributed photoreceptors from marine bacteria. Their discovery revealed a high degree of evolutionary adaptation to ambient light, resulting in blue- and green-absorbing variants that correlate with a conserved glutamine/leucine at position 105. On the basis of an integrated approach combining sensitivity-enhanced solid-state nuclear magnetic resonance (ssNMR) spectroscopy and linear-scaling quantum mechanics/molecular mechanics (QM/MM) methods, this single residue is shown to be responsible for a variety of synergistically coupled structural and electrostatic changes along the retinal polyene chain, ionone ring, and within the binding pocket. They collectively explain the observed color shift. Furthermore, analysis of the differences in chemical shift between nuclei within the same residues in green and blue proteorhodopsins also reveals a correlation with the respective degree of conservation. Our data show that the highly conserved color change mainly affects other highly conserved residues, illustrating a high degree of robustness of the color phenotype to sequence variation.

Iron Heme Enzyme-Catalyzed Cyclopropanations with Diazirines as Carbene Precursors: Computational Explorations of Diazirine Activation and Cyclopropanation Mechanism.

The mechanism of cyclopropanations with diazirines as air-stable and user-friendly alternatives to commonly employed diazo compounds within iron heme enzyme-catalyzed carbene transfer reactions has been studied by means of density functional theory (DFT) calculations of model systems, quantum mechanics/molecular mechanics (QM/MM) calculations, and molecular dynamics (MD) simulations of the iron carbene and the cyclopropanation transition state in the enzyme active site. The reaction is initiated by a direct diazirine-diazo isomerization occurring in the active site of the enzyme. In contrast, an isomerization mechanism proceeding via the formation of a free carbene intermediate in lieu of a direct, one-step isomerization process was observed for model systems. Subsequent reaction with benzyl acrylate takes place through stepwise C-C bond formation via a diradical intermediate, delivering the cyclopropane product. The origin of the observed diastereo- and enantioselectivity in the enzyme was investigated through MD simulations, which indicate a preferred formation of the cis-cyclopropane by steric control.

Discovery and Mechanistic Understanding of a Lipase from Rhizorhabdus dicambivorans for Efficient Ester Aminolysis in Aromatic Amines.

The formation of amide bonds via aminolysis of esters by lipases generates a diverse range of amide frameworks in biosynthetic chemistry. Few lipases have satisfactory activity towards bulky aromatic amines despite numerous attempts to improve the efficiency of this transformation. Here, we report the discovery of a new intracellular lipase (Ndbn) with a broad substrate scope. Ndbn turns over a range of esters and aromatic amines in the presence of water (2 %; v/v), producing a high yield of multiple valuable amides. Remarkably, a higher conversion rate was observed for the synthesis of amides from substrates with aromatic amine rather than aliphatic amines. Molecular dynamics (MD) and quantum mechanical/molecular mechanical (QM/MM) studies showcase the mechanism for the preference for aromatic amines, including a more suitable orientation, shorter catalytic distances in the active site pocket and a lower reaction barrier for aromatic than for aliphatic amines. This unique lipase is thus a promising biocatalyst for the efficient synthesis of aromatic amides.

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV.

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

How Aqueous Solvation Impacts the Frequencies and Intensities of Infrared Absorption Bands in Flavin: The Quest for a Suitable Solvent Model.

FTIR spectroscopy accompanied by quantum chemical simulations can reveal important information about molecular structure and intermolecular interactions in the condensed phase. Simulations typically account for the solvent either through cluster quantum mechanical (QM) models, polarizable continuum models (PCM), or hybrid quantum mechanical/molecular mechanical (QM/MM) models. Recently, we studied the effect of aqueous solvent interactions on the vibrational frequencies of lumiflavin, a minimal flavin model, using cluster QM and PCM models. Those models successfully reproduced the relative frequencies of four prominent stretching modes of flavin's isoalloxazine ring in the diagnostic 1450-1750 cm-1 range but poorly reproduced the relative band intensities. Here, we extend our studies on this system and account for solvation through a series of increasingly sophisticated models. Only by combining elements of QM clusters, QM/MM, and PCM approaches do we obtain an improved agreement with the experiment. The study sheds light more generally on factors that can impact the computed frequencies and intensities of IR bands in solution.

Investigation of the effects of solvent on oxygen evolution reactions on the surface of magnesium oxide

In this research, due to the importance of oxygen evolution reactions (OER) in the production of hydrogen gas, the effects of the solvent on the OER steps on the magnesium oxide surface have been studied. For this purpose, the difference in free energy value in OER steps on the Magnesium oxide surface for two systems (in the presence and the absence of solvent) have been investigated using molecular dynamics (MD) simulation and thermodynamic integration (TI). These results are in perfect agreement with the results obtained from ab initio molecular dynamics simulations and quantum mechanics calculations. Our findings show that the presence of water solvent around the surface of magnesium oxide has a crucial role on OER and leads to an increase more than twofold in the free energy of all steps. Also, the presence of the solvent has the most effect on the third step of OER, and its ΔG increases by about 1.9 eV. This fact can be attributed to the approach of the water molecules to the substrate spontaneously, which increases the probability of forming hydrogen bonds and the number of contacts, leading to more favorable thermodynamic reactions. Close inspection of the calculated binding energies between the substrate and the intermediates confirms that the binding energy of OER's second step is significantly higher than the other steps and is the velocity determining step.

Tailoring parameters for QM/MM simulations: accurate modeling of adsorption and catalysis in zirconium-based metal-organic frameworks.

Quantum mechanics/molecular mechanics (QM/MM) simulations offer an efficient way to model reactions occurring in complex environments. This study introduces a specialized set of charge and Lennard-Jones parameters tailored for electrostatically embedded QM/MM calculations, aiming to accurately model both adsorption processes and catalytic reactions in zirconium-based metal-organic frameworks (Zr-MOFs). To validate our approach, we compare adsorption energies derived from QM/MM simulations against experimental results and Monte Carlo simulation outcomes. The developed parameters showcase the ability of QM/MM simulations to represent long-range electrostatic and van der Waals interactions faithfully. This capability is evidenced by the prediction of adsorption energies with a low root mean square error of 1.1 kcal mol-1 across a wide range of adsorbates. The practical applicability of our QM/MM model is further illustrated through the study of glucose isomerization and epimerization reactions catalyzed by two structurally distinct Zr-MOF catalysts, UiO-66 and MOF-808. Our QM/MM calculations closely align with experimental activation energies. Importantly, the parameter set introduced here is compatible with the widely used universal force field (UFF). Moreover, we thoroughly explore how the size of the cluster model and the choice of density functional theory (DFT) methodologies influence the simulation outcomes. This work provides an accurate and computationally efficient framework for modeling complex catalytic reactions within Zr-MOFs, contributing valuable insights into their mechanistic behaviors and facilitating further advancements in this dynamic area of research.

EXPLORANDO O MECANISMO DE DEGRADAÇÃO ENZIMÁTICA DO POLIURETANO (PU) POR MEIO DAS POLIURETANASES ESTERASES A E B (PueA E PueB): UMA ABORDAGEM INTEGRADA DE DOCKING, DINÂMICA MOLECULAR (MD) E QM/MM

EXPLORING THE MECHANISM OF ENZYMATIC DEGRADATION OF POLYURETHANE (PU) BY MEANS OF POLYURETHANASE ESTERASES A AND B (PueA AND PueB): AN INTEGRATED APPROACH OF DOCKING, MOLECULAR DYNAMICS (MD) AND QM/MM. The continuous global production of polyurethane (PU) and its environmental persistence generate negative impacts on ecosystems, driving the search for sustainable handling solutions, specifically the use of enzymes. Therefore, this study seeks to elucidate the understanding of the enzymatic degradation mechanism of polyurethanase esterases A and B (PueA and PueB) applied to PU-MDI (methylene diphenyl diisocyanate). In this sense, previously built structures were used and dockings were carried out for the PueA and PueB complexes. The results indicate binding scores of –10.743 and –9.587 kcal mol-1 for PueA and PueB, respectively. Furthermore, the complexes were then subjected to molecular dynamics simulations (MD) for 500 ns in triplicate, showing a good interaction between the enzymes and the PU-MDI. Finally, a QM/MM (quantum mechanics/molecular mechanics) simulation, starting with an adequate snapshot from classical MD was carried out, and a relaxed scan was performed to obtain the occurrence coordinates, followed by umbrella sampling to obtain the free energy profile of the systems. The results indicate that the PueA and PueB reaction mechanisms occur in two steps, with determining energy barriers of 11.95 and 15.10 kcal mol-1, respectively, indicating that these enzymes may be used to degrade PU.

Histidine-assisted reduction of arylnitrenes upon photo-activation of phenyl azide chromophores in GFP-like fluorescent proteins.

The photochemically active sites of the proteins sfGFP66azF and Venus66azF, members of the green fluorescent protein (GFP) family, contain a non-canonical amino acid residue p-azidophenylalanine (azF) instead of Tyr66. The light-induced decomposition of azF at these sites leads to the formation of reactive arylnitrene (nF) intermediates followed by the formation of phenylamine-containing chromophores. We report the first study of the reaction mechanism of the reduction of the arylnitrene intermediates in sfGFP66nF and Venus66nF using molecular modeling methods. The Gibbs energy profiles for the elementary steps of the chemical reaction in sfGFP66nF are computed using molecular dynamics simulations with quantum mechanics/molecular mechanics (QM/MM) potentials. Structures and energies along the reaction pathway in Venus66nF are evaluated using a QM/MM approach. According to the results of the simulations, arylnitrene reduction is coupled with oxidation of the histidine side chain on the His148 residue located near the chromophore.

Mechanism of a novel metal-free carbonic anhydrase.

The recently discovered metal-free carbonic anhydrase (CA) enzyme may significantly impact the global carbon dioxide (CO2) cycle, as it can irreversibly perform the CO2 hydration reaction. In this study, we investigated several key aspects of metal-free CA, including the identification of the catalytic site, the determination of the CO2 binding site, and the mechanism of catalysis. This is achieved through classical molecular dynamics (MD) simulations, quantum chemical density functional theory (DFT), and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Our study indicates that the experimental structure based on X-ray crystallography, which shows the 'bicarbonate (HCO3-) product' trapped in the hydrophilic region of metal-free CA, might not accurately depict the actual enzyme-substrate interaction. Instead, the simulation reveals that CO2 prefers the hydrophobic zone, which serves as the primary catalytic site. It also highlights the strategic role of a gatekeeper residue (Phe504), which assists in regulating the transportation of CO2 by tilting its aromatic plane. Additionally, the hybrid QM/MM calculations establish that CO2 hydration is catalyzed within the hydrophobic zone by a deprotonated tyrosine with the help of an organized water chain.

Ultrafast dynamics in spatially confined photoisomerization: accelerated simulations through machine learning models.

This study sheds light on the exploration of photoresponsive host-guest systems, highlighting the intricate interplay between confined spaces and photosensitive guest molecules. Conducting nonadiabatic molecular dynamics (NAMD) simulations based on electronic structure calculations for such large systems remains a formidable challenge. By leveraging machine learning (ML) as an accelerator for NAMD simulations, we analytically constructed excited-state potential energy surfaces along relevant collective variables to investigate photoisomerization processes efficiently. Combining the quantum mechanics/molecular mechanics (QM/MM) methodology with ML-based NAMD simulations, we elucidated the reaction pathways and identified the key degrees of freedom as reaction coordinates leading to conical intersections. A machine learning-based nonadiabatic dynamics model has been developed to compare the excited-state dynamics of the guest molecule, benzopyran, in both the gas phase and its behavior within the confined space of cucurbit[5]uril. This comparative analysis was designed to determine the influence of the environment on the photoisomerization rate of the guest molecule. The results underscore the effectiveness of ML models in simulating trajectory evolution in a cost-effective manner. This research offers a practical approach to accelerate NAMD simulations in large-scale systems of photochemical reactions, with potential applications in other host-guest complex systems.

Mechanistic insights into the electron attachment process to guanosine in the presence of arginine.

The attachment of low-energy electrons (LEEs) to DNA biomolecules leads to irreversible damage. However, the behavior of this interaction can be influenced by the presence of amino acids. Herein, we have delved into the mechanism of electron attachment to the guanosine in the presence of arginine. This study used combined molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) approaches to collect and optimize the geometries having hydrogen-bonds (H-bonds) between guanosine and arginine, respectively, followed by atom centered density matrix propagation (ADMP) simulations to assess the electron attachment ability of guanine with and without arginine. The vertical detachment energy (VDE) and natural population analysis (NPA) suggest that the electron attached to guanosine occurs more readily due to the H-bonds between guanosine and arginine. The singly occupied molecular orbitals (SOMOs), VDE, and NPA from ADMP results corroborated the idea that in the presence of arginine, the electron effectively attached to the guanosine moiety while the auto detachment process becomes less probable in the case of arg-guanosine (two-H bonds). However, in the presence of arginine, the dissociative electron attachment (DEA) process for guanosine is exothermic, while in the absence of arginine, it is endothermic. This study provides new insight into the process of radiation damaging biological systems by elucidating the DEA to DNA subunits in the presence of amino acids, paving the way for a deeper understanding of radiation-induced damage in biological systems.

Switching the proton-coupled electron transfer mechanism for non-canonical tyrosine residues in a de novo protein.

The proton-coupled electron transfer (PCET) reactions of tyrosine (Y) are instrumental to many redox reactions in nature. This study investigates how the local environment and the thermodynamic properties of Y influence its PCET characteristics. Herein, 2- and 4-mercaptophenol (MP) are placed in the well-folded α3C protein (forming 2MP-α3C and 4MP-α3C) and oxidized by external light-generated [Ru(L)3]3+ complexes. The resulting neutral radicals are long-lived (>100 s) with distinct optical and EPR spectra. Calculated spin-density distributions are similar to canonical Y˙ and display very little spin on the S-S bridge that ligates the MPs to C32 inside the protein. With 2MP-α3C and 4MP-α3C we probe how proton transfer (PT) affects the PCET rate constants and mechanisms by varying the degree of solvent exposure or the potential to form an internal hydrogen bond. Solution NMR ensemble structures confirmed our intended design by displaying a major difference in the phenol OH solvent accessible surface area (≤∼2% for 2MP and 30-40% for 4MP). Additionally, 2MP-C32 is within hydrogen bonding distance to a nearby glutamate (average O-O distance is 3.2 ± 0.5 Å), which is suggested also by quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations. Neither increased exposure of the phenol OH to solvent (buffered water), nor the internal hydrogen bond, was found to significantly affect the PCET rates. However, the lower phenol pKa values associated with the MP-α3C proteins compared to α3Y provided a sufficient change in PT driving force to alter the PCET mechanism. The PCET mechanism for 2MP-α3C and 4MP-α3C with moderately strong oxidants was predominantly step-wise PTET for pH values, but changed to concerted PCET at neutral pH values and below when a stronger oxidant was used, as found previously for α3Y. This shows how the balance of ET and PT driving forces is critical for controlling PCET mechanisms. The presented results improve our general understanding of amino-acid based PCET in enzymes.