Enzyme catalytic nitration of aromatic compounds

Nitroaromatic compounds are important intermediates in organic synthesis. The classic method used to synthesize them is chemical nitration, which involves the use of nitric acid diluted in water or acetic acid, both harmful to the environment. With the development of green chemistry, environmental friendly enzyme catalysis is increasingly employed in chemical processes. In this work, we adopted a non-aqueous horseradish peroxidase (HRP)/NaNO₂/H₂O₂ reaction system to study the structural characteristics of aromatic compounds potentially nitrated by enzyme catalysis, as well as the relationship between the charges on carbon atoms in benzene ring and the nitro product distribution. Investigation of various reaction parameters showed that mild reaction conditions (ambient temperature and neutral pH), plus appropriate use of H₂O₂ and NaNO₂ could prevent inactivation of HRP and polymerization of the substrates. Compared to aqueous–organic co-solvent reaction media, the aqueous–organic two-liquid phase system had great advantages in increasing the dissolved concentration of substrate and alleviating substrate inhibition. Analysis of the aromatic compounds’ structural characteristics indicated that substrates containing substituents of NH₂ or OH were readily catalyzed. Furthermore, analysis of the relationship between natural bond orbital (NBO) charges on carbon atoms in benzene ring, as calculated by the density functional method, and the nitro product distribution characteristics, demonstrated that the favored nitration sites were the ortho and para positions of substituents in benzene ring, similar to the selectivity of chemical nitration.     

Computational characterization of the chemical step in the GTP hydrolysis by Ras‐GAP for the wild‐type and G13V mutated Ras

The free energy profiles for the chemical reaction of the guanosine triphosphate hydrolysis GTP + H₂O → GDP + Pi by Ras‐GAP for the wild‐type and G13V mutated Ras were computed by using molecular dynamics protocols with the QM(ab initio)/MM potentials. The results are consistent with the recent measurements of reaction kinetics in Ras‐GAP showing about two‐order reduction of the rate constant upon G13V mutation in Ras: the computed activation barrier on the free energy profile is increased by 3 kcal/mol upon the G13V replacement. The major reason for a higher energy barrier is a shift of the “arginine finger” (R789 from GAP) from the favorable position in the active site. The results of simulations provide support for the mechanism of the reference reaction according to which the Q61 side chain directly participates in chemical transformations at the proton transfer stage. Proteins 2015; 83:1046–1053. © 2015 Wiley Periodicals, Inc.     

Proton reduction at surface of transition metal nanocatalysts

Catalytic activities of neutral and charged palladium (Pd) nanoparticles are compared for hydrogen reduction half-reaction. In this work the sequential H₂ dissociation from the surface of a Pd₁₃H₂₄ cluster is systematically studied by ab initio molecular dynamics (AIMD) at the density functional theory level. AIMD simulation is launched by preparing initial values of momenta of all nuclei in the model corresponding to a temperature range of 0–1700 K. AIMD simulation provides the trajectories of all the atoms in the cluster. A sequential H₂ desorption up to seven molecules is observed from the cluster surface due to thermal motion of nuclei. Modifications of total charge on the neutral Pd₁₃H₂₄ cluster model are found to affect surface H₂ desorption behaviour. A desorption rate of H₂ molecule on both neutral and charged Pd₁₃H₂₄ clusters is compared to the data of Pt₁₃H₂₄ cluster reported previously. The H₂ desorption energy on all the investigated clusters is also determined. The results reveal that Pd₁₃ cluster presents a higher catalytic activity than Pt₁₃ cluster.     

Crystal structure of human insulin‐regulated aminopeptidase with specificity for cyclic peptides

Insulin‐regulated aminopeptidase (IRAP or oxytocinase) is a membrane‐bound zinc‐metallopeptidase that cleaves neuroactive peptides in the brain and produces memory enhancing effects when inhibited. We have determined the crystal structure of human IRAP revealing a closed, four domain arrangement with a large, mostly buried cavity abutting the active site. The structure reveals that the GAMEN exopeptidase loop adopts a very different conformation from other aminopeptidases, thus explaining IRAP's unique specificity for cyclic peptides such as oxytocin and vasopressin. Computational docking of a series of IRAP‐specific cognitive enhancers into the crystal structure provides a molecular basis for their structure–activity relationships and demonstrates that the structure will be a powerful tool in the development of new classes of cognitive enhancers for treating a variety of memory disorders such as Alzheimer's disease.     

Quantitative functional characterization of conserved molecular interactions in the active site of mannitol 2-dehydrogenase

Enzyme active site residues are often highly conserved, indicating a significant role in function. In this study we quantitate the functional contribution for all conserved molecular interactions occurring within a Michaelis complex for mannitol 2-dehydrogenase derived from Pseudomonas fluorescens (pfMDH). Through systematic mutagenesis of active site residues, we reveal that the molecular interactions in pfMDH mediated by highly conserved residues not directly involved in reaction chemistry can be as important to catalysis as those directly involved in the reaction chemistry. This quantitative analysis of the molecular interactions within the pfMDH active site provides direct insight into the functional role of each molecular interaction, several of which were unexpected based on canonical sequence conservation and structural analyses.     

Two-step Ligand Binding in a (βα)8 Barrel Enzyme: SUBSTRATE-BOUND STRUCTURES SHED NEW LIGHT ON THE CATALYTIC CYCLE OF HisA*

HisA is a (βα)₈ barrel enzyme that catalyzes the Amadori rearrangement of N′-[(5′-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (ProFAR) to N′-((5′-phosphoribulosyl) formimino)-5-aminoimidazole-4-carboxamide-ribonucleotide (PRFAR) in the histidine biosynthesis pathway, and it is a paradigm for the study of enzyme evolution. Still, its exact catalytic mechanism has remained unclear. Here, we present crystal structures of wild type Salmonella enterica HisA (SeHisA) in its apo-state and of mutants D7N and D7N/D176A in complex with two different conformations of the labile substrate ProFAR, which was structurally visualized for the first time. Site-directed mutagenesis and kinetics demonstrated that Asp-7 acts as the catalytic base, and Asp-176 acts as the catalytic acid. The SeHisA structures with ProFAR display two different states of the long loops on the catalytic face of the structure and demonstrate that initial binding of ProFAR to the active site is independent of loop interactions. When the long loops enclose the substrate, ProFAR adopts an extended conformation where its non-reacting half is in a product-like conformation. This change is associated with shifts in a hydrogen bond network including His-47, Asp-129, Thr-171, and Ser-202, all shown to be functionally important. The closed conformation structure is highly similar to the bifunctional HisA homologue PriA in complex with PRFAR, thus proving that structure and mechanism are conserved between HisA and PriA. This study clarifies the mechanistic cycle of HisA and provides a striking example of how an enzyme and its substrate can undergo coordinated conformational changes before catalysis.     

Engineered and Native Coenzyme B12-dependent Isovaleryl-CoA/Pivalyl-CoA Mutase

Adenosylcobalamin-dependent isomerases catalyze carbon skeleton rearrangements using radical chemistry. We have recently demonstrated that an isobutyryl-CoA mutase variant, IcmF, a member of this enzyme family that catalyzes the interconversion of isobutyryl-CoA and n-butyryl-CoA also catalyzes the interconversion between isovaleryl-CoA and pivalyl-CoA, albeit with low efficiency and high susceptibility to inactivation. Given the biotechnological potential of the isovaleryl-CoA/pivalyl-CoA mutase (PCM) reaction, we initially attempted to engineer IcmF to be a more proficient PCM by targeting two active site residues predicted based on sequence alignments and crystal structures, to be key to substrate selectivity. Of the eight mutants tested, the F598A mutation was the most robust, resulting in an ∼17-fold increase in the catalytic efficiency of the PCM activity and a concomitant ∼240-fold decrease in the isobutyryl-CoA mutase activity compared with wild-type IcmF. Hence, mutation of a single residue in IcmF tuned substrate specificity yielding an ∼4000-fold increase in the specificity for an unnatural substrate. However, the F598A mutant was even more susceptible to inactivation than wild-type IcmF. To circumvent this limitation, we used bioinformatics analysis to identify an authentic PCM in genomic databases. Cloning and expression of the putative AdoCbl-dependent PCM with an α₂β₂ heterotetrameric organization similar to that of isobutyryl-CoA mutase and a recently characterized archaeal methylmalonyl-CoA mutase, allowed demonstration of its robust PCM activity. To simplify kinetic analysis and handling, a variant PCM-F was generated in which the αβ subunits were fused into a single polypeptide via a short 11-amino acid linker. The fusion protein, PCM-F, retained high PCM activity and like PCM, was resistant to inactivation. Neither PCM nor PCM-F displayed detectable isobutyryl-CoA mutase activity, demonstrating that PCM represents a novel 5′-deoxyadenosylcobalamin-dependent acyl-CoA mutase. The newly discovered PCM and the derivative PCM-F, have potential applications in bioremediation of pivalic acid found in sludge, in stereospecific synthesis of C5 carboxylic acids and alcohols, and in the production of potential commodity and specialty chemicals.     

Reconstitution of Formylglycine-generating Enzyme with Copper(II) for Aldehyde Tag Conversion

To further our aim of synthesizing aldehyde-tagged proteins for research and biotechnology applications, we developed methods for recombinant production of aerobic formylglycine-generating enzyme (FGE) in good yield. We then optimized the FGE biocatalytic reaction conditions for conversion of cysteine to formylglycine in aldehyde tags on intact monoclonal antibodies. During the development of these conditions, we discovered that pretreating FGE with copper(II) is required for high turnover rates and yields. After further investigation, we confirmed that both aerobic prokaryotic (Streptomyces coelicolor) and eukaryotic (Homo sapiens) FGEs contain a copper cofactor. The complete kinetic parameters for both forms of FGE are described, along with a proposed mechanism for FGE catalysis that accounts for the copper-dependent activity.     

Prechemistry versus preorganization in DNA replication fidelity

The molecular origin of nucleotide insertion catalysis and fidelity of DNA polymerases is explored by means of computational simulations. Special attention is paid to the examination of the validity of proposals that invoke prechemistry effects, checkpoints concepts, and dynamical effects. The simulations reproduce the observed fidelity in Pol β, starting with the relevant observed X-ray structures of the complex with the right (R) and wrong (W) nucleotides. The generation of free energy surfaces for the R and W systems also allowed us to analyze different proposals about the origin of the fidelity and to reach several important conclusions. It is found that the potential of mean force (PMF) obtained by proper sampling does not support QM/MM-based proposals of a large barrier before the prechemistry state. Furthermore, examination of dynamical proposals by the renormalization approach indicates that the motions from open to close configurations do not contribute to catalysis or fidelity. Finally we discuss and analyze the induced fit concept and show that, despite its importance, it does not explain fidelity. That is, the fidelity is apparently due to the change in the preorganization of the chemical site, as a result of the relaxation of the binding site upon binding of the incorrect nucleotide. Finally and importantly, since the issue is the barrier associated with the enzyme-substrate (ES)/DNA complex at the chemical transition state and not the path to this complex formation (unless this path involves rate determining steps), it is also not useful to invoke checkpoints while discussing fidelity.     

Calcium-independent calmodulin binding and two-metal-ion catalytic mechanism of anthrax edema factor

Edema factor (EF), a key anthrax exotoxin, has an anthrax protective antigen-binding domain (PABD) and a calmodulin (CaM)-activated adenylyl cyclase domain. Here, we report the crystal structures of CaM-bound EF, revealing the architecture of EF PABD. CaM has N- and C-terminal domains and each domain can bind two calcium ions. Calcium binding induces the conformational change of CaM from closed to open. Structures of the EF-CaM complex show how EF locks the N-terminal domain of CaM into a closed conformation regardless of its calcium-loading state. This represents a mechanism of how CaM effector alters the calcium affinity of CaM and uncouples the conformational change of CaM from calcium loading. Furthermore, structures of EF-CaM complexed with nucleotides show that EF uses two-metal-ion catalysis, a prevalent mechanism in DNA and RNA polymerases. A histidine (H351) further facilitates the catalysis of EF by activating a water to deprotonate 3'OH of ATP. Mammalian adenylyl cyclases share no structural similarity with EF and they also use two-metal-ion catalysis, suggesting the catalytic mechanism-driven convergent evolution of two structurally diverse adenylyl cyclases.