A Computational Physics-based Approach to Predict Unbound Brain-to-Plasma Partition Coefficient, Kp,uu.

The blood-brain barrier (BBB) plays a critical role in preventing harmful endogenous and exogenous substances from penetrating the brain. Optimal brain penetration of small-molecule central nervous system (CNS) drugs is characterized by a high unbound brain/plasma ratio (Kp,uu). While various medicinal chemistry strategies and in silico models have been reported to improve BBB penetration, they have limited application in predicting Kp,uu directly. We describe a physics-based computational approach, a quantum mechanics (QM)-based energy of solvation (E-sol), to predict Kp,uu. Prospective application of this method in internal CNS drug discovery programs highlights the utility and accuracy of this new method, which showed a categorical accuracy of 79% and an R2 of 0.61 from a linear regression model.

A network of molecular switches controls the activation of the two-component response regulator NtrC.

Recent successes in simulating protein structure and folding dynamics have demonstrated the power of molecular dynamics to predict the long timescale behaviour of proteins. Here, we extend and improve these methods to predict molecular switches that characterize conformational change pathways between the active and inactive state of nitrogen regulatory protein C (NtrC). By employing unbiased Markov state model-based molecular dynamics simulations, we construct a dynamic picture of the activation pathways of this key bacterial signalling protein that is consistent with experimental observations and predicts new mutants that could be used for validation of the mechanism. Moreover, these results suggest a novel mechanistic paradigm for conformational switching.

Elucidating Ligand-Modulated Conformational Landscape of GPCRs Using Cloud-Computing Approaches.

G-protein-coupled receptors (GPCRs) are a versatile family of membrane-bound signaling proteins. Despite the recent successes in obtaining crystal structures of GPCRs, much needs to be learned about the conformational changes associated with their activation. Furthermore, the mechanism by which ligands modulate the activation of GPCRs has remained elusive. Molecular simulations provide a way of obtaining detailed an atomistic description of GPCR activation dynamics. However, simulating GPCR activation is challenging due to the long timescales involved and the associated challenge of gaining insights from the "Big" simulation datasets. Here, we demonstrate how cloud-computing approaches have been used to tackle these challenges and obtain insights into the activation mechanism of GPCRs. In particular, we review the use of Markov state model (MSM)-based sampling algorithms for sampling milliseconds of dynamics of a major drug target, the G-protein-coupled receptor ß2-AR. MSMs of agonist and inverse agonist-bound ß2-AR reveal multiple activation pathways and how ligands function via modulation of the ensemble of activation pathways. We target this ensemble of conformations with computer-aided drug design approaches, with the goal of designing drugs that interact more closely with diverse receptor states, for overall increased efficacy and specificity. We conclude by discussing how cloud-based approaches present a powerful and broadly available tool for studying the complex biological systems routinely.

Cloud computing approaches for prediction of ligand binding poses and pathways.

We describe an innovative protocol for ab initio prediction of ligand crystallographic binding poses and highly effective analysis of large datasets generated for protein-ligand dynamics. We include a procedure for setup and performance of distributed molecular dynamics simulations on cloud computing architectures, a model for efficient analysis of simulation data, and a metric for evaluation of model convergence. We give accurate binding pose predictions for five ligands ranging in affinity from 7 nM to > 200 µM for the immunophilin protein FKBP12, for expedited results in cases where experimental structures are difficult to produce. Our approach goes beyond single, low energy ligand poses to give quantitative kinetic information that can inform protein engineering and ligand design.

Cloud-based simulations on Google Exacycle reveal ligand modulation of GPCR activation pathways.

Simulations can provide tremendous insight into the atomistic details of biological mechanisms, but micro- to millisecond timescales are historically only accessible on dedicated supercomputers. We demonstrate that cloud computing is a viable alternative that brings long-timescale processes within reach of a broader community. We used Google's Exacycle cloud-computing platform to simulate two milliseconds of dynamics of a major drug target, the G-protein-coupled receptor ß2AR. Markov state models aggregate independent simulations into a single statistical model that is validated by previous computational and experimental results. Moreover, our models provide an atomistic description of the activation of a G-protein-coupled receptor and reveal multiple activation pathways. Agonists and inverse agonists interact differentially with these pathways, with profound implications for drug design.

Thermodynamic integration to predict host-guest binding affinities.

An alchemical free energy method with explicit solvent molecular dynamics simulations was applied as part of the blind prediction contest SAMPL3 to calculate binding free energies for seven guests to an acyclic cucurbit-[n]uril host. The predictions included determination of protonation states for both host and guests, docking pose generation, and binding free energy calculations using thermodynamic integration. We found a root mean square error (RMSE) of 3.6 kcal mol(-1) from the reference experimental results, with an R(2) correlation of 0.51. The agreement with experiment for the largest contributor to this error, guest 6, is improved by 1.7 kcal mol(-1) when a periodicity-induced free energy correction is applied. The corrections for the other ligands were significantly smaller, and altogether the RMSE was reduced by 0.4 kcal mol(-1). We link properties of the host-guest systems during simulation to errors in the computed free energies. Overall, we show that charged host-guest systems studied here, initialized in unconfirmed docking poses, present a challenge to accurate alchemical simulation methods.

Independent-Trajectory Thermodynamic Integration: a practical guide to protein-drug binding free energy calculations using distributed computing.

The Independent-Trajectory Thermodynamic Integration (IT-TI) approach for free energy calculation with distributed computing is described. IT-TI utilizes diverse conformational sampling obtained from multiple, independent simulations to obtain more reliable free energy estimates compared to single TI predictions. The latter may significantly under- or over-estimate the binding free energy due to finite sampling. We exemplify the advantages of the IT-TI approach using two distinct cases of protein-ligand binding. In both cases, IT-TI yields distributions of absolute binding free energy estimates that are remarkably centered on the target experimental values. Alternative protocols for the practical and general application of IT-TI calculations are investigated. We highlight a protocol that maximizes predictive power and computational efficiency.

Effects of Biomolecular Flexibility on Alchemical Calculations of Absolute Binding Free Energies.

The independent trajectory thermodynamic integration (IT-TI) approach (Lawrenz et. al J. Chem. Theory. Comput. 2009, 5:1106-1116(1)) for free energy calculations with distributed computing is employed to study two distinct cases of protein-ligand binding: first, the influenza surface protein N1 neuraminidase bound to the inhibitor oseltamivir, and second, the M. tuberculosis enzyme RmlC complexed with the molecule CID 77074. For both systems, finite molecular dynamics (MD) sampling and varied molecular flexibility give rise to IT-TI free energy distributions that are remarkably centered on the target experimental values, with a spread directly related to protein, ligand, and solvent dynamics. Using over 2 µs of total MD simulation, alternative protocols for the practical, general implementation of IT-TI are investigated, including the optimal use of distributed computing, the total number of alchemical intermediates, and the procedure to perturb electrostatics and van der Waals interactions. A protocol that maximizes predictive power and computational efficiency is proposed. IT-TI outperforms traditional TI predictions and allows a straightforward evaluation of the reliability of free energy estimates. Our study has broad implications for the use of distributed computing in free energy calculations of macromolecular systems.

Porous protein frameworks with unsaturated metal centers in sterically encumbered coordination sites.

Described is an engineered metal-binding protein, MBPPhen2, which forms porous crystalline frameworks that feature coordinatively unsaturated Zn- and Ni-centers.

Impact of calcium on N1 influenza neuraminidase dynamics and binding free energy.

The highly pathogenic influenza strains H5N1 and H1N1 are currently treated with inhibitors of the viral surface protein neuraminidase (N1). Crystal structures of N1 indicate a conserved, high affinity calcium binding site located near the active site. The specific role of this calcium in the enzyme mechanism is unknown, though it has been shown to be important for enzymatic activity and thermostability. We report molecular dynamics (MD) simulations of calcium-bound and calcium-free N1 complexes with the inhibitor oseltamivir (marketed as the drug Tamiflu), independently using both the AMBER FF99SB and GROMOS96 force fields, to give structural insight into calcium stabilization of key framework residues. Y347, which demonstrates similar sampling patterns in the simulations of both force fields, is implicated as an important N1 residue that can "clamp" the ligand into a favorable binding pose. Free energy perturbation and thermodynamic integration calculations, using two different force fields, support the importance of Y347 and indicate a +3 to +5 kcal/mol change in the binding free energy of oseltamivir in the absence of calcium. With the important role of structure-based drug design for neuraminidase inhibitors and the growing literature on emerging strains and subtypes, inclusion of this calcium for active site stability is particularly crucial for computational efforts such as homology modeling, virtual screening, and free energy methods.

Independent-Trajectories Thermodynamic-Integration Free-Energy Changes for Biomolecular Systems: Determinants of H5N1 Avian Influenza Virus Neuraminidase Inhibition by Peramivir.

Free-energy changes are essential physicochemical quantities for understanding most biochemical processes. Yet, the application of accurate thermodynamic-integration (TI) computation to biological and macromolecular systems is limited by finite-sampling artifacts. In this paper, we employ independent-trajectories thermodynamic-integration (IT-TI) computation to estimate improved free-energy changes and their uncertainties for (bio)molecular systems. IT-TI aids sampling statistics of the thermodynamic macrostates for flexible associating partners by ensemble averaging of multiple, independent simulation trajectories. We study peramivir (PVR) inhibition of the H5N1 avian influenza virus neuraminidase flexible receptor (N1). Binding site loops 150 and 119 are highly mobile, as revealed by N1-PVR 20-ns molecular dynamics. Due to such heterogeneous sampling, standard TI binding free-energy estimates span a rather large free-energy range, from a 19% underestimation to a 29% overestimation of the experimental reference value (-62.2 +/- 1.8 kJ mol(-1)). Remarkably, our IT-TI binding free-energy estimate (-61.1 +/- 5.4 kJ mol(-1)) agrees with a 2% relative difference. In addition, IT-TI runs provide a statistics-based free-energy uncertainty for the process of interest. Using approximately 800 ns of overall sampling, we investigate N1-PVR binding determinants by IT-TI alchemical modifications of PVR moieties. These results emphasize the dominant electrostatic contribution, particularly through the N1 E277-PVR guanidinium interaction. Future drug development may be also guided by properly tuning ligand flexibility and hydrophobicity. IT-TI will allow estimation of relative free energies for systems of increasing size, with improved reliability by employing large-scale distributed computing.