Vienna-PTM web server: a toolkit for MD simulations of protein post-translational modifications

S-nitrosylation of EZH2 alters PRC2 assembly, methyltransferase activity, and EZH2 stability to maintain endothelial homeostasis

Nitric oxide (NO), a versatile bio-active molecule modulates cellular functions through diverse mechanisms including S-nitrosylation of proteins. Herein, we report S-nitrosylation of selected cysteine residues of EZH2 in endothelial cells, which interplays with its stability and functions. We detect a significant reduction in H3K27me3 upon S-nitrosylation of EZH2 as contributed by the early dissociation of SUZ12 from the PRC2. Moreover, S-nitrosylation of EZH2 causes its cytosolic translocation, ubiquitination, and degradation. Further analysis reveal S-nitrosylation of cysteine 329 induces EZH2 instability, whereas S-nitrosylation of cysteine 700 abrogates its catalytic activity. We further show that S-nitrosylation-dependent regulation of EZH2 maintains endothelial homeostasis in both physiological and pathological settings. Molecular dynamics simulation reveals the inability of SUZ12 to efficiently bind to the SAL domain of EZH2 upon S-nitrosylation. Taken together, our study reports S-nitrosylation-dependent regulation of EZH2 and its associated PRC2 complex, thereby influencing the epigenetics of endothelial homeostasis.

Protein editing using a coordinated transposition reaction

Protein engineering through the ligation of polypeptide fragments has proven enormously powerful for studying biochemical processes. In general, this strategy necessitates a final protein-folding step, constraining the types of systems amenable to the approach. Here, we report a method that allows internal regions of target proteins to be replaced in a single operation. Conceptually, our system is analogous to a DNA transposition reaction but uses orthogonal pairs of engineered split inteins to mediate the editing process. This "protein transposition" reaction is applied to several systems, including folded protein complexes, allowing the efficient introduction of a variety of noncoded elements. By carrying out a molecular "cut and paste" under native protein-folding conditions, our approach substantially expands the scope of protein semisynthesis.

Identification and inhibition of PIN1-NRF2 protein-protein interactions through computational and biophysical approaches

NRF2 is a transcription factor responsible for coordinating the expression of over a thousand cytoprotective genes. Although NRF2 is constitutively expressed, its stability is modulated by the redox-sensitive protein KEAP1 and other conditional binding partner regulators. The new era of NRF2 research has highlighted the cooperation between NRF2 and PIN1 in modifying its cytoprotective effect. Despite numerous studies, the understanding of the PIN1-NRF2 interaction remains limited. Herein, we described the binding interaction of PIN1 and three different 14-mer long phospho-peptides mimicking NRF2 protein using computer-based, biophysical, and biochemical approaches. According to our computational analyses, the residues positioned in the WW domain of PIN1 (Ser16, Arg17, Ser18, Tyr23, Ser32, Gln33, and Trp34) were found to be crucial for PIN1-NRF2 interactions. Biophysical FP assays were used to verify the computational prediction. The data demonstrated that Pintide, a peptide predominantly interacting with the PIN1 WW-domain, led to a significant reduction in the binding affinity of the NRF2 mimicking peptides. Moreover, we evaluated the impact of known PIN1 inhibitors (juglone, KPT-6566, and EGCG) on the PIN1-NRF2 interaction. Among the inhibitors, KPT-6566 showed the most potent inhibitory effect on PIN1-NRF2 interaction within an IC50 range of 0.3-1.4 µM. Furthermore, our mass spectrometry analyses showed that KPT-6566 appeared to covalently modify PIN1 via conjugate addition, rather than disulfide exchange of the sulfonyl-acetate moiety. Altogether, such inhibitors would also be highly valuable molecular probes for further investigation of PIN1 regulation of NRF2 in the cellular context and potentially pave the way for drug molecules that specifically inhibit the cytoprotective effects of NRF2 in cancer.

Structural Analysis of Plasma-Induced Oxidation and Electric Field Effect on the Heat Shock Protein (Hsp60) Structure: A Computational Viewpoint

In recent years, there has been an increase in the study of the mechanisms behind plasma oncology. For this, many wet lab experiments and computational studies were conducted. Computational studies give an advantage in examining protein structures that are costly to extract in enough amounts to analyze the biophysical properties following plasma treatment. Therefore, in this work, we studied the effect of plasma oxidation and electric field on the human mitochondrial heat shock protein (mHsp60). Hsp60, alias chaperonin, is one of the most conserved proteins expressed across all species. Hence, we performed molecular dynamic simulations to calculate the root-mean-square deviation, root-mean-square fluctuation, and solvent-accessible surface area of mHsp60 with and without oxidation. In addition to the oxidation state, we also applied an electric field (0.003 and 2.0 V/nm) to check the changes in the mHsp60 protein. Through simulations, we observed that the electric field strongly affects the structure of mHsp60 protein compared with the oxidation. The combination of oxidation and electric field effect increases the destabilization of the mHsp60 structure compared with their respective control states.

Navigating Uncertainty: Assessing Variants of Uncertain Significance in the CDKL5 Gene for Developmental and Epileptic Encephalopathy Using In Silico Prediction Tools and Computational Analysis

Mutations in the CDKL5 gene are associated with developmental and epileptic encephalopathy (DEE), a severe disorder characterized by developmental delay and epileptic activity. In genetic analyses of DEEs, variants classified as pathogenic confirm the diagnosis of the disease while Variants of Uncertain Significance (VUS) remain in a gray area due to insufficient evidence. This study aimed to optimize the interpretation of VUS in the CDKL5 gene by evaluating the performance of 22 in silico prediction tools using 186 known pathogenic or benign missense variants from the ClinVar database. The best-performing tools were then applied to analyze CDKL5 VUS variants, complemented by the evaluation of evolutionary conservation, structural analyses, and molecular dynamics simulations to assess their impact on protein structure and function. The results identified SNPred as the most reliable tool, achieving 100% accuracy, sensitivity, and specificity. Other high-performing tools, including ESM-1v, AlphaMissense, EVE, and ClinPred, demonstrated over 98% accuracy. Among 44 CDKL5 VUS variants evaluated, 20 were initially classified as pathogenic by these tools. However, further evaluation using stringent criteria-incorporating conservation scores, structural disruptions identified by Missense3D and PyMol, and molecular dynamics simulation results-led to the reclassification of 8 VUS variants as "potentially pathogenic" and the remaining 12 as "variants with conflicting data". This comprehensive approach provides a robust framework for the classification of VUS in the CDKL5 gene, offering critical insights for accurate diagnosis and treatment strategies in DEE. These findings will serve as a valuable resource for clinicians and geneticists in resolving the diagnostic ambiguity associated with VUS.

Effect of plasma-induced oxidation on NK cell immune checkpoint ligands: A computational-experimental approach

Non-thermal plasma (NTP) shows promise as a potent anti-cancer therapy with both cytotoxic and immunomodulatory effects. In this study, we investigate the chemical and biological effects of NTP-induced oxidation on several key, determinant immune checkpoints of natural killer (NK) cell function. We used molecular dynamics (MD) and umbrella sampling simulations to investigate the effect of NTP-induced oxidative changes on the MHC-I complexes HLA-Cw4 and HLA-E. Our simulations indicate that these chemical alterations do not significantly affect the binding affinity of these markers to their corresponding NK cell receptor, which is supported with experimental read-outs of ligand expression on human head and neck squamous cell carcinoma cells after NTP application. Broadening our scope to other key ligands for NK cell reactivity, we demonstrate rapid reduction in CD155 and CD112, target ligands of the inhibitory TIGIT axis, and in immune checkpoint CD73 immediately after treatment. Besides these transient chemical alterations, the reactive species in NTP cause a cascade of downstream cellular reactions. This is underlined by the upregulation of the stress proteins MICA/B, potent ligands for NK cell activation, 24 h post treatment. Taken together, this work corroborates the immunomodulatory potential of NTP, and sheds light on the interaction mechanisms between NTP and cancer cells.

Molecular Dynamics Reveals Altered Interactions between Belzutifan and HIF-2 with Natural Variant G323E or Proximal Phosphorylation at T324

In patients with von-Hippel Lindau (VHL) disease, hypoxia-independent accumulation of HIF-2α leads to increased transcriptional activity of HIF-2α:ARNT that drives cancers such as renal cell carcinoma. Belzutifan, a recently FDA-approved drug, is designed to prevent the transcriptional activity of HIF-2α:ARNT, thereby overcoming the consequences of its unnatural accumulation in VHL-dependent cancers. Emerging evidence suggests that the naturally occurring variant G323E located in the HIF-2α drug binding pocket prevents inhibitory activity of belzutifan analogs, though the mechanism of inhibition remains unclear. Interestingly, proximal phosphorylation at neighboring T324, previously shown to regulate HIF-2 protein interactions, has also been proposed to affect HIF-2 drug binding. Here, we used molecular dynamics (MD) simulations to understand and compare the molecular-level effects of G323E and phospho-T324 (pT324) on the belzutifan bound-HIF-2α:ARNT complex. We find that both G323E and pT324 increase structural flexibility within the drug binding site and reduce the apparent binding affinity for belzutifan. Whereas the effects of G323E are concentrated in the binding pocket Fα helix within the HIF-2α PAS-B domain, pT324 decreased the belzutifan binding affinity and stabilized the HIF-2 heterodimer through an alternate mechanism involving polar interactions between the HIF-2α PAS-B and PAS-A domains. Further analysis via ensemble machine learning uncovered important and distinct interchain residue interactions modified by G323E and pT324. These findings reveal a molecular mechanism of G323E-induced drug resistance and suggest that pT324 may also affect the efficacy of HIF-2 drug binding interactions via allosteric effects.

An Integrated Analytical Approach for Screening Functional Post-Translational Modification Sites in Metabolic Enzymes

Post-translational modifications (PTMs) are pivotal in the orchestration of diverse physiological and pathological processes. Despite this, the identification of functional PTM sites within the vast amount of data remains challenging. Conventionally, those PTM sites are discerned through labor-intensive and time-consuming experiments. Here, we developed an integrated analytical approach for the identification of functional PTM sites on metabolic enzymes via a screening process. Through gene ontology (GO) analysis, we identified 269 enzymes with lysine 2-hydroxyisobutyrylation (Khib) from our proteomics data set of Escherichia coli. The first round of screening was performed based on the enzyme structures/predicted structures using the TM-score engineer, a tool designed to evaluate the impact of PTM on the protein structure. Subsequently, we examined the influence of Khib on the enzyme-substrate interactions through both static and dynamic analyses, molecular docking, and molecular dynamics simulation. Ultimately, we identified NfsB K181hib and ThiF K83hib as potential functional sites. This work has established a novel analytical approach for the identification of functional protein PTM sites, thereby contributing to the understanding of Khib functions.

Effect of α-tubulin acetylation on the doublet microtubule structure

Acetylation of α-tubulin at the lysine 40 residue (αK40) by αTAT1/MEC-17 acetyltransferase modulates microtubule properties and occurs in most eukaryotic cells. Previous literatures suggest that acetylated microtubules are more stable and damage resistant. αK40 acetylation is the only known microtubule luminal post-translational modification site. The luminal location suggests that the modification tunes the lateral interaction of protofilaments inside the microtubule. In this study, we examined the effect of tubulin acetylation on the doublet microtubule (DMT) in the cilia of Tetrahymena thermophila using a combination of cryo-electron microscopy, molecular dynamics, and mass spectrometry. We found that αK40 acetylation exerts a small-scale effect on the DMT structure and stability by influencing the lateral rotational angle. In addition, comparative mass spectrometry revealed a link between αK40 acetylation and phosphorylation in cilia.

Trioxidized cysteine in the aging proteome mimics the structural dynamics and interactome of phosphorylated serine

Aging is the primary risk factor for the development of numerous human chronic diseases. On a molecular level, it significantly impacts the regulation of protein modifications, leading to the accumulation of degenerative protein modifications (DPMs) such as aberrant serine phosphorylation (p-Ser) and trioxidized cysteine (t-Cys) within the proteome. The altered p-Ser is linked to abnormal cell signaling, while the accumulation of t-Cys is associated with chronic diseases induced by oxidative stress. Despite this, the potential cross-effects and functional interplay between these two critical molecular factors of aging remain undisclosed. This study analyzes the aging proteome of wild-type C57BL/6NTac mice over 2 years using advanced proteomics and bioinformatics. Our objective is to provide a comprehensive analysis of how t-Cys affects cell signaling and protein structure in the aging process. The results obtained indicate that t-Cys residues accumulate in the aging proteome, interact with p-Ser interacting enzymes, as validated in vitro, and alter their structures similarly to p-Ser. These findings have significant implications for understanding the interplay of oxidative stress and phosphorylation in the aging process. Additionally, they open new venues for further research on the role(s) of these protein modifications in various human chronic diseases and aging, wherein exacerbated oxidation and aberrant phosphorylation are implicated.

Tyrosine phosphorylation of recombinant hirudin increases affinity to thrombin and antithrombotic activity

Thrombin is one of the key enzymes of the blood coagulation system and a promising target for the development of anticoagulants. One of the most specific natural thrombin inhibitors is hirudin, contained in the salivary glands of medicinal leeches. The medicinal use of recombinant hirudin is limited because of the lack of sulfation on Tyr63, resulting in a 10-fold decrease in activity compared to native (sulfated) hirudin. In the present work, a set of hirudin derivatives was tested for affinity to thrombin: phospho-Tyr63, Tyr63(carboxymethyl)Phe, and Tyr63Glu mutants, which mimic Tyr63 sulfation and Gln65Glu mutant and lysine-succinylated hirudin, which enhance the overall negative charge of hirudin, as well as sulfo-hirudin and desulfo-hirudin as references. Using steered molecular dynamics simulations with subsequent umbrella sampling, phospho-hirudin was shown to exhibit the highest affinity to thrombin among all hirudin analogs, including native sulfo-hirudin; succinylated hirudin was also prospective. Phospho-hirudin exhibited the highest antithrombotic activity in in vitro assay in human plasma. Taking into account the modern methods for obtaining phospho-hirudin and succinylated hirudin, they are prospective as anticoagulants in clinical practice.

Prediction and biological analysis of yeast VDAC1 phosphorylation

The mitochondrial outer membrane protein porin 1 (Por1), the yeast orthologue of mammalian voltage-dependent anion channel (VDAC), is the major permeability pathway for the flux of metabolites and ions between cytosol and mitochondria. In yeast, several Por1 phosphorylation sites have been identified. Protein phosphorylation is a major modification regulating a variety of biological activities, but the potential biological roles of Por1 phosphorylation remains unaddressed. In this work, we analysed 10 experimentally observed phosphorylation sites in yeast Por1 using bioinformatics tools. Two of the residues, T100 and S133, predicted to reduce and increase pore permeability, respectively, were validated using biological assays. In accordance, Por1T100D reduced mitochondrial respiration, while Por1S133E phosphomimetic mutant increased it. Por1T100A expression also improved respiratory growth, while Por1S133A caused defects in all growth conditions tested, notably in fermenting media. In conclusion, we found phosphorylation has the potential to modulate Por1, causing a marked effect on mitochondrial function. It can also impact on cell morphology and growth both in respiratory and, unpredictably, also in fermenting conditions, expanding our knowledge on the role of Por1 in cell physiology.

A Model of the Full-Length Cytokinin Receptor: New Insights and Prospects

Cytokinins (CK) are one of the most important classes of phytohormones that regulate a wide range of processes in plants. A CK receptor, a sensor hybrid histidine kinase, was discovered more than 20 years ago, but the structural basis for its signaling is still a challenge for plant biologists. To date, only two fragments of the CK receptor structure, the sensory module and the receiver domain, were experimentally resolved. Some other regions were built up by molecular modeling based on structures of proteins homologous to CK receptors. However, in the long term, these data have proven insufficient for solving the structure of the full-sized CK receptor. The functional unit of CK receptor is the receptor dimer. In this article, a molecular structure of the dimeric form of the full-length CK receptor based on AlphaFold Multimer and ColabFold modeling is presented for the first time. Structural changes of the receptor upon interacting with phosphotransfer protein are visualized. According to mathematical simulation and available data, both types of dimeric receptor complexes with hormones, either half- or fully liganded, appear to be active in triggering signals. In addition, the prospects of using this and similar models to address remaining fundamental problems of CK signaling were outlined.

In silico designed microtubule-stabilizer drugs against tauopathy in Alzheimer's disease

Microtubules are the main building blocks of the cytoskeleton that maintain the shape of the cell. Microtubule-associated proteins, such as Tau protein, facilitate their plasticity in cells. Highly phosphorylated Tau has weak affinity to microtubule and, hence, high probability of aggregation into neurofibrillary tangles (tauopathy). Alzheimer's disease evolves when Tau proteins are abnormally phosphorylated. To prevent tauopathy in Alzheimer's disease, we designed drugs de novo targeting them in silico to the phosphorylated Tau-microtubule complexes. Our molecular docking (AutoDock, MOE, GOLD) and molecular dynamics (GROMACS, 2019.6) simulation results revealed compound 23 (C12H28N4O5) as a potential drug candidate, since it can bind (-11.1 kcal/mol by AutoDock) and fix not only phosphorylated Tau on the surface of microtubules, but also prevent their aggregation into bundles. In addition, compound 23 has shown its ability to de-bundle already grouped phosphorylated peptides into single pieces.Communicated by Ramaswamy H. Sarma.

Tuning transcription factor DegU for developing extracellular protease overproducer in Bacillus pumilus

BACKGROUND

Global transcription machinery engineering (gTME) is an effective approach employed in strain engineering to rewire gene expression and reshape cellular metabolic fluxes at the transcriptional level.

RESULTS

In this study, we utilized gTME to engineer the positive transcription factor, DegU, in the regulation network of major alkaline protease, AprE, in Bacillus pumilus. To validate its functionality when incorporated into the chromosome, we performed several experiments. First, three negative transcription factors, SinR, Hpr, and AbrB, were deleted to promote AprE synthesis. Second, several hyper-active DegU mutants, designated as DegU(hy), were selected using the fluorescence colorimetric method with the host of the Bacillus subtilis ΔdegSU mutant. Third, we integrated a screened degU(L113F) sequence into the chromosome of the Δhpr mutant of B. pumilus SCU11 to replace the original degU gene using a CRISPR/Cas9 system. Finally, based on transcriptomic and molecular dynamic analysis, we interpreted the possible mechanism of high-yielding and found that the strain produced alkaline proteases 2.7 times higher than that of the control strain (B. pumilus SCU11) in LB medium.

CONCLUSION

Our findings serve as a proof-of-concept that tuning the global regulator is feasible and crucial for improving the production performance of B. pumilus. Additionally, our study established a paradigm for gene function research in strains that are difficult to handle.

Adnp-mutant mice with cognitive inflexibility, CaMKIIα hyperactivity, and synaptic plasticity deficits

ADNP syndrome, involving the ADNP transcription factor of the SWI/SNF chromatin-remodeling complex, is characterized by developmental delay, intellectual disability, and autism spectrum disorders (ASD). Although Adnp-haploinsufficient (Adnp-HT) mice display various phenotypic deficits, whether these mice display abnormal synaptic functions remain poorly understood. Here, we report synaptic plasticity deficits associated with cognitive inflexibility and CaMKIIα hyperactivity in Adnp-HT mice. These mice show impaired and inflexible contextual learning and memory, additional to social deficits, long after the juvenile-stage decrease of ADNP protein levels to ~10% of the newborn level. The adult Adnp-HT hippocampus shows hyperphosphorylated CaMKIIα and its substrates, including SynGAP1, and excessive long-term potentiation that is normalized by CaMKIIα inhibition. Therefore, Adnp haploinsufficiency in mice leads to cognitive inflexibility involving CaMKIIα hyperphosphorylation and excessive LTP in adults long after its marked expressional decrease in juveniles.

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.

Phosphorylation Modification Force Field FB18CMAP Improving Conformation Sampling of Phosphoproteins

Phosphorylation of proteins plays an important regulatory role at almost all levels of cellular organization. Molecular dynamics (MD) simulation is a promising tool to reveal the mechanism of how phosphorylation regulates many key biological processes at the atomistic level. MD simulation accuracy depends on force field precision, while the current force fields for phospho-amino acids have resulted in notable inconsistency with experimental data. Here, a new force field parameter (named FB18CMAP) is generated by fitting against quantum mechanics (QM) energy in aqueous solution with φ/ψ dihedral potential-energy surfaces optimized using CMAP parameters. MD simulations of phosphorylated dipeptides, intrinsically disordered proteins (IDPs), and ordered (folded) proteins show that FB18CMAP can mimic NMR observables and structural characteristics of phosphorylated dipeptides and proteins more accurately than the FB18 force field. These findings suggest that FB18CMAP performs well in both the simulation of ordered and disordered states of phosphorylated proteins.

Structural basis of ZAP-70 activation upon phosphorylation of tyrosines 315, 319 and 493

ZAP-70 (Zeta-chain-Associated Protein kinase 70) is a key kinase in the regulation of the adaptive immune response. Zap-70 acts by binding its SH2-domains to the T-cell-associated CD3ζ protein, thus transmitting a T-cell activation signal induced by the interaction of Major Histocompatibility Complex with T-cell Receptor. It has been established that for ZAP-70 kinase activation, the phosphorylation of Tyr315, Tyr319, and Tyr493 is required, however the mechanisms are unclear. In the present study, we use the tools of structural modeling to elucidate the ZAP-70 activation mechanisms.

Tyrosine phosphorylation controlled poly(A) polymerase I activity regulates general stress response in bacteria

RNA 3'-end polyadenylation that marks transcripts for degradation is implicated in general stress response in Escherichia coli Yet, the mechanism and regulation of poly(A) polymerase I (PAPI) in stress response are obscure. We show that pcnB (that encodes PAPI)-null mutation widely stabilises stress response mRNAs and imparts cellular tolerance to multiple stresses, whereas PAPI ectopic expression renders cells stress-sensitive. We demonstrate that there is a substantial loss of PAPI activity on stress exposure that functionally phenocopies pcnB-null mutation stabilising target mRNAs. We identify PAPI tyrosine phosphorylation at the 202 residue (Y202) that is enormously enhanced on stress exposure. This phosphorylation inhibits PAPI polyadenylation activity under stress. Consequentially, PAPI phosphodeficient mutation (tyrosine 202 to phenylalanine, Y202F) fails to stimulate mRNA expression rendering cells stress-sensitive. Bacterial tyrosine kinase Wzc phosphorylates PAPI-Y202 residue, and that wzc-null mutation renders cells stress-sensitive. Accordingly, wzc-null mutation has no effect on stress sensitivity in the presence of pcnB-null or pcnB-Y202F mutation. We also establish that PAPI phosphorylation-dependent stress tolerance mechanism is distinct and operates downstream of the primary stress regulator RpoS.

Facilitating CG Simulations with MAD: The MArtini Database Server

The MArtini Database (MAD - https://mad.ibcp.fr) is a web server designed for the sharing of structures and topologies of molecules parametrized with the Martini coarse-grained (CG) force field. MAD can also convert atomistic structures into CG structures and prepare complex systems (including proteins, lipids, etc.) for molecular dynamics (MD) simulations at the CG level. It is dedicated to the generation of input files for Martini 3, the most recent version of this popular CG force field. Specifically, the MAD server currently includes tools to submit or retrieve CG models of a wide range of molecules (lipids, carbohydrates, nanoparticles, etc.), transform atomistic protein structures into CG structures and topologies, with fine control on the process and assemble biomolecules into large systems, and deliver all files necessary to start simulations in the GROMACS MD engine.

The origin of pigment-binding differences in CP29 and LHCII: the role of protein structure and dynamics

The first step of photosynthesis in plants is performed by the light-harvesting complexes (LHC), a large family of pigment-binding proteins embedded in the photosynthetic membranes. These complexes are conserved across species, suggesting that each has a distinct role. However, they display a high degree of sequence homology and their static structures are almost identical. What are then the structural features that determine their different properties? In this work, we compared the two best-characterized LHCs of plants: LHCII and CP29. Using molecular dynamics simulations, we could rationalize the difference between them in terms of pigment-binding properties. The data also show that while the loops between the helices are very flexible, the structure of the transmembrane regions remains very similar in the crystal and the membranes. However, the small structural differences significantly affect the excitonic coupling between some pigment pairs. Finally, we analyzed in detail the structure of the long N-terminus of CP29, showing that it is structurally stable and it remains on top of the membrane even in the absence of other proteins. Although the structural changes upon phosphorylation are minor, they can explain the differences in the absorption properties of the pigments observed experimentally.

Acetylation and phosphorylation processes modulate Tau's binding to microtubules: A molecular dynamics study

The microtubule-associated protein Tau has its normal function impaired when undergoing post-translational modifications. In this work, molecular modelling techniques were used to infer the effects of acetylation and phosphorylation in Tau's overall conformation, electrostatics, and interactions, but mostly in Tau's ability to bind microtubules. Reported harmful Lys sites were mutated by its acetylated form, generating eight different acetylated Tau (aTau) analogues. Similarly, phosphorylation sites found in normal brains and in Alzheimer's lesioned brains were considered to design phosphorylated Tau (pTau) analogues. All these designed variants were evaluated in intracellular fluid and near a microtubule (MT) model. Our in silico findings demonstrated that the electrostatic changes, due to the absence of positive Lys' charges in acetylation cases, or the increasingly negative charge in the phosphorylated forms, hamper the association to the MT tubulins in most cases. Post-translational modifications also pose very distinct conformations to the ones described for native Tau, which hinders the microtubule-binding region (MTBR) and turns difficult the expected binding. Our study elucidates important molecular processes behind Tau abnormal function which can inspire novel therapeutics to address Alzheimer's disease.

Identification of peptide binding sequence of TRIM25 on 14-3-3σ by bioinformatics and biophysical techniques

14-3-3σ protein is one of the seven isoforms from the highly conserved eukaryotic 14-3-3 protein family. Downregulation of 14-3-3σ expression has been observed in various tumors. TRIM25 is responsible for the proteolytic degradation of 14-3-3σ, in which abrogation of TRIM25 suppressed tumor growth through 14-3-3σ upregulation. However, to date, the exact 14-3-3σ interacting residues of TRIM25 have yet to be resolved. Thus, this study attempts to identify the peptide binding sequence of TRIM25 on 14-3-3σ via both bioinformatics and biophysical techniques. Multiple sequence alignment of the CC domain of TRIM25 revealed five potential peptide binding sequences (Peptide 1-5). Nuclear magnetic resonance (NMR) assay (1H CPMG) identified Peptide 1 as an important sequence for binding to 14-3-3σ. Competition NMR assay suggested that Peptide 1 binds to the amphipathic pocket of 14-3-3σ with an estimated KD of 116.4 µM by isothermal titration calorimetry. Further in silico docking and molecular dynamics simulations studies proposed that Peptide 1 is likely to interact with Lys49, Arg56, Arg129, and Tyr130 residues at the amphipathic pocket of 14-3-3σ. These results suggest that Peptide 1 may serve as a biological probe or a template to design inhibitors of TRIM25-14-3-3σ interaction as a potentially novel class of anticancer agents.Communicated by Ramaswamy H. Sarma.

Molecular cloning, expression, mRNA secondary structure and immunological characterization of mussel foot proteins (Mfps) (Mollusca: Bivalvia)

The macroscale production of mussel foot proteins (Mfps) in the expression system has not succeeded to date. The principal reasons for this are low levels of expression and yield of Mfps, lack of post-translational modifications (PTMs), and immunological toxic effects on the host system. Identification of post-translational modification sites, suitable expression hosts, and immunological responses through an experimental approach is very costly and time-consuming. However, in the present study, in silico post-translation modification, antigenicity, allergenicity, and the immunological reaction of all available Mfps were characterized. Furthermore, all Mfps were codon optimized in three different expression systems to determine the best expression host. Finally, we performed the in-silico cloning of all codon-optimized Mfps in a suitable host (E. coli K12, pET28a(+) vector) and analyzed the secondary structure of mRNA and its structural stability. Among the 78 Mfps, six fps are considered potential allergenic proteins, six fps are considered non-allergenic proteins, and all other fps are probably allergenic. High antigenicity was observed in bacterial cells as compared to yeast and tumor cells. Nevertheless, the predicted expression of Mfps in a bacterial host is higher than in other expression hosts. Important to note that all Mfps showed significant immunological activity in the human system, and we concluded that these antigenic, allergenic, and immunological properties are directly correlated with their amino acid composition. The study's major goal is to provide a comprehensive understanding of Mfps and aid in the future genetic engineering and expression of Mfps and its diverse applications in different fields.Communicated by Ramaswamy H. Sarma.

New simulation insights on the structural transition mechanism of bovine rhodopsin activation

Inactive rhodopsin can absorb photons, which induces different structural transitions that finally activate rhodopsin. We have examined the change in spatial configurations and physicochemical factors that result during the transition mechanism from the inactive to the active rhodopsin state via intermediates. During the activation process, many existing atomic contacts are disrupted, and new ones are formed. This is related to the movement of Helix 5, which tilts away from Helix 3 in the intermediate state in lumirhodopsin and moves closer to Helix 3 again in the active state. Similar patterns of changing atomic contacts are observed between Helices 3 and 5 of the adenosine and neurotensin receptors. In addition, residues 220-238 of rhodopsin, which are disordered in the inactive state, fold in the active state before binding to the Gα, where it catalyzes GDP/GTP exchange on the Gα subunit. Finally, molecular dynamics simulations in the membrane environment revealed that the arrestin binding region adopts a more flexible extended conformation upon phosphorylation, likely promoting arrestin binding and inactivation. In summary, our results provide additional structural understanding of specific rhodopsin activation which might be relevant to other Class A G protein-coupled receptor proteins.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Balanced Force Field ff03CMAP Improving the Dynamics Conformation Sampling of Phosphorylation Site

Phosphorylation plays a key role in plant biology, such as the accumulation of plant cells to form the observed proteome. Statistical analysis found that many phosphorylation sites are located in disordered regions. However, current force fields are mainly trained for structural proteins, which might not have the capacity to perfectly capture the dynamic conformation of the phosphorylated proteins. Therefore, we evaluated the performance of ff03CMAP, a balanced force field between structural and disordered proteins, for the sampling of the phosphorylated proteins. The test results of 11 different phosphorylated systems, including dipeptides, disordered proteins, folded proteins, and their complex, indicate that the ff03CMAP force field can better sample the conformations of phosphorylation sites for disordered proteins and disordered regions than ff03. For the solvent model, the results strongly suggest that the ff03CMAP force field with the TIP4PD water model is the best combination for the conformer sampling. Additional tests of CHARMM36m and FB18 force fields on two phosphorylated systems suggest that the overall performance of ff03CMAP is similar to that of FB18 and better than that of CHARMM36m. These results can help other researchers to choose suitable force field and solvent models to investigate the dynamic properties of phosphorylation proteins.

Electrostatic and steric effects underlie acetylation-induced changes in ubiquitin structure and function

Covalent attachment of ubiquitin (Ub) to proteins is a highly versatile posttranslational modification. Moreover, Ub is not only a modifier but itself is modified by phosphorylation and lysine acetylation. However, the functional consequences of Ub acetylation are poorly understood. By generation and comprehensive characterization of all seven possible mono-acetylated Ub variants, we show that each acetylation site has a particular impact on Ub structure. This is reflected in selective usage of the acetylated variants by different E3 ligases and overlapping but distinct interactomes, linking different acetylated variants to different cellular pathways. Notably, not only electrostatic but also steric effects contribute to acetylation-induced changes in Ub structure and, thus, function. Finally, we provide evidence that p300 acts as a position-specific Ub acetyltransferase and HDAC6 as a general Ub deacetylase. Our findings provide intimate insights into the structural and functional consequences of Ub acetylation and highlight the general importance of Ub acetylation.

Ubiquitin is not only a posttranslational modifier but itself is subject to modifications, such as acetylation. Characterization of distinct acetylated ubiquitin variants reveals that each acetylation site has a particular impact on ubiquitin structure and its protein-protein interaction properties.

Structural insights into the pSer/pThr dependent regulation of the SHP2 tyrosine phosphatase in insulin and CD28 signaling

Serine/threonine phosphorylation of insulin receptor substrate (IRS) proteins is well known to modulate insulin signaling. However, the molecular details of this process have mostly been elusive. While exploring the role of phosphoserines, we have detected a direct link between Tyr-flanking Ser/Thr phosphorylation sites and regulation of specific phosphotyrosine phosphatases. Here we present a concise structural study on how the activity of SHP2 phosphatase is controlled by an asymmetric, dual phosphorylation of its substrates. The structure of SHP2 has been determined with three different substrate peptides, unveiling the versatile and highly dynamic nature of substrate recruitment. What is more, the relatively stable pre-catalytic state of SHP2 could potentially be useful for inhibitor design. Our findings not only show an unusual dependence of SHP2 catalytic activity on Ser/Thr phosphorylation sites in IRS1 and CD28, but also suggest a negative regulatory mechanism that may also apply to other tyrosine kinase pathways as well.

SHP2 is an important human tyrosine phosphatase with key roles in cancer, immune responses and insulin signaling. Here, the authors explore its substrate recognition mechanism in molecular detail and uncover a complex regulatory mechanism for this enzyme that marks specific target sites for dephosphorylation.

Using Atomistic Simulations to Explore the Role of Methylation and ATP in Chemotaxis Signal Transduction

A bacterial chemotaxis mechanism is activated when nutrients bind to surface receptors. The sequence of intra- and interprotein events in this signal cascade from the receptors to the eventual molecular motors has been clearly identified. However, the atomistic details remain elusive, as in general may be expected of intraprotein signal transduction pathways, especially when fibrillar proteins are involved. We performed atomistic calculations of the methyl accepting chemoprotein (MCP)-CheA-CheW multidomain complex from Escherichia coli, simulating the methylated and unmethylated conditions in the chemoreceptors and the ATP-bound and apo conditions of the CheA. Our results indicate that these atomistic simulations, especially with one of the two force fields we tried, capture several relevant features of the downstream effects, such as the methylation favoring an intermediate structure that is more toward a dipped state and increases the chance of ATP hydrolysis. The results thus suggest the sensitivity of the model to reflect the nutrient signal response, a nontrivial validation considering the complexity of the system, encouraging even more detailed studies on the thermodynamic quantification of the effects and the identification of the signaling networks.

O-GlcNAc modification of leucyl-tRNA synthetase 1 integrates leucine and glucose availability to regulate mTORC1 and the metabolic fate of leucine

All living organisms have the ability to sense nutrient levels to coordinate cellular metabolism. Despite the importance of nutrient-sensing pathways that detect the levels of amino acids and glucose, how the availability of these two types of nutrients is integrated is unclear. Here, we show that glucose availability regulates the central nutrient effector mTORC1 through intracellular leucine sensor leucyl-tRNA synthetase 1 (LARS1). Glucose starvation results in O-GlcNAcylation of LARS1 on residue S1042. This modification inhibits the interaction of LARS1 with RagD GTPase and reduces the affinity of LARS1 for leucine by promoting phosphorylation of its leucine-binding site by the autophagy-activating kinase ULK1, decreasing mTORC1 activity. The lack of LARS1 O-GlcNAcylation constitutively activates mTORC1, supporting its ability to sense leucine, and deregulates protein synthesis and leucine catabolism under glucose starvation. This work demonstrates that LARS1 integrates leucine and glucose availability to regulate mTORC1 and the metabolic fate of leucine.

Leucyl-tRNA synthetase 1 (LARS1) is a leucine sensor for mTORC1 signaling and regulates leucine utilization depending on glucose availability. Here, the author show that O-GlcNAcylation of LARS1 is crucial for its ability to regulate mTORC1 activity and leucine metabolism upon glucose starvation.

Ferrochelatase: Mapping the Intersection of Iron and Porphyrin Metabolism in the Mitochondria

Porphyrin and iron are ubiquitous and essential for sustaining life in virtually all living organisms. Unlike iron, which exists in many forms, porphyrin macrocycles are mostly functional as metal complexes. The iron-containing porphyrin, heme, serves as a prosthetic group in a wide array of metabolic pathways; including respiratory cytochromes, hemoglobin, cytochrome P450s, catalases, and other hemoproteins. Despite playing crucial roles in many biological processes, heme, iron, and porphyrin intermediates are potentially cytotoxic. Thus, the intersection of porphyrin and iron metabolism at heme synthesis, and intracellular trafficking of heme and its porphyrin precursors are tightly regulated processes. In this review, we discuss recent advances in understanding the physiological dynamics of eukaryotic ferrochelatase, a mitochondrially localized metalloenzyme. Ferrochelatase catalyzes the terminal step of heme biosynthesis, the insertion of ferrous iron into protoporphyrin IX to produce heme. In most eukaryotes, except plants, ferrochelatase is localized to the mitochondrial matrix, where substrates are delivered and heme is synthesized for trafficking to multiple cellular locales. Herein, we delve into the structural and functional features of ferrochelatase, as well as its metabolic regulation in the mitochondria. We discuss the regulation of ferrochelatase via post-translational modifications, transportation of substrates and product across the mitochondrial membrane, protein-protein interactions, inhibition by small-molecule inhibitors, and ferrochelatase in protozoal parasites. Overall, this review presents insight on mitochondrial heme homeostasis from the perspective of ferrochelatase.

Structures and interactions of insulin-like peptides from cone snail venom

The venomous insulin-like peptides released by certain cone snails stimulate hypoglycemic shock to immobilize fish and catch the prey. Compared to human insulin (hIns), the cone snail insulins (Con-Ins) are typically monomeric and shorter in sequence, yet they exhibit moderate hIns-like biological activity. We have modeled six variants of Con-Ins (G3, K1, K2, T1A, T1B, and T2) and carried out explicit-solvent molecular dynamics (MD) simulations of eight types of insulins, two with known structures (hIns and Con-Ins-G1) and six Con-Ins with modeled structures, to characterize key residues of each insulin that interact with the truncated human insulin receptor (μIR). We show that each insulin/μIR complex is stable during explicit-solvent MD simulations and hIns interactions indicate the highest affinity for the "site 1" of IR. The residue contact maps reveal that each insulin preferably interacts with the αCT peptide than the L1 domain of IR. Through analysis of the average nonbonded interaction energy contribution of every residue of each insulin for the μIR, we probe the residues establishing favorable interactions with the receptor. We compared the interaction energy of each residue of every Con-Ins to the μIR and observed that γ-carboxylated glutamate (Gla), His, Thr, Tyr, Tyr/His, and Asn in Con-Ins are favorable substitutions for GluA4, AsnA21, ValB12, LeuB15, GlyB20, and ArgB22 in hIns, respectively. The identified insulin analogs, although lacking the last eight residues of the B-chain of hIns, bind strongly to μIR. Our findings are potentially useful in designing potent fast-acting therapeutic insulin.

S-Nitrosation of E3 Ubiquitin Ligase Complex Components Regulates Hormonal Signalings in Arabidopsis

E3 ubiquitin ligases mediate the last step of the ubiquitination pathway in the ubiquitin-proteasome system (UPS). By targeting transcriptional regulators for their turnover, E3s play a crucial role in every aspect of plant biology. In plants, SKP1/CULLIN1/F-BOX PROTEIN (SCF)-type E3 ubiquitin ligases are essential for the perception and signaling of several key hormones including auxins and jasmonates (JAs). F-box proteins, TRANSPORT INHIBITOR RESPONSE 1 (TIR1) and CORONATINE INSENSITIVE 1 (COI1), bind directly transcriptional repressors AUXIN/INDOLE-3-ACETIC ACID (AUX/IAA) and JASMONATE ZIM-DOMAIN (JAZ) in auxin- and JAs-depending manner, respectively, which permits the perception of the hormones and transcriptional activation of signaling pathways. Redox modification of proteins mainly by S-nitrosation of cysteines (Cys) residues via nitric oxide (NO) has emerged as a valued regulatory mechanism in physiological processes requiring its rapid and versatile integration. Previously, we demonstrated that TIR1 and Arabidopsis thaliana SKP1 (ASK1) are targets of S-nitrosation, and these NO-dependent posttranslational modifications enhance protein-protein interactions and positively regulate SCFTIR1 complex assembly and expression of auxin response genes. In this work, we confirmed S-nitrosation of Cys140 in TIR1, which was associated in planta to auxin-dependent developmental and stress-associated responses. In addition, we provide evidence on the modulation of the SCFCOI1 complex by different S-nitrosation events. We demonstrated that S-nitrosation of ASK1 Cys118 enhanced ASK1-COI1 protein-protein interaction. Overexpression of non-nitrosable ask1 mutant protein impaired the activation of JA-responsive genes mediated by SCFCOI1 illustrating the functional relevance of this redox-mediated regulation in planta. In silico analysis positions COI1 as a promising S-nitrosation target, and demonstrated that plants treated with methyl JA (MeJA) or S-nitrosocysteine (NO-Cys, S-nitrosation agent) develop shared responses at a genome-wide level. The regulation of SCF components involved in hormonal perception by S-nitrosation may represent a key strategy to determine the precise time and site-dependent activation of each hormonal signaling pathway and highlights NO as a pivotal molecular player in these scenarios.

Current status of PTMs structural databases: applications, limitations and prospects

Protein 3D structures, determined by their amino acid sequences, are the support of major crucial biological functions. Post-translational modifications (PTMs) play an essential role in regulating these functions by altering the physicochemical properties of proteins. By virtue of their importance, several PTM databases have been developed and released in decades, but very few of these databases incorporate real 3D structural data. Since PTMs influence the function of the protein and their aberrant states are frequently implicated in human diseases, providing structural insights to understand the influence and dynamics of PTMs is crucial for unraveling the underlying processes. This review is dedicated to the current status of databases providing 3D structural data on PTM sites in proteins. Some of these databases are general, covering multiple types of PTMs in different organisms, while others are specific to one particular type of PTM, class of proteins or organism. The importance of these databases is illustrated with two major types of in silico applications: predicting PTM sites in proteins using machine learning approaches and investigating protein structure-function relationships involving PTMs. Finally, these databases suffer from multiple problems and care must be taken when analyzing the PTMs data.

Charge neutralization of lysine via carbamylation reveals hidden aggregation hot-spots in tau protein flanking regions

Tau protein is found abundantly in neurofibrillary tangles in Alzheimer's disease (AD). The longest human tau isoform (2N4R) has 44 lysine residues. Several lysine-based post-translational modifications (PTMs) such as glycation, acetylation, ubiquitination, and sumoylation have been implicated not only in AD, but also in other tauopathies. Carbamylation is one such lysine neutralizing age-related nonenzymatic PTM which can modulate the aggregation propensity of tau. In this work, we have studied the aggregation potential of lysine-rich regions of tau upon carbamylation which do not aggregate in their native form. Using an array of biophysical and microscopic analyses, such as ThT kinetic assay, fluorescence microscopy, Congo red staining, and scanning electron microscopy, we demonstrate that peptides derived from four of five such regions exhibit robust fibrillar amyloid formation. These regions are found in the N-terminal projection domain that encompasses proline-rich domain (148-153 and 223-230), repeat domain R1 (253-260), as well as fibrillary core region (368-378), and can be described as hidden aggregation hot-spots which become activated upon carbamylation. We have further compared the impact of carbamylation with acetylation on the aggregation propensity of lysine-rich peptide (²⁵⁴ KKVAVV²⁵⁹ ) using biophysical experiments and molecular dynamics simulations and deduced that carbamylation is a much stronger driver of aggregation than acetylation. Our findings may offer more insight into amyloid fibrils' interaction with hidden aggregation-prone nucleating sequences that act as hot-spots for inducing tau fibrillation.

Phosphorylation of eukaryotic initiation factor eIFiso4E enhances the binding rates to VPg of turnip mosaic virus

Binding of phosphorylated eIFiso4E with viral genome-linked protein (VPg) of turnip mosaic virus was examined by stopped-flow, fluorescence, circular dichroism (CD) spectroscopy, and molecular docking analysis. Phosphorylation of eIFiso4E increased (4-fold) the binding rates as compared to unphosphorylated eIFiso4E with VPg. Stopped-flow kinetic studies of phosphorylated eIFiso4E with VPg showed a concentration-independent conformational change. The dissociation rate was about 3-fold slower for eIFiso4E∙VPg complex upon phosphorylation. Phosphorylation enhanced the association rates and lowered the dissociation rates for the eIFiso4E∙VPg binding, with having higher preferential binding to eIFiso4Ep. Binding rates for the interaction of eIFiso4Ep with VPg increased (6-fold) with an increase in temperature, 278 K to 298 K. The activation energies for binding of eIFiso4Ep and eIFiso4E with VPg were 37.2 ± 2.8 and 52.6 ± 3.6 kJ/mol, respectively. Phosphorylation decreased the activation energy for the binding of eIFiso4E to VPg. The reduced energy barrier suggests more stable platform for eIFiso4Ep∙VPg initiation complex formation, which was further supported by molecular docking analysis. Moreover, far-UV CD studies revealed that VPg formed complex with eIFiso4Ep with substantial change in the secondary structure. These results suggested that phosphorylation, not only reduced the energy barrier and dissociation rate but also enhanced binding rate, and an overall conformational change, which provides a more stable platform for efficient viral translation.

Development and Validation of AMBER-FB15-Compatible Force Field Parameters for Phosphorylated Amino Acids

Phosphorylation of select amino acid residues is one of the most common biological mechanisms for regulating protein structures and functions. While computational modeling can be used to explore the detailed structural changes associated with phosphorylation, most molecular mechanics force fields developed for the simulation of phosphoproteins have been noted to be inconsistent with experimental data. In this work, we parameterize force fields for the phosphorylated forms of the amino acids serine, threonine, and tyrosine using the ForceBalance software package with the goal of improving agreement with experiments for these residues. Our optimized force field, denoted as FB18, is parameterized using high-quality ab initio potential energy scans and is designed to be fully compatible with the AMBER-FB15 protein force field. When utilized in MD simulations together with the TIP3P-FB water model, we find that FB18 consistently enhances the prediction of experimental quantities such as ³J NMR couplings and intramolecular hydrogen-bonding propensities in comparison to previously published models. As was reported with AMBER-FB15, we also see improved agreement with the reference QM calculations in regions at and away from local minima. We thus believe that the FB18 parameter set provides a promising route for the further investigation of the varied effects of protein phosphorylation.

PTMdyna: exploring the influence of post-translation modifications on protein conformational dynamics

Protein post-translational modifications (PTM) play vital roles in cellular regulation, modulating functions by driving changes in protein structure and dynamics. Exploring comprehensively the influence of PTM on conformational dynamics can facilitate the understanding of the related biological function and molecular mechanism. Currently, a series of excellent computation tools have been designed to analyze the time-dependent structural properties of proteins. However, the protocol aimed to explore conformational dynamics of post-translational modified protein is still a blank. To fill this gap, we present PTMdyna to visually predict the conformational dynamics differences between unmodified and modified proteins, thus indicating the influence of specific PTM. PTMdyna exhibits an AUC of 0.884 tested on 220 protein-protein complex structures. The case of heterochromatin protein 1α complexed with lysine 9-methylated histone H3, which is critical for genomic stability and cell differentiation, was used to demonstrate its applicability. PTMdyna provides a reliable platform to predict the influence of PTM on protein dynamics, making it easier to interpret PTM functionality at the structure level. The web server is freely available at http://ccbportal.com/PTMdyna.

Changes in Protein Structural Motifs upon Post-Translational Modification in Kidney Cancer

Post-translational modification (PTM) leads to conformational changes in protein structure, modulates the biological function of proteins, and, consequently, changes the signature of metabolic transformations and the immune response in the body. Common PTMs are reversible and serve as a mechanism for modulating metabolic trans-formations in cells. It is likely that dysregulation of post-translational cellular signaling leads to abnormal proliferation and oncogenesis. We examined protein PTMs in the blood samples from patients with kidney cancer. Conformational changes in proteins after modification were analyzed. The proteins were analyzed using ultra-high resolution HPLC-MS/MS and structural analysis was performed with the AMBER and GROMACS software packages. Fifteen proteins containing PTMs were identified in blood samples from patients with kidney cancer. For proteins with PDB structures, a comparative analysis of the structural changes accompanying the modifications was performed. Results revealed that PTMs are localized in stable and compact space protein globule motifs that are exposed to a solvent. The phenomenon of modification is accompanied, as a rule, by an increase in the area available for the solvent of the modified amino acid residue and its active environment.

Towards an Understanding of Oxidative Damage in an α-L-Arabinofuranosidase of Trichoderma reesei: a Molecular Dynamics Approach

Trichoderma reesei is a "workhorse" fungus that produces glycosyl hydrolases (e.g., cellulases) at high titers for use in industrial bioprocessing. In this study, we focused on α-L-arabinofuranosidase, an enzyme important for the treatment of lignocellulosic biomass, but susceptible to oxidative damage that can occur during industrial processing. The molecular details that render this enzyme inactive have not yet been identified. To approach this issue, we used proteomics to identify amino acid residues that were oxidized after a relevant oxidative treatment (Fenton reaction). These oxidative modifications were included in the 3D protein structures, and using molecular dynamics simulations, we then studied the behaviors of non-modified and oxidized enzymes. These simulations showed significant alterations of the conformational stability of the protein when oxidized, as evidenced by changes in root mean square deviation (RMSD) and principal component analyses (PCA) trajectories. Likewise, enzyme-ligand interactions such as hydrogen bonds were greatly reduced in quantity and quality in the oxidized protein. Finally, free energy landscape plots showed that there was a more rugged energy surface in the oxidized protein, implying a less favorable reaction pathway. These results reveal the basis for loss of function in this carbohydrate active enzyme (CAZY) in the commercially relevant fungus T. reesei.

In Silico Molecular Docking and Molecular Dynamic Simulation Analysis of Phytochemicals From Indian Foods as Potential Inhibitors of SARS-CoV-2 RdRp and 3CLpro

With the current pandemic of the novel coronavirus disease 2019 (COVID-19) in hand, researchers around the globe are dexterously working to find the best suitable drug candidates and overcome vaccination-related challenges, to achieve efficient control over the second surge of COVID-19. The medical consultants time and again have been reiterating the need to abide by the precautionary steps to prevent the spread of the coronavirus by maintaining social distancing when outside, sanitizing hands regularly, and wearing masks and gloves. They also suggest taking a good and hygienic meal so as to boost immunity. Indians have an inborn nature of using natural spices, food, and medicines in their daily lives. Indian researchers have paid heed to deploy compounds from natural sources to explore potential antiviral agents against COVID-19 as the chances of acquiring side effects are perceived as less, and the efficacy of phytochemicals from medicinal plants is sometimes greater when compared to their synthetic counterparts. In the present study, we performed an in silico molecular docking and molecular dynamic simulation analysis of screened phytochemicals from a comprehensive list of Ayurvedic herbs/functional foods that are present in natural food products against key receptor proteins of severe acute respiratory syndrome coronavirus 2. We found that Aegle marmelos, Vetiveria zizanoides, Moringaolifera, and Punica granatum have antiviral potential to prevent coronavirus infection in the populace.

Perturbation Free-Energy Toolkit: An Automated Alchemical Topology Builder

Free-energy calculations play an important role in the application of computational chemistry to a range of fields, including protein biochemistry, rational drug design, or materials science. Importantly, the free-energy difference is directly related to experimentally measurable quantities such as partition and adsorption coefficients, water activity, and binding affinities. Among several techniques aimed at predicting free-energy differences, perturbation approaches, involving the alchemical transformation of one molecule into another through intermediate states, stand out as rigorous methods based on statistical mechanics. However, despite the importance of free-energy calculations, the applicability of the perturbation approaches is still largely impeded by a number of challenges, including the definition of the perturbation path, i.e., alchemical changes leading to the transformation of one molecule to the other. To address this, an automatic perturbation topology builder based on a graph-matching algorithm is developed, which can identify the maximum common substructure (MCS) of two or multiple molecules and provide the perturbation topologies suitable for free-energy calculations using the GROMOS and the GROMACS simulation packages. Various MCS search options are presented leading to alternative definitions of the perturbation pathway. Moreover, perturbation topologies generated using the default multistate MCS search are used to calculate the changes in free energy between lysine and its two post-translational modifications, 3-methyllysine and acetyllysine. The pairwise free-energy calculations performed on this test system led to a cycle closure of 0.5 ± 0.3 and 0.2 ± 0.2 kJ mol-1, with GROMOS and GROMACS simulation packages, respectively. The same relative free energies between the three states are obtained by employing the enveloping distribution sampling (EDS) approach when compared to the pairwise perturbations. Importantly, this toolkit is made available online as an open-source Python package (https://github.com/drazen-petrov/SMArt).

Systematic optimization of a fragment-based force field against experimental pure-liquid properties considering large compound families: application to oxygen and nitrogen compounds

The CombiFF approach is a workflow for the automated refinement of force-field parameters against experimental condensed-phase data, considering entire classes of organic molecules constructed using a fragment library via combinatorial isomer enumeration. One peculiarity of this approach is that it relies on an electronegativity-equalization scheme to account for induction effects within molecules, with values of the atomic hardness and electronegativity as electrostatic parameters, rather than the partial charges themselves. In a previous article [M. P. Oliveira, M. Andrey, S. R. Rieder, L. Kern, D. F. Hahn, S. Riniker, B. A. C. Horta and P. H. Hünenberger, J. Chem. Theory. Comput. 2020, 16, 7525], CombiFF was introduced and applied to calibrate a GROMOS-compatible united-atom force field for the saturated acyclic (halo-)alkane family. Here, this scheme is employed for the construction of a corresponding force field for saturated acyclic compounds encompassing eight common chemical functional groups involving oxygen and/or nitrogen atoms, namely: ether, aldehyde, ketone, ester, alcohol, carboxylic acid, amine, and amide. Monofunctional as well as homo-polyfunctional compounds are considered. A total of 1712 experimental liquid densities ρliq and vaporization enthalpies ΔHvap concerning 1175 molecules are used for the calibration (339 molecules) and validation (836 molecules) of the 102 non-bonded interaction parameters of the force field. Using initial parameter values based on the GROMOS 2016H66 parameter set, convergence is reached after five iterations. Given access to one processor per simulated system, this operation only requires a few days of wall-clock computing time. After optimization, the root-mean-square deviations from experiment are 29.9 (22.4) kg m-3 for ρliq and 4.1 (5.5) kJ mol-1 for ΔHvap for the calibration (validation) set. Thus, a very good level of agreement with experiment is achieved in terms of these two properties, although the errors are inhomogeneously distributed across the different chemical functional groups.

Impact of Reactive Oxygen and Nitrogen Species Produced by Plasma on Mdm2-p53 Complex

The study of protein–protein interactions is of great interest. Several early studies focused on the murine double minute 2 (Mdm2)–tumor suppressor protein p53 interactions. However, the effect of plasma treatment on Mdm2 and p53 is still absent from the literature. This study investigated the structural changes in Mdm2, p53, and the Mdm2–p53 complex before and after possible plasma oxidation through molecular dynamic (MD) simulations. MD calculation revealed that the oxidized Mdm2 bounded or unbounded showed high flexibility that might increase the availability of tumor suppressor protein p53 in plasma-treated cells. This study provides insight into Mdm2 and p53 for a better understanding of plasma oncology.

Structural basis for the hyperthermostability of an archaeal enzyme induced by succinimide formation

Stability of proteins from hyperthermophiles (organisms existing under boiling water conditions) enabled by a reduction of conformational flexibility is realized through various mechanisms. A succinimide (SNN) arising from the post-translational cyclization of the side chains of aspartyl/asparaginyl residues with the backbone amide -NH of the succeeding residue would restrain the torsion angle Ψ and can serve as a new route for hyperthermostability. However, such a succinimide is typically prone to hydrolysis, transforming to either an aspartyl or β-isoaspartyl residue. Here, we present the crystal structure of Methanocaldococcus jannaschii glutamine amidotransferase and, using enhanced sampling molecular dynamics simulations, address the mechanism of its increased thermostability, up to 100°C, imparted by an unexpectedly stable succinimidyl residue at position 109. The stability of SNN109 to hydrolysis is seen to arise from its electrostatic shielding by the side-chain carboxylate group of its succeeding residue Asp110, as well as through n → π^∗ interactions between SNN109 and its preceding residue Glu108, both of which prevent water access to SNN. The stable succinimidyl residue induces the formation of an α-turn structure involving 13-atom hydrogen bonding, which locks the local conformation, reducing protein flexibility. The destabilization of the protein upon replacement of SNN with a Φ-restricted prolyl residue highlights the specificity of the succinimidyl residue in imparting hyperthermostability to the enzyme. The conservation of the succinimide-forming tripeptide sequence (E(N/D)(E/D)) in several archaeal GATases strongly suggests an adaptation of this otherwise detrimental post-translational modification as a harbinger of thermostability.

In Silico and Electrochemical Studies for a ZnO-CuO-Based Immunosensor for Sensitive and Selective Detection of E. coli

Escherichia coli is a harmful Gram-negative bacterium commonly found in the gut of warm-blooded organisms and affects millions of people annually worldwide. In this study, we have synthesized a ZnO-CuO nanocomposite (NC) by a co-precipitation method and characterized the as-synthesized NC using FTIR spectroscopy, XRD, Raman spectroscopy, and FESEM techniques. To fabricate the immunosensor, the ZnO-CuO NC composite was screen-printed on gold-plated electrodes followed by physisorption of the anti-LPS E. coli antibody. The biosensor was optimized for higher specificity and sensitivity. The immunosensor exhibited a high sensitivity (11.04 μA CFU mL-1) with a low detection limit of 2 CFU mL-1 with a redox couple. The improved performance of the immunosensor is attributed to the synergistic effect of the NC and the antilipopolysaccharide antibody against E. coli. The selectivity studies were also carried out with Staphylococcus aureus to assess the specificity of the immunosensor. Testing in milk samples was done by spiking the milk samples with different concentrations of E. coli to check the potential of this immunosensor. We further checked the affinity between ZnO-CuO NC with E. coli LPS and the anti-LPS antibody using molecular docking studies. Atomic charge computation and interaction analyses were performed to support our hypothesis. Our results discern that there is a strong correlation between molecular docking studies and electrochemical characterization. The interaction analysis further displays the strong affinity between the antibody-LPS complex when immobilized with a nanoparticle composite (ZnO-CuO).

Additive CHARMM36 Force Field for Nonstandard Amino Acids

Nonstandard amino acids are both abundant in nature, where they play a key role in various cellular processes, and can be synthesized in laboratories, for example, for the manufacture of a range of pharmaceutical agents. In this work, we have extended the additive all-atom CHARMM36 and CHARMM General force field (CGenFF) to a large set of 333 nonstandard amino acids. These include both amino acids with nonstandard side chains, such as post-translationally modified and artificial amino acids, as well as amino acids with modified backbone groups, such as chromophores composed of several amino acids. Model compounds representative of the nonstandard amino acids were parametrized for protonation states that are likely at the physiological pH of 7 and, for some more common residues, in both d- and l-stereoisomers. Considering all protonation, tautomeric, and stereoisomeric forms, a total of 406 nonstandard amino acids were parametrized. Emphasis was placed on the quality of both intra- and intermolecular parameters. Partial charges were derived using quantum mechanical (QM) data on model compound dipole moments, electrostatic potentials, and interactions with water. Optimization of all intramolecular parameters, including torsion angle parameters, was performed against information from QM adiabatic potential energy surface (PES) scans. Special emphasis was put on the quality of terms corresponding to PES around rotatable dihedral angles. Validation of the force field was based on molecular dynamics simulations of 20 protein complexes containing different nonstandard amino acids. Overall, the presented parameters will allow for computational studies of a wide range of proteins containing nonstandard amino acids, including natural and artificial residues.