Glycan flexibility: insights into nanosecond dynamics from a microsecond molecular dynamics simulation explaining an unusual nuclear Overhauser effect

Theoretical models of staurosporine and analogs uncover detailed structural information in biological solution

Staurosporine and its analogs (STA-analogs) are indolocarbazoles (ICZs) compounds able to inhibit kinase proteins in a non-specific way, while present antimicrobial and cytostatic properties. The knowledge of molecular features associated to the complexation, including the ligand shape in solution and thermodynamics of complexation, is substantial to the development of new bioactive ICZs with improved therapeutic properties. In this context, the empirical approach of GROMOS force field is able to accurately reproduce condensed phase physicochemical properties of molecular systems after parameterization. Hence, through parameterization under GROMOS force field and molecular simulations, we assessed STA-analogs dynamics in aqueous solution, as well as its interaction with water to probe conformational and structural features involved in complexation to therapeutic targets. The coexistence of multiple conformers observed in simulations, and confirmed by metadynamics calculations, expanding the conformational space knowledge of these ligands with potential implications in understanding the ligand conformational selection during complexation. Also, changes in availability to H-bonding concerning the different substituents and water can reflect on effects at complexation free energy due to variation at the desolvation energetic costs. Based on these results, we expect the obtained structural data provide systemic framework for rational chemical modification of STA-analogs.

Multiscale Simulations Examining Glycan Shield Effects on Drug Binding to Influenza Neuraminidase

Influenza neuraminidase is an important drug target. Glycans are present on neuraminidase and are generally considered to inhibit antibody binding via their glycan shield. In this work, we studied the effect of glycans on the binding kinetics of antiviral drugs to the influenza neuraminidase. We created all-atom in silico systems of influenza neuraminidase with experimentally derived glycoprofiles consisting of four systems with different glycan conformations and one system without glycans. Using Brownian dynamics simulations, we observe a two- to eightfold decrease in the rate of ligand binding to the primary binding site of neuraminidase due to the presence of glycans. These glycans are capable of covering much of the surface area of neuraminidase, and the ligand binding inhibition is derived from glycans sterically occluding the primary binding site on a neighboring monomer. Our work also indicates that drugs preferentially bind to the primary binding site (i.e., the active site) over the secondary binding site, and we propose a binding mechanism illustrating this. These results help illuminate the complex interplay between glycans and ligand binding on the influenza membrane protein neuraminidase.

Enantioselective Binding of Propranolol and Analogues Thereof to Cellobiohydrolase Cel7A

At the catalytic site for the hydrolysis of cellulose the enzyme cellobiohydrolase Cel7A binds the enantiomers of the adrenergic beta-blocker propranolol with different selectivity. Methyl-to-hydroxymethyl group modifications of propranolol, which result in higher affinity and improved selectivity, were herein studied by ¹ H,¹ H and ¹ H,¹³ C scalar spin-spin coupling constants as well as utilizing the nuclear Overhauser effect (NOE) in conjunction with molecular dynamics simulations of the ligands per se, which showed the presence of all-antiperiplanar conformations, except for the one containing a vicinal oxygen-oxygen arrangement governed by the gauche effect. For the ligand-protein complexes investigated by NMR spectroscopy using, inter alia, transferred NOESY and saturation-transfer difference (STD) NMR experiments the S-isomers were shown to bind with a higher affinity and a conformation similar to that preferred in solution, in contrast to the R-isomer. The fact that the S-form of the propranolol enantiomer is pre-arranged for binding to the protein is also observed for a crystal structure of dihydroxy-(S)-propranolol and Cel7A presented herein. Whereas the binding of propranolol is entropy driven, the complexation with the dihydroxy analogue is anticipated to be favored also by an enthalpic term, such as for its enantiomer, that is, dihydroxy-(R)-propranolol, because hydrogen-bond donation replaces the corresponding bonding from hydroxyl groups in glucosyl residues of the natural substrate. In addition to a favorable entropy component, albeit lesser in magnitude, this represents an effect of enthalpy-to-entropy compensation in ligand-protein interactions.

Predicting the Structures of Glycans, Glycoproteins, and Their Complexes

Complex carbohydrates are ubiquitous in nature, and together with proteins and nucleic acids they comprise the building blocks of life. But unlike proteins and nucleic acids, carbohydrates form nonlinear polymers, and they are not characterized by robust secondary or tertiary structures but rather by distributions of well-defined conformational states. Their molecular flexibility means that oligosaccharides are often refractory to crystallization, and nuclear magnetic resonance (NMR) spectroscopy augmented by molecular dynamics (MD) simulation is the leading method for their characterization in solution. The biological importance of carbohydrate-protein interactions, in organismal development as well as in disease, places urgency on the creation of innovative experimental and theoretical methods that can predict the specificity of such interactions and quantify their strengths. Additionally, the emerging realization that protein glycosylation impacts protein function and immunogenicity places the ability to define the mechanisms by which glycosylation impacts these features at the forefront of carbohydrate modeling. This review will discuss the relevant theoretical approaches to studying the three-dimensional structures of this fascinating class of molecules and interactions, with reference to the relevant experimental data and techniques that are key for validation of the theoretical predictions.

Conformation and dynamics at a flexible glycosidic linkage revealed by NMR spectroscopy and molecular dynamics simulations: analysis of β-L-Fucp-(1→6)-α-D-Glcp-OMe in water solution

The intrinsic flexibility of carbohydrates facilitates different 3D structures in response to altered environments. At glycosidic (1→6)-linkages, three torsion angles are variable, and herein the conformation and dynamics of β-L-Fucp-(1→6)-α-D-Glcp-OMe are investigated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. The disaccharide shows evidence of conformational averaging for the ψ and ω torsion angles, best explained by a four-state conformational distribution. Notably, there is a significant population of conformations having ψ = 85° (clinal) in addition to those having ψ = 180° (antiperiplanar). Moderate differences in (13)C R1 relaxation rates are found to be best explained by axially symmetric tumbling in combination with minor differences in librational motion for the two residues, whereas the isomerization motions are occurring too slowly to be contributing significantly to the observed relaxation rates. The MD simulation was found to give a reasonably good agreement with experiment, especially with respect to diffusive properties, among which the rotational anisotropy, D∥/D⊥, is found to be 2.35. The force field employed showed too narrow ω torsion angles in the gauche-trans and gauche-gauche states as well as overestimating the population of the gauche-trans conformer. This information can subsequently be used in directing parameter developments and emphasizes the need for refinement of force fields for (1→6)-linked carbohydrates.

Shaping up for structural glycomics: a predictive protocol for oligosaccharide conformational analysis applied to N-linked glycans

•

Aqueous 10 μs simulations of N-linked mannosyl cores and sialyl Lewis (sLe) antennae are validated.

•

Sequence dependent glycosidic linkage and pyranose ring μs motions are implicated in bioactivity.

•

Stacked pyranoses in sLe^a and sLe^x are predicted to be atypically rigid on μs timescales.

•

In a 25 μs simulation of sLe^x, all known conformers were sampled within the initial 10 μs of dynamics.

•

Unbiased 10 μs simulations are proposed as a route to systematic and accurate glycomic 3D-analysis.

The human glycome comprises a vast untapped repository of 3D-structural information that holds the key to glycan recognition and a new era of rationally designed mimetic chemical probes, drugs, and biomaterials. Toward routine prediction of oligosaccharide conformational populations and exchange rates at thermodynamic equilibrium, we apply hardware-accelerated aqueous molecular dynamics to model μs motions in N-glycans that underpin inflammation and immunity. In 10 μs simulations, conformational equilibria of mannosyl cores, sialyl Lewis (sLe) antennae, and constituent sub-sequences agreed with prior refinements (X-ray and NMR). Glycosidic linkage and pyranose ring flexing were affected by branching, linkage position, and secondary structure, implicating sequence dependent motions in glycomic functional diversity. Linkage and ring conformational transitions that have eluded precise quantification by experiment and conventional (ns) simulations were predicted to occur on μs timescales. All rings populated non-chair shapes and the stacked galactose and fucose pyranoses of sLe^a and sLe^x were rigidified, suggesting an exploitable 3D-signature of cell adhesion protein binding. Analyses of sLe^x dynamics over 25 μs revealed that only 10 μs were sufficient to explore all aqueous conformers. This simulation protocol, which yields conformational ensembles that are independent of initial 3D-structure, is proposed as a route to understanding oligosaccharide recognition and structure–activity relationships, toward development of carbohydrate-based novel chemical entities.

A perspective on the primary and three-dimensional structures of carbohydrates

Carbohydrates, in more biologically oriented areas referred to as glycans, constitute one of the four groups of biomolecules. The glycans, often present as glycoproteins or glycolipids, form highly complex structures. In mammals ten monosaccharides are utilized in building glycoconjugates in the form of oligo- (up to about a dozen monomers) and polysaccharides. Subsequent modifications and additions create a large number of different compounds. In bacteria, more than a hundred monosaccharides have been reported to be constituents of lipopolysaccharides, capsular polysaccharides, and exopolysaccharides. Thus, the number of polysaccharide structures possible to create is huge. NMR spectroscopy plays an essential part in elucidating the primary structure, that is, monosaccharide identity and ring size, anomeric configuration, linkage position, and sequence, of the sugar residues. The structural studies may also employ computational approaches for NMR chemical shift predictions (CASPER program). Once the components and sequence of sugar residues have been unraveled, the three-dimensional arrangement of the sugar residues relative to each other (conformation), their flexibility (transitions between and populations of conformational states), together with the dynamics (timescales) should be addressed. To shed light on these aspects we have utilized a combination of experimental liquid state NMR techniques together with molecular dynamics simulations. For the latter a molecular mechanics force field such as our CHARMM-based PARM22/SU01 has been used. The experimental NMR parameters acquired are typically (1)H,(1)H cross-relaxation rates (related to NOEs), (3)JCH and (3)JCCtrans-glycosidic coupling constants and (1)H,(13)C- and (1)H,(1)H-residual dipolar couplings. At a glycosidic linkage two torsion angles ϕ and ψ are defined and for 6-substituted residues also the ω torsion angle is required. Major conformers can be identified for which highly populated states are present. Thus, in many cases a well-defined albeit not rigid structure can be identified. However, on longer timescales, oligosaccharides must be considered as highly flexible molecules since also anti-conformations have been shown to exist with H-C-O-C torsion angles of ∼180°, compared to syn-conformations in which the protons at the carbon atoms forming the glycosidic linkage are in close proximity. The accessible conformational space governs possible interactions with proteins and both minor changes and significant alterations occur for the oligosaccharides in these interaction processes. Transferred NOE NMR experiments give information on the conformation of the glycan ligand when bound to the proteins whereas saturation transfer difference NMR experiments report on the carbohydrate part in contact with the protein. It is anticipated that the subtle differences in conformational preferences for glycan structures facilitate a means to regulate biochemical processes in different environments. Further developments in the analysis of glycan structure and in particular its role in interactions with other molecules, will lead to clarifications of the importance of structure in biochemical regulation processes essential to health and disease.

Density functional conformational study of 2-O-sulfated 3,6 anhydro-α-D-galactose and of neo-κ- and ι-carrabiose molecules in gas phase and water

We examined the conformational preferences of the 2-O-sulfated-3,6-α-D-anhydrogalactose (compound I) and two 1,3 linked disaccharides constituting-κ or ι-carrageenans using density functional and ab initio methods in gas phase and aqueous solution. Systematic modifications of two torsion angles leading to 324 and 144 starting geometries for the compound I and each disaccharide were used to generate adiabatic maps using B3LYP/6-31G(d). The lower energy conformers were then fully optimized using B3LYP, B3PW91 and MP2 with several basis sets. Overall, we discuss the impact of full relaxation on the energy and structure of the dominant conformations, present the performance comparison with previous molecular mechanics calculations if available, and determine whether our results are impacted, when polarization and diffuse functions are added to the 6-31G(d) basis set, or when the MP2 level of theory is used.

Conformational flexibility of N-glycans in solution studied by REMD simulations

Protein-glycan recognition regulates a wide range of biological and pathogenic processes. Conformational diversity of glycans in solution is apparently incompatible with specific binding to their receptor proteins. One possibility is that among the different conformational states of a glycan, only one conformer is utilized for specific binding to a protein. However, the labile nature of glycans makes characterizing their conformational states a challenging issue. All-atom molecular dynamics (MD) simulations provide the atomic details of glycan structures in solution, but fairly extensive sampling is required for simulating the transitions between rotameric states. This difficulty limits application of conventional MD simulations to small fragments like di- and tri-saccharides. Replica-exchange molecular dynamics (REMD) simulation, with extensive sampling of structures in solution, provides a valuable way to identify a family of glycan conformers. This article reviews recent REMD simulations of glycans carried out by us or other research groups and provides new insights into the conformational equilibria of N-glycans and their alteration by chemical modification. We also emphasize the importance of statistical averaging over the multiple conformers of glycans for comparing simulation results with experimental observables. The results support the concept of "conformer selection" in protein-glycan recognition.

Marine polysaccharides in microencapsulation and application to aquaculture: "from sea to sea"

This review's main objective is to discuss some physico-chemical features of polysaccharides as intrinsic determinants for the supramolecular structures that can efficiently provide encapsulation of drugs and other biological entities. Thus, the general characteristics of some basic polysaccharides are outlined in terms of their conformational, dynamic and thermodynamic properties. The analysis of some polysaccharide gelling properties is also provided, including the peculiarity of the charged polysaccharides. Then, the way the basic physical chemistry of polymer self-assembly is made in practice through the laboratory methods is highlighted. A description of the several literature procedures used to influence molecular interactions into the macroscopic goal of the encapsulation is given with an attempt at classification. Finally, a practical case study of specific interest, the use of marine polysaccharide matrices for encapsulation of vaccines in aquaculture, is reported.

Dynamic aspects of antibody:oligosaccharide complexes characterized by molecular dynamics simulations and saturation transfer difference nuclear magnetic resonance

Carbohydrates are likely to maintain significant conformational flexibility in antibody (Ab):carbohydrate complexes. As demonstrated herein for the protective monoclonal Ab (mAb) F22-4 recognizing the Shigella flexneri 2a O-antigen (O-Ag) and numerous synthetic oligosaccharide fragments thereof, the combination of molecular dynamics simulations and nuclear magnetic resonance saturation transfer difference experiments, supported by physicochemical analysis, allows us to determine the binding epitope and its various contributions to affinity without using any modified oligosaccharides. Moreover, the methods used provide insights into ligand flexibility in the complex, thus enabling a better understanding of the Ab affinities observed for a representative set of synthetic O-Ag fragments. Additionally, these complementary pieces of information give evidence to the ability of the studied mAb to recognize internal as well as terminal epitopes of its cognate polysaccharide antigen. Hence, we show that an appropriate combination of computational and experimental methods provides a basis to explore carbohydrate functional mimicry and receptor binding. The strategy may facilitate the design of either ligands or carbohydrate recognition domains, according to needed improvements of the natural carbohydrate:receptor properties.

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.