The adaptability of the ion-binding site by the Ag(I)/Cu(I) periplasmic chaperone SilF.

The periplasmic chaperone SilF has been identified as part of an Ag(I) detoxification system in Gram-negative bacteria. Sil proteins also bind Cu(I) but with reported weaker affinity, therefore leading to the designation of a specific detoxification system for Ag(I). Using isothermal titration calorimetry, we show that binding of both ions is not only tighter than previously thought but of very similar affinities. We investigated the structural origins of ion binding using molecular dynamics and QM/MM simulations underpinned by structural and biophysical experiments. The results of this analysis showed that the binding site adapts to accommodate either ion, with key interactions with the solvent in the case of Cu(I). The implications of this are that Gram-negative bacteria do not appear to have evolved a specific Ag(I) efflux system but take advantage of the existing Cu(I) detoxification system. Therefore, there are consequences for how we define a particular metal resistance mechanism and understand its evolution in the environment.

If It Is Hard, It Is Worth Doing: Engineering Radical Enzymes from Anaerobes.

With a pressing need for sustainable chemistries, radical enzymes from anaerobes offer a shortcut for many chemical transformations and deliver highly sought-after functionalizations such as late-stage C-H functionalization, C-C bond formation, and carbon-skeleton rearrangements, among others. The challenges in handling these oxygen-sensitive enzymes are reflected in their limited industrial exploitation, despite what they may deliver. With an influx of structures and mechanistic understanding, the scope for designed radical enzymes to deliver wanted processes becomes ever closer. Combined with new advances in computational methods and workflows for these complex systems, the outlook for an increased use of radical enzymes in future processes is exciting.

Application of a Synthetic Ferredoxin-Inspired [4Fe4S]-Peptide Maquette as the Redox Partner for an [FeFe]-Hydrogenase.

'Bacterial-type' ferredoxins host a cubane [4Fe4S]2+/+ cluster that enables these proteins to mediate electron transfer and facilitate a broad range of biological processes. Peptide maquettes based on the conserved cluster-forming motif have previously been reported and used to model the ferredoxins. Herein we explore the integration of a [4Fe4S]-peptide maquette into a H2 -powered electron transport chain. While routinely formed under anaerobic conditions, we illustrate by electron paramagnetic resonance (EPR) analysis that these maquettes can be reconstituted under aerobic conditions by using photoactivated NADH to reduce the cluster at 240 K. Attempts to tune the redox properties of the iron-sulfur cluster by introducing an Fe-coordinating selenocysteine residue were also explored. To demonstrate the integration of these artificial metalloproteins into a semi-synthetic electron transport chain, we utilize a ferredoxin-inspired [4Fe4S]-peptide maquette as the redox partner in the hydrogenase-mediated oxidation of H2 .

Activation of Glycyl Radical Enzymes─Multiscale Modeling Insights into Catalysis and Radical Control in a Pyruvate Formate-Lyase-Activating Enzyme.

Pyruvate formate-lyase (PFL) is a glycyl radical enzyme (GRE) playing a pivotal role in the metabolism of strict and facultative anaerobes. Its activation is carried out by a PFL-activating enzyme, a member of the radical S-adenosylmethionine (rSAM) superfamily of metalloenzymes, which introduces a glycyl radical into the Gly radical domain of PFL. The activation mechanism is still not fully understood and is structurally based on a complex with a short model peptide of PFL. Here, we present extensive molecular dynamics simulations in combination with quantum mechanics/molecular mechanics (QM/MM)-based kinetic and thermodynamic reaction evaluations of a more complete activation model comprising the 49 amino acid long C-terminus region of PFL. We reveal the benefits and pitfalls of the current activation model, providing evidence that the bound peptide conformation does not resemble the bound protein-protein complex conformation with PFL, with implications for the activation process. Substitution of the central glycine with (S)- and (R)-alanine showed excellent binding of (R)-alanine over unstable binding of (S)-alanine. Radical stabilization calculations indicate that a higher radical stability of the glycyl radical might not be the sole origin of the evolutionary development of GREs. QM/MM-derived radical formation kinetics further demonstrate feasible activation barriers for both peptide and C-terminus activation, demonstrating why the crystalized model peptide system is an excellent inhibitory system for natural activation. This new evidence supports the theory that GREs converged on glycyl radical formation due to the better conformational accessibility of the glycine radical loop, rather than the highest radical stability of the formed peptide radicals.

Effect of Oriented Electric Fields on Biologically Relevant Iron-Sulfur Clusters: Tuning Redox Reactivity for Catalysis.

Enzyme-based iron-sulfur clusters, exemplified in families such as hydrogenases, nitrogenases, and radical S-adenosylmethionine enzymes, feature in many essential biological processes. The functionality of biological iron-sulfur clusters extends beyond simple electron transfer, relying primarily on the redox activity of the clusters, with a remarkable diversity for different enzymes. The active-site structure and the electrostatic environment in which the cluster resides direct this redox reactivity. Oriented electric fields in enzymatic active sites can be significantly strong, and understanding the extent of their effect on iron-sulfur cluster reactivity can inform first steps toward rationally engineering their reactivity. An extensive systematic density functional theory-based screening approach using OPBE/TZP has afforded a simple electric field-effect representation. The results demonstrate that the orientation of an external electric field of strength 28.8 MV cm-1 at the center of the cluster can have a significant effect on its relative stability in the order of 35 kJ mol-1. This shows clear implications for the reactivity of iron-sulfur clusters in enzymes. The results also demonstrate that the orientation of the electric field can alter the most stable broken-symmetry state, which further has implications on the directionality of initiated electron-transfer reactions. These insights open the path for manipulating the enzymatic redox reactivity of iron-sulfur cluster-containing enzymes by rationally engineering oriented electric fields within the enzymes.

Controlling the outcome of SN2 reactions in ionic liquids: from rational data set design to predictive linear regression models.

Rate constants for a bimolecular nucleophilic substitution (SN2) process in a range of ionic liquids are correlated with calculated parameters associated with the charge localisation on the cation of the ionic liquid (including the molecular electrostatic potential). Simple linear regression models proved effective, though the interdependency of the descriptors needs to be taken into account when considering generality. A series of ionic liquids were then prepared and evaluated as solvents for the same process; this data set was rationally chosen to incorporate homologous series (to evaluate systematic variation) and functionalities not available in the original data set. These new data were used to evaluate and refine the original models, which were expanded to include simple artificial neural networks. Along with showing the importance of an appropriate data set and the perils of overfitting, the work demonstrates that such models can be used to reliably predict ionic liquid solvent effects on an organic process, within the limits of the data set.

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 Mg2+, Ca2+, 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.

Structural elements of an NRPS cyclization domain and its intermodule docking domain.

Epothilones are thiazole-containing natural products with anticancer activity that are biosynthesized by polyketide synthase (PKS)-nonribosomal peptide synthetase (NRPS) enzymes EpoA-F. A cyclization domain of EpoB (Cy) assembles the thiazole functionality from an acetyl group and l-cysteine via condensation, cyclization, and dehydration. The PKS carrier protein of EpoA contributes the acetyl moiety, guided by a docking domain, whereas an NRPS EpoB carrier protein contributes l-cysteine. To visualize the structure of a cyclization domain with an accompanying docking domain, we solved a 2.03-Å resolution structure of this bidomain EpoB unit, comprising residues M1-Q497 (62 kDa) of the 160-kDa EpoB protein. We find that the N-terminal docking domain is connected to the V-shaped Cy domain by a 20-residue linker but otherwise makes no contacts to Cy. Molecular dynamic simulations and additional crystal structures reveal a high degree of flexibility for this docking domain, emphasizing the modular nature of the components of PKS-NRPS hybrid systems. These structures further reveal two 20-Å-long channels that run from distant sites on the Cy domain to the active site at the core of the enzyme, allowing two carrier proteins to dock with Cy and deliver their substrates simultaneously. Through mutagenesis and activity assays, catalytic residues N335 and D449 have been identified. Surprisingly, these residues do not map to the location of the conserved HHxxxDG motif in the structurally homologous NRPS condensation (C) domain. Thus, although both C and Cy domains have the same basic fold, their active sites appear distinct.

Ion-Reagent Interactions Contributing to Ionic Liquid Solvent Effects on a Condensation Reaction.

Molecular dynamics simulations of solutions of hexan-1-amine or 4-methoxybenzaldehyde in acetonitrile, an ionic liquid/acetonitrile mixture (χIL =0.2), and a number of different (neat) ionic liquids were performed, to further understand the solvent effects on the condensation reaction of these species. This work indicates that, in the presence of an ionic liquid, the amine group of hexan-1-amine is exclusively solvated by the components of the ionic liquid, and not by acetonitrile, and that the anion interacts with the aldehyde group of 4-methoxybenzaldehyde. These interactions showed little dependence on the proportion of the ionic liquid present. When varying the cation of the ionic liquid there were changes in the cation-amine interaction, and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethanesulfonyl)imide ([Bm2 im][N(CF3 SO2 )2 ]) was found to order more than expected about the amine. This ordering is likely the origin of the large rate constant values determined in [Bm2 im][N(CF3 SO2 )2 ] for this condensation reaction and explains an anomaly seen previously. When changing the anion, changes were seen in the interactions between both the cation and anion with hexan-1-amine, and the anion with 4-methoxybenzaldehyde. The differing magnitude of these interactions likely causes subtle changes in the activation parameters for this condensation reaction, and provides an explanation for the anomalous rate constant values previously determined when varying the anion.

Improved NOE fitting for flexible molecules based on molecular mechanics data - a case study with S-adenosylmethionine.

The use of molecular dynamics (MD) calculations to derive relative populations of conformers is highly sensitive to both timescale and parameterisation of the MD. Where these calculations are coupled with NOE data to determine the dynamics of a molecular system, this can present issues if these populations are thus relied upon. We present an approach that refines the highly accurate PANIC NMR methodology combined with clustering approaches to generate conformers, but without restraining the simulations or considering the relative population distributions generated by MD. Combining this structural sampling with NOE fitting, we demonstrate, for S-adenosylmethionine (aqueous solution at pH 7.0), significant improvements are made to the fit of populations to the experimental data, revealing a strong overall preference for the syn conformation of the adenosyl group relative to the ribose ring, but with less discrimination for the conformation of the ribose ring itself.

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.

Developing energy efficient lignin biomass processing - towards understanding mediator behaviour in ionic liquids.

Environmental concerns have brought attention to the requirement for more efficient and renewable processes for chemicals production. Lignin is the second most abundant natural polymer, and might serve as a sustainable resource for manufacturing fuels and aromatic derivatives for the chemicals industry after being depolymerised. In this work, the mediator 2,2'-azino-bis(3-ethylbenthiazoline-6-sulfonic acid) diammonium salt (ABTS), commonly used with enzyme degradation systems, has been evaluated by means of cyclic voltammetry (CV) for enhancing the oxidation of the non-phenolic lignin model compound veratryl alcohol and three types of lignin (organosolv, Kraft and lignosulfonate) in the ionic liquid 1-ethyl-3-methylimidazolium ethyl sulfate, ([C2mim][C2SO4]). The presence of either veratryl alcohol or organosolv lignin increased the second oxidation peak of ABTS under select conditions, indicating the ABTS-mediated oxidation of these molecules at high potentials in [C2mim][C2SO4]. Furthermore, CV was applied as a quick and efficient way to explore the impact of water in the ABTS-mediated oxidation of both organosolv and lignosulfonate lignin. Higher catalytic efficiencies of ABTS were observed for lignosulfonate solutions either in sodium acetate buffer or when [C2mim][C2SO4] (15 v/v%) was present in the buffer solution, whilst there was no change found in the catalytic efficiency of ABTS in [C2mim][C2SO4]-lignosulfonate mixtures relative to ABTS alone. In contrast, organosolv showed an initial increase in oxidation, followed by a significant decrease on increasing the water content of a [C2mim][C2SO4] solution.

Three dimensional ink-jet printing of biomaterials using ionic liquids and co-solvents.

1-Ethyl-3-methylimidazolium acetate ([C2C1Im][OAc]) and 1-butyl-3-methylimidazolium acetate ([C4C1Im][OAc]) have been used as solvents for the dissolution and ink-jet printing of cellulose from 1.0 to 4.8 wt%, mixed with the co-solvents 1-butanol and DMSO. 1-Butanol and DMSO were used as rheological modifiers to ensure consistent printing, with DMSO in the range of 41-47 wt% producing samples within the printable range of a DIMATIX print-head used (printability parameter < 10) at 55 °C, whilst maintaining cellulose solubility. Regeneration of cellulose from printed samples using water was demonstrated, with the resulting structural changes to the cellulose sample assessed by scanning electron microscopy (SEM) and white light interferometry (WLI). These results indicate the potential of biorenewable materials to be used in the 3D additive manufacture process to generate single-component and composite materials.

Ion-combination specific effects driving the enzymatic activity of halophilic alcohol dehydrogenase 2 from Haloferax volcanii in aqueous ionic liquid solvent mixtures.

Biocatalysis in ionic liquids enables novel routes for bioprocessing. Enzymes derived from extremophiles promise greater stability and activity under ionic liquid (IL) influence. Here, we probe the enzyme alcohol dehydrogenase 2 from the halophilic archaeon Haloferax volcanii in thirteen different ion combinations for relative activity and analyse the results against molecular dynamics (MD) simulations of the same IL systems. We probe the ionic liquid property space based on ion polarizability and molecular electrostatic potential. Using the radial distribution functions, survival probabilities and spatial distribution functions of ions, we show that cooperative ion-ion interactions determine ion-protein interactions, and specifically, strong ion-ion interactions equate to higher enzymatic activity if neither of the ions interact strongly with the protein surface. We further demonstrate a tendency for cations interacting with the protein surface to be least detrimental to enzymatic activity if they show a low polarizability when combined with small hydrophilic anions. We also find that the IL ion influence is not mitigated by the surplus of negatively charged residues of the halophilic enzyme. This is shown by free energy landscape analysis in root mean square deviation and distance variation plots of active site gating residues (Trp43 and His273) demonstrating no protection of specific structural elements relevant to preserving enzymatic activity. On the other hand, we observe a general effect across all IL systems that a tight binding of water at acidic residues is preferentially interrupted at these residues through the increased presence of potassium ions. Overall, this study demonstrates a co-ion interaction dependent influence on allosteric surface residues controlling the active/inactive conformation of halophilic alcohol dehydrogenase 2 and the necessity to engineer ionic liquid systems for enzymes that rely on the integrity of functional surface residues regardless of their halophilicity or thermophilicity for use in bioprocessing.

Structural diversity in the AdoMet radical enzyme superfamily.

AdoMet radical enzymes are involved in processes such as cofactor biosynthesis, anaerobic metabolism, and natural product biosynthesis. These enzymes utilize the reductive cleavage of S-adenosylmethionine (AdoMet) to afford l-methionine and a transient 5'-deoxyadenosyl radical, which subsequently generates a substrate radical species. By harnessing radical reactivity, the AdoMet radical enzyme superfamily is responsible for an incredible diversity of chemical transformations. Structural analysis reveals that family members adopt a full or partial Triose-phosphate Isomerase Mutase (TIM) barrel protein fold, containing core motifs responsible for binding a catalytic [4Fe-4S] cluster and AdoMet. Here we evaluate over twenty structures of AdoMet radical enzymes and classify them into two categories: 'traditional' and 'ThiC-like' (named for the structure of 4-amino-5-hydroxymethyl-2-methylpyrimidine phosphate synthase (ThiC)). In light of new structural data, we reexamine the 'traditional' structural motifs responsible for binding the [4Fe-4S] cluster and AdoMet, and compare and contrast these motifs with the ThiC case. We also review how structural data combine with biochemical, spectroscopic, and computational data to help us understand key features of this enzyme superfamily, such as the energetics, the triggering, and the molecular mechanisms of AdoMet reductive cleavage. This article is part of a Special Issue entitled: Radical SAM Enzymes and Radical Enzymology.

Does the cation really matter? The effect of modifying an ionic liquid cation on an SN2 process.

The rate of reaction of a Menschutkin process in a range of ionic liquids with different cations was investigated, with temperature-dependent kinetic data giving access to activation parameters for the process in each solvent. These data, along with molecular dynamics simulations, demonstrate the importance of accessibility of the charged centre on the cation and that the key interactions are of a generalised electrostatic nature.

Probing the importance of ionic liquid structure: a general ionic liquid effect on an S(N)Ar process.

The effect of a range of ionic liquids, with systematic variations in the cation and anion, on the rate constant of an aromatic substitution process was investigated. Temperature-dependent kinetic data allowed calculation of activation parameters for the process in each solvent. These data demonstrate a generalised ionic liquid effect, with an increase in rate constant observed in each ionic solvent, though the microscopic origins of the rate constant enhancement differ with the nature of the ionic liquid.

Anion-π interactions influence pK(a) values.

Five 8-(4-R-phenyl)-1-naphthol derivatives were prepared by PdCl(2)-catalysed electrophilic aromatic substitution. The pK(a)' values for these 1,8-disubstituted arene naphthols have been measured in acetonitrile/water (R = NO(2), 8.42; R = Cl, 8.52; R = H, 8.56; R = Me 8.68; and R = OMe, 8.71) and indicate a correlation with the electronic nature of the arene substituent, as determined through LFER analysis. Contributions to the relative pK(a)' values have been interpreted, using M06-2X DFT calculations, as consisting of two components: A small contribution from initial OH-π bonding in the starting materials and a larger contribution from anion-π interactions in the products. Such effects have implications for a range of other systems.

Ionic liquids through the looking glass: theory mirrors experiment and provides further insight into aromatic substitution processes.

Aromatic substitution of 1-fluoro-2,4-dinitrobenzene has been shown to proceed more rapidly in an ionic liquid than in ethanol, primarily due to entropic, rather than enthalpic, effects. By using molecular dynamics simulations, an anion-pi interaction appears to contribute to strong ordering within the starting material complex that accounts for the relative differences in the S(N)Ar k.

Solvent reorganisation as the driving force for rate changes of Menschutkin reactions in an ionic liquid.

The effect on the rate of reaction of each of a series of Menschutkin processes on changing from a molecular solvent to an ionic liquid was investigated. In each case, the rate acceleration observed at room temperature could be attributed to the change in the entropy of the system on reaching the transition state, offsetting any enthalpic cost.