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.

Magnesium ions (Mg2+) are crucial in class II terpene cyclases that utilize substrates with diphosphate groups. Interestingly, these enzymes catalyze reactions without cleaving the diphosphate group, instead initiating the reaction through protonation. In our recent research, we discovered a novel class II sesquiterpene cyclase in Streptomyces showdoensis. Notably, we determined its crystal structure and identified Mg2+ within its active site. This finding has shed light on the previously elusive question of Mg2+ binding in class II terpene cyclases. In this chapter, we outline our methods for discovering this novel enzyme, including steps for its purification, crystallization, and kinetic analysis.

Terpenes are a diverse class of natural products which have long been sought after for their chemical properties as medicine, perfumes, and for food flavoring. Computational docking studies of terpene mechanisms have been a challenge due to the lack of strong directing groups which many docking programs rely on. In this chapter, we dive into our computational method Terdockin (Terpene-Docking) as a successful methodology in modeling terpene synthase mechanisms. This method could also be used as inspiration for any multi-ligand docking project.

Proteases represent common targets in combating infectious diseases, including COVID-19. The 3-chymotrypsin-like protease (3CLpro) is a validated molecular target for COVID-19, and it is key for developing potent and selective inhibitors for inhibiting viral replication of SARS-CoV-2. In this review, we discuss structural relationships and diverse subsites of 3CLpro, shedding light on the pivotal role of dimerization and active site architecture in substrate recognition and catalysis. Our analysis of bioinformatics and other published studies motivated us to investigate a novel catalytic mechanism for the SARS-CoV-2 polyprotein cleavage by 3CLpro, centering on the triad mechanism involving His41-Cys145-Asp187 and its indispensable role in viral replication. Our hypothesis is that Asp187 may participate in modulating the pKa of the His41, in which catalytic histidine may act as an acid and/or a base in the catalytic mechanism. Recognizing Asp187 as a crucial component in the catalytic process underscores its significance as a fundamental pharmacophoric element in drug design. Next, we provide an overview of both covalent and non-covalent inhibitors, elucidating advancements in drug development observed in preclinical and clinical trials. By highlighting various chemical classes and their pharmacokinetic profiles, our review aims to guide future research directions toward the development of highly selective inhibitors, underscore the significance of 3CLpro as a validated therapeutic target, and propel the progression of drug candidates through preclinical and clinical phases.

Hepatitis B virus (HBV) remains a global health threat. Ribonuclease H (RNase H), part of the virus polymerase protein, cleaves the pgRNA template during viral genome replication. Inhibition of RNase H activity prevents (+) DNA strand synthesis and results in the accumulation of non-functional genomes, terminating the viral replication cycle. RNase H, though promising, remains an under-explored drug target against HBV. We previously reported the identification of a series of N-hydroxypyridinedione (HPD) imines that effectively inhibit the HBV RNase H. In our effort to further explore the HPD scaffold, we designed, synthesized, and evaluated 18 novel HPD oximes, as well as 4 structurally related minoxidil derivatives and 2 barbituric acid counterparts. The new analogs were docked on the RNase H active site and all proved able to coordinate the two Mg2+ ions in the catalytic site. All of the new HPDs effectively inhibited the viral replication in cell assays exhibiting EC50 values in the low μM range (1.1-7.7 μM) with low cytotoxicity, resulting in selectivity indexes (SI) of up to 92, one of the highest reported to date among HBV RNase H inhibitors. Our findings expand the structure-activity relationships on the HPD scaffold, facilitating the development of even more potent anti-HBV agents.

The development of new compounds to treat Chagas disease is imperative due to the adverse effects of current drugs and their low efficacy in the chronic phase. This study aims to investigate nitroisoxazole derivatives that produce oxidative stress while modifying the compounds\' lipophilicity, affecting their ability to fight trypanosomes. The results indicate that these compounds are more effective against the epimastigote form of T. cruzi, with a 52 ± 4% trypanocidal effect for compound 9. However, they are less effective against the trypomastigote form, with a 15 ± 3% trypanocidal effect. Additionally, compound 11 interacts with a higher number of amino acid residues within the active site of the enzyme cruzipain. Furthermore, it was also found that the presence of a nitro group allows for the generation of free radicals; likewise, the large size of the compound enables increased interaction with aminoacidic residues in the active site of cruzipain, contributing to trypanocidal activity. This activity depends on the size and lipophilicity of the compounds. The study recommends exploring new compounds based on the nitroisoxazole skeleton, with larger substituents and lipophilicity to enhance their trypanocidal activity.

Bacterial nitroreductase enzymes capable of activating imaging probes and prodrugs are valuable tools for gene-directed enzyme prodrug therapies and targeted cell ablation models. We recently engineered a nitroreductase (E. coli NfsB F70A/F108Y) for the substantially enhanced reduction of the 5-nitroimidazole PET-capable probe, SN33623, which permits the theranostic imaging of vectors labeled with oxygen-insensitive bacterial nitroreductases. This mutant enzyme also shows improved activation of the DNA-alkylation prodrugs CB1954 and metronidazole. To elucidate the mechanism behind these enhancements, we resolved the crystal structure of the mutant enzyme to 1.98 Å and compared it to the wild-type enzyme. Structural analysis revealed an expanded substrate access channel and new hydrogen bonding interactions. Additionally, computational modeling of SN33623, CB1954, and metronidazole binding in the active sites of both the mutant and wild-type enzymes revealed key differences in substrate orientations and interactions, with improvements in activity being mirrored by reduced distances between the N5-H of isoalloxazine and the substrate nitro group oxygen in the mutant models. These findings deepen our understanding of nitroreductase substrate specificity and catalytic mechanisms and have potential implications for developing more effective theranostic imaging strategies in cancer treatment.

Ferrochelatase (FECH) is the terminal enzyme in human heme biosynthesis, catalyzing the insertion of ferrous iron into protoporphyrin IX (PPIX) to form protoheme IX (Heme). Phosphorylation increases the activity of FECH, and it has been confirmed that the activity of FECH phosphorylated at T116 increases. However, it remains unclear whether the T116 site and other potential phosphorylation modification sites collaboratively regulate the activity of FECH. In this study, we identified a new phosphorylation site, T218, and explored the allosteric effects of unphosphorylated (UP), PT116, PT218, and PT116 + PT218 states on FECH in the presence and absence of substrates (PPIX and Heme) using molecular dynamics (MD) simulations. Binding free energies were evaluated with the MM/PBSA method. Our findings indicate that the PT116 + PT218 state exhibits the lowest binding free energy with PPIX, suggesting the strongest binding affinity. Additionally, this state showed a higher binding free energy with Heme compared to UP, which facilitates Heme release. Moreover, employing multiple analysis methods, including free energy landscape (FEL), principal component analysis (PCA), dynamic cross-correlation matrix (DCCM), and hydrogen bond interaction analysis, we demonstrated that phosphorylation significantly affects the dynamic behavior and binding patterns of substrates to FECH. Insights from this study provide valuable theoretical guidance for treating conditions related to disrupted heme metabolism, such as various porphyrias and iron-related disorders.

Amino acid deprivation therapy (AADT) is a novel anticancer therapy, considered nontoxic and selective. Thermophilic L-asparaginase enzymes display high stability and activity at elevated temperatures. However, they are of limited use in clinical applications because of their low substrate affinity and reduced activity under physiological conditions, which may necessitate an improved dosage, leading to side effects and greater costs. Thus, in an attempt to improve the activity of L-Asn at 37 °C, with the use of a semi-rational design, eight active-site mutants of Thermococcus litoralis DSM 5473 L-asparaginase Tli10209 were developed. T70A exhibited a 5.11-fold increase compared with the wild enzyme in physiological conditions. Double-mutant enzymes were created by combining mutants with higher hydrolysis activity. T70A/F36Y, T70A/K48L, and T70A/D50G were enhanced by 5.59-, 6.38-, and 5.58-fold. The immobilized enzyme applied in MCF-7 breast cancer cells only required one-seventh of the dose of the free enzyme to achieve the same inhibition rate under near-infrared irradiation. This provides a proof of concept that it is possible to reduce the consumption of L-Asn by improving its activity, thus providing a method to manage side effects.

Insertion sequence (IS) elements are the simplest autonomous transposable elements found in prokaryotic genomes1. We recently discovered that IS110 family elements encode a recombinase and a non-coding bridge RNA (bRNA) that confers modular specificity for target DNA and donor DNA through two programmable loops2. Here we report the cryo-electron microscopy structures of the IS110 recombinase in complex with its bRNA, target DNA and donor DNA in three different stages of the recombination reaction cycle. The IS110 synaptic complex comprises two recombinase dimers, one of which houses the target-binding loop of the bRNA and binds to target DNA, whereas the other coordinates the bRNA donor-binding loop and donor DNA. We uncovered the formation of a composite RuvC-Tnp active site that spans the two dimers, positioning the catalytic serine residues adjacent to the recombination sites in both target and donor DNA. A comparison of the three structures revealed that (1) the top strands of target and donor DNA are cleaved at the composite active sites to form covalent 5\'-phosphoserine intermediates, (2) the cleaved DNA strands are exchanged and religated to create a Holliday junction intermediate, and (3) this intermediate is subsequently resolved by cleavage of the bottom strands. Overall, this study reveals the mechanism by which a bispecific RNA confers target and donor DNA specificity to IS110 recombinases for programmable DNA recombination.