Search

Controlling the selectivity of a reaction is critical for target-oriented synthesis. Accessing complementary selectivity profiles enables divergent synthetic strategies, but is challenging to achieve in biocatalytic reactions given enzymes’ innate preferences of a single selectivity. Thus, it is critical to understand the structural features that control selectivity in biocatalytic reactions to achieve tunable selectivity. Here, we investigate the structural features that control the stereoselectivity in an oxidative dearomatization reaction that is key to making azaphilone natural products. Crystal structures of enantiocomplementary biocatalysts guided the development of multiple hypotheses centered on the structural features that control the stereochemical outcome of the reaction; however, in many cases, direct substitutions of active site residues in natural proteins led to inactive enzymes. Ancestral sequence reconstruction (ASR) and resurrection were employed as an alternative strategy to probe the impact of each residue on the stereochemical outcome of the dearomatization reaction. These studies suggest that two mechanisms are active in controlling the stereochemical outcome of the oxidative dearomatization reaction: one involving multiple active site residues in AzaH and the other dominated by a single Phe to Tyr switch in TropB and AfoD. Moreover, this study suggests that the flavin-dependent monooxygenases (FDMOs) adopt simple and flexible strategies to control stereoselectivity, which has led to stereocomplementary azaphilone natural products produced by fungi. This paradigm of combining ASR and resurrection with mutational and computational studies showcases sets of tools for understanding enzyme mechanisms and provides a solid foundation for future protein engineering efforts.

Acid–base catalysis, which involves one or more proton transfer reactions, is a chemical mechanism commonly employed by many enzymes. The molecular basis for catalysis is often derived from structures determined at the optimal pH for enzyme activity. However, direct observation of protons from experimental structures is quite difficult; thus, a complete mechanistic description for most enzymes remains lacking. Dihydrofolate reductase (DHFR) exemplifies general acid–base catalysis, requiring hydride transfer and protonation of its substrate, DHF, to form the product, tetrahydrofolate (THF). Previous X-ray and neutron crystal structures coupled with theoretical calculations have proposed that solvent mediates the protonation step. However, visualization of a proton transfer has been elusive. Based on a 2.1 Å resolution neutron structure of a pseudo-Michaelis complex of E. coli DHFR determined at acidic pH, we report the direct observation of the catalytic proton and its parent solvent molecule. Comparison of X-ray and neutron structures elucidated at acidic and neutral pH reveals dampened dynamics at acidic pH, even for the regulatory Met20 loop. Guided by the structures and calculations, we propose a mechanism where dynamics are crucial for solvent entry and protonation of substrate. This mechanism invokes the release of a sole proton from a hydronium (H3O+) ion, its pathway through a narrow channel that sterically hinders the passage of water, and the ultimate protonation of DHF at the N5 atom.

HIV-1 protease is indispensable for virus propagation and an important therapeutic target for antiviral inhibitors to treat AIDS. As such inhibitors are transition-state mimics, a detailed understanding of the enzyme mechanism is crucial for the development of better anti-HIV drugs. Here, we used room-temperature joint X-ray/neutron crystallography to directly visualize hydrogen atoms and map hydrogen bonding interactions in a protease complex with peptidomimetic inhibitor KVS-1 containing a reactive nonhydrolyzable ketomethylene isostere, which, upon reacting with the catalytic water molecule, is converted into a tetrahedral intermediate state, KVS-1TI. We unambiguously determined that the resulting tetrahedral intermediate is an oxyanion, rather than the gem-diol, and both catalytic aspartic acid residues are protonated. The oxyanion tetrahedral intermediate appears to be unstable, even though the negative charge on the oxyanion is delocalized through a strong n → π* hyperconjugative interaction into the nearby peptidic carbonyl group of the inhibitor. To better understand the influence of the ketomethylene isostere as a protease inhibitor, we have also examined the protease structure and binding affinity with keto-darunavir (keto-DRV), which similar to KVS-1 includes the ketomethylene isostere. We show that keto-DRV is a significantly less potent protease inhibitor than DRV. These findings shed light on the reaction mechanism of peptide hydrolysis catalyzed by HIV-1 protease and provide valuable insights into further improvements in the design of protease inhibitors.

Non-heme iron (NHI) enzymes perform a variety of oxidative rearrangements to advance simple building blocks toward complex molecular scaffolds within secondary metabolite pathways. Many of these transformations occur with selectivity that is unprecedented in small molecule catalysis, spurring an interest in the enzymatic processes which lead to a particular rearrangement. In-depth investigations of NHI mechanisms examine the source of this selectivity and can offer inspiration for the development of novel synthetic transformations. However, the mechanistic details of many NHI-catalyzed rearrangements remain underexplored, hindering full characterization of the chemistry accessible to this functionally diverse class of enzymes. For NHI-catalyzed rearrangements which have been investigated, mechanistic proposals often describe one-electron processes, followed by single electron oxidation from the substrate to the iron(III)-hydroxyl active site species. Here, we examine the ring expansion mechanism employed in fungal tropolone biosynthesis. TropC, an α-ketoglutarate- dependent NHI dioxygenase, catalyzes a ring expansion in the biosynthesis of tropolone natural product stipitatic acid through an under-studied mechanism. Investigation of both polar and radical mechanistic proposals suggests tropolones are constructed through a radical ring expansion. This biosynthetic route to tropolones is supported by X-ray crystal structure data combined with molecular dynamics simulations, alanine-scanning of active site residues, assessed reactivity of putative biosynthetic intermediates, and quantum mechanical (QM) calculations. These studies support a radical ring expansion in fungal tropolone biosynthesis.

Class A flavin-dependent hydroxylases (FdHs) catalyze the hydroxylation of organic compounds in a site and stereoselective manner. In stark contrast, conventional synthetic routes require environmentally hazardous reagents and give modest yields. Thus, understanding the detailed mechanism of this class of enzymes is essential to their rational manipulation for applications in green chemistry and pharmaceutical production. Both electrophilic substitution and radical intermediate mechanisms have been proposed as interpretations of FdH hydroxylation rates and optical spectra. While radical mechanistic steps are often difficult to examine directly, modern quantum chemistry calculations combined with statistical mechanical approaches can yield detailed mechanistic models providing insight that can be used to differentiate reaction pathways. In the current work, we report quantum mechanical/molecular mechanical (QM/MM) calculations on the fungal TropB enzyme that show an alternative reaction pathway in which hydroxylation through a hydroxyl radical coupled electron transfer (HRC-ET) mechanism is significantly favored over electrophilic substitution. Furthermore, QM/MM calculations on several modified flavins provide a more consistent interpretation of the experimental trends in the reaction rates seen experimentally for a related enzyme, para-hydroxybenzoate hydroxylase (PHBH). These calculations should guide future enzyme and substrate design strategies and broaden the scope of biological spin chemistry.

Myoglobin (Mb) binds oxygen with high affinity as a low spin singlet complex and thus functions as an oxygen storage protein. Yet, hybrid Density Functional Theory/Molecular Mechanical (DFT/MM) calculations of oxy-Mb models predict that the O(2) bond is much less resistant to breaking in the presence of hydrogen sulfide (H(2)S) compared with water. Specifically, a hydrogen atom from H(2)S can be transferred to the distal oxygen atom through homolytic cleavage of the S–H bond to form the intermediate Compound (Cpd) 0 structure and a thiyl radical. In the presence of a neutral His64 (Nε protonation, His64-ε) and H(2)S, only a metastable Cpd 0 would be formed as the active site is devoid of any additional proton donor to fully break the O(2) bond. In contrast, the calculations predict that the triplet state is significantly favored over the open shell singlet diradical state throughout the entire reaction coordinate in the presence of H(2)S and a positively charged His64. Furthermore, a positively charged His64 can readily donate a proton to Cpd 0 to fully break the O(2) bond resulting in a configuration analogous to reported reaction models of a hemoglobin mutant bound to H(2)O(2) with H(2)S present. Typically, exotic techniques are required to generate Cpd 0 but under the conditions just described the intermediate is readily detected in UV–Vis spectra at room temperature. The effect is observed as a 2 nm red shift of the Soret band from 414 nm to 416 nm (pH 5.0, His64-εδ) and from 416 nm to 418 nm (pH 6.6, His64-ε).

Biocatalytic reactions embody many features of ideal chemical transformations, including the potential for impeccable selectivity, high catalytic efficiency, mild reaction conditions and the use of environmentally benign reagents. These advantages have created a demand for biocatalysts that expand the portfolio of complexity-generating reactions available to synthetic chemists. However, the tradeoff that often exists between the substrate scope of a biocatalyst and its selectivity limits the application of enzymes in synthesis. We recently demonstrated that a flavin-dependent monooxygenase, TropB, maintains high levels of site- and stereoselectivity across a range of structurally diverse substrates. Herein, we disclose the structural basis for substrate binding in TropB, which performs a synthetically challenging asymmetric oxidative dearomatization reaction with exquisite site- and stereoselectivity across a range of phenol substrates, providing a foundation for future protein engineering and reaction development efforts. Our hypothesis for substrate binding is informed by a crystal structure of TropB and molecular dynamics simulations with the corresponding computational TropB model and is supported by experimental data. In contrast to canonical class A FAD-dependent monooxygenases in which substrates bind in a protonated form, our data indicate that the phenolate form of the substrate binds in the active site. Furthermore, the substrate position is controlled through twopoint binding of the phenolate oxygen to Arg206 and Tyr239, which are shown to have distinct and essential roles in catalysis. Arg206 is involved in the reduction of the flavin cofactor, suggesting a role in flavin dynamics. Further, QM/MM simulations reveal the interactions that govern the facial selectivity that leads to a highly enantioselective transformation. Thus, the structural origins of the high levels of site-and stereoselectivity observed in reactions of TropB across a range of substrates are elucidated, providing a foundation for future protein engineering and reaction development efforts.

Exposure to organophosphorus compounds (OPs) may be fatal if untreated, and a clear and present danger posed by nerve agent OPs has become palpable in recent years. OPs inactivate acetylcholinesterase (AChE) by covalently modifying its catalytic serine. Inhibited AChE cannot hydrolyze the neurotransmitter acetylcholine leading to its build-up at the cholinergic synapses and creating an acute cholinergic crisis. Current antidotes, including oxime reactivators that attack the OP-AChE conjugate to free the active enzyme, are inefficient. Better reactivators are sought, but their design is hampered by a conformationally rigid portrait of AChE extracted exclusively from 100K X-ray crystallography and scarcity of structural knowledge on human AChE (hAChE). Here, we present room temperature X-ray structures of native and VX-phosphonylated hAChE with an imidazole-based oxime reactivator, RS-170B. We discovered that inhibition with VX triggers substantial conformational changes in bound RS-170B from a “nonproductive” pose (the reactive aldoxime group points away from the VX-bound serine) in the reactivator-only complex to a “semi-productive” orientation in the VX-modified complex. This observation, supported by concurrent molecular simulations, suggested that the narrow active-site gorge of hAChE may be significantly more dynamic than previously thought, allowing RS-170B to reorient inside the gorge. Furthermore, we found that small molecules can bind in the choline-binding site hindering approach to the phosphorous of VX-bound serine. Our results provide structural and mechanistic perspectives on the reactivation of OP-inhibited hAChE and demonstrate that structural studies at physiologically relevant temperatures can deliver previously overlooked insights applicable for designing next-generation antidotes.

International audience; The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Since the 1863 discovery of a new green hemoglobin derivative called “sulfhemoglobin”, the nature of the characteristic 618 nm absorption band has been the subject of several hypotheses. The experimental spectra are a function of the observation time and interplay between two major sulfheme isomer concentrations (a three- and five-membered ring adduct), with the latter being the dominant isomer at longer times. Thus, time-dependent density functional theory (TDDFT) was used to calculate the sulfheme excited states and visualize the highest occupied molecular orbitals (HOMOs) and lowest unoccupied MOs (LUMOs) of both isomers in order to interpret the transitions between them. These two isomers have distinguishable a1u and a2u HOMO energies. Formation of the three-membered ring SA isomeric structure decreases the energy of the HOMO a1u and a2u orbitals compared to the unmodified heme due to the electron-withdrawing, sulfur-containing, three-membered ring. Conversely, formation of the SC isomeric structure d...

The triterpene squalene is a key metabolic intermediate for sterols, hopanoids and various other triterpenoids. The biosynthesis of squalene is catalyzed by squalene synthase (SQS), which converts two molecules of farnesyl diphosphate to squalene. In this study, a lycophyte Selaginella moellendorffii was found to emit squalene as a volatile compound under a number of conditions that mimic biotic stresses. Searching the genome sequence of S. moellendorffii led to the identification of a putative squalene synthase gene. It was designated as SmSQS. SmSQS is homologous to known squalene synthases from other plants and animals at both the amino acid level and structural level. Recombinant SmSQS expressed in Escherichia coli catalyzed the formation of squalene using farnesyl diphosphate as substrate. The expression of SmSQS was significantly induced by the same set of stress factors that induced the emission of volatile squalene from S. moellendorffii plants. Taken together, these results support that SmSQS is responsible for the biosynthesis of volatile squalene and volatile squalene may have a role in the defense of S. moellendorffii plants against biotic stresses.

The carbapenem-hydrolyzing class D β-lactamases OXA-23 and OXA-24/40 have emerged worldwide as causative agents for β-lactam antibiotic resistance in Acinetobacter species. Many variants of these enzymes have appeared clinically, including OXA-160 and OXA-225, which both contain a P → S substitution at homologous positions in the OXA-24/40 and OXA-23 backgrounds, respectively. We purified OXA-160 and OXA-225 and used steady-state kinetic analysis to compare the substrate profiles of these variants to their parental enzymes, OXA-24/40 and OXA-23. OXA-160 and OXA-225 possess greatly enhanced hydrolytic activities against aztreonam, ceftazidime, cefotaxime, and ceftriaxone when compared to OXA-24/40 and OXA-23. These enhanced activities are the result of much lower Km values, suggesting that the P → S substitution enhances the binding affinity of these drugs. We have determined the structures of the acylated forms of OXA-160 (with ceftazidime and aztreonam) and OXA-225 (ceftazidime). These structures show that the R1 oxyimino side-chain of these drugs occupies a space near the β5-β6 loop and the omega loop of the enzymes. The P → S substitution found in OXA-160 and OXA-225 results in a deviation of the β5-β6 loop, relieving the steric clash with the R1 side-chain carboxypropyl group of aztreonam and ceftazidime. These results reveal worrying trends in the enhancement of substrate spectrum of class D β-lactamases but may also provide a map for β-lactam improvement.

The inactivation of acetylcholinesterase (AChE) by organophosphorus agent (OP) compounds is a serious problem regardless of how the individual was exposed. The reactivation of OP-inactivated AChE is dependent on the OP conjugate, and commonly a specific oxime is better at reactivating a specific OP conjugate than several diverse OP conjugates. The presented research explores the physicochemical properties needed for the reactivation of OP-inactivated AChE. Four different OPs, cyclosarin, sarin, tabun, and VX, were analyzed using the same set of oxime reactivators. A trial descriptor pool of semiempirical, traditional, and molecular interaction field descriptors was used to construct an ensemble of QSAR models for each OP-conjugate pair. Based on the molecular information and the cross-validation ability, individual QSAR models were selected to be part of an OP-conjugate consensus model. The OP-conjugate specific models provide important insight into the physicochemical properties required to reactivate the OP conjugates of interest. The reactivation of AChE inactivated with either cyclosarin or tabun requires the oxime therapeutic to possess an overall polar-positive surface area. Oxime therapeutics for the reactivation of sarin-inactivated AChE are conformationally dependent while oxime reverse therapeutics for VX require a compact region with a highly hydrophilic region and two positively charged pyridine rings.

Hybrid quantum chemical (QC)/molecular mechanical (MM) potentials are very powerful tools for molecular simulation. They are especially useful for studying processes in condensed phase systems, such as chemical reactions, that involve a relatively localized change in electronic structure and where the surrounding environment contributes to these changes but can be represented with more computationally efficient functional forms. Despite their utility, however, these potentials are not always straightforward to apply since the extent of significant electronic structure changes occurring in the condensed phase process may not be intuitively obvious. To facilitate their use we have developed an open-source graphical plug-in, GTKDynamo, that links the PyMOL visualization program and the pDynamo QC/MM simulation library. This article describes the implementation of GTKDynamo and its capabilities and illustrates its application to QC/MM simulations.

The anesthetic propofol inhibits the currents of the homo-pentameric ligand-gated ion channel GLIC, yet the crystal structure of GLIC with five propofol molecules bound symmetrically shows an open-channel conformation. To address this dilemma and determine if symmetry of propofol binding sites affects the channel conformational transition, we performed a total of 1.5 As of molecular dynamics simulations for different GLIC systems with propofol occupancies of 0, 1, 2, 3, and 5. GLIC without propofol binding or with five propofol molecules bound symmetrically showed similar channel conformation and hydration status over multiple replicates of 100-ns simulations. In contrast, asymmetric binding to one, two or three equivalent sites in different subunits accelerated the channel dehydration, which was accompanied by increased conformational heterogeneity of the pore and shifted the lateral and radial tilting angles of the pore-lining TM2 towards a closed-channel conformation. The results differentiate two groups of systems based on the propofol binding symmetry. The difference between symmetric and asymmetric groups is correlated with the variance in the propofol-binding cavity adjacent to the hydrophobic gate and the force imposed by the bound propofol. Asymmetrically bound propofol produced greater variance in the cavity size that could further elevate the conformation heterogeneity. The force trajectory generated by propofol in each subunit over the course of a simulation exhibits an ellipsoidal shape, which has the larger component tangential to the pore. Asymmetric propofol binding creates an unbalanced force that expedites the channel conformation transitions. The findings from this study not only suggest that asymmetric binding underlies the propofol functional inhibition of GLIC, but also advocate for the role of symmetry breaking in facilitating channel conformational transitions.

Widespread antibiotic resistance, particularly when mediated by broad-spectrum β-lactamases, has major implications for public health. Substitutions in the active site often allow broad-spectrum enzymes to accommodate diverse types of β-lactams. Substitutions observed outside the active site are thought to compensate for the loss of thermal stability. The OXA-1 clade of class D β-lactamases contains a pair of conserved cysteines located outside the active site that forms a disulfide bond in the periplasm. Here, the effect of the distal disulfide bond on the structure and dynamics of OXA-1 was investigated via 4 μs molecular dynamics simulations. The results reveal that the disulfide promotes the preorganized orientation of the catalytic residues and affects the conformation of the functionally important Ω loop. Furthermore, principal component analysis reveals differences in the global dynamics between the oxidized and reduced forms, especially in the motions involving the Ω loop. A dynamical network analysis indicates that, in the oxidized form, in addition to its role in ligand binding, the KTG family motif is a central hub of the global dynamics. As activity of OXA-1 has been measured only in the reduced form, we suggest that accurate assessment of its functional profile would require oxidative conditions mimicking periplasm.

Many heme-containing proteins with a histidine in the distal E7 (HisE7) position can form sulfheme in the presence of hydrogen sulfide (H2S) and a reactive oxygen species such as hydrogen peroxide. For reasons unknown, sulfheme derivatives are formed specifically on solvent-excluded heme pyrrole B. Sulfhemes severely decrease the oxygen-binding affinity in hemoglobin (Hb) and myoglobin (Mb). Here, use of hybrid quantum mechanical/molecular mechanical methods has permitted characterization of the entire process of sulfheme formation in the HisE7 mutant of hemoglobin I (HbI) from Lucina pectinata. This process includes a mechanism for H2S to enter the solvent-excluded active site through a hydrophobic channel to ultimately form a hydrogen bond with H2O2 bound to Fe(III). Proton transfer from H2O2 to His64 to form compound (Cpd) 0, followed by hydrogen transfer from H2S to the Fe(III)-H2O2 complex, results in homolytic cleavage of the O-O and S-H bonds to form a reactive thiyl radical (HS(•)), ferryl heme Cpd II, and a water molecule. Subsequently, the addition of HS(•) to Cpd II, followed by three proton transfer reactions, results in the formation of a three-membered ring ferric sulfheme that avoids migration of the radical to the protein matrix, in contrast to that in other peroxidative reactions. The transformation of this three-membered episulfide ring structure to the five-membered thiochlorin ring structure occurs through a significant potential energy barrier, although both structures are nearly isoenergetic. Both three- and five-membered ring structures reveal longer NB-Fe(III) bonds compared with other pyrrole nitrogen-Fe(III) bonds, which would lead to decreased oxygen binding. Overall, these results are in agreement with a wide range of experimental data and provide fertile ground for further investigations of sulfheme formation in other heme proteins and additional effects of H2S on cell signaling and reactivity.

Class D β-lactamases, a major source of bacterial resistance to β-lactam antibiotic therapies, represent a distinct subset of the β-lactamase superfamily. They share a serine hydrolase mechanism with Classes A/C vs. Class B. Further understanding of their sequence–structure–function relationships would benefit efforts to design a new generation of antibiotics as well as to predict evolutionary mechanisms in response to such therapies. Here we describe analyses based on our high-resolution multiple sequence alignment and phylogenetic tree of ∼80 Class D β-lactamases that leverage several 3D structures of these enzymes. We observe several sequence clusters on the phylogenetic tree, some that are species specific while others include several species from α-, β- and γ-proteobacteria. Residues characteristic of a specific cluster were identified and shown to be located just outside the active site, possibly modulating the function of the catalytic residues to facilitate reactions with specific types of β-lactams. Most significant was the discovery of a likely disulfide bond in a large group composed of α-, β- and γ-proteobacteria that would contribute to enzyme stability and hence bacterial viability under antibiotic assault. A network of co-evolving residues was identified which suggested the importance of maintaining a surface for binding a highly conserved Phe69.

Our computer simulations of gamma glutamyl semialdehyde dehydrogenase (GGSALDH, pyrroline 5-carboxylate dehydrogenase, ALDH4) were initiated from the Thermus thermophilus crystal structures in an effort to understand the effects of a seemingly subtle mutation. In humans, a natural S352L mutation gives rise to type II hyperprolinemia (mental retardation). The mutation occurs in what might be a priori considered the outer shell of the active site, affecting a residue of no obvious significance. In another member of the superfamily (ALDH3) this serine residue is an aspartate, which tethers the “distal” Lys. It has been our hypothesis that in ALDH3 this is a beneficial interaction, enabling the “proximal” Lys to interact with the carbonyl oxygen of the peptide bond with the catalytic Cys, allowing the Cys amide N to transiently protonate the tetrahedral intermediate. That the role of this Asp is significant is proved by a natural Asp-to-Asn mutation that abolishes activity. The Ser-to-Leu exchange in GGSALDH might be expected to alter the water structure at the site of mutation, and the MM simulations clearly support this. It was our hypothesis, based on initial static models of the mutation, that the leucyl side chain would block the direct or indirect interaction of the distal Lys with the active site. Our simulations indicate that this lysine residue is indeed important in explaining the molecular pathology of the mutation. Through small rotations of its C–C bonds, the Lys ɛ-amino group comes into H-bonding distance with Ser-326, the equivalent of human Ser-352. In the S326L mutant, this interaction is not possible, while the water network from this residue to the target main-chain carbonyl oxygen is disturbed as well.

A reaction path including transition states is generated for the Dowd mechanism [P. Dowd, R. Hershlne, S.W. Ham, S. Naganathan. Vitamin K and energy transduction: a base strength amplification mechanism. Science 269 (2005) 1684-1691] of action for Vitamin K carboxylase (VKC) using quantum chemical methods (B3LYP/6-311G**). VKC, an essential enzyme in mammalian systems, catalyzes the conversion of hydroquinone form of Vitamin K to the epoxide form in the presence of oxygen. An intermediate species of the oxidation of Vitamin K, an alkoxide, acts apparently to abstract the gamma hydrogen from specifically located glutamate residues. We are able to follow the Dowd proposed path to generate this alkoxide species. The geometries of the proposed model intermediates and transition states in the mechanism are energy optimized. We find that the most energetic step in the mechanism is the uni-deprotonation of the hydroquinone - once this occurs, there is only a small barrier of 3.5kcal/mol for the interaction of oxygen with the carbon to be attacked - and then the reaction proceeds downhill in free energy to form the critical alkoxide species. The results are consistent with the idea that the enzyme probably acts to facilitate the formation of the epoxide by reducing the energy required to deprotonate the hydroquinone form.

A reaction path including transition states is generated for the Silverman mechanism [R.B. Silverman, Chemical model studies for the mechanism of Vitamin K epoxide reductase, J. Am. Chem. Soc. 103 (1981) 5939–5941] of action for Vitamin K epoxide reductase (VKOR) using quantum mechanical methods (B3LYP/6-311G**). VKOR, an essential enzyme in mammalian systems, acts to convert Vitamin K epoxide, formed by Vitamin K carboxylase, to its (initial) quinone form for cellular reuse. This study elaborates on a prior work that focused on the thermodynamics of VKOR [D.W. Deerfield II, C.H. Davis, T. Wymore, D.W. Stafford, L.G. Pedersen, Int. J. Quant. Chem. 106 (2006) 2944–2952]. The geometries of proposed model intermediates and transition states in the mechanism are energy optimized. We find that once a key disulfide bond is broken, the reaction proceeds largely downhill. An important step in the conversion of the epoxide back to the quinone form involves initial protonation of the epoxide oxygen. We find that the source of this proton is likely a free mercapto group rather than a water molecule. The results are consistent with the current view that the widely used drug Warfarin likely acts by blocking binding of Vitamin K at the VKOR active site and thereby effectively blocking the initiating step. These results will be useful for designing more complete QM/MM studies of the enzymatic pathway once three-dimensional structural data is determined and available for VKOR.

Recent computer simulations of the cysteine nucleophilic attack on propanal in human mitochondrial Aldehyde Dehydrogenase (ALDH2) yielded an unexpected result; the chemically reasonable formation of a dead-end cysteine-cofactor adduct when NAD+ was in the “hydride transfer” position. More recently, this adduct found independent crystallographic support in work on formyltetrahydrofolate dehydrogenase, work which further found evidence of the same adduct on re-examination of deposited electron densities of ALDH2. Although the experimental data showed this adduct was reversible, several mechanistic questions arise from the fact that it forms at all. Here, we present results from further Quantum Mechanical/Molecular Mechanical (QM/MM) simulations toward understanding the mechanistic implications of adduct formation. These simulations revealed that formation of the oxyanion thiohemiacetal intermediate only when the nicotinamide ring of NAD+ is oriented away from the active site, contrary to prior arguments. In contrast, and in seeming paradox, when NAD is oriented to receive the hydride, disassociation of the oxyanion intermediate to form the dead-end adduct is more thermodynamically-favored than maintaining the oxyanion intermediate necessary for catalysis to proceed. However, this disassociation to the adduct could be avoided through proton transfer from the enzyme to the intermediate. Our results continue to indicate that the unlikely source of this proton is the Cys302 main chain amide.

Proton transfer is integral to many biochemical processes. However, the direct observation and structural characterization of biological proton transfer has hitherto not been possible. Here, neutron crystallography directly locates two protons before and after a pH-induced two-proton transfer between the catalytic aspartic acid residues and the hydroxyl group of the bound clinical drug, darunavir, in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by QM/MM calculations. We propose that the low-pH proton configuration in the catalytic site is critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.

Possible model, but simplistic, mechanisms for the action of vitamin K epoxide reductase (VKOR) are investigated with quantum mechanical methods (B3LYP/6-311G**). The geometries of proposed model intermediates in the mechanisms are energy optimized. Finally, the energetics of the proposed (pseudo-enzymatic) pathways are compared. We find that the several pathways are all energetically feasible. These results will be useful for designing quantum mechanical/molecular mechanical method (QM/MM) studies of the enzymatic pathway once three-dimensional structural data are determined and available for VKOR. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006

Experimental structural data on the state of substrates bound to class 3 Aldehyde Dehydrogenases (ALDH3A1) is currently unknown. We have utilized molecular mechanics (MM) simulations, in conjunction with new force field parameters for aldehydes, to study the atomic details of benzaldehyde binding to ALDH3A1. Our results indicate that while the nucleophilic Cys243 must be in the neutral state to form what are commonly called near-attack conformers (NACs), these structures do not correlate with increased complexation energy calculated with the MM-Generalized Born Molecular Volume (GBMV) method. The negatively charged Cys243 (thiolate form) of ALDH3A1 also binds benzaldehyde in a stable conformation but in this complex the sulfur of Cys243 is oriented away from benzaldehyde yet yields the most favorable MM-GBMV complexation energy. The identity of the general base, Glu209 or Glu333, in ALDHs remains uncertain. The MM simulations reveal structural and possible functional roles for both Glu209 and Glu333. Structures from the MM simulations that would support either glutamate residue as the general base were further examined with Hybrid Quantum Mechanical (QM)/MM simulations. These simulations show that, with the PM3/OPLS potential, Glu209 must go through a step-wise mechanism to activate Cys243 through an intervening water molecule while Glu333 can go through a more favorable concerted mechanism for the same activation process.

Ab initio and empirical methods were combined to optimize the partial atomic charges and Lennard-Jones parameters for two halogenated compounds, halothane (CF3CHClBr, a potent volatile anesthetic) and hexafluoroethane (CF3CF3, a nonanesthetic). Charge optimization was achieved using empirical calculations by systematically adjusting the charge assignments to fit minimum interaction energies and geometries between a TIP3 water molecule and the halogenated compounds to the corresponding values from the ab initio calculations, which were carried out at the HF/6-311+G(2d,p) and HF/6-31G(d) levels for halothane and hexafluoroethane, respectively. To optimize the Lennard-Jones parameters, the initial estimates were obtained from scaling the values from the ab initio minimum interaction energies and geometries between neon and the halogenated compounds calculated at the MP3/6-311++G(3d,3p) level. The Lennard-Jones parameters were further refined by fitting the empirical interaction energies to the corresponding ...

SummaryD-xylose isomerase (XI) is capable of sugar isomerization and slow conversion of some monosaccharides into their C2-epimers. We present X-ray and neutron crystallographic studies to locate H and D atoms during the respective isomerization and epimerization of L-arabinose to L-ribulose and L-ribose, respectively. Neutron structures in complex with cyclic and linear L-arabinose have demonstrated that the mechanism of ring-opening is the same as for the reaction with D-xylose. Structural evidence and QM/MM calculations show that in the reactive Michaelis complex L-arabinose is distorted to the high-energy 5S1 conformation; this may explain the apparent high KM for this sugar. MD-FEP simulations indicate that amino acid substitutions in a hydrophobic pocket near C5 of L-arabinose can enhance sugar binding. L-ribulose and L-ribose were found in furanose forms when bound to XI. We propose that these complexes containing Ni2+ cofactors are Michaelis-like and the isomerization between these two sugars proceeds via a cis-ene-diol mechanism.

Molecular dynamics (MD) simulation of the rat class 3 aldehyde dehydrogenase (ALDH) with nicotinamide dinucleotide (NAD) cofactors and explicit water molecules are reported. Our results demonstrate that MD simulation using the latest methodologies can maintain the crystal structure of the enzyme, as well as closely reproduce the short timescale dynamics of the enzyme. Furthermore, the examination of the distance between the nucleophilic Cys-243 and the NAD cofactor reveal important fluctuations that could be linked to ALDH catalysis. Finally, our quantum mechanical model of benzaldehyde in the active site of ALDH demonstrates that the enzyme requires only minor conformational changes to be poised for nucleophilic attack on the substrate.

A molecular dynamics (MD) simulation study of a dodecylphosphocholine (DPC) micelle in water is presented. This system contains 60 DPC molecules in 5294 water molecules. The structure, shape, hydrocarbon chain fluidity, the hydration of the head group and the hydrocarbon chain, and the dynamics of the micelle were analyzed from the 1.2 ns constant pressure MD simulation. The micelle was found to be slightly prolate, with the ratio of the moments of inertia 1:24:1.11:1. The penetration of water into the interior of the micelle is limited. The interaction of water with the micelle mainly comes from the head group, and the modes of interaction with the positively charged choline group and the negatively charged phosphate group can be deciphered from the radial distribution functions between these groups and both the hydrogen and the oxygen atoms of water. From the ratio of the trans/gauche conformers on the hydrocarbon chain, and the density distribution of the various carbons with respect to the center of the mass of the micelle, the conformational properties of the hydrocarbon chains have been analyzed. The dynamics of the micelle was probed by the trans–gauche conformational transition rates and from the time correlation function of the C–H bonds. From the latter, the order parameters and the correlation times for the internal motions have been obtained. These parameters were found to be, for the most part, in excellent agreement with those obtained from NMR relaxation.

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in an explicit sodium dodecylsulfate (SDS) micelle. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) during the equilibration stage of the simulation. It was shown that when SP-Y8 was initially placed in an insertion (perpendicular) configuration, the peptide equilibrated to a surface-bound (parallel) configuration in approximately 450 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained from the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbones of all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions come from the hydrophobic interaction between the side chains of Lys-3, Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, from Lys-3, Phe-7, Leu-10, and Met-11 in SP-Y8, and the micellar interior. Significant interactions, electrostatic and hydrogen bonding, between the N-terminal residues, Arg-Pro-Lys, and the micellar headgroups were observed. These latter interactions served to affect both the structure and, especially, the flexibility, of the N-terminus. The results from simulation of the same peptides in a water/CCl4 biphasic cell were compared with the results of the present study, and the validity of using the biphasic system as an approximation for peptide-micelle or peptide-bilayer systems is discussed.

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in a TIP3P water/CCl4 biphasic solvent system as a mimic for the water-membrane system. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) in the equilibration stage of the simulation. The starting orientation/position of the peptides for the MD simulation was either parallel to the water/CCl4 interface or in a perpendicular/insertion mode. In both cases the peptides equilibrated and adopted a near-parallel orientation within approximately 250 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained in the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbone of nearly all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions with the hydrophobic phase come from the hydrophobic interactions of the side chains of Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, and Phe-7, Leu-10, Met-11 and, to a lesser extent, Tyr-8 in SP-Y8. Concerted conformational transitions took place in the time frame of hundreds of picoseconds. The concertedness of the transition was due to the tendency of the peptide to maintain the necessary secondary structure to position the peptide properly with respect to the water/CCl4 interface.

Amorphous hydrogenated carbon (a‐C:H) films prepared by plasma polymerization can contain chemical species due not only to the process gas species, but species from prior processes in the deposition chamber. Our studies showed that standard O2 etching/ cleaning processes done prior to plasma film formation will affect future film deposition rate measurements and carbonyl group formation. Additionally, substrate etching by polymer forming gas species is shown through the 956 cm−1 negative peak in the IR spectrum.

Plant sesquiterpene synthases, a subset of the terpene synthase superfamily, are a mechanistically diverse family of enzymes capable of synthesizing hundreds of complex compounds with high regio- and stereospecificity and are of biological importance due to their role in plant defense mechanisms. In the current report we describe a large-scale, high-resolution phylogenetic analysis of ∼200 plant sesquiterpene synthases integrated with structural and experimental data that address these issues. We observe that all sequences that cluster together on the phylogenetic tree into well-defined groups share at least the first reaction in the catalytic mechanism subsequent to the initial ionization step and many share steps beyond this down to proton transfers between the enzyme and substrate. Most significant is the previously unreported high conservation of an Asp-Tyr-Asp triad. Due to its high conservation, patterns in the phylogenetic tree as well as experimental and modeling results, we suggest that this Asp-Tyr-Asp triad is an important functional element responsible for many proton transfers to and from the substrate and intermediates along the plant sesquiterpene synthase catalytic cycle and whose position can be tuned by residues outside the active site that can lead to the evolution of novel enzyme function.

Despite recent advances in fold recognition algorithms that identify template structures with distant homology to the target sequence, the quality of the target-template alignment can be a major problem for distantly related proteins in comparative modeling. Here we report for the first time on the use of ensembles of pairwise alignments obtained by stochastic backtracking as a means to improve three-dimensional comparative protein models. In every one of the 35 cases, the ensemble produced by the program probA resulted in alignments that were closer to the structural alignment than those obtained from the optimal alignment. In addition, we examined the lowest energy structure among these ensembles from four different structural assessment methods and compared these with the optimal and structural alignment model. The structural assessment methods consisted of the DFIRE, DOPE, and ProsaII statistical potential energies and the potential energy from the CHARMM protein force field coupled to a Generalized Born implicit solvent model. The results demonstrate that the generation of alignment ensembles through stochastic backtracking using probA combined with one of the statistical potentials for assessing three-dimensional structures can be used to improve comparative models.

A reaction path including transition states is generated for the Dowd mechanism [P. Dowd, R. Hershlne, S.W. Ham, S. Naganathan. Vitamin K and energy transduction: a base strength amplification mechanism. Science 269 (2005) 1684-1691] of action for Vitamin K carboxylase (VKC) using quantum chemical methods (B3LYP/6-311G**). VKC, an essential enzyme in mammalian systems, catalyzes the conversion of hydroquinone form of Vitamin K to the epoxide form in the presence of oxygen. An intermediate species of the oxidation of Vitamin K, an alkoxide, acts apparently to abstract the gamma hydrogen from specifically located glutamate residues. We are able to follow the Dowd proposed path to generate this alkoxide species. The geometries of the proposed model intermediates and transition states in the mechanism are energy optimized. We find that the most energetic step in the mechanism is the uni-deprotonation of the hydroquinone - once this occurs, there is only a small barrier of 3.5kcal/mol for the interaction of oxygen with the carbon to be attacked - and then the reaction proceeds downhill in free energy to form the critical alkoxide species. The results are consistent with the idea that the enzyme probably acts to facilitate the formation of the epoxide by reducing the energy required to deprotonate the hydroquinone form.

A reaction path including transition states is generated for the Silverman mechanism [R.B. Silverman, Chemical model studies for the mechanism of Vitamin K epoxide reductase, J. Am. Chem. Soc. 103 (1981) 5939-5941] of action for Vitamin K epoxide reductase (VKOR) using quantum mechanical methods (B3LYP/6-311G**). VKOR, an essential enzyme in mammalian systems, acts to convert Vitamin K epoxide, formed by Vitamin K carboxylase, to its (initial) quinone form for cellular reuse. This study elaborates on a prior work that focused on the thermodynamics of VKOR [D.W. Deerfield II, C.H. Davis, T. Wymore, D.W. Stafford, L.G. Pedersen, Int. J. Quant. Chem. 106 (2006) 2944-2952]. The geometries of proposed model intermediates and transition states in the mechanism are energy optimized. We find that once a key disulfide bond is broken, the reaction proceeds largely downhill. An important step in the conversion of the epoxide back to the quinone form involves initial protonation of the epoxide oxygen. We find that the source of this proton is likely a free mercapto group rather than a water molecule. The results are consistent with the current view that the widely used drug Warfarin likely acts by blocking binding of Vitamin K at the VKOR active site and thereby effectively blocking the initiating step. These results will be useful for designing more complete QM/MM studies of the enzymatic pathway once three-dimensional structural data is determined and available for VKOR.

Experimental structural data on the state of substrates bound to class 3 Aldehyde Dehydrogenases (ALDH3A1) is currently unknown. We have utilized molecular mechanics (MM) simulations, in conjunction with new force field parameters for aldehydes, to study the atomic details of benzaldehyde binding to ALDH3A1. Our results indicate that while the nucleophilic Cys243 must be in the neutral state to form what are commonly called near-attack conformers (NACs), these structures do not correlate with increased complexation energy calculated with the MM-Generalized Born Molecular Volume (GBMV) method. The negatively charged Cys243 (thiolate form) of ALDH3A1 also binds benzaldehyde in a stable conformation but in this complex the sulfur of Cys243 is oriented away from benzaldehyde yet yields the most favorable MM-GBMV complexation energy. The identity of the general base, Glu209 or Glu333, in ALDHs remains uncertain. The MM simulations reveal structural and possible functional roles for both Glu209 and Glu333. Structures from the MM simulations that would support either glutamate residue as the general base were further examined with Hybrid Quantum Mechanical (QM)/MM simulations. These simulations show that, with the PM3/OPLS potential, Glu209 must go through a step-wise mechanism to activate Cys243 through an intervening water molecule while Glu333 can go through a more favorable concerted mechanism for the same activation process.

An algorithm for detecting amino acid residues characteristic of individual protein families from within aligned collections of paralogous sequences, and its application to the ALDH3 family versus the rest of the ALDH extended family is described. Residues illuminated by this analysis include a key intramolecular tether, a lysine that makes an intersubunit contact at the dimer interface, three residues in close association with the substrate-binding funnel, and a pair of residues suggested to participate in proton relay during the catalytic cycle.

The 2013 Nobel Prize in Chemistry was awarded for the development of multiscale models for complex chemical systems, whereas the 2013 Peace Prize was given to the Organisation for the Prohibition of Chemical Weapons for their efforts to eliminate chemical warfare agents. This review relates the two by introducing the field of multiscale modeling and highlighting its application to the study of the biological mechanisms by which selected chemical weapon agents exert their effects at an atomic level.

Enzymes continue to expand their role in industry as a "green" option for the synthesis of value-added products. They are targeted for the design of drugs in pharmaceutical applications and also for protein engineering in industry to improve their efficiency, stability, and specificity. Knowledge of the exact mechanisms of enzymatic reactions may provide essential information for more effective drug design and enzyme engineering. For the first time, we are employing a joint X-ray/neutron (XN) protein crystallographic technique in combination with high-performance computing, including QM and QM/MM calculations, MD and Rosetta simulations, to investigate the mechanisms of several enzymes that are important to renewable energy and chemical synthesis. D-xylose isomerase (XI) is an enzyme which can be used to increase the production of biofuels from lignocellulosic biomass and also to synthesize rare sugars for pharmaceutical industry. XI catalyzes the reversible multi-stage sugar inter-conversion reaction facilitated by the presence of two divalent metal cations in its active site. It primarily catalyzes the isomerization of the aldo-sugar D-xylose to the keto-isomer D-xylulose, but can also epimerize L-arabinose into L-ribose, albeit much less efficiently. The reaction involves moving hydrogen atoms between the protein residues, sugar and water molecules, and can only be understood if hydrogen atoms are visualized at each reaction stage. We have obtained a number of joint XN structures of XI complexes representing snapshots along the reaction path with D-glucose, D-xylose and L-arabinose. The suggested reaction mechanism has been verified by QM calculations using the novel O(N) methodology. We are using this structural and mechanistic information to re-design XI to be more efficient on D-xylose and L-arabinose for biofuels and biomedical applications by employing QM/MM, MD, and Rosetta methodologies.

Hydrogen sulfide (H2S) in the presence of O2 or H2O2, reacts with hemoglobin(Hb) and myoglobin(Mb), to generate specifically a covalent modification of the heme pyrrole ring bearing the4 − vinyl group, to produce sulfmyoglobin and sulfhemoglobin derivatives. These species have lower O2 affinity affecting the heme protein functionality and generating sulfhemoglobinemia. Peroxidases and catalases, with distal heme pocket, similar to Hb and Mb, also produce sulfheme. The reaction has been explored, upon interaction of H2S, RSH species, with O2 and H2O2 and several hemeproteins, using UV–vis, resonance Raman, EPR spectroscopy, and theoretical calculations. The results show that the formation of the sulfheme derivative requires the presence of a: (a) proton donor from the HSR species; (b) HisE7 residue in the heme distal site with an adequate orientation form an active ternary complex; and (c) the intermediate involving the HisE7, the hydro peroxo {HemeFe (III)-OOH}, and H2S species. Homolysis of the moiety generates oxoheme ferryl compound II and HS radical. Since the rate of sulfMb formation decreases with the concentration of hydrosulfide HS− ion, ferryl compound II and the HS species are the predominant intermediates in the reaction of sulfheme formation. The HS attacks the ferryl compound II, producing the final sulfheme derivative. Under this state, H2S avoids the presence of oxoferryl heme species compound I and limits the further generation protein radicals and ROS species. The finding opens a window to further explore H2S as an antioxidant agent and unravels the arena for H2S activation by other proteins.

ACTH (1-10), an adrenocorticotropin hormone fragment, was studied by molecular dynamics (MD) simulation in the NPT ensemble in an explicit sodium dodecylsulfate (SDS) micelle. Initially, distance restraints derived from NMR nuclear Overhauser enhancements were incorporated during the equilibration stage of the simulation. The analyses of the trajectories from the subsequent unrestrained MD showed that ACTH (1-10) does not conform to a helical structure at the micelle-water interface; however, the structure is amphipathic. The loss of the helical structure is due to decreased intramolecular hydrogen bonding accompanied by an increase of hydrogen bonding between the amide hydrogens of the peptide and the micelle head-groups. ACTH (1-10) was found to lie on the surface of the SDS micelle. Most of the hydrophobic interactions came from the side-chains of Met-4, Phe-7 and Trp-9. The peptide bonds were either hydrated or involved in intramolecular hydrogen bonding. Decreased hydration for the backbone of His-6 and Phe-7 was due to intermolecular hydrogen bonding with the SDS head-groups. The time correlation functions of the N-H bonds of the peptide in water and in the micelle showed that the motions of the peptide, except for the N- and C-termini, are significantly reduced when partitioned in the micelle.

For a homo-pentameric ligand gated ion channel, agonist occupancy of only two or three of the five possible binding sites has been found to produce the maximal channel response. Does asymmetric ligand binding also apply for channel antagonism? The anesthetic propofol inhibits the homo-pentameric GLIC, but the crystal structure of GLIC bound symmetrically with propofol to each of the five binding sites shows a virtually identical open channel as the one without drug binding. We hypothesize that symmetric ligand occupancy of all five binding sites may not occur under physiological conditions, and asymmetric ligand binding facilitates channel conformational changes? To test this hypothesis, we performed multiple simulations on multiple GLIC systems with propofol occupying 0, 1, 2, 3, or 5 of the potential sites. We found that systems with symmetric ligand occupancy (0 or 5 bound propofols) showed similar channel conformation and hydration statuses. However, systems breaking the five-fold symmetry (1, 2, or 3 bound propofols) showed accelerated channel dehydration, increased conformational entropy of the pore-lining TM2 helix, and shifted tilting angles of the TM2 helix towards the closed-channel conformation. The trajectory of the force generated by propofol within each subunit is ellipsoidal with the primary force component tangential to the pore. Asymmetric ligand occupancy induced an unbalanced force around the channel, perturbed global stability, and facilitated conformational change. Our study suggests that asymmetry induced by ligand binding plays an important role in eliciting conformational changes related to protein function. Supported by NIH (R01GM066358, R01GM056257, R37GM049202, and T32GM075770) as well as the NSF via TeraGrid resources.