Multiple Conformational States Contribute to the 3D Structure of a Glucan Decasaccharide: A Combined SAXS and MD Simulation Study

Conformational analysis and interaction of the Staphylococcus aureus transmembrane peptidase AgrB with its AgrD propeptide substrate

Virulence gene expression in the human pathogen, S. aureus is regulated by the agr (accessory gene regulator) quorum sensing (QS) system which is conserved in diverse Gram-positive bacteria. The agr QS signal molecule is an autoinducing peptide (AIP) generated via the initial processing of the AgrD pro-peptide by the transmembrane peptidase AgrB. Since structural information for AgrB and AgrBD interactions are lacking, we used homology modelling and molecular dynamics (MD) annealing to characterise the conformations of AgrB and AgrD in model membranes and in solution. These revealed a six helical transmembrane domain (6TMD) topology for AgrB. In solution, AgrD behaves as a disordered peptide, which binds N-terminally to membranes in the absence and in the presence of AgrB. In silico, membrane complexes of AgrD and dimeric AgrB show non-equivalent AgrB monomers responsible for initial binding and for processing, respectively. By exploiting split luciferase assays in Staphylococcus aureus, we provide experimental evidence that AgrB interacts directly with itself and with AgrD. We confirmed the in vitro formation of an AgrBD complex and AIP production after Western blotting using either membranes from Escherichia coli expressing AgrB or with purified AgrB and T7-tagged AgrD. AgrB and AgrD formed stable complexes in detergent micelles revealed using synchrotron radiation CD (SRCD) and Landau analysis consistent with the enhanced thermal stability of AgrB in the presence of AgrD. Conformational alteration of AgrB following provision of AgrD was observed by small angle X-ray scattering from proteodetergent micelles. An atomistic description of AgrB and AgrD has been obtained together with confirmation of the AgrB 6TMD membrane topology and existence of AgrBD molecular complexes in vitro and in vivo.

Assessing the sequence dependence of pyrimidine-pyrimidone (6-4) photoproduct in a duplex double-stranded DNA: A pitfall for microsecond range simulation

Sequence dependence of the (6-4) photoproduct conformational landscape when embedded in six 25-bp duplexes is evaluated along extensive unbiased and enhanced (replica exchange with solute tempering, REST2) molecular dynamics simulations. The structural reorganization as the central pyrimidines become covalently tethered is traced back in terms of non-covalent interactions, DNA bending, and extrusion of adenines of the opposite strands. The close sequence pattern impacts the conformational landscape around the lesion, inducing different upstream and downstream flexibilities. Moreover, REST2 simulations allow us to probe structures possibly important for damaged DNA recognition.

Relation between glycosidic linkage, structure and dynamics of α- and β-glucans in water

In a molecular dynamics simulation study of several oligosaccharides comprising of the very basic building block of carbohydrate, the α- or β-d glucopyranose units, linked by any one of the 1-3/1-4 or 1-6 glycosidic linkages, we compare and contrast their structural and dynamical properties. Results indicate that the litheness of the oligosaccharide chain is noticeably controlled by the composition, anomeric nature and glycosidic linkage type of the units. In mixed β 1-4/1-3 d-glucopyranosides, as those found in oats and barley, the ratio of the β 1-4 and β 1-3 linked residues is crucial in determining the structural and dynamical attributes. Principal component analysis (PCA) using the internal coordinates of torsion angles subtended by glycosidic oxygen atoms and subsequent K-means clustering of the dynamical space spanned by PC1 to PC2 point to the dynamical and structural disparity in the various types of oligosaccharides studied. The properties simulated in this work are meant to provide a systematic yet comparative understanding of the importance of linkage and anomericity on the oligosaccharide chain properties and are in line with some experimental structural attributes.

Modeling the conformations of Neisseria meningitidis serogroup a CPS and a carba-analogue: Implications for vaccine development

Neisseria meningitidis is a major cause of bacterial meningitis worldwide, especially in Africa. The capsular polysaccharide is the main virulence factor and the target antigen for polysaccharide- and conjugate vaccines. Three tetravalent conjugate vaccines against serogroups A, C, Y and W have been licensed and the monovalent MenAfriVac® was introduced to address the high burden of serogroup A disease in the Meningitis Belt of sub-Saharan Africa. Three of these four vaccines are lyophilized due to the instability of the serogroup A antigen (MenA) in aqueous solution, resulting in a two vial presentation with concomitant additional costs for storage and distribution. Replacement of the saccharide ring oxygen with a methylene group is a promising approach to preparing a stable oligosaccharide MenA analogue (Carba-MenA) vaccine suitable for a liquid formulation. However, to be effective, Carba-MenA must elicit an immune response that is cross-reactive to the native MenA. Here we employ microsecond molecular dynamics simulations of ten repeats of MenA and Carba-MenA to establish that there are significant differences in the conformation and dynamics of these antigens in solution. Carba-MenA has a more random extended, conformation than MenA; MenA has a significant population of compact S-bend conformations that are absent in the analogue. We also find that the disaccharides are poor models of the conformational behaviour of longer chains. This information is relevant for the rational design of optimal analogues for conjugate vaccines.

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.

Inter-residual Hydrogen Bonding in Carbohydrates Unraveled by NMR Spectroscopy and Molecular Dynamics Simulations

Carbohydrates, also known as glycans in biological systems, are omnipresent in nature where they as glycoconjugates occur as oligo- and polysaccharides linked to lipids and proteins. Their three-dimensional structure is defined by two or three torsion angles at each glycosidic linkage. In addition, transglycosidic hydrogen bonding between sugar residues may be important. Herein we investigate the presence of these inter-residue interactions by NMR spectroscopy in D₂ O/[D₆ ]DMSO (70:30) or D₂ O and by molecular dynamics (MD) simulations with explicit water as solvent for disaccharides with structural elements α-d-Manp-(1→2)-d-Manp, β-d-GlcpNAc-(1→2)-d-Manp, and α-d-Glcp-(1→4)-β-d-Glcp, all of which have been suggested to exhibit inter-residue hydrogen bonding. For the disaccharide β-d-GlcpNAc-(1→2)-β-d-Manp-OMe, the large extent of O5'⋅⋅⋅HO3 hydrogen bonding as seen from the MD simulation is implicitly supported by the ¹ H NMR chemical shift and ³ J_HO3,H3 value of the hydroxy proton. In the case of α-d-Glcp-(1→4)-β-d-Glcp-OMe, the existence of a transglycosidic hydrogen bond O2'⋅⋅⋅HO3 was proven by the presence of a cross-peak in ¹ H,¹³ C HSQC-TOCSY experiments as a result of direct TOCSY transfer between HO3 of the reducing end residue and H2' (detected at C2') of the terminal residue. The occurrence of inter-residue hydrogen bonding, albeit transient, is judged important for the stabilization of three-dimensional structures, which may be essential in maintaining a conformational state for carbohydrate-protein interactions of glycans to take place in biologically important environments.