Andrew Almond

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.

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.

N-Acetylated amino sugars: the dependence of NMR 3J(HNH2)-couplings on conformation, dynamics and solvent.

Density functional theory (DFT) calculations are used to determine the angular dependence of scalar couplings (Karplus equation) for N-acetylated amino sugars and an equation derived that can estimate bond libration from experimental measurements.

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

Hyaluronan: The absence of amide-carboxylate hydrogen bonds and the chain conformation in aqueous solution are incompatible with stable secondary and tertiary structure models

Contradictory descriptions for the aqueous solution conformation of the glycosaminoglycan hyaluronan (HA) exist in the literature. According to hydrodynamic and simulation data, HA molecules are stiffened by a rapidly interchanging network of transient hydrogen bonds at the local level and do not significantly associate at the global level. In marked contrast, models derived from NMR data suggest that the secondary structure involves persistent hydrogen bonds and that strong associations between chains can occur to form vast stable tertiary structures. These models require an extended 2-fold helical conformation of the HA chain and specific hydrogen bonds between amide and carboxylate groups. To test these descriptions, we have used 15N-labelled oligosaccharides and high-field NMR to measure pertinent properties of the acetamido group. The amide proton chemical shift perturbation and carboxylate group pK(a) value are inconsistent with a highly populated hydrogen bond between the amide and carboxylate groups. Amide proton temperature coefficients and chemical exchange rates confirm this conclusion. Comparison of oligomer properties with polymeric HA indicates that there is no discernible difference in amide proton environment between the centre of octasaccharides and the polymer, inconsistent with the formation of tertiary structures. A [1H-1H-15N] NOESY-HSQC (heteronuclear single-quantum correlation) spectrum recorded on an HA octasaccharide revealed that amide groups in the centre are in a trans orientation and that the average solution conformation is not an extended 2-fold helix. Therefore the two key aspects of the secondary and tertiary structure models are unlikely to be correct. Rather, these new NMR data agree with descriptions from hydrodynamic and simulations data.

A new route to carbohydrate secondary and tertiary structure using Raman spectroscopy and Raman optical activity.

The structural characterization of carbohydrate polymers is important for understanding their functions and behavior. However, mainstream structural biology tools are not applicable to many carbohydrate polymers, particularly at physiological concentrations. We report Raman and Raman optical activity spectra of hyaluronan polymer, the hyaluronan tetramer building block, and the two monosaccharide components glucuronic acid and N-acetylglucosamine and identify marker bands corresponding to primary and secondary structure in glycosaminoglycans. Furthermore, we show that the hyaluronan polymer does not adopt tertiary structure under near-physiological conditions, confirming a proposed model of hyaluronan structural organization.

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.

Pd(II)-mediated assembly of porphyrin channels in bilayer membranes.

A membrane-spanning bis(meso-3-pyridyl) porphyrin 1 has been synthesized, embedded in EYPC vesicles, and upon Pd(II) addition has been shown to form ionophores that allow the passage of anionic 5/6-carboxyfluorescein through membranes. The geometric matching of bis(meso-3-pyridyl) porphyrin 1 and trans-Pd(II) was designed to give a cyclic porphyrin trimer [PdCl(2)(1)](3). However, solution-phase studies showed that PdCl(2)(PhCN)(2) cross linked 1 into linear oligomers at porphyrin concentrations above 10 mM, although the formation of cyclic species was inferred from studies at concentrations below 2 μM. Fluorescence titrations showed that embedding porphyrin 1 in bilayers greatly reduced its affinity for Pd(II), but the combination of porphyrin 1 and Pd(II) gave an ionophoric species that increased the rate of 5/6-carboxyfluorescein (5/6-CF) transit through the phospholipid bilayer 12-fold. A maximum in the 5/6-CF release rate was observed at a Pd(II) concentration of 4 μM, and the application of a solution-phase binding model to the membrane phase showed that this peak in ionophoric activity corresponded to the greatest extent of porphyrin oligomerization. Further studies suggested these Pd(II)/porphyrin oligomers transported 5/6-CF via a channel mechanism.

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.