Blind prediction of solvation free energies from the SAMPL4 challenge

Here, we give an overview of the small molecule hydration portion of the SAMPL4 challenge, which focused on predicting hydration free energies for a series of 47 small molecules. These gas-to-water transfer free energies have in the past proven a valuable test of a variety of computational methods and force fields. Here, in contrast to some previous SAMPL challenges, we find a relatively wide range of methods perform quite well on this test set, with RMS errors in the 1.2 kcal/mol range for several of the best performing methods. Top-performers included a quantum mechanical approach with continuum solvent models and functional group corrections, alchemical molecular dynamics simulations with a classical all-atom force field, and a single-conformation Poisson–Boltzmann approach. While 1.2 kcal/mol is still a significant error, experimental hydration free energies covered a range of nearly 20 kcal/mol, so methods typically showed substantial predictive power. Here, a substantial new focus was on evaluation of error estimates, as predicting when a computational prediction is reliable versus unreliable has considerable practical value. We found, however, that in many cases errors are substantially underestimated, and that typically little effort has been invested in estimating likely error. We believe this is an important area for further research.     

Overview of the SAMPL5 host–guest challenge: Are we doing better?

The ability to computationally predict protein-small molecule binding affinities with high accuracy would accelerate drug discovery and reduce its cost by eliminating rounds of trial-and-error synthesis and experimental evaluation of candidate ligands. As academic and industrial groups work toward this capability, there is an ongoing need for datasets that can be used to rigorously test new computational methods. Although protein–ligand data are clearly important for this purpose, their size and complexity make it difficult to obtain well-converged results and to troubleshoot computational methods. Host–guest systems offer a valuable alternative class of test cases, as they exemplify noncovalent molecular recognition but are far smaller and simpler. As a consequence, host–guest systems have been part of the prior two rounds of SAMPL prediction exercises, and they also figure in the present SAMPL5 round. In addition to being blinded, and thus avoiding biases that may arise in retrospective studies, the SAMPL challenges have the merit of focusing multiple researchers on a common set of molecular systems, so that methods may be compared and ideas exchanged. The present paper provides an overview of the host–guest component of SAMPL5, which centers on three different hosts, two octa-acids and a glycoluril-based molecular clip, and two different sets of guest molecules, in aqueous solution. A range of methods were applied, including electronic structure calculations with implicit solvent models; methods that combine empirical force fields with implicit solvent models; and explicit solvent free energy simulations. The most reliable methods tend to fall in the latter class, consistent with results in prior SAMPL rounds, but the level of accuracy is still below that sought for reliable computer-aided drug design. Advances in force field accuracy, modeling of protonation equilibria, electronic structure methods, and solvent models, hold promise for future improvements.     

Comparison of charge models for fixed-charge force fields: small-molecule hydration free energies in explicit solvent

In molecular simulations with fixed-charge force fields, the choice of partial atomic charges influences numerous computed physical properties, including binding free energies. Many molecular mechanics force fields specify how nonbonded parameters should be determined, but various choices are often available for how these charges are to be determined for arbitrary small molecules. Here, we compute hydration free energies for a set of 44 small, neutral molecules in two different explicit water models (TIP3P and TIP4P-Ew) to examine the influence of charge model on agreement with experiment. Using the AMBER GAFF force field for nonbonded parameters, we test several different methods for obtaining partial atomic charges, including two fast methods exploiting semiempirical quantum calculations and methods deriving charges from the electrostatic potentials computed with several different levels of ab initio quantum calculations with and without a continuum reaction field treatment of solvent. We find that the best charge sets give a root-mean-square error from experiment of roughly 1 kcal/mol. Surprisingly, agreement with experimental hydration free energies does not increase substantially with increasing level of quantum theory, even when the quantum calculations are performed with a reaction field treatment to better model the aqueous phase. We also find that the semiempirical AM1-BCC method for computing charges works almost as well as any of the more computationally expensive ab initio methods and that the root-mean-square error reported here is similar to that for implicit solvent models reported in the literature. Further, we find that the discrepancy with experimental hydration free energies grows substantially with the polarity of the compound, as does its variation across theory levels.     

The Confine-and-Release Method: Obtaining Correct Binding Free Energies in the Presence of Protein Conformational Change

Free energy calculations are increasingly being used to estimate absolute and relative binding free energies of ligands to proteins. However, computed free energies often appear to depend on the initial protein conformation, indicating incomplete sampling. This is especially true when proteins can change conformation on ligand binding, as free energies associated with these conformational changes are either ignored or assumed to be included by virtue of the sampling performed in the calculation. Here, we show that, in a model protein system (a designed binding site in T4 Lysozyme), conformational changes can make a difference of several kcal/mol in computed binding free energies, and that they are neglected in computed binding free energies if the system remains kinetically trapped in a particular metastable state on simulation timescales. We introduce a general "confine-and-release" framework for free energy calculations that accounts for these free energies of conformational change. We illustrate its use in this model system by demonstrating that an umbrella sampling protocol can obtain converged binding free energies that are independent of the starting protein structure and include these conformational change free energies.     

Alchemical prediction of hydration free energies for SAMPL

Hydration free energy calculations have become important tests of force fields. Alchemical free energy calculations based on molecular dynamics simulations provide a rigorous way to calculate these free energies for a particular force field, given sufficient sampling. Here, we report results of alchemical hydration free energy calculations for the set of small molecules comprising the 2011 Statistical Assessment of Modeling of Proteins and Ligands challenge. Our calculations are largely based on the Generalized Amber Force Field with several different charge models, and we achieved RMS errors in the 1.4-2.2 kcal/mol range depending on charge model, marginally higher than what we typically observed in previous studies (Mobley et al. in J Phys Chem B 111(9):2242-2254, 2007, J Chem Theory Comput 5(2):350-358, 2009, J Phys Chem B 115:1329-1332, 2011; Nicholls et al. in J Med Chem 51:769-779, 2008; Klimovich and Mobley in J Comput Aided Mol Design 24(4):307-316, 2010). The test set consists of ethane, biphenyl, and a dibenzyl dioxin, as well as a series of chlorinated derivatives of each. We found that, for this set, using high-quality partial charges from MP2/cc-PVTZ SCRF RESP fits provided marginally improved agreement with experiment over using AM1-BCC partial charges as we have more typically done, in keeping with our recent findings (Mobley et al. in J Phys Chem B 115:1329-1332, 2011). Switching to OPLS Lennard-Jones parameters with AM1-BCC charges also improves agreement with experiment. We also find a number of chemical trends within each molecular series which we can explain, but there are also some surprises, including some that are captured by the calculations and some that are not.     

The effect of static pressure on the inertial cavitation threshold

The amplitude of the acoustic pressure required to nucleate a gas or vapor bubble in a fluid, and to have that bubble undergo an inertial collapse, is termed the inertial cavitation threshold. The magnitude of the inertial cavitation threshold is typically limited by mechanisms other than homogeneous nucleation such that the theoretical maximum is never achieved. However, the onset of inertial cavitation can be suppressed by increasing the static pressure of the fluid. The inertial cavitation threshold was measured in ultrapure water at static pressures up to 30?MPa (300 bars) by exciting a radially symmetric standing wave field in a spherical resonator driven at a resonant frequency of 25.5 kHz. The threshold was found to increase linearly with the static pressure; an exponentially decaying temperature dependence was also found. The nature and properties of the nucleating mechanisms were investigated by comparing the measured thresholds to an independent analysis of the particulate content and available models for nucleation.     

Integrated circuit microchip system with multiplex capillary electrophoresis module for DNA analysis

In this paper, we describe the use of an integrated circuit (IC) microchip system as a detector in multiplex capillary electrophoresis (CE). This combination of multiplex capillary gel electrophoresis and the IC microchip technology represents a novel approach to DNA analysis on the microchip platform. Separation of DNA ladders using a multiplex CE microsystem of four capillaries was monitored simultaneously using the IC microchip system. The IC microchip-CE system has advantages such as low cost, rapid analysis, compactness, and multiplex capability, and has great potential as an alternative system to conventional capillary array gel electrophoresis systems based on charge-coupled device (CCD) detection.     

Intensified biochip system using chemiluminescence for the detection of Bacillus globigii spores

This paper reports the first intensified biochip system for chemiluminescence detection and the feasibility of using this system for the analysis of biological warfare agents is demonstrated. An enzyme-linked immunosorbent assay targeting Bacillus globigii spores, a surrogate species for Bacillus anthracis, using a chemiluminescent alkaline phosphatase substrate is combined with a compact intensified biochip detection system. The enzymatic amplification was found to be an attractive method for detection of low spore concentrations when combined with the intensified biochip device. This system was capable of detecting approximately 1 × 10⁵ Bacillus globigii spores. Moreover, the chemiluminescence method, combined with the self-contained biochip design, allows for a simple, compact system that does not require laser excitation and is readily adaptable to field use. Figure Schematic diagram of the miniature biochip detection system     

Review of Chemometrics Applied to Spectroscopy: 1985-95, Part 2

INTRODUCTION

Before beginning Part 2 of this review, a caveat noted by Deming and Palasota is brought to the reader's attention: [1] “Press et al. [2] have emphasized that data ‘consist of numbers, of course. But these numbers are fed into the computer, not produced by it. These are numbers to be treated with considerable respect, never to be tampered with, nor subjected to a numerical process whose character you do not completely understand. You are well advised to acquire reverence for data that is rather diferent fiom the “sporty” attitude which is sometimes allowable, or even commendable, in other numerical tasks.’ Yet by and large within chemometrics, preprocessing often seems to be carried out with little understanding of its fundamental efect on the structure of the data.”