Stability of domain structures in multi-domain proteins

Fibrinogen-binding M-related proteins facilitate the recruitment of plasminogen by Streptococcus pyogenes

Group A Streptococcus (GAS) M-related proteins (Mrp) are dimeric α-helical coiled-coil cell-wall-attached proteins. During infection, Mrp recruit human fibrinogen (Fg) to the bacterial surface, enhancing phagocytosis resistance and promoting growth in human blood. However, Mrp exhibit a high degree of sequence diversity, clustering into four evolutionarily distinct groups. It is currently unknown whether this diversity affects the host-pathogen interactions mediated by Mrp. In this study, nine Mrp sequences from the four major evolutionary groups were selected to examine the effect of sequence diversity on protein-protein interactions with Fg. Negative staining transmission electron microscopy confirmed that Mrp are fibrillar proteins measuring between 45.4 and 47.3 nm in length, and mass photometry confirmed the ability of Mrp to form dimers. Surface plasmon resonance was used to evaluate the affinity of each Mrp for Fg. All Mrp studied bound to Fg via Fragment D (FgD) with nanomolar affinity. Previous studies have linked the acquisition of plasminogen (Plg) by GAS Fg-binding M proteins to tissue destruction and excessive stimulation of the human inflammatory response during infection. Our findings show that Mrp provide an alternative mechanism for Plg recruitment, as Plg binding by Mrp was significantly enhanced following pre-incubation with Fg. These data suggest that Mrp play an important role in GAS host-pathogen interactions. However, further studies are necessary to investigate the relevance of these findings in vivo.

Comparative analysis of permanent and transient domain-domain interactions in multi-domain proteins

Protein domains are structural, functional, and evolutionary units. These domains bring out the diversity of functionality by means of interactions with other co-existing domains and provide stability. Hence, it is important to study intra-protein inter-domain interactions from the perspective of types of interactions. Domains within a chain could interact over short timeframes or permanently, rather like protein-protein interactions (PPIs). However, no systematic study has been carried out between two classes, namely permanent and transient domain-domain interactions. In this work, we studied 263 two-domain proteins, belonging to either of these classes and their interfaces on the basis of several factors, such as interface area and details of interactions (number, strength, and types of interactions). We also characterized them based on residue conservation at the interface, correlation of residue motions across domains, its involvement in repeat formation, and their involvement in particular molecular processes. Finally, we could analyze the interactions arising from domains in two-domain monomeric proteins, and we observed significant differences between these two classes of domain interactions and a few similarities. This study will help to obtain a better understanding of structure-function and folding principles of multi-domain proteins.

Bio-conversion of whey lactose using enzymatic hydrolysis with β-galactosidase: an experimental and kinetic study

Lactose in cheese whey is increasingly challenging to metabolise under normal conditions. The hydrolysis of whey lactose into glucose and galactose using enzymatic methods has been acclaimed to confer benefits like enhanced substrate availability for better degradation in anaerobic digestion. In the present study, whey lactose was subjected to hydrolysis using the enzyme β-galactosidase derived from Aspergillus oryzae fungus to reduce the difficulty of lipid and fat transformation in anaerobic digestion. The individual and combined effects of hydrolysis parameters, pH, enzyme load, reaction time and temperature were studied using Response Surface Methodology by Central Composite Design. The optimum conditions were determined based on variance analyses and surface plots; pH 4.63, temperature 40.47°C, reaction time 25.96 min and enzyme load 0.49%. Results showed a maximum lactose hydrolysis value of 86.21%, while the predicted value was 87.44%. Indeed, enzyme hydrolysis induced a change of soluble chemical oxygen demand around 24.6% and 75.8% reduction in volatile fatty acid concentration. Upon anaerobic digestion, the pre-hydrolysed whey revealed a 3.6-fold higher bio-methane production than that of raw hey, and a visible decrease in volatile fatty acid concentrations. The resultant data agreed with the Gompertz model, and lag phase times were significantly reduced for hydrolysed whey.

Inter-Domain Repulsion of Dumbbell-Shaped Calmodulin during Electrospray Ionization Revealed by Molecular Dynamics Simulations

The mechanisms whereby protein ions are released from nanodroplets at the liquid-gas interface have continued to be controversial since electrospray ionization (ESI) mass spectrometry was widely applied in biomolecular structure analysis in solution. Several viable pathways have been proposed and verified for single-domain proteins. However, the ESI mechanism of multi-domain proteins with more complicated and flexible structures remains unclear. Herein, dumbbell-shaped calmodulin was chosen as a multi-domain protein model to perform molecular dynamics simulations to investigate the structural evolution during the ESI process. For [Ca₄CAM], the protein followed the classical charge residue model. As the inter-domain electrostatic repulsion increased, the droplet was found to split into two sub-droplets, while stronger-repulsive apo-calmodulin unfolded during the early evaporation stage. We designated this novel ESI mechanism as the domain repulsion model, which provides new mechanistic insights into further exploration of proteins containing more domains. Our results suggest that greater attention should be paid to the effect of domain-domain interactions on structure retention during liquid-gas interface transfer when mass spectrometry is used as the developing technique in gas phase structural biology.

Domain-wise dissection of thermal stability enhancement in multidomain proteins

Stability is critical for the proper functioning of all proteins. Optimization of protein thermostability is a key step in the development of industrial enzymes and biologics. Herein, we demonstrate that multidomain proteins can be stabilized significantly using domain-based engineering followed by the recombination of the optimized domains. Domain-level analysis of designed protein variants with similar structures but different thermal profiles showed that the independent enhancement of the thermostability of a constituent domain improves the overall stability of the whole multidomain protein. The crystal structure and AlphaFold-predicted model of the designed proteins via domain-recombination provided a molecular explanation for domain-based stepwise stabilization. Our study suggests that domain-based modular engineering can minimize the sequence space for calculations in computational design and experimental errors, thereby offering useful guidance for multidomain protein engineering.

No evidence for widespread positive selection on double substitutions within codons in primates and yeasts

Nucleotide substitutions in protein-coding genes can be divided into synonymous (S) and non-synonymous (N) ones that alter amino acids (including nonsense mutations causing stop codons). The S substitutions are expected to have little effect on function. The N substitutions almost always are affected by strong purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases can modulate the deleterious effect of single N substitutions and, thus, could be subjected to the positive selection. This effect has been demonstrated for mutations in the serine codons, stop codons and double N substitutions in prokaryotes. In all abovementioned cases, a novel technique was applied that allows elucidating the effects of selection on double substitutions considering mutational biases. Here, we applied the same technique to study double N substitutions in eukaryotic lineages of primates and yeast. We identified markedly fewer cases of purifying selection relative to prokaryotes and no evidence of codon double substitutions under positive selection. This is consistent with previous studies of serine codons in primates and yeast. In general, the obtained results strongly suggest that there are major differences between studied pro- and eukaryotes; double substitutions in primates and yeasts largely reflect mutational biases and are not hallmarks of selection. This is especially important in the context of detection of positive selection in codons because it has been suggested that multiple mutations in codons cause false inferences of lineage-specific site positive selection. It is likely that this concern is applicable to previously studied prokaryotes but not to primates and yeasts where markedly fewer double substitutions are affected by positive selection.

Identification of 121 variants of honey bee Vitellogenin protein sequences with structural differences at functional sites

Proteins are under selection to maintain central functions and to accommodate needs that arise in ever-changing environments. The positive selection and neutral drift that preserve functions result in a diversity of protein variants. The amount of diversity differs between proteins: multifunctional or disease-related proteins tend to have fewer variants than proteins involved in some aspects of immunity. Our work focuses on the extensively studied protein Vitellogenin (Vg), which in honey bees (Apis mellifera) is multifunctional and highly expressed and plays roles in immunity. Yet, almost nothing is known about the natural variation in the coding sequences of this protein or how amino acid-altering variants might impact structure-function relationships. Here, we map out allelic variation in honey bee Vg using biological samples from 15 countries. The successful barcoded amplicon Nanopore sequencing of 543 bees revealed 121 protein variants, indicating a high level of diversity in Vg. We find that the distribution of non-synonymous single nucleotide polymorphisms (nsSNPs) differs between protein regions with different functions; domains involved in DNA and protein-protein interactions contain fewer nsSNPs than the protein's lipid binding cavities. We outline how the central functions of the protein can be maintained in different variants and how the variation pattern may inform about selection from pathogens and nutrition.

Structural insights into the role of the WW2 domain on tandem WW/PPxY-motif interactions of oxidoreductase WWOX

Class I WW domains are present in many proteins of various functions and mediate protein interactions by binding to short linear PPxY motifs. Tandem WW domains often bind peptides with multiple PPxY motifs, but the interplay of WW–peptide interactions is not always intuitive. The WW domain–containing oxidoreductase (WWOX) harbors two WW domains: an unstable WW1 capable of PPxY binding and stable WW2 that cannot bind PPxY. The WW2 domain has been suggested to act as a WW1 domain chaperone, but the underlying mechanism of its chaperone activity remains to be revealed. Here, we combined NMR, isothermal calorimetry, and structural modeling to elucidate the roles of both WW domains in WWOX binding to its PPxY-containing substrate ErbB4. Using NMR, we identified an interaction surface between these two domains that supports a WWOX conformation compatible with peptide substrate binding. Isothermal calorimetry and NMR measurements also indicated that while binding affinity to a single PPxY motif is marginally increased in the presence of WW2, affinity to a dual-motif peptide increases 10-fold. Furthermore, we found WW2 can directly bind double-motif peptides using its canonical binding site. Finally, differential binding of peptides in mutagenesis experiments was consistent with a parallel N- to C-terminal PPxY tandem motif orientation in binding to the WW1–WW2 tandem domain, validating structural models of the interaction. Taken together, our results reveal the complex nature of tandem WW-domain organization and substrate binding, highlighting the contribution of WWOX WW2 to both protein stability and target binding.

Structural and dynamic studies of the peptidase domain from Clostridium thermocellum PCAT1

The export of antimicrobial peptides is mediated by diverse mechanisms in bacterial quorum sensing pathways. One such binary system employed by gram-positive bacteria is the PCAT1 ABC transporter coupled to a cysteine protease. The focus of this study is the N-terminal C39 peptidase (PEP) domain from Clostridium thermocellum PCAT1 that processes its natural substrate CtA by cleaving a conserved -GG- motif to separate the cargo from the leader peptide prior to secretion. In this study, we are primarily interested in elucidating the dynamic and structural determinants of CtA binding and how it is coupled to cleavage efficiency in the PCAT1 PEP domain. To this end, we have characterized CtA interactions with PEP domain and PCAT1 transporter in detergent micelles using solution nuclear magnetic resonance spectroscopy. The bound CtA structure revealed the disordered C-terminal cargo peptide is linked by a sterically hindered cleavage site to a helix docked within a hydrophobic cavity in the PEP domain. The wide range of internal motions detected by amide nitrogen (N¹⁵ ) relaxation measurements in the free enzyme and substrate-bound complex suggests the binding site is relatively floppy. This flexibility plays a key role in the structural rearrangement necessary to relax steric inhibition in the bound substrate. In conjunction with previously reported PCAT1 structures, we offer fresh insight into the ATP-mediated association between PEP and transmembrane domains as a putative mechanism to optimize peptide cleavage by regulating the width and flexibility of the enzyme active site.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

PPIDomainMiner: Inferring domain-domain interactions from multiple sources of protein-protein interactions

Many biological processes are mediated by protein-protein interactions (PPIs). Because protein domains are the building blocks of proteins, PPIs likely rely on domain-domain interactions (DDIs). Several attempts exist to infer DDIs from PPI networks but the produced datasets are heterogeneous and sometimes not accessible, while the PPI interactome data keeps growing.

We describe a new computational approach called “PPIDM” (Protein-Protein Interactions Domain Miner) for inferring DDIs using multiple sources of PPIs. The approach is an extension of our previously described “CODAC” (Computational Discovery of Direct Associations using Common neighbors) method for inferring new edges in a tripartite graph. The PPIDM method has been applied to seven widely used PPI resources, using as “Gold-Standard” a set of DDIs extracted from 3D structural databases. Overall, PPIDM has produced a dataset of 84,552 non-redundant DDIs. Statistical significance (p-value) is calculated for each source of PPI and used to classify the PPIDM DDIs in Gold (9,175 DDIs), Silver (24,934 DDIs) and Bronze (50,443 DDIs) categories. Dataset comparison reveals that PPIDM has inferred from the 2017 releases of PPI sources about 46% of the DDIs present in the 2020 release of the 3did database, not counting the DDIs present in the Gold-Standard. The PPIDM dataset contains 10,229 DDIs that are consistent with more than 13,300 PPIs extracted from the IMEx database, and nearly 23,300 DDIs (27.5%) that are consistent with more than 214,000 human PPIs extracted from the STRING database. Examples of newly inferred DDIs covering more than 10 PPIs in the IMEx database are provided.

Further exploitation of the PPIDM DDI reservoir includes the inventory of possible partners of a protein of interest and characterization of protein interactions at the domain level in combination with other methods. The result is publicly available at http://ppidm.loria.fr/.

We revisit at a large scale the question of inferring DDIs from PPIs. Compared to previous studies, we take a unified approach accross multiple sources of PPIs. This approach is a method for inferring new edges in a tripartite graph setting and can be compared to link prediction approaches in knowledge graphs. Aggregation of several sources is performed using an optimized weighted average of the individual scores calculated in each source. A huge dataset of over 84K DDIs is produced which far exceeds the previous datasets. We show that a significant portion of the PPIDM dataset covers a large number of PPIs from curated (IMEx) or non curated (STRING) databases. Such a reservoir of DDIs deserves further exploration and can be combined with high-throughput methods such as cross-linking mass spectrometry to identify plausible protein partners of proteins of interest.

TopDomain: Exhaustive Protein Domain Boundary Metaprediction Combining Multisource Information and Deep Learning

Protein domains are independent, functional, and stable structural units of proteins. Accurate protein domain boundary prediction plays an important role in understanding protein structure and evolution, as well as for protein structure prediction. Current domain boundary prediction methods differ in terms of boundary definition, methodology, and training databases resulting in disparate performance for different proteins. We developed TopDomain, an exhaustive metapredictor, that uses deep neural networks to combine multisource information from sequence- and homology-based features of over 50 primary predictors. For this purpose, we developed a new domain boundary data set termed the TopDomain data set, in which the true annotations are informed by SCOPe annotations, structural domain parsers, human inspection, and deep learning. We benchmark TopDomain against 2484 targets with 3354 boundaries from the TopDomain test set and achieve F1 scores of 78.4% and 73.8% for multidomain boundary prediction within ±20 residues and ±10 residues of the true boundary, respectively. When examined on targets from CASP11-13 competitions, TopDomain achieves F1 scores of 47.5% and 42.8% for multidomain proteins. TopDomain significantly outperforms 15 widely used, state-of-the-art ab initio and homology-based domain boundary predictors. Finally, we implemented TopDomain_TMC, which accurately predicts whether domain parsing is necessary for the target protein.

Computational-Driven Epitope Verification and Affinity Maturation of TLR4-Targeting Antibodies

Toll-like receptor (TLR) signaling plays a critical role in the induction and progression of autoimmune diseases such as rheumatoid arthritis, systemic lupus erythematous, experimental autoimmune encephalitis, type 1 diabetes mellitus and neurodegenerative diseases. Deciphering antigen recognition by antibodies provides insights and defines the mechanism of action into the progression of immune responses. Multiple strategies, including phage display and hybridoma technologies, have been used to enhance the affinity of antibodies for their respective epitopes. Here, we investigate the TLR4 antibody-binding epitope by computational-driven approach. We demonstrate that three important residues, i.e., Y328, N329, and K349 of TLR4 antibody binding epitope identified upon in silico mutagenesis, affect not only the interaction and binding affinity of antibody but also influence the structural integrity of TLR4. Furthermore, we predict a novel epitope at the TLR4-MD2 interface which can be targeted and explored for therapeutic antibodies and small molecules. This technique provides an in-depth insight into antibody-antigen interactions at the resolution and will be beneficial for the development of new monoclonal antibodies. Computational techniques, if coupled with experimental methods, will shorten the duration of rational design and development of antibody therapeutics.

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

Identification of Ku70 Domain-Specific Interactors Using BioID2

Since its inception, proximity-dependent biotin identification (BioID), an in vivo biochemical screening method to identify proximal protein interactors, has seen extensive developments. Improvements and variants of the original BioID technique are being reported regularly, each expanding upon the existing potential of the original technique. While this is advancing our capabilities to study protein interactions under different contexts, we have yet to explore the full potential of the existing BioID variants already at our disposal. Here, we used BioID2 in an innovative manner to identify and map domain-specific protein interactions for the human Ku70 protein. Four HEK293 cell lines were created, each stably expressing various BioID2-tagged Ku70 segments designed to collectively identify factors that interact with different regions of Ku70. Historically, although many interactions have been mapped to the C-terminus of the Ku70 protein, few have been mapped to the N-terminal von Willebrand A-like domain, a canonical protein-binding domain ideally situated as a site for protein interaction. Using this segmented approach, we were able to identify domain-specific interactors as well as evaluate advantages and drawbacks of the BioID2 technique. Our study identifies several potential new Ku70 interactors and validates RNF113A and Spindly as proteins that contact or co-localize with Ku in a Ku70 vWA domain-specific manner.

Biophysical and Dynamic Characterization of Fine-Tuned Binding of the Human Respiratory Syncytial Virus M2-1 Core Domain to Long RNAs

The human respiratory syncytial virus (hRSV) M2-1 protein functions as a processivity and antitermination factor of the viral polymerase complex. Here, the first evidence that the hRSV M2-1 core domain (cdM2-1) alone has an unfolding activity for long RNAs is presented and the biophysical and dynamic characterization of the cdM2-1/RNA complex is provided. The main contact region of cdM2-1 with RNA was the α1-α2-α5-α6 helix bundle, which suffered local conformational changes and promoted the RNA unfolding activity. This activity may be triggered by base-pairing recognition. RNA molecules wrap around the whole cdM2-1, protruding their termini over the domain. The α2-α3 and α3-α4 loops of cdM2-1 were marked by an increase in picosecond internal motions upon RNA binding, even though they are not directly involved in the interaction. The results revealed that the cdM2-1/RNA complex originates from a fine-tuned binding, contributing to the unraveling interaction aspects necessary for M2-1 activity.IMPORTANCE The main outcome is the molecular description of the fine-tuned binding of the cdM2-1/RNA complex and the provision of evidence that the domain alone has unfolding activity for long RNAs. This binding mode is essential in the understanding of the function in the full-length protein. Human respiratory syncytial virus (hRSV), an orthopneumovirus, stands out for the unique role of its M2-1 protein as a transcriptional antitermination factor able to increase RNA polymerase processivity.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.

Functional Disassociation Between the Protein Domains of MSMEG_4305 of Mycolicibacterium smegmatis (Mycobacterium smegmatis) in vivo

MSMEG_4305 is a two-domain protein of Mycolicibacterium smegmatis (Mycobacterium smegmatis) (Mycolicibacterium smegmatis). The N-terminal domain of MSMEG_4305 encodes an RNase H type I. The C-terminal domain is a presumed CobC, predicted to be involved in the aerobic synthesis of vitamin B12. Both domains reach their maximum at distinct pH, approximately 8.5 and 4.5, respectively. The presence of the CobC domain influenced RNase activity in vitro in homolog Rv2228c. Here, we analyzed the role of MSMEG_4305 in vitamin B12 synthesis and the functional association between both domains in vivo in M. smegmatis. We used knock-out mutant of M. smegmatis, deficient in MSMEG_4305. Whole-cell lysates of the mutants strain contained a lower concentration of vitamin B12, as it determined with immunoenzimatic assay. We observed growth deficits, related to vitamin B12 production, on media containing sulfamethazine and propionate. Removal of the CobC domain of MSMEG_4305 in ΔrnhA background hardly affected the growth rate of M. smegmatis in vivo. The strain carrying truncation showed no fitness deficit in the competitive assay and it did not show increased level of RNA/DNA hybrids in its genome. We show that homologs of MSMEG_4305 are present only in the Actinomycetales phylogenetic branch (according to the old classification system). The domains of MSMEG_4305 homologs accumulate mutations at a different rate, while the linker region is highly variable. We conclude that MSMEG_4305 is a multidomain protein that most probably was fixed in the phylogenetic tree of life due to genetic drift.

A Change in Domain Cooperativity Drives the Function of Calnuc

With the increasing incidence of neurodegenerative disorders, there is an urgent need to understand the protein folding process. Examining the folding process of multidomain proteins remains a prime challenge, as their complex conformational dynamics make them highly susceptible to misfolding and/or aggregation. The presence of multiple domains in a protein can lead to interaction between the partially folded domains, thereby driving misfolding and/or aggregation. Calnuc is one such multidomain protein for which Ca²⁺ binding plays a pivotal role in governing its structural dynamics and stability and, presumably, in directing its interactions with other proteins. We demonstrate differential structural dynamics between the Ca²⁺-free and Ca²⁺-bound forms of calnuc. In the absence of Ca²⁺, full-length calnuc displays equilibrium structural transitions with four intermediate states, reporting a sum of the behavioral properties of its individual domains. Fragment-based studies illustrate the sequential events of structure adoption proceeding in the following order: EF domain followed by the NT and LZ domains in the apo state. On the other hand, Ca²⁺ binding increases domain cooperativity and enables the protein to fold as a single unit. Single-tryptophan mutant proteins, designed in a domain-dependent manner, confirm an increase in the number of interdomain interactions in the Ca²⁺-bound form as compared to the Ca²⁺-free state of the protein, thereby providing insight into its folding process. The attenuated domain crosstalk in apo-calnuc is likely to influence and regulate its physiologically important intermolecular interactions.

MetaGeneHunt for protein domain annotation in short-read metagenomes

The annotation of short-reads metagenomes is an essential process to understand the functional potential of sequenced microbial communities. Annotation techniques based solely on the identification of local matches tend to confound local sequence similarity and overall protein homology and thus don't mirror the complex multidomain architecture and the shuffling of functional domains in many protein families. Here, we present MetaGeneHunt to identify specific protein domains and to normalize the hit-counts based on the domain length. We used MetaGeneHunt to investigate the potential for carbohydrate processing in the mouse gastrointestinal tract. We sampled, sequenced, and analyzed the microbial communities associated with the bolus in the stomach, intestine, cecum, and colon of five captive mice. Focusing on Glycoside Hydrolases (GHs) we found that, across samples, 58.3% of the 4,726,023 short-read sequences matching with a GH domain-containing protein were located outside the domain of interest. Next, before comparing the samples, the counts of localized hits matching the domains of interest were normalized to account for the corresponding domain length. Microbial communities in the intestine and cecum displayed characteristic GH profiles matching distinct microbial assemblages. Conversely, the stomach and colon were associated with structurally and functionally more diverse and variable microbial communities. Across samples, despite fluctuations, changes in the functional potential for carbohydrate processing correlated with changes in community composition. Overall MetaGeneHunt is a new way to quickly and precisely identify discrete protein domains in sequenced metagenomes processed with MG-RAST. In addition, using the sister program "GeneHunt" to create custom Reference Annotation Table, MetaGeneHunt provides an unprecedented way to (re)investigate the precise distribution of any protein domain in short-reads metagenomes.

Ubiquitin Signaling Regulates RNA Biogenesis, Processing, and Metabolism

The fate of eukaryotic proteins, from their synthesis to destruction, is supervised by the ubiquitin-proteasome system (UPS). The UPS is the primary pathway responsible for selective proteolysis of intracellular proteins, which is guided by covalent attachment of ubiquitin to target proteins by E1 (activating), E2 (conjugating), and E3 (ligating) enzymes in a process known as ubiquitylation. The UPS can also regulate protein synthesis by influencing multiple steps of RNA (ribonucleic acid) metabolism. Here, recent publications concerning the interplay between the UPS and different types of RNA are reviewed. This interplay mainly involves specific RNA-binding E3 ligases that link RNA-dependent processes with protein ubiquitylation. The emerging understanding of their modes of RNA binding, their RNA targets, and their molecular and cellular functions are primarily focused on. It is discussed how the UPS adapted to interact with different types of RNA and how RNA molecules influence the ubiquitin signaling components.

Energetic dependencies dictate folding mechanism in a complex protein

Large proteins composed of multiple domains are abundant in all proteomes, but their folding and structural dynamics remain poorly understood. Using single-molecule force spectroscopy, we have defined how stabilizing interfaces among the domains of elongation factor G (EF-G) shape its folding pathway. Contrary to the expectation that multidomain proteins fold sequentially as they emerge from the ribosome, we find that folding cannot be completed until the full protein has been synthesized. This posttranslational folding mechanism results in a propensity for misfolding. It is dictated by an energetic coupling among domains that enables conformational flexibility crucial for EF-G function. EF-G thus provides an example of how distinct biological ends—robust folding and functionally important flexibility—come into conflict during protein biogenesis.

Large proteins with multiple domains are thought to fold cotranslationally to minimize interdomain misfolding. Once folded, domains interact with each other through the formation of extensive interfaces that are important for protein stability and function. However, multidomain protein folding and the energetics of domain interactions remain poorly understood. In elongation factor G (EF-G), a highly conserved protein composed of 5 domains, the 2 N-terminal domains form a stably structured unit cotranslationally. Using single-molecule optical tweezers, we have defined the steps leading to fully folded EF-G. We find that the central domain III of EF-G is highly dynamic and does not fold upon emerging from the ribosome. Surprisingly, a large interface with the N-terminal domains does not contribute to the stability of domain III. Instead, it requires interactions with its folded C-terminal neighbors to be stably structured. Because of the directionality of protein synthesis, this energetic dependency of domain III on its C-terminal neighbors disrupts cotranslational folding and imposes a posttranslational mechanism on the folding of the C-terminal part of EF-G. As a consequence, unfolded domains accumulate during synthesis, leading to the extensive population of misfolded species that interfere with productive folding. Domain III flexibility enables large-scale conformational transitions that are part of the EF-G functional cycle during ribosome translocation. Our results suggest that energetic tuning of domain stabilities, which is likely crucial for EF-G function, complicates the folding of this large multidomain protein.

Unfolding of polymers tethered to viscoelastic substrates

The problem of globular polymer unfolding under applied force is a widely-studied fundamental topic in biological and chemical physics, with important applications in cell biology. Much of the existing theoretical and experimental literature focuses on the case where force is applied while fixing the opposite end of the polymer chain in space. However, in a realistic biological microenvironment, forces will be applied against viscoelastic references, and the deformation of the folded polymer chain will be combined with the deformation of viscoelastic substrate. In this paper, we consider several simple viscoelastic models for the substrate, and show that its relaxation properties determine the unfolding kinetics. In particular, for low pulling forces, substrates with longer relaxation times cause lower unfolding rates for the pulled polymer chain, whereas for high forces, those substrates with longer relaxation times instead produce higher unfolding rates.

Structural studies of the Hsp70/Hsp90 organizing protein of Plasmodium falciparum and its modulation of Hsp70 and Hsp90 ATPase activities

HOP is a cochaperone belonging to the foldosome, a system formed by the cytoplasmic Hsp70 and Hsp90 chaperones. HOP acts as an adapter protein capable of transferring client proteins from the first to the second molecular chaperone. HOP is a modular protein that regulates the ATPase activity of Hsp70 and Hsp90 to perform its function. To obtain more detailed information on the structure and function of this protein, we produced the recombinant HOP of Plasmodium falciparum (PfHOP). The protein was obtained in a folded form, with a high content of α-helix secondary structure. Unfolding experiments showed that PfHOP unfolds through two transitions, suggesting the presence of at least two domains with different stabilities. In addition, PfHOP primarily behaved as an elongated dimer in equilibrium with the monomer. Small-angle X-ray scattering data corroborated this interpretation and led to the reconstruction of a PfHOP ab initio model as a dimer. Finally, the PfHOP protein was able to inhibit and to stimulate the ATPase activity of the recombinant Hsp90 and Hsp70-1, respectively, of P. falciparum. Our results deepened the knowledge of the structure and function of PfHOP and further clarified its participation in the P. falciparum foldosome.

Unraveling the structural landscape of intra-chain domain interfaces: Implication in the evolution of domain-domain interactions

Intra-chain domain interactions are known to play a significant role in the function and stability of multidomain proteins. These interactions are mediated through a physical interaction at domain-domain interfaces (DDIs). With a motivation to understand evolution of interfaces, we have investigated similarities among DDIs. Even though interfaces of protein-protein interactions (PPIs) have been previously studied by structurally aligning interfaces, similar analyses have not yet been performed on DDIs of either multidomain proteins or PPIs. For studying the structural landscape of DDIs, we have used iAlign to structurally align intra-chain domain interfaces of domains. The interface alignment of spatially constrained domains (due to inter-domain linkers) showed that ~88% of these could identify a structural matching interface having similar C-alpha geometry and contact pattern despite that aligned domain pairs are not structurally related. Moreover, the mean interface similarity score (IS-score) is 0.307, which is higher compared to the average random IS-score (0.207) suggesting domain interfaces are not random. The structural space of DDIs is highly connected as ~84% of all possible directed edges among interfaces are found to have at most path length of 8 when 0.26 is IS-score threshold. At this threshold, ~83% of interfaces form the largest strongly connected component. Thus, suggesting that structural space of intra-chain domain interfaces is degenerate and highly connected, as has been found in PPI interfaces. Interestingly, searching for structural neighbors of inter-chain interfaces among intra-chain interfaces showed that ~86% could find a statistically significant match to intra-chain interface with a mean IS-score of 0.311. This implies that domain interfaces are degenerate whether formed within a protein or between proteins. The interface degeneracy is most likely due to limited possible ways of packing secondary structures. In principle, interface similarities can be exploited to accurately model domain interfaces in structure prediction of multidomain proteins.

Fold combinations in multi-domain proteins

Domain-domain interactions in multi-domain proteins play an important role in the combined function of individual domains for the overall biological activity of the protein. The functions of the tethered domains are often coupled and hence, limited numbers of domain architectures with defined folds are known in nature. Therefore, it is of interest to document the available fold-fold combinations and their preference in multi-domain proteins. Hence, we analyzed all multi-domain proteins with known structures in the protein databank and observed that only about 860 fold-fold combinations are present among them. Analyses of multi-domain proteins represented in sequence database result in recognition of 29,860 fold-fold combinations and it accounts for only 2.8% of the theoretically possible 1,036,080 (1439C2) fold-fold combinations. The observed preference for fold-fold combinations in multi-domain proteins is interesting in the context of multiple functions through structural adaptation by gene fusion.

Directed evolution of the metalloproteinase inhibitor TIMP-1 reveals that its N- and C-terminal domains cooperate in matrix metalloproteinase recognition

Tissue inhibitors of metalloproteinases (TIMPs) are natural inhibitors of matrix metalloproteinases (MMPs), enzymes that contribute to cancer and many inflammatory and degenerative diseases. The TIMP N-terminal domain binds and inhibits an MMP catalytic domain, but the role of the TIMP C-terminal domain in MMP inhibition is poorly understood. Here, we employed yeast surface display for directed evolution of full-length human TIMP-1 to develop MMP-3-targeting ultrabinders. By simultaneously incorporating diversity into both domains, we identified TIMP-1 variants that were up to 10-fold improved in binding MMP-3 compared with WT TIMP-1, with inhibition constants (K_i ) in the low picomolar range. Analysis of individual and paired mutations from the selected TIMP-1 variants revealed cooperative effects between distant residues located on the N- and C-terminal TIMP domains, positioned on opposite sides of the interaction interface with MMP-3. Crystal structures of MMP-3 complexes with TIMP-1 variants revealed conformational changes in TIMP-1 near the cooperative mutation sites. Affinity was strengthened by cinching of a reciprocal "tyrosine clasp" formed between the N-terminal domain of TIMP-1 and proximal MMP-3 interface and by changes in secondary structure within the TIMP-1 C-terminal domain that stabilize interdomain interactions and improve complementarity to MMP-3. Our protein engineering and structural studies provide critical insight into the cooperative function of TIMP domains and the significance of peripheral TIMP epitopes in MMP recognition. Our findings suggest new strategies to engineer TIMP proteins for therapeutic applications, and our directed evolution approach may also enable exploration of functional domain interactions in other protein systems.

Global computational mutagenesis of domain structures associated with inherited eye disease

Multidomain proteins account for 70% of the eukaryotic proteome. In genetic disease, multidomain proteins are often affected by numerous mutations, but the effects of these mutations on protein stability and their roles in genetic disease are not well understood. Here, we analyzed protein globular domains to understand how genetic mutations affect the stability of multidomain proteins in inherited disease. In total, 291 domain atomic structures from nine multidomain proteins were modeled by homology, equilibrated using molecular dynamics in water, and subjected to global computational mutagenesis. The domains were separated into 7 groups based on protein fold homology. Mutation propensities within each group of domains were then averaged to select residues critical for domain fold stability. The consensus derived from the sequence alignment shows that the critical residues determined by global mutagenesis are conserved within each group. From this analysis, we concluded that 80% of known disease-related genetic variants are associated with critical residues and are expected to have significant destabilizing effects on domain structure. Our work provides an in silico quantification of protein stability and could help to analyze the complex relationship among missense mutations, multidomain protein stability, and disease phenotypes in inherited eye disease.

Same but not alike: Structure, flexibility and energetics of domains in multi-domain proteins are influenced by the presence of other domains

The majority of the proteins encoded in the genomes of eukaryotes contain more than one domain. Reasons for high prevalence of multi-domain proteins in various organisms have been attributed to higher stability and functional and folding advantages over single-domain proteins. Despite these advantages, many proteins are composed of only one domain while their homologous domains are part of multi-domain proteins. In the study presented here, differences in the properties of protein domains in single-domain and multi-domain systems and their influence on functions are discussed. We studied 20 pairs of identical protein domains, which were crystallized in two forms (a) tethered to other proteins domains and (b) tethered to fewer protein domains than (a) or not tethered to any protein domain. Results suggest that tethering of domains in multi-domain proteins influences the structural, dynamic and energetic properties of the constituent protein domains. 50% of the protein domain pairs show significant structural deviations while 90% of the protein domain pairs show differences in dynamics and 12% of the residues show differences in the energetics. To gain further insights on the influence of tethering on the function of the domains, 4 pairs of homologous protein domains, where one of them is a full-length single-domain protein and the other protein domain is a part of a multi-domain protein, were studied. Analyses showed that identical and structurally equivalent functional residues show differential dynamics in homologous protein domains; though comparable dynamics between in-silico generated chimera protein and multi-domain proteins were observed. From these observations, the differences observed in the functions of homologous proteins could be attributed to the presence of tethered domain. Overall, we conclude that tethered domains in multi-domain proteins not only provide stability or folding advantages but also influence pathways resulting in differences in function or regulatory properties.

High prevalence of multi-domain proteins in proteomes has been attributed to higher stability and functional and folding advantages of the multi-domain proteins. Influence of tethering of domains on the overall properties of proteins has been well studied but its influence on the properties of the constituent domains is largely unaddressed. Here, we investigate the influence of tethering of domains in multi-domain proteins on the structural, dynamics and energetics properties of the constituent domains and its implications on the functions of proteins. To this end, comparative analyses were carried out for identical protein domains crystallized in tethered and untethered forms. Also, comparative analyses of single-domain proteins and their homologous multi-domain proteins were performed. The analyses suggest that tethering influences the structural, dynamic and energetic properties of constituent protein domains. Our observations hint at regulation of protein domains by tethered domains in multi-domain systems, which may manifest at the differential function observed between single-domain and homologous multi-domain proteins.

FoldX as Protein Engineering Tool: Better Than Random Based Approaches?

Improving protein stability is an important goal for basic research as well as for clinical and industrial applications but no commonly accepted and widely used strategy for efficient engineering is known. Beside random approaches like error prone PCR or physical techniques to stabilize proteins, e.g. by immobilization, in silico approaches are gaining more attention to apply target-oriented mutagenesis. In this review different algorithms for the prediction of beneficial mutation sites to enhance protein stability are summarized and the advantages and disadvantages of FoldX are highlighted. The question whether the prediction of mutation sites by the algorithm FoldX is more accurate than random based approaches is addressed.

Cooperative mechanosensitivity and allostery of focal adhesion clusters

We analyse the role of cooperative interaction between neighbouring adhesion-mechanosensor complexes by constructing an Ising-like Hamiltonian describing the free energy of cell adhesion on a substrate as a lattice of 3-state mechanosensing sites involving focal adhesion kinase (FAK). We use a Monte Carlo stochastic algorithm to find equilibrium configurations of these mechanosensors in two representative geometries: on a 1D ring representing the rim of a cell on a flat surface, and a 2D bounded surface representing the whole area of cell contact with a flat surface. The level of FAK activation depends on the pulling force applied to the individual FAK-integrin via actin-myosin contractile networks, and the details of the coupling between individual sensors in a cluster. Strong coupling is shown to make the FAK sensors experience a sharp on-off behaviour in their activation, while at low coupling the activation/autoinhibition transition occurs over a broad range of pulling force. We find that the activation/autoinhibition transition of FAK in the 2D system with strong coupling occurs with a hysteresis, the width of which depends on the rate of change of force. The effect of introducing a regulating protein (such as Src) in a limited quantity to control FAK activation is explored, and visualizations of clustering in both topologies are presented. In particular the results on the bounded 2D surface indicate that clustering of active FAK occurs preferentially at the boundary, in agreement with experimental observations of focal adhesions in cells.

Protein Assembly and Building Blocks: Beyond the Limits of the LEGO Brick Metaphor

Proteins, like other biomolecules, have a modular and hierarchical structure. Various building blocks are used to construct proteins of high structural complexity and diverse functionality. In multidomain proteins, for example, domains are fused to each other in different combinations to achieve different functions. Although the LEGO brick metaphor is justified as a means of simplifying the complexity of three-dimensional protein structures, several fundamental properties (such as allostery or the induced-fit mechanism) make deviation from it necessary to respect the plasticity, softness, and cross-talk that are essential to protein function. In this work, we illustrate recently reported protein behavior in multidomain proteins that deviates from the LEGO brick analogy. While earlier studies showed that a protein domain is often unaffected by being fused to another domain or becomes more stable following the formation of a new interface between the tethered domains, destabilization due to tethering has been reported for several systems. We illustrate that tethering may sometimes result in a multidomain protein behaving as "less than the sum of its parts". We survey these cases for which structure additivity does not guarantee thermodynamic additivity. Protein destabilization due to fusion to other domains may be linked in some cases to biological function and should be taken into account when designing large assemblies.

Focal Adhesion Kinase: The Reversible Molecular Mechanosensor

Sensors are the first element of the pathways that control the response of cells to their environment. Protein complexes that produce or enable a chemical signal in response to a mechanical stimulus are called "mechanosensors". In this work, we develop a theoretical model describing the physical mechanism of a reversible single-molecule stiffness sensor. Although this has the potential for general application, here we apply the model to focal adhesion kinase, which initiates the chemical signal in its active phosphorylated conformation, but can spontaneously return to its closed folded conformation. We find how the rates of conformation changes depend on the substrate stiffness and the pulling force applied from the cell cytoskeleton. We find the sensor is homeostatic, spontaneously self-adjusting to reach a state where its range of maximum sensitivity matches the substrate stiffness. The results compare well with the phenotype observations of cells on different substrates.

Multifactorial level of extremostability of proteins: can they be exploited for protein engineering?

Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.

Concatemers of Outer Membrane Protein A Take Detours in the Folding Landscape

Outer membrane protein A (OmpA) is the most abundant protein in the outer membrane of Escherichia coli. The N-terminal domain forms an eight-stranded membrane-embedded β-barrel that is widely used as a model protein for in vitro folding into the membrane and into surfactant micelles. Under conditions that include a low surfactant concentration, OmpA can form stable higher-order structures by intermolecular association. Other β-barrel membrane proteins also associate to form noncovalently linked trimers in vivo. This inspired us to test how topological constraints imposed by intramolecular links between individual OmpA molecules affect this process. Here we report on the properties of concatemers consisting of two and three copies of the transmembrane part of OmpA. Both concatemers could be folded to a native state in surfactant micelles according to spectroscopy and electrophoretic band shifts. This native state had the same thermodynamic stability against chemical denaturation as the original OmpA. Above 1.5 M GdmCl, concatemerization increased both refolding and unfolding rates, which we attribute to entropic effects. However, below 1.5 M GdmCl, folding kinetics were 2-3 orders of magnitude slower and more complex, involving a greater degree of parallel folding steps and species that could be classified as off-pathway. Only OmpA2 could quantitatively be folded into vesicles (though to an extent lower than that of OmpA), while OmpA3 formed three species with different levels of folding. Thus, close spatial and sequential proximity of OmpA domains on the same polypeptide chain have a strong tendency to trap the protein in different misfolded states.

Ubiquitylation Directly Induces Fold Destabilization of Proteins

Ubiquitin is a common post-translational modifier and its conjugation is a key signal for proteolysis by the proteasome. Because the molecular mass of ubiquitin is larger than that of other modifiers such as phosphate, acetyl, or methyl groups, ubiquitylation not only influences biochemical signaling, but also may exert physical effects on its substrate proteins by increasing molecular volume and altering shape anisotropy. Here we show that ubiquitylation destabilizes the fold of two proteins, FKBP12 and FABP4, and that elongation of the conjugated ubiquitin chains further enhances this destabilization effect. Moreover, NMR relaxation analysis shows that ubiquitylation induces characteristic structural fluctuations in the backbone of both proteins. These results suggest that the ubiquitylation-driven structural fluctuations lead to fold destabilization of its substrate proteins. Thus, physical destabilization by ubiquitylation may facilitate protein degradation by the proteasome.

Interdomain salt-bridges in the Ebola virus protein VP40 and their role in domain association and plasma membrane localization

The Ebola virus protein VP40 is a transformer protein that possesses an extraordinary ability to accomplish multiple functions by transforming into various oligomeric conformations. The disengagement of the C-terminal domain (CTD) from the N-terminal domain (NTD) is a crucial step in the conformational transformations of VP40 from the dimeric form to the hexameric form or octameric ring structure. Here, we use various molecular dynamics (MD) simulations to investigate the dynamics of the VP40 protein and the roles of interdomain interactions that are important for the domain-domain association and dissociation, and report on experimental results of the behavior of mutant variants of VP40. The MD studies find that various salt-bridge interactions modulate the VP40 domain dynamics by providing conformational specificity through interdomain interactions. The MD simulations reveal a novel salt-bridge between D45-K326 when the CTD participates in a latch-like interaction with the NTD. The D45-K326 salt-bridge interaction is proposed to help domain-domain association, whereas the E76-K291 interaction is important for stabilizing the closed-form structure. The effects of the removal of important VP40 salt-bridges on plasma membrane (PM) localization, VP40 oligomerization, and virus like particle (VLP) budding assays were investigated experimentally by live cell imaging using an EGFP-tagged VP40 system. It is found that the mutations K291E and D45K show enhanced PM localization but D45K significantly reduced VLP formation.

Natural diversity of cellulases, xylanases, and chitinases in bacteria

BACKGROUND

Glycoside hydrolases (GH) targeting cellulose, xylan, and chitin are common in the bacterial genomes that have been sequenced. Little is known, however, about the architecture of multi-domain and multi-activity glycoside hydrolases. In these enzymes, combined catalytic domains act synergistically and thus display overall improved catalytic efficiency, making these proteins of high interest for the biofuel technology industry.

RESULTS

Here, we identify the domain organization in 40,946 proteins targeting cellulose, xylan, and chitin derived from 11,953 sequenced bacterial genomes. These bacteria are known to be capable, or to have the potential, to degrade polysaccharides, or are newly identified potential degraders (e.g., Actinospica, Hamadaea, Cystobacter, and Microbispora). Most of the proteins we identified contain a single catalytic domain that is frequently associated with an accessory non-catalytic domain. Regarding multi-domain proteins, we found that many bacterial strains have unique GH protein architectures and that the overall protein organization is not conserved across most genera. We identified 217 multi-activity proteins with at least two GH domains for cellulose, xylan, and chitin. Of these proteins, 211 have GH domains targeting similar or associated substrates (i.e., cellulose and xylan), whereas only six proteins target both cellulose and chitin. Fifty-two percent of multi-activity GHs are hetero-GHs. Finally, GH6, -10, -44 and -48 domains were mostly C-terminal; GH9, -11, -12, and -18 were mostly N-terminal; and GH5 domains were either N- or C-terminal.

CONCLUSION

We identified 40,946 multi-domain/multi-activity proteins targeting cellulase, chitinase, and xylanase in bacterial genomes and proposed new candidate lineages and protein architectures for carbohydrate processing that may play a role in biofuel production.

Interpreting functional effects of coding variants: challenges in proteome-scale prediction, annotation and assessment

Accurate assessment of genetic variation in human DNA sequencing studies remains a nontrivial challenge in clinical genomics and genome informatics. Ascribing functional roles and/or clinical significances to single nucleotide variants identified from a next-generation sequencing study is an important step in genome interpretation. Experimental characterization of all the observed functional variants is yet impractical; thus, the prediction of functional and/or regulatory impacts of the various mutations using in silico approaches is an important step toward the identification of functionally significant or clinically actionable variants. The relationships between genotypes and the expressed phenotypes are multilayered and biologically complex; such relationships present numerous challenges and at the same time offer various opportunities for the design of in silico variant assessment strategies. Over the past decade, many bioinformatics algorithms have been developed to predict functional consequences of single nucleotide variants in the protein coding regions. In this review, we provide an overview of the bioinformatics resources for the prediction, annotation and visualization of coding single nucleotide variants. We discuss the currently available approaches and major challenges from the perspective of protein sequence, structure, function and interactions that require consideration when interpreting the impact of putatively functional variants. We also discuss the relevance of incorporating integrated workflows for predicting the biomedical impact of the functionally important variations encoded in a genome, exome or transcriptome. Finally, we propose a framework to classify variant assessment approaches and strategies for incorporation of variant assessment within electronic health records.

Structural and functional stabilization of protein entities: state-of-the-art

Within the context of biomedicine and pharmaceutical sciences, the issue of (therapeutic) protein stabilization assumes particular relevance. Stabilization of protein and protein-like molecules translates into preservation of both structure and functionality during storage and/or targeting, and such stabilization is mostly attained through establishment of a thermodynamic equilibrium with the (micro)environment. The basic thermodynamic principles that govern protein structural transitions and the interactions of the protein molecule with its (micro)environment are, therefore, tackled in a systematic fashion. Highlights are given to the major classes of (bio)therapeutic molecules, viz. enzymes, recombinant proteins, (macro)peptides, (monoclonal) antibodies and bacteriophages. Modification of the microenvironment of the biomolecule via multipoint covalent attachment onto a solid surface followed by hydrophilic polymer co-immobilization, or physical containment within nanocarriers, are some of the (latest) strategies discussed aiming at full structural and functional stabilization of said biomolecules.

The organization of domains in proteins obeys Menzerath-Altmann's law of language

BACKGROUND

The combination of domains in multidomain proteins enhances their function and structure but lengthens the molecules and increases their cost at cellular level.

METHODS

The dependence of domain length on the number of domains a protein holds was surveyed for a set of 60 proteomes representing free-living organisms from all kingdoms of life. Distributions were fitted using non-linear functions and fitted parameters interpreted with a formulation of decreasing returns.

RESULTS

We find that domain length decreases with increasing number of domains in proteins, following the Menzerath-Altmann (MA) law of language. Highly significant negative correlations exist for the set of proteomes examined. Mathematically, the MA law expresses as a power law relationship that unfolds when molecular persistence P is a function of domain accretion. P holds two terms, one reflecting the matter-energy cost of adding domains and extending their length, the other reflecting how domain length and number impinges on information and biophysics. The pattern of diminishing returns can therefore be explained as a frustrated interplay between the strategies of economy, flexibility and robustness, matching previously observed trade-offs in the domain makeup of proteomes. Proteomes of Archaea, Fungi and to a lesser degree Plants show the largest push towards molecular economy, each at their own economic stratum. Fungi increase domain size in single domain proteins while reinforcing the pattern of diminishing returns. In contrast, Metazoa, and to lesser degrees Protista and Bacteria, relax economy. Metazoa achieves maximum flexibility and robustness by harboring compact molecules and complex domain organization, offering a new functional vocabulary for molecular biology.

CONCLUSIONS

The tendency of parts to decrease their size when systems enlarge is universal for language and music, and now for parts of macromolecules, extending the MA law to natural systems.

Computing the origin and evolution of the ribosome from its structure - Uncovering processes of macromolecular accretion benefiting synthetic biology

Accretion occurs pervasively in nature at widely different timeframes. The process also manifests in the evolution of macromolecules. Here we review recent computational and structural biology studies of evolutionary accretion that make use of the ideographic (historical, retrodictive) and nomothetic (universal, predictive) scientific frameworks. Computational studies uncover explicit timelines of accretion of structural parts in molecular repertoires and molecules. Phylogenetic trees of protein structural domains and proteomes and their molecular functions were built from a genomic census of millions of encoded proteins and associated terminal Gene Ontology terms. Trees reveal a ‘metabolic-first’ origin of proteins, the late development of translation, and a patchwork distribution of proteins in biological networks mediated by molecular recruitment. Similarly, the natural history of ancient RNA molecules inferred from trees of molecular substructures built from a census of molecular features shows patchwork-like accretion patterns. Ideographic analyses of ribosomal history uncover the early appearance of structures supporting mRNA decoding and tRNA translocation, the coevolution of ribosomal proteins and RNA, and a first evolutionary transition that brings ribosomal subunits together into a processive protein biosynthetic complex. Nomothetic structural biology studies of tertiary interactions and ancient insertions in rRNA complement these findings, once concentric layering assumptions are removed. Patterns of coaxial helical stacking reveal a frustrated dynamics of outward and inward ribosomal growth possibly mediated by structural grafting. The early rise of the ribosomal ‘turnstile’ suggests an evolutionary transition in natural biological computation. Results make explicit the need to understand processes of molecular growth and information transfer of macromolecules.