Benedict M. Sattelle

Free energy landscapes of iduronic acid and related monosaccharides

The pyranose ring of L-iduronic acid (IdoA), a major constituent of the anticoagulant heparin, is an equilibrium of multiple ring puckers that have evaded quantification by experiment or computation. In order to resolve this enigma, we have calculated the free energy landscape of IdoA and two related monosaccharides from extensive microsecond simulations. After establishing that the simulated puckers had reached equilibrium, hypotheses were confirmed that (a) IdoA (1)C(4)- and (4)C(1)-chair conformations exchange on the microsecond time scale, (b) C5 epimerization leads to a (4)C(1)-chair, and (c) IdoA 2-O-sulfation (IdoA2S) stabilizes the (1)C(4) conformer. The IdoA and IdoA2S (1)C(4) conformers were isoenergetic and computed to be 0.9 and 2.6 kcal mol(-1) lower in free energy than their respective (4)C(1)-chair conformations. The simulations also predicted that the IdoA (2)S(O)-skew-boat was less populated than previously thought. Novel chemical synthesis and ultra-high-field NMR supported these observations, but slight discrepancies in observed and predicted NMR vicinal couplings implied that the simulation overestimated the population of the IdoA (4)C(1)-chair with respect to (1)C(4)-chair due to small force field inaccuracies that only manifest in long simulations. These free-energy calculations drive improvements in computational methods and provide a novel route to carbohydrate mimetic biomaterials and pharmaceuticals.

Splice-Variant- and Stage-Specific RNA Editing of the Drosophila GABA Receptor Modulates Agonist Potency

The molecular diversity of many gene products functioning in the nervous system is enhanced by alternative splicing and adenosine-to-inosine editing of pre-mRNA. Using RDL, a Drosophila melanogaster GABA-gated ion channel, we examined the functional impact of RNA editing at several sites along with alternative splicing of more than one exon. We show that alternative splicing and RNA editing have a combined influence on the potency of the neurotransmitter GABA, and the editing isoforms detected in vivo span the entire functional range of potencies seen for all possible edit variants expressed in Xenopus laevis oocytes. The extent of RNA editing is developmentally regulated and can also be linked to the choice of alternative exons. These results provide insights into how the rich diversity of signaling necessary for complex brain function can be achieved by relatively few genes.

Comparative pharmacology and computational modelling yield insights into allosteric modulation of human α7 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.

Is N-acetyl-d-glucosamine a rigid 4C1 chair?

Understanding microsecond-timescale dynamics is crucial to establish three-dimensional (3D) structure–activity relationships in sugars but has been intractable to experiments and simulations. As a consequence, whether arguably the most important chemical scaffold in glycobiology, N-acetyl-d-glucosamine (GlcNAc), deviates from a rigid 4C1 chair is unknown. Here, conformer populations and exchange kinetics were quantified from the longest aqueous carbohydrate simulations to date (0.2 ms total) of GlcNAc, four derivatives from heparan sulfate and their methylglycosides. Unmodified GlcNAc took 3–5 μs to reach a conformational equilibrium, which comprised a metastable 4C1 chair that underwent 4C1 ↔ 1C4 transitions at a predicted forward rate of 0.8 μs−1 with an average 1C4-chair lifetime of 3 ns. These predictions agree with high-resolution crystallography and nuclear magnetic resonance but not with the hypothesis that GlcNAc is a rigid 4C1 chair, concluded from previous experimental analyses and non-aqueous modeling. The methylglycoside was calculated to have a slower forward rate (0.3 μs−1) and a more stable 4C1 conformer (0.2 kcal mol−1), suggesting that pivotal 3D intermediates (particularly 2SO, 1S5 and B2,5) increased in energy, and water was implicated as a major cause. Sulfonation (N-, 3-O and 6-O) significantly augmented this effect by blocking pseudorotation, but did not alter the rotational preferences of hydroyxl or hydroxymethyl groups. We therefore propose that GlcNAc undergoes puckering exchange that is dependent on polymerization and sulfo substituents. Our analyses, and 3D model of the equilibrium GlcNAc conformer in water, can be used as dictionary data and present new opportunities to rationally modify puckering and carbohydrate bioactivity, with diverse applications from improving crop yields to disease amelioration.

A 3D-structural model of unsulfated chondroitin from high-field NMR: 4-sulfation has little effect on backbone conformation

Graphical abstract

Does microsecond sugar ring flexing encode 3D-shape and bioactivity in the heparanome?

The biological information encoded in carbohydrate sequences dwarfs that of proteins and nucleic acids. Deciphering structure-function relationships in heparin and heparan sulfate (the heparanome) is further compounded by extreme sequence diversity, experimental difficulties, and the computational cost of rigorous modeling. Here we perform unbiased microsecond dynamics simulations of 11 heparanome oligosaccharides (55 microseconds total) to investigate the effect of sequence on 3D-structure and to underpin a coarse-grained model that is consistent with long-time scale experimentally validated atomic motions in water. Pyranose ring flexing (puckering) in 2-O-sulfo-α-l-iduronic acid, which underlies heparin-mediated anticoagulation, was modulated by polymerization (chain position and adjacent residues), which is supported by previous experiments. Furthermore, in coarse-grained simulations, inclusion of puckering was essential to predict macroscopic hydrodynamic properties of heparan sulfate chains containing hundreds of monosaccharaides. Our structural findings and model enable rational molecular design, and we propose that, in the heparanome, puckering, polymer 3D-shape, and bioactivity are inextricably linked.

Less is more when simulating unsulfated glycosaminoglycan 3D-structure: comparison of GLYCAM06/TIP3P, PM3-CARB1/TIP3P, and SCC-DFTB-D/TIP3P predictions with experiment.

The 3D‐structure of extracellular matrix glycosaminoglycans is central to function, but is currently poorly understood. Resolving this will provide insight into their heterogeneous biological roles and help to realize their significant therapeutic potential. Glycosaminoglycan chemical isoforms are too numerous to study experimentally and simulation provides a tractable alternative. However, best practice for accurate calculation of glycosaminoglycan 3D‐structure within biologically relevant nanosecond timescales is uncertain. Here, we evaluate the ability of three potentials to reproduce experimentally observed glycosaminoglycan monosaccharide puckering, disaccharide 3D‐conformation, and characteristic solvent interactions. Temporal dynamics of unsulfated chondroitin, chondroitin‐4‐sulfate, and hyaluronan β(1→3) disaccharides were simulated within TIP3P explicit solvent unrestrained for 20 ns using the GLYCAM06 force‐field and two semi‐empirical quantum mechanics methods, PM3‐CARB1 and SCC‐DFTB‐D (both within a hybrid QM/MM formalism). Comparison of calculated and experimental properties (vicinal couplings, nuclear Overhauser enhancements, and glycosidic linkage geometries) showed that the carbohydrate‐specific parameterization of PM3‐CARB1 imparted quantifiable benefits on monosaccharide puckering and that the SCC‐DFTB‐D method (including an empirical correction for dispersion) best modeled the effects of hexosamine 4‐sulfation. However, paradoxically, the most approximate approach (GLYCAM06/TIP3P) was the best at predicting monosaccharide puckering, 3D‐conformation, and solvent interactions. Our data contribute to the debate and emerging consensus on the relative performance of these levels of theory for biological molecules.

A Refined Model for the TSG-6 LINK Module in Complex with Hyaluronan: use of Defined Oligosaccharides to Probe Structure and Function.

Background: The polysaccharide hyaluronan is organized through interactions with the protein TSG-6 during inflammation and ovulation. Results: NMR spectroscopy on TSG-6 in the presence of defined sugars provided restraints that allowed modeling of a refined hyaluronan/TSG-6 complex. Conclusion: TSG-6 binding causes bending of hyaluronan that explains its condensation of this polysaccharide. Significance: This provides novel structural insights into protein-hyaluronan interactions. 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 13C-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.

Dependence of pyranose ring puckering on anomeric configuration: Methyl idopyranosides

In the aldohexopyranose idose, the unique presence of three axial ring hydroxyl groups causes considerable conformational flexibility, rendering it challenging to study experimentally and an excellent model for rationalizing the relationship between puckering and anomeric configuration. Puckering in methyl α- and β-L-idopyranosides was predicted from kinetically rigorous 10 μs simulations using GLYCAM11 and three explicit water models (TIP3P, TIP4P, and TIP4P-EW). In each case, computed pyranose ring three-bond (vicinal) (1)H-(1)H spin couplings ((3)J(H,H)) trended with NMR measurements. These values, calculated puckering exchange rates and free energies, were independent of the water model. The α- and β-anomers were (1)C(4) chairs for 85 and >99% of their respective trajectories and underwent (1)C(4)→(4)C(1) exchange at rates of 20 μs(-1) and 1 μs(-1). Computed α-anomer (1)C(4)↔(4)C(1) puckering rates depended on the exocyclic C6 substituent, comparing hydroxymethyl with carboxyl from previous work. The slower kinetics and restricted pseudorotational profile of the β-anomer were caused by water occupying a cavity bounded by the anomeric 1-O-methyl and the C6 hydroxymethyl groups. This finding rationalizes the different methyl α- and β-L-idopyranoside (3)J(H,H) values. Identifying a relationship between idopyranose anomeric configuration, microsecond puckering, and water structure facilitates engineering of biologically and commercially important derivatives and underpins deciphering presently elusive structure-function relationships in the glycome.

Calculating chemically accurate redox potentials for engineered flavoproteins from classical molecular dynamics free energy simulations.

The tricyclic isoalloxazine nucleus of the redox cofactors flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) acts as an electron sink in life-sustaining biological electron transfer (eT). The functional diversity of flavin-containing proteins (flavoproteins) transcends that of free flavins. A large body of experimental evidence attributes natural control of flavoprotein-mediated eT to tuning of the thermodynamic driving force by the protein environment. Understanding and engineering such modulation by the protein environment of the flavin redox potential (DeltaE(o)) is valuable in biotechnology and device design. In this study we employed classical molecular dynamics free energy simulations (MDFES), within a thermodynamic integration (TI) formalism, to calculate the change in FMN first reduction potential (DeltaDeltaE(o)(ox/sq)) imparted by 6 flavoprotein active site mutations. The combined performance of the AMBER ff03 (protein) and GAFF (cofactor) force fields was benchmarked against experimental data for mutations close to the isoalloxazine re- and si-faces that perturb the wild-type DeltaE(o)(ox/sq) value in Anabaena flavodoxin. The classical alchemical approach used in this study overestimates the magnitude of DeltaE(o) values, in common with other studies. Nevertheless, chemically accurate DeltaDeltaE(o) values--calculated to within 1 kcal mol(-1) of the experimental value--were obtained for five of the six mutations studied. We have shown that this approach is practical for quantitative in silico screening of the effect of mutations on the first reduction potential where experimental values and structural data are available for the wild-type flavoprotein. This approach promises to be useful as an integral part of future interdisciplinary strategies to engineer desired thermodynamic properties in flavoproteins of biotechnological interest.