The power of computational proteomics platforms to decipher protein-protein interactions.

Adopting computational tools for analyzing extensive biological datasets has profoundly transformed our understanding and interpretation of biological phenomena. Innovative platforms have emerged, providing automated analysis to unravel essential insights about proteins and the complexities of their interactions. These computational advancements align with traditional studies, which employ experimental techniques to discern and quantify physical and functional protein-protein interactions (PPIs). Among these techniques, tandem mass spectrometry is notably recognized for its precision and sensitivity in identifying PPIs. These approaches might serve as important information enabling the identification of PPIs with potential pharmacological significance. This review aims to convey our experience using computational tools for detecting PPI networks and offer an analysis of platforms that facilitate predictions derived from experimental data.

Understanding the Differences of Danusertib's Residence Time in Aurora Kinases A/B: Dissociation Paths and Key Residues Identified using Conventional and Enhanced Molecular Dynamics Simulations.

The accurate experimental estimation of protein-ligand systems' residence time (τ) has become very relevant in drug design projects due to its importance in the last stages of refinement of the drug's pharmacodynamics and pharmacokinetics. It is now well-known that it is not sufficient to estimate the affinity of a protein-drug complex in the thermodynamic equilibrium process in in vitro experiments (closed systems), where the concentrations of the drug and protein remain constant. On the contrary, it is mandatory to consider the conformational dynamics of the system in terms of the binding and unbinding processes between protein and drugs in in vivo experiments (open systems), where their concentrations are in constant flux. This last model has been proven to dictate much of several drugs' pharmacological activities in vivo. At the atomistic level, molecular dynamics simulations can explain why some drugs are more effective than others or unveil the molecular aspects that make some drugs work better in one molecular target. Here, the protein kinases Aurora A/B, complexed with its inhibitor Danusertib, were studied using conventional and enhanced molecular dynamics (MD) simulations to estimate the dissociation paths and, therefore, the computational τ values and their comparison with experimental ones. Using classical molecular dynamics (cMD), three differential residues within the Aurora A/B active site, which seems to play an essential role in the observed experimental Danusertib's residence time against these kinases, were characterized. Then, using WT-MetaD, the relative Danusertib's residence times against Aurora A/B kinases were measured in a nanosecond time scale and were compared to those τ values observed experimentally. In addition, the potential dissociation paths of Danusertib in Aurora A and B were characterized, and differences that might be explained by the differential residues in the enzyme's active sites were found. In perspective, it is expected that this computational protocol can be applied to other protein-ligand complexes to understand, at the molecular level, the differences in residence times and amino acids that may contribute to it.

New NADPH Oxidase 2 Inhibitors Display Potent Activity against Oxidative Stress by Targeting p22phox-p47phox Interactions.

NADPH oxidase (NOX2) is responsible for reactive oxygen species (ROS) production in neutrophils and has been recognized as a key mediator in inflammatory and cardiovascular pathologies. Nevertheless, there is a lack of specific NOX2 pharmacological inhibitors. In medicinal chemistry, heterocyclic compounds are essential scaffolds for drug design, and among them, indole is a very versatile pharmacophore. We tested the hypothesis that indole heteroaryl-acrylonitrile derivatives may serve as NOX2 inhibitors by evaluating the capacity of 19 of these molecules to inhibit NOX2-derived ROS production in human neutrophils (HL-60 cells). Of these compounds, C6 and C14 exhibited concentration-dependent inhibition of NOX2 (IC50~1 µM). These molecules also reduced NOX2-derived oxidative stress in cardiomyocytes and prevented cardiac damage induced by ischemia-reperfusion. Compound C6 significantly reduced the membrane translocation of p47phox, a cytosolic subunit that is required for NOX2 activation. Molecular docking analyses of the binding modes of these molecules with p47phox indicated that C6 and C14 interact with specific residues in the inner part of the groove of p47phox, the binding cavity for p22phox. This combination of methods showed that novel indole heteroaryl acrylonitriles represent interesting lead compounds for developing specific and potent NOX2 inhibitors.

Modes of action of lysophospholipids as endogenous activators of the TRPV4 ion channel.

The Transient Receptor Potential Vanilloid 4 (TRPV4) channel has been shown to function in many physiological and pathophysiological processes. Despite abundant information on its importance in physiology, very few endogenous agonists for this channel have been described, and very few underlying mechanisms for its activation have been clarified. TRPV4 is expressed by several types of cells, such as vascular endothelial, and skin and lung epithelial cells, where it plays pivotal roles in their function. In the present study, we show that TRPV4 is activated by lysophosphatidic acid (LPA) in both endogenous and heterologous expression systems, pinpointing this molecule as one of the few known endogenous agonists for TRPV4. Importantly, LPA is a bioactive glycerophospholipid, relevant in several physiological conditions, including inflammation and vascular function, where TRPV4 has also been found to be essential. Here we also provide mechanistic details of the activation of TRPV4 by LPA and another glycerophospholipid, lysophosphatidylcholine (LPC), and show that LPA directly interacts with both the N- and C-terminal regions of TRPV4 to activate this channel. Moreover, we show that LPC activates TRPV4 by producing an open state with a different single-channel conductance to that observed with LPA. Our data suggest that the activation of TRPV4 can be finely tuned in response to different endogenous lipids, highlighting this phenomenon as a regulator of cell and organismal physiology. KEY POINTS: The Transient Receptor Potential Vaniloid (TRPV) 4 ion channel is a widely distributed protein with important roles in normal and disease physiology for which few endogenous ligands are known. TRPV4 is activated by a bioactive lipid, lysophosphatidic acid (LPA) 18:1, in a dose-dependent manner, in both a primary and a heterologous expression system. Activation of TRPV4 by LPA18:1 requires residues in the N- and C-termini of the ion channel. Single-channel recordings show that TRPV4 is activated with a decreased current amplitude (conductance) in the presence of lysophosphatidylcholine (LPC) 18:1, while LPA18:1 and GSK101 activate the channel with a larger single-channel amplitude. Distinct single-channel amplitudes produced by LPA18:1 and LPC18:1 could differentially modulate the responses of the cells expressing TRPV4 under different physiological conditions.

Unconventional interactions of the TRPV4 ion channel with beta-adrenergic receptor ligands.

The transient receptor potential vanilloid 4 (TRPV4) ion channel is present in different tissues including those of the airways. This channel is activated in response to stimuli such as changes in temperature, hypoosmotic conditions, mechanical stress, and chemicals from plants, lipids, and others. TRPV4's overactivity and/or dysfunction has been associated with several diseases, such as skeletal dysplasias, neuromuscular disorders, and lung pathologies such as asthma and cardiogenic lung edema and COVID-19-related respiratory malfunction. TRPV4 antagonists and blockers have been described; nonetheless, the mechanisms involved in achieving inhibition of the channel remain scarce, and the search for safe use of these molecules in humans continues. Here, we show that the widely used bronchodilator salbutamol and other ligands of ß-adrenergic receptors inhibit TRPV4's activation. We also demonstrate that inhibition of TRPV4 by salbutamol is achieved through interaction with two residues located in the outer region of the pore and that salbutamol leads to channel closing, consistent with an allosteric mechanism. Our study provides molecular insights into the mechanisms that regulate the activity of this physiopathologically important ion channel.

Structural Insights into the Inhibition Site in the Phosphorylcholine Phosphatase Enzyme of Pseudomonas aeruginosa.

Pseudomonas aeruginosa is a highly pathogenic Gram-negative microorganism associated with high mortality levels in burned or immunosuppressed patients or individuals affected by cystic fibrosis. Studies support a colonization mechanism whereby P. aeruginosa can breakdown the host cell membrane phospholipids through the sequential action of two enzymes: (I) hemolytic phospholipase C acting upon phosphatidylcholine or sphingomyelin to produce phosphorylcholine (Pcho) and (II) phosphorylcholine phosphatase (PchP) that hydrolyzes Pcho to generate choline and inorganic phosphate. This coordinated action provides the bacteria with carbon, nitrogen, and inorganic phosphate to support growth. Furthermore, PchP exhibits a distinctive inhibition mechanism by high substrate concentration. Here, we combine kinetic assays and computational approaches such as molecular docking, molecular dynamics, and free-energy calculations to describe the inhibitory site of PchP, which shares specific residues with the enzyme's active site. Our study provides insights into a coupled inhibition mechanism by the substrate, allowing us to postulate that the integrity of the inhibition site is needed to the correct functioning of the active site. Our results allow us to gain a better understanding of PchP function and provide the basis for a rational drug design that might contribute to the treatment of infections caused by this important opportunistic pathogen.

Structural Insights into the Substrate Transport Mechanisms in GTR Transporters through Ensemble Docking.

Glucosinolate transporters (GTRs) are part of the nitrate/peptide transporter (NPF) family, members of which also transport specialized secondary metabolites as substrates. Glucosinolates are defense compounds derived from amino acids. We selected 4-methylthiobutyl (4MTB) and indol-3-ylmethyl (I3M) glucosinolates to study how GTR1 from Arabidopsis thaliana transports these substrates in computational simulation approaches. The designed pipeline reported here includes massive docking of 4MTB and I3M in an ensemble of GTR1 conformations (in both inward and outward conformations) extracted from molecular dynamics simulations, followed by clustered and substrate-protein interactions profiling. The identified key residues were mutated, and their role in substrate transport was tested. We were able to identify key residues that integrate a major binding site of these substrates, which is critical for transport activity. In silico approaches employed here represent a breakthrough in the plant transportomics field, as the identification of key residues usually takes a long time if performed from a purely wet-lab experimental perspective. The inclusion of structural bioinformatics in the analyses of plant transporters significantly speeds up the knowledge-gaining process and optimizes valuable time and resources.

Bitopic Sigma 1 Receptor Modulators to Shed Light on Molecular Mechanisms Underpinning Ligand Binding and Receptor Oligomerization.

The sigma 1 receptor (S1R) is an enigmatic ligand-operated chaperone involved in many important biological processes, and its functions are not fully understood yet. Herein, we developed a novel series of bitopic S1R ligands as versatile tools to investigate binding processes, allosteric modulation, and the oligomerization mechanism. These molecules have been prepared in the enantiopure form and subjected to a preliminary biological evaluation, while in silico investigations helped to rationalize the results. Compound 7 emerged as the first bitopic S1R ligand endowed with low nanomolar affinity (Ki = 2.6 nM) reported thus far. Computational analyses suggested that 7 may stabilize the open conformation of the S1R by simultaneously binding the occluded primary binding site and a peripheral site on the cytosol-exposed surface. These findings pave the way to new S1R ligands with enhanced activity and/or selectivity, which could also be used as probes for the identification of a potential allosteric site.

PPI-MASS: An Interactive Web Server to Identify Protein-Protein Interactions From Mass Spectrometry-Based Proteomics Data.

Mass spectrometry-based proteomics methods are widely used to identify and quantify protein complexes involved in diverse biological processes. Specifically, tandem mass spectrometry methods represent an accurate and sensitive strategy for identifying protein-protein interactions. However, most of these approaches provide only lists of peptide fragments associated with a target protein, without performing further analyses to discriminate physical or functional protein-protein interactions. Here, we present the PPI-MASS web server, which provides an interactive analytics platform to identify protein-protein interactions with pharmacological potential by filtering a large protein set according to different biological features. Starting from a list of proteins detected by MS-based methods, PPI-MASS integrates an automatized pipeline to obtain information of each protein from freely accessible databases. The collected data include protein sequence, functional and structural properties, associated pathologies and drugs, as well as location and expression in human tissues. Based on this information, users can manipulate different filters in the web platform to identify candidate proteins to establish physical contacts with a target protein. Thus, our server offers a simple but powerful tool to detect novel protein-protein interactions, avoiding tedious and time-consuming data postprocessing. To test the web server, we employed the interactome of the TRPM4 and TMPRSS11a proteins as a use case. From these data, protein-protein interactions were identified, which have been validated through biochemical and bioinformatic studies. Accordingly, our web platform provides a comprehensive and complementary tool for identifying protein-protein complexes assisting the future design of associated therapies.

Beta-1,3 Oligoglucans Specifically Bind to Immune Receptor CD28 and May Enhance T Cell Activation.

Beta glucans are known to have immunomodulatory effects that mediated by a variety of mechanisms. In this article, we describe experiments and simulations suggesting that beta-1,3 glucans may promote activation of T cells by a previously unknown mechanism. First, we find that treatment of a T lymphoblast cell line with beta-1,3 oligoglucan significantly increases mRNA levels of T cell activation-associated cytokines, especially in the presence of the agonistic anti-CD3 antibody. This immunostimulatory activity was observed in the absence of dectin-1, a known receptor for beta-1,3 glucans. To clarify the molecular mechanism underlying this activity, we performed a series of molecular dynamics simulations and free-energy calculations to explore the interaction of beta-1,3 oligoglucans with potential immune receptors. While the simulations reveal little association between beta-1,3 oligoglucan and the immune receptor CD3, we find that beta-1,3 oligoglucans bind to CD28 near the region identified as the binding site for its natural ligands CD80 and CD86. Using a rigorous absolute binding free-energy technique, we calculate a dissociation constant in the low millimolar range for binding of 8-mer beta-1,3 oligoglucan to this site on CD28. The simulations show this binding to be specific, as no such association is computed for alpha-1,4 oligoglucan. This study suggests that beta-1,3 glucans bind to CD28 and may stimulate T cell activation collaboratively with T cell receptor activation, thereby stimulating immune function.

TMPRSS11a is a novel age-altered, tissue specific regulator of migration and wound healing.

Aging is a gradual biological process characterized by a decrease in cellular and organism functions. Aging-related processes involve changes in the expression and activity of several proteins. Here, we identified the transmembrane protease serine 11a (TMPRSS11a) as a new age-specific protein that plays an important role in skin wound healing. TMPRSS11a levels increased with age in rodent and human skin and gingival samples. Strikingly, overexpression of TMPRSS11a decreased cell migration and spreading, and inducing cellular senescence. Mass spectrometry, bioinformatics, and functional analyses revealed that TMPRSS11a interacts with integrin ß1 through an RGD sequence contained within the C-terminal domain and that this motif was relevant for cell migration. Moreover, TMPRSS11a was associated with cellular senescence, as shown by overexpression and downregulation experiments. In agreement with tissue-specific expression of TMPRSS11a, shRNA-mediated downregulation of this protein improved wound healing in the skin, but not in the skeletal muscle of old mice, where TMPRSS11a is undetectable. Collectively, these findings indicate that TMPRSS11a is a tissue-specific factor relevant for wound healing, which becomes elevated with aging, promoting cellular senescence and inhibiting cell migration and skin repair.

Structural Determinants of the Dopamine Transporter Regulation Mediated by G Proteins.

Dopamine clearance in the brain is controlled by the dopamine transporter (DAT), a protein residing in the plasma membrane, which drives reuptake of extracellular dopamine into presynaptic neurons. Studies have revealed that the ßγ subunits of heterotrimeric G proteins modulate DAT function through a physical association with the C-terminal region of the transporter. Regulation of neurotransmitter transporters by Gßγ subunits is unprecedented in the literature; therefore, it is interesting to investigate the structural details of this particular protein-protein interaction. Here, we refined the crystal structure of the Drosophila melanogaster DAT (dDAT), modeling de novo the N- and C-terminal domains; subsequently, we used the full-length dDAT structure to generate a comparative model of human DAT (hDAT). Both proteins were assembled with Gß1γ2 subunits employing protein-protein docking, and subsequent molecular dynamics simulations were run to identify the specific interactions governing the formation of the hDAT:Gßγ and dDAT:Gßγ complexes. A [L/F]R[Q/E]R sequence motif containing the residues R588 in hDAT and R587 in dDAT was found as key to bind the Gßγ subunits through electrostatic interactions with a cluster of negatively charged residues located at the top face of the Gß subunit. Alterations of DAT function have been associated with multiple devastating neuropathological conditions; therefore, this work represents a step toward better understanding DAT regulation by signaling proteins, allowing us to predict therapeutic target regions.

KCTD5, a novel TRPM4-regulatory protein required for cell migration as a new predictor for breast cancer prognosis.

Transient receptor potential melastatin 4 (TRPM4) is a Ca2+ -activated nonselective cationic channel that regulates cell migration and contractility. Increased TRPM4 expression has been related to pathologies, in which cytoskeletal rearrangement and cell migration are altered, such as metastatic cancer. Here, we identify the K+ channel tetramerization domain 5 (KCTD5) protein, a putative adaptor of cullin3 E3 ubiquitin ligase, as a novel TRPM4-interacting protein. We demonstrate that KCTD5 is a positive regulator of TRPM4 activity by enhancing its Ca2+ sensitivity. We show that through its effects on TRPM4 that KCTD5 promotes cell migration and contractility. Finally, we observed that both TRPM4 and KCTD5 expression are increased in distinct patterns in different classes of breast cancer tumor samples. Together, these data support that TRPM4 activity can be regulated through expression levels of either TRPM4 or KCTD5, not only contributing to increased understanding of the molecular mechanisms involved on the regulation of these important ion channels, but also providing information that could inform treatments based on targeting these distinct molecules that define TRPM4 activity.

A transient amphipathic helix in the prodomain of PCSK9 facilitates binding to low-density lipoprotein particles.

Proprotein convertase subtilisin/kexin type-9 (PCSK9) is a ligand of low-density lipoprotein (LDL) receptor (LDLR) that promotes LDLR degradation in late endosomes/lysosomes. In human plasma, 30-40% of PCSK9 is bound to LDL particles; however, the physiological significance of this interaction remains unknown. LDL binding in vitro requires a disordered N-terminal region in PCSK9's prodomain. Here, we report that peptides corresponding to a predicted amphipathic α-helix in the prodomain N terminus adopt helical structure in a membrane-mimetic environment. This effect was greatly enhanced by an R46L substitution representing an atheroprotective PCSK9 loss-of-function mutation. A helix-disrupting proline substitution within the putative α-helical motif in full-length PCSK9 lowered LDL binding affinity >5-fold. Modeling studies suggested that the transient α-helix aligns multiple polar residues to interact with positively charged residues in the C-terminal domain. Gain-of-function PCSK9 mutations associated with familial hypercholesterolemia (FH) and clustered at the predicted interdomain interface (R469W, R496W, and F515L) inhibited LDL binding, which was completely abolished in the case of the R496W variant. These findings shed light on allosteric conformational changes in PCSK9 required for high-affinity binding to LDL particles. Moreover, the initial identification of FH-associated mutations that diminish PCSK9's ability to bind LDL reported here supports the notion that PCSK9-LDL association in the circulation inhibits PCSK9 activity.

New Insights into the Opening of the Occluded Ligand-Binding Pocket of Sigma1 Receptor: Binding of a Novel Bivalent RC-33 Derivative.

Significant progresses have been made to understand the molecular basis of the Sigma1 receptor (S1R) operating in normal and pathological conditions. S1R is a transmembrane protein that participates in a wide variety of processes at the central nervous system; hence, its function has been associated with mental and neurological disorders. Several ligands have been proposed to regulate the function of S1R revealing a high plasticity of the ligand-binding pocket. Previous drug-design studies have been mainly based on pharmacophore models; however, the recently revealed crystal structure of S1R provides an excellent opportunity for verifying previous predictions and for evaluating the binding of novel compounds. Interestingly, the crystal structure shows that the binding pocket of S1R is highly occluded from solvent; therefore, it is not clear how ligands access this site. In the present work, we applied steered molecular dynamics (SMD) simulations to open the occluded ligand-binding pocket in the S1R crystal structure and to determine the preferred ligand pathway to enter and exit the binding site. The intracellular surface of the ß-barrel ligand-binding region was found the most favorable route to accommodate ligands. This route supports the binding of RC-33 (our in-house-developed S1R modulator) and a new bivalent derivative that constitutes the first divalent structure shown to interact with S1R. Free energy calculations of these compounds associated with S1R agree with experimental Ki values and provide molecular insights of the binding mode of modulators that could access the S1R ligand-binding pocket through the cytoplasmic region.

A Structural Model of the Inactivation Gate of Voltage-Activated Potassium Channels.

After opening, the Shaker voltage-gated potassium (KV) channel rapidly inactivates when one of its four N-termini enters and occludes the channel pore. Although it is known that the tip of the N-terminus reaches deep into the central cavity, the conformation adopted by this domain during inactivation and the nature of its interactions with the rest of the channel remain unclear. Here, we use molecular dynamics simulations coupled with electrophysiology experiments to reveal the atomic-scale mechanisms of inactivation. We find that the first six amino acids of the N-terminus spontaneously enter the central cavity in an extended conformation, establishing hydrophobic contacts with residues lining the pore. A second portion of the N-terminus, consisting of a long 24 amino acid α-helix, forms numerous polar contacts with residues in the intracellular entryway of the T1 domain. Double mutant cycle analysis revealed a strong relationship between predicted interatomic distances and empirically observed thermodynamic coupling, establishing a plausible model of the transition of KV channels to the inactivated state.

An extracellular cation coordination site influences ion conduction of OsHKT2;2.

BACKGROUND: HKT channels mediate sodium uniport or sodium and potassium symport in plants. Monocotyledons express a higher number of HKT proteins than dicotyledons, and it is only within this clade of HKT channels that cation symport mechanisms are found. The prevailing ion composition in the extracellular medium affects the transport abilities of various HKT channels by changing their selectivity or ion transport rates. How this mutual effect is achieved at the molecular level is still unknown. Here, we built a homology model of the monocotyledonous OsHKT2;2, which shows sodium and potassium symport activity. We performed molecular dynamics simulations in the presence of sodium and potassium ions to investigate the mutual effect of cation species. RESULTS: By analyzing ion-protein interactions, we identified a cation coordination site on the extracellular protein surface, which is formed by residues P71, D75, D501 and K504. Proline and the two aspartate residues coordinate cations, while K504 forms salt bridges with D75 and D501 and may be involved in the forwarding of cations towards the pore entrance. Functional validation via electrophysiological experiments confirmed the biological relevance of the predicted ion coordination site and identified K504 as a central key residue. Mutation of the cation coordinating residues affected the functionality of HKT only slightly. Additional in silico mutants and simulations of K504 supported experimental results. CONCLUSION: We identified an extracellular cation coordination site, which is involved in ion coordination and influences the conduction of OsHKT2;2. This finding proposes a new viewpoint in the discussion of how the mutual effect of variable ion species may be achieved in HKT channels.

EB1- and EB2-dependent anterograde trafficking of TRPM4 regulates focal adhesion turnover and cell invasion.

Transient receptor potential melastatin 4 (TRPM4) is a Ca2+-activated nonselective cationic channel involved in a wide variety of physiologic and pathophysiological processes. Bioinformatics analyses of the primary sequence of TRPM4 allowed us to identify a putative motif for interaction with end-binding (EB) proteins, which are microtubule plus-end tracking proteins. Here, we provide novel data suggesting that TRPM4 interacts with EB proteins. We show that mutations of the putative EB binding motif abolish the TRPM4-EB interaction, leading to a reduced expression of the mature population of the plasma membrane channel and instead display an endoplasmic reticulum-associated distribution. Furthermore, we demonstrate that EB1 and EB2 proteins are required for TRPM4 trafficking and functional activity. Finally, we demonstrated that the expression of a soluble fragment containing the EB binding motif of TRPM4 reduces the plasma membrane expression of the channel and affects TRPM4-dependent focal adhesion disassembly and cell invasion processes.-Blanco, C., Morales, D., Mogollones, I., Vergara-Jaque, A., Vargas, C., Álvarez, A., Riquelme, D., Leiva-Salcedo, E., González, W., Morales, D., Maureira, D., Aldunate, I., Cáceres, M., Varela, D., Cerda, O. EB1- and EB2-dependent anterograde trafficking of TRPM4 regulates focal adhesion turnover and cell invasion.

Polymeric substances for the removal of ochratoxin A from red wine followed by computational modeling of the complexes formed.

Ochratoxin A (OTA) is a mycotoxin produced by filamentous-type fungi that contaminates a wide variety of foods and beverages such as wines. In these trials, we evaluated the capacity of the following polymers for the removal of OTA from acidic model solutions and red wine: polyvinylpolypyrrolidone (PVPP), resin of N-vinyl-2-pyrrolidinone with ethylene glycol dimethacrylate and triallyl isocyanurate (PVP-DEGMA-TAIC), and poly(acrylamide-co-ethylene glycol-dimethacrylate) (PA-EGDMA). In acidic model solution, PVP-DEGMA-TAIC and PA-EGDMA polymers removed up to 99.9% of OTA, but their trapping capacity was highly reduced by the presence of competing phenolic substances (i.e. gallic acid and 4-methylcathecol). In real red wine, PA-EGDMA polymer showed the most promising results, with more than 68.0% OTA removal and less than 14.0% reduction in total phenolic. Finally, computational chemistry analyses showed that the affinity between OTA and the polymers studied would be due to Van der Waals interactions.