Search

Neutron diffraction with isotopic substitution (NDIS) experiments were done on both natural abundance potassium and isotopically labeled 41KCl heavy water solutions to characterize the solvent structuring around the potassium ion in water. Preliminary measurements suggested that the literature value for the coherent neutron scattering length (2.69 fm) for 41K was significantly in error. This value was remeasured using a neutron powder diffractometer and found to be 2.40 fm. This revision increases significantly the contrast between the natural abundance K and 41K by about 30% (from 1.0 to 1.3 fm). The experimentally determined structure factor of the potassium ion was then compared to that calculated from molecular dynamics (MD) simulations. Previous neutron scattering measurements of potassium gave a solvation number of 5.5 (see below). In this study, the NDIS and MD results are in good agreement and allowed us to derive a coordination number of 6.1 for water molecules and 0.8 for chloride ions around each K+ ion in 4 molal aqueous KCl solution.

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.

Cel48A of Thermobifida fusca (TfCel48A) is a processive exocellulase that contains an active site tunnel and digests lignocellulosic biomass via synergistic interactions between different cellulases. Cel48A possesses a number of aromatic amino acids lining the tunnel entrance, which are highly conserved across a diverse number of microbial species and appear to play a role in the selection and threading of individual strands of cellulose from highly recalcitrant substrates. In this study, we sought to further elucidate the roles of these tunnel entrance aromatic amino acids by creating a series of double mutants and examining their effect on TfCel48A activity, processivity, and synergistic interactions with the well-studied processive endocellulase TfCel9A. Our results provide further insight concerning the mechanism of Cel48A kinetics with soluble and insoluble substrates and could play an influential role in the application of Cel48A and other exocellulases for industrial purposes.

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.

Insight into the molecular interactions of homotactic and heterotactic association of caffeine and pyridine in aqueous solution is given on the basis of both experimental and simulation studies. Caffeine is about 5 times more soluble in a 3 m aqueous pyridine solution than it is in pure water (an increase from ∼0.1 m to 0.5 m). At this elevated concentration the system becomes suitable for neutron scattering study. Caffeine-pyridine interactions were studied by neutron scattering and molecular dynamics simulations, allowing a detailed characterization of the spatial and orientational structure of the solution. It was found that while pyridine-caffeine interactions are not as strong as caffeine-caffeine interactions, the pyridine-caffeine interactions still significantly disrupted caffeine-caffeine stacking. The alteration of the caffeine-caffeine stacking, occasioned by the presence of pyridine molecules in solution and the consequent formation of heterotactic interactions, leads to the experimentally detected increase in caffeine solubility.

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.

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.

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 dynamical and structural properties of caffeine solutions at the solubility limit have been investigated as a function of temperature by means of MD simulations, static and dynamic light scattering, and small angle neutron scattering experiments. A clear picture unambiguously supported by both experiment and simulation emerges: caffeine self-aggregation promotes the formation of two distinct types of clusters: linear aggregates of stacked molecules, formed by 2-14 caffeine molecules depending on the thermodynamic conditions and disordered branched aggregates with a size in the range 1000-3000 Å. While the first type of association is well-known to occur under room temperature conditions for both caffeine and other purine systems, such as nucleotides, the presence of the supramolecular aggregates has not been reported previously. MD simulations indicate that branched structures are formed by caffeine molecules in a T-shaped arrangement. An increase of the solubility limit (higher temperature but also higher concentration) broadens the distribution of cluster sizes, promoting the formation of stacked aggregates composed by a larger number of caffeine molecules. Surprisingly, the effect on the branched aggregates is rather limited. Their internal structure and size do not change considerably in the range of solubility limits investigated.

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.

As water/trehalose structures seem to gain more biological relevance in anhydrobiosis, the intriguing mechanism of bioprotection may involve water channels in the crystalline dihydrate trehalose. This work aims at understanding the water state(s) in the bioprotectant sugar and its removal under “controlled conditions”. The results of some simple calorimetric experiments are presented to corroborate the above findings and to support the role of molecular nanostructures of trehalose possessing the ability to transform reversibly from the dihydrate. In addition to updating some of our previously results, new calorimetric data and MD analysis are provided for the peculiar transitions of TRE-h and for the energetics of water molecules flowing along the TRE-h channels.

Glycoside hydrolase (GH) family 48 is an understudied and increasingly important exoglucanase family found in the majority of bacterial cellulase systems. Moreover, many thermophilic enzyme systems contain GH48 enzymes. Deletion of GH48 enzymes in these microorganisms results in drastic reduction in biomass deconstruction. Surprisingly, given their importance for these microorganisms, GH48s have intrinsically low cellulolytic activity but even in low ratios synergize greatly with GH9 endoglucanases. In this study, we explore the structural and enzymatic diversity of these enzymes across a wide range of temperature optima. We have crystallized one new GH48 module from Bacillus pumilus in a complex with cellobiose and cellohexaose (BpumGH48). We compare this structure to other known GH48 enzymes in an attempt to understand GH48 structure/function relationships and draw general rules correlating amino acid sequences and secondary structures to thermostability in this GH family. Electronic supplementary material The online version of this article (10.1186/s13068-017-0951-5) contains supplementary material, which is available to authorized users.

Many important biological solutes possess not only polar and hydrogen-bonding functionalities but also weakly hydrating, or hydrophobic, surfaces. While the aggregation of these hydrophobic surfaces has been shown to play an important role in the aggregation of individual chains of cellulose, it is not known whether the water structuring imposed by these hydrophobic surfaces more closely resembles that associated with small hydrophobic solutes like methane and fats or more closely resembles that associated with extended hydrophobic surfaces like mica or waxy planes. By using molecular dynamics simulations to characterize the water molecule orientations over different regions of the 100 surface of cellulose in contact with water, it was found that the hydrophobic strips of the cellulose crystal are sufficiently narrow that they hydrate like a fatty acid chain, rather than like a more extended surface, suggesting that their aggregation would be dominated by entropy rather than enthalpy.

Room temperature ionic liquids (ILs) of the imidazolium family have attracted much attention during the past decade for their capability to dissolve biomass. Besides experimental work, numerous compuational studies have been concerned with the physical properties of both neat ILs and their interactions with different solutes, in particular, carbohydrates. Many classical force fields designed specifically for ILs have been found to yield viscosities that are too high for the liquid state, which has been attributed to the fact that the effective charge densities are too high due to the lack of electronic polarizability. One solution to this problem has been uniform scaling of the partial charges by a scale factor in the range 0.6-0.9, depending on model. This procedure has been shown to improve the viscosity of the models, and also to positively affect other properties, such as diffusion constants and ionic conductivity. However, less attention has been paid to how this affects the overall thermodynamics of the system, and the problems it might create when the IL models are combined with other force fields (e.g., for solutes). In the present work, we employ three widely used IL force fields to simulate 1-n-buty1-3-methyl-imidazolium chloride in both the crystal and the liquid state, as well as its binary mixture with ethanol. Two approaches are used: one in which the ionic charge is retained at its full integer value and one in which the partial charges are uniformly reduced to 85%. We investigate and calculate crystal and liquid structures, molar heat capacities, heats of fusion, self-diffusion constants, ionic conductivity, and viscosity for the neat IL, and ethanol activity as a function of ethanol concentration for the binary mixture. We show that properties of the crystal are less affected by charge scaling compared to the liquid. In the liquid state, transport properties of the neat IL are generally improved by scaling, whereas values for the heat of fusion are unaffected, and results for the heat capacity are ambiguous. Neither full nor reduced charges could reproduce experimental ethanol activities for the whole range of compositions. AuthorCount:5

Molecular dynamics simulations were used to study the possible catalytic role of an unusual conserved water-filled pore structure in the family 48 cellulase enzyme Cel48A from Thermobifida fusca. It was hypothesized that this pore serves as the pathway for the water molecules consumed in the hydrolysis catalyzed by the enzyme to reach the active site in a continuous stream to participate in the processive reactions. Theoretical mutants of this enzyme were created in which all of the residues lining the pore were made hydrophobic, which had the effect in molecular dynamics simulations of emptying the pore of water molecules and preventing any from passing through the pore on the simulation time scale. Mutants with smaller numbers of substitutions of this nature, which could be created experimentally by site-directed mutagenesis, were also identified from simulations, and these proteins were subsequently produced in Escherichia coli, expressed and purified, but were found to not fold in a manner similar to the wild type protein, preventing the determination of the importance of the water pore for activity. It is possible that the presence of a small vacuum in the pore was responsible for the instability of the mutants. In addition, alternate pathways were observed in the simulations that would allow water molecules to reach the active site of the enzyme, suggesting that the hypothesis that the pore has functional significance might be incorrect.

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.

Molecular dynamics simulations were carried out on a system of caffeine interacting with the sugar alcohol sorbitol. The system examined had a caffeine concentration 0.083 m and a sugar concentration 1.08 m. The trajectories of all molecules in the system were collected over a period of 80 ns and analyzed to determine whether there is any tendency for sorbitol to bind to caffeine, and if so, by what mechanism. The results show that the sorbitol molecules have an affinity for the caffeine molecules and that the binding occurred by the interaction of the aliphatic hydrophobic protons of the sugar with the caffeine face. This intermolecular association via face-to-face stacking, as suggested by simulation studies, is similar to that found for sucrose and for D-glucose, which overwhelmingly exists in the pyranose ring chair form in aqueous solution, as well as for caffeine-caffeine association. The sorbitol molecules, however, exist as relatively extended chains and are, therefore, topologically quite different from the sugars sucrose and glucose. The comparison of the average conformation of sorbitol molecules bound to caffeine with that of molecules in the free state shows a substantial similarity.

During the last few decades it has been ascertained that base stacking is one of the major contributions stabilizing nucleic acid conformations. However, the understanding of the nature of the interactions involved in the stacking process remains under debate and it is a subject of theoretical and experimental studies. Structural similarity between purine bases (guanine and adenine) in DNA and the caffeine molecule makes caffeine an excellent model for the purine bases. The present study clearly shows that dipolar interactions play a fundamental role in determining stacking of purine molecules in solution. In order to reach this achievement, polarized ultraviolet Raman resonant scattering experiments have been carried out on caffeine aqueous solutions as a function of concentration and temperature. The investigation pointed out at the aggregation and solvation properties, particularly at elevated temperatures. Kubo-Anderson theory was used as a framework to investigate the non-coincidence effect (NCE) occurring in the totally symmetric breathing modes of the purine rings, and in the bending modes of the methyl groups of caffeine. The NCE concentration dependence shows that caffeine aggregation at 80 °C occurs by planar stacking of the hydrophobic faces. The data clearly indicate that dipolar interactions determine the reorientational motion of the molecules in solution and are the driving force for the stacking of caffeine. In parallel, the observed dephasing times imply a change in caffeine interactions as a function of temperature and concentration. A decrease, at low water content, of the dephasing time for the ring breathing vibration mode indicates that self-association alters the solvation structure that is detectable at low concentration. These results are in agreement with simulation predictions and serve as an important validation of the models used in those calculations.

Molecular dynamics simulations were carried out on several systems of caffeine interacting with simple sugars. These included a single caffeine molecule in a 3 molal solution of α-D-glucopyranose, at a caffeine concentration of 0.083 molal; a single caffeine in a 3 molal solution of β-D-glucopyranose, and a single caffeine molecule in a 1.08 molal solution of sucrose (table sugar). Parallel Nuclear Magnetic Resonance titration experiments were carried out on the same solutions under similar conditions. Consistent with previous thermodynamic experiments, the sugars were found to have an affinity for the caffeine molecules in both the simulations and experiments, and that the binding in these complexes occurs by face-to-face stacking of the hydrophobic triad of protons of the pyranose rings against the caffeine face, rather than by hydrogen bonding. For the disaccharide, the binding occurs via stacking of the glucose ring against the caffeine, with a lesser affinity for the fructose observed. These findings are consistent with the association being driven by hydrophobic hydration, and are similar to the previously observed binding of glucose rings to various other planar molecules, including indole, serotonin, and phenol.

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.

A combined solid-state NMR and Molecular Dynamics simulation study of cellulose in urea aqueous solution and in pure water was conducted. It was found that the local concentration of urea is signif ...

Molecular dynamics simulations were carried out on a system of eight independent caffeine molecules in a periodic box of water at 300 K, representing a solution near the solubility limit for caffeine at room temperature, using a newly-developed CHARMM-type force field for caffeine in water. Simulations were also conducted for single caffeine molecules in water using two different water models (TIP3P and TIP4P). Water was found to structure in a complex fashion around the planar caffeine molecules, which was not sensitive to the water model used. As expected, extensive aggregation of the caffeine molecules was observed, with the molecules stacking their flat faces against one another like coins, with their methylene groups staggered to avoid steric clashes. A dynamic equilibrum was observed between large n-mers, including stacks with all eight solute molecules, and smaller clusters, with the calculated osmotic coefficient being in acceptable agreement with the experimental value. The insensitivity of the results to water model and the congruence with experimental thermodynamic data suggest that the observed stacking interactions are a realistic representation of the actual association mechanism in aqueous caffeine solutions.

International audience; 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 kBT). 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.

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.

Neutron diffraction with isotopic substitution (NDIS) experiments have been used to examine the structuring of aqueous solutions of pyridine. A new method is described for extracting the structure factors relating to intermolecular correlations from neutron scattering experiments on liquid solutions of complex molecular species. This approach performs experiments at different concentrations and exploits the intramolecular coordination number concentration invariance (ICNCI) to separate the internal and intermolecular contributions to the total intensities. The ability of this method to deconvolute molecular and intermolecular correlations is tested and demonstrated using simulated neutron scattering results predicted from molecular dynamics simulations of aqueous solutions of the polyatomic solute pyridine in which the inter- and intramolecular terms are known. The method is then implemented using neutron scattering measurements on solutions of pyridine. The results confirm that pyridine shows a significant propensity to aggregate in solution and demonstrate the prospects for the future application of the ICNCI approach to the study of large polyatomic solutes using next-generation neutron sources and detectors.

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 non-polar 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 non-polar groups of the solutes, with pyridine aggregation producing longer-range structuring than isopropanol in 3 molal solutions. Addition of GdmCl at 3 molal 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.

Cellobiohydrolases are the dominant components of the commercially relevant Trichoderma reesei cellulase system. Although natural cellulases can totally hydrolyze crystalline cellulose to soluble sugars, the current enzyme loadings and long digestion times required render these enzymes less than cost effective for biomass conversion processes. It is clear that cellobiohydrolases must be improved via protein engineering to reduce processing costs. To better understand cellobiohydrolase function, new simulations have been conducted using CHARMM of cellobiohydrolase I (CBH I) from T. reesei interacting with a model segment (cellodextrin) of a cellulose microfibril in which one chain from the substrate has been placed into the active site tunnel mimicking the hypothesized configuration prior to final substrate docking (i.e., the +1 and +2 sites are unoccupied), which is also the structure following a catalytic bond scission. No tendency was found for the protein to dissociate from or translate along the substrate surface during this initial simulation, nor to align with the direction of the cellulose chains. However, a tendency for the decrystallized cellodextrin to partially re-anneal into the cellulose surface hints that the arbitrary starting configuration selected was not ideal.

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.

We present an all-atom additive empirical force field for the hexopyranose monosaccharide form of glucose and its diastereomers allose, altrose, galactose, gulose, idose, mannose, and talose. The model is developed to be consistent with the CHARMM all-atom biomolecular force fields, and the same parameters are used for all diastereomers, including both the alpha- and beta-anomers of each monosaccharide. The force field is developed in a hierarchical manner and reproduces the gas-phase and condensed-phase properties of small-molecule model compounds corresponding to fragments of pyranose monosaccharides. The resultant parameters are transferred to the full pyranose monosaccharides, and additional parameter development is done to achieve a complete hexopyranose monosaccharide force field. Parametrization target data include vibrational frequencies, crystal geometries, solute-water interaction energies, molecular volumes, heats of vaporization, and conformational energies, including those for over 1800 monosaccharide conformations at the MP2/cc-pVTZ//MP2/6-31G(d) level of theory. Although not targeted during parametrization, free energies of aqueous solvation for the model compounds compare favorably with experimental values. Also well-reproduced are monosaccharide crystal unit cell dimensions and ring pucker, densities of concentrated aqueous glucose systems, and the thermodynamic and dynamic properties of the exocyclic torsion in dilute aqueous systems. The new parameter set expands the CHARMM additive force field to allow for simulation of heterogeneous systems that include hexopyranose monosaccharides in addition to proteins, nucleic acids, and lipids.

Guanidinium (Gdm+) chloride is a powerful protein denaturant, whereas the sulfate dianion (SO42-) is a strong stabilizer of folded protein states; Gdm2SO4 is effectively neutral in its effects on protein stability. While the "neutralizing" effects of protein-stabilizing solutes on the activity of denaturants can be broadly interpreted in terms of additive effects of the solutes, recent experimental and simulation studies support a role for hetero-ion interactions in the effect of sulfate on Gdm+ denaturation [Mason, P. E.; et al. J. Phys. Chem. B 2005, 109, 24185-24196]. Here we describe an experimental strategy for testing this mechanism that involves spectroscopic analysis of the separate effects of alkali metal sulfates (Na2SO4, Rb2SO4), GdmCl, and Gdm2SO4 on the folded populations of several peptides chosen to dissect specific noncovalent contributions to the conformational stability of proteins [alanine-based helical peptides stabilized by hydrogen bonds, tryptophan zipper (trpzip) peptides stabilized largely by cross-strand indole-indole interactions]. While the trpzip peptides are highly sensitive to GdmCl denaturation, they are unaffected by NaCl, Na2SO4, or Gdm2SO4, indicating that the reversal of the denaturant activity of Gdm+ by sulfate in this case is not due to competing stabilizing (sulfate) and destabilizing (Gdm+) interactions. Gdm2SO4 was found to retain considerable denaturant activity against alanine-based alpha-helical peptides. The differences in the effects of Gdm2SO4 on the two peptide types can be understood in terms of the different mechanisms of Gdm+ denaturation of trpzip peptides and helical peptides, respectively, and the specific nature of Gdm+ and SO42- ionic "clustering" that differentially affects the ability of Gdm+ to make the molecular interactions with the peptides that underlie its denaturant activity.

We describe the construction of a model complex of the cellobiohydrolase I (CBH I) cellulase from Trichoderma reesei bound to a cellulose microfibril in an aqueous environment for use in molecular dynamics (MD) simulations. Preliminary characterization from the initial phases of an MD simulation of this complex is also described. The linker sequence between the two globular domains was found to be quite flexible, and the oligosaccharides bound to this linker were found to prefer to be splayed like the spokes in a wheel due to their hydration requirements. The overall conformations of the two globular domains remained stable in the simulations, although both underwent changes in their orientations.

Molecular dynamics simulations of pure liquid water under ambient conditions using four common empirical water models have been analyzed to determine how well the oxygen-oxygen radial distribution function, g(OO)(r), used as the sole criterion of congruence with experiment, captures variations in the actual anisotropic collective structuring for these models. The largest systematic deviations from tetrahedrality were found to be due to deformations of the angle between the two closet hydrogen bond donor neighbors, but for intrinsic geometric reasons, these were found to contribute less to g(OO)(r) than deformations of the angles between one hydrogen bond donor neighbor and one hydrogen bond acceptor neighbor. Relying exclusively on a qualitative characterization of the second peak in g(OO)(r) seems to overemphasize the differences between the structuring in some of these models.

Lignocellulosic biomass has long been recognized as a potential sustainable source of mixed sugars for fermentation to biofuels and other biomaterials. Several technologies have been developed during the past 80 years that allow this conversion process to occur, and the clear objective now is to make this process cost-competitive in today's markets. Here, we consider the natural resistance of plant cell walls to microbial and enzymatic deconstruction, collectively known as “biomass recalcitrance.” It is this property of plants that is largely responsible for the high cost of lignocellulose conversion. To achieve sustainable energy production, it will be necessary to overcome the chemical and structural properties that have evolved in biomass to prevent its disassembly.

Molecular dynamics simulations of cellulose have recently become more prevalent due to increased interest in renewable energy applications, and many atomistic and coarse-grained force fields exist that can be applied to cellulose. However, to date no systematic comparison between carbohydrate force fields has been conducted for this important system. To that end, we present a molecular dynamics simulation study of hydrated, 36-chain cellulose Iβ microfibrils at room temperature with three carbohydrate force fields (CHARMM35, GLYCAM06, and Gromos 45a4) up to the near-microsecond time scale. Our results indicate that each of these simulated microfibrils diverge from the cellulose Iβ crystal structure to varying degrees under the conditions tested. The CHARMM35 and GLYCAM06 force fields eventually result in structures similar to those observed at 500 K with the same force fields, which are consistent with the experimentally observed high-temperature behavior of cellulose I. The third force field, Gromos 45a4, produces behavior significantly different from experiment, from the other two force fields, and from previously reported simulations with this force field using shorter simulation times and constrained periodic boundary conditions. For the GLYCAM06 force field, initial hydrogen-bond conformations and choice of electrostatic scaling factors significantly affect the rate of structural divergence. Our results suggest dramatically different time scales for convergence of properties of interest, which is important in the design of computational studies and comparisons to experimental data. This study highlights that further experimental and theoretical work is required to understand the structure of small diameter cellulose microfibrils typical of plant cellulose.

Three independently developed molecular mechanics force fields for carbohydrates have been used to simulate a suite of small molecule analogs of cellulose for which crystal structures have been reported, as a test to determine which might be best for simulations of cellulose itself. Such evaluation is necessary since the reported cellulose crystal structure is not stable in molecular dynamics simulations with any available force field. The present simulations found that all three resulted in small deviations from the reported crystal structures, but that all were reasonably accurate and none was clearly superior to the others for the entire suite of structures examined.

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.

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to study the structuring in aqueous solution of two cesium salts, cesium carbonate, and cesium nitrate. As was previously found for guanidinium salts of carbonate, mesoscopic-scale clusters were seen to form in the Cs2CO3 solution both in the MD simulations and in the diffraction experiments. No such large scale ion clusters were found in the CsNO3 solutions in either the modeling or experiments. The results are dominated by the strength and geometry of the direct first-neighbor interactions, which explain the differences in the clustering behavior between the two solutions without need to refer to longer-range water-water structuring.

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to examine the pentose D-xylose in aqueous solution. By specifically labeling D-xylose molecules with a deuterium atom at the nonexchangeable hydrogen position on C4, it was possible to extract information about the atomic structuring around just that specific position. The MD simulations were found to give satisfactory agreement with the experimental NDIS results and could be used to help interpret the scattering data in terms of the solvent structuring as well as the intramolecular hydroxyl conformations. Although the experiment is challenging and on the limit of modern instrumentation, it is possible by careful analysis, in conjunction with MD studies, to show that the conformation trans to H4 at 180 degrees is strongly disfavored, in excellent agreement with the MD results. This is the first attempt to use NDIS experiments to determine the rotameric conformation of a hydroxyl group.

The effects of the solution of glucose molecules on the dynamics of solvent water have been studied by quasielastic neutron scattering (QENS) measurements on solutions of selectively deuterated glucose in natural water. The data are fitted to two Lorentzians ascribed to pure translational and mixed translational and rotational character, respectively. The addition of the glucose to the water causes a substantial slowing down, by a factor 10 for the translational diffusion and 3−4 for the rotational motion at the highest concentration studied, 1:11 C6H12O6:H2O. The values obtained for water diffusion constants are consistent with previous QENS and NMR experiments on monosaccharide solutions but an order of magnitude higher than those derived from a recent molecular dynamics simulation.

A series of molecular mechanics calculations were used to analyze the energetics for moving a single polysaccharide chain from the surface of microcrystalline cellulose into the binding cleft of the Cel5A cellulase from Acidothermus cellulolyticus. A build-up procedure was used to model the placement of a 12-residue oligosaccharide chain along the surface of the enzyme, using as a guide the four residues of the tetrasaccharide substrate co-crystallized with the protein in the crystallographic structure determination. The position of this 12-residue oligosaccharide was used to orient the enzyme properly above two different surfaces of cellulose 1beta, the (1,0,0) and the (1,1,0) faces of the crystal. Constrained molecular dynamics simulations were then used to pull a target chain directly below the enzyme up out of the crystal surface and into the binding groove. The energetics for this process were favorable for both faces, with the step face being more favorable than the planar face, implying that this surface could be hydrolyzed more readily.

Neutron diffraction experiments were carried out at ambient temperature on aqueous solutions of glucose at concentrations of 1, 3, and 5 molal. The difference methods of neutron diffraction and isotopic substitution were used to determine various structure factors associated with all the exchangeable hydrogen atoms of the solutions. The results imply that the hydration of the glucose molecule is not concentration sensitive. It is also found that the average water structure itself is relatively unperturbed by the presence of glucose.

International audience; Mutagenesis experiments suggest that Asp79 in cellulase Cel6A (E2) from Thermobifida fusca has a catalytic role, in spite of the fact that this residue is more than 13 Å from the scissile bond in models of the enzyme–substrate complex built upon the crystal structure of the protein. This suggests that there is a substantial conformational shift in the protein upon substrate binding. Molecular mechanics simulations were used to investigate possible alternate conformations of the protein bound to a tetrasaccharide substrate, primarily involving shifts of the loop containing Asp79, and to model the role of water in the active site complex for both the native conformation and alternative low-energy conformations. Several alternative conformations of reasonable energy have been identified, including one in which the overall energy of the enzyme–substrate complex in solution is lower than that of the conformation in the crystal structure. This conformation was found to be stable in molecular dynamics simulations with a cellotetraose substrate and water. In simulations of the substrate complexed with the native protein conformation, the sugar ring in the –1 binding site was observed to make a spontaneous transition from the 4C1 conformation to a twist-boat conformer, consistent with generally accepted glycosidase mechanisms. Also, from these simulations Tyr73 and Arg78 were found to have important roles in the active site. Based on the results of these various MD simulations, a new catalytic mechanism is proposed. Using this mechanism, predictions about the effects of changes in Arg78 were made which were confirmed by site-directed mutagenesis.

The solution conformation of neocarrabiose 41-sulfate, which is a building block of κ-carrageenan, was investigated by a combination of molecular dynamics simulations and 1H NMR experiments. The calculated Ramachandran-type potential energy map for the glycosidic torsion angles φ and ψ of this dimer showed that there are two principal minimum-energy conformations, designated here as “a” and “b”. Comparison of the depths of these minima and the MD trajectories started from each minimum indicated that in a vacuum the molecular conformation in the “a” well is more stable than that in the “b” well. On the other hand, MD simulations in water showed that trajectories started from each minimum position stayed within their respective wells throughout the time courses of the simulations. However, analysis of the hydration water around the molecule indicated that the stability of the molecular conformation in the “a” well would be favored by solvation compared to the “b” well. Analysis of the experimental NOE−NMR d...

Molecular dynamics simulations employing umbrella sampling techniques have been used to calculate the potential of mean force (pmf) for the binding of a β-d-glucopyranose molecule to a benzene molecule in aqueous solution, as a model for the binding of sugar substrates to phenylalanine residues in proteins. The interaction of these two molecules was found to be strongly affected by hydration, as expected, with their nonpolar faces pairing by hydrophobic association to minimize the exposure of apolar groups to water. The pmf for the approach of these molecules is oscillatory in character, with two primary low-energy minima separated by a high free energy barrier, and with the net binding energy being approximately 1.1 kcal/mol and the intervening barrier being almost 2 kcal/mol. A third, weaker minimum was observed between 11 and 12 A. The oscillatory nature of the pmf results from the successive removal of water layers between the two molecules as they approach one another.

Molecular dynamics simulations employing adaptive umbrella sampling have been used to calculate the Ramachandran conformational potential of mean force in aqueous (TIP3P) solution for the α(1→4)-linked dimer of d-xylopyranose (4-O-α-d-xylopyranosyl-α-d-xylopyranose), a pentose analogue of maltose and a useful general model for the effects of solvent structuring upon biopolymer hydration. The vacuum adiabatic energy map for this molecule closely resembles that for maltose, but the solution pmf is quite different, with one of the principal vacuum minima almost completely disappearing in solution and with the global minimum-energy conformation being a new minimum which does not occur at all on the vacuum surface. This conformation is apparently stabilized by a water molecule which hydrogen bonds to a hydroxyl group on each ring, bridging between the two rings. The new conformation also places the two hydrophobic methylene groups almost in van der Waals contact, reducing their exposed surface area. Unfortunat...

Molecular mechanics calculations have been performed for the disaccharide neocarrabiose, one of the repeat units of beta-carrageenan, as a general model for the (1-->3)-linkage in the carrageenans. An adiabatic conformational energy map for this molecule has been prepared by constrained energy minimization and compared to previously reported relaxed maps. Neither the experimentally determined crystal structure of neocarrabiose nor the fiber diffraction conformation of iota-carrageenan is a low energy conformation on the relaxed Ramachandran map. Molecular dynamics simulations in vacuum produced trajectories consistent with this relaxed vacuum surface. However, a simulation with explicitly included solvent water molecules produced a trajectory that remained in the region of the two experimental structures. This dramatic solvation effect is apparently the result of the breaking of an interring hydrogen bond between the O2 hydroxyl groups of neocarrabiose as both groups hydrogen bond to solvent.

Recently, it was shown using structural neutron diffraction with isotopic substitutions (NDIS) measurements, combined with molecular dynamics simulations, that in an aqueous solution of D-xylose the hydroxyl group on the C4 position does not significantly occupy the position trans to the H4 atom. Here, a similar combination of NDIS and MD studies is described which uses D-xylose deuterated at the C5 position to further characterize this hydroxyl conformation as being trans to the C5 atom, as predicted by constrained MD simulations, confirming the previous study.

Limitations in experimental measurements complicate the investigations of biologically important disaccharides. Previous computational studies of disaccharides have been very CPU intensive. The application of simulated annealing to conformational space studies of disaccharides is shown to decrease the computational effort by more than an order of magnitude. At the same time this approach produces a more accurate description of the conformational space particularly in the high energy regions. The method is demonstrated on the biologically important GlcNac-β-(1-4)-GlcNac disaccharide and can be easily translated for application to other carbohydrates or related molecules. The energy surface E ( Π , ψ ) for the GlcNac-β-(1-4)-GlcNac disaccharide which is a function of the glycosidic linkage is analysed and compared to a similar surface generated by the standard computationally more intense exhaustive search method.

Molecular dynamics simulations have been used to model the flexibility of the seven-sugar oligosaccharide of the lectin from Erythrina corallodendron in three separate simulations: one of the isolated oligosaccharide in vacuo, one of the oligosaccharide in solution and one of the oligosaccharide linked to the protein in solution. Adiabatic conformational energy maps were prepared for each of the disaccharide linkages as a means of interpreting the observed dynamics and conformational averages in terms of intramolecular energy. The inclusion of aqueous solvent molecules appears to be necessary to reproduce the experimental conformational behavior, which also cannot be predicted well from conformational energy maps for the disaccharide linkages alone. The crystallographically determined conformation does not appear to be induced by the crystal dimerization, but is rather stable in solution. The build-up of fluctuations along the successive linkages of the oligosaccharide is significant and would be sufficient to prevent branch residues from being located in most crystal structure determinations. Good general agreement between the calculated solution structure and the average structure determined by NMR was found for most of the oligosaccharide linkages.

The structure of the lectin from Erythrina corallodendron (EcorL) has recently been solved crystallographically, with the oligosaccharide portion of this glycoprotein unusually well resolved. In general, the structure and dynamics of the oligosaccharides in glycoproteins are poorly understood. The crystal structure of EcorL provides an opportunity to probe the conformational properties of the oligosaccharide using computer modeling. Two 300 ps molecular dynamics simulations are reported of the oligosaccharide in solution isolated from the protein and bound to the protein and surrounded by water. The motions of the oligosaccharide are similar both free in solution and bonded to the protein. However, different parts of the glycosidic linkage conformational space have been explored by key component disaccharides in the two environments. The disaccharide of the “core” region shows greater flexibility in the isolated oligosaccharide than in the same region in the bound oligosaccharide. Because of the short length of these simulations, it could not be determined whether these differences are due to random chance or are a true reflection of the influence of the protein on the conformation of the oligosaccharide.