Comparison of the Glucocorticoid Receptor Binding and Agonist Activities of Typical Glucocorticoids: Insights into Their Endocrine Disrupting Effects

Over the past decades, the synthetic glucocorticoids (GCs) have been widely used in clinical practice and animal husbandry. Given the health hazard of these toxic residues in food, it is necessary to explore the detailed interaction mechanisms of typical GCs and their main target glucocorticoid receptor (GR). Hence, this work compared the GR binding and agonist activities of typical GCs. Fluorescence polarization assay showed that these GCs were potent ligands of GR. Their GR binding affinities were in the order of methylprednisolone>betamethasone≈prednisolone>dexamethasone, with IC50 values of 1.67, 2.94, 2.95, and 5.58 nM. Additionally, the limits of detection of dexamethasone, betamethasone, prednisolone, and methylprednisolone were 0.32, 0.14, 0.19, and 0.09 μg/kg in fluorescence polarization assay. Reporter gene assay showed that these GCs induced GR transactivation in a dose-dependent manner, confirming their GR agonist activities. Among which, dexamethasone at the concentration of 100 nM produced a maximal induction of more than 11-fold over the blank control. Molecular docking and molecular dynamics simulations suggested that hydrogen-bonding and hydrophobic interactions played an important role in stabilizing the GC-GR-LBD complexes. In summary, this work might help to understand the GR-mediated endocrine disrupting effects of typical GCs.

The Impact of Data on Structure-Based Binding Affinity Predictions Using Deep Neural Networks

Artificial intelligence (AI) has gained significant traction in the field of drug discovery, with deep learning (DL) algorithms playing a crucial role in predicting protein-ligand binding affinities. Despite advancements in neural network architectures, system representation, and training techniques, the performance of DL affinity prediction has reached a plateau, prompting the question of whether it is truly solved or if the current performance is overly optimistic and reliant on biased, easily predictable data. Like other DL-related problems, this issue seems to stem from the training and test sets used when building the models. In this work, we investigate the impact of several parameters related to the input data on the performance of neural network affinity prediction models. Notably, we identify the size of the binding pocket as a critical factor influencing the performance of our statistical models; furthermore, it is more important to train a model with as much data as possible than to restrict the training to only high-quality datasets. Finally, we also confirm the bias in the typically used current test sets. Therefore, several types of evaluation and benchmarking are required to understand models' decision-making processes and accurately compare the performance of models.

Differences between Omicron SARS-CoV-2 RBD and other variants in their ability to interact with cell receptors and monoclonal antibodies

SARS-CoV-2 remains a health threat with the continuous emergence of new variants. This work aims to expand the knowledge about the SARS-CoV-2 receptor-binding domain (RBD) interactions with cell receptors and monoclonal antibodies (mAbs). By using constant-pH Monte Carlo simulations, the free energy of interactions between the RBD from different variants and several partners (Angiotensin-Converting Enzyme-2 (ACE2) polymorphisms and various mAbs) were predicted. Computed RBD-ACE2-binding affinities were higher for two ACE2 polymorphisms (rs142984500 and rs4646116) typically found in Europeans which indicates a genetic susceptibility. This is amplified for Omicron (BA.1) and its sublineages BA.2 and BA.3. The antibody landscape was computationally investigated with the largest set of mAbs so far in the literature. From the 32 studied binders, groups of mAbs were identified from weak to strong binding affinities (e.g. S2K146). These mAbs with strong binding capacity and especially their combination are amenable to experimentation and clinical trials because of their high predicted binding affinities and possible neutralization potential for current known virus mutations and a universal coronavirus.Communicated by Ramaswamy H. Sarma.

Additive energetic contributions of multiple peptide positions determine the relative promiscuity of viral and human sequences for PDZ domain targets

Protein-protein interactions that involve recognition of short peptides are critical in cellular processes. Protein-peptide interaction surface areas are relatively small and shallow, and there are often overlapping specificities in families of peptide-binding domains. Therefore, dissecting selectivity determinants can be challenging. PDZ domains are a family of peptide-binding domains located in several intracellular signaling and trafficking pathways. These domains are also directly targeted by pathogens, and a hallmark of many oncogenic viral proteins is a PDZ-binding motif. However, amidst sequences that target PDZ domains, there is a wide spectrum in relative promiscuity. For example, the viral HPV16 E6 oncoprotein recognizes over double the number of PDZ domain-containing proteins as the cystic fibrosis transmembrane conductance regulator (CFTR) in the cell, despite similar PDZ targeting-sequences and identical motif residues. Here, we determine binding affinities for PDZ domains known to bind either HPV16 E6 alone or both CFTR and HPV16 E6, using peptides matching WT and hybrid sequences. We also use energy minimization to model PDZ-peptide complexes and use sequence analyses to investigate this difference. We find that while the majority of single mutations had marginal effects on overall affinity, the additive effect on the free energy of binding accurately describes the selectivity observed. Taken together, our results describe how complex and differing PDZ interactomes can be programmed in the cell.

Network Pharmacology and Molecular Modeling to Elucidate the Potential Mechanism of Neem Oil against Acne vulgaris

Acne vulgaris is a common skin disorder with a complicated etiology. Papules, lesions, comedones, blackheads, and other skin lesions are common physical manifestations of Acne vulgaris, but the individual who has it also regularly has psychological repercussions. Natural oils are being utilized more and more to treat skin conditions since they have fewer negative effects and are expected to provide benefits. Using network pharmacology, this study aims to ascertain if neem oil has any anti-acne benefits and, if so, to speculate on probable mechanisms of action for such effects. The neem leaves (Azadirachta indica) were collected, verified, authenticated, and assigned a voucher number. After steam distillation was used to extract the neem oil, the phytochemical components of the oil were examined using gas chromatography-mass spectrometry (GC-MS). The components of the oil were computationally examined for drug-likeness using Lipinski's criteria. The Pharm Mapper service was used to anticipate the targets. Prior to pathway and protein-protein interaction investigations, molecular docking was performed to predict binding affinity. Neem oil was discovered to be a potential target for STAT1, CSK, CRABP2, and SYK genes in the treatment of Acne vulgaris. In conclusion, it was discovered that the neem oil components with PubChem IDs: ID_610088 (2-(1-adamantyl)-N-methylacetamide), ID_600826 (N-benzyl-2-(2-methyl-5-phenyl-3H-1,3,4-thiadiazol-2-yl)acetamide), and ID_16451547 (N-(3-methoxyphenyl)-2-(1-phenyltetrazol-5-yl)sulfanylpropanamide) have strong affinities for these drug targets and may thus be used as therapeutic agents in the treatment of acne.

A Highly Expressed Antennae Odorant-Binding Protein Involved in Recognition of Herbivore-Induced Plant Volatiles in Dastarcus helophoroides

Natural enemies such as parasitoids and parasites depend on sensitive olfactory to search for their specific hosts. Herbivore-induced plant volatiles (HIPVs) are vital components in providing host information for many natural enemies of herbivores. However, the olfactory-related proteins involved in the recognition of HIPVs are rarely reported. In this study, we established an exhaustive tissue and developmental expression profile of odorant-binding proteins (OBPs) from Dastarcus helophoroides, an essential natural enemy in the forestry ecosystem. Twenty DhelOBPs displayed various expression patterns in different organs and adult physiological states, suggesting a potential involvement in olfactory perception. In silico AlphaFold2-based modeling and molecular docking showed similar binding energies between six DhelOBPs (DhelOBP4, 5, 6, 14, 18, and 20) and HIPVs from Pinus massoniana. While in vitro fluorescence competitive binding assays showed only recombinant DhelOBP4, the most highly expressed in the antennae of emerging adults could bind to HIPVs with high binding affinities. RNAi-mediated behavioral assays indicated that DhelOBP4 was an essential functional protein for D. helophoroides adults recognizing two behaviorally attractive substances: p-cymene and γ-terpinene. Further binding conformation analyses revealed that Phe 54, Val 56, and Phe 71 might be the key binding sites for DhelOBP4 interacting with HIPVs. In conclusion, our results provide an essential molecular basis for the olfactory perception of D. helophoroides and reliable evidence for recognizing the HIPVs of natural enemies from insect OBPs' perspective.

Additive energetic contributions of multiple peptide positions determine the relative promiscuity of viral and human sequences for PDZ domain targets

Protein-protein interactions that include recognition of short sequences of amino acids, or peptides, are critical in cellular processes. Protein-peptide interaction surface areas are relatively small and shallow, and there are often overlapping specificities in families of peptide-binding domains. Therefore, dissecting selectivity determinants can be challenging. PDZ domains are an example of a peptide-binding domain located in several intracellular signaling and trafficking pathways, which form interactions critical for the regulation of receptor endocytic trafficking, tight junction formation, organization of supramolecular complexes in neurons, and other biological systems. These domains are also directly targeted by pathogens, and a hallmark of many oncogenic viral proteins is a PDZ-binding motif. However, amidst sequences that target PDZ domains, there is a wide spectrum in relative promiscuity. For example, the viral HPV16 E6 oncoprotein recognizes over double the number of PDZ domain-containing proteins as the cystic fibrosis transmembrane conductance regulator (CFTR) in the cell, despite similar PDZ targeting-sequences and identical motif residues. Here, we determine binding affinities for PDZ domains known to bind either HPV16 E6 alone or both CFTR and HPV16 E6, using peptides matching WT and hybrid sequences. We also use energy minimization to model PDZ-peptide complexes and use sequence analyses to investigate this difference. We find that while the majority of single mutations had a marginal effect on overall affinity, the additive effect on the free energy of binding accurately describes the selectivity observed. Taken together, our results describe how complex and differing PDZ interactomes can be programmed in the cell.

Comparative Genomic Characterization of Buffalo Fibronectin Type III Domain Proteins: Exploring the Novel Role of FNDC5/Irisin as a Ligand of Gonadal Receptors

FN-III proteins are widely distributed in mammals and are usually involved in cellular growth, differentiation, and adhesion. The FNDC5/irisin regulates energy metabolism and is present in different tissues (liver, brain, etc.). The present study aimed to investigate the physiochemical characteristics and the evolution of FN-III proteins and FNDC5/irisin as a ligand targeting the gonadal receptors including androgen (AR), DDB1 and CUL4 associated factor 6 (DCAF6), estrogen-related receptor β (ERR-β), estrogen-related receptor γ (ERR-γ), Krüppel-like factor 15 (KLF15), and nuclear receptor subfamily 3 group C member 1 (NR3C1). Moreover, the putative role of irisin in folliculogenesis and spermatogenesis was also elucidated. We presented the molecular structure and function of 29 FN-III genes widely distributed in the buffalo genome. Phylogenetic analysis, motif, and conserved domain pattern demonstrated the evolutionary well-conserved nature of FN-III proteins with a variety of stable to unstable, hydrophobic to hydrophilic, and thermostable to thermo-unstable properties. The comparative structural configuration of FNDC5 revealed amino acid variations but still the FNDC5 structure of humans, buffalo, and cattle was quite similar to each other. For the first time, we predicted the binding scores and interface residues of FNDC5/irisin as a ligand for six representative receptors having a functional role in energy homeostasis, and a significant involvement in folliculogenesis and spermatogenesis in buffalo.

MRBLE-pep Measurements Reveal Accurate Binding Affinities for B56, a PP2A Regulatory Subunit

Signal transduction pathways rely on dynamic interactions between protein globular domains and short linear motifs (SLiMs). The weak affinities of these interactions are essential to allow fast rewiring of signaling pathways and downstream responses but also pose technical challenges for interaction detection and measurement. We recently developed a technique (MRBLE-pep) that leverages spectrally encoded hydrogel beads to measure binding affinities between a single protein of interest and 48 different peptide sequences in a single small volume. In prior work, we applied it to map the binding specificity landscape between calcineurin and the PxIxIT SLiM (Nguyen, H. Q. et al. Elife 2019, 8). Here, using peptide sequences known to bind the PP2A regulatory subunit B56α, we systematically compare affinities measured by MRBLE-pep or isothermal calorimetry (ITC) and confirm that MRBLE-pep accurately quantifies relative affinity over a wide dynamic range while using a fraction of the material required for traditional methods such as ITC.

Boundary values of binding constants determined by affinity capillary electrophoresis

This study shows that the upper limit of binding (stability) constants determined by mobility shift affinity capillary electrophoresis can be increased from 10⁴ to 10⁶ -10⁹ L/mol if the lowest possible analyte concentration in samples is used (for example, the concentration that gives electrophoretic peaks with a signal-to-noise ratio of 10) and the effective electrophoretic mobility of the analyte is calculated via the parameter a₁ of the Haarhoff-Van der Linde function. The equation to calculate the boundary values of binding constants for 1:1 complexes was derived for the case when the constants cannot be calculated in the usual way. These values are obtained from the inequality: the difference between the ionic mobility of the analyte-ligand complex and the effective electrophoretic mobility of the analyte determined at the lowest ligand concentration in the background electrolyte at which the analyte appears as an undistorted peak in electropherograms is less than or equal to the absolute error in mobility measurements. The application of the equation was illustrated by the example of electrophoretic data for a complex between betulin 3,28-diphthalate and (2-hydroxypropyl)-γ-cyclodextrin. An algorithm to determine the binding constants for strong complexation by mobility shift affinity capillary electrophoresis was suggested.

Structural characterization and computational analysis of PDZ domains in Monosiga brevicollis

Identification of the molecular networks that facilitated the evolution of multicellular animals from their unicellular ancestors is a fundamental problem in evolutionary cellular biology. Choanoflagellates are recognized as the closest extant nonmetazoan ancestors to animals. These unicellular eukaryotes can adopt a multicellular-like "rosette" state. Therefore, they are compelling models for the study of early multicellularity. Comparative studies revealed that a number of putative human orthologs are present in choanoflagellate genomes, suggesting that a subset of these genes were necessary for the emergence of multicellularity. However, previous work is largely based on sequence alignments alone, which does not confirm structural nor functional similarity. Here, we focus on the PDZ domain, a peptide-binding domain which plays critical roles in myriad cellular signaling networks and which underwent a gene family expansion in metazoan lineages. Using a customized sequence similarity search algorithm, we identified 178 PDZ domains in the Monosiga brevicollis proteome. This includes 11 previously unidentified sequences, which we analyzed using Rosetta and homology modeling. To assess conservation of protein structure, we solved high-resolution crystal structures of representative M. brevicollis PDZ domains that are homologous to human Dlg1 PDZ2, Dlg1 PDZ3, GIPC, and SHANK1 PDZ domains. To assess functional conservation, we calculated binding affinities for mbGIPC, mbSHANK1, mbSNX27, and mbDLG-3 PDZ domains from M. brevicollis. Overall, we find that peptide selectivity is generally conserved between these two disparate organisms, with one possible exception, mbDLG-3. Overall, our results provide novel insight into signaling pathways in a choanoflagellate model of primitive multicellularity.

Tomato Chlorosis Virus Minor Coat Protein as a Novel Target To Screen Antiviral Drugs

Minor coat protein (mCP), an important component of tomato chlorosis virus (ToCV), plays a significant role in the process of virus assembly and movement and is directly related to the virus-insect transmission. Therefore, ToCV mCP could be considered as a potent target for anti-ToCV drugs. In this study, ToCV mCP was first cloned, expressed, purified, and a novel target to screen the antiviral agents. The results showed that some antiviral compounds bound to ToCV mCP with strongly affinities in vitro, including quinazoline derivatives 4a and 4b, Ningnanmycin, and Ribavirin. Subsequently, three-dimensional-quantitative structure-activity relationship (3D-QSAR) analysis was performed based on the binding affinities, and the model indicated that 4a and 4b had indeed stronger binding effects on ToCV mCP than other quinazoline derivatives. Finally, the anti-ToCV activities of compounds 4a and 4b were evaluated by quantitative real-time polymerase chain reaction in vivo. Compounds 4a and 4b inhibited infection of ToCV in the host and as well as reduced the level of ToCV mCP gene expression. Thus, ToCV mCP can be used as a novel drug target for screening anti-ToCV agents, and the ligand-based 3D-QSAR analysis of quinazoline derivatives provided new insights into the design and optimization of novel anti-ToCV drug molecules based on ToCV mCP.

Updates to Binding MOAD (Mother of All Databases): Polypharmacology Tools and Their Utility in Drug Repurposing

The goal of Binding MOAD is to provide users with a data set focused on high-quality x-ray crystal structures that have been solved with biologically relevant ligands bound. Where available, experimental binding affinities (K_a, K_d, K_i, IC₅₀) are provided from the primary literature of the crystal structure. The database has been updated regularly since 2005, and this most recent update has added nearly 7000 new structures (growth of 21%). MOAD currently contains 32,747 structures, composed of 9117 protein families and 16,044 unique ligands. The data are freely available on www.BindingMOAD.org. This paper outlines updates to the data in Binding MOAD as well as improvements made to both the website and its contents. The NGL viewer has been added to improve visualization of the ligands and protein structures. MarvinJS has been implemented, over the outdated MarvinView, to work with JChem for small molecule searching in the database. To add tools for predicting polypharmacology, we have added information about sequence, binding-site, and ligand similarity between entries in the database. A main premise behind polypharmacology is that similar binding sites will bind similar ligands. The large amount of protein-ligand information available in Binding MOAD allows us to compute pairwise ligand and binding-site similarities. Lists of similar ligands and similar binding sites have been added to allow users to identify potential polypharmacology pairs. To show the utility of the polypharmacology data, we detail a few examples from Binding MOAD of drug repurposing targets with their respective similarities.

Analysis of α-Dystroglycan/LG Domain Binding Modes: Investigating Protein Motifs That Regulate the Affinity of Isolated LG Domains

Dystroglycan (DG) is an adhesion complex that links the cytoskeleton to the surrounding extracellular matrix in skeletal muscle and a wide variety of other tissues. It is composed of a highly glycosylated extracellular α-DG associated noncovalently with a transmembrane β-DG whose cytodomain interacts with dystrophin and its isoforms. Alpha-dystroglycan (α-DG) binds tightly and in a calcium-dependent fashion to multiple extracellular proteins and proteoglycans, each of which harbors at least one, or, more frequently, tandem arrays of laminin-globular (LG) domains. Considerable biochemical and structural work has accumulated on the α-DG-binding LG domains, highlighting a significant heterogeneity in ligand-binding properties of domains from different proteins as well as between single and multiple LG domains within the same protein. Here we review biochemical, structural, and functional information on the LG domains reported to bind α-dystroglycan. In addition, we have incorporated bioinformatics and modeling to explore whether specific motifs responsible for α-dystroglycan recognition can be identified within isolated LG domains. In particular, we analyzed the LG domains of slits and agrin as well as those of paradigmatic α-DG non-binders such as laminin-α3. While some stretches of basic residues may be important, no universally conserved motifs could be identified. However, the data confirm that the coordinated calcium atom within the LG domain is needed to establish an interaction with the sugars of α-DG, although it appears that this alone is insufficient to mediate significant α-DG binding. We develop a scenario involving different binding modes of a single LG domain unit, or tandemly repeated units, with α-DG. A variability of binding modes might be important to generate a range of affinities to allow physiological regulation of this interaction, reflecting its crucial biological importance.

"Methylene Bridge" to 5-HT3 Receptor Antagonists: Conformationally Constrained Phenylguanidines

Arylguanidines, depending upon their aromatic substitution pattern, display varying actions at 5-HT₃ receptors (e.g., partial agonist, agonist, superagonist). Here, we demonstrate that conformational constraint of these agents as dihydroquinazolines (such as A6CDQ; 1) results in their conversion to 5-HT₃ receptor antagonists. We examined the structure-activity relationships of 1. Replacement/removal of any of the guanidinium nitrogen atoms of 1 resulted in decreased affinity. All three nitrogen atoms of 1 are necessary for optimal binding affinity at 5-HT₃ receptors. Introduction of substituents as small as an N2-methyl group abolishes affinity. The results are consistent with homology modeling/docking studies and binding data from site-directed mutagenesis studies. Introducing a "methylene bridge" to the arylguanidine structure, regardless of its functional activity, results in a 5-HT₃ receptor antagonist.

N1-Substituted Quinoxaline-2,3-diones as Kainate Receptor Antagonists: X-ray Crystallography, Structure-Affinity Relationships, and in Vitro Pharmacology

Among the ionotropic glutamate receptors, the physiological role of kainate receptors is less well understood. Although ligands with selectivity toward the kainate receptor subtype GluK1 are available, tool compounds with selectivity at the remaining kainate receptor subtypes are sparse. Here, we have synthesized a series of quinoxaline-2,3-diones with substitutions in the N1-, 6-, and 7-position to investigate the structure-activity relationship (SAR) at GluK1-3 and GluK5. Pharmacological characterization at native and recombinant kainate and AMPA receptors revealed that compound 37 had a GluK3-binding affinity ( K_i) of 0.142 μM and 8-fold preference for GluK3 over GluK1. Despite lower binding affinity of 22 at GluK3 ( K_i = 2.91 μM), its preference for GluK3 over GluK1 and GluK2 was >30-fold. Compound 37 was crystallized with the GluK1 ligand-binding domain to understand the SAR. The X-ray structure showed that 37 stabilized the protein in an open conformation, consistent with an antagonist binding mode.

Ranking Power of the SQM/COSMO Scoring Function on Carbonic Anhydrase II-Inhibitor Complexes

Accurate prediction of protein-ligand binding affinities is essential for hit-to-lead optimization and virtual screening. The reliability of scoring functions can be improved by including quantum effects. Here, we demonstrate the ranking power of the semiempirical quantum mechanics (SQM)/implicit solvent (COSMO) scoring function by using a challenging set of 10 inhibitors binding to carbonic anhydrase II through Zn²⁺ in the active site. This new dataset consists of the high-resolution (1.1-1.4 Å) crystal structures and experimentally determined inhibitory constant (K_i ) values. It allows for evaluation of the common approximations, such as representing the solvent implicitly or by using a single target conformation combined with a set of ligand docking poses. SQM/COSMO attained a good correlation of R² of 0.56-0.77 with the experimental inhibitory activities, benefiting from careful handling of both noncovalent interactions (e.g. charge transfer) and solvation. This proof-of-concept study of SQM/COSMO ranking for metalloprotein-ligand systems demonstrates its potential for hit-to-lead applications.

Systematic study on metabolism and activity evaluation of Radix Scutellaria extract in rat plasma using UHPLC with quadrupole time-of-flight mass spectrometry and microdialysis intensity-fading mass spectrometry

Radix Scutellaria is a widely used traditional Chinese medicine in the treatment of various diseases. However, the activities of the absorbed components and metabolites of its main flavones in rat plasma need further investigation. In this study, a systematic method based on ultra-high performance liquid chromatography with quadruple time-of-flight mass spectrometry was developed to speculate the absorbed components and metabolites of the main flavonoids in Radix Scutellaria extract in rat plasma sample after oral administration of the extract. Twelve compounds, including four prototype components and eight metabolites, were confirmed in drug-containing plasma. In these metabolites, five were originally detected in rat plasma. The possible metabolic pathways of these polyhydroxy flavones in vivo were described and clarified. Microdialysis with intensity-fading mass spectrometry was originally employed to investigate the binding affinities of the absorbed components and metabolites with α-glucosidase. The order of their binding affinities was P4 > P3 > P2 > P1≥M5 > M3 > M1. The research result is helpful to deepen the understanding of the absorbed components and metabolic pathways of main flavones from Radix Scutellaria, and provide a new approach to screen potential inhibitors from in vivo components originated from Chinese herb.

Structure and Affinity of Two Bicyclic Glutamate Analogues at AMPA and Kainate Receptors

Ionotropic glutamate receptors (iGluRs) are involved in most of the fast excitatory synaptic transmission in the central nervous system. These receptors are important for learning and memory formation, but are also involved in the development of diseases such as Alzheimer's disease, epilepsy and depression. To understand the function of different types of iGluRs, selective agonists are invaluable as pharmacological tool compounds. Here, we report binding affinities of two bicyclic, conformationally restricted analogues of glutamate (CIP-AS and LM-12b) at AMPA (GluA2 and GluA3) and kainate receptor subunits (GluK1-3 and GluK5). Both CIP-AS and LM-12b were found to be GluK3-preferring agonists, with K_i of 6 and 22 nM, respectively, at recombinant GluK3 receptors. The detailed binding mode of CIP-AS and LM-12b in the ligand-binding domains of the AMPA receptor subunit GluA2 (GluA2-LBD) and the kainate receptor subunits GluK1 (GluK1-LBD) and GluK3 (GluK3-LBD) was investigated by X-ray crystallography. CIP-AS stabilized all three receptor constructs in conformations similar to those with kainate. Remarkably, whereas LM-12b bound in a similar manner to CIP-AS in GluA2-LBD and GluK3-LBD, it introduced full closure of the ligand-binding domain in GluK1-LBD and formation of a D1-D2 interlobe hydrogen bond between Glu441 and Ser721, as also observed with glutamate. As the binding affinity of LM-12b at GluK1 is ∼8-fold better than that for CIP-AS (K_i of 85 and 656 nM, respectively), it shows that small changes in agonist structure can lead to prominent differences in structure and function.

Superagonist, Full Agonist, Partial Agonist, and Antagonist Actions of Arylguanidines at 5-Hydroxytryptamine-3 (5-HT3) Subunit A Receptors

Introduction of minor variations to the substitution pattern of arylguanidine 5-hydroxytryptamine-3 (5-HT₃) receptor ligands resulted in a broad spectrum of functionally-active ligands from antagonist to superagonist. For example, (i) introduction of an additional Cl-substituent(s) to our lead full agonist N-(3-chlorophenyl)guanidine (mCPG, 2; efficacy % = 106) yielded superagonists 7-9 (efficacy % = 186, 139, and 129, respectively), (ii) a positional isomer of 2, p-Cl analog 11, displayed partial agonist actions (efficacy % = 12), and (iii) replacing the halogen atom at the meta or para position with an electron donating OCH₃ group or a stronger electron withdrawing (i.e., CF₃) group resulted in antagonists 13-16. We posit based on combined mutagenesis, crystallographic, and computational analyses that for the 5-HT₃ receptor, the arylguanidines that are better able to simultaneously engage the primary and complementary subunits, thus keeping them in close proximity, have greater agonist character while those that are deficient in this ability are antagonists.

Is It Reliable to Use Common Molecular Docking Methods for Comparing the Binding Affinities of Enantiomer Pairs for Their Protein Target?

Molecular docking is a computational chemistry method which has become essential for the rational drug design process. In this context, it has had great impact as a successful tool for the study of ligand–receptor interaction modes, and for the exploration of large chemical datasets through virtual screening experiments. Despite their unquestionable merits, docking methods are not reliable for predicting binding energies due to the simple scoring functions they use. However, comparisons between two or three complexes using the predicted binding energies as a criterion are commonly found in the literature. In the present work we tested how wise is it to trust the docking energies when two complexes between a target protein and enantiomer pairs are compared. For this purpose, a ligand library composed by 141 enantiomeric pairs was used, including compounds with biological activities reported against seven protein targets. Docking results using the software Glide (considering extra precision (XP), standard precision (SP), and high-throughput virtual screening (HTVS) modes) and AutoDock Vina were compared with the reported biological activities using a classification scheme. Our test failed for all modes and targets, demonstrating that an accurate prediction when binding energies of enantiomers are compared using docking may be due to chance. We also compared pairs of compounds with different molecular weights and found the same results.

Enzymatic Synthesis of Galactosylated Serine/Threonine Derivatives by β-Galactosidase from Escherichia coli

The transgalactosylations of serine/threonine derivatives were investigated using β-galactosidase from Escherichia coli as biocatalyst. Using ortho-nitrophenyl-β-d-galactoside as donor, the highest bioconversion yield of transgalactosylated N-carboxy benzyl l-serine benzyl ester (23.2%) was achieved in heptane:buffer medium (70:30), whereas with the lactose, the highest bioconversion yield (3.94%) was obtained in the buffer reaction system. The structures of most abundant galactosylated serine products were characterized by MS/MS. The molecular docking simulation revealed that the binding of serine/threonine derivatives to the enzyme’s active site was stronger (−4.6~−7.9 kcal/mol) than that of the natural acceptor, glucose, and mainly occurred through interactions with aromatic residues. For N-tert-butoxycarbonyl serine methyl ester (6.8%) and N-carboxybenzyl serine benzyl ester (3.4%), their binding affinities and the distances between their hydroxyl side chain and the 1′-OH group of galactose moiety were in good accordance with the quantified bioconversion yields. Despite its lower predicted bioconversion yield, the high experimental bioconversion yield obtained with N-carboxybenzyl serine methyl ester (23.2%) demonstrated the importance of the thermodynamically-driven nature of the transgalactosylation reaction.

Mechanically Induced Trapping of Molecular Interactions and Its Applications

Measuring binding affinities and association/dissociation rates of molecular interactions is important for a quantitative understanding of cellular mechanisms. Many low-throughput methods have been developed throughout the years to obtain these parameters. Acquiring data with higher accuracy and throughput is, however, necessary to characterize complex biological networks. Here, we provide an overview of a high-throughput microfluidic method based on mechanically induced trapping of molecular interactions (MITOMI). MITOMI can be used to obtain affinity constants and kinetic rates of hundreds of protein-ligand interactions in parallel. It has been used in dozens of studies to measure binding affinities of transcription factors, map protein interaction networks, identify pharmacological inhibitors, and perform high-throughput, low-cost molecular diagnostics. This article covers the technological aspects of MITOMI and its applications.

Anticholinergics in palliative medicine: an update

Anticholinergics, or antimuscarinic drugs, are drugs that competitively inhibit the action of acetylcholine at muscarinic receptors, leading to a blockade of the actions of the parasympathetic nervous system at sites where overactivity can lead to increased symptom burden. Successful blockade of the parasympathetic nervous system ultimately leads to decreased production of secretions in the salivary, bronchial, and gastrointestinal tracts. These effects are often used for several symptoms that originate due to parasympathetic nervous system overactivity, such as the "death rattle" and malignant bowel obstruction. Anticholinergic agents are divided into either tertiary amines or quaternary ammonium compounds, which differ in their ability to cross into the central nervous system. Quaternary compounds do not cross into the central nervous system and have a different adverse effect profile than the tertiary amines. The purpose of this review is to highlight anticholinergic agents, their pharmacology, and an evidence-based assessment of their role in palliative care.

Iron-binding activity in yeast frataxin entails a trade off with stability in the alpha1/beta1 acidic ridge region

Frataxin is a highly conserved mitochondrial protein whose deficiency in humans results in Friedreich's ataxia (FRDA), an autosomal recessive disorder characterized by progressive ataxia and cardiomyopathy. Although its cellular function is still not fully clear, the fact that frataxin plays a crucial role in Fe-S assembly on the scaffold protein Isu is well accepted. In the present paper, we report the characterization of eight frataxin variants having alterations on two putative functional regions: the alpha1/beta1 acidic ridge and the conserved beta-sheet surface. We report that frataxin iron-binding capacity is quite robust: even when five of the most conserved residues from the putative iron-binding region are altered, at least two iron atoms per monomer can be bound, although with decreased affinity. Furthermore, we conclude that the acidic ridge is designed to favour function over stability. The negative charges have a functional role, but at the same time significantly impair frataxin's stability. Removing five of those charges results in a thermal stabilization of approximately 24 degrees C and reduces the inherent conformational plasticity. Alterations on the conserved beta-sheet residues have only a modest impact on the protein stability, highlighting the functional importance of residues 122-124.

Site-specific Eu(III) binding affinities to a Datura innoxia biosorbent

The binding of Eu(III) to a biosorbent derived from cultured cells of the plant Datura innoxia, have been investigated through elucidation of apparent affinity constants associated with different chemical environments present on the cell wall. Adsorption isotherms for separate types of binding sites were generated using metal ion luminescence measurements. Application of regularized regression analysis to these isotherm data for four chemically distinguishable sites revealed the presence of sites exhibiting both low (mean log K(app)=-0.3 to 0.6) and higher (mean log K(app)=3.2-3.5) apparent affinities for pH conditions of 2.0, 4.0, and 5.0. Low affinity sites were observed for all pH conditions and attributed to non-specific binding of the metal ions to the negatively charged biomaterial. The pH-dependent higher affinity sites are ascribed to specific sites involving either an ion-exchange mechanism or formation of weak surface-metal ion complexes. These results differed significantly from a similar analysis of total metal binding isotherms that indicated mean log K(app) values of -0.5 to 0.25 (low affinity) and 5.6-6.0 (high affinity).

The role of multiple hydrogen-bonding groups in specific alcohol binding sites in proteins: insights from structural studies of LUSH

It is now generally accepted that many of the physiological effects of alcohol consumption are a direct result of binding to specific sites in neuronal proteins such as ion channels or other components of neuronal signaling cascades. Binding to these targets generally occurs in water-filled pockets and leads to alterations in protein structure and dynamics. However, the precise interactions required to confer alcohol sensitivity to a particular protein remain undefined. Using information from the previously solved crystal structures of the Drosophila melanogaster protein LUSH in complexes with short-chain alcohols, we have designed and tested the effects of specific amino acid substitutions on alcohol binding. The effects of these substitutions, specifically S52A, T57S, and T57A, were examined using a combination of molecular dynamics, X-ray crystallography, fluorescence spectroscopy, and thermal unfolding. These studies reveal that the binding of ethanol is highly sensitive to small changes in the composition of the alcohol binding site. We find that T57 is the most critical residue for binding alcohols; the T57A substitution completely abolishes binding, while the T57S substitution differentially affects ethanol binding compared to longer-chain alcohols. The additional requirement for a potential hydrogen-bond acceptor at position 52 suggests that both the presence of multiple hydrogen-bonding groups and the identity of the hydrogen-bonding residues are critical for defining an ethanol binding site. These results provide new insights into the detailed chemistry of alcohol's interactions with proteins.

Quantitative autoradiographic studies of relaxin binding in rat atria, uterus and cerebral cortex: characterization and effects of oestrogen treatment

The binding characteristics of the relaxin receptor in rat atria, uterus and cortex were studied using a [33P]-labelled human gene 2 relaxin (B33) and quantitative receptor autoradiography. The binding kinetics of [33P]-human gene 2 relaxin (B33) were investigated in slide-mounted rat atrial sections. The binding achieved equilibrium after 60 min incubation at room temperature (23+/-1 degrees C) and dissociated slowly. The association and dissociation rate constants were 4.31+/-0.34x10(8) M(-1) x min(-1) and 1.55+/-0.38x10(-3) min(-1) respectively. Thus, the kinetic dissociation constant was 3.46+/-0.59 pM. Binding was saturable to a single population of non-interacting sites throughout atria, in uterine myometrium and the 5th layer of cerebral cortex. The binding affinities (pK(D)) of [33P]-human gene 2 relaxin (B33) were 8.92+/-0.09 in atrial myocardium and 8.79+/-0.04 in cerebral cortex of male rats, and 8.79+/-0.10 in uterine myometrium. Receptor densities in the cerebral cortex and atria were higher than in uterine myometrium, indicating that relaxin also has important roles in non-reproductive tissues. In male rats, treatment with 17beta-oestradiol (20 microg in 0.1 ml sesame oil s.c., 18-24 h) significantly decreased the density of relaxin receptors in atria and cerebral cortex. Identical treatment in female rats had no significant effect in atria and cerebral cortex, but it significantly increased the density of relaxin receptors in uterine myometrium. Relaxin binding was competitively displaced by porcine and rat native relaxins. Porcine native relaxin binds to the relaxin receptor in male rat atria (8.90+/-0.02), and cerebral cortex (8.90+/-0.03) and uterine myometrium (8.89+/-0.03) with affinities not significantly different from human gene 2 (B33) relaxin. Nevertheless, rat relaxin binds to the receptors with affinities (8.35+/-0.09 in atria, 8.22+/-0.07 in cerebral cortex and 8.48+/-0.06 in uterine myometrium) significantly less than human gene 2 (B33) and porcine relaxins. Quantitative receptor autoradiography is the method of choice for measurement of affinities and densities of relaxin receptor in atria, uterine myometrium and cerebral cortex. High densities were found in all these tissues. 17beta-oestradiol treatment produced complex effects where it increased the densities of relaxin receptors in uterus but decreased those in atria and cerebral cortex of the male rats, and had no effect on the atria and cerebral cortex of the female rats.