In silico assessment of biocompatibility and toxicity: molecular docking and dynamics simulation of PMMA-based dental materials for interim prosthetic restorations.

AIM: This study aimed to comprehensively assess the biocompatibility and toxicity profiles of poly(methyl methacrylate) (PMMA) and its monomeric unit, methyl methacrylate (MMA), crucial components in dental materials for interim prosthetic restorations. METHODOLOGY: Molecular docking was employed to predict the binding affinities, energetics, and steric features of MMA and PMMA with selected receptors involved in bone metabolism and tissue development, including RANKL, Fibronectin, BMP9, NOTCH2, and other related receptors. The HADDOCK standalone version was utilized for docking calculations, employing a Lamarckian genetic algorithm to explore the conformational space of ligand-receptor interactions. Furthermore, molecular dynamics (MD) simulations over 100 nanoseconds were conducted using the GROMACS package to evaluate dynamic actions and structural stability. The LigandScout was utilized for pharmacophore modeling, which employs a shape-based screening approach to identify potential ligand binding sites on protein targets. RESULTS: The molecular docking studies elucidated promising interactions between PMMA and MMA with key biomolecular targets relevant to dental applications. MD simulation results provided strong evidence supporting the structural stability of PMMA complexes over time. Pharmacophore modeling highlighted the significance of carbonyl and hydroxyl groups as pharmacophoric features, indicating compounds with favorable biocompatibility profiles. CONCLUSION: This study underscores the potential of PMMA in dental applications, emphasizing its structural stability, molecular interactions, and safety considerations. These findings lay a foundation for future advancements in dental biomaterials, guiding the design and optimization of materials for enhanced biocompatibility. Future directions include experimental validation of computational findings and the development of PMMA-based dental materials with improved biocompatibility and clinical performance.

Watermelon: setup and validation of an in silico fragment-based approach.

We present a new computational approach, named Watermelon, designed for the development of pharmacophore models based on receptor structures. The methodology involves the sampling of potential hotspots for ligand interactions within a protein target's binding site, utilising molecular fragments as probes. By employing docking and molecular dynamics (MD) simulations, the most significant interactions formed by these probes within distinct regions of the binding site are identified. These interactions are subsequently transformed into pharmacophore features that delineates key anchoring sites for potential ligands. The reliability of the approach was experimentally validated using the monoacylglycerol lipase (MAGL) enzyme. The generated pharmacophore model captured features representing ligand-MAGL interactions observed in various X-ray co-crystal structures and was employed to screen a database of commercially available compounds, in combination with consensus docking and MD simulations. The screening successfully identified two new MAGL inhibitors with micromolar potency, thus confirming the reliability of the Watermelon approach.

Water Ligands Regulate the Redox Leveling Mechanism of the Oxygen-Evolving Complex of the Photosystem II.

Understanding how water ligands regulate the conformational changes and functionality of the oxygen-evolving complex (OEC) in photosystem II (PSII) throughout the catalytic cycle of oxygen evolution remains a highly intriguing and unresolved challenge. In this study, we investigate the effect of water insertion (WI) on the redox state of the OEC by using the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) hybrid methods. We find that water binding significantly reduces the free energy change for proton-coupled electron transfer (PCET) from Mn to YZ•, underscoring the important regulatory role of water binding, which is essential for enabling the OEC redox-leveling mechanism along the catalytic cycle. We propose a water binding mechanism in which WI is thermodynamically favored by the closed-cubane form of the OEC, with water delivery mediated by Ca2+ ligand exchange. Isomerization from the closed- to open-cubane conformation at three post-WI states highlights the importance of the location of the MnIII center in the OEC and the orientation of its Jahn-Teller axis to conformational changes of the OEC, which might be critical for the formation of the O-O bond. These findings reveal a complex interplay between conformational changes in the OEC and the ligand environment during the activation of the OEC by YZ•. Analogous regulatory effects due to water ligand binding are expected to be important for a wide range of catalysts activated by redox state transitions in aqueous environments.

Atomistic description of the OCTN1 recognition mechanism via in silico methods.

The Organic Cation Transporter Novel 1 (OCTN1), also known as SLC22A4, is widely expressed in various human tissues, and involved in numerous physiological and pathological processes remains. It facilitates the transport of organic cations, zwitterions, with selectivity for positively charged solutes. Ergothioneine, an antioxidant compound, and acetylcholine (Ach) are among its substrates. Given the lack of experimentally solved structures of this protein, this study aimed at generating a reliable 3D model of OCTN1 to shed light on its substrate-binding preferences and the role of sodium in substrate recognition and transport. A chimeric model was built by grafting the large extracellular loop 1 (EL1) from an AlphaFold-generated model onto a homology model. Molecular dynamics simulations revealed domain-specific mobility, with EL1 exhibiting the highest impact on overall stability. Molecular docking simulations identified cytarabine and verapamil as highest affinity ligands, consistent with their known inhibitory effects on OCTN1. Furthermore, MM/GBSA analysis allowed the categorization of substrates into weak, good, and strong binders, with molecular weight strongly correlating with binding affinity to the recognition site. Key recognition residues, including Tyr211, Glu381, and Arg469, were identified through interaction analysis. Ach demonstrated a low interaction energy, supporting the hypothesis of its one-directional transport towards to outside of the membrane. Regarding the role of sodium, our model suggested the involvement of Glu381 in sodium binding. Molecular dynamics simulations of systems at increasing levels of Na+ concentrations revealed increased sodium occupancy around Glu381, supporting experimental data associating Na+ concentration to molecule transport. In conclusion, this study provides valuable insights into the 3D structure of OCTN1, its substrate-binding preferences, and the role of sodium in the recognition. These findings contribute to the understanding of OCTN1 involvement in various physiological and pathological processes and may have implications for drug development and disease management.

Thermostability-assisted limited proteolysis-coupled mass spectrometry for capturing drug target proteins and sites.

BACKGROUND: Identifying drug-binding targets and their corresponding sites is crucial for drug discovery and mechanism studies. Limited proteolysis-coupled mass spectrometry (LiP-MS) is a sophisticated method used for the detection of compound and protein interactions. However, in some cases, LiP-MS cannot identify the target proteins due to the small structure changes or the lack of enrichment of low-abundant protein. To overcome this drawback, we developed a thermostability-assisted limited proteolysis-coupled mass spectrometry (TALiP-MS) approach for efficient drug target discovery. RESULTS: We proved that the novel strategy, TALiP-MS, could efficiently identify target proteins of various ligands, including cyclosporin A (a calcineurin inhibitor), geldanamycin (an HSP90 inhibitor), and staurosporine (a kinase inhibitor), with accurately recognizing drug-binding domains. The TALiP protocol increased the number of target peptides detected in LiP-MS experiments by 2- to 8-fold. Meanwhile, the TALiP-MS approach can not only identify both ligand-binding stability and destabilization proteins but also shows high complementarity with the thermal proteome profiling (TPP) and machine learning-based limited proteolysis (LiP-Quant) methods. The developed TALiP-MS approach was applied to identify the target proteins of celastrol (CEL), a natural product known for its strong antioxidant and anti-cancer angiogenesis effect. Among them, four proteins, MTHFD1, UBA1, ACLY, and SND1 were further validated for their strong affinity to CEL by using cellular thermal shift assay. Additionally, the destabilized proteins induced by CEL such as TAGLN2 and CFL1 were also validated. SIGNIFICANCE: Collectively, these findings underscore the efficacy of the TALiP-MS method for identifying drug targets, elucidating binding sites, and even detecting drug-induced conformational changes in target proteins in complex proteomes.

Rational Design of Drugs Targeting G-Protein-Coupled Receptors: A Structural Biology Perspective.

G protein-coupled receptors (GPCRs) play a key role in the transduction of extracellular signals to cells and regulation of many biological processes, which makes these membrane proteins one of the most important targets for pharmacological agents. A significant increase in the number of resolved atomic structures of GPCRs has opened the possibility of developing pharmaceuticals targeting these receptors via structure-based drug design (SBDD). SBDD employs information on the structure of receptor-ligand complexes to search for selective ligands without the need for an extensive high-throughput experimental ligand screening and can significantly expand the chemical space for ligand search. In this review, we describe the process of deciphering GPCR structures using X-ray diffraction analysis and cryoelectron microscopy as an important stage in the rational design of drugs targeting this receptor class. Our main goal was to present modern developments and key features of experimental methods used in SBDD of GPCR-targeting agents to a wide range of specialists.

Highly biased agonism for GPCR ligands via nanobody tethering.

Ligand-induced activation of G protein-coupled receptors (GPCRs) can initiate signaling through multiple distinct pathways with differing biological and physiological outcomes. There is intense interest in understanding how variation in GPCR ligand structure can be used to promote pathway selective signaling ("biased agonism") with the goal of promoting desirable responses and avoiding deleterious side effects. Here we present an approach in which a conventional peptide ligand for the type 1 parathyroid hormone receptor (PTHR1) is converted from an agonist which induces signaling through all relevant pathways to a compound that is highly selective for a single pathway. This is achieved not through variation in the core structure of the agonist, but rather by linking it to a nanobody tethering agent that binds with high affinity to a separate site on the receptor not involved in signal transduction. The resulting conjugate represents the most biased agonist of PTHR1 reported to date. This approach holds promise for facile generation of pathway selective ligands for other GPCRs.

Use of Surface Plasmon Resonance Technique for Studies of Inter-domain Interactions in Ion Channels.

Ion channels are transmembrane proteins essential for cellular functions and are important drug targets. Surface plasmon resonance (SPR) is a powerful technique for investigating protein-protein and protein-small molecule ligand interactions. SPR has been underutilized for studies of ion channels, even though it could provide a wealth of information on the mechanisms of ion channel regulation and aid in ion channel drug discovery. Here we provide a detailed description of the use of SPR technology for investigating inter-domain interactions in KCNH potassium-selective and voltage-gated ion channels.

Isothermal Titration Calorimetry for Fragment-Based Analysis of Ion Channel Interactions.

Ion channels are membrane proteins that may also have intracellular and extracellular domains that interact with other ligands. In many cases, these interaction sites are highly mobile and may undergo changes in the configuration upon binding with regulatory signaling molecules. Isothermal titration calorimetry (ITC) is a powerful technique to quantify protein-ligand interactions of purified samples in solution. This chapter describes a fragment-based analysis method using ITC to quantify the interactions between a domain of the voltage-gated Kv7 channel and the calcium-regulated protein calmodulin. This example can be used to quantify the interactions between specific domains of other ion channels and their regulatory signaling proteins.

[Microorganism Immobilization Device Using Artificial Siderophores].

Inspired by the mechanism by which microorganisms utilize siderophores to ingest iron, four different FeIII complexes of typical artificial siderophore ligands containing catecholate and/or hydroxamate groups, K3[FeIII-LC3], K2[FeIII-LC2H1], K[FeIII-LC1H2], and [FeIII-LH3], were prepared. They were modified on an Au substrate surface (Fe-L/Au) and applied as microorganism immobilization devices for fast, sensitive, selective detection of microorganisms, where H6LC3, H5LC2H1, H4LC1H2, and H3LH3 denote the tri-catecholate, biscatecholate-monohydroxamate, monocatecholate-bishydroxamate, and tri-hydroxamate type of artificial siderophores, respectively. Their adsorption properties for the several microorganisms were investigated using scanning electron microscopy (SEM), quartz crystal microbalance (QCM), and electric impedance spectroscopy (EIS) methods. The artificial siderophore-iron complexes modified on the Au substrates Fe-LC3/Au, Fe-LC2H1/Au, Fe-LC1H2/Au, and Fe-LH3/Au showed specific microorganism immobilization behavior with selectivity based on the structure of the artificial siderophores. Their specificities corresponded well with the structural characteristics of natural siderophores that microorganisms release from the cell and/or use to take up an iron. These findings suggest that release and uptake are achieved through specific interactions between the artificial siderophore-FeIII complexes and receptors on the cell surfaces of microorganisms. This study revealed that Fe-L/Au systems have specific potential to serve as effective immobilization probes of microorganisms for rapid, selective detection and identification of a variety of microorganisms.

Identification of structural fingerprints among natural inhibitors of HDAC1 to accelerate nature-inspired drug discovery in cancer epigenetics.

Histone deacetylase 1 (HDAC1), a class I HDAC enzyme, is crucial for histone modification. Currently, it is emerged as one of the important biological targets for designing small molecule drugs through cancer epigenetics. Along with synthetic inhibitors different natural inhibitors are showing potential HDAC1 inhibitions. In order to gain insights into the relationship between the molecular structures of the natural inhibitors and HDAC1, different molecular modelling techniques (Bayesian classification, recursive partitioning, molecular docking and molecular dynamics simulations) have been applied on a dataset of 155 HDAC1 nature-inspired inhibitors with diverse scaffolds. The Bayesian study showed acceptable ROC values for both the training set and test sets. The Recursive partitioning study produced decision tree 1 with 6 leaves. Further, molecular docking study was processed for generating the protein ligand complex which identified some potential amino acid residues such as F205, H28, L271, P29, F150, Y204 for the binding interactions in case of natural inhibitors. Stability of these HDAC1-natutal inhibitors complexes has been also evaluated by molecular dynamics simulation study. The current modelling study is an attempt to get a deep insight into the different important structural fingerprints among different natural compounds modulating HDAC1 inhibition.Communicated by Ramaswamy H. Sarma.

Revolutionizing GPCR-ligand predictions: DeepGPCR with experimental validation for high-precision drug discovery.

G-protein coupled receptors (GPCRs), crucial in various diseases, are targeted of over 40% of approved drugs. However, the reliable acquisition of experimental GPCRs structures is hindered by their lipid-embedded conformations. Traditional protein-ligand interaction models falter in GPCR-drug interactions, caused by limited and low-quality structures. Generalized models, trained on soluble protein-ligand pairs, are also inadequate. To address these issues, we developed two models, DeepGPCR_BC for binary classification and DeepGPCR_RG for affinity prediction. These models use non-structural GPCR-ligand interaction data, leveraging graph convolutional networks and mol2vec techniques to represent binding pockets and ligands as graphs. This approach significantly speeds up predictions while preserving critical physical-chemical and spatial information. In independent tests, DeepGPCR_BC surpassed Autodock Vina and Schrödinger Dock with an area under the curve of 0.72, accuracy of 0.68 and true positive rate of 0.73, whereas DeepGPCR_RG demonstrated a Pearson correlation of 0.39 and root mean squared error of 1.34. We applied these models to screen drug candidates for GPR35 (Q9HC97), yielding promising results with three (F545-1970, K297-0698, S948-0241) out of eight candidates. Furthermore, we also successfully obtained six active inhibitors for GLP-1R. Our GPCR-specific models pave the way for efficient and accurate large-scale virtual screening, potentially revolutionizing drug discovery in the GPCR field.

Kinetic and dynamical properties of truncated hemoglobins of the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125.

Due to the low temperature, the Antarctic marine environment is challenging for protein functioning. Cold-adapted organisms have evolved proteins endowed with higher flexibility and lower stability in comparison to their thermophilic homologs, resulting in enhanced reaction rates at low temperatures. The Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 (PhTAC125) genome is one of the few examples of coexistence of multiple hemoglobin genes encoding, among others, two constitutively transcribed 2/2 hemoglobins (2/2Hbs), also named truncated Hbs (TrHbs), belonging to the Group II (or O), annotated as PSHAa0030 and PSHAa2217. In this work, we describe the ligand binding kinetics and their interrelationship with the dynamical properties of globin Ph-2/2HbO-2217 by combining experimental and computational approaches and implementing a new computational method to retrieve information from molecular dynamic trajectories. We show that our approach allows us to identify docking sites within the protein matrix that are potentially able to transiently accommodate ligands and migration pathways connecting them. Consistently with ligand rebinding studies, our modeling suggests that the distal heme pocket is connected to the solvent through a low energy barrier, while inner cavities play only a minor role in modulating rebinding kinetics.

Biased signalling in analgesic research and development.

Ligand bias offers a novel means to improve the therapeutic profile of drugs. With regard to G protein-coupled receptors involved in analgesia, it could be advantageous to develop such drugs if the analgesic effect is mediated by a different cellular signalling pathway than the adverse effects associated with the drug. Whilst this has been explored over a number of years for the µ receptor, it remains unclear whether this approach offers significant benefit for the treatment of pain. Nevertheless, the development of biased ligands at other G protein-coupled receptors in the CNS does offer some promise for the development of novel analgesic drugs in the future. Here we summarise and discuss the recent evidence to support this.

Characterizing Modulators of Protease-activated Receptors with a Calcium Mobilization Assay Using a Plate Reader.

Changes in calcium concentration in cells are rapidly monitored in a high-throughput fashion with the use of intracellular, fluorescent, calcium-binding dyes and imaging instruments that can measure fluorescent emissions from up to 1,536 wells simultaneously. However, these instruments are much more expensive and can be challenging to maintain relative to widely available plate readers that scan wells individually. Described here is an optimized plate reader assay for use with an endothelial cell line (EA.hy926) to measure the protease-activated receptor (PAR)-driven activation of Gαq signaling and subsequent calcium mobilization using the calcium-binding dye Fluo-4. This assay has been used to characterize a range of PAR ligands, including the allosteric PAR1-targeting anti-inflammatory "parmodulin" ligands identified in the Dockendorff lab. This protocol obviates the need for an automated liquid handler and permits the medium-throughput screening of PAR ligands in 96-well plates and should be applicable to the study of other receptors that initiate calcium mobilization.

Q-BioLiP: A Comprehensive Resource for Quaternary Structure-based Protein-ligand Interactions.

Since its establishment in 2013, BioLiP has become one of the widely used resources for protein-ligand interactions. Nevertheless, several known issues occurred with it over the past decade. For example, the protein-ligand interactions are represented in the form of single chain-based tertiary structures, which may be inappropriate as many interactions involve multiple protein chains (known as quaternary structures). We sought to address these issues, resulting in Q-BioLiP, a comprehensive resource for quaternary structure-based protein-ligand interactions. The major features of Q-BioLiP include: (1) representing protein structures in the form of quaternary structures rather than single chain-based tertiary structures; (2) pairing DNA/RNA chains properly rather than separation; (3) providing both experimental and predicted binding affinities; (4) retaining both biologically relevant and irrelevant interactions to alleviate the wrong justification of ligands' biological relevance; and (5) developing a new quaternary structure-based algorithm for the modelling of protein-ligand complex structure. With these new features, Q-BioLiP is expected to be a valuable resource for studying biomolecule interactions, including protein-small molecule interaction, protein-metal ion interaction, protein-peptide interaction, protein-protein interaction, protein-DNA/RNA interaction, and RNA-small molecule interaction. Q-BioLiP is freely available at https://yanglab.qd.sdu.edu.cn/Q-BioLiP/.

Trace amine associated receptor 1: predicted effects of single nucleotide variants on structure-function in geographically diverse populations.

Trace Amine Associated Receptor 1 (TAAR1) is a novel pharmaceutical target under investigation for the treatment of several neuropsychiatric conditions. TAAR1 single nucleotide variants (SNV) have been found in patients with schizophrenia and metabolic disorders. However, the frequency of variants in geographically diverse populations and the functional effects of such variants are unknown. In this study, we aimed to characterise the distribution of TAAR1 SNVs in five different WHO regions using the Database of Genotypes and Phenotypes (dbGaP) and conducted a critical computational analysis using available TAAR1 structural data to identify SNVs affecting ligand binding and/or functional regions. Our analysis shows 19 orthosteric, 9 signalling and 16 micro-switch SNVs hypothesised to critically influence the agonist induced TAAR1 activation. These SNVs may non-proportionally influence populations from discrete regions and differentially influence the activity of TAAR1-targeting therapeutics in genetically and geographically diverse populations. Notably, our dataset presented with orthosteric SNVs D1033.32N (found only in the South-East Asian Region and Western Pacific Region) and T1945.42A (found only in South-East Asian Region), and 2 signalling SNVs (V1253.54A/T2526.36A, found in African Region and commonly, respectively), all of which have previously demonstrated to influence ligand induced functions of TAAR1. Furthermore, bioinformatics analysis using SIFT4G, MutationTaster 2, PROVEAN and MutationAssessor predicted all 16 micro-switch SNVs are damaging and may further influence the agonist activation of TAAR1, thereby possibly impacting upon clinical outcomes. Understanding the genetic basis of TAAR1 function and the impact of common mutations within clinical populations is important for the safe and effective utilisation of novel and existing pharmacotherapies.

Ligand-based design of [18F]OXD-2314 for PET imaging in non-Alzheimer's disease tauopathies.

Positron emission tomography (PET) imaging of tau aggregation in Alzheimer's disease (AD) is helping to map and quantify the in vivo progression of AD pathology. To date, no high-affinity tau-PET radiopharmaceutical has been optimized for imaging non-AD tauopathies. Here we show the properties of analogues of a first-in-class 4R-tau lead, [18F]OXD-2115, using ligand-based design. Over 150 analogues of OXD-2115 were synthesized and screened in post-mortem brain tissue for tau affinity against [3H]OXD-2115, and in silico models were used to predict brain uptake. [18F]OXD-2314 was identified as a selective, high-affinity non-AD tau PET radiotracer with favorable brain uptake, dosimetry, and radiometabolite profiles in rats and non-human primate and is being translated for first-in-human PET studies.

Rational Design of Drugs Targeting G-Protein-Coupled Receptors: Ligand Search and Screening.

G protein-coupled receptors (GPCRs) are transmembrane proteins that participate in many physiological processes and represent major pharmacological targets. Recent advances in structural biology of GPCRs have enabled the development of drugs based on the receptor structure (structure-based drug design, SBDD). SBDD utilizes information about the receptor-ligand complex to search for suitable compounds, thus expanding the chemical space of possible receptor ligands without the need for experimental screening. The review describes the use of structure-based virtual screening (SBVS) for GPCR ligands and approaches for the functional testing of potential drug compounds, as well as discusses recent advances and successful examples in the application of SBDD for the identification of GPCR ligands.

Construction and application of H1R ligand screening materials based on SMA stabilization and his-tag covalent immobilization of membrane proteins.

The histamine H1 receptor (H1R) plays a pivotal role in allergy initiation and undergoes the necessity of devising a high-throughput screening approach centered on H1R to screen novel ligands effectively. This study suggests a method employing styrene maleic acid (SMA) extraction and His-tag covalent bonding to immobilize H1R membrane proteins, minimizing the interference of nonspecific proteins interference while preserving native protein structure and maximizing target exposure. This approach was utilized to develop a novel material for high-throughput ligand screening and implemented in cell membrane chromatography (CMC). An H1R-His-SMALPs/CMC model was established and its chromatographic performance (selectivity, specificity and lifespan) validated, demonstrating a significant enhancement in lifespan compared to previous CMC models. Subsequently, this model facilitated high-throughput screening of H1R ligands in the compound library and preliminary activity verification of potential H1R antagonists. Identification of a novel H1R antagonist laid the foundation for further development in this area.