Protein interaction evolution from promiscuity to specificity with reduced flexibility in an increasingly complex network

Key interaction networks: Identifying evolutionarily conserved non-covalent interaction networks across protein families

Protein structure (and thus function) is dictated by non-covalent interaction networks. These can be highly evolutionarily conserved across protein families, the members of which can diverge in sequence and evolutionary history. Here we present KIN, a tool to identify and analyze conserved non-covalent interaction networks across evolutionarily related groups of proteins. KIN is available for download under a GNU General Public License, version 2, from https://www.github.com/kamerlinlab/KIN. KIN can operate on experimentally determined structures, predicted structures, or molecular dynamics trajectories, providing insight into both conserved and missing interactions across evolutionarily related proteins. This provides useful insight both into protein evolution, as well as a tool that can be exploited for protein engineering efforts. As a showcase system, we demonstrate applications of this tool to understanding the evolutionary-relevant conserved interaction networks across the class A β-lactamases.

Investigating the Cytotoxicity of Ru(II) Polypyridyl Complexes by Changing the Electronic Structure of Salicylaldehyde Ligands

A novel class of Ru(II)-based polypyridyl complexes with an auxiliary salicylaldehyde ligand [Ru(phen)2(X-Sal)]BF4 {X: H (1), 5-Cl (2), 5-Br (3), 3,5-Cl2 (4), 3,5-Br2 (5), 3-Br,5-Cl (6), 3,5-I2 (7), 5-NO2 (8), 5-Me (9), 4-Me (10), 4-OMe (11), and 4-DEA (12), has been synthesized and characterized by elemental analysis, FT-IR, and 1H/13C NMR spectroscopy. The molecular structure of 4, 6, 9, 10, and 11 was determined by single-crystal X-ray diffraction analysis which revealed structural similarities. DFT and TD-DFT calculations showed that they also possess similar electronic structures. Absorption/emission spectra were recorded for 2, 3, 10, and 11. All Ru-complexes, unlike the pure ligands and the complex lacking the salicylaldehyde component, displayed outstanding antiproliferative activity in the screening test (10 μM) against CCRF-CEM leukemia cells underlining the crucial role of the presence of the auxiliary ligand for the biological activity. The two most active derivatives, namely 7 and 10, were selected for continuous assays showing IC50 values in the submicromolar and micromolar range against drug-sensitive CCRF-CEM and multidrug-resistant CEM/ADR5000 leukemia cells, respectively. These two compounds were investigated in silico for their potential binding to duplex DNA well-matched and mismatched base pairs, since they showed remarkable selectivity indexes (2.2 and 19.5 respectively) on PBMC cells.

Protein-protein interaction network module changes associated with the vertebrate fin-to-limb transition

Evolutionary phenotypic transitions, such as the fin-to-limb transition in vertebrates, result from modifications in related proteins and their interactions, often in response to changing environment. Identifying these alterations in protein networks is crucial for a more comprehensive understanding of these transitions. However, previous research has not attempted to compare protein-protein interaction (PPI) networks associated with evolutionary transitions, and most experimental studies concentrate on a limited set of proteins. Therefore, the goal of this work was to develop a network-based platform for investigating the fin-to-limb transition using PPI networks. Quality-enhanced protein networks, constructed by integrating PPI networks with anatomy ontology data, were leveraged to compare protein modules for paired fins (pectoral fin and pelvic fin) of fishes (zebrafish) to those of the paired limbs (forelimb and hindlimb) of mammals (mouse). This also included prediction of novel protein candidates and their validation by enrichment and homology analyses. Hub proteins such as shh and bmp4, which are crucial for module stability, were identified, and their changing roles throughout the transition were examined. Proteins with preserved roles during the fin-to-limb transition were more likely to be hub proteins. This study also addressed hypotheses regarding the role of non-preserved proteins associated with the transition.

Transcription factor expression levels and environmental signals constrain transcription factor innovation

Evolutionary innovation of transcription factors frequently drives phenotypic diversification and adaptation to environmental change. Transcription factors can gain or lose connections to target genes, resulting in novel regulatory responses and phenotypes. However the frequency of functional adaptation varies between different regulators, even when they are closely related. To identify factors influencing propensity for innovation, we utilise a Pseudomonas fluorescens SBW25 strain rendered incapable of flagellar mediated motility in soft-agar plates via deletion of the flagellar master regulator (fleQ ). This bacterium can evolve to rescue flagellar motility via gene regulatory network rewiring of an alternative transcription factor to rescue activity of FleQ. Previously, we have identified two members (out of 22) of the RpoN-dependent enhancer binding protein (RpoN-EBP) family of transcription factors (NtrC and PFLU1132) that are capable of innovating in this way. These two transcription factors rescue motility repeatably and reliably in a strict hierarchy – with NtrC the only route in a ∆fleQ background, and PFLU1132 the only route in a ∆fleQ ∆ntrC background. However, why other members in the same transcription factor family have not been observed to rescue flagellar activity is unclear. Previous work shows that protein homology cannot explain this pattern within the protein family (RpoN-EBPs), and mutations in strains that rescued motility suggested high levels of transcription factor expression and activation drive innovation. We predict that mutations that increase expression of the transcription factor are vital to unlock evolutionary potential for innovation. Here, we construct titratable expression mutant lines for 11 of the RpoN-EBPs in P. fluorescens . We show that in five additional RpoN-EBPs (FleR, HbcR, GcsR, DctD, AauR and PFLU2209), high expression levels result in different mutations conferring motility rescue, suggesting alternative rewiring pathways. Our results indicate that expression levels (and not protein homology) of RpoN-EBPs are a key constraining factor in determining evolutionary potential for innovation. This suggests that transcription factors that can achieve high expression through few mutational changes, or transcription factors that are active in the selective environment, are more likely to innovate and contribute to adaptive gene regulatory network evolution.

Cracking the Floral Quartet Code: How Do Multimers of MIKCC-Type MADS-Domain Transcription Factors Recognize Their Target Genes?

MADS-domain transcription factors (MTFs) are involved in the control of many important processes in eukaryotes. They are defined by the presence of a unique and highly conserved DNA-binding domain, the MADS domain. MTFs bind to double-stranded DNA as dimers and recognize specific sequences termed CArG boxes (such as 5'-CC(A/T)₆GG-3') and similar sequences that occur hundreds of thousands of times in a typical flowering plant genome. The number of MTF-encoding genes increased by around two orders of magnitude during land plant evolution, resulting in roughly 100 genes in flowering plant genomes. This raises the question as to how dozens of different but highly similar MTFs accurately recognize the cis-regulatory elements of diverse target genes when the core binding sequence (CArG box) occurs at such a high frequency. Besides the usual processes, such as the base and shape readout of individual DNA sequences by dimers of MTFs, an important sublineage of MTFs in plants, termed MIKC^C-type MTFs (M^C-MTFs), has evolved an additional mechanism to increase the accurate recognition of target genes: the formation of heterotetramers of closely related proteins that bind to two CArG boxes on the same DNA strand involving DNA looping. M^C-MTFs control important developmental processes in flowering plants, ranging from root and shoot to flower, fruit and seed development. The way in which M^C-MTFs bind to DNA and select their target genes is hence not only of high biological interest, but also of great agronomic and economic importance. In this article, we review the interplay of the different mechanisms of target gene recognition, from the ordinary (base readout) via the extravagant (shape readout) to the idiosyncratic (recognition of the distance and orientation of two CArG boxes by heterotetramers of M^C-MTFs). A special focus of our review is on the structural prerequisites of M^C-MTFs that enable the specific recognition of target genes.

Synergistic Combination of Citrus Flavanones as Strong Antioxidant and COX-Inhibitor Agent

Recently, we demonstrated that a Citrus flavanone mix (FM) shows antioxidant and anti-inflammatory activity, even after gastro-duodenal digestion (DFM). The aim of this study was to investigate the possible involvement of the cyclooxygenases (COXs) in the anti-inflammatory activity previously detected, using a human COX inhibitor screening assay, molecular modeling studies, and PGE2 release by Caco-2 cells stimulated with IL-1β and arachidonic acid. Furthermore, the ability to counteract pro-oxidative processes induced by IL-1β was evaluated by measuring four oxidative stress markers, namely, carbonylated proteins, thiobarbituric acid-reactive substances, reactive oxygen species, and reduced glutathione/oxidized glutathione ratio in Caco-2 cells. All flavonoids showed a strong inhibitory activity on COXs, confirmed by molecular modeling studies, with DFM, which showed the best and most synergistic activity on COX-2 (82.45% vs. 87.93% of nimesulide). These results were also corroborated by the cell-based assays. Indeed, DFM proves to be the most powerful anti-inflammatory and antioxidant agent reducing, synergistically and in a statistically significant manner (p < 0.05), PGE2 release than the oxidative stress markers, also with respect to the nimesulide and trolox used as reference compounds. This leads to the hypothesis that FM could be an excellent antioxidant and COX inhibitor candidate to counteract intestinal inflammation.

Antileukemia Activity and Mechanism of Platinum(II)-Based Metal Complexes

Transition metal complexes have continued to constitute an appealing class of medicinal compounds since the exceptional discovery of cisplatin in the late 1960s. Pt(II)-based complexes are endowed with a broad range of biological properties, which are mainly exerted by targeting DNA. In this study, we report a significant biological investigation into and computation analyses of four Pt(II)-complexes, namely, LDP-1-4, synthesized and characterized according to previously reported procedures. Molecular-modelling studies highlighted that the top two LDP compounds (i.e., LDP-1 and LDP-4) might bind to both matched and mismatched base pair sites of the oligonucleotide 5'-(dCGGAAATTACCG)₂-3', supporting their anticancer potential. These two complexes displayed noteworthy cytotoxicity in vitro (sub-micromolar-micromolar range) against two leukaemia cell lines, i.e., CCRF-CEM and its multi-drug-resistant counterpart CEM/ADR5000, and remarkable anti-angiogenic properties (in the sub-micromolar range) evaluated in an in vivo model, i.e., a chick embryo chorioallantoic membrane (CAM) assay.

Bicalutamide Anticancer Activity Enhancement by Formulation of Soluble Inclusion Complexes with Cyclodextrins

Bicalutamide (BCL) is a nonsteroidal antiandrogen drug that represents an alternative to castration in the treatment of prostate cancer, due to its relatively long half-life and tolerable side effects. However, it possesses a very low water solubility that can affect its oral bioavailability. In this work, we developed inclusion complexes of BCL with the highly soluble hydroxypropyl-β-cyclodextrin (HP-β-CyD) and sulfobutylether-β-cyclodextrin (SBE-β-CyD) to increase the water solubility and anticancer activity of BCL. The inclusion complexes were prepared using the freeze-drying method and were then characterized in a solid state via differential scanning calorimetry and X-ray analysis and in solution via phase-solubility studies and UV-vis and NMR spectroscopy. The BCL/HP-β-CyD and BCL/SBE-β-CyD inclusion complexes were amorphous and rapidly dissolved in water. Both the 1H-NMR spectra and molecular modeling studies confirmed the penetration of the 2-(trifluoromethyl)benzonitrile ring of BCL within the cavity of both cyclodextrins (CyDs). Due to the consistent improvement of the water solubility of BCL, the inclusion complexes showed higher antiproliferative activity toward the human prostate androgen-independent cell lines, DU-145 and PC-3, with respect to free BCL. These results demonstrate the ability of HP-β-CyD and SBE-β-CyD to complex BCL, permitting the realization of liquid formulations with potentially high oral bioavailability and/or possible parenteral administration.

αα-hub coregulator structure and flexibility determine transcription factor binding and selection in regulatory interactomes

Formation of transcription factor (TF)-coregulator complexes is a key step in transcriptional regulation, with coregulators having essential functions as hub nodes in molecular networks. How specificity and selectivity are maintained in these nodes remain open questions. In this work, we addressed specificity in transcriptional networks using complexes formed between TFs and αα-hubs, which are defined by a common αα-hairpin secondary structure motif, as a model. Using NMR spectroscopy and binding thermodynamics, we analyzed the structure, dynamics, stability, and ligand-binding properties of the Arabidopsis thaliana RST domains from TAF4 and known binding partner RCD1, and the TAFH domain from human TAF4, allowing comparison across species, functions, and architectural contexts. While these αα-hubs shared the αα-hairpin motif, they differed in length and orientation of accessory helices as well as in their thermodynamic profiles of ligand binding. Whereas biologically relevant RCD1-ligand pairs displayed high affinity driven by enthalpy, TAF4-ligand interactions were entropy driven and exhibited less binding-induced structuring. We in addition identified a thermal unfolding state with a structured core for all three domains, although the temperature sensitivity differed. Thermal stability studies suggested that initial unfolding of the RCD1-RST domain localized around helix 1, lending this region structural malleability, while effects in TAF4-RST were more stochastic, suggesting variability in structural adaptability upon binding. Collectively, our results support a model in which hub structure, flexibility, and binding thermodynamics contribute to αα-hub-TF binding specificity, a finding of general relevance to the understanding of coregulator-ligand interactions and interactome sizes.

Stop CRYing! Inhibition of cryptochrome function by small proteins

Plants can detect the presence of light using specialised photoreceptor proteins. These photoreceptors measure the intensity of light, but they can also respond to different spectra of light and thus ‘see' different colours. Cryptochromes, which are also present in animals, are flavin-based photoreceptors that enable plants to detect blue and ultraviolet-A (UV-A) light. In Arabidopsis, there are two cryptochromes, CRYPTOCHROME 1 (CRY1) and CRYPTOCHROME 2 (CRY2) with known sensory roles. They function in various processes such as blue-light mediated inhibition of hypocotyl elongation, photoperiodic promotion of floral initiation, cotyledon expansion, anthocyanin production, and magnetoreception, to name a few. In the dark, the cryptochromes are in an inactive monomeric state and undergo photochemical and conformational change in response to illumination. This results in flavin reduction, oligomerisation, and the formation of the ‘cryptochrome complexome'. Mechanisms of cryptochrome activation and signalling have been extensively studied and found to be conserved across phylogenetic lines. In this review, we will therefore focus on a far lesser-known mechanism of regulation that is unique to plant cryptochromes. This involves inhibition of cryptochrome activity by small proteins that prevent its dimerisation in response to light. The resulting inhibition of function cause profound alterations in economically important traits such as plant growth, flowering, and fruit production. This review will describe the known mechanisms of cryptochrome activation and signalling in the context of their modulation by these endogenous and artificial small inhibitor proteins. Promising new applications for biotechnological and agricultural applications will be discussed.

Development of isoquinolinone derivatives as immunoproteasome inhibitors

The inhibition of immunoproteasome is considered nowadays a promising strategy for the treatment of hematologic malignancies. In this paper we report the design, synthesis, and biological evaluation as immunoproteasome inhibitors of a new series of isoquinolinone derivatives characterized by a (E)-prop-1-ene fragment that connects the heterocycle to a distal amide functionality. Among all the synthesized compounds, we identified an inhibitor with K_i values in the low micromolar or submicromolar range towards the chymotrypsin-like activities of both proteasome and immunoproteasome (β5c, β5i and β1i subunits). Molecular modeling studies suggest that the most potent compound of the series may act a single-site binder. In particular, through its isopentyl group, it might dock into P1 site in the case of the β1i catalytic subunit, while in the case of β5c and β5i subunits, the P3 site might be the preferred binding site.

Analysis of the structure and function of the tomato Solanum lycopersicum L. MADS-box gene SlMADS5

At all stages of f lowering, a decisive role is played by the family of MADS-domain transcription factors, the combinatorial action of which is described by the ABCDE-model of f lower development. The current volume of data suggests a high conservatism of ABCDE genes in angiosperms. The E-proteins SEPALLATA are the central hub of the MADS-complexes, which determine the identity of the f loral organs. The only representative of the SEPALLATA3 clade in tomato Solanum lycopersicum L., SlMADS5, is involved in determining the identity of petals, stamens, and carpels; however, data on the functions of the gene are limited. The study was focused on the SlMADS5 functional characterization. Structural and phylogenetic analyses of SlMADS5 conf irmed its belonging to the SEP3 clade. An in silico expression analysis revealed the absence of gene transcripts in roots, leaves, and shoot apical meristem, and their presence in f lowers, fruits, and seeds at different stages of development. Two-hybrid analysis showed the ability of SlMADS5 to activate transcription of the target gene and interact with TAGL1. Transgenic plants Nicotiana tabacum L. with constitutive overexpression of SlMADS5 cDNA f lowered 2.2 times later than the control; plants formed thickened leaves, 2.5–3.0 times thicker stems, 1.5–2.7 times shortened internodes, and 1.9 times fewer f lowers and capsules than non-transgenic plants. The f lower structure did not differ from the control; however, the corolla petals changed color from light pink to magenta. Analysis of the expression of SlMADS5 and the tobacco genes NtLFY, NtAP1, NtWUS, NtAG, NtPLE, NtSEP1, NtSEP2, and NtSEP3 in leaves and apexes of transgenic and control plants showed that SlMADS5 mRNA is present only in tissues of transgenic lines. The other genes analyzed were highly expressed in the reproductive meristem of control plants. Gene transcripts were absent or were imperceptibly present in the leaves and vegetative apex of the control, as well as in the leaves and apexes of transgenic lines. The results obtained indicate the possible involvement of SlMADS5 in the regulation of f lower meristem development and the pathway of anthocyanin biosynthesis in petals.

Genome-Wide Identification and Analysis of the MADS-Box Gene Family in American Beautyberry (Callicarpa americana)

The MADS-box gene family encodes a number of transcription factors that play key roles in various plant growth and development processes from response to environmental cues to cell differentiation and organ identity, especially the floral organogenesis, as in the prominent ABCDE model of flower development. Recently, the genome of American beautyberry (Callicarpa americana) has been sequenced. It is a shrub native to the southern region of United States with edible purple-colored berries; it is a member of the Lamiaceae family, a family of medical and agricultural importance. Seventy-eight MADS-box genes were identified from 17 chromosomes of the C. americana assembled genome. Peptide sequences blast and analysis of phylogenetic relationships with MADS-box genes of Sesame indicum, Solanum lycopersicum, Arabidopsis thaliana, and Amborella trichopoda were performed. Genes were separated into 32 type I and 46 type II MADS-box genes. C. americana MADS-box genes were clustered into four groups: MIKC^C, MIKC*, Mα-type, and Mγ-type, while the Mβ-type group was absent. Analysis of the gene structure revealed that from 1 to 15 exons exist in C. americana MADS-box genes. The number of exons in type II MADS-box genes (5–15) greatly exceeded the number in type I genes (1–9). The motif distribution analysis of the two types of MADS-box genes showed that type II MADS-box genes contained more motifs than type I genes. These results suggested that C. americana MADS-box genes type II had more complex structures and might have more diverse functions. The role of MIKC-type MADS-box genes in flower and fruit development was highlighted when the expression profile was analyzed in different organs transcriptomes. This study is the first genome-wide analysis of the C. americana MADS-box gene family, and the results will further support any functional and evolutionary studies of C. americana MADS-box genes and serve as a reference for related studies of other plants in the medically important Lamiaceae family.

Were Ancestral Proteins Less Specific?

Some have hypothesized that ancestral proteins were, on average, less specific than their descendants. If true, this would provide a universal axis along which to organize protein evolution and suggests that reconstructed ancestral proteins may be uniquely powerful tools for protein engineering. Ancestral sequence reconstruction studies are one line of evidence used to support this hypothesis. Previously, we performed such a study, investigating the evolution of peptide-binding specificity for the paralogs S100A5 and S100A6. The modern proteins appeared more specific than their last common ancestor (ancA5/A6), as each paralog bound a subset of the peptides bound by ancA5/A6. In this study, we revisit this transition, using quantitative phage display to measure the interactions of 30,533 random peptides with human S100A5, S100A6, and ancA5/A6. This unbiased screen reveals a different picture. While S100A5 and S100A6 do indeed bind to a subset of the peptides recognized by ancA5/A6, they also acquired new peptide partners outside of the set recognized by ancA5/A6. Our previous work showed that ancA5/A6 had lower specificity than its descendants when measured against biological targets; our new work shows that ancA5/A6 has similar specificity to the modern proteins when measured against a random set of peptide targets. This demonstrates that altered biological specificity does not necessarily indicate altered intrinsic specificity, and sounds a cautionary note for using ancestral reconstruction studies with biological targets as a means to infer global evolutionary trends in specificity.

Interaction of selected terpenoids with two SARS-CoV-2 key therapeutic targets: An in silico study through molecular docking and dynamics simulations

The outbreak of COVID-19 disease caused by SARS-CoV-2, along with the lack of targeted medicaments, forced the scientific world to search for new antiviral formulations. In the current emergent situation, drug repurposing of well-known traditional and/or approved drugs could be the most effective strategy. Herein, through computational approaches, we aimed to screen 14 natural compounds from limonoids and terpenoids class for their ability to inhibit the key therapeutic target proteins of SARS-CoV-2. Among these, some limonoids, namely deacetylnomilin, ichangin and nomilin, and the terpenoid β-amyrin provided good interaction energies with SARS-CoV-2 3CL hydrolase (Mpro) in molecular dynamic simulation. Interestingly, deacetylnomilin and ichangin showed direct interaction with the catalytic dyad of the enzyme so supporting their potential role in preventing SARS-CoV-2 replication and growth. On the contrary, despite the good affinity with the spike protein RBD site, all the selected phytochemicals lose contact with the amino acid residues over the course of 120ns-long molecular dynamics simulations therefore suggesting they scarcely can interfere in SARS-CoV-2 binding to the ACE2 receptor. The in silico analyses of docking score and binding energies, along with predicted pharmacokinetic profiles, indicate that these triterpenoids might have potential as inhibitors of SARS-CoV-2 Mpro, recommending further in vitro and in vivo investigations for a complete understanding and confirmation of their inhibitory potential.

On the specificity of protein-protein interactions in the context of disorder

With the increased focus on intrinsically disordered proteins (IDPs) and their large interactomes, the question about their specificity — or more so on their multispecificity — arise. Here we recapitulate how specificity and multispecificity are quantified and address through examples if IDPs in this respect differ from globular proteins. The conclusion is that quantitatively, globular proteins and IDPs are similar when it comes to specificity. However, compared with globular proteins, IDPs have larger interactome sizes, a phenomenon that is further enabled by their flexibility, repetitive binding motifs and propensity to adapt to different binding partners. For IDPs, this adaptability, interactome size and a higher degree of multivalency opens for new interaction mechanisms such as facilitated exchange through trimer formation and ultra-sensitivity via threshold effects and ensemble redistribution. IDPs and their interactions, thus, do not compromise the definition of specificity. Instead, it is the sheer size of their interactomes that complicates its calculation. More importantly, it is this size that challenges how we conceptually envision, interpret and speak about their specificity.

αα-Hub domains and intrinsically disordered proteins: A decisive combo

Hub proteins are central nodes in protein-protein interaction networks with critical importance to all living organisms. Recently, a new group of folded hub domains, the αα-hubs, was defined based on a shared αα-hairpin supersecondary structural foundation. The members PAH, RST, TAFH, NCBD, and HHD are found in large proteins such as Sin3, RCD1, TAF4, CBP, and harmonin, which organize disordered transcriptional regulators and membrane scaffolds in interactomes of importance to human diseases and plant quality. In this review, studies of structures, functions, and complexes across the αα-hubs are described and compared to provide a unified description of the group. This analysis expands the associated molecular concepts of "one domain-one binding site", motif-based ligand binding, and coupled folding and binding of intrinsically disordered ligands to additional concepts of importance to signal fidelity. These include context, motif reversibility, multivalency, complex heterogeneity, synergistic αα-hub:ligand folding, accessory binding sites, and supramodules. We propose that these multifaceted protein-protein interaction properties are made possible by the characteristics of the αα-hub fold, including supersite properties, dynamics, variable topologies, accessory helices, and malleability and abetted by adaptability of the disordered ligands. Critically, these features provide additional filters for specificity. With the presentations of new concepts, this review opens for new research questions addressing properties across the group, which are driven from concepts discovered in studies of the individual members. Combined, the members of the αα-hubs are ideal models for deconvoluting signal fidelity maintained by folded hubs and their interactions with intrinsically disordered ligands.

A Fifth of the Protein World: Rossmann-like Proteins as an Evolutionarily Successful Structural unit

The Rossmann-like fold is the most prevalent and diversified doubly-wound superfold of ancient evolutionary origin. Rossmann-like domains are present in a variety of metabolic enzymes and are capable of binding diverse ligands. Discerning evolutionary relationships among these domains is challenging because of their diverse functions and ancient origin. We defined a minimal Rossmann-like structural motif (RLM), identified RLM-containing domains among known 3D structures (20%) and classified them according to their homologous relationships. New classifications were incorporated into our Evolutionary Classification of protein Domains (ECOD) database. We defined 156 homology groups (H-groups), which were further clustered into 123 possible homology groups (X-groups). Our analysis revealed that RLM-containing proteins constitute approximately 15% of the human proteome. We found that disease-causing mutations are more frequent within RLM domains than within non-RLM domains of these proteins, highlighting the importance of RLM-containing proteins for human health.

Exposing the distinctive modular behavior of β-strands and α-helices in folded proteins

Although folded proteins are commonly depicted as simplistic combinations of β-strands and α-helices, the actual properties and functions of these secondary-structure elements in their native contexts are just partly understood. The principal reason is that the behavior of individual β- and α-elements is obscured by the global folding cooperativity. In this study, we have circumvented this problem by designing frustrated variants of the mixed α/β-protein S6, which allow the structural behavior of individual β-strands and α-helices to be targeted selectively by stopped-flow kinetics, X-ray crystallography, and solution-state NMR. Essentially, our approach is based on provoking intramolecular "domain swap." The results show that the α- and β-elements have quite different characteristics: The swaps of β-strands proceed via global unfolding, whereas the α-helices are free to swap locally in the native basin. Moreover, the α-helices tend to hybridize and to promote protein association by gliding over to neighboring molecules. This difference in structural behavior follows directly from hydrogen-bonding restrictions and suggests that the protein secondary structure defines not only tertiary geometry, but also maintains control in function and structural evolution. Finally, our alternative approach to protein folding and native-state dynamics presents a generally applicable strategy for in silico design of protein models that are computationally testable in the microsecond-millisecond regime.

Evolutionary Variation in MADS Box Dimerization Affects Floral Development and Protein Abundance in Maize

Interactions between MADS box transcription factors are critical in the regulation of floral development, and shifting MADS box protein-protein interactions are predicted to have influenced floral evolution. However, precisely how evolutionary variation in protein-protein interactions affects MADS box protein function remains unknown. To assess the impact of changing MADS box protein-protein interactions on transcription factor function, we turned to the grasses, where interactions between B-class MADS box proteins vary. We tested the functional consequences of this evolutionary variability using maize (Zea mays) as an experimental system. We found that differential B-class dimerization was associated with subtle, quantitative differences in stamen shape. In contrast, differential dimerization resulted in large-scale changes to downstream gene expression. Differential dimerization also affected B-class complex composition and abundance, independent of transcript levels. This indicates that differential B-class dimerization affects protein degradation, revealing an important consequence for evolutionary variability in MADS box interactions. Our results highlight complexity in the evolution of developmental gene networks: changing protein-protein interactions could affect not only the composition of transcription factor complexes but also their degradation and persistence in developing flowers. Our results also show how coding change in a pleiotropic master regulator could have small, quantitative effects on development.

Antibody Clustering Using a Machine Learning Pipeline that Fuses Genetic, Structural, and Physicochemical Properties

Antibody V domain clustering is of paramount importance to a repertoire of immunology-related areas. Although several approaches have been proposed for antibody clustering, still no consensus has been reached. Numerous attempts use information from genes, protein sequences, 3D structures, and 3D surfaces in an effort to elucidate unknown action mechanisms directly related to their function and to either link them directly to diseases or drive the discovery of new medicines, such as antibody drug conjugates (ADC). Herein, we describe a new V domain antibody clustering method based on the comparison of the interaction sites between each antibody and its antigen. A more specific clustering analysis of the antibody's V domain was provided using deep learning and data mining techniques. The multidimensional information was extracted from the structural resolved antibodies when they were captured to interact with other proteins. The available 3D structures of protein antigen-antibody (Ag-Ab) interfaces contain information about how antibody V domains recognize antigens as well as about which amino acids are involved in the recognition. As such, the antibody surface holds information about antigens' folding that reside with the Ab-Ag interface residues and how they interact. In order to gain insight into the nature of such interactions, we propose a new simple philosophy to transform the conserved framework (fragment regions, complementarity-determining regions) of antibody V domain in a binary form using structural features of antibody-antigen interactions, toward identifying new antibody signatures in V domain binding activity. Finally, an advanced three-level hybrid classification scheme has been set for clustering antibodies in subgroups, which can combine the information from the protein sequences, the three-dimensional structures, and specific "key patterns" of recognized interactions. The clusters provide multilevel information about antibodies and antibody-antigen complexes.

Involvement of Local, Rapid Conformational Dynamics in Binding of Flexible Recognition Motifs

Flexible protein sequences populate ensembles of rapidly interconverting states differentiated by small-scale fluctuations; however, elucidating whether and how the ensembles determine function experimentally is challenged by the combined high spatial and temporal resolution needed to capture the states. We used carbon-deuterium (C-D) bond vibrations incorporated as infrared probes to characterize with residue-specific detail the heterogeneity of states adopted by proline-rich (PR) sequences and assess their involvement in recognition of Src homology 3 domains. The C-D absorption envelopes provided evidence for two or three sub-populations at all proline residues. The changes in the subpopulations induced by binding generally reflected recognition by conformational selection but depended on the residue and the state of the ligand to illuminate distinct mechanisms among the PR ligands. Notably, the spectral data indicate that greater adaptability among the states is associated with reduced recognition specificity and that perturbation to the ensemble populations contributes to differences in binding entropy. Broadly, the study quantifies rapidly interconverting ensembles with residue-specific detail and implicates them in function.

Structural Basis for Plant MADS Transcription Factor Oligomerization

MADS transcription factors (TFs) are DNA binding proteins found in almost all eukaryotes that play essential roles in diverse biological processes. While present in animals and fungi as a small TF family, the family has dramatically expanded in plants over the course of evolution, with the model flowering plant, Arabidopsis thaliana, possessing over 100 type I and type II MADS TFs. All MADS TFs contain a core and highly conserved DNA binding domain called the MADS or M domain. Plant MADS TFs have diversified this domain with plant-specific auxiliary domains. Plant type I MADS TFs have a highly diverse and largely unstructured Carboxy-terminal (C domain), whereas type II MADS have added oligomerization domains, called Intervening (I domain) and Keratin-like (K domain), in addition to the C domain. In this mini review, we describe the overall structure of the type II "MIKC" type MADS TFs in plants, with a focus on the K domain, a critical oligomerization module. We summarize the determining factors for oligomerization and provide mechanistic insights on how secondary structural elements are required for oligomerization capability and specificity. Using MADS TFs that are involved in flower organ specification as an example, we provide case studies and homology modeling of MADS TFs complex formation. Finally, we highlight outstanding questions in the field.

Gene duplication and relaxation from selective constraints of GCYC genes correlated with various floral symmetry patterns in Asiatic Gesneriaceae tribe Trichosporeae

Floral bilateral symmetry is one of the most important acquisitions in flower shape evolution in angiosperms. Members of Gesneriaceae possess predominantly zygomorphic flowers yet natural reversal to actinomorphy have independently evolved multiple times. The development of floral bilateral symmetry relies greatly on the gene CYCLOIDEA (CYC). Our reconstructed GCYC phylogeny indicated at least five GCYC duplication events occurred over the evolutionary history of Gesneriaceae. However, the patterns of GCYC expression following the duplications and the role of natural selection on GCYC copies in relation to floral symmetry remained largely unstudied. The Asiatic tribe Trichosporeae contains most reversals to actinomorphy. We thus investigated shifts in GCYC gene expression among selected zygomorphic species (Hemiboea bicornuta and Lysionotus pauciflorus) and species with reversals to actinomorphy (Conandron ramondioides) by RT-PCR. In the actinomorphic C. ramondioides, none of the three copies of GCYC was found expressed in petals implying that the reversal was a loss-of-function event. On the other hand, both zygomorphic species retained one GCYC1 copy that was expressed in the dorsal petals but each species utilized a different copy (GCYC1C for H. bicornuta and GCYC1D for L. pauciflorus). Together with previously published data, it appeared that GCYC1C and GCYC1D copies diversified their expression in a distinct species-specific pattern. To detect whether the selection signal (ω) changed before and after the duplication of GCYC1 in Asiatic Trichosporeae, we reconstructed a GCYC phylogeny using maximum likelihood and Bayesian inference algorithms and examined selection signals using PAML. The PAML analysis detected relaxation from selection right after the GCYC1 duplication (ωpre-duplication = 0.2819, ωpost-duplication = 0.3985) among Asiatic Trichosporeae species. We propose that the selection relaxation after the GCYC1 duplication created an "evolutionary window of flexibility" in which multiple copies were retained with randomly diverged roles for dorsal-specific expressions in association with floral symmetry changes.

A two-level model for the role of complex and young genes in the formation of organism complexity and new insights into the relationship between evolution and development

Background

How genome complexity affects organismal phenotypic complexity is a fundamental question in evolutionary developmental biology. Previous studies proposed various contributing factors of genome complexity and tried to find the connection between genomic complexity and organism complexity. However, a general model to answer this question is lacking. Here, we introduce a ‘two-level’ model for the realization of genome complexity at phenotypic level.

Results

Five representative species across Protostomia and Deuterostomia were involved in this study. The intrinsic gene properties contributing to genome complexity were classified into two generalized groups: the complexity and age degree of both protein-coding and noncoding genes. We found that young genes tend to be simpler; however, the mid-age genes, rather than the oldest genes, show the highest proportion of high complexity. Complex genes tend to be utilized preferentially in each stage of embryonic development, with maximum representation during the late stage of organogenesis. This trend is mainly attributed to mid-age complex genes. In contrast, young genes tend to be expressed in specific spatiotemporal states. An obvious correlation between the time point of the change in over- and under-representation and the order of gene age was observed, which supports the funnel-like model of the conservation pattern of development. In addition, we found some probable causes for the seemingly contradictory ‘funnel-like’ or ‘hourglass’ model.

Conclusions

These results indicate that complex and young genes contribute to organismal complexity at two different levels: Complex genes contribute to the complexity of individual proteomes in certain states, whereas young genes contribute to the diversity of proteomes in different spatiotemporal states. This conclusion is valid across the five species investigated, indicating it is a conserved model across Protostomia and Deuterostomia. The results in this study also support ‘funnel-like model’ from a new viewpoint and explain why there are different evo–devo relation models.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0111-4) contains supplementary material, which is available to authorized users.

Resurrected Protein Interaction Networks Reveal the Innovation Potential of Ancient Whole-Genome Duplication

The evolution of plants is characterized by whole-genome duplications, sometimes closely associated with the origin of large groups of species. The gamma (γ) genome triplication occurred at the origin of the core eudicots, which comprise ∼75% of flowering plants. To better understand the impact of whole-genome duplication, we studied the protein interaction network of MADS domain transcription factors, which are key regulators of reproductive development. We reconstructed, synthesized, and tested the interactions of ancestral proteins immediately before and closely after the triplication and directly compared these ancestral networks to the extant networks of Arabidopsis thaliana and tomato (Solanum lycopersicum). We found that gamma expanded the MADS domain interaction network more strongly than subsequent genomic events. This event strongly rewired MADS domain interactions and allowed for the evolution of new functions and installed robustness through new redundancy. Despite extensive rewiring, the organization of the network was maintained through gamma. New interactions and protein retention compensated for its potentially destructive impact on network organization. Post gamma, the network evolved from an organization around the single hub SEP3 to a network organized around multiple hubs and well-connected proteins lost, rather than gained, interactions. The data provide a resource for comparative developmental biology in flowering plants.

Sporadic cases of lumpy skin disease among cattle in Sharkia province, Egypt: Genetic characterization of lumpy skin disease virus isolates and pathological findings

BACKGROUND AND AIM

Lumpy skin disease (LSD) is a highly infectious viral disease upsetting cattle, caused by LSD virus (LSDV) within the family Poxviridae. Sporadic cases of LSD have been observed in cattle previously vaccinated with the Romanian sheep poxvirus (SPPV) vaccine during the summer of 2016 in Sharkia province, Egypt. The present study was undertaken to perform molecular characterization of LSDV strains which circulated in this period as well as investigate their phylogenetic relatedness with published reference capripoxvirus genome sequences.

MATERIALS AND METHODS

A total of 82 skin nodules, as well as 5 lymph nodes, were collected from suspect LSD cases, and the virus was isolated in embryonated chicken eggs (ECEs). LSD was confirmed by polymerase chain reactions amplification of the partial and full-length sequences of the attachment and G-protein-coupled chemokine receptor (GPCR) genes, respectively, as well as a histopathological examination of the lesions. Molecular characterization of the LSDV isolates was conducted by sequencing the GPCR gene.

RESULTS

Characteristic skin nodules that covered the whole intact skin, as well as lymphadenopathy, were significant clinical signs in all suspected cases. LSDV isolation in ECEs revealed the characteristic focal white pock lesions dispersed on the chorioallantoic membranes. Histopathologic examination showed characteristic eosinophilic intracytoplasmic inclusion bodies within inflammatory cell infiltration. Phylogenetic analysis revealed that the LSDV isolates were clustered together with other African and European LSDV strains. In addition, the LSDV isolates have a unique signature of LSDVs (A11, T12, T34, S99, and P199).

CONCLUSION

LSDV infections have been detected in cattle previously vaccinated with Romanian SPPV vaccine during the summer of 2016 and making the evaluation of vaccine efficacy under field conditions necessary.

MADS-box genes and crop domestication: the jack of all traits

MADS-box genes are key regulators of virtually every aspect of plant reproductive development. They play especially prominent roles in flowering time control, inflorescence architecture, floral organ identity determination, and seed development. The developmental and evolutionary importance of MADS-box genes is widely acknowledged. However, their role during flowering plant domestication is less well recognized. Here, we provide an overview illustrating that MADS-box genes have been important targets of selection during crop domestication and improvement. Numerous examples from a diversity of crop plants show that various developmental processes have been shaped by allelic variations in MADS-box genes. We propose that new genomic and genome editing resources provide an excellent starting point for further harnessing the potential of MADS-box genes to improve a variety of reproductive traits in crops. We also suggest that the biophysics of MADS-domain protein-protein and protein-DNA interactions, which is becoming increasingly well characterized, makes them especially suited to exploit coding sequence variations for targeted breeding approaches.

Real Talk: The Inter-play Between the mTOR, AMPK, and Hexosamine Biosynthetic Pathways in Cell Signaling

O-linked N-acetylglucosamine, better known as O-GlcNAc, is a sugar post-translational modification participating in a diverse range of cell functions. Disruptions in the cycling of O-GlcNAc mediated by O-GlcNAc transferase (OGT) and O-GlcNAcase (OGA), respectively, is a driving force for aberrant cell signaling in disease pathologies, such as diabetes, obesity, Alzheimer's disease, and cancer. Production of UDP-GlcNAc, the metabolic substrate for OGT, by the Hexosamine Biosynthetic Pathway (HBP) is controlled by the input of amino acids, fats, and nucleic acids, making O-GlcNAc a key nutrient-sensor for fluctuations in these macromolecules. The mammalian target of rapamycin (mTOR) and AMP-activated protein kinase (AMPK) pathways also participate in nutrient-sensing as a means of controlling cell activity and are significant factors in a variety of pathologies. Research into the individual nutrient-sensitivities of the HBP, AMPK, and mTOR pathways has revealed a complex regulatory dynamic, where their unique responses to macromolecule levels coordinate cell behavior. Importantly, cross-talk between these pathways fine-tunes the cellular response to nutrients. Strong evidence demonstrates that AMPK negatively regulates the mTOR pathway, but O-GlcNAcylation of AMPK lowers enzymatic activity and promotes growth. On the other hand, AMPK can phosphorylate OGT leading to changes in OGT function. Complex sets of interactions between the HBP, AMPK, and mTOR pathways integrate nutritional signals to respond to changes in the environment. In particular, examining these relationships using systems biology approaches might prove a useful method of exploring the complex nature of cell signaling. Overall, understanding the complex interactions of these nutrient pathways will provide novel mechanistic information into how nutrients influence health and disease.

Conservation of Specificity in Two Low-Specificity Proteins

Many regulatory proteins bind peptide regions of target proteins and modulate their activity. Such regulatory proteins can often interact with highly diverse target peptides. In many instances, it is not known if the peptide-binding interface discriminates targets in a biological context, or whether biological specificity is achieved exclusively through external factors such as subcellular localization. We used an evolutionary biochemical approach to distinguish these possibilities for two such low-specificity proteins: S100A5 and S100A6. We used isothermal titration calorimetry to study the binding of peptides with diverse sequence and biochemistry to human S100A5 and S100A6. These proteins bound distinct, but overlapping, sets of peptide targets. We then studied the peptide binding properties of orthologs sampled from across five amniote species. Binding specificity was conserved along all lineages, for the last 320 million years, despite the low specificity of each protein. We used ancestral sequence reconstruction to determine the binding specificity of the last common ancestor of the paralogs. The ancestor bound the entire set of peptides bound by modern S100A5 and S100A6 proteins, suggesting that paralog specificity evolved via subfunctionalization. To rule out the possibility that specificity is conserved because it is difficult to modify, we identified a single historical mutation that, when reverted in human S100A5, gave it the ability to bind an S100A6-specific peptide. These results reveal strong evolutionary constraints on peptide binding specificity. Despite being able to bind a large number of targets, the specificity of S100 peptide interfaces is likely important for the biology of these proteins.