Protein Stability Analysis Research Articles

The Effects of SARS CoV-2 nsp13 Mutations on the Structure and Stability of Helicase in Chinese Isolates

Objective: Coronavirus Disease 2019 (COVID19) is a viral disease caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS CoV-2). The high mutation propensity of the SARS CoV-2 genome is one of the biggest threats to the long-term validity of treatment options. Helicases are anti-viral targets because of the vital role they play in the viral life cycle. In this study, changes in the protein structure caused by SARS CoV-2 nsp13 mutations were investigated to contribute to the development of effective antiviral drugs. Materials and Methods: Genome data of 298 individuals located in the China location were examined. The mutant model was built using deep learning algorithms. Model quality assessment was done with QMEAN. Protein stability analyses were performed with DynaMut2 and Cutoff Scanning Matrix stability. Changes in substrate affinity were performed with Haddock v2.4. Results: In this study, twenty-eight mutations in nsp13 were identified (23 sense, 5 missense). The changes in protein structure caused by the five missense mutations (Leu14Phe, Arg15Ser, Arg21Ser, Leu235Phe, Ala454Thr) were modeled. The mutations caused a decrease in the stability of SARS CoV-2 helicase (-0.99, -1.66, -1.15, -0.54, and -0.73 for Leu14Phe, Arg15Ser, Arg21Ser, Leu235Phe, Ala454Thr, respectively). The mutations reduced the helicase's affinity to the substrate. The docking scores for wild-type and mutant helicase were -84.4±1.4 kcal.mol-1 and -71.1±6.7 kcal.mol-1, respectively. Conclusion: Helicase mutations caused a decrease in the protein stability and nucleic acid affinity of the SARS CoV-2 helicase. The results provide important data on the development of potential antivirals and the effect of mutation on the functions of viral proteins.

In Silico Analysis Identified Putative Pathogenic Missense nsSNPs in Human SLITRK1 Gene.

Human DNA contains several variations, which can affect the structure and normal functioning of a protein. These variations could be single nucleotide polymorphisms (SNPs) or insertion-deletions (InDels). SNPs, as opposed to InDels, are more commonly present in DNA and may cause genetic disorders. In the current study, several bioinformatic tools were used to prioritize the pathogenic variants in the SLITRK1 gene. Out of all of the variants, 16 were commonly predicted to be pathogenic by these tools. All the variants had very low frequency, i.e., <0.0001 in the global population. The secondary structure of all filtered variants was predicted, but no structural change was observed at the site of variation in any variant. Protein stability analysis of these variants was then performed, which determined a decrease in protein stability of 10 of the variants. Amino acid conservation analysis revealed that all the amino acids were highly conserved, indicating their structural and functional importance. Protein 3D structure of wildtype SLITRK1 and all of its variants was predicted using I-TASSER, and the effect of variation on 3D structure of the protein was observed using the Missense3D tool, which presented the probable structural loss in three variants, i.e., Asn529Lys, Leu496Pro and Leu94Phe. The wildtype SLITRK1 protein and these three variants were independently docked with their close interactor protein PTPRD, and remarkable differences were observed in the docking sites of normal and variants, which will ultimately affect the functional activity of the SLITRK1 protein. Previous studies have shown that mutations in SLITRK1 are involved in Tourette syndrome. The present study may assist a molecular geneticist in interpreting the variant pathogenicity in research as well as diagnostic setup.

Landscape-Based Protein Stability Analysis and Network Modeling of Multiple Conformational States of the SARS-CoV-2 Spike D614G Mutant: Conformational Plasticity and Frustration-Induced Allostery as Energetic Drivers of Highly Transmissible Spike Variants.

The structural and functional studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across many variants of concern (VOCs), suggesting the effect of this mutation on the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies provided important evidence about multiple conformational substates of the D614G spike protein. The development of a plausible mechanistic model that can explain the experimental observations from a more unified thermodynamic perspective is an important objective of the current work. In this study, we employed efficient and accurate coarse-grained simulations of multiple structural substates of the D614G spike trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration profiling of the conformational states with residue-based mutational scanning of protein stability and network analysis of allosteric interactions and communications, we determine the patterns of mutational sensitivity in the functional regions and sites of variants. We found that the D614G mutation may induce a considerable conformational adaptability of the open states in the SARS-CoV-2 spike protein without compromising the folding stability and integrity of the spike protein. The results suggest that the D614G mutant may employ a hinge-shift mechanism in which the dynamic couplings between the site of mutation and the interprotomer hinge modulate the interdomain interactions, global mobility change, and the increased stability of the open form. This study proposes that mutation-induced modulation of the conformational flexibility and energetic frustration at the interprotomer interfaces may serve as an efficient mechanism for allosteric regulation of the SARS-CoV-2 spike proteins.

ProSAP: a GUI software tool for statistical analysis and assessment of thermal stability data.

The Cellular Thermal Shift Assay (CETSA) plays an important role in drug-target identification, and statistical analysis is a crucial step significantly affecting conclusion. We put forward ProSAP (Protein Stability Analysis Pod), an open-source, cross-platform and user-friendly software tool, which provides multiple methods for thermal proteome profiling (TPP) analysis, nonparametric analysis (NPA), proteome integral solubility alteration and isothermal shift assay (iTSA). For testing the performance of ProSAP, we processed several datasets and compare the performance of different algorithms. Overall, TPP analysis is more accurate with fewer false positive targets, but NPA methods are flexible and free from parameters. For iTSA, edgeR and DESeq2 identify more true targets than t-test and Limma, but when it comes to ranking, the four methods show not much difference. ProSAP software is available at https://github.com/hcji/ProSAP and https://zenodo.org/record/5763315.

Predicting pathogenicity for novel hearing loss mutations based on genetic and protein structure approaches

Hearing loss is a heterogeneous disorder. Identification of causative mutations is demanding due to genetic heterogeneity. In this study, we investigated the genetic cause of sensorineural hearing loss in patients with severe/profound deafness. After the exclusion of GJB2-GJB6 mutations, we performed whole exome sequencing in 32 unrelated Argentinean families. Mutations were detected in 16 known deafness genes in 20 patients: ACTG1, ADGRV1 (GPR98), CDH23, COL4A3, COL4A5, DFNA5 (GSDDE), EYA4, LARS2, LOXHD1, MITF, MYO6, MYO7A, TECTA, TMPRSS3, USH2A and WSF1. Notably, 11 variants affecting 9 different non-GJB2 genes resulted novel: c.12829C > T, p.(Arg4277*) in ADGRV1; c.337del, p.(Asp109*) and c.3352del, p.(Gly1118Alafs*7) in CDH23; c.3500G > A, p.(Gly1167Glu) in COL4A3; c.1183C > T, p.(Pro395Ser) and c.1759C > T, p.(Pro587Ser) in COL4A5; c.580 + 2 T > C in EYA4; c.1481dup, p.(Leu495Profs*31) in LARS2; c.1939 T > C, p.(Phe647Leu), in MYO6; c.733C > T, p.(Gln245*) in MYO7A and c.242C > G, p.(Ser81*) in TMPRSS3 genes. To predict the effect of these variants, novel protein modeling and protein stability analysis were employed. These results highlight the value of whole exome sequencing to identify candidate variants, as well as bioinformatic strategies to infer their pathogenicity.

Investigation of changes in protein stability and substrate affinity of 3CL-protease of SARS-CoV-2 caused by mutations.

3CLpro of SARS-CoV-2 is one of the enzymes required for the replication process of the virus responsible for the COVID-19 pandemic. In this study, changes in protein stability and substrate affinity caused by mutations were investigated to stir the development of potent inhibitors. Sequence data of samples were obtained from the NCBI Virus database. Mutation analyses were performed with RDP4 and MegaX. 3CLpro tertiary models were created using Robetta. Molecular docking for peptidomimetic substrate and inhibitor ligand was done with Autodock v4.2 and Haddock v2.4. Protein stability analysis was performed using mCSM stability and DynaMut2. Twenty-four missense mutations in 3CLpro were identified in this study. Changes in the 3CLpro structure induced by the mutations Met49Thr, Leu167Ser, and Val202Ala resulted in significant levels of instability (-2.029,-2.612,-2.177 kcal.mol-1, respectively). The lowest interaction energy for substrate was -58.7 kcal.mol-1 and -62.6 kcal.mol-1 in wild-type and mutant, respectively. The lowest docking energy for ligand was -6.19 and -9.52 kcal.mol-1 for wild-type and mutant, respectively. This study reports for the first time that mutations cause increased substrate affinity of 3CLpro from SARS-CoV-2. This research provides important data for the development of potent peptidomimetic inhibitors for the treatment of COVID-19.

Quantitative analysis of biomolecule release from polystyrene-block-polyethylene oxide thin films.

Block copolymers have garnered recent attention due to their ability to contain molecular cargo within nanoscale domains and release said cargo in aqueous environments. However, the release kinetics of cargo from these thin-films has not yet been reported. Knowledge of the release quantities and release profiles of these systems is paramount for applications of these systems. Here, Polystyrene-block-poly(ethylene oxide) (PS-b-PEO) was co-assembled with fluorescein isothiocyanate isomer I-lysozyme (FITC-LZ) and fluorescein isothiocyanate isomer I-TAT (FITC-TAT), such that these molecular cargos arrange within the PEO domains of the thin films. We show that high loading ratios of cargo/PS-b-PEO do not significantly impact the nanostructure of the films; however, a loading limit appears to be present with aggregates of protein forming at the microscale with higher loading ratios. The presence of lysozyme (LZ) within the films was confirmed qualitatively after aqueous exposure through photo-induced force microscopy (PiFM) imaging at the Amide I characteristic peak (∼1650 cm-1). Furthermore, we demonstrate that LZ maintains activity and structure after exposure to the polymer solvent (benzene/methanol/water mix). Finally, we demonstrate quantitatively 20-80 ng cm-2 of cargo is released from these films, depending on the cargo incorporated. We show that the larger molecule lysozyme is released over a longer time than the smaller TAT peptide. Finally, we demonstrate the ability to tune the quantity of cargo released by altering the thickness of the PS-b-PEO thin-films during fabrication.

Effect of different plants’ aromatic essential oils on frozen Awassi lamb meat’s chemical and physical characteristics

The effect of drenching Awassi lambs with three aromatic essential oils from sage (Salvia officinalis L.), clove (Syzygium aromaticum L.), and laurel (Laurus nobilis L.) was investigated on meat chemical and physical characteristics, and oxidative and deterioration measurements. Twenty-four Awassi lambs, five to six months old, were divided into four groups. A concentrated diet was provided to the lambs at a rate of 3% of the body weight. The treatments were as follows: T1 was served as the untreated control, while T2, T3, and T4 were drenched with oils of sage, clove, and laurel, respectively. Drenching was carried out using water-soluble capsules containing 500 mg oil/capsule/day. Treatments lasted 90 days. At the end of the treatment period, the animals were fasted overnight and slaughtered. The carcasses were cleaned and kept at 4°C for 24 h. The longissimus dorsi (LD) muscle was then separated and preserved in a plastic bag for three preservation periods: no freezing and 30 days and 60 days freezing at −18°C. Several physical, fat, and protein stability analyses of meat were done after the preservation periods. The results indicated no significant effect of drenching Awassi lambs with different aromatic essential oils on the meat’s physical and chemical characteristics. However, these oils, especially clove oil, affected fat and protein stability with increasing preservation period by freezing.

Functional and structural analysis of non-synonymous single nucleotide polymorphisms (nsSNPs) in the MYB oncoproteins associated with human cancer

MYB proteins are highly conserved DNA-binding domains (DBD) and mutations in MYB oncoproteins have been reported to cause aberrant and augmented cancer progression. Identification of MYB molecular biomarkers predictive of cancer progression can be used for improving cancer management. To address this, a biomarker discovery pipeline was employed in investigating deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in predicting damaging and potential alterations on the properties of proteins. The nsSNP of the MYB family; MYB, MYBL1, and MYBL2 was extracted from the NCBI database. Five in silico tools (PROVEAN, SIFT, PolyPhen-2, SNPs&GO and PhD-SNP) were utilized to investigate the outcomes of nsSNPs. A total of 45 nsSNPs were predicted as high-risk and damaging, and were subjected to PMut and I-Mutant 2.0 for protein stability analysis. This resulted in 32 nsSNPs with decreased stability with a DDG score lower than − 0.5, indicating damaging effect. G111S, N183S, G122S, and S178C located within the helix-turn-helix (HTH) domain were predicted to be conserved, further posttranslational modifications and 3-D protein analysis indicated these nsSNPs to shift DNA-binding specificity of the protein thus altering the protein function. Findings from this study would help in the field of pharmacogenomic and cancer therapy towards better intervention and management of cancer.

Whole Exome Sequencing Identifies a Novel Homozygous Missense Mutation in the CSB Protein-Encoding ERCC6 Gene in a Taiwanese Boy with Cockayne Syndrome.

Background: Cockayne syndrome (CS) is a rare form of dwarfism that is characterized by progressive premature aging. CS is typically caused by mutations in the excision repair cross-complementing protein group 6 (ERCC6) gene that encodes the CS group B (CSB) protein. Using whole exome sequencing, we recently identified a novel homozygous missense mutation (Leu536Trp) in CSB in a Taiwanese boy with CS. Since the current database (Varsome) interprets this variant as likely pathogenic, we utilized a bioinformatic tool to investigate the impact of Leu536Trp as well as two other variants (Arg453Ter, Asp532Gly) in similar articles on the CSB protein structure stability. Methods: We used iterative threading assembly refinement (I-TASSER) to generate a predictive 3D structure of CSB. We calculated the change of mutation energy after residues substitution on the protein stability using I-TASSER as well as the artificial intelligence program Alphafold. Results: The Asp532Gly variant destabilized both modeled structures, while the Leu536Trp variant showed no effect on I-TASSER’s model but destabilized the Alphafold’s modeled structure. Conclusions: We propose here the first case of CS associated with a novel homozygous missense mutation (Leu536Trp) in CSB. Furthermore, we suggest that the Asp532Gly and Leu536Trp variants are both pathogenic after bioinformatic analysis of protein stability.

Novel MYO1D Missense Variant Identified Through Whole Exome Sequencing and Computational Biology Analysis Expands the Spectrum of Causal Genes of Laterality Defects.

Laterality defects (LDs) or asymmetrically positioned organs are a group of rare developmental disorders caused by environmental and/or genetic factors. However, the exact molecular pathophysiology of LD is not yet fully characterised. In this context, studying Arab population presents an ideal opportunity to discover the novel molecular basis of diseases owing to the high rate of consanguinity and genetic disorders. Therefore, in the present study, we studied the molecular basis of LD in Arab patients, using next-generation sequencing method. We discovered an extremely rare novel missense variant in MYO1D gene (Pro765Ser) presenting with visceral heterotaxy and left isomerism with polysplenia syndrome. The proband in this index family has inherited this homozygous variant from her heterozygous parents following the autosomal recessive pattern. This is the first report to show MYO1D genetic variant causing left–right axis defects in humans, besides previous known evidence from zebrafish, frog and Drosophila models. Moreover, our multilevel bioinformatics-based structural (protein variant structural modelling, divergence, and stability) analysis has suggested that Ser765 causes minor structural drifts and stability changes, potentially affecting the biophysical and functional properties of MYO1D protein like calmodulin binding and microfilament motor activities. Functional bioinformatics analysis has shown that MYO1D is ubiquitously expressed across several human tissues and is reported to induce severe phenotypes in knockout mouse models. In conclusion, our findings show the expanded genetic spectrum of LD, which could potentially pave way for the novel drug target identification and development of personalised medicine for high-risk families.

Clinically significant findings of high-risk mutations in human SLC29A4 gene associated with diabetes mellitus type 2 in Pakistani population

This study conducted an in-depth analysis combining computational and experimental verifications of the deleterious missense mutations associated with the SLC29A4 protein. The functional annotation of the non-synonymous single nucleotide polymorphism (nsSNPs), followed by structure-function analysis, revealed 13 single nucleotide polymorphisms (SNP) as the most damaging. Among these, six mutants P429T/S, L144S, M108V, N86H, and V79E, were predicted as structurally and functionally damaging by protein stability analysis. Also, these variants are located at evolutionary conserved regions, either buried, contributing to the structural damage, or exposed, causing functional changes in the protein. These mutants were further taken for molecular docking studies. When verified via experimental analysis, the SNPs M108V (rs149798710), N86H (rs151039853), and V79E (rs17854505) showed an association with type 2 diabetes mellitus (T2DM). Minor allele frequency for rs149798710 (A > G) was 0.23 in controls, 0.29 in metformin responders, 0.37 in metformin non-responder, for rs151039853 (A > C) was 0.21 in controls, 0.28 in metformin responders, 0.36 in metformin non-responder and for rs17854505 (T > A) was 0.20 in controls, 0.25 in metformin responders, 0.37 in metformin non-responder. Hence, this study concludes that SLC29A4 M108V (rs149798710), N86H (rs151039853), and V79E (rs17854505) polymorphisms were associated with the increased risk of T2DM as well as with the increased risk towards the failure of metformin therapeutic response in T2DM patients of Pakistan. Communicated by Ramaswamy H. Sarma

The investigation of nonsynonymous SNPs of human SLC6A4 gene associated with depression: An in silico approach

Genetic polymorphism of the SLC6A4 gene is associated with several behavioral disorders, including depression. Since studying the total nonsynonymous single nucleotide polymorphisms (nsSNPs) of the SLC6A4 gene at the population level is a difficult task, we aim to utilize in silico approach to detect the most deleterious nsSNPs of the SLC6A4 gene. In our study, 7 computational tools were used in the initial stage, including SIFT, Polyphen-2, PROVEAN, SNAP2, PhD-SNP, PANTHER, and SNPs&GO to find out the most damaging nsSNPs. In the second phase, we performed structural, functional, and stability analysis of SLC6A4 protein by popular computation tools, including I-Mutant 2.0 and MutPred2. Also, the ConSurf server was utilized to find the conserved region of the SLC6A4 protein to determine the relationship between these conserved regions with high-risk nsSNPs. Based on these analyses, 5 high-risk mutations of the SLC6A4 protein were selected. Then, we carried out comparative modeling by using the Robetta server and aligned the mutant protein model with the native protein structure. Later, we performed the post-translational modification and functional domain analysis of the SLC6A4 protein. This study concludes that Arginine → Tryptophan at position 79 and Arginine → Cysteine at position 104 are the two significant mutations in SLC6A4 protein which might play an essential role in causing diseases. Future studies should take these high-risk nsSNPs (rs1221448303, rs200953188) into consideration while exploring diseases related to the SLC6A4 gene. Besides, our research is the first-ever comprehensive in silico investigation of the SLC6A4 gene. Thus, the findings of this study could be beneficial for developing precision medicines against diseases caused by SLC6A4 malfunction. Furthermore, extensive wet-lab research and experiments on various model organisms might be helpful to investigate the precise role of these damaging nsSNPs of the SLC6A4 gene.

Analysis of δ-globin gene alleles in Tunisians: description of three new delta-thalassemia mutations.

Thalassemia is one of the most prevalent worldwide autosomal recessive disorders characterized by a great molecular and clinical expression heterogeneity. Alpha and beta-thalassemia are the main two types observed in case of mutations affecting alpha and beta-globin genes respectively. Delta-thalassemia is noted when mutations occur on the delta-globin gene. In Tunisia, β-thalassemia prevalence isestimated at 2.21% of carriers. However, few reports investigated the delta-globin gene. In this work, we aimed to perform a molecular study to help define the molecular spectrum of δ-thalassemia mutations in Tunisia. The study involved 7558 patients among whom we selected 179 individuals with abnormal HbA2 values or fractions. Hemoglobin analysis was performed using Capillary electrophoresis (CE) and high-performance liquid chromatography (HPLC). DNA sequencing was performed on ABI prism 310 Genetic Analyzer Applied Biosystems. CUPSAT (Cologne University Protein Stability Analysis Tool) was used for the prediction of protein stability changes upon missense mutations and mutants were modeled via DeepView-SwissPdbViewer and POV-Ray softwares for molecular dynamics simulation studies. We identified four mutations: HbA2-Yialousa described for the first time in Tunisia ( in 72.72% of cases) and 3 mutations reported for the first time in the world: (i) c.442T > C Stop147Arg ext 15aa-stop observed in 18.18% of cases, (ii) c.187 G > C (Ala62Pro) noted in 4.54% of casesand (iii) c.93-1G > C found in 4.54% of cases. Our data provide genetic basis that would be especially useful in screening for beta-thalassemia trait during delta-beta thalassemia associations.

Lactucin induces apoptosis through reactive oxygen species-mediated BCL-2 and CFLARL downregulation in Caki-1 cells.

Lactucin, a naturally occurring active sesquiterpene lactone, is abundantly found in chicory and romaine lettuce. A recent study reported that lactucin could induce apoptosis in leukemia cells. However, its cytotoxicity and potential molecular mechanisms underlying cancer cell death remain unclear. Therefore, in this study, we aimed to investigate the direct effect and underlying mechanism of action of lactucin on renal cancer cells. MTT assay and flow cytometry were performed to evaluate the rate of cell proliferation and apoptosis, respectively. Western blotting, reverse transcription polymerase chain reaction, and protein stability analyses were performed to analyze the effect of lactucin on the expression of apoptosis-related proteins such as B-cell lymphoma 2 (BCL-2) and CFLAR (CASP8 and FADD like apoptosis regulator) long isoform (CFLARL) in Caki-1 human renal cancer cells. In addition, reactive oxygen species (ROS) generation was evaluated using flow cytometry. Lactucin treatment induced apoptosis in Caki-1 cells in a dose-dependent manner via activation of the caspase pathway. It downregulated BCL-2 and CFLARL expression levels by suppressing BCL-2 transcription and CFLARL protein stability, respectively. Pretreatment with N-acetyl-1-cysteine, a ROS scavenger, attenuated the lactucin-induced apoptosis and restored the BCL-2 and CFLARL expression to basal levels. Lactucin-facilitated BCL-2 downregulation was regulated at the transcriptional level through the inactivation of the NF-κB pathway. Our study is the first to demonstrate that lactucin-induced apoptosis is mediated by ROS production, which in turn activates the caspase-dependent apoptotic pathway by inhibiting BCL-2 and CFLARL expression in Caki-1 cells.

Computational analysis of protein stability and allosteric interaction networks in distinct conformational forms of the SARS-CoV-2 spike D614G mutant: reconciling functional mechanisms through allosteric model of spike regulation

In this study, we used an integrative computational approach to examine molecular mechanisms underlying functional effects of the D614G mutation by exploring atomistic modeling of the SARS-CoV-2 spike proteins as allosteric regulatory machines. We combined coarse-grained simulations, protein stability and dynamic fluctuation communication analysis with network-based community analysis to examine structures of the native and mutant SARS-CoV-2 spike proteins in different functional states. Through distance fluctuations communication analysis, we probed stability and allosteric communication propensities of protein residues in the native and mutant SARS-CoV-2 spike proteins, providing evidence that the D614G mutation can enhance long-range signaling of the allosteric spike engine. By combining functional dynamics analysis and ensemble-based alanine scanning of the SARS-CoV-2 spike proteins we found that the D614G mutation can improve stability of the spike protein in both closed and open forms, but shifting thermodynamic preferences towards the open mutant form. Our results revealed that the D614G mutation can promote the increased number of stable communities and allosteric hub centers in the open form by reorganizing and enhancing the stability of the S1-S2 inter-domain interactions and restricting mobility of the S1 regions. This study provides atomistic-based view of allosteric communications in the SARS-CoV-2 spike proteins, suggesting that the D614G mutation can exert its primary effect through allosterically induced changes on stability and communications in the residue interaction networks. Communicated by Ramaswamy H. Sarma

Dynamic Network Modeling of Allosteric Interactions and Communication Pathways in the SARS-CoV-2 Spike Trimer Mutants: Differential Modulation of Conformational Landscapes and Signal Transmission via Cascades of Regulatory Switches.

The rapidly growing body of structural and biochemical studies of the SARS-CoV-2 spike glycoprotein has revealed a variety of distinct functional states with radically different arrangements of the receptor-binding domain, highlighting a remarkable function-driven conformational plasticity and adaptability of the spike proteins. In this study, we examined molecular mechanisms underlying conformational and dynamic changes in the SARS-CoV-2 spike mutant trimers through the lens of dynamic analysis of allosteric interaction networks and atomistic modeling of signal transmission. Using an integrated approach that combined coarse-grained molecular simulations, protein stability analysis, and perturbation-based modeling of residue interaction networks, we examined how mutations in the regulatory regions of the SARS-CoV-2 spike protein can differentially affect dynamics and allosteric signaling in distinct functional states. The results of this study revealed key functional regions and regulatory centers that govern collective dynamics, allosteric interactions, and control signal transmission in the SARS-CoV-2 spike proteins. We found that the experimentally confirmed regulatory hotspots that dictate dynamic switching between conformational states of the SARS-CoV-2 spike protein correspond to the key hinge sites and global mediating centers of the allosteric interaction networks. The results of this study provide a novel insight into allosteric regulatory mechanisms of SARS-CoV-2 spike proteins showing that mutations at the key regulatory positions can differentially modulate distribution of states and determine topography of signal communication pathways operating through state-specific cascades of control switch points. This analysis provides a plausible strategy for allosteric probing of the conformational equilibrium and therapeutic intervention by targeting specific hotspots of allosteric interactions and communications in the SARS-CoV-2 spike proteins.

Cataract-causing mutations L45P and Y46D impair the thermal stability of γC-crystallin

Crystallin gene mutations are responsible for about half of the congenital cataract caused by genetic disorders. L45P and Y46D mutations of γC-crystallin have been reported in patients with nuclear congenital cataract. In this study, we explored the thermal stability of wild type (WT), L45P, and Y46D mutants of γC-crystallin at low and high concentrations, as well as the effect of αA-crystallin on the thermal stability of mutants. Spectroscopic experiments were used to monitor the structural changes on temperature-gradient and time-course heating process. Intermediate morphologies were determined through cryo-electron microscopy. The thermal stability of WT and mutants at concentrations ranging up to hundreds of milligrams were assessed via the UNcle multifunctional protein stability analysis system. The results showed that L45P and Y46D mutations impaired the thermal stability of γC-crystallin at low (0.2 mg/mL) and high concentrations (up to 200 mg/mL). Notably, with increase in protein concentration, the thermal stability of L45P and Y46D mutants of γC-crystallin simultaneously decreased. Thermal stability of L45P and Y46D mutants could be rescued by αA-crystallin in a concentration-dependent manner. The dramatic decrease in thermal stability of γC-crystallin caused by L45P and Y46D mutations contributed to congenital cataract in the mature human lens.

Molecular basis for arginine C-terminal degron recognition by Cul2FEM1 E3 ligase.

Degrons are elements within protein substrates that mediate the interaction with specific degradation machineries to control proteolysis. Recently, a few classes of C-terminal degrons (C-degrons) that are recognized by dedicated cullin-RING ligases (CRLs) have been identified. Specifically, CRL2 using the related substrate adapters FEM1A/B/C was found to recognize C degrons ending with arginine (Arg/C-degron). Here, we uncover the molecular mechanism of Arg/C-degron recognition by solving a subset of structures of FEM1 proteins in complex with Arg/C-degron-bearing substrates. Our structural research, complemented by binding assays and global protein stability (GPS) analyses, demonstrates that FEM1A/C and FEM1B selectively target distinct classes of Arg/C-degrons. Overall, our study not only sheds light on the molecular mechanism underlying Arg/C-degron recognition for precise control of substrate turnover, but also provides valuable information for development of chemical probes for selectively regulating proteostasis.

Analysis of Protein Stability by Synthesis Shutoff.

In this protocol, we describe the analysis of protein stability over time, using synthesis shutoff. As an example, we express HA-tagged yeast mitofusin Fzo1 in Saccharomyces cerevisiae and inhibit translation via cycloheximide (CHX). Proteasomal inhibition with MG132 is performed, as an optional step, before the addition of CHX. Proteins are extracted via trichloroacetic acid (TCA) precipitation and subsequently separated via SDS-PAGE. Immunoblotting and antibody-decoration are performed to detect Fzo1 using HA-specific antibodies. We have adapted the method of blocking protein translation with cycloheximide to analyze the stability of high molecular weight proteins, including post-translational modifications and their impact on protein turnover.