Ranking ligand affinity for the DNA minor groove by experiment and simulation

The structural and thermodynamic basis for the strength and selectivity of the interactions of minor-groove binders (MGBs) with DNA is not fully understood. In 2003 we reported the first example of a thiazole containing MGB that bound in a phase shifted pattern that spanned 6 base-pairs rather than the usual 4 (for tricyclic distamycin-like compounds). Since then, using DNA footprinting, nuclear magnetic resonance spectroscopy, isothermal titration calorimetry and molecular dynamics, we have established that the flanking bases around the central 4 being read by the ligand have subtle effects on recognition. We have investigated the effect of these flanking sequences on binding and the reasons for the differences and established a computational method to rank ligand affinity against varying DNA sequences

Strength of Hydrogen Bond Network Takes Crucial Roles in the Dissociation Process of Inhibitors from the HIV-1 Protease Binding Pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions