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.

Water Thermodynamics of Peptide Toxin Binding Sites on Ion Channels.

Peptide toxins isolated from venomous creatures, long prized as research tools due to their innate potency for ion channels, are emerging as drugs as well. However, it remains challenging to understand why peptide toxins bind with high potency to ion channels, to identify residues that are key for activity, and to improve their affinities via mutagenesis. We use WaterMap, a molecular dynamics simulation-based method, to gain computational insight into these three questions by calculating the locations and thermodynamic properties of water molecules in the peptide toxin binding sites of five ion channels.

High throughput evaluation of macrocyclization strategies for conformer stabilization.

While macrocyclization of a linear compound to stabilize a known bioactive conformation can be a useful strategy to increase binding potency, the difficulty of macrocycle synthesis can limit the throughput of such strategies. Thus computational techniques may offer the higher throughput required to screen large numbers of compounds. Here we introduce a method for evaluating the propensity of a macrocyclic compound to adopt a conformation similar that of a known active linear compound in the binding site.

Improving Accuracy, Diversity, and Speed with Prime Macrocycle Conformational Sampling.

A novel method for exploring macrocycle conformational space, Prime macrocycle conformational sampling (Prime-MCS), is introduced and evaluated in the context of other available algorithms (Molecular Dynamics, LowModeMD in MOE, and MacroModel Baseline Search). The algorithms were benchmarked on a data set of 208 macrocycles which was curated for diversity from the Cambridge Structural Database, the Protein Data Bank, and the Biologically Interesting Molecule Reference Dictionary. The algorithms were evaluated in terms of accuracy (ability to reproduce the crystal structure), diversity (coverage of conformational space), and computational speed.

Prediction of Protein-Ligand Binding Poses via a Combination of Induced Fit Docking and Metadynamics Simulations.

Ligand docking is a widely used tool for lead discovery and binding mode prediction based drug discovery. The greatest challenges in docking occur when the receptor significantly reorganizes upon small molecule binding, thereby requiring an induced fit docking (IFD) approach in which the receptor is allowed to move in order to bind to the ligand optimally. IFD methods have had some success but suffer from a lack of reliability.