Structural insight into tetrameric hTRPV1 from homology modeling, molecular docking, molecular dynamics simulation, virtual screening, and bioassay validations

Parameterization of General Organic Polymers within the Open Force Field Framework

Polymer and chemically modified biopolymer systems present unique challenges to traditional molecular simulation preparation workflows. First, typical polymer and biomolecular input formats, such as Protein Data Bank (PDB) files, lack adequate chemical information needed for the parameterization of new chemistries. Second, polymers are typically too large for accurate partial charge generation methods. In this work, we employ direct chemical perception through the Open Force Field toolkit to create a flexible polymer simulation workflow for organic polymers, encompassing everything from biopolymers to soft materials. We propose and test a new input specification for monomer information that can, along with a 3D conformational geometry, parametrize and simulate most soft-material systems within the same workflow used for smaller ligands. The monomer format encompasses a subset of the SMIRKS substructure query language to uniquely identify chemical information and repeating charges in underspecified systems through matching atomic connectivity. This workflow is combined with several different approaches for automatic partial-charge generation for larger systems. As an initial proof of concept, a variety of diverse polymeric systems were parametrized with the Open Force Field toolkit, including functionalized proteins, DNA, homopolymers, cross-linked systems, and sugars. Additionally, shape properties and radial distribution functions were computed from molecular dynamics simulations of poly(ethylene glycol), polyacrylamide, and poly(N-isopropylacrylamide) homopolymers in aqueous solution and compared to previous simulation results in order to demonstrate a start-to-finish workflow for simulation and property prediction. We expect that these tools will greatly expedite the day-to-day computational research of soft-matter simulations and create a robust atomic-scale polymer specification in conjunction with existing polymer structural notations.

Insights from molecular dynamics simulations of TRPV1 channel modulators in pain

TRPV1 is a nonselective cation channel vital for detecting noxious stimuli (heat, acid, capsaicin). Its role in pain makes it a potential drug target for chronic pain management, migraines, and related disorders. This review updates molecular dynamics (MD) simulation studies on the TRPV1 channel, focusing on its gating mechanism, ligand-binding sites, and implications for drug design. The article also explores challenges in developing modulators, SAR optimization, and clinical trial studies. Efforts have been undertaken to concisely present MD simulation findings, with a focus on their relevance to drug discovery.

Identification of a Partial and Selective TRPV1 Agonist CPIPC for Alleviation of Inflammatory Pain

Transient receptor potential vanilloid 1 (TRPV1) is a non-selective cation channel, predominantly expressed in a subset of peripheral sensory neurons for pain signaling. Topical application of agonist capsaicin for desensitizing TRPV1 currents has been approved for relief of chronic pain. However, the potent TRPV1 capsaicin is not ingestible and even topical capsaicin causes common side effects such as skin irritation, swelling, erythema and pruritus, suggesting that a mild TRPV1 agonist might be helpful for reducing side effects while reliving pain. In this study, we reported on a partial and selective TRPV1 agonist 4-(5-chloropyridin-2-yl)-N-(1H-indazol-6-yl)piperazine-1-carboxamide named CPIPC that was modified based on targeting the residue Arg557, important for conversion between the channel antagonism and agonism. Whole-cell patch clamp recordings indicated a concentration-dependent activation of TRPV1 currents by CPIPC with an EC₅₀ of 1.56 ± 0.13 μM. The maximum efficacy of CPIPC (30 μM) was about 60% of saturated capsaicin (10 μM). Repetitive additions of CPIPC caused TRPV1 current desensitization in both TRPV1-expressing HEK293 cells and dorsal root ganglion (DRG) sensory neurons. Oral administration of CPIPC dose-dependently alleviated inflammatory pain in mice. Further site-directed mutagenesis combined with molecular docking revealed that residue Arg557 is critical for TRPV1 activation by CPIPC. Taken together, we identified a novel partial and selective TRPV1 agonist CPIPC that exhibits antinociceptive activity in mice.

Discovery of Potent and Selective Transient Receptor Potential Vanilloid 1 (TRPV1) Agonists with Analgesic Effects In Vivo Based on the Functional Conversion Induced by Altering the Orientation of the Indazole Core

Transient receptor potential vanilloid 1 (TRPV1) is a promising target for developing antinociceptive agents. Here, we report the synthesis of N-indazole-4-aryl piperazine carboxamide analogues as TRPV1 modulators. The structure-activity relationship (SAR) reveals that substituting indazole at the 5-/6-position leads to TRPV1 agonism, whereas the 4- and 7-positions of indazole obtain mild antagonism and loss of activity, respectively. The whole-cell clamp patch assay shows that 28 is a potent and selective TRPV1 agonist and it relieves inflammatory and thermal pain by desensitizing the native TRPV1 current in the dorsal root ganglion (DRG) in mice. Additionally, site-directed mutagenesis combined with molecular docking shows an important hydrogen interaction between Arg557 and the indazole of 28. Taken together, our findings provide insight into TRPV1 agonism-antagonism conversion based on the interaction between indazole and Arg557, which provides a strategy to obtain new TRPV1 agonists by structural modification of antagonists. Compound 28 may be used as a lead compound for further optimization.

The antimicrobial cyclic peptide B2 combats multidrug resistant Acinetobacter baumannii infection

In silico methods were employed for the development of antimicrobial peptides against MDR A. baumannii by binding to BamA.

A molecular perspective on identifying TRPV1 thermosensitive regions and disentangling polymodal activation

TRPV1 is a polymodal receptor ion channel that is best known to function as a molecular thermometer. It is activated in diverse ways, including by heat, protons (low pH), and vanilloid compounds, such as capsaicin. In this review, we summarize molecular studies of TRPV1 thermosensing, focusing on the cross-talk between heat and other activation modes. Additional insights from TRPV1 isoforms and non-rodent/non-human TRPV1 ortholog studies are also discussed in this context. While the molecular mechanism of heat activation is still emerging, it is clear that TRPV1 thermosensing is modulated allosterically, i.e., at a distance, with contributions from many distinct regions of the channel. Similarly, current studies identify cross-talk between heat and other TRPV1 activation modes, such as protons and capsaicin, and that these modes can generally be selectively disentangled. In aggregate, this suggests that future TRPV1 molecular studies should define allosteric pathways and provide mechanistic insight, thereby enabling mode-selective manipulation of the polymodal receptor. These advances are anticipated to have significant implications in both basic and applied biomedical sciences.

Discovery of modulators for the PD-1/PD-L1 interaction by molecular simulation and bioassay

PD-1-targeted discovery of modulators by virtual screening, molecular docking, surface plasmon resonance binding assay and T-cell activation assay.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Molecular docking studies of tea (Thea sinensis Linn.) polyphenols inhibition pattern with Rat P-glycoprotein

Since 3000 B.C., evergreen plant Thea sinensis (Theaceae) is used both as a social and medicinal beverage. Leaves of T. sinensis contain amino acids, vitamins, caffeine, polysaccharides and polyphenols. Most of the natural medicinal actions of tea are due to the availability and abundance of polyphenols mainly catechins. It has also been stated that some catechins were absorbed more rapidly than other compounds after the oral administration of tea and could increase the bio-enhancing activities of anticancer drugs by inhibiting P-glycoprotein (P-gp). The results of the molecular docking showed that polyphenols bind easily to the active P-gp site. All compounds exhibited fluctuating binding affinity ranged from −11.67 to −8.36 kcal/mol. Observed binding energy required for theaflavin to bind to P-gp was lowest (−11.67 kcal/mol). The obtained data that supports all the selected polyphenols inhibited P-gp and therefore may enhance the bioavailability of drugs. This study may play a vital role in finding hotspots in P-gp and eventually may be proved useful in designing compounds with high affinity and specificity to the protein.

A single TRPV1 amino acid controls species sensitivity to capsaicin

Chili peppers produce capsaicin (a vanilloid) that activates the transient receptor potential cation channel subfamily V member 1 (TRPV1) on sensory neurons to alter their membrane potential and induce pain. To identify residues responsible for differential TRPV1 capsaicin sensitivity among species, we used intracellular Ca²⁺ imaging to characterize chimeras composed of capsaicin-sensitive rat TRPV1 (rTRPV1) and capsaicin-insensitive chicken TRPV1 (cTRPV1) exposed to a series of capsaicinoids. We found that chimeras containing rat E570-V686 swapped into chicken receptors displayed capsaicin sensitivity, and that simply changing the alanine at position 578 in the S4-S5 helix of the chicken receptor to a glutamic acid was sufficient to endow it with capsaicin sensitivity in the micromolar range. Moreover, introduction of lysine, glutamine or proline at residue A578 also elicited capsaicin sensitivity in cTRPV1. Similarly, replacing corresponding rTRPV1 residue E570 with lysine or glutamine retained capsaicin sensitivity. The hydrophilic capsaicin analog Cap-EA activated a cTRPV1-A578E mutant, suggesting that A578 may participate in vanilloid binding. The hydrophilic vanilloid agonist zingerone did not activate any A578 mutants with capsaicin sensitivity, suggesting that the vanilloid group alone is not sufficient for receptor activation. Our study demonstrates that a subtle modification of TRPV1 in different species globally alters capsaicin responses.

Structure-Based Design of Novel Biphenyl Amide Antagonists of Human Transient Receptor Potential Cation Channel Subfamily M Member 8 Channels with Potential Implications in the Treatment of Sensory Neuropathies

Structure-activity relationship studies of a reported menthol-based transient receptor potential cation channel subfamily M member 8 channel (TRPM8) antagonist, guided by computational simulations and structure-based design, uncovers a novel series of TRPM8 antagonists with >10-fold selectivity versus related TRP subtypes. Spiro[4.5]decan-8-yl analogue 14 inhibits icilin-evoked Ca²⁺ entry in HEK-293 cells stably expressing human TRPM8 (hTRPM8) with an IC₅₀ of 2.4 ± 1.0 nM, while in whole-cell patch-clamp recordings this analogue inhibits menthol-evoked currents with a hTRPM8 IC₅₀ of 64 ± 2 nM. Molecular dynamics (MD) simulations of compound 14 in our homology model of hTRPM8 suggest that this antagonist forms extensive hydrophobic contacts within the orthosteric site. In the wet dog shakes (WDS) assay, compound 14 dose-dependently blocks icilin-triggered shaking behaviors in mice. Upon local administration, compound 14 dose dependently inhibits cold allodynia evoked by the chemotherapy oxaliplatin in a murine model of peripheral neuropathy at microgram doses. Our findings suggest that 14 and other biphenyl amide analogues within our series can find utility as potent antagonist chemical probes derived from (-)-menthol as well as small molecule therapeutic scaffolds for chemotherapy-induced peripheral neuropathy (CIPN) and other sensory neuropathies.

Cryo-EM Structure of the Human Cannabinoid Receptor CB2-Gi Signaling Complex

Drugs selectively targeting CB2 hold promise for treating neurodegenerative disorders, inflammation, and pain while avoiding psychotropic side effects mediated by CB1. The mechanisms underlying CB2 activation and signaling are poorly understood but critical for drug design. Here we report the cryo-EM structure of the human CB2-G_i signaling complex bound to the agonist WIN 55,212-2. The 3D structure reveals the binding mode of WIN 55,212-2 and structural determinants for distinguishing CB2 agonists from antagonists, which are supported by a pair of rationally designed agonist and antagonist. Further structural analyses with computational docking results uncover the differences between CB2 and CB1 in receptor activation, ligand recognition, and G_i coupling. These findings are expected to facilitate rational structure-based discovery of drugs targeting the cannabinoid system.

Combined approach of homology modeling, molecular dynamics, and docking: computer-aided drug discovery

With the continuous development in software, algorithms, and increase in computer speed, the field of computer-aided drug design has been witnessing reduction in the time and cost of the drug designing process. Structure based drug design (SBDD), which is based on the 3D structure of the enzyme, is helping in proposing novel inhibitors. Although a number of crystal structures are available in various repositories, there are various proteins whose experimental crystallization is difficult. In such cases, homology modeling, along with the combined application of MD and docking, helps in establishing a reliable 3D structure that can be used for SBDD. In this review, we have reported recent works, which have employed these three techniques for generating structures and further proposing novel inhibitors, for cytoplasmic proteins, membrane proteins, and metal containing proteins. Also, we have discussed these techniques in brief in terms of the theory involved and the various software employed. Hence, this review can give a brief idea about using these tools specifically for a particular problem.

PD-1-Targeted Discovery of Peptide Inhibitors by Virtual Screening, Molecular Dynamics Simulation, and Surface Plasmon Resonance

The blockade of the programmed cell death protein 1/programmed cell death ligand 1 (PD-1/PD-L1) pathway plays a critical role in cancer immunotherapy by reducing the immune escape. Five monoclonal antibodies that antagonized PD-1/PD-L1 interaction have been approved by the Food and Drug Administration (FDA) and marketed as immunotherapy for cancer treatment. However, some weaknesses of antibodies, such as high cost, low stability, poor amenability for oral administration, and immunogenicity, should not be overlooked. To overcome these disadvantages, small-molecule inhibitors targeting PD-L1 were developed. In the present work, we applied in silico and in vitro approaches to develop short peptides targeting PD-1 as chemical probes for the inhibition of PD-1-PD-L1 interaction. We first predicted the potential binding pocket on PD-1/PD-L1 protein-protein interface (PPI). Sequentially, we carried out virtual screening against our in-house peptide library to identify potential ligands. WANG-003, WANG-004, and WANG-005, three of our in-house peptides, were predicted to bind to PD-1 with promising docking scores. Next, we conducted molecular docking and molecular dynamics (MD) simulation for the further analysis of interactions between our peptides and PD-1. Finally, we evaluated the affinity between peptides and PD-1 by surface plasmon resonance (SPR) binding technology. The present study provides a new perspective for the development of PD-1 inhibitors that disrupt PD-1-PD-L1 interactions. These promising peptides have the potential to be utilized as a novel chemical probe for further studies, as well as providing a foundation for further designs of potent small-molecule inhibitors targeting PD-1.

Trilobatin, a Novel SGLT1/2 Inhibitor, Selectively Induces the Proliferation of Human Hepatoblastoma Cells

Studies have indicated that Na⁺-d-glucose co-transporter (SGLT) inhibitors had anti-proliferative activity by attenuating the uptake of glucose in several tumor cell lines. In this study, the molecular docking showed that, trilobatin, one of the dihydrochalcones from leaves of Lithocarpus polystachyus Rehd., might be a novel inhibitor of SGLT1 and SGLT2, which evidently attenuated the uptake of glucose in vitro and in vivo. To our surprise, we observed that trilobatin did not inhibit, but promoted the proliferation of human hepatoblastoma HepG2 and Huh 7 cells when it was present at high concentrations. At the same time, incubation with high concentrations of trilobatin arrested the cell cycle at S phase in HepG2 cells. We also found that treatment with trilobatin had no significant effect on the expression of hepatitis B x-interacting protein (HBXIP) and hepatocyte nuclear factor (HNF)-4α, the two key regulators of hepatocyte proliferation. Taken together, although trilobatin worked as a novel inhibitor of SGLTs to attenuate the uptake of glucose, it also selectively induced the cell proliferation of HepG2 cells, suggesting that not all the SGLT inhibitors inhibited the proliferation of tumor cells, and further studies are needed to assess the anti-cancer potentials of new glucose-lowering agents.

Key Topics in Molecular Docking for Drug Design

Molecular docking has been widely employed as a fast and inexpensive technique in the past decades, both in academic and industrial settings. Although this discipline has now had enough time to consolidate, many aspects remain challenging and there is still not a straightforward and accurate route to readily pinpoint true ligands among a set of molecules, nor to identify with precision the correct ligand conformation within the binding pocket of a given target molecule. Nevertheless, new approaches continue to be developed and the volume of published works grows at a rapid pace. In this review, we present an overview of the method and attempt to summarise recent developments regarding four main aspects of molecular docking approaches: (i) the available benchmarking sets, highlighting their advantages and caveats, (ii) the advances in consensus methods, (iii) recent algorithms and applications using fragment-based approaches, and (iv) the use of machine learning algorithms in molecular docking. These recent developments incrementally contribute to an increase in accuracy and are expected, given time, and together with advances in computing power and hardware capability, to eventually accomplish the full potential of this area.

Combination of a Rapidly Penetrating Agonist and a Slowly Penetrating Antagonist Affords Agonist Action of Limited Duration at the Cellular Level

The capsaicin receptor TRPV1 (transient receptor potential vanilloid 1) has been an object of intense interest for pharmacological development on account of its critical role in nociception. In the course of structure activity analysis, it has become apparent that TRPV1 ligands may vary dramatically in the rates at which they interact with TRPV1, presumably reflecting differences in their abilities to penetrate into the cell. Using a fast penetrating agonist together with an excess of a slower penetrating antagonist, we find that we can induce an agonist response of limited duration and, moreover, the duration of the agonist response remains largely independent of the absolute dose of agonist, as long as the ratio of antagonist to agonist is held constant. This general approach for limiting agonist duration under conditions in which absolute agonist dose is variable should have more general applicability.

Computational Systems Pharmacology-Target Mapping for Fentanyl-Laced Cocaine Overdose

The United States of America is fighting against one of its worst-ever drug crises. Over 900 people a week die from opioid- or heroin-related overdoses, while millions more suffer from opioid prescription addiction. Recently, drug overdoses caused by fentanyl-laced cocaine specifically are on the rise. Due to drug synergy and an increase in side effects, polydrug addiction can cause more risk than addiction to a single drug. In the present work, we systematically analyzed the overdose and addiction mechanism of cocaine and fentanyl. First, we applied our established chemogenomics knowledgebase and machine-learning-based methods to map out the potential and known proteins, transporters, and metabolic enzymes and the potential therapeutic target(s) for cocaine and fentanyl. Sequentially, we looked into the detail of (1) the addiction to cocaine and fentanyl by binding to the dopamine transporter and the μ opioid receptor (DAT and μOR, respectively), (2) the potential drug-drug interaction of cocaine and fentanyl via p-glycoprotein (P-gp) efflux, (3) the metabolism of cocaine and fentanyl in CYP3A4, and (4) the physiologically based pharmacokinetic (PBPK) model for two drugs and their drug-drug interaction at the absorption, distribution, metabolism, and excretion (ADME) level. Finally, we looked into the detail of JWH133, an agonist of cannabinoid 2-receptor (CB2) with potential as a therapy for cocaine and fentanyl overdose. All these results provide a better understanding of fentanyl and cocaine polydrug addiction and future drug abuse prevention.

Structural dynamics behind variants in pyrazinamidase and pyrazinamide resistance

Pyrazinamide (PZA) is an important component of first-line anti-tuberculosis (anti-TB) drugs. The anti-TB agent is activated into an active form, pyrazinoic acid (POA), by Mycobacterium tuberculosis (MTB) pncA gene encoding pyrazinamidase (PZase). The major cause of PZA-resistance has been associated with mutations in the pncA gene. We have detected several novel mutations including V131F, Q141P, R154T, A170P, and V180F (GeneBank Accession No. MH461111) in the pncA gene of PZA-resistant isolates during PZA drug susceptibility testing followed by pncA gene sequencing. Here, we investigated molecular mechanism of PZA-resistance by comparing the results of experimental and molecular dynamics. The mutants (MTs) and wild type (WT) PZase structures in apo and complex with PZA were subjected to molecular dynamic simulations (MD) at the 40 ns. Multiple factors, including root mean square deviations (RMSD), binding pocket, total energy, dynamic cross correlation, and root mean square fluctuations (RMSF) of MTs and WT were compared. The MTs attained a high deviation and fluctuation compared to WT. Binding pocket volumes of the MTs, were found, lower than the WT, and the docking scores were high than WT while shape complementarity scores were lower than that of the WT. Residual motion in MTs are seemed to be dominant in anti-correlated motion. Mutations at locations, V131F, Q141P, R154T, A170P, and V180F, might be involved in the structural changes, possibly affecting the catalytic property of PZase to convert PZA into POA. Our study provides useful information that will enhance the understanding for better management of TB. AbbreviationsDSTdrug susceptibility testingΔelecelectrostatic energyLJLowenstein-Jensen mediumMGITmycobacterium growth indicator tubesMTsmutantsMDmolecular dynamic simulationsMTBMycobacterium tuberculosisNALC-NaOHN-acetyl-l-cysteine-sodium hydroxideNIHNational Institutes of HealthNPTamount of substance (N), pressure (P) temperature (T)NVTmoles (N), volume (V) temperature (T)PZasepyrazinamidaseΔpspolar solvation energyPTRLProvincial Tuberculosis Reference LaboratoryRMSDroot mean square deviationsRMSFroot mean square fluctuationsΔSASAsolvent accessible surface area energyTBtuberculosisGTotaltotal binding free energyΔvdWVan der Waals energyWTwild typeCommunicated by Ramaswamy H. Sarma.

DAKB-GPCRs: An Integrated Computational Platform for Drug Abuse Related GPCRs

Drug abuse (DA) or drug addiction is a complicated brain disorder which is commonly considered as neurobiological impairments caused by both genetic factors and environmental effects. Among DA-related targets, G protein-coupled receptors (GPCRs) play an important role in DA therapy. However, only 52 GPCRs have been published with crystal structures in the recent two decades. In the effort to overcome the limitations of crystal structure and conformational diversity of GPCRs, we built homology models and performed conformational searches by molecular dynamics (MD) simulation. To accelerate and facilitate the drug abuse research, we construct a DA-related GPCR-specific chemogenomics knowledgebase (KB) (DAKB-GPCRs) for its research that can be implemented with our established and novel chemogenomics tools as well as algorithms for data analysis and visualization. Our established TargetHunter and HTDocking tools, as well as our novel tools that include target classification and Spider Plot, are compiled into the platform. Our DAKB-GPCRs provides the following results for a query compound: (1) blood-brain barrier (BBB) plot via our BBB predictor, (2) docking scores via HTDocking, (3) similarity score via TargetHunter, (4) target classification via machine learning methods that utilize both docking scores and similarity scores, and (5) a drug-target interaction network via Spider Plot.

Insight of Captagon Abuse by Chemogenomics Knowledgebase-guided Systems Pharmacology Target Mapping Analyses

Captagon, known by its genetic name Fenethylline, is an addictive drug that complicates the War on Drugs. Captagon has a strong CNS stimulating effect than its primary metabolite, Amphetamine. However, multi-targets issues associated with the drug and metabolites as well as its underlying mechanisms have not been fully defined. In the present work, we applied our established drug-abuse chemogenomics-knowledgebase systems pharmacology approach to conduct targets/off-targets mapping (SP-Targets) investigation of Captagon and its metabolites for hallucination addiction, and also analyzed the cell signaling pathways for both Amphetamine and Theophylline with data mining of available literature. Of note, Amphetamine, an agonist for trace amine-associated receptor 1 (TAAR1) with enhancing dopamine signaling (increase of irritability, aggression, etc.), is the main cause of Captagon addiction; Theophylline, an antagonist that blocks adenosine receptors (e.g. A2aR) in the brain responsible for restlessness and painlessness, may attenuate the behavioral sensitization caused by Amphetamine. We uncovered that Theophylline's metabolism and elimination could be retarded due to competition and/or blockage of the CYP2D6 enzyme by Amphetamine; We also found that the synergies between these two metabolites cause Captagon's psychoactive effects to act faster and far more potently than those of Amphetamine alone. We carried out further molecular docking modeling and molecular dynamics simulation to explore the molecular interactions between Amphetamine and Theophylline and their important GPCRs targets, including TAAR1 and adenosine receptors. All of the systems pharmacology analyses and results will shed light insight into a better understanding of Captagon addiction and future drug abuse prevention.

Heat activation mechanism of TRPV1: New insights from molecular dynamics simulation

As a member of the transient receptor potential (TRP) channels superfamily, the TRPV1 channel undergoes a closed-to-open gating transition in response to various physical and chemical stimuli including heat. Thanks to recent progress in cryo-electron microscopy, high-resolution structures are becoming available for various TRP channels including TRPV1. This has enabled us to study the molecular mechanism of TRPV1 channel gating by using molecular simulation. Here we review recent progress in molecular simulations of TRPV1 channel by us and others, with focus on our molecular dynamics (MD) simulations of TRPV1 at different temperatures. While no consensus has been reached on the heat activation mechanism of TRPV1, the simulations have offered specific predictions and models for future experimental studies to test.

Development and Testing of Druglike Screening Libraries

Although significant advances in experimental high throughput screening (HTS) have been made for drug lead identification, in silico virtual screening (VS) is indispensable owing to its unique advantage over experimental HTS, target-focused, cheap, and efficient, albeit its disadvantage of producing false positive hits. For both experimental HTS and VS, the quality of screening libraries is crucial and determines the outcome of those studies. In this paper, we first reviewed the recent progress on screening library construction. We realized the urgent need for compiling high-quality screening libraries in drug discovery. Then we compiled a set of screening libraries from about 20 million druglike ZINC molecules by running fingerprint-based similarity searches against known drug molecules. Lastly, the screening libraries were objectively evaluated using 5847 external actives covering more than 2000 drug targets. The result of the assessment is very encouraging. For example, with the Tanimoto coefficient being set to 0.75, 36% of external actives were retrieved and the enrichment factor was 13. Additionally, drug target family specific screening libraries were also constructed and evaluated. The druglike screening libraries are available for download from https://mulan.pharmacy.pitt.edu .

Inhibitory Reactivity of Capsaicin with α-Amylase and α-Glucosidase Related to Antidiabetes using Molecular Docking and Quantum Calculation Methods

This work aims to investigate the inhibitory activity of capsaicin, which is one of capsaicinoid compounds, on these enzymes using a molecular docking and quantum calculation. Acarbose, a commercial diabetes drug, was also investigated for comparison. The docking results revealed that acarbose yields better inhibition efficiency with binding free energy (ΔGbinding) of about -8.2 to -11.9 kcal/mol, and inhibition constant (Ki) of about 0.0002 to 0.4 µM, whereas capsaicin provided the ΔGbinding of -5.8 to -6.1 kcal/mol and Ki of 23.7 to 45.9 µM. The total binding energy (ΔEbinding) between each inhibitor and amino acids in active site of enzyme obtained from quantum calculation with MP2/6-31G(d,p) level is in agreement with the ΔGbinding, i.e. the ΔEbinding of acarbose was larger negative than that of capsaicin. The amino acids interacting with inhibitor as hydrogen bond mainly contribute to the total binding energy. Nevertheless, it could be concluded that capsaicinoids have high potential to be developed as an alternative drug for diabetes disease.

An insight into paracetamol and its metabolites using molecular docking and molecular dynamics simulation

Paracetamol is a relatively safe analgesia/antipyretic drug without the risks of addiction, dependence, tolerance, and withdrawal when used alone. However, when administrated in an opioid/paracetamol combination product, which often contains a large quantity of paracetamol, it can be potentially dangerous due to the risk of hepatotoxicity. Paracetamol is known to be metabolized into N-(4-hydroxyphenyl)-arachidonamide (AM404) via fatty acid amide hydrolase (FAAH) and into N-acetyl-p-benzoquinone imine (NAPQI) via cytochrome P450 (CYP) enzymes. However, the underlying mechanism of paracetamol is still unclear. In addition, paracetamol has the potential to interact with other drugs that are also involved with CYP family enzymes (inducer/inhibitor/substrate), an example being illicit drugs. In our present work, we looked into the relationship between paracetamol and its metabolites (AM404 and NAPQI) using molecular docking and molecular dynamics (MD) simulations. We first carried out a series of molecular docking studies between paracetamol/AM404/NAQPI and their reported targets, including CYP 2E1, FAAH, TRPA1, CB1, and TRPV1. Subsequently, we performed MD simulations and energy decomposition for CB1-AM404, TRPV1-AM404, and TRPV1-NAPQI for further investigation of the dynamics interactions. Finally, we summarized and discussed the reported drug-drug interactions between paracetamol and central nervous system drugs, especially illicit drugs. Overall, we are able to provide new insights into the structural and functional roles of paracetamol and its metabolites that can inform the potential prevention and treatment of paracetamol overdose. Graphical abstract Paracetamol and its metabolites.

The efficacy and safety of cilostazol as an alternative to aspirin in Chinese patients with aspirin intolerance after coronary stent implantation: a combined clinical study and computational system pharmacology analysis

Dual antiplatelet therapy (DAT) with aspirin and clopidogrel is the standard regimen to achieve rapid platelet inhibition and prevent thrombotic events. Currently, little information is available regarding alternative antiplatelet therapy in patients with an allergy or intolerance to aspirin. Although cilostazol is already a common alternative to aspirin in clinical practice in China, its efficacy and safety remain to be determined. We retrospectively analyzed 613 Chinese patients who had undergone primary percutaneous coronary intervention (PCI). Among them, 405 patients received standard DAT (aspirin plus clopidogrel) and 205 patients were identified with intolerance to aspirin and received alternative DAT (cilostazol plus clopidogrel). There were no significant differences between the two groups in their baseline clinical characteristics. The main outcomes of the study included major adverse cardiac events (MACEs) and bleeding events during 12 months of follow-up. The MACEs endpoint was reached in 10 of 205 patients treated with cilostazol (4.9%) and in 34 of 408 patients treated with aspirin (8.3%). No statistically significant difference was observed in MACEs between the two groups. However, patients in the cilostazol group had less restenosis than did patients in the aspirin group (1.5% vs 4.9%, P=0.035). The occurrence of bleeding events tended to be lower in the cilostazol group (0.49% vs 2.7%, P=0.063). These clinical observations were further analyzed using network system pharmacology analysis, and the outcomes were consistent with clinical observations and preclinical data reports. We conclude that in Chinese patients with aspirin intolerance undergoing coronary stent implantation, the combination of clopidogrel with cilostazol may be an efficacious and safe alternative to the standard DAT regimen.

Allosteric Modulation of Intact γ-Secretase Structural Dynamics

As a protease complex involved in the cleavage of amyloid precursor proteins that lead to the formation of amyloid β fibrils implicated in Alzheimer's disease, γ-secretase is an important target for developing therapeutics against Alzheimer's disease. γ-secretase is composed of four subunits: nicastrin (NCT) in the extracellular (EC) domain, presenilin-1 (PS1), anterior pharynx defective 1, and presenilin enhancer 2 in the transmembrane (TM) domain. NCT and PS1 play important roles in binding amyloid β precursor proteins and modulating PS1 catalytic activity. Yet, the molecular mechanisms of coupling between substrate/modulator binding and catalytic activity remain to be elucidated. Recent determination of intact human γ-secretase cryo-electron microscopy structure has opened the way for a detailed investigation of the structural dynamics of this complex. Our analysis, based on a membrane-coupled anisotropic network model, reveals two types of NCT motions, bending and twisting, with respect to PS1. These underlie the fluctuations between the "open" and "closed" states of the lid-like NCT with respect to a hydrophilic loop 1 (HL1) on PS1, thus allowing or blocking access of the substrate peptide (EC portion) to HL1 and to the neighboring helix TM2. In addition to this alternating access mechanism, fluctuations in the volume of the PS1 central cavity facilitate the exposure of the catalytic site for substrate cleavage. Druggability simulations show that γ-secretase presents several hot spots for either orthosteric or allosteric inhibition of catalytic activity, consistent with experimental data. In particular, a hinge region at the interface between the EC and TM domains, near the interlobe groove of NCT, emerges as an allo-targeting site that would impact the coupling between HL1/TM2 and the catalytic pocket, opening, to our knowledge, new avenues for structure-based design of novel allosteric modulators of γ-secretase protease activity.

Structural basis of TRPV5 channel inhibition by econazole revealed by cryo-EM

The transient receptor potential vanilloid 5 (TRPV5) channel is a member of the transient receptor potential (TRP) channel family, which is highly selective for Ca2+, that is present primarily at the apical membrane of distal tubule epithelial cells in the kidney and plays a key role in Ca2+ reabsorption. Here we present the structure of the full-length rabbit TRPV5 channel as determined using cryo-EM in complex with its inhibitor econazole. This structure reveals that econazole resides in a hydrophobic pocket analogous to that occupied by phosphatidylinositides and vanilloids in TRPV1, thus suggesting conserved mechanisms for ligand recognition and lipid binding among TRPV channels. The econazole-bound TRPV5 structure adopts a closed conformation with a distinct lower gate that occludes Ca2+ permeation through the channel. Structural comparisons between TRPV5 and other TRPV channels, complemented with molecular dynamics (MD) simulations of the econazole-bound TRPV5 structure, allowed us to gain mechanistic insight into TRPV5 channel inhibition by small molecules.

Novel Radiolabeled Vanilloid with Enhanced Specificity for Human Transient Receptor Potential Vanilloid 1 (TRPV1)

Transient receptor potential vanilloid 1 (TRPV1) has emerged as a promising therapeutic target. While radiolabeled resiniferatoxin (RTX) has provided a powerful tool for characterization of vanilloid binding to TRPV1, TRPV1 shows 20-fold weaker binding to the human TRPV1 than to the rodent TRPV1. We now describe a tritium radiolabeled synthetic vanilloid antagonist, 1-((2-(4-(methyl-[³H])piperidin-1-yl-4-[³H])-6-(trifluoromethyl)pyridin-3-yl)methyl)-3-(3-oxo-3,4-dihydro-2H-benzo[b][1,4]oxazin-8-yl)urea ([³H]MPOU), that embodies improved absolute affinity for human TRPV1 and improved synthetic accessibility.

Pyrazole C-region analogues of 2-(3-fluoro-4-methylsulfonylaminophenyl)propanamides as potent TRPV1 antagonists

A series of 1-substituted 3-(t-butyl/trifluoromethyl)pyrazole C-region analogues of 2-(3-fluoro-4-methylsulfonamidophenyl)propanamides were investigated for hTRPV1 antagonism. The structure activity relationship indicated that the 3-chlorophenyl group at the 1-position of pyrazole was the optimized hydrophobic group for antagonistic potency and the activity was stereospecific to the S-configuration, providing exceptionally potent antagonists 13S and 16S with K_i(CAP)=0.1nM. Particularly significant, 13S exhibited antagonism selective for capsaicin and NADA and not for low pH or elevated temperature. Both compounds also proved to be very potent antagonists for rTRPV1, blocking in vivo the hypothermic action of capsaicin, consistent with their in vitro mechanism. The docking study of compounds 13S and 16S in our hTRPV1 homology model indicated that the binding modes differed somewhat, with that of 13S more closely resembling that of GRT12360.

Discovery of Novel Ligands for TNF-α and TNF Receptor-1 through Structure-Based Virtual Screening and Biological Assay

Tumor necrosis factor α (TNF-α) is overexpressed in various diseases, and it has been a validated therapeutic target for autoimmune diseases. All therapeutics currently used to target TNF-α are biomacromolecules, and limited numbers of TNF-α chemical inhibitors have been reported, which makes the identification of small-molecule alternatives an urgent need. Recent studies have mainly focused on identifying small molecules that directly bind to TNF-α or TNF receptor-1 (TNFR1), inhibit the interaction between TNF-α and TNFR1, and/or regulate related signaling pathways. In this study, we combined in silico methods with biophysical and cell-based assays to identify novel antagonists that bind to TNF-α or TNFR1. Pharmacophore model filtering and molecular docking were applied to identify potential TNF-α antagonists. In regard to TNFR1, we constructed a three-dimensional model of the TNF-α-TNFR1 complex and carried out molecular dynamics simulations to sample the conformations. The residues in TNF-α that have been reported to play important roles in the TNF-α-TNFR1 complex were removed to form a pocket for further virtual screening of TNFR1-binding ligands. We obtained 20 virtual hits and tested them using surface plasmon resonance-based assays, which resulted in one ligand that binds to TNFR1 and four ligands with different scaffolds that bind to TNF-α. T1 and R1, the two most active compounds with Kd values of 11 and 16 μM for TNF-α and TNFR1, respectively, showed activities similar to those of known antagonists. Further cell-based assays also demonstrated that T1 and R1 have similar activities compared to the known TNF-α antagonist C87. Our work has not only produced several TNF-α and TNFR1 antagonists with novel scaffolds for further structural optimization but also showcases the power of our in silico methods for TNF-α- and TNFR1-based drug discovery.

t-Butyl pyridine and phenyl C-region analogues of 2-(3-fluoro-4-methylsulfonylaminophenyl)propanamides as potent TRPV1 antagonists

A series of 2-substituted 6-t-butylpyridine and 4-t-butylphenyl C-region analogues of 2-(3-fluoro-4-methylsulfonamidophenyl)propanamides were investigated for hTRPV1 antagonism. The analysis of structure activity relationships indicated that the pyridine derivatives generally exhibited a little better antagonism than did the corresponding phenyl surrogates for most of the series. Among the compounds, compound 7 showed excellent antagonism toward capsaicin activation with Ki=0.1nM and compound 60S demonstrated a strong antiallodynic effect with 83% MPE at 10mg/kg in the neuropathic pain model. The docking study of 7S in our hTRPV1 homology model indicated that the interactions between the A/B-regions of 7S with Tyr511 and the interactions between the t-butyl and ethyl groups in the C-region of 7S with the two hydrophobic binding pockets of hTRPV1 contributed to the high potency.

TRPV1 antagonism by piperazinyl-aryl compounds: A Topomer-CoMFA study and its use in virtual screening for identification of novel antagonists

Transient Receptor Potential Vanilloid, member 1 (TRPV1), is a non-selective cation channel belonging to the transient receptor potential (TRP) family of ion channels. It occurs in the peripheral and central nervous system, activated by a variety of exogenous and endogenous stimuli, thus playing a key role in transmission of pain. This has been a target for chronic pain since more than a decade and a number of antagonists that progressed into clinical trials have failed due to the unexpected side effect of core body temperature rise, thus halting progress in this field. Of late, there has been an upsurge in research on this target, with the rat TRPV1 structure being determined, many new antagonists discovered that are temperature-neutral and many new therapeutic avenues being discovered for TRPV1, including diseases of respiratory and digestive systems, skin and bladder. Towards identifying diverse compounds to decipher the role of this target in various indications, here we report a 3D-QSAR model built using the new topomer-CoMFA methodology on a series of piperazinyl-aryl TRPV1 antagonists and the use of this model, along with a pharmacophore model and the shape of one of the potent compounds of this series, to virtually screen a subset of the ZINC database to find novel and diverse hits. These can serve as starting points to develop modality-selective antagonists for chronic pain and to elucidate the critical role of TRPV1 in the various new therapeutic areas.

Cardiovascular Disease Chemogenomics Knowledgebase-guided Target Identification and Drug Synergy Mechanism Study of an Herbal Formula

Combination therapy is a popular treatment for various diseases in the clinic. Among the successful cases, Traditional Chinese Medicinal (TCM) formulae can achieve synergistic effects in therapeutics and antagonistic effects in toxicity. However, characterizing the underlying molecular synergisms for the combination of drugs remains a challenging task due to high experimental expenses and complication of multicomponent herbal medicines. To understand the rationale of combination therapy, we investigated Sini Decoction, a well-known TCM consisting of three herbs, as a model. We applied our established diseases-specific chemogenomics databases and our systems pharmacology approach TargetHunter to explore synergistic mechanisms of Sini Decoction in the treatment of cardiovascular diseases. (1) We constructed a cardiovascular diseases-specific chemogenomics database, including drugs, target proteins, chemicals, and associated pathways. (2) Using our implemented chemoinformatics tools, we mapped out the interaction networks between active ingredients of Sini Decoction and their targets. (3) We also in silico predicted and experimentally confirmed that the side effects can be alleviated by the combination of the components. Overall, our results demonstrated that our cardiovascular disease-specific database was successfully applied for systems pharmacology analysis of a complicated herbal formula in predicting molecular synergetic mechanisms, and led to better understanding of a combinational therapy.

Temperature-sensitive gating of TRPV1 channel as probed by atomistic simulations of its trans- and juxtamembrane domains

Heat-activated transient receptor potential channel TRPV1 is one of the most studied eukaryotic proteins involved in temperature sensation. Upon heating, it exhibits rapid reversible pore gating, which depolarizes neurons and generates action potentials. Underlying molecular details of such effects in the pore region of TRPV1 is of a crucial importance to control temperature responses of the organism. Despite the spatial structure of the channel in both open (O) and closed (C) states is known, microscopic nature of channel gating and mechanism of thermal sensitivity are still poorly understood. In this work, we used unrestrained atomistic molecular dynamics simulations of TRPV1 (without N- and C-terminal cytoplasmic domains) embedded into explicit lipid bilayer in its O- and C-states. We found that the pore domain with its neighboring loops undergoes large temperature-dependent conformational transitions in an asymmetric way, when fragments of only one monomer move with large amplitude, freeing the pore upon heating. Such an asymmetrical gating looks rather biologically relevant because it is faster and more reliable than traditionally proposed “iris-like” symmetric scheme of channel opening. Analysis of structural, dynamic, and hydrophobic organization of the pore domain revealed entropy growth upon TRPV1 gating, which is in line with current concepts of thermal sensitivity.

TRPV1: A Target for Rational Drug Design

Transient Receptor Potential Vanilloid 1 (TRPV1) is a non-selective, Ca²⁺ permeable cation channel activated by noxious heat, and chemical ligands, such as capsaicin and resiniferatoxin (RTX). Many compounds have been developed that either activate or inhibit TRPV1, but none of them are in routine clinical practice. This review will discuss the rationale for antagonists and agonists of TRPV1 for pain relief and other conditions, and strategies to develop new, better drugs to target this ion channel, using the newly available high-resolution structures.

Multi-Functional Diarylurea Small Molecule Inhibitors of TRPV1 with Therapeutic Potential for Neuroinflammation

Transient receptor potential vanilloid type 1 (TRPV1), a heat-sensitive calcium channel protein, contributes to inflammation as well as to acute and persistent pain. Since TRPV1 occupies a central position in pathways of neuronal inflammatory signaling, it represents a highly attractive potential therapeutic target for neuroinflammation. In the present work, we have in silico identified a series of diarylurea analogues for hTRPV1, of which 11 compounds showed activity in the nanomolar to micromolar range as validated by in vitro biological assays. Then, we utilized molecular docking to explore the detailed interactions between TRPV1 and the compounds to understand the contributions of the different substituent groups. Tyr511, Leu518, Leu547, Thr550, Asn551, Arg557, and Leu670 were important for the recognition of the small molecules by TRPV1. A hydrophobic group in R2 or a polar/hydrophilic group in R1 contributed significantly to the activities of the antagonists at TRPV1. In addition, the subtle different binding pose of meta-chloro in place of para-fluoro in the R2 group converted antagonism into partial agonism, as was predicted by our short-term molecular dynamics (MD) simulation and validated by bioassay. Importantly, compound 15, one of our best TRPV1 inhibitors, also showed potential binding affinity (1.39 μM) at cannabinoid receptor 2 (CB2), which is another attractive target for immuno-inflammation diseases. Furthermore, compound 1 and its diarylurea analogues were predicted to target the C-X-C chemokine receptor 2 (CXCR2), although bioassay validation of CXCR2 with these compounds still needs to be performed. This prediction from the modeling is of interest, since CXCR2 is also a potential therapeutic target for chronic inflammatory diseases. Our findings provide novel strategies to develop a small molecule inhibitor to simultaneously target two or more inflammation-related proteins for the treatment of a wide range of inflammatory disorders including neuroinflammation and neurodegenerative diseases with potential synergistic effect.

Difference and Influence of Inactive and Active States of Cannabinoid Receptor Subtype CB2: From Conformation to Drug Discovery

Cannabinoid receptor 2 (CB2), a G protein-coupled receptor (GPCR), is a promising target for the treatment of neuropathic pain, osteoporosis, immune system, cancer, and drug abuse. The lack of an experimental three-dimensional CB2 structure has hindered not only the development of studies of conformational differences between the inactive and active CB2 but also the rational discovery of novel functional compounds targeting CB2. In this work, we constructed models of both inactive and active CB2 by homology modeling. Then we conducted two comparative 100 ns molecular dynamics (MD) simulations on the two systems-the active CB2 bound with both the agonist and G protein and the inactive CB2 bound with inverse agonist-to analyze the conformational difference of CB2 proteins and the key residues involved in molecular recognition. Our results showed that the inactive CB2 and the inverse agonist remained stable during the MD simulation. However, during the MD simulations, we observed dynamical details about the breakdown of the "ionic lock" between R131(3.50) and D240(6.30) as well as the outward/inward movements of transmembrane domains of the active CB2 that bind with G proteins and agonist (TM5, TM6, and TM7). All of these results are congruent with the experimental data and recent reports. Moreover, our results indicate that W258(6.48) in TM6 and residues in TM4 (V164(4.56)-L169(4.61)) contribute greatly to the binding of the agonist on the basis of the binding energy decomposition, while residues S180-F183 in extracellular loop 2 (ECL2) may be of importance in recognition of the inverse agonist. Furthermore, pharmacophore modeling and virtual screening were carried out for the inactive and active CB2 models in parallel. Among all 10 hits, two compounds exhibited novel scaffolds and can be used as novel chemical probes for future studies of CB2. Importantly, our studies show that the hits obtained from the inactive CB2 model mainly act as inverse agonist(s) or neutral antagonist(s) at low concentration. Moreover, the hit from the active CB2 model also behaves as a neutral antagonist at low concentration. Our studies provide new insight leading to a better understanding of the structural and conformational differences between two states of CB2 and illuminate the effects of structure on virtual screening and drug design.

Structure-Driven Pharmacology of Transient Receptor Potential Channel Vanilloid 1

The transient receptor potential vanilloid 1 (TRPV1) ion channel is a polymodal receptor that mediates the flux of cations across the membrane in response to several stimuli, including heat, voltage, and ligands. The best known agonist of TRPV1 channels is capsaicin, the pungent component of "hot" chili peppers. In addition, peptides found in the venom of poisonous animals, along with the lipids phosphatidylinositol 4,5-biphosphate, lysophosphatidic acid, and cholesterol, bind to TRPV1 with high affinity to modulate channel gating. Here, we discuss the functional evidence regarding ligand-dependent activation of TRPV1 channels in light of structural data recently obtained by cryoelectron microscopy. This review focuses on the mechanistic insights into ligand binding and allosteric gating of TRPV1 channels and the relevance of accurate polymodal receptor biophysical characterization for drug design in novel pain therapies.

Metal binding mediated conformational change of XPA protein:a potential cytotoxic mechanism of nickel in the nucleotide excision repair

Nucleotide excision repair (NER) is a pivotal life process for repairing DNA nucleotide mismatch caused by chemicals, metal ions, radiation, and other factors. As the initiation step of NER, the xeroderma pigmentosum complementation group A protein (XPA) recognizes damaged DNA molecules, and recruits the replication protein A (RPA), another important player in the NER process. The stability of the Zn(2+)-chelated Zn-finger domain of XPA center core portion (i.e., XPA98-210) is the foundation of its biological functionality, while the displacement of the Zn(2+) by toxic metal ions (such as Ni(2+), a known human carcinogen and allergen) may impair the effectiveness of NER and hence elevate the chance of carcinogenesis. In this study, we first calculated the force field parameters for the bonded model in the metal center of the XPA98-210 system, showing that the calculated results, including charges, bonds, angles etc., are congruent with previously reported results measured by spectrometry experiments and quantum chemistry computation. Then, comparative molecular dynamics simulations using these parameters revealed the changes in the conformation and motion mode of XPA98-210 Zn-finger after the substitution of Zn(2+) by Ni(2+). The results showed that Ni(2+) dramatically disrupted the relative positions of the four Cys residues in the Zn-finger structure, forcing them to collapse from a tetrahedron into an almost planar structure. Finally, we acquired the binding mode of XPA98-210 with its ligands RPA70N and DNA based on molecular docking and structural alignment. We found that XPA98-210's Zn-finger domain primarily binds to a V-shaped cleft in RPA70N, while the cationic band in its C-terminal subdomain participates in the recognition of damaged DNA. In addition, this article sheds light on the multi-component interaction pattern among XPA, DNA, and other NER-related proteins (i.e., RPA70N, RPA70A, RPA70B, RPA70C, RPA32, and RPA14) based on previously reported structural biology information. Thus, we derived a putative cytotoxic mechanism associated with the nickel ion, where the Ni(2+) disrupts the conformation of the XPA Zn-finger, directly weakening its interaction with RPA70N, and thus lowering the effectiveness of the NER process. In sum, this work not only provides a theoretical insight into the multi-protein interactions involved in the NER process and potential cytotoxic mechanism associated with Ni(2+) binding in XPA, but may also facilitate rational anti-cancer drug design based on the NER mechanism.

Design and Synthesis of New Transient Receptor Potential Vanilloid Type-1 (TRPV1) Channel Modulators: Identification, Molecular Modeling Analysis, and Pharmacological Characterization of the N-(4-Hydroxy-3-methoxybenzyl)-4-(thiophen-2-yl)butanamide, a Small Molecule Endowed with Agonist TRPV1 Activity and Protective Effects against Oxidative Stress

4-(Thiophen-2-yl)butanoic acid was identified as a cyclic substitute of the unsaturated alkyl chain of the natural ligand, capsaicin. Accordingly, a new class of amides was synthesized in good yield and high purity and their molecular recognition against the target was investigated by means of docking experiments followed by molecular dynamics simulations, in order to rationalize their geometrical and thermodynamic profiles. The pharmacological properties of these new compounds were expressed as activation (EC50) and desensitization (IC50) potencies. Several compounds were found to activate TRPV1 channels, and in particular, derivatives 1 and 10 behaved as TRPV1 agonists endowed with good efficacy as compared to capsaicin. The most promising compound 1 was also evaluated for its protective role against oxidative stress on keratinocytes and differentiated human neuroblastoma cell lines expressing the TRPV1 receptor as well as for its cytotoxicity and analgesic activity in vivo.

Discovery of N-(3-fluoro-4-methylsulfonamidomethylphenyl)urea as a potent TRPV1 antagonistic template

A series of homologous analogues of prototype antagonist 1 and its urea surrogate were investigated as hTRPV1 ligands. Through one-carbon elongation in the respective pharmacophoric regions, N-(3-fluoro-4-methylsulfonamidomethylphenyl)urea was identified as a novel and potent TRPV1 antagonistic template. Its representative compound 27 showed a potency comparable to that of lead compound 1. Docking analysis of compound 27 in our hTRPV1 homology model indicated that its binding mode was similar with that of 1S.

Engineering vanilloid-sensitivity into the rat TRPV2 channel

The TRPV1 channel is a detector of noxious stimuli, including heat, acidosis, vanilloid compounds and lipids. The gating mechanisms of the related TRPV2 channel are poorly understood because selective high affinity ligands are not available, and the threshold for heat activation is extremely high (>50°C). Cryo-EM structures of TRPV1 and TRPV2 reveal that they adopt similar structures, and identify a putative vanilloid binding pocket near the internal side of TRPV1. Here we use biochemical and electrophysiological approaches to investigate the resiniferatoxin(RTx) binding site in TRPV1 and to explore the functional relationships between TRPV1 and TRPV2. Collectively, our results support the interaction of vanilloids with the proposed RTx binding pocket, and demonstrate an allosteric influence of a tarantula toxin on vanilloid binding. Moreover, we show that sensitivity to RTx can be engineered into TRPV2, demonstrating that the gating and permeation properties of this channel are similar to TRPV1.

2-Sulfonamidopyridine C-region analogs of 2-(3-fluoro-4-methylsulfonamidophenyl)propanamides as potent TRPV1 antagonists

A series of 2-sulfonamidopyridine C-region derivatives of 2-(3-fluoro-4-methylsulfonamidophenyl)propanamide were investigated as hTRPV1 ligands. Systematic modification on the 2-sulfonamido group provided highly potent TRPV1 antagonists. The N-benzyl phenylsulfonamide derivatives 12 and 23 in particular showed higher affinities than that of lead compound 1. Compound 12 exhibited strong analgesic activity in the formalin pain model. Docking analysis of its chiral S-form 12S in our hTRPV1 homology model indicated that its high affinity might arise from additional hydrophobic interactions not present in lead compound 1S.

Cloning, Expression, Mutagenesis Library Construction of Glycerol Dehydratase, and Binding Mode Simulation of Its Reactivase with Ligands

The production of 1, 3-propanediol (1, 3-PD) and 3-hydroxypropionaldehyde (3-HPA) by enzyme reaction has been a hot field, and glycerol dehydratase (GDHt) is the key and rate-limiting enzyme involved in their biosynthesis. The gldABC gene encoding GDHt was cloned from Klebsiella pneumoniae, and the activity of the corresponding proteins expressed extracellularly and intracellularly was 6.8 and 3.2 U/mg, respectively, about six and three times higher than that of the wild strain. The change of amino acids for the β subunit can adjust the length of the Co-N bond and affect the homolysis rate of the Co-C bond to change GDHt activity. The expression plasmid, pET-32a-gldAC (containing no gldB which encodes the β subunit of GDHt), was constructed to build the mutagenesis library to improve the GDHt activity. The binding models of glycerol dehydratase reactivation factor (GDHtR) with ATP, CTP, or GTP were simulated by semi-flexible docking, respectively, and there was almost no difference between them. This research provided the basis for studying the quantitative structure-activity relationships between GDHtR and its ligands, as well as searching inexpensive ligands to replace ATP. These results and methods are of great use in economical and highly efficient production of 3-HPA and 1, 3-PD by the enzyme method.