Search

A new superacid, H(CB11H6X6) (where X = chlorine or bromine), whose conjugate base is the exceptionally inert CB11H6X6- carborane anion, separates Bronsted acidity from oxidizing capacity and anion nucleophilicity in a manner not previously achieved. Reaction of this superacid with C60 gives HC60+ as a stable ion in solution and in the solid state. In a separate experiment, an oxidant was developed such that the long-sought C60.+ ion can be synthesized in solution. The preparation of these two fullerene carbocations is a notable departure from the prevalent chemistry of C60, which is dominated by the formation of anions or the addition of nucleophiles. The H(CB11H6X6) superacid overcomes the major limitations of presently known superacids and has potentially wide application.

A nearly universal component of NMR crystallography is the ranking of candidate structures based on how well their first-principles-predicted NMR parameters align with the results of solid-state NMR experiments. Here, a novel approach for assigning probabilities to candidate models is proposed that quantifies the likelihood that each model is the correct experimental structure. This method employs hierarchical Bayesian inference and leverages explicit prior probabilities derived from a uniform distribution of potential candidate structures with respect to chi-squared values. The resulting uniform chi-squared (UC) model provides a more cautious estimate of candidate probabilities compared to previous approaches, assigning decreased likelihood to the best-fit structure and increased likelihood to alternate candidates. Although developed here within the context of NMR crystallography, the UC model represents a general method for assigning likelihoods based on chi-squared goodness-of-fit assessments.

Ab initio chemical shift prediction plays a central role in nuclear magnetic resonance (NMR) crystallography, and the accuracy with which chemical shifts can be predicted relative to experiment impacts the confidence with which structures can be assigned. For organic crystals, periodic density functional theory calculations with the gauge-including projector augmented wave (GIPAW) approximation and the PBE functional are widely used at present. Many previous studies have examined how using more advanced density functionals can increase the accuracy of predicted chemical shifts relative to experiment, but nearly all of those studies employed crystal structures that were optimized with generalized-gradient approximation (GGA) functionals. Here, we investigate how the accuracy of the predicted chemical shifts in organic crystals is affected by replacing GGA-level PBE-D3(BJ) crystal geometries with more accurate hybrid functional PBE0-D3(BJ) ones. Based on benchmark data sets containing 132 13 C and 35 15 N chemical shifts, plus case studies on testosterone, acetaminophen, and phenobarbital, we find that switching from GGA-level geometries and chemical shifts to hybrid-functional ones reduces 13 C and 15 N chemical shift errors by ∼40-60% versus experiment. However, most of the improvement stems from the use of the hybrid functional for the chemical shift calculations, rather than from the refined geometries. In addition, even with the improved geometries, we find that double-hybrid functionals still do not systematically increase chemical shift agreement with experiment beyond what hybrid functionals provide. In the end, these results suggest that the combination of GGA-level crystal structures and hybrid-functional chemical shifts represents a particularly cost-effective combination for NMR crystallography in organic systems.

Solution NMR is typically applied to biological systems with molecular weights < 40 kDa whereas magic-angle-spinning (MAS) solid-state NMR traditionally targets very large, oligomeric proteins and complexes exceeding 500 kDa in mass, including fibrils and crystalline protein preparations. Here, we propose that the gap between these size regimes can be filled by the approach presented that enables investigation of large, soluble and fully protonated proteins in the range of 40-140 kDa. As a key step, ultracentrifugation produces a highly concentrated, gel-like state, resembling a dense phase in spontaneous liquid-liquid phase separation (LLPS). By means of three examples, a Sulfolobus acidocaldarius bifurcating electron transfer flavoprotein (SaETF), tryptophan synthases from Salmonella typhimurium (StTS) and their dimeric β-subunits from Pyrococcus furiosus (PfTrpB), we show that such samples yield well-resolved proton-detected 2D and 3D NMR spectra at 100 kHz MAS without heterogeneous broadening, similar to diluted liquids. Herein, we provide practical guidance on centrifugation conditions and tools, sample behavior, and line widths expected. We demonstrate that the observed chemical shifts correspond to those obtained from µM/low mM solutions or crystalline samples, indicating structural integrity. Nitrogen line widths as low as 20-30 Hz are observed. The presented approach is advantageous for proteins or nucleic acids that cannot be deuterated due to the expression system used, or where relevant protons cannot be re-incorporated after expression in deuterated medium, and it circumvents crystallization. Importantly, it allows the use of low-glycerol buffers in dynamic nuclear polarization (DNP) NMR of proteins as demonstrated with the cyanobacterial phytochrome Cph1., (© 2024. The Author(s).)

TensorView for MATLAB is a GUI-based visualization tool for depicting second-rank Cartesian tensors as surfaces on three-dimensional molecular models. Both ellipsoid and ovaloid tensor display formats are supported, and the software allows for easy conversion of Euler angles from common rotation schemes (active, passive, ZXZ, and ZYZ conventions) with visual feedback. In addition, the software displays all four orientation-equivalent Euler angle solutions for the placement of a single tensor in the molecular frame and can report relative orientations of two tensors with all 16 orientation-equivalent Euler angle sets that relate them. The salient relations are derived and illustrated through several examples. TensorView for MATLAB expands and complements the earlier implementation of TensorView within the Mathematica programming environment and can be run without a MATLAB license. TensorView for MATLAB is available through github at https://github.com/LeoSvenningsson/TensorViewforMatlab, and can also be accessed directly via the NMRbox resource., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper, (Copyright © 2022 Elsevier Inc. All rights reserved.)

Photomechanical molecular crystals have garnered attention for their ability to transform light into mechanical work, but difficulties in characterizing the structural changes and mechanical responses experimentally have hindered the development of practical organic crystal engines. This study proposes a new computational framework for predicting the solid-state crystal-to-crystal photochemical transformations entirely from first principles, and it establishes a photomechanical engine cycle that quantifies the anisotropic mechanical performance resulting from the transformation. The approach relies on crystal structure prediction, solid-state topochemical principles, and high-quality electronic structure methods. After validating the framework on the well-studied [4 + 4] cycloadditions in 9-methyl anthracene and 9- tert -butyl anthracene ester, the experimentally-unknown solid-state transformation of 9-carboxylic acid anthracene is predicted for the first time. The results illustrate how the mechanical work is done by relaxation of the crystal lattice to accommodate the photoproduct, rather than by the photochemistry itself. The large ∼10 7 J m -3 work densities computed for all three systems highlight the promise of photomechanical crystal engines. This study demonstrates the importance of crystal packing in determining molecular crystal engine performance and provides tools and insights to design improved materials in silico ., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

The reaction of W(NAr)( 13 C 4 H 8 )(OSiPh 3 ) 2 ( 1 ) (NAr = 2,6-diisopropylphenylimido) with silica partially dehydroxylated at 700 °C (SiO 2-700 ) is highly dependent on the reaction conditions. The primary product of this reaction is W(NAr)( 13 C 4 H 8 )(OSiPh 3 )(OSi(O-) 3 ) ( 2 ) when the reaction is carried out in the dark. Grafting 1 onto SiO 2-700 in ambient lab light results in the formation of 2 , W(NAr)( 13 CH 2 13 CH 2 )(OSiPh 3 )(OSi(O-) 3 ) ( 4 ), and one isomer of square-pyramidal W(NAr)( 13 CH 2 13 CH( 13 Me) 13 CH 2 )(OSiPh 3 )(OSi(O-) 3 ) ( 3 ). Heating 2 to 85 °C for 6 h results in the formation of 3 , 4 , W(NAr)( 13 CH( 13 Me) 13 CH 2 13 CH 2 )(OSiPh 3 )(OSi(O-) 3 ) ( 5 ), and W(NAr)(( 13 CH 2 ) 2 13 CH( 13 Me)( 13 CH 2 ) 2 )(OSiPh 3 )(OSi(O-) 3 ) ( 6 ). Photolysis of 2 with blue LEDs (λ max = 450 nm) produces 4, both isomers of 3, 5 , and free ethylene. In the presence of excess ethylene and blue LED irradiation at 85 °C, 1 /SiO 2-700 catalyzes the direct conversion of ethylene to propylene.

The regulation of the synthesis of L-tryptophan (L-Trp) in enteric bacteria begins at the level of gene expression where the cellular concentration of L-Trp tightly controls expression of the five enzymes of the Trp operon responsible for the synthesis of L-Trp. Two of these enzymes, trpA and trpB, form an αββα bienzyme complex, designated as tryptophan synthase (TS). TS carries out the last two enzymatic processes comprising the synthesis of L-Trp. The TS α-subunits catalyze the cleavage of 3-indole D-glyceraldehyde 3'-phosphate to indole and D-glyceraldehyde 3-phosphate; the pyridoxal phosphate-requiring β-subunits catalyze a nine-step reaction sequence to replace the L-Ser hydroxyl by indole giving L-Trp and a water molecule. Within αβ dimeric units of the αββα bienzyme complex, the common intermediate indole is channeled from the α site to the β site via an interconnecting 25 Å-long tunnel. The TS system provides an unusual example of allosteric control wherein the structures of the nine different covalent intermediates along the β-reaction catalytic path and substrate binding to the α-site provide the allosteric triggers for switching the αββα system between the open (T) and closed (R) allosteric states. This triggering provides a linkage that couples the allosteric conformational coordinate to the covalent chemical reaction coordinates at the α- and β-sites. This coupling drives the α- and β-sites between T and R conformations to achieve regulation of substrate binding and/or product release, modulation of the α- and β-site catalytic activities, prevention of indole escape from the confines of the active sites and the interconnecting tunnel, and synchronization of the α- and β-site catalytic activities. Here we review recent advances in the understanding of the relationships between structure, function, and allosteric regulation of the complex found in Salmonella typhimurium., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Ghosh, Hilario, Chang, Mueller and Dunn.)

A self-assembled Fe II 4 L 6 cage was synthesized with 12 internal amines in the cavity. The cage forms as the dodeca-ammonium salt, despite the cage carrying an overall 8+ charge at the metal centers, extracting protons from displaced water in the reaction. Despite this, the basicity of the internal amines is lower than their counterparts in free solution. The 12 amines have a sliding scale of basicity, with a ≈6 pK a unit difference between the first and last protons to be removed. This moderation of side-chain basicity in an active site is a hallmark of enzymatic catalysis., (© 2022 Wiley-VCH GmbH.)

Antibiotic resistance is a continually growing challenge in the treatment of various bacterial infections worldwide. New drugs and new drug targets are necessary to curb the threat of infectious diseases caused by multidrug-resistant pathogens. The tryptophan biosynthesis pathway is essential for bacterial growth but is absent in higher animals and humans. Drugs that can inhibit the bacterial biosynthesis of tryptophan offer a new class of antibiotics. In this work, we combined a structure-based strategy using in silico docking screening and molecular dynamics (MD) simulations to identify compounds targeting the α subunit of tryptophan synthase with experimental methods involving the whole-cell minimum inhibitory concentration (MIC) test, solution state NMR, and crystallography to confirm the inhibition of L-tryptophan biosynthesis. Screening 1,800 compounds from the National Cancer Institute Diversity Set I against α subunit revealed 28 compounds for experimental validation; four of the 28 hit compounds showed promising activity in MIC testing. We performed solution state NMR experiments to demonstrate that a one successful inhibitor, 3-amino-3-imino-2-phenyldiazenylpropanamide (Compound 1) binds to the α subunit. We also report a crystal structure of Salmonella enterica serotype Typhimurium tryptophan synthase in complex with Compound 1 which revealed a binding site at the αβ interface of the dimeric enzyme. MD simulations were carried out to examine two binding sites for the compound. Our results show that this small molecule inhibitor could be a promising lead for future drug development., (© 2021 The Protein Society.)

NMR chemical shifts provide detailed information on the chemical properties of molecules, thereby complementing structural data from techniques like X-ray crystallography and electron microscopy. Detailed analysis of protein NMR data, however, often hinges on comprehensive, site-specific assignment of backbone resonances, which becomes a bottleneck for molecular weights beyond 40 to 45 kDa. Here, we show that assignments for the (2x)72-kDa protein tryptophan synthase (665 amino acids per asymmetric unit) can be achieved via higher-dimensional, proton-detected, solid-state NMR using a single, 1-mg, uniformly labeled, microcrystalline sample. This framework grants access to atom-specific characterization of chemical properties and relaxation for the backbone and side chains, including those residues important for the catalytic turnover. Combined with first-principles calculations, the chemical shifts in the β-subunit active site suggest a connection between active-site chemistry, the electrostatic environment, and catalytically important dynamics of the portal to the β-subunit from solution., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

NMR-assisted crystallography-the integrated application of solid-state NMR, X-ray crystallography, and first-principles computational chemistry-holds significant promise for mechanistic enzymology: by providing atomic-resolution characterization of stable intermediates in enzyme active sites, including hydrogen atom locations and tautomeric equilibria, NMR crystallography offers insight into both structure and chemical dynamics. Here, this integrated approach is used to characterize the tryptophan synthase α-aminoacrylate intermediate, a defining species for pyridoxal-5'-phosphate-dependent enzymes that catalyze β-elimination and replacement reactions. For this intermediate, NMR-assisted crystallography is able to identify the protonation states of the ionizable sites on the cofactor, substrate, and catalytic side chains as well as the location and orientation of crystallographic waters within the active site. Most notable is the water molecule immediately adjacent to the substrate β-carbon, which serves as a hydrogen bond donor to the ε-amino group of the acid-base catalytic residue βLys87. From this analysis, a detailed three-dimensional picture of structure and reactivity emerges, highlighting the fate of the L-serine hydroxyl leaving group and the reaction pathway back to the preceding transition state. Reaction of the α-aminoacrylate intermediate with benzimidazole, an isostere of the natural substrate indole, shows benzimidazole bound in the active site and poised for, but unable to initiate, the subsequent bond formation step. When modeled into the benzimidazole position, indole is positioned with C3 in contact with the α-aminoacrylate C β and aligned for nucleophilic attack. Here, the chemically detailed, three-dimensional structure from NMR-assisted crystallography is key to understanding why benzimidazole does not react, while indole does., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

Photomechanical molecular crystals that expand under illumination could potentially be used as photon-powered actuators. In this study, we find that the use of high-quality single crystals of 9-methylanthracene (9MA) leads to more homogeneous reaction kinetics than that previously seen for polycrystalline samples, presumably due to a lower concentration of defects. Furthermore, simultaneous observation of absorbance and shape changes in single crystals revealed that the dimensional change mirrors the reaction progress, resulting in a smooth expansion of 7 % along the c-axis that is linearly correlated with reaction progress. The same expansion dynamics are highly reproducible across different single crystal samples. Organic single crystals exhibit well-defined linear expansions during 100 % photoconversion, suggesting that this class of solid-state phase change material could be used for actuation., (© 2021 Wiley-VCH GmbH.)

The tryptophan synthase (TS) bienzyme complexes found in bacteria, yeasts, and molds are pyridoxal 5'-phosphate (PLP)-requiring enzymes that synthesize l-Trp. In the TS catalytic cycle, switching between the open and closed states of the α- and β-subunits via allosteric interactions is key to the efficient conversion of 3-indole-d-glycerol-3'-phosphate and l-Ser to l-Trp. In this process, the roles played by β-site residues proximal to the PLP cofactor have not yet been fully established. βGln114 is one such residue. To explore the roles played by βQ114, we conducted a detailed investigation of the βQ114A mutation on the structure and function of tryptophan synthase. Initial steady-state kinetic and static ultraviolet-visible spectroscopic analyses showed the Q to A mutation impairs catalytic activity and alters the stabilities of intermediates in the β-reaction. Therefore, we conducted X-ray structural and solid-state nuclear magnetic resonance spectroscopic studies to compare the wild-type and βQ114A mutant enzymes. These comparisons establish that the protein structural changes are limited to the Gln to Ala replacement, the loss of hydrogen bonds among the side chains of βGln114, βAsn145, and βArg148, and the inclusion of waters in the cavity created by substitution of the smaller Ala side chain. Because the conformations of the open and closed allosteric states are not changed by the mutation, we hypothesize that the altered properties arise from the lost hydrogen bonds that alter the relative stabilities of the open (β T state) and closed (β R state) conformations of the β-subunit and consequently alter the distribution of intermediates along the β-subunit catalytic path.

Backbone chemical shift assignments for the Toho-1 β-lactamase (263 amino acids, 28.9 kDa) are reported based on triple resonance solution-state NMR experiments performed on a uniformly 2 H, 13 C, 15 N-labeled sample. These assignments allow for subsequent site-specific characterization at the chemical, structural, and dynamical levels. At the chemical level, titration with the non-β-lactam β-lactamase inhibitor avibactam is found to give chemical shift perturbations indicative of tight covalent binding that allow for mapping of the inhibitor binding site. At the structural level, protein secondary structure is predicted based on the backbone chemical shifts and protein residue sequence using TALOS-N and found to agree well with structural characterization from X-ray crystallography. At the dynamical level, model-free analysis of 15 N relaxation data at a single field of 16.4 T reveals well-ordered structures for the ligand-free and avibactam-bound enzymes with generalized order parameters of ~ 0.85. Complementary relaxation dispersion experiments indicate that there is an escalation in motions on the millisecond timescale in the vicinity of the active site upon substrate binding. The combination of high rigidity on short timescales and active site flexibility on longer timescales is consistent with hypotheses for achieving both high catalytic efficiency and broad substrate specificity: the induced active site dynamics allows variously sized substrates to be accommodated and increases the probability that the optimal conformation for catalysis will be sampled., (© 2021. The Author(s), under exclusive licence to Springer Nature B.V.)

A spacious Fe(ii)-iminopyridine self-assembled cage complex can catalyze the oxidative dimerization of alkanethiols, with air as stoichiometric oxidant. The reaction is aided by selective molecular recognition of the reactants, and the active catalyst is derived from the Fe(ii) centers that provide the structural vertices of the host. The host is even capable of size-selective oxidation and can discriminate between alkanethiols of identical reactivity, based solely on size.

Crystals composed of photoreactive molecules represent a new class of photomechanical materials with the potential to generate large forces on fast timescales. An example is the photodimerization of 9- tert -butyl-anthracene ester ( 9TBAE ) in molecular crystal nanorods that leads to an average elongation of 8%. Previous work showed that this expansion results from the formation of a metastable crystalline product. In this article, it is shown how a novel combination of ensemble oriented-crystal solid-state NMR, X-ray diffraction, and first principles computational modeling can be used to establish the absolute unit cell orientations relative to the shape change, revealing the atomic-resolution mechanism for the photomechanical response and enabling the construction of a model that predicts an elongation of 7.4%, in good agreement with the experimental value. According to this model, the nanorod expansion does not result from an overall change in the volume of the unit cell, but rather from an anisotropic rearrangement of the molecular contents. The ability to understand quantitatively how molecular-level photochemistry generates mechanical displacements allows us to predict that the expansion could be tuned from +9% to -9.5% by controlling the initial orientation of the unit cell with respect to the nanorod axis. This application of NMR-assisted crystallography provides a new tool capable of tying the atomic-level structural rearrangement of the reacting molecular species to the mechanical response of a nanostructured sample., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Structural studies with tryptophan synthase (TS) bienzyme complex (α2β2 TS) from Salmonella typhimurium have been performed to better understand its catalytic mechanism, allosteric behavior, and details of the enzymatic transformation of substrate to product in PLP-dependent enzymes. In this work, a novel expression system to produce the isolated α- and isolated β-subunit allowed the purification of high amounts of pure subunits and α2β2 StTS complex from the isolated subunits within 2 days. Purification was carried out by affinity chromatography followed by cleavage of the affinity tag, ammonium sulfate precipitation, and size exclusion chromatography (SEC). To better understand the role of key residues at the enzyme β-site, site-direct mutagenesis was performed in prior structural studies. Another protocol was created to purify the wild type and mutant α2β2 StTS complexes. A simple, fast and efficient protocol using ammonium sulfate fractionation and SEC allowed purification of α2β2 StTS complex in a single day. Both purification protocols described in this work have considerable advantages when compared with previous protocols to purify the same complex using PEG 8000 and spermine to crystalize the α2β2 StTS complex along the purification protocol. Crystallization of wild type and some mutant forms occurs under slightly different conditions, impairing the purification of some mutants using PEG 8000 and spermine. To prepare crystals suitable for x-ray crystallographic studies several efforts were made to optimize crystallization, crystal quality and cryoprotection. The methods presented here should be generally applicable for purification of tryptophan synthase subunits and wild type and mutant α2β2 StTS complexes.

The increased sensitivity under weighted non-uniform sampling (NUS) is demonstrated and quantified using Monte Carlo simulations of nuclear magnetic resonance (NMR) time- and frequency-domain signals. The concept of spectral knowledge is introduced and shown to be superior to the frequency-domain signal-to-noise ratio for assessing the quality of NMR data. Two methods for rigorously preserving spectral knowledge and the time-domain NUS knowledge enhancement upon transformation to the frequency domain are demonstrated, both theoretically and numerically. The first, non-uniform weighted sampling using consistent root-mean-square noise, is applicable to data sampled on the Nyquist grid, whereas the second, the block Fourier transform using consistent root-mean-square noise, can be used to transform time-domain data acquired with arbitrary, off-grid NUS.

Backbone assignments for the isolated α-subunit of Salmonella typhimurium tryptophan synthase (TS) are reported based on triple resonance solution-state NMR experiments on a uniformly 2 H, 13 C, 15 N-labeled sample. From the backbone chemical shifts, secondary structure and random coil index order parameters (RCI-S 2 ) are predicted. Titration with the 3-indole-D-glycerol 3'-phosphate analog, N-(4'-trifluoromethoxybenzenesulfonyl)-2-aminoethyl phosphate (F9), leads to chemical shift perturbations indicative of conformational changes from which an estimate of the dissociation constant is obtained. Comparisons of the backbone chemical-shifts, RCI-S 2 values, and site-specific relaxation times with and without F9 reveal allosteric changes including modulation in secondary structures and loop rigidity induced upon ligand binding. A comparison is made to the X-ray crystal structure of the α-subunit in the full TS αββα bi-enzyme complex and to two new X-ray crystal structures of the isolated TS α-subunit reported in this work.

A self-assembled Fe 4 L 6 cage is capable of co-encapsulating multiple carboxylic acid containing guests in its cavity, and these acids can act as cofactors for cage-catalyzed nucleophilic substitutions. The kinetics of the substitution reaction depend on the size, shape, and binding affinity of each of the components, and small structural changes in guest size can have large effects on the reaction. The host is quite promiscuous and is capable of binding multiple guests with micromolar binding affinities while retaining the ability to effect turnover and catalysis. Substrate binding modes vary widely, from simple 1:1 complexes to 1:2 complexes that can show either negative or positive cooperativity, depending on the guest. The molecularity of the dissociative substitution reaction varies, depending on the electrophile leaving group, acid cofactor, and nucleophile size: small changes in the nature of substrate can have large effects on reaction kinetics, all controlled by selective molecular recognition in the cage interior.

One-dimensional 1 H NMR experiments were conducted for aqueous solutions of glycosaminoglycan oligosaccharides to measure the amide proton temperature coefficients and activation energy barriers for solvent exchange and evaluate the effect of pH on the solvent exchange properties. A library of mono- and oligosaccharides was prepared by enzymatic depolymerization of amide-containing polysaccharides and by chemical modification of heparin and heparan sulfate saccharides including members that contain a 3- O-sulfated glucosamine residue. The systematic evaluation of this saccharide library facilitated assessment of the effects of structural characteristics, such as size, sulfation number and site, and glycosidic linkage, on amide proton solvent exchange rates. Charge repulsion by neighboring negatively charged sulfate and carboxylate groups was found to have a significant impact on the catalysis of amide proton solvent exchange by hydroxide. This observation leads to the conclusion that solvent exchange rates must be interpreted within the context of a given chemical environment. On their own, slow exchange rates do not conclusively establish the involvement of a labile proton in a hydrogen bond, and additional supporting experimental evidence such as reduced temperature coefficients is required.

The representation of nuclear magnetic resonance (NMR) tensors as surfaces on three-dimensional molecular models is an information-rich presentation that highlights the geometric relationship between tensor principal components and the underlying molecular and electronic structure. Here, we describe a new computational tool, TensorView, for depicting NMR tensors on the molecular framework. This package makes use of the graphical interface and built-in molecular display functionality present within the Mathematica programming environment and is robust for displaying tensor properties from a broad range of commercial and user-specific computational chemistry packages. Two mathematical forms for representing tensor interaction surfaces are presented, the popular ellipsoidal construct and the more technically correct "ovaloid" form. Examples are provided for chemical shielding and shift tensors, dipole-dipole and quadrupolar couplings, and atomic anisotropic displacement parameters (thermal ellipsoids) derived from NMR crystallography., (© 2018 John Wiley & Sons, Ltd.)

The hydrogen-bonding environments at the COOH moiety in eight polycrystalline polymorphs of palmitic acid are explored using solid-state NMR. Although most phases have no previously reported crystal structure, measured 13C chemical shift tensors for COOH moieties, combined with DFT modeling establish that all phases crystallize with a cyclic dimer (R22(8)) hydrogen bonding arrangement. Phases A2, Bm and Em have localized OH hydrogens while phase C has a dynamically disordered OH hydrogen. The phase designated As is a mix of five forms, including 27.4% of Bm and four novel phases not fully characterized here due to insufficient sample mass. For phases A2, Bm, Em, and C the anisotropic uncertainties in the COOH hydrogen atom positions are established using a Monte Carlo sampling scheme. Sampled points are retained or rejected at the ±1σ level based upon agreement of DFT computed 13COOH tensors with experimental values. The collection of retained hydrogen positions bear a remarkable resemblance to the anisotropic displacement parameters (i.e. thermal ellipsoids) from diffraction studies. We posit that this similarity is no mere coincidence and that the two are fundamentally related. The volumes of NMR-derived anisotropic displacement ellipsoids for phases with localized OH hydrogens are 4.1 times smaller than those derived from single crystal X-ray diffraction and 1.8 times smaller than the volume of benchmark single crystal neutron diffraction values.

The effects on the molecular recognition properties of water-soluble deep cavitand hosts upon embedding them in phosphocholine lipid bilayer environments have been studied by 2D NMR experiments. By employing suitable guests containing 19 F or 13 C nuclei that can be encapsulated inside the host, 2D EXSY NMR experiments can be used to analyze and compare the in/out guest exchange rates in aqueous solution, isotropically tumbling micelles, or magnetically ordered bicelles. These analyses show that embedding the deep cavitands in lipid bilayers slows the guest exchange rate, due to the lipids acting as a "compression sleeve" around the host, restricting guest egress. This effect also enhances guest conformations in the host that are not observed in free solution, such as axial cyclohexane conformers and ketone hydrates.

NMR studies measuring chemical shift tensors are increasingly being employed to assign structure in difficult-to-crystallize solids. For small organic molecules, such studies usually focus on 13 C sites, but proteins and peptides are more commonly described using 15 N amide sites. An important and often neglected consideration when measuring shift tensors is the evaluation of their accuracy against benchmark standards, where available. Here we measure 15 N tensors in the dipeptide glycylglycine at natural abundance using the slow-spinning FIREMAT method with SPINAL-64 decoupling. The accuracy of these 15 N tensors is evaluated by comparing to benchmark single crystal NMR 15 N measurements and found to be statistically indistinguishable. These FIREMAT experimental results are further used to evaluate the accuracy of theoretical predictions of tensors from four different density functional theory (DFT) methods that include lattice effects. The best theoretical approach provides a root mean square (rms) difference of ±3.9 ppm and is obtained from a fragment-based method and the PBE0 density functional., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Carbanionic intermediates play a central role in the catalytic transformations of amino acids performed by pyridoxal-5'-phosphate (PLP)-dependent enzymes. Here, we make use of NMR crystallography-the synergistic combination of solid-state nuclear magnetic resonance, X-ray crystallography, and computational chemistry-to interrogate a carbanionic/quinonoid intermediate analogue in the β-subunit active site of the PLP-requiring enzyme tryptophan synthase. The solid-state NMR chemical shifts of the PLP pyridine ring nitrogen and additional sites, coupled with first-principles computational models, allow a detailed model of protonation states for ionizable groups on the cofactor, substrates, and nearby catalytic residues to be established. Most significantly, we find that a deprotonated pyridine nitrogen on PLP precludes formation of a true quinonoid species and that there is an equilibrium between the phenolic and protonated Schiff base tautomeric forms of this intermediate. Natural bond orbital analysis indicates that the latter builds up negative charge at the substrate C α and positive charge at C4' of the cofactor, consistent with its role as the catalytic tautomer. These findings support the hypothesis that the specificity for β-elimination/replacement versus transamination is dictated in part by the protonation states of ionizable groups on PLP and the reacting substrates and underscore the essential role that NMR crystallography can play in characterizing both chemical structure and dynamics within functioning enzyme active sites.

The performance of fragment-based ab initio(1)H, (13)C, (15)N and (17)O chemical shift predictions is assessed against experimental NMR chemical shift data in four benchmark sets of molecular crystals. Employing a variety of commonly used density functionals (PBE0, B3LYP, TPSSh, OPBE, PBE, TPSS), we explore the relative performance of cluster, two-body fragment, and combined cluster/fragment models. The hybrid density functionals (PBE0, B3LYP and TPSSh) generally out-perform their generalized gradient approximation (GGA)-based counterparts. (1)H, (13)C, (15)N, and (17)O isotropic chemical shifts can be predicted with root-mean-square errors of 0.3, 1.5, 4.2, and 9.8 ppm, respectively, using a computationally inexpensive electrostatically embedded two-body PBE0 fragment model. Oxygen chemical shieldings prove particularly sensitive to local many-body effects, and using a combined cluster/fragment model instead of the simple two-body fragment model decreases the root-mean-square errors to 7.6 ppm. These fragment-based model errors compare favorably with GIPAW PBE ones of 0.4, 2.2, 5.4, and 7.2 ppm for the same (1)H, (13)C, (15)N, and (17)O test sets. Using these benchmark calculations, a set of recommended linear regression parameters for mapping between calculated chemical shieldings and observed chemical shifts are provided and their robustness assessed using statistical cross-validation. We demonstrate the utility of these approaches and the reported scaling parameters on applications to 9-tert-butyl anthracene, several histidine co-crystals, benzoic acid and the C-nitrosoarene SnCl2(CH3)2(NODMA)2.

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of β-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations., Competing Interests: The authors declare no conflicts of interest.

Four new X-ray structures of tryptophan synthase (TS) crystallized with varying numbers of the amphipathic N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) molecule are presented. These structures show one of the F6 ligands threaded into the tunnel from the β-site and reveal a distinct hydrophobic region. Over this expanse, the interactions between F6 and the tunnel are primarily nonpolar, while the F6 phosphoryl group fits into a polar pocket of the β-subunit active site. Further examination of TS structures reveals that one portion of the tunnel (T1) binds clusters of water molecules, whereas waters are not observed in the nonpolar F6 binding region of the tunnel (T2). MD simulation of another TS structure with an unobstructed tunnel also indicates the T2 region of the tunnel excludes water, consistent with a dewetted state that presents a significant barrier to the transfer of water into the closed β-site. We conclude that hydrophobic molecules can freely diffuse between the α- and β-sites via the tunnel, while water does not. We propose that exclusion of water serves to inhibit reaction of water with the α-aminoacrylate intermediate to form ammonium ion and pyruvate, a deleterious side reaction in the αβ-catalytic cycle. Finally, while most TS structures show βPhe280 partially blocking the tunnel between the α- and β-sites, new structures show an open tunnel, suggesting the flexibility of the βPhe280 side chain. Flexible docking studies and MD simulations confirm that the dynamic behavior of βPhe280 allows unhindered transfer of indole through the tunnel, therefore excluding a gating role for this residue., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Oxygen is an essential participant in the acid-base chemistry that takes place within many enzyme active sites, yet has remained virtually silent as a probe in NMR spectroscopy. Here, we demonstrate the first use of solution-state (17)O quadrupole central-transition NMR spectroscopy to characterize enzymatic intermediates under conditions of active catalysis. In the 143 kDa pyridoxal-5'-phosphate-dependent enzyme tryptophan synthase, reactions of the α-aminoacrylate intermediate with the nucleophiles indoline and 2-aminophenol correlate with an upfield shift of the substrate carboxylate oxygen resonances. First principles calculations suggest that the increased shieldings for these quinonoid intermediates result from the net increase in the charge density of the substrate-cofactor π-bonding network, particularly at the adjacent α-carbon site., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The importance of protonation states and proton transfer in pyridoxal 5'-phosphate (PLP)-chemistry can hardly be overstated. Although experimental approaches to investigate pKa values can provide general guidance for assigning proton locations, only static pictures of the chemical species are available. To obtain the overall protein dynamics for the interpretation of detailed enzyme catalysis in this study, guided by information from solid-state NMR, we performed molecular dynamics (MD) simulations for the PLP-dependent enzyme tryptophan synthase (TRPS), whose catalytic mechanism features multiple quasi-stable intermediates. The primary objective of this work is to elucidate how the position of a single proton on the reacting substrate affects local and global protein dynamics during the catalytic cycle. In general, proteins create a chemical environment and an ensemble of conformational motions to recognize different substrates with different protonations. The study of these interactions in TRPS shows that functional groups on the reacting substrate, such as the phosphoryl group, pyridine nitrogen, phenolic oxygen and carboxyl group, of each PLP-bound intermediate play a crucial role in constructing an appropriate molecular interface with TRPS. In particular, the protonation states of the ionizable groups on the PLP cofactor may enhance or weaken the attractions between the enzyme and substrate. In addition, remodulation of the charge distribution for the intermediates may help generate a suitable environment for chemical reactions. The results of our study enhance knowledge of protonation states for several PLP intermediates and help to elucidate their effects on protein dynamics in the function of TRPS and other PLP-dependent enzymes., (© 2015 The Protein Society.)

The proposed mechanism for tryptophan synthase shows βLys87 playing multiple catalytic roles: it bonds to the PLP cofactor, activates C4' for nucleophilic attack via a protonated Schiff base nitrogen, and abstracts and returns protons to PLP-bound substrates (i.e. acid-base catalysis). ε-¹⁵N-lysine TS was prepared to access the protonation state of βLys87 using ¹⁵N solid-state nuclear magnetic resonance (SSNMR) spectroscopy for three quasi-stable intermediates along the reaction pathway. These experiments establish that the protonation state of the ε-amino group switches between protonated and neutral states as the β-site undergoes conversion from one intermediate to the next during catalysis, corresponding to mechanistic steps where this lysine residue has been anticipated to play alternating acid and base catalytic roles that help steer reaction specificity in tryptophan synthase catalysis. This article is part of a Special Issue entitled: Cofactor-dependent proteins: evolution, chemical diversity and bio-applications. Guest Editors: Andrea Mozzarelli and Loredano Pollegioni., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Ab initio chemical shielding calculations greatly facilitate the interpretation of nuclear magnetic resonance (NMR) chemical shifts in biological systems, but the large sizes of these systems requires approximations in the chemical models used to represent them. Achieving good convergence in the predicted chemical shieldings is necessary before one can unravel how other complex structural and dynamical factors affect the NMR measurements. Here, we investigate how to balance trade-offs between using a better basis set or a larger cluster model for predicting the chemical shieldings of the substrates in two representative examples of protein-substrate systems involving different domains in tryptophan synthase: the N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F9) ligand which binds in the α active site, and the 2-aminophenol quinonoid intermediate formed in the β active site. We first demonstrate that a chemically intuitive three-layer, locally dense basis model that uses a large basis on the substrate, a medium triple-zeta basis to describe its hydrogen-bonding partners and/or surrounding van der Waals cavity, and a crude basis set for more distant atoms provides chemical shieldings in good agreement with much more expensive large basis calculations. Second, long-range quantum mechanical interactions are important, and one can accurately estimate them as a small-basis correction to larger-basis calculations on a smaller cluster. The combination of these approaches enables one to perform density functional theory NMR chemical shift calculations in protein systems that are well-converged with respect to both basis set and cluster size.

The acid-base chemistry that drives catalysis in pyridoxal-5'-phosphate (PLP)-dependent enzymes has been the subject of intense interest and investigation since the initial identification of PLP's role as a coenzyme in this extensive class of enzymes. It was first proposed over 50 years ago that the initial step in the catalytic cycle is facilitated by a protonated Schiff base form of the holoenzyme in which the linking lysine ε-imine nitrogen, which covalently binds the coenzyme, is protonated. Here we provide the first (15)N NMR chemical shift measurements of such a Schiff base linkage in the resting holoenzyme form, the internal aldimine state of tryptophan synthase. Double-resonance experiments confirm the assignment of the Schiff base nitrogen, and additional (13)C, (15)N, and (31)P chemical shift measurements of sites on the PLP coenzyme allow a detailed model of coenzyme protonation states to be established.

Metal-organic frameworks (MOFs) with cationic frameworks and mobile anions have many applications from sensing, anion exchange and separation, to fast ion conductivity. Despite recent progress, the vast majority of MOFs have neutral frameworks. A common mechanism for the formation of neutral frameworks is the attachment of anionic species such as F(-) or OH(-) to the framework metal sites, neutralizing an otherwise cationic scaffolding. Here, we report a general method capable of converting such neutral frameworks directly into cationic ones with concurrent generation of mobile anions. Our method is based on the differential affinity between distinct metal ions with framework anionic species. Specifically, Al(3+) is used to strip F(-) anions away from framework Cr(3+) sites, leading to cationic frameworks with mobile Cl(-) anions. The subsequent anion exchange with OH(-) further leads to a porous network with mobile OH(-) anions. New materials prepared by anion stripping can undergo ion exchange with anionic organic dyes and also exhibit much improved ionic conductivity compared to the original unmodified MOFs.

The solid-state photodimerization of 9-methylanthracene is used as a model system to investigate how crystal morphology and reaction dynamics affect photomechanical deformations of single microcrystals. By varying the crystallization conditions, two different crystal shapes, microneedles and microribbons, are grown on a clean water surface. The microribbons twist under irradiation, while the microneedles bend. In both shapes, the maximum deformation occurs at roughly the midpoint of the reaction, while further dimerization causes the crystals return to their original shapes. Powder X-ray diffraction patterns establish that the needles and ribbons have the same crystal orientation and that the photoreaction proceeds in a crystal-to-crystal manner. NMR spin-lattice relaxation measurements are consistent with the rapid formation of large (>100 nm) dimer crystal domains. Simultaneous measurement of the needle bending and monomer fluorescence signal allows us to correlate the bending with the reaction progress. The behavior is qualitatively reproduced by a model in which the motion is driven by strain between spatially distinct reactant and product domains, also called heterometry. We consider several different mechanisms that could give rise to these spatially distinct domains. The ability to control the photoinduced crystal deformation by manipulating crystal shape and solid-state reaction kinetics suggests that photoreactive molecular crystals may be useful for generating well-defined motions on small length scales.

Heparin is best known for its anticoagulant activity, which is mediated by the binding of a specific pentasaccharide sequence to the protease inhibitor antithrombin-III (AT-III). Although heparin oligosaccharides are thought to be flexible in aqueous solution, the recent discovery of a hydrogen bond between the sulfamate (NHSO3(-)) proton and the adjacent 3-O-sulfo group of the 3,6-O-sulfated N-sulfoglucosamine residue of the Arixtra (fondaparinux sodium) pentasaccharide demonstrates that definable elements of local structure are accessed. Molecular dynamics simulations of Arixtra suggest the presence of additional hydrogen bonds involving the C3-OH groups of the glucuronic acid and 2-O-sulfo-iduronic acid residues. NMR measurements of temperature coefficients, chemical shift differences, and solvent exchange rate constants provide experimental confirmation of these hydrogen bonds. We note that the extraction of rate constants from cross-peak buildup curves in 2D exchange spectroscopy is complicated by the presence of radiation damping in aqueous solution. A straightforward model is presented that explicitly takes into account the effects of radiation damping on the water proton relaxation and is sufficiently robust to provide an accurate measure of the proton exchange rate between the analyte hydroxyl protons and water.

Two-dimensional (15)N chemical shift/(1)H chemical shift and three-dimensional (1)H-(15)N dipolar coupling/(15)N chemical shift/(1)H chemical shift MAS solid-state NMR correlation spectra of the filamentous bacteriophage Pf1 major coat protein show single-site resolution in noncrystalline, intact-phage preparations. The high sensitivity and resolution result from (1)H detection at 600MHz under 50kHz magic angle spinning using ∼0.5mg of perdeuterated and uniformly (15)N-labeled protein in which the exchangeable amide sites are partially or completely back-exchanged (reprotonated). Notably, the heteronuclear (1)H-(15)N dipolar coupling frequency dimension is shown to select among (15)N resonances, which will be useful in structural studies of larger proteins where the resonances exhibit a high degree of overlap in multidimensional chemical shift correlation spectra., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Vibrational spectroscopy and NMR demonstrate that the proton-bound dimer of 1-methylcytosine, 1, has an unsymmetrical structure at room temperature. In the gas phase, investigation of isolated homodimer 1 reveals five fundamental NH vibrations by IR Multiple Photon Dissociation (IRMPD) action spectroscopy. The NH···N stretching vibration between the two ring nitrogens exhibits a frequency of 1570 cm(-1), as confirmed by examination of the proton-bound homodimers of 5-fluoro-1-methycytosine, 2, and of 1,5-dimethylcytosine, 3, which display absorptions in the same region that disappear upon deuterium substitution. (13)C, and (15)N NMR of the solid iodide salt of 1 confirm the nonequivalence of the two rings in the anhydrous proton-bound homodimer at room temperature. IRMPD spectra of the three possible heterodimers also show NH···N stretches in the same domain, and at least one of the heterodimers, the proton-bound dimer of 1,5-dimethylcytosine with 1-methylcytosine, exhibits two bands suggestive of the presence of two tautomers close in energy.

Sensitivity and resolution are the two fundamental obstacles to extending solid-state nuclear magnetic resonance to even larger protein systems. Here, a novel long-observation-window band-selective homonuclear decoupling (LOW BASHD) scheme is introduced that increases resolution up to a factor of 3 and sensitivity up to 1.8 by decoupling backbone alpha-carbon (C(α)) and carbonyl (C') nuclei in U-(13)C-labeled proteins during direct (13)C acquisition. This approach introduces short (<200 μs) pulse breaks into much longer (~8 ms) sampling windows to efficiently refocus the J-coupling interaction during detection while avoiding the deleterious effects on sensitivity inherent in rapid stroboscopic band-selective homonuclear decoupling techniques. A significant advantage of LOW-BASHD detection is that it can be directly incorporated into existing correlation methods, as illustrated here for 2D CACO, NCO, and NCA correlation spectroscopy applied to the β1 immunoglobulin binding domain of protein G and 3D CBCACO correlation spectroscopy applied to the α-subunit of tryptophan synthase., (Copyright © 2013 Elsevier Inc. All rights reserved.)

NMR crystallography--the synergistic combination of X-ray diffraction, solid-state NMR spectroscopy, and computational chemistry--offers unprecedented insight into three-dimensional, chemically detailed structure. Initially, researchers used NMR crystallography to refine diffraction data from organic and inorganic solids. Now we are applying this technique to explore active sites in biomolecules, where it reveals chemically rich detail concerning the interactions between enzyme site residues and the reacting substrate. Researchers cannot achieve this level of detail from X-ray, NMR,or computational methodologies in isolation. For example, typical X-ray crystal structures (1.5-2.5 Å resolution) of enzyme-bound intermediates identify possible hydrogen-bonding interactions between site residues and substrate but do not directly identify the protonation states. Solid-state NMR can provide chemical shifts for selected atoms of enzyme-substrate complexes, but without a larger structural framework in which to interpret them only empirical correlations with local chemical structure are possible. Ab initio calculations and molecular mechanics can build models for enzymatic processes, but they rely on researcher-specified chemical details. Together, however, X-ray diffraction, solid-state NMR spectroscopy, and computational chemistry can provide consistent and testable models for structure and function of enzyme active sites: X-ray crystallography provides a coarse framework upon which scientists can develop models of the active site using computational chemistry; they can then distinguish these models by comparing calculated NMR chemical shifts with the results of solid-state NMR spectroscopy experiments. Conceptually, each technique is a puzzle piece offering a generous view of the big picture. Only when correctly pieced together, however, can they reveal the big picture at the highest possible resolution. In this Account, we detail our first steps in the development of NMR crystallography applied to enzyme catalysis. We begin with a brief introduction to NMR crystallography and then define the process that we have employed to probe the active site in the β-subunit of tryptophan synthase with unprecedented atomic-level resolution. This approach has resulted in a novel structural hypothesis for the protonation state of the quinonoid intermediate in tryptophan synthase and its surprising role in directing the next step in the catalysis of L-Trp formation.

The allosteric regulation of substrate channeling in tryptophan synthase involves ligand-mediated allosteric signaling that switches the α- and β-subunits between open (low activity) and closed (high activity) conformations. This switching prevents the escape of the common intermediate, indole, and synchronizes the α- and β-catalytic cycles. (19)F NMR studies of bound α-site substrate analogues, N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) and N-(4'-trifluoromethoxybenzenesulfonyl)-2-aminoethyl phosphate (F9), were found to be sensitive NMR probes of β-subunit conformation. Both the internal and external aldimine F6 complexes gave a single bound peak at the same chemical shift, while α-aminoacrylate and quinonoid F6 complexes all gave a different bound peak shifted by +1.07 ppm. The F9 complexes exhibited similar behavior, but with a corresponding shift of -0.12 ppm. X-ray crystal structures show the F6 and F9 CF3 groups located at the α-β subunit interface and report changes in both the ligand conformation and the surrounding protein microenvironment. Ab initio computational modeling suggests that the change in (19)F chemical shift results primarily from changes in the α-site ligand conformation. Structures of α-aminoacrylate F6 and F9 complexes and quinonoid F6 and F9 complexes show the α- and β-subunits have closed conformations wherein access of ligands into the α- and β-sites from solution is blocked. Internal and external aldimine structures show the α- and β-subunits with closed and open global conformations, respectively. These results establish that β-subunits exist in two global conformational states, designated open, where the β-sites are freely accessible to substrates, and closed, where the β-site portal into solution is blocked. Switching between these conformations is critically important for the αβ-catalytic cycle.

Sulfamate groups (NHSO(3)(-)) are important structural elements in the glycosaminoglycans (GAGs) heparin and heparan sulfate (HS). In this work, proton nuclear magnetic resonance (NMR) line-shape analysis is used to explore the solvent exchange properties of the sulfamate NH groups within heparin-related mono-, di-, tetra- and pentasaccharides as a function of pH and temperature. The results of these experiments identified a persistent hydrogen bond within the Arixtra (fondaparinux sodium) pentasaccharide between the internal glucosamine sulfamate NH and the adjacent 3-O-sulfo group. This discovery provides new insights into the solution structure of the Arixtra pentasaccharide and suggests that 3-O-sulfation of the heparin N-sulfoglucosamine (GlcNS) residues pre-organize the secondary structure in a way that facilitates binding to antithrombin-III. NMR studies of the GlcNS NH groups can provide important information about heparin structure complementary to that available from NMR spectral analysis of the carbon-bound protons.