Molecular dynamics simulations of pea (Pisum sativum) lectin structure with octyl glucoside detergents: the ligand interactions and dynamics

Binding structures of tri-N-acetyl-β-glucosamine in hen egg white lysozyme using molecular dynamics with a polarizable force field

Lysozyme is a well-studied enzyme that hydrolyzes the β-(1,4)-glycosidic linkage of N-acetyl-β-glucosamine (NAG)(n) oligomers. The active site of hen egg-white lysozyme (HEWL) is believed to consist of six subsites, A-F that can accommodate six sugar residues. We present studies exploring the use of polarizable force fields in conjunction with all-atom molecular dynamics (MD) simulations to analyze binding structures of complexes of lysozyme and NAG trisaccharide, (NAG)(3). MD trajectories are applied to analyze structures and conformation of the complex as well as protein-ligand interactions, including the hydrogen-bonding network in the binding pocket. Two binding modes (ABC and BCD) of (NAG)(3) are investigated independently based on a fixed-charge model and a polarizable model. We also apply molecular mechanics with generalized born and surface area (MM-GBSA) methods based on MD using both nonpolarizable and polarizable force fields to compute binding free energies. We also study the correlation between root-mean-squared deviation and binding free energies of the wildtype and W62Y mutant; we find that for this prototypical system, approaches using the MD trajectories coupled with implicit solvent models are equivalent for polarizable and fixed-charge models.

A molecular simulation based assessment of binding of metal ions on micelles

An assessment of the ability of a micellar surface to bind different metal ions using molecular simulation is presented in this study. Sodium dodecyl sulfate (SDS) is considered as the anionic surfactant. Various relevant characteristics of SDS-metal ion systems are estimated to quantify preferential binding of metal ions. These are electrostatic energy, total potential energy of the system, radial distribution function, and entropy and free energy change of the system. By examining these parameters, the relative extents of binding of different metal ions to the micellar surface are assessed.

Molecular simulations of dodecyl-β-maltoside micelles in water: influence of the headgroup conformation and force field parameters

This paper deals with the development and validation of new potential parameter sets, based on the CHARMM36 and GLYCAM06 force fields, to simulate micelles of the two anomeric forms (α and β) of N-dodecyl-β-maltoside (C(12)G(2)), a surfactant widely used in the extraction and purification of membrane proteins. In this context, properties such as size, shape, internal structure, and hydration of the C(12)G(2) anomer micelles were thoroughly investigated by molecular dynamics simulations and the results compared with experiments. Additional simulations were also performed with the older CHARMM22 force field for carbohydrates (Kuttel, M.; et al. J. Comput. Chem. 2002, 23, 1236-1243). We find that our CHARMM and GLYCAM parameter sets yield similar results in the case of properties related to the micelle structure but differ for other properties such as the headgroup conformation or the micelle hydration. In agreement with experiments, our results show that for all model potentials the β-C(12)G(2) micelles have a more pronounced ellipsoidal shape than those containing α anomers. The computed radius of gyration is 20.2 and 25.4 Å for the α- and β-anomer micelles, respectively. Finally, we show that depending on the potential the water translational diffusion of the interfacial water is 7-11.5 times slower than that of bulk water due to the entrapment of the water in the micelle crevices. This retardation is independent of the headgroup in α- or β-anomers.

Structure, dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X-ray studies

Molecular dynamics simulations have been carried out on all the jacalin-carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X-ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin-carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate-binding residues are consistent with the known thermodynamic parameters of jacalin-carbohydrate interactions. The simulations, along with X-ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin-carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit.

High processing tolerances of immunomodulatory proteins in Enoki and Reishi mushrooms

This study investigated the processing tolerances of two mushroom proteins with immunomodulatory activities, including FVE from Enoki ( Flammulia velutipes ) and LZ8 from Reishi ( Ganoderma lucidum ) mushrooms, under food processing treatments such as heating, sterilization, frozen storage, extraction in acid/alkaline conditions, and dehydration. Results showed that the capability of these two proteins to induce IFN-gamma secretion by murine splenocytes remained after 100 degrees C heating for 30 min, 121 degrees C autoclaving for 15 min, and -80 degrees C freezing. The retained activities of both proteins on cell proliferation and IFN-gamma production did not decrease at 0.6 M hydrochloric acid (at pH 2) but strikingly dropped at 5 M sodium hydrate (at pH 13). After vacuum dehydration, FVE and LZ8 retained most of their activities on cell proliferation; nevertheless, the IFN-gamma secretion decreased to about half of the initial values. These findings suggest that these two mushroom proteins have a good thermal/freezing resistance, acid tolerance, and dehydration stability and are candidates for processing in food and nutraceutical utilization.