Protein charge and mass contribute to the spatio-temporal dynamics of protein-protein interactions in a minimal proteome

Profiling genetic mutations in the DNA damage repair genes of oral squamous cell carcinoma patients from Pakistan

Herein, we reported mutations in five DNA Damage Repair (DDR) i.e., TP53, ATR, ATM, CHEK1 and CHEK2 involved in OSCC using NG-WES and their analysis using bioinformatics tools. Out of 42 identified mutations, 16.7% are reported for the 1st time. A total of 28 nonsynonymous SNVs are identified. TP53 harbored the highest number of mutations followed by ATM, ATR, CHEK1 and CHEK2. Nine mutations (TP53p.R43H, TP53p.L125Q, TP53p.R116Q, TP53p.C110Y, TP53p.L62F, ATRp.H120Y, ATMp.P1054R, ATMp.D1853V, ATMp.T2934N) were predicted highly pathogenic. SAAFEQ-SEQ predicted destabilizing effects for all mutations, while ISPRED-SEQ identified 09 IS mutations, 07 on TP53, 01 in ATR and 01 in CHEK1 with no IS mutations predicted for ATM and CHEK2. Among the IS mutations, only SNVs were used in MDS simulations. The gyration radius for all IS SNVs was larger for mutant as compared to the wild type indicating perturbed folding behavior of the mutant proteins. Structural deviations across the carbon back bone were noted by RMSD for mutant and wild type. The TP53 IS mutations include TP53p.R116Q, TP53 p.C110Y, TP53p.R43H, TP53p.E214X, TP53p.R210X, TP53 p.C110Afs*5 and TP53 p,S108Ffs*23 whereas ATR and CHEK1 IS mutations consist of ATRp.M1932T and CHEK1p.E76Kfs*21. ConSurf analysis revealed four SNVs with a high conservation score (9) on TP53 and ATM. TP53p.P33R was predominantly associated with moderately differentiated tumors (84.60%), naswar users (86.60%) and positive family history of cancer (91.60%). The TP53p.P33R, ATRp.M211T and CHEK1p.I437V mutations were found recurrently in 21/27 (77.7%), 20/27 (74.04%), and 27/27 (100%) patients, suggesting its potential biomarker applications in local screening.

Computational exploration of SLC14A1 genetic variants through structure modeling, protein-ligand docking, and molecular dynamics simulation

The urea transporter UT-B1, encoded by the SLC14A1 gene, has been hypothesized to be a significant protein whose deficiency and dysfunction contribute to the pathogenesis of bladder cancer and many other diseases. Several studies reported the association of genetic alterations in the SLC14A1 (UT-B1) gene with bladder carcinogenesis, suggesting a need for thorough characterization of the UT-B1 protein's coding and non-coding variants. This study used various computational techniques to investigate the commonly occurring germ-line missense and non-coding SNPs (ncSNPs) of the SLC14A1 gene (UT-B1) for their structural, functional, and molecular implications for disease susceptibility and dysfunctionality. SLC14A1 missense variants, primarily identified from the ENSEMBL genome browser, were screened through twelve functionality prediction tools leading to two variants D280Y (predicted detrimental by maximum tools) and D280N (high global MAF) for rs1058396. Subsequently, the ConSurf and NetSurf tools revealed the D280 residue to be in a variable site and exposed on the protein surface. According to I-Mutant2.0 and MUpro, both variants are predicted to cause a significant effect on protein stability. Analysis of molecular docking anticipated these two variants to decrease the binding affinity of UT-B1 protein for the examined ligands to a significant extent. Molecular dynamics also disclosed the possible destabilization of the UT-B1 protein due to single nucleotide polymorphism compared to wild-type protein which may result in impaired protein function. Furthermore, several non-coding SNPs were estimated to affect transcription factor binding and regulation of SLC14A1 gene expression. Additionally, two ncSNPs were found to affect miRNA-based post-transcriptional regulation by creating new seed regions for miRNA binding. This comprehensive in-silico study of SLC14A1 gene variants may serve as a springboard for future large-scale investigations examining SLC14A1 polymorphisms.

Prediction of the most deleterious non-synonymous SNPs in the human IL1B gene: evidence from bioinformatics analyses

BACKGROUND

Polymorphisms in IL1B play a significant role in depression, multiple inflammatory-associated disorders, and susceptibility to infection. Functional non-synonymous SNPs (nsSNPs) result in changes in the encoded amino acids, potentially leading to structural and functional alterations in the mutant proteins. So far, most genetic studies have concentrated on SNPs located in the IL1B promoter region, without addressing nsSNPs and their association with multifactorial diseases. Therefore, this study aimed to explore the impact of deleterious nsSNPs retrieved from the dbSNP database on the structure and functions of the IL1B protein.

RESULTS

Six web servers (SIFT, PolyPhen-2, PROVEAN, SNPs&GO, PHD-SNP, PANTHER) were used to analyze the impact of 222 missense SNPs on the function and structure of IL1B protein. Five novel nsSNPs (E100K, T240I, S53Y, D128Y, and F228S) were found to be deleterious and had a mutational impact on the structure and function of the IL1B protein. The I-mutant v2.0 and MUPro servers predicted that these mutations decreased the stability of the IL1B protein. Additionally, these five mutations were found to be conserved, underscoring their significance in protein structure and function. Three of them (T240I, D128Y, and F228S) were predicted to be cancer-causing nsSNPs. To analyze the behavior of the mutant structures under physiological conditions, we conducted a 50 ns molecular dynamics simulation using the WebGro online tool. Our findings indicate that the mutant values differ from those of the IL1B wild type in terms of RMSD, RMSF, Rg, SASA, and the number of hydrogen bonds.

CONCLUSIONS

This study provides valuable insights into nsSNPs located in the coding regions of IL1B, which lead to direct deleterious effects on the functional and structural aspects of the IL1B protein. Thus, these nsSNPs could be considered significant candidates in the pathogenesis of disorders caused by IL1B dysfunction, contributing to effective drug discovery and the development of precision medications. Thorough research and wet lab experiments are required to verify our findings. Moreover, bioinformatic tools were found valuable in the prediction of deleterious nsSNPs.

The Rapid Evolution of De Novo Proteins in Structure and Complex

Recent studies in the rice genome-wide have established that de novo genes, evolving from noncoding sequences, enhance protein diversity through a stepwise process. However, the pattern and rate of their evolution in protein structure over time remain unclear. Here, we addressed these issues within a surprisingly short evolutionary timescale (<1 million years for 97% of Oryza de novo genes) with comparative approaches to gene duplicates. We found that de novo genes evolve faster than gene duplicates in the intrinsically disordered regions (such as random coils), secondary structure elements (such as α helix and β strand), hydrophobicity, and molecular recognition features. In de novo proteins, specifically, we observed an 8% to 14% decay in random coils and intrinsically disordered region lengths and a 2.3% to 6.5% increase in structured elements, hydrophobicity, and molecular recognition features, per million years on average. These patterns of structural evolution align with changes in amino acid composition over time as well. We also revealed higher positive charges but smaller molecular weights for de novo proteins than duplicates. Tertiary structure predictions showed that most de novo proteins, though not typically well folded on their own, readily form low-energy and compact complexes with other proteins facilitated by extensive residue contacts and conformational flexibility, suggesting a faster-binding scenario in de novo proteins to promote interaction. These analyses illuminate a rapid evolution of protein structure in de novo genes in rice genomes, originating from noncoding sequences, highlighting their quick transformation into active, protein complex-forming components within a remarkably short evolutionary timeframe.

Fractionation and characterisation of sialylated-mucin glycoprotein from edible birds' nest hydrolysates through anion exchange chromatography

Edible bird's nest (EBN) is made up of sialylated-mucin glycoprotein with various health benefits due to its high antioxidative activity. However, as a macromolecule with distinct charged sialic acid and amino acids, fractions with different charges would have varied physicochemical properties and antioxidant activity, which have not been studied. Therefore, this study aimed to fractionate and purify the enzymatic hydrolysed of cleaned EBN (EBNhc) and EBN by-product (EBNhbyp) through anion exchange chromatography (AEC), and determine their molecular weights, physicochemical properties, and antioxidative activities. Overall, 26 fractionates were collected from enzymatic hydrolysate by AEC, which were classified into 5 fractions. It was found that the positively charged fraction of EBNhc (CF 1) and EBNhbyp (DF 1) showed the significantly highest (p < 0.05) soluble protein contents (22.86 and 18.40 mg/g), total peptide contents (511.13 and 800.47 mg/g) and ferric reducing antioxidant power (17.44 and 6.96 mg/g) among the fractionates. In conclusion, a positively charged fraction (CF 1 and DF 1) showed more desired physicochemical properties and antioxidative activities. This research suggests the potential of AEC fractionation as a technology to purify EBN and produce positively charged EBN fractionates with antioxidative potential that could be applied as food components to provide health benefits.

A computational analysis to evaluate deleterious SNPs of GSK3β, a multifunctional and regulatory protein, for metabolism, wound healing, and migratory processes

This study focused on GSK-3β, a critical serine/threonine kinase with diverse cellular functions. However, there is limited understanding of the impact of non-synonymous single nucleotide polymorphisms (nsSNPs) on its structure and function. Through an exhaustive in-silico investigation 12 harmful nsSNPs were predicted from a pool of 172 acquired from the NCBI dbSNP database using 12 established tools that detects deleterious SNPs. Consistently, these nsSNPs were discovered in locations with high levels of conservation. Notably, the three harmful nsSNPs F67C, A83T, and T138I were situated in the active/binding site of GSK-3β, which may affect the protein's capacity to bind to substrates and other proteins. Molecular dynamics simulations revealed that the F67C and T138I mutants had stable structures, indicating rigidness, whereas the A83T mutant was unstable. Analysis of secondary structures revealed different modifications in all mutant forms, which may affect the stability, functioning, and interactions of the protein. These mutations appear to alter the structural dynamics of GSK-3β, which may have functional ramifications, such as the formation of novel secondary structures and variations in coil-to-helix transitions. In conclusion, this study illuminates the possible structural and functional ramifications of these GSK-3 nsSNPs, revealing how protein compactness, stiffness, and interactions may affect biological activities.

Mutational Landscape and In-Silico Analysis of TP53, PIK3CA, and PTEN in Patients with Breast Cancer from Khyber Pakhtunkhwa

Herein, we report the mutational spectrum of three breast cancer candidate genes (TP53, PIK3CA, and PTEN) using WES for identifying potential biomarkers. The WES data were thoroughly analyzed using SAMtools for variant calling and identification of the mutations. Various bioinformatic tools (SIFT, PolyPhen-2, Mutation Taster, ISPRED-SEQ, SAAFEQ-SEQ, ConSurf, PROCHECK etc.) were used to determine the pathogenicity and nature of the SNVs. Selected interaction site (IS) mutations were visualized in PyMOL after building 3D structures in Swiss-Model. Ramachandran plots were generated by using the PROCHECK server. The selected IS mutations were subjected to molecular dynamic simulation (MDS) studies using Gromacs 4.5. STRING and GeneMANIA were used for the prediction of gene-gene interactions and pathways. Our results revealed that the luminal A molecular subtype of the breast cancer was most common, whereas a high percentage of was Her2 negatives. Moreover, the somatic mutations were more common as compared to the germline mutations in TP53, PIK3CA, and PTEN. 20% of the identified mutations are reported for the first time from Khyber Pakhtunkhwa. In the enrolled cohort, 23 mutations were nonsynonymous SNVs. The frequency of mutations was the highest in PIK3CA, followed by TP53 and PTEN. A total of 13 mutations were found to be highly pathogenic. Four novel mutations were identified on PIK3CA and one each on PTEN and TP53. SAAFEQ-SEQ predicted the destabilizing effect for all mutations. ISPRED-SEQ predicted 9 IS mutations (6 on TP53 and 3 on PIK3CA), whereas no IS mutation was predicted on PTEN. The TP53 IS mutations were TP53R43H, TP53Y73X, TP53K93Q, TP53K93R, TP53D149E, and TP53Q199X; whereas for PIK3CA, the IS mutations were PIK3CAL156V, PIK3CAM610K, and PIK3CAH1047R. Analysis from the ConSurf Web server revealed five SNVs with a highly conserved status (conservation score 9) across TP53 and PTEN. TP53P33R was found predominant in the grade 3 tumors, whereas PTENp.C65S was distributed on ER+, ER-, PR+, PR-, Her2+, and Her2- patients. TP53p.P33R mutation was found to be recurring in the 14/19 (73.6%) patients and, therefore, can be considered as a potential biomarker. Finally, these mutations were studied in the context of their potential association with different hormonal and social factors.

Deleterious single nucleotide polymorphisms (SNPs) of human IFNAR2 gene facilitate COVID-19 severity in patients: a comprehensive in silico approach

In humans, the dimeric receptor complex IFNAR2-IFNAR1 accelerates cellular response triggered by type I interferon (IFN) family proteins in response to viral infection including Coronavirus infection. Studies have revealed the association of the IFNAR2 gene with severe illness in Coronavirus infection and indicated the association of genomic variants, i.e. single nucleotide polymorphisms (SNPs). However, comprehensive analysis of SNPs of the IFNAR2 gene has not been performed in both coding and non-coding region to find the causes of loss of function of IFNAR2 in COVID-19 patients. In this study, we have characterized coding SNPs (nsSNPs) of IFNAR2 gene using different bioinformatics tools and identified deleterious SNPs. We found 9 nsSNPs as pathogenic and disease-causing along with a decrease in protein stability. We employed molecular docking analysis that showed 5 nsSNPs to decrease binding affinity to IFN. Later, MD simulations showed that P136R mutant may destabilize crucial binding with the IFN molecule in response to COVID-19. Thus, P136R is likely to have a high impact on disrupting the structure of the IFNAR2 protein. GTEx portal analysis predicted 14 sQTLs and 5 eQTLs SNPs in lung tissues hampering the post-transcriptional modification (splicing) and altering the expression of the IFNAR2 gene. sQTLs and eQTLs SNPs potentially explain the reduced IFNAR2 production leading to severe diseases. These mutants in the coding and non-coding region of the IFNAR2 gene can help to recognize severe illness due to COVID 19 and consequently assist to develop an effective drug against the infection.Communicated by Ramaswamy H. Sarma.

Comprehensive Characterization of the Coding and Non-Coding Single Nucleotide Polymorphisms in the Tumor Protein p63 (TP63) Gene Using In Silico Tools

Single nucleotide polymorphisms (SNPs) help to understand the phenotypic variations in humans. Genome-wide association studies (GWAS) have identified SNPs located in the tumor protein 63 (TP63) locus to be associated with the genetic susceptibility of cancers. However, there is a lack of in-depth characterization of the structural and functional impacts of the SNPs located at the TP63 gene. The current study was designed for the comprehensive characterization of the coding and non-coding SNPs in the human TP63 gene for their functional and structural significance. The functional and structural effects of the SNPs were investigated using a wide variety of computational tools and approaches, including molecular dynamics (MD) simulation. The deleterious impact of eight nonsynonymous SNPs (nsSNPs) affecting protein stability, structure, and functions was measured by using 13 bioinformatics tools. These eight nsSNPs are in highly conserved positions in protein and were predicted to decrease protein stability and have a deleterious impact on the TP63 protein function. Molecular docking analysis showed five nsSNPs to reduce the binding affinity of TP63 protein to DNA with significant results for three SNPs (R319H, G349E, and C347F). Further, MD simulations revealed the possible disruption of TP63 and DNA binding, hampering the essential protein function. PolymiRTS study found five non-coding SNPs in miRNA binding sites, and the GTEx portal recognized five eQTLs SNPs in single tissue of the lung, heart (LV), and cerebral hemisphere (brain). Characterized nsSNPs and non-coding SNPs will help researchers to focus on TP63 gene loci and ascertain their association with certain diseases.

Pathogenic nsSNPs that increase the risks of cancers among the Orang Asli and Malays

Single-nucleotide polymorphisms (SNPs) are the most common genetic variations for various complex human diseases, including cancers. Genome-wide association studies (GWAS) have identified numerous SNPs that increase cancer risks, such as breast cancer, colorectal cancer, and leukemia. These SNPs were cataloged for scientific use. However, GWAS are often conducted on certain populations in which the Orang Asli and Malays were not included. Therefore, we have developed a bioinformatic pipeline to mine the whole-genome sequence databases of the Orang Asli and Malays to determine the presence of pathogenic SNPs that might increase the risks of cancers among them. Five different in silico tools, SIFT, PROVEAN, Poly-Phen-2, Condel, and PANTHER, were used to predict and assess the functional impacts of the SNPs. Out of the 80 cancer-related nsSNPs from the GWAS dataset, 52 nsSNPs were found among the Orang Asli and Malays. They were further analyzed using the bioinformatic pipeline to identify the pathogenic variants. Three nsSNPs; rs1126809 (TYR), rs10936600 (LRRC34), and rs757978 (FARP2), were found as the most damaging cancer pathogenic variants. These mutations alter the protein interface and change the allosteric sites of the respective proteins. As TYR, LRRC34, and FARP2 genes play important roles in numerous cellular processes such as cell proliferation, differentiation, growth, and cell survival; therefore, any impairment on the protein function could be involved in the development of cancer. rs1126809, rs10936600, and rs757978 are the important pathogenic variants that increase the risks of cancers among the Orang Asli and Malays. The roles and impacts of these variants in cancers will require further investigations using in vitro cancer models.

Charge-driven condensation of RNA and proteins suggests broad role of phase separation in cytoplasmic environments

Phase separation processes are increasingly being recognized as important organizing mechanisms of biological macromolecules in cellular environments. Well-established drivers of phase separation are multi-valency and intrinsic disorder. Here, we show that globular macromolecules may condense simply based on electrostatic complementarity. More specifically, phase separation of mixtures between RNA and positively charged proteins is described from a combination of multiscale computer simulations with microscopy and spectroscopy experiments. Phase diagrams were mapped out as a function of molecular concentrations in experiment and as a function of molecular size and temperature via simulations. The resulting condensates were found to retain at least some degree of internal dynamics varying as a function of the molecular composition. The results suggest a more general principle for phase separation that is based primarily on electrostatic complementarity without invoking polymer properties as in most previous studies. Simulation results furthermore suggest that such phase separation may occur widely in heterogenous cellular environment between nucleic acid and protein components.

Molecular signatures of beef tenderness: Underlying mechanisms based on integromics of protein biomarkers from multi-platform proteomics studies

Over the last two decades, proteomics have been employed to decipher the underlying factors contributing to variation in the quality of muscle foods, including beef tenderness. One such approach is the application of high-throughput protein analytical platforms in the identification of meat quality biomarkers. To broaden our understanding about the biological mechanisms underpinning meat tenderization across a large number of studies, an integromics study was performed to review the current status of protein biomarker discovery targeting beef tenderness. This meta-analysis is the first to gather and propose a comprehensive list of 124 putative protein biomarkers derived from 28 independent proteomics-based experiments, from which 33 robust candidates were identified worthy of evaluation using targeted or untargeted data-independent acquisition proteomic methods. We further provide an overview of the interconnectedness of the main biological pathways impacting tenderness determination after multistep analyses including Gene Ontology annotations, pathway and process enrichment and literature mining, and specifically discuss the major proteins and pathways most often reported in proteomics research.

Prediction of Deleterious Non-synonymous SNPs of Human STK11 Gene by Combining Algorithms, Molecular Docking, and Molecular Dynamics Simulation

Serine-threonine kinase11 (STK11) is a tumor suppressor gene which plays a key role in regulating cell growth and apoptosis. It is widely known as a multitasking kinase and engaged in cell polarity, cell cycle arrest, chromatin remodeling, energy metabolism, and Wnt signaling. The substitutions of single amino acids in highly conserved regions of the STK11 protein are associated with Peutz-Jeghers syndrome (PJS), which is an autosomal dominant inherited disorder. The abnormal function of the STK11 protein is still not well understood. In this study, we classified disease susceptible single nucleotide polymorphisms (SNPs) in STK11 by using different computational algorithms. We identified the deleterious nsSNPs, constructed mutant protein structures, and evaluated the impact of mutation by employing molecular docking and molecular dynamics analysis. Our results show that W239R and W308C variants are likely to be highly deleterious mutations found in the catalytic kinase domain, which may destabilize structure and disrupt the activation of the STK11 protein as well as reduce its catalytic efficiency. The W239R mutant is likely to have a greater impact on destabilizing the protein structure compared to the W308C mutant. In conclusion, these mutants can help to further realize the large pool of disease susceptibilities linked with catalytic kinase domain activation of STK11 and assist to develop an effective drug for associated diseases.

Local and global anatomy of antibody-protein antigen recognition

Deciphering antibody-protein antigen recognition is of fundamental and practical significance. We constructed an antibody structural dataset, partitioned it into human and murine subgroups, and compared it with nonantibody protein-protein complexes. We investigated the physicochemical properties of regions on and away from the antibody-antigen interfaces, including net charge, overall antibody charge distributions, and their potential role in antigen interaction. We observed that amino acid preference in antibody-protein antigen recognition is entropy driven, with residues having low side-chain entropy appearing to compensate for the high backbone entropy in interaction with protein antigens. Antibodies prefer charged and polar antigen residues and bridging water molecules. They also prefer positive net charge, presumably to promote interaction with negatively charged protein antigens, which are common in proteomes. Antibody-antigen interfaces have large percentages of Tyr, Ser, and Asp, but little Lys. Electrostatic and hydrophobic interactions in the Ag binding sites might be coupled with Fab domains through organized charge and residue distributions away from the binding interfaces. Here we describe some features of antibody-antigen interfaces and of Fab domains as compared with nonantibody protein-protein interactions. The distributions of interface residues in human and murine antibodies do not differ significantly. Overall, our results provide not only a local but also a global anatomy of antibody structures.

Protein charge determination and implications for interactions in cell extracts

Decades of dilute-solution studies have revealed the influence of charged residues on protein stability, solubility and stickiness. Similar characterizations are now required in physiological solutions to understand the effect of charge on protein behavior under native conditions. Toward this end, we used free boundary and native gel electrophoresis to explore the charge of cytochrome c in buffer and in Escherichia coli extracts. We find that the charge of cytochrome c was ∼2-fold lower than predicted from primary structure analysis. Cytochrome c charge was tuned by sulfate binding and was rendered anionic in E. coli extracts due to interactions with macroanions. Mutants in which three or four cationic residues were replaced with glutamate were charge-neutral and "inert" in extracts. A comparison of the interaction propensities of cytochrome c and the mutants emphasizes the role of negative charge in stabilizing physiological environments. Charge-charge repulsion and preferential hydration appear to prevent aggregation. The implications for molecular organization in vivo are discussed.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Exploring proteome-wide occurrence of clusters of charged residues in eukaryotes

Clusters of charged residues are one of the key features of protein primary structure since they have been associated to important functions of proteins. Here, we present a proteome wide scan for the occurrence of Charge Clusters in Protein sequences using a new search tool (FCCP) based on a score-based methodology. The FCCP was run to search charge clusters in seven eukaryotic proteomes: Arabidopsis thaliana, Caenorhabditis elegans, Danio rerio, Drosophila melanogaster, Homo sapiens, Mus musculus, and Saccharomyces cerevisiae. We found that negative charge clusters (NCCs) are three to four times more frequent than positive charge clusters (PCCs). The Drosophila proteome is on average the most charged, whereas the human proteome is the least charged. Only 3 to 8% of the studied protein sequences have negative charge clusters, while 1.6 to 3% having PCCs and only 0.07 to 0.6% have both types of clusters. NCCs are localized predominantly in the N-terminal and C-terminal domains, while PCCs tend to be localized within the functional domains of the protein sequences. Furthermore, the gene ontology classification revealed that the protein sequences with negative and PCCs are mainly binding proteins.

Structured crowding and its effects on enzyme catalysis

Macromolecular crowding decreases the diffusion rate, shifts the equilibrium of protein-protein and protein-substrate interactions, and changes protein conformational dynamics. Collectively, these effects contribute to enzyme catalysis. Here we describe how crowding may bias the conformational change and dynamics of enzyme populations and in this way affect catalysis. Crowding effects have been studied using artificial crowding agents and in vivo-like environments. These studies revealed a correlation between protein dynamics and function in the crowded environment. We suggest that crowded environments be classified into uniform crowding and structured crowding. Uniform crowding represents random crowding conditions created by synthetic particles with a narrow size distribution. Structured crowding refers to the highly coordinated cellular environment, where proteins and other macromolecules are clustered and organized. In structured crowded environments the perturbation of protein thermal stability may be lower; however, it may still be able to modulate functions effectively and dynamically. Dynamic, allosteric enzymes could be more sensitive to cellular perturbations if their free energy landscape is flatter around the native state; on the other hand, if their free energy landscape is rougher, with high kinetic barriers separating deep minima, they could be more robust. Above all, cells are structured; and this holds both for the cytosol and for the membrane environment. The crowded environment is organized, which limits the search, and the crowders are not necessarily inert. More likely, they too transmit allosteric effects, and as such play important functional roles. Overall, structured cellular crowding may lead to higher enzyme efficiency and specificity.