Solution structure of the N-terminal extension domain of a Schistosoma japonicum asparaginyl-tRNA synthetase

Several secreted proteins from helminths (parasitic worms) have been shown to have immunomodulatory activities. Asparaginyl-tRNA synthetases are abundantly secreted in the filarial nematode Brugia malayi (BmAsnRS) and the parasitic flatworm Schistosoma japonicum (SjAsnRS), indicating a possible immune function. The suggestion is supported by BmAsnRS alleviating disease symptoms in a T-cell transfer mouse model of colitis. This immunomodulatory function is potentially related to an N-terminal extension domain present in eukaryotic AsnRS proteins but few structure/function studies have been done on this domain. Here we have determined the three-dimensional solution structure of the N-terminal extension domain of SjAsnRS. A protein containing the 114 N-terminal amino acids of SjAsnRS was recombinantly expressed with isotopic labelling to allow structure determination using 3D NMR spectroscopy, and analysis of dynamics using NMR relaxation experiments. Structural comparisons of the N-terminal extension domain of SjAsnRS with filarial and human homologues highlight a high degree of variability in the β-hairpin region of these eukaryotic N-AsnRS proteins, but similarities in the disorder of the C-terminal regions. Limitations in PrDOS-based intrinsically disordered region (IDR) model predictions were also evident in this comparison. Empirical structural data such as that presented in our study for N-SjAsnRS will enhance the prediction of sequence-homology based structure modelling and prediction of IDRs in the future.Communicated by Ramaswamy H. Sarma.

Control of charge transport in electronically active systems towards integrated biomolecular circuits (IbC)

The miniaturization of traditional silicon-based electronics will soon reach its limitation as quantum tunneling and heat become serious problems at the several-nanometer scale. Crafting integrated circuits via self-assembly of electronically active molecules using a "bottom-up" paradigm provides a potential solution to these technological challenges. In particular, integrated biomolecular circuits (IbC) offer promising advantages to achieve this goal, as nature offers countless examples of functionalities entailed by self-assembly and examples of controlling charge transport at the molecular level within the self-assembled structures. To this end, the review summarizes the progress in understanding how charge transport is regulated in biosystems and the key redox-active amino acids that enable the charge transport. In addition, charge transport mechanisms at different length scales are also reviewed, offering key insights for controlling charge transport in IbC in the future.

Infrared Spectroscopy of (Benzene-H2S-Xn)+, X = H2O (n = 1 and 2) and CH3OH (n = 1), Radical Cation Clusters: Microsolvation Effects on the S-π Hemibond

An unconventional covalent bond in which three electrons are shared by two centers is called hemibond. Hemibond formation frequently competes with proton transfer (or ionic hydrogen bond formation), but there have been a few experimental reports on such competition. In the present study, we focus on the (benzene-H₂S)⁺ radical cation cluster, which is a model system of the S-π hemibond. The stability of the S-π hemibond to the microsolvation by water and methanol is explored with infrared spectroscopy of (benzene-H₂S-X_n)⁺, X = H₂O (n = 1 and 2) and CH₃OH (n = 1), clusters. We also perform energy-optimization and vibrational simulations of (benzene-H₂S-X_n)⁺. By comparison among the observed and simulated spectra, we determine the intermolecular binding motifs in (benzene-H₂S-X_n)⁺. While the S-π hemibonded isomer is exclusively populated in (benzene-H₂S-H₂O)⁺, both the hemibonded and proton-transferred isomers coexist in [benzene-H₂S-(H₂O)₂]⁺ and (benzene-H₂S-CH₃OH)⁺. Breaking of the S-π hemibond by the microsolvation is observed, and its solvent and cluster size dependence is interpreted by the proton affinity and the coordination property of the solvent moiety.

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.

Exploring amino acid functions in a deep mutational landscape

Amino acids fulfil a diverse range of roles in proteins, each utilising its chemical properties in different ways in different contexts to create required functions. For example, cysteines form disulphide or hydrogen bonds in different circumstances and charged amino acids do not always make use of their charge. The repertoire of amino acid functions and the frequency at which they occur in proteins remains understudied. Measuring large numbers of mutational consequences, which can elucidate the role an amino acid plays, was prohibitively time-consuming until recent developments in deep mutational scanning. In this study, we gathered data from 28 deep mutational scanning studies, covering 6,291 positions in 30 proteins, and used the consequences of mutation at each position to define a mutational landscape. We demonstrated rich relationships between this landscape and biophysical or evolutionary properties. Finally, we identified 100 functional amino acid subtypes with a data-driven clustering analysis and studied their features, including their frequencies and chemical properties such as tolerating polarity, hydrophobicity or being intolerant of charge or specific amino acids. The mutational landscape and amino acid subtypes provide a foundational catalogue of amino acid functional diversity, which will be refined as the number of studied protein positions increases.

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.

Computational analysis of vibrational frequencies and rovibrational spectroscopic constants of hydrogen sulfide dimer using MP2 and CCSD(T)

Previous studies have shown that the weakly bonded H₂S dimer demands high level quantum chemical calculations to reproduce experimental values. We investigated the hydrogen bonding of H₂S dimer using MP2 and CCSD(T) levels of theory in combination with aug-cc-pV(D,T,Q)Z basis sets. More precisely, the binding energies, potential energy curves, rovibrational spectroscopic constants, decomposition lifetime, and normal vibrational frequencies were calculated. In addition, we introduced the local mode analysis of Konkoli-Cremer to quantify the hydrogen bonding in the H₂S dimer as well as providing for the first time the comprehensive decomposition of normal vibrational modes into local modes contributions, and a decomposition lifetime based on rate constant. The local mode force constant of the H₂S dimer hydrogen bond is smaller than that of the water dimer, in accordance with the weaker hydrogen bonding in the H₂S dimer.

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.

Sulfur-centered hemi-bond radicals as active intermediates in S-DNA phosphorothioate oxidation

Phosphorothioate (PS) modifications naturally appear in bacteria and archaea genome and are widely used as antisense strategy in gene therapy. But the chemical effects of PS introduction as a redox active site into DNA (S-DNA) is still poorly understood. Herein, we perform time-resolved spectroscopy to examine the underlying mechanisms and dynamics of the PS oxidation by potent radicals in free model, in dinucleotide, and in S-oligomer. The crucial sulphur-centered hemi-bonded intermediates -P–S∴S–P- were observed and found to play critical roles leading to the stable adducts of -P–S–S–P-, which are backbone DNA lesion products. Moreover, the oxidation of the PS moiety in dinucleotides d[G_PSG], d[A_PSA], d[G_PSA], d[A_PSG] and in S-oligomers was monitored in real-time, showing that PS oxidation can compete with adenine but not with guanine. Significantly, hole transfer process from A^+• to PS and concomitant -P–S∴S–P- formation was observed, demonstrating the base-to-backbone hole transfer unique to S-DNA, which is different from the normally adopted backbone-to-base hole transfer in native DNA. These findings reveal the distinct backbone lesion pathway brought by the PS modification and also imply an alternative -P–S∴S–P-/-P–S–S–P- pathway accounting for the interesting protective role of PS as an oxidation sacrifice in bacterial genome.

A crosslinked, low pH-stable, mixed-mode cation-exchange like stationary phase made using the thiol-yne click reaction

Mixed-mode cation-exchange stationary phases are useful for the separation of mixtures containing hydrophobic, acidic, and basic molecules. To ensure that weak organic bases are protonated and carboxylic acids are neutral low pH mobile phases are required. Mixed-mode stationary phases that are stable at pH < 3 are needed. We synthesized a crosslinked structure along the surface of thiol functionalized silica gel particles using the thiol-yne click reaction. The alkyne, 1,7-octadiyne, was added to the 3-mercaptopropyl silica gel, then crosslinked using 1,6-hexanedithiol. Elemental analysis showed low octadiyne ligand surface coverage, but, stoichiometrically, three sulfurs were added to each octadiyne ligand during the crosslinking step, indicating that crosslinking occurred. The effect of the crosslinking on the stability was tested with a 50:50 (v/v) pH 0.50 5% TFA aqueous:acetonitrile mobile phase at 70 °C for six days, over 35,000 column volumes. The stationary phase showed good stability with the retention of triphenylene decreasing only 20% during that time. The Tanaka test showed that the phase has a methylene selectivity of 1.20 ± 0.04, a high shape selectivity of 2.71 ± 0.03, and a 3.98 ± 0.05 cation-exchange factor at pH 2.70. The phase has a selectivity factor for nitrobenzene and benzene of 1.41 ± 0.01, indicating the electron donating charge transfer characteristic of the phase. The mixed-mode characteristics of the phase were investigated using a mixture of the monoamine neurotransmitters norepinephrine, dopamine, and serotonin. Baseline resolution of the monoamines could be obtained using a simple 20 mM potassium phosphate (pH 2.70)/methanol mobile phase. Altering both the methanol content and the potassium ion concentration altered the retention of the monoamines indicating mixed-mode cation exchange characteristic of the crosslinked stationary phase.

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.

The S∴π hemibond and its competition with the S∴S hemibond in the simplest model system: infrared spectroscopy of the [benzene-(H2S) n ]+ (n = 1-4) radical cation clusters

The S∴π hemibond (two-center three-electron, 2c-3e, bond) is an attractive interaction between a sulfur atom and π electrons. The S∴π hemibond is of essential importance in understanding chemistry of sulfur radical cations, and its roles in biochemistry have recently attracted much interest. In the present study, we observe the S∴π hemibond in the simplest model system in the gas phase. Infrared spectroscopy is applied to the [benzene-(H2S) n ]+ (n = 1-4) radical cation clusters. In n = 1, the CH stretch and SH stretch bands of the benzene and H2S moieties, respectively, are clearly different from those of the neutral molecules but similar to those of the ionic species. These vibrational features show that the positive charge is delocalized over the cluster due to the S∴π hemibond formation. In n = 2-4, the S∴S hemibond and S-π-S multicenter hemibond (three-center five-electron, 3c-5e, bond) can compete with the S∴π hemibond. The observed vibrational features clearly indicate that the S∴S hemibond formation is superior to the S∴π hemibond and S-π-S multicenter hemibond. Calculations of several dispersion-corrected density functionals are compared with the observations. While all the tested functionals qualitatively catch the feature of the S∴π hemibond, the energy order among the isomers of the different hemibond motifs strongly depends on the functionals. These results demonstrate that the [benzene-(H2S) n ]+ clusters can be a benchmark of density functionals to evaluate the sulfur hemibonds.

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.

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

S-aromatic motifs are important noncovalent forces for protein stability and function but remain poorly understood. Hence, we performed quantum calculations at the MP2(full)/6-311++G(d,p) level on complexes between Cys (H₂S, MeSH) and Met (Me₂S) models with models of Phe (benzene, toluene), Trp (indole, 3-methylindole), Tyr (phenol, 4-methylphenol), and His (imidazole, 4-methylimidazole). The most stable gas-phase conformers exhibit binding energies of -2 to -6 kcal/mol, and the S atom lies perpendicular to the ring plane. This reveals preferential interaction with the ring π-system, except in the imidazoles where S binds edge-on to an N atom. Complexation tunes the gas-phase vertical ionization potentials of the ligands over as much as 1 eV, and strong σ- or π-type H-bonding supports charge transfer to the H-bond donor, rendering it more oxidizable. When the S atom acts as an H-bond acceptor (N/O-H_ar···S), calibration of the CHARMM36 force field (by optimizing pair-specific Lennard-Jones parameters) is required. Implementing the optimized parameters in molecular dynamics simulations in bulk water, we find stable S-aromatic complexes with binding free energies of -0.6 to -1.1 kcal/mol at ligand separations up to 8 Å. The aqueous S-aromatics exhibit flexible binding conformations, but edge-on conformers are less stable in water. Reflecting this, only 0.3 to 10% of the S-indole, S-phenol, and S-imidazole structures are stabilized by N/O-H_ar···S or S-H···O_ar/N_ar σ-type H-bonding. The wide range of energies and geometries found for S-aromatic interactions and their tunable redox properties expose the versatility and variability of the S-aromatic motif in proteins and allow us to predict a number of their reported properties.

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.

Insights into Thiol-Aromatic Interactions: A Stereoelectronic Basis for S-H/π Interactions

Thiols can engage favorably with aromatic rings in S-H/π interactions, within abiological systems and within proteins. However, the underlying bases for S-H/π interactions are not well understood. The crystal structure of Boc-l-4-thiolphenylalanine tert-butyl ester revealed crystal organization centered on the interaction of the thiol S-H with the aromatic ring of an adjacent molecule, with a through-space Hthiol···Caromatic distance of 2.71 Å, below the 2.90 Å sum of the van der Waals radii of H and C. The nature of this interaction was further examined by DFT calculations, IR spectroscopy, solid-state NMR spectroscopy, and analysis of the Cambridge Structural Database. The S-H/π interaction was found to be driven significantly by favorable molecular orbital interactions, between an aromatic π donor orbital and the S-H σ* acceptor orbital (a π → σ* interaction). For comparison, a structural analysis of O-H/π interactions and of cation/π interactions of alkali metal cations with aromatic rings was conducted. Na+ and K+ exhibit a significant preference for the centroid of the aromatic ring and distances near the sum of the van der Waals and ionic radii, as expected for predominantly electrostatic interactions. Li+ deviates substantially from Na+ and K+. The S-H/π interaction differs from classical cation/π interactions by the preferential alignment of the S-H σ* toward the ring carbons and an aromatic π orbital rather than toward the aromatic centroid. These results describe a potentially broadly applicable approach to understanding the interactions of weakly polar bonds with π systems.