Rigid‐body methods, particularly Fourier correlation techniques, are very efficient for docking bound (co‐crystallized) protein conformations using measures of surface complementarity as the target function. However, when docking unbound (separately crystallized) conformations, the method generally yields hundreds of false positive structures with good scores but high root mean square deviations (RMSDs). This paper describes a two‐step scoring algorithm that can discriminate near‐native conformations (with less than 5 Å RMSD) from other structures. The first step includes two rigid‐body filters that use the desolvation free energy and the electrostatic energy to select a manageable number of conformations for further processing, but are unable to eliminate all false positives. Complete discrimination is achieved in the second step that minimizes the molecular mechanics energy of the retained structures, and re‐ranks them with a combined free‐energy function which includes electrostatic, solvation, and van der Waals energy terms. After minimization, the improved fit in near‐native complex conformations provides the free‐energy gap required for discrimination. The algorithm has been developed and tested using docking decoys, i.e., docked conformations generated by Fourier correlation techniques. The decoy sets are available on the web for testing other discrimination procedures. Proteins 2000;40:525–537. © 2000 Wiley‐Liss, Inc.

The role of desolvation in protein binding kinetics is investigated using Brownian dynamics simulations in complexes in which the electrostatic interactions are relatively weak. We find that partial desolvation, modeled by a short-range atomic contact potential, is not only a major contributor to the binding free energy but also substantially increases the diffusion-limited rate for complexes in which long-range electrostatics is weak. This rate enhancement is mostly due to weakly specific pathways leading to a low free-energy attractor, i.e., a precursor state before docking. For alpha-chymotrypsin and human leukocyte elastase, both interacting with turkey ovomucoid third domain, we find that the forward rate constant associated with a collision within a solid angle phi around their corresponding attractor approaches 10(7) and 10(6) M(-1)s(-1), respectively, in the limit phi approximately 2 degrees. Because these estimates agree well with experiments, we conclude that the final bound conformation must be preceded by a small set of well-defined diffusion-accessible precursor states. The inclusion of the otherwise repulsive desolvation interaction also explains the lack of aggregation in proteins by restricting nonspecific association times to approximately 4 ns. Under the same reaction conditions but without short range forces, the association rate would be only approximately 10(3) M(-1)s(-1). Although desolvation increases these rates by three orders of magnitude, desolvation-mediated association is still at least 100-fold slower than the electrostatically assisted binding in complexes such as barnase and barstar.

In recent years, molecular dynamics simulations of proteins in explicit mixed solvents have been applied to various problems in protein biophysics and drug discovery, including protein folding, protein surface characterization, fragment screening, allostery, and druggability assessment. In this study, we perform a systematic study on how mixtures of organic solvent probes in water can reveal cryptic ligand binding pockets that are not evident in crystal structures of apo proteins. We examine a diverse set of eight PDB proteins that show pocket opening induced by ligand binding and investigate whether solvent MD simulations on the apo structures can induce the binding site observed in the holo structures. The cosolvent simulations were found to induce conformational changes on the protein surface, which were characterized and compared with the holo structures. Analyses of the biological systems, choice of probes and concentrations, druggability of the resulting induced pockets, and application to drug discovery are discussed here.

When a complex is constructed from the separately determined rigid structures of a receptor and its ligand, some key side chains are usually in wrong positions. These distortions of the interface yield an apparent loss in affinity and would unfavorably affect the kinetics of association. It is generally assumed that the interacting proteins should drive the appropriate conformational changes, leading to their complementarity, but this hypothesis does not explain their fast association rates. However, nanosecond explicit solvent molecular dynamics simulations of misfolded surface side chains from the independently solved structures of barstar, bovine pancreatic trypsin inhibitor, and lysozyme show that even before any receptor-ligand interaction, key side chains frequently visit the rotamer conformations seen in the complex. We show that these simple structural motifs can reconcile most of the binding affinity required for a rapid and highly specific association process. Side chains amenable to induced fit are also identified. These results corroborate that solvent-side chain interactions play a critical role in the recognition process. Our findings are also supported by crystallographic data.

HIV-1 gp120 undergoes multiple conformational changes both before and after binding to the host CD4 receptor. BMS-626529 is an attachment inhibitor (AI) in clinical development (administered as prodrug BMS-663068) that binds to HIV-1 gp120. To investigate the mechanism of action of this new class of antiretroviral compounds, we constructed homology models of unliganded HIV-1 gp120 (UNLIG), a pre-CD4 binding-intermediate conformation (pCD4), a CD4 bound-intermediate conformation (bCD4), and a CD4/co-receptor-bound gp120 (LIG) from a series of partial structures. We also describe a simple pathway illustrating the transition between these four states. Guided by the positions of BMS-626529 resistance substitutions and structure–activity relationship data for the AI series, putative binding sites for BMS-626529 were identified, supported by biochemical and biophysical data. BMS-626529 was docked into the UNLIG model and molecular dynamics simulations were used to demonstrate the thermodynamic stability of the different gp120 UNLIG/BMS-626529 models. We propose that BMS-626529 binds to the UNLIG conformation of gp120 within the structurally conserved outer domain, under the antiparallel β20–β21 sheet, and adjacent to the CD4 binding loop. Through this binding mode, BMS-626529 can inhibit both CD4-induced and CD4-independent formation of the “open state” four-stranded gp120 bridging sheet, and the subsequent formation and exposure of the chemokine co-receptor binding site. This unique mechanism of action prevents the initial interaction of HIV-1 with the host CD4+ T cell, and subsequent HIV-1 binding and entry. Our findings clarify the novel mechanism of BMS-626529, supporting its ongoing clinical development. Proteins 2015; 83:331–350. © 2014 Wiley Periodicals, Inc.

Long-chain fatty acids (FAs) with low water solubility require fatty-acid-binding proteins (FABPs) to transport them from cytoplasm to the mitochondria for energy production. However, the precise mechanism by which these proteins recognize the various lengths of simple alkyl chains of FAs with similar high affinity remains unknown. To address this question, we employed a newly developed calorimetric method for comprehensively evaluating the affinity of FAs, sub-Angstrom X-ray crystallography to accurately determine their 3D structure, and energy calculations of the coexisting water molecules using the computer program WaterMap. Our results clearly showed that the heart-type FABP (FABP3) preferentially incorporates a U-shaped FA of C10–C18 using a lipid-compatible water cluster, and excludes longer FAs using a chain-length-limiting water cluster. These mechanisms could help us gain a general understanding of how proteins recognize diverse lipids with different chain lengths.

Homology modeling of G protein-coupled receptors is becoming a widely used tool in drug discovery. However, unrefined models built using the bovine rhodopsin crystal structure as the template, often have binding sites that are too small to accommodate known ligands. Here, we present a novel systematic method to refine model active sites based on a pressure-guided molecular dynamics simulation. A distinct advantage of this approach is the ability to introduce systematic perturbations in model backbone atoms in addition to side chain adjustments. The method is validated on two test cases: (1) docking of retinal into an MD-relaxed structure of opsin and (2) docking of known ligands into a homology model of the CCR2 receptor. In both cases, we show that the MD expansion algorithm makes it possible to dock the ligands in poses that agree with the crystal structure or mutagenesis data.

No abstract