miPEP31 alleviates Ang II-induced hypertension in mice by occupying Cebpα binding sites in the pri-miR-31 promoter.

BACKGROUND: Previous studies have shown that peptides encoded by noncoding RNAs (ncRNAs) can be used as peptide drugs to alleviate diseases. We found that microRNA-31 (miR-31) is involved in the regulation of hypertension and that the peptide miPEP31, which is encoded by the primary transcript of miR-31 (pri-miR-31), can inhibit miR-31 expression. However, the role and mechanism of miPEP31 in hypertension have not been elucidated. METHODS: miPEP31 expression was determined by western blot analysis. miPEP31-deficient mice (miPEP31-/-) were used, and synthetic miPEP31 was injected into Ang II-induced hypertensive mice. Blood pressure was monitored through the tail-cuff method. Histological staining was used to evaluate renal damage. Regulatory T (Treg) cells were assessed by flow cytometry. Differentially expressed genes were analysed through RNA sequencing. The transcription factors were predicted by JASPAR. Luciferase reporter and electrophoretic mobility shift assays (EMSAs) were used to determine the effect of pri-miR-31 on the promoter activity of miPEP31. Images were taken to track the entry of miPEP31 into the cell. RESULTS: miPEP31 is endogenously expressed in target organs and cells related to hypertension. miPEP31 deficiency exacerbated but exogenous miPEP31 administration mitigated the Ang II-induced systolic blood pressure (SBP) elevation, renal impairment and Treg cell decreases in the kidney. Moreover, miPEP31 deletion increased the expression of genes related to Ang II-induced renal fibrosis. miPEP31 inhibited the transcription of miR-31 and promoted Treg differentiation by occupying the Cebpα binding site. The minimal functional domain of miPEP31 was identified and shown to regulate miR-31. CONCLUSION: miPEP31 was identified as a potential therapeutic peptide for treating hypertension by promoting Treg cell differentiation in vivo. Mechanistically, we found that miPEP31 acted as a transcriptional repressor to specifically inhibit miR-31 transcription by competitively occupying the Cebpα binding site in the pri-miR-31 promoter. Our study highlights the significant therapeutic effect of miPEP31 on hypertension and provides novel insight into the role and mechanism of miPEPs.

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction.

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

Computational Investigation of BMAA and Its Carbamate Adducts as Potential GluR2 Modulators.

Beta-N-methylamino-l-alanine (BMAA) is a potential neurotoxic nonprotein amino acid, which can reach the human body through the food chain. When BMAA interacts with bicarbonate in the human body, carbamate adducts are produced, which share a high structural similarity with the neurotransmitter glutamate. It is believed that BMAA and its l-carbamate adducts bind in the glutamate binding site of ionotropic glutamate receptor 2 (GluR2). Chronic exposure to BMAA and its adducts could cause neurological illness such as neurodegenerative diseases. However, the mechanism of BMAA action and its carbamate adducts bound to GluR2 has not yet been elucidated. Here, we investigate the binding modes and the affinity of BMAA and its carbamate adducts to GluR2 in comparison to the natural agonist, glutamate, to understand whether these can act as GluR2 modulators. Initially, we perform molecular dynamics simulations of BMAA and its carbamate adducts bound to GluR2 to examine the stability of the ligands in the S1/S2 ligand-binding core of the receptor. In addition, we utilize alchemical free energy calculations to compute the difference in the free energy of binding of the beta-carbamate adduct of BMAA to GluR2 compared to that of glutamate. Our findings indicate that carbamate adducts of BMAA and glutamate remain stable in the binding site of the GluR2 compared to BMAA. Additionally, alchemical free energy results reveal that glutamate and the beta-carbamate adduct of BMAA have comparable binding affinity to the GluR2. These results provide a rationale that BMAA carbamate adducts may be, in fact, the modulators of GluR2 and not BMAA itself.

Bitopic Ligands Support the Presence of a Metastable Binding Site at the ß2 Adrenergic Receptor.

Metastable binding sites (MBS) have been observed in a multitude of molecular dynamics simulations and can be considered low affinity allosteric binding sites (ABS) that function as stepping stones as the ligand moves toward the orthosteric binding site (OBS). Herein, we show that MBS can be utilized as ABS in ligand design, resulting in ligands with improved binding kinetics. Four homobivalent bitopic ligands (1-4) were designed by molecular docking of (S)-alprenolol ((S)-ALP) in the cocrystal structure of the ß2 adrenergic receptor (ß2AR) bound to the antagonist ALP. Ligand 4 displayed a potency and affinity similar to (S)-ALP, but with a >4-fold increase in residence time. The proposed binding mode was confirmed by X-ray crystallography of ligand 4 in complex with the ß2AR. This ligand design principle can find applications beyond the ß2AR and G protein-coupled receptors (GPCRs) as a general approach for improving the pharmacological profile of orthosteric ligands by targeting the OBS and an MBS simultaneously.

Crystal Structure, Theoretical Analysis, and Protein/DNA Binding Activity of Iron(III) Complex Containing Differently Protonated Pyridoxal-S-Methyl-Isothiosemicarbazone Ligands.

Pyridoxal-S-methyl-isothiosemicarbazone (PLITSC) is a member of an important group of ligands characterized by different complexation modes to various transition metals. In this contribution, a new complex containing two differently protonated PLITSC ligands ([Fe(PLITSC-H)(PLITSC)]SO4)∙2.5H2O was obtained. The crystal structure was solved by the X-ray analysis and used further for the optimization at B3LYP/6-311++G(d,p)(H,C,N,O,S)/def2-TZVP(Fe) level of theory. Changes in the interaction strength and bond distance due to protonation were observed upon examination by the Quantum Theory of Atoms in Molecules. The protein binding affinity of [Fe(PLITSC-H)(PLITSC)]SO4 towards transport proteins (Bovine Serum Albumin (BSA) and Human Serum Albumin (HSA)) was investigated by the spectrofluorimetric titration and molecular docking. The interactions with the active pocket containing fluorescent amino acids were examined in detail, which explained the fluorescence quenching. The interactions between complex and DNA were followed by the ethidium-bromide displacement titration and molecular docking. The binding along the minor groove was the dominant process involving complex in the proximity of DNA.

Binding Affinity Determination in Drug Design: Insights from Lock and Key, Induced Fit, Conformational Selection, and Inhibitor Trapping Models.

Binding affinity is a fundamental parameter in drug design, describing the strength of the interaction between a molecule and its target protein. Accurately predicting binding affinity is crucial for the rapid development of novel therapeutics, the prioritization of promising candidates, and the optimization of their properties through rational design strategies. Binding affinity is determined by the mechanism of recognition between proteins and ligands. Various models, including the lock and key, induced fit, and conformational selection, have been proposed to explain this recognition process. However, current computational strategies to predict binding affinity, which are based on these models, have yet to produce satisfactory results. This article explores the connection between binding affinity and these protein-ligand interaction models, highlighting that they offer an incomplete picture of the mechanism governing binding affinity. Specifically, current models primarily center on the binding of the ligand and do not address its dissociation. In this context, the concept of ligand trapping is introduced, which models the mechanisms of dissociation. When combined with the current models, this concept can provide a unified theoretical framework that may allow for the accurate determination of the ligands' binding affinity.

The impact of von Willebrand factor on fibrin formation and structure unveiled with type 3 von Willebrand disease plasma.

Normally, von Willebrand factor (VWF) remains inactive unless its A1A2 domains undergo a shear stress-triggered conformational change. We demonstrated the capacity of a recombinant A2 domain of VWF to bind and to affect fibrin formation, altering the fibrin clot structure. The data indicated that VWF contains an additional binding site for fibrin in the A2 domain that plays a role in the incorporation of VWF to the polymerizing fibrin. This study is to examine the hypothesis that active plasma VWF directly influence fibrin polymerization and the structure of fibrin clots. The study used healthy and type 3 von Willebrand disease (VWD) plasma, purified plasma VWF, fibrin polymerization assays, confocal microscopy and scanning electron microscopy. The exposed A2 domain in active VWF harbors additional binding sites for fibrinogen, and significantly potentiates fibrin formation (P < 0.02). Antibody against the A2 domain of VWF significantly decreased the initial rate of change of fibrin formation (P < 0.002). Clot analyses revealed a significant difference in porosity between normal and type 3 VWD plasma (P < 0.008), further supported by scanning electron microscopy, which demonstrated thicker fibrin fibers in the presence of plasma VWF (P < 0.0003). Confocal immunofluorescence microscopy showed punctate VWF staining along fibrin fibrils, providing visual evidence of the integration of plasma VWF into the fibrin matrix. The study with type 3 VWD plasma supports the hypothesis that plasma VWF directly influences fibrin polymerization and clot structure. In addition, a conformational change in the A1A2 domains exposes a hidden fibrin(ogen) binding site, indicating that plasma VWF determines the fibrin clot structure.

From Multiple Protein Docking to Protein-Protein Docking at Interactome Level.

Molecular docking is used to anticipate the optimal orientation of a particular molecule to a target to form a stable complex. It makes predictions about the 3D structure of any complex based on the binding characteristics of the ligand and the target receptor usually a protein. It is an exceptionally useful tool, which is used as a model to study how ligands attach to proteins. Docking can also be used for studying the interaction of ligands and proteins to analyze inhibitory efficacy. The ligand may also be a protein, making it possible to study interactions between two different proteins using the numerous docking tools available for basic research on protein interactions. The protein-protein docking is a crucial approach to understanding the protein interactions and predicting the structure of protein complexes that have not yet been experimentally determined. Moreover, the protein-protein interactions can predict the function of target proteins and the drug-like properties of molecules. Therefore, protein docking assists in uncovering insights into protein interactions and also aids in a better understanding of molecular pathways/mechanisms. This chapter comprehends the various tools for protein-protein docking (pairwise and multiple), including their methodologies and analysis of output as results.

TFscope: systematic analysis of the sequence features involved in the binding preferences of transcription factors.

Characterizing the binding preferences of transcription factors (TFs) in different cell types and conditions is key to understand how they orchestrate gene expression. Here, we develop TFscope, a machine learning approach that identifies sequence features explaining the binding differences observed between two ChIP-seq experiments targeting either the same TF in two conditions or two TFs with similar motifs (paralogous TFs). TFscope systematically investigates differences in the core motif, nucleotide environment and co-factor motifs, and provides the contribution of each key feature in the two experiments. TFscope was applied to > 305 ChIP-seq pairs, and several examples are discussed.

Identification of a magnesium-binding site at the primary allosteric calcium sensor of the sodium-calcium exchanger: Implications for physiological regulation.

Sodium-calcium exchanger (NCX) proteins are ubiquitously expressed and play a pivotal role in cellular calcium homeostasis by mediating uphill calcium efflux across the cell membrane. Intracellular calcium allosterically regulates the exchange activity by binding to two cytoplasmic calcium-binding domains, CBD1 and CBD2. However, the calcium-binding affinities of these domains are seemingly inadequate to sense physiological calcium oscillations. Previously, magnesium binding to either domain was shown to tune their affinity for calcium, bringing it into the physiological range. However, while the magnesium-binding site of CBD2 was identified, the identity of the CBD1 magnesium site remains elusive. Here, using molecular dynamics in combination with differential scanning fluorimetry and mutational analysis, we pinpoint the magnesium-binding site in CBD1. Specifically, among four calcium-binding sites (Ca1-Ca4) in this domain, only Ca1 can accommodate magnesium with an affinity similar to its free intracellular concentration. Moreover, our results provide mechanistic insights into the modulation of the regulatory calcium affinity by magnesium, which allows an adequate NCX activity level throughout varying physiological needs.

GAPS: a geometric attention-based network for peptide binding site identification by the transfer learning approach.

Protein-peptide interactions (PPepIs) are vital to understanding cellular functions, which can facilitate the design of novel drugs. As an essential component in forming a PPepI, protein-peptide binding sites are the basis for understanding the mechanisms involved in PPepIs. Therefore, accurately identifying protein-peptide binding sites becomes a critical task. The traditional experimental methods for researching these binding sites are labor-intensive and time-consuming, and some computational tools have been invented to supplement it. However, these computational tools have limitations in generality or accuracy due to the need for ligand information, complex feature construction, or their reliance on modeling based on amino acid residues. To deal with the drawbacks of these computational algorithms, we describe a geometric attention-based network for peptide binding site identification (GAPS) in this work. The proposed model utilizes geometric feature engineering to construct atom representations and incorporates multiple attention mechanisms to update relevant biological features. In addition, the transfer learning strategy is implemented for leveraging the protein-protein binding sites information to enhance the protein-peptide binding sites recognition capability, taking into account the common structure and biological bias between proteins and peptides. Consequently, GAPS demonstrates the state-of-the-art performance and excellent robustness in this task. Moreover, our model exhibits exceptional performance across several extended experiments including predicting the apo protein-peptide, protein-cyclic peptide and the AlphaFold-predicted protein-peptide binding sites. These results confirm that the GAPS model is a powerful, versatile, stable method suitable for diverse binding site predictions.

Identifying small-molecules binding sites in RNA conformational ensembles with SHAMAN.

The rational targeting of RNA with small molecules is hampered by our still limited understanding of RNA structural and dynamic properties. Most in silico tools for binding site identification rely on static structures and therefore cannot face the challenges posed by the dynamic nature of RNA molecules. Here, we present SHAMAN, a computational technique to identify potential small-molecule binding sites in RNA structural ensembles. SHAMAN enables exploring the conformational landscape of RNA with atomistic molecular dynamics simulations and at the same time identifying RNA pockets in an efficient way with the aid of probes and enhanced-sampling techniques. In our benchmark composed of large, structured riboswitches as well as small, flexible viral RNAs, SHAMAN successfully identifies all the experimentally resolved pockets and ranks them among the most favorite probe hotspots. Overall, SHAMAN sets a solid foundation for future drug design efforts targeting RNA with small molecules, effectively addressing the long-standing challenges in the field.

Metal-ligand dual-site single-atom nanozyme mimicking urate oxidase with high substrates specificity.

In nature, coenzyme-independent oxidases have evolved in selective catalysis using isolated substrate-binding pockets. Single-atom nanozymes (SAzymes), an emerging type of non-protein artificial enzymes, are promising to simulate enzyme active centers, but owing to the lack of recognition sites, realizing substrate specificity is a formidable task. Here we report a metal-ligand dual-site SAzyme (Ni-DAB) that exhibited selectivity in uric acid (UA) oxidation. Ni-DAB mimics the dual-site catalytic mechanism of urate oxidase, in which the Ni metal center and the C atom in the ligand serve as the specific UA and O2 binding sites, respectively, characterized by synchrotron soft X-ray absorption spectroscopy, in situ near ambient pressure X-ray photoelectron spectroscopy, and isotope labeling. The theoretical calculations reveal the high catalytic specificity is derived from not only the delicate interaction between UA and the Ni center but also the complementary oxygen reduction at the beta C site in the ligand. As a potential application, a Ni-DAB-based biofuel cell using human urine is constructed. This work unlocks an approach of enzyme-like isolated dual sites in boosting the selectivity of non-protein artificial enzymes.

Deviated binding of anti-HBV nucleoside analog E-CFCP-TP to the reverse transcriptase active site attenuates the effect of drug-resistant mutations.

While certain human hepatitis B virus-targeting nucleoside analogs (NAs) serve as crucial anti-HBV drugs, HBV yet remains to be a major global health threat. E-CFCP is a 4'-modified and fluoromethylenated NA that exhibits potent antiviral activity against both wild-type and drug-resistant HBVs but less potent against human immunodeficiency virus type-1 (HIV-1). Here, we show that HIV-1 with HBV-associated amino acid substitutions introduced into the RT's dNTP-binding site (N-site) is highly susceptible to E-CFCP. We determined the X-ray structures of HBV-associated HIV-1 RT mutants complexed with DNA:E-CFCP-triphosphate (E-CFCP-TP). The structures revealed that exocyclic fluoromethylene pushes the Met184 sidechain backward, and the resultant enlarged hydrophobic pocket accommodates both the fluoromethylene and 4'-cyano moiety of E-CFCP. Structural comparison with the DNA:dGTP/entecavir-triphosphate complex also indicated that the cyclopentene moiety of the bound E-CFCP-TP is slightly skewed and deviated. This positioning partly corresponds to that of the bound dNTP observed in the HIV-1 RT mutant with drug-resistant mutations F160M/M184V, resulting in the attenuation of the structural effects of F160M/M184V substitutions. These results expand our knowledge of the interactions between NAs and the RT N-site and should help further design antiviral NAs against both HIV-1 and HBV.

Binding sites and design strategies for small molecule GLP-1R agonists.

Glucagon-like peptide-1 receptor (GLP-1R) is a pivotal receptor involved in blood glucose regulation and influencing feeding behavior. It has received significant attention in the treatment of obesity and diabetes due to its potent incretin effect. Peptide GLP-1 receptor agonists (GLP-1RAs) have achieved tremendous success in the market, driving the vigorous development of small molecule GLP-1RAs. Currently, several small molecules have entered the clinical research stage. Additionally, recent discoveries of GLP-1R positive allosteric modulators (PAMs) are also unveiling new regulatory patterns and treatment methods. This article reviews the structure and functional mechanisms of GLP-1R, recent reports on small molecule GLP-1RAs and PAMs, as well as the optimization process. Furthermore, it combines computer simulations to analyze structure-activity relationships (SAR) studies, providing a foundation for exploring new strategies for designing small molecule GLP-1RAs.

Identification of heparin-binding amino acid residues in antibody HS4C3 with the potential to design antibodies against heparan sulfate domains.

Heparan sulfate (HS) is a linear polysaccharide with high structural and functional diversity. Detection and localization of HS in tissues can be performed using single chain variable fragment (scFv) antibodies. Although several anti-HS antibodies recognizing different sulfation motifs have been identified, little is known about their interaction with HS. In this study the interaction between the scFv antibody HS4C3 and heparin was investigated. Heparin-binding lysine and arginine residues were identified using a protect and label methodology. Site-directed mutagenesis was applied to further identify critical heparin-binding lysine/arginine residues using immunohistochemical and biochemical assays. In addition, computational docking of a heparin tetrasaccharide towards a 3-D homology model of HS4C3 was applied to identify potential heparin-binding sites. Of the 12 lysine and 15 arginine residues within the HS4C3 antibody, 6 and 9, respectively, were identified as heparin-binding. Most of these residues are located within one of the complementarity determining regions (CDR) or in their proximity. All basic amino acid residues in the CDR3 region of the heavy chain were involved in binding. Computational docking showed a heparin tetrasaccharide close to these regions. Mutagenesis of heparin-binding residues reduced or altered reactivity towards HS and heparin. Identification of heparin-binding arginine and lysine residues in HS4C3 allows for better understanding of the interaction with HS and creates a framework to rationally design antibodies targeting specific HS motifs.

Ligand-Receptor Interactions and Structure-Function Relationships in Off-Target Binding of the ß3-Adrenergic Agonist Mirabegron to α1A-Adrenergic Receptors.

The ß3-adrenoceptor agonist mirabegron is available for the treatment of storage symptoms of overactive bladder, including frequency, urgency, and incontinence. The off-target effects of mirabegron include binding to α1-adrenoceptors, which are central in the treatment of voiding symptoms. Here, we examined the structure-function relationships in the binding of mirabegron to a cryo-electron microscopy structure of α1A. The binding was simulated by docking mirabegron to a 3D structure of a human α1A-adrenoceptor (7YMH) using Autodock Vina. The simulations identified two binding states: slope orientation involving 10 positions and horizontal binding to the receptor surface involving 4 positions. No interactions occurred with positions constituting the α1A binding pocket, including Asp-106, Ser-188, or Phe-312, despite the positioning of the phenylethanolamine moiety in transmembrane regions close to the binding pocket by contact with Phe-288, -289, and Val-107. Contact with the unique positions of α1A included the transmembrane Met-292 during slope binding and exosite Phe-86 during horizontal binding. Exosite binding in slope orientation involved contact of the anilino part, rather than the aminothiazol end, to Ile-178, Ala-103, and Asn-179. In conclusion, contact with Met-292 and Phe-86, which are unique positions of α1A, accounts for mirabegron binding to α1A. Because of its lack of interactions with the binding pocket, mirabegron has lower affinity compared to α1A-blockers and no effects on voiding symptoms.

Conformational dynamics underlying atypical chemokine receptor 3 activation.

Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon ß-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of ß-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins.

Mapping the Escherichia coli DnaA-binding landscape reveals a preference for binding pairs of closely spaced DNA sites.

DnaA is a widely conserved DNA-binding protein that is essential for the initiation of DNA replication in many bacterial species, including Escherichia coli. Cooperative binding of ATP-bound DnaA to multiple 9mer sites ('DnaA boxes') at the origin of replication results in local unwinding of the DNA and recruitment of the replication machinery. DnaA also functions as a transcription regulator by binding to DNA sites upstream of target genes. Previous studies have identified many sites of direct positive and negative regulation by E. coli DnaA. Here, we use a ChIP-seq to map the E. coli DnaA-binding landscape. Our data reveal a compact regulon for DnaA that coordinates the initiation of DNA replication with expression of genes associated with nucleotide synthesis, replication, DNA repair and RNA metabolism. We also show that DnaA binds preferentially to pairs of DnaA boxes spaced 2 or 3 bp apart. Mutation of either the upstream or downstream site in a pair disrupts DnaA binding, as does altering the spacing between sites. We conclude that binding of DnaA at almost all target sites requires a dimer of DnaA, with each subunit making critical contacts with a DnaA box.

Summary of ChIP-Seq Methods and Description of an Optimized ChIP-Seq Protocol.

Chromatin immunoprecipitation (ChIP) is an invaluable method to characterize interactions between proteins and genomic DNA, such as the genomic localization of transcription factors and post-translational modification of histones. DNA and proteins are reversibly and covalently crosslinked using formaldehyde. Then the cells are lysed to release the chromatin. The chromatin is fragmented into smaller sizes either by micrococcal nuclease (MN) or sonication and then purified from other cellular components. The protein-DNA complexes are enriched by immunoprecipitation (IP) with antibodies that target the epitope of interest. The DNA is released from the proteins by heat and protease treatment, followed by degradation of contaminating RNAs with RNase. The resulting DNA is analyzed using various methods, including polymerase chain reaction (PCR), quantitative PCR (qPCR), or sequencing. This protocol outlines each of these steps for both yeast and human cells. This chapter includes a contextual discussion of the combination of ChIP with DNA analysis methods such as ChIP-on-Chip, ChIP-qPCR, and ChIP-Seq, recent updates on ChIP-Seq data analysis pipelines, complementary methods for identification of binding sites of DNA binding proteins, and additional protocol information about ChIP-qPCR and ChIP-Seq.