On the Role of Molecular Conformation of the 8-Oxoguanine Lesion in Damaged DNA Processing by Polymerases.

A common and insidious DNA damage is 8-oxoguanine (8OG), bypassed with low catalytic efficiency and high error frequency by polymerases (Pols) during DNA replication. This is a fundamental process with far-reaching implications in cell function and diseases. However, the molecular determinants of how 8OG exactly affects the catalytic efficiency of Pols remain largely unclear. By examining ternary deoxycytidine triphosphate/DNA/Pol complexes containing the 8OG damage, we found that 8OG consistently adopts different conformations when bound to Pols, compared to when in isolated DNA. Equilibrium molecular dynamics and metadynamics free energy calculations quantified that 8OG is in the lowest energy conformation in isolated DNA. In contrast, 8OG adopts high-energy conformations often characterized by intramolecular steric repulsion when bound to Pols. We show that the 8OG conformation can be regulated by mutating Pol residues interacting with the 8OG phosphate group. These findings propose the 8OG conformation as a factor in Pol-mediated processing of damaged DNA.

Kinetic and molecular dynamics study of inhibition and transglycosylation in Hypocrea jecorina family 3 ß-glucosidases.

ß-Glucosidases enhance enzymatic biomass conversion by relieving cellobiose inhibition of endoglucanases and cellobiohydrolases. However, the susceptibility of these enzymes to inhibition and transglycosylation at high glucose or cellobiose concentrations severely limits their activity and, consequently, the overall efficiency of enzyme mixtures. We determined the impact of these two processes on the hydrolytic activity of the industrially relevant family 3 ß-glucosidases from Hypocrea jecorina, HjCel3A and HjCel3B, and investigated the underlying molecular mechanisms through kinetic studies, binding free energy calculations, and molecular dynamics (MD) simulations. HjCel3B had a 7-fold higher specificity for cellobiose than HjCel3A but greater tendency for glucose inhibition. Energy decomposition analysis indicated that cellobiose has relatively weak electrostatic interactions with binding site residues, allowing it to be easily displaced by glucose and free to inhibit other hydrolytic enzymes. HjCel3A is, thus, preferable as an industrial ß-glucosidase despite its lower activity caused by transglycosylation. This competing pathway to hydrolysis arises from binding of glucose or cellobiose at the product site after formation of the glycosyl-enzyme intermediate. MD simulations revealed that binding is facilitated by hydrophobic interactions with Trp-37, Phe-260, and Tyr-443. Targeting these aromatic residues for mutation to reduce substrate affinity at the product site would therefore potentially mitigate transglycosidic activity. Engineering improved variants of HjCel3A and other structurally similar ß-glucosidases would have a significant economic effect on enzymatic biomass conversion in terms of yield and production cost as the process can be consequently conducted at higher substrate loadings.

Molecular Determinants of Substrate Affinity and Enzyme Activity of a Cytochrome P450BM3 Variant.

Cytochrome P450BM3 catalyzes the hydroxylation and/or epoxidation of fatty acids, fatty amides, and alcohols. Protein engineering has produced P450BM3 variants capable of accepting drug molecules normally metabolized by human P450 enzymes. The enhanced substrate promiscuity has been attributed to the greater flexibility of the lid of the substrate channel. However, it is not well understood how structurally different and highly polar drug molecules can stably bind in the active site nor how the activity and coupling efficiency of the enzyme may be affected by the lack of enzyme-substrate complementarity. To address these important aspects of non-native small molecule binding, this study investigated the binding of drug molecules with different size, charge, polar surface area, and human P450 affinity on the promiscuous R47L/F87V/L188Q/E267V/F81I pentuple mutant of P450BM3. Binding free energy data and energy decomposition analysis showed that pentuple mutant P450BM3 stably binds (i.e., negative ΔGb°) a broad range of substrate and inhibitor types because dispersion interactions with active site residues overcome unfavorable repulsive and electrostatic effects. Molecular dynamics simulations revealed that 1) acidic substrates tend to disrupt the heme propionate A-K69 salt bridge, which may reduce heme oxidizing ability, and 2) the lack of complementarity leads to high substrate mobility and water density in the active site, which may lead to uncoupling. These factors must be considered in future developments of P450BM3 as a biocatalyst in the large-scale production of drug metabolites.

The role of catalytic residue pKa on the hydrolysis/transglycosylation partition in family 3 ß-glucosidases.

ß-Glucosidases (ßgls) primarily catalyze the hydrolysis of the terminal glycosidic bond at the non-reducing end of ß-glucosides, although glycosidic bond synthesis (called transglycosylation) can also occur in the presence of another acceptor. In the final reaction step, the glucose product or another substrate competes with water for transfer to the glycosyl-enzyme intermediate. The factors governing the balance between the two pathways are not fully known; however, the involvement of ionizable residues in binding and catalysis suggests that their pKa may play a role. Through constant pH molecular dynamics simulations of a glycoside hydrolase Family 3 (GH3) ßgl, we showed that the pKa of the catalytic acid/base residue, E441, is low (∼2) during either reaction due to E441-R125-E128 and E441-R125-E166 hydrogen bond networks. The low basicity of E441 would reduce its ability to deprotonate the acceptor. This may be less critical for transglycosylation because sugars have a lower deprotonation enthalpy than water. Moreover, their acidity would be increased by hydrogen bonding with R169 at the acceptor binding site. In contrast, no such interaction was observed for catalytic water. The results are likely applicable to other GH3 ßgls because R125, E128, E166, and R169 are conserved. As these enzymes are commonly used in biomass degradation, there is interest in developing variants with enhanced hydrolytic activity. This may be accomplished by elevating the acid/base residue pKa by disrupting its hydrogen bond networks and reducing the affinity and reactivity of a sugar acceptor by mutating R169.

Effect of Mutation and Substrate Binding on the Stability of Cytochrome P450BM3 Variants.

Cytochrome P450BM3 is a heme-containing enzyme from Bacillus megaterium that exhibits high monooxygenase activity and has a self-sufficient electron transfer system in the full-length enzyme. Its potential synthetic applications drive protein engineering efforts to produce variants capable of oxidizing nonnative substrates such as pharmaceuticals and aromatic pollutants. However, promiscuous P450BM3 mutants often exhibit lower stability, thereby hindering their industrial application. This study demonstrated that the heme domain R47L/F87V/L188Q/E267V/F81I pentuple mutant (PM) is destabilized because of the disruption of hydrophobic contacts and salt bridge interactions. This was directly observed from crystal structures of PM in the presence and absence of ligands (palmitic acid and metyrapone). The instability of the tertiary structure and heme environment of substrate-free PM was confirmed by pulse proteolysis and circular dichroism, respectively. Binding of the inhibitor, metyrapone, significantly stabilized PM, but the presence of the native substrate, palmitic acid, had no effect. On the basis of high-temperature molecular dynamics simulations, the lid domain, ß-sheet 1, and Cys ligand loop (a ß-bulge segment connected to the heme) are the most labile regions and, thus, potential sites for stabilizing mutations. Possible approaches to stabilization include improvement of hydrophobic packing interactions in the lid domain and introduction of new salt bridges into ß-sheet 1 and the heme region. An understanding of the molecular factors behind the loss of stability of P450BM3 variants therefore expedites site-directed mutagenesis studies aimed at developing thermostability.

A DFT and ONIOM study of C-H hydroxylation catalyzed by nitrobenzene 1,2-dioxygenase.

A detailed description of the mechanism of C-H hydroxylation by Rieske non-heme iron dioxygenases remains elusive, as the nature of the oxidizing species is not definitively known. DFT calculations on cluster models of nitrobenzene 1,2-dioxygenase were done to explore possible mechanisms arising from oxidation by either the experimentally observed Fe(III)-OOH complex or the putative high-valent HO-Fe(V)=O intermediate formed through a heterolytic O-O bond cleavage. Hydrogen abstraction by HO-Fe(V)=O, followed by oxygen rebound, was found to be consistent with experimental studies. The findings from the quantum mechanical cluster approach were verified by accounting for the effect of the protein environment on transition state geometries and reaction barriers through ONIOM calculations.

Computer-aided identification, synthesis, and biological evaluation of DNA polymerase η inhibitors for the treatment of cancer.

In cancer cells, Pol η allows DNA replication and cell proliferation even in the presence of chemotherapeutic drug-induced damages, like in the case of platinum-containing drugs. Inhibition of Pol η thus represents a promising strategy to overcome drug resistance and preserve the effectiveness of chemotherapeutic drugs. Here, we report the discovery of a novel class of Pol ƞ inhibitors, with 35 active close analogs. Compound 21 (ARN24964) stands out as the best inhibitor, with an IC50 value of 14.7 µM against Pol η and a good antiproliferative activity when used in combination with cisplatin - with a synergistic effect in three different cancer cell lines (A375, A549, OVCAR3). Moreover, it is characterized by a favorable drug-like profile in terms of its aqueous kinetic solubility, plasma and metabolic stability. Thus, ARN24964 is a promising compound for further structure-based drug design efforts toward developing drugs to solve or limit the issue of drug resistance to platinum-containing drugs in cancer patients.

Fluorogenic sensing of CH3CO2- and H2PO4- by ditopic receptor through conformational change.

Cyclo-bis-(urea-3,6-dichlorocarbazole) (1) forms a 1 : 2 complex with CH(3)CO(2)(-) and H(2)PO(4)(-) through hydrogen bonding with the two urea moieties, resulting in fluorescence enhancement via a combined photoinduced electron transfer (PET) and energy transfer mechanism. The binding mechanism involves a conformational change of the two urea receptors to a trans orientation after binding of the first anion, which facilitates the second interaction.

Alchemical Free-Energy Calculations of Watson-Crick and Hoogsteen Base Pairing Interconversion in DNA.

Hoogsteen (HG) base pairs have a transient nature and can be structurally similar to Watson-Crick (WC) base pairs, making their occurrence and thermodynamic stability difficult to determine experimentally. Herein, we employed the restrain-free-energy perturbation-release (R-FEP-R) method to calculate the relative free energy of the WC and HG base pairing modes in isolated and bound DNA systems and predict the glycosyl torsion conformational preference of purine bases. Notably, this method does not require prior knowledge of the transition pathway between the two end states. Remarkably, relatively fast convergence was reached, with results in excellent agreement with experimental data for all the examined DNA systems. The R-REP-R method successfully determined the stability of HG base pairing and more generally, the conformational preference of purine bases, in these systems. Therefore, this computational approach can help to understand the dynamic equilibrium between the WC and HG base pairing modes in DNA.

Nature of anion-templated π(+)-π(+) interactions.

Interaction between positively charged aromatic groups (π(+)-π(+)) is characterized by anti-parallel, displaced-stacked structures in the presence of counteranions. Binding energies of pyridinium, N-methylpyridinium and N-methylimidazolium dimers are much larger than that of benzene-pyridine (π-π) and pyridinium-benzene (π(+)-π). Stabilization is attributed to attractive electrostatic interaction with significant dispersion contribution.

Hydrolysis and Transglycosylation Transition States of Glycoside Hydrolase Family 3 ß-Glucosidases Differ in Charge and Puckering Conformation.

ß-Glucosidases (ßgls) from glycoside hydrolase family 3 play an important role in biomass degradation by catalyzing cellobiose hydrolysis. However, the hydrolysis rate decreases when the glucose product or another cellobiose competes with water to form oligosaccharides in a reaction called transglycosylation. Both reactions involve proton transfer to the acid/base residue and nucleophilic attack on the glycosyl-enzyme intermediate. To gain a deeper understanding of these competing reactions, quantum mechanics/molecular mechanics calculations were performed. Although both reactions are exothermic and have similar free-energy barriers (∼18 kcal/mol), the transition-state (TS) characteristics are different. The glycosyl-water bond is nearly formed in the hydrolysis TS, leading to reduced ionic character and a 4C1 chair conformation. The transglycosylation TS is more positively charged and adopts the 4H3 half-chair conformation because bond formation is less advanced. Water interacts solely with acid/base residue E441, though the long distance between them (2.1 Å) suggests that E441 does not activate water for nucleophilic attack. In comparison, a glucose acceptor has a lower deprotonation enthalpy and hydrogen bonds to E441 (1.6 Å) as well as to Y204, R169, and R67. Knowledge of these factors that are relevant to TS formation and stability is valuable for engineering ßgls with enhanced hydrolytic activity for industrial applications.

Desulfination by 2'-hydroxybiphenyl-2-sulfinate desulfinase proceeds via electrophilic aromatic substitution by the cysteine-27 proton.

Biodesulfurization is an attractive option for enzymatically removing sulfur from the recalcitrant thiophenic derivatives that comprise the majority of organosulfur compounds remaining in hydrotreated petroleum products. Desulfurization in the bacteria Rhodococcus erythropolis follows a four-step pathway culminating in C-S bond cleavage in the 2'-hydroxybiphenyl-2-sulfinate (HBPS) intermediate to yield 2-hydroxybiphenyl and bisulfite. The reaction, catalyzed by 2'-hydroxybiphenyl-2-sulfinate desulfinase (DszB), is the rate-limiting step and also the least understood, as experimental evidence points to a mechanism unlike that of other desulfinases. On the basis of structural and biochemical evidence, two possible mechanisms have been proposed: nucleophilic addition and electrophilic aromatic substitution. Density functional theory calculations showed that electrophilic substitution by a proton is the lower energy pathway and is consistent with previous kinetic and site-directed mutagenesis studies. C27 transfers its proton to HBPS, leading directly to the release of SO2 without the formation of a carbocation intermediate. The H60-S25 dyad stabilizes the transition state by withdrawing the developing negative charge on cysteine. Establishing the desulfination mechanism and specific role of active site residues, accomplished in this study, is essential to protein engineering efforts to increase DszB catalytic activity, which is currently too low for industrial-scale application.

Analyzing sites of OH radical attack (ring vs. side chain) in oxidation of substituted benzenes via dual stable isotope analysis (δ(13)C and δ(2)H).

OH radicals generated by the photolysis of H2O2 can degrade aromatic contaminants by either attacking the aromatic ring to form phenolic products or oxidizing the substituent. We characterized these competing pathways by analyzing the carbon and hydrogen isotope fractionation (εC and εH) of various substituted benzenes. For benzene and halobenzenes that only undergo ring addition, low values of εC (-0.7‰ to -1.0‰) were observed compared with theoretical values (-7.2‰ to -8‰), possibly owing to masking effect caused by pre-equilibrium between the substrate and OH radical preceding the rate-limiting step. In contrast, the addition of OH radicals to nitrobenzene ring showed a higher εC (-3.9‰), probably due to the lower reactivity. Xylene isomers, anisole, aniline, N,N-dimethylaniline, and benzonitrile yielded normal εH values (-2.8‰ to -29‰) indicating the occurrence of side-chain reactions, in contrast to the inverse εH (11.7‰ to 30‰) observed for ring addition due to an sp(2) to sp(3) hybridization change at the reacting carbon. Inverse εH values for toluene (14‰) and ethylbenzene (30‰) were observed despite the formation of side-chain oxidation products, suggesting that ring addition has a larger contribution to isotope fractionation. Dual element isotope slopes (∆δ(2)H/∆δ(13)C) therefore allow identification of significant degradation pathways of aromatic compounds by photochemically induced OH radicals. Issues that should be addressed in future studies include quantitative determination of the contribution of each competing pathway to the observed isotope fractionation and characterization of physical processes preceding the reaction that could affect isotope fractionation.

A DFT study of the cis-dihydroxylation of nitroaromatic compounds catalyzed by nitrobenzene dioxygenase.

The mechanism of cis-dihydroxylation of nitrobenzene and 2-nitrotoluene catalyzed by nitrobenzene 1,2-dioxygenase (NBDO), a member of the naphthalene family of Rieske non-heme iron dioxygenases, was studied by means of the density functional theory method using four models of the enzyme active site. Different possible reaction pathways for the substrate dioxygenation initiated either by the Fe(III)-OOH or HO-Fe(V)═O attack on the aromatic ring were considered and the computed activation barriers compared with the Gibbs free energy of activation for the oxygen-oxygen cleavage leading to the formation of the iron(V)-oxo species from its ferric hydroperoxo precursor. The mechanism of the substrate cis-dihydroxylation leading to the formation of a cis-dihydrodiol was then investigated, and the most feasible mechanism was found to be starting with the attack of the high-valent iron-oxo species on the substrate ring yielding a radical intermediate, which further evolves toward the final product.

Molecular Dynamics Simulation of Nitrobenzene Dioxygenase Using AMBER Force Field.

Molecular dynamics simulation of the oxygenase component of nitrobenzene dioxygenase (NBDO) system, a member of the naphthalene family of Rieske nonheme iron dioxygenases, has been carried out using the AMBER force field combined with a new set of parameters for the description of the mononuclear nonheme iron center and iron-sulfur Rieske cluster. Simulation results provide information on the structure and dynamics of nitrobenzene dioxygenase in an aqueous environment and shed light on specific interactions that occur in its catalytic center. The results suggest that the architecture of the active site is stabilized by key hydrogen bonds, and Asn258 positions the substrate for oxidation. Analysis of protein-water interactions reveal the presence of a network of solvent molecules at the entrance to the active site, which could be of potential catalytic importance.

Cytochrome P450-catalyzed dealkylation of atrazine by Rhodococcus sp. strain NI86/21 involves hydrogen atom transfer rather than single electron transfer.

Cytochrome P450 enzymes are responsible for a multitude of natural transformation reactions. For oxidative N-dealkylation, single electron (SET) and hydrogen atom abstraction (HAT) have been debated as underlying mechanisms. Combined evidence from (i) product distribution and (ii) isotope effects indicate that HAT, rather than SET, initiates N-dealkylation of atrazine to desethyl- and desisopropylatrazine by the microorganism Rhodococcus sp. strain NI86/21. (i) Product analysis revealed a non-selective oxidation at both the αC and ßC-atom of the alkyl chain, which is expected for a radical reaction, but not SET. (ii) Normal (13)C and (15)N as well as pronounced (2)H isotope effects (εcarbon: -4.0‰ ± 0.2‰; εnitrogen: -1.4‰ ± 0.3‰, KIEH: 3.6 ± 0.8) agree qualitatively with calculated values for HAT, whereas inverse (13)C and (15)N isotope effects are predicted for SET. Analogous results are observed with the Fe(iv)[double bond, length as m-dash]O model system [5,10,15,20-tetrakis(pentafluorophenyl)porphyrin-iron(iii)-chloride + NaIO4], but not with permanganate. These results emphasize the relevance of the HAT mechanism for N-dealkylation by P450.

Anion Binding by Electron-Deficient Arenes Based on Complementary Geometry and Charge Distribution.

Extended electron-deficient arenes are investigated as potential neutral receptors for polyanions. Anion binds via σ interaction with extended arenes, which are composed solely of C and N ring atoms and CN substituents. As a result, the positive charge on the aromatic C is enhanced, consequently maximizing binding strength. Selectivity is achieved because different charge distributions can be obtained for target anions of a particular geometry. The halides F(-) and Cl(-) form the most stable complex with 6, while the linear N3(-) interacts most favorably with 7. The trigonal NO3(-) and tetrahedral ClO4(-) fit the 3-fold rotational axis of 6 but do not form stable complexes with 5 and 7. The Y-shaped HCOO(-) forms complexes with 4, 5, and 7, with the latter being the most stable. Thus, the anion complexes exhibit strong binding and the best geometrical fit between guest and host, reminiscent of Lego blocks.

Fluorescent imidazolium-based cyclophane for detection of guanosine-5'-triphosphate and I(-) in aqueous solution of physiological pH.

A new water-soluble and fluorescent imidazolium-anthracene cyclophane (1) effectively recognizes the biologically important GTP and I(-) over other anions in a 100% aqueous solution of physiological pH 7.4. Fluorescence and (1)H NMR spectra and ab initio calculations demonstrate that emission arises from the formation of an excimer state and quenching occurs upon GTP/I(-) binding through (C-H)(+)···A(-) hydrogen bond interactions.

Can Electron-Rich π Systems Bind Anions?

In general, anion-π interactions exist between anions and aromatics with a positive quadrupole moment. The interaction between anions and aromatics with a negative quadrupole moment is expected to be unstable due to Coulombic repulsion. However, here we investigated the cases of aromatics with a negative quadrupole moment such as electron-rich alkyl/alkenyl/alkynyl-substituted benzenes and triphenylene, which interact with halides. Favorable binding was demonstrated with coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)] at the complete basis set (CBS) limit. Stability increases with chain length, unsaturation, and halogenation. Energy decomposition analysis based on symmetry adapted perturbation theory (SAPT) shows that electrostatic repulsion is overcome by induction effects arising from the alkyl substituents.