Search

Liquid electrolytes are one of the most important components of Li-ion batteries, which are a critical technology of the modern world. However, we still lack the computational tools required to accurately calculate key properties of these materials (viscosity and ionic diffusivity) from first principles necessary to support improved designs. In this work, we report a machine learning-based force field for liquid electrolyte simulations, which bridges the gap between the accuracy of range-separated hybrid density functional theory and the efficiency of classical force fields. Predictions of material properties made with this force field are quantitatively accurate compared to experimental data. Our model uses the QRNN deep neural network architecture, which includes both long-range interactions and global charge equilibration. The training data set is composed solely of non-periodic density functional theory (DFT), allowing the practical use of an accurate theory (here, ωB97X-D3BJ/def2-TZVPD), which would be prohibitively expensive for generating large data sets with periodic DFT. In this report, we focus on seven common carbonates and LiPF 6 , but this methodology has very few assumptions and can be readily applied to any liquid electrolyte system. This provides a promising path forward for large-scale atomistic modeling of many important battery chemistries.

Pyridoxal 5'-phosphate (PLP)-dependent enzymes have been extensively studied for their ability to fine-tune PLP cofactor electronics to promote a wide array of chemistries. Neutron crystallography offers a straightforward approach to studying the electronic states of PLP and the electrostatics of enzyme active sites, responsible for the reaction specificities, by enabling direct visualization of hydrogen atom positions. Here we report a room-temperature joint X-ray/neutron structure of aspartate aminotransferase (AAT) with pyridoxamine 5'-phosphate (PMP), the cofactor product of the first half reaction catalyzed by the enzyme. Between PMP N SB and catalytic Lys258 Nζ amino groups an equally shared deuterium is observed in an apparent low-barrier hydrogen bond (LBHB). Density functional theory calculations were performed to provide further evidence of this LBHB interaction. The structural arrangement and the juxtaposition of PMP and Lys258, facilitated by the LBHB, suggests active site preorganization for the incoming ketoacid substrate that initiates the second half-reaction., Competing Interests: There are no conflicts of interest to declare., (This journal is © The Royal Society of Chemistry.)

To address some of the inherent challenges in modeling metalloenzymes, we here report an extension to the functional form of the OPLS3e force field to include terms adopted from the ligand field molecular mechanics (LFMM) model, including the angular overlap and Morse potential terms. The integration of these terms with OPLS3e, herein referred to as OPLS3e+M, improves the description of metal-ligand interactions and provides accurate relative binding energies and geometric preferences of transition-metal complexes by training to gas-phase density functional theory (DFT) energies. For [Cu(H 2 O) 4 ] 2+ , OPLS3e+M significantly improves H 2 O binding energies and the geometric preference of the tetra-aqua Cu 2+ complex. In addition, we conduct free-energy perturbation calculations on two pharmaceutically relevant metalloenzyme targets, which include chemical modifications at varying proximity to the binding-site metals, including changes to the metal-binding moiety of the ligand itself. The extensions made to OPLS3e lead to accurate predicted relative binding free energies for these series (mean unsigned error of 1.29 kcal mol -1 ). Our results provide evidence that integration of the LFMM model with OPLS3e can be utilized to predict thermodynamic quantities for such systems near chemical accuracy. With these improvements, we anticipate that robust free-energy perturbation calculations can be employed to accelerate the drug development efforts for metalloenzyme targets.

Protein dynamics on various time scales from femtoseconds to milliseconds impacts biological function by driving proteins to conformations conducive to ligand binding and creating functional states in enzyme catalysis. Neutron vibrational spectroscopy carried out by measuring inelastic neutron scattering from protein molecules in combination with molecular simulations has the unique ability of detecting and visualizing changes in the picosecond protein vibrational dynamics due to ligand binding. Here we present neutron vibrational spectra of a homodimeric pyridoxal 5'-phosphate-dependent enzyme, aspartate aminotransferase, obtained from the open internal aldimine and closed external aldimine conformational states. We observe that in the external aldimine state the protein structure stiffens relative to the internal aldimine state, indicating rigidified vibrational dynamics on the picosecond time scale in the low-frequency regime of 5-50 cm -1 . Our molecular dynamics simulations indicate substantial changes in the picosecond dynamics of the enzyme secondary structure elements upon substrate binding, with the largest contributions from just two helices and the β-sheet., Competing Interests: The authors declare no competing financial interest., (Copyright © 2020 American Chemical Society.)

Enzyme catalysis is the primary activity in energy and information metabolism and enzyme cofactors are key to the catalytic ability of most enzymes. Pyridoxal 5'-phosphate (PLP) cofactor, derived from Vitamin B6, is widely distributed in nature and has significant latitude in catalytic diversity. X-ray crystallography has revealed the structures of diverse PLP dependent enzymes from multiple families. But these structures are incomplete, lacking the positions of protons essential for understanding enzymatic mechanisms. Here, we review the diversity of PLP and discuss the use of neutron crystallography and joint X-ray/neutron refinement of Fold Type I aspartate aminotransferase to visualize the positions of protons in both the internal and external aldimine forms. Strategies used to prepare extremely large crystals required for neutron diffraction and the approach to data refinement including the PLP cofactor are discussed. The observed positions of protons, including one located in a previously unknown low-barrier hydrogen bond, have been used to create more accurate models for computational analysis. The results revealed a new mechanism for the transaminase reaction where hyperconjugation is key to reducing the energy barrier which finally provides a clear explanation of the Dunathan alignment., (© 2020 Elsevier Inc. All rights reserved.)

[This retracts the article DOI: 10.1107/S2052252519005761.].

Mycobacterium tuberculosis antigen 85 (Ag85) enzymes catalyze the transfer of mycolic acid (MA) from trehalose monomycolate to produce the mycolyl arabinogalactan (mAG) or trehalose dimycolate (TDM). These lipids define the protective mycomembrane of mycobacteria. The current model of substrate binding within the active sites of Ag85s for the production of TDM is not sterically and geometrically feasible; additionally, this model does not account for the production of mAG. Furthermore, this model does not address how Ag85s limit the hydrolysis of the acyl-enzyme intermediate while catalyzing acyl transfer. To inform an updated model, we obtained an Ag85 acyl-enzyme intermediate structure that resembles the mycolated form. Here, we present a 1.45-Å X-ray crystal structure of M. tuberculosis Ag85C covalently modified by tetrahydrolipstatin (THL), an esterase inhibitor that suppresses M. tuberculosis growth and mimics structural attributes of MAs. The mode of covalent inhibition differs from that observed in the reversible inhibition of the human fatty-acid synthase by THL. Similarities between the Ag85-THL structure and previously determined Ag85C structures suggest that the enzyme undergoes structural changes upon acylation, and positioning of the peptidyl arm of THL limits hydrolysis of the acyl-enzyme adduct. Molecular dynamics simulations of the modeled mycolated-enzyme form corroborate the structural analysis. From these findings, we propose an alternative arrangement of substrates that rectifies issues with the previous model and suggest a direct role for the β-hydroxy of MA in the second half-reaction of Ag85 catalysis. This information affords the visualization of a complete mycolyltransferase catalytic cycle.

The Mycobacterium tuberculosis rv3802c gene encodes an essential enzyme with thioesterase and phospholipase A activity. Overexpression of Rv3802 orthologs in Mycobacterium smegmatis and Corynebacterium glutamicum increases mycolate content and decreases glycerophospholipids. Although a role in modulating the lipid composition of the unique mycomembrane has been proposed, the true biological function of Rv3802 remains uncertain. In this study, we present the first M. tuberculosis Rv3802 X-ray crystal structure, solved to 1.7 Å resolution. On the basis of the binding of PEG molecules to Rv3802, we identified its lipid-binding site and the structural basis for phosphatidyl-based substrate binding and phospholipase A activity. We found that movement of the α8-helix affords lipid binding and is required for catalytic turnover through covalent tethering. We gained insights into the mechanism of acyl hydrolysis by observing differing arrangements of PEG and water molecules within the active site. This study provides structural insights into biological function and facilitates future structure-based drug design toward Rv3802.

Enzymes dependent on pyridoxal 5'-phosphate (PLP, the active form of vitamin B 6 ) perform a myriad of diverse chemical transformations. They promote various reactions by modulating the electronic states of PLP through weak interactions in the active site. Neutron crystallography has the unique ability of visualizing the nuclear positions of hydrogen atoms in macromolecules. Here we present a room-temperature neutron structure of a homodimeric PLP-dependent enzyme, aspartate aminotransferase, which was reacted in situ with α-methylaspartate. In one monomer, the PLP remained as an internal aldimine with a deprotonated Schiff base. In the second monomer, the external aldimine formed with the substrate analog. We observe a deuterium equidistant between the Schiff base and the C-terminal carboxylate of the substrate, a position indicative of a low-barrier hydrogen bond. Quantum chemical calculations and a low-pH room-temperature X-ray structure provide insight into the physical phenomena that control the electronic modulation in aspartate aminotransferase.Pyridoxal 5'-phosphate (PLP) is a ubiquitous co factor for diverse enzymes, among them aspartate aminotransferase. Here the authors use neutron crystallography, which allows the visualization of the positions of hydrogen atoms, and computation to characterize the catalytic mechanism of the enzyme.

Previous studies identified ebselen as a potent in vitro and in vivo inhibitor of the Mycobacterium tuberculosis (Mtb) antigen 85 (Ag85) complex, comprising three homologous enzymes required for the biosynthesis of the mycobacterial cell wall. In this study, the Mtb Ag85C enzyme was cocrystallized with azido and adamantyl ebselen derivatives, resulting in two crystallographic structures of 2.01 and 1.30 Å resolution, respectively. Both structures displayed the anticipated covalent modification of the solvent accessible, noncatalytic Cys209 residue forming a selenenylsulfide bond. Continuous difference density for both thiol modifiers allowed for the assessment of interactions that influence ebselen binding and inhibitor orientation that were unobserved in previous Ag85C ebselen structures. The k inact /K I values for ebselen, adamantyl ebselen, and azido ebselen support the importance of observed constructive chemical interactions with Arg239 for increased in vitro efficacy toward Ag85C. To better understand the in vitro kinetic properties of these ebselen derivatives, the energetics of specific protein-inhibitor interactions and relative reaction free energies were calculated for ebselen and both derivatives using density functional theory. These studies further support the different in vitro properties of ebselen and two select ebselen derivatives from our previously published ebselen library with respect to kinetics and protein-inhibitor interactions. In both structures, the α9 helix was displaced farther from the enzyme active site than the previous Ag85C ebselen structure, resulting in the restructuring of a connecting loop and imparting a conformational change to residues believed to play a role in substrate binding specific to Ag85C. These notable structural changes directly affect protein stability, reducing the overall melting temperature by up to 14.5 °C, resulting in the unfolding of protein at physiological temperatures. Additionally, this structural rearrangement due to covalent allosteric modification creates a sizable solvent network that encompasses the active site and extends to the modified Cys209 residue. In all, this study outlines factors that influence enzyme inhibition by ebselen and its derivatives while further highlighting the effects of the covalent modification of Cys209 by said inhibitors on the structure and stability of Ag85C. Furthermore, the results suggest a strategy for developing new classes of Ag85 inhibitors with increased specificity and potency.

Pyridoxal 5'-phosphate (PLP) is a fundamental, multifunctional enzyme cofactor used to catalyze a wide variety of chemical reactions involved in amino acid metabolism. PLP-dependent enzymes optimize specific chemical reactions by modulating the electronic states of PLP through distinct active site environments. In aspartate aminotransferase (AAT), an extended hydrogen bond network is coupled to the pyridinyl nitrogen of the PLP, influencing the electrophilicity of the cofactor. This network, which involves residues Asp-222, His-143, Thr-139, His-189, and structural waters, is located at the edge of PLP opposite the reactive Schiff base. We demonstrate that this hydrogen bond network directly influences the protonation state of the pyridine nitrogen of PLP, which affects the rates of catalysis. We analyzed perturbations caused by single- and double-mutant variants using steady-state kinetics, high resolution X-ray crystallography, and quantum chemical calculations. Protonation of the pyridinyl nitrogen to form a pyridinium cation induces electronic delocalization in the PLP, which correlates with the enhancement in catalytic rate in AAT. Thus, PLP activation is controlled by the proximity of the pyridinyl nitrogen to the hydrogen bond microenvironment. Quantum chemical calculations indicate that Asp-222, which is directly coupled to the pyridinyl nitrogen, increases the p K a of the pyridine nitrogen and stabilizes the pyridinium cation. His-143 and His-189 also increase the p K a of the pyridine nitrogen but, more significantly, influence the position of the proton that resides between Asp-222 and the pyridinyl nitrogen. These findings indicate that the second shell residues directly enhance the rate of catalysis in AAT., (© 2017 by The American Society for Biochemistry and Molecular Biology, Inc.)

Neutron crystallography provides direct visual evidence of the atomic positions of deuterium-exchanged H atoms, enabling the accurate determination of the protonation/deuteration state of hydrated biomolecules. Comparison of two neutron structures of hemoglobins, human deoxyhemoglobin (T state) and equine cyanomethemoglobin (R state), offers a direct observation of histidine residues that are likely to contribute to the Bohr effect. Previous studies have shown that the T-state N-terminal and C-terminal salt bridges appear to have a partial instead of a primary overall contribution. Four conserved histidine residues [αHis72(EF1), αHis103(G10), αHis89(FG1), αHis112(G19) and βHis97(FG4)] can become protonated/deuterated from the R to the T state, while two histidine residues [αHis20(B1) and βHis117(G19)] can lose a proton/deuteron. αHis103(G10), located in the α1:β1 dimer interface, appears to be a Bohr group that undergoes structural changes: in the R state it is singly protonated/deuterated and hydrogen-bonded through a water network to βAsn108(G10) and in the T state it is doubly protonated/deuterated with the network uncoupled. The very long-term H/D exchange of the amide protons identifies regions that are accessible to exchange as well as regions that are impermeable to exchange. The liganded relaxed state (R state) has comparable levels of exchange (17.1% non-exchanged) compared with the deoxy tense state (T state; 11.8% non-exchanged). Interestingly, the regions of non-exchanged protons shift from the tetramer interfaces in the T-state interface (α1:β2 and α2:β1) to the cores of the individual monomers and to the dimer interfaces (α1:β1 and α2:β2) in the R state. The comparison of regions of stability in the two states allows a visualization of the conservation of fold energy necessary for ligand binding and release.