Shahila Mehboob

To validate a method for predicting the binding affinities of FabI inhibitors, three implicit solvent methods, MM-PBSA, MM-GBSA, and QM/MM-GBSA were carefully compared using 16 benzimidazole inhibitors in complex with Francisella tularensis FabI. The data suggests that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length, if only one simulation trajectory is used for each ligand. In this case, QM/MM-GBSA using 6 ns MD simulation trajectories together with GB(neck2) , PM3, and the mbondi2 radii set, generate the closest agreement with experimental values (r(2) = 0.88). However, if the three implicit solvent methods are averaged from six 1 ns MD simulations for each ligand (called "multiple independent sampling"), the prediction results are relatively insensitive to all the tested parameters. Moreover, MM/GBSA together with GB(HCT) and mbondi, using 600 frames extracted evenly from six 0.25 ns MD simulations, can also provide accurate prediction to experimental values (r(2) = 0.84). Therefore, the multiple independent sampling method can be more efficient than a single, long simulation method. Since future scaffold expansions may significantly change the benzimidazole's physiochemical properties (charges, etc.) and possibly binding modes, which may affect the sensitivities of various parameters, the relatively insensitive "multiple independent sampling method" may avoid the need of an entirely new validation study. Moreover, due to large fluctuating entropy values, (QM/)MM-P(G)BSA were limited to inhibitors' relative affinity prediction, but not the absolute affinity. The developed protocol will support an ongoing benzimidazole lead optimization program.

Enoyl-acyl carrier protein (ACP) reductase, FabI, is a key enzyme in the bacterial fatty acid biosynthesis pathway (FAS II). FabI is an NADH-dependent oxidoreductase that acts to reduce enoyl-ACP substrates in a final step of the pathway. The absence of this enzyme in humans makes it an attractive target for the development of new antibacterial agents. FabI is known to be unresponsive to structure-based design efforts due to a high degree of induced fit and a mobile flexible loop encompassing the active site. Here we discuss the development, validation, and careful application of a ligand-based virtual screen used for the identification of novel inhibitors of the Francisella tularensis FabI target. In this study, four known classes of FabI inhibitors were used as templates for virtual screens that involved molecular shape and electrostatic matching. The program ROCS was used to search a high-throughput screening library for compounds that matched any of the four molecular shape queries. Matching compounds were further refined using the program EON, which compares and scores compounds by matching electrostatic properties. Using these techniques, 50 compounds were selected, ordered, and tested. The tested compounds possessed novel chemical scaffolds when compared to the input query compounds. Several hits with low micromolar activity were identified and follow-up scaffold-based searches resulted in the identification of a lead series with submicromolar enzyme inhibition, high ligand efficiency, and a novel scaffold. Additionally, one of the most active compounds showed promising whole-cell antibacterial activity against several Gram-positive and Gram-negative species, including the target pathogen. The results of a preliminary structure-activity relationship analysis are presented.

On the basis of sequence homology studies, it has been suggested that the association of human erythrocytes alpha and beta spectrin at the tetramerization site involves interactions between helices. However, no empirical details are available, presumably due to the experimental difficulties in studying spectrin molecules because of its size and/or its structural flexibility. It has been speculated that erythrocyte tetramerization involves helical bundling rather than coiled coil association. We have used recombinant spectrin peptides to model alpha and beta spectrin to study their association at the tetramerization site. Two alpha peptides, Sp alpha 1-156 and Sp alpha 1-368, and one beta peptide, Sp beta 1898-2083, were used as model peptides to demonstrate the formation of the alpha beta complex. We also found that the replacement of R28 in Sp alpha 1-368 to give Sp alpha 1-368R28C abolished complex formation with the beta peptide. Circular dichroism techniques were used to monitor the secondary structures of the individual peptides and of the complex, and the results showed that both Sp alpha 1-156 and Sp beta 1898-2083 peptides in solution, separately, included helices that were not paired with other helices in the absence of their binding partners. However, in a mixture of Sp alpha 1-156 and Sp beta 1898-2083 and formation of the alpha beta complex, the unpaired helices associated to form coiled coils. Since the sequences of these two peptides that are involved in the coiled coil association are derived from a native protein, the information obtained from this study also provides insight toward a better understanding of naturally occurring coiled coil subunit-subunit association.