Molecular dynamics simulation of methyl group relaxation in water

Increasing the Affinity of an O-Antigen Polysaccharide Binding Site in Shigella flexneri Bacteriophage Sf6 Tailspike Protein

Broad and unspecific use of antibiotics accelerates spread of resistances. Sensitive and robust pathogen detection is thus important for a more targeted application. Bacteriophages contain a large repertoire of pathogen-binding proteins. These tailspike proteins (TSP) often bind surface glycans and represent a promising design platform for specific pathogen sensors. We analysed bacteriophage Sf6 TSP that recognizes the O-polysaccharide of dysentery-causing Shigella flexneri to develop variants with increased sensitivity for sensor applications. Ligand polyrhamnose backbone conformations were obtained from 2D 1 H,1 H-trNOESY NMR utilizing methine-methine and methine-methyl correlations. They agreed well with conformations obtained from molecular dynamics (MD), validating the method for further predictions. In a set of mutants, MD predicted ligand flexibilities that were in good correlation with binding strength as confirmed on immobilized S. flexneri O-polysaccharide (PS) with surface plasmon resonance. In silico approaches combined with rapid screening on PS surfaces hence provide valuable strategies for TSP-based pathogen sensor design.

Glycans Confer Specificity to the Recognition of Ganglioside Receptors by Botulinum Neurotoxin A

The highly poisonous botulinum neurotoxins, produced by the bacterium Clostridium botulinum, act on their hosts by a high-affinity association to two receptors on neuronal cell surfaces as the first step of invasion. The glycan motifs of gangliosides serve as initial coreceptors for these protein complexes, whereby a membrane protein receptor is bound. Herein we set out to characterize the carbohydrate minimal binding epitope of the botulinum neurotoxin serotype A. By means of ligand-based NMR spectroscopy, X-ray crystallography, computer simulations, and isothermal titration calorimetry, a screening of ganglioside analogues together with a detailed characterization of various carbohydrate ligand complexes with the toxin were accomplished. We show that the representation of the glycan epitope to the protein affects the details of binding. Notably, both branches of the oligosaccharide GD1a can associate to botulinum neurotoxin serotype A when expressed as individual trisaccharides. It is, however, the terminal branch of GD1a as well as this trisaccharide motif alone, corresponding to the sialyl-Thomsen-Friedenreich antigen, that represents the active ligand epitope, and these compounds bind to the neurotoxin with a high degree of predisposition but with low affinities. This finding does not correlate with the oligosaccharide moieties having a strong contribution to the total affinity, which was expected to be the case. We here propose that the glycan part of the ganglioside receptors mainly provides abundance and specificity, whereas the interaction with the membrane itself and protein receptor brings about the strong total binding of the toxin to the neuronal membrane.

Dynamics of exocyclic groups in the Escherichia coli O91 O-antigen polysaccharide in solution studied by carbon-13 NMR relaxation

Carbon-13 relaxation data are reported for exocyclic groups of hexopyranosyl sugar residues in the repeating unit within the Escherichia coli O91 O-antigen polysaccharide in a dilute D2O solution. The measurements of T 1, T 2 and heteronuclear nuclear Overhauser enhancements were carried out at 310 K at two magnetic fields (16.4 T, 21.1 T). The data were analyzed using the standard and extended Lipari-Szabo models, as well as a conformational jump model. The extended version of the Lipari-Szabo and the two-site jump models were most successful for the hydroxymethyl groups of Gal and GlcNAc sugar residues. Different dynamics was found for the hydroxymethyl groups associated with different configurations (D-gluco, D-galacto) of the sugar residues, the latter being faster than the former.

Steric and electronic influences on the torsional energy landscape of retinal

We have performed quantum mechanical calculations for retinal model compounds to establish the rotational energy barriers for the C5-, C9-, and C13-methyl groups known to play an essential role in rhodopsin activation. Intraretinal steric interactions as well as electronic effects lower the rotational barriers of both the C9- and C13-methyl groups, consistent with experimental (2)H NMR data. Each retinal methyl group has a unique rotational behavior which must be treated individually. These results are highly relevant for the parameterization of molecular mechanics force fields which form the basis of molecular dynamics simulations of retinal proteins such as rhodopsin.

Comparing simulated and experimental translation and rotation constants: range of validity for viscosity scaling

Proper simulation of dynamic properties, including molecular diffusion, is an important goal of empirical force fields. However, the widely used TIP3P water model does not reproduce the experimental viscosity of water. Consequently, scaling of simulated diffusion constants of solutes in aqueous solutions is required to effectively compare them with experiment. It is proposed that scaling by the ratio of viscosities of model and real water is appropriate in the regime where the concentration dependence of simulated and experimental solution viscosities is parallel. With this ansatz, viscosity scaling can be carried out for glucose and trehalose up to 20 wt % for simulations carried out with the CHARMM additive carbohydrate force field C35 and TIP3P water; above this value, the concentration dependence of simulated viscosities lags that of experiment, and scaling is not advised. Scaled translational diffusion constants for glucose and the disaccharides trehalose, maltose, and melibiose at low concentration agree nearly quantitatively with experiment, as do NMR (13)C T(1)'s for glucose, trehalose, and maltose; these results support the use of C35 for simulations of sugar transport properties at low concentration. At high concentrations the scaled diffusion constants for glucose and trehalose underestimate and overestimate experiment, respectively. Hydrodynamic bead model calculations indicate a hydration level of approximately 1 water/hydroxyl for glucose. Patterns for the disaccharides are more complicated, though trehalose binds 0.5 to 1 more water than does maltose depending on the analysis.

Retinal dynamics during light activation of rhodopsin revealed by solid-state NMR spectroscopy

Rhodopsin is a canonical member of class A of the G protein-coupled receptors (GPCRs) that are implicated in many of the drug interventions in humans and are of great pharmaceutical interest. The molecular mechanism of rhodopsin activation remains unknown as atomistic structural information for the active metarhodopsin II state is currently lacking. Solid-state (2)H NMR constitutes a powerful approach to study atomic-level dynamics of membrane proteins. In the present application, we describe how information is obtained about interactions of the retinal cofactor with rhodopsin that change with light activation of the photoreceptor. The retinal methyl groups play an important role in rhodopsin function by directing conformational changes upon transition into the active state. Site-specific (2)H labels have been introduced into the methyl groups of retinal and solid-state (2)H NMR methods applied to obtain order parameters and correlation times that quantify the mobility of the cofactor in the inactive dark state, as well as the cryotrapped metarhodopsin I and metarhodopsin II states. Analysis of the angular-dependent (2)H NMR line shapes for selectively deuterated methyl groups of rhodopsin in aligned membranes enables determination of the average ligand conformation within the binding pocket. The relaxation data suggest that the beta-ionone ring is not expelled from its hydrophobic pocket in the transition from the pre-activated metarhodopsin I to the active metarhodopsin II state. Rather, the major structural changes of the retinal cofactor occur already at the metarhodopsin I state in the activation process. The metarhodopsin I to metarhodopsin II transition involves mainly conformational changes of the protein within the membrane lipid bilayer rather than the ligand. The dynamics of the retinylidene methyl groups upon isomerization are explained by an activation mechanism involving cooperative rearrangements of extracellular loop E2 together with transmembrane helices H5 and H6. These activating movements are triggered by steric clashes of the isomerized all-trans retinal with the beta4 strand of the E2 loop and the side chains of Glu(122) and Trp(265) within the binding pocket. The solid-state (2)H NMR data are discussed with regard to the pathway of the energy flow in the receptor activation mechanism.