Next Article in Journal
CaDHN4, a Salt and Cold Stress-Responsive Dehydrin Gene from Pepper Decreases Abscisic Acid Sensitivity in Arabidopsis
Previous Article in Journal
Cell-Based Mechanosensation, Epigenetics, and Non-Coding RNAs in Progression of Cardiac Fibrosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Conformations of Di-, Tri-, and Tetra-α-(2→8)-Linked Sialic Acid from NMR Spectroscopy and MD Simulations

1
Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences Vienna, 1190 Vienna, Austria
2
Department of Chemistry, University of Natural Resources and Life Sciences Vienna, 1190 Vienna, Austria
*
Author to whom correspondence should be addressed.
Current address: Department of Chemistry, Yale University, New Haven, CT 06520-8107, USA.
Int. J. Mol. Sci. 2020, 21(1), 30; https://doi.org/10.3390/ijms21010030
Submission received: 12 November 2019 / Revised: 12 December 2019 / Accepted: 14 December 2019 / Published: 19 December 2019
(This article belongs to the Section Molecular Biophysics)

Abstract

:
By using molecular dynamics simulations with an efficient enhanced sampling technique and in combination with nuclear magnetic resonance (NMR) spectroscopy quantitative structural information on α -2,8-linked sialic acids is presented. We used a bottom-up approach to obtain a set of larger ensembles for tetra- and deca-sialic acid from model dimer and trimer systems that are in agreement with the available J-coupling constants and nuclear Overhauser effects. The molecular dynamic (MD) simulations with enhanced sampling are used to validate the force field used in this study for its further use. This empowered us to couple NMR observables in the MD framework via J-coupling and distance restraining simulations to obtain conformations that are supported by experimental data. We used these conformations in thermodynamic integration and one-step perturbation simulations to calculate the free-energy of suggested helical conformations. This study brings most of the available NMR experiments together and supplies information to resolve the conflict on the structures of poly- α -2,8-linked sialic acid.

1. Introduction

Neuraminic (sialic) acids are commonly found in the carbohydrate moiety of glycoproteins and glycolipids, typically at the outermost non-reducing ends. The chemical nature of this sugar, being negatively charged and hydrophilic, makes it suitable to perform key roles in molecular recognition events. In addition, sialic acids can form polymers with degrees of polymerization varying from 8 to 200, called polysialic acids (polySia) [1].
Polysialic acids play crucial roles in various steps of neural development, such as cell survival and cell migration of neural precursors, neuronal guidance, and synapse formation [2]. Modification of therapeutic recombinant proteins by polySia is of growing pharmaceutical relevance due to the significant increase in their serum lifetime [3,4]. However, the structural properties of these polymeric units are still under debate. First, Jennings et al. [5] found that at least a decamer of N-acetylated neuraminic acid (Neu5Ac) 10 was needed for binding to the group B polysaccharide specific antibody. As the size of the antibody is limited to the accommodation of 6–7 residues, α -2,8-linked polySia was considered to form a structural epitope. After that, structural studies on the conformation of oligomers of different degrees of polymerization suggested a flexible helical structure by early NMR [6,7]. One model was built from 3 J H H coupling constants by setting the appropriate torsional angles and using the antibody binding site as a target. This resulted in a helical model with 6 residues per turn [8], contradicting another helical polymer model, which was built by setting torsional angles satisfying coupling constants and suggests 3–4 residues per turn [9]. Yongye et al. [10] suggested another helical type by taking trisialic acid as a model through a combination of NMR and MD studies. Battistel et al. observed transient H-bonds in α -2,8-linked tetra-sialic acid from 3 J N C correlations at 263 K [11]. A contrasting NMR study based on 13 C-{ 1 H} NMR relaxation times and steady-state NOE’s suggests a random coil rather than a helical structure [12]. The actual conformation(s) are still unclear since there is no consensus on either the type of the helical conformation or on the random coil.
In the current work, we perform molecular simulations of di-, tri-, and tetra- α -2,8-linked sialic acid to study its conformational freedom and preferences. The simulations will be compared to and biased by experimental data derived from NMR experiments which were both conducted in this study and available in the literature, to characterize the conformational flexibility, within the limits set by the measured data. This study allowed us to further comment on the conformation on polymeric sialic acid with a decamer model not by using indirect information through setting torsional angles from coupling constants but by direct use of NOEs and coupling constants in the simulation leading to an ensemble satisfying the experimental findings.
Initially, we have studied the conformational landscapes of α -2,8-linked sialic acids in three different dimers, one trimer, and one tetramer (see Figure 1). The analyzed dimers were subsequently used to model larger sialic acid systems from tri-, tetra- to decasialic acid. In contrast to the α -linkages in polysialic acid, the reducing end, residue a, is present as the β -anomer which has a major impact on the NMR-characteristics of this unit and as well as of the preceding residue b. Therefore, for building tetra-sialic acid, we have used the information from its exact disaccharide units with respective anomericity. In the trimer system, the reducing end was kept in α -anomeric form to study the effect of chain prolongation on the dimeric linkage.

2. Results and Discussion

α -2,8-linked polySia has gained interest in recent drug design, especially in the development of specific antibodies and for glycoengineering of therapeutic proteins. However, the conformation of this polymer remains unsolved because of incomplete information on its structure. To shed light on the structure of α -2,8-linked tetra-sialic acid and the conformational preferences of its dimeric units we have performed various molecular simulations and validated these against NMR experiments. We previously used this approach to obtain 3D structures of various glycan units where the structure was unknown [13,14]. The 53A6 G L Y C [15] force field was used both to validate the sialic acid conformation and use it to propose conformational preferences of higher units. This force field was shown to reproduce experimental observables (NOE’s, 3 J coupling constants, etc.) well [16,17].
We have already shown in our recent work [17] that the LEUS simulations offer a more complete sampling than unbiased MD simulations. Therefore, we chose to compare the NMR findings with the LEUS simulations. For the three dimers, one trimer, and one tetramer, LEUS simulations were analyzed in terms of 3 J H 8 H 9 R / S , 3 J H 7 H 8 , and 3 J H 6 H 7 coupling constants, intramolecular hydrogen bonding occurrences, and glycosidic dihedral angle distributions. All ensemble averages were unbiased using Equation (1) (see Section 3).

2.1. Free-Energy Landscape

Inspection of the glycosidic free-energy maps (Figure 2) reveals that the α -2,8-linked systems show three significant minima regions (termed as A, B, and C) unlike other 1 n (where n is 2, 3, 4, and 6) linked disaccharides displaying four regions (see Ref. [18]). Although this behavior was shared among the studied dimer systems, the population of the states differs, depending on the stereochemistry of the reducing end. The corresponding states (A, B, C) have free-energy values of 0.0, 23.3, 13.67 kJ mol 1 for dimer1 α H and 6.6, 12.1, 0.0 kJ mol 1 for dimer2 β H, respectively. In region A, the minimum for the ϕ dihedral angle slightly shifts towards 30 in dimer2 and dimer3 compared to dimer1, for which it is centered around 60 (for ψ = 120 ). The additional degree of freedom, ω 8 in these linkages is almost freely rotating while ω 7 is restricted to 60 as can be seen in Figure 2. This predominance of the g + conformation of ω 7 is in agreement with NMR studies [6,7]. Although the right upper part of the ϕ , ψ free-energy landscape became more favorable compared to dimer1 α H, as expected, the effect of the end-group in dimer3 α C is relatively minor.

2.2. 3J Coupling Constants

For systems dimer1 α H and dimer2 β H, Table 1 shows the experimental 3 JHH coupling constants and the computed values from the LEUS simulations for each residue. For the trimer and dimer α C, the corresponding values are given in Table 2.
3 JH6H7 coupling constants for all systems was calculated to be lower than 1.5 Hz which were also reported in NMR. 3 JH7H8 couplings showed the maximum deviation for residue b in dimer1 β H system with 5.8 Hz deviation. This can be inspected from Figure 3. Although the energetically most favored region was calculated to be around 180 , the conformation around 60 was also highly populated and this resulted in a significant decrease in the average 3 JH7H8 coupling down to 3.7 Hz. In our LEUS simulations, the nonglycosylated 3 JH7H8 coupling of residue b in the dimers and residue c in the trimer is calculated as 3.7 and 2.5 Hz for α ending dimers (dimer1 α H, dimer3 α C) and 4.7 Hz for trimer1 α C. However, in the β ending dimer2 β H, it was calculated as 8.4 Hz. As can be seen from the shape of the Karplus relations of the 3 JH7H8 coupling in Figure 3 and Figure 4, the two extrema at 60 and 180 contribute values of 3 JH7H8 coupling lower than 2 Hz and higher than 9 Hz, respectively for the α ending systems, resulting in an average value of the two. This is in contrast to the high value of nonglycosylated 3 JH7H8 coupling in dimer2 β H which had a lower population around 60 , leading to a dominance of the 180 region, giving a high J-value. However, reported values from NMR suggest that 3 JH7H8 coupling on the nonglycosylated chain is ∼ 9.5 Hz in all α and β ending systems. Possibly, the nonglycosylated chains are slightly too flexible in our simulations of the α ending systems.
If we turn our attention to the 3 JH7H8 coupling constant on the glycosylated chain (residue a), all systems are in agreement with NMR values with a maximum of 0.3 Hz deviation. While α dimer1 α H and α dimer3 α C give 3 JH7H8 values of 1.8 and 1.0 Hz, β dimer2 β H gives 7.4 Hz, in close agreement with the experimental values. This 3 J coupling constant shows the main difference between having an α and β terminus, emphasizing the effect of the stereochemistry of the reducing end on the glycosidic dihedral angle preference. In Figure 2 and the colors of Figure 3 and Figure 4, this is reflected by the increased preference for conformations with ω 8 ≈ 180 .
Systems were largely in agreement with the experimental 3 J H 8 H 9 R / S values with a maximum deviation of 2.5 Hz for residue a in dimer2 β H. The average absolute deviation over eight 3 J H 8 H 9 R / S coupling constants in dimer1 α H and dimer2 β H amounts to 1.1 Hz. No significant difference between α and β -terminated systems in 3 J H 8 H 9 R / S values were found in the LEUS simulations which were also reported in NMR experiments.
To complement the J-coupling data of the tetramer, coupling constants of tetrasialic acid were derived from a published 850 MHz proton spectrum as well as from J-resolved 600 MHz experiments (see Figures S1 and S2 in the supplementary material). A high-order spin system of H9c H8c overlapping with H9b was cross-checked by spin simulation [20].
The full 2 J and 3 J coupling assignments from NMR experiments of the tetramer are reported in Table S1 in the supplementary material. The comparison of the 3 J H 6 H 7 , 3 J H 7 H 8 , and 3 J H 8 H 9 R / S coupling constants with LEUS simulations of the tetramer are represented in Table 3 and in Figures S3 and S4. 3 J H 6 H 7 couplings show maximum deviation in residue d with 1.1 Hz. 3 J H 7 H 8 coupling at the free, non-glycosylated end in residue d is calculated as 9.8 which is in agreement with the NMR findings. Only residue a with a β -terminus showed a 3.6 Hz deviation. If it is assumed to be similar to dimer2 β H, 3 J H 7 H 8 would have to be bigger than 7 Hz. The reason for not capturing the higher value might be due to strong interactions with the other residues, resulting in a different conformational preference of the tetramer or it can be due to a too pronounced sampling of the lower extreme of the Karplus curve for ω 8 around 60 . For the 3 J H 8 H 9 R / S coupling constants in the tetramer, the highest deviation is seen at the second residue (c) with 6.5 Hz deviation. NMR showed values of 5.9 and 4.1 for 3 J H 8 H 9 R and 3 J H 8 H 9 S while LEUS calculations gave a value at the two extrema of the Karplus relation. This might be an indication of poor sampling of one of the two conformations.

2.3. Hydrogen Bonding

The occurrence of H-bonds was reweighed to the unbiased ensemble for the LEUS simulations. H-bonds with occurrences larger than 2% are reported in Table 4 for dimer and trimer systems.
Hydrogen bonds between the oxygen of the carboxyl group and the hydroxyl group of either the free or the glycosylated side chain were found to be highly populated in these systems. Dimer2 β H showed a more prominent intermolecular hydrogen-bonding pattern, as was reported by Azurmendi et al. as well [21].

2.4. MD Simulations of Tetramer with NOE Restraining

The CBCANH and SOLEXSY experimental studies of Battistel et al. indicated that the tetramer shows an H-bond for each of the first three residues of the molecule between HN5-O8 at 263 K [11]. They suggested that conformations with this transient H-bond pattern can be in equilibrium with other conformations at higher temperatures. In our LEUS simulations of the tetramer, this inter-residue hydrogen bond was not prominent. To obtain the related conformation, we have performed two additional simulations in which we imposed instantaneous distance restraints derived from the NOE data at 263 K and 300 K. In addition to NOE’s, Battistel et al. reported 4 JH7-C2 coupling constants (also called W-couplings) via HSQMBC with 0.8 and 1.5 Hz for the second and the third residue, respectively. This spin coupling is measurable when C2 and H7 are in a W-shaped arrangement (C2-O6-C6-C7-H7 is quasi-coplanar) such that a C6-C7 conformation orienting H6 approximately anti to O7 is highly populated in solution. This long-range coupling was also reported as 0.8 Hz in the monomeric form of β -sialic acid by Klepach et al. [22]. We have used two additional dihedral restraining potentials for these two residues to sustain 4 JH7-C2 couplings by restraining the H7-C7-C6-H6 dihedral angle to approximately 90 .
Battistel et al. have generated structures based on three H-bond restraints and setting two H7-C7-C6-H6 dihedral angles at residues c and b at approximately 90 . By using only these observations they created models resembling a 2 4 and 1 4 helix. Their NOE’s were complementing the static 2 4 helix model with a higher percentage (71%) in which the exo-anomericity was not satisfied. Therefore, we aimed to build conformation(s) of oligosialic acids from the experimental data (instead of imposing experimental data on static models). For this aim, we have combined all experimental findings for a model with 2 dihedral angle restraints for the 4 J -couplings, 3 distance restraints for hydrogen bonding and 82 NOE-derived distance restraints. Simulations were performed at 263 K and 300 K.
The energetic barriers between possible ring conformers of sialic acid are relatively shallow. The free energy penalty of the ring conformational change from the most stable conformer 2 C 5 to the closest alternative conformations 4 , O B/ O S 3 and 4 S 2 was calculated to be less than 10 k B T [23]. As imposing restraints might disrupt the chair conformation, we have checked the ring puckering parameters and all sialic acid residues stayed at their energetically favored 2 C 5 chair conformation. Ring puckering parameters were calculated and are reported as ϕ vs. θ for all residues in Figure S5. Idealized ring puckering coordinate parameters are given in Table S4.
Our dynamical model resulted in compact conformations satisfying their NMR findings. These conformations are consistent with the exo-anomeric effect. The obtained set of conformations satisfied all NOE’s and H-bonds to within 0.1 Å (see Figure S6) and the two dihedral angles stayed at 70 with a total of 880 kJ/mol restraining energy in both the 263 K and 300 K simulations. We have clustered the resulting ensemble by using an RMSD matrix with a 0.15 nm cutoff which resulted in six clusters [24]. 97% of the ensemble was clustered into a single cluster and the hydrogen bonding analysis of this cluster reveals a similar inter-residue pattern as in our LEUS simulation. Inter-residue hydrogen bonds between successive residues were observed between HN5-O1A/B, HO7-O1A/B, and HO9-O6 for 95%, 90%, and 37% of simulation time, respectively. Among the rest of the clusters, only one of them, with 15 structures out of 8000, showed the successive HN5-O8 H-bonding pattern as suggested by Battistel et al. (see Figure S7). No significant differences were observed between 263 K and 300 K. Note that the experiments were performed at 263 K, but also at this lower temperature, the HN5-O8 H-bond was only rarely observed in our simulations, in spite of the restraints.

2.5. Helical Conformation?

Thermodynamic integration (TI) was used to compute the free-energy difference between conformation satisfying NOE’s and 4 J -couplings (state A) and conformations with three consecutive H-bonds that are observed in the NMR experiment [11] (state B). The initial structure for state A was taken from the end of the simulation with NOE restraining. Then, with the coupling parameter λ , three hydrogen bond restraints were added. By integrating over the derivative of the Hamiltonian with respect to λ , we obtain the free-energy difference ( Δ G A B ) as 81.5 kJ/mol. Then, with one step perturbation, we calculated the free energy of releasing the H-bond restraints for the conformations of state B by using Equation (8) as -15.2 kJ/mol (see Materials and Methods and Figure S8). Based on the free-energy calculations, the proposed helical structure with three consecutive hydrogen bonds between O8 and HN5 is unfavorable by about 65 kJ/mol. This 25 k B T shows that the conformation with three consecutive hydrogen bonds is not a thermodynamically accessible state for the tetramer.
To check if the energetic gain of the Hydrogen bonds may counteract the entropic loss of conformational freedom for larger oligomers, we have constructed a decameric structure as well. By using our free-energy landscapes, we have constructed an decamer of α -2,8-linked-polysialic acid by setting the glycosidic dihedral angles to their respective lowest free-energy states (A, B and C) in Figure 2. For each lowest energetic state, a set of structures was generated (81 decameric structures) by setting ϕ and ψ dihedral angles for the four linkages in the pentamer system and repeating this to a decamer. Among those 81 structures, the two non-clashing lowest energetic ones were selected and equilibrated to be used as the initial structure for 100 ns simulations with plain MD and with hydrogen bond restraints on the O8-HN5 distance. In the plain MD simulations, the O8-HN5 distances fluctuate in the range of 3.7–5.0 nm and conformations largely resemble random coils. In contrast, the restraint simulations showed different helical patterns. We have jointly clustered the configurations from plain MD and H-bond restrained simulations with a 0.25 nm cutoff. No overlap of the most observed conformations was observed. Of the first nine clusters, clusters 1–3 and 5–9 contained conformations that originated from the restrained simulations. Cluster 4 originated from the plain MD simulations, with all remaining conformations forming smaller undefined clusters. Figure 5 shows the central member structures of the first nine clusters. Helical patterns are observed in clusters 1, 3, 7 and 9 amounting to 37% of the restrained simulations. The plain MD simulations, however, confirm that helical structures with O8-HN5 hydrogen bonds are also not thermodynamically likely in decasialic acid. Along with the relaxation NMR study of Henderson et al., our results support a flexible conformation for polysialic acid. Another supporting NMR study is from Hanashima et al. [25] who found significant differences in the interresidue transferred NOE correlations of bound and free conformations of trisialic and hexasialic acid units suggesting distinct conformations. This shows that highly unfavorable conformations may still be stabilized when polysialic acid is bound to an antibody or within a glycoprotein environment.
Finally, we have modeled an antibody-decasialic acid complex by using the crystal structure of the single-chain variable fragment of an anti-polysialic acid antibody, which was in a complex with octasialic acid (see Figure S9). We created a trajectory from all the structures that belong to cluster 4 (which is the most populated cluster originating from plain MD simulations) and aligned this to the octasialic acid in the crystal structure. It shows that the alternative ensemble of decasialic acid obtained in this study can fit well to the crystal structure of the anti-polysialic acid antibody.

3. Materials and Methods

3.1. NMR Experiments

Tetrasialic acid was obtained from Hycultec (Cat No C-7005-10) and measured in D 2 O at 300 K at a concentration of 3.5 mM. For NMR experiments the Bruker Topspin 3.5 PL 6 Software Suite was used on a Bruker Avance III console with a 600 MHz magnet. Assignments of proton and carbon chemical shifts and the 850 MHz 1 H NMR spectrum were available from previous work [26].
For JRES experiments, a predefined Bruker parameter set (PROF_JRES) that employs the jresgpprqf pulse program with a a time domain for F2 of 32k and F1 200 data points (FIDRES F2 = 0.3 Hz and F1 = 0.8 Hz) was used. For HSQC, the pulseprogram hsqcedetgp and for COSY cosygpqf was employed.

3.2. MD Simulation Settings

In this study, we have used the six sialic acid systems represented in Figure 1. These involved three dimeric systems with different anomericity and end-groups at the reducing end, one trimer, one tetramer, and one decamer. All MD simulations were performed using the GROMOS11 biomolecular simulation package (http://www.gromos.net) [27] and the 53A6 G L Y C carbohydrate parameter set [15]. We have previously introduced minor modifications to the original 53A6 G L Y C carbohydrate parameter set for the sialic acid building block to ensure consistency with the rest of the force field [18]. Initial structures of the studied units were modeled in the molecular operating environment (MOE) [28] by setting their glycosidic dihedral angles to their respective free-energy minima which have been previously reported [18]. Short energy minimization was performed using the steepest-descent algorithm in a vacuum. The compounds were placed in a periodic cubic water box with simple point charge (SPC) water [29] molecules and initialized with a 1.4 nm minimum distance of the solute to the box walls. With position restraints on the solute atoms, the system was further relaxed by a steepest descent minimization. Then, the systems were equilibrated with initial random velocities generated from a Maxwell–Boltzmann distribution at 60 K and heated up to 300 K in five discrete steps. While heating up the system, position restraints on the solute atoms were reduced from 2.5 × 10 4 to 0.0 kJ mol 1 nm 2 .
The production simulations were performed at a constant temperature of 300 K and a constant pressure of 1 atm using a weak coupling scheme [30] for both temperature and pressure with coupling times τ T = 0.1 ps and τ P = 0.5 ps, respectively with isothermal compressibility of 4.575 × 10 4 kJ 1 mol nm 3 . Newton’s equations of motion were integrated using the leapfrog scheme [31] with a time step of 2 fs. The SHAKE algorithm [32] was used to maintain the bond lengths at their optimal values. Long-range electrostatic interactions beyond a cutoff of 1.4 nm were truncated and approximated by a generalized reaction field [33] with a relative dielectric permittivity of 61 [34]. Nonbonded interactions up to a distance of 0.8 nm, were computed at every time step using a pairlist [35] that was updated every 5 steps. Interactions up to 1.4 nm, were computed at pairlist updates and kept constant in between.
The GROMOS++ software [36] is used for time series analysis. A geometrical criterion was used to identify hydrogen bonds if a hydrogen-acceptor distance is smaller than 0.25 nm and the donor-hydrogen-acceptor angle is larger than 135 .

3.3. Creating Biased Potentials with Local Elevation and Sampling with Umbrella Sampling

For each system, unbiased MD simulations were carried out for 100 ns after equilibration. In addition to unbiased simulations, an enhanced sampling method, local elevation with umbrella sampling (LEUS) [37,38] was applied. For the local elevation potential build-up, the glycosidic dihedral angles of the systems were used. For a detailed description of the methodology, see our previous work in Ref. [18]. In short, in the LEUS method dihedral angles are binned in N g = 36 bins, and a biasing potential width of σ = 360 / N g was used with a force constant increment of c = 0.005 kJ mol 1 . In the current work, the studied systems have four conformational degrees of freedom along their 2→8 glycosidic linkage. To ensure a near-to-complete sampling along the glycosidic linkage, two different two-dimensional potentials with ϕ , ψ and ω 8 , ω 7 were used. For definitions of the dihedral angles see the caption of Figure 1. Local elevation potentials were only built up for dimer systems (dimers 1–3). Dimer2 and dimer3 potentials were used for umbrella sampling for the tetramer and dimer3 for the trimer. Two different end groups used in the dimer systems; OH group (dimer 1 and dimer 2) and OMe group (dimer 3), to compare the effect of the end group. In addition, different combinations of 3D and 4D potentials for the ϕ , ψ , ω 8 and ω 7 were checked but those attempts did not offer additional sampling. Although in principle a long four-dimensional build up along all four degrees of freedom ( ϕ , ψ , ω 8 , and ω 7 ) offers a complete sampling, it was not as efficient as using two different 2D potentials. Figure 6 shows the free energy landscape obtained after a four-dimensional build up of 200 and 400 ns for dimer3. Only after 400 ns the same coverage is achieved as with two 2D LE potentials (Figure 2C).
Therefore, after the optimization of the creation of the LE potentials two 2D LE potentials are built for the glycosidic linkage with t L E = 100 ns. These potentials are created for all the three dimer systems as they differ in type and anomericity of the end groups. In the US phase, the LE biased potentials were frozen and sampling was applied by using both 2D potentials. Then, umbrella sampling was applied to those potentials by saving trajectories every 0.1 ps for 100 ns to achieve statistical efficiency as discussed in Ref. [37].
The unbiased probability of any property Q can be obtained from the LEUS (biased) simulations through reweighing:
P ( Q ) = δ ( Q Q ) exp [ U L E U S ( Q ) / k B T ] exp [ U L E U S ( Q ) / k B T ] ,
where indicates an ensemble average of the biased LEUS simulation, U LEUS (Q) is the biasing energy at a particular value of Q, δ is the Kronecker delta function, kB is the Boltzmann constant and T is the absolute temperature. The corresponding free energies can be obtained from the calculated probabilities,
G ( Q ) = k B T ln P ( Q ) .

3.4. Analysis

For dimer systems (1–3) two free-energy maps, G( ϕ , ψ ) and G( ω 8 , ω 7 ) were created from the LEUS simulations after reweighing of the biased energy with Equations (1) and (2). The global minimum of each map represents the lowest free-energy conformation which is set to 0 kJ/mol and the colormap is drawn using a 5 kJ/mol contour. A more detailed explanation for the construction of the free-energy maps can be found in the methods section of Ref. [18].

3.5. 3J-coupling Constants and NOE Calculations

Simulations are compared with NOE data and 3 J-coupling constants. Aliphatic carbon atoms are treated as united atoms in the GROMOS force field. Therefore, virtual atomic positions for prochiral CH2 (C3 and C9), for CH (C4, C5, C6, C7, and C8) and pseudo atomic positions for CH3 (C5A) were used to calculate interproton distances. Upper bound corrections are applied by centre averaging approach [39] and reported in Tables S2 and S3 (using for prochiral CH2 0.09 nm and for CH3 0.1 nm; no corrections are added for CH and explicit HN hydrogen atoms). For the NOE analysis, averaging is performed as < r 6 > 1 / 6 [40].
3 J-coupling constants can be related to a torsional angle through the Karplus relation. Coupling constants related to glycosidic dihedral angles are 3 J H 8 H 9 R / S , 3 J H 7 H 8 , and 3 J H 6 H 7 corresponding to ω 9, ω 8, and ω 7 dihedral angles, respectively (Figure 1). 3 J H 6 H 7 and 3 J H 7 H 8 were calculated from MD and LEUS simulations using the Haasnoot Equation (3) which involves a modification to the Karplus relation, to include the effect of electronegativity of neighboring groups [41].
3 J H C C H = P 1 c o s 2 θ + P 2 c o s θ + P 3 + Δ χ i [ P 4 + P 5 c o s 2 ( s i θ + P 6 | Δ χ i | ) ] .
Here, the sum runs over the four substituents, Δ χ i = Δ χ α + P 7 Δ χ β and s i is the direction coefficient which is +1 for “positive” and −1 for “negative” substituents defined by their orientation relative to the attached α or β atom. Parameters P1–P7 used in Haasnoot equation are: 13.7, 0.73 , 0, 0.56, 2.47 , 16.9, and 0.14, respectively.
For prochiral protons, 3 J H 8 H 9 R and 3 J H 8 H 9 S values were calculated from the following equation which was derived by DFT methods [42]:
3 J H 8 H 9 R = 5.08 + 0.47 c o s ω + 0.9 s i n ω 0.12 c o s 2 ω + 4.86 s i n 2 ω , ω = ω 9 3 J H 8 H 9 S = 4.92 1.29 c o s ω + 0.05 s i n ω + 4.58 c o s 2 ω + 0.07 s i n 2 ω , ω = ω 9 120 .

3.6. Simulations with Distance Restraints

To study the possible conformations under the condition that available experimental data is fulfilled, we performed additional simulations with full harmonic distance restraints. In the restraint tetramer simulations, 82 instantaneous NOE restraints were applied [43] with a force constant of 500 kJ mol 1 nm 2 and 3 instantaneous restraints on HN5-O8 pairs with a force constant of 2500 kJ mol 1 nm 2 to satisfy hydrogen bonding. The NOE distance restraints used are listed in Tables S2 and S3 in the Supplementary Materials. Additionally, 2 instantaneous dihedral angle restraints were included with a force constant of 100 kJ mol 1 degrees 2 for residues c and b at H7-C7-C6-H6 at 90 to fulfill 4 J H 7 C 2 coupling. Since the hydrogen bonding evidence and the NOE data experiments [11] were conducted at 263 K we simulated the tetramer at 263 K. An additional simulation with the same parameters at a temperature of 300 K was run to check the system at a temperature closer to physiological conditions. The conformations that are observed in the restraint simulations were further used as an initial structure for the free-energy calculation. Thereby, to the best of our knowledge, this is the first time where a quantitative conformational analysis was made for proposed polysialic acid structures and models.

3.7. Thermodynamic Integration and One-Step Perturbation

The change from a conformation satisfying experimental NOEs and 4 J-coupling constants (state A) to a state B in which additionally three H-bonds are enforced, is performed with thermodynamic integration (TI). The path from state A to B is defined by a scaling parameter λ , where the Hamiltonian ( H ) at λ = 0 describes state A and at λ = 1 describes state B. The free energy difference between state A and B can be extracted by numerical integration of H / λ λ as a function of λ with the following equation [44]:
Δ G A B = G B G A = 0 1 H λ λ d λ ,
where the angular brackets denote the ensemble average of the derivative of H ( λ ) with respect to λ . In state A, the force constant on the three hydrogen-bond restraints is 0, but the derivative of H still depends on the unrestrained distance between the atoms. This leads to instabilities in λ -derivatives during the perturbation. This is avoided by using a soft bond potential energy term which is introduced with a softness term (S( r , λ )) [45]:
U ( r , λ ) = 1 2 ( 1 λ ) k A S A ( r , λ ) + λ k B S B ( r , 1 λ ) r ( 1 λ ) r 0 A + λ r 0 B 2 ,
where the softness term for state X is
S X ( r , λ ) = 1 + α λ ( r r 0 X ) 2 ,
with α a unitless softness parameter set to 250. The force constants k A and k B were set to be k A = 0 and k B = 5 × 10 4 kJ mol 1 nm 2 and the optimal bond length was set to r 0 A = r 0 B = 0.25   nm . First equidistant λ values were used for the change between states A and B. Then, to achieve a smooth transformation, additional λ points were introduced. In total 12 λ points were first equilibrated for 50 ps, followed by a production run of 2 ns for each λ point sequentially. Initial configuration at a certain λ point were taken from the final configuration of the previous λ value. The coordinates and energies were stored every 0.2 ps.
Then, one-step perturbation was applied on the restrained ensemble to remove the hydrogen bond restraints (state C). State C represents an unrestrained ensemble, which only covers the phase space in which the hydrogen bonds are observed. The free energy difference between state B and state C is calculated using the Zwanzig formula [46]:
Δ G B C = G C G B = k B T ln e ( H C H B ) / k B T B ,
where k B is the Boltzmann constant, T is the temperature, and the difference between the Hamiltonian of two states is the restraining energy ( H B H C = U r e s ). B represents the ensemble average over all configurations generated during the simulation at state B which was run for 100 ns.

4. Conclusions

We have studied the conformational freedom of the α -2,8-linkage between two subsequent sialic acid monomers with an enhanced conformational searching and sampling method, LEUS. Our simulations showed that the NMR observations can be satisfied with a set of conformations that are very similar to the ones observed in our LEUS simulations. The largest deviations in terms of 3 J -coupling constants are observed for the free non-glycosylated tail of the terminal residue in the C8-C9 sidearm, while the parameters determining the glycosidic linkage between two sialic acid units are excellently reproduced. Due to the extra degrees of freedom in the α -2,8-linkage, multiple conformations are accessible, as observed in the free-energy landscape of the dimers. However, the hydrogen bond that was previously suggested to support the helical conformation was found to be very transient in our restrained tetra-sialic acid simulations. The thermodynamic penalty of enforcing all three of these hydrogen bonds simultaneously was computed to be highly unfavorable, with about 25 k B T . While enforcing the hydrogen bonds in decasialic acid does lead to helical structures (up to 37%), such conformations are not observed in plain MD simulations of the same molecule. From the current work, we conclude that helical conformations, while possible, are unlikely to play a dominant role in free polysialic acid, but would need to be induced upon binding to a partner, requiring a considerable enthalpy-entropy compensation.

Supplementary Materials

Supplementary materials are available online at https://www.mdpi.com/1422-0067/21/1/30/s1, Figures S1 and S2: Expansion plot of the 600 MHz J-resolved spectrum of tetrasialic acid, Figures S3 and S4: θ vs. calculated 3 JH7H8, 3 JH8H9R and 3 JH8H9S couplings from LEUS simulations for tetramer β H, Figure S5: Ring conformation of the residues in the tetramer, Figure S6: NOE violations from restrained simulation of tetramer, Figure S7: Representative structures from clusters to illustrate hydrogen bonding patterns from NOE-restrained simulation of the tetramer, Figure S8: TI and OSP plot, Figure S9: Representation of the anti-polysialic acid crystal structure of single-chain variable fragment in complex with octasialic acid, Table S1: 2D NMR assignments for tetramer β H, Tables S2 and S3: Intra and Inter-residue NOE restraining values, Table S4: Cremer-Pople parameters for the idealized ring conformations.

Author Contributions

This study was initiated and conceptually drafted by A.T. and C.O. NMR Experiments were performed and analyzed by M.B. and P.K. Simulations were performed and analyzed by A.T. Interpretation of the data was performed by all authors. The first draft of the manuscript was written by A.T., the manuscript was finalized by all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Austrian Science Fund in the context of the Doctoral Program Biomolecular Technology of Proteins (BioToP; Grant W 1224).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

MDMolecular Dynamics
NMRNuclear Magnetic Resonance
NOENuclear Overhauser Effect
LEUSLocal Elevation Umbrella Sampling

References

  1. Sato, C.; Kitajima, K. Structural Analysis of Polysialic Acid. In Experimental Glycoscience; Springer: Tokyo, Japan, 2008; pp. 77–81. [Google Scholar]
  2. Rutishauser, U. Polysialic acid in the plasticity of the developing and adult vertebrate nervous system. Nat. Rev. Neurosci. 2008, 9, 26–35. [Google Scholar] [CrossRef] [PubMed]
  3. Lindhout, T.; Iqbal, U.; Willis, L.M.; Reid, A.N.; Li, J.; Liu, X.; Moreno, M.; Wakarchuk, W.W. Site-specific enzymatic polysialylation of therapeutic proteins using bacterial enzymes. Proc. Natl. Acad. Sci. USA 2011, 108, 7397–7402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Chen, C.; Constantinou, A.; Chester, K.A.; Vyas, B.; Canis, K.; Haslam, S.M.; Dell, A.; Epenetos, A.A.; Deonarain, M.P. Glycoengineering Approach to Half-Life Extension of Recombinant Biotherapeutics. Bioconjug. Chem. 2012, 23, 1524–1533. [Google Scholar] [CrossRef] [PubMed]
  5. Jennings, H.J.; Roy, R.; Michon, F. Determinant specificities of the groups B and C polysaccharides of Neisseria meningitidis. J. Immunol. 1985, 134, 2651–2657. [Google Scholar]
  6. Michon, F.; Brisson, J.R.; Jennings, H.J. Conformational Differences Between Linear α(2→8)-linked Homosialooligosaccharides and the Epitope of the Group B Meningococcal Polysaccharide. Biochemistry 1987, 26, 8399–8405. [Google Scholar] [CrossRef]
  7. Brisson, J.R.; Baumann, H.; Imberty, A.; Perez, S.; Jennings, H.J. Helical Epitope of the Group B Meningococcal α(2-8)-linked Sialic acid Polysaccharide. Biochemistry 1992, 31, 4996–5004. [Google Scholar] [CrossRef]
  8. Evans, S.V.; Sigurskjold, B.W.; Jennings, H.J.; Brisson, J.R.; To, R.; Tse, W.C.; Altman, E.; Weisgerber, C.; Young, N.M.; Frosch, M.; et al. Evidence for the Extended Helical Nature of Polysaccharide Epitopes. The 2.8 Å Resolution Structure and Thermodynamics of Ligand Binding of an Antigen Binding Fragment Specific for α-(2→8)-Polysialic Acid. Biochemistry 1995, 34, 6737–6744. [Google Scholar] [CrossRef]
  9. Yamasaki, R.; Bacon, B. Three-Dimensional Structural Analysis of the Group B Polysaccharide of Neisseria meningitidis 6275 by Two-Dimensional NMR: The Polysaccharide Is Suggested to Exist in Helical Conformations in Solution. Biochemistry 1991, 30, 851–857. [Google Scholar] [CrossRef]
  10. Yongye, A.B.; Gonzalez-Outeiriño, J.; Glushka, J.; Schultheis, V.; Woods, R.J. The conformational properties of methyl alpha-(2,8)-di/trisialosides and their N-acyl analogues: Implications for anti-Neisseria meningitidis B vaccine design. Biochemistry 2008, 47, 12493–12514. [Google Scholar] [CrossRef] [Green Version]
  11. Battistel, M.D.; Shangold, M.; Trinh, L.; Shiloach, J.; Freedberg, D.I. Evidence for helical structure in a tetramer of α2-8 sialic acid: Unveiling a structural antigen. J. Am. Chem. Soc. 2012, 134, 10717–10720. [Google Scholar] [CrossRef] [Green Version]
  12. Henderson, T.J.; Venable, R.M.; Egan, W. Conformational flexibility of the Group B meningococcal polysaccharide in solution. J. Am. Chem. Soc. 2003, 125, 2930–2939. [Google Scholar] [CrossRef] [PubMed]
  13. Montero-Morales, L.; Maresch, D.; Castilho, A.; Turupcu, A.; Ilieva, K.M.K.; Crescioli, S.; Karagiannis, S.S.N.; Lupinek, C.; Oostenbrink, C.; Altmann, F.; et al. Recombinant plant-derived human IgE glycoproteomics. J. Proteom. 2017, 161, 81–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Göritzer, K.; Turupcu, A.; Maresch, D.; Novak, J.; Altmann, F.; Oostenbrink, C.; Obinger, C.; Strasser, R. Distinct Fc alpha receptor N-glycans modulate the binding affinity to immunoglobulin A (IgA) antibodies. J. Biol. Chem. 2019, 294, 13995–14008. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Pol-Fachin, L.; Rusu, V.H.; Verli, H.; Lins, R.D. GROMOS 53A6 GLYC, an Improved GROMOS Force Field for Hexopyranose-Based Carbohydrates. J. Chem. Theory Comput. 2012, 8, 4681–4690. [Google Scholar] [CrossRef] [PubMed]
  16. Gebhardt, J.; Kleist, C.; Jakobtorweihen, S.; Hansen, N. Validation and Comparison of Force Fields for Native Cyclodextrins in Aqueous Solution. J. Phys. Chem. B 2018, 122, 1608–1626. [Google Scholar] [CrossRef] [PubMed]
  17. Turupcu, A.; Diem, M.; Smith, L.J.; Oostenbrink, C. Structural Aspects of the O-glycosylation Linkage in Glycopeptides via MD Simulations and Comparison with NMR Experiments. ChemPhysChem 2019, 20, 1527–1537. [Google Scholar] [CrossRef]
  18. Turupcu, A.; Oostenbrink, C. Modeling of Oligosaccharides within Glycoproteins from Free-Energy Landscapes. J. Chem. Inf. Model. 2017, 57, 2222–2236. [Google Scholar] [CrossRef] [Green Version]
  19. Baumann, H.; Brisson, J.R.; Michon, F.; Pon, R.; Jennings, H.J. Comparison of the conformation of the epitope of α(2→8) polysialic acid with its reduced and N-acyl derivatives. Biochemistry 1993, 32, 4007–4013. [Google Scholar] [CrossRef]
  20. Castillo, A.M.; Patiny, L.; Wist, J. Fast and accurate algorithm for the simulation of NMR spectra of large spin systems. J. Magn. Res. 2011, 209, 123–130. [Google Scholar] [CrossRef]
  21. Azurmendi, H.F.; Battistel, M.D.; Zarb, J.; Lichaa, F.; Virgen, A.N.; Shiloach, J.; Freedberg, D.I. The β-reducing end in α(2–8)-polysialic acid constitutes a unique structural motif. Glycobiology 2017, 27, 909–911. [Google Scholar] [CrossRef] [Green Version]
  22. Klepach, T.; Zhang, W.; Carmichael, I.; Serianni, A.S. 13C-1H and 13C-13C NMR J-Couplings in 13C-Labeled N-Acetyl-neuraminic Acid: Correlations with Molecular Structure. J. Org. Chem. 2008, 73, 4376–4387. [Google Scholar] [CrossRef] [PubMed]
  23. Spiwok, V.; Tvaroska, I. Conformational Free Energy Surface of α-N-Acetylneuraminic Acid: An Interplay Between Hydrogen Bonding and Solvation. J. Phys. Chem. B 2009, 113, 9589–9594. [Google Scholar] [CrossRef] [PubMed]
  24. Daura, X.; Gunsteren, W.F.V.; Mark, A.E. Folding-Unfolding Thermodynamics of a β-Heptapeptide From Equilibrium Simulations. Proteins 1999, 280, 269–280. [Google Scholar] [CrossRef]
  25. Hanashima, S.; Sato, C.; Tanaka, H.; Takahashi, T.; Ken, K.; Yamaguchi, Y. NMR study into the mechanism of recognition of the degree of polymerization by oligo/polysialic acid antibodies. Bioorg. Med. Chem. 2013, 21, 6069–6076. [Google Scholar] [CrossRef] [PubMed]
  26. Ray, G.J.; Ravenscroft, N.; Siekmann, J.J.; Zhang, Z.; Sanders, P.; Shaligram, U.; Szabo, C.M.; Kosma, P. Complete Structural Elucidation of an Oxidized Polysialic Acid Drug Intermediate by Nuclear Magnetic Resonance Spectroscopy. Bioconjug. Chem. 2014, 25, 665–676. [Google Scholar] [CrossRef]
  27. Schmid, N.; Christ, C.D.; Christen, M.; Eichenberger, A.P.; van Gunsteren, W.F. Architecture, implementation and parallelisation of the GROMOS software for biomolecular simulation. Comput. Phys. Commun. 2012, 183, 890–903. [Google Scholar] [CrossRef]
  28. Molecular Operating Environment (MOE); Chemical Computing Group: Montreal, QC, Canada, 2013.
  29. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; Hermans, J. Interaction models for water in relation to protein hydration. In Intermolecular Forces; Pullman, B., Ed.; Reidel: Dordrecht, The Netherlands, 1981; pp. 331–342. [Google Scholar]
  30. Berendsen, H.J.; Postma, J.P.; Van Gunsteren, W.F.; Dinola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  31. Hockney, R.W. The potential calculation and some applications. Methods Comput. Phys. 1970, 9, 136–211. [Google Scholar]
  32. Ryckaert, J.P.; Ciccotti, G.; Berendsen, H.J. Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 1977, 23, 327–341. [Google Scholar] [CrossRef] [Green Version]
  33. Tironi, I.G.; Sperb, R.; Smith, P.E.; van Gunsteren, W.F. A generalized reaction field method for molecular dynamics simulations. J. Chem. Phys. 1995, 102, 5451–5459. [Google Scholar] [CrossRef]
  34. Heinz, T.N.; van Gunsteren, W.F.; Hünenberger, P.H. Comparison of four methods to compute the dielectric permittivity of liquids from molecular dynamics simulations. J. Chem. Phys. 2001, 115, 1125–1136. [Google Scholar] [CrossRef]
  35. Heinz, T.N.; Hunenberger, P.H. A fast pairlist-construction algorithm for molecular simulations under periodic boundary conditions. J. Comput. Chem. 2004, 25, 1474–1486. [Google Scholar] [CrossRef] [PubMed]
  36. Eichenberger, A.P.; Allison, J.R.; Dolenc, J.; Geerke, D.P.; Horta, B.A.; Meier, K.; Oostenbrink, C.; Schmid, N.; Steiner, D.; Wang, D.; et al. GROMOS++ software for the analysis of biomolecular simulation trajectories. J. Chem. Theory Comput. 2011, 7, 3379–3390. [Google Scholar] [CrossRef] [PubMed]
  37. Hansen, H.S.; Hünenberger, P.H. Using the local elevation method to construct optimized umbrella sampling potentials: Calculation of the relative free energies and interconversion barriers of glucopyranose ring conformers in water. J. Comput. Chem. 2010, 31, 1–23. [Google Scholar] [CrossRef]
  38. Huber, T.; Torda, A.E.; van Gunsteren, W.F. Local elevation: A method for improving the searching properties of molecular dynamics simulation. J. Comput. Aided Mol. Des. 1994, 8, 695–708. [Google Scholar] [CrossRef]
  39. Van Gunsteren, W.F.; Billeter, S.; Eising, A.; Hünenberger, P.; Krüger, P.; Mark, A.; Scott, W.; Tironi, I. Biomolecular Simulation: The GROMOS96 Manual and User Guide; Vdf Hochschulverlag AG an der ETH Zürich: Groningen, Zürich, 1996. [Google Scholar]
  40. van Gunsteren, W.F.; Allison, J.R.; Daura, X.; Dolenc, J.; Hansen, N.; Mark, A.E.; Oostenbrink, C.; Rusu, V.H.; Smith, L.J. Deriving Structural Information from Experimentally Measured Data on Biomolecules. Angew. Chem. Int. Ed. 2016, 55, 15990–16010. [Google Scholar] [CrossRef]
  41. Haasnoot, C.A.; de Leeuw, F.A.; Altona, C. The relationship between proton-proton NMR coupling constants and substituent electronegativities-I. An empirical generalization of the karplus equation. Tetrahedron 1980, 36, 2783–2792. [Google Scholar] [CrossRef]
  42. Stenutz, R.; Carmichael, I.; Widmalm, G.; Serianni, A.S. Hydroxymethyl group conformation in saccharides: Structural dependencies of 2JHH, 3JHH, and 1JCH spin-spin coupling constants. J. Org. Chem. 2002, 67, 949–958. [Google Scholar] [CrossRef]
  43. Schmid, N.; Allison, J.R.; Dolenc, J.; Eichenberger, A.P.; Kunz, A.P.E.; van Gunsteren, W.F. Biomolecular structure refinement using the GROMOS simulation software. J. Biomol. NMR 2011, 51, 265–281. [Google Scholar] [CrossRef]
  44. Kirkwood, J.G. Statistical Mechanics of Fluid Mixtures. J. Chem. Phys. 1935, 3, 300–313. [Google Scholar] [CrossRef]
  45. Wang, L.; Deng, Y.; Wu, Y.; Kim, B.; LeBard, D.N.; Wandschneider, D.; Beachy, M.; Friesner, R.A.; Abel, R. Accurate Modeling of Scaffold Hopping Transformations in Drug Discovery. J. Chem. Theory Comput. 2017, 13, 42–54. [Google Scholar] [CrossRef] [PubMed]
  46. Zwanzig, R.W. High-Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases. J. Chem. Phys. 1954, 22, 1420–1426. [Google Scholar] [CrossRef]
Figure 1. Studied systems. (A) Schematic representation of dimer systems R=H or C, with the anomericity indicated in the superscript (B); trimer with residue labels (C); tetramer with residue labels; (D) decamer with n = 10 . Definition of the dihedral angles are ϕ =O6-C2-O8’-C8’, ψ =C2-O8’-C8’-C7’, ω 7=O7-C7-C6-O6, ω 8=O8-C8-C7-O7, and ω 9=O9-C9-C8-O8.
Figure 1. Studied systems. (A) Schematic representation of dimer systems R=H or C, with the anomericity indicated in the superscript (B); trimer with residue labels (C); tetramer with residue labels; (D) decamer with n = 10 . Definition of the dihedral angles are ϕ =O6-C2-O8’-C8’, ψ =C2-O8’-C8’-C7’, ω 7=O7-C7-C6-O6, ω 8=O8-C8-C7-O7, and ω 9=O9-C9-C8-O8.
Ijms 21 00030 g001
Figure 2. Free-energy maps G ( ϕ , ψ ) and G ( ω 7 , ω 8 ) of the glycosidic dihedral angles from the local elevation of umbrella sampling (LEUS) simulations. Contour maps are drawn with 5 kJ/mol spacing starting from the global minimum energy which is set to 0 kJ/mol. The regions that were never visited are shown in dark blue and the corresponding unbiased free energies are represented in the color maps. (A) dimer1 α H; (B) dimer2 β H (C); dimer3 α C.
Figure 2. Free-energy maps G ( ϕ , ψ ) and G ( ω 7 , ω 8 ) of the glycosidic dihedral angles from the local elevation of umbrella sampling (LEUS) simulations. Contour maps are drawn with 5 kJ/mol spacing starting from the global minimum energy which is set to 0 kJ/mol. The regions that were never visited are shown in dark blue and the corresponding unbiased free energies are represented in the color maps. (A) dimer1 α H; (B) dimer2 β H (C); dimer3 α C.
Ijms 21 00030 g002
Figure 3. θ vs. calculated 3 JH7H8, 3 JH8H9R, and 3 JH8H9S couplings from LEUS simulations for dimer1 α H. The first column represents the first non-reducing residue b where the ω 7 and ω 8 are free and the second column is for residue a where they are part of the glycosidic linkage. Experimental and calculated 3 J values are represented with solid and dashed horizontal lines, respectively. The colors on this Karplus curve indicate the preferred sampling after unbiasing of the LEUS simulations. In the unbiasing procedure, LEUS occurrences (P) are binned with 6 grid spacing. Negative values of l n P are set to zero. The definition of the θ for each J value is given in Figure 1 and Equation (4).
Figure 3. θ vs. calculated 3 JH7H8, 3 JH8H9R, and 3 JH8H9S couplings from LEUS simulations for dimer1 α H. The first column represents the first non-reducing residue b where the ω 7 and ω 8 are free and the second column is for residue a where they are part of the glycosidic linkage. Experimental and calculated 3 J values are represented with solid and dashed horizontal lines, respectively. The colors on this Karplus curve indicate the preferred sampling after unbiasing of the LEUS simulations. In the unbiasing procedure, LEUS occurrences (P) are binned with 6 grid spacing. Negative values of l n P are set to zero. The definition of the θ for each J value is given in Figure 1 and Equation (4).
Ijms 21 00030 g003
Figure 4. θ vs. calculated 3 JH7H8, 3 JH8H9R, and 3 JH8H9S couplings from LEUS simulations for dimer2 β H. The first column represents the non-reducing residue b, where the ω 7 and ω 8 are free and the second column is residue a where they are part of the glycosidic linkage. Experimental and calculated 3 J values are represented with solid and dashed horizontal lines, respectively. The colors on this Karplus curve indicate the preferred sampling after unbiasing of the LEUS simulations. In the unbiasing procedure, LEUS occurrences (P) are binned with 6 grid spacing. Negative values of l n P set to zero. The definition of the θ for each J value is given in Figure 1.
Figure 4. θ vs. calculated 3 JH7H8, 3 JH8H9R, and 3 JH8H9S couplings from LEUS simulations for dimer2 β H. The first column represents the non-reducing residue b, where the ω 7 and ω 8 are free and the second column is residue a where they are part of the glycosidic linkage. Experimental and calculated 3 J values are represented with solid and dashed horizontal lines, respectively. The colors on this Karplus curve indicate the preferred sampling after unbiasing of the LEUS simulations. In the unbiasing procedure, LEUS occurrences (P) are binned with 6 grid spacing. Negative values of l n P set to zero. The definition of the θ for each J value is given in Figure 1.
Ijms 21 00030 g004
Figure 5. Central members of the most populated clusters of the jointly clustered configurations of α -2,8-decasialic acid from plain MD and H-bond restrained simulations. Clusters 1–3 and 5–9 contained conformations that originated from the restrained simulations. Cluster 4 originated from the plain MD simulations.
Figure 5. Central members of the most populated clusters of the jointly clustered configurations of α -2,8-decasialic acid from plain MD and H-bond restrained simulations. Clusters 1–3 and 5–9 contained conformations that originated from the restrained simulations. Cluster 4 originated from the plain MD simulations.
Ijms 21 00030 g005
Figure 6. Conformational sampling for dimer3 α C in a 4D local elevation buildup simulation after 200 ns (A) and 400 ns (B).
Figure 6. Conformational sampling for dimer3 α C in a 4D local elevation buildup simulation after 200 ns (A) and 400 ns (B).
Ijms 21 00030 g006
Table 1. Experimental and calculated 3 JHH coupling constants for α -Neu5Ac-(2-8)- α -Neu5Ac-OH (dimer1 α H) and α -Neu5Ac-(2-8)- β -Neu5Ac-OH (dimer2 β H). Calculated values from LEUS simulations after unbiasing.
Table 1. Experimental and calculated 3 JHH coupling constants for α -Neu5Ac-(2-8)- α -Neu5Ac-OH (dimer1 α H) and α -Neu5Ac-(2-8)- β -Neu5Ac-OH (dimer2 β H). Calculated values from LEUS simulations after unbiasing.
dimer1 α Hdimer2 β H
ResidueNMR aLEUSResidueNMR bLEUS
3 JH6H7 [Hz]b1.4 ± 0.11.4 ± 0.2b< 1.50.9 ± 0.1
a<1.00.9 ± 0.3a< 1.50.8 ± 0.1
3 JH7H8 [Hz]b9.5 ± 0.93.7 ± 0.6b9.6 ± 0.58.4 ± 0.1
a1.5 ± 0.21.8 ± 0.1a7.3 ± 0.57.4 ± 0.1
3 JH8H9R [Hz]b6.1 ± 0.66.1 ± 0.1b6.0 ± 0.56.6 ± 0.1
a6.1 ± 0.66.3 ± 0.1a4.0 ± 0.56.5 ± 0.1
3 JH8H9S [Hz]b2.4 ± 0.21.1 ± 0.1b2.6 ± 0.51.5 ± 0.1
a4.1 ± 0.42.2 ± 0.5a3.3 ± 0.51.4 ± 0.1
a Experimental values from Ref. [10]; b from Ref. [19].
Table 2. Experimental and calculated 3 JHH coupling constants for α -Neu5Ac-(2-8)- α -Neu5Ac- α -(2-8)- α - Neu5Ac-OMe (trimer1) and α -Neu5Ac-(2-8)- α -Neu5Ac-OMe (dimer 3 α C). Calculated values from LEUS simulations after reweighing.
Table 2. Experimental and calculated 3 JHH coupling constants for α -Neu5Ac-(2-8)- α -Neu5Ac- α -(2-8)- α - Neu5Ac-OMe (trimer1) and α -Neu5Ac-(2-8)- α -Neu5Ac-OMe (dimer 3 α C). Calculated values from LEUS simulations after reweighing.
trimer1 α Cdimer3 α C
ResidueNMR aLEUSResidueLEUS
3 JH6H7 [Hz]c1.5 ± 0.21.0 ± 0.2b1.0 ± 0.1
b< 1.01.0 ± 0.3a0.9 ± 0.1
a< 1.01.0 ± 0.2
3 JH7H8 [Hz]c9.6 ± 1.04.7 ± 0.5b2.5 ± 0.5
b< 4.03.3 ± 0.7a1.0 ± 0.5
a< 4.06.4 ± 0.7
3 JH8H9R [Hz]cn.d.7.1 ± 0.1b5.9 ± 0.5
bn.d.6.0 ± 0.1a5.9 ± 0.5
an.d.6.6 ± 0.1
3 JH8H9S [Hz]cn.d.2.8 ± 0.2b1.8 ± 0.5
bn.d.1.2 ± 0.1a1.4 ± 0.5
an.d.0.9 ± 0.1
a Experimental values from Ref. [10].
Table 3. Experimental and calculated 3 JHH coupling constants for α -Neu5Ac-(2-8)- α -Neu5Ac-(2-8)- α - Neu5Ac-(2-8)- β -Neu5Ac-OH (tetramer β H). Calculated values from LEUS simulations after reweighing.
Table 3. Experimental and calculated 3 JHH coupling constants for α -Neu5Ac-(2-8)- α -Neu5Ac-(2-8)- α - Neu5Ac-(2-8)- β -Neu5Ac-OH (tetramer β H). Calculated values from LEUS simulations after reweighing.
3 Jtetramer β H
ResidueNMR aLEUS
3 JH6H7d2.0 ± 0.10.9 ± 0.1
c< 1.00.8 ± 0.1
b< 1.00.6 ± 0.1
a< 1.00.6 ± 0.1
3 JH7H8d8.9 ± 0.19.8 ± 0.0
c4.5 ± 0.11.7 ± 0.1
b2.7 ± 0.11.5 ± 0.4
a6.3 ± 0.12.7 ± 0.6
3 JH8H9Rd6.1 ± 0.16.9 ± 0.1
c5.9 ± 0.10.3 ± 0.6
b6.1 ± 0.15.1 ± 0.1
a4.3 ± 0.16.1 ± 0.1
3 JH8H9Sd2.5 ± 0.10.7 ± 0.3
c4.1 ± 0.110.6 ± 0.0
b5.4 ± 0.11.6 ± 0.1
a2.8 ± 0.11.2 ± 0.4
a NMR experiments were conducted at 600 MHz and 800 MHz.
Table 4. Intra- and Inter-residue hydrogen bond occurrences along with the water bridges calculated from LEUS simulations for dimers and trimer. Reweighted populations higher than 2% are reported. ’-’ indicates that the hydrogen bond is not observed. ’n.a.’ indicates that the hydrogen bond not applicable for this system.
Table 4. Intra- and Inter-residue hydrogen bond occurrences along with the water bridges calculated from LEUS simulations for dimers and trimer. Reweighted populations higher than 2% are reported. ’-’ indicates that the hydrogen bond is not observed. ’n.a.’ indicates that the hydrogen bond not applicable for this system.
H-bondSystem
Typedimer1 α Hdimer2 β Hdimer3 α Ctrimer1 α C
Intra-residuebHO8-bO1A/B-5.9%5.3%-
aHO2-aO1A/B-25.2%--
Inter-residueaHO9-bO1A/B38.4%3.2%19.4%1.7%
aHO7-bO1A/B-46.7 %8.5%14.9%
aHO7-bO64.3%--14.7%
aHO9-bO6-16.7%8.2%-
aHN5-bO1A/B--3.1%-
aHO4-bO1A/B-7.7%--
bHO9-aO1A/B-40.0%--
bHO7-aO1A/B-2.4%--
cHO9-bO1A/Bn.an.a.n.a.3.8%
cHO7-bO6n.a.n.a.n.a.6.3%
Water-bridgeaO6-bO8-11.5%--
bO6-bO2-56.5%14.9%-
bHN5-bO2-17.9%--

Share and Cite

MDPI and ACS Style

Turupcu, A.; Blaukopf, M.; Kosma, P.; Oostenbrink, C. Molecular Conformations of Di-, Tri-, and Tetra-α-(2→8)-Linked Sialic Acid from NMR Spectroscopy and MD Simulations. Int. J. Mol. Sci. 2020, 21, 30. https://doi.org/10.3390/ijms21010030

AMA Style

Turupcu A, Blaukopf M, Kosma P, Oostenbrink C. Molecular Conformations of Di-, Tri-, and Tetra-α-(2→8)-Linked Sialic Acid from NMR Spectroscopy and MD Simulations. International Journal of Molecular Sciences. 2020; 21(1):30. https://doi.org/10.3390/ijms21010030

Chicago/Turabian Style

Turupcu, Aysegül, Markus Blaukopf, Paul Kosma, and Chris Oostenbrink. 2020. "Molecular Conformations of Di-, Tri-, and Tetra-α-(2→8)-Linked Sialic Acid from NMR Spectroscopy and MD Simulations" International Journal of Molecular Sciences 21, no. 1: 30. https://doi.org/10.3390/ijms21010030

APA Style

Turupcu, A., Blaukopf, M., Kosma, P., & Oostenbrink, C. (2020). Molecular Conformations of Di-, Tri-, and Tetra-α-(2→8)-Linked Sialic Acid from NMR Spectroscopy and MD Simulations. International Journal of Molecular Sciences, 21(1), 30. https://doi.org/10.3390/ijms21010030

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop