Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Effect of Knotting on the Shape of Polymers

Momentary configurations of long polymers at thermal equilibrium usually deviate from spherical symmetry and can be better described, on average, by a prolate ellipsoid. The asphericity and nature of asphericity (or prolateness) that describe these momentary ellipsoidal shapes of a polymer are determined by specific expressions involving the three principal moments of inertia calculated for configurations of the polymer. Earlier theoretical studies and numerical simulations have established that as the length of the polymer increases, the average shape for the statistical ensemble of random configurations asymptotically approaches a characteristic universal shape that depends on the solvent quality. It has been established, however, that these universal shapes differ for linear, circular, and branched chains. We investigate here the effect of knotting on the shape of cyclic polymers modeled as random isosegmental polygons. We observe that random polygons forming different knot types reach asymptotic shapes that are distinct from the ensemble average shape. For the same chain length, more complex knots are, on average, more spherical than less complex knots.