Dataset Information

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

ABSTRACT: 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.

SUBMITTER: Mobley DL

PROVIDER: S-EPMC2562444 | biostudies-literature | 2007

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

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

Mobley David L DL Chodera John D JD Dill Ken A KA

Journal of chemical theory and computation 20070101 4

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. ...[more]

PMID: 18843379

Similar Datasets

Project description:The thermodynamic reasons why membrane proteins form stable complexes inside the hydrophobic lipid bilayer remain poorly understood. This is largely because of a lack of membrane-protein systems amenable for equilibrium studies and a limited number of methods for measuring these reactions. Recently, we reported the equilibrium dimerization of the CLC-ec1 Cl-/H+ transporter in lipid bilayers (Chadda et al. 2016. eLife https://doi.org/10.7554/eLife.17438), which provided a new type of model system for studying protein association in membranes. The measurement was conducted using the subunit-capture approach, involving passive dilution of the protein in large multilamellar vesicles, followed by single-molecule photobleaching analysis of the Poisson distribution describing protein encapsulation into extruded liposomes. To estimate the fraction of dimers (FDimer ) as a function of protein density, the photobleaching distributions for the nonreactive, ideal monomer and dimer species must be known so that random co-capture probabilities can be accounted for. Previously, this was done by simulating the Poisson process of protein reconstitution into a known size distribution of liposomes composed of Escherichia coli polar lipids (EPLs). In the present study, we investigate the dependency of FDimer and ?G° on the modeling through a comparison of different liposome size distributions (EPL versus 2:1 POPE/POPG). The results show that the estimated FDimer values are comparable, except at higher densities when liposomes become saturated with protein. We then develop empirical controls to directly measure the photobleaching distributions of the nonreactive monomer (CLC-ec1 I201W/I422W) and ideal dimer (WT CLC-ec1 cross-linked by glutaraldehyde or CLC-ec1 R230C/L249C cross-linked by a disulfide bond). The measured equilibrium constants do not depend on the correction method used, indicating the robustness of the subunit-capture approach. This strategy therefore presents a model-free way to quantify protein dimerization in lipid bilayers, offering a simplified strategy in the ongoing effort to characterize equilibrium membrane-protein reactions in membranes.

Project description:We present here a greatly updated version of an earlier study on the conformational energies of protein-ligand complexes in the Protein Data Bank (PDB) [Nicklaus et al. Bioorg. Med. Chem. 1995, 3, 411-428], with the goal of improving on all possible aspects such as number and selection of ligand instances, energy calculations performed, and additional analyses conducted. Starting from about 357,000 ligand instances deposited in the 2008 version of the Ligand Expo database of the experimental 3D coordinates of all small-molecule instances in the PDB, we created a "high-quality" subset of ligand instances by various filtering steps including application of crystallographic quality criteria and structural unambiguousness. Submission of 640 Gaussian 03 jobs yielded a set of about 415 successfully concluded runs. We used a stepwise optimization of internal degrees of freedom at the DFT level of theory with the B3LYP/6-31G(d) basis set and a single-point energy calculation at B3LYP/6-311++G(3df,2p) after each round of (partial) optimization to separate energy changes due to bond length stretches vs bond angle changes vs torsion changes. Even for the most "conservative" choice of all the possible conformational energies-the energy difference between the conformation in which all internal degrees of freedom except torsions have been optimized and the fully optimized conformer-significant energy values were found. The range of 0 to ~25 kcal/mol was populated quite evenly and independently of the crystallographic resolution. A smaller number of "outliers" of yet higher energies were seen only at resolutions above 1.3 Å. The energies showed some correlation with molecular size and flexibility but not with crystallographic quality metrics such as the Cruickshank diffraction-component precision index (DPI) and R(free)-R, or with the ligand instance-specific metrics such as occupancy-weighted B-factor (OWAB), real-space R factor (RSR), and real-space correlation coefficient (RSCC). We repeated these calculations with the solvent model IEFPCM, which yielded energy differences that were generally somewhat lower than the corresponding vacuum results but did not produce a qualitatively different picture. Torsional sampling around the crystal conformation at the molecular mechanics level using the MMFF94s force field typically led to an increase in energy.

Dataset Information

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

Publications

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

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets