CORCEMA refinement of the bound ligand conformation within the protein binding pocket in reversibly forming weak complexes using STD-NMR intensities

We describe an intensity-restrained optimization procedure for refining approximate structures of ligands within the protein binding pockets using STD-NMR intensity data on reversibly forming weak complexes. In this approach, the global minimum for the bound-ligand conformation is obtained by a hybrid structure refinement method involving CORCEMA calculation of intensities and simulated annealing optimization of torsion angles of the bound ligand using STD-NMR intensities as experimental constraints and the NOE R-factor as the pseudo-energy function to be minimized. This method is illustrated using simulated STD data sets for typical carbohydrate and peptide ligands. Our procedure also allows for the optimization of side chain torsion angles of protein residues within the binding pocket. This procedure is useful in refining and improving initial models based on crystallography or computer docking or other algorithms to generate models for the bound ligand (e.g., a lead compound) within the protein binding pocket compatible with solution STD-NMR data. This method may facilitate structure-based drug design efforts.