RINspector: a Cytoscape app for centrality analyses and DynaMine flexibility prediction

What is the potential impact of genetic divergence of plastid ribosomal genes between Silene nutans lineages in hybrids? An in silico approach using the 3D structure of the plastid ribosome

Introduction

Following the integration of cyanobacteria into the eukaryotic cells, many genes were transferred from the plastid to the nucleus. As a result, plastid complexes are encoded both by plastid and nuclear genes. Tight co-adaptation is required between these genes as plastid and nuclear genomes differ in several characteristics, such as mutation rate and inheritance patterns. Among these are complexes from the plastid ribosome, composed of two main subunits: a large and a small one, both composed of nuclear and plastid gene products. This complex has been identified as a potential candidate for sheltering plastid-nuclear incompatibilities in a Caryophyllaceae species, Silene nutans. This species is composed of four genetically differentiated lineages, which exhibit hybrid breakdown when interlineage crosses are conducted. As this complex is composed of numerous interacting plastid-nuclear gene pairs, in the present study, the goal was to reduce the number of gene pairs that could induce such incompatibilities.

Method

We used the previously published 3D structure of the spinach ribosome to further elucidate which of the potential gene pairs might disrupt plastid-nuclear interactions within this complex. After modeling the impact of the identified mutations on the 3D structure, we further focused on one strongly mutated plastid-nuclear gene pair: rps11-rps21. We used the centrality measure of the mutated residues to further understand if the modified interactions and associated modified centralities might be correlated with hybrid breakdown.

Results and discussion

This study highlights that lineage-specific mutations in essential plastid and nuclear genes might disrupt plastid-nuclear protein interactions of the plastid ribosome and that reproductive isolation correlates with changes in residue centrality values. Because of this, the plastid ribosome might be involved in hybrid breakdown in this system.

Direct interaction of TrkA/CD44v3 is essential for NGF-promoted aggressiveness of breast cancer cells

Background

CD44 is a multifunctional membrane glycoprotein. Through its heparan sulfate chain, CD44 presents growth factors to their receptors. We have shown that CD44 and Tropomyosin kinase A (TrkA) form a complex following nerve growth factor (NGF) induction. Our study aimed to understand how CD44 and TrkA interact and the consequences of inhibiting this interaction regarding the pro-tumoral effect of NGF in breast cancer.

Methods

After determining which CD44 isoforms (variants) are involved in forming the TrkA/CD44 complex using proximity ligation assays, we investigated the molecular determinants of this interaction. By molecular modeling, we isolated the amino acids involved and confirmed their involvement using mutations. A CD44v3 mimetic peptide was then synthesized to block the TrkA/CD44v3 interaction. The effects of this peptide on the growth, migration and invasion of xenografted triple-negative breast cancer cells were assessed. Finally, we investigated the correlations between the expression of the TrkA/CD44v3 complex in tumors and histo-pronostic parameters.

Results

We demonstrated that isoform v3 (CD44v3), but not v6, binds to TrkA in response to NGF stimulation. The final 10 amino acids of exon v3 and the TrkA H112 residue are necessary for the association of CD44v3 with TrkA. Functionally, the CD44v3 mimetic peptide impairs not only NGF-induced RhoA activation, clonogenicity, and migration/invasion of breast cancer cells in vitro but also tumor growth and metastasis in a xenograft mouse model. We also detected TrkA/CD44v3 only in cancerous cells, not in normal adjacent tissues.

Conclusion

Collectively, our results suggest that blocking the CD44v3/TrkA interaction can be a new therapeutic option for triple-negative breast cancers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02314-4.

Targeting SARS-CoV-2 nsp13 Helicase and Assessment of Druggability Pockets: Identification of Two Potent Inhibitors by a Multi-Site In Silico Drug Repurposing Approach

The SARS-CoV-2 non-structural protein 13 (nsp13) helicase is an essential enzyme for viral replication and has been identified as an attractive target for the development of new antiviral drugs. In detail, the helicase catalyzes the unwinding of double-stranded DNA or RNA in a 5' to 3' direction and acts in concert with the replication-transcription complex (nsp7/nsp8/nsp12). In this work, bioinformatics and computational tools allowed us to perform a detailed conservation analysis of the SARS-CoV-2 helicase genome and to further predict the druggable enzyme's binding pockets. Thus, a structure-based virtual screening was used to identify valuable compounds that are capable of recognizing multiple nsp13 pockets. Starting from a database of around 4000 drugs already approved by the Food and Drug Administration (FDA), we chose 14 shared compounds capable of recognizing three out of four sites. Finally, by means of visual inspection analysis and based on their commercial availability, five promising compounds were submitted to in vitro assays. Among them, PF-03715455 was able to block both the unwinding and NTPase activities of nsp13 in a micromolar range.

Design, construction and in vivo functional assessment of a hinge truncated sFLT01

Gene therapy for the treatment of ocular neovascularization has reached clinical trial phases. The AAV2-sFLT01 construct was already evaluated in a phase 1 open-label trial administered intravitreally to patients with advanced neovascular age-related macular degeneration. SFLT01 protein functions by binding to VEGF and PlGF molecules and inhibiting their activities simultaneously. It consists of human VEGFR1/Flt-1 (hVEGFR1), a polyglycine linker, and the Fc region of human IgG1. The IgG1 upper hinge region of the sFLT01 molecule makes it vulnerable to radical attacks and prone to causing immune reactions. This study pursued two goals: (i) minimizing the immunogenicity and vulnerability of the molecule by designing a truncated molecule called htsFLT01 (hinge truncated sFLT01) that lacked the IgG1 upper hinge and lacked 2 amino acids from the core hinge region; and (ii) investigating the structural and functional properties of the aforesaid chimeric molecule at different levels (in silico, in vitro, and in vivo). Molecular dynamics simulations and molecular mechanics energies combined with Poisson-Boltzmann and surface area continuum solvation calculations revealed comparable free energy of binding and binding affinity for sFLT01 and htsFLT01 to their cognate ligands. Conditioned media from human retinal pigment epithelial (hRPE) cells that expressed htsFLT01 significantly reduced tube formation in HUVECs. The AAV2-htsFLT01 virus suppressed vascular development in the eyes of newborn mice. The htsFLT01 gene construct is a novel anti-angiogenic tool with promising improvements compared to existing treatments.

Binding Networks Identify Targetable Protein Pockets for Mechanism-Based Drug Design

The human genome codes only a few thousand druggable proteins, mainly receptors and enzymes. While this pool of available drug targets is limited, there is an untapped potential for discovering new drug-binding mechanisms and modes. For example, enzymes with long binding cavities offer numerous prerequisite binding sites that may be visited by an inhibitor during migration from a bulk solution to the destination site. Drug design can use these prerequisite sites as new structural targets. However, identifying these ephemeral sites is challenging. Here, we introduce a new method called NetBinder for the systematic identification and classification of prerequisite binding sites at atomic resolution. NetBinder is based on atomistic simulations of the full inhibitor binding process and provides a networking framework on which to select the most important binding modes and uncover the entire binding mechanism, including previously undiscovered events. NetBinder was validated by a study of the binding mechanism of blebbistatin (a potent inhibitor) to myosin 2 (a promising target for cancer chemotherapy). Myosin 2 is a good test enzyme because, like other potential targets, it has a long internal binding cavity that provides blebbistatin with numerous potential prerequisite binding sites. The mechanism proposed by NetBinder of myosin 2 structural changes during blebbistatin binding shows excellent agreement with experimentally determined binding sites and structural changes. While NetBinder was tested on myosin 2, it may easily be adopted to other proteins with long internal cavities, such as G-protein-coupled receptors or ion channels, the most popular current drug targets. NetBinder provides a new paradigm for drug design by a network-based elucidation of binding mechanisms at an atomic resolution.

TRPM8-Rap1A Interaction Sites as Critical Determinants for Adhesion and Migration of Prostate and Other Epithelial Cancer Cells

Emerging evidence indicates that the TRPM8 channel plays an important role in prostate cancer (PCa) progression, by impairing the motility of these cancer cells. Here, we reveal a novel facet of PCa motility control via direct protein-protein interaction (PPI) of the channel with the small GTPase Rap1A. The functional interaction of the two proteins was assessed by active Rap1 pull-down assays and live-cell imaging experiments. Molecular modeling analysis allowed the identification of four putative residues involved in TRPM8-Rap1A interaction. Point mutations of these sites impaired PPI as shown by GST-pull-down, co-immunoprecipitation, and PLA experiments and revealed their key functional role in the adhesion and migration of PC3 prostate cancer cells. More precisely, TRPM8 inhibits cell migration and adhesion by trapping Rap1A in its GDP-bound inactive form, thus preventing its activation at the plasma membrane. In particular, residues E207 and Y240 in the sequence of TRPM8 and Y32 in that of Rap1A are critical for the interaction between the two proteins not only in PC3 cells but also in cervical (HeLa) and breast (MCF-7) cancer cells. This study deepens our knowledge of the mechanism through which TRPM8 would exert a protective role in cancer progression and provides new insights into the possible use of TRPM8 as a new therapeutic target in cancer treatment.

Landscape-Based Protein Stability Analysis and Network Modeling of Multiple Conformational States of the SARS-CoV-2 Spike D614G Mutant: Conformational Plasticity and Frustration-Induced Allostery as Energetic Drivers of Highly Transmissible Spike Variants

The structural and functional studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across many variants of concern (VOCs), suggesting the effect of this mutation on the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies provided important evidence about multiple conformational substates of the D614G spike protein. The development of a plausible mechanistic model that can explain the experimental observations from a more unified thermodynamic perspective is an important objective of the current work. In this study, we employed efficient and accurate coarse-grained simulations of multiple structural substates of the D614G spike trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration profiling of the conformational states with residue-based mutational scanning of protein stability and network analysis of allosteric interactions and communications, we determine the patterns of mutational sensitivity in the functional regions and sites of variants. We found that the D614G mutation may induce a considerable conformational adaptability of the open states in the SARS-CoV-2 spike protein without compromising the folding stability and integrity of the spike protein. The results suggest that the D614G mutant may employ a hinge-shift mechanism in which the dynamic couplings between the site of mutation and the interprotomer hinge modulate the interdomain interactions, global mobility change, and the increased stability of the open form. This study proposes that mutation-induced modulation of the conformational flexibility and energetic frustration at the interprotomer interfaces may serve as an efficient mechanism for allosteric regulation of the SARS-CoV-2 spike proteins.

Electrophysiological, Morphologic, and Transcriptomic Profiling of the Ogura-CMS, DGMS and Maintainer Broccoli Lines

To better serve breeding of broccoli, the electrophysiological, morphological and transcriptomic profiling of the isogenic Ogura-CMS, DGMS and their maintainer fertile lines, were carried out by scanning electron microscopy, investigation of agronomic traits and RNA-sequencing analysis. The agronomic traits of plant height, length of the largest leaf, plant spread angle, single head weight, head width and stem diameter showed stronger performance in Ogura-CMS broccoli than in DGMS line or maintainer fertile line. However, the Ogura-CMS broccoli was poorer in the seed yield and seed germination than in the DGMS line and maintainer fertile line. Additionally, the DGMS broccoli had longer maturation and flowering periods than the Ogura-CMS and maintainer fertile lines. There were obvious differences in the honey gland, happening in the male sterility and fertile lines of broccoli. Additionally, the mechanism regulating Ogura-CMS and DGMS in broccoli was investigated using florets transcriptome analyses of the Ogura-CMS, DGMS and maintainer fertile lines. As a result, a total of 2670 differentially expressed genes (DEGs) were detected, including 1054 up- and 1616 downregulated genes in the Ogura-CMS and DGMS lines compared to the maintainer fertile line. A number of functionally known genes involved in plant hormones (auxin, salicylic acid and brassinosteroid), five Mitochondrial Oxidative Phosphorylation (OXPHOS) genes of atp8, LOC106319879, LOC106324734, LOC106314622 and LOC106298585, and three upregulated genes (Lhcb1, Lhcb3 and Lhcb5) associated with the photosynthesis-antenna protein pathway, were obviously detected to be highly associated with reproductive development including flowering time, maturity and reproductive period in the Ogura-CMS and DGMS broccoli comparing to their maintainer fertile line. Our research would provide a comprehensive foundation for understanding the differences of electrophysiological, morphological and transcriptomic profiles in the Ogura-CMS, DGMS and maintainer broccoli, and as well as being beneficial to exploring the mechanism of male sterility in Brassica crops.

Theta-Defensins to Counter COVID-19 as Furin Inhibitors: In Silico Efficiency Prediction and Novel Compound Design

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was characterized as a pandemic by the World Health Organization (WHO) in Dec. 2019. SARS-CoV-2 binds to the cell membrane through spike proteins on its surface and infects the cell. Furin, a host-cell enzyme, possesses a binding site for the spike protein. Thus, molecules that block furin could potentially be a therapeutic solution. Defensins are antimicrobial peptides that can hypothetically inhibit furin because of their arginine-rich structure. Theta-defensins, a subclass of defensins, have attracted attention as drug candidates due to their small size, unique structure, and involvement in several defense mechanisms. Theta-defensins could be a potential treatment for COVID-19 through furin inhibition and an anti-inflammatory mechanism. Note that inflammatory events are a significant and deadly condition that could happen at the later stages of COVID-19 infection. Here, the potential of theta-defensins against SARS-CoV-2 infection was investigated through in silico approaches. Based on docking analysis results, theta-defensins can function as furin inhibitors. Additionally, a novel candidate peptide against COVID-19 with optimal properties regarding antigenicity, stability, electrostatic potential, and binding strength was proposed. Further in vitro/in vivo investigations could verify the efficiency of the designed novel peptide.

Hotspots for mutations in the SARS-CoV-2 spike glycoprotein: a correspondence analysis

Spike glycoprotein (Sgp) is liable for binding of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to the host receptors. Since Sgp is the main target for vaccine and drug designing, elucidating its mutation pattern could help in this regard. This study is aimed at investigating the correspondence of specific residues to the SgpSARS-CoV-2 functionality by explorative interpretation of sequence alignments. Centrality analysis of the Sgp dissects the importance of these residues in the interaction network of the RBD-ACE2 (receptor-binding domain) complex and furin cleavage site. Correspondence of RBD to threonine500 and asparagine501 and furin cleavage site to glutamine675, glutamine677, threonine678, and alanine684 was observed; all residues are exactly located at the interaction interfaces. The harmonious location of residues dictates the RBD binding property and the flexibility, hydrophobicity, and accessibility of the furin cleavage site. These species-specific residues can be assumed as real targets of evolution, while other substitutions tend to support them. Moreover, all these residues are parts of experimentally identified epitopes. Therefore, their substitution may affect vaccine efficacy. Higher rate of RBD maintenance than furin cleavage site was predicted. The accumulation of substitutions reinforces the probability of the multi-host circulation of the virus and emphasizes the enduring evolutionary events.

Residue interaction networks of K-Ras protein with water molecules identifies the potential role of switch II and P-loop

The mutant K-Ras with aberrant signaling is the primary cause of several cancers. The proposed study investigated the influence of water molecules in K-Ras crystal structure, where they have a significant function by understanding their residue interaction networks (RINs). We analyzed the RINs of K-Ras with and without water molecules and determined their interaction properties. RINs were developed with the help of StructureViz2 and RINspector; further, the changes in K-Ras backbone flexibility were predicted with the DynaMine. We found that the residues K42, I142, and L159 are the hotspots from water, including the K-Ras-GTP complex with the highest residue centrality analysis (RCA) Z-score. The DynaMine prediction calculated the NMR S² value for the frequently mutated positions G12, G13, and Q61 showing a minor shift in flexibility, which make up the P-Loop and switch II of the K-Ras protein. This flexibility shift can account for changes in conformational activity and the protein's GTPase activity, making it difficult to recognize by the effectors and exchange factors. Taken together, our study helps in understanding the functional importance of the water molecules in K-Ras protein and the impact of mutation that modulate the conformational state of the protein.

Centrality Measures in Residue Interaction Networks to Highlight Amino Acids in Protein–Protein Binding

Residue interaction networks (RINs) describe a protein structure as a network of interacting residues. Central nodes in these networks, identified by centrality analyses, highlight those residues that play a role in the structure and function of the protein. However, little is known about the capability of such analyses to identify residues involved in the formation of macromolecular complexes. Here, we performed six different centrality measures on the RINs generated from the complexes of the SKEMPI 2 database of changes in protein–protein binding upon mutation in order to evaluate the capability of each of these measures to identify major binding residues. The analyses were performed with and without the crystallographic water molecules, in addition to the protein residues. We also investigated the use of a weight factor based on the inter-residue distances to improve the detection of these residues. We show that for the identification of major binding residues, closeness, degree, and PageRank result in good precision, whereas betweenness, eigenvector, and residue centrality analyses give a higher sensitivity. Including water in the analysis improves the sensitivity of all measures without losing precision. Applying weights only slightly raises the sensitivity of eigenvector centrality analysis. We finally show that a combination of multiple centrality analyses is the optimal approach to identify residues that play a role in protein–protein interaction.

Design of an engineered ACE2 as a novel therapeutics against COVID-19

The interaction between the angiotensin-converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of the spike protein from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plays a pivotal role in virus entry into the host cells. Since recombinant ACE2 protein has been suggested as an anti-SARS-CoV-2 therapeutic agent, this study was conducted to design an ACE2 protein with more desirable properties. In this regard, the amino acids with central roles in enzymatic activity of the ACE2 were substituted. Moreover, saturation mutagenesis at the interaction interface between the ACE2 and RBD was performed to increase their interaction affinity. The best mutations to increase the structural and thermal stability of the ACE2 were also selected based on B factors and mutation effects. The obtained resulted revealed that the Arg273Gln and Thr445Gly mutation have drastically reduced the binding affinity of the angiotensin-II into the active site of ACE2. The Thr27Arg mutation was determined to be the most potent mutation to increase the binding affinity. The Asp427Arg mutation was done to decrease the flexibility of the region with high B factor. The Pro451Met mutation along with the Gly448Trp mutation was predicted to increase the thermodynamic stability and thermostability of the ACE2. The designed therapeutic ACE2 would have no enzymatic activity while it could bear stronger interaction with Spike glycoprotein of the SARS-CoV-2. Moreover, decreased in vivo enzymatic degradation would be anticipated due to increased thermostability. This engineered ACE2 could be exploited as a novel therapeutic agent against COVID-19 after necessary evaluations.

Identifying Synergistic Mechanisms of Multiple Ingredients in Shuangbai Tablets against Proteinuria by Virtual Screening and a Network Pharmacology Approach

Shuangbai Tablets (SBT), a traditional herbal mixture, has shown substantial clinical efficacy. However, a systematic mechanism of its active ingredients and pharmacological mechanisms of action against proteinuria continues being lacking. A network pharmacology approach was effectual in discovering the relationship of multiple ingredients and targets of the herbal mixture. This study aimed to identify key targets, major active ingredients, and pathways of SBT against proteinuria by network pharmacology approach combined with thin layer chromatography (TLC). Human phenotype (HP) disease analysis, gene ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and molecular docking were used in this study. To this end, a total of 48 candidate targets of 118 active ingredients of SBT were identified. Network analysis showed PTGS2, ESR1, and NOS2 to be the three key targets, and beta-sitosterol, quercetin, and berberine were the three major active ingredients; among them one of the major active ingredients, quercetin, was discriminated by TLC. These results of the functional enrichment analysis indicated that the most relevant disease including these 48 candidate proteins is proteinuria, SBT treated proteinuria by sympathetically regulating multiple biological pathways, such as the HIF-1, RAS, AGE-RAGE, and VEGF signaling pathways. Additionally, molecular docking validation suggested that major active ingredients of SBT were capable of binding to HIF-1A and VEGFA of the main pathways. Consequently, key targets, major active ingredients, and pathways based on data analysis of SBT against proteinuria were systematically identified confirming its utility and providing a new drug against proteinuria.

Identification of Novel Interaction Partners of Ets-1: Focus on DNA Repair

The transcription factor Ets-1 (ETS proto-oncogene 1) shows low expression levels except in specific biological processes like haematopoiesis or angiogenesis. Elevated levels of expression are observed in tumor progression, resulting in Ets-1 being named an oncoprotein. It has recently been shown that Ets-1 interacts with two DNA repair enzymes, PARP-1 (poly(ADP-ribose) polymerase 1) and DNA-PK (DNA-dependent protein kinase), through two different domains and that these interactions play a role in cancer. Considering that Ets-1 can bind to distinctly different domains of two DNA repair enzymes, we hypothesized that the interaction can be transposed onto homologs of the respective domains. We have searched for sequence and structure homologs of the interacting ETS(Ets-1), BRCT(PARP-1) and SAP(DNA-PK) domains, and have identified several candidate binding pairs that are currently not annotated as such. Many of the Ets-1 partners are associated to DNA repair mechanisms. We have applied protein-protein docking to establish putative interaction poses and investigated these using centrality analyses at the protein residue level. Most of the identified poses are virtually similar to our recently established interaction model for Ets-1/PARP-1 and Ets-1/DNA-PK. Our work illustrates the potentially high number of interactors of Ets-1, in particular involved in DNA repair mechanisms, which shows the oncoprotein as a potential important regulator of the mechanism.

The Inclusion of Water Molecules in Residue Interaction Networks Identifies Additional Central Residues

The relevance of water molecules for the recognition and the interaction of biomolecules is widely appreciated. In this paper we address the role that water molecules associated to protein complexes play for the functional relevance of residues by considering their residue interaction networks (RINs). These are commonly defined on the basis of the amino acid composition of the proteins themselves, disregarding the solvation state of the protein. We determine properties of the RINs of two protein complexes, colicin E2/Im2 and barnase/barstar, with and without associated water molecules, using a previously developed methodology and its associated application RINspector. We find that the inclusion of water molecules in RINs leads to an increase in the number of central residues which adds a novel mechanism to the relevance of water molecules for protein function.

Computational characterization of the binding mode between oncoprotein Ets-1 and DNA-repair enzymes

The Ets-1 oncoprotein is a transcription factor that promotes target gene expression in specific biological processes. Typically, Ets-1 activity is low in healthy cells, but elevated levels of expression have been found in cancerous cells, specifically related to tumor progression. Like the vast majority of the cellular effectors, Ets-1 does not act alone but in association with partners. Given the important role that is attributed to Ets-1 in major human diseases, it is crucial to identify its partners and characterize their interactions. In this context, two DNA-repair enzymes, PARP-1 and DNA-PK, have been identified recently as interaction partners of Ets-1. We here identify their binding mode by means of protein docking. The results identify the interacting surface between Ets-1 and the two DNA-repair enzymes centered on the α-helix H1 of the ETS domain, leaving α-helix H3 available to bind DNA. The models highlight a hydrophobic patch on Ets-1 at the center of the interaction interface that includes three tryptophans (Trp338, Trp356, and Trp361). We rationalize the binding mode using a series of computational analyses, including alanine scanning, molecular dynamics simulation, and residue centrality analysis. Our study constitutes a first but important step in the characterization, at the molecular level, of the interaction between an oncoprotein and DNA-repair enzymes.

Dynamic perspectives into the mechanisms of mutation-induced p53-DNA binding loss and inactivation using active perturbation theory: Structural and molecular insights toward the design of potent reactivators in cancer therapy

The DNA-binding ability of p53 represents the crux of its tumor suppressive activities, which involves transcriptional activation of target genes responsible for apoptosis and cell-cycle arrest. Mutational occurrences within or in close proximity to the DNA-binding surface of p53 have accounted for the loss of direct DNA-binding ability and inactivation implicated in many cases of cancer. Moreover, the design of therapeutic compounds that can restore DNA-binding ability in p53 mutants has been identified as a way forward in curtailing their oncogenic activities. However, there is still the need for more insights into evaluate the perturbations that occur at the DNA-binding interface of mp53 relative to DNA-binding loss, inactivation, and design of potent reactivators, hence the purpose of this study. Therefore, we evaluated p53-structural (R175H) and contact (R273C) mutational effects using tunnel perturbation analysis and other computational tools. We identified significant perturbations in the active tunnels of p53, which resulted in altered geometry and loss, unlike in the wild-type p53. This corroborated with structural, DNA-binding, and interaction network analysis, which showed that loss of flexibility, repulsion of DNA-interactive residues, and instability occurred at the binding interface of both mutants. Also, these mutations altered bonding interactions and network topology at the DNA-binding interface, resulting in the reduction of p53-DNA binding proximity and affinity. Therefore, these findings would aid the structure-based design of novel chemical entities capable of restoring p53-DNA binding and activation.

Structure and dynamics of a human myelin protein P2 portal region mutant indicate opening of the β barrel in fatty acid binding proteins

Background

Myelin is a multilayered proteolipid sheath wrapped around selected axons in the nervous system. Its constituent proteins play major roles in forming of the highly regular membrane structure. P2 is a myelin-specific protein of the fatty acid binding protein (FABP) superfamily, which is able to stack lipid bilayers together, and it is a target for mutations in the human inherited neuropathy Charcot-Marie-Tooth disease. A conserved residue that has been proposed to participate in membrane and fatty acid binding and conformational changes in FABPs is Phe57. This residue is thought to be a gatekeeper for the opening of the portal region upon ligand entry and egress.

Results

We performed a structural characterization of the F57A mutant of human P2. The mutant protein was crystallized in three crystal forms, all of which showed changes in the portal region and helix α2. In addition, the behaviour of the mutant protein upon lipid bilayer binding suggested more unfolding than previously observed for wild-type P2. On the other hand, membrane binding rendered F57A heat-stable, similarly to wild-type P2. Atomistic molecular dynamics simulations showed opening of the side of the discontinuous β barrel, giving important indications on the mechanism of portal region opening and ligand entry into FABPs. The results suggest a central role for Phe57 in regulating the opening of the portal region in human P2 and other FABPs, and the F57A mutation disturbs dynamic cross-correlation networks in the portal region of P2.

Conclusions

Overall, the F57A variant presents similar properties to the P2 patient mutations recently linked to Charcot-Marie-Tooth disease. Our results identify Phe57 as a residue regulating conformational changes that may accompany membrane surface binding and ligand exchange in P2 and other FABPs.

Electronic supplementary material

The online version of this article (10.1186/s12900-018-0087-2) contains supplementary material, which is available to authorized users.

Comparing protein structures with RINspector automation in Cytoscape

Residue interaction networks (RINs) have been shown to be relevant representations of the tertiary or quaternary structures of proteins, in particular thanks to network centrality analyses. We recently developed the RINspector 1.0.0 Cytoscape app, which couples centrality analyses with backbone flexibility predictions. This combined approach permits the identification of crucial residues for the folding or function of the protein that can constitute good targets for mutagenesis experiments. Here we present an application programming interface (API) for RINspector 1.1.0 that enables interplay between Cytoscape, RINspector and external languages, such as R or Python. This API provides easy access to batch centrality calculations and flexibility predictions, and allows for the easy comparison of results between different structures. These comparisons can lead to the identification of specific and conserved central residues, and show the impact of mutations to these and other residues on the flexibility of the proteins. We give two use cases to demonstrate the interest of these functionalities and provide the corresponding scripts: the first concerns NMR conformers, the second focuses on mutations in a structure.