Abundance and Temperature Dependency of Protein-Protein Interaction Revealed by Interface Structure Analysis and Stability Evolution

Application of Machine Learning in the Quantitative Analysis of the Surface Characteristics of Highly Abundant Cytoplasmic Proteins: Toward AI-Based Biomimetics

Proteins in the crowded environment of human cells have often been studied regarding nonspecific interactions, misfolding, and aggregation, which may cause cellular malfunction and disease. Specifically, proteins with high abundance are more susceptible to these issues due to the law of mass action. Therefore, the surfaces of highly abundant cytoplasmic (HAC) proteins directly exposed to the environment can exhibit specific physicochemical, structural, and geometrical characteristics that reduce nonspecific interactions and adapt to the environment. However, the quantitative relationships between the overall surface descriptors still need clarification. Here, we used machine learning to identify HAC proteins using hydrophobicity, charge, roughness, secondary structures, and B-factor from the protein surfaces and quantified the contribution of each descriptor. First, several supervised learning algorithms were compared to solve binary classification problems for the surfaces of HAC and extracellular proteins. Then, logistic regression was used for the feature importance analysis of descriptors considering model performance (80.2% accuracy and 87.6% AUC) and interpretability. The HAC proteins showed positive correlations with negatively and positively charged areas but negative correlations with hydrophobicity, the B-factor, the proportion of beta structures, roughness, and the proportion of disordered regions. Finally, the details of each descriptor could be explained concerning adaptative surface strategies of HAC proteins to regulate nonspecific interactions, protein folding, flexibility, stability, and adsorption. This study presented a novel approach using various surface descriptors to identify HAC proteins and provided quantitative design rules for the surfaces well-suited to human cellular crowded environments.

Underlying features for the enhanced electrostatic strength of the extremophilic malate dehydrogenase interface salt-bridge compared to the mesophilic one

Malate dehydrogenase (MDH) exists in multimeric form in normal and extreme solvent conditions where residues of the interface are involved in specific interactions. The interface salt-bridge (ISB) and its microenvironment (ME) residues may have a crucial role in the stability and specificity of the interface. To gain insight into this, we have analyzed 218 ISBs from 42 interfaces of 15 crystal structures along with their sequences. Comparative analyses demonstrate that the ISB strength is ∼30 times greater in extremophilic cases than that of the normal one. To this end, the interface residue propensity, ISB design and pair selection, and ME-residue's types, i.e., type-I and type-II, are seen to be intrinsically involved. Although Type-I is a common type, Type-II appears to be extremophile-specific, where the net ME-residue count is much lower with an excessive net ME-energy contribution, which seems to be a novel interface compaction strategy. Furthermore, the interface strength can be enhanced by selecting the desired mutant from the net-energy profile of all possible mutations of an unfavorable ME-residue. The study that applies to other similar systems finds applications in protein-protein interaction and protein engineering.Communicated by Ramaswamy H. Sarma.

Simultaneous Improvement in the Thermostability and Catalytic Activity of Epoxidase Lsd18 for the Synthesis of Lasalocid A

Enzymes used in the synthesis of natural products are potent catalysts, capable of efficient and stereoselective chemical transformations. Lsd18 catalyzes two sequential epoxidations during the biosynthesis of lasalocid A, a polyether polyketide natural product. We performed protein engineering on Lsd18 to improve its thermostability and catalytic activity. Utilizing structure-guided methods of FoldX and Rosetta-ddG, we designed 15 mutants of Lsd18. Screening of these mutants using thermal shift assay identified stabilized variants Lsd18-T189M, Lsd18-S195M, and the double mutant Lsd18-T189M-S195M. Trypsin digestion, molecular dynamic simulation, circular dichroism (CD) spectroscopy, and X-ray crystallography provided insights into the molecular basis for the improved enzyme properties. Notably, enhanced hydrophobic interaction within the enzyme core and interaction of the protein with the FAD cofactor appear to be responsible for its better thermostability.

Structure of full-length cobalamin-dependent methionine synthase and cofactor loading captured in crystallo

Cobalamin-dependent methionine synthase (MS) is a key enzyme in methionine and folate one-carbon metabolism. MS is a large multi-domain protein capable of binding and activating three substrates: homocysteine, folate, and S-adenosylmethionine for methylation. Achieving three chemically distinct methylations necessitates significant domain rearrangements to facilitate substrate access to the cobalamin cofactor at the right time. The distinct conformations required for each reaction have eluded structural characterization as its inherently dynamic nature renders structural studies difficult. Here, we use a thermophilic MS homolog (tMS) as a functional MS model. Its exceptional stability enabled characterization of MS in the absence of cobalamin, marking the only studies of a cobalamin-binding protein in its apoenzyme state. More importantly, we report the high-resolution full-length MS structure, ending a multi-decade quest. We also capture cobalamin loading in crystallo, providing structural insights into holoenzyme formation. Our work paves the way for unraveling how MS orchestrates large-scale domain rearrangements crucial for achieving challenging chemistries.

Effect of temperature and nucleotide on the binding of BiP chaperone to a protein substrate

BiP (immunoglobulin heavy-chain binding protein) is a Hsp70 monomeric ATPase motor that plays broad and crucial roles in maintaining proteostasis inside the cell. Structurally, BiP is formed by two domains, a nucleotide-binding domain (NBD) with ATPase activity connected by a flexible hydrophobic linker to the substrate-binding domain. While the ATPase and substrate binding activities of BiP are allosterically coupled, the latter is also dependent on nucleotide binding. Recent structural studies have provided new insights into BiP's allostery; however, the influence of temperature on the coupling between substrate and nucleotide binding to BiP remains unexplored. Here, we study BiP's binding to its substrate at the single molecule level using thermo-regulated optical tweezers which allows us to mechanically unfold the client protein and explore the effect of temperature and different nucleotides on BiP binding. Our results confirm that the affinity of BiP for its protein substrate relies on nucleotide binding, by mainly regulating the binding kinetics between BiP and its substrate. Interestingly, our findings also showed that the apparent affinity of BiP for its protein substrate in the presence of nucleotides remains invariable over a wide range of temperatures, suggesting that BiP may interact with its client proteins with similar affinities even when the temperature is not optimal. Thus, BiP could play a role as a "thermal buffer" in proteostasis.

A novel EDAR variant identified in non-syndromic tooth agenesis: Insights from molecular dynamics

OBJECTIVE

This study aims to investigate a novel pathogenic variant in a Chinese family of non-syndromic tooth agenesis (NSTA) and study the impact of the variant on related protein and pathway.

DESIGN

One NSTA family was collected. Whole exome sequencing and Sanger sequencing were performed on the proband with NSTA and his 5 family members. The pathogenic influence of the mutant is evaluated by bioinformatics analyses including evolutionary conservation analysis and secondary structure prediction. Molecular dynamics (MD) simulations and binding free energy calculations were then performed to explore changes in the tertiary structure and binding ability of the protein.

RESULTS

We found a novel missense ectodysplasin A receptor (EDAR) variant (c .1292 T > G; p.Ile431Arg) in all affected family members. The results of bioinformatics analyses revealed that the EDAR had harmful changes after mutation. MD simulations and the binding free energy calculations results showed that the mutant EDAR protein and EDAR/ectodysplasin-A receptor-associated adapter (EDARADD) complex displayed tertiary structural change, and EDAR possessed a lower affinity to EDARADD after mutation.

CONCLUSIONS

We found a novel EDAR variant (c.1292 T > G; p.Ile431Arg) in one NSTA family, which affects the binding of EDAR and EDARADD.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.