Molecular dynamics studies of side chain effect on the beta-1,3-D-glucan triple helix in aqueous solution

beta-1,3-D-glucans have been isolated from fungi as right-handed 6(1) triple helices. They are categorized by the side chains bound to the main triple helix through beta-(1-->6)-D-glycosyl linkage. Indeed, since a glucose-based side chain is water soluble, the presence and frequency of glucose-based side chains give rise to significant variation in the physical properties of the glucan family. Curdlan has no side chains and self-assembles to form an water-insoluble triple helical structure, while schizophyllan, which has a 1,6-D-glucose side chain on every third glucose unit along the main chain, is completely water soluble. A thermal fluctuation in the optical rotatory dispersion is observed for the side chain, indicating probable co-operative interaction between the side chains and water molecules. This paper documents molecular dynamics simulations in aqueous solution for three models of the beta-1,3-D-glucan series: curdlan (no side chain), schizophyllan (a beta-(1-->6)-D-glycosyl side-chain at every third position), and a hypothetical triple helix with a side chain at every sixth main-chain glucose unit. A decrease was observed in the helical pitch as the population of the side chain increased. Two types of hydrogen bonding via water molecules, the side chain/main chain and the side chain/side chain hydrogen bonding, play an important role in determination of the triple helix conformation. The formation of a one-dimensional cavity of diameter about 3.5 A was observed in the schizophyllan triple helix, while curdlan showed no such cavity. The side chain/side chain hydrogen bonding in schizophyllan and the hypothetical beta-1,3-D-glucan triple helix could cause the tilt of the main-chain glucose residues to the helix.     

How hydrogen bonds impact P-glycoprotein transport and permeability

The requirement to cross a biological membrane can be a complex process especially if multidrug transporters such as P-gp must be considered. Drug partitioning into the lipid membrane and efflux by P-gp are tightly coupled processes wherein H-bonding interactions play a key role. All H-bond donors and acceptors are not equal in terms of the strength of the H-bonds that they form, hence it is important to consider their relative strength. Using various examples from literature, we illustrate the benefits of considering the relative strengths of individual H-bonds and introducing intramolecular H-bonds to increase membrane permeability and/or decrease P-gp efflux.     

Side-Chain to Main-Chain Hydrogen Bonding Controls the Intrinsic Backbone Dynamics of the Amyloid Precursor Protein Transmembrane Helix

Many transmembrane helices contain serine and/or threonine residues whose side chains form intrahelical H-bonds with upstream carbonyl oxygens. Here, we investigated the impact of threonine side-chain/main-chain backbonding on the backbone dynamics of the amyloid precursor protein transmembrane helix. This helix consists of a N-terminal dimerization region and a C-terminal cleavage region, which is processed by γ-secretase to a series of products. Threonine mutations within this transmembrane helix are known to alter the cleavage pattern, which can lead to early-onset Alzheimer’s disease. Circular dichroism spectroscopy and amide exchange experiments of synthetic transmembrane domain peptides reveal that mutating threonine enhances the flexibility of this helix. Molecular dynamics simulations show that the mutations reduce intrahelical amide H-bonding and H-bond lifetimes. In addition, the removal of side-chain/main-chain backbonding distorts the helix, which alters bending and rotation at a diglycine hinge connecting the dimerization and cleavage regions. We propose that the backbone dynamics of the substrate profoundly affects the way by which the substrate is presented to the catalytic site within the enzyme. Changing this conformational flexibility may thus change the pattern of proteolytic processing.     

Hydrogen-bonding interactions in adrenaline–water complexes: DFT and QTAIM studies of structures, properties, and topologies

ωB97XD/6-311++G(d,p) calculations were carried out to investigate the hydrogen-bonding interactions between adrenaline (Ad) and water. Six Ad-H(2)O complexes possessing various types of hydrogen bonds (H-bonds) were characterized in terms of their geometries, energies, vibrational frequencies, and electron-density topology. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses were performed to elucidate the nature of the hydrogen-bonding interactions in these complexes. The intramolecular H-bond between the amino and carboxyl oxygen atom of Ad was retained in most of the complexes, and cooperativity between the intra- and intermolecular H-bonds was present in some of the complexes. H-bonds in which hydroxyls of Ad/water acted as proton donors were stronger than other H-bonds. Both hydrogen-bonding interactions and structural deformation play important roles in the relative stabilities of the complexes. The intramolecular H-bond was broken during the formation of the most stable complex, which indicates that Ad tends to break the intramolecular H-bond and form two new intermolecular H-bonds with the first water molecule.     

Crystal structures of substrate-free and nitrosyl cytochrome p450cin: implications for o(2) activation

The crystal structure of the P450cin substrate-bound nitric oxide complex and the substrate-free form have been determined revealing a substrate-free structure that adopts an open conformation relative to the substrate-bound structure. The region of the I helix that forms part of the O(2) binding pocket shifts from an α helix in the substrate-free form to a π helix in the substrate-bound form. Unique to P450cin is an active site residue, Asn242, in the I helix that H-bonds with the substrate. In most other P450s this residue is a Thr and plays an important role in O(2) activation by participating in an H-bonding network required for O(2) activation. The π/α I helix transition results in the carbonyl O atom of Gly238 moving in to form an H-bond with the water/hydroxide ligand in the substrate-free form. The corresponding residue, Gly248, in the substrate-free P450cam structure experiences a similar motion. Most significantly, in the oxy-P450cam complex Gly248 adopts a position midway between the substrate-free and -bound states. A comparison between these P450cam and the new P450cin structures provides insights into differences in how the two P450s activate O(2). The structure of P450cin complexed with nitric oxide, a close mimic of the O(2) complex, shows that Gly238 is likely to form tighter interactions with ligands than the corresponding Gly248 in P450cam. Having a close interaction between an H-bond acceptor, the Gly238 carbonyl O atom, and the distal oxygen atom of O(2) will promote protonation and hence further reduction of the oxy complex to the hydroperoxy intermediate resulting in heterolytic cleavage of the peroxide O-O bond and formation of the active ferryl intermediate required for substrate hydroxylation.     

Chiral recognition ability of curdlan triacetate: solvent and temperature effects

The chiral recognition ability of the chiral stationary phase (CSP) consisting of curdlan (beta-1,3-glucan) triacetate coated on silica gel was clearly changed by the contacting solvents and heat treatment. The chiral recognition ability significantly decreased, particularly at temperatures above 45 degrees C, depending on the racemates. The molecular weight of the curdlan triacetate slightly influenced its ability. The recognition abilities of curdlan tricetate that was lost by heat treatment were partially recovered by contact with methanol. However, when it was contacted with ethanol a different selectivity was observed. The labile chiral recognition ability of curdlan triacetate is in striking contrast to the very stable chiral recognition of cellulose (beta-1,4-glucan) triacetate (Chiralcel OA). This difference may be ascribed to the conformational stability of the acetates consisting of curdlan (beta-1,3-glucan) and cellulose (beta-1,4-glucan) with different sugar linkages.     

Structural role of (bacterio)chlorophyll ligated in the energetically unfavorable beta-position

Chlorophyll is attached to apoprotein in diastereotopically distinct ways, by beta- and alpha-ligation. Both the beta- and alpha-ligated chlorophylls of photosystem I are shown to have ample contacts to apoprotein within their proteinaceous binding sites, in particular, at C-13 of the isocyclic ring. The H-bonding patterns for the C-13(1) oxo groups, however, are clearly distinct for the beta-ligated and alpha-ligated chlorophylls. The beta-ligated chlorophylls frequently employ their C-13(1) oxo in H-bonds to neighboring helices and subunits. In contrast, the C-13(1) oxo of alpha-ligated chlorophylls are significantly less involved in H-bonding interactions, particularly to neighboring helices. Remarkably, in the peripheral antenna, light harvesting complex (LH2) from Rhodobacter sphaeroides, a single mutation in the alpha-subunit, introduced to eliminate H-bonding to the beta-bacteriochlorophyll-B850, which is ligated in the "beta-position," results in significant thermal destabilization of the LH2 in the membrane. In addition, in comparison with wild type LH2, the expression level of the LH2 lacking this H-bond is significantly reduced. These findings show that H-bonding to the C-13(1) keto group ofbeta-ligated (bacterio)-chlorophyll is a key structural motif and significantly contributes to the stability of bacteriochlorophyll proteins in the native membrane. Our analysis of photosystem I and II suggests that this hitherto unrecognized motif involving H-bonding to beta-ligated chlorophylls may be equally critical for the stable assembly of the inner core antenna of these multicomponent chlorophyll proteins.     

Side-Chain to Main-Chain Hydrogen Bonding Controls the Intrinsic Backbone Dynamics of the Amyloid Precursor Protein Transmembrane Helix

Many transmembrane helices contain serine and/or threonine residues whose side chains form intrahelical H-bonds with upstream carbonyl oxygens. Here, we investigated the impact of threonine side-chain/main-chain backbonding on the backbone dynamics of the amyloid precursor protein transmembrane helix. This helix consists of a N-terminal dimerization region and a C-terminal cleavage region, which is processed by γ-secretase to a series of products. Threonine mutations within this transmembrane helix are known to alter the cleavage pattern, which can lead to early-onset Alzheimer’s disease. Circular dichroism spectroscopy and amide exchange experiments of synthetic transmembrane domain peptides reveal that mutating threonine enhances the flexibility of this helix. Molecular dynamics simulations show that the mutations reduce intrahelical amide H-bonding and H-bond lifetimes. In addition, the removal of side-chain/main-chain backbonding distorts the helix, which alters bending and rotation at a diglycine hinge connecting the dimerization and cleavage regions. We propose that the backbone dynamics of the substrate profoundly affects the way by which the substrate is presented to the catalytic site within the enzyme. Changing this conformational flexibility may thus change the pattern of proteolytic processing.     

Representation of DNA sequences by use of helical wheels to identify pattern formation with protein binding potential

Similar to the protein segments with helical potential the three-dimensional structures of DNA sequences can be represented by a two-dimensional projection along the axis of the double helix. We will call this projection according to Schiffer and Edmundson "helical wheel". In the first place the wheels are employed here as graphical restatements of known sequences which represent binding sites of restriction-modification enzymes. Specific recognition of these sites has to be based on sequence related multi contact interactions of weak bonds. Because of the spatial and directional characteristics of H-bonds we will utilize the pattern formation of H-bond donor--and acceptor groups protruding into the major groove of the helix to identify the particular sequences. The characteristics of these patterns can be readily visualized and compared by examination of these wheels.     

Novel carbohydrate-binding module of beta-1,3-xylanase from a marine bacterium, Alcaligenes sp. strain XY-234

A beta-1,3-xylanase gene (txyA) from a marine bacterium, Alcaligenes sp. strain XY-234, has been cloned and sequenced. txyA consists of a 1,410-bp open reading frame that encodes 469 amino acid residues with a calculated molecular mass of 52,256 Da. The domain structure of the beta-1,3-xylanase (TxyA) consists of a signal peptide of 22 amino acid residues, followed by a catalytic domain which belongs to family 26 of the glycosyl hydrolases, a linker region with one array of DGG and six repeats of DNGG, and a novel carbohydrate-binding module (CBM) at the C terminus. The recombinant TxyA hydrolyzed beta-1,3-xylan but not other polysaccharides such as beta-1,4-xylan, carboxymethylcellulose, curdlan, glucomannan, or beta-1,4-mannan. TxyA was capable of binding specifically to beta-1,3-xylan. The analysis using truncated TxyA lacking either the N- or C-terminal region indicated that the region encoding the CBM was located between residues 376 and 469. Binding studies on the CBM revealed that the K(d) and the maximum amount of protein bound to beta-1,3-xylan were 4.2 microM and 18.2 micromol/g of beta-1,3-xylan, respectively. Furthermore, comparison of the enzymatic properties between proteins with and without the CBM strongly indicated that the CBM of TxyA plays an important role in the hydrolysis of beta-1,3-xylan.     

Interaction of nucleic acid bases with water molecules and formation of mismatched nucleotide pairs

The possibility of the inclusion of water molecules in the formation of mismatched nucleotide pairs was considered in relation to the mechanisms of point errors in template directed biosynthesis. Calculations of the intermolecular interaction energy for systems containing two bases and one water molecule were carried out by the method of atom-atom potential functions. There exist energy minima for each base pair, corresponding to a single N--H...O or N--H...N H-bond between the bases and H-bonding of the water molecule with both bases. The relative positions of glycosyl bonds in some of these minima are closer to those for Watson--Crick pairs, than the positions of minima for these pairs without water. For other minima, the H-bond formation between the water molecule and the two bases additionally stabilizes the relative base position in wobble-pairs with two H-bonds between the bases. The base and water positions in energy minima are compared with the positions in some pairs proposed on the basis of NMR and X-ray data for double helical oligonucleotides.     

NMR determination of the global structure of the 113Cd derivative of desulforedoxin: investigation of the hydrogen bonding pattern at the metal center.

Desulforedoxin (Dx) is a simple homodimeric protein isolated from Desulfovibrio gigas (Dg) containing a distorted rubredoxin-like center with one iron coordinated by four cysteinyl residues (7.9 kDa with 36 amino acids per monomer). In order to probe the geometry and the H-bonding at the active site of Dx, the protein was reconstituted with 113Cd and the solution structure determined using 2D NMR methods. The structure of this derivative was initially compared with the NMR solution structure of the Zn form (Goodfellow BJ et al., 1996, J Biol Inorg Chem 1:341-353). Backbone amide protons for G4, D5, G13, L11 NH, and the Q14 NH side-chain protons, H-bonded in the X-ray structure, were readily exchanged with solvent. Chemical shift differences observed for amide protons near the metal center confirm the H-bonding pattern seen in the X-ray model (Archer M et al., 1995, J Mol Biol 251:690-702) and also suggest that H-bond lengths may vary between the Fe, Zn, and 113Cd forms. The H-bonding pattern was further probed using a heteronuclear spin echo difference (HSED) experiment; the results confirm the presence of NH-S H-bonds inferred from D2O exchange data and observed in the NMR family of structures. The presence of "H-bond mediated" coupling in Dx indicates that the NH-S H-bonds at the metal center have significant covalent character. The HSED experiment also identified an intermonomer "through space" coupling for one of the L26 methyl groups, indicating its proximity to the 113Cd center in the opposing monomer. This is the first example of an intermonomer "through space" coupling. Initial structure calculations produced subsets of NMR families with the S of C28 pointing away from or toward the L26 methyl: only the subset with the C28 sulfur pointing toward the L26 methyl could result in a "through space" coupling. The HSED result was therefore included in the structure calculations. Comparison of the Fe, Zn, and 113Cd forms of Dx suggests that the geometry of the metal center and the global fold of the protein does not vary to any great extent, although the H-bond network varies slightly when Cd is introduced. The similarity between the H-bonding pattern seen at the metal center in Dx, Rd (including H-bonded and through space-mediated coupling), and many zinc-finger proteins suggests that these H-bonds are structurally vital for stabilization of the metal centers in these proteins.     

Conformational studies of sphingolipids by NMR spectroscopy. II. Sphingomyelin

Sphingomyelin (SM) is the most prevalent sphingolipid in the majority of mammalian membranes. Proton and 31P nuclear magnetic resonance spectral data were acquired to establish the nature of intra- and intermolecular H-bonds in the monomeric and aggregated forms of SM and to assess possible differences between this lipid and dihydrosphingomyelin (DHSM), which lacks the double bond between carbons 4 and 5 of the sphingoid base. The spectral trends suggest the formation of an intramolecular H-bond between the OH group of the sphingosine moiety and the phosphate ester oxygen of the head group. The narrower linewidth and the downfield shift of the resonance corresponding to OH proton in SM suggest that this H-bond is stronger in SM than in DHSM. The NH group appears to be involved predominantly in intramolecular H-bonding in the monomer. As the concentration of SM increases and the molecules come in closer proximity, these intramolecular bonds are partially disrupted and the NH group becomes involved in lipid-water interactions. The difference between the SM and DHSM appears to be not in the nature of these interactions but rather in the degree to which these intermolecular interactions prevail. As SM molecules cannot come as close together as DHSM molecules can, both the NH and OH moieties remain, on average, more intramolecularly bonded as compared to DHSM.     

Observation of the closing of individual hydrogen bonds during TFE-induced helix formation in a peptide

Helix formation of an S-peptide analog, comprising the first 20 residues of Ribonuclease A and two additional N-terminal residues, was studied by measuring hydrogen bond (H-bond) (h3)J(NC') scalar couplings as a function of 2,2,2-trifluoroethanol (TFE) concentration. The (h3)J(NC') couplings give direct evidence for the closing of individual backbone N-H***O = C H-bonds during the TFE-induced formation of secondary structure. Whereas no (h3)J(NC') correlations could be detected without TFE, alpha-helical (i,i +4) H-bond correlations were observed for the amides of residues A5 to M15 in the presence of TFE. The analysis of individual coupling constants indicates that alpha-helix formation starts at the center of the S-peptide around residue E11 and proceeds gradually from there to both peptide ends as the TFE concentration is increased. At 60% to 90% TFE, well-formed alpha-helical H-bonds were observed for the amides hydrogens of residues K9 to Q13, whereas H-bonds of residues T5 to A8, H14, and M15 are affected by fraying. No intramolecular backbone H-bonds are present at and beyond the putative helix stop signal D16. As the (h3)J(NC') constants represent ensemble averages and the dependence of (h3)J(NC') on H-bond lengths is very steep, the size of the individual (h3)J(NC') coupling constants can be used as a measure for the population of a closed H-bond. These individual populations are in agreement with results derived from the Lifson-Roig theory for coil-to-helix transitions. The present work shows that the closing of individual H-bonds during TFE-induced helix formation can be monitored by changes in the size of H-bond scalar couplings.     

Beta-1,3-glucan receptor and peptidoglycan receptor are present as separate entities within insect prophenoloxidase activating system

Silkworm plasma was passed over a peptidoglycan-Sepharose 4B column or a CPB column [CPB, curdlan type polysaccharide (beta-1,3-glucan) bead] in the absence of divalent cation and the effluents from the columns were named plasma-PG and plasma-CPB, respectively. Prophenoloxidase activating system in plasma-PG was triggered by beta-1,3-glucan but not by peptidoglycan and the system in plasma-CPB was triggered by peptidoglycan but not by beta-1,3-glucan, suggesting that the peptidoglycan-Sepharose 4B column and the CPB column remove peptidoglycan-receptor and beta-1,3-glucan-receptor, respectively, from plasma. This result indicates that both receptors exist as separate entities in silkworm plasma. It is suggested that plasma-PG and plasma-CPB may be used as specific reagents to detect minute amounts of beta-1,3-glucan and peptidoglycan.     

100ps Molecular Dynamic Simulation of d(TATCACC)2

Abstract

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage X is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H- bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.     

Conformational studies of sphingolipids by NMR spectroscopy. I. Dihydrosphingomyelin

The conformational features of dihydrosphingomyelin (DHSM), the major phospholipid of human lens membranes, were investigated by 1H and 31P nuclear magnetic resonance spectroscopy. Several postulates emerge from the observed trends: (a) in partially hydrated samples of DHSM in CDCl3 above 13 mM, at which lipid-lipid interactions prevail, the amide proton is mostly involved in intermolecular H-bonds that link neighboring phospholipids through bridging water molecules. In the absence of water, the NH group is involved in an intramolecular H-bond that restricts the mobility of the phosphate group. (b) In the monomeric form of the lipid molecule, the amide proton of the major conformer is bound intramolecularly with one of the anionic and/or ester oxygens of the phosphate group. A minor conformer may also be present in which the NH proton participates in an intramolecular H-bond linking to the OH group of the sphingoid base. (c) Complete hydration leads to an extension of the head group as water molecules bind to the phosphate and NH groups via H-bonds, thus disrupting the intramolecular H-bonds prevalent at low concentrations.     

Disordered Single Crystal Evidence for a Quadruple Helix Formed by Guanosine 5′-Monophosphate

Abstract

The disodium salt of guanosine 5′-monophosphate (5′-GMP) has been crystallized earlier in an orthorhombic array. We have obtained a new crystal form of 5′-GMP at pH 8 which reveals a clear helical nature, with guanine bases stacked perpendicular to the helix axis. Although the X-ray pictures show partial disorder, they can be indexed on a hexagonal net with a = b = 28.6 Å,c = 9.8 Å, V= 6942ų(1Å = 0.1 nm). The probable space group is P6₄, and past experience with ca. 600 ų per base in oligonucleotide crystals suggests that the cell contains 12 GMP molecules. The crystal packing parameters and the intensity distribution agree with a model of three hydrogen-bonded guanine tetrads in the unit cell, stacked so as to build a quadruple helix similar to that proposed earlier from fiber studies (Zimmerman, S.B., J. Mol. Biol. 106, 663–672 (1976)).     

A relationship between amide hydrogen bond strength and quinone reduction potential: implications for photosystem I and bacterial reaction center quinone function

A series of 11 simple phylloquinone derivatives, each lacking the extended phytyl side chain but featuring H-bond donor amides at one or both peri positions, were prepared and some salient physical properties were measured. A correlation between both IR frequency and NMR peak position, as indicators of internal H-bond strength, and the quinone half-wave reduction potential, was observed. These data are consistent with the prevailing hypothesis that quinone carbonyl H-bonding in general, and stronger H-bonds in particular, favorably bias the endogenous quinone's electrochemical potential toward easier reduction.     

Control of pH and PIP2 gating in heteromeric Kir4.1/Kir5.1 channels by H-Bonding at the helix-bundle crossing

Inhibition by intracellular H(+) (pH gating) and activation by phosphoinositides such as PIP(2) (PIP(2)-gating) are key regulatory mechanisms in the physiology of inwardly-rectifying potassium (Kir) channels. Our recent findings suggest that PIP(2) gating and pH gating are controlled by an intra-subunit H-bond at the helix-bundle crossing between a lysine in TM1 and a backbone carbonyl group in TM2. This interaction only occurs in the closed state and channel opening requires this H-bond to be broken, thereby influencing the kinetics of PIP(2) and pH gating in Kir channels. In this addendum, we explore the role of H-bonding in heteromeric Kir4.1/Kir5.1 channels. Kir5.1 subunits do not possess a TM1 lysine. However, homology modelling and molecular dynamics simulations demonstrate that the TM1 lysine in Kir4.1 is capable of H-bonding at the helix-bundle crossing. Consistent with this, the rates of pH and PIP2 gating in Kir4.1/Kir5.1 channels (two H-bonds) were intermediate between those of wild-type homomeric Kir4.1 (four H-bonds) and Kir4.1(K67M) channels (no H-bonds) suggesting that the number of H-bonds in the tetrameric channel complex determines the gating kinetics. Furthermore, in heteromeric Kir4.1(K67M)/Kir5.1 channels, where the two remaining H-bonds are disrupted, we found that the gating kinetics were similar to Kir4.1(K67M) homomeric channels despite the fact that these two channels differ considerably in their PIP(2) affinities. This indicates that Kir channel PIP(2) affinity has little impact on either the PIP(2) or pH gating kinetics.