Relationship between Structure, Entropy and Diffusivity in Water and Water-like Liquids

Anomalous behaviour of the excess entropy ($S_e$) and the associated scaling relationship with diffusivity are compared in liquids with very different underlying interactions but similar water-like anomalies: water (SPC/E and TIP3P models), tetrahedral ionic melts (SiO$_2$ and BeF$_2$) and a fluid with core-softened, two-scale ramp (2SRP) interactions. We demonstrate the presence of an excess entropy anomaly in the two water models. Using length and energy scales appropriate for onset of anomalous behaviour, the density range of the excess entropy anomaly is shown to be much narrower in water than in ionic melts or the 2SRP fluid. While the reduced diffusivities ($D^*$) conform to the excess entropy scaling relation, $D^* =A\exp (\alpha S_e)$ for all the systems (Y. Rosenfeld, Phys. Rev. A {\bf 1977}, {\it 15}, 2545), the exponential scaling parameter, $\alpha$, shows a small isochore-dependence in the case of water. Replacing $S_e$ by pair correlation-based approximants accentuates the isochore-dependence of the diffusivity scaling. Isochores with similar diffusivity scaling parameters are shown to have the temperature dependence of the corresponding entropic contribution. The relationship between diffusivity, excess entropy and pair correlation approximants to the excess entropy are very similar in all the tetrahedral liquids.Comment: 24 pages, 4 figures, to be published in Journal of Physical Chemistry

Charting Evolution’s Trajectory: Using Molluscan Eye Diversity to Understand Parallel and Convergent Evolution

For over 100 years, molluscan eyes have been used as an example of convergent evolution and, more recently, as a textbook example of stepwise evolution of a complex lens eye via natural selection. Yet, little is known about the underlying mechanisms that create the eye and generate different morphologies. Assessing molluscan eye diversity and understanding how this diversity came about will be important to developing meaningful interpretations of evolutionary processes. This paper provides an introduction to the myriad of eye types found in molluscs, focusing on some of the more unusual structures. We discuss how molluscan eyes can be applied to the study of evolution by examining patterns of convergent and parallel evolution and provide several examples, including the putative convergence of the camera-type eyes of cephalopods and vertebrates

Rapid and Accurate Prediction and Scoring of Water Molecules in Protein Binding Sites

Water plays a critical role in ligand-protein interactions. However, it is still challenging to predict accurately not only where water molecules prefer to bind, but also which of those water molecules might be displaceable. The latter is often seen as a route to optimizing affinity of potential drug candidates. Using a protocol we call WaterDock, we show that the freely available AutoDock Vina tool can be used to predict accurately the binding sites of water molecules. WaterDock was validated using data from X-ray crystallography, neutron diffraction and molecular dynamics simulations and correctly predicted 97% of the water molecules in the test set. In addition, we combined data-mining, heuristic and machine learning techniques to develop probabilistic water molecule classifiers. When applied to WaterDock predictions in the Astex Diverse Set of protein ligand complexes, we could identify whether a water molecule was conserved or displaced to an accuracy of 75%. A second model predicted whether water molecules were displaced by polar groups or by non-polar groups to an accuracy of 80%. These results should prove useful for anyone wishing to undertake rational design of new compounds where the displacement of water molecules is being considered as a route to improved affinity

The Dynamics of Ligand Barrier Crossing inside the Acetylcholinesterase Gorge

The dynamics of ligand movement through the constricted region of the acetylcholinesterase gorge is important in understanding how the ligand gains access to and is released from the active site of the enzyme. Molecular dynamics simulations of the simple ligand, tetramethylammonium, crossing this bottleneck region are conducted using umbrella potential sampling and activated flux techniques. The low potential of mean force obtained is consistent with the fast reaction rate of acetylcholinesterase observed experimentally. From the results of the activated dynamics simulations, local conformational fluctuations of the gorge residues and larger scale collective motions of the protein are found to correlate highly with the ligand crossing

Transition-state theory model for the diffusion coefficients of small penetrants in glassy polymers

Previous molecular dynamics simulations have shown that the diffusion of a penetrant in a glassy polymer involves occasional jumps between cavities through the opening of a "neck", and thus, because this is a rare event, the diffusion coefficient can be estimated using transition-state theory. We treat this process as a unimolecular rearrangement and develop semiempirical means of estimating the activation energy, frequency factor and jump length. The activation energy is obtained by treating the polymer as a continuous solid and calculating the energy required to expand a neck in that continuum. The model for the frequency factor uses the result from simulations that the distribution of frequencies of the modes in the transition is very similar to that distribution for the reactant state. The frequency factor is estimated by considering only the motion of the penetrant. These motions are treated as harmonic oscillators. The jump length is obtained from simple geometric considerations of the polymer chain. The ibility. The model reproduces experimental trends semiquantitatively and could be used to interpolate and extrapolate experimental diffusion data

Revisiting Free Energy Calculations: A Theoretical Connection to MM/PBSA and Direct Calculation of the Association Free Energy

The prediction of absolute ligand-receptor binding affinities is essential in a wide range of biophysical queries, from the study of protein-protein interactions to structure-based drug design. End-point free energy methods, such as the Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) model, have received much attention and widespread application in recent literature. These methods benefit from computational efficiency as only the initial and final states of the system are evaluated, yet there remains a need for strengthening their theoretical foundation. Here a clear connection between statistical thermodynamics and end-point free energy models is presented. The importance of the association free energy, arising from one molecule's loss of translational and rotational freedom from the standard state concentration, is addressed. A novel method for calculating this quantity directly from a molecular dynamics simulation is described. The challenges of accounting for changes in the protein conformation and its fluctuations from separate simulations are discussed. A simple first-order approximation of the configuration integral is presented to lay the groundwork for future efforts. This model has been applied to FKBP12, a small immunophilin that has been widely studied in the drug industry for its potential immunosuppressive and neuroregenerative effects