Ligand binding to G-quadruplex DNA: new insights from ultraviolet resonance Raman spectroscopy.

G-Quadruplexes (G4s) are noncanonical nucleic acid structures involved in the regulation of several biological processes of many organisms. The rational design of G4-targeting molecules developed as potential anticancer and antiviral therapeutics is a complex problem intrinsically due to the structural polymorphism of these peculiar DNA structures. The aim of the present work is to show how Ultraviolet Resonance Raman (UVRR) spectroscopy can complement other techniques in providing valuable information about ligand/G4 interactions in solution. Here, the binding of BRACO-19 and Pyridostatin - two of the most potent ligands - to selected biologically relevant G4s was investigated by polarized UVRR scattering at 266 nm. The results give new insights into the binding mode of these ligands to G4s having different sequences and topologies by performing an accurate analysis of peaks assigned to specific groups and their changes upon binding. Indeed, the UVRR data not only show that BRACO-19 and Pyridostatin interact with different G4 sites, but also shed light on the ligand and G4 chemical groups really involved in the interaction. In addition, UVRR results complemented by circular dichroism data clearly indicate that the binding mode of a ligand can also depend on the conformation(s) of the target G4. Overall, these findings demonstrate the utility of using UVRR spectroscopy in the investigation of G4s and G4-ligand interactions in solution.

Burkholderia cenocepacia H111 Produces a Water-Insoluble Exopolysaccharide in Biofilm: Structural Determination and Molecular Modelling.

Biofilms are a multicellular way of life, where bacterial cells are close together and embedded in a hydrated macromolecular matrix which offers a number of advantages to the cells. Extracellular polysaccharides play an important role in matrix setup and maintenance. A water-insoluble polysaccharide was isolated and purified from the biofilm produced by Burkholderia cenocepacia strain H111, a cystic fibrosis pathogen. Its composition and glycosidic linkages were determined using Gas-Liquid Chromatography-Mass Spectrometry (GLC-MS) on appropriate carbohydrate derivatives while its complete structure was unraveled by 1D and 2D NMR spectroscopy in deuterated sodium hydroxide (NaOD) aqueous solutions. All the collected data demonstrated the following repeating unit for the water-insoluble B. cenocepacia biofilm polysaccharide: [3)-α-d-Galp-(1→3)-α-d-Glcp-(1→3)-α-d-Galp-(1→3)-α-d-Manp-(1→]n Molecular modelling was used, coupled with NMR Nuclear Overhauser Effect (NOE) data, to obtain information about local structural motifs which could give hints about the polysaccharide insolubility. Both modelling and NMR data pointed at restricted dynamics of local conformations which were ascribed to the presence of inter-residue hydrogen bonds and to steric restrictions. In addition, the good correlation between NOE data and calculated interatomic distances by molecular dynamics simulations validated potential energy functions used for calculations.

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 D2 O/[D6 ]DMSO (70:30) or D2 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 1 H NMR chemical shift and 3 JHO3,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 1 H,13 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.

Crowding and conformation interplay on human DNA G-quadruplex by ultraviolet resonant Raman scattering.

The G-quadruplex-forming telomeric sequence (TTAGGG)4TT was investigated by polarized Ultraviolet Resonance Raman Scattering (UVRR) at 266 nm. The presence of 40% poly(ethylene glycol) and the so-called "self-crowding" condition were used to induce the hybrid-to-parallel topology transition. Analysis of frequency shifts with temperature showed the role of several functional groups in the topological transitions and provides structural dynamical information. Circular dichroism under similar conditions was used as a reference. UVRR shed light on the effect of intramolecular interactions and of local and environmental dynamics in promoting different G-quadruplex topologies, induced by solution conditions or by temperature changes. Overall, these findings showed the enormous potential of this spectroscopy for G-quadruplex conformational studies.

Water structuring above solutes with planar hydrophobic surfaces.

Many important biological solutes possess not only polar and hydrogen bonding functionalities, but also weakly-hydrating, or hydrophobic, surfaces. Theories of the hydration of such surfaces predict that their solvent interactions will change from a wetting type interaction to a dewetting regime as a function of the solute size, with a gradual transition in behavior taking place around characteristic lengths of ∼1 nm. Aggregations of non-polar species over this size range will undergo a transition from being dominated by entropy to being dominated by enthalpy. These transitions can be understood in part in terms of the geometries required of the solvating water molecules. We report here a series of simulations in aqueous solution of organic molecules with planar faces of increasing size, ranging from cyclopropane to circumcircumcoronene, in order to explore the transition in behavior for such solutes as their size increases. For this series, the dewetting transition occurred gradually, converging asymptotically to a limiting separation value for first layer water molecules of around 3.3 Å, while the transition in hydrogen bonding orientational structure occurred between cyclopropane and cyclopentadene. Water immediately adjacent to the largest planar hydrophobic surfaces oriented in ways that resembled on average the structural organization of the basal planes of ice.

Strategies to reduce end-product inhibition in family 48 glycoside hydrolases.

Family 48 cellobiohydrolases are some of the most abundant glycoside hydrolases in nature. They are able to degrade cellulosic biomass and therefore serve as good enzyme candidates for biofuel production. Family 48 cellulases hydrolyze cellulose chains via a processive mechanism, and produce end products composed primarily of cellobiose as well as other cellooligomers (dp ≤ 4). The challenge of utilizing cellulases in biofuel production lies in their extremely slow turnover rate. A factor contributing to the low enzyme activity is suggested to be product binding to enzyme and the resulting performance inhibition. In this study, we quantitatively evaluated the product inhibitory effect of four family 48 glycoside hydrolases using molecular dynamics simulations and product expulsion free-energy calculations. We also suggested a series of single mutants of the four family 48 glycoside hydrolases with theoretically reduced level of product inhibition. The theoretical calculations provide a guide for future experimental studies designed to produce mutant cellulases with enhanced activity.

Simulation studies of substrate recognition by the exocellulase CelF from Clostridium cellulolyticum.

Molecular dynamics (MD) simulations were used to study substrate recognition by the family 48 exocellulase CelF from Clostridium cellulolyticum. It was hypothesized that residues around the entrance of the active site tunnel of this enzyme might serve to recognize and bind the substrate through an affinity for the cellulose monomer repeat unit, ß-d-glucopyranose. Simulations were conducted of the catalytic domain of this enzyme surrounded by a concentrated solution of ß-d-glucopyranose, and the full three-dimensional probability distribution for finding sugar molecules adjacent to the enzyme was calculated from the trajectory. A significant probability of finding the sugar stacked against the planar faces of Trp 310 and Trp 312 at the entrance of the active site tunnel was observed. Biotechnol. Bioeng. 2016;113: 1433-1440. © 2015 Wiley Periodicals, Inc.

Myelography iodinated contrast media. I. Unraveling the atropisomerism properties in solution.

The present work reports a thorough conformational analysis of iodinated contrast media: iomeprol, iopamidol (the world's most utilized contrast agent), and iopromide. Its main aim is the understanding of the complex structural features of these atropisomeric molecules, characterized by the presence of many conformers with hindered rotations, and of the role of atropisomerism in the physicochemical properties of their aqueous solutions. The problem was tackled by using an extensive analysis of (13)C NMR data on the solutions of whole molecules and of simple precursors in addition to FT-IR investigation and molecular simulations. This analysis demonstrated that out of the many possible atropisomers, only a few are significantly populated, and their relative population is provided. The conformational analysis also indicated that the presence of a sterically hindered amidic bond, allowing a significant population of cis forms (E in iopromide and exo in iomeprol), may be the basis for an increased thermodynamic solubility of concentrated solutions of iomeprol.

The conformation of a ribose derivative in aqueous solution: a neutron-scattering and molecular dynamics study.

The structure of aqueous solutions of methyl ß-D-ribofuranoside was investigated by coupling molecular dynamics (MD) simulations and neutron scattering measurements with isotopic substitution. Using a sample of the sugar isotopically-labeled at a single unique position, neutron scattering structure factors and radial distribution functions can be compared with MD simulations constrained to different conformations to determine which conformer best fits the experimental results. Three different simulations were performed with the methyl ether group of the sugar unconstrained and constrained in each of its staggered orientations. The results of the unconstrained simulation showed that the methyl ester group occupied predominantly the 300° position, which is in agreement with the diffraction experimental results. This result suggests that the molecular mechanics force field used in the simulation adequately describes the conformation of the 1-methyl ether group in the methyl ß-D-ribofuranoside.

Potential of mean force conformational energy maps for disaccharide linkages of the Burkholderia multivorans exopolysaccharide C1576 in aqueous solution.

Potential of Mean Force Ramachandran energy maps in aqueous solution have been prepared for all of the glycosidic linkages found in the C1576 exopolysaccharide from the biofilms of the bacterial species Burkholderia multivorans, a member of the Burkholderia cepacian complex that was isolated from a cystic fibrosis patient. C1576 is a rhamnomannan with a tetrasaccharide repeat unit. In general, for the four linkage types in this polymer, hydration did not produce dramatic changes in the Ramachandran energy surfaces, with the 3-methyl-α-d-rhamnopyranose-(1→3)-α-d-rhamnopyranose case exhibiting the greatest hydration change, with the global minimum energy conformation shifting by more than 80° in ψ. However, hydration did reduce the rigidity of all the linkages, increasing the overall flexibility of this polysaccharide.

Microscopy and modelling investigations on the morphology of the biofilm exopolysaccharide produced by Burkholderia multivorans strain C1576.

Bacteria form very often biofilms where they embed in a self-synthesized matrix exhibiting a gel-like appearance. Matrices offer several advantages, including defence against external threats and the easiness of intercellular communication. In infections, biofilm formation enhances bacteria resistance against antimicrobials, causing serious clinical problems for patients' treatments. Biofilm matrices are composed of proteins, extracellular DNA, and polysaccharides, the latter being the major responsible for matrix architecture. The repeating unit of the biofilm polysaccharide synthesized by Burkholderia multivorans strain C1576 contains two mannoses and two sequentially linked rhamnoses, one of them 50 % methylated on C-3. Rhamnose, a 6-deoxysugar, has lower polarity than other common monosaccharides and its methylation further reduces polarity. This suggests a possible role of this polysaccharide in the biofilm matrix; in fact, computer modelling and atomic force microscopy studies evidenced intra- and inter-molecular non-polar interactions both within polysaccharides and with aliphatic molecules. In this paper, the polysaccharide three-dimensional morphology was investigated using atomic force microscopy in both solid and solution states. Independent evidence of the polymer conformation was obtained by transmission electron microscopy which confirmed the formation of globular compact structures. Finally, data from computer dynamic simulations were used to model the three-dimensional structure.

Weakly hydrated surfaces and the binding interactions of small biological solutes.

Extended planar hydrophobic surfaces, such as are found in the side chains of the amino acids histidine, phenylalanine, tyrosine, and tryptophan, exhibit an affinity for the weakly hydrated faces of glucopyranose. In addition, molecular species such as these, including indole, caffeine, and imidazole, exhibit a weak tendency to pair together by hydrophobic stacking in aqueous solution. These interactions can be partially understood in terms of recent models for the hydration of extended hydrophobic faces and should provide insight into the architecture of sugar-binding sites in proteins.

Glucose interactions with a model peptide.

Molecular dynamics simulations have been conducted of the helical polypeptide melittin, in concentrated aqueous solutions of the alpha and beta anomers of D-glucopyranose. Glucose is an osmolyte, and it is expected to be preferentially excluded from the surfaces of proteins. This was indeed found to be the case in the simulations. The results indicate that the observed exclusion may have a contribution from an under-representation of hydrogen bonding interactions between glucose groups and exposed side chains, compared to water. However, glucose was found to bind quite specifically to melittin by stacking its hydrophobic face, consisting of aliphatic protons, against the flat hydrophobic face of the indole group of the tryptophan-19 side chain. Although the binding site for this interaction is localized, the binding is weak for both anomers, with a binding free energy estimated as only ∼0.5 kcal/mol (i.e. near k(B)T). The face of the sugar stacked against the Trp indole ring is different for the two anomers of glucose, due to the disruption of the H1-H3-H5 hydrophobic triad of the beta anomer by the axial C1 hydroxyl group in the alpha anomer. The measurable affinity of the sugar for the Trp side chain is consistent with the very frequent occurrence of this group in the binding sites of proteins that complex with sugars.

The biofilm of Burkholderia cenocepacia H111 contains an exopolysaccharide composed of l-rhamnose and l-mannose: Structural characterization and molecular modelling.

Burkholderia cenocepacia belongs to the Burkholderia Cepacia Complex, a group of 22 closely related species both of clinical and environmental origin, infecting cystic fibrosis patients. B. cenocepacia accounts for the majority of the clinical isolates, comprising the most virulent and transmissible strains. The capacity to form biofilms is among the many virulence determinants of B. cenocepacia, a characteristic that confers enhanced tolerance to some antibiotics, desiccation, oxidizing agents, and host defenses. Exopolysaccharides are a major component of biofilm matrices, particularly providing mechanical stability to biofilms. Recently, a water-insoluble exopolysaccharide produced by B. cenocepacia H111 in biofilm was characterized. In the present study, a water-soluble exopolysaccharide was extracted from B. cenocepacia H111 biofilm, and its structure was determined by GLC-MS, NMR and ESI-MS. The repeating unit is a linear rhamno-tetrasaccharide with 50% replacement of a 3-α-L-Rha with a α-3-L-Man. [2)-α-L-Rhap-(1→3)-α-L-[Rhap or Manp]-(1→3)-α-L-Rhap-(1→2)-α-L-Rhap-(1→]n Molecular modelling was used to obtain information about local structural motifs which could give information about the polysaccharide conformation.

The O-glycosylated linker from the Trichoderma reesei Family 7 cellulase is a flexible, disordered protein.

Fungi and bacteria secrete glycoprotein cocktails to deconstruct cellulose. Cellulose-degrading enzymes (cellulases) are often modular, with catalytic domains for cellulose hydrolysis and carbohydrate-binding modules connected by linkers rich in serine and threonine with O-glycosylation. Few studies have probed the role that the linker and O-glycans play in catalysis. Since different expression and growth conditions produce different glycosylation patterns that affect enzyme activity, the structure-function relationships that glycosylation imparts to linkers are relevant for understanding cellulase mechanisms. Here, the linker of the Trichoderma reesei Family 7 cellobiohydrolase (Cel7A) is examined by simulation. Our results suggest that the Cel7A linker is an intrinsically disordered protein with and without glycosylation. Contrary to the predominant view, the O-glycosylation does not change the stiffness of the linker, as measured by the relative fluctuations in the end-to-end distance; rather, it provides a 16 Å extension, thus expanding the operating range of Cel7A. We explain observations from previous biochemical experiments in the light of results obtained here, and compare the Cel7A linker with linkers from other cellulases with sequence-based tools to predict disorder. This preliminary screen indicates that linkers from Family 7 enzymes from other genera and other cellulases within T. reesei may not be as disordered, warranting further study.

Free energy surfaces for the interaction of D-glucose with planar aromatic groups in aqueous solution.

Multidimensional potentials of mean force for the interactions in aqueous solution of both anomers of D-glucopyranose with two planar aromatic molecules, indole and para-methyl-phenol, have been calculated using molecular dynamics simulations with umbrella sampling and were subsequently used to estimate binding free energies. Indole and para-methyl-phenol serve as models for the side chains of the amino acids tryptophan and tyrosine, respectively. In all cases, a weak affinity between the glucose molecules and the flat aromatic surfaces was found. The global minimum for these interactions was found to be for the case when the pseudoplanar face of ß-D-glucopyranose is stacked against the planar surfaces of the aromatic residues. The calculated binding free energies are in good agreement with both experiment and previous simulations. The multidimensional free energy maps suggest a mechanism that could lend kinetic stability to the complexes formed by sugars bound to sugar-binding proteins.

Ramachandran conformational energy maps for disaccharide linkages found in Burkholderia multivorans biofilm polysaccharides.

Ramachandran conformational energy maps have been prepared for all of the glycosidic linkages found in the C1576 exopolysaccharide that constitutes the biofilms of the bacterial species Burkholderia multivorans, a member of the Burkholderia cepacian complex that was isolated from a cystic fibrosis patient. This polysaccharide is a rhamnomannan with a tetrasaccharide repeat unit containing two mannose residues and two rhamnose residues, -[3-α-d-Man-(1→2)-α-d-Man-(1→2)-α-d-Rha-(1→3)-α-d-Rha-(1→]n-, where approximately 50% of the rhamnoses are randomly methylated on their O3 hydroxyl groups, further increasing the overall hydrophobicity of the chains. Because of the methylation, the tetrasaccharide repeat unit actually contains six possible linkages. The conformational energy maps are fully adiabatic relaxed maps in which the energy for each (ϕ,ψ) grid point on the map represents the lowest possible energy for the molecule in that conformation, considering all the combinations of the other degrees of freedom, such as hydroxyl orientations. Molecular dynamics simulations were used to verify that these maps indeed describe the conformational dynamics of these linkages. All six linkages were found to be quite restricted in possible ϕ angles, but to exhibit several possible low-energy ψ angles, suggesting that these chains could be quite flexible.

Preferential interactions of guanidinum ions with aromatic groups over aliphatic groups.

Small angle neutron scattering (SANS) and molecular dynamics (MD) simulations were used to characterize the long-range structuring (aggregation) of aqueous solutions of isopropanol (IPA) and pyridine and the effect on structuring of guanidinium chloride (GdmCl). These solutes serve as highly soluble analogs of the nonpolar aliphatic (IPA) and aromatic (pyridine) side chains of proteins. SANS data showed that isopropanol and pyridine both form clusters in water resulting from interaction between nonpolar groups of the solutes, with pyridine aggregation producing longer-range structuring than isopropanol in 3 m solutions. Addition of GdmCl at 3 m concentration considerably reduced pyridine aggregation but had no effect on isopropanol aggregation. MD simulations of these solutions support the conclusion that long-range structuring involves hydrophobic solute interactions and that Gdm(+) interacts with the planar pyridine group to suppress pyridine-pyridine interactions in solution. Hydrophobic interactions involving the aliphatic groups of isopropanol were unaffected by GdmCl, indicating that the planar and weakly hydrated Gdm(+) cation cannot make productive interactions with the highly curved or "lumpy" aliphatic groups of this solute. These observations support the conclusion that the effects of Gdm(+) ions on protein-stabilizing interactions involving aromatic amino acid side chains make significant contributions to the denaturant activity of GdmCl, whereas interactions with the "lumpy" aliphatic side chains are likely to be less important.

Specificity of ion-protein interactions: complementary and competitive effects of tetrapropylammonium, guanidinium, sulfate, and chloride ions.

The interactions of ions with a model peptide (a single melittin alpha-helix) in solutions of tetrapropylammonium sulfate or guanidinium chloride were examined by molecular dynamics simulations. The tetrapropylammonium cation shares the geometrical property of essentially flat faces with the previously examined guanidinium cation, and it was found that that this geometry leads to a strong preference for tetrapropylammonium to interact in a similar stacking-type fashion with flat nonpolar groups such as the indole side chain of tryptophan. In contrast to guanidinium, however, tetrapropylammonium does not exhibit strong ion pairing or clustering with sulfate counterions in the solution. Sulfate was found to interact almost exclusively and strongly with the cationic groups of the peptide, such that, already in a 0.1 m solution of tetrapropylammonium sulfate, the 6+ charge of the peptide is effectively locally neutralized. In combination with previous simulations, neutron scattering studies, and experiments on the conformational stability of model peptides, the present results suggest that the Hofmeister series can be explained in higher detail by splitting ions according to the effect they have on hydrogen bonding, salt bridges, and hydrophobic interactions in the protein and how these effects are altered by the counterion.

The polysaccharide extracted from the biofilm of Burkholderia multivorans strain C1576 binds hydrophobic species and exhibits a compact 3D-structure.

Microorganisms often grow in communities called biofilms where cells are imbedded in a complex self-produced biopolymeric matrix composed mainly of polysaccharides, proteins, and DNA. This matrix, together with cell proximity, confers many advantages to these microbial communities, but also constitutes a serious concern when biofilms develop in human tissues or on implanted prostheses. Although polysaccharides are considered the main constituents of the matrices, their specific role needs to be clarified. We have investigated the chemical and morphological properties of the polysaccharide extracted from biofilms produced by the C1576 reference strain of the opportunistic pathogen Burkholderia multivorans, which causes lung infections in cystic fibrosis patients. The aim of the present study is the definition of possible interactions of the polysaccharide and the three-dimensional conformation of its chain within the biofilm matrix. Surface plasmon resonance experiments confirmed the ability of the polysaccharide to bind hydrophobic molecules, due to the presence of rhamnose dimers in its primary structure. In addition, atomic force microscopy studies evidenced an extremely compact three-dimensional structure of the polysaccharide which may form aggregates, suggesting a novel view of its structural role into the biofilm matrix.