Glycosidic α-linked mannopyranose disaccharides: an NMR spectroscopy and molecular dynamics simulation study employing additive and Drude polarizable force fields.

D-Mannose is a structural component in N-linked glycoproteins from viruses and mammals as well as in polysaccharides from fungi and bacteria. Structural components often consist of D-Manp residues joined via α-(1→2)-, α-(1→3)-, α-(1→4)- or α-(1→6)-linkages. As models for these oligo- and polysaccharides, a series of mannose-containing disaccharides have been investigated with respect to conformation and dynamics. Translational diffusion NMR experiments were performed to deduce rotational correlation times for the molecules, 1D 1H,1H-NOESY and 1D 1H,1H-T-ROESY NMR experiments were carried out to obtain inter-residue proton-proton distances and one-dimensional long-range and 2D J-HMBC experiments were acquired to gain information about conformationally dependent heteronuclear coupling constants across glycosidic linkages. To attain further spectroscopic data, the doubly 13C-isotope labeled α-D-[1,2-13C2]Manp-(1→4)-α-D-Manp-OMe was synthesized thereby facilitating conformational analysis based on 13C,13C coupling constants as interpreted by Karplus-type relationships. Molecular dynamics simulations were carried out for the disaccharides with explicit water as solvent using the additive CHARMM36 and Drude polarizable force fields for carbohydrates, where the latter showed broader population distributions. Both simulations sampled conformational space in such a way that inter-glycosidic proton-proton distances were very well described whereas in some cases deviations were observed between calculated conformationally dependent NMR scalar coupling constants and those determined from experiment, with closely similar root-mean-square differences for the two force fields. However, analyses of dipole moments and radial distribution functions with water of the hydroxyl groups indicate differences in the underlying physical forces dictating the wider conformational sampling with the Drude polarizable versus additive C36 force field and indicate the improved utility of the Drude polarizable model in investigating complex carbohydrates.

Extension of the CHARMM Classical Drude Polarizable Force Field to N- and O-Linked Glycopeptides and Glycoproteins.

Molecular dynamic simulations are an effective tool to study complex molecular systems and are contingent upon the availability of an accurate and reliable molecular mechanics force field. The Drude polarizable force field, which allows for the explicit treatment of electronic polarization in a computationally efficient fashion, has been shown to reproduce experimental properties that were difficult or impossible to reproduce with the CHARMM additive force field, including peptide folding cooperativity, RNA hairpin structures, and DNA base flipping. Glycoproteins are essential components of glycoconjugate vaccines, antibodies, and many pharmaceutically important molecules, and an accurate polarizable force field that includes compatibility between the protein and carbohydrate aspect of the force field is essential to study these types of systems. In this work, we present an extension of the Drude polarizable force field to glycoproteins, including both N- and O-linked species. Parameter optimization focused on the dihedral terms using a reweighting protocol targeting NMR solution J-coupling data for model glycopeptides. Validation of the model include eight model glycopeptides and four glycoproteins with multiple N- and O-linked glycosylations. The new glycoprotein carbohydrate force field can be used in conjunction with the remainder of Drude polarizable force field through a variety of MD simulation programs including GROMACS, OPENMM, NAMD, and CHARMM and may be accessed through the Drude Prepper module in the CHARMM-GUI.

Development of CHARMM Additive Potential Energy Parameters for α-Methyl Amino Acids.

Potential energy parameters for α-methyl amino acids were generated with ab initio calculations on α-methyl-N-acetylalanyl-N'-methylamide (the α-methyl "alanine dipeptide") which served as an input to a grid-based correction to the backbone torsional potential (known as CMAP) consistent with the CHARMM36m additive protein force field. The new parameters were validated by comparison with experimentally determined helicities of the 22 residue C-terminal peptide (H10) from apolipoprotein A1 and five α-methylated variants in water and 0.3:0.7 trifluoroethanol (TFE)/water. Conventional molecular dynamics simulation totaling 30 µs for each peptide is in overall good agreement with the experiment, including the increased helicity in 30% TFE. An additional 500 ns of simulation using two-dimensional dihedral biasing (bpCMAP) replica exchange reduced left-handed conformations, increased right-handed helices, and thereby mostly decreased agreement with the experiment. Analysis of side chain-side chain salt bridges suggests that the overestimation of the helical content may be, in part, due to such interactions. The increased helicity of the peptides in 30% TFE arises from decreased hydrogen bonding of the backbone atoms to water and a concomitant increase in intramolecular backbone hydrogen bonds.

Insights into substrate recognition and specificity for IgG by Endoglycosidase S2.

Antibodies bind foreign antigens with high affinity and specificity leading to their neutralization and/or clearance by the immune system. The conserved N-glycan on IgG has significant impact on antibody effector function, with the endoglycosidases of Streptococcus pyogenes deglycosylating the IgG to evade the immune system, a process catalyzed by the endoglycosidase EndoS2. Studies have shown that two of the four domains of EndoS2, the carbohydrate binding module (CBM) and the glycoside hydrolase (GH) domain are critical for catalytic activity. To yield structural insights into contributions of the CBM and the GH domains as well as the overall flexibility of EndoS2 to the proteins' catalytic activity, models of EndoS2-Fc complexes were generated through enhanced-sampling molecular-dynamics (MD) simulations and site-identification by ligand competitive saturation (SILCS) docking followed by reconstruction and multi-microsecond MD simulations. Modeling results predict that EndoS2 initially interacts with the IgG through its CBM followed by interactions with the GH yielding catalytically competent states. These may involve the CBM and GH of EndoS2 simultaneously interacting with either the same Fc CH2/CH3 domain or individually with the two Fc CH2/CH3 domains, with EndoS2 predicted to assume closed conformations in the former case and open conformations in the latter. Apo EndoS2 is predicted to sample both the open and closed states, suggesting that either complex can directly form following initial IgG-EndoS2 encounter. Interactions of the CBM and GH domains with the IgG are predicted to occur through both its glycan and protein regions. Simulations also predict that the Fc glycan can directly transfer from the CBM to the GH, facilitating formation of catalytically competent complexes and how the 734 to 751 loop on the CBM can facilitate extraction of the glycan away from the Fc CH2/CH3 domain. The predicted models are compared and consistent with Hydrogen/Deuterium Exchange data. In addition, the complex models are consistent with the high specificity of EndoS2 for the glycans on IgG supporting the validity of the predicted models.

Site-Selective Chemoenzymatic Modification on the Core Fucose of an Antibody Enhances Its Fcγ Receptor Affinity and ADCC Activity.

Fc glycosylation profoundly impacts the effector functions of antibodies and often dictates an antibody's pro- or anti-inflammatory activities. It is well established that core fucosylation of the Fc domain N-glycans of an antibody significantly reduces its affinity for FcγRIIIa receptors and antibody-dependent cellular cytotoxicity (ADCC). Previous structural studies have suggested that the presence of a core fucose remarkably decreases the unique and favorable carbohydrate-carbohydrate interactions between the Fc and the receptor N-glycans, leading to reduced affinity. We report here that in contrast to natural core fucose, special site-specific modification on the core fucose could dramatically enhance the affinity of an antibody for FcγRIIIa. The site-selective modification was achieved through an enzymatic transfucosylation with a novel fucosidase mutant, which was shown to be able to use modified α-fucosyl fluoride as the donor substrate. We found that replacement of the core l-fucose with 6-azide- or 6-hydroxy-l-fucose (l-galactose) significantly enhanced the antibody's affinity for FcγRIIIa receptors and substantially increased the ADCC activity. To understand the mechanism of the modified fucose-mediated affinity enhancement, we performed molecular dynamics simulations. Our data revealed that the number of glycan contacts between the Fc and the Fc receptor was increased by the selective core-fucose modifications, showing the importance of unique carbohydrate-carbohydrate interactions in achieving high FcγRIIIa affinity and ADCC activity of antibodies. Thus, the direct site-selective modification turns the adverse effect of the core fucose into a favorable force to promote the carbohydrate-carbohydrate interactions.

Balanced polarizable Drude force field parameters for molecular anions: phosphates, sulfates, sulfamates, and oxides.

Polarizable force fields are emerging as a more accurate alternative to additive force fields in terms of modeling and simulations of a variety of chemicals including biomolecules. Explicit treatment of induced polarization in charged species such as phosphates and sulfates offers the potential for achieving an improved atomistic understanding of the physical forces driving their interactions with their environments. To help achieve this, in this study we present balanced Drude polarizable force field parameters for molecular ions including phosphates, sulfates, sulfamates, and oxides. Better balance was primarily achieved in the relative values of minimum interaction energies and distances of the anionic model compounds with water at the Drude and quantum mechanical (QM) model chemistries. Parametrization involved reoptimizing available parameters as well as extending the force field to new molecules with the goal of achieving self-consistency with respect to the Lennard-Jones and electrostatic parameters targeting QM and experimental hydration free energies. The resulting force field parameters achieve consistent treatment across the studied anions, facilitating more balanced simulations of biomolecules and small organic molecules in the context of the classical Drude polarizable force field. Graphical abstract.

Impact of electronic polarizability on protein-functional group interactions.

Interactions of proteins with functional groups are key to their biological functions, making it essential that they be accurately modeled. To investigate the impact of the inclusion of explicit treatment of electronic polarizability in force fields on protein-functional group interactions, the additive CHARMM and Drude polarizable force field are compared in the context of the Site-Identification by Ligand Competitive Saturation (SILCS) simulation methodology from which functional group interaction patterns with five proteins for which experimental binding affinities of multiple ligands are available, were obtained. The explicit treatment of polarizability produces significant differences in the functional group interactions in the ligand binding sites including overall enhanced binding of functional groups to the proteins. This is associated with variations of the dipole moments of solutes representative of functional groups in the binding sites relative to aqueous solution with higher dipole moments systematically occurring in the latter, though exceptions occur with positively charged methylammonium. Such variation indicates the complex, heterogeneous nature of the electronic environments of ligand binding sites and emphasizes the inherent limitation of fixed charged, additive force fields for modeling ligand-protein interactions. These effects yield more defined orientation of the functional groups in the binding pockets and a small, but systematic improvement in the ability of the SILCS method to predict the binding orientation and relative affinities of ligands to their target proteins. Overall, these results indicate that the physical model associated with the explicit treatment of polarizability along with the presence of lone pairs in a force field leads to changes in the nature of the interactions of functional groups with proteins versus that occurring with additive force fields, suggesting the utility of polarizable force fields in obtaining a more realistic understanding of protein-ligand interactions.

Drude Polarizable Force Field Parametrization of Carboxylate and N-Acetyl Amine Carbohydrate Derivatives.

In this work, we report the development of Drude polarizable force field parameters for the carboxylate and N-acetyl amine derivatives, extending the functionality of the existing Drude polarizable carbohydrate force field. The force field parameters have been developed in a hierarchical manner, reproducing the quantum mechanical gas-phase properties of small model compounds representing the key functional group in the carbohydrate derivatives, including optimization of the electrostatic and bonded parameters. The optimized parameters were then used to generate the models for carboxylate and N-acetyl amine carbohydrate derivatives. The transferred parameters were further tested and optimized to reproduce crystal geometries and J-coupling data from nuclear magnetic resonance experiments. The parameter development resulted in the incorporation of d-glucuronate, l-iduronate, N-acetyl-d-glucosamine (GlcNAc), and N-acetyl-d-galactosamine (GalNAc) sugars into the Drude polarizable force field. The parameters developed in this study were then applied to study the conformational properties of glycosaminoglycan polymer hyaluronan, composed of d-glucuronate and N-acetyl-d-glucosamine, in aqueous solution. Upon comparing the results from the additive and polarizable simulations, it was found that the inclusion of polarization improved the description of the electrostatic interactions observed in hyaluronan, resulting in enhanced conformational flexibility. The developed Drude polarizable force field parameters in conjunction with the remainder of the Drude polarizable force field parameters can be used for future studies involving carbohydrates and their conjugates in complex, heterogeneous systems.

Impact of branching on the conformational heterogeneity of the lipopolysaccharide from Klebsiella pneumoniae: Implications for vaccine design.

Resistance of Klebsiella pneumoniae (KP) to antibiotics has motivated the development of an efficacious KP human vaccine that would not be subject to antibiotic resistance. Klebsiella lipopolysaccharide (LPS) associated O polysaccharide (OPS) types have provoked broad interest as a vaccine antigen as there are only 4 that predominate worldwide (O1, O2a, O3, O5). Klebsiella O1 and O2 OPS are polygalactans that share a common D-Gal-I structure, for which a variant D-Gal-III was recently discovered. To understand the potential impact of this variability on antigenicity, a detailed molecular picture of the conformational differences associated with the addition of the D-Gal-III (1 → 4)-α-Galp branch is presented using enhanced-sampling molecular dynamics simulations. In D-Gal-I two major conformational states are observed while the presence of the 1 → 4 branch in D-Gal-III resulted in only a single dominant extended state. Stabilization of the more folded states in D-Gal-I is due to a O4-H⋯O2 hydrogen bond in the linear backbone that cannot occur in D-Gal-III as the O4 is in the Galp(1 → 4)Galp glycosidic linkage. The impact of branching in D-Gal-III also significantly decreases the accessibility of the monosaccharides in the linear backbone region of D-Gal-I, while the accessibility of the terminal D-Gal-II region of the OPS is not substantially altered. The present results suggest that a vaccine that targets both the D-Gal-I and D-Gal-III LPS can be developed by using D-Gal-III as the antigen combined with cross-reactivity experiments using the Gal-II polysaccharide to assure that this region of the LPS is the primary epitope of the antigen.

CHARMM Drude Polarizable Force Field for Glycosidic Linkages Involving Pyranoses and Furanoses.

We present an extension of the CHARMM Drude polarizable force field to enable modeling of polysaccharides containing pyranose and furanose monosaccharides. The new force field parameters encompass 1↔2, 1→3, 1→4, and 1→6 pyranose-furanose linkages, 2→1 and 2→6 furanose-furanose linkages, 2→2, 2→3, and 2→4 furanose-pyranose, and 1↔1, 1→2, 1→3, 1→4, and 1→6 pyranose-pyranose linkages. For the glycosidic linkages, both simple model compounds and the full disaccharides with methylation at the reducing end were used for parameter optimization. The model compounds were chosen to be monomers or glycosidic-linked dimers of tetrahydropyran (THP) and tetrahydrofuran (THF). Target data for optimization included one- and two-dimensional potential energy scans of ω and the Φ/Ψ glycosidic dihedral angles in the model compounds and full disaccharides computed by quantum mechanical (QM) RIMP2/cc-pVQZ single point energies on MP2/6-31G(d) optimized structures. Also included in the target data are extensive sets of QM gas phase monohydrate water-saccharide interactions, dipole moments, and molecular polarizabilities for both model compounds and full disaccharides. The resulting polarizable model is shown to be in good agreement with a range of QM data, offering a significant improvement over the additive CHARMM36 carbohydrate force field, as well as experimental data including crystal structures and conformational properties of disaccharides and a trisaccharide in aqueous solution.

Proper balance of solvent-solute and solute-solute interactions in the treatment of the diffusion of glucose using the Drude polarizable force field.

Motivated by underestimation of the diffusion constant of glucose in aqueous solution at high glucose concentrations we performed additional optimization of the Drude polarizable hexopyranose monosaccharide force field. This indicated aggregation of the glucose at higher concentrations, which is a concern for studies of complex glycan systems such as the HIV Envelope where high effective concentrations of sugars are present. High-level quantum mechanical calculations were undertaken on water monohydrate-glucose interactions, on water cluster-glucose interactions and on glucose-glucose dimers in stacked (parallel) and perpendicular orientations. Optimization of the nonbond and dihedral parameters targeting these data yielded a revised model that showed improved agreement with experimental aqueous diffusion data. However, limitations in the diffusion constants were still present. These were due to the SWM4-NDP inherently overestimating the diffusion constant of water, a problem that was validated by calculation of the aqueous diffusion constants using the SWM6-NDP water model. In addition, results show the water diffusion constant to be significantly overestimated at high glucose concentrations though the glucose diffusion is in satisfactory agreement with experiment. These results indicate the subtle balance of water-sugar, water-water and sugar-sugar interactions that needs to be properly modeled to account for the full range of aqueous behavior of sugars in aqueous solution.

Drude polarizable force field for aliphatic ketones and aldehydes, and their associated acyclic carbohydrates.

The majority of computer simulations exploring biomolecular function employ Class I additive force fields (FF), which do not treat polarization explicitly. Accordingly, much effort has been made into developing models that go beyond the additive approximation. Development and optimization of the Drude polarizable FF has yielded parameters for selected lipids, proteins, DNA and a limited number of carbohydrates. The work presented here details parametrization of aliphatic aldehydes and ketones (viz. acetaldehyde, propionaldehyde, butaryaldehyde, isobutaryaldehyde, acetone, and butanone) as well as their associated acyclic sugars (D-allose and D-psicose). LJ parameters are optimized targeting experimental heats of vaporization and molecular volumes, while the electrostatic parameters are optimized targeting QM water interactions, dipole moments, and molecular polarizabilities. Bonded parameters are targeted to both QM and crystal survey values, with the models for ketones and aldehydes shown to be in good agreement with QM and experimental target data. The reported heats of vaporization and molecular volumes represent a compromise between the studied model compounds. Simulations of the model compounds show an increase in the magnitude and the fluctuations of the dipole moments in moving from gas phase to condensed phases, which is a phenomenon that the additive FF is intrinsically unable to reproduce. The result is a polarizable model for aliphatic ketones and aldehydes including the acyclic sugars D-allose and D-psicose, thereby extending the available biomolecules in the Drude polarizable FF.

Revised RNA Dihedral Parameters for the Amber Force Field Improve RNA Molecular Dynamics.

The backbone dihedral parameters of the Amber RNA force field were improved by fitting using multiple linear regression to potential energies determined by quantum chemistry calculations. Five backbone and four glycosidic dihedral parameters were fit simultaneously to reproduce the potential energies determined by a high-level density functional theory calculation (B97D3 functional with the AUG-CC-PVTZ basis set). Umbrella sampling was used to determine conformational free energies along the dihedral angles, and these better agree with the population of conformations observed in the protein data bank for the new parameters than for the conventional parameters. Molecular dynamics simulations performed on a set of hairpin loops, duplexes and tetramers with the new parameter set show improved modeling for the structures of tetramers CCCC, CAAU, and GACC, and an RNA internal loop of noncanonical pairs, as compared to the conventional parameters. For the tetramers, the new parameters largely avoid the incorrect intercalated structures that dominate the conformational samples from the conventional parameters. For the internal loop, the major conformation solved by NMR is stable with the new parameters, but not with the conventional parameters. The new force field performs similarly to the conventional parameters for the UUCG and GCAA hairpin loops and the [U(UA)6A]2 duplex.

Molecular mechanism for preQ1-II riboswitch function revealed by molecular dynamics.

Riboswitches are RNA molecules that regulate gene expression using conformational change, affected by binding of small molecule ligands. A crystal structure of a ligand-bound class II preQ1 riboswitch has been determined in a previous structural study. To gain insight into the dynamics of this riboswitch in solution, eight total molecular dynamic simulations, four with and four without ligand, were performed using the Amber force field. In the presence of ligand, all four of the simulations demonstrated rearranged base pairs at the 3' end, consistent with expected base-pairing from comparative sequence analysis in a prior bioinformatic analysis; this suggests the pairing in this region was altered by crystallization. Additionally, in the absence of ligand, three of the simulations demonstrated similar changes in base-pairing at the ligand binding site. Significantly, although most of the riboswitch architecture remained intact in the respective trajectories, the P3 stem was destabilized in the ligand-free simulations in a way that exposed the Shine-Dalgarno sequence. This work illustrates how destabilization of two major groove base triples can influence a nearby H-type pseudoknot and provides a mechanism for control of gene expression by a fold that is frequently found in bacterial riboswitches.

Structural analysis of a class III preQ1 riboswitch reveals an aptamer distant from a ribosome-binding site regulated by fast dynamics.

PreQ1-III riboswitches are newly identified RNA elements that control bacterial genes in response to preQ1 (7-aminomethyl-7-deazaguanine), a precursor to the essential hypermodified tRNA base queuosine. Although numerous riboswitches fold as H-type or HLout-type pseudoknots that integrate ligand-binding and regulatory sequences within a single folded domain, the preQ1-III riboswitch aptamer forms a HLout-type pseudoknot that does not appear to incorporate its ribosome-binding site (RBS). To understand how this unusual organization confers function, we determined the crystal structure of the class III preQ1 riboswitch from Faecalibacterium prausnitzii at 2.75 Å resolution. PreQ1 binds tightly (KD,app 6.5 ± 0.5 nM) between helices P1 and P2 of a three-way helical junction wherein the third helix, P4, projects orthogonally from the ligand-binding pocket, exposing its stem-loop to base pair with the 3' RBS. Biochemical analysis, computational modeling, and single-molecule FRET imaging demonstrated that preQ1 enhances P4 reorientation toward P1-P2, promoting a partially nested, H-type pseudoknot in which the RBS undergoes rapid docking (kdock ∼ 0.6 s(-1)) and undocking (kundock ∼ 1.1 s(-1)). Discovery of such dynamic conformational switching provides insight into how a riboswitch with bipartite architecture uses dynamics to modulate expression platform accessibility, thus expanding the known repertoire of gene control strategies used by regulatory RNAs.

Modified Amber Force Field Correctly Models the Conformational Preference for Tandem GA pairs in RNA.

Molecular mechanics with all-atom models was used to understand the conformational preference of tandem guanine-adenine (GA) noncanonical pairs in RNA. These tandem GA pairs play important roles in determining stability, flexibility, and structural dynamics of RNA tertiary structures. Previous solution structures showed that these tandem GA pairs adopt either imino (cis Watson-Crick/Watson-Crick A-G) or sheared (trans Hoogsteen/sugar edge A-G) conformations depending on the sequence and orientation of the adjacent closing base pairs. The solution structures (GCGGACGC)2 [Biochemistry, 1996, 35, 9677-9689] and (GCGGAUGC)2 [Biochemistry, 2007, 46, 1511-1522] demonstrate imino and sheared conformations for the two central GA pairs, respectively. These systems were studied using molecular dynamics and free energy change calculations for conformational changes, using umbrella sampling. For the structures to maintain their native conformations during molecular dynamics simulations, a modification to the standard Amber ff10 force field was required, which allowed the amino group of guanine to leave the plane of the base [J. Chem. Theory Comput., 2009, 5, 2088-2100] and form out-of-plane hydrogen bonds with a cross-strand cytosine or uracil. The requirement for this modification suggests the importance of out-of-plane hydrogen bonds in stabilizing the native structures. Free energy change calculations for each sequence demonstrated the correct conformational preference when the force field modification was used, but the extent of the preference is underestimated.