A systematic framework for molecular dynamics simulations of protein post-translational modifications

PCSK9 inhibitors as safer therapeutics for atherosclerotic cardiovascular disease (ASCVD): Pharmacophore design and molecular dynamics analysis

Cardiovascular disorders are still challenging and are among the deadly diseases. As a major risk factor for atherosclerotic cardiovascular disease, dyslipidemia, and high low-density lipoprotein cholesterol in particular, can be prevented primary and secondary by lipid-lowering medications. Therefore, insights are still needed into designing new drugs with minimal side effects. Proprotein convertase subtilisin/kexin 9 (PCSK9) enzyme catalyses protein-protein interactions with low-density lipoprotein, making it a critical target for designing promising inhibitors compared to statins. Therefore, we screened for potential compounds using a redesigned PCSK9 conformational behaviour to search for a significantly extensive chemical library and investigated the inhibitory mechanisms of the final compounds using integrated computational methods, from ligand essential functional group screening to all-atoms MD simulations and MMGBSA-based binding free energy. The inhibitory mechanisms of the screened compounds compared with the standard inhibitor. K31 and K34 molecules showed stronger interactions for PCSK9, having binding energy (kcal/mol) of -33.39 and -63.51, respectively, against -27.97 of control. The final molecules showed suitable drug-likeness, non-mutagenesis, permeability, and high solubility values. The C-α atoms root mean square deviation and root mean square fluctuation of the bound-PCSK9 complexes showed stable and lower fluctuations compared to apo PCSK9. The findings present a model that unravels the mechanism by which the final molecules proposedly inhibit the PCSK9 function and could further improve the design of novel drugs against cardiovascular diseases.

Protein Oxidative Modifications in Neurodegenerative Diseases: From Advances in Detection and Modelling to Their Use as Disease Biomarkers

Oxidation-reduction post-translational modifications (redox-PTMs) are chemical alterations to amino acids of proteins. Redox-PTMs participate in the regulation of protein conformation, localization and function, acting as signalling effectors that impact many essential biochemical processes in the cells. Crucially, the dysregulation of redox-PTMs of proteins has been implicated in the pathophysiology of numerous human diseases, including neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease. This review aims to highlight the current gaps in knowledge in the field of redox-PTMs biology and to explore new methodological advances in proteomics and computational modelling that will pave the way for a better understanding of the role and therapeutic potential of redox-PTMs of proteins in neurodegenerative diseases. Here, we summarize the main types of redox-PTMs of proteins while providing examples of their occurrence in neurodegenerative diseases and an overview of the state-of-the-art methods used for their detection. We explore the potential of novel computational modelling approaches as essential tools to obtain insights into the precise role of redox-PTMs in regulating protein structure and function. We also discuss the complex crosstalk between various PTMs that occur in living cells. Finally, we argue that redox-PTMs of proteins could be used in the future as diagnosis and prognosis biomarkers for neurodegenerative diseases.

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.

[Modern approaches to studying the molecular mechanisms of lung functioning in normal and pathological conditions]

Problems with breathing and lung function are caused by the development of various lung diseases associated with lifestyle, harmful environmental factors and genetic predisposition. Knowledge of the molecular mechanisms of the development of the pathological process will allow on time identification of the disease or the development of targeted therapy. The article provides an overview of modern methods that make it possible to most accurately reproduce the structural, functional and mechanical properties of the lung (organ-on-a-chip), to perform non-invasive molecular studies of biomarkers of bronchopulmonary pathology using saliva diagnostics, as well as using DNA and RNA aptamers, verify tumor markers in biological samples of human tissue. Analysis of alterations in the pattern of protein glycosylation using glycodiagnostic methods makes it possible to detect lung cancer in the early stages.

Top-Down Machine Learning of Coarse-Grained Protein Force Fields

Developing accurate and efficient coarse-grained representations of proteins is crucial for understanding their folding, function, and interactions over extended time scales. Our methodology involves simulating proteins with molecular dynamics and utilizing the resulting trajectories to train a neural network potential through differentiable trajectory reweighting. Remarkably, this method requires only the native conformation of proteins, eliminating the need for labeled data derived from extensive simulations or memory-intensive end-to-end differentiable simulations. Once trained, the model can be employed to run parallel molecular dynamics simulations and sample folding events for proteins both within and beyond the training distribution, showcasing its extrapolation capabilities. By applying Markov state models, native-like conformations of the simulated proteins can be predicted from the coarse-grained simulations. Owing to its theoretical transferability and ability to use solely experimental static structures as training data, we anticipate that this approach will prove advantageous for developing new protein force fields and further advancing the study of protein dynamics, folding, and interactions.

Cysteine carbonylation with reactive carbonyl species from lipid peroxidation induce local structural changes on thioredoxin active site

Protein oxidative modifications with reactive carbonyl species (RCS) is directly linked to metabolic processes in premature aging, cancer, neurodegenerative and infectious diseases. RCS as 4-Hydroxy-2-nonal (HNE), 4-Hydroxy-2-hexenal (HHE), 4-Oxo-2-nonenal (ONE) and Malondialdehyde (MDA) attack nucleophilic amino acids residues forming irreversible adducts with proteins as Thioredoxins (Trx). This is a class of small thiol oxide-reductases playing a key role in redox signaling and oxidative stress responses in mammals. Although proteomic studies have identified to Cys-32 residue as a target of HNE attack that inhibit its enzymatic activity, how this carbonylation affects its structure and dynamic behavior at the atomic level is unknown. Even more, the molecular bases for the atomistic behavior of these modified proteins have not been completely understood. We present molecular dynamics simulations of Trx-modified with four different RCS to analyze its global and local structural effects. For this, parameters supported in the AMBER force fields were built and validated for three non-natural cysteine residues modified with HHE, ONE and MDA. Results obtained showed a slight change in the global conformational stability of Trx modified with HNE and MDA, establishing that all modified proteins presented local regions of high mobility in the modified catalytic site and some regions far from the modification area. In addition, essential domain movement modes evidences that proteins modified with the RCS assayed induce changes in conformational flexibility. Finally, these data showed that the given conformational changes did not caused global changes in proteins but rather localized changes in particular regions.

Towards targeting EGFR and COX-2 inhibitors: comprehensive computational studies on the role of chlorine group in novel thienyl-pyrazoline derivative

To enhance the effectiveness of chemotherapy and overcome resistance, scientists must develop novel drugs or scaffolds that have a combined effect, such as the inhibition of EGFR and COX-2. This research employed virtual screening techniques, such as docking, and dynamics simulation, to predict chlorinated thienyl-pyrazoline derivatives that inhibit these proteins. The study proposed eleven (11) ligands with binding energies ranging from -7.8 kcal/mol to -8.7 kcal/mol for EGFR and -6.4 kcal/mol to -8.4 kcal/mol for COX-2. Ligands P1 and P11 exhibited the highest binding affinity for both proteins. The results of RMSD, RMSF, RoG, SASA the number of hydrogen bonds, and BAR free binding energy demonstrated the good stability of ligands P1 and P11 when binding to both proteins over 180 ns simulations. In addition, the absorption, distribution, metabolism, excretion, and toxicity properties of the selected ligands were assessed to predict their toxicity and drug likeliness. Based on the results, these compounds can be proposed for further synthesis and in vitro studies.Communicated by Ramaswamy H. Sarma.

The protein dynamics of bovine and caprine β-lactoglobulin differ as a function of pH

The properties of milk proteins differ between mammalian species. β-Lactoglobulin (βlg) proteins from caprine and bovine milk are sequentially and structurally highly similar, yet their physicochemical properties differ, particularly in response to pH. To resolve this conundrum, we compared the dynamics of both the monomeric and dimeric states for each homologue at pH 6.9 and 7.5 using hydrogen/deuterium exchange experiments. At pH 7.5, the rate of exchange is similar across both homologues, but at pH 6.9 the dimeric states of the bovine βlg B variant homologue have significantly more conformational flexibility compared with caprine βlg. Molecular dynamics simulations provide a mechanistic rationale for the experimental observations, revealing that variant-specific substitutions encode different conformational ensembles with different dynamic properties consistent with the hydrogen/deuterium exchange experiments. Understanding the dynamic differences across βlg homologues is essential to understand the different responses of these milks to processing, human digestion, and differences in immunogenicity.

TRAP1 S-nitrosylation as a model of population-shift mechanism to study the effects of nitric oxide on redox-sensitive oncoproteins

S-nitrosylation is a post-translational modification in which nitric oxide (NO) binds to the thiol group of cysteine, generating an S-nitrosothiol (SNO) adduct. S-nitrosylation has different physiological roles, and its alteration has also been linked to a growing list of pathologies, including cancer. SNO can affect the function and stability of different proteins, such as the mitochondrial chaperone TRAP1. Interestingly, the SNO site (C501) of TRAP1 is in the proximity of another cysteine (C527). This feature suggests that the S-nitrosylated C501 could engage in a disulfide bridge with C527 in TRAP1, resembling the well-known ability of S-nitrosylated cysteines to resolve in disulfide bridge with vicinal cysteines. We used enhanced sampling simulations and in-vitro biochemical assays to address the structural mechanisms induced by TRAP1 S-nitrosylation. We showed that the SNO site induces conformational changes in the proximal cysteine and favors conformations suitable for disulfide bridge formation. We explored 4172 known S-nitrosylated proteins using high-throughput structural analyses. Furthermore, we used a coarse-grained model for 44 protein targets to account for protein flexibility. This resulted in the identification of up to 1248 proximal cysteines, which could sense the redox state of the SNO site, opening new perspectives on the biological effects of redox switches. In addition, we devised two bioinformatic workflows ( https://github.com/ELELAB/SNO_investigation_pipelines ) to identify proximal or vicinal cysteines for a SNO site with accompanying structural annotations. Finally, we analyzed mutations in tumor suppressors or oncogenes in connection with the conformational switch induced by S-nitrosylation. We classified the variants as neutral, stabilizing, or destabilizing for the propensity to be S-nitrosylated and undergo the population-shift mechanism. The methods applied here provide a comprehensive toolkit for future high-throughput studies of new protein candidates, variant classification, and a rich data source for the research community in the NO field.

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.

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.

Bioactive Nonthermal Biocompatible Plasma Enhances Migration on Human Gingival Fibroblasts

This study hypothesizes that the application of low-dose nonthermal biocompatible dielectric barrier discharge plasma (DBD-NBP) to human gingival fibroblasts (HGFs) will inhibit colony formation but not cell death and induce matrix metalloproteinase (MMP) expression, extracellular matrix (ECM) degradation, and subsequent cell migration, which can result in enhanced wound healing. HGFs treated with plasma for 3 min migrate to each other across the gap faster than those in the control and 5-min treatment groups on days 1 and 3. The plasma-treated HGFs show significantly high expression levels of the cell cycle arrest-related p21 gene and enhanced MMP activity. Focal adhesion kinase (FAK) mediated attenuation of wound healing or actin cytoskeleton rearrangement, and plasma-mediated reversal of this attenuation support the migratory effect of DBD-NBP. Further, this work performs computer simulations to investigate the effect of oxidation on the stability and conformation of the catalytic kinase domain (KD) of FAK. It is found that the oxidation of highly reactive amino acids (AAs) Cys427, Met442, Cys559, Met571, Met617, and Met643 changes the conformation and increases the structural flexibility of the FAK protein and thus modulates its function and activity. Low-dose DBD-NBP-induces host cell cycle arrest, ECM breakdown, and subsequent migration, thus contributing to the enhanced wound healing process.

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.

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.

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.

Engineering of enzymes using non-natural amino acids

In enzyme engineering, the main targets for enhancing properties are enzyme activity, stereoselective specificity, stability, substrate range, and the development of unique functions. With the advent of genetic code extension technology, non-natural amino acids (nnAAs) are able to be incorporated into proteins in a site-specific or residue-specific manner, which breaks the limit of 20 natural amino acids for protein engineering. Benefitting from this approach, numerous enzymes have been engineered with nnAAs for improved properties or extended functionality. In this review, we focus on applications and strategies for using nnAAs in enzyme engineering. Notably, approaches to computational modelling of enzymes with nnAAs are also addressed. Finally, we discuss the bottlenecks that currently need to be addressed in order to realise the broader prospects of this genetic code extension technique.

Leveraging Protein Dynamics to Identify Functional Phosphorylation Sites using Deep Learning Models

Accurate prediction of post-translational modifications (PTMs) is of great significance in understanding cellular processes, by modulating protein structure and dynamics. Nowadays, with the rapid growth of protein data at different "omics" levels, machine learning models largely enriched the prediction of PTMs. However, most machine learning models only rely on protein sequence and little structural information. The lack of the systematic dynamics analysis underlying PTMs largely limits the PTM functional predictions. In this research, we present two dynamics-centric deep learning models, namely, cDL-PAU and cDL-FuncPhos, by incorporating sequence, structure, and dynamics-based features to elucidate the molecular basis and underlying functional landscape of PTMs. cDL-PAU achieved satisfactory area under the curve (AUC) scores of 0.804-0.888 for predicting phosphorylation, acetylation, and ubiquitination (PAU) sites, while cDL-FuncPhos achieved an AUC value of 0.771 for predicting functional phosphorylation (FuncPhos) sites, displaying reliable improvements. Through a feature selection, the dynamics-based coupling and commute ability show large contributions in discovering PAU sites and FuncPhos sites, suggesting the allosteric propensity for important PTMs. The application of cDL-FuncPhos in three oncoproteins not only corroborates its strong performance in FuncPhos prioritization but also gains insight into the physical basis for the functions. The source code and data set of cDL-PAU and cDL-FuncPhos are available at https://github.com/ComputeSuda/PTM_ML.

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.

Molecular Dynamics Study of Citrullinated Proteins Associated with the Development of Rheumatoid Arthritis

Biological activity regulation by protein post-translational modification (PTM) is critical for cell function, development, differentiation, and survival. Dysregulation of PTM proteins is present in various pathological conditions, including rheumatoid arthritis (RA). RA is a systemic autoimmune disease that primarily affects joints, and there are three main types of protein PTMs associated with the development of this disease, namely, glycosylation, citrullination, and carbamylation. Glycosylation is important for the processing and presentation of antigen fragments on the cell surface and can modulate immunoglobulin activity. The citrullination of autoantigens is closely associated with RA, as evidenced by the presence of antibodies specific to citrullinated proteins in the serum of patients. Carbamylation and dysregulation have recently been associated with RA development in humans.In this study, we performed an overview analysis of proteins with post-translational modifications associated with the development of RA adverted in peer-reviewed scientific papers for the past 20 years. As a result of the search, a list of target proteins and corresponding amino acid sequences with PTM in RA was formed. Structural characteristics of the listed modified proteins were extracted from the Protein Data Bank. Then, molecular dynamics experiments of intact protein structures and corresponding structures with PTMs were performed regarding structures in the list announced in the ProtDB service. This study aimed to conduct a molecular dynamics study of intact proteins and proteins, including post-translational modification and protein citrullination, likely associated with RA development. We observed another exhibition of the fundamental physics concept, symmetry, at the submolecular level, unveiled as the autonomous repetitions of outside the protein structural motif performance globule corresponding to those in the whole protein molecule.

Molecular Modeling Insights into Upadacitinib Selectivity upon Binding to JAK Protein Family

Rheumatoid arthritis (RA) is a chronic disease characterized by bone joint damage and incapacitation. The mechanism underlying RA pathogenesis is autoimmunity in the connective tissue. Cytokines play an important role in the human immune system for signal transduction and in the development of inflammatory responses. Janus kinases (JAK) participate in the JAK/STAT pathway, which mediates cytokine effects, in particular interleukin 6 and IFNγ. The discovery of small molecule inhibitors of the JAK protein family has led to a revolution in RA therapy. The novel JAK inhibitor upadacitinib (Rinvoq^TM) has a higher selectivity for JAK1 compared to JAK2 and JAK3 in vivo. Currently, details on the molecular recognition of JAK1 by upadacitinib are not available. We found that characteristics of hydrogen bond formation with the glycine loop and hinge in JAKs define the selectivity. Our molecular modeling study could provide insight into the drug action mechanism and pharmacophore model differences in JAK isoforms.

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.

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).

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.

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.

Lipid Oxidation: Role of Membrane Phase-Separated Domains

Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are "hot spots" for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

KEAP1 Cancer Mutants: A Large-Scale Molecular Dynamics Study of Protein Stability

We have performed 280 μs of unbiased molecular dynamics (MD) simulations to investigate the effects of 12 different cancer mutations on Kelch-like ECH-associated protein 1 (KEAP1) (G333C, G350S, G364C, G379D, R413L, R415G, A427V, G430C, R470C, R470H, R470S and G476R), one of the frequently mutated proteins in lung cancer. The aim was to provide structural insight into the effects of these mutants, including a new class of ANCHOR (additionally NRF2-complexed hypomorph) mutant variants. Our work provides additional insight into the structural dynamics of mutants that could not be analyzed experimentally, painting a more complete picture of their mutagenic effects. Notably, blade-wise analysis of the Kelch domain points to stability as a possible target of cancer in KEAP1. Interestingly, structural analysis of the R470C ANCHOR mutant, the most prevalent missense mutation in KEAP1, revealed no significant change in structural stability or NRF2 binding site dynamics, possibly indicating an covalent modification as this mutant's mode of action.

Molecular dynamics shows complex interplay and long-range effects of post-translational modifications in yeast protein interactions

Post-translational modifications (PTMs) play a vital, yet often overlooked role in the living cells through modulation of protein properties, such as localization and affinity towards their interactors, thereby enabling quick adaptation to changing environmental conditions. We have previously benchmarked a computational framework for the prediction of PTMs’ effects on the stability of protein-protein interactions, which has molecular dynamics simulations followed by free energy calculations at its core. In the present work, we apply this framework to publicly available data on Saccharomyces cerevisiae protein structures and PTM sites, identified in both normal and stress conditions. We predict proteome-wide effects of acetylations and phosphorylations on protein-protein interactions and find that acetylations more frequently have locally stabilizing roles in protein interactions, while the opposite is true for phosphorylations. However, the overall impact of PTMs on protein-protein interactions is more complex than a simple sum of local changes caused by the introduction of PTMs and adds to our understanding of PTM cross-talk. We further use the obtained data to calculate the conformational changes brought about by PTMs. Finally, conservation of the analyzed PTM residues in orthologues shows that some predictions for yeast proteins will be mirrored to other organisms, including human. This work, therefore, contributes to our overall understanding of the modulation of the cellular protein interaction networks in yeast and beyond.

Proteins are a diverse set of biological molecules responsible for numerous functions within cells, such as obtaining energy from food or transport of small molecules, and many processes rely on interactions of specific proteins. Moreover, a single protein may acquire different roles depending on cellular requirements and as a response to changes in the environment. A commonly used way to quickly change protein’s function or activity is by introducing small chemical modifications on specific locations within the protein. These modifications can cause the protein to interact in a more or less stable way with other proteins. We have previously developed a computational pipeline for predicting the effect of modifications on interactions of proteins, and in this work we apply it to all yeast proteins with known structures. We find differences in effects on the binding for different types of modifications. Importantly, we demonstrate that the modifications far from the interaction interface also significantly contribute to binding due to their impact on protein’s shape, which is often neglected by other methods. This work contributes to our understanding of the modulation of protein interactions in yeast due to modifications, while our widely applicable method will allow similar investigations in other organisms.

Molecular Dynamic Simulations Suggest That Metabolite-Induced Post-Translational Modifications Alter the Behavior of the Fibrinogen Coiled-Coil Domain

Fibrinogen is an abundant blood plasma protein that, inter alia, participates in blood coagulation. It polymerizes to form a fibrin clot that is among the major components of the thrombus. Fibrinogen reactions with various reactive metabolites may induce post-translational modifications (PTMs) into the protein structure that affect the architecture and properties of fibrin clots. We reviewed the previous literature to find the positions of PTMs of fibrinogen. For 7 out of 307 reported PTMs, we used molecular dynamics simulations to characterize their effect on the behavior of the fibrinogen coiled-coil domain. Interactions of the γ-coil with adjacent chains give rise to π-helices in Aα and Bβ chains of even unmodified fibrinogen. The examined PTMs suppress fluctuations of the γ-coil, which may affect the fibrinolysis and stiffness of the fibrin fibers. Citrullination of AαR104 and oxidations of γP70 and γP76 to glutamic semialdehyde unfold the α-helical structure of Aα and Bβ chains. Oxidation of γM78 to methionine sulfoxide induces the formation of an α-helix in the γ-coil region. Our findings suggest that certain PTMs alter the protein secondary structure. Thus, the altered protein structure may indicate the presence of PTMs in the molecule and consequently of certain metabolites within the system.

Oxidative damage to hyaluronan-CD44 interactions as an underlying mechanism of action of oxidative stress-inducing cancer therapy

Multiple cancer therapies nowadays rely on oxidative stress to damage cancer cells. Here we investigated the biological and molecular effect of oxidative stress on the interaction between CD44 and hyaluronan (HA), as interrupting their binding can hinder cancer progression. Our experiments demonstrated that the oxidation of HA decreased its recognition by CD44, which was further enhanced when both CD44 and HA were oxidized. The reduction of CD44-HA binding negatively affected the proliferative state of cancer cells. Our multi-level atomistic simulations revealed that the binding free energy of HA to CD44 decreased upon oxidation. The effect of HA and CD44 oxidation on CD44-HA binding was similar, but when both HA and CD44 were oxidized, the effect was much larger, in agreement with our experiments. Hence, our experiments and computations support our hypothesis on the role of oxidation in the disturbance of CD44-HA interaction, which can lead to the inhibition of proliferative signaling pathways inside the tumor cell to induce cell death.

Cell division requires RNA eviction from condensing chromosomes

During mitosis, the genome is transformed from a decondensed, transcriptionally active state to a highly condensed, transcriptionally inactive state. Mitotic chromosome reorganization is marked by the general attenuation of transcription on chromosome arms, yet how the cell regulates nuclear and chromatin-associated RNAs after chromosome condensation and nuclear envelope breakdown is unknown. SAF-A/hnRNPU is an abundant nuclear protein with RNA-to-DNA tethering activity, coordinated by two spatially distinct nucleic acid–binding domains. Here we show that RNA is evicted from prophase chromosomes through Aurora-B–dependent phosphorylation of the SAF-A DNA-binding domain; failure to execute this pathway leads to accumulation of SAF-A–RNA complexes on mitotic chromosomes, defects in metaphase chromosome alignment, and elevated rates of chromosome missegregation in anaphase. This work reveals a role for Aurora-B in removing chromatin-associated RNAs during prophase and demonstrates that Aurora-B–dependent relocalization of SAF-A during cell division contributes to the fidelity of chromosome segregation.

Spectroscopy, Dynamics, and Hydration of S-Nitrosylated Myoglobin

S-Nitrosylation, the covalent addition of NO to the thiol side chain of cysteine, is an important post-transitional modification that can alter the function of various proteins. The structural dynamics and vibrational spectroscopy of S-nitrosylation in the condensed phase are investigated for the methyl-capped cysteine model system and for myoglobin. Using conventional point charge and physically more realistic multipolar force fields for the -SNO group, it is found that the SN- and NO-stretch and the SNO-bend vibrations can be located and distinguished from the other protein modes for simulations of MbSNO at 50 K. The finding of stable cis- and trans-MbSNO agrees with experimental findings on other proteins as is the observation of buried -SNO. For MbSNO the observed relocation of the EF loop in the simulations by ∼3 Å is consistent with the available X-ray structure, and the conformations adopted by the -SNO label are in good overall agreement with the X-ray structure. Despite the larger size of the -SNO group compared with -SH, MbSNO recruits more water molecules in the first two hydration shells due to stronger electrostatic interactions. Similarly, when comparing the hydration between the A- and H-helices, they differ by up to 30% between WT and MbSNO. This suggests that local hydration can also be significantly modulated through nitrosylation.

Dimeric phosphorylation of glyoxalase I alters its symmetry and substrate binding mechanism: simulation studies

Glyoxalase I (GLO1) is a dimeric esterase of the glyoxalase system. Phosphorylation of the residue T106 has been found to inhibit GLO1 activity, and contribute to the onset of oxidative stress and cellular damage. This research uses multiple molecular dynamics simulations and automated docking of both GLO1 and dimerically phosphorylated GLO1 (p2-GLO1) to predict the initial structural differences induced by phosphorylation, and their interaction with the intermediate substrate Hemimercaptal. This research indicates that immediately following phosphorylation, GLO1 exhibits reduced sphericity, partly caused by outward splaying of the loop region surrounding T106. Phosphorylation induces enhanced concerted motions in the loop composed of residues immediately surrounding T106, which are correlated with motions at the active site pocket at the distant, opposite end of the dimer. These T106 region loop motions result in the distortion of the shape of the active site, and potentially alter its accessibility. Phosphorylation alters the manner in which GLO1 interacts with Hemimercaptal. For GLO1, Hemimercaptal is predicted to bind to T106, which we propose constitutes a novel, highly accessible 'capture site' responsible for initial contact with the substrate. In contrast, for p2-GLO1, Hemimercaptal is unable to bind favourably to (phosphorylated) position T106, suggesting that this proposed transient 'capture site' is abolished upon phosphorylation of GLO1. Hence, a novel physiological role is here proposed for the known essential GLO1 residue T106. These results may further contribute to understanding the inhibition mechanism of GLO1 upon phosphorylation.Communicated by Ramaswamy H. Sarma.

Development and Testing of Force Field Parameters for Phenylalanine and Tyrosine Derivatives

Theoretical analyses are valuable for the exploration of the effects of unnatural amino acids on enzyme functions; however, many necessary parameters for unnatural amino acids remain lacking. In this study, we developed and tested force field parameters compatible with Amber ff14SB for 18 phenylalanine and tyrosine derivatives. The charge parameters were derived from ab initio calculations using the RESP fitting approach and then adjusted to reproduce the benchmark relative energies (at the MP2/TZ level) of the α- and β-backbones for each unnatural amino acid dipeptide. The structures optimized under the proposed force field parameters for the 18 unnatural amino acid dipeptides in both the α- and β-backbone forms were in good agreement with their QM structures, as the average RMSD was as small as 0.1 Å. The force field parameters were then tested in their application to seven proteins containing unnatural amino acids. The RMSDs of the simulated configurations of these unnatural amino acids were approximately 1.0 Å compared with those of the crystal structures. The vital interactions between proteins and unnatural amino acids in five protein–ligand complexes were also predicted using MM/PBSA analysis, and they were largely consistent with experimental observations. This work will provide theoretical aid for drug design involving unnatural amino acids.

Exploring the structural significance of molecular switch mechanism alias motif phosphorylation in Wnt/β-catenin and their crucial role in triple-negative breast cancer

β-Catenin, a key transcriptional factor involved in the canonical Wnt signaling pathway, is regulated by a cascade of phosphorylations and plays a major role in the progression of triple-negative breast cancer (TNBC). However, the phosphorylation induced conformational changes in a β-Catenin is still poorly understood. Hence, we adopted a conventional molecular dynamics approach to study phosphorylations present in a sequence motif Ser ^552 675 and Tyr⁶⁷⁰ of the β-Catenin domain and analyzed in terms of structural transitions, bond formation, and folding-misfolding conformations. Our results unveil the β-Catenin linear motif 549-555 (RRTSMGG) of armadillo repeats domain prefers order to disorder state. In contrast, helix C associated with 670-678 (YKKRLSVEL) motif prefers disorder to order upon phosphorylation of Ser ^552 675 and Tyr⁶⁷⁰. In addition, the crucial secondary structural transition from α-helix to coil induced by phospho Ser⁵⁵² and phospho Tyr⁶⁷⁰ of β-Catenin ARM domain connecting helix C modifies conformational diversity and binding affinities of the complex interaction in functional regulation significantly. Moreover, the post phosphorylation disrupted the hydrogen bond interactions (Ser552-Arg549, Arg550-Asp546 and Ser675-Lys672) and abolished the residual alliance with hydrophobic interactions (Tyr670-Leu674) that easily interrupt in secondary structure packing as well as folding conformations connecting ARM and helix C (R10, 12 & R1C) compared to unphosphorylation. Our integrated computational analysis may help in shedding light on understanding the induced folding and unfolding pattern due to motif phosphorylations. Overall, our results provide an atomistic structural description of the way phosphorylation facilitates conformational and dynamic changes in β-Catenin, a fundamental molecular switch mechanism in triple-negative breast cancer pathogenesis.

Improved GROMOS 54A7 Charge Sets for Phosphorylated Tyr, Ser, and Thr to Deal with pH-Dependent Binding Phenomena

Phosphorylation is a ubiquitous post-translational modification in proteins, and the phosphate group is present constitutively or transiently in most biological building blocks. These phosphorylated biomolecules are involved in many high-affinity binding/unbinding events that rely predominantly on electrostatic interactions. To build accurate models of these molecules, we need an improved description of the atomic partial charges for all relevant protonation states. In this work, we showed that the commonly used protocols to derive atomic partial charges using well-solvated molecules are inadequate to model the protonation equilibria in binding events. We introduced a protocol based on PB/MC calculations with a single representative conformation (of both protonation states) and used the resulting pK_a estimations to help manually curate the atomic partial charges. The final charge set, which is fully compatible with the GROMOS 54A7 force field, proved to be very effective in modeling the protonation equilibrium in different phosphorylated peptides in the free (tetrapeptides, pentapeptides, and pY1021) and protein-complexed forms (pY1021/PLC-γ1 complex). This was particularly important in the case of the pY1021 bound to the SH2 domain of PLC-γ1, where only our curated charge set captured the correct protonation equilibrium at the neutral to slightly acidic pH range. The binding/unbinding phenomena in that pH range are biologically relevant, and to improve our models, we need to go beyond the commonly used protocols and obtain revised force field parameters for these molecules.

Molecular modelling of sorption processes of a range of diverse small organic molecules in Leonardite humic acid

Soil organic matter (SOM) is abundant in the environment and plays an important role in several biogeochemical processes, including microbial activity, soil aggregation, plant growth and carbon storage. One of its key functions is the retention and release of various chemical compounds, primarily governed by the sorption process, which strongly affects the environmental fate of nutrients and pollutants. Sorption largely depends on the composition of SOM, as well as its structure, dynamics and the thermodynamic conditions. Although several approaches are available, experimental characterization of sorption mechanisms is not easy. Computational models for predicting sorption coefficients often require a wealth of experimental data for training and are only applicable to compounds and conditions related to the training dataset. Here, we use molecular dynamics (MD) simulations to study the sorption of a range of small organic compounds. As a model SOM system we use the standard Leonardite humic acid (LHA) sample, which physicochemical properties have recently been characterized computationally in detail. This model allowed us to estimate sorption propensities of the systems at two different hydration levels (water activities close to 0 and 1), showing a remarkable correlation with experimental data. Importantly, this molecular modelling approach based on perturbation free-energy calculations is rigorously derived from statistical thermodynamics and requires no experimental sorption data for training. It is therefore in principle applicable to any SOM model or thermodynamic condition. Moreover, the power of MD simulations to provide high-resolution insight into atomistic and molecular interactions was employed to explore how sorbate molecules associate with the LHA matrix and which contacts they form. The heteroatoms of both sorbate and sorbent play an important role and water molecules are identified as further key players in facilitating the sorption process.

HIGHLIGHTS

Modelling of the sorption processes in soil organic matter at atomistic level.Rigorous, physics-based approach applicable to a range of SOM systems and conditions.Remarkable level of matching with experimental data with additional insight into the molecular mechanism.Interactions between the sorbate and local environment strongly affects the sorption process.

Evidence for NAD+-dependent histone dynamics and tunneling associated conformational transitions in circadian deacetylase SIRT1

Circadian rhythm is a biological cycle that is involved in all processes over 24 h day and night period. Sirtuin 1 (SIRT1) is a 747 amino acid-long class III Nicotinamide adenine dinucleotide (NAD⁺)-dependent histone that acts as a circadian deacetylase. Here we present a detailed in-silico analysis to address comparative structure-function relationship and interaction pattern of SIRT1-NAD⁺/Zn⁺² and SIRT1^NAD+/Zn+2-acetylated histone H4 (H4^KAC16) complexes. MD-based ensemble analysis suggested an overall loss of helical content (21.144-17.230%) in H4^KAC16-bound SIRT1^NAD+/Zn+2 due to conformational readjustments of 32 residues, as compared to SIRT1^NAD+/Zn+2. Due to increased flexibility, SIRT1-specific SER275, SER442 and ARG466 residues involved in NAD⁺ association facilitated in the formation of a transient tunnel (17.77 Å) that was further elongated to 19.25 Å upon SIRT1^NAD+/Zn+2 binding to H4^KAC16. A close conformation of SIRT1^NAD+/Zn+2 was achieved due to binding of H4^KAC16 that results in the movement of helical module towards Zn⁺² binding module together with Rossmann fold at NAD⁺ binding region. Furthermore, a 2-fold increase (4.31-8.82 Å) in the measured inter-atomic distance between imidazole nitrogen of conserved HIS363 and NAD⁺-specific 2'-hydroxyl group of ribose ring was evident in SIRT1^NAD+/Zn+2-H4^KAC16 complex. At 90 ns time scale, the distance between C6A of adenine ring and C2N of nicotinamide ring was more extended (19.32 Å) in SIRT1^NAD+/Zn+2-H4^KAC16 as compared to SIRT1^NAD+/Zn+2 (11.54 Å). These data suggest that H4^KAC16 binding to SIRT1 may coordinate an unusual conformational readjustment of nicotinamide ring at site-b and reposition of HIS363 to facilitate SIRT1-dependent deacetylase activity. Taken together, our findings will help in understanding the precise structural changes occurring in response to SIRT1 deacetylase activity of core histone.

S-nitrosylation affects TRAP1 structure and ATPase activity and modulates cell response to apoptotic stimuli

The mitochondrial chaperone TRAP1 has been involved in several mitochondrial functions, and modulation of its expression/activity has been suggested to play a role in the metabolic reprogramming distinctive of cancer cells. TRAP1 posttranslational modifications, i.e. phosphorylation, can modify its capability to bind to different client proteins and modulate its oncogenic activity. Recently, it has been also demonstrated that TRAP1 is S-nitrosylated at Cys501, a redox modification associated with its degradation via the proteasome. Here we report molecular dynamics simulations of TRAP1, together with analysis of long-range structural communication, providing a model according to which Cys501 S-nitrosylation induces conformational changes to distal sites in the structure of the protein. The modification is also predicted to alter open and closing motions for the chaperone function. By means of colorimetric assays and site directed mutagenesis aimed at generating C501S variant, we also experimentally confirmed that selective S-nitrosylation of Cys501 decreases ATPase activity of recombinant TRAP1. Coherently, C501S mutant was more active and conferred protection to cell death induced by staurosporine. Overall, our results provide the first in silico, in vitro and cellular evidence of the relevance of Cys501 S-nitrosylation in TRAP1 biology.

Impact of posttranslational modifications on atomistic structure of fibrinogen

Oxidative stress in humans is related to various pathophysiological processes, which can manifest in numerous diseases including cancer, cardiovascular diseases, and Alzheimer’s disease. On the atomistic level, oxidative stress causes posttranslational modifications, thus inducing structural and functional changes into the proteins structure. This study focuses on fibrinogen, a blood plasma protein that is frequently targeted by reagents causing posttranslational modifications in proteins. Fibrinogen was in vitro modified by three reagents, namely sodium hypochlorite, malondialdehyde, and 3-morpholinosydnonimine that mimic the oxidative stress in diseases. Newly induced posttranslational modifications were detected via mass spectrometry. Electron microscopy was used to visualize changes in the fibrin networks, which highlight the extent of disturbances in fibrinogen behavior after exposure to reagents. We used molecular dynamics simulations to observe the impact of selected posttranslational modifications on the fibrinogen structure at the atomistic level. In total, 154 posttranslational modifications were identified, 84 of them were in fibrinogen treated with hypochlorite, 51 resulted from a reaction of fibrinogen with malondialdehyde, and 19 were caused by 3-morpholinosydnonimine. Our data reveal that the stronger reagents induce more posttranslational modifications in the fibrinogen structure than the weaker ones, and they extensively alter the architecture of the fibrin network. Molecular dynamics simulations revealed that the effect of posttranslational modifications on fibrinogen secondary structure varies from negligible alternations to serious disruptions. Among the serious disruptions is the oxidation of γR375 resulting in the release of Ca²⁺ ion that is necessary for appropriate fibrin fiber formation. Folding of amino acids γE72–γN77 into a short α-helix is a result of oxidation of γP76 to glutamic acid. The study describes behaviour of fibrinogen coiled-coil connecter in the vicinity of plasmin and hementin cleavage sites.

Lipid Head Group Parameterization for GROMOS 54A8: A Consistent Approach with Protein Force Field Description

Membranes are a crucial component of both bacterial and mammalian cells, being involved in signaling, transport, and compartmentalization. This versatility requires a variety of lipid species to tailor the membrane’s behavior as needed, increasing the complexity of the system. Molecular dynamics simulations have been successfully applied to study model membranes and their interactions with proteins, elucidating some crucial mechanisms at the atomistic detail and thus complementing experimental techniques. An accurate description of the functional interplay of the diverse membrane components crucially depends on the selected parameters that define the adopted force field. A coherent parameterization for lipids and proteins is therefore needed. In this work, we propose and validate new lipid head group parameters for the GROMOS 54A8 force field, making use of recently published parametrizations for key chemical moieties present in lipids. We make use additionally of a new canonical set of partial charges for lipids, chosen to be consistent with the parameterization of soluble molecules such as proteins. We test the derived parameters on five phosphocholine model bilayers, composed of lipid patches four times larger than the ones used in previous studies, and run 500 ns long simulations of each system. Reproduction of experimental data like area per lipid and deuterium order parameters is good and comparable with previous parameterizations, as well as the description of liquid crystal to gel-phase transition. On the other hand, the orientational behavior of the head groups is more realistic for this new parameter set, and this can be crucial in the description of interactions with other polar molecules. For that reason, we tested the interaction of the antimicrobial peptide lactoferricin with two model membranes showing that the new parameters lead to a weaker peptide–membrane binding and give a more realistic outcome in comparing binding to antimicrobial versus mammal membranes.

Effect of in vivo post-translational modifications of the HMGB1 protein upon binding to platinated DNA: a molecular simulation study

Cisplatin is one of the most widely used anticancer drugs. Its efficiency is unfortunately severely hampered by resistance. The High Mobility Group Box (HMGB) proteins may sensitize tumor cells to cisplatin by specifically binding to platinated DNA (PtDNA) lesions. In vivo, the HMGB/PtDNA binding is regulated by multisite post-translational modifications (PTMs). The impact of PTMs on the HMGB/PtDNA complex at atomistic level is here investigated by enhanced sampling molecular simulations. The PTMs turn out to affect the structure of the complex, the mobility of several regions (including the platinated site), and the nature of the protein/PtDNA non-covalent interactions. Overall, the multisite PTMs increase significantly the apparent synchrony of all the contacts between the protein and PtDNA. Consequently, the hydrophobic anchoring of the side chain of F37 between the two cross-linked guanines at the platinated site-a key element of the complexes formation - is more stable than in the complex without PTM. These differences can account for the experimentally measured greater affinity for PtDNA of the protein isoforms with PTMs. The collective behavior of multisite PTMs, as revealed here by the synchrony of contacts, may have a general significance for the modulation of intermolecular recognitions occurring in vivo.

Transport of Reactive Oxygen and Nitrogen Species across Aquaporin: A Molecular Level Picture

Aquaporins (AQPs) are transmembrane proteins that conduct not only water molecules across the cell membrane but also other solutes, such as reactive oxygen and nitrogen species (RONS), produced (among others) by cold atmospheric plasma (CAP). These RONS may induce oxidative stress in the cell interior, which plays a role in cancer treatment. The underlying mechanisms of the transport of RONS across AQPs, however, still remain obscure. We apply molecular dynamics simulations to investigate the permeation of both hydrophilic (H₂O₂ and OH) and hydrophobic (NO₂ and NO) RONS through AQP1. Our simulations show that these RONS can all penetrate across the pores of AQP1. The permeation free energy barrier of OH and NO is lower than that of H₂O₂ and NO₂, indicating that these radicals may have easier access to the pore interior and interact with the amino acid residues of AQP1. We also study the effect of RONS-induced oxidation of both the phospholipids and AQP1 (i.e., sulfenylation of Cys₁₉₁) on the transport of the above-mentioned RONS across AQP1. Both lipid and protein oxidation seem to slightly increase the free energy barrier for H₂O₂ and NO₂ permeation, while for OH and NO, we do not observe a strong effect of oxidation. The simulation results help to gain insight in the underlying mechanisms of the noticeable rise of CAP-induced RONS in cancer cells, thereby improving our understanding on the role of AQPs in the selective anticancer capacity of CAP.

Oxidation destabilizes toxic amyloid beta peptide aggregation

The aggregation of insoluble amyloid beta (Aβ) peptides in the brain is known to trigger the onset of neurodegenerative diseases, such as Alzheimer’s disease. In spite of the massive number of investigations, the underlying mechanisms to destabilize the Aβ aggregates are still poorly understood. Some studies indicate the importance of oxidation to destabilize the Aβ aggregates. In particular, oxidation induced by cold atmospheric plasma (CAP) has demonstrated promising results in eliminating these toxic aggregates. In this paper, we investigate the effect of oxidation on the stability of an Aβ pentamer. By means of molecular dynamics simulations and umbrella sampling, we elucidate the conformational changes of Aβ pentamer in the presence of oxidized residues, and we estimate the dissociation free energy of the terminal peptide out of the pentamer form. The calculated dissociation free energy of the terminal peptide is also found to decrease with increasing oxidation. This indicates that Aβ pentamer aggregation becomes less favorable upon oxidation. Our study contributes to a better insight in one of the potential mechanisms for inhibition of toxic Aβ peptide aggregation, which is considered to be the main culprit to Alzheimer’s disease.

Effects of phosphorylation on the NLRP3 inflammasome

The pyrin domain containing Nod-like receptors (NLRPs) are a family of pattern recognition receptors known to regulate an array of immune signaling pathways. Emergent studies demonstrate the potential for regulatory control of inflammasome assembly by phosphorylation, notably NLRP3. Over a dozen phosphorylation sites have been identified for NLRP3 with many more suggested by phosphoproteomic studies of the NLRP family. Well characterized NLRP3 phosphorylation events include Ser198 by c-Jun terminal kinase (JNK), Ser295 by protein kinase D (PKD) and/or protein kinase A (PKA), and Tyr861 by an unknown kinase but is dephosphorylated by protein tyrosine phosphatase non-receptor 22 (PTPN22). Since the PKA- and PKD-dependent phosphorylation of NLRP3 at Ser295 is best characterized, we provide detailed review of this aspect of NLRP3 regulation. Phosphorylation of Ser295 can attenuate ATPase activity as compared to its dephosphorylated counterpart, and this event is likely unique to NLRP3. In silico modeling of NLRP3 is useful in predicting how Ser295 phosphorylation might impact upon the structural topology of the ATP-binding domain to influence catalytic activity. It is important to gain as complete understanding as possible of the complex phosphorylation-mediated mechanisms of regulation for NLRP3 in part because of its involvement in many pathological processes.

The Effect of FG-Nup Phosphorylation on NPC Selectivity: A One-Bead-Per-Amino-Acid Molecular Dynamics Study

Nuclear pore complexes (NPCs) are large protein complexes embedded in the nuclear envelope separating the cytoplasm from the nucleoplasm in eukaryotic cells. They function as selective gates for the transport of molecules in and out of the nucleus. The inner wall of the NPC is coated with intrinsically disordered proteins rich in phenylalanine-glycine repeats (FG-repeats), which are responsible for the intriguing selectivity of NPCs. The phosphorylation state of the FG-Nups is controlled by kinases and phosphatases. In the current study, we extended our one-bead-per-amino-acid (1BPA) model for intrinsically disordered proteins to account for phosphorylation. With this, we performed molecular dynamics simulations to probe the effect of phosphorylation on the Stokes radius of isolated FG-Nups, and on the structure and transport properties of the NPC. Our results indicate that phosphorylation causes a reduced attraction between the residues, leading to an extension of the FG-Nups and the formation of a significantly less dense FG-network inside the NPC. Furthermore, our simulations show that upon phosphorylation, the transport rate of inert molecules increases, while that of nuclear transport receptors decreases, which can be rationalized in terms of modified hydrophobic, electrostatic, and steric interactions. Altogether, our models provide a molecular framework to explain how extensive phosphorylation of FG-Nups decreases the selectivity of the NPC.