(15g) Understanding the Antagonistic Mechanism of a De Novo Designed Peptide With Experimental HIV Fusion Inhibition Via Molecular Simulation

Authors: 
Khoury, G. A., Pennsylvania State University-University Park
Tamamis, P., Texas A&M University
Smadbeck, J., Princeton University
Szekely, Z., Rutgers University
Sinko, P. J., Rutgers University
Floudas, C. A., Princeton University



A noble goal in protein design is to predict the potency and understand in detail the mechanism of inhibition of a candidate antagonist. Previous de novo design has produced several HIV fusion inhibitors targeting HIV-gp411. These inhibitors were experimentally confirmed to inhibit viral entry.  Here, we probe, through computational studies based on molecular dynamics simulations the mechanism of inhibition of SQ175, a design among the best in antagonistic activity.

Multiple independent explicit solvent (TIP3) molecular dynamics simulations of SQ175 bound to the trimeric coiled-coil are performed in the AMBER2 MD package. Using the snapshots from the trajectories, interaction maps were constructed identifying key SQ175:gp41 hydrophobic and electrostatic interactions. Binding free energy calculations of SQ175 confirm its affinity for the trimeric coiled-coil of gp41. Subsequently, in CHARMM3 and using generalized Born implicit solvent4, simulations of the C-terminus binding fragment bound to the trimeric coiled-coil of gp41 were performed. These simulations aimed at evaluating the native interactions between the C-terminus binding fragment:trimeric core and to construct a map denoting the minimum Cα-Cα distances the C-terminus binding fragment residues experience with respect to the gp41 trimeric coiled-coil. Next, multiple independent docking simulations of the C-terminus binding fragment starting 10-12 Å away from its fusion intermediate structure were performed, where after the MD, the C-terminus binding fragment was found to spontaneously dock within 1.5 Å of its conformation in the crystal structure in one simulation, indicating the system behaves well. Then, in the presence of SQ175, 12 independent docking simulations were performed starting from 12 different initial coordinates to mimic the C-terminus binding fragment approaching gp41. No restraints were implemented to drive them together. The SQ175 peptide was observed to dramatically reduce the number of native C-terminus binding fragment:gp41 contacts of which several are necessary for HIV fusion and subsequent entry, as observed by minimum Cα-Cα distances. The antagonist SQ175 was observed to stably bind to the trimeric coiled-coil of gp41 and physically block a large fraction of the native C-terminus binding fragment:gp41 contacts. Based on these results, new insights for future generations of fusion inhibitors are described.

 

References

1.         Bellows, M. L.; Taylor, M. S.; Cole, P. A.; Shen, L.; Siliciano, R. F.; Fung, H. K.; Floudas, C. A., Discovery of Entry Inhibitors for HIV-1 via a New De Novo Protein Design Framework. Biophysical Journal 2010, 99 (10), 3445-3453.

2.         Weiner, P. K.; Kollman, P. A., AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions. Journal of Computational Chemistry 1981, 2 (3), 287-303.

3.         MacKerell, J., A. D., Brooks, B., Brooks, III, C.L., Nilsson, L., Roux, B., Won, Y., and Karplus, M., CHARMM: The Energy Function and Its Parameterization with an Overview of the Program. In The Encyclopedia of Computational Chemistry, al., P. v. R. S. e., Ed. John Wiley & Sons: Chichester: 1998; Vol. 1, pp 271-277.

4.         (a) Chen, J.; Im, W.; Brooks, C. L., 3rd, Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field. Journal of the American Chemical Society 2006, 128 (11), 3728-36; (b) Im, W.; Lee, M. S.; Brooks, C. L., 3rd, Generalized born model with a simple smoothing function. J Comput Chem 2003, 24 (14), 1691-702.