The ConSurf-DB: pre-calculated evolutionary conservation profiles of protein structures

Identification of deleterious non-synonymous single nucleotide polymorphisms in the mRNA decay activator ZFP36L2

More than 4,000 single nucleotide polymorphisms (SNP) variants have been identified in the human ZFP36L2 gene, however only a few have been studied in the context of protein function. The tandem zinc finger domain of ZFP36L2, an RNA binding protein, is the functional domain that binds to its target mRNAs. This protein/RNA interaction triggers mRNA degradation, controlling gene expression. We identified 32 non-synonymous SNPs (nsSNPs) in the tandem zinc finger domain of ZFP36L2 that could have possible deleterious impacts in humans. Using different bioinformatic strategies, we prioritized five among these 32 nsSNPs, namely rs375096815, rs1183688047, rs1214015428, rs1215671792 and rs920398592 to be validated. When we experimentally tested the functionality of these protein variants using gel shift assays, all five (Y154H, R160W, R184C, G204D, and C206F) resulted in a dramatic reduction in RNA binding compared to the WT protein. To understand the mechanistic effect of these variants on the protein/RNA interaction, we employed DUET, DynaMut and PyMOL to investigate structural changes in the protein. Additionally, we conducted Molecular Docking and Molecular Dynamics Simulations to fine tune the active behaviour of this biomolecular system at an atomic level. Our results propose atomic explanations for the impact of each of these five genetic variants identified.

The enterobactin biosynthetic intermediate 2,3-dihydroxybenzoic acid is a competitive inhibitor of the Escherichia coli isochorismatase EntB

The Escherichia coli enterobactin biosynthetic protein EntB is a bifunctional enzyme that catalyzes hydrolysis of isochorismate via its N-terminal isochorismatase (IC) domain, and then transfers phosphopantetheinylated 2,3-DHB to EntF via the EntB C-terminal aryl carrier protein (ArCP) domain. Here we used a fluorescence anisotropy binding assay to investigate the ability of 2,3-DHB to bind to enzymes in the DHB synthetic arm of the pathway. We found that 2,3-DHB binds to EntE as a natural substrate with high affinity (KD = 0.54 μM). Furthermore, apo-EntB was found to bind to 2,3-DHB with moderate affinity (KD = 8.95 μM), despite the fact that this intermediate is neither a substrate nor a product of EntB. Molecular docking simulations predicted a top-ranked ensemble in which 2,3-DHB is bound at the isochorismatase active site of apo-EntB. Steady-state coupled enzymatic assays revealed that 2,3-DHB is a competitive inhibitor of apo-EntB isochorismatase activity (Ki ~ 200 μM), consistent with modeling predictions. Monitoring the EntC-EntB coupled reaction in real time via isothermal titration microcalorimetry confirmed that EntB was required to drive the EntC reaction toward isochorismate formation. Furthermore, addition of 2,3-DHB to the ITC-monitored reaction resulted in a suppression of integrated reaction heats, consistent with our observation that the molecule acts as a competitive inhibitor of EntB. Finally, we found that 2,3-DHB lowered the efficiency of EntC-EntB isochorismate channeling by approximately 70%, consistent with steric blockage of the isochorismatase active site by bound 2,3-DHB. Given its inhibitory properties, we hypothesize that 2,3-DHB plays a regulatory role in feedback inhibition in order to maintain iron homeostasis upon intracellular accumulation of sufficient ferric enterobactin.

Clinical Exome Sequencing Identifies, Two Homozygous LOXHD1 Variants in Two Inbred Families With Pre-Lingual Hearing Loss From South India

In recent years, numerous genetic variants have been linked with prelingual hearing loss (HL). Variants in the LOXHD1 gene (lipoxygenase homology domain-1) associated with DFNB77 are highly heterogeneous, with different auditory characteristics varying from stable to progressive and mild to profound. To date, 168 DFNB77 cases have been recorded worldwide. Forty-one hearing-impaired (HI) probands, who were previously excluded for a set of four common deafness-causing genes (viz., GJB2, GJB6, SLC26A4, and CDH23) from 33 HI families, were subjected to clinical exome sequencing (CES) involving 285 genes associated with HL. This was followed by a segregation analysis of the available members in the family. We identified two pathogenic LOXHD1 variants in two unrelated inbred families. One is a novel homozygous pathogenic nonsense variant (c.3999C > A; p.C1333X), whereas the other is a likely pathogenic missense variant (c.6046G > T; p.E2046K). In silico tools such as SIFT, PolyPhen-2, Mutation Taster, CADD, and REVEL scores were used to predict variant pathogenicity. Furthermore, American College of Medical Genetics and Genomics guidelines specific to HL were applied to finally classify a variant as pathogenic or otherwise. The frequency of LOXHD1 variants identified in our study is 4.88% (2/41). This is the first LOXHD1 report associated with non-syndromic HL in South Indian families.

Dynamic Allostery: Evolution's Double-Edged Sword in Protein Function and Disease

Allostery is a core mechanism in biology that allows proteins to communicate and regulate activity over long structural distances. While classical models of allostery focus on conformational changes triggered by ligand binding, dynamic allostery-where protein function is modulated through alterations in thermal fluctuations without major conformational shifts-has emerged as a critical evolutionary mechanism. This review explores how evolution leverages dynamic allostery to fine-tune protein function through subtle mutations at distal sites, preserving core structural architecture while dramatically altering functional properties. Using a combination of computational approaches including Dynamic Flexibility Index (DFI), Dynamic Coupling Index (DCI), and vibrational density of states (VDOS) analysis, we demonstrate that functional adaptations in proteins often involve "hinge-shift" mechanisms, where redistribution of rigid and flexible regions modulates collective motions without changing the overall fold. This evolutionary principle is a double-edged sword: the same mechanisms that enable functional innovation also create vulnerabilities that can be exploited in disease states. Disease-associated variants frequently occur at positions highly coupled to functional sites despite being physically distant, forming Dynamic Allosteric Residue Couples (DARC sites). We demonstrate applications of these principles in understanding viral evolution, drug resistance, and capsid assembly dynamics. Understanding dynamic allostery provides critical insights into protein evolution and offers new avenues for therapeutic interventions targeting allosteric regulation.

Adoption of an in-silico analysis approach to assess the functional and structural impacts of rpoB-encoded protein mutations on Chlamydia pneumoniae sensitivity to antibiotics

BACKGROUND

Antibiotics are frequently used to treat infections caused by Chlamydia pneumoniae; an obligate intracellular gram-negative bacterium commonly associated with respiratory diseases. However, improper or overuse of these drugs has raised concerns about the development of antibiotic resistance, which poses a significant global health challenge. Previous studies have revealed a link between mutations in the rpoB-encoded protein of C. pneumoniae and antibiotic resistance. This study assessed these mutations via various bioinformatics tools to predict their impact on function, structural stability, antibiotic binding, and, ultimately, their effect on bacterial sensitivity to antibiotics.

RESULTS

Eight mutations in the rpoB-encoded protein (R421S, F450S, L456I, S454F, D461E, S476F, L478S, and S519Y) are associated with resistance to rifampin and rifalazil. These mutations occur in conserved regions of the protein, leading to decreased stability and affecting essential functional sites of RNA polymerase, the target of these antibiotics. Although the structural differences between the native and mutant proteins are minimal, notable changes in local hydrogen bonding have been observed. Despite similar binding energies, variations in hydrogen bonds and hydrophobic interactions in certain mutants (for instance, D461E for rifalazil and S476F for rifampin) indicate that these changes may diminish ligand affinity and specificity. Furthermore, protein-protein network analysis demonstrated a strong correlation between wild-type rpoB and ten C. pneumoniae proteins, each fulfilling specific functional roles. Consequently, some of these mutations can reduce the bacterium's sensitivity to rifampin and rifalazil, thereby contributing to antibiotic resistance.

CONCLUSION

The findings of this study indicate that mutations in the rpoB gene, which encodes the beta subunit of RNA polymerase, are pivotal in the resistance of C. pneumoniae to rifampin and rifalazil. Some of these mutations may result in reduced protein stability and changes in the structure, function, and antibiotic binding. As a consequence, the efficacy of these drugs in inhibiting RNA polymerase is compromised, allowing the bacteria to persist in transcription and replication even in the presence of antibiotics. Overall, these insights enhance our understanding of the resistance mechanisms in C. pneumoniae and could guide the development of strategies to address this challenge.

CLINICAL TRIAL NUMBER

Not applicable.

Enhanced thermostability of nattokinase by rational design of disulfide bond

Nattokinase, the thrombolytically active substance in the health food natto, nevertheless, its lower thermostability restricts its use in food and pharmaceutical applications. In this study, two heat-resistant variants of nattokinase, designated 50-109 (M1) and 15-271 (M2), were successfully obtained by introducing a disulfide bonding strategy. Their half-lives at 55℃ were found to be 2.50-fold and 5.17-fold higher, respectively, than that of the wild type. Furthermore, the specific enzyme activities of the variants, M1 and M2, were also increased by 2.37 and 1.66-fold, respectively. Meanwhile, the combination of two mutants increased the thermostability of nattokinase by 8.0-fold. Bioinformatics analyses indicated that the enhanced thermostability of the M1 and M2 variants was due to the increased rigidity and structural contraction of the overall structure. Finally, the fermentation process of mutant M1 was optimized to increase the expression of nattokinase. Study provides substantial molecular and theoretical support for the industrial production and application of nattokinase.

An ADPRS variant disrupts ARH3 stability and subcellular localization in children with neurodegeneration and respiratory failure

ADP-ribosylation is a post-translational modification involving the transfer of one or more ADP-ribose units from NAD+ to target proteins. Dysregulation of ADP-ribosylation is implicated in neurodegenerative diseases. In this study, genetic testing via exome sequencing was used to identify the underlying disease in two siblings with developmental delay, seizures, progressive muscle weakness, and respiratory failure following an episodic course. This identified a novel homozygous variant in the ADPRS gene (c.545A>G, p.His182Arg) encoding the mono(ADP-ribosyl) hydrolase ARH3, confirming the diagnosis of childhood-onset neurodegeneration with stress-induced ataxia and seizures (CONDSIAS) in these two children. Mechanistically, the ARH3H182R variant affects a highly conserved residue in the active site of ARH3, leading to protein instability, degradation, and, subsequently, reduced protein expression. The ARH3H182R mutant additionally fails to localize to the nucleus, which further resulted in accumulated mono-ADP ribosylated species in cells. The children's clinical course combined with the biochemical characterization of their genetic variant develops our understanding of the pathogenic mechanisms driving CONDSIAS and highlights a critical role for ARH3-regulated ADP-ribosylation in nervous system integrity.

Computer-Aided Design to Improve the Thermal Stability of Rhizomucor miehei Lipase

Lipase, a green biocatalyst, finds extensive application in the food sector. Enhancing the thermal stability of lipase presents both challenges and opportunities within the food industry. This research employed multiple rounds of cross-screening using tools like FoldX and I-Mutant 3.0 to strategically design mutations for Rhizomucor miehei lipase (RML), resulting in eight unique single-point mutation designs. E230I, N120M, and N264M have been confirmed experimentally to be potential combination mutation candidates. The resulting triple mutant N120M/E230I/N264M showed a higher thermal stability, with an optimum temperature of 55 °C, 10 °C higher than that of the wild-type RML. The half-life was extended from 46 to 462 min at 50 °C. Furthermore, the catalytic activity of N120M/E230I/N264M on camphor tree seed oil increased by 140% to 600 U/mg, which aids in the production of novel structured lipids. Using molecular docking and molecular dynamics simulations, we analyzed the molecular mechanism of enhanced thermal stability. This study validated the efficacy and dependability of computer-aided design to generate heat-resistant RML mutants and indicated that RML N120M/E230I/N264M lipase can be used as an effective biocatalyst for fat processing in the food industry.

Modification of the active centre of nattokinase to enhance its thermostability using a strategy based on molecular dynamics simulation, steered dynamics simulation, and conservative prediction

Introduction

The poor thermostability of nattokinase represents a significant limitation in its potential applications. Additionally, there is a notable absence of studies focused on modifying residues within the active site region of nattokinase with the aim of enhancing its catalytic properties. Furthermore, the direct utilisation of directed evolution often yields unfavourable outcomes, with a considerable workload being a common consequence.

Methods

In order to solve the above problems, a new method based on molecular dynamics simulation, steered dynamics simulation and conservative analysis with site-directed mutagenesis was proposed to screen nattokinase mutants with improved thermal stability. Molecular dynamics simulation was used to explain the mechanism of catalytic performance improvement of positive mutants. Finally, the fermentation process of the positive mutant was optimized.

Results and discussion

Based on these findings, the mutant A216K was selected for a 5.7-fold increase in half-life at 55°C with a small increase in activity, which further enhanced the mutation library of the thermal stability enhancement site in the enzyme's active centre. The results of the molecular dynamics simulation indicated that the enhancement of the number of hydrogen bonds within the protein and between the protein and the solvent, as well as the augmentation of the rigidity around the calcium ion binding site and the mutation site, were the primary factors contributing to the improvement of the thermal stability of A216K. It is anticipated that this strategy will provide novel insights into enzyme engineering research.

Characterization of atypical BAR domain-containing proteins coded by Toxoplasma gondii

Toxoplasma gondii, the causative agent of toxoplasmosis, infects cells and replicates inside via the secretion of factors stored in specialized organelles (rhoptries, micronemes, and dense granules) and the capture of host materials. The genesis of the secretory organelles and the processes of secretion and endocytosis depend on vesicular trafficking events whose molecular bases remain poorly known. Notably, there is no characterization of the BAR (Bin/Amphiphysin/Rvs) domain-containing proteins expressed by T. gondii and other apicomplexans, although such proteins are known to play critical roles in vesicular trafficking in other eukaryotes. Here, by combining structural analyses with in vitro assays and cellular observations, we have characterized TgREMIND (regulators of membrane interacting domains), involved in the genesis of rhoptries and dense granules, and TgBAR2 found at the parasite cortex. We establish that TgREMIND comprises an F-BAR domain that can bind curved neutral membranes with no strict phosphoinositide requirement and exert a membrane remodeling activity. Next, we establish that TgREMIND contains a new structural domain called REMIND, which negatively regulates the membrane-binding capacities of the F-BAR domain. In parallel, we report that TgBAR2 contains a BAR domain with an extremely basic membrane-binding interface able to deform anionic membranes into very narrow tubules. Our data show that T. gondii codes for two atypical BAR domain-containing proteins with very contrasting membrane-binding properties, allowing them to function in two distinct regions of the parasite trafficking system.

Post-translational modification sites are present in hydrophilic cavities of alpha-synuclein, tau, FUS, and TDP-43 fibrils: A molecular dynamics study

Hydration plays a crucial role in the refolding of intrinsically disordered proteins into amyloid fibrils; however, the specific interactions between water and protein that may contribute to this process are still unknown. In our previous studies of alpha-synuclein (aSyn), we have shown that waters confined in fibril cavities are stabilizing features of this pathological fold; and that amino acids that hydrogen bond with these confined waters modulate primary and seeded aggregation. Here, we extend our aSyn molecular dynamics (MD) simulations with three new polymorphs and correlate MD trajectory information with known post-translational modifications (PTMs) and experimental data. We show that cavity residues are more evolutionarily conserved than non-cavity residues and are enriched with PTM sites. As expected, the confinement within hydrophilic cavities results in more stably hydrated amino acids. Interestingly, cavity PTM sites display the longest protein-water hydrogen bond lifetimes, three-fold greater than non-PTM cavity sites. Utilizing the deep mutational screen dataset by Newberry et al. and the Thioflavin T aggregation review by Pancoe et al. parsed using a fibril cavity/non-cavity definition, we show that hydrophobic changes to amino acids in cavities have a larger effect on fitness and aggregation rate than residues outside cavities, supporting our hypothesis that these sites are involved in the inhibition of aSyn toxic fibrillization. Finally, we expand our study to include analysis of fibril structures of tau, FUS, TDP-43, prion, and hnRNPA1; all of which contained hydrated cavities, with tau, FUS, and TDP-43 recapitulating our PTM results in aSyn fibril cavities.

Enhanced Thermostability of Nattokinase by Computation-Based Rational Redesign of Flexible Regions

Nattokinase is a nutrient in healthy food natto that has the function of preventing and treating blood thrombus. However, its low thermostability and fibrinolytic activity limit its application in food and pharmaceuticals. In this study, we used bioinformatics analysis to identify two loops (loop10 and loop12) in the flexible region of nattokinase rAprY. Using this basis, we screened the G131S-S161T variant, which showed a 2.38-fold increase in half-life at 55 °C, and the M3 variant, which showed a 2.01-fold increase in activity, by using a thermostability prediction algorithm. Bioinformatics analysis revealed that the enhanced thermostability of the G131S-S161T variant was due to the increased rigidity and structural shrinkage of the overall structure. Additionally, the increased rigidity of the local region surrounding the active center and its mutated sites helps maintain its normal conformation in high-temperature environments. The increased catalytic activity of the M3 variant may be due to its more efficient substrate binding mechanism. We investigated strategies to improve the thermostability and fibrinolytic activity of nattokinase, and the resulting variants show promise for industrial production and application.

Modeling allosteric mechanisms of eukaryotic type II topoisomerases

Type II topoisomerases (TopoIIs) are ubiquitous enzymes that are involved in crucial nuclear processes such as genome organization, chromosome segregation, and other DNA metabolic processes. These enzymes function as large, homodimeric complexes that undergo a complex cycle of binding and hydrolysis of two ATP molecules in their ATPase domains, which regulates the capture and passage of one DNA double-helix through a second, cleaved DNA molecule. This process requires the transmission of information about the state of the bound nucleotide over vast ranges in the TopoII complex. How this information is transmitted at the molecular level to regulate TopoII functions and how protein substitutions disrupt these mechanisms remains largely unknown. Here, we employed extensive microsecond-scale molecular dynamics simulations of the yeast TopoII enzyme in multiple nucleotide-bound states and with amino acid substitutions near both the N and C termini of the complex. Simulation results indicate that the ATPase domains are remarkably flexible on the sub-microsecond timescale and that these dynamics are modulated by the identity of the bound nucleotides and both local and distant amino acid substitutions. Network analyses point toward specific allosteric networks that transmit information about the hydrolysis cycle throughout the complex, which include residues in both the protein and the bound DNA molecule. Amino acid substitutions weaken many of these pathways. Together, our results provide molecular level details on how the TopoII catalytic cycle is controlled through nucleotide binding and hydrolysis and how mutations may disrupt this process.

A novel variant in ADPRS disrupts ARH3 stability and subcellular localization in children with neurodegeneration and respiratory failure

Purpose

ADP-ribosylation is a post-translational modification involving the transfer of one or more ADP-ribose units from NAD+ to target proteins. Dysregulation of ADP-ribosylation is implicated in neurodegenerative diseases. Here we report a novel homozygous variant in the ADPRS gene (c.545A>G, p.His182Arg) encoding the mono(ADP-ribosyl) hydrolase ARH3 found in 2 patients with childhood-onset neurodegeneration with stress-induced ataxia and seizures (CONDSIAS).

Methods

Genetic testing via exome sequencing was used to identify the underlying disease cause in two siblings with developmental delay, seizures, progressive muscle weakness, and respiratory failure following an episodic course. Studies in a cell culture model uncover biochemical and cellular consequences of the identified genetic change.

Results

The ARH3 H182R variant affects a highly conserved residue in the active site of ARH3, leading to protein instability, degradation, and reduced expression. ARH3 H182R additionally fails to localize to the nucleus. The combination of reduced expression and mislocalization of ARH3 H182R resulted in accumulation of mono-ADP ribosylated species in cells.

Conclusions

The children's clinical course combined with the biochemical characterization of their genetic variant develops our understanding of the pathogenic mechanisms driving CONDSIAS and highlights a critical role for ARH3-regulated ADP ribosylation in nervous system integrity.

Unraveling the Genetic Basis of Combined Deafness and Male Infertility Phenotypes through High-Throughput Sequencing in a Unique Cohort from South India

The co-occurrence of sensorineural hearing loss and male infertility has been reported in several instances, suggesting potential shared genetic underpinnings. One such example is the contiguous gene deletion of CATSPER2 and STRC genes, previously associated with deafness-infertility syndrome (DIS) in males. Fifteen males with both hearing loss and infertility from southern India after exclusion for the DIS contiguous gene deletion and the FOXI1 gene mutations are subjected to exome sequencing. This resolves the genetic etiology in four probands for both the phenotypes; In the remaining 11 probands, two each conclusively accounted for deafness and male infertility etiologies. Genetic heterogeneity is well reflected in both phenotypes. Four recessive (TRIOBP, SLC26A4, GJB2, COL4A3) and one dominant (SOX10) for the deafness; six recessive genes (LRGUK, DNAH9, ARMC4, DNAH2, RSPH6A, and ACE) for male infertility can be conclusively ascribed. LRGUK and RSPH6A genes are implicated earlier only in mice models, while the ARMC4 gene is implicated in chronic destructive airway diseases due to primary ciliary dyskinesia. This study would be the first to document the role of these genes in the male infertility phenotype in humans. The result suggests that deafness and infertility are independent events and do not segregate together among the probands.

Modeling Allosteric Mechanisms of Eukaryotic Type II Topoisomerases

Type II topoisomerases (TopoIIs) are ubiquitous enzymes that are involved in crucial nuclear processes such as genome organization, chromosome segregation, and other DNA metabolic processes. These enzymes function as large, homodimeric complexes that undergo a complex cycle of binding and hydrolysis of two ATP molecules in their ATPase domains, which regulates the capture and passage of one DNA double-helix through a second, cleaved DNA molecule. This process requires the transmission of information about the state of the bound nucleotide over vast ranges in the TopoII complex. How this information is transmitted at the molecular level to regulate TopoII functions and how protein substitutions disrupt these mechanisms remains largely unknown. Here, we employed extensive microsecond scale molecular dynamics simulations of the yeast TopoII enzyme in multiple nucleotide-bound states and with amino acid substitutions near both the N- and C-terminals of the complex. Simulation results indicate that the ATPase domains are remarkably flexible on the sub-microsecond timescale and that these dynamics are modulated by the identity of the bound nucleotides and both local and distant amino acid substitutions. Network analyses point towards specific allosteric networks that transmit information about the hydrolysis cycle throughout the complex, which include residues in both the protein and the bound DNA molecule. Amino acid substitutions weaken many of these pathways. Together, our results provide molecular-level details on how the TopoII catalytic cycle is controlled through nucleotide binding and hydrolysis and how mutations may disrupt this process.

Empowering Protein Engineering through Recombination of Beneficial Substitutions

Directed evolution stands as a seminal technology for generating novel protein functionalities, a cornerstone in biocatalysis, metabolic engineering, and synthetic biology. Today, with the development of various mutagenesis methods and advanced analytical machines, the challenge of diversity generation and high-throughput screening platforms is largely solved, and one of the remaining challenges is: how to empower the potential of single beneficial substitutions with recombination to achieve the epistatic effect. This review overviews experimental and computer-assisted recombination methods in protein engineering campaigns. In addition, integrated and machine learning-guided strategies were highlighted to discuss how these recombination approaches contribute to generating the screening library with better diversity, coverage, and size. A decision tree was finally summarized to guide the further selection of proper recombination strategies in practice, which was beneficial for accelerating protein engineering.

The catalytic mechanism of the RNA methyltransferase METTL3

The complex of methyltransferase-like proteins 3 and 14 (METTL3-14) is the major enzyme that deposits N6-methyladenosine (m6A) modifications on messenger RNA (mRNA) in humans. METTL3-14 plays key roles in various biological processes through its methyltransferase (MTase) activity. However, little is known about its substrate recognition and methyl transfer mechanism from its cofactor and methyl donor S-adenosylmethionine (SAM). Here, we study the MTase mechanism of METTL3-14 by a combined experimental and multiscale simulation approach using bisubstrate analogues (BAs), conjugates of a SAM-like moiety connected to the N6-atom of adenosine. Molecular dynamics simulations based on crystal structures of METTL3-14 with BAs suggest that the Y406 side chain of METTL3 is involved in the recruitment of adenosine and release of m6A. A crystal structure with a BA representing the transition state of methyl transfer shows a direct involvement of the METTL3 side chains E481 and K513 in adenosine binding which is supported by mutational analysis. Quantum mechanics/molecular mechanics (QM/MM) free energy calculations indicate that methyl transfer occurs without prior deprotonation of adenosine-N6. Furthermore, the QM/MM calculations provide further support for the role of electrostatic contributions of E481 and K513 to catalysis. The multidisciplinary approach used here sheds light on the (co)substrate binding mechanism, catalytic step, and (co)product release, and suggests that the latter step is rate-limiting for METTL3. The atomistic information on the substrate binding and methyl transfer reaction of METTL3 can be useful for understanding the mechanisms of other RNA MTases and for the design of transition state analogues as their inhibitors.

Bi-allelic variants in CELSR3 are implicated in central nervous system and urinary tract anomalies

CELSR3 codes for a planar cell polarity protein. We describe twelve affected individuals from eleven independent families with bi-allelic variants in CELSR3. Affected individuals presented with an overlapping phenotypic spectrum comprising central nervous system (CNS) anomalies (7/12), combined CNS anomalies and congenital anomalies of the kidneys and urinary tract (CAKUT) (3/12) and CAKUT only (2/12). Computational simulation of the 3D protein structure suggests the position of the identified variants to be implicated in penetrance and phenotype expression. CELSR3 immunolocalization in human embryonic urinary tract and transient suppression and rescue experiments of Celsr3 in fluorescent zebrafish reporter lines further support an embryonic role of CELSR3 in CNS and urinary tract formation.

A cryptic plasmid is among the most numerous genetic elements in the human gut

Plasmids are extrachromosomal genetic elements that often encode fitness-enhancing features. However, many bacteria carry "cryptic" plasmids that do not confer clear beneficial functions. We identified one such cryptic plasmid, pBI143, which is ubiquitous across industrialized gut microbiomes and is 14 times as numerous as crAssphage, currently established as the most abundant extrachromosomal genetic element in the human gut. The majority of mutations in pBI143 accumulate in specific positions across thousands of metagenomes, indicating strong purifying selection. pBI143 is monoclonal in most individuals, likely due to the priority effect of the version first acquired, often from one's mother. pBI143 can transfer between Bacteroidales, and although it does not appear to impact bacterial host fitness in vivo, it can transiently acquire additional genetic content. We identified important practical applications of pBI143, including its use in identifying human fecal contamination and its potential as an alternative approach to track human colonic inflammatory states.

To kill a microRNA: emerging concepts in target-directed microRNA degradation

MicroRNAs (miRNAs) guide Argonaute (AGO) proteins to bind mRNA targets. Although most targets are destabilized by miRNA-AGO binding, some targets induce degradation of the miRNA instead. These special targets are also referred to as trigger RNAs. All triggers identified thus far have binding sites with greater complementarity to the miRNA than typical target sites. Target-directed miRNA degradation (TDMD) occurs when trigger RNAs bind the miRNA-AGO complex and recruit the ZSWIM8 E3 ubiquitin ligase, leading to AGO ubiquitination and proteolysis and subsequent miRNA destruction. More than 100 different miRNAs are regulated by ZSWIM8 in bilaterian animals, and hundreds of trigger RNAs have been predicted computationally. Disruption of individual trigger RNAs or ZSWIM8 has uncovered important developmental and physiologic roles for TDMD across a variety of model organisms and cell types. In this review, we highlight recent progress in understanding the mechanistic basis and functions of TDMD, describe common features of trigger RNAs, outline best practices for validating trigger RNAs, and discuss outstanding questions in the field.

Phylogenomics analysis of velvet regulators in the fungal kingdom

All life forms have evolved to respond appropriately to various environmental and internal cues. In the animal kingdom, the prototypical regulator class of such cellular responses is the Rel homology domain proteins including nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB). Fungi, the close relatives of animals, have also evolved with their own NF-κB-like regulators called velvet family proteins to govern cellular and chemical development. Here, we conducted a detailed investigation of the taxonomic broad presence of velvet proteins. We observed that velvet proteins are widely distributed in the fungal kingdom. Moreover, we have identified and characterized 21 major velvet clades in fungi. We have further revealed that the highly conserved velvet domain is composed of three distinct motifs and acts as an evolutionarily independent domain, which can be shuffled with various functional domains. Such rearrangements of the velvet domain have resulted in the functional and type diversity of the present velvet regulators. Importantly, our in-deep analyses of the primary and 3D structures of the various velvet domains showed that the fungal velvet domains can be divided into two major clans: the VelB and the VosA clans. The 3D structure comparisons revealed a close similarity of the velvet domain with many other eukaryotic DNA-binding proteins, including those of the Rel, Runt, and signal transducer and activator of transcription families, sharing a common β-sandwich fold. Altogether, this study improves our understanding of velvet regulators in the fungal kingdom.IMPORTANCEFungi are the relatives of animals in Opisthokonta and closely associated with human life by interactive ways such as pathogenicity, food, and secondary metabolites including beneficial ones like penicillin and harmful ones like the carcinogenic aflatoxins. Similar to animals, fungi have also evolved with NF-κB-like velvet family regulators. The velvet proteins constitute a large protein family of fungal transcription factors sharing a common velvet domain and play a key role in coordinating fungal secondary metabolism, developmental and differentiation processes. Our current understanding on velvet regulators is mostly from Ascomycota fungi; however, they remain largely unknown outside Ascomycota. Therefore, this study performed a taxonomic broad investigation of velvet proteins across the fungal kingdom and conducted a detailed analysis on velvet distribution, structure, diversity, and evolution. The results provide a holistic view of velvet regulatory system in the fungal kingdom.

Improvement of the fibrinolytic activity, acid resistance and thermostability of nattokinase by surface charge engineering

Nattokinase is a promising thrombolytic drug due to its powerful fibrinolytic effect and few side effects. However, the low fibrinolytic activity and stability of nattokinase have limited its industrial production and oral application. In this study, the basic and neutral amino acid residues on the surface of recombinant nattokinase AprY from Bacillus mojavensis LY-06 (rAprY) were mutated to acidic amino acid residues by surface charge engineering strategy, and two variants K12D and N109D with 92.6 % and 8.4 % increased fibrinolytic activity were obtained. The R45E variant with enhanced acid stability and thermostability was also screened, its acid stability at pH 4 and t1/2 at 55 °C were 3.7-fold and 1.8-fold higher than that of wild type rAprY, respectively. Bioinformatics analysis showed that the increased activities of K12D and N109D variants were related to the increased flexibility of the region around their active centers. The increased rigidity of 97-103 amino acid residues around the active center of R45E may be the reason for its enhanced stability and reduced catalytic activity. The multipoint mutation K12D-N109D (M2)'s catalytic activity did not increase cumulatively, but its pH stability did. The nattokinase variants generated in this study have potential for industrial production and application.

Constructing and interpreting a large-scale variant effect map for an ultrarare disease gene: Comprehensive prediction of the functional impact of PSAT1 genotypes

Reduced activity of the enzymes encoded by PHGDH, PSAT1, and PSPH causes a set of ultrarare, autosomal recessive diseases known as serine biosynthesis defects. These diseases present in a broad phenotypic spectrum: at the severe end is Neu-Laxova syndrome, in the intermediate range are infantile serine biosynthesis defects with severe neurological manifestations and growth deficiency, and at the mild end is childhood disease with intellectual disability. However, L-serine supplementation, especially if started early, can ameliorate and in some cases even prevent symptoms. Therefore, knowledge of pathogenic variants can improve clinical outcomes. Here, we use a yeast-based assay to individually measure the functional impact of 1,914 SNV-accessible amino acid substitutions in PSAT. Results of our assay agree well with clinical interpretations and protein structure-function relationships, supporting the inclusion of our data as functional evidence as part of the ACMG variant interpretation guidelines. We use existing ClinVar variants, disease alleles reported in the literature and variants present as homozygotes in the primAD database to define assay ranges that could aid clinical variant interpretation for up to 98% of the tested variants. In addition to measuring the functional impact of individual variants in yeast haploid cells, we also assay pairwise combinations of PSAT1 alleles that recapitulate human genotypes, including compound heterozygotes, in yeast diploids. Results from our diploid assay successfully distinguish the genotypes of affected individuals from those of healthy carriers and agree well with disease severity. Finally, we present a linear model that uses individual allele measurements to predict the biallelic function of ~1.8 million allele combinations corresponding to potential human genotypes. Taken together, our work provides an example of how large-scale functional assays in model systems can be powerfully applied to the study of ultrarare diseases.

Allosteric regulatory control in dihydrofolate reductase is revealed by dynamic asymmetry

We investigated the relationship between mutations and dynamics in Escherichia coli dihydrofolate reductase (DHFR) using computational methods. Our study focused on the M20 and FG loops, which are known to be functionally important and affected by mutations distal to the loops. We used molecular dynamics simulations and developed position-specific metrics, including the dynamic flexibility index (DFI) and dynamic coupling index (DCI), to analyze the dynamics of wild-type DHFR and compared our results with existing deep mutational scanning data. Our analysis showed a statistically significant association between DFI and mutational tolerance of the DHFR positions, indicating that DFI can predict functionally beneficial or detrimental substitutions. We also applied an asymmetric version of our DCI metric (DCIasym ) to DHFR and found that certain distal residues control the dynamics of the M20 and FG loops, whereas others are controlled by them. Residues that are suggested to control the M20 and FG loops by our DCIasym metric are evolutionarily nonconserved; mutations at these sites can enhance enzyme activity. On the other hand, residues controlled by the loops are mostly deleterious to function when mutated and are also evolutionary conserved. Our results suggest that dynamics-based metrics can identify residues that explain the relationship between mutation and protein function or can be targeted to rationally engineer enzymes with enhanced activity.

EvoRator2: Predicting Site-specific Amino Acid Substitutions Based on Protein Structural Information Using Deep Learning

Multiple sequence alignments (MSAs) are the workhorse of molecular evolution and structural biology research. From MSAs, the amino acids that are tolerated at each site during protein evolution can be inferred. However, little is known regarding the repertoire of tolerated amino acids in proteins when only a few or no sequence homologs are available, such as orphan and de novo designed proteins. Here we present EvoRator2, a deep-learning algorithm trained on over 15,000 protein structures that can predict which amino acids are tolerated at any given site, based exclusively on protein structural information mined from atomic coordinate files. We show that EvoRator2 obtained satisfying results for the prediction of position-weighted scoring matrices (PSSM). We further show that EvoRator2 obtained near state-of-the-art performance on proteins with high quality structures in predicting the effect of mutations in deep mutation scanning (DMS) experiments and that for certain DMS targets, EvoRator2 outperformed state-of-the-art methods. We also show that by combining EvoRator2's predictions with those obtained by a state-of-the-art deep-learning method that accounts for the information in the MSA, the prediction of the effect of mutation in DMS experiments was improved in terms of both accuracy and stability. EvoRator2 is designed to predict which amino-acid substitutions are tolerated in such proteins without many homologous sequences, including orphan or de novo designed proteins. We implemented our approach in the EvoRator web server (https://evorator.tau.ac.il).

Evolutionarily conserved amino acids in MHC-II mediate bat influenza A virus entry into human cells

The viral hemagglutinins of conventional influenza A viruses (IAVs) bind to sialylated glycans on host cell surfaces for attachment and subsequent infection. In contrast, hemagglutinins of bat-derived IAVs target major histocompatibility complex class II (MHC-II) for cell entry. MHC-II proteins from various vertebrate species can facilitate infection with the bat IAV H18N11. Yet, it has been difficult to biochemically determine the H18:MHC-II binding. Here, we followed a different approach and generated MHC-II chimeras from the human leukocyte antigen DR (HLA-DR), which supports H18-mediated entry, and the nonclassical MHC-II molecule HLA-DM, which does not. In this context, viral entry was supported only by a chimera containing the HLA-DR α1, α2, and β1 domains. Subsequent modeling of the H18:HLA-DR interaction identified the α2 domain as central for this interaction. Further mutational analyses revealed highly conserved amino acids within loop 4 (N149) and β-sheet 6 (V190) of the α2 domain as critical for virus entry. This suggests that conserved residues in the α1, α2, and β1 domains of MHC-II mediate H18-binding and virus propagation. The conservation of MHC-II amino acids, which are critical for H18N11 binding, may explain the broad species specificity of this virus.

Mitochondrial DNA Changes in Respiratory Complex I Genes in Brain Gliomas

Mitochondria are organelles necessary for oxidative phosphorylation. The interest in the role of mitochondria in the process of carcinogenesis results from the fact that a respiratory deficit is found in dividing cells, especially in cells with accelerated proliferation. The study included tumor and blood material from 30 patients diagnosed with glioma grade II, III and IV according to WHO (World Health Organization). DNA was isolated from the collected material and next-generation sequencing was performed on the MiSeqFGx apparatus (Illumina). The study searched for a possible relationship between the occurrence of specific mitochondrial DNA polymorphisms in the respiratory complex I genes and brain gliomas of grade II, III and IV. The impact of missense changes on the biochemical properties, structure and functioning of the encoded protein, as well as their potential harmfulness, were assessed in silico along with their belonging to a given mitochondrial subgroup. The A3505G, C3992T, A4024G, T4216C, G5046A, G7444A, T11253C, G12406A and G13604C polymorphisms were assessed as deleterious changes in silico, indicating their association with carcinogenesis.

A highly conserved and globally prevalent cryptic plasmid is among the most numerous mobile genetic elements in the human gut

Plasmids are extrachromosomal genetic elements that often encode fitness enhancing features. However, many bacteria carry 'cryptic' plasmids that do not confer clear beneficial functions. We identified one such cryptic plasmid, pBI143, which is ubiquitous across industrialized gut microbiomes, and is 14 times as numerous as crAssphage, currently established as the most abundant genetic element in the human gut. The majority of mutations in pBI143 accumulate in specific positions across thousands of metagenomes, indicating strong purifying selection. pBI143 is monoclonal in most individuals, likely due to the priority effect of the version first acquired, often from one's mother. pBI143 can transfer between Bacteroidales and although it does not appear to impact bacterial host fitness in vivo, can transiently acquire additional genetic content. We identified important practical applications of pBI143, including its use in identifying human fecal contamination and its potential as an inexpensive alternative for detecting human colonic inflammatory states.

Domain Architecture of the Nonreceptor Tyrosine Kinase Ack1

The nonreceptor tyrosine kinase (NRTK) Ack1 comprises a distinct arrangement of non-catalytic modules. Its SH3 domain has a C-terminal to the kinase domain (SH1), in contrast to the typical SH3-SH2-SH1 layout in NRTKs. The Ack1 is the only protein that shares a region of high homology to the tumor suppressor protein Mig6, a modulator of EGFR. The vertebrate Acks make up the only tyrosine kinase (TK) family known to carry a UBA domain. The GTPase binding and SAM domains are also uncommon in the NRTKs. In addition to being a downstream effector of receptor tyrosine kinases (RTKs) and integrins, Ack1 can act as an epigenetic regulator, modulate the degradation of the epidermal growth factor receptor (EGFR), confer drug resistance, and mediate the progression of hormone-sensitive tumors. In this review, we discuss the domain architecture of Ack1 in relation to other protein kinases that possess such defined regulatory domains.

Molecular dynamics simulations reveal the impact of NUDT15 R139C and R139H variants in structural conformation and dynamics

NUDT15, also known as MTH2, is a member of the NUDIX protein family that catalyzes the hydrolysis of nucleotides and deoxynucleotides, as well as thioguanine analogues. NUDT15 has been reported as a DNA sanitizer in humans, and more recent studies have shown that some genetic variants are related to a poor prognosis in neoplastic and immunologic diseases treated with thioguanine drugs. Despite this, the role of NUDT15 in physiology and molecular biology is quite unclear, as is the mechanism of action of this enzyme. The existence of clinically relevant variants has prompted the study of these enzymes, whose capacity to bind and hydrolyze thioguanine nucleotides is still poorly understood. By using a combination of biomolecular modeling techniques and molecular dynamics, we have studied the monomeric wild type NUDT15 as well as two important variants, R139C and R139H. Our findings reveal not only how nucleotide binding stabilizes the enzyme but also how two loops are responsible for keeping the enzyme in a packed, close conformation. Mutations in α2 helix affect a network of hydrophobic and π-interactions that enclose the active site. This knowledge contributes to the understanding of NUDT15 structural dynamics and will be valuable for the design of new chemical probes and drugs targeting this protein.Communicated by Ramaswamy H. Sarma.

Predicting the functional effect of compound heterozygous genotypes from large scale variant effect maps

Background

Pathogenic variants in PHGDH, PSAT1 , and PSPH cause a set of rare, autosomal recessive diseases known as serine biosynthesis defects. Serine biosynthesis defects present in a broad phenotypic spectrum that includes, at the severe end, Neu-Laxova syndrome, a lethal multiple congenital anomaly disease, intermediately in the form of infantile serine biosynthesis defects with severe neurological manifestations and growth deficiency, and at the mild end, as childhood disease with intellectual disability. However, because L-serine supplementation, especially if started early, can ameliorate and in some cases even prevent symptoms, knowledge of pathogenic variants is highly actionable.

Methods

Recently, our laboratory established a yeast-based assay for human PSAT1 function. We have now applied it at scale to assay the functional impact of 1,914 SNV-accessible amino acid substitutions. In addition to assaying the functional impact of individual variants in yeast haploid cells, we can assay pairwise combinations of PSAT1 alleles that recapitulate human genotypes, including compound heterozygotes, in yeast diploids.

Results

Results of our assays of individual variants (in haploid yeast cells) agree well with clinical interpretations and protein structure-function relationships, supporting the use of our data as functional evidence under the ACMG interpretation guidelines. Results from our diploid assay successfully distinguish patient genotypes from those of healthy carriers and agree well with disease severity. Finally, we present a linear model that uses individual allele measurements (in haploid yeast cells) to accurately predict the biallelic function (in diploid yeast cells) of ~ 1.8 million allele combinations corresponding to potential human genotypes.

Conclusions

Taken together, our work provides an example of how large-scale functional assays in model systems can be powerfully applied to the study of a rare disease.

RLBind: a deep learning method to predict RNA-ligand binding sites

Identification of RNA-small molecule binding sites plays an essential role in RNA-targeted drug discovery and development. These small molecules are expected to be leading compounds to guide the development of new types of RNA-targeted therapeutics compared with regular therapeutics targeting proteins. RNAs can provide many potential drug targets with diverse structures and functions. However, up to now, only a few methods have been proposed. Predicting RNA-small molecule binding sites still remains a big challenge. New computational model is required to better extract the features and predict RNA-small molecule binding sites more accurately. In this paper, a deep learning model, RLBind, was proposed to predict RNA-small molecule binding sites from sequence-dependent and structure-dependent properties by combining global RNA sequence channel and local neighbor nucleotides channel. To our best knowledge, this research was the first to develop a convolutional neural network for RNA-small molecule binding sites prediction. Furthermore, RLBind also can be used as a potential tool when the RNA experimental tertiary structure is not available. The experimental results show that RLBind outperforms other state-of-the-art methods in predicting binding sites. Therefore, our study demonstrates that the combination of global information for full-length sequences and local information for limited local neighbor nucleotides in RNAs can improve the model's predictive performance for binding sites prediction. All datasets and resource codes are available at https://github.com/KailiWang1/RLBind.

Targeting RNA:protein interactions with an integrative approach leads to the identification of potent YBX1 inhibitors

RNA-protein interactions (RPIs) are promising targets for developing new molecules of therapeutic interest. Nevertheless, challenges arise from the lack of methods and feedback between computational and experimental techniques during the drug discovery process. Here, we tackle these challenges by developing a drug screening approach that integrates chemical, structural and cellular data from both advanced computational techniques and a method to score RPIs in cells for the development of small RPI inhibitors; and we demonstrate its robustness by targeting Y-box binding protein 1 (YB-1), a messenger RNA-binding protein involved in cancer progression and resistance to chemotherapy. This approach led to the identification of 22 hits validated by molecular dynamics (MD) simulations and nuclear magnetic resonance (NMR) spectroscopy of which 11 were found to significantly interfere with the binding of messenger RNA (mRNA) to YB-1 in cells. One of our leads is an FDA-approved poly(ADP-ribose) polymerase 1 (PARP-1) inhibitor. This work shows the potential of our integrative approach and paves the way for the rational development of RPI inhibitors.

Molecular dynamics simulations of the calmodulin-induced α-helix in the SK2 calcium-gated potassium ion channel

The family of small-conductance Ca²⁺-activated potassium ion channels (SK channels) is composed of four members (SK1, SK2, SK3, and SK4) involved in neuron-firing regulation. The gating of these channels depends on the intracellular Ca²⁺ concentration, and their sensitivity to this ion is provided by calmodulin (CaM). This protein binds to a specific region in SK channels known as the calmodulin-binding domain (CaMBD), an event which is essential for their gating. While CaMBDs are typically disordered in the absence of CaM, the SK2 channel subtype displays a small prefolded α-helical region in its CaMBD even if CaM is not present. This small helix is known to turn into a full α-helix upon CaM binding, although the molecular-level details for this conversion are not fully understood yet. In this work, we offer new insights on this physiologically relevant process by means of enhanced sampling, atomistic Hamiltonian replica exchange molecular dynamics simulations, providing a more detailed understanding of CaM binding to this target. Our results show that CaM is necessary for inducing a full α-helix along the SK2 CaMBD through hydrophobic interactions with V426 and L427. However, it is also necessary that W431 does not compete for these interactions; the role of the small prefolded α-helix in the SK2 CaMBD would be to stabilize W431 so that this is the case. In conclusion, our findings provide further insight into a key interaction between CaM and SK channels that is important for channel sensitivity to Ca²⁺.

Structure of a TOC-TIC supercomplex spanning two chloroplast envelope membranes

The TOC and TIC complexes are essential translocons that facilitate the import of the nuclear genome-encoded preproteins across the two envelope membranes of chloroplast, but their exact molecular identities and assembly remain unclear. Here, we report a cryoelectron microscopy structure of TOC-TIC supercomplex from Chlamydomonas, containing a total of 14 identified components. The preprotein-conducting pore of TOC is a hybrid β-barrel co-assembled by Toc120 and Toc75, while the potential translocation path of TIC is formed by transmembrane helices from Tic20 and YlmG, rather than a classic model of Tic110. A rigid intermembrane space (IMS) scaffold bridges two chloroplast membranes, and a large hydrophilic cleft on the IMS scaffold connects TOC and TIC, forming a pathway for preprotein translocation. Our study provides structural insights into the TOC-TIC supercomplex composition, assembly, and preprotein translocation mechanism, and lays a foundation to interpret the evolutionary conservation and diversity of this fundamental translocon machinery.

rs12329760 Polymorphism in Transmembrane Serine Protease 2 Gene and Risk of Coronavirus Disease 2019 Mortality

The protease produced by the transmembrane serine protease 2 (TMPRSS2) gene enhances viral infections and has been linked to severe acute respiratory syndrome coronavirus 2 pathogenesis. Therefore, this study evaluated the association between TMPRSS2 and coronavirus disease 2019 (COVID-19) mortality. TMPRSS2 rs12329760 polymorphism was genotyped using the tetraprimer amplification refractory mutation system-polymerase chain reaction method in 592 dead and 693 improved patients. In the current study, the frequency of TMPRSS2 rs12329760 CC than TT genotypes was significantly lower in improved patients than in dead patients. According to the findings of the multivariate logistic regression test, higher levels of mean age, creatinine, erythrocyte sedimentation rate, C-reactive protein, aspartate aminotransferase, lower levels of 25-hydroxyvitamin D, uric acid, and real-time PCR Ct values and TMPRSS2 rs12329760 CC genotype were observed to be associated with increased COVID-19 mortality rates. In conclusion, the TMPRSS2 rs12329760 CC genotype was a polymorphism linked to a significantly higher incidence of severe COVID-19. Further studies are required to corroborate the obtained findings.

Mitochondrial DNA Changes in Genes of Respiratory Complexes III, IV and V Could Be Related to Brain Tumours in Humans

Mitochondrial DNA changes can contribute to both an increased and decreased likelihood of cancer. This process is complex and not fully understood. Polymorphisms and mutations, especially those of the missense type, can affect mitochondrial functions, particularly if the conservative domain of the protein is concerned. This study aimed to identify the possible relationships between brain gliomas and the occurrence of specific mitochondrial DNA polymorphisms and mutations in respiratory complexes III, IV and V. The investigated material included blood and tumour material collected from 30 Caucasian patients diagnosed with WHO grade II, III or IV glioma. The mitochondrial genetic variants were investigated across the mitochondrial genome using next-generation sequencing (MiSeq/FGx system—Illumina). The study investigated, in silico, the effects of missense mutations on the biochemical properties, structure and functioning of the encoded protein, as well as their potential harmfulness. The A14793G (MTCYB), A15758G, (MT-CYB), A15218G (MT-CYB), G7444A (MT-CO1) polymorphisms, and the T15663C (MT-CYB) and G8959A (ATP6) mutations were assessed in silico as harmful alterations that could be involved in oncogenesis. The G8959A (E145K) ATP6 missense mutation has not been described in the literature so far. In light of these results, further research into the role of mtDNA changes in brain tumours should be conducted.

Screening of Pathogenic Missense Single Nucleotide Variants From LHPP Gene Associated With the Hepatocellular Carcinoma: An In silico Approach

LHPP gene encodes a phospholysine phosphohistidine inorganic pyrophosphate phosphatase, which functions as a tumor-suppressor protein. The tumor suppression by this protein has been confirmed in various cancers, including hepatocellular carcinoma (HCC). LHPP downregulation promotes cell growth and proliferation by modulating the PI3K/AKT signaling pathway. This study identifies potentially deleterious missense single nucleotide variants (SNVs) associated with the LHPP gene using multiple computational tools based on different algorithms. A total of 4 destabilizing mutants are identified as L22P, I212T, G227R, and G236R, from the conserved region of the phosphatase. The 3-dimensional (3D) modeling and structural comparison of variants with the native protein reveals significant structural and conformational variations after mutations, suggesting disruption in the function of phospholysine phosphohistidine inorganic pyrophosphate phosphatase. The identified mutations might, therefore, participate in the cause of HCC.

The Inhibitors of CDK4/6 from a Library of Marine Compound Database: A Pharmacophore, ADMET, Molecular Docking and Molecular Dynamics Study

BACKGROUND

CDK4/6 (Cyclin-dependent kinases 4/6) are the key promoters of cell cycle transition from G1 phase to S phase. Thus, selective inhibition of CDK4/6 is a promising cancer treatment.

METHODS

A total of 52,765 marine natural products were screened for CDK4/6. To screen out better natural compounds, pharmacophore models were first generated, then the absorption, distribution, metabolism, elimination, and toxicity (ADMET) were tested, followed by molecular docking. Finally, molecular dynamics simulation was carried out to verify the binding characteristics of the selected compounds.

RESULTS

Eighty-seven marine small molecules were screened based on the pharmacophore model. Then, compounds 41369 and 50843 were selected according to the ADMET and molecular docking score for further kinetic simulation evaluation. Finally, through molecular dynamics analysis, it was confirmed that compound 50843 maintained a stable conformation with the target protein, so it has the opportunity to become an inhibitor of CDK4/6.

CONCLUSION

Through structure-based pharmacophore modeling, ADMET, the molecular docking method and molecular dynamics (MD) simulation, marine natural compound 50843 was proposed as a promising marine inhibitor of CDK4/6.

Association of TMPRSS2 Gene Polymorphisms with COVID-19 Severity and Mortality: a Case-Control Study with Computational Analyses

Coronavirus disease 2019 (COVID-19) is a severe disease caused by a new variant of beta-coronavirus that first appeared in China. Human genetic factors, including polymorphisms, serve pivotal roles in the high transmission of SARS-CoV-2 and the stubbornly progressing sickness seen in a small but significant percentage of infected people; however, but these factors remain ill-defined. A total of 288 COVID-19 patients and 288 controls were genotyped for TMPRSS2 polymorphisms using both restriction fragment length polymorphism polymerase chain reaction (RFLP-PCR) and amplification refractory mutation system (ARMS)-PCR techniques. Different genotypes of TMPRSS2 polymorphisms were compared in terms of disease susceptibility and mortality. The statistical analysis showed that minor alleles of all studied variants statistically increased the risk of COVID-19, except for the rs75603675 C > A variant. The T allele of rs12329760 conferred an increased risk of COVID-19. Moreover, the AG/AC/TT/AG combination of genotypes significantly enhanced the risk of COVID-19 in our population. Different haplotypes of rs17854725/rs75603675/rs12329760/rs4303795 polymorphisms, including GACA, GACG, GATG, GATA, AATA, ACCG, ACTG, ACTA, GCCA, and GCTG, were found to be associated with increased risk of the disease (odds ratio > 1). Regarding the clinical and paraclinical characteristics, a statistically significant difference was found between non-severe and severe forms except for gender, platelet, C-reactive protein (CRP), erythrocyte sedimentation rate (ESR), and underlying diseases. In addition, case genotypes of TMPRSS2 rs17854725 A > G, rs12329760 C > T, and rs4303795 A > G were significantly different regarding severe and non-severe forms of the disease (P-value < 0.001). Specifically, death was more frequent in carriers of the AG genotype of rs17854725 A > G (P-value = 0.022). Patients who carry the minor alleles of the four studied TMPRSS2 variants were rather vulnerable to COVID-19 infection. Our findings indicated that rs17854725 A > G (AA vs. AG and AA vs. GG), rs12329760 C > T (CC vs. CT and CC vs. TT), and rs4303795 A > G (AA vs. AG) genotypes of TMPRSS2 variations are associated with a more invasive disorder pattern. More studies on larger populations are needed to confirm our results.

Dynamic coupling of residues within proteins as a mechanistic foundation of many enigmatic pathogenic missense variants

Many pathogenic missense mutations are found in protein positions that are neither well-conserved nor fall in any known functional domains. Consequently, we lack any mechanistic underpinning of dysfunction caused by such mutations. We explored the disruption of allosteric dynamic coupling between these positions and the known functional sites as a possible mechanism for pathogenesis. In this study, we present an analysis of 591 pathogenic missense variants in 144 human enzymes that suggests that allosteric dynamic coupling of mutated positions with known active sites is a plausible biophysical mechanism and evidence of their functional importance. We illustrate this mechanism in a case study of β-Glucocerebrosidase (GCase) in which a vast majority of 94 sites harboring Gaucher disease-associated missense variants are located some distance away from the active site. An analysis of the conformational dynamics of GCase suggests that mutations on these distal sites cause changes in the flexibility of active site residues despite their distance, indicating a dynamic communication network throughout the protein. The disruption of the long-distance dynamic coupling caused by missense mutations may provide a plausible general mechanistic explanation for biological dysfunction and disease.

Structural and Biochemical Characterization of Endo-β-1,4-glucanase from Dictyoglomus thermophilum, a Hyperthermostable and Halotolerant Cellulase

Enzymatic conversion of polysaccharides in the lignocellulosic biomass is currently the subject of intensive research and will be a key technology in future biorefineries. Using a bioinformatics approach, we previously identified a putative endo-β-1,4-glucanase (DtCel5A) from Dictyoglomus thermophilum, a chemoorganotrophic and thermophilic bacterium. Here, we structurally and functionally characterize DtCel5A and show that it is endowed with remarkable thermal and chemical stability. The structural features of DtCel5A and of its complex with cellobiose have been investigated by combining X-ray crystallography and other biophysical studies. Importantly, biochemical assays show that DtCel5A retains its activity on cellulose at high temperatures and at elevated salt concentrations. These features make DtCel5A an enzyme with interesting biotechnological applications for biomass degradation.

Structural dynamics in the evolution of a bilobed protein scaffold

Proteins conduct numerous complex biological functions by use of tailored structural dynamics. The molecular details of how these emerged from ancestral peptides remains mysterious. How does nature utilize the same repertoire of folds to diversify function? To shed light on this, we analyzed bilobed proteins with a common structural core, which is spread throughout the tree of life and is involved in diverse biological functions such as transcription, enzymatic catalysis, membrane transport, and signaling. We show here that the structural dynamics of the structural core differentiate predominantly via terminal additions during a long-period evolution. This diversifies substrate specificity and, ultimately, biological function.

Novel biophysical tools allow the structural dynamics of proteins and the regulation of such dynamics by binding partners to be explored in unprecedented detail. Although this has provided critical insights into protein function, the means by which structural dynamics direct protein evolution remain poorly understood. Here, we investigated how proteins with a bilobed structure, composed of two related domains from the periplasmic-binding protein–like II domain family, have undergone divergent evolution, leading to adaptation of their structural dynamics. We performed a structural analysis on ∼600 bilobed proteins with a common primordial structural core, which we complemented with biophysical studies to explore the structural dynamics of selected examples by single-molecule Förster resonance energy transfer and Hydrogen–Deuterium exchange mass spectrometry. We show that evolutionary modifications of the structural core, largely at its termini, enable distinct structural dynamics, allowing the diversification of these proteins into transcription factors, enzymes, and extracytoplasmic transport-related proteins. Structural embellishments of the core created interdomain interactions that stabilized structural states, reshaping the active site geometry, and ultimately altered substrate specificity. Our findings reveal an as-yet-unrecognized mechanism for the emergence of functional promiscuity during long periods of evolution and are applicable to a large number of domain architectures.

Structure of a dimer of the Sulfolobus solfataricus MCM N-terminal domain reveals a potential role in MCM ring opening

Cells strongly regulate DNA replication to ensure genomic stability and prevent several diseases, including cancers. Eukaryotes and archaea strictly control DNA-replication initiation by the regulated loading of hexameric minichromosome maintenance (MCM) rings to encircle both strands of the DNA double helix followed by regulated activation of the loaded rings such that they then encircle one DNA strand while excluding the other. Both steps involve an open/closed ring transformation, allowing DNA strands to enter or exit. Here, the crystal structure of a dimer of the N-terminal domain of Sulfolobus solfataricus MCM with an intersubunit interface that is more extensive than in closed-ring structures, while including common interactions to enable facile interconversion, is presented. It is shown that the identified interface could stabilize open MCM rings by compensating for lost interactions at an open neighbor interface and that the prior open-ring cryo-EM structure of MCM loading has a similar extended interface adjacent to its open interface.

Quantifying evolutionary importance of protein sites: A Tale of two measures

A key challenge in evolutionary biology is the accurate quantification of selective pressure on proteins and other biological macromolecules at single-site resolution. The evolutionary importance of a protein site under purifying selection is typically measured by the degree of conservation of the protein site itself. A possible alternative measure is the strength of the site-induced conservation gradient in the rest of the protein structure. However, the quantitative relationship between these two measures remains unknown. Here, we show that despite major differences, there is a strong linear relationship between the two measures such that more conserved protein sites also induce stronger conservation gradient in the rest of the protein. This linear relationship is universal as it holds for different types of proteins and functional sites in proteins. Our results show that the strong selective pressure acting on the functional site in general percolates through the rest of the protein via residue-residue contacts. Surprisingly however, catalytic sites in enzymes are the principal exception to this rule. Catalytic sites induce significantly stronger conservation gradients in the rest of the protein than expected from the degree of conservation of the site alone. The unique requirement for the active site to selectively stabilize the transition state of the catalyzed chemical reaction imposes additional selective constraints on the rest of the enzyme.

Sites within proteins which are important for stability or function are under stronger selective pressure and evolve more slowly than other sites. Catalytic sites in enzymes are such highly conserved sites with relatively low evolutionary rates. Recently, catalytic sites were shown to induce a strong gradient of conservation such that the closer a residue is to the catalytic site, the more conserved it is. Here we show that there is a universal linear relationship between the degree of evolutionary conservation of a protein site and the conservation gradient it induces in the protein tertiary structure, applicable to all types of sites. Our findings suggest that selective pressure acting on a protein site generally percolates through the rest of the protein via residue-residue contacts. Remarkably however, catalytic sites induce significantly stronger conservation gradients than expected from their degree of conservation alone. Our results indicate that the strong conservation gradient induced by catalytic sites is driven by the unique function of enzyme catalysis, which requires the participation of many residues beyond the few key catalytic residues. Our results provide insights into evolutionary conservation patterns of and surrounding proteins functional sites, with implications for functional site prediction and protein design.

Simplified geometric representations of protein structures identify complementary interaction interfaces

Protein-protein interactions are critical to protein function, but three-dimensional (3D) arrangements of interacting proteins have proven hard to predict, even given the identities and 3D structures of the interacting partners. Specifically, identifying the relevant pairwise interaction surfaces remains difficult, often relying on shape complementarity with molecular docking while accounting for molecular motions to optimize rigid 3D translations and rotations. However, such approaches can be computationally expensive, and faster, less accurate approximations may prove useful for large-scale prediction and assembly of 3D structures of multi-protein complexes. We asked if a reduced representation of protein geometry retains enough information about molecular properties to predict pairwise protein interaction interfaces that are tolerant of limited structural rearrangements. Here, we describe a reduced representation of 3D protein accessible surfaces on which molecular properties such as charge, hydrophobicity, and evolutionary rate can be easily mapped, implemented in the MorphProt package. Pairs of surfaces are compared to rapidly assess partner-specific potential surface complementarity. On two available benchmarks of 185 overall known protein complexes, we observe predictions comparable to other structure-based tools at correctly identifying protein interaction surfaces. Furthermore, we examined the effect of molecular motion through normal mode simulation on a benchmark receptor-ligand pair and observed no marked loss of predictive accuracy for distortions of up to 6 Å Cα-RMSD. Thus, a shape reduction of protein surfaces retains considerable information about surface complementarity, offers enhanced speed of comparison relative to more complex geometric representations, and exhibits tolerance to conformational changes.

Proteostasis of Islet Amyloid Polypeptide: A Molecular Perspective of Risk Factors and Protective Strategies for Type II Diabetes

The possible link between hIAPP accumulation and β-cell death in diabetic patients has inspired numerous studies focusing on amyloid structures and aggregation pathways of this hormone. Recent studies have reported on the importance of early oligomeric intermediates, the many roles of their interactions with lipid membrane, pH, insulin, and zinc on the mechanism of aggregation of hIAPP. The challenges posed by the transient nature of amyloid oligomers, their structural heterogeneity, and the complex nature of their interaction with lipid membranes have resulted in the development of a wide range of biophysical and chemical approaches to characterize the aggregation process. While the cellular processes and factors activating hIAPP-mediated cytotoxicity are still not clear, it has recently been suggested that its impaired turnover and cellular processing by proteasome and autophagy may contribute significantly toward toxic hIAPP accumulation and, eventually, β-cell death. Therefore, studies focusing on the restoration of hIAPP proteostasis may represent a promising arena for the design of effective therapies. In this review we discuss the current knowledge of the structures and pathology associated with hIAPP self-assembly and point out the opportunities for therapy that a detailed biochemical, biophysical, and cellular understanding of its aggregation may unveil.