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.

Quantitatively Accounting for Protein Reorganization in Computer-Aided Drug Design.

Structure-based drug design frequently operates under the assumption that a single holo structure is relevant. However, a large number of crystallographic examples clearly show that multiple conformations are possible. In those cases, the protein reorganization free energy must be known to accurately predict binding free energies for ligands.

Structural and dynamic effects of paraoxon binding to human acetylcholinesterase by X-ray crystallography and inelastic neutron scattering.

Organophosphorus (OP) compounds, including nerve agents and some pesticides, covalently bind to the catalytic serine of human acetylcholinesterase (hAChE), thereby inhibiting acetylcholine hydrolysis necessary for efficient neurotransmission. Oxime antidotes can reactivate the OP-conjugated hAChE, but reactivation efficiency can be low for pesticides, such as paraoxon (POX). Understanding structural and dynamic determinants of OP inhibition and reactivation can provide insights to design improved reactivators.

Computational study of non-conductive selectivity filter conformations and C-type inactivation in a voltage-dependent potassium channel.

C-type inactivation is a time-dependent process of great physiological significance that is observed in a large class of K+ channels. Experimental and computational studies of the pH-activated KcsA channel show that the functional C-type inactivated state, for this channel, is associated with a structural constriction of the selectivity filter at the level of the central glycine residue in the signature sequence, TTV(G)YGD. The structural constriction is allosterically promoted by the wide opening of the intracellular activation gate.

Covalent inhibition of hAChE by organophosphates causes homodimer dissociation through long-range allosteric effects.

Acetylcholinesterase (EC 3.1.1.

The Activation of c-Src Tyrosine Kinase: Conformational Transition Pathway and Free Energy Landscape.

Tyrosine kinases are important cellular signaling allosteric enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Their activity must be tightly controlled, and malfunction can lead to a variety of diseases, particularly cancer. The nonreceptor tyrosine kinase c-Src, a prototypical model system and a representative member of the Src-family, functions as complex multidomain allosteric molecular switches comprising SH2 and SH3 domains modulating the activity of the catalytic domain.

Full Atom Simulations of Spin Label Conformations.

Interpretation of EPR measurables from spin labels in terms of structure and dynamics requires knowledge of label behavior. General strategies were developed for simulation of labels used in EPR of proteins. The criteria for those simulations are (a) exhaustive sampling of rotamer space, (b) consensus of results independent of starting points, and (c) inclusion of entropy.

Generalized Scalable Multiple Copy Algorithms for Molecular Dynamics Simulations in NAMD.

Computational methodologies that couple the dynamical evolution of a set of replicated copies of a system of interest offer powerful and flexible approaches to characterize complex molecular processes. Such multiple copy algorithms (MCAs) can be used to enhance sampling, compute reversible work and free energies, as well as refine transition pathways. Widely used examples of MCAs include temperature and Hamiltonian-tempering replica-exchange molecular dynamics (-REMD and -REMD), alchemical free energy perturbation with lambda replica-exchange (FEP/λ-REMD), umbrella sampling with Hamiltonian replica exchange (US/-REMD), and string method with swarms-of-trajectories conformational transition pathways.

On the Role of Dewetting Transitions in Host-Guest Binding Free Energy Calculations.

We use thermodynamic integration (TI) and explicit solvent molecular dynamics (MD) simulation to estimate the absolute free energy of host-guest binding. In the unbound state, water molecules visit all of the internally accessible volume of the host, which is fully hydrated on all sides. Upon binding of an apolar guest, the toroidal host cavity is fully dehydrated; thus, during the intermediate λ stages along the integration, the hydration of the host fluctuates between hydrated and dehydrated states.

w-REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations.

Free energy governs the equilibrium extent of many biological processes. High barriers separating free energy minima often limit the sampling in molecular dynamics (MD) simulations, leading to inaccurate free energies. Here, we demonstrate enhanced sampling and improved free energy calculations, relative to conventional MD, using windowed accelerated MD within a Hamiltonian replica exchange framework (w-REXAMD).

Replica-Exchange Accelerated Molecular Dynamics (REXAMD) Applied to Thermodynamic Integration.

Accelerated molecular dynamics (AMD) is an efficient strategy for accelerating the sampling of molecular dynamics simulations, and observable quantities such as free energies derived on the biased AMD potential can be reweighted to yield results consistent with the original, unmodified potential. In conventional AMD the reweighting procedure has an inherent statistical problem in systems with large acceleration, where the points with the largest biases will dominate the reweighted result and reduce the effective number of data points. We propose a replica exchange of various degrees of acceleration (REXAMD) to retain good statistics while achieving enhanced sampling.

Using multistate free energy techniques to improve the efficiency of replica exchange accelerated molecular dynamics.

Replica exchange accelerated molecular dynamics (REXAMD) is a method that enhances conformational sampling while retaining at least one replica on the original potential, thus avoiding the statistical problems of exponential reweighting. In this article, we study three methods that can combine the data from the accelerated replicas to enhance the estimate of properties on the original potential: weighted histogram analysis method (WHAM), pairwise multistate Bennett acceptance ratio (PBAR), and multistate Bennett acceptance ratio (MBAR). We show that the method that makes the most efficient use of equilibrium data from REXAMD simulations is the MBAR method.

Intrinsic Free Energy of the Conformational Transition of the KcsA Signature Peptide from Conducting to Nonconducting State.

We explore a conformational transition of the TATTVGYG signature peptide of the KcsA ion selectivity filter and its GYG to AYA mutant from the conducting α-strand state into the nonconducting pII-like state using a novel technique for multidimensional optimization of transition path ensembles and free energy calculations. We find that the wild type peptide, unlike the mutant, intrinsically favors the conducting state due to G77 backbone propensities and additional hydrophobic interaction between the V76 and Y78 side chains in water. The molecular mechanical free energy profiles in explicit water are in very good agreement with the corresponding adiabatic energies from the Generalized Born Molecular Volume (GBMV) implicit solvent model.

Mapping electron paramagnetic resonance spin label conformations by the simulated scaling method.

In order to efficiently simulate spin label behavior when attached to the protein backbone we developed a novel approach that enhances local conformational sampling. The simulated scaling (SS) approach (Li, H., et al.

Simulated scaling method for localized enhanced sampling and simultaneous "alchemical" free energy simulations: a general method for molecular mechanical, quantum mechanical, and quantum mechanical/molecular mechanical simulations.

A potential scaling version of simulated tempering is presented to efficiently sample configuration space in a localized region. The present "simulated scaling" method is developed with a Wang-Landau type of updating scheme in order to quickly flatten the distributions in the scaling parameter lambdam space. This proposal is meaningful for a broad range of biophysical problems, in which localized sampling is required.