Potential Energy Curves for Cation−π Interactions: Off-Axis Configurations Are Also Attractive

Cation-π Interactions between Quaternary Ammonium Ions and Amino Acid Aromatic Groups in Aqueous Solution

Cation-π interactions play important roles in the stabilization of protein structures and protein-ligand complexes. They contribute to the binding of quaternary ammonium ligands (mainly RNH₃⁺ and RN(CH₃)₃⁺) to various protein receptors and are likely involved in the blockage of potassium channels by tetramethylammonium (TMA⁺) and tetraethylammonium (TEA⁺). Polarizable molecular models are calibrated for NH₄⁺, TMA⁺, and TEA⁺ interacting with benzene, toluene, 4-methylphenol, and 3-methylindole (representing aromatic amino acid side chains) based on the ab initio MP2(full)/6-311++G(d,p) properties of the complexes. Whereas the gas-phase affinity of the ions with a given aromatic follows the trend NH₄⁺ > TMA⁺ > TEA⁺, molecular dynamics simulations using the polarizable models show a reverse trend in water, likely due to a contribution from the hydrophobic effect. This reversed trend follows the solubility of aromatic hydrocarbons in quaternary ammonium salt solutions, which suggests a role for cation-π interactions in the salting-in of aromatic compounds in solution. Simulations in water show that the complexes possess binding free energies ranging from -1.3 to -3.3 kcal/mol (compared to gas-phase binding energies between -8.5 and -25.0 kcal/mol). Interestingly, whereas the most stable complexes involve TEA⁺ (the largest ion), the most stable solvent-separated complexes involve TMA⁺ (the intermediate-size ion).

Improving the Force Field Description of Tyrosine-Choline Cation-π Interactions: QM Investigation of Phenol-N(Me)4+ Interactions

Cation-π interactions between tyrosine amino acids and compounds containing N,N,N-trimethylethanolammonium (N(CH3)3) are involved in the recognition of histone tails by chromodomains and in the recognition of phosphatidylcholine (PC) phospholipids by membrane-binding proteins. Yet, the lack of explicit polarization or charge transfer effects in molecular mechanics force fields raises questions about the reliability of the representation of these interactions in biomolecular simulations. Here, we investigate the nature of phenol-tetramethylammonium (TMA) interactions using quantum mechanical (QM) calculations, which we also use to evaluate the accuracy of the additive CHARMM36 and Drude polarizable force fields in modeling tyrosine-choline interactions. We show that the potential energy surface (PES) obtained using SAPT2+/aug-cc-pVDZ compares well with the large basis-set CCSD(T) PES when TMA approaches the phenol ring perpendicularly. Furthermore, the SAPT energy decomposition reveals comparable contributions from electrostatics and dispersion in phenol-TMA interactions. We then compared the SAPT2+/aug-cc-pVDZ PES obtained along various approach directions to the corresponding PES obtained with CHARMM, and we show that the force field accurately reproduces the minimum distances while the interaction energies are underestimated. The use of the Drude polarizable force field significantly improves the interaction energies but decreases the agreement on distances at energy minima. The best agreement between force field and QM PES is obtained by modifying the Lennard-Jones terms for atom pairs involved in the phenol-TMA cation-π interactions. This is further shown to improve the correlation between the occupancy of tyrosine-choline cation-π interactions obtained from molecular dynamics simulations of a bilayer-bound bacterial phospholipase and experimental affinity data of the wild-type protein and selected mutants.

Ground-State Charge Transfer: Lithium-Benzene and the Role of Hartree-Fock Exchange

Most approximations to the exchange-correlation functional of Kohn-Sham density functional theory lead to delocalization errors that undermine the description of charge-transfer phenomena. We explore how various approximate functionals and charge-distribution schemes describe ground-state atomic-charge distributions in the lithium-benzene complex, a model system of relevance to carbon-based supercapacitors. To understand the trends, we compare Hartree-Fock (HF) and correlated post-HF calculations, confirming that the HOMO-LUMO gap is narrower in semilocal functionals but widened by hybrid functionals with large fractions of HF exchange. For semilocal functionals, natural bond orbital (NBO) and Mulliken schemes yield opposite pictures of how charge transfer occurs. In PBE, for example, when lithium and benzene are <1.5 Å apart, NBO yields a positive charge on the lithium atom, but the Mulliken scheme yields a negative charge. Furthermore, the partial charges in conjugated materials depend on the interplay between the charge-distribution scheme employed and the underlying exchange-correlation functional, being critically sensitive to the admixture of HF exchange. We analyze and explain why this happens, discuss implications, and conclude that hybrid functionals with an admixture of about one-fourth of HF exchange are particularly useful in describing charge transfer in the lithium-benzene model.

Generation of self-clusters of galectin-1 in the farnesyl-bound form

Ras protein is involved in a signal transduction cascade in cell growth, and cluster formation of H-Ras and human galectin-1 (Gal-1) complex is considered to be crucial to achieve its physiological roles. It is considered that the complex is formed through interactions between Gal-1 and the farnesyl group (farnesyl-dependent model), post-translationally modified to the C-terminal Cys, of H-Ras. We investigated the role of farnesyl-bound Gal-1 in the cluster formation by analyzing the structure and properties of Gal-1 bound to farnesyl thiosalicylic acid (FTS), a competitive inhibitor of the binding of H-Ras to Gal-1. Gal-1 exhibited self-cluster formation upon interaction with FTS, and small- and large-size clusters were formed depending on FTS concentration. The galactoside-binding pocket of Gal-1 in the FTS-bound form was found to play an important role in small-size cluster formation. Large-size clusters were likely formed by the interaction among the hydrophobic sites of Gal-1 in the FTS-bound form. The present results indicate that Gal-1 in the FTS-bound form has the ability to form self-clusters as well as intrinsic lectin activity. Relevance of the self-clustering of FTS-bound Gal-1 to the cluster formation of the H-Ras–Gal-1complex was discussed by taking account of the farnesyl-dependent model and another (Raf-dependent) model.

Gas-Phase Spectroscopic Signatures of Carboxylate-Li(+) Contact Ion Pairs: New Benchmarks For Characterizing Ion Pairing in Solution

The coexistence of several types of ion pairs in solution together with their elusive nature hampers their experimental characterization, which relies in practice on theoretical models resorting to numerous approximations. In this context, a series of isolated contact ion pairs between a lithium cation and phenyl-tagged carboxylate anions of various lengths (Ph-(CH2)n-COO(-), n = 1-3) has been investigated in a conformer-selective manner by IR and UV laser spectroscopy, in conjunction with quantum chemistry calculations. The typical gas-phase IR signature of the bidentate structure formed between the carboxylate moiety and Li(+) has thus been obtained in the CO2(-) stretch region. In addition to the cation-anion interaction, a cation-π interaction occurs simultaneously in the largest system investigated (n = 3). The resulting distorted ion pair structure has been evidenced from both the IR signature of the CO2(-) stretches and the unique vibrationally resolved UV spectroscopy of a phenyl ring interacting with a cation. Such specific spectroscopic signatures of contact ion pairs provide experimental benchmarks, alternative to theoretical predictions, that can assist the assignment of vibrational spectra in solution.

Comprehensive Structural Characterization of the Bacterial Homospermidine Synthase-an Essential Enzyme of the Polyamine Metabolism

The highly conserved bacterial homospermidine synthase (HSS) is a key enzyme of the polyamine metabolism of many proteobacteria including pathogenic strains such as Legionella pneumophila and Pseudomonas aeruginosa; The unique usage of NAD(H) as a prosthetic group is a common feature of bacterial HSS, eukaryotic HSS and deoxyhypusine synthase (DHS). The structure of the bacterial enzyme does not possess a lysine residue in the active center and thus does not form an enzyme-substrate Schiff base intermediate as observed for the DHS. In contrast to the DHS the active site is not formed by the interface of two subunits but resides within one subunit of the bacterial HSS. Crystal structures of Blastochloris viridis HSS (BvHSS) reveal two distinct substrate binding sites, one of which is highly specific for putrescine. BvHSS features a side pocket in the direct vicinity of the active site formed by conserved amino acids and a potential substrate discrimination, guiding, and sensing mechanism. The proposed reaction steps for the catalysis of BvHSS emphasize cation-π interaction through a conserved Trp residue as a key stabilizer of high energetic transition states.

Structure of a TCR-Mimic Antibody with Target Predicts Pharmacogenetics

Antibody therapies currently target only extracellular antigens. A strategy to recognize intracellular antigens is to target peptides presented by immune HLA receptors. ESK1 is a human, T-cell receptor (TCR)-mimic antibody that binds with subnanomolar affinity to the RMF peptide from the intracellular Wilms tumor oncoprotein WT1 in complex with HLA-A*02:01. ESK1 is therapeutically effective in mouse models of WT1(+) human cancers. TCR-based therapies have been presumed to be restricted to one HLA subtype. The mechanism for the specificity and high affinity of ESK1 is unknown. We show in a crystal structure that ESK1 Fab binds to RMF/HLA-A*02:01 in a mode different from that of TCRs. From the structure, we predict and then experimentally confirm high-affinity binding with multiple other HLA-A*02 subtypes, broadening the potential patient pool for ESK1 therapy. Using the crystal structure, we also predict potential off-target binding that we experimentally confirm. Our results demonstrate how protein structure information can contribute to personalized immunotherapy.

An additional role for the Brønsted acid-base catalysts of mandelate racemase in transition state stabilization

Mandelate racemase (MR) catalyzes the interconversion of the enantiomers of mandelate and serves as a paradigm for understanding the enzyme-catalyzed abstraction of an α-proton from a carbon acid substrate with a high pKa. The enzyme utilizes a two-base mechanism with Lys 166 and His 297 acting as Brønsted acid and base catalysts, respectively, in the R → S reaction direction. In the S → R reaction direction, their roles are reversed. Using isothermal titration calorimetry (ITC), MR is shown to bind the intermediate/transition state (TS) analogue inhibitor benzohydroxamate (BzH) in an entropy-driven process with a value of ΔCp equal to -358 ± 3 cal mol(-1) K(-1), consistent with an increased number of hydrophobic interactions. However, MR binds BzH with an affinity that is ∼2 orders of magnitude greater than that predicted solely on the basis of hydrophobic interactions [St. Maurice, M., and Bearne, S. L. (2004) Biochemistry 43, 2524], suggesting that additional specific interactions contribute to binding. To test the hypothesis that cation-π/NH-π interactions between the side chains of Lys 166 and His 297 and the aromatic ring and/or the hydroxamate/hydroximate moiety of BzH contribute to the binding of BzH, site-directed mutagenesis was used to generate the MR variants K166M, K166C, H297N, and K166M/H297N and their binding affinity for various ligands determined using ITC. Comparison of the binding affinities of these MR variants with the intermediate/TS analogues BzH and cyclohexanecarbohydroxamate revealed that cation-π/NH-π interactions between His 297 and the hydroxamate/hydroximate moiety and the phenyl ring of BzH contribute approximately 0.26 and 0.91 kcal/mol to binding, respectively, while interactions with Lys 166 contribute approximately 1.74 and 1.74 kcal/mol, respectively. Similarly, comparison of the binding affinities of these mutants with substrate analogues revealed that Lys 166 contributes >2.93 kcal/mol to the binding of (R)-atrolactate, and His 297 contributes 2.46 kcal/mol to the binding of (S)-atrolactate. These results are consistent with Lys 166 and His 297 playing dual roles in catalysis: they act as Brønsted acid-base catalysts, and they stabilize both the enolate moiety and phenyl ring of the altered substrate in the TS.

Blockage of Water Flow in Carbon Nanotubes by Ions Due to Interactions between Cations and Aromatic Rings

Combining classical molecular dynamics simulations and density functional theory calculations, we find that cations block water flow through narrow (6,6)-type carbon nanotubes (CNTs) because of interactions between cations and aromatic rings in CNTs. In wide CNTs, these interactions trap the cations in the interior of the CNT, inducing unexpected open or closed state switching of ion transfer under a strong electric field, which is consistent with experiments. These findings will help to develop new methods to facilitate water and ion transport across CNTs.

Exploring Strong Interactions in Proteins with Quantum Chemistry and Examples of Their Applications in Drug Design

OBJECTIVES

Three strong interactions between amino acid side chains (salt bridge, cation-π, and amide bridge) are studied that are stronger than (or comparable to) the common hydrogen bond interactions, and play important roles in protein-protein interactions.

METHODS

Quantum chemical methods MP2 and CCSD(T) are used in calculations of interaction energies and structural optimizations.

RESULTS

The energies of three types of amino acid side chain interactions in gaseous phase and in aqueous solutions are calculated using high level quantum chemical methods and basis sets. Typical examples of amino acid salt bridge, cation-π, and amide bridge interactions are analyzed, including the inhibitor design targeting neuraminidase (NA) enzyme of influenza A virus, and the ligand binding interactions in the HCV p7 ion channel. The inhibition mechanism of the M2 proton channel in the influenza A virus is analyzed based on strong amino acid interactions.

CONCLUSION

(1) The salt bridge interactions between acidic amino acids (Glu- and Asp-) and alkaline amino acids (Arg+, Lys+ and His+) are the strongest residue-residue interactions. However, this type of interaction may be weakened by solvation effects and broken by lower pH conditions. (2) The cation- interactions between protonated amino acids (Arg+, Lys+ and His+) and aromatic amino acids (Phe, Tyr, Trp and His) are 2.5 to 5-fold stronger than common hydrogen bond interactions and are less affected by the solvation environment. (3) The amide bridge interactions between the two amide-containing amino acids (Asn and Gln) are three times stronger than hydrogen bond interactions, which are less influenced by the pH of the solution. (4) Ten of the twenty natural amino acids are involved in salt bridge, or cation-, or amide bridge interactions that often play important roles in protein-protein, protein-peptide, protein-ligand, and protein-DNA interactions.

Allosterism in human complement component 5a ((h)C5a): a damper of C5a receptor (C5aR) signaling

The phenomena of allosterism continues to advance the field of drug discovery, by illuminating gainful insights for many key processes, related to the structure-function relationships in proteins and enzymes, including the transmembrane G-protein coupled receptors (GPCRs), both in normal as well as in the disease states. However, allosterism is completely unexplored in the native protein ligands, especially when a small covalent change significantly modulates the pharmacology of the protein ligands toward the signaling axes of the GPCRs. One such example is the human C5a ((h)C5a), the potent cationic anaphylatoxin that engages C5aR and C5L2 to elicit numerous immunological and non-immunological responses in humans. From the recently available structure-function data, it is clear that unlike the mouse C5a ((m)C5a), the (h)C5a displays conformational heterogeneity. However, the molecular basis of such conformational heterogeneity, otherwise allosterism in (h)C5a and its precise contribution toward the overall C5aR signaling is not known. This study attempts to decipher the functional role of allosterism in (h)C5a, by exploring the inherent conformational dynamics in (m)C5a, (h)C5a and in its point mutants, including the proteolytic mutant des-Arg(74)-(h)C5a. Prima facie, the comparative molecular dynamics study, over total 500 ns, identifies Arg(74)-Tyr(23) and Arg(37)-Phe(51) "cation-π" pairs as the molecular "allosteric switches" on (h)C5a that potentially functions as a damper of C5aR signaling.

Ortho-7 bound to the active-site gorge of free and OP-conjugated acetylcholinesterase: cation-π interactions

Despite the immense importance of cation-π interactions prevailing in bispyridinium drug acetylcholinesterase (AChE) complexes, a precise description of cation-π interactions at molecular level has remained elusive. Here, we consider a bispyridinium drug, namely, ortho-7 in three different structures of AChE, with and without complexation with organophosphorus (OP) compounds for detailed investigation using all atom molecular dynamics simulation. By quantum mechanical calculations, Y72, W86, Y124, W286, Y337, and Y341 aromatic residues of the enzyme are investigated for possible cation-π interactions with ortho-7. The cation-π interactions in each of the protein-drug complexes are studied using distance, angle, a suitable functional form of them, and electrostatic criteria. The variation of cation-π functional is remarkably consistent with that of the Columbic variation. It is clearly observed that cation-π interactions for some of the residues in the catalytic active site (CAS) and peripheral anionic site (PAS) of the enzyme are either enhanced or reduced based on the nature of OP conjugation (i.e., nerve gas, tabun or pesticide, fenamiphos) when compared with the OP-free enzyme. The strength of cation-π interaction is strongly dependent on the type OP conjugation. The effect of conjugation at CAS is also seen to influence the cation-π interaction at the PAS region. The variation of cation-π interactions on the type of conjugating OP compounds might be suggestive of a reason as to why wide spectrum drug against any OP poisoning is yet to arrive in the market.

Spatial assignment of symmetry adapted perturbation theory interaction energy components: The atomic SAPT partition

We develop a physically-motivated assignment of symmetry adapted perturbation theory for intermolecular interactions (SAPT) into atom-pairwise contributions (the A-SAPT partition). The basic precept of A-SAPT is that the many-body interaction energy components are computed normally under the formalism of SAPT, following which a spatially-localized two-body quasiparticle interaction is extracted from the many-body interaction terms. For electrostatics and induction source terms, the relevant quasiparticles are atoms, which are obtained in this work through the iterative stockholder analysis (ISA) procedure. For the exchange, induction response, and dispersion terms, the relevant quasiparticles are local occupied orbitals, which are obtained in this work through the Pipek-Mezey procedure. The local orbital atomic charges obtained from ISA additionally allow the terms involving local orbitals to be assigned in an atom-pairwise manner. Further summation over the atoms of one or the other monomer allows for a chemically intuitive visualization of the contribution of each atom and interaction component to the overall noncovalent interaction strength. Herein, we present the intuitive development and mathematical form for A-SAPT applied in the SAPT0 approximation (the A-SAPT0 partition). We also provide an efficient series of algorithms for the computation of the A-SAPT0 partition with essentially the same computational cost as the corresponding SAPT0 decomposition. We probe the sensitivity of the A-SAPT0 partition to the ISA grid and convergence parameter, orbital localization metric, and induction coupling treatment, and recommend a set of practical choices which closes the definition of the A-SAPT0 partition. We demonstrate the utility and computational tractability of the A-SAPT0 partition in the context of side-on cation-π interactions and the intercalation of DNA by proflavine. A-SAPT0 clearly shows the key processes in these complicated noncovalent interactions, in systems with up to 220 atoms and 2845 basis functions.

Model structures of inactive and peptide agonist bound C5aR: Insights into agonist binding, selectivity and activation

C5a receptor (C5aR) is one of the major chemoattractant receptors of the druggable proteome that binds C5a, the proinflammatory polypeptide of complement cascade, triggering inflammation and SEPSIS. Here, we report the model structures of C5aR in both inactive and peptide agonist (YSFKPMPLaR; a=D-Ala) bound meta-active state. Assembled in CYANA and evolved over molecular dynamics (MD) in POPC bilayer, the inactive C5aR demonstrates a topologically unique compact heptahelical bundle topology harboring a β-hairpin in extracellular loop 2 (ECL2), derived from the atomistic folding simulations. The peptide agonist bound meta-active C5aR deciphers the “site2” at an atomistic resolution in the extracellular surface (ECS), in contrast to the previously hypothesized inter-helical crevice. With estimated Ki≈2.75 μM, the meta-active C5aR excellently rationalizes the IC₅₀ (0.1–13 μM) and EC₅₀ (0.01–6 μM) values, displayed by the peptide agonist in several signaling studies. Moreover, with Ki≈5.3×10⁵ μM, the “site2” also illustrates selectivity, by discriminating the stereochemical mutant peptide (YSFkPMPLaR; k=D-Lys), known to be inert toward C5aR, up to 1 mM concentration. Topologically juxtaposed between the structures of rhodopsin and CXCR1, the C5aR models also display excellent structural correlations with the other G-protein coupled receptors (GPCRs). The models elaborated in the current study unravel many important structural insights previously not known for regulating the agonist binding and activation mechanism of C5aR.

•

Topologically unique inactive and meta-active atomistic models of C5aR.

•

Model demonstrates excellent structural correlation with the other reported GPCRs.

•

Model deciphers the “site2” in the ECS and also demonstrates agonist selectivity.

•

Agonist binding and activation requires “cation–π” interaction with F275 of C5aR.

•

Inactive to meta-active transition involves TM3–TM6 movements (ΔΘ≈+11.1°) in C5aR.

Diarylferrocene tweezers for cation binding

Diarylferrocenes can act as molecular tweezers of cations. Their unique molecular shape and low torsional potentials allow for strong binding of small cations in the gas phase.

The leading role of cation–π interactions in polymer chemistry: the control of the helical sense in solution

Cation–π interactions determine the helical sense adopted by a polyphenylacetylene bearing (R)-α-methoxy-α-phenylacetamide as a pendant group (poly-1).

Cation-π interactions in protonated phenylalkylamines

Phenylalkylamines of the general formula C6H5(CH2)nNH2 (n = 1-4) have been delivered to the gas phase as protonated species using electrospray ionization. The ions thus formed have been assayed by IRMPD spectroscopy in two different spectroscopic domains, namely, the 600-1800 and the 3000-3500 cm(-1) regions using either an IR free electron laser or a tabletop OPO/OPA laser source. The interpretation of the experimental spectra is aided by density functional theory calculations of candidate species and vibrational frequency analyses. Protonated benzylamine presents a relatively straightforward instance of a single stable conformer, providing a trial case for the adopted approach. Turning to the higher homologues, C6H5(CH2)nNH3(+) (n = 2-4), more conformations become accessible. For each C6H5(CH2)nNH3(+) ion (n = 2-4), the most stable geometry is characterized by cation-π interactions between the positively charged ammonium group and the aromatic π-electronic system, permitted by the folding of the polymethylene chain. The IRMPD spectra of the sampled ions confirm the presence of the folded structures by comparison with the calculated IR spectra of the various possible conformers. An inspection of the NH stretching region is helpful in this regard.

Effect of substituent and solvent on cation-π interactions in benzene and borazine: a computational study

A DFT and ab initio quantum chemical study has been carried out at different theoretical levels to delve into the role of the cation-π interaction within the main group metal cations (Li(+), Na(+) and K(+)), substituted benzene and borazine. The effects of electron withdrawing and electron donating groups on these non-covalent forces of interaction were also studied. The excellent correlation between Hammett constants and binding energy values indicates that the cation-π interaction is influenced by both inductive and resonance effects. Electron donating groups (EDG) such as -CH3 and -NH2 attached to benzene at the 1, 3 and 5 position and the three boron atoms of borazine were found to strengthen these interactions, while electron withdrawing groups (EWG) such as -NO2 did the reverse. These results were further substantiated by topological analysis using the quantum theory of atoms in molecules (QTAIM). The polarized continuum model (PCM) and the discrete solvation model were used to elucidate the effect of solvation on the cation-π interaction. The size of the cations and the nature of the substituents were found to influence the enthalpy and binding energy of the systems (or complex). In the gas phase, the cation-π interaction was found to be exothermic, whereas in the presence of a polar solvent the interaction was highly endothermic. Thermochemical analysis predicts the presence of thermodynamic driving forces for borazine and benzene substituted with EDG. DFT based reactivity descriptors, such as global hardness (η), chemical potential (μ) and the electrophilicity index (ω) were used to elucidate the effect of the substituent on the reactivity of the cation-π complexes.

Cation-π interactions of methylated ammonium ions: a quantum mechanical study

Cation-π interactions of methylated ammonium ions play a key role in a broad range of biochemical systems. These include methyl-lysine binding proteins which bind to methylated sites on histone proteins, lysine demethylase enzymes which demethylate these sites, and neurotransmitter receptor complexes which bind choline-derived ligands. Recognition in these systems is achieved through an 'aromatic cage' motif in the binding site. Here we use high-level quantum mechanical calculations to address how cation-π interactions of methylated ammonium ions are modulated by a change in methylation state and interaction geometry. We survey methyl-lysine and choline-derived complexes in the Protein Databank to validate our results against available structural data. A quantitative description of cation-π interactions of methylated ammonium systems is critical to structure-based efforts to target methyl-lysine binding proteins and demethylase enzymes in the treatment of unregulated transcriptional control, and neurotransmitter receptors in the treatment of neurological disease. It is our hope that our work will serve as a benchmark for the development of physical chemistry based force fields that can accurately model the contribution of cation-π interactions to binding and specificity in these systems.

A mixed valent heterometallic CuII/NaI coordination polymer with sodium–phenyl bonds

A mixed valent heterometallic Cu^II/Na^I coordination polymer (1) is generated by the reaction of a Schiff base ligand, (6,6′-(1E,1′E)-(2-hydroxypropane-1,3-diyl)bis(azan-1-yl-1-ylidene)bis(methan-1-yl-1-ylidene)bis(2-methoxyphenol)) with copper(ii) acetate and sodium perchlorate.

Chemical synthesis and biological activities of novel pleuromutilin derivatives with substituted amino moiety

Novel pleuromutilin derivatives designed based on the structure of valnemulin were synthesized and evaluated for their in vitro antibacterial activities. These pleuromutilin derivatives with substituted amino moiety exhibited excellent activities against methicillin-resistant Staphylococcus aureus, methicillin-resistant Staphylococcus epidermidis, Escherichia coli, and Streptococcus agalactiae. Compound 5b showed the highest antibacterial activities and even exceeded tiamulin. Moreover, the docking experiments provided information about the binding model between the synthesized compounds and peptidyl transferase center (PTC) of 23S rRNA.

Polycation-π interactions are a driving force for molecular recognition by an intrinsically disordered oncoprotein family

Molecular recognition by intrinsically disordered proteins (IDPs) commonly involves specific localized contacts and target-induced disorder to order transitions. However, some IDPs remain disordered in the bound state, a phenomenon coined “fuzziness”, often characterized by IDP polyvalency, sequence-insensitivity and a dynamic ensemble of disordered bound-state conformations. Besides the above general features, specific biophysical models for fuzzy interactions are mostly lacking. The transcriptional activation domain of the Ewing's Sarcoma oncoprotein family (EAD) is an IDP that exhibits many features of fuzziness, with multiple EAD aromatic side chains driving molecular recognition. Considering the prevalent role of cation-π interactions at various protein-protein interfaces, we hypothesized that EAD-target binding involves polycation- π contacts between a disordered EAD and basic residues on the target. Herein we evaluated the polycation-π hypothesis via functional and theoretical interrogation of EAD variants. The experimental effects of a range of EAD sequence variations, including aromatic number, aromatic density and charge perturbations, all support the cation-π model. Moreover, the activity trends observed are well captured by a coarse-grained EAD chain model and a corresponding analytical model based on interaction between EAD aromatics and surface cations of a generic globular target. EAD-target binding, in the context of pathological Ewing's Sarcoma oncoproteins, is thus seen to be driven by a balance between EAD conformational entropy and favorable EAD-target cation-π contacts. Such a highly versatile mode of molecular recognition offers a general conceptual framework for promiscuous target recognition by polyvalent IDPs.

Understanding how proteins recognize each other is central to deciphering the inner workings of living things and for biomedical research. It has long been known that the sequence of a protein, which is a string of different amino acids, can dictate how a protein molecule folds into a well-defined shape required for biological tasks. Many folded proteins recognize and bind with each other by a tight geometric fit similar to that between a lock and its key. Recently, however, it has become clear that some proteins function as a flexible string, in constant motion, without forming a stable shape. Understanding how such “disordered” proteins work is challenging. To gain insight, we studied a disordered protein region that causes a large family of human cancers. Employing an innovative combination of experimental and theoretical techniques, we describe a new mode of protein interaction based on multiple simple contacts between one type of amino acid (aromatic) in the disordered protein and another type (positively charged) on the partner protein. Because this mechanism also underlies the ability of the disordered protein to cause cancer, further investigation of this unprecedented mode of protein-protein interaction may open up new avenues for cancer therapy.

Energy component analysis of π interactions

Fundamental features of biomolecules, such as their structure, solvation, and crystal packing and even the docking of drugs, rely on noncovalent interactions. Theory can help elucidate the nature of these interactions, and energy component analysis reveals the contributions from the various intermolecular forces: electrostatics, London dispersion terms, induction (polarization), and short-range exchange-repulsion. Symmetry-adapted perturbation theory (SAPT) provides one method for this type of analysis. In this Account, we show several examples of how SAPT provides insight into the nature of noncovalent π-interactions. In cation-π interactions, the cation strongly polarizes electrons in π-orbitals, leading to substantially attractive induction terms. This polarization is so important that a cation and a benzene attract each other when placed in the same plane, even though a consideration of the electrostatic interactions alone would suggest otherwise. SAPT analysis can also support an understanding of substituent effects in π-π interactions. Trends in face-to-face sandwich benzene dimers cannot be understood solely in terms of electrostatic effects, especially for multiply substituted dimers, but SAPT analysis demonstrates the importance of London dispersion forces. Moreover, detailed SAPT studies also reveal the critical importance of charge penetration effects in π-stacking interactions. These effects arise in cases with substantial orbital overlap, such as in π-stacking in DNA or in crystal structures of π-conjugated materials. These charge penetration effects lead to attractive electrostatic terms where a simpler analysis based on atom-centered charges, electrostatic potential plots, or even distributed multipole analysis would incorrectly predict repulsive electrostatics. SAPT analysis of sandwich benzene, benzene-pyridine, and pyridine dimers indicates that dipole/induced-dipole terms present in benzene-pyridine but not in benzene dimer are relatively unimportant. In general, a nitrogen heteroatom contracts the electron density, reducing the magnitude of both the London dispersion and the exchange-repulsion terms, but with an overall net increase in attraction. Finally, using recent advances in SAPT algorithms, researchers can now perform SAPT computations on systems with 200 atoms or more. We discuss a recent study of the intercalation complex of proflavine with a trinucleotide duplex of DNA. Here, London dispersion forces are the strongest contributors to binding, as is typical for π-π interactions. However, the electrostatic terms are larger than usual on a fractional basis, which likely results from the positive charge on the intercalator and its location between two electron-rich base pairs. These cation-π interactions also increase the induction term beyond those of typical noncovalent π-interactions.

Effects of the aromatic substitution pattern in cation-π sandwich complexes

A computational study investigating the effects of the aromatic substitution pattern on the structure and binding energies of cation-π sandwich complexes is reported. The correlation between the binding energies (Ebind) and Hammett substituent constants is approximately the same as what is observed for cation-π half-sandwich complexes. For cation-π sandwich complexes where both aromatics contain substituents the issue of relative conformation is a possible factor in the strength of the binding; however, the work presented here shows the Ebind values are approximately the same regardless of the relative conformation of the two substituted aromatics. Finally, recent computational work has shown conflicting results on whether cation-π sandwich Ebind values (Ebind,S) are approximately equal to twice the respective half-sandwich Ebind values (Ebind,HS), or if cation-π sandwich Ebind,S values are less than double the respective half-sandwich Ebind,HS values. The work presented here shows that for cation-π sandwich complexes involving substituted aromatics the Ebind,S values are less than twice the respective half-sandwich Ebind,HS values, and this is termed nonadditive. The extent to which the cation-π sandwich complexes investigated here are nonadditive is greater for B3LYP calculated values than for MP2 calculated values and for sandwich complexes with electron-donating substituents than those with electron-withdrawing groups.

Exploring Adsorption of Water and Ions on Carbon Surfaces using a Polarizable Force Field

Graphene, carbon nanotubes, and fullerenes are of great interest due to their unique properties and diverse applications in biology, molecular electronics, and materials science. Therefore, there is demand for methods that can accurately model the interface between carbon surfaces and their environment. In this letter we compare results for complexes of water, potassium ion, and chloride ion with graphene, carbon nanotube, and fullerene surfaces using a standard non-polarizable force field (OPLS-AA), a polarizable force field (OPLS-AAP), DFT, and ab initio theory. For interactions with water, OPLS-AA with the TIP3P or TIP4P water models describes the interactions with benzene (C(6)H(6)) and coronene (C(24)H(12)) well; however, for acenes larger than circumcoronene (C(54)H(18)) and especially for C(60), the interaction energies are somewhat too weak and polarization is needed. For ions interacting with carbon surfaces, inclusion of polarization is essential, and OPLS-AAP is found to perform well in comparison to the highest-level quantum mechanical methods. Overall, OPLS-AAP provides an accurate and computationally efficient force field for modeling condensed-phase systems featuring carbon surfaces.

JACOB: a dynamic database for computational chemistry benchmarking

JACOB (just a collection of benchmarks) is a database that contains four diverse benchmark studies, which in-turn included 72 data sets, with a total of 122,356 individual results. The database is constructed upon a dynamic web framework that allows users to retrieve data from the database via predefined categories. Additional flexibility is made available via user-defined text-based queries. Requested sets of results are then automatically presented as bar graphs, with parameters of the graphs being controllable via the URL. JACOB is currently available at www.wallerlab.org/jacob.