Multi-Level High-Throughput Screening for Discovery of Ligands That Inhibit Insulin Aggregation

Abstract

We have developed a multi-level virtual screening protocol to identify lead molecules from the FDA inactives database that can inhibit insulin aggregation. The method is based on the presence of structural and interaction specificity in non-native aggregation pathway protein−protein interactions. Some key challenges specific to the present problem, when compared with native protein association, include structural heterogeneity of the protein species involved, multiple association pathways, and relatively higher probability of conformational rearrangement of the association complex. In this multi-step method, the inactives database was first screened using the dominant pharmacophore features of previously identified molecules shown to significantly inhibit insulin aggregation nucleation by binding to its aggregation-prone conformers. We then performed ensemble docking of several low-energy ligand conformations on these aggregation-prone conformers followed by molecular dynamics simulations and binding affinity calculations on a subset of docked complexes to identify a final set of five potential lead molecules to inhibit insulin aggregation nucleation. Their effect on aggregation inhibition was extensively investigated by incubating insulin under aggregation-prone aqueous buffer conditions (low pH, high temperature). Aggregation kinetics were characterized using size exclusion chromatography and Thioflavin T fluorescence assay, and the secondary structure was determined using circular dichroism spectroscopy. Riboflavin provided the best aggregation inhibition, with 85% native monomer retention after 48 h incubation under aggregation-prone conditions, whereas the no-ligand formulation showed complete monomer loss after 36 h. Further, insulin incubated with two of the screened inactives (aspartame, riboflavin) had the characteristic α-helical dip in CD spectra, while the no-ligand formulation showed a change to β-sheet rich conformations.