Binding free energy calculations using MMPB/GBSA approaches for PAMAM-G4-drug complexes at neutral, basic and acid pH conditions

Dendrimers are synthetic macromolecules with a highly-branched structure and high concentration of surface groups. Among dendrimers, Poly(amidoamine) (PAMAM) has received substantial attention as a novel drug carrier and delivery system. Depending on the generation and type of terminal groups, dendrimer toxicity could change and include cytotoxicity. Although PAMAM is water soluble, molecular modeling of the dendrimer-drug complex is considered challenging for exploring the conformational mobility of dendrimers and atomic specific interactions during the dendrimer-drug association. However, conventional protocols for predicting binding affinities have been designed for small protein molecules or protein-protein complexes that can be applied to study the dendrimer-drug association. In this work, we performed docking calculations for a set of 94 previously reported compounds on PAMAM of fourth generation (G4-PAMAM) to select six compounds, cromoglicic acid (CRO) − a mast cell stabilizer, Fusidic acid (FUS) − a bacteriostatic antibiotic, and Methotrexate (MTX) − a chemotherapy agent and immune system suppressant, which have the highest affinities for G4-PAMAM, and Lidocaine (LDC) − used to numb tissue in a specific area and for ventricular tachycardia treatment, Metoprolol (MET) − a β₁ receptor blocker, and Pindolol (PIN) − a β blocker, which have the lowest affinities for the G4-PAMAM dendrimer, to perform MD simulations combined with the molecular mechanics generalized/Poisson-Boltzmann surface area MMGBSA/MMPBSA approach to investigate the interactions of generating 4 charge-neutral, charge-basic and charge-acid G4-PAMAM dendrimers. In addition, to validate these theoretical G4-PAMAM-drug complexes, the complexes were experimentally conjugated to determine their stability in aqueous solubility studies immediately and over one year. Our results show that among the different commercial drugs, both charged and neutral PAMAM have the most favorable binding free energies for CRO, MTX, and FUS, which appears to be due to a complex counterbalance of electrostatics and van der Waals interactions. These theoretical and aqueous solubility studies supported the high affinity of methotrexate for the G4-PAMAM-drug due to its carboxyl and aryl moieties that favor its accommodation by noncovalent interactions.