Polymorph Identification for Flexible Molecules: Linear Regression Analysis of Experimental and Calculated Solution- and Solid-State NMR Data.

The Δδ regression approach of Blade et al. [ J. Phys. Chem. A 2020, 124(43), 8959-8977] for accurately discriminating between solid forms using a combination of experimental solution- and solid-state NMR data with density functional theory (DFT) calculation is here extended to molecules with multiple conformational degrees of freedom, using furosemide polymorphs as an exemplar. As before, the differences in measured 1H and 13C chemical shifts between solution-state NMR and solid-state magic-angle spinning (MAS) NMR (Δδexperimental) are compared to those determined by gauge-including projector augmented wave (GIPAW) calculations (Δδcalculated) by regression analysis and a t-test, allowing the correct furosemide polymorph to be precisely identified. Monte Carlo random sampling is used to calculate solution-state NMR chemical shifts, reducing computation times by avoiding the need to systematically sample the multidimensional conformational landscape that furosemide occupies in solution. The solvent conditions should be chosen to match the molecule's charge state between the solution and solid states. The Δδ regression approach indicates whether or not correlations between Δδexperimental and Δδcalculated are statistically significant; the approach is differently sensitive to the popular root mean squared error (RMSE) method, being shown to exhibit a much greater dynamic range. An alternative method for estimating solution-state NMR chemical shifts by approximating the measured solution-state dynamic 3D behavior with an ensemble of 54 furosemide crystal structures (polymorphs and cocrystals) from the Cambridge Structural Database (CSD) was also successful in this case, suggesting new avenues for this method that may overcome its current dependency on the prior determination of solution dynamic 3D structures.

Crystal structures of tolfenamic acid polymorphic forms I and II with precise hydrogen-atom positions for nuclear magnetic resonance studies.

The structures of tolfenamic acid [TFA; 2-(3-chloro-2-methyl-anilino)benzoic acid, C14H12ClNO2] polymorph forms I and II have been redetermined [compare Andersen et al. (1989 ▸). J. Chem. Soc., Perkin Trans. 2, pp. 1443-1447] with improved precision using high-resolution X-ray diffraction data and Hirshfield atom refinement in order to better define both hydrogen-atom locations and their associated bond lengths. Covalent bond lengths to hydrogen were found to be significantly longer throughout both structures, especially for the anilino H atom, which is involved in an important intra-molecular N-H⋯O hydrogen bond to the carb-oxy-lic acid group. This hydrogen bond is shown to clearly perturb the electron density around both oxygen atoms in the latter group. The extended structures of both polymorphs feature carb-oxy-lic acid inversion dimers. These structures provide an improved foundation for nuclear magnetic resonance studies in both solution and the solid state.

Conformations in Solution and in Solid-State Polymorphs: Correlating Experimental and Calculated Nuclear Magnetic Resonance Chemical Shifts for Tolfenamic Acid.

A new approach for quantitively assessing putative crystal structures with applications in crystal structure prediction (CSP) is introduced that is based upon experimental solution- and magic-angle spinning (MAS) solid-state NMR data and density functional theory (DFT) calculation. For the specific case of tolfenamic acid (TFA), we consider experimental solution-state NMR for a range of solvents, experimental MAS NMR of polymorphs I and II, and DFT calculations for four polymorphs. The change in NMR chemical shift observed in passing from the solution state to the solid state (ΔδExperimental) is calculated as the difference between 1H and 13C experimental solid-state chemical shifts for each polymorphic form (δSolid expt) and the corresponding solution-state NMR chemical shifts (δSolution expt). Separately, we use the gauge-included projector augmented wave (GIPAW) method to calculate the NMR chemical shifts for each form (δSolid calc) and for TFA in solution (δSolution calc) using the dynamic 3D solution conformational ensemble determined from NMR spectroscopy. The calculated change in passing from the solution state to the solid state (ΔδCalculated) is then calculated as the difference of δSolid calc and δSolution calc. Regression analysis for ΔδCalculated against ΔδExperimental followed by a t-test for statistical significance provides a robust quantitative assessment. We show that this assessment clearly identifies the correct polymorph, i.e., when comparing ΔδExperimental based on the experimental MAS NMR chemical shifts of form I or II with ΔδCalculated based on calculated chemical shifts for polymorphs I, II, III, and IV. Complementarity to the established approach of comparing δSolid expt to δSolid calc is explored. We further show that our approach is applicable if there are no solid-state crystal structure data. Specifically, δSolid calc in ΔδCalculated is replaced by the chemical shift for an isolated molecule with a specific conformation. Sampling conformations at specific 15° angle values and comparing them against experimental 13C chemical shift data for forms I and II identifies matching narrow ranges of conformations, successfully predicting the conformation of tolfenamic acid in each form. This methodology can therefore be used in crystal structure prediction to both reduce the initial conformational search space and also quantitatively assess subsequent putative structures to reliably and unambiguously identify the correct structure.

Pushing the Limits of Molecular Crystal Structure Determination From Powder Diffraction Data in High-Throughput Chemical Environments.

Crystal structure determination from powder diffraction data (SDPD) using the DASH software package is evaluated for data recorded using transmission capillary, transmission flat plate, and reflection flat plate geometries on a selection of pharmaceutical compounds. We show that transmission capillary geometry remains the best option when crystal structure determination is the primary consideration and, as expected, reflection flat plate geometry is not recommended for SDPD because of preferred orientation effects. However, the quality of crystal structures obtained from transmission plate instruments can be excellent, and the convenience factor for sample preparation, throughput, and retrieval is higher than that of transmission capillary instruments. Indeed, it is possible to solve crystal structures within an hour of a polycrystalline sample arriving in the laboratory, which has clear implications for making small-molecule crystal structures more routinely available to the practicing laboratory medicinal chemist. With appropriate modifications to crystal structure determination software, it can be imagined that SDPD could become a rapid turn-around walk-up analytical service in high-throughput chemical environments.

Accurate measurement of long range proton-carbon scalar coupling constants.

The accuracy and practicality of measuring heteronuclear scalar coupling constants, nJCH, from modern NMR experimental methods is examined, based on F1 or F2 evolution of nJCH in HSQMBC (including EXSIDE) and HMBC experiments. The results from these methods are compared to both robust experimental data (derived from coupled 13C spectra), computed (Density Functional Theory) and literature values where available. We report on the accuracy, ease of use and time efficiency of these multi-dimensional methods and highlight their extent and limitations.

GPR103 antagonists demonstrating anorexigenic activity in vivo: design and development of pyrrolo[2,3-c]pyridines that mimic the C-terminal Arg-Phe motif of QRFP26.

GPR103, a G-protein coupled receptor, has been reported to have orexigenic properties through activation by the endogenous neuropeptide ligands QRFP26 and QRFP43. Recognizing that central administration of QRFP26 and QRFP43 increases high fat food intake in rats, we decided to investigate if antagonists of GPR103 could play a role in managing feeding behaviors. Here we present the development of a new series of pyrrolo[2,3-c]pyridines as GPR103 small molecule antagonists with GPR103 affinity, drug metabolism and pharmacokinetics and safety parameters suitable for drug development. In a preclinical obesity model measuring food intake, the anorexigenic effect of a pyrrolo[2,3-c]pyridine GPR103 antagonist was demonstrated. In addition, the dynamic 3D solution structure of the C-terminal heptapeptide of the endogenous agonist QRFP26(20-26) was determined using NMR. The synthetic pyrrolo[2,3-c]pyridine antagonists were compared to this experimental structure, which displayed a possible overlay of pharmacophore features supportive for further design of GPR103 antagonists.

A refined model for the TSG-6 link module in complex with hyaluronan: use of defined oligosaccharides to probe structure and function.

Tumor necrosis factor-stimulated gene-6 (TSG-6) is an inflammation-associated hyaluronan (HA)-binding protein that contributes to remodeling of HA-rich extracellular matrices during inflammatory processes and ovulation. The HA-binding domain of TSG-6 consists solely of a Link module, making it a prototypical member of the superfamily of proteins that interacts with this high molecular weight polysaccharide composed of repeating disaccharides of D-glucuronic acid and N-acetyl-D-glucosamine (GlcNAc). Previously we modeled a complex of the TSG-6 Link module in association with an HA octasaccharide based on the structure of the domain in its HA-bound conformation. Here we have generated a refined model for a HA/Link module complex using novel restraints identified from NMR spectroscopy of the protein in the presence of 10 distinct HA oligosaccharides (from 4- to 8-mers); the model was then tested using unique sugar reagents, i.e. chondroitin/HA hybrid oligomers and an octasaccharide in which a single sugar ring was (13)C-labeled. The HA chain was found to make more extensive contacts with the TSG-6 surface than thought previously, such that a D-glucuronic acid ring makes stacking and ionic interactions with a histidine and lysine, respectively. Importantly, this causes the HA to bend around two faces of the Link module (resembling the way that HA binds to CD44), potentially providing a mechanism for how TSG-6 can reorganize HA during inflammation. However, the HA-binding site defined here may not play a role in TSG-6-mediated transfer of heavy chains from inter-α-inhibitor onto HA, a process known to be essential for ovulation.

Quantification of free ligand conformational preferences by NMR and their relationship to the bioactive conformation.

Accurate unbound solution 3D-structures of ligands provide unique opportunities for medicinal chemistry and, in particular, a context to understand binding thermodynamics and kinetics. Previous methods of deriving these 3D-structures have had neither the accuracy nor resolution needed for drug design and have not yet realized their potential. Here, we describe and apply a NMR methodology to the aminoglycoside streptomycin that can accurately quantify accessible 3D-space and rank the occupancy of observed conformers to a resolution that enables medicinal chemistry understanding and design. Importantly, it is based upon conventional small molecule NMR techniques and can be performed in physiologically-relevant solvents. The methodology uses multiple datasets, an order of magnitude more experimental data than previous NMR approaches and a dynamic model during refinement, is independent of computational chemistry and avoids the problem of virtual conformations. The refined set of solution 3D-shapes for streptomycin can be grouped into two major families, of which the most populated is almost identical to the 30S ribosomal subunit bioactive shape. We therefore propose that accurate unbound ligand solution conformations may, in some cases, provide a subsidiary route to bioactive shape without crystallography. This experimental technique opens up new opportunities for drug design and more so when complemented with protein co-crystal structures, SAR data and pharmacophore modeling.

Comparative pharmacology and computational modelling yield insights into allosteric modulation of human alpha7 nicotinic acetylcholine receptors.

The human alpha7 nicotinic acetylcholine receptor (nAChR) subunit and its Caenorhabditis elegans homolog, ACR-16, can generate functional recombinant homomeric receptors when expressed in Xenopus laevis oocytes. Both nAChRs express robustly in the presence of the co-injected chaperone, RIC-3, and show striking differences in the actions of a type I positive allosteric modulator (PAM), ivermectin (IVM). Type I PAMs are characterised by an increase in amplitude only of the response to acetylcholine (ACh), whereas type II PAMs exhibit, in addition, changes in time-course/desensitization of the ACh response. The type I PAMs, ivermectin, 5-hydroxyindole (5-HI), NS-1738 and genistein and the type II PAM, PNU-120596, are all active on human alpha7 but are without PAM activity on ACR-16, where they attenuate the amplitude of the ACh response. We used the published structure of avermectin B1a to generate a model of IVM, which was then docked into the candidate transmembrane allosteric binding site on alpha7 and ACR-16 in an attempt to gain insights into the observed differences in IVM actions. The new pharmacological findings and computational approaches being developed may inform the design of novel PAM drugs targeting major neurological disorders.

Investigating the molecular basis for the virulence of Escherichia coli K5 by nuclear magnetic resonance analysis of the capsule polysaccharide.

The capsular polysaccharide of Escherichia coli K5 has been hypothesised to promote virulence through its molecular mimicry of host heparan sulphate. To test this hypothesis, we have produced pure oligosaccharides from K5 capsular polysaccharide and investigated their conformational properties with ultra-high-field nuclear magnetic resonance (NMR) (900 MHz). Ultra-high-field affords a significant resolution enhancement over previous studies and allowed a full-atomic assignment of the K5 hexasaccharide for the first time. All carbohydrate rings adopt a (4)C(1) conformation, the amide sidechains have a trans orientation and the hydroxymethyl group is freely exposed to bulk solvent. Initial models of the glycosidic linkage conformation based upon simple interpretation of NOE cross-peaks suggests that the beta1-->4 linkage adopts a 3D geometry of phi approximately 60 degrees , psi approximately 0 degrees and the alpha1-->4 linkage prefers phi approximately -30 degrees , psi approximately -30 degrees (phi and psi being defined by dihedral angles involving linkage protons). In this conformation the overall molecular geometries of K5 polysaccharide, heparan sulphate and even fully-sulphated heparin are remarkably similar. These results substantiate the hypothesis that the K5 capsular polysaccharide confers virulence to E. coli K5 by being a 3D molecular mimetic of host heparan sulphate, helping it to evade detection by the mammalian immune system.

Plasticity of the TSG-6 HA-binding loop and mobility in the TSG-6-HA complex revealed by NMR and X-ray crystallography.

Tumour necrosis factor-stimulated gene-6 (TSG-6) is a glycosaminoglycan-binding protein expressed during inflammatory and inflammation-like processes. Previously NMR structures were calculated for the Link module of TSG-6 (Link_TSG6) in its free state and when bound to an octasaccharide of hyaluronan (HA(8)). Heparin was found to compete for HA binding even though it interacts at a site that is distinct from the HA-binding surface. Here we present crystallography data on the free protein, and (15)N NMR relaxation data for the uncomplexed and HA(8)-bound forms of Link_TSG6. Although the Link module is comparatively rigid overall, the free protein shows a high degree of mobility in the beta4/beta5 loop and at the Cys47-Cys68 disulfide bond, both of which are regions involved in HA binding. When bound to HA(8), this dynamic behaviour is dampened, but not eliminated, suggesting a degree of dynamic matching between the protein and sugar that may decrease the entropic penalty of complex formation. A further highly dynamic residue is Lys54, which is distant from the HA-binding site, but was previously shown to be involved in heparin binding. When HA is bound, Lys54 becomes less mobile, providing evidence for an allosteric effect linking the HA and heparin-binding sites. A mechanism is suggested involving the beta2-strand and alpha2-helix. The crystal structure of free Link_TSG6 contains five molecules in the asymmetric unit that are highly similar to the NMR structure and support the dynamic behaviour seen near the HA-binding site: they show little or no electron density for the beta4/beta5 loop and display multiple conformations for the Cys47-Cys68 disulfide bond. The crystal structures were used in docking calculations with heparin. An extended interface between a Link_TSG6 dimer and heparin 11-mer was identified that is in excellent agreement with previous mutagenesis and calorimetric data, providing the basis for further investigation of this interaction.

Temperature dependencies of amide 1H- and 15N-chemical shifts in hyaluronan oligosaccharides.

Temperature coefficients (Deltadelta/DeltaT) of amide chemical shifts of N-acetylglucosamine residues have been measured in a range of oligosaccharides of the important vertebrate polysaccharide hyaluronan. Odd- and even-numbered oligosaccharides with glucuronic acid, Delta-4,5-unsaturated glucuronic acid and N-acetylglucosamine at the termini were investigated. All amide proton temperature coefficients were only slightly less negative (-6.9 to - 9.1 ppb/ degrees C) than those of amide protons in free exchange with water (approximately equal to -11 ppb/ degrees C), indicating an absence of persistent intramolecular hydrogen bonds. With the exception of amide groups in reducing-terminal N-acetylglucosamine rings, all amide proton environments have the same temperature coefficient (-6.9 ppb/ degrees C), irrespective of differences in amide group chemical shifts and (3)J(HH) coupling constants, i.e. they do not sense subtle differences in orientation of the amide group. Amide nitrogen temperature coefficients report the same phenomena but with greater sensitivity. These data provide a set of reference values for temperature coefficients measured in other carbohydrates with acetamido sugars.

Determining the molecular basis for the pH-dependent interaction between the link module of human TSG-6 and hyaluronan.

TSG-6 is an inflammation-associated hyaluronan (HA)-binding protein that has anti-inflammatory and protective functions in arthritis and asthma as well as a critical role in mammalian ovulation. The interaction between TSG-6 and HA is pH-dependent, with a marked reduction in affinity on increasing the pH from 6.0 to 8.0. Here we have investigated the mechanism underlying this pH dependence using a combined approach of site-directed mutagenesis, NMR, isothermal titration calorimetry and microtiter plate assays. Analysis of single-site mutants of the TSG-6 Link module indicated that the loss in affinity above pH 6.0 is mediated by the change in ionization state of a histidine residue (His(4)) that is not within the HA-binding site. To understand this in molecular terms, the pH-dependent folding profile and the pK(a) values of charged residues within the Link module were determined using NMR. These data indicated that His(4) makes a salt bridge to one side-chain oxygen atom of a buried aspartate residue (Asp(89)), whereas the other oxygen is simultaneously hydrogen-bonded to a key HA-binding residue (Tyr(12)). This molecular network transmits the change in ionization state of His(4) to the HA-binding site, which explains the loss of affinity at high pH. In contrast, simulations of the pH affinity curves indicate that another histidine residue, His(45), is largely responsible for the gain in affinity for HA between pH 3.5 and 6.0. The pH-dependent interaction of TSG-6 with HA (and other ligands) provides a means of differentially regulating the functional activity of this protein in different tissue microenvironments.

Using Molecular Dynamics Simulations To Provide New Insights into Protein Structure on the Nanosecond Timescale: Comparison with Experimental Data and Biological Inferences for the Hyaluronan-Binding Link Module of TSG-6.

Link module domains play an essential role in extracellular matrix assembly and remodeling by binding to the flexible glycosaminoglycan hyaluronan. A high-resolution NMR-structure of the Link module from the protein product of tumor necrosis factor-stimulated gene-6 (Link_TSG6) has been determined, but a fuller appreciation of protein dynamics may be necessary to understand its hyaluronan-binding. Therefore, we have performed a 0.25 µs MD simulation, starting from the lowest-energy NMR-derived solution structure of Link_TSG6, with explicit water and ions, using the CHARMM22 protein force field. The simulation was as good a fit to the NMR data as the ensemble from simulated annealing, except in the ß5-ß6 loop. Furthermore, analysis revealed that secondary structure elements extended further than previously reported and underwent fast picosecond time scale dynamics, whereas nanosecond dynamics was found in certain loops. In particular, surface side chains proposed to interact with glycosaminoglycans were predicted to be highly mobile and be directed away from the protein surface. Furthermore, the hyaluronan-binding ß4-ß5 loop remained in a closed conformation, favoring an allosteric interaction mechanism. This enhanced view of the Link module provides general insight into protein dynamics and may be helpful for understanding the dynamic molecular basis of tissue assembly, remodeling, and disease processes.

Complete assignment of hyaluronan oligosaccharides up to hexasaccharides.

The glycosaminoglycan hyaluronan is involved in a diverse range of physiological and diseases processes and comprises repeated disaccharide units of N-acetyl-d-glucosamine (GlcNAc) and d-glucuronic acid (GlcA). A molecular description of the solution conformation of HA is required to account for this biology, which is best attained using nuclear magnetic resonance (NMR). NMR studies of the polymer, however, are frustrated by resonance overlap arising from the highly degenerate structure. In contrast, end-effects in oligosaccharides can produce some chemical shift dispersion, giving the possibility that their conformational properties can be measured and extrapolated to models of the polymer. We report the complete resolution and assignment of (1)H, (13)C and (15)N nuclei in hyaluronan oligosaccharides with seven different naturally occurring terminal rings. At 900MHz, all (1)H nuclei in the hexasaccharide GlcA-beta-(1-->3)-GlcNAc-beta-(1-->4)-GlcA-beta-(1-->3)-GlcNAc-beta-(1-->4)-GlcA-beta-(1-->3)-GlcNAc-OH were uniquely resolved and the two central rings were found to be a good model for the polymer environment. These assignments now allow resolved, unambiguous structural restraints to be acquired on this oligosaccharide and extrapolated to models for the solution conformation of the polymer.

NMR spectra of oligosaccharides at ultra-high field (900 MHz) have better resolution than expected due to favourable molecular tumbling.

Nuclear magnetic resonance (NMR) remains the most promising technique for acquiring atomic-resolution information in complex carbohydrates. Significant obstacles to the acquisition of such data are the poor chemical-shift dispersion and artifacts resultant from their degenerate chemical structures. The recent development of ultra-high-field NMR (at 900 MHz and beyond) gives new potential to overcome these problems, as we demonstrate on a hexasaccharide of the highly repetitive glycosaminoglycan hyaluronan. At 900 MHz, the expected increase in spectral dispersion due to higher resonance frequencies and reduction in strong coupling-associated distortions are observed. In addition, the fortuitous molecular tumbling rate of oligosaccharides results in longer T2-values that further significantly enhances resolution, an effect not available to proteins. Combined, the resolution enhancement can be as much as twofold relative to 600 MHz, allowing all 1H-resonances in the hexasaccharide to be unambiguously assigned using standard natural-abundance experiments. The use of ultra-high-field spectrometers is clearly advantageous and promises a new and exciting era in carbohydrate structural biology.

Enzymatic and chemical methods for the generation of pure hyaluronan oligosaccharides with both odd and even numbers of monosaccharide units.

Hyaluronan oligosaccharides display physiological activities not associated with the polymer and are widely used to characterize hyaluronan-binding proteins. They can also be used as biocompatible starting blocks for chemical derivatization. Here we present methods for generating milligram quantities of unusual odd- and even-numbered oligosaccharides, greatly increasing the diversity of reagents for use in such studies. These methods are based upon protocols from the 1960s, at which time it was very difficult to assess the stereochemical purity of the products. To address this, products were analyzed with modern high-field nuclear magnetic resonance spectroscopy. Alkaline beta-elimination conditions previously used to remove reducing-terminal N-acetylglucosamine residues in fact introduce a significant ( approximately 30%) level of stereoisomerism in the products by alkali-catalyzed keto-enol tautomerizations. Milder alkaline conditions were used to overcome this problem, reducing the contamination to <5%. The elimination by-products from this reaction were isolated and characterized, allowing the mechanism of alkaline degradation of hyaluronan to be investigated for the first time. beta-Glucuronidase was used to remove nonreducing-terminal glucuronic acid residues from oligosaccharides. Odd-numbered oligosaccharides with terminal glucuronic acid residues isolated from hyaluronidase digests are shown to originate from acid-catalyzed acetal hydrolysis during boiling denaturation and also have significant levels of stereochemical impurities.

Hyaluronan: the local solution conformation determined by NMR and computer modeling is close to a contracted left-handed 4-fold helix.

The polysaccharide hyaluronan (HA) is a ubiquitous component of the vertebrate extracellular matrix with diverse physiological roles from space-filling to acting as a scaffold for other macromolecules. The molecular interactions responsible for these solution properties have been the subject of much debate and, primarily due to the lack of residue-specific experimental data, no consensus model for the three-dimensional conformation nor dynamics of HA in solution has emerged. Here, the solution conformation of HA is investigated using molecular dynamics (MD) simulations and high-field nuclear magnetic resonance (NMR). In contrast to previous studies, MD simulations incorporated explicit water molecules and sodium ions, while NMR experiments utilized (15)N-enriched oligosaccharides to allow residue-specific information to be obtained. The resultant average conformation is predicted to be almost a contracted left-handed 4-fold helix; i.e. similar to that observed for sodium hyaluronate fibers by X-ray diffraction, but with the acetamido side-chain trans to H(2). The glycosidic linkages and acetamido side-chains are predicted to have standard deviation rotations of 13 degrees and 18 degrees around their mean conformations in free solution, respectively, and are not observed to be stabilized by strong intramolecular hydrogen bonds as X-ray fiber diffraction refinements describe for the solid-state. Rather, weak and transient hydrogen bonds that are in rapid interchange with solvent molecules are predicted. These predictions are quantitatively consistent with demanding residue-specific NMR data and correspond to an HA molecule that is rod-like as an oligosaccharide and behaves as a stiffened random coil at large molecular mass, in close agreement with previous hydrodynamic observations. This new description of the solution conformation of HA is consistent with all available experimental data and accounts for its viscoelastic space-filling properties. This representation can be used as a basis for modeling the association between HA and proteins, which will elucidate important aspects of extracellular matrix assembly.