Modeling Protein S-Aromatic Motifs Reveals Their Structural and Redox Flexibility

Proximal Methionine Amino Acid Residue Affects the Properties of Redox-Active Tryptophan in an Artificial Model Protein

Redox-active amino acid residues are at the heart of biological electron-transfer reactions. They play important roles in natural protein functions and are implicated in disease states (e.g., oxidative-stress-associated disorders). Tryptophan (Trp) is one such redox-active amino acid residue, and it has long been known to serve a functional role in proteins. Broadly speaking, there is still much to learn about the local features that make some Trp redox active and others inactive. Herein, we describe a new protein model system where we investigate how a methionine (Met) residue proximal to a redox-active Trp affects its reactivity and spectroscopy. We use an artificial variant of azurin from Pseudomonas aeruginosa to produce these models. We employ a series of UV-visible spectroscopy, electrochemistry, electron paramagnetic resonance, and density functional theory experiments to demonstrate the effect that placing Met near Trp radicals has in the context of redox proteins. The introduction of Met proximal to Trp lowers its reduction potential by ca. 30 mV and causes clear shifts in the optical spectra of the corresponding radicals. While the effect may be small, it is significant enough to be a way for natural systems to tune Trp reactivity.

The Potential of a Protein Model Synthesized Absent of Methionine

Methionine is an amino acid long thought to be essential, but only in the case of protein synthesis initiation. In more recent years, methionine has been found to play an important role in antioxidant defense, stability, and modulation of cell and protein activity. Though these findings have expanded the previously held sentiment of methionine having a singular purpose within cells and proteins, the essential nature of methionine can still be challenged. Many of the features that give methionine its newfound functions are shared by the other sulfur-containing amino acid: cysteine. While the antioxidant, stabilizing, and cell/protein modulatory functions of cysteine have already been well established, recent findings have shown a similar hydrophobicity to methionine which suggests cysteine may be able to replace methionine in all functions outside of protein synthesis initiation with little effect on cell and protein function. Furthermore, a number of novel mechanisms for alternative initiation of protein synthesis have been identified that suggest a potential to bypass the traditional methionine-dependent initiation during times of stress. In this review, these findings are discussed with a number of examples that demonstrate a potential model for synthesizing a protein in the absence of methionine.

On the roles of methionine and the importance of its microenvironments in redox metalloproteins

The amino acid residue methionine (Met) is commonly thought of as a ligand in redox metalloproteins, for example in cytochromes c and in blue copper proteins. However, the roles of Met can go beyond a simple ligand. The thioether functional group of Met allows it to be considered as a hydrophobic residue as well as one that is capable of weak dipolar interactions. In addition, the lone pairs on sulphur allow Met to interact with other groups, inluding the aforementioned metal ions. Because of its properties, Met can play diverse roles in metal coordination, fine tuning of redox reactions, or supporting protein structures. These roles are strongly influenced by the nature of the surrounding medium. Herein, we describe several common interactions between Met and surrounding aromatic amino acids and how they affect the physical properties of both copper and iron metalloproteins. While the importance of interactions between Met and other groups is established in biological systems, less is known about their roles in redox metalloproteins and our view is that this is an area that is ready for greater attention.

Comparative Analysis of Sulfonium-π, Ammonium-π, and Sulfur-π Interactions and Relevance to SAM-Dependent Methyltransferases

We report the measurement and analysis of sulfonium-π, thioether-π, and ammonium-π interactions in a β-hairpin peptide model system, coupled with computational investigation and PDB analysis. These studies indicated that the sulfonium-π interaction is the strongest and that polarizability contributes to the stronger interaction with sulfonium relative to ammonium. Computational studies demonstrate that differences in solvation of the trimethylsulfonium versus the trimethylammonium group also contribute to the stronger sulfonium-π interaction. In comparing sulfonium-π versus sulfur-π interactions in proteins, analysis of SAM- and SAH-bound enzymes in the PDB suggests that aromatic residues are enriched in close proximity to the sulfur of both SAM and SAH, but the populations of aromatic interactions of the two cofactors are not significantly different, with the exception of the Me-π interactions in SAM, which are the most prevalent interaction in SAM but are not possible for SAH. This suggests that the weaker interaction energies due to loss of the cation-π interaction in going from SAM to SAH may contribute to turnover of the cofactor.

Corrections in the CHARMM36 Parametrization of Chloride Interactions with Proteins, Lipids, and Alkali Cations, and Extension to Other Halide Anions

The nonpolarizable CHARMM force field is one of the most widely used energy functions for all-atom biomolecular simulations. Chloride is the only halide ion included in the latest version, CHARMM36m, and is used widely in simulation studies, often as an electrolyte ion but also as the biological substrate of transport proteins and enzymes. Here, we find that existing parameters systematically underestimate the interaction of Cl^- with proteins and lipids. Accordingly, when examined in solution, little to no Cl^-association can be observed with most components of the protein, including backbone, polar side chains and aromatic rings. The strength of the interaction with cationic side chains and with alkali ions is also incongruent with experimental measurements, specifically osmotic coefficients of concentrated solutions. Consistent with these findings, a 4-μs trajectory of the Cl^--specific transport protein CLC-ec1 shows irreversible Cl^- dissociation from the so-called S_cen binding site, even in a 150 mM NaCl buffer. To correct for these deficiencies, we formulate a series of pair-specific Lennard-Jones parameters that override those resulting from the conventional Lorentz-Berthelot combination rules. These parameters, referred to as NBFIX, are systematically calibrated against available experimental data as well as ab initio geometry optimizations and energy evaluations, for a wide set of binary and ternary Cl^- complexes with protein and lipid analogs and alkali cations. Analogously, we also formulate parameter sets for the other three biological halide ions, namely, fluoride, bromide, and iodide. The resulting parameters are used to calculate the potential of mean force defining the interaction of each anion and each of the protein and lipid analogues in bulk water, revealing association free energies in the range of -0.3 to -3.3 kcal/mol, with the F^- complexes being the least stable. The NBFIX corrections also preserve the Cl^- occupancy of CLC-ec1 in a second 4-μs trajectory. We posit that these optimized molecular-mechanics models provide a more realistic foundation for all-atom simulation studies of processes entailing changes in hydration, recognition, or transport of halide anions.

The Role of Methionine Residues in the Regulation of Liquid-Liquid Phase Separation

Membraneless organelles are non-stoichiometric supramolecular structures in the micron scale. These structures can be quickly assembled/disassembled in a regulated fashion in response to specific stimuli. Membraneless organelles contribute to the spatiotemporal compartmentalization of the cell, and they are involved in diverse cellular processes often, but not exclusively, related to RNA metabolism. Liquid-liquid phase separation, a reversible event involving demixing into two distinct liquid phases, provides a physical framework to gain insights concerning the molecular forces underlying the process and how they can be tuned according to the cellular needs. Proteins able to undergo phase separation usually present a modular architecture, which favors a multivalency-driven demixing. We discuss the role of low complexity regions in establishing networks of intra- and intermolecular interactions that collectively control the phase regime. Post-translational modifications of the residues present in these domains provide a convenient strategy to reshape the residue-residue interaction networks that determine the dynamics of phase separation. Focus will be placed on those proteins with low complexity domains exhibiting a biased composition towards the amino acid methionine and the prominent role that reversible methionine sulfoxidation plays in the assembly/disassembly of biomolecular condensates.

Molecular basis for lipid recognition by the prostaglandin D2 receptor CRTH2

Prostaglandin D₂ (PGD₂) signals through the G protein-coupled receptor (GPCR) CRTH2 to mediate various inflammatory responses. CRTH2 is the only member of the prostanoid receptor family that is phylogenetically distant from others, implying a nonconserved mechanism of lipid action on CRTH2. Here, we report a crystal structure of human CRTH2 bound to a PGD₂ derivative, 15R-methyl-PGD₂ (15mPGD₂), by serial femtosecond crystallography. The structure revealed a "polar group in"-binding mode of 15mPGD₂ contrasting the "polar group out"-binding mode of PGE₂ in its receptor EP3. Structural comparison analysis suggested that these two lipid-binding modes, associated with distinct charge distributions of ligand-binding pockets, may apply to other lipid GPCRs. Molecular dynamics simulations together with mutagenesis studies also identified charged residues at the ligand entry port that function to capture lipid ligands of CRTH2 from the lipid bilayer. Together, our studies suggest critical roles of charge environment in lipid recognition by GPCRs.

An "up" oriented methionine-aromatic structural motif in SUMO is critical for its stability and activity

Protein structural bioinformatic analyses suggest preferential associations between methionine and aromatic amino acid residues in proteins. Ab initio energy calculations highlight a conformation-dependent stabilizing interaction between the interacting sulfur-aromatic molecular pair. However, the relevance of buried methionine-aromatic motifs to protein folding and function is relatively unexplored. The Small Ubiquitin-Like Modifier (SUMO) is a β-grasp fold protein and a common posttranslational modifier that affects diverse cellular processes, including transcriptional regulation, chromatin remodeling, metabolic regulation, mitosis, and meiosis. SUMO is a member of the Ubiquitin-Like (UBL) protein family. Herein, we report that a highly conserved and buried methionine-phenylalanine motif is a unique signature of SUMO proteins but absent in other homologous UBL proteins. We also detect that a specific "up" conformation between the methionine-phenylalanine pair of interacting residues in SUMO is critical to its β-grasp fold. The noncovalent interactions of SUMO with its ligands are dependent on the methionine-phenylalanine pair. MD simulations, NMR, and biophysical and biochemical studies suggest that perturbation of the methionine-aromatic motif disrupts native contacts, modulates noncovalent interactions, and attenuates SUMOylation activity. Our results highlight the importance of conserved orientations of Met-aromatic structural motifs inside a protein core for its structure and function.

Modeling Shows that Rotation about the Peroxide O-O Bond Assists Protein and Lipid Functional Groups in Discriminating between H2O2 and H2O

Long associated with cell death, hydrogen peroxide (H₂O₂) is now known to perform many physiological roles. Unraveling its biological mechanisms of action requires atomic-level knowledge of its association with proteins and lipids, which we address here. High-level [MP2(full)/6-311++G(3df,3pd)] ab initio calculations reveal skew rotamers as the lowest-energy states of isolated H₂O₂ (ϕ_HOOH ∼ 112°) with minimum and maximum electrostatic potentials (kcal/mol) of -24.8 (V_s,min) and 36.5 (V_s,max), respectively. Transition-state, nonpolar trans rotamers (ϕ_HOOH ∼ 180°) at 1.2 kcal/mol higher in energy are poorer H-bond acceptors (V_s,min = -16.6) than the skew rotamers, while highly polar cis rotamers (ϕ_HOOH ∼ 0°) at 7.8 kcal/mol are much better H-bond donors (V_s,max = 52.7). Modeling H₂O₂ association with neutral and charged analogs of protein residues and lipid groups (e.g., ester, phosphate, choline) reveals that skew rotamers (ϕ_HOOH = 84-122°) are favored in the neutral and cationic complexes, which display gas-phase interaction energies (E^CP, kcal/mol) of -1.5 to -18. The neutral and cationic complexes of H₂O exhibit a similar range of stabilities (E^CP ∼ -1 to -18). However, considerably higher energies (E^CP ∼ -14 to -36) are found for the H₂O₂ complexes of the anionic ligands, which are stabilized by charge-assisted H-bond donation from cis and distorted cis rotamers (ϕ_HOOH = 0-60°). H₂O is a much poorer H-bond donor (V_s,max = 33.4) than cis-H₂O₂, so its anionic complexes are significantly weaker (E^CP ∼ -11 to -20). Thus, by dictating the rotamer preference of H₂O₂, functional groups in biomolecules can discriminate between H₂O₂ and H₂O. Finally, exploiting the present ab initio data, we calibrated and validated our published molecular mechanics model for H₂O₂ (Orabi, E. A.; English, A. M. J. Chem. Theory Comput. 2018, 14, 2808-2821) to provide an important tool for simulating H₂O₂ in biology.

New Molecular-Mechanics Model for Simulations of Hydrogen Fluoride in Chemistry and Biology

Hydrogen fluoride (HF) is the most polar diatomic molecule and one of the simplest molecules capable of hydrogen-bonding. HF deviates from ideality both in the gas phase and in solution and is thus of great interest from a fundamental standpoint. Pure and aqueous HF solutions are broadly used in chemical and industrial processes, despite their high toxicity. HF is a stable species also in some biological conditions, because it does not readily dissociate in water unlike other hydrogen halides; yet, little is known about how HF interacts with biomolecules. Here, we set out to develop a molecular-mechanics model to enable computer simulations of HF in chemical and biological applications. This model is based on a comprehensive high-level ab initio quantum chemical investigation of the structure and energetics of the HF monomer and dimer; (HF)_n clusters, for n = 3-7; various clusters of HF and H₂O; and complexes of HF with analogs of all 20 amino acids and of several commonly occurring lipids, both neutral and ionized. This systematic analysis explains the unique properties of this molecule: for example, that interacting HF molecules favor nonlinear geometries despite being diatomic and that HF is a strong H-bond donor but a poor acceptor. The ab initio data also enables us to calibrate a three-site molecular-mechanics model, with which we investigate the structure and thermodynamic properties of gaseous, liquid, and supercritical HF in a wide range of temperatures and pressures; the solvation structure of HF in water and of H₂O in liquid HF; and the free diffusion of HF across a lipid bilayer, a key process underlying the high cytotoxicity of HF. Despite its inherent simplifications, the model presented significantly improves upon previous efforts to capture the properties of pure and aqueous HF fluids by molecular-mechanics methods and to our knowledge constitutes the first parameter set calibrated for biomolecular simulations.

Hole Hopping through Cytochrome P450

High-potential iron-oxo species are intermediates in the catalytic cycles of oxygenase enzymes. They can cause heme degradation and irreversible oxidation of nearby amino acids. We have proposed that there are protective mechanisms in which hole hopping from oxidized hemes through tryptophan/tyrosine chains generates a surface-exposed amino-acid oxidant that could be rapidly disarmed by reaction with cellular reductants. In investigations of cytochrome P450_BM3, we identified Trp96 as a critical residue that could play such a protective role. This Trp is cation-π paired with Arg398 in 81% of mammalian P450s. Here we report on the effect of the Trp/Arg cation-π interaction on Trp96 formal potentials as well as on electronic coupling strengths between Trp96 and the heme both for wild type cytochrome P450 and selected mutants. Mutation of Arg398 to His, which decreases the Trp96 formal potential, increases Trp-heme electronic coupling; however, surprisingly, the rate of phototriggered electron transfer from a Ru-sensitizer (through Trp96) to the P450_BM3 heme was unaffected by the Arg398His mutation. We conclude that Trp96 has moved away from Arg398, suggesting that the protective mechanism for P450s with this Trp-Arg pair is conformationally gated.

77Se NMR Probes the Protein Environment of Selenomethionine

Sulfur is critical for the correct structure and proper function of proteins. Yet, lacking a sensitive enough isotope, nuclear magnetic resonance (NMR) experiments are unable to deliver for sulfur in proteins the usual wealth of chemical, dynamic, and structural information. This limitation can be circumvented by substituting sulfur with selenium, which has similar physicochemical properties and minimal impact on protein structures but possesses an NMR compatible isotope (⁷⁷Se). Here we exploit the sensitivity of ⁷⁷Se NMR to the nucleus' chemical milieu and use selenomethionine as a probe for its proteinaceous environment. However, such selenium NMR spectra of proteins currently resist a reliable interpretation because systematic connections between variations of system variables and changes in ⁷⁷Se NMR parameters are still lacking. To start narrowing this knowledge gap, we report here on a biological ⁷⁷Se magnetic resonance data bank based on a systematically designed library of GB1 variants in which a single selenomethionine was introduced at different locations within the protein. We recorded the resulting isotropic ⁷⁷Se chemical shifts and relaxation times for six GB1 variants by solution-state ⁷⁷Se NMR. For four of the GB1 variants we were also able to determine the chemical shift anisotropy tensor of SeM by solid-state ⁷⁷Se NMR. To enable interpretation of the NMR data, the structures of five of the GB1 variants were solved by X-ray crystallography to a resolution of 1.2 Å, allowing us to unambiguously determine the conformation of the selenomethionine. Finally, we combine our solution- and solid-state NMR data with the structural information to arrive at general insights regarding the execution and interpretation of ⁷⁷Se NMR experiments that exploit selenomethionine to probe proteins.

Methionine in proteins: The Cinderella of the proteinogenic amino acids

Methionine in proteins, apart from its role in the initiation of translation, is assumed to play a simple structural role in the hydrophobic core, in a similar way to other hydrophobic amino acids such as leucine, isoleucine, and valine. However, research from a number of laboratories supports the concept that methionine serves as an important cellular antioxidant, stabilizes the structure of proteins, participates in the sequence-independent recognition of protein surfaces, and can act as a regulatory switch through reversible oxidation and reduction. Despite all these evidences, the role of methionine in protein structure and function is largely overlooked by most biochemists. Thus, the main aim of the current article is not so much to carry out an exhaustive review of the many and diverse processes in which methionine residues are involved, but to review some illustrative examples that may help the nonspecialized reader to form a richer and more precise insight regarding the role-played by methionine residues in such processes.

Computational insight into hydrogen persulfide and a new additive model for chemical and biological simulations

S-Sulfhydration of cysteine to the Cys-SSH persulfide is an oxidative post-translational modification that plays an important regulatory role in many physiological systems. Though hydrogen persulfide (H2S2) has recently been established as a signaling and cellular sulfhydration reagent, the chemistry and chemical biology of persulfides remain poorly explored. We first report an extensive high-level ab initio quantum chemical investigation of (H2S2)n, (H2S2)m·H2O, and (H2O)m·H2S2 clusters (n = 1-3 and m = 1, 2) and of H2S2 complexes with 19 compounds that model the side chains of naturally-occurring amino acids. The high polarizability of S necessitates the use of large, very diffuse, basis sets for proper description of H2S2 and its complexes. H2S2 possesses a skewed equilibrium geometry, with nonpolar trans and more polar cis conformers 6 and 8 kcal mol-1 higher in energy, respectively; the skewed conformation is preserved in all neutral and cationic complexes while a cis geometry prevails in some anionic complexes. H2S2 is found to be a better H-bond donor and a poorer acceptor than H2S, and that in complexes with H2O, alcohols and amines, H2S2 is a better H-bond donor. Radical delocalization on both S atoms stabilizes the perthiyl (HSS˙) over the thiyl (HS˙) radical and results in a ∼20 kcal mol-1 lower S-H homolytic bond dissociation in H2S2, making it a potential antioxidant. A simple additive model is optimized for H2S2 and used together with the TIP3P model and the CHARMM36 all-atom force field (FF) to investigate the structure and thermodynamic properties of liquid H2S2 and the solubility of H2S2 in water, and to model H2S2-protein interactions (for which new FF parameters are further developed). Very weak H-bonding characterizes liquid H2S2 and it is found immiscible in liquid water with a trend in H-bonding strengths between H2S2 and H2O in the order O-HO ≫ S-HO > O-HS. This work does not only provide a thorough description of the structure and energetics of H2S2 and its various complexes, but also yields a reliable FF for investigating H2S2 in chemistry and biology.

Expanding the range of binding energies and oxidizability of biologically relevant S-aromatic interactions: imidazolium and phenolate binding to sulfoxide and sulfone

Oxidation and protonation/deprotonation strongly impact intermolecular noncovalent interactions. For example, S-aromatic interactions are stabilized up to three-fold in the gas phase on oxidation of the sulfur ligand or protonation/deprotonation of the aromatic. To probe if such stabilizing effects are additive and to model interactions of oxidized methionine (MetOn) with protonated histidine and deprotonated tyrosine residues in proteins, we examined Me2SOn (n = 1, 2) binding to imidazolium, phenolate and their 4-methylated forms. Ab initio MP2(full)/6-311++G(d,p) gas-phase calculations reveal that the Me2SOn-imidazolium complexes adopt edge-on geometry with σ-type (N/C-HarO) H-bonding and interaction energies of -17.2 to -31.1 kcal mol-1. The less stable (-13.8 to -21.0 kcal mol-1) Me2SOn-phenolates possess en-face geometry stabilized by π-type (C-Hπar) H-bonding. Comparing these energies with those reported for the Me2S-neutral aromatics affirms the additive effects of ligand protonation/deprotonation and oxidation on gas-phase stability. However, this is not the case in water although the aqueous complexes retain their preferred gas-phase σ- and π-type H-bonded structures. Binding free energies (kcal mol-1) calculated from molecular dynamics simulations in bulk water (preceded by CHARMM36 force field calibration where necessary) reveal that Me2SO-imidazolium (-4.4) is more stable than Me2SO-phenolate (-2.4) but Me2SO2-imidazolium (-0.6) is less stable than Me2SO2-phenolate (-3.8). Vertical ionization potentials (IPV) calculated for the gas-phase complexes indicate that the Me2SOn-phenolates, but not the Me2SOn-imidazoles, are oxidizable under biological conditions. Charge transfer from the phenolate increases its IPV by ∼20%, decreasing its susceptibility to oxidation. Overall, this work provides fundamental data to predict the behaviour of protein-based MetOn-aromatic-ion interactions.

LC-MS/MS Proteoform Profiling Exposes Cytochrome c Peroxidase Self-Oxidation in Mitochondria and Functionally Important Hole Hopping from Its Heme

LC-MS/MS profiling reveals that the proteoforms of cytochrome c peroxidase (Ccp1) isolated from respiring yeast mitochondria are oxidized at numerous Met, Trp, and Tyr residues. In vitro oxidation of recombinant Ccp1 by H₂O₂ in the absence of its reducing substrate, ferrocytochrome c, gives rise to similar proteoforms, indicating uncoupling of Ccp1 oxidation and reduction in mitochondria. The oxidative modifications found in the Ccp1 proteoforms are consistent with radical transfer (hole hopping) from the heme along several chains of redox-active residues (Trp, Met, Tyr). These modifications delineate likely hole-hopping pathways to novel substrate-binding sites. Moreover, a decrease in recombinant Ccp1 oxidation by H₂O₂ in vitro in the presence of glutathione supports a protective role for hole hopping to this antioxidant. Isolation and characterization of extramitochondrial Ccp1 proteoforms reveals that hole hopping from the heme in these proteoforms results in selective oxidation of the proximal heme ligand (H175) and heme labilization. Previously, we demonstrated that this labilized heme is recruited for catalase maturation (Kathiresan, M.; Martins, D.; English, A. M. Respiration triggers heme transfer from cytochrome c peroxidase to catalase in yeast mitochondria. Proc. Natl. Acad. Sci. U. S. A. 2014, 111, 17468-17473; DOI: 10.1073/pnas.1409692111 ). Following heme release, apoCcp1 exits mitochondria, yielding the extramitochondrial proteoforms that we characterize here. The targeting of Ccp1 for selective H175 oxidation may be linked to the phosphorylation status of Y153 close to the heme since pY153 is abundant in certain proteoforms. In sum, when insufficient electrons from ferrocytochrome c are available to Ccp1 in mitochondria, hole hopping from its heme expands its physiological functions. Specifically, we observe an unprecedented hole-hopping sequence for heme labilization and identify hole-hopping pathways from the heme to novel substrates and to glutathione at Ccp1's surface. Furthermore, our results underscore the power of proteoform profiling by LC-MS/MS in exploring the cellular roles of oxidoreductases.

Predicting structural and energetic changes in Met–aromatic motifs on methionine oxidation to the sulfoxide and sulfone

Methionine oxidation increases its affinity for aromatics in the gas phase but lowers it for most complexes in water.