Eric Pellegrini

Vibrational Density of States of Hydration Water at Biomolecular Sites: Hydrophobicity Promotes Low Density Amorphous Ice Behavior

Inelastic neutron scattering experiments and molecular dynamics simulations have been used to investigate the low frequency modes, in the region between 0 and 100 meV, of hydration water in selected hydrophilic and hydrophobic biomolecules. The results show changes in the plasticity of the hydrogen-bond network of hydration water molecules depending on the biomolecular site. At 200 K, the measured low frequency density of states of hydration water molecules of hydrophilic peptides is remarkably similar to that of high density amorphous ice, whereas, for hydrophobic biomolecules, it is comparable to that of low density amorphous ice behavior. In both hydrophilic and hydrophobic biomolecules, the high frequency modes show a blue shift of the libration mode as compared to the room temperature data. These results can be related to the density of water molecules around the biological interface, suggesting that the apparent local density of water is larger in a hydrophilic environment.

Evidence of dynamical constraints imposed by water organization around a bio-hydrophobic interface.

Molecular dynamics simulations and elastic neutron scattering experiments have been used to highlight how the structural organization of hydration water is able in some cases to locally constrain atomic movements at biologic interfaces. Using fully hydrated small peptides as models of protein interfaces, we show that the length of the side chains and the hydrophilic backbone have specific signatures. The dynamics of the side chain, which is part of biomolecules, have not only a crucial role in the whole flexibility as compared to the backbone, but also modify the values of transition temperatures. The analysis of the activation energies of methyl group dynamics suggests that the interaction between hydrophobic side chain and surrounding water plays an important role in the whole flexibility as well. We suggest that the progressive water cluster organization, around hydrophobic interfaces increases the activation energy and that a plateau regime is reached only when an extended hydrogen-bond network is established. The cluster size corresponds to a single layer of water molecules.

Development and testing of a de novo drug-design algorithm.

In this article we present an implementation of a de novo drug-design algorithm. The algorithm starts with a molecule placed in the binding site of a protein and then modifies it using a sequential growth approach. This involves successive cycles of suppression of randomly picked groups in the molecule and their replacement by other groups chosen from databanks of linear or cyclic fragments. The algorithm has been coupled with the Dynamo library which allows the simulation of macromolecules using molecular mechanical and quantum chemical methods. The main body of the article describes the methodologies we use to create, characterize and evaluate putative ligands. We also consider briefly an application of the algorithm to a protein of pharmacological interest, the neuraminidase of the influenza virus, and discuss the strengths and weaknesses of our approach.

Direct comparison of elastic incoherent neutron scattering experiments with molecular dynamics simulations of DMPC phase transitions.

Abstract.Neutron scattering techniques have been employed to investigate 1,2-dimyristoyl-sn -glycero-3-phosphocholine (DMPC) membranes in the form of multilamellar vesicles (MLVs) and deposited, stacked multilamellar-bilayers (MLBs), covering transitions from the gel to the liquid phase. Neutron diffraction was used to characterise the samples in terms of transition temperatures, whereas elastic incoherent neutron scattering (EINS) demonstrates that the dynamics on the sub-macromolecular length-scale and pico- to nano-second time-scale are correlated with the structural transitions through a discontinuity in the observed elastic intensities and the derived mean square displacements. Molecular dynamics simulations have been performed in parallel focussing on the length-, time- and temperature-scales of the neutron experiments. They correctly reproduce the structural features of the main gel-liquid phase transition. Particular emphasis is placed on the dynamical amplitudes derived from experiment and simulations. Two methods are used to analyse the experimental data and mean square displacements. They agree within a factor of 2 irrespective of the probed time-scale, i.e. the instrument utilized. Mean square displacements computed from simulations show a comparable level of agreement with the experimental values, albeit, the best match with the two methods varies for the two instruments. Consequently, experiments and simulations together give a consistent picture of the structural and dynamical aspects of the main lipid transition and provide a basis for future, theoretical modelling of dynamics and phase behaviour in membranes. The need for more detailed analytical models is pointed out by the remaining variation of the dynamical amplitudes derived in two different ways from experiments on the one hand and simulations on the other.Graphical abstract

Perspectives in biological physics: the nDDB project for a neutron Dynamics Data Bank for biological macromolecules.

Neutron spectroscopy provides experimental data on time-dependent trajectories, which can be directly compared to molecular dynamics simulations. Its importance in helping us to understand biological macromolecules at a molecular level is demonstrated by the results of a literature survey over the last two to three decades. Around 300 articles in refereed journals relate to neutron scattering studies of biological macromolecular dynamics, and the results of the survey are presented here. The scope of the publications ranges from the general physics of protein and solvent dynamics, to the biologically relevant dynamics-function relationships in live cells. As a result of the survey we are currently setting up a neutron Dynamics Data Bank (n DDB) with the aim to make the neutron data on biological systems widely available. This will benefit, in particular, the MD simulation community to validate and improve their force fields. The aim of the database is to expose and give easy access to a body of experimental data to the scientific community. The database will be populated with as much of the existing data as possible. In the future it will give value, as part of a bigger whole, to high throughput data, as well as more detailed studies. A range and volume of experimental data will be of interest in determining how quantitatively MD simulations can reproduce trends across a range of systems and to what extent such trends may depend on sample preparation and data reduction and analysis methods. In this context, we strongly encourage researchers in the field to deposit their data in the n DDB.Graphical abstract

Computational Tools to Understand Inelastic and Quasielastic Neutron Scattering Data

The availability of user-friendly software and affordable multi-processor computers opens the door to the world of simulations to experimentalists, for “advanced data analysis”. Neutron scattering (NS), which explores length and time scales and probes the relative positions and motions of atoms as in simulations, constitutes the ideal partner for atomistic simulations. On the experimental side, the ever-increasing complexity of samples and therefore data requires more elaborate and realistic models. This chapter therefore describes, in practical terms, the simulation methods that can be used to interpret quasielastic and inelastic NS data, namely molecular dynamics and lattice dynamics. Both of these methods are based on the knowledge of inter-atomic interactions and total energy for which density functional theory and classical, force field-based methods are presented as the most viable.