Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Computational approaches to investigate fluoride binding, selectivity and transport across the membrane

The use of molecular dynamics (MD) simulations to study biomolecular systems has proven reliable in elucidating atomic-level details of structure and function. In this chapter, MD simulations were used to uncover new insights into two phylogenetically unrelated bacterial fluoride (F-) exporters: the CLCF F-/H+ antiporter and the Fluc F- channel. The CLCF antiporter, a member of the broader CLC family, has previously revealed unique stoichiometry, anion-coordinating residues, and the absence of an internal glutamate crucial for proton import in the CLCs. Through MD simulations enhanced with umbrella sampling, we provide insights into the energetics and mechanism of the CLCF transport process, including its selectivity for F- over HF. In contrast, the Fluc F- channel presents a novel architecture as a dual topology dimer, featuring two pores for F- export and a central non-transported sodium ion. Using computational electrophysiology, we simulate the electrochemical gradient necessary for F- export in Fluc and reveal details about the coordination and hydration of both F- and the central sodium ion. The procedures described here delineate the specifics of these advanced techniques and can also be adapted to investigate other membrane protein systems.

Targeting reversible post-translational modifications with PROTACs: a focus on enzymes modifying protein lysine and arginine residues

PROTACs represent an emerging field in medicinal chemistry, which has already led to the development of compounds that reached clinical studies. Posttranslational modifications contribute to the complexity of proteomes, with 2846 disease-associated sites. PROTAC field is very advanced in targeting kinases, while its use for enzymes mediating posttranslational modifications of the basic amino acid residues, started to be developed recently. Therefore, we bring together this less popular class of PROTACs, targeting lysine acetyltransferases/deacetylases, lysine and arginine methyltransferases, ADP-ribosyltransferases, E3 ligases, and ubiquitin-specific proteases. We put special emphasis on structural aspects of PROTAC elements to facilitate the lengthy experimental endeavours directed towards developing PROTACs. We will cover the period from the inception of the field, 2017, to April 2023.

Robust AMBER Force Field Parameters for Glutathionylated Cysteines

S-glutathionylation is an oxidative post-translational modification, which is involved in the regulation of many cell signaling pathways. Increasing amounts of studies show that it is crucial in cell homeostasis and deregulated in several pathologies. However, the effect of S-glutathionylation on proteins' structure and activity is poorly understood, and a drastic lack of structural information at the atomic scale remains. Studies based on the use of molecular dynamics simulations, which can provide important information about modification-induced modulation of proteins' structure and function, are also sparse, and there is no benchmarked force field parameters for this modified cysteine. In this contribution, we provide robust AMBER parameters for S-glutathionylation, which we tested extensively against experimental data through a total of 33 μs molecular dynamics simulations. We show that our parameter set efficiently describes the global and local structural properties of S-glutathionylated proteins. These data provide the community with an important tool to foster new investigations into the effect of S-glutathionylation on protein dynamics and function, in a common effort to unravel the structural mechanisms underlying its critical role in cellular processes.

Effect of histidine covalent modification on strigolactone receptor activation and selectivity

The parasitic weed Striga has led to billions of dollars' worth of agricultural productivity loss worldwide. Striga detects host plants using compounds of the strigolactone class of phytohormones. Early steps in the strigolactone signaling pathway involve substrate binding and hydrolysis followed by a conformational change to an "active" or "closed" state, after which it associates with a MAX2-family downstream signaling partner. The structures of the inactive and active states of strigolactone receptors are known through X-ray crystallography, and the transition pathway from the inactive to active state in apo receptors has previously been characterized using molecular dynamics simulations. However, it also has been suggested that a covalent butenolide modification of the receptor on the catalytic histidine through substrate hydrolysis promotes formation of the active state. Using molecular dynamics simulations, we show that the presence of the covalent butenolide enhances activation in both AtD14, a receptor found in Arabidopsis, and ShHTL7, a receptor found in Striga, but the enhancement is ∼50 times greater in ShHTL7. We also show that several conserved interactions with the covalent butenolide modification promote transition to the active state in both AtD14 (non-parasite) and ShHTL7 (parasite). Finally, we demonstrate that the enhanced activation of ShHTL7 likely results from disruption of ShHTL7-specific histidine interactions that inhibited activation in the apo case. These results provide a possible explanation for difference in strigolactone sensitivity seen between different strigolactone-sensitive proteins and can be used to aid the design of selective modulators to control Striga parasites.