Search

The diastereo- and enantioselectivity obtained experimentally by List on the proline-catalized intramolecular aldol reaction of substituted 1,7-dicarbonylic compounds was accurately predicted using density functional theory methods at the B3LYP/6-31 ++G** level. A polarizable continuum model was used to describe solvent effects. The theoretical data agree in good extension with List's experimental results, both in enantioselectivity and diastereoselectivity, and allow for the confirmation of our previous rationalization of the main factors contributing to the reaction selectivity. While the enantioselectivity results from an important electrostatic contact between the forming alkoxyde group and the proline moiety, the calculated diasteroselectivity results from several steric contacts that can be established between the different substituents and from their relative orientation in respect to the ring conformation. However, for dialdehydes that can originate two diastereomeric enamine intermediates, the proline attack and the immonium formation steps can also be of major importance in the rationalization of the final reaction selectivity, as is the case in two of the six studied systems. The obtained data allows for a full rationalization of the known experimental systems as well as for the extrapolation to new ones with variable substitution at the carbonylic chain. [ABSTRACT FROM AUTHOR]

Metal ions play a central, functional, and structural role in many molecular structures, from small catalysts to metal-organic frameworks (MOFs) and proteins. Computational studies of these systems typically employ classical or quantum mechanical approaches or a combination of both. Among classical models, only the covalent metal model reproduces both geometries and charge transfer effects but requires time-consuming parameterization, especially for supramolecular systems containing repetitive units. To streamline this process, we introduce metallicious , a Python tool designed for efficient force-field parameterization of supramolecular structures. Metallicious has been tested on diverse systems including supramolecular cages, knots, and MOFs. Our benchmarks demonstrate that parameters accurately reproduce the reference properties obtained from quantum calculations and crystal structures. Molecular dynamics simulations of the generated structures consistently yield stable simulations in explicit solvent, in contrast to similar simulations performed with nonbonded and cationic dummy models. Overall, metallicious facilitates the atomistic modeling of supramolecular systems, key for understanding their dynamic properties and host-guest interactions. The tool is freely available on GitHub (https://github.com/duartegroup/metallicious).

Extended π-systems often form supramolecular aggregates, drastically changing their optical and electronic properties. However, aggregation processes can be difficult to characterize or predict. Here, we show that butadiyne-linked 8- and 12-porphyrin nanorings form stable and well-defined bimolecular aggregates with remarkably sharp NMR spectra, despite their dynamic structures and high molecular weights (12.7 to 26.0 kDa). Pyridine breaks up the aggregates into their constituent rings, which are in slow exchange with the aggregates on the NMR time scale. All the aggregates have the same general two-layer sandwich structure, as deduced from NMR spectroscopy experiments, including 1 H DOSY, 1 H- 1 H COSY, TOCSY, NOESY, and 1 H- 13 C HSQC. This structure was confirmed by analysis of residual dipolar couplings from 13 C-coupled 1 H- 13 C HSQC experiments on one of the 12-ring aggregates. Variable-temperature NMR spectroscopy revealed an internal ring-on-ring rotation process by which two π-π stacked conformers interconvert via a staggered conformation. A slower dynamic process, involving rotation of individual porphyrin units, was also detected by exchange spectroscopy in the 8-ring aggregates, implying partial disaggregation and reassociation. Molecular dynamics simulations indicate that the 8-ring aggregates are bowl-shaped and highly fluxional, compared to the 12-ring aggregates, which are cylindrical. This work demonstrates that large π-systems can form surprisingly well-defined aggregates and may inspire the design of other noncovalent assemblies.

Organic(porous) and metal-organic cages are promising biomimetic platforms with diverse applications spanning recognition, sensing, and catalysis. The key to the emergence of these functions is the presence of well-defined inner cavities capable of binding a wide range of guest molecules and modulating their properties. However, despite the myriad cage architectures currently available, the rational design of structurally diverse and functional cages with specific host-guest properties remains challenging. Efficiently predicting such properties is critical for accelerating the discovery of novel functional cages. Herein, we introduce CageCavityCalc ( C 3), a Python-based tool for calculating the cavity size of molecular cages. The code is available on GitHub at https://github.com/VicenteMartiCentelles/CageCavityCalc. C3 utilizes a novel algorithm that enables the rapid calculation of cavity sizes for a wide range of molecular structures and porous systems. Moreover, C 3 facilitates easy visualization of the computed cavity size alongside hydrophobic and electrostatic potentials, providing insights into host-guest interactions within the cage. Furthermore, the calculated cavity can be visualized using widely available visualization software, such as PyMol, VMD, or ChimeraX. To enhance user accessibility, a PyMol plugin has been created, allowing nonspecialists to use this tool without requiring computer programming expertise. We anticipate that the deployment of this computational tool will significantly streamline cage cavity calculations, thereby accelerating the discovery of functional cages.

The release of strain energy is a fundamental driving force for organic reactions. However, absolute strain energy alone is an insufficient predictor of reactivity, evidenced by the similar ring strain but disparate reactivity of cyclopropanes and cyclobutanes. In this work, we demonstrate that electronic delocalization is a key factor that operates alongside strain release to boost, or even dominate, reactivity. This delocalization principle extends across a wide range of molecules containing three-membered rings such as epoxides, aziridines, and propellanes and also applies to strain-driven cycloaddition reactions. Our findings lead to a "rule of thumb" for the accurate prediction of activation barriers in such systems, which can be easily applied to reactions involving many of the strained building blocks commonly encountered in organic synthesis, medicinal chemistry, polymer science, and bioconjugation. Given the significance of electronic delocalization in organic chemistry, for example in aromatic π-systems and hyperconjugation, we anticipate that this concept will serve as a versatile tool to understand and predict organic reactivity.

Cage catalysis continues to create significant interest, yet catalyst function remains poorly understood. Herein, we report mechanistic insights into coordination-cage-catalyzed Michael addition using kinetic and computational methods. The study has been enabled by the detection of identifiable catalyst intermediates, which allow the evolution of different cage species to be monitored and modeled alongside reactants and products. The investigations show that the overall acceleration results from two distinct effects. First, the cage reaction shows a thousand-fold increase in the rate constant for the turnover-limiting C-C bond-forming step compared to a reference state. Computational modeling and experimental analysis of activation parameters indicate that this stems from a significant reduction in entropy, suggesting substrate coencapsulation. Second, the cage markedly acidifies the bound pronucleophile, shifting this equilibrium by up to 6 orders of magnitude. The combination of these two factors results in accelerations up to 10 9 relative to bulk-phase reference reactions. We also show that the catalyst can fundamentally alter the reaction mechanism, leading to intermediates and products that are not observable outside of the cage. Collectively, the results show that cage catalysis can proceed with very high activity and unique selectivity by harnessing a series of individually weak noncovalent interactions.

Aldehyde dehydrogenase enzymes (ALDHs) are widely studied for their roles in disease propagation and cell metabolism. Their use in biocatalysis applications, for the conversion of aldehydes to carboxylic acids, has also been recognized. Understanding the structural features and functions of both prokaryotic and eukaryotic ALDHs is key to uncovering novel applications of the enzyme and probing its role in disease propagation. The thermostable enzyme ALDH Tt originating from Thermus thermophilus , strain HB27, possesses a unique extension of its C-terminus, which has been evolutionarily excluded from mesophilic counterparts and other thermophilic enzymes in the same genus. In this work, the thermophilic adaptation is studied by the expression and optimized purification of mutant ALDH Tt- 508, with a 22-amino acid truncation of the C-terminus. The mutant shows increased activity throughout production compared to native ALDH Tt , indicating an opening of the active site upon C-terminus truncation and giving rationale into the evolutionary exclusion of the C-terminal extension from similar thermophilic and mesophilic ALDH proteins. Additionally, the C-terminus is shown to play a role in controlling substrate specificity of native ALDH, particularly in excluding catalysis of certain large and certain aromatic ortho-substituted aldehydes, as well as modulating the protein's pH tolerance by increasing surface charge. Dynamic light scattering and size-exclusion HPLC methods are used to show the role of the C-terminus in ALDH Tt oligomeric stability at the cost of catalytic efficiency. Studying the aggregation rate of ALDH Tt with and without a C-terminal extension leads to the conclusion that ALDH Tt follows a monomolecular reaction aggregation mechanism.

While human mobility plays a crucial role in determining ambient air pollution exposures and health risks, research to date has assessed risks on the basis of almost solely residential location. Here, we leveraged a database of ∼128-144 million workers in the United States and published ambient PM 2.5 data between 2011 and 2018 to explore how incorporating information on both workplace and residential location changes our understanding of disparities in air pollution exposure. In general, we observed higher workplace exposures relative to home exposures, as well as increased exposures for nonwhite and less educated workers relative to the national average. Workplace exposure disparities were higher among racial and ethnic groups and job types than by income, education, age, and sex. Not considering workplace exposures can lead to systematic underestimations in disparities in exposure among these subpopulations. We also quantified the error in assigning workers home instead of a weighted home-and-work exposure. We observed that biases in associations between PM 2.5 and health impacts by using home instead of home-and-work exposure were the highest among urban, younger populations.

Target validation remains a challenge in drug discovery, which leads to a high attrition rate in the drug discovery process, particularly in Phase II clinical trials. Consequently, new approaches to enhance target validation are valuable tools to improve the drug discovery process. Here, we report the combination of site-directed mutagenesis and electrophilic fragments to enable the rapid identification of small molecules that selectively inhibit the mutant protein. Using the bromodomain-containing protein BRD4 as an example, we employed a structure-based approach to identify the L94C mutation in the first bromodomain of BRD4 [BRD4(1)] as having a minimal effect on BRD4(1) function. We then screened a focused, KAc mimic-containing fragment set and a diverse fragment library against the mutant and wild-type proteins and identified a series of fragments that showed high selectivity for the mutant protein. These compounds were elaborated to include an alkyne click tag to enable the attachment of a fluorescent dye. These clickable compounds were then assessed in HEK293T cells, transiently expressing BRD4(1) WT or BRD4(1) L94C , to determine their selectivity for BRD4(1) L94C over other possible cellular targets. One compound was identified that shows very high selectivity for BRD4(1) L94C over all other proteins. This work provides a proof-of-concept that the combination of site-directed mutagenesis and electrophilic fragments, in a mutate and conjugate approach, can enable rapid identification of small molecule inhibitors for an appropriately mutated protein of interest. This technology can be used to assess the cellular phenotype of inhibiting the protein of interest, and the electrophilic ligand provides a starting point for noncovalent ligand development.

Posttranslational modifications of proteins can impact their therapeutic efficacy, stability, and potential for pharmaceutical development. The Group A Streptococcus pyogenes C5a peptidase (ScpA) is a multi-domain protein composed of an N-terminal signal peptide, a catalytic domain (including propeptide), three fibronectin domains, and cell membrane-associated domains. It is one of several proteins produced by Group A S. pyogenes known to cleave components of the human complement system. After signal peptide removal, ScpA undergoes autoproteolysis and cleaves its propeptide for full maturation. The exact location and mechanism of the propeptide cleavage, and the impact of this cleavage on stability and activity, are not clearly understood, and the exact primary sequence of the final enzyme is not known. A form of ScpA with no autoproteolysis fragments of propeptide present may be more desirable for pharmaceutical development from a regulatory and a biocompatibility in the body perspective. The current study describes an in-depth structural and functional characterization of propeptide truncated variants of ScpA expressed in Escherichia coli cells. All three purified ScpA variants, ScpA, 79ΔPro, and 92ΔPro, starting with N32, D79, and A92 positions, respectively, showed similar activity against C5a, which suggests a propeptide-independent activity profile of ScpA. CE-SDS and MALDI top-down sequencing analyses highlight a time-dependent propeptide autoproteolysis of ScpA at 37 °C with a distinct end point at A92 and/or D93. In comparison, all three variants of ScpA exhibit similar stability, melting temperatures, and secondary structure orientation. In summary, this work not only highlights propeptide localization but also provides a strategy to recombinantly produce a final mature and active form of ScpA without any propeptide-related fragments.

Mobile ambient air quality monitoring is rapidly changing the current paradigm of air quality monitoring and growing as an important tool to address air quality and climate data gaps across the globe. This review seeks to provide a systematic understanding of the current landscape of advances and applications in this field. We observe a rapidly growing number of air quality studies employing mobile monitoring, with low-cost sensor usage drastically increasing in recent years. A prominent research gap was revealed, highlighting the double burden of severe air pollution and poor air quality monitoring in low- and middle-income regions. Experiment-design-wise, the advances in low-cost monitoring technology show great potential in bridging this gap while bringing unique opportunities for real-time personal exposure, large-scale deployment, and diversified monitoring strategies. The median value of unique observations at the same location in spatial regression studies is ten, which can be used as a rule-of-thumb for future experiment design. Data-analysis-wise, even though data mining techniques have been extensively employed in air quality analysis and modeling, future research can benefit from exploring air quality information from nontabular data, such as images and natural language.

We applied a bottom-up approach to develop biofunctional supramolecular hydrogels from an aromatic glycodipeptide. The self-assembly of the glycopeptide was induced by either temperature manipulation (heating-cooling cycle) or solvent (DMSO to water) switch. The sol-gel transition was salt-triggered in cell culture media and resulted in gels with the same chemical compositions but different mechanical properties. Human adipose derived stem cells (hASCs) cultured on these gels under basal conditions (i.e., without differentiation factors) overexpressed neural markers, such as GFAP, Nestin, MAP2, and βIII-tubulin, confirming the differentiation into neural lineages. The mechanical properties of the gels influenced the number and distribution of the adhered cells. A comparison with gels obtained from the nonglycosylated peptide showed that glycosylation is crucial for the biofunctionality of the hydrogels by capturing and preserving essential growth factors, e.g., FGF-2.

The SARS-CoV-2 main protease (M pro ) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. M pro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how M pro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of M pro to the presence and absence of a substrate. The results highlight communication between the M pro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of M pro . More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

Four-membered heterocycles offer exciting potential as small polar motifs in medicinal chemistry but require further methods for incorporation. Photoredox catalysis is a powerful method for the mild generation of alkyl radicals for C-C bond formation. The effect of ring strain on radical reactivity is not well understood, with no studies that address this question systematically. Examples of reactions that involve benzylic radicals are rare, and their reactivity is challenging to harness. This work develops a radical functionalization of benzylic oxetanes and azetidines using visible light photoredox catalysis to prepare 3-aryl-3-alkyl substituted derivatives and assesses the influence of ring strain and heterosubstitution on the reactivity of small-ring radicals. 3-Aryl-3-carboxylic acid oxetanes and azetidines are suitable precursors to tertiary benzylic oxetane/azetidine radicals which undergo conjugate addition into activated alkenes. We compare the reactivity of oxetane radicals to other benzylic systems. Computational studies indicate that Giese additions of unstrained benzylic radicals into acrylates are reversible and result in low yields and radical dimerization. Benzylic radicals as part of a strained ring, however, are less stable and more π-delocalized, decreasing dimer and increasing Giese product formation. Oxetanes show high product yields due to ring strain and Bent's rule rendering the Giese addition irreversible.

TRIM33 is a member of the tripartite motif (TRIM) family of proteins, some of which possess E3 ligase activity and are involved in the ubiquitin-dependent degradation of proteins. Four of the TRIM family proteins, TRIM24 (TIF1α), TRIM28 (TIF1β), TRIM33 (TIF1γ) and TRIM66, contain C-terminal plant homeodomain (PHD) and bromodomain (BRD) modules, which bind to methylated lysine (KMe n ) and acetylated lysine (KAc), respectively. Here we investigate the differences between the two isoforms of TRIM33, TRIM33α and TRIM33β, using structural and biophysical approaches. We show that the N1039 residue, which is equivalent to N140 in BRD4(1) and which is conserved in most BRDs, has a different orientation in each isoform. In TRIM33β, this residue coordinates KAc, but this is not the case in TRIM33α. Despite these differences, both isoforms show similar affinities for H3 1-27 K18Ac, and bind preferentially to H3 1-27 K9Me 3 K18Ac. We used this information to develop an AlphaScreen assay, with which we have identified four new ligands for the TRIM33 PHD-BRD cassette. These findings provide fundamental new information regarding which histone marks are recognized by both isoforms of TRIM33 and suggest starting points for the development of chemical probes to investigate the cellular function of TRIM33.

Thiophene S,S- dioxides are underutilized tools for the de novo construction of benzene rings in organic synthesis. We report a collective synthesis of nine illudalane sesquiterpenes using bicyclic thiophene S,S- dioxides as generalized precursors to the indane core of the natural products. Exploiting furans as unusual dienophiles in this inverse electron demand Diels-Alder cascade, this concise and convergent approach enables the synthesis of these targets in as little as five steps. Theoretical studies rationalize the reactivity of thiophene S,S- dioxides with both electron-poor and electron-rich dienophiles and reveal reaction pathways involving either nonpolar pericyclic or bifurcating ambimodal cycloadditions. Overall, this work demonstrates the wider potential of thiophene S,S -dioxides as convenient and flexible precursors to polysubstituted arenes.

Predictive models based on mobile measurements have been increasingly used to understand the spatiotemporal variations of intraurban air quality. However, the effects of meteorological factors, which significantly affect the dispersion of air pollution, on the urban-form-air-quality relationship have not been understood on a granular level. We attempt to fill this gap by developing predictive models of particulate matter (PM) in the Bronx (New York City) using meteorological and urban form parameters. The granular PM data was collected by mobile low-cost sensors as the ground truth. To evaluate the effects of meteorological factors, we compared the performance of models using the urban form within fixed and wind-sensitive buffers, respectively. We find better predictive power in the wind-sensitive group ( R = 0.85) for NC 10 (number concentration for particles with diameters of 1 μm-10 μm) than the control group ( R = 0.01), and modest improvements for PM 2.5 ( R = 0.84 for the wind sensitive group, R = 0.77 for the control group), indicating that incorporating meteorological factors improved the predictive power of our models. We also found that urban form factors account for 62.95% of feature importance for NC 10 and 14.90% for PM 2.5 (9.99% and 4.91% for 3-D and 2-D urban form factors, respectively) in our Random Forest models. It suggests the importance of incorporating urban form factors, especially for the uncommonly used 3-D characteristics, in estimating intraurban PM. Our method can be applied in other cities to better capture the influence of urban context on PM levels.

Self-assembled metallo-organic cages have emerged as promising biomimetic platforms that can encapsulate whole substrates akin to an enzyme active site. Extensive experimental work has enabled access to a variety of structures, with a few notable examples showing catalytic behavior. However, computational investigations of metallo-organic cages are scarce, not least due to the challenges associated with their modeling and the lack of accurate and efficient protocols to evaluate these systems. In this review, we discuss key molecular principles governing the design of functional metallo-organic cages, from the assembly of building blocks through binding and catalysis. For each of these processes, computational protocols will be reviewed, considering their inherent strengths and weaknesses. We will demonstrate that while each approach may have its own specific pitfalls, they can be a powerful tool for rationalizing experimental observables and to guide synthetic efforts. To illustrate this point, we present several examples where modeling has helped to elucidate fundamental principles behind molecular recognition and reactivity. We highlight the importance of combining computational and experimental efforts to speed up supramolecular catalyst design while reducing time and resources., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

Carbon-carbon bond-forming processes that involve the deprotonation of a weakly acidic C-H pro-nucleophile using a strong Brønsted base are central to synthetic methodology. Enzymes also catalyze C-C bond formation from weakly C-H acidic substrates; however, they accomplish this at pH 7 using only collections of noncovalent interactions. Here, we show that a simple, bioinspired synthetic cage catalyzes Michael addition reactions using only Coulombic and other weak interactions to activate various pro-nucleophiles and electrophiles. The anion-stabilizing property of the cage promotes spontaneous pro-nucleophile deprotonation, suggesting acidity enhancement equivalent to several p K a units. Using a second noncovalent reagent-commercially available 18-crown-6-facilitates catalytic base-free addition of several challenging Michael partners. The cage's microenvironment also promotes high diastereoselectivity compared to a conventional base-catalyzed reaction.

Sulfoximines and sulfonimidamides are promising compounds for medicinal and agrochemistry. As monoaza analogues of sulfones and sulfonamides, respectively, they combine good physicochemical properties, high stability, and the ability to build complexity from a three-dimensional core. However, a lack of quick and efficient methods to prepare these compounds has hindered their uptake in molecule discovery programmes. Herein, we describe a unified, one-pot approach to both sulfoximines and sulfonimidamides, which exploits the high electrophilicity of sulfinyl nitrenes. We generate these rare reactive intermediates from a novel sulfinylhydroxylamine (R-O-N=S=O) reagent through an N-O bond fragmentation process. Combining sulfinyl nitrenes with carbon and nitrogen nucleophiles enables the synthesis of sulfoximines and sulfonimidamides in a reaction time of just 15 min. Alkyl, (hetero)aryl, and alkenyl organometallic reagents can all be used as the first or second component in the reaction, while primary and secondary amines, and anilines, all react with high efficiency as the second nucleophile. The tolerance of the reaction to steric and electronic factors has allowed for the synthesis of the most diverse set of sulfoximines and sulfonimidamides yet described. Experimental and computational investigations support the intermediacy of sulfinyl nitrenes, with nitrene formation proceeding via a transient triplet intermediate before reaching a planar singlet species.

Metallocages offer a diverse and underexplored region of chemical space in which to search for novel catalysts and substrate hosts. However, the ability to tailor such structures toward applications in binding and catalysis is a challenging task. Here, we present an open-source computational toolkit, cgbind , that facilitates the construction, characterization, and prediction of functional metallocages. It employs known structural scaffolds as starting points and computationally efficient approaches for property evaluation. We demonstrate the ability of cgbind to construct libraries of cages with varied topologies and linker functionalities, generate accurate geometries (RMSD < 1.5 Å to crystal structures), and predict substrate binding with accuracy on par with semiempirical QM, all in seconds. The cgbind code presented here is freely available at github.com/duartegroup/cgbind and also via a web-based graphical user interface at cgbind.chem.ox.ac.uk. The protocol described here paves the way for high-throughput virtual screening of potential supramolecular structures, accelerating the search for new hosts and catalysts.

Modifying the reactivity of substrates by encapsulation is a fundamental principle of capsule catalysis. Here we show an alternative strategy, wherein catalytic activation of otherwise inactive quinone "co-factors" by a simple Pd 2 L 4 capsule promotes a range of bulk-phase, radical-cation cycloadditions. Solution electron-transfer experiments and cyclic voltammetry show that the cage anodically shifts the redox potential of the encapsulated quinone by a significant 1 V. Moreover, the capsule also protects the reduced semiquinone from protonation, thus transforming the role of quinones from stoichiometric oxidants into catalytic single-electron acceptors. We envisage that the host-guest-induced release of an "electron hole" will translate to various forms of non-encapsulated catalysis that involve other difficult-to-handle, highly reactive species.

Self-assembled cages have emerged as novel platforms to explore bioinspired catalysis. While many different size and shape supramolecular structures are now readily accessible, only a few are known to accelerate chemical reactions under substoichiometric conditions. These limited examples point to a poor understanding of cage catalysis in general, limiting the ability to design new systems. Here we show that a simple and efficient density-functional-theory-based methodology, informed by explicitly solvated molecular dynamics and coupled cluster calculations, is sufficient to accurately reproduce experimental guest binding affinities (MAD = 1.9 kcal mol -1 ) and identify the catalytic Diels-Alder proficiencies (>80% accuracy) of two homologous Pd 2 L 4 metallocages with a variety of substrates. This analysis reveals how subtle structural differences in the cage framework affect binding and catalysis. These effects manifest in a smaller distortion and more favorable interaction energy for the catalytic cage compared to the inactive structure. This study gives detailed insight that would otherwise be difficult to obtain from experiments, providing new opportunities in the design of catalytically active supramolecular cages.

Transcriptional repressor EthR from Mycobacterium tuberculosis is a valuable target for antibiotic booster drugs. We previously reported a virtual screening campaign to identify EthR inhibitors for development. Two ligand binding orientations were often proposed, though only the top scoring pose was utilized for filtering of the large data set. We obtained biophysically validated hits, some of which yielded complex crystal structures. In some cases, the crystallized binding mode and top scoring mode agree, while for others an alternate ligand binding orientation was found. In this contribution, we combine rigid docking, molecular dynamics simulations, and the linear interaction energy method to calculate binding free energies and derive relative binding energies for a number of EthR inhibitors in both modes. This strategy allowed us to correctly predict the most favorable orientation. Therefore, this widely applicable approach will be suitable to triage multiple binding modes within EthR and other potential drug targets with similar characteristics.

The base-catalyzed rearrangement of arylindenols is a rare example of a suprafacial [1,3]-hydrogen atom transfer. The mechanism has been proposed to proceed via sequential [1,5]-sigmatropic shifts, which occur in a selective sense and avoid an achiral intermediate. A computational analysis using quantum chemistry casts serious doubt on these suggestions: These pathways have enormous activation barriers, and in constrast to what is observed experimentally, they overwhelmingly favor a racemic product. Instead we propose that a suprafacial [1,3]-prototopic shift occurs in a two-step deprotonation/reprotonation sequence. This mechanism is favored by 15 kcal mol -1 over that previously proposed. Most importantly, this is also consistent with stereospecificity since reprotonation occurs rapidly on the same π-face. We have used explicitly solvated molecular dynamics studies to study the persistence and condensed-phase dynamics of the intermediate ion-pair formed in this reaction. Chirality transfer is the result of a particularly resilient contact ion-pair, held together by electrostatic attraction and a critical NH···π interaction which ensures that this species has an appreciable lifetime even in polar solvents such as DMSO and MeOH.

We describe the first theoretical study of a landmark example of chiral anion phase-transfer catalysis: the enantioselective ring-opening of meso-aziridinium and episulfonium cations promoted by asymmetric counteranion-directed catalysis (ACDC). The mechanism of ion-pairing, ring-opening, and catalyst deactivation have been studied in the condensed phase with both classical and quantum methods using explicitly and implicitly solvated models. We find that the stability of chiral ion-pairs, a prerequisite for asymmetric catalysis, is dominated by electrostatic interactions at long range and by CH···O interactions at short range. The decisive role of solvent upon ion-pair formation and of nonbonding interactions upon enantioselectivity are quantified by complementary computational approaches. The major enantiomer is favored by a smaller distortion of the substrate, demonstrated by a distortion/interaction analysis. Our computational results rationalize the stereoselectivity for several experimental results and demonstrate a combined classical/quantum approach to perform realistic modeling of chiral counterion catalysis in solution.

Despite the numerous experimental and theoretical studies on phosphate monoester hydrolysis, significant questions remain concerning the mechanistic details of these biologically critical reactions. In the present work we construct a linear free energy relationship for phosphate monoester hydrolysis to explore the effect of modulating leaving group pKa on the competition between solvent- and substrate-assisted pathways for the hydrolysis of these compounds. Through detailed comparative electronic-structure studies of methyl phosphate and a series of substituted aryl phosphate monoesters, we demonstrate that the preferred mechanism is dependent on the nature of the leaving group. For good leaving groups, a strong preference is observed for a more dissociative solvent-assisted pathway. However, the energy difference between the two pathways gradually reduces as the leaving group pKa increases and creates mechanistic ambiguity for reactions involving relatively poor alkoxy leaving groups. Our calculations show that the transition-state structures vary smoothly across the range of pKas studied and that the pathways remain discrete mechanistic alternatives. Therefore, while not impossible, a biological catalyst would have to surmount a significantly higher activation barrier to facilitate a substrate-assisted pathway than for the solvent-assisted pathway when phosphate is bonded to good leaving groups. For poor leaving groups, this intrinsic preference disappears.

The enzymes that facilitate phosphate and sulfate hydrolysis are among the most proficient natural catalysts known to date. Interestingly, a large number of these enzymes are promiscuous catalysts that exhibit both phosphatase and sulfatase activities in the same active site and, on top of that, have also been demonstrated to efficiently catalyze the hydrolysis of other additional substrates with varying degrees of efficiency. Understanding the factors that underlie such multifunctionality is crucial both for understanding functional evolution in enzyme superfamilies and for the development of artificial enzymes. In this Current Topic, we have primarily focused on the structural and mechanistic basis for catalytic promiscuity among enzymes that facilitate both phosphoryl and sulfuryl transfer in the same active site, while comparing this to how catalytic promiscuity manifests in other promiscuous phosphatases. We have also drawn on the large number of experimental and computational studies of selected model systems in the literature to explore the different features driving the catalytic promiscuity of such enzymes. Finally, on the basis of this comparative analysis, we probe the plausible origins and determinants of catalytic promiscuity in enzymes that catalyze phosphoryl and sulfuryl transfer.

Potato epoxide hydrolase 1 exhibits rich enantio- and regioselectivity in the hydrolysis of a broad range of substrates. The enzyme can be engineered to increase the yield of optically pure products as a result of changes in both enantio- and regioselectivity. It is thus highly attractive in biocatalysis, particularly for the generation of enantiopure fine chemicals and pharmaceuticals. The present work aims to establish the principles underlying the activity and selectivity of the enzyme through a combined computational, structural, and kinetic study using the substrate trans -stilbene oxide as a model system. Extensive empirical valence bond simulations have been performed on the wild-type enzyme together with several experimentally characterized mutants. We are able to computationally reproduce the differences between the activities of different stereoisomers of the substrate and the effects of mutations of several active-site residues. In addition, our results indicate the involvement of a previously neglected residue, H104, which is electrostatically linked to the general base H300. We find that this residue, which is highly conserved in epoxide hydrolases and related hydrolytic enzymes, needs to be in its protonated form in order to provide charge balance in an otherwise negatively charged active site. Our data show that unless the active-site charge balance is correctly treated in simulations, it is not possible to generate a physically meaningful model for the enzyme that can accurately reproduce activity and selectivity trends. We also expand our understanding of other catalytic residues, demonstrating in particular the role of a noncanonical residue, E35, as a "backup base" in the absence of H300. Our results provide a detailed view of the main factors driving catalysis and regioselectivity in this enzyme and identify targets for subsequent enzyme design efforts.

It is becoming widely accepted that catalytic promiscuity, i.e., the ability of a single enzyme to catalyze the turnover of multiple, chemically distinct substrates, plays a key role in the evolution of new enzyme functions. In this context, the members of the alkaline phosphatase superfamily have been extensively studied as model systems in order to understand the phenomenon of enzyme multifunctionality. In the present work, we model the selectivity of two multiply promiscuous members of this superfamily, namely the phosphonate monoester hydrolases from Burkholderia caryophylli and Rhizobium leguminosarum. We have performed extensive simulations of the enzymatic reaction of both wild-type enzymes and several experimentally characterized mutants. Our computational models are in agreement with key experimental observables, such as the observed activities of the wild-type enzymes, qualitative interpretations of experimental pH-rate profiles, and activity trends among several active site mutants. In all cases the substrates of interest bind to the enzyme in similar conformations, with largely unperturbed transition states from their corresponding analogues in aqueous solution. Examination of transition-state geometries and the contribution of individual residues to the calculated activation barriers suggest that the broad promiscuity of these enzymes arises from cooperative electrostatic interactions in the active site, allowing each enzyme to adapt to the electrostatic needs of different substrates. By comparing the structural and electrostatic features of several alkaline phosphatases, we suggest that this phenomenon is a generalized feature driving selectivity and promiscuity within this superfamily and can be in turn used for artificial enzyme design.

Understanding phosphoryl and sulfuryl transfer is central to many biochemical processes. However, despite decades of experimental and computational studies, a consensus concerning the precise mechanistic details of these reactions has yet to be reached. In this work we perform a detailed comparative theoretical study of the hydrolysis of p-nitrophenyl phosphate, methyl phosphate and p-nitrophenyl sulfate, all of which have served as key model systems for understanding phosphoryl and sulfuryl transfer reactions, respectively. We demonstrate the existence of energetically similar but mechanistically distinct possibilities for phosphate monoester hydrolysis. The calculated kinetic isotope effects for p-nitrophenyl phosphate provide a means to discriminate between substrate- and solvent-assisted pathways of phosphate monoester hydrolysis, and show that the solvent-assisted pathway dominates in solution. This preferred mechanism for p-nitrophenyl phosphate hydrolysis is difficult to find computationally due to the limitations of compressing multiple bonding changes onto a 2-dimensional energy surface. This problem is compounded by the need to include implicit solvation to at least microsolvate the system and stabilize the highly charged species. In contrast, methyl phosphate hydrolysis shows a preference for a substrate-assisted mechanism. For p-nitrophenyl sulfate hydrolysis there is only one viable reaction pathway, which is similar to the solvent-assisted pathway for phosphate hydrolysis, and the substrate-assisted pathway is not accessible. Overall, our results provide a unifying mechanistic framework that is consistent with the experimentally measured kinetic isotope effects and reconciles the discrepancies between theoretical and experimental models for these biochemically ubiquitous classes of reaction.

The cationic dummy atom approach provides a powerful nonbonded description for a range of alkaline-earth and transition-metal centers, capturing both structural and electrostatic effects. In this work we refine existing literature parameters for octahedrally coordinated Mn(2+), Zn(2+), Mg(2+), and Ca(2+), as well as providing new parameters for Ni(2+), Co(2+), and Fe(2+). In all the cases, we are able to reproduce both M(2+)-O distances and experimental solvation free energies, which has not been achieved to date for transition metals using any other model. The parameters have also been tested using two different water models and show consistent performance. Therefore, our parameters are easily transferable to any force field that describes nonbonded interactions using Coulomb and Lennard-Jones potentials. Finally, we demonstrate the stability of our parameters in both the human and Escherichia coli variants of the enzyme glyoxalase I as showcase systems, as both enzymes are active with a range of transition metals. The parameters presented in this work provide a valuable resource for the molecular simulation community, as they extend the range of metal ions that can be studied using classical approaches, while also providing a starting point for subsequent parametrization of new metal centers.

Sulfonate ester hydrolysis has been the subject of recent debate, with experimental evidence interpreted in terms of both stepwise and concerted mechanisms. In particular, a recent study of the alkaline hydrolysis of a series of benzene arylsulfonates (Babtie et al., Org. Biomol. Chem. 10, 2012, 8095) presented a nonlinear Brønsted plot, which was explained in terms of a change from a stepwise mechanism involving a pentavalent intermediate for poorer leaving groups to a fully concerted mechanism for good leaving groups and supported by a theoretical study. In the present work, we have performed a detailed computational study of the hydrolysis of these compounds and find no computational evidence for a thermodynamically stable intermediate for any of these compounds. Additionally, we have extended the experimental data to include pyridine-3-yl benzene sulfonate and its N-oxide and N-methylpyridinium derivatives. Inclusion of these compounds converts the Brønsted plot to a moderately scattered but linear correlation and gives a very good Hammett correlation. These data suggest a concerted pathway for this reaction that proceeds via an early transition state with little bond cleavage to the leaving group, highlighting the care that needs to be taken with the interpretation of experimental and especially theoretical data.

The base-catalyzed dehydration of benzene cis-1,2-dihydrodiols is driven by formation of an aromatic product as well as intermediates potentially stabilized by hyperaromaticity. Experiments exhibit surprising shifts in isotope effects, indicating an unusual mechanistic balance on the E2-E1cB continuum. In this study, both 1- and 2-dimensional free energy surfaces are generated for these compounds with various substituents, using density functional theory and a mixed implicit/explicit solvation model. The computational data help unravel hidden intermediates along the reaction coordinate and provide a novel conceptual framework for distinguishing between competing pathways in this and any other system with borderline reaction mechanisms.

Pin1 represents an enzyme that specifically catalyzes the isomerization of peptide bonds between phosphorylated threonine or serine residues and proline. Despite its relevance as molecular timer in a number of biological processes related to cancer and Alzheimer disease, a detailed understanding of the factors contributing to the catalysis is still missing. In this study, we employ extensive QM/MM molecular dynamics simulations in combination with the mean reaction force (MRF) to discern the influence of the enzyme on the reaction mechanism and the origin of the catalysis. As a recently introduced method, the MRF separates the activation free energy barrier to reach the transition state into structural and electronic contributions providing a more detailed description of the enzyme's function. As a reference, we first study the isomerization starting from the cis form in solution and obtain a free energy barrier and a reaction free energy, which are in agreement with previous studies and experiment. With the new mean reaction force method, intramolecular hydrogen bonds in the peptide were identified that stabilize the transition state and reduce the electronic contribution to the free energy barrier. To elucidate the mechanism of catalysis of Pin1, the reaction in solution and in the catalytic cavity of the enzyme were compared. Both yield the same free energy barrier for the isomerization of the cis form, but with different decomposition in structural and electronic contributions by the mean reaction force. The enzyme reduces the energy required for structural rearrangements to reach the transition state, pointing to a destabilization of the reactant, but increases the electronic contribution to the barrier through specific enzyme-peptide hydrogen bonds. In the reverse reaction, the isomerization of the trans form, the enzyme alters the energetics and the mechanism of the reaction considerably. Unfavorable enzyme-peptide interactions in the catalytic cavity during the isomerization change the reaction coordinate, resulting in two minima with small energy differences to the transition state. These small free energy barriers should in principle make the reaction feasible at room temperature once the conformer is bound in the right conformation.

We have computationally investigated the mechanism of the H(2) activation reaction by (amino)carbenes compounds. Describing the electronic activity taking place during the reaction through the Reaction Electronic Flux, it has been possible to elucidate the mechanism of the hydrogen activation process and assign the energetic cost associated to every chemical event that drives the process along the reaction coordinate; this is crucial information to rationalize the reported experimental results. It has been observed that the substituent effect may induce early charge-transfer phenomena that increases the energy barrier and lowers the exothermicity of the reaction. Reversibility of the process is discussed in light of specific interactions defining the components of the reverse activation energy.

On the basis of the lead compound 4-(N-acetylamino)-3-guanidinobenzoic acid (BANA 113), which inhibits influenza A sialidase with a Ki of 2.5 microM, several novel aromatic inhibitors of influenza sialidases were designed. In this study the N-acetyl group of BANA 113 was replaced with a 2-pyrrolidinone ring, which was designed in part to offer opportunities for introduction of spatially directed side chains that could potentially interact with the 4-, 5-, and/or 6-subsites of sialidase. While the parent structure 1-(4-carboxy-2-guanidinophenyl)pyrrolidin-2-one (8) was only a modest inhibitor of sialidase, the introduction of a hydroxymethyl or bis(hydroxymethyl) substituent at the C5' position of the 2-pyrrolidinone ring resulted in inhibitors (9 and 12, respectively) with low micromolar activity. Crystal structures of these inhibitors in complex with sialidase demonstrated that the substituents at the 5'-position of the 2-pyrrolidinone ring interact in the 4- and/or 5-subsites of the enzyme. Replacement of the guanidine in 12 with a hydrophobic 3-pentylamino group resulted in a large enhancement in binding to produce an inhibitor (14) with an IC50 of about 50 nM against influenza A sialidase, although the inhibition of influenza B sialidase was 2000-fold less. This represents the first reported example of a simple, achiral benzoic acid with potent (low nanomolar) activity as an inhibitor of influenza sialidase.