Publications

BACKGROUND: The nuclear factor erythroid 2-related factor 2 (Nrf2) is a potential molecular target for cancer chemoprevention. Si-Wu-Tang (SWT), a popular traditional Chinese medicine for women's health, was reported with a novel activity of cancer prevention. PURPOSE: The present study was aimed to identify the bioactive constituents in SWT responsible for the Nrf2 activating and cancer preventive activity and explore the pharmacological mechanisms. METHODS: Nine compounds detectable from various batches of SWT were ranked using in silico molecular docking based on their ability to interfere the forming of Nrf2-Keap1 complex. The predicted Nrf2 activating effect was validated using the antioxidant response element (ARE) luciferase reporter assay and quantitative RT-PCR analysis for select Nrf2 regulated genes Hmox1, Nqo1 and Slc7a11. The antimutagenic activity of the compounds were determined by the Ames test. The chemopreventive activity of these compounds were assessed on EGF-induced neoplastic transformation of JB6 P+ cells, an established non-cancerous murine epidermal model for studying tumor promotion and identifying cancer preventive agents. These compounds were further characterized using luciferase reporter assay on EGF-induced activation of AP-1, a known transcription factor mediating carcinogenesis. RESULTS: Three of the nine compounds predicted as Nrf2 activators by molecular docking, gallic acid (GA), Z-liguistilide (LIG), and senkyunolide A (SA), were confirmed with highest potency of increasing the Nrf2/ARE promoter activity and upregulating the expression of Hmox1, Nqo1 and Slc7a11. In addition, GA, LIG and SA exhibited an antimutagenic activity against the direct mutagen 2-nitrofluorene while no mutagenic effects were observed at the same time in Ames test. At nontoxic concentrations, GA, LIG, and SA inhibited EGF-induced neoplastic transformation of JB6 P+ cells. Combined treatment of GA, LIG and SA, in the same ratio as detected in SWT, showed enhanced effect against JB6 transformation compared with that of the single compound alone. GA, LIG and SA, alone or in combination, suppressed EGF-induced activation of AP-1. CONCLUSION: We identified three bioactive constituents in SWT responsible for the Nrf2 activating and cancer preventive activity. This study provides evidence supporting novel molecular basis of SWT in cancer prevention.

Piezo proteins are transmembrane ion channels which transduce many forms of mechanical stimuli into electrochemical signals. Their pore, formed by the assembly of three identical subunits, opens by an unknown mechanism. Here, to probe this mechanism, we investigate the interaction of Piezo1 with the small molecule agonist Yoda1. By engineering chimeras between mouse Piezo1 and its Yoda1-insensitive paralog Piezo2, we first identify a minimal protein region required for Yoda1 sensitivity. We next study the effect of Yoda1 on heterotrimeric Piezo1 channels harboring wild type subunits and Yoda1-insensitive mutant subunits. Using calcium imaging and patch-clamp electrophysiology, we show that hybrid channels harboring as few as one Yoda1-sensitive subunit exhibit Yoda1 sensitivity undistinguishable from homotrimeric wild type channels. Our results show that the Piezo1 pore remains fully open if only one subunit remains activated. This study sheds light on the gating and pharmacological mechanisms of a member of the Piezo channel family.

Type 1 Serine/Threonine Kinase Receptors (STKR1) transduce a wide spectrum of biological signals mediated by TGF-beta superfamily members. The STKR1 activity is tightly controlled by their regulatory glycine-serine rich (GS) domain adjacent to the kinase domain. Despite decades of studies, it remains unknown how physiological or pathological GS domain modifications are coupled to STKR1 kinase activity. Here, by performing molecular dynamics simulations and free energy calculation of Activin-Like Kinase 2 (ALK2), we found that GS domain phosphorylation, FKBP12 dissociation, and disease mutations all destabilize a D354-R375 salt-bridge, which normally acts as an electrostatic lock to prevent coordination of adenosine triphosphate (ATP) to the catalytic site. We developed a WAFEX-guided principal analysis and unraveled how phosphorylation destabilizes this highly conserved salt-bridge in temporal and physical space. Using current-flow betweenness scores, we identified an allosteric network of residue-residue contacts between the GS domain and the catalytic site that controls the formation and disruption of this salt bridge. Importantly, our novel network analysis approach revealed how certain disease-causing mutations bypass FKBP12-mediated kinase inhibition to produce leaky signaling. We further provide experimental evidence that this salt-bridge lock exists in other STKR1s, and acts as a general safety mechanism in STKR1 to prevent pathological leaky signaling. In summary, our study provides a compelling and unifying allosteric activation mechanism in STKR1 kinases that reconciles a large number of experimental studies and sheds light on a novel therapeutic avenue to target disease-related STKR1 mutants.

Reversible covalent inhibitors have many clinical advantages over noncovalent or irreversible covalent drugs. However, apart from selecting a warhead, substantial efforts in design and synthesis are needed to optimize noncovalent interactions to improve target-selective binding. Computational prediction of binding affinity for reversible covalent inhibitors presents a unique challenge since the binding process consists of multiple steps, which are not necessarily independent of each other. In this study, we lay out the relation between relative binding free energy and the overall reversible covalent binding affinity using a two-state binding model. To prove the concept, we employed free energy perturbation (FEP) coupled with lambda-exchange molecular dynamics method to calculate the binding free energy of a series of alpha-ketoamide analogues relative to a common warhead scaffold, in both noncovalent and covalent binding states, and for two highly homologous proteases, calpain-1 and calpain-2. We conclude that covalent binding state alone, in general, can be used to predict reversible covalent binding selectivity. However, exceptions may exist. Therefore, we also discuss the conditions under which the noncovalent binding step is no longer negligible and propose to combine the relative FEP calculations with a single QM/MM calculation of warhead to predict the binding affinity and binding kinetics. Our FEP calculations also revealed that covalent and noncovalent binding states of an inhibitor do not necessarily exhibit the same selectivity. Thus, investigating both binding states, as well as the kinetics will provide extremely useful information for optimizing reversible covalent inhibitors.

Aberrant opening of nonjunctional connexin hemichannels at the plasma membrane is associated with many diseases, including ischemia and muscular dystrophy. Proper control of hemichannel opening is essential to maintain cell viability and is achieved by physiological levels of extracellular Ca(2+), which drastically reduce hemichannel activity. Here we examined the role of conserved charged residues that form electrostatic networks near the extracellular entrance of the connexin pore, a region thought to be involved in gating rearrangements of hemichannels. Molecular dynamics simulations indicate discrete sites for Ca(2+) interaction and consequent disruption of salt bridges in the open hemichannels. Experimentally, we found that disruption of these salt bridges by mutations facilitates hemichannel closing. Two negatively charged residues in these networks are putative Ca(2+) binding sites, forming a Ca(2+)-gating ring near the extracellular entrance of the pore. Accessibility studies showed that this Ca(2+)-bound gating ring does not prevent access of ions or small molecules to positions deeper into the pore, indicating that the physical gate is below the Ca(2+)-gating ring. We conclude that intra- and intersubunit electrostatic networks at the extracellular entrance of the hemichannel pore play critical roles in hemichannel gating reactions and are tightly controlled by extracellular Ca(2+) Our findings provide a general mechanism for Ca(2+) gating among different connexin hemichannel isoforms.

A signal property of connexin channels is the ability to mediate selective diffusive movement of molecules through plasma membrane(s), but the energetics and determinants of molecular movement through these channels have yet to be understood. Different connexin channels have distinct molecular selectivities that cannot be explained simply on the basis of size or charge of the permeants. To gain insight into the forces and interactions that underlie selective molecular permeation, we investigated the energetics of two uncharged derivatized sugars, one permeable and one impermeable, through a validated connexin26 (Cx26) channel structural model, using molecular dynamics and associated analytic tools. The system is a Cx26 channel equilibrated in explicit membrane/solvent, shown by Brownian dynamics to reproduce key conductance characteristics of the native channel. The results are consistent with the known difference in permeability to each molecule. The energetic barriers extend through most of the pore length, rather than being highly localized as in ion-specific channels. There is little evidence for binding within the pore. Force decomposition reveals how, for each tested molecule, interactions with water and the Cx26 protein vary over the length of the pore and reveals a significant contribution from hydrogen bonding and interaction with K(+). The flexibility of the pore width varies along its length, and the tested molecules have differential effects on pore width as they pass through. Potential sites of interaction within the pore are defined for each molecule. The results suggest that for the tested molecules, differences in hydrogen bonding and entropic factors arising from permeant flexibility substantially contribute to the energetics of permeation. This work highlights factors involved in selective molecular permeation that differ from those that define selectivity among atomic ions.

Gated solely by activity-induced changes in intracellular calcium, small-conductance potassium channels (SKs) are critical for a variety of functions in the CNS, from learning and memory to rhythmic activity and sleep. While there is a wealth of information on SK2 gating, kinetics, and Ca(2+) sensitivity, little is known regarding the regulation of SK2 subcellular localization. We report here that synaptic SK2 levels are regulated by the E3 ubiquitin ligase UBE3A, whose deficiency results in Angelman syndrome and overexpression in increased risk of autistic spectrum disorder. UBE3A directly ubiquitinates SK2 in the C-terminal domain, which facilitates endocytosis. In UBE3A-deficient mice, increased postsynaptic SK2 levels result in decreased NMDA receptor activation, thereby impairing hippocampal long-term synaptic plasticity. Impairments in both synaptic plasticity and fear conditioning memory in UBE3A-deficient mice are significantly ameliorated by blocking SK2. These results elucidate a mechanism by which UBE3A directly influences cognitive function.

Abnormal alteration of bone morphogenetic protein (BMP) signaling is implicated in many types of diseases including cancer and heterotopic ossifications. Hence, small molecules targeting BMP type I receptors (BMPRI) to interrupt BMP signaling are believed to be an effective approach to treat these diseases. However, lack of understanding of the molecular determinants responsible for the binding selectivity of current BMP inhibitors has been a big hindrance to the development of BMP inhibitors for clinical use. To address this issue, we carried out in silico experiments to test whether computational methods can reproduce and explain the high selectivity of a small molecule BMP inhibitor DMH1 on BMPRI kinase ALK2 vs. the closely related TGF-beta type I receptor kinase ALK5 and vascular endothelial growth factor receptor type 2 (VEGFR2) tyrosine kinase. We found that, while the rigid docking method used here gave nearly identical binding affinity scores among the three kinases; free energy perturbation coupled with Hamiltonian replica-exchange molecular dynamics (FEP/H-REMD) simulations reproduced the absolute binding free energies in excellent agreement with experimental data. Furthermore, the binding poses identified by FEP/H-REMD led to a quantitative analysis of physical/chemical determinants governing DMH1 selectivity. The current work illustrates that small changes in the binding site residue type (e.g. pre-hinge region in ALK2 vs. ALK5) or side chain orientation (e.g. Tyr219 in caALK2 vs. wtALK2), as well as a subtle structural modification on the ligand (e.g. DMH1 vs. LDN193189) will cause distinct binding profiles and selectivity among BMP inhibitors. Therefore, the current computational approach represents a new way of investigating BMP inhibitors. Our results provide critical information for designing exclusively selective BMP inhibitors for the development of effective pharmacotherapy for diseases caused by aberrant BMP signaling.

Computational methodologies that couple the dynamical evolution of a set of replicated copies of a system of interest offer powerful and flexible approaches to characterize complex molecular processes. Such multiple copy algorithms (MCAs) can be used to enhance sampling, compute reversible work and free energies, as well as refine transition pathways. Widely used examples of MCAs include temperature and Hamiltonian-tempering replica-exchange molecular dynamics (T-REMD and H-REMD), alchemical free energy perturbation with lambda replica-exchange (FEP/lambda-REMD), umbrella sampling with Hamiltonian replica exchange (US/H-REMD), and string method with swarms-of-trajectories conformational transition pathways. Here, we report a robust and general implementation of MCAs for molecular dynamics (MD) simulations in the highly scalable program NAMD built upon the parallel programming system Charm++. Multiple concurrent NAMD instances are launched with internal partitions of Charm++ and located continuously within a single communication world. Messages between NAMD instances are passed by low-level point-to-point communication functions, which are accessible through NAMD's Tcl scripting interface. The communication-enabled Tcl scripting provides a sustainable application interface for end users to realize generalized MCAs without modifying the source code. Illustrative applications of MCAs with fine-grained inter-copy communication structure, including global lambda exchange in FEP/lambda-REMD, window swapping US/H-REMD in multidimensional order parameter space, and string method with swarms-of-trajectories were carried out on IBM Blue Gene/Q to demonstrate the versatility and massive scalability of the present implementation. Published by Elsevier B.V.

Overbinding of ions to lipid head groups is a potentially serious artifact in simulations of charged lipid bilayers. In this study, the Lennard-Jones radii in the CHARMM force field for interactions of Na+ and lipid oxygen atoms of carboxyl, phosphate, and ester groups were revised to match osmotic pressure data on sodium acetate and electrophoresis data on palmitoyloleoyl phosphatidylcholine (POPC) vesicles. The new parameters were then validated by successfully reproducing previously published experimental NMR deuterium order parameters for dimyristoyl phosphatidylglycerol (DMPG) and newly obtained values for palmitoyloleoyl phosphatidylserine (POPS). Although the increases in Lennard-Jones diameters are only 0.02-0.12 angstrom, they are sufficient of reduce Na+ binding, and thereby increase surface areas per lipid by 5-10% compared with the unmodified parameters.