Leveraging the Thermodynamics of Protein Conformations in Drug Discovery.

As the name implies, structure-based drug design requires confidence in the holo complex structure. The ability to clarify which protein conformation to use when ambiguity arises would be incredibly useful. We present a large scale validation of the computational method Protein Reorganization Free Energy Perturbation (PReorg-FEP) and demonstrate its quantitative accuracy in selecting the correct protein conformation among candidate models in apo or ligand induced states for 14 different systems.

Is the Functional Response of a Receptor Determined by the Thermodynamics of Ligand Binding?

For an effective drug, strong binding to the target protein is a prerequisite, but it is not enough. To produce a particular functional response, drugs need to either block the proteins' functions or modulate their activities by changing their conformational equilibrium. The binding free energy of a compound to its target is routinely calculated, but the timescales for the protein conformational changes are prohibitively long to be efficiently modeled via physics-based simulations.

Protein Loop Conformational Free Energy Changes via an Alchemical Path without Reaction Coordinates.

We introduce a method called restrain-free energy perturbation-release 2.0 (R-FEP-R 2.0) to estimate conformational free energy changes of protein loops via an alchemical path.

Solvation Thermodynamics from the Perspective of Endpoints DFT.

Solvation thermodynamics is concerned with the evaluation and physical interpretation of solvation free energies. Endpoints DFT provides a framework for computing solvation free energies by combining molecular simulations with a version of the classical density-functional theory of solutions which focuses on ω, the indirect (solvent-mediated) part of the solute-solvent potential of mean force (indirect PMF). The simulations are performed at the endpoints of a hypothetical charging process which transforms the solvent density from the pure liquid state to that of the solution state.

Cavity Particle in Aqueous Solution with a Hydrophobic Solute: Structure, Energetics, and Functionals.

Endpoints density functional theory (DFT) provides a framework for calculating the excess chemical potential of a solute in solution using solvent distribution functions obtained from both physical endpoints of a hypothetical charging process which transforms the solvent density from that of the pure liquid to the solution state. In this work, the endpoints DFT equations are formulated in terms of the indirect (solvent-mediated) contribution ω() to the solute-solvent potential of mean force, and their connections are established with the conventional DFT expressions which are based on the use of direct correlation functions. ω actually corresponds to the free-energy cost to move a cavity particle (a tagged solvent molecule which interacts with the other solvent molecules but not the solute) from the bulk to the configuration of a solvent molecule relative to the solute and is a suitable variable to describe the solvent effects on the solute-solvent interactions.

Ligand Binding Thermodynamic Cycles: Hysteresis, the Locally Weighted Histogram Analysis Method, and the Overlapping States Matrix.

Free energy perturbation (FEP) simulations have been widely applied to obtain predictions of the relative binding free energy for a series of congeneric ligands binding to the same receptor, which is an essential component for the lead optimization process in computer-aided drug discovery. In the case of several congeneric ligands forming a perturbation map involving a closed thermodynamic cycle, the summation of the estimated free energy change along each edge in the cycle using Bennett acceptance ratio (BAR) usually will deviate from zero due to systematic and random errors, which is the hysteresis of cycle closure. In this work, the advanced reweighting techniques binless weighted histogram analysis method (UWHAM) and locally weighted histogram analysis method (LWHAM) are applied to provide statistical estimators of the free energy change along each edge in order to eliminate the hysteresis effect.

Absolute Protein Binding Free Energy Simulations for Ligands with Multiple Poses, a Thermodynamic Path That Avoids Exhaustive Enumeration of the Poses.

We propose a free energy calculation method for receptor-ligand binding, which have multiple binding poses that avoids exhaustive enumeration of the poses. For systems with multiple binding poses, the standard procedure is to enumerate orientations of the binding poses, restrain the ligand to each orientation, and then, calculate the binding free energies for each binding pose. In this study, we modify a part of the thermodynamic cycle in order to sample a broader conformational space of the ligand in the binding site.

Spatially-Decomposed Free Energy of Solvation Based on the Endpoint Density-Functional Method.

A spatially resolved version of the density-functional method for solvation thermodynamics is presented by extending the free-energy functional previously established in the one-dimensional, energy representation and formulating a new expression in a mixed four-dimensional representation (three dimensions for position and one dimension for energy). The space was further divided into a set of discrete regions with respect to the relative position of a solvent molecule from the solute, and the spatially decomposed energetics of solvation were analyzed for small molecules with a methyl, amine, or hydroxyl group and alanine dipeptide in solvent water. It was observed that the density of the solvation free energy is weakly dependent on the solute site in the excluded-volume region and is distinctively favorable in the first shells of the solute atoms that can readily form hydrogen bonds with water.

The UWHAM and SWHAM Software Package.

We introduce the UWHAM (binless weighted histogram analysis method) and SWHAM (stochastic UWHAM) software package that can be used to estimate the density of states and free energy differences based on the data generated by multi-state simulations. The programs used to solve the UWHAM equations are written in the C++ language and operated via the command line interface. In this paper, first we review the theoretical bases of UWHAM, its stochastic solver RE-SWHAM (replica exchange-like SWHAM)and ST-SWHAM (serial tempering-like SWHAM).

Conformational Free Energy Changes via an Alchemical Path without Reaction Coordinates.

We introduce a novel method called restrain-free energy perturbation-release (R-FEP-R) to estimate conformational free energy changes via an alchemical path, which for some conformational landscapes like those associated with cellular signaling proteins in the kinase family is more direct and readily converged than the corresponding free energy changes along the physical path. The R-FEP-R method was developed from the dual topology free energy perturbation method that is widely applied to estimate the binding free energy difference between two ligands. In R-FEP-R, the free energy change between two conformational basins is calculated by free energy perturbations that remove those atoms involved in the conformational change from their initial conformational basin while simultaneously growing them back according to the final conformational basin.

Comparing alchemical and physical pathway methods for computing the absolute binding free energy of charged ligands.

Accurately predicting absolute binding free energies of protein-ligand complexes is important as a fundamental problem in both computational biophysics and pharmaceutical discovery. Calculating binding free energies for charged ligands is generally considered to be challenging because of the strong electrostatic interactions between the ligand and its environment in aqueous solution. In this work, we compare the performance of the potential of mean force (PMF) method and the double decoupling method (DDM) for computing absolute binding free energies for charged ligands.

The Excess Chemical Potential of Water at the Interface with a Protein from End Point Simulations.

We use end point simulations to estimate the excess chemical potential of water in the homogeneous liquid and at the interface with a protein in solution. When the pure liquid is taken as the reference, the excess chemical potential of interfacial water is the difference between the solvation free energy of a water molecule at the interface and in the bulk. Using the homogeneous liquid as an example, we show that the solvation free energy for growing a water molecule can be estimated by applying UWHAM to the simulation data generated from the initial and final states (i.

The Role of Interfacial Water in Protein-Ligand Binding: Insights from the Indirect Solvent Mediated Potential of Mean Force.

Classical density functional theory (DFT) can be used to relate the thermodynamic properties of solutions to the indirect solvent mediated part of the solute-solvent potential of mean force (PMF). Standard, but powerful numerical methods can be used to estimate the solute-solvent PMF from which the indirect part can be extracted. In this work we show how knowledge of the direct and indirect parts of the solute-solvent PMF for water at the interface of a protein receptor can be used to gain insights about how to design tighter binding ligands.

Stratified UWHAM and Its Stochastic Approximation for Multicanonical Simulations Which Are Far from Equilibrium.

We describe a new analysis tool called Stratified unbinned Weighted Histogram Analysis Method (Stratified-UWHAM), which can be used to compute free energies and expectations from a multicanonical ensemble when a subset of the parallel simulations is far from being equilibrated because of barriers between free energy basins which are only rarely (or never) crossed at some states. The Stratified-UWHAM equations can be obtained in the form of UWHAM equations but with an expanded set of states. We also provide a stochastic solver, Stratified RE-SWHAM, for Stratified-UWHAM to remove its computational bottleneck.

Relationship between Solvation Thermodynamics from IST and DFT Perspectives.

Inhomogeneous solvation theory (IST) and classical density functional theory (DFT) each provide a framework for relating distribution functions of solutions to their thermodynamic properties. As reviewed in this work, both IST and DFT can be formulated in a way that use two "end point" simulations, one of the pure solvent and the other of the solution, to determine the solute chemical potential and other thermodynamic properties of the solution and of subvolumes in regions local to the solute containing hydrating waters. In contrast to IST, where expressions for the excess energy and entropy of solution are the object of analysis, in the DFT end point formulation of the problem, the solute-solvent potential of mean force (PMF) plays a central role.

Simulating Replica Exchange: Markov State Models, Proposal Schemes, and the Infinite Swapping Limit.

Replica exchange molecular dynamics is a multicanonical simulation technique commonly used to enhance the sampling of solvated biomolecules on rugged free energy landscapes. While replica exchange is relatively easy to implement, there are many unanswered questions about how to use this technique most efficiently, especially because it is frequently the case in practice that replica exchange simulations are not fully converged. A replica exchange cycle consists of a series of molecular dynamics steps of a set of replicas moving under different Hamiltonians or at different thermodynamic states followed by one or more replica exchange attempts to swap replicas among the different states.

Locally weighted histogram analysis and stochastic solution for large-scale multi-state free energy estimation.

The weighted histogram analysis method (WHAM) including its binless extension has been developed independently in several different contexts, and widely used in chemistry, physics, and statistics, for computing free energies and expectations from multiple ensembles. However, this method, while statistically efficient, is computationally costly or even infeasible when a large number, hundreds or more, of distributions are studied. We develop a locally WHAM (local WHAM) from the perspective of simulations of simulations (SOS), using generalized serial tempering (GST) to resample simulated data from multiple ensembles.

A Stochastic Solution to the Unbinned WHAM Equations.

The weighted histogram analysis method (WHAM) and unbinned versions such as the multistate Bennett acceptance ratio (MBAR) and unbinned WHAM (UWHAM) are widely used to compute free energies and expectations from data generated by independent or coupled parallel simulations. Here we introduce a replica exchange-like algorithm (RE-SWHAM) that can be used to solve the UWHAM equations stochastically. This method is capable of analyzing large data sets generated by hundreds or even thousands of parallel simulations that are too large to be "WHAMMED" using standard methods.

Connecting free energy surfaces in implicit and explicit solvent: an efficient method to compute conformational and solvation free energies.

The ability to accurately model solvent effects on free energy surfaces is important for understanding many biophysical processes including protein folding and misfolding, allosteric transitions, and protein–ligand binding. Although all-atom simulations in explicit solvent can provide an accurate model for biomolecules in solution, explicit solvent simulations are hampered by the slow equilibration on rugged landscapes containing multiple basins separated by barriers. In many cases, implicit solvent models can be used to significantly speed up the conformational sampling; however, implicit solvent simulations do not fully capture the effects of a molecular solvent, and this can lead to loss of accuracy in the estimated free energies.

Large-scale asynchronous and distributed multidimensional replica exchange molecular simulations and efficiency analysis.

We describe methods to perform replica exchange molecular dynamics (REMD) simulations asynchronously (ASyncRE). The methods are designed to facilitate large scale REMD simulations on grid computing networks consisting of heterogeneous and distributed computing environments as well as on homogeneous high-performance clusters. We have implemented these methods on NSF (National Science Foundation) XSEDE (Extreme Science and Engineering Discovery Environment) clusters and BOINC (Berkeley Open Infrastructure for Network Computing) distributed computing networks at Temple University and Brooklyn College at CUNY (the City University of New York).

Chaperone activation by unfolding.

Conditionally disordered proteins can alternate between highly ordered and less ordered configurations under physiological conditions. Whereas protein function is often associated with the ordered conformation, for some of these conditionally unstructured proteins, the opposite applies: Their activation is associated with their unfolding. An example is the small periplasmic chaperone HdeA, which is critical for the ability of enteric bacterial pathogens like Escherichia coli to survive passage through extremely acidic environments, such as the human stomach.

pH-sensitive residues in the p19 RNA silencing suppressor protein from carnation Italian ringspot virus affect siRNA binding stability.

Tombusviruses, such as Carnation Italian ringspot virus (CIRV), encode a protein homodimer called p19 that is capable of suppressing RNA silencing in their infected hosts by binding to and sequestering short-interfering RNA (siRNA) away from the RNA silencing pathway. P19 binding stability has been shown to be sensitive to changes in pH but the specific amino acid residues involved have remained unclear. Using constant pH molecular dynamics simulations, we have identified key pH-dependent residues that affect CIRV p19-siRNA binding stability at various pH ranges based on calculated changes in the free energy contribution from each titratable residue.

Probing pH-dependent dissociation of HdeA dimers.

HdeA protein is a small, ATP-independent, acid stress chaperone that undergoes a dimer-to-monomer transition in acidic environments. The HdeA monomer binds a broad range of proteins to prevent their acid-induced aggregation. To understand better HdeA's function and mechanism, we perform constant-pH molecular dynamics simulations (CPHMD) to elucidate the details of the HdeA dimer dissociation process.

Steady-state simulations using weighted ensemble path sampling.

We extend the weighted ensemble (WE) path sampling method to perform rigorous statistical sampling for systems at steady state. A straightforward steady-state implementation of WE is directly practical for simple landscapes, but not when significant metastable intermediates states are present. We therefore develop an enhanced WE scheme, building on existing ideas, which accelerates attainment of steady state in complex systems.

The "weighted ensemble" path sampling method is statistically exact for a broad class of stochastic processes and binning procedures.

The "weighted ensemble" method, introduced by Huber and Kim [Biophys. J. 70, 97 (1996)], is one of a handful of rigorous approaches to path sampling of rare events.