Search

Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da that weakly bind to the protein of interest. This technique is challenging for computational docking methods as binding is determined by only a few specific interactions. Inaccuracies in the energy function or slight deviations in the docking pose can lead to the prediction of incorrect binding or difficulties in ranking fragments inin silicoscreening. Here we test RosettaLigand by docking a series of fragments to a cysteine-depleted variant of the TIM-barrel protein, HisF. We compare the computational results with experimental NMR spectroscopy screens. NMR spectroscopy gives details on binding affinities of individual ligands, which allows assessment of the ligand-ranking ability by RosettaLigand, and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigand’s ability to identify plausible binding poses. From a library screen of 3456 fragments, we identified a set of 31 ligands with intrinsic affinities to HisF with dissociation constants as low as 400 µM. The same library of fragments was blindly screenedin silico. RosettaLigand was able to rank binders before non-binders with an area under the curve (AUC) of the receiver operating characteristics (ROC) of 0.74. The docking poses observed for binders agreed with the binding pocket identified by NMR chemical shift perturbations for all fragments. Taken together, these results provide a baseline performance of RosettaLigand in a fragment-based drug discovery setting.

The trans- and juxtamembrane domains (TMD/JMD) of human epidermal growth factor receptors (HER/EGFR/ErbB) regulate the helix-helix dimerization interface to control receptor tyrosine kinase (RTK) activation. Patient-identified non-small cell lung cancer (NSCLC) mutations occur all along the EGFR (HER1) and HER2 TMD and JMD, but our structural understanding of how these patient-derived TMD/JMD mutations hijack wild-type RTK function to drive cancer is lacking. In this work, we integrate tools in experimental biophysics and computational modeling to evaluate the structural mechanism(s) driving these interactions in atomic detail. We use nuclear magnetic resonance (NMR) spectroscopy to assess how oncogenic mutations within the TMD/JMD perturb the architecture/interactions of the dimer interface. NMR-derived models are used to seed molecular dynamics (MD) simulations to probe how mutations affect TMD association and activation/stabilization via dimerization in a realistic membrane environment. We apply this strategy and focus our initial efforts on the EGFRA647T variant. We computed the free energy landscape from our simulations and show that the active state for the EGFRA647T-HER2WT heterodimer is stabilized by 10 kcal/mol. Conversely, the wild-type EGFR-HER2 heterodimer is stabilized by 5 kcal/mol for the inactive state suggesting the propensity of the EGFRA647T mutant favors the active state. These results are part of a multidisciplinary collaboration that share the common goal of constructing holistic models that detail how TMD/JMD variants are aberrantly activated in cancer. Citation Format: Christina Mercado, Abdullah Al Mamun, Younghee Park, Christine M. Lovly, Jens Meiler, Kaitlyn Ledwitch. Integrative structural models of oncogenic mutants within the trans- and juxta-membrane region of human epidermal growth factor receptors (HERs) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3173.

SummarySite-directed spin labeling of proteins via non-canonical amino acids (ncAAs) is a non-traditional method for the measurement of pseudocontact shifts (PCSs) by nuclear magnetic resonance (NMR) spectroscopy. PCSs provide long-range distance and orientational information between a paramagnetic center and protein nuclei that can be used as restraints for computational structural modeling techniques. Here, we present the first experimental structure of an ncAA chemically linked to a lanthanide tag conjugated to the protein, T4-Lysozyme (T4L). T4L was crystallized with a cyclen-based C3 tag coordinated to the paramagnetic ion terbium (Tb3+). The paramagnetic C3-lanthanide tag generated PCSs measured at four different ncAA sites. We show that the addition of these restraints improves structure prediction protocols for T4L using the RosettaNMR framework. Generated models provide insight into T4L conformational flexibility sampled in solution. This integrative modeling protocol is readily transferable to larger proteins. Methods to predict protein structures are advancing into an exciting arena such that reliable experimental data will play important roles for evaluating the biophysical relevance of predicted structural models. Our contribution here caters to the growing interest in using ncAAs for a range of biophysical studies, and these methods can be readily transferred to larger protein systems of interest.

The translocator protein (TSPO), previously known as the peripheral benzodiazepine receptor (PBR), is a membrane protein located on the outer mitochondrial membrane. Experimentally-derived structures of mouse TSPO (mTSPO) and its homologs from bacterial species have been determined by NMR spectroscopy and X-ray crystallography, respectively. These structures and ligand interactions within the TSPO binding pocket display distinct differences. Here, we leverage experimental and computational studies to derive a unified structural model of mTSPO in the presence and absence of the TSPO ligand, PK11195, and study the effects of DPC detergent micelles on the TSPO structure and ligand binding. From this work, we conclude that that the lipid-mimetic system used to solubilize mTSPO for NMR studies thermodynamically destabilizes the protein, introduces structural perturbations, and alters the characteristics of ligand binding. Furthermore, we used Rosetta to construct a unified mTSPO model that reconciles deviating features of the mammalian and bacterial TSPO. These deviating features are likely a consequence of the detergent system used for structure determination of mTSPO by NMR. The unified mTSPO model agrees with available experimental NMR data, appears to be physically realistic (i.e. thermodynamically not frustrated as judged by the Rosetta energy function), and simultaneously shares the structural features observed in sequence-conserved regions of the bacterial proteins. Finally, we identified the binding site for an imaging ligand VUIIS8310 that is currently positioned for clinical translation using NMR spectroscopy and propose a computational model of the VUIIS8310-mTSPO complex.

Parkinson disease (PD) is a progressive, neurodegenerative disorder affecting over 6.1 million people worldwide. Although the cause of PD remains unclear, studies of highly-penetrant mutations identified in early-onset familial parkinsonism have contributed to our understanding of the molecular mechanisms underlying disease pathology. Dopamine (DA) transporter (DAT) deficiency syndrome (DTDS) is a distinct type of infantile parkinsonism-dystonia that shares key clinical features with PD, including motor deficits (progressive bradykinesia, tremor, hypomimia) and altered DA neurotransmission. Here, we define structural, functional, and behavioral consequences of a Cys substitution at R445 in human DAT (hDAT R445C), identified in a patient with DTDS. We found that this R445 substitution disrupts a phylogenetically conserved intracellular (IC) network of interactions that compromise the hDAT IC gate. This is demonstrated by both Rosetta molecular modeling and fine-grained simulations using hDAT R445C, as well as EPR analysis and X-ray crystallography of the bacterial homolog leucine transporter. Notably, the disruption of this IC network of interactions supported a channel-like intermediate of hDAT and compromised hDAT function. We demonstrate thatDrosophila melanogasterexpressing hDAT R445C show impaired hDAT activity, which is associated with DA dysfunction in isolated brains and with abnormal behaviors monitored at high-speed time resolution.We show that hDAT R445CDrosophilaexhibit motor deficits, lack of motor coordination (i.e. flight coordination) and phenotypic heterogeneity in these behaviors that is typically associated with DTDS and PD. These behaviors are linked with altered dopaminergic signaling stemming from loss of DA neurons and decreased DA availability. We rescued flight coordination through enhanced DAT surface expressionviathe lysosomal inhibitor chloroquine. Together, these studies shed light on how a DTDS-linked DAT mutation underlies DA dysfunction and, more broadly, the clinical phenotypes shared by DTDS and PD.

P-glycoprotein (Pgp) is a multidrug resistance transporter that limits the penetration of a wide range of neurotherapeutics into the brain including opioids. The diphenylpropylamine opioids methadone and loperamide are structurally similar, but loperamide has about a 4-fold higher Pgp-mediated transport rate. In addition to these differences, they showed significant differences in their effects on Pgp-mediated adenosine triphosphate (ATP) hydrolysis. The activation of Pgp-mediated ATP hydrolysis by methadone was monophasic, whereas loperamide activation of ATP hydrolysis was biphasic implying methadone has a single binding site and loperamide has 2 binding sites on Pgp. Quenching of tryptophan fluorescence with these drugs and digoxin showed competition between the opioids and that loperamide does not compete for the digoxin-binding site. Acrylamide quenching of tryptophan fluorescence to probe Pgp conformational changes revealed that methadone- and loperamide-induced conformational changes were distinct. These results were used to develop a model for Pgp-mediated transport of methadone and loperamide where opioid binding and conformational changes are used to explain the differences in the opioid transport rates between methadone and loperamide.

Background Cytochrome P450 (P450) BM3, from Bacillus megaterium, catalyzes a wide range of chemical reactions and is routinely used as a model system to study mammalian P450 reactions and structure. Methods The metabolism of 2,6-di-tert-butyl-4-hydroperoxy-4-methyl-2,5-cyclohexadienone (BHTOOH) and 2-tert-butyl-4-hydroperoxy-4-methyl-2,5-cyclohexadien-1-one (BMPOOH) was examined with P450 BM3 and with the conserved T268 and F87 residues mutated to investigate their effects on organic hydroperoxide metabolism. To determine the effects of the mutations on the active site volume and architecture, the X-ray crystal structure of the F87A/T268A P450 BM3 heme domain (BMP) was determined and compared to previous structures. To investigate the interactions of the substrates with the F87 and T268 residues, BHTOOH and BMPOOH were docked into the BMP X-ray crystal structures. Results Lower metabolism of BHTOOH and BMPOOH was observed in the WT P450 BM3 and the T268A P450 BM3 mutant than in the F87A and F87A/T268A P450 BM3 mutants. Large differences were found in the F–G loop regions and active site cavity volumes for the F87A mutated structures. Conclusions Analysis of the metabolism, X-ray crystal structures, and molecular docking simulations suggests that P450 BM3 activity toward BHTOOH and BMPOOH is mediated through substrate recognition by T268 and F87, and the active site cavity volume. Based on this information, a simplified representation is presented with the relative orientation of organic hydroperoxides in the P450 BM3 active site. General significance The metabolism results and structural analysis of this model P450 allowed us to rationalize the structural factors that influence organic hydroperoxide metabolism.

ATP binding cassette type-C (ABCC) transporters move molecules across cell membranes upon hydrolysis of ATP; however, their coupling of ATP hydrolysis to substrate transport remains elusive. Drosophila multidrug resistance-associated protein (DMRP) is the functional ortholog of human long ABCC transporters, with similar substrate and inhibitor specificity, but higher activity. Exploiting its high activity, we kinetically dissected the catalytic mechanism of DMRP by using E2-d-glucuronide (E2G), the physiologic substrate of human ABCC. We examined the DMRP-mediated interdependence of ATP and E2G in biochemical assays. We observed E2G-dependent ATPase activity to be biphasic at subsaturating ATP concentrations, which implies at least 2 E2G binding sites on DMRP. Furthermore, transport measurements indicated strong nonreciprocal cooperativity between ATP and E2G. In addition to confirming these findings, our kinetic modeling with the Complex Pathway Simulator indicated a 10-fold decrease in the E2G-mediated activation of ATP hydrolysis upon saturation of the second E2G binding site. Surprisingly, the binding of the second E2G allowed for substrate transport with a constant rate, which tightly coupled ATP hydrolysis to transport. In summary, we show that the second E2G binding-similar to human ABCC2-allosterically stimulates transport activity of DMRP. Our data suggest that this is achieved by a significant increase in the coupling of ATP hydrolysis to transport.-Karasik, A., Ledwitch, K. V., Aranyi, T., Varadi, A., Roberts, A., Szeri, F. Boosted coupling of ATP hydrolysis to substrate transport upon cooperative estradiol-17-β-D-glucuronide binding in a Drosophila ATP binding cassette type-C transporter.

P-glycoprotein (Pgp) is an ATP-binding cassette (ABC) transporter that plays a major role in cardiovascular drug disposition by effluxing a chemically and structurally diverse range of cardiovascular therapeutics. Unfortunately, drug-drug interactions (DDIs) with the transporter have become a major roadblock to effective cardiovascular drug administration because they can cause adverse drug reactions (ADRs) or reduce the efficacy of drugs. Cardiovascular ion channel inhibitors are particularly susceptible to DDIs and ADRs with Pgp because they often have low therapeutic indexes and are commonly coadministered with other drugs that are also Pgp substrates. DDIs from cardiovascular ion channel inhibitors with the transporter occur because of inhibition or induction of the transporter and the transporter’s tissue and cellular localization. Inhibiting Pgp can increase absorption and reduce excretion of drugs, leading to elevated drug plasma concentrations and drug toxicity. In contrast, inducing Pgp can have the opposite effect by reducing the drug plasma concentration and its efficacy. A number of in vitro and in vivo studies have already demonstrated DDIs from several cardiovascular ion channel inhibitors with human Pgp and its animal analogs, including verapamil, digoxin, and amiodarone. In this review, Pgp-mediated DDIs and their effects on pharmacokinetics for different categories of cardiovascular ion channel inhibitors are discussed. This information is essential for improving pharmacokinetic predictions of cardiovascular therapeutics, for safer cardiovascular drug administration and for mitigating ADRs emanating from Pgp.

In medicinal chemistry, one of the most studied molecules in recent history is taxol. Taxol is a versatile natural product that is used in various cancer treatment regimens. It is administered to patients with breast, lung, and ovarian cancers, and is currently being studied for the treatment of squamous cell carcinoma of the oral cavity and tongue. Taxol has been tested in a number of research and clinical phase trials to determine feasibility, toxicity, and cytotoxicity against oral squamous cell carcinoma as a single drug regimen and as a contributing drug component in treatment plans. This paper reviews over forty articles that examine cell lines, murine models, and human results for the response of taxol against squamous cell carcinoma (SCC) of the oral cavity and the tongue.

The P-glycoprotein (Pgp) transporter plays a central role in drug disposition by effluxing a chemically diverse range of drugs from cells through conformational changes and ATP hydrolysis. A number of drugs are known to activate ATP hydrolysis of Pgp, but coupling between ATP and drug binding is not well understood. The cardiovascular drug verapamil is one of the most widely studied Pgp substrates and therefore, represents an ideal drug to investigate the drug-induced ATPase activation of Pgp. As previously noted, verapamil-induced Pgp-mediated ATP hydrolysis kinetics was biphasic at saturating ATP concentrations. However, at subsaturating ATP concentrations, verapamil-induced ATPase activation kinetics became monophasic. To further understand this switch in kinetic behavior, the Pgp-coupled ATPase activity kinetics was checked with a panel of verapamil and ATP concentrations and fit with the substrate inhibition equation and the kinetic fitting software COPASI. The fits suggested that cooperativity between ATP and verapamil switched between low and high verapamil concentration. Fluorescence spectroscopy of Pgp revealed that cooperativity between verapamil and a non-hydrolyzable ATP analog leads to distinct global conformational changes of Pgp. NMR of Pgp reconstituted in liposomes showed that cooperativity between verapamil and the non-hydrolyzable ATP analog modulate each other's interactions. This information was used to produce a conformationally-gated model of drug-induced activation of Pgp-mediated ATP hydrolysis.

P-glycoprotein (Pgp) plays a major role in promoting drug–drug interactions (DDIs) with verapamil and digoxin. In the present study, we present a comprehensive molecular and mechanistic model of Pgp DDIs encompassing drug binding, ATP hydrolysis, transport and conformational changes., Drug–drug interactions (DDIs) and associated toxicity from cardiovascular drugs represents a major problem for effective co-administration of cardiovascular therapeutics. A significant amount of drug toxicity from DDIs occurs because of drug interactions and multiple cardiovascular drug binding to the efflux transporter P-glycoprotein (Pgp), which is particularly problematic for cardiovascular drugs because of their relatively low therapeutic indexes. The calcium channel antagonist, verapamil and the cardiac glycoside, digoxin, exhibit DDIs with Pgp through non-competitive inhibition of digoxin transport, which leads to elevated digoxin plasma concentrations and digoxin toxicity. In the present study, verapamil-induced ATPase activation kinetics were biphasic implying at least two verapamil-binding sites on Pgp, whereas monophasic digoxin activation of Pgp-coupled ATPase kinetics suggested a single digoxin-binding site. Using intrinsic protein fluorescence and the saturation transfer double difference (STDD) NMR techniques to probe drug–Pgp interactions, verapamil was found to have little effect on digoxin–Pgp interactions at low concentrations of verapamil, which is consistent with simultaneous binding of the drugs and non-competitive inhibition. Higher concentrations of verapamil caused significant disruption of digoxin–Pgp interactions that suggested overlapping and competing drug-binding sites. These interactions correlated to drug-induced conformational changes deduced from acrylamide quenching of Pgp tryptophan fluorescence. Also, Pgp-coupled ATPase activity kinetics measured with a range of verapamil and digoxin concentrations fit well to a DDI model encompassing non-competitive and competitive inhibition of digoxin by verapamil. The results and previous transport studies were combined into a comprehensive model of verapamil–digoxin DDIs encompassing drug binding, ATP hydrolysis, transport and conformational changes.

A number of delivery agents, such as proteins, liposomes, micelles, and nanoparticles, are utilized for transporting pharmaceutical agents in a physiological environment. This Letter focuses on the use of the copper(II) ion and its potential role as a delivery agent for the taxanes and taxol couple to a malaria drug. Nuclear magnetic resonance (NMR, (1)H, (13)C, (15)N), Mass Spectrometry (LC-MS, MALDI-TOF, FT-ICR) and computational methods are used to examine the structure of the complex. The National Cancer Institute's benchmark 60 cell line panel is used to compare the efficacy of the copper-taxol and copper-taxol-hydroxychloroquin complexes to that of iron-taxol and pure taxol.