RESUMEN

This study examined the antiradical activity of three synthesized coumarin derivatives: (E)-3-(1-((2-hydroxyphenyl)amino)ethylidene)-2,4-dioxochroman-7-yl acetate (A1-OH), (E)-3-(1-((3-hydroxyphenyl)amino)ethylidene)-2,4-dioxochroman-7-yl acetate (A2-OH), and (E)-3-(1-((4-hydroxyphenyl)amino)ethylidene)-2,4-dioxochroman-7-yl acetate (A3-OH) against HOO•/O2•- radical species. The investigation included electron spin resonance (ESR) measurements and a DFT kinetic study. Thermodynamic and kinetic parameters of antiradical mechanisms-Formal Hydrogen Atom Transfer (f-HAT), Radical Adduct Formation (RAF), Sequential Proton Loss followed by Electron Transfer (SPLET), and Single-Electron Transfer followed by Proton Transfer (SET-PT)-were evaluated using the Quantum Mechanics-based test for Overall Free Radical Scavenging Activity (QM-ORSA) under physiological conditions. ESR results indicated antiradical activity decreased in the sequence A1-OH (58.7%) > A2-OH (57.5%) > A3-OH (53.1%). Kinetic analysis revealed the f-HAT mechanism dominated HOO• inactivation. A newly formulated Sequential Proton Loss followed by Radical Adduct Formation (SPL-RAF) mechanism described interactions with O2•-. The activity toward O2•- was A2-OH (1.26 × 106 M-1s-1) > A3-OH (7.71 × 105 M-1s-1) > A1-OH (4.22 × 105 M-1s-1). Molecular docking and dynamics studies tested inhibitory capability against enzymes producing reactive species: Lipoxygenase (LOX), Myeloperoxidase (MPO), NAD(P)H oxidase (NOX), and Xanthine Oxidase (XOD). Affinity to enzymes decreased in the order: XOD > LOX > NOX > MPO.

RESUMEN

The cellular environment plays a significant role in low energy electron-mediated radiation damage to genetic materials. In this study, we have modeled the effect of the bulk medium on electron attachment to nucleobases in diethylene glycol (DEG) using uracil as a test case, in accordance with recent experimental work on the observation of dissociative quasi-free electron attachment to nucleoside via excited anion radical in solution (in DEG). Our EOM-CCSD-based quantum mechanical/molecular mechanical (QM/MM) simulations indicate that the electron scavenging by uracil in DEG is much slower than that observed in the aqueous medium due to its viscosity. This work also establishes that a doorway mechanism exists in uracil microsolvated and bulk solvated with DEG, with the dipole-bound state and solvent-bound state acting as doorway states, respectively.

RESUMEN

Epidermal growth factor receptor (EGFR) kinase has been implicated in the uncontrolled cell growth associated with non-small cell lung cancer (NSCLC). This has prompted the development of 3 generations of EGFR inhibitors over the last 2 decades due to the rapid development of drug resistance issues caused by clinical mutations, including T790M, L858R and the double mutant T790M & L858R. In this work we report the design, preparation and biological assessment of new irreversible 2,4-diaminopyrimidine-based inhibitors of EGFR kinase. Twenty new compounds have been prepared and evaluated which incorporate a range of electrophilic moieties. These include acrylamide, 2-chloroacetamide and (2E)-3-phenylprop-2-enamide, to allow reaction with residue Cys797. In addition, more polar groups have been incorporated to provide a better balance of physical properties than clinical candidate Rociletinib. Inhibitory activities against EGFR wildtype (WT) and EGFR T790M & L858R have been evaluated along with cytotoxicity against EGFR-overexpressing (A549, A431) and normal cell lines (HepG2). Selectivity against JAK3 kinase as well as physicochemical properties determination (logD7.4 and phosphate buffer solubility) have been used to profile the compounds. We have identified 20, 21 and 23 as potent mutant EGFR inhibitors (≤20 nM), with comparable or better selectivity over WT EGFR, and lower activity at JAK3, than Osimertinib or Rociletinib. Compounds 21 displayed the best combination of EGFR mutant activity, JAK3 selectivity, cellular activity and physicochemical properties. Finally, kinetic studies on 21 were performed, confirming a covalent mechanism of action at EGFR.

RESUMEN

Among azobenzene derivatives, azothiophenes represent a relatively recent family of compounds that exhibit similar characteristics as dyes and photoreactive systems. Their technological applications are extensive thanks to the additional design flexibility conferred by the heteroaromatic ring. In this study, we present a comprehensive investigation of the structural and electronic properties of novel dyes derived from 3-thiophenamine, utilizing a multilevel approach. We thoroughly examined the potential energy surfaces of the E and Z isomers for three molecules, each bearing different substituents on the phenyl ring at the para position relative to the diazo group. This exploration was conducted through quantum chemistry calculations at various levels of theory, employing a continuum solvent model. Subsequently, we incorporated an explicit solvent (a dimethyl sulfoxide-water mixture) to simulate the most stable isomers using classical molecular dynamics, delivering a clear picture of the local solvation structure and intermolecular interactions. Finally, a hybrid quantum mechanics/molecular mechanics (QM/MM) approach was employed to accurately describe the evolution of the solutes' properties within their environment, accounting for finite temperature effects.

RESUMEN

Melanopsin is a photopigment belonging to the G Protein-Coupled Receptor (GPCR) family expressed in a subset of intrinsically photosensitive retinal ganglion cells (ipRGCs) and responsible for a variety of processes. The bistability and, thus, the possibility to function under low retinal availability would make melanopsin a powerful optogenetic tool. Here, we aim to utilize mouse melanopsin to trigger macrophage migration by its subcellular optical activation with localized blue light, while simultaneously imaging the migration with red light. To reduce melanopsin's red light sensitivity, we employ a combination of in silico structure prediction and automated quantum mechanics/molecular mechanics modeling to predict minimally invasive mutations to shift its absorption spectrum towards the shorter wavelength region of the visible spectrum without compromising the signaling efficiency. The results demonstrate that it is possible to achieve melanopsin mutants that resist red light-induced activation but are activated by blue light and display properties indicating preserved bistability. Using the A333T mutant, we show that the blue light-induced subcellular melanopsin activation triggers localized PIP3 generation and macrophage migration, which we imaged using red light, demonstrating the optogenetic utility of minimally engineered melanopsins.

Asunto(s)

RESUMEN

Industrial fermentation applications typically require enzymes that exhibit high stability and activity at high temperatures. However, efforts to simultaneously improve these properties are usually limited by a trade-off between stability and activity. This report describes a computational strategy to enhance both activity and thermal stability of the mesophilic organophosphate-degrading enzyme, methyl parathion hydrolase (MPH). To predict hotspot mutation sites, we assembled a library of features associated with the target properties for each residue and then prioritized candidate sites by hierarchical clustering. Subsequent in silico screening with multiple algorithms to simulate selective pressures yielded a subset of 23 candidate mutations. Iterative parallel screening of mutations that improved thermal stability and activity yielded, MPHase-m5b, which exhibited 13.3 °C higher Tm and 4.2 times higher catalytic activity than wild-type (WT) MPH over a wide temperature range. Systematic analysis of crystal structures, molecular dynamics (MD) simulations, and Quantum Mechanics/Molecular Mechanics (QM/MM) calculations revealed a wider entrance to the active site that increased substrate access with an extensive network of interactions outside the active site that reinforced αß/ßα sandwich architecture to improve thermal stability. This study thus provides an advanced, rational design framework to improve efficiency in engineering highly active, thermostable biocatalysts for industrial applications.

RESUMEN

INTRODUCTION: Molecular Dynamics (MD) simulations can support mechanism-based drug design. Indeed, MD simulations by capturing biomolecule motions at finite temperatures can reveal hidden binding sites, accurately predict drug-binding poses, and estimate the thermodynamics and kinetics, crucial information for drug discovery campaigns. Small-Guanosine Triphosphate Phosphohydrolases (GTPases) regulate a cascade of signaling events, that affect most cellular processes. Their deregulation is linked to several diseases, making them appealing drug targets. The broad roles of small-GTPases in cellular processes and the recent approval of a covalent KRas inhibitor as an anticancer agent renewed the interest in targeting small-GTPase with small molecules. AREA COVERED: This review emphasizes the role of MD simulations in elucidating small-GTPase mechanisms, assessing the impact of cancer-related variants, and discovering novel inhibitors. EXPERT OPINION: The application of MD simulations to small-GTPases exemplifies the role of MD simulations in the structure-based drug design process for challenging biomolecular targets. Furthermore, AI and machine learning-enhanced MD simulations, coupled with the upcoming power of quantum computing, are promising instruments to target elusive small-GTPases mutations and splice variants. This powerful synergy will aid in developing innovative therapeutic strategies associated to small-GTPases deregulation, which could potentially be used for personalized therapies and in a tissue-agnostic manner to treat tumors with mutations in small-GTPases.

RESUMEN

Introduction: In clinical practice, phosphodiesterase 5 (PDE5) inhibitors are commonly used to treat erectile dysfunction and pulmonary arterial hypertension. However, due to the high structural similarity between PDE5 and Phosphodiesterase 6 (PDE6), there is a risk that existing drugs will cause off-target effects on PDE6 resulting in visual disorders such as low visual acuity and color blindness. Previous research on the selectivity of PDE5 inhibitors focused on marketed drugs such as sildenafil and tadalafil. Methods: In this study, a highly selective PDE5 inhibitor, ligand3, was used as the subject, and molecular docking, molecular dynamics simulations, MM-GBSA, alanine scanning, and independent gradient model analysis were employed to investigate the biological mechanism underlying the selectivity of PDE5 inhibitors. Results and Discussion: The present work revealed that the binding mode of ligand3 to the PDE5A and PDE6C targets was distinctly different. Ligand3 exhibited stronger coulombic forces when binding to PDE5A, while showing stronger van der waals forces when binding to PDE6C. Ligand3 binds more deeply at the active site of PDE5A than at PDE6C, allowing its side chains to effectively bind to the critical TYR612, whereas in the case of the shallow binding to PDE6C, ligand3 lacks a similar effect. Mechanism investigations of highly selective inhibitors through computational simulation might provide an insight into potent treatment of drugs.

RESUMEN

Objective: The primary objective of this study was to investigate the role and regulatory mechanisms of Platelet-derived Growth Factor Subunit B (PDGFB) in muscle differentiation. Methods: In this study, a vector for PDGFB was designed and transfected into quail muscle cells to investigate its role and regulatory mechanism during muscle formation. To investigate the inhibitory mechanisms of PDGFB on myogenic differentiation, the mRNA expression levels of various genes and the phosphorylation of extracellular signal-regulated kinase 1/2 (ERK 1/2), both known to regulate muscle development and differentiation were compared. Results: PDGFB-overexpressed (OE) cells formed morphologically shorter and thinner myotubes and demonstrated a smaller total myotube area than did the control cells. This result was also confirmed at the molecular level by a reduced amount of myosin heavy chain protein in the PDGFB-OE cells. Therefore, PDGFB inhibits the differentiation of muscle cells. Additionally, the expression of myogenin (MYOG) significantly decreased in the PDGFB-OE cells on days 2 and 4 compared with that in the control cells. The phosphorylation of ERK 1/2, an upstream protein that inhibits MYOG expression, increased in the PDGFB-OE cells on day 4 compared with that in the control cells. The decreased expression of MYOG in the PDGFB-OE cells increased by inhibition ERK 1/2 phosphorylation. Conclusion: PDGFB may suppress myogenesis by reducing MYOG expression through ERK 1/2 phosphorylation. These findings can help understand muscle differentiation and potentially improve poultry meat production.

RESUMEN

The increase in the resistance of mutant strains of Neisseria gonorrhoeae to the antibiotic ceftriaxone is pronounced in the decrease in the second-order acylation rate constant, k2/KS, by penicillin-binding protein 2 (PBP2). These changes can be caused by both the decrease in the acylation rate constant, k2, and the weakening of the binding affinity, i.e., an increase in the substrate constant, KS. A501X mutations in PBP2 affect second-order acylation rate constants. The PBP2A501V variant exhibits a higher k2/KS value, whereas for PBP2A501R and PBP2A501P variants, these values are lower. We performed molecular dynamic simulations with both classical and QM/MM potentials to model both acylation energy profiles and conformational dynamics of four PBP2 variants to explain the origin of k2/KS changes. The acylation reaction occurs in two elementary steps, specifically, a nucleophilic attack by the oxygen atom of the Ser310 residue and C-N bond cleavage in the ß-lactam ring accompanied by the elimination of the leaving group of ceftriaxone. The energy barrier of the first step increases for PBP2 variants with a decrease in the observed k2/KS value. Submicrosecond classic molecular dynamic trajectories with subsequent cluster analysis reveal that the conformation of the ß3-ß4 loop switches from open to closed and its flexibility decreases for PBP2 variants with a lower k2/KS value. Thus, the experimentally observed decrease in the k2/KS in A501X variants of PBP2 occurs due to both the decrease in the acylation rate constant, k2, and the increase in KS.

Asunto(s)

RESUMEN

Formate oxidase (FOx), which contains 8-formyl flavin adenine dinucleotide (FAD), exhibits a distinct advantage in utilizing ambient oxygen molecules for the oxidation of formic acid compared to other glucose-methanol-choline (GMC) oxidoreductase enzymes that contain only the standard FAD cofactor. The FOx-mediated conversion of FAD to 8-formyl FAD results in an approximate 10-fold increase in formate oxidase activity. However, the mechanistic details underlying the autocatalytic formation of 8-formyl FAD are still not well understood, which impedes further utilization of FOx. In this study, we employ molecular dynamics simulation, QM/MM umbrella sampling simulation, enzyme activity assay, site-directed mutagenesis, and spectroscopic analysis to elucidate the oxidation mechanism of FAD to 8-formyl FAD. Our results reveal that a catalytic water molecule, rather than any catalytic amino acids, serves as a general base to deprotonate the C8 methyl group on FAD, thus facilitating the formation of a quinone-methide tautomer intermediate. An oxygen molecule subsequently oxidizes this intermediate, resulting in a C8 methyl hydroperoxide anion that is protonated and dissociated to form OHC-RP and OH-. During the oxidation of FAD to 8-formyl FAD, the energy barrier for the rate-limiting step is calculated to be 22.8 kcal/mol, which corresponds to the required 14-hour transformation time observed experimentally. Further, the elucidated oxidation mechanism reveals that the autocatalytic formation of 8-formyl FAD depends on the proximal arginine and serine residues, R87 and S94, respectively. Enzymatic activity assay validates that the mutation of R87 to lysine reduces the kcat value to 75% of the wild-type, while the mutation to histidine results in a complete loss of activity. Similarly, the mutant S94I also leads to the deactivation of enzyme. This dependency arises because the nucleophilic OH- group and the quinone-methide tautomer intermediate are stabilized through the noncovalent interaction provided by R87 and S94. These findings not only explain the mechanistic details of each reaction step but also clarify the functional role of R87 and S94 during the oxidative maturation of 8-formyl FAD, thereby providing crucial theoretical support for the development of novel flavoenzymes with enhanced redox properties.

RESUMEN

Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.

Asunto(s)

RESUMEN

In the green alga Chlamydomonas reinhardtii, hydrogenase HydA1 converts protons and electrons to H2 at the H-cluster, which includes a [4Fe-4S] cluster linked to a [2Fe] cluster. The yield of H2 is limited by the electron transfer to HydA1, mediated by the iron-sulfur unit of a photosynthetic electron transfer ferredoxin (PetF). In this study, I have investigated by molecular dynamics and the hybrid quantum mechanics/molecular mechanics method two canonical iron-sulfur peptides (PM1 and FBM) that hold potential as PetF replacements. Using a docking approach, I predict that the distance between the two iron-sulfur clusters in FBM/HydA1 is shorter than in PM1/HydA1, ensuring a greater electron transfer rate. This finding is in line with the reported higher H2 production rates for FBM/HydA1. I also show that the redox potential of these peptides, and therefore their electron transfer properties, can be changed by single-residue mutations in the secondary coordination sphere of their cluster. In particular, I have designed a PM1 variant that disrupts the hydrogen-bonding network between water and the cluster, shifting the redox potential negatively compared to PM1. These results will guide experiments aimed at replacing PetF with peptides that can unlock the biotechnological potential of the alga.

RESUMEN

Nucleic acid processing enzymes use a two-Mg2+-ion motif to promote the formation and cleavage of phosphodiester bonds. Yet, recent evidence demonstrates the presence of spatially conserved second-shell cations surrounding the catalytic architecture of proteinaceous and RNA-dependent enzymes. The RNase mitochondrial RNA processing (MRP) complex, which cleaves the ribosomal RNA (rRNA) precursor at the A3 cleavage site to yield mature 5'-end of 5.8S rRNA, hosts in the catalytic core one atypically-located Mg2+ ion, in addition to the ions forming the canonical catalytic motif. Here, we employ biased quantum classical molecular dynamics simulations of RNase MRP to discover that the third Mg2+ ion inhibits the catalytic process. Instead, its displacement in favour of a second-shell monovalent K+ ion propels phosphodiester bond cleavage by enabling the formation of a specific hydrogen bonding network that mediates the essential proton transfer step. This study points to a direct involvement of a transient K+ ion in the catalytic cleavage of the phosphodiester bond and implicates cation trafficking as a general mechanism in nucleic acid processing enzymes and ribozymes.

RESUMEN

This paper investigates the esterase activity of minimalist amyloid fibers composed of short seven-residue peptides, IHIHIHI (IH7) and IHIHIQI (IH7Q), with a particular focus on the role of the sixth residue position within the peptide sequence. Through computational simulations and analyses, we explore the molecular mechanisms underlying catalysis in these amyloid-based enzymes. Contrary to initial hypotheses, our study reveals that the twist angle of the fiber, and thus the catalytic site's environment, is not notably affected by the sixth residue. Instead, the sixth residue interacts with the p-nitrophenylacetate (pNPA) substrate, particularly through its -NO2 group, potentially enhancing catalysis. Quantum mechanics/molecular mechanics (QM/MM) simulations of the reaction mechanism suggest that the polarizing effect of glutamine enhances catalytic activity by forming a stabilizing network of hydrogen bonds with pNPA, leading to lower energy barriers and a more exergonic reaction. Our findings provide valuable insights into the intricate interplay between peptide sequence, structural arrangement, and catalytic function in amyloid-based enzymes, offering potentially valuable information for the design and optimization of biomimetic catalysts.

RESUMEN

This study investigates the application of quantum mechanical (QM) and multiscale computational methods in understanding the reaction mechanisms and kinetics of SN2 reactions involving methyl iodide with NH2OH and NH2O-, as well as the Claisen rearrangement of 8-(vinyloxy)dec-9-enoate. Our aim is to evaluate the accuracy and effectiveness of these methods in predicting experimental outcomes for these organic reactions. We achieve this by employing QM-only calculations and several hybrids of QM and molecular mechanics (MM) methods, namely QM/MM, QM1/QM2, and QM1/QM2/MM methodologies. For the SN2 reactions, our results demonstrate the importance of explicitly including solvent effects in the calculations to accurately reproduce the transition state geometry and energetics. The multiscale methods, particularly QM/MM and QM1/QM2, show promising performance in predicting activation energies. Moreover, we observe that the size of the MM active region significantly affects the accuracy of calculated activation energies, highlighting the need for careful consideration during the setup of multiscale calculations. In the case of the Claisen rearrangement, both QM-only and multiscale methods successfully reproduce the proposed reaction mechanism. However, the activation free energies calculated using a continuum solvation model, based on single-point calculations of QM-only structures, fail to account for solvent effects. On the other hand, multiscale methods more accurately capture the impact of solvents on activation free energies, with systematic error correction enhancing the accuracy of the results. Furthermore, we introduce a Python code for setting up multiscale calculations with ORCA, which is available on GitHub at https://github.com/iranimehdi/pdbtoORCA .

RESUMEN

We investigate the inhibition mechanism between pomotrelvir and the SARS-CoV-2 main protease using molecular mechanics and quantum mechanics/molecular mechanics simulations. Alchemical transformations where each Pi group of pomotrelvir was transformed into its counterpart in nirmatrelvir were performed to unravel the individual contribution of each group to the binding and reaction processes. We have shown that while a γ-lactam ring is preferred at position P1, a δ-lactam ring is a good alternative for the design of inhibitors for variants presenting mutations at position 166. For the P2 position, tertiary amines are preferred with respect to secondary amines. Flexible side chains at the P2 position can disrupt the preorganization of the active site, favouring the exploration of non-reactive conformations. The substitution of the P2 group of pomotrelvir by that of nirmatrelvir resulted in a compound, named as C2, that presents a better binding free energy and a higher population of reactive conformations in the Michaelis complex. Analysis of the chemical reaction to form the covalent complex has shown a similar reaction mechanism and activation free energies for pomotrelvir, nirmatrelvir and C2. We hope that these findings could be useful to design better inhibitors to fight present and future variants of the SARS-CoV-2 virus.

RESUMEN

Stimulus-responsive organic room temperature phosphorescence (RTP) materials exhibit variations in their luminescent characteristics (lifetime and efficiency) upon exposure to external stimuli, including force, heat, light and acid-base conditions, the development of stimulus-responsive RTP molecules becomes imperative. However, the inner responsive mechanism is unclear, theoretical investigations to reveal the relationship among hydrostatic pressures, molecular structures and photophysical properties are highly desired. Herein, taking the Se-containing RTP molecule (SeAN) as a model, based on the dispersion corrected density functional theory (DFT-D), the combined quantum mechanics and molecular dynamics (QM/MM) method and thermal vibration correlation function (TVCF) theory, the influences of hydrostatic pressure on molecular structures, transition properties as well as lifetimes and efficiencies of RTP molecule are theoretically studied. Results show that extended lifetime and enhanced efficiency are observed at 2 Gpa compared with molecule at normal pressure, and this is related with the small reorganization energy and large oscillator strength. Moreover, due to the small energy gap (0.34 eV) and remarkable spin-orbit coupling (SOC) constant (8.56 cm-1) between first singlet excited state and triplet state, fast intersystem crossing (ISC) process is determined for molecule at 6 Gpa. Furthermore, the intermolecular interactions are visualized using independent gradient model based on Hirshfeld partition (IGMH) and the changes of molecular packing modes, SOC values, lifetimes and efficiencies with pressures are detected. These results reveal the relationship between molecular structures and RTP properties. Our work provides theoretical insights into the hydrostatic pressure response mechanism and could promote the development new efficient stimulus-responsive molecules.

RESUMEN

7H-Dibenzo[c,g]carbazole (DBC) is a prevalent environmental contaminant that induces tumorigenesis in several experimental animals. Recently, it has been utilized to develop high-performance solar cells and organic phosphorescent materials. It is imperative to strengthen investigations of DBC metabolism to understand its potential risks to human health. In this study, human CYP1A1 was employed as the metabolic enzyme to investigate the metabolic mechanism of DBC by molecular docking, molecular dynamics (MD) simulation, and quantum mechanical (QM) calculation. The results indicate that DBC binds to CYP1A1 in two modes (mode 1 and mode 2) mainly through nonpolar solvation energies (ΔGnonpolar). The formation of the two binding modes is attributed to the anchoring effect of the hydrogen bond formed by DBC with Asp320 (mode 1) or Ser116 (mode 2). Mode 1 is a "reactive" conformation, while mode 2 is not considered a "reactive" conformation. C5 is identified as the dominant site, and the pyrrole nitrogen cannot participate in the metabolism. DBC is metabolized mainly by a distinct electrophilic addition-rearrangement mechanism, with an energy barrier of 21.74 kcal/mol. The results provide meaningful insights into the biometabolic process of DBC and contribute to understanding its environmental effects and health risks.

Asunto(s)

RESUMEN

Terpene Synthases (TPS) catalyze the formation of multicyclic, complex terpenes and terpenoids from linear substrates. Molecular docking is an important research tool that can further our understanding of TPS multistep mechanisms and guide enzyme design. Standard docking programs are not well suited to tackle the unique challenges of TPS, like the many chemical steps which form multiple stereo-centers, the weak dispersion interactions between the isoprenoid chain and the hydrophobic region of the active site, description of carbocation intermediates, and finding mechanistically meaningful sets of docked poses. To address these and other unique challenges, we developed the multistate, multiscale docking program EnzyDock and used it to study many TPS and other enzymes. In this review we discuss the unique challenges of TPS, the special features of EnzyDock developed to address these challenges and demonstrate its successful use in ongoing research on the bacterial TPS CotB2.

Asunto(s)