Definition and testing of the GROMOS force-field versions 54A7 and 54B7

New parameter sets of the GROMOS biomolecular force field, 54A7 and 54B7, are introduced. These parameter sets summarise some previously published force field modifications: The 53A6 helical propensities are corrected through new φ/ψ torsional angle terms and a modification of the N-H, C=O repulsion, a new atom type for a charged −CH3 in the choline moiety is added, the Na+ and Cl− ions are modified to reproduce the free energy of hydration, and additional improper torsional angle types for free energy calculations involving a chirality change are introduced. The new helical propensity modification is tested using the benchmark proteins hen egg-white lysozyme, fox1 RNA binding domain, chorismate mutase and the GCN4-p1 peptide. The stability of the proteins is improved in comparison with the 53A6 force field, and good agreement with a range of primary experimental data is obtaine

Exposing the flexibility of human parainfluenza virus hemagglutinin- neuraminidase

Human parainfluenza virus type 3 (hPIV-3) is a clinically significant pathogen and is the causative agent of pneumonia and bronchiolitis in children. In this study the solution dynamics of human parainfluenza type 3 hemagglutinin-neuraminidase (HN) have been investigated. A flexible loop around Asp216 that adopts an open conformation in direct vicinity of the active site of the apo-form of the protein and closes upon inhibitor binding has been identified. To date, no available X-ray crystal structure has shown the molecular dynamics simulation-derived predominant loop-conformation states found in the present study. The outcomes of this study provide additional insight into the dynamical properties of hPIV-3 HN and may have important implications in defining HN glycan recognition events, receptor specificity, and antiparainfluenza virus drug discovery

Exposing the Flexibility of Human Parainfluenza Virus Hemagglutinin-neuraminidase

Human parainfluenza virus type 3 (hPIV-3) is a clinically significant pathogen and is the causative agent of pneumonia and bronchiolitis in children. In this study the solution dynamics of human parainfluenza type 3 hemagglutinin-neuraminidase (HN) have been investigated. A flexible loop around Asp216 that adopts an open conformation in direct vicinity of the active site of the <i>apo</i>-form of the protein and closes upon inhibitor binding has been identified. To date, no available X-ray crystal structure has shown the molecular dynamics simulation-derived predominant loop-conformation states found in the present study. The outcomes of this study provide additional insight into the dynamical properties of hPIV-3 HN and may have important implications in defining HN glycan recognition events, receptor specificity, and antiparainfluenza virus drug discovery

On the conformational properties of amylose and cellulose oligomers in solution

Molecular dynamics (MD) simulations were used to monitor the stability and conformation of double-stranded and single-stranded amyloses and single-stranded cellulose oligomers containing 9 sugar moieties in solution as a function of solvent composition, ionic strength, temperature, and methylation state. This study along with other previous studies suggests that hydrogen bonds are crucial for guaranteeing the stability of the amylose double helix. Single-stranded amylose forms a helical structure as well, and cellulose stays highly elongated throughout the simulation time, a behavior that was also observed experimentally. In terms of coordination of solute hydroxyl groups with ions, amylose shows entropy-driven coordination of calcium and sulfate ions, whereas cellulose-ion coordination seems to be enthalpy-dominated. This indicates that entropy considerations cannot be neglected when explaining the structural differences between amyloses and cellulose

Force-field dependence of the conformational properties of a,x-dimethoxypolyethylene glycol

A molecular dynamics (MD) study of ,-dimethoxypolyethylene glycol has been carried out under various conditions with respect to solvent composition, ionic strength, chain length, force field and temperature. A previous MD study on a 15-mer of polyethyleneglycol (PEG) suggested a helical equilibrium structure that was stabilised by hydrogen bonding and bridging water molecules. Experiments show that PEG is highly soluble in water, and indicate that clustering is not favoured. In the present study using different force fields, the GROMOS force fields 45A3 and 53A6, a variation on the latter 53A6_OE, and a force field by Smith et al. produced different results. For the GROMOS force fields 45A3 and 53A6 no helical structure was found, but formation of more or less compact random coils in aqueous solution due to hydrophobic interactions was observed. For the other two force fields used, ,-dimethoxypolyethylene glycol stayed flexible and more or less elongated in aqueous solution, more in agreement with experimental observations and the previous MD study

On using a too large integration time step in molecular dynamics simulations of coarse-grained molecular models

The use of a coarse-grained (CG) model that is widely used in molecular dynamics simulations of biomolecular systems is investigated with respect to the dependence of a variety of quantities upon the size of the used integration time step and cutoff radius. The results suggest that when using a non-bonded interaction-cutoff radius of 1.4 nm a time step of maximally 10 fs should be used, in order not to produce energy sinks or wells. Using a too-large time step, e.g. 50 fs with a cutoff of 1.2 nm, as is done in the coarse-grained model of Marrink et al. (J. Phys. Chem. B, 2004, 108, 250 and 2007, 111, 7812), induces errors due to the linear approximation of the integrators that are commonly used to integrate the equations of motion. As a spin-off of the investigation of the mentioned CG models, we found that the parameters of the CG water model place it at physiological temperatures well into the solid phase of the phase diagram. Copyright the Owner Societies 200

Definition and testing of the GROMOS force-field versions 54A7 and 54B7

New parameter sets of the GROMOS biomolecular force field, 54A7 and 54B7, are introduced. These parameter sets summarise some previously published force field modifications: The 53A6 helical propensities are corrected through new phi/psi torsional angle terms and a modification of the N-H, C=O repulsion, a new atom type for a charged -CH3 in the choline moiety is added, the Na+ and Cl- ions are modified to reproduce the free energy of hydration, and additional improper torsional angle types for free energy calculations involving a chirality change are introduced. The new helical propensity modification is tested using the benchmark proteins hen egg-white lysozyme, fox1 RNA binding domain, chorismate mutase and the GCN4-p1 peptide. The stability of the proteins is improved in comparison with the 53A6 force field, and good agreement with a range of primary experimental data is obtained

Investigation into the Feasibility of Thioditaloside as a Novel Scaffold for Galectin-3-Specific Inhibitors

Galectin-3 is extensively involved in metabolic and disease processes, such as cancer metastasis, thus giving impetus for the design of specific inhibitors targeting this -galactose-binding protein. Thiodigalactoside (TDG) presents a scaffold for construction of galectin inhibitors, and its inhibition of galectin-1 has already demonstrated beneficial effects as an adjuvant with vaccine immunotherapy, thereby improving the survival outcome of tumour-challenged mice. A novel approachreplacing galactose with its C2 epimer, taloseoffers an alternative framework, as extensions at C2 permit exploitation of a galectin-3-specific binding groove, thereby facilitating the design of selective inhibitors. We report the synthesis of thioditaloside (TDT) and crystal structures of the galectin-3 carbohydrate recognition domain in complexes with TDT and TDG. The different abilities of galactose and talose to anchor to the protein correlate with molecular dynamics studies, likely explaining the relative disaccharide binding affinities. The feasibility of a TDT scaffold to enable access to a particular galectin-3 binding groove and the need for modifications to optimise such a scaffold for use in the design of potent and selective inhibitors are assessed