In this paper we propose a nonlinear elasticity model of macromolecular conformational change (deformation) induced by electrostatic forces generated by an implicit solvation model. The Poisson-Boltzmann equation for the electrostatic potential is analyzed in a domain varying with the elastic deformation of molecules, and a new continuous model of the electrostatic forces is developed to ensure solvability of the nonlinear elasticity equations. We derive the estimates of electrostatic forces corresponding to four types of perturbations to an electrostatic potential field, and establish the existance of an equilibrium configuration using a fixed-point argument, under the assumption that the change in the ionic strength and charges due to the additional molecules causing the deformation are sufficiently small. The results are valid for elastic models with arbitrarily complex dielectric interfaces and cavities, and can be generalized to large elastic deformation caused by high ionic strength, large charges, and strong external fields by using continuation methods.

With the rise in antibiotic resistance, there is interest in discovering new drugs active against new targets. Here, we investigate the dynamic structures of three isoprenoid synthases from Mycobacterium tuberculosis using molecular dynamics (MD) methods with a view to discovering new drug leads. Two of the enzymes, cis-farnesyl diphosphate synthase (cis-FPPS) and cis-decaprenyl diphosphate synthase (cis-DPPS), are involved in bacterial cell wall biosynthesis, while the third, tuberculosinyl adenosine synthase (Rv3378c), is involved in virulence factor formation. The MD results for these three enzymes were then compared with previous results on undecaprenyl diphosphate synthase (UPPS) by means of active site volume fluctuation and principal component analyses. In addition, an analysis of the binding of prenyl diphosphates to cis-FPPS, cis-DPPS, and UPPS utilizing the new MD results is reported. We also screened libraries of inhibitors against cis-DPPS, finding ~1 μm inhibitors, and used the receiver operating characteristic-area under the curve (ROC-AUC) method to test the predictive power of X-ray and MD-derived cis-DPPS receptors. We found that one compound with potent M. tuberculosis cell growth inhibition activity was an IC(50) ~0.5- to 20-μm inhibitor (depending on substrate) of cis-DPPS, a ~660-nm inhibitor of Rv3378c as well as a 4.8-μm inhibitor of cis-FPPS, opening up the possibility of multitarget inhibition involving both cell wall biosynthesis and virulence factor formation.