Analysis of the interactions of sulfur-containing amino acids in membrane proteins

Genetic and structural insights into the functional importance of the conserved gly-met-rich C-terminal tails in bacterial chaperonins

E. coli chaperonin GroEL forms nano-cages for protein folding. Although the chaperonin-mediated protein folding mechanism is well understood, the role of the conserved glycine and methionine-rich carboxy-terminal residues remains unclear. Bacteria with multiple chaperonins always retain at least one paralogue having the gly-met-rich C-terminus, indicating an essential conserved function. Here, we observed a stronger selection pressure on the paralogues with gly-met-rich C-termini, consistent with their ancestral functional importance. E. coli GroEL variants having mutations in their C-termini failed to functionally replace GroEL, suggesting the functional significance of the gly-met-rich C-termini. Further, our structural modelling and normal mode analysis showed that the C-terminal region shuttles between two cavity-specific conformations that correlate with the client-protein-binding apical domains, supporting C-termini's role in client protein encapsulation. Therefore, employing phylogenetic, genetic, and structural tools, we demonstrate that the gly-met-rich C-termini are functionally significant in chaperonin-mediated protein folding function. Owing to the pathogenic roles of the chaperonins having non-canonical C-termini, future investigations on the client protein selectivity will enable understanding the disease-specific client protein folding pathways and treatment options.

Design of Small Non-Peptidic Ligands That Alter Heteromerization between Cannabinoid CB1 and Serotonin 5HT2A Receptors

Activation of cannabinoid CB1 receptors (CB1R) by agonists induces analgesia but also induces cognitive impairment through the heteromer formed between CB1R and the serotonin 5HT2A receptor (5HT2AR). This side effect poses a serious drawback in the therapeutic use of cannabis for pain alleviation. Peptides designed from the transmembrane helices of CB1R, which are predicted to bind 5HT2AR and alter the stability of the CB1R-5HT2AR heteromer, have been shown to avert CB1R agonist-induced cognitive impairment while preserving analgesia. Using these peptides as templates, we have now designed nonpeptidic small molecules that prevent CB1R-5HT2AR heteromerization in bimolecular fluorescence complementation assays and the heteromerization-dependent allosteric modulations in cell signaling experiments. These results provide proof-of-principle for the design of optimized ligand-based disruptors of the CB1R-5HT2AR heteromer, opening new perspectives for in vivo studies.

Homodimerization of CB2 cannabinoid receptor triggered by a bivalent ligand enhances cellular signaling

G protein-coupled receptors (GPCRs) exist within a landscape of interconvertible conformational states and in dynamic equilibrium between monomers and higher-order oligomers, both influenced by ligand binding. Here, we show that a homobivalent ligand formed by equal chromenopyrazole moieties as pharmacophores, connected by 14 methylene units, can modulate the dynamics of the cannabinoid CB2 receptor (CB2R) homodimerization by simultaneously binding both protomers of the CB2R-CB2R homodimer. Computational and pharmacological experiments showed that one of the ligand pharmacophores binds to the orthosteric site of one protomer, and the other pharmacophore to a membrane-oriented pocket between transmembranes 1 and 7 of the partner protomer. This results in unique pharmacological properties, including increased potency in Gi-mediated signaling and enhanced recruitment of β-arrestin. Thus, by modulating dimerization dynamics, it may be possible to fine-tune CB2R activity, potentially leading to improved therapeutic outcomes.

PhI(OAc)2-Promoted 1,2-Transfer Reaction between 1,1-Disubstituted Allylic Alcohols and Thiophenols

The PhI(OAc)2-promoted 1,2-transfer reaction between allylic alcohols and thiophenols, conducted in an argon atmosphere, has proven to be effective in producing β-carbonyl sulfides from 1,1-disubstituted allylic alcohols in high yields. This method offers a fast and efficient way to synthesize β-carbonyl sulfides, which are valuable intermediates in organic synthesis. This discussion focuses on the effects of the oxidizer, temperature, and solvent on the reaction. A proposed tentative mechanism for this reaction is also discussed.

Discovering immunoreceptor coupling and organization motifs

The recently determined cryo-EM structures of the T cell antigen receptor (TCR) and B cell antigen receptor (BCR) show in molecular details the interactions of the ligand-binding part with the signaling subunits but they do not reveal the signaling mechanism of these antigen receptors. Without knowing the molecular basis of antigen sensing by these receptors, a rational design of optimal vaccines is not possible. The existence of conserved amino acids (AAs) that are not involved in the subunit interaction suggests that antigen receptors form higher complexes and/or have lateral interactors that control their activity. Here, I describe evolutionary conserved leucine zipper (LZ) motifs within the transmembrane domains (TMD) of antigen and coreceptor components that are likely to be involved in the oligomerization and lateral interaction of antigen receptor complexes on T and B cells. These immunoreceptor coupling and organization motifs (ICOMs) are also found within the TMDs of other important receptor types and viral envelope proteins. This discovery suggests that antigen receptors do not function as isolated entities but rather as part of an ICOM-based interactome that controls their nanoscale organization on resting cells and their dynamic remodeling on activated lymphocytes.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Peptide-mediated inhibition of the transcriptional regulator Elongin BC induces apoptosis in cancer cells

Inhibition of protein-protein interactions (PPIs) via designed peptides is an effective strategy to perturb their biological functions. The Elongin BC heterodimer (ELOB/C) binds to a BC-box motif and is essential for cancer cell growth. Here, we report a peptide that mimics the high-affinity BC-box of the PRC2-associated protein EPOP. This peptide tightly binds to the ELOB/C dimer (k_D = 0.46 ± 0.02 nM) and blocks the association of ELOB/C with its interaction partners, both in vitro and in the cellular environment. Cancer cells treated with our peptide inhibitor showed decreased cell viability, increased apoptosis, and perturbed gene expression. Therefore, our work proposes that blocking the BC-box-binding pocket of ELOB/C is a feasible strategy to impair its function and inhibit cancer cell growth. Our peptide inhibitor promises novel mechanistic insights into the biological function of the ELOB/C dimer and offers a starting point for therapeutics linked to ELOB/C dysfunction.

Bioderived Lipoic Acid-Based Dynamic Covalent Nanonetworks of Poly(disulfide)s: Enhanced Encapsulation Stability and Cancer Cell-Selective Delivery of Drugs

Dynamic covalent poly(disulfide)-based cross-linked nanoaggregates, termed nanonetworks (NNs), endowed with pH- and redox-responsive degradation features have been fabricated for stable noncovalent encapsulation and triggered cargo release in a controlled fashion. A bioderived lipoic acid-based Gemini surfactant-like amphiphilic molecule was synthesized for the preparation of nanoaggregates. It self-assembles by a entropy-driven self-assembly process in aqueous milieu. To further stabilize the self-assembled nanostructure, the core was cross-linked by ring-opening disulfide exchange polymerization (RODEP) of 1,2-dithiolane rings situated inside the core of the nanoaggregates. The cross-linked nanoaggregates, i.e., nanonetwork, are found to be stable in the presence of blood serum, and also, they maintain the self-assembled structure even below the critical aggregation concentration (CAC) as probed by dynamic light scattering (DLS) experiments. The nanonetwork showed almost 50% reduction in guest leakage compared to that of the nanoaggregates as shown by the release profile in the absence of stimuli, suggesting high encapsulation stability as evidenced by the fluorescence resonance energy transfer (FRET) experiment. The decross-linking of the nanonetwork occurs in response to redox and pH stimuli due to disulfide reduction and β-thioester hydrolysis, respectively, thus empowering disassembly-mediated controlled cargo release up to ∼87% for 55 h of incubation. The biological evaluation of the doxorubicin (DOX)-loaded nanonetwork revealed environment-specific surface charge modulation-mediated cancer cell-selective cellular uptake and cytotoxicity. The benign nature of the nanonetwork toward normal cells makes the system very promising in targeted drug delivery applications. Thus, the ease of synthesis, nanonetwork fabrication reproducibility, robust stability, triggered drug release in a controlled fashion, and cell-selective cytotoxicity behavior, we believe, will make the system a potential candidate in the development of robust materials for chemotherapeutic applications.

In silico investigation of uncoupling protein function in avian genomes

Introduction

The uncoupling proteins (UCPs) are involved in lipid metabolism and belong to a family of mitochondrial anionic transporters. In poultry, only one UCP homologue has been identified and experimentally shown to be associated with growth, feed conversion ratio, and abdominal fat according to its predominant expression in bird muscles. In endotherm birds, cell metabolic efficiency can be tuned by the rate of mitochondrial coupling. Thus, avUCP may be a key contributor to controlling metabolic rate during particular environmental changes.

Methods

This study aimed to perform a set of in-silico investigations primarily focused on the structural, biological, and biomimetic functions of avUCP. Thereby, using in silico genome analyses among 8 avian species (chicken, turkey, swallow, manakin, sparrow, wagtail, pigeon, and mallard) and a series of bioinformatic approaches, we provide phylogenetic inference and comparative genomics of avUCPs and investigate whether sequence variation can alter coding sequence characteristics, the protein structure, and its biological features. Complementarily, a combination of literature mining and prediction approaches was also applied to predict the gene networks of avUCP to identify genes, pathways, and biological crosstalk associated with avUCP function.

Results

The results showed the evolutionary alteration of UCP proteins in different avian species. Uncoupling proteins in avian species are highly conserved trans membrane proteins as seen by sequence alignment, physio-chemical parameters, and predicted protein structures. Taken together, avUCP has the potential to be considered a functional marker for the identification of cell metabolic state, thermogenesis, and oxidative stress caused by cold, heat, fasting, transfer, and other chemical stimuli stresses in birds. It can also be deduced that avUCP, in migrant or domestic birds, may increase heat stress resistance by reducing fatty acid transport/b-oxidation and thermoregulation alongside antioxidant defense mechanisms. The predicted gene network for avUCP highlighted a cluster of 21 genes involved in response to stress and 28 genes related to lipid metabolism and the proton buffering system. Finally, among 11 enriched pathways, crosstalk of 5 signaling pathways including MAPK, adipocytokine, mTOR, insulin, ErbB, and GnRH was predicted, indicating a possible combination of positive or negative feedback among pathways to regulate avUCP functions.

Discussion

Genetic selection for fast-growing commercial poultry has unintentionally increased susceptibility to many kinds of oxidative stress, and so avUCP could be considered as a potential candidate gene for balancing energy expenditure and reactive oxygen species production, especially in breeding programs. In conclusion, avUCP can be introduced as a pleiotropic gene that requires the contribution of regulatory genes, hormones, pathways, and genetic crosstalk to allow its finely-tuned function.

Impacts of noncovalent interactions involving sulfur atoms on protein stability, structure, folding, and bioactivity

This review discusses the various types of noncovalent interactions in which sulfur atoms participate and their effects on protein stability, structure, folding and bioactivity. Current approaches and recommendations for modelling these noncovalent interactions (in terms of both geometries and interaction energies) are highlighted.

A Single Point Mutation Blocks the Entrance of Ligands to the Cannabinoid CB2 Receptor via the Lipid Bilayer

Molecular dynamic (MD) simulations have become a common tool to study the pathway of ligand entry to the orthosteric binding site of G protein-coupled receptors. Here, we have combined MD simulations and site-directed mutagenesis to study the binding process of the potent JWH-133 agonist to the cannabinoid CB₂ receptor (CB₂R). In CB₂R, the N-terminus and extracellular loop 2 fold over the ligand binding pocket, blocking access to the binding cavity from the extracellular environment. We, thus, hypothesized that the binding pathway is a multistage process consisting of the hydrophobic ligand diffusing in the lipid bilayer to contact a lipid-facing vestibule, from which the ligand enters an allosteric site inside the transmembrane bundle through a tunnel formed between TMs 1 and 7 and finally moving from the allosteric to the orthosteric binding cavity. This pathway was experimentally validated by the Ala282^7.36Phe mutation that blocks the entrance of the ligand, as JWH-133 was not able to decrease the forskolin-induced cAMP levels in cells expressing the mutant receptor. This proposed ligand entry pathway defines transient binding sites that are potential cavities for the design of synthetic modulators.

A Census of Human Methionine-Rich Prion-like Domain-Containing Proteins

Methionine-rich prion-like proteins can regulate liquid–liquid phase separation processes in response to stresses. To date, however, very few proteins have been identified as methionine-rich prion-like. Herein, we have performed a computational survey of the human proteome to search for methionine-rich prion-like domains. We present a census of 51 manually curated methionine-rich prion-like proteins. Our results show that these proteins tend to be modular in nature, with molecular sizes significantly greater than those we would expect due to random sampling effects. These proteins also exhibit a remarkably high degree of spatial compaction when compared to average human proteins, even when protein size is accounted for. Computational evidence suggests that such a high degree of compactness might be due to the aggregation of methionine residues, pointing to a potential redox regulation of compactness. Gene ontology and network analyses, performed to shed light on the biological processes in which these proteins might participate, indicate that methionine-rich and non-methionine-rich prion-like proteins share gene ontology terms related to the regulation of transcription and translation but, more interestingly, these analyses also reveal that proteins from the methionine-rich group tend to share more gene ontology terms among them than they do with their non-methionine-rich prion-like counterparts.

Cofactor Dynamics Couples the Protein Surface to the Heme in Cytochrome c, Facilitating Electron Transfer

Electron transport through biomolecules and in biological transport networks is of great importance to bioenergetics and biocatalysis. More generally, it is of crucial importance to understand how the pathways that connect buried metallocofactors to other cofactors, and to protein surfaces, affect the biological chemistry of metalloproteins. In terms of electron transfer (ET), the strongest coupling pathways usually comprise covalent and hydrogen bonded networks, with a limited number of through-space contacts. Herein, we set out to determine the relative roles of hydrogen bonds involved in ET via an established heme-to-surface tunneling pathway in cytochrome (cyt) c (i.e., heme-W59-D60-E61-N62). A series of cyt c variants were produced where a ruthenium tris(diimine) photooxidant was placed at position 62 via covalent modification of the N62C residue. Surprisingly, variants where the H-bonding residues W59 and D60 were replaced (i.e., W59F and D60A) showed no change in ET rate from the ferrous heme to Ru(III). In contrast, changing the composition of an alternative tunneling pathway (i.e., heme-M64-N63-C62) with the M64L substitution shows a factor of 2 decrease in the rate of heme-to-Ru ET. This pathway involves a through-space tunneling step between the heme and M64 residue, and such steps are usually disfavored. To rationalize why the heme-M64-N63-C62 is preferred, molecular dynamics (MD) simulations and Pathways analysis were employed. These simulations show that the change in heme-Ru ET rates is attributed to different conformations with compressed donor-acceptor distances, by ∼2 Å in pathway distance, in the M64-containing protein as compared to the M64L protein. The change in distance is correlated with changes in the electronic coupling that are in accord with the experimentally observed heme-Ru ET rates. Remarkably, the M64L variation at the core of the protein translates to changes in cofactor dynamics at the protein surface. The surface changes identified by MD simulations include dynamic anion-π and dipole-dipole interactions. These interactions influence the strength of tunneling pathways and ET rates by facilitating decreases in through-space tunneling distances in key coupling pathways.

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.

Analysis of edible characteristics, antioxidant capacities, and phenolic pigment monomers in Lilium bulbs native to China

Lilium is cherished for its health-promoting properties in China. The bulbs of Lilium are rich in phenolic compounds, which are associated with antioxidant capacity. However, no systematic evaluation on phenolic compositions and antioxidant capacities for the edible Lilium native to China has been conducted. Herein, bulbs of 56 wild populations and three cultivars were collected. Their edible characteristics, antioxidant capacities, and pigments have been investigated and analyzed. The results showed that phenolic compounds contributed to the major colors (red, yellow and white) in Lilium bulbs. The seven phenolic pigment monomers responsible for the color of bulbs-cyanidin-3-O-rutinoside, isoquercitrin, regaloside B, regaloside C, regaloside H, regaloside A and regaloside D-were identified by the combination of HPLC-MS and NMR analysis. The population Lilium regale E. H. Wilson (Maoxian County, Sichuan Province) had the highest antioxidant capacity. According to the quantification results, Lilium bulbs with darker and redder colors possessed larger biomass, better nutrient compositions, significantly higher bioactive constituents, and higher antioxidant capacities than the three currently consumed cultivars of edible lily bulbs. Overall, these findings suggest that the mountainous area of southwest China could be the fourth source of edible lilies with the bulb-colored Lilium species.

Clustering of Aromatic Amino Acid Residues around Methionine in Proteins

Short-range, non-covalent interactions between amino acid residues determine protein structures and contribute to protein functions in diverse ways. The interactions of the thioether of methionine with the aromatic rings of tyrosine, tryptophan, and/or phenylalanine has long been discussed and such interactions are favorable on the order of 1–3 kcal mol⁻¹. Here, we carry out a new bioinformatics survey of known protein structures where we assay the propensity of three aromatic residues to localize around the [-CH₂-S-CH₃] of methionine. We term these groups “3-bridge clusters”. A dataset consisting of 33,819 proteins with less than 90% sequence identity was analyzed and such clusters were found in 4093 structures (or 12% of the non-redundant dataset). All sub-classes of enzymes were represented. A 3D coordinate analysis shows that most aromatic groups localize near the CH₂ and CH₃ of methionine. Quantum chemical calculations support that the 3-bridge clusters involve a network of interactions that involve the Met-S, Met-CH₂, Met-CH₃, and the π systems of nearby aromatic amino acid residues. Selected examples of proposed functions of 3-bridge clusters are discussed.

Evolutionary plasticity of SH3 domain binding by Nef proteins of the HIV-1/SIVcpz lentiviral lineage

The accessory protein Nef of human and simian immunodeficiency viruses (HIV and SIV) is an important pathogenicity factor known to interact with cellular protein kinases and other signaling proteins. A canonical SH3 domain binding motif in Nef is required for most of these interactions. For example, HIV-1 Nef activates the tyrosine kinase Hck by tightly binding to its SH3 domain. An archetypal contact between a negatively charged SH3 residue and a highly conserved arginine in Nef (Arg77) plays a key role here. Combining structural analyses with functional assays, we here show that Nef proteins have also developed a distinct structural strategy—termed the "R-clamp”—that favors the formation of this salt bridge via buttressing Arg77. Comparison of evolutionarily diverse Nef proteins revealed that several distinct R-clamps have evolved that are functionally equivalent but differ in the side chain compositions of Nef residues 83 and 120. Whereas a similar R-clamp design is shared by Nef proteins of HIV-1 groups M, O, and P, as well as SIVgor, the Nef proteins of SIV from the Eastern chimpanzee subspecies (SIVcpz^P.t.s.) exclusively utilize another type of R-clamp. By contrast, SIV of Central chimpanzees (SIVcpz^P.t.t.) and HIV-1 group N strains show more heterogenous R-clamp design principles, including a non-functional evolutionary intermediate of the aforementioned two classes. These data add to our understanding of the structural basis of SH3 binding and kinase deregulation by Nef, and provide an interesting example of primate lentiviral protein evolution.

Viral replication depends on interactions with a plethora of host cell proteins. Cellular protein interactions are typically mediated by specialized binding modules, such as the SH3 domain. To gain access to host cell regulation viruses have evolved to contain SH3 domain binding sites in their proteins, a notable example of which is the HIV-1 Nef protein. Here we show that during the primate lentivirus evolution the structural strategy that underlies the avid binding of Nef to cellular SH3 domains, which we have dubbed the R-clamp, has been generated via alternative but functionally interchangeable molecular designs. These patterns of SH3 recognition depend on the amino acid combinations at the positions corresponding to residues 83 and 120 in the consensus HIV-1 Nef sequence, and are distinctly different in Nef proteins from SIVs of Eastern and Central chimpanzees, gorillas, and the four groups of HIV-1 that have independently originated from the latter two. These results highlight the evolutionary plasticity of viral proteins, and have implications on therapeutic development aiming to interfere with SH3 binding of Nef.

Altered Peptide Self-Assembly and Co-Assembly with DNA by Modification of Aromatic Residues

Aromatic residues are widely used as building blocks for driving self-assemblies in natural and designer biomaterials. The noncovalent interactions involving aromatic rings determine proteins' structure and biofunction. Here, we studied the effects of changes in the proximity of the aromatic rings in a self-assembling peptide for modulating interactions involving the aromatic residues. By changing the distance between the aromatic ring and peptide backbone and replacing the side chain with a sulfur atom, we altered the nanostructures and gene transfection efficiency of peptide-DNA co-assemblies. This study demonstrates the significance of subtle alterations in aromatic interactions and facilitates deeper understanding of the aromatic-involving interactions.

Contributions of methionine to recognition of trimethyllysine in aromatic cage of PHD domains: implications of polarizability, hydrophobicity, and charge on binding

Recognition of trimethyllysine (Kme3) by reader proteins is an important regulator of gene expression. This recognition event is mediated by an aromatic cage made up of 2-4 aromatic residues in the reader proteins that bind Kme3 via cation-π interactions. A small subset of reader proteins contain a methionine (Met) residue in place of an aromatic sidechain in the binding pocket. The unique role of sulfur in molecular recognition has been demonstrated in a number of noncovalent interactions recently, including interactions of thiols, thioethers, and sulfoxides with aromatic rings. However, the interaction of a thioether with an ammonium ion has not previously been investigated and the role of Met in binding Kme3 has not yet been explored. Herein, we systematically vary the Met in two reader proteins, DIDO1 and TAF3, and the ligand, Kme3 or its neutral analog tert-butyl norleucine (tBuNle), to determine the role of Met in the recognition of the cationic Kme3. Our studies demonstrate that Met contributes to binding via dispersion forces, with about an equal contribution to binding Kme3 and tBuNle, indicating that electrostatic interactions do not play a role. During the course of these studies, we also discovered that DIDO1 exhibits equivalent binding to tBuNle and Kme3 through a change in the mechanism of binding.

Effect of Methionine Sulfoxide on the Synthesis and Purification of Aggregation-Prone Peptides

A two-step synthesis for methionine-containing hydrophobic and/or aggregation-prone peptides is presented that takes advantage of the reversibility of methionine oxidation. The use of polar methionine sulfoxide as a building block in solid-phase peptide synthesis improves the synthesis quality and yields the crude peptide, with significantly improved solubility compared to the reduced species. This facilitates the otherwise often laborious peptide purification by high-performance liquid chromatography. The subsequent reduction proceeds quantitatively. This approach has been optimised with the methionine-rich Tar-DNA-binding protein 43 (307-347), but is also more generally applicable, as demonstrated by the syntheses of human calcitonin and two aggregation-prone peptides from the human prion protein.

Interrogating permeation and gating of Orai channels using chemical modification of cysteine residues

Chemical modification of ion channels using the substituted cysteine accessibility method has a rich and successful history in elucidating the structural basis of ion channel function. In this approach, cysteine residues are introduced in regions of interest into the protein and their accessibility to water soluble thiol-reactive reagents is determined by monitoring ion channel activity. Because a wide range of these reagents are available with differing size, charge, and membrane solubility, the physio-chemical environment of the introduced cysteine residue and therefore the protein domain of interest can be probed with great precision. The approach has been widely employed for determining the secondary structure of specific ion channel domains, the location and nature of the channel gate, and the conformational rearrangements in the channel pore that underlie the opening/closing of the pore. In this chapter, we describe the use of these and related approaches to probe the functional architecture and gating of store-operated Orai1 channels.

The 1.3-Å resolution structure of bovine cytochrome c oxidase suggests a dimerization mechanism

Cytochrome c oxidase (CcO) in the respiratory chain catalyzes oxygen reduction by coupling electron and proton transfer through the enzyme and proton pumping across the membrane. Although the functional unit of CcO is monomeric, mitochondrial CcO forms a monomer and a dimer, as well as a supercomplex with respiratory complexes I and III. A recent study showed that dimeric CcO has lower activity than monomeric CcO and proposed that dimeric CcO is a standby form for enzymatic activation in the mitochondrial membrane. Other studies have suggested that the dimerization is dependent on specifically arranged lipid molecules, peptide segments, and post-translationally modified amino acid residues. To re-examine the structural basis of dimerization, we improved the resolution of the crystallographic structure to 1.3 Å by optimizing the method for cryoprotectant soaking. The observed electron density map revealed many weakly bound detergent and lipid molecules at the interface of the dimer. The dimer interface also contained hydrogen bonds with tightly bound cholate molecules, hydrophobic interactions between the transmembrane helices, and a Met-Met interaction between the extramembrane regions. These results imply that binding of physiological ligands structurally similar to cholate could trigger dimerization in the mitochondrial membrane and that non-specifically bound lipid molecules at the transmembrane surface between monomers support the stabilization of the dimer. The weak interactions involving the transmembrane helices and extramembrane regions may play a role in positioning each monomer at the correct orientation in the dimer.

A sulfur-aromatic gate latch is essential for opening of the Orai1 channel pore

Sulfur-aromatic interactions occur in the majority of protein structures, yet little is known about their functional roles in ion channels. Here, we describe a novel molecular motif, the M101 gate latch, which is essential for gating of human Orai1 channels via its sulfur-aromatic interactions with the F99 hydrophobic gate. Molecular dynamics simulations of different Orai variants reveal that the gate latch is mostly engaged in open but not closed channels. In experimental studies, we use metal-ion bridges to show that promoting an M101-F99 bond directly activates Orai1, whereas disrupting this interaction triggers channel closure. Mutational analysis demonstrates that the methionine residue at this position has a unique combination of length, flexibility, and chemistry to act as an effective latch for the phenylalanine gate. Because sulfur-aromatic interactions provide additional stabilization compared to purely hydrophobic interactions, we infer that the six M101-F99 pairs in the hexameric channel provide a substantial energetic contribution to Orai1 activation.

Experimental and computational analysis of biased agonism on full-length and a C-terminally truncated adenosine A2A receptor

Biased agonism, the ability of agonists to differentially activate downstream signaling pathways by stabilizing specific receptor conformations, is a key issue for G protein-coupled receptor (GPCR) signaling. The C-terminal domain might influence this functional selectivity of GPCRs as it engages G proteins, GPCR kinases, β-arrestins, and several other proteins. Thus, the aim of this paper is to compare the agonist-dependent selectivity for intracellular pathways in a heterologous system expressing the full-length (A_2AR) and a C-tail truncated (A_2A^Δ40R lacking the last 40 amino acids) adenosine A_2A receptor, a GPCR that is already targeted in Parkinson’s disease using a first-in-class drug. Experimental data such as ligand binding, cAMP production, β-arrestin recruitment, ERK1/2 phosphorylation and dynamic mass redistribution assays, which correspond to different aspects of signal transduction, were measured upon the action of structurally diverse compounds (the agonists adenosine, NECA, CGS-21680, PSB-0777 and LUF-5834 and the SCH-58261 antagonist) in cells expressing A_2AR and A_2A^Δ40R. The results show that taking cAMP levels and the endogenous adenosine agonist as references, the main difference in bias was obtained with PSB-0777 and LUF-5834. The C-terminus is dispensable for both G-protein and β-arrestin recruitment and also for MAPK activation. Unrestrained molecular dynamics simulations, at the μs timescale, were used to understand the structural arrangements of the binding cavity, triggered by these chemically different agonists, facilitating G protein binding with different efficacy.

The Prodigious Hydrogen Bonds with Sulfur and Selenium in Molecular Assemblies, Structural Biology, and Functional Materials

Hydrogen bonds (H-bonds) play important roles in imparting functionality to the basic molecules of life by stabilizing their structures and directing their interactions. Numerous studies have been devoted to understanding H-bonds involving highly electronegative atoms like nitrogen, oxygen, and halogens and consequences of those H-bonds in chemical reactions, catalysis, and structure and function of biomolecules; but the involvement of less electronegative atoms like sulfur and selenium in H-bond formation establishes the concept of noncanonical H-bonds. Initially belittled for the "weak" nature of their interactions, these perceptions have gradually evolved over time through dedicated efforts by several research groups. This has been facilitated by advancements in experimental methods for their detection through gas-phase laser spectroscopy and solution NMR spectroscopy, as well as through theoretical predictions from high level quantum chemical calculations.In this Account, we present insights into the versatility of the sulfur and selenium centered H-bonds (S/SeCHBs) by highlighting their multifarious applications in various fields from chemical reactions to optoelectronic properties to structural biology. Our group has highlighted the significance and strength of such H-bonds in natural and modified biomolecules. Here, we have reviewed several molecular assemblies, biomolecules, and functional materials, where the role of these H-bonds is pivotal in influencing biological functions. It is worth mentioning here that the precise experimental data obtained from gas-phase laser spectroscopy have contributed considerably to changing the existing perceptions toward S/SeCHBs. Thus, molecular beam experiments, though difficult to perform on smaller model thio- or seleno-substituted Molecules, etc. (amides, nucleobases, drug molecules), are inevitable to gather elementary knowledge and convincing concepts on S/SeCHBs that can be extended from a small four-atom sulfanyl dimer to a large 14 kDa iron-sulfur protein, ferredoxin. These H-bonds can also tailor a fascinating array of molecular frameworks and design supramolecular assemblies by inter- and intralinking of individual "molecular Lego-like" units.The discussion is indeed intriguing when it turns to the usage of S/SeCHBs in facile synthetic strategies like tuning regioselectivity in reactions, as well as invoking phenomena like dual phosphorescence and chemiluminescence. This is in addition to our investigations of the dispersive nature of the hydrogen bond between metal hydrides and sulfur or selenium as acceptor, which we anticipate would lead to progress in the areas of proton and hydride transfer, as well as force-field design. This Account demonstrates how ease of fabrication, enhanced efficiency, and alteration of physicochemical properties of several functional materials is facilitated owing to the presence of S/SeCHBs. Our efforts have been instrumental in the evaluation of various S/SeCHBs in flue gas capture, as well as design of organic energy harvesting materials, where dipole moment and polarizability have important roles to play. We hope this Account invokes newer perspectives with regard to how H-bonds with sulfur and selenium can be adequately adopted for crystal engineering, for more photo- and biophysical studies with different spectroscopic methods, and for developing next-generation field-effect transistors, batteries, superconductors, and organic thin-film transistors, among many other multifunctional materials for the future.

Two-stage electro-mechanical coupling of a KV channel in voltage-dependent activation

In voltage-gated potassium (K_V) channels, the voltage-sensing domain (VSD) undergoes sequential activation from the resting state to the intermediate state and activated state to trigger pore opening via electro-mechanical (E-M) coupling. However, the spatial and temporal details underlying E-M coupling remain elusive. Here, utilizing K_V7.1's unique two open states, we report a two-stage E-M coupling mechanism in voltage-dependent gating of K_V7.1 as triggered by VSD activations to the intermediate and then activated state. When the S4 segment transitions to the intermediate state, the hand-like C-terminus of the VSD-pore linker (S4-S5L) interacts with the pore in the same subunit. When S4 then proceeds to the fully-activated state, the elbow-like hinge between S4 and S4-S5L engages with the pore of the neighboring subunit to activate conductance. This two-stage hand-and-elbow gating mechanism elucidates distinct tissue-specific modulations, pharmacology, and disease pathogenesis of K_V7.1, and likely applies to numerous domain-swapped K_V channels.

Identification of seven novel variants in the β-globin gene in transfusion-dependent and normal patients

INTRODUCTION

Abnormality in HBB results in an inherited recessive blood disorder, which can be caused by variants at the transcriptional or translational level affecting the stability and the production of the HBB chain. The severity of the disease relies on the variant's characteristics. This study aimed to identify the common β-globin HBB variants in the population of the Eastern Province, which has the highest prevalence of blood diseases in Saudi Arabia.

MATERIAL AND METHODS

Direct sequence of β-globin HBB gene, and alpha-globin HBA1 and HBA2 genes was performed on a total of 545 blood samples (transfusion-dependent: 215, 106 men and 109 women; normal healthy subjects: 330, 197 men and 133 women) collected from Saudi Arabian participants in the Eastern region.

RESULTS

A total of 36 variants in HBB gene were revealed with 11 variants that have been reported for the first time in Saudi Arabia, including 7 novel variants that have been identified for the first time in HBB gene. The novel variants consisted of two exonic (HBB:c.252C>T; HBB:c.281G>T) and five intronic variants (c.316-183_316-168del; c.315+241T>A; c.315+376T>C; c.316-114C>G; c.315+208T>G) at HBB gene. The novel exonic variants and three (c.316-183_316-168del; c.315+241T>A; c.315+376T>C) intronic variants were co-inherited with α deletion.

CONCLUSIONS

This current study updated the HBB gene variations with newly identified variants of HBB gene and co-inheritance with α-globin deletions. The identified β-globin mutations will strengthen the genetic reference that could aid in characterizing mutations that are associated with phenotype of thalassemia in a specific region.

The functionally relevant site for paxilline inhibition of BK channels

The tremorgenic fungal alkaloid paxilline (PAX) is a commonly used specific inhibitor of the large-conductance, voltage- and Ca²⁺-dependent BK-type K⁺ channel. PAX inhibits BK channels by selective interaction with closed states. BK inhibition by PAX is best characterized by the idea that PAX gains access to the channel through the central cavity of the BK channel, and that only a single PAX molecule can interact with the BK channel at a time. The notion that PAX reaches its binding site via the central cavity and involves only a single PAX molecule would be consistent with binding on the axis of the permeation pathway, similar to classical open channel block and inconsistent with the observation that PAX selectively inhibits closed channels. To explore the potential sites of interaction of PAX with the BK channel, we undertook a computational analysis of the interaction of PAX with the BK channel pore gate domain guided by recently available liganded (open) and metal-free (closed) Aplysia BK channel structures. The analysis unambiguously identified a preferred position of PAX occupancy that accounts for all previously described features of PAX inhibition, including state dependence, G311 sensitivity, stoichiometry, and central cavity accessibility. This PAX-binding pose in closed BK channels is supported by additional functional results.

Strong Enrichment of Aromatic and Sulfur-Containing Residues in Ligand-Protein Binding Sites

While certain residues have clear involvement in determining the 3D structure of a macromolecule because they affect the folding topology or the overall protein stability, the role of different residues in ligand accommodation and binding has attracted less attention. On the basis of the assumption that drug-binding sites on target molecules have specific amino acid compositions, the incidence of each standard amino acid at the binding sites of small molecules and their correlations are calculated for an unprecedented large set of high-quality X-ray structures. Results show, for the first time, strong and highly correlated enrichments of aromatic and sulfur-containing residues, which play an important role in ligand binding and shape the nature of the chemical interactions.

Sequence specificity despite intrinsic disorder: How a disease-associated Val/Met polymorphism rearranges tertiary interactions in a long disordered protein

The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) in precursor brain-derived neurotrophic factor (BDNF) is one of the earliest SNPs to be associated with neuropsychiatric disorders, and the underlying molecular mechanism is unknown. Here we report on over 250 μs of fully-atomistic, explicit solvent, temperature replica-exchange molecular dynamics (MD) simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence. The simulations were able to correctly reproduce the location of both local and non-local secondary structure changes due to the Val66Met mutation, when compared with NMR spectroscopy. We find that the change in local structure is mediated via entropic and sequence specific effects. We developed a hierarchical sequence-based framework for analysis and conceptualization, which first identifies “blobs” of 4-15 residues representing local globular regions or linkers. We use this framework within a novel test for enrichment of higher-order (tertiary) structure in disordered proteins; the size and shape of each blob is extracted from MD simulation of the real protein (RP), and used to parameterize a self-avoiding heterogenous polymer (SAHP). The SAHP version of the BDNF prodomain suggested a protein segmented into three regions, with a central long, highly disordered polyampholyte linker separating two globular regions. This effective segmentation was also observed in full simulations of the RP, but the Val66Met substitution significantly increased interactions across the linker, as well as the number of participating residues. The Val66Met substitution replaces β-bridging between V66 and V94 (on either side of the linker) with specific side-chain interactions between M66 and M95. The protein backbone in the vicinity of M95 is then free to form β-bridges with residues 31-41 near the N-terminus, which condenses the protein. A significant role for Met/Met interactions is consistent with previously-observed non-local effects of the Val66Met SNP, as well as established interactions between the Met66 sequence and a Met-rich receptor that initiates neuronal growth cone retraction.

Intrinsically disordered proteins are proteins that have no well-defined structure in at least one functional form. Mutations in one amino acid may still affect their function significantly, especially in subtle ways with cumulative adverse effects on health. Here we report on molecular dynamics simulations of a protein that is critical for neuronal health throughout adulthood (brain-derived neurotrophic factor). We investigate the effects of a mutation carried by 30% of human population, which has been widely studied for its association with aging-related and stress-related disorders, reduced volume of the hippocampus, and variations in episodic memory. We identify a molecular mechanism in which the mutation may change the global conformations of the protein and its ability to bind to receptors.

A comprehensive computational study of amino acid interactions in membrane proteins

Transmembrane proteins play a fundamental role in a wide series of biological processes but, despite their importance, they are less studied than globular proteins, essentially because their embedding in lipid membranes hampers their experimental characterization. In this paper, we improved our understanding of their structural stability through the development of new knowledge-based energy functions describing amino acid pair interactions that prevail in the transmembrane and extramembrane regions of membrane proteins. The comparison of these potentials and those derived from globular proteins yields an objective view of the relative strength of amino acid interactions in the different protein environments, and their role in protein stabilization. Separate potentials were also derived from α-helical and β-barrel transmembrane regions to investigate possible dissimilarities. We found that, in extramembrane regions, hydrophobic residues are less frequent but interactions between aromatic and aliphatic amino acids as well as aromatic-sulfur interactions contribute more to stability. In transmembrane regions, polar residues are less abundant but interactions between residues of equal or opposite charges or non-charged polar residues as well as anion-π interactions appear stronger. This shows indirectly the preference of the water and lipid molecules to interact with polar and hydrophobic residues, respectively. We applied these new energy functions to predict whether a residue is located in the trans- or extramembrane region, and obtained an AUC score of 83% in cross validation, which demonstrates their accuracy. As their application is, moreover, extremely fast, they are optimal instruments for membrane protein design and large-scale investigations of membrane protein stability.

Inter-residue interactions in alpha-helical transmembrane proteins

MOTIVATION

The number of available membrane protein structures has markedly increased in the last years and, in parallel, the reliability of the methods to detect transmembrane (TM) segments. In the present report, we characterized inter-residue interactions in α-helical membrane proteins using a dataset of 3462 TM helices from 430 proteins. This is by far the largest analysis published to date.

RESULTS

Our analysis of residue-residue interactions in TM segments of membrane proteins shows that almost all interactions involve aliphatic residues and Phe. There is lack of polar-polar, polar-charged and charged-charged interactions except for those between Thr or Ser sidechains and the backbone carbonyl of aliphatic and Phe residues. The results are discussed in the context of the preferences of amino acids to be in the protein core or exposed to the lipid bilayer and to occupy specific positions along the TM segment. Comparison to datasets of β-barrel membrane proteins and of α-helical globular proteins unveils the specific patterns of interactions and residue composition characteristic of α-helical membrane proteins that are the clue to understanding their structure.

AVAILABILITY AND IMPLEMENTATION

Results data and datasets used are available at http://lmc.uab.cat/TMalphaDB/interactions.php.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Structure-function relationships of nonviral gene vectors: Lessons from antimicrobial polymers

In recent years, substantial advances have been achieved in the design and synthesis of nonviral gene vectors. However, lack of effective and biocompatible vectors still remains a major challenge that hinders their application in clinical settings. In the past decade, there has been a rapid expansion of cationic antimicrobial polymers, due to their potent, rapid, and broad-spectrum biocidal activity against resistant microbes, and biocompatible features. Given that antimicrobial polymers share common features with nonviral gene vectors in various aspects, such as membrane affinity, functional groups, physicochemical characteristics, and unique macromolecular architectures, these polymers may provide us with inspirations to overcome challenges in the design of novel vectors toward more safe and efficient gene delivery in clinic. Building off these observations, we provide here an overview of the structure-function relationships of polymers for both antimicrobial applications and gene delivery by elaborating some key structural parameters, including functional groups, charge density, hydrophobic/hydrophilic balance, MW, and macromolecular architectures. By borrowing a leaf from antimicrobial agents, great advancement in the development of newer nonviral gene vectors with high transfection efficiency and biocompatibility will be more promising. STATEMENT OF SIGNIFICANCE: The development of gene delivery is still in the preclinical stage for the lack of effective and biocompatible vectors. Given that antimicrobial polymers share common features with gene vectors in various aspects, such as membrane affinity, functional groups, physicochemical characteristics, and unique macromolecular architectures, these polymers may provide us with inspirations to overcome challenges in the design of novel vectors toward more safe and efficient gene delivery in clinic. In this review, we systematically summarized the structure-function relationships of antimicrobial polymers and gene vectors, with which the design of more advanced nonviral gene vectors is anticipated to be further boosted in the future.

Tipping the Balance between S-π and O-π Interactions

A comprehensive experimental survey consisting of 36 molecular balances was conducted to compare 18 pairs of S-π versus O-π interactions over a wide range of structural, geometric, and solvent parameters. A strong linear correlation was observed between the folding energies of the sulfur and oxygen balances across the entire library of balance pairs. The more stable interaction systematically switched from the O-π to S-π interaction. Computational studies of bimolecular PhSCH₃-arene and PhOCH₃-arene complexes were able to replicate the experimental trends in the molecular balances. The change in preference for the O-π to S-π interaction was due to the interplay of stabilizing (dispersion and solvophobic) and destabilizing (exchange-repulsion) terms arising from the differences in size and polarizability of the oxygen and sulfur atoms.

Silver nanoparticles as antimicrobial therapeutics: current perspectives and future challenges

Utility of silver metal in antimicrobial therapy is an accepted practice since ages that faded with time because of the identification of a few silver resistant strains in the contemporary era. A successive development of antibiotics soon followed. However, due to an indiscriminate and unregulated use coupled with poor legal control measures and a dearth of expertise in handling the critical episodes, the antibiotics era has already seen a steep decline in the past decades due to the evolution of multi-drug resistant 'superbugs' which pose a sizeable challenge to manage with. Due to limited options in the pipeline and no clear strategy in the forefront, the aspirations for novel, MDR focused drug discovery to target the 'superbugs' arose which once again led to the rise of AgNPs in antimicrobial research. In this review, we have focused on the green routes for the synthesis of AgNPs, the mode of microbial inhibition by AgNPs, synergistic effect of AgNPs with antibiotics and future challenges for the development of nano-silver-based therapeutics.

2-(4'-Pyridyl-N-oxide)-Substituted Hemithioindigos as Photoresponsive Guests for a Super Aryl-Extended Calix[4]pyrrole Receptor

We report the synthesis of two 2-(4'-pyridyl-N-oxide)-substituted hemithioindigos (HTIs). We probed their photoisomerization by using UV/Vis and ¹ H NMR spectroscopy techniques. Light irradiation at λ=450 nm provoked the isomerization of the HTI Z isomer to the E counterpart to a large extent (≈80 % at the photostationary state). ¹ H NMR titration experiments revealed the formation of thermodynamically and kinetically stable 1:1 inclusion complexes of the (Z)-HTI isomers with a super aryl-extended host (association constant>10⁴  m^-1 ). Photoirradiation at λ=450 nm of the inclusion complexes induced the isomerization of the bound HTI N-oxide to afford the (E)-HTI⊂calix[4]pyrrole complex. We determined accurate association constant values for the 1:1 inclusion complexes of the (Z)- and (E)-HTI isomers by using isothermal titration calorimetry experiments. The results showed that the stability constants of the (E)-HTI complexes were 2.2-2.8-fold lower than those of the (Z)-HTI counterparts, which explains the lack of light-induced release of the former to the bulk solution.

Rosetta Broker for membrane protein structure prediction: concentrative nucleoside transporter 3 and corticotropin-releasing factor receptor 1 test cases

Background

Membrane proteins are difficult targets for structure prediction due to the limited structural data deposited in Protein Data Bank. Most computational methods for membrane protein structure prediction are based on the comparative modeling. There are only few de novo methods targeting that distinct protein family. In this work an example of such de novo method was used to structurally and functionally characterize two representatives of distinct membrane proteins families of solute carrier transporters and G protein-coupled receptors. The well-known Rosetta program and one of its protocols named Broker was used in two test cases. The first case was de novo structure prediction of three N-terminal transmembrane helices of the human concentrative nucleoside transporter 3 (hCNT3) homotrimer belonging to the solute carrier 28 family of transporters (SLC28). The second case concerned the large scale refinement of transmembrane helices of a homology model of the corticotropin-releasing factor receptor 1 (CRFR1) belonging to the G protein-coupled receptors family.

Results

The inward-facing model of the hCNT3 homotrimer was used to propose the functional impact of its single nucleotide polymorphisms. Additionally, the 100 ns molecular dynamics simulation of the unliganded hCNT3 model confirmed its validity and revealed mobility of the selected binding site and homotrimer interface residues. The large scale refinement of transmembrane helices of the CRFR1 homology model resulted in the significant improvement of its accuracy with respect to the crystal structure of CRFR1, especially in the binding site area. Consequently, the antagonist CP-376395 could be docked with Autodock VINA to the CRFR1 model without any steric clashes.

Conclusions

The presented work demonstrated that Rosetta Broker can be a versatile tool for solving various issues referring to protein biology. Two distinct examples of de novo membrane protein structure prediction presented here provided important insights into three major areas of protein biology. Namely, the dynamics of the inward-facing hCNT3 homotrimer system, the structural changes of the CRFR1 receptor upon the antagonist binding and finally, the role of single nucleotide polymorphisms in both, hCNT3 and CRFR1 proteins, were investigated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0078-8) contains supplementary material, which is available to authorized users.

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.