CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed.

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Mechanistic docking in terpene synthases using EnzyDock.

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.

Fluorinated Co-Solvents Enable Excellent Performances of Practical Cells Comprising LixSiO-Graphite Composite Anodes and LiNi0.89Co0.05Mn0.05Al0.01O2 (NCMA) Cathodes.

Li-ion batteries based on high specific capacity LixSiO-Graphite anodes and LiNi0.89Co0.05 Mn0.05Al0.01O2 (NCMA) cathodes may have numerous practical applications owing to high energy density without a necessary compromise on safety. SiO, which is an attractive Li insertion anode material, offers more cycling stability than Si and a higher capacity than graphite. Therefore, a new trend has emerged for developing composite C-Si anodes, possessing the excellent cyclability of graphite coupled with high capacity SiO. The composite structure described herein prevents the volume expansion of SiO and maintains the structural integrity during prolonged cycling. However, graphite electrodes suffer from exfoliation in propylene carbonate (PC) based electrolyte solutions, which avoids well known safety benefits related to a possible use of PC based electrolyte solutions in all kinds of Li batteries. Herein, it is reported that trifluoro propylene carbonate (TFPC) is compatible with graphite anodes. New electrolyte formulations are developed and tested containing fluorinated co-solvents and compared the performance of several electrolyte solutions, including conventional alkyl carbonates-based solutions in full Li-ion cells, which included LixSiO-Graphite anodes and LiNi0.89Co0.05Mn0.05Al0.01O2 (NCMA) cathodes. Cells with new electrolyte solutions developed herein demonstrated nearly twice capacity retention in prolonged cycling experiments compared to similar reference cells containing conventional electrolyte solutions.

Differential Substrate Sensing in Terpene Synthases from Plants and Microorganisms: Insight from Structural, Bioinformatic, and EnzyDock Analyses.

Terpene synthases (TPSs) catalyze the first step in the formation of terpenoids, which comprise the largest class of natural products in nature. TPSs employ a family of universal natural substrates, composed of isoprenoid units bound to a diphosphate moiety. The intricate structures generated by TPSs are the result of substrate binding and folding in the active site, enzyme-controlled carbocation reaction cascades, and final reaction quenching. A key unaddressed question in class I TPSs is the asymmetric nature of the diphosphate-(Mg2+)3 cluster, which forms a critical part of the active site. In this asymmetric ion cluster, two diphosphate oxygen atoms protrude into the active site pocket. The substrate hydrocarbon tail, which is eventually molded into terpenes, can bind to either of these oxygen atoms, yet to which is unknown. Herein, we employ structural, bioinformatics, and EnzyDock docking tools to address this enigma. We bring initial data suggesting that this difference is rooted in evolutionary differences between TPSs. We hypothesize that this alteration in binding, and subsequent chemistry, is due to TPSs originating from plants or microorganisms. We further suggest that this difference can cast light on the frequent observation that the chiral products or intermediates of plant and bacterial terpene synthases represent opposite enantiomers.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development.

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Monte Carlo-Simulated Annealing and Machine Learning-Based Funneled Approach for Finding the Global Minimum Structure of Molecular Clusters.

Understanding the physical underpinnings and geometry of molecular clusters is of great importance in many fields, ranging from studying the beginning of the universe to the formation of atmospheric particles. To this end, several approaches have been suggested, yet identifying the most stable cluster geometry (i.e., global potential energy minimum) remains a challenge, especially for highly symmetric clusters. Here, we suggest a new funneled Monte Carlo-based simulated annealing (SA) approach, which includes two key steps: generation of symmetrical clusters and classification of the clusters according to their geometry using machine learning (MCSA-ML). We demonstrate the merits of the MCSA-ML method in comparison to other approaches on several Lennard-Jones (LJ) clusters and four molecular clusters-Ser8(Cl-)2, H+(H2O)6, Ag+(CO2)8, and Bet4Cl-. For the latter of these clusters, the correct structure is unknown, and hence, we compare the experimental and simulated fragmentation patterns, and the fragmentation of the proposed global minimum matches experiments closely. Additionally, based on the fragmentation of the predicted betaine cluster, we were able to identify hitherto unknown neutral fragmentation channels. In comparison to results obtained with other methods, we demonstrated a superior ability of MCSA-ML to predict clusters with high symmetry and similar abilities to predict clusters with asymmetrical structures.

Understanding the competing pathways leading to hydropyrene and isoelisabethatriene.

Terpene synthases are responsible for the biosynthesis of terpenes, the largest family of natural products. Hydropyrene synthase generates hydropyrene and hydropyrenol as its main products along with two byproducts, isoelisabethatrienes A and B. Fascinatingly, a single active site mutation (M75L) diverts the product distribution towards isoelisabethatrienes A and B. In the current work, we study the competing pathways leading to these products using quantum chemical calculations in the gas phase. We show that there is a great thermodynamic preference for hydropyrene and hydropyrenol formation, and hence most likely in the synthesis of the isoelisabethatriene products kinetic control is at play.

Can Anions Be Inserted into MXene?

Despite the continuous progress in the research and development of Ti3C2Tx (MXene) electrodes for high-power batteries and supercapacitor applications, the role of the anions in the electrochemical energy storage and their ability to intercalate between the MXene sheets upon application of positive voltage have not been clarified. A decade after the discovery of MXenes, the information about the possibility of anion insertion into the restacked MXene electrode is still being questioned. Since the positive potential stability range in diluted aqueous electrolytes is severely limited by anodic oxidation of the Ti, the possibility of anion insertion was evaluated in concentrated aqueous electrolyte solutions and aprotic electrolytes as well. To address this issue, we have conducted in situ gravimetric electrochemical quartz crystal microbalance with dissipation monitoring (EQCM-D) measurements in highly concentrated LiCl and LiBr electrolytes, which enable a significant extension of the operation range of the MXene electrodes toward positive potentials. Also, halogens are among the smallest anions and should be easier to intercalate between MXene layers, in comparison to multiatomic anions. On the basis of mass change variations in the positive voltage range and complementary density functional theory calculations, it was demonstrated that insertion of anionic species into MXene, within the range of potentials of interest for capacitive energy storage, is not likely to occur. This can be explained by the strong negative charge on Ti3C2Tx sheets terminated by functional groups.

Benchmarking the Ability of Common Docking Programs to Correctly Reproduce and Score Binding Modes in SARS-CoV-2 Protease Mpro.

The coronavirus SARS-CoV-2 main protease, Mpro, is conserved among coronaviruses with no human homolog and has therefore attracted significant attention as an enzyme drug target for COVID-19. The number of studies targeting Mpro for in silico screening has grown rapidly, and it would be of great interest to know in advance how well docking methods can reproduce the correct ligand binding modes and rank these correctly. Clearly, current attempts at designing drugs targeting Mpro with the aid of computational docking would benefit from a priori knowledge of the ability of docking programs to predict correct binding modes and score these correctly. In the current work, we tested the ability of several leading docking programs, namely, Glide, DOCK, AutoDock, AutoDock Vina, FRED, and EnzyDock, to correctly identify and score the binding mode of Mpro ligands in 193 crystal structures. None of the codes were able to correctly identify the crystal structure binding mode (lowest energy pose with root-mean-square deviation < 2 Å) in more than 26% of the cases for noncovalently bound ligands (Glide: top performer), whereas for covalently bound ligands the top score was 45% (EnzyDock). These results suggest that one should perform in silico campaigns of Mpro with care and that more comprehensive strategies including ligand free energy perturbation might be necessary in conjunction with virtual screening and docking.

Concentration-Dependent Self-Assembly of an Unusually Large Hexameric Hydrogen-Bonded Molecular Cage.

The sizes of available self-assembled hydrogen-bond-based supramolecular capsules and cages are rather limited. The largest systems have volumes of approximately 1400-2300 Å3 . Herein, we report a large, hexameric cage based on intermolecular amide-amide dimerization. The unusual structure with openings, reminiscent of covalently linked cages, is held together by 24 hydrogen bonds. With a diameter of 2.3 nm and a cavity volume of ∼2800 Å3 , the assembly is larger than any previously known capsule/cage structure relying exclusively on hydrogen bonds. The self-assembly process in chlorinated, organic solvents was found to be strongly concentration dependent, with the monomeric form prevailing at low concentrations. Additionally, the formation of host-guest complexes with fullerenes (C60 and C70 ) was observed.

Four-Step Access to the Sesquiterpene Natural Product Presilphiperfolan-1ß-ol and Unnatural Derivatives via Supramolecular Catalysis.

Terpenes constitute one of the most structurally varied classes of natural products. A wide range of these structures are produced in nature by type I terpene cyclase enzymes from one single substrate. However, such reactivity has proven difficult to reproduce in solution with man-made systems. Herein we report the shortest synthesis of the tricyclic sesquiterpene presilphiperfolan-1ß-ol to date, utilizing the supramolecular resorcinarene capsule as catalyst for the key step. This synthetic approach also allows access to unnatural derivatives of the natural product, which would not be accessible through the biosynthetic machinery. Additionally, this study provides useful insight into the biosynthesis of the presilphiperfolanol natural products, including the first experimental evidence consistent with the proposed biosynthetic connection between caryophyllene and the presilphiperfolanols.

Understanding the role of active site residues in CotB2 catalysis using a cluster model.

Terpene cyclases are responsible for the initial cyclization cascade in the multistep synthesis of a large number of terpenes. CotB2 is a diterpene cyclase from Streptomyces melanosporofaciens, which catalyzes the formation of cycloocta-9-en-7-ol, a precursor to the next-generation anti-inflammatory drug cyclooctatin. In this work, we present evidence for the significant role of the active site's residues in CotB2 on the reaction energetics using quantum mechanical calculations in an active site cluster model. The results revealed the significant effect of the active site residues on the relative electronic energy of the intermediates and transition state structures with respect to gas phase data. A detailed understanding of the role of the enzyme environment on the CotB2 reaction cascade can provide important information towards a biosynthetic strategy for cyclooctatin and the biomanufacturing of related terpene structures.

First principles study of electrocatalytic behavior of olivine phosphates with mixed alkali and mixed transition metal atoms.

Lithium transition metal olivine phosphates are well known Li-ion battery cathode materials, but these materials can also be used as electrocatalyst. Recent experimental studies showed that olivine phosphates with mixed alkali metals (Li and Na) and mixed transition metals (Ni and Fe) provide better electrocatalytic activity compared to single alkali and transition metal alternatives. In the current work, we analyzed the role of alkali metals, transition metals and vacancies on the reactivity of a series of olivine phosphates with different stoichiometries using first principles calculations. To this end, we investigated the adsorption of water at the surface of these materials. We found that water binds preferably at Ni surface sites for materials devoid of alkali ion vacancies. We further found correlation between the calculated adsorption energy with experimentally measured overpotentials for a series of olivine phosphates. Additionally, we found correlation between the adsorption energy of the systems with the total charge polarization of surface and adsorbate. To explain the computed trends, we analyzed the occupancies of the partial density of states of the Ni and Fe 3d states and Bader atomic charges.

EnzyDock: Protein-Ligand Docking of Multiple Reactive States along a Reaction Coordinate in Enzymes.

Enzymes play a pivotal role in all biological systems. These biomachines are the most effective catalysts known, dramatically enhancing the rate of reactions by more than 10 orders of magnitude relative to the uncatalyzed reactions in solution. Predicting the correct, mechanistically appropriate binding modes for substrate and product, as well as all reaction intermediates and transition states, along a reaction pathway is immensely challenging and remains an unsolved problem. In the present work, we developed an effective methodology for identifying probable binding modes of multiple ligand states along a reaction coordinate in an enzyme active site. The program is called EnzyDock and is a CHARMM-based multistate consensus docking program that includes a series of protocols to predict the chemically relevant orientation of substrate, reaction intermediates, transition states, product, and inhibitors. EnzyDock is based on simulated annealing molecular dynamics and Monte Carlo sampling and allows ligand, as well as enzyme side-chain and backbone flexibility. The program can employ many user-defined constraints and restraints and classical force field potentials, as well as a range of hybrid quantum mechanics-molecular mechanics potentials. Herein, we apply EnzyDock to several different kinds of problems. First, we study two terpene synthase reactions, namely bornyl diphosphate synthase and the bacterial diterpene synthase CotB2. Second, we use EnzyDock to predict reaction coordinate states in a pair of Diels-Alder reactions in the enzymes spirotetronate AbyU and LepI. Third, we study a couple of racemases: the cofactor-dependent serine racemase and the cofactor independent proline racemase. Finally, we study several cases of covalent docking involving the Michael addition reaction. For all systems we predict binding modes that are consistent with available experimental observations, as well as with theoretical modeling studies from the literature. EnzyDock provides a platform for generating mechanistic insight into enzyme reactions, useful and reliable starting points for in-depth multiscale modeling projects, and rational design of noncovalent and covalent enzyme inhibitors.

Discovering Monoterpene Catalysis Inside Nanocapsules with Multiscale Modeling and Experiments.

Large-scale production of natural products, such as terpenes, presents a significant scientific and technological challenge. One promising approach to tackle this problem is chemical synthesis inside nanocapsules, although enzyme-like control of such chemistry has not yet been achieved. In order to better understand the complex chemistry inside nanocapsules, we design a multiscale nanoreactor simulation approach. The nanoreactor simulation protocol consists of hybrid quantum mechanics-molecular mechanics-based high temperature Langevin molecular dynamics simulations. Using this approach we model the tail-to-head formation of monoterpenes inside a resorcin[4]arene-based capsule (capsule I). We provide a rationale for the experimentally observed kinetics of monoterpene product formation and product distribution using capsule I, and we explain why additional stable monoterpenes, like camphene, are not observed. On the basis of the in-capsule I simulations, and mechanistic insights, we propose that feeding the capsule with pinene can yield camphene, and this proposal is verified experimentally. This suggests that the capsule may direct the dynamic reaction cascades by virtue of π-cation interactions.

Oscillatory Active-site Motions Correlate with Kinetic Isotope Effects in Formate Dehydrogenase.

Thermal motions of enzymes have been invoked to explain the temperature dependence of kinetic isotope effects (KIE) in enzyme-catalyzed hydride transfers. Formate dehydrogenase (FDH) from Candida boidinii exhibits a temperature independent KIE that becomes temperature dependent upon mutation of hydrophobic residues in the active site. Ternary complexes of FDH that mimic the transition state structure allow investigation of how these mutations influence active-site dynamics. A combination of X-ray crystallography, two-dimensional infrared (2D IR) spectroscopy, and molecular dynamic simulations characterize the structure and dynamics of the active site. FDH exhibits oscillatory frequency fluctuations on the picosecond timescale, and the amplitude of these fluctuations correlates with the temperature dependence of the KIE. Both the kinetic and dynamic phenomena can be reproduced computationally. These results provide experimental evidence for a connection between the temperature dependence of KIEs and motions of the active site in an enzyme-catalyzed reaction consistent with activated tunneling models of the hydride transfer reaction.

Understanding the Role of Minor Molybdenum Doping in LiNi0.5Co0.2Mn0.3O2 Electrodes: from Structural and Surface Analyses and Theoretical Modeling to Practical Electrochemical Cells.

Doping LiNi0.5Co0.2Mn0.3O2 (NCM523) cathode material by small amount of Mo6+ ions, around 1 mol %, affects pronouncedly its structure, surface properties, and electronic and electrochemical behavior. Cathodes comprising Mo6+-doped NCM523 exhibited in Li cells higher specific capacities, higher rate capabilities, lower capacity fading, and lower charge-transfer resistance that relates to a more stable electrode/solution interface due to doping. This, in turn, is ascribed to the fact that the Mo6+ ions tend to concentrate more at the surface, as a result of a synthesis that always includes a necessary calcination, high-temperature stage. This phenomenon of the Mo dopant segregation at the surface in NCM523 material was discovered in the present work for the first time. It appears that Mo doping reduces the reactivity of the Ni-rich NCM cathode materials toward the standard electrolyte solutions of Li-ion batteries. Using density functional theory (DFT) calculations, we showed that Mo6+ ions are preferably incorporated at Ni sites and that the doping increases the amount of Ni2+ ions at the expense of Ni3+ ions, due to charge compensation, in accord with X-ray absorption fine structure (XAFS) spectroscopy measurements. Furthermore, DFT calculations predicted Ni-O bond length distributions in good agreement with the XAFS results, supporting a model of partial substitution of Ni sites by molybdenum.

A promising drug candidate for the treatment of glaucoma based on a P2Y6-receptor agonist.

Extracellular nucleotides can regulate the production/drainage of the aqueous humor via activation of P2 receptors, thus affecting the intraocular pressure (IOP). We evaluated 5-OMe-UDP(α-B), 1A, a potent P2Y6-receptor agonist, for reducing IOP and treating glaucoma. Cell viability in the presence of 1A was measured using [3-(4, 5-dimethyl-thiazol-2-yl) 2, 5-diphenyl-tetrazolium bromide] (MTT) assay in rabbit NPE ciliary non-pigmented and corneal epithelial cells, human retinoblastoma, and liver Huh7 cells. The effect of 1A on IOP was determined in acute glaucomatous rabbit hyaluronate model and phenol-induced chronic glaucomatous rabbit model. The origin of activity of 1A was investigated by generation of a homology model of hP2Y6-R and docking studies. 1A did not exert cytotoxic effects up to 100 mM vs. trusopt and timolol in MTT assay in ocular and liver cells. In normotensive rabbits, 100 µM 1A vs. xalatan, trusopt, and pilocarpine reduced IOP by 45 vs. 20-30%, respectively. In the phenol animal model, 1A (100 µM) showed reduction of IOP by 40 and 20%, following early and late administration, respectively. Docking results suggest that the high activity and selectivity of 1A is due to intramolecular interaction between Pα-BH3 and C5-OMe which positions 1A in a most favorable site inside the receptor. P2Y6-receptor agonist 1A effectively and safely reduces IOP in normotense, acute, and chronic glaucomatous rabbits, and hence may be suggested as a novel approach for the treatment of glaucoma.

Comment on "Substrate Folding Modes in Trichodiene Synthase: A Determinant of Chemo- and Stereoselectivity".

Wang et al. recently reported an in silico study of the trichodiene synthase (TDS) conversion of farnesyl diphosphate (FPP) to trichodiene (TD) (Wang et al., ACS Catal. 2017, 7, 5841-5846). Although the methods and level of theory used in that work are nearly identical to our own recent work on this system (Dixit et al., ACS Catal. 2017, 7, 812-818), Wang et al. reach rather different conclusions. The authors claimed to obtain a "very credible" mechanism for the biosynthesis of TD and optimized the optimal folding mode of FPP in the 1,6-ring closure in TDS. However, the folding mode of the FPP substrate that was presented contradicts well-established NMR and mass spectrometry data. Moreover, the authors make numerous incorrect statements regarding our earlier work.

A surprising substituent effect in corroles on the electrochemical activation of oxygen reduction.

Four ß-pyrrole-substituted cobalt(iii) corroles were studied as electrocatalysts for the oxygen reduction reaction. The results disclose high dependence of the corrole's performance on its substituents, but once adsorbed on a high surface area carbon, this effect vanishes, resulting in a better catalytic performance than most well-defined molecular electrocatalysts for this reaction.