Contribution of cation-pi interactions to the stability of protein-DNA complexes

AQcalc: A web server that identifies weak molecular interactions in protein structures

Weak molecular interactions play an important role in protein structure and function. Computational tools that identify weak molecular interactions are, therefore, valuable for the study of proteins. Here, we present AQcalc, a web server (https://aqcalcbiocomputing.com/) that can be used to identify anion-quadrupole (AQ) interactions, which are weak interactions involving aromatic residue (Trp, Tyr, and Phe) ring edges and anions (Asp, Glu, and phosphate ion) both within proteins and at their interfaces (protein-protein, protein-nucleic acids, and protein-lipid bilayer). AQcalc identifies AQ interactions as well as clusters involving AQ, cation-π, and salt bridges, among others. Utilizing AQcalc we analyzed weak interactions in protein models, even in the absence of experimental structures, to understand the contributions of weak interactions to deleterious structural changes, including those associated with oncogenic and germline disease variants. We identified several deleterious variants with disrupted AQ interactions (comparable in frequency to cation-π disruptions). Amyloid fibrils utilize AQ to bury anions at frequencies that far exceed those observed for globular proteins. AQ interactions were detected three and five times more frequently than the hydrogen-bonded AQ (HBAQ) in fibril structures and protein-lipid bilayer interfaces, respectively. By contrast, AQ and HBAQ interactions were detected with similar frequencies in globular proteins. Collectively, these findings suggest AQcalc will be effective in facilitating fine structural analysis. As other web utilities designed to identify protein residue interaction networks do not report AQ interactions, wide use of AQcalc will enrich our understanding of residue interaction networks and facilitate hypothesis testing by identifying and experimentally characterizing these comparably weak but important interactions.

Nanoporous Crystalline Materials for the Recognition and Applications of Nucleic Acids

Nucleic acid plays a crucial role in countless biological processes. Hence, there is great interest in its detection and analysis in various fields from chemistry, biology, to medicine. Nanoporous crystalline materials exhibit enormous potential as an effective platform for nucleic acid recognition and application. These materials have highly ordered and uniform pore structures, as well as adjustable surface chemistry and pore size, making them good carriers for nucleic acid extraction, detection, and delivery. In this review, the latest developments in nanoporous crystalline materials, including metal organic frameworks (MOFs), covalent organic frameworks (COFs), and supramolecular organic frameworks (SOFs) for nucleic acid recognition and applications are discussed. Different strategies for functionalizing these materials are explored to specifically identify nucleic acid targets. Their applications in selective separation and detection of nucleic acids are highlighted. They can also be used as DNA/RNA sensors, gene delivery agents, host DNAzymes, and in DNA-based computing. Other applications include catalysis, data storage, and biomimetics. The development of novel nanoporous crystalline materials with enhanced biocompatibility has opened up new avenues in the fields of nucleic acid analysis and therapy, paving the way for the development of sensitive, selective, and cost-effective diagnostic and therapeutic tools with widespread applications.

Algebraic Diagrammatic Construction Theory for Simulating Charged Excited States and Photoelectron Spectra

Charged excitations are electronic transitions that involve a change in the total charge of a molecule or material. Understanding the properties and reactivity of charged species requires insights from theoretical calculations that can accurately describe orbital relaxation and electron correlation effects in open-shell electronic states. In this Review, we describe the current state of algebraic diagrammatic construction (ADC) theory for simulating charged excitations and its recent developments. We start with a short overview of ADC formalism for the one-particle Green's function, including its single- and multireference formulations and extension to periodic systems. Next, we focus on the capabilities of ADC methods and discuss recent findings about their accuracy for calculating a wide range of excited-state properties. We conclude our Review by outlining possible directions for future developments of this theoretical approach.

Screening, characterization and specific binding mechanism of aptamers against human plasminogen Kringle 5

Plasminogen Kringle 5 is one of the most potent cytokines identified to inhibit the proliferation and migration of vascular endothelial cells. Herein, six aptamer candidates that specifically bind to Kringle 5 were generated by the systematic evolution of ligands by exponential enrichment (SELEX). After 10 rounds of screening against Kringle 5, a highly enriched ssDNA pool was sequenced and the representative aptamers were subjected to binding assays to evaluate their affinity and specificity. The preferred aptamer KG-4, which demonstrated a low dissociation constant (K_d) of ∼ 432 nM and excellent selectivity for Kringle 5. A conserved "motif" of eight bases located at the stem-loop intersection, common to the aptamer, was further confirmed as the recognition element for binding with Kringle 5. The bulge formed by the motif and depression on the lysine binding site of Kringle 5 were both located at the binding interface, and the "induced fit" between their structures played a central role in the recognition process. Kringle 5 interacts KG-4 primarily through enthalpy-driven van der Waals forces and hydrogen bond. The key nucleotides A34 and C35 at motif on KG-4 and the positively charged amino acids in the loop 1 and loop 4 regions on Kringle 5 play a major role in the interaction. Furthermore, KG-4 dose-dependently reduced the proliferation inhibition of vascular endothelial cells by Kringle 5 and had a blocking effect on the function of Kringle 5 in inhibiting migration and promoting apoptosis of vascular endothelial cells in vitro. This study put a new light on protein-aptamer binding mechanism and may provide insight into the treatment of ischemic diseases by target depletion of Kringle 5.

Identification of a broadly conserved family of enzymes that hydrolyze (p)ppApp

Bacteria produce a variety of nucleotide second messengers to adapt to their surroundings. Although chemically similar, the nucleotides guanosine penta- and tetraphosphate [(p)ppGpp] and adenosine penta- and tetraphosphate [(p)ppApp] have distinct functions in bacteria. (p)ppGpp mediates survival under nutrient-limiting conditions and its intracellular levels are regulated by synthetases and hydrolases belonging to the RelA-SpoT homolog (RSH) family of enzymes. By contrast, (p)ppApp is not known to be involved in nutrient stress responses and is synthesized by RSH-resembling toxins that inhibit the growth of bacterial cells. However, it remains unclear whether there exists a family of hydrolases that specifically act on (p)ppApp to reverse its toxic effects. Here, we present the structure and biochemical characterization of adenosine 3'-pyrophosphohydrolase 1 (Aph1), the founding member of a monofunctional (p)ppApp hydrolase family of enzymes. Our work reveals that Aph1 adopts a histidine-aspartate (HD)-domain fold characteristic of phosphohydrolase metalloenzymes and its activity mitigates the growth inhibitory effects of (p)ppApp-synthesizing toxins. Using an informatic approach, we identify over 2,000 putative (p)ppApp hydrolases that are widely distributed across bacterial phyla and found in diverse genomic contexts, and we demonstrate that 12 representative members hydrolyze ppApp. In addition, our in silico analyses reveal a unique molecular signature that is specific to (p)ppApp hydrolases, and we show that mutation of two residues within this signature broadens the specificity of Aph1 to promiscuously hydrolyze (p)ppGpp in vitro. Overall, our findings indicate that like (p)ppGpp hydrolases, (p)ppApp hydrolases are widespread in bacteria and may play important and underappreciated role(s) in bacterial physiology.

Using Macrocyclic G-Quadruplex Ligands to Decipher the Interactions Between Small Molecules and G-Quadruplex DNA

This study aims to deepen the knowledge of the current state of rational G4-ligand design through the design and synthesis of a novel set of compounds based on indoles, quinolines, and benzofurans and their comparisons with well-known G4-ligands. This resulted in novel synthetic methods and G4-ligands that bind and stabilize G4 DNA with high selectivity. Furthermore, the study corroborates previous studies on the design of G4-ligands and adds deeper explanations to why a) macrocycles offer advantages in terms of G4-binding and -selectivity, b) molecular pre-organization is of key importance in the development of strong novel binders, c) an electron-deficient aromatic core is essential to engage in strong arene-arene interactions with the G4-surface, and d) aliphatic amines can strengthen interactions indirectly through changing the arene electrostatic nature of the compound. Finally, fundamental physicochemical properties of selected G4-binders are evaluated, underscoring the complexity of aligning the properties required for efficient G4 binding and stabilization with feasible pharmacokinetic properties.

Influence of cation-π interactions to the structural stability of phycocyanin proteins: A computational study

The influences of cation-π interactions in phycocyanin proteins and their environmental preferences were analyzed. The number of interactions formed by arginine showed to be higher than those formed by the lysine in the cationic group, while histidine is comparatively higher than phenylalanine and N-terminal residue in the π group. Arg-Tyr and Arg-Phe interacting pairs are predominant among the various pairs analyzed. Cation-π interactions are distance-dependent and can be realized above a wider area above the π ring. We analyzed the energy contribution resulting from cation-π interactions using ab initio calculations. The energy contribution resulting from the most frequent cation-π interactions was in the lower range of strong hydrogen bonds. The results showed that, while most of their interaction energies lay ranged from - 2 to - 8 kcal/mol, those energies could be up to -12- 12 kcal/mol. Stabilization centers for these proteins showed that all residues found in cation-π interactions are important in locating one or more of such centers. In the cation-π interacting residues, 54% of the amino acid residues involved in these interactions might be conserved in phycocyanins. From this study, we infer that cation-π forming residues play an important role in the stability of the multiply commercially used phycocyanin proteins and could help structural biologists and medicinal chemists to design better and safer drugs.

Intramitochondrial co-assembly between ATP and nucleopeptides induces cancer cell apoptosis

Mitochondria are essential intracellular organelles involved in many cellular processes, especially adenosine triphosphate (ATP) production. Since cancer cells require high ATP levels for proliferation, ATP elimination can be a unique target for cancer growth inhibition. We describe a newly developed mitochondria-targeting nucleopeptide (MNP) that sequesters ATP by self-assembling with ATP inside mitochondria. MNP interacts strongly with ATP through electrostatic and hydrogen bonding interactions. MNP exhibits higher binding affinity for ATP (−637.5 kJ mol⁻¹) than for adenosine diphosphate (ADP) (−578.2 kJ mol⁻¹). To improve anticancer efficacy, the small-sized MNP/ADP complex formed large assemblies with ATP inside cancer cell mitochondria. ATP sequestration and formation of large assemblies of the MNP/ADP–ATP complex inside mitochondria caused physical stress by large structures and metabolic disorders in cancer cells, leading to apoptosis. This work illustrates a facile approach to developing cancer therapeutics that relies on molecular assemblies.

Mitochondria-targeting nucleopeptide (MNP) can sequester ATP by self-assembling with ATP. A small nanosized MNP/ADP complex forms a large assembly with ATP. Thus, intramitochondrial co-assembly causes stress by large structures and apoptosis.

Investigations on the role of cation-π interactions in active centers of superoxide dismutase

In this study, we have analysed the influence of cation-? interactions on stability and properties of superoxide dismutase (SOD) active centres. The number of interactions formed by arginine is higher than lysine in the cationic group, while histidine is comparatively higher in the ? group. The energy contribution resulting from most frequent cation-? interactions was in the lower range of strong hydrogen bonds. The cation-? interactions involving transition metal ions as cation have energy more negative than -418.4 kJ mol-1. Stabilization centres for these proteins showed that all residues involved in cation-? interactions were important in locating one or more of such centres. The majority of the residues involved in cation-p interactions were evolutionarily conserved and might have a significant contribution towards the stability of SOD proteins. The results presented in this work can be very useful for understanding the contribution of cation-? interactions to the stability of SOD active canters.

Inhibition of aquaporin-1 prevents myocardial remodeling by blocking the transmembrane transport of hydrogen peroxide

Pathological remodeling of the myocardium has long been known to involve oxidant signaling, but strategies using systemic antioxidants have generally failed to prevent it. We sought to identify key regulators of oxidant-mediated cardiac hypertrophy amenable to targeted pharmacological therapy. Specific isoforms of the aquaporin water channels have been implicated in oxidant sensing, but their role in heart muscle is unknown. RNA sequencing from human cardiac myocytes revealed that the archetypal AQP1 is a major isoform. AQP1 expression correlates with the severity of hypertrophic remodeling in patients with aortic stenosis. The AQP1 channel was detected at the plasma membrane of human and mouse cardiac myocytes from hypertrophic hearts, where it colocalized with NADPH oxidase-2 and caveolin-3. We show that hydrogen peroxide (H₂O₂), produced extracellularly, is necessary for the hypertrophic response of isolated cardiac myocytes and that AQP1 facilitates the transmembrane transport of H₂O₂ through its water pore, resulting in activation of oxidant-sensitive kinases in cardiac myocytes. Structural analysis of the amino acid residues lining the water pore of AQP1 supports its permeation by H₂O₂ Deletion of Aqp1 or selective blockade of the AQP1 intrasubunit pore inhibited H₂O₂ transport in mouse and human cells and rescued the myocyte hypertrophy in human induced pluripotent stem cell-derived engineered heart muscle. Treatment of mice with a clinically approved AQP1 inhibitor, Bacopaside, attenuated cardiac hypertrophy. We conclude that cardiac hypertrophy is mediated by the transmembrane transport of H₂O₂ by the water channel AQP1 and that inhibitors of AQP1 represent new possibilities for treating hypertrophic cardiomyopathies.

In Silico Analysis of the Molecular-Level Impact of SMPD1 Variants on Niemann-Pick Disease Severity

Sphingomyelin phosphodiesterase (SMPD1) is a key enzyme in the sphingolipid metabolism. Genetic SMPD1 variants have been related to the Niemann-Pick lysosomal storage disorder, which has different degrees of phenotypic severity ranging from severe symptomatology involving the central nervous system (type A) to milder ones (type B). They have also been linked to neurodegenerative disorders such as Parkinson and Alzheimer. In this paper, we leveraged structural, evolutionary and stability information on SMPD1 to predict and analyze the impact of variants at the molecular level. We developed the SMPD1-ZooM algorithm, which is able to predict with good accuracy whether variants cause Niemann-Pick disease and its phenotypic severity; the predictor is freely available for download. We performed a large-scale analysis of all possible SMPD1 variants, which led us to identify protein regions that are either robust or fragile with respect to amino acid variations, and show the importance of aromatic-involving interactions in SMPD1 function and stability. Our study also revealed a good correlation between SMPD1-ZooM scores and in vitro loss of SMPD1 activity. The understanding of the molecular effects of SMPD1 variants is of crucial importance to improve genetic screening of SMPD1-related disorders and to develop personalized treatments that restore SMPD1 functionality.

Structural Biology for the Molecular Insight between Aptamers and Target Proteins

Aptamers are promising therapeutic and diagnostic agents for various diseases due to their high affinity and specificity against target proteins. Structural determination in combination with multiple biochemical and biophysical methods could help to explore the interacting mechanism between aptamers and their targets. Regrettably, structural studies for aptamer–target interactions are still the bottleneck in this field, which are facing various difficulties. In this review, we first reviewed the methods for resolving structures of aptamer–protein complexes and for analyzing the interactions between aptamers and target proteins. We summarized the general features of the interacting nucleotides and residues involved in the interactions between aptamers and proteins. Challenges and perspectives in current methodologies were discussed. Approaches for determining the binding affinity between aptamers and target proteins as well as modification strategies for stabilizing the binding affinity of aptamers to target proteins were also reviewed. The review could help to understand how aptamers interact with their targets and how alterations such as chemical modifications in the structures affect the affinity and function of aptamers, which could facilitate the optimization and translation of aptamers-based theranostics.

Cross-Talk between Overlap Interactions in Biomolecules: A Case Study of the β-Turn Motif

Noncovalent interactions play a pivotal role in regulating protein conformation, stability and dynamics. Among the quantum mechanical (QM) overlap-based noncovalent interactions, n→π* is the best understood with studies ranging from small molecules to β-turns of model proteins such as GB1. However, these investigations do not explore the interplay between multiple overlap interactions in contributing to local structure and stability. In this work, we identify and characterize all noncovalent overlap interactions in the β-turn, an important secondary structural element that facilitates the folding of a polypeptide chain. Invoking a QM framework of natural bond orbitals, we demonstrate the role of several additional interactions such as n→σ* and π→π* that are energetically comparable to or larger than n→π*. We find that these interactions are sensitive to changes in the side chain of the residues in the β-turn of GB1, suggesting that the n→π* may not be the only component in dictating β-turn conformation and stability. Furthermore, a database search of n→σ* and π→π* in the PDB reveals that they are prevalent in most proteins and have significant interaction energies (∼1 kcal/mol). This indicates that all overlap interactions must be taken into account to obtain a comprehensive picture of their contributions to protein structure and energetics. Lastly, based on the extent of QM overlaps and interaction energies, we propose geometric criteria using which these additional interactions can be efficiently tracked in broad database searches.

The effect of cation-π interactions on the stability and electronic properties of anticancer drug Altretamine: a theoretical study

The present work utilizes density functional theory (DFT) calculations to study the influence of cation-π interactions on the electronic properties of the complexes formed by Altretamine [2,4,6-tris(dimethylamino)-1,3,5-triazine], an anticancer drug, with mono- and divalent (Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺ and Ca²⁺) metal cations. The structures were optimized with the M06-2X method and the 6-311++G(d,p) basis set in the gas phase and in solution. The theory of `Atoms in Molecules' (AIM) was applied to study the nature of the interactions by calculating the electron density ρ(r) and its Laplacian at the bond critical points. The charge-transfer process during complexation was evaluated using natural bond orbital (NBO) analysis. The results of DFT calculations demonstrate that the strongest/weakest interactions belong to Be²⁺/K⁺ complexes. There are good correlations between the achieved densities and the amounts of charge transfer with the interaction energies. Finally, the stability and reactivity of the cation-π interactions can be determined by quantum chemical computation based on the molecular orbital (MO) theory.

Dynamics of Ionic Interactions at Protein-Nucleic Acid Interfaces

Molecular association of proteins with nucleic acids is required for many biological processes essential to life. Electrostatic interactions via ion pairs (salt bridges) of nucleic acid phosphates and protein side chains are crucial for proteins to bind to DNA or RNA. Counterions around the macromolecules are also key constituents for the thermodynamics of protein–nucleic acid association. Until recently, there had been only a limited amount of experiment-based information about how ions and ionic moieties behave in biological macromolecular processes. In the past decade, there has been significant progress in quantitative experimental research on ionic interactions with nucleic acids and their complexes with proteins. The highly negatively charged surfaces of DNA and RNA electrostatically attract and condense cations, creating a zone called the ion atmosphere. Recent experimental studies were able to examine and validate theoretical models on ions and their mobility and interactions with macromolecules. The ionic interactions are highly dynamic. The counterions rapidly diffuse within the ion atmosphere. Some of the ions are released from the ion atmosphere when proteins bind to nucleic acids, balancing the charge via intermolecular ion pairs of positively charged side chains and negatively charged backbone phosphates. Previously, the release of counterions had been implicated indirectly by the salt-concentration dependence of the association constant.

Recently, direct detection of counterion release by NMR spectroscopy has become possible and enabled more accurate and quantitative analysis of the counterion release and its entropic impact on the thermodynamics of protein–nucleic acid association. Recent studies also revealed the dynamic nature of ion pairs of protein side chains and nucleic acid phosphates. These ion pairs undergo transitions between two major states. In one of the major states, the cation and the anion are in direct contact and form hydrogen bonds. In the other major state, the cation and the anion are separated by water. Transitions between these states rapidly occur on a picosecond to nanosecond time scale. When proteins interact with nucleic acids, interfacial arginine (Arg) and lysine (Lys) side chains exhibit considerably different behaviors. Arg side chains show a higher propensity to form rigid contacts with nucleotide bases, whereas Lys side chains tend to be more mobile at the molecular interfaces. The dynamic ionic interactions may facilitate adaptive molecular recognition and play both thermodynamic and kinetic roles in protein–nucleic acid interactions.

Racemic phosphorothioate as a tool for NMR investigations of protein-DNA complexes

A major driving force for protein-nucleic acid association is electrostatic interactions via ion pairs of the positively charged basic side chains and negatively charged phosphates. For a better understanding of how proteins scan DNA and recognize particular signatures, it is important to gain atomic-level insight into the behavior of basic side chains at the protein-DNA interfaces. NMR spectroscopy is a powerful tool for investigating the structural, dynamic, and kinetic aspects of protein-DNA interactions. However, resonance assignment of basic side-chain cationic moieties at the molecular interfaces remains to be a major challenge. Here, we propose a fast, robust, and inexpensive approach that greatly facilitates resonance assignment of interfacial moieties and also allows for kinetic measurements of protein translocation between two DNA duplexes. This approach utilizes site-specific incorporation of racemic phosphorothioate at the position of a phosphate that interacts with a protein side chain. This modification retains the electric charge of phosphate and therefore is mild, but causes significant chemical shift perturbations for the proximal protein side chains, which facilitates resonance assignment. Due to the racemic nature of the modification, two different chemical shifts are observed for the species with different diastereomers R_P and S_P of the incorporated phosphorothioate group. Kinetic information on the exchange of the protein molecule between R_P and S_P DNA duplexes can be obtained by ¹⁵N_z exchange spectroscopy. We demonstrate the applications of this approach to the Antennapedia homeodomain-DNA complex and the CREB1 basic leucine-zipper (bZIP)-DNA complex.

New insights into protein-DNA binding specificity from hydrogen bond based comparative study

Knowledge of protein-DNA binding specificity has important implications in understanding DNA metabolism, transcriptional regulation and developing therapeutic drugs. Previous studies demonstrated hydrogen bonds between amino acid side chains and DNA bases play major roles in specific protein-DNA interactions. In this paper, we investigated the roles of individual DNA strands and protein secondary structure types in specific protein-DNA recognition based on side chain-base hydrogen bonds. By comparing the contribution of each DNA strand to the overall binding specificity between DNA-binding proteins with different degrees of binding specificity, we found that highly specific DNA-binding proteins show balanced hydrogen bonding with each of the two DNA strands while multi-specific DNA binding proteins are generally biased towards one strand. Protein-base pair hydrogen bonds, in which both bases of a base pair are involved in forming hydrogen bonds with amino acid side chains, are more prevalent in the highly specific protein-DNA complexes than those in the multi-specific group. Amino acids involved in side chain-base hydrogen bonds favor strand and coil secondary structure types in highly specific DNA-binding proteins while multi-specific DNA-binding proteins prefer helices.

Improved Modeling of Cation-π and Anion-Ring Interactions Using the Drude Polarizable Empirical Force Field for Proteins

Cation-π interactions are noncovalent interactions between a π-electron system and a positively charged ion that are regarded as a strong noncovalent interaction and are ubiquitous in biological systems. Similarly, though less studied, anion-ring interactions are present in proteins along with in-plane interactions of anions with aromatic rings. As these interactions are between a polarizing ion and a polarizable π system, the accuracy of the treatment of these interactions in molecular dynamics (MD) simulations using additive force fields (FFs) may be limited. In the present work, to allow for a better description of ion-π interactions in proteins in the Drude-2013 protein polarizable FF, we systematically optimized the parameters for these interactions targeting model compound quantum mechanical (QM) interaction energies with atom pair-specific Lennard-Jones parameters along with virtual particles as selected ring centroids introduced to target the QM interaction energies and geometries. Subsequently, MD simulations were performed on a series of protein structures where ion-π pairs occur to evaluate the optimized parameters in the context of the Drude-2013 FF. The resulting FF leads to a significant improvement in reproducing the ion-π pair distances observed in experimental protein structures, as well as a smaller root-mean-square differences and fluctuations of the overall protein structures from experimental structures. Accordingly, the optimized Drude-2013 protein polarizable FF is suggested for use in MD simulations of proteins where cation-π and anion-ring interactions are critical. © 2019 Wiley Periodicals, Inc.

Interplay between Short-Range Attraction and Long-Range Repulsion Controls Reentrant Liquid Condensation of Ribonucleoprotein-RNA Complexes

In eukaryotic cells, ribonucleoproteins (RNPs) form mesoscale condensates by liquid-liquid phase separation that play essential roles in subcellular dynamic compartmentalization. The formation and dissolution of many RNP condensates are finely dependent on the RNA-to-RNP ratio, giving rise to a windowlike phase separation behavior. This is commonly referred to as reentrant liquid condensation (RLC). Here, using ribonucleoprotein-inspired polypeptides with low-complexity RNA-binding sequences as well as an archetypal disordered RNP, fused in sarcoma, as model systems, we investigate the molecular driving forces underlying this nonmonotonous phase transition. We show that an interplay between short-range cation-π attractions and long-range electrostatic forces governs the heterotypic RLC behavior of RNP-RNA complexes. Short-range attractions, which can be encoded by both polypeptide chain primary sequence and nucleic acid base sequence, control the two-phase coexistence regime, regulate material properties of polypeptide-RNA condensates, and oppose condensate reentrant dissolution. In the presence of excess RNA, a competition between short-range attraction and long-range electrostatic repulsion drives the formation of a colloidlike cluster phase. With increasing short-range attraction, the fluid dynamics of the cluster phase is arrested, leading to the formation of a colloidal gel. Our results reveal that phase behavior, supramolecular organization, and material states of RNP-RNA assemblies are controlled by a dynamic interplay between molecular interactions at different length scales.

Novel Deazaflavin Analogues Potently Inhibited Tyrosyl DNA Phosphodiesterase 2 (TDP2) and Strongly Sensitized Cancer Cells toward Treatment with Topoisomerase II (TOP2) Poison Etoposide

Topoisomerase II (TOP2) poisons as anticancer drugs work by trapping TOP2 cleavage complexes (TOP2cc) to generate DNA damage. Repair of such damage by tyrosyl DNA phosphodiesterase 2 (TDP2) could render cancer cells resistant to TOP2 poisons. Inhibiting TDP2, thus, represents an attractive mechanism-based chemosensitization approach. Currently known TDP2 inhibitors lack cellular potency and/or permeability. We report herein two novel subtypes of the deazaflavin TDP2 inhibitor core. By introducing an additional phenyl ring to the N-10 phenyl ring (subtype 11) or to the N-3 site of the deazaflavin scaffold (subtype 12), we have generated novel analogues with considerably improved biochemical potency and/or permeability. Importantly, many analogues of both subtypes, particularly compounds 11a, 11e, 12a, 12b, and 12h, exhibited much stronger cancer cell sensitizing effect than the best previous analogue 4a toward the treatment with etoposide, suggesting that these analogues could serve as effective cellular probes.

Tryptophan usage by Helicobacter pylori differs among strains

Because of its association with severe gastric pathologies, including gastric cancer, Helicobacter pylori has been subject of research for more than 30 years. Its capacity to adapt and survive in the human stomach can be attributed to its genetic flexibility. Its natural competence and its capacity to turn genes on and off allows H. pylori to adapt rapidly to the changing conditions of its host. Because of its genetic variability, it is difficult to establish the uniqueness of each strain obtained from a human host. The methods considered to-date to deliver the best result for differentiation of strains are Rapid Amplification of Polymorphic DNA (RAPD), Multilocus Sequence Typing (MLST) and Whole Genome Sequencing (WGS) analysis. While RAPD analysis is cost-effective, it requires a stable genome for its reliability. MLST and WGS are optimal for strain identification, however, they require analysis of data at the bioinformatics level. Using the StainFree method, which modifies tryptophan residues on proteins using 2, 2, 2, - trichloroethanol (TCE), we observed a strain specific pattern of tryptophan in 1D acrylamide gels. In order to establish the effectiveness of tryptophan fingerprinting for strain identification, we compared the graphic analysis of tryptophan-labelled bands in the gel images with MLST results. Based on this, we find that tryptophan banding patterns can be used as an alternative method for the differentiation of H. pylori strains. Furthermore, investigating the origin for these differences, we found that H. pylori strains alters the number and/or position of tryptophan present in several proteins at the genetic code level, with most exchanges taking place in membrane- and cation-binding proteins, which could be part of a novel response of H. pylori to host adaptation.

Computational analysis of the amino acid interactions that promote or decrease protein solubility

The solubility of globular proteins is a basic biophysical property that is usually a prerequisite for their functioning. In this study, we probed the solubility of globular proteins with the help of the statistical potential formalism, in view of objectifying the connection of solubility with structural and energetic properties and of the solubility-dependence of specific amino acid interactions. We started by setting up two independent datasets containing either soluble or aggregation-prone proteins with known structures. From these two datasets, we computed solubility-dependent distance potentials that are by construction biased towards the solubility of the proteins from which they are derived. Their analysis showed the clear preference of amino acid interactions such as Lys-containing salt bridges and aliphatic interactions to promote protein solubility, whereas others such as aromatic, His-π, cation-π, amino-π and anion-π interactions rather tend to reduce it. These results indicate that interactions involving delocalized π-electrons favor aggregation, unlike those involving no (or few) dispersion forces. Furthermore, using our potentials derived from either highly or weakly soluble proteins to compute protein folding free energies, we found that the difference between these two energies correlates better with solubility than other properties analyzed before such as protein length, isoelectric point and aliphatic index. This is, to the best of our knowledge, the first comprehensive in silico study of the impact of residue-residue interactions on protein solubility properties.The results of this analysis provide new insights that will facilitate future rational protein design applications aimed at modulating the solubility of targeted proteins.

Influenza D virus M2 protein exhibits ion channel activity in Xenopus laevis oocytes

BACKGROUND

A new type of influenza virus, known as type D, has recently been identified in cattle and pigs. Influenza D virus infection in cattle is typically asymptomatic; however, its infection in swine can result in clinical disease. Swine can also be infected with all other types of influenza viruses, namely A, B, and C. Consequently, swine can serve as a "mixing vessel" for highly pathogenic influenza viruses, including those with zoonotic potential. Currently, the only antiviral drug available targets influenza M2 protein ion channel is not completely effective. Thus, it is necessary to develop an M2 ion channel blocker capable of suppressing the induction of resistance to the genetic shift. To provide a basis for developing novel ion channel-blocking compounds, we investigated the properties of influenza D virus M2 protein (DM2) as a drug target.

RESULTS

To test the ion channel activity of DM2, the DNA corresponding to DM2 with cMyc-tag conjugated to its carboxyl end was cloned into the shuttle vector pNCB1. The mRNA of the DM2-cMyc gene was synthesized and injected into Xenopus oocytes. The translation products of DM2-cMyc mRNA were confirmed by immunofluorescence and mass spectrometry analyses. The DM2-cMyc mRNA-injected oocytes were subjected to the two-electrode voltage-clamp (TEVC) method, and the induced inward current was observed. The midpoint (Vmid) values in Boltzmann modeling for oocytes injected with DM2-cMyc RNA or a buffer were -152 and -200 mV, respectively. Assuming the same expression level in the Xenopus oocytes, DM2 without tag and influenza C virus M2 protein (CM2) were subjected to the TEVC method. DM2 exhibited ion channel activity under the condition that CM2 ion channel activity was reproduced. The gating voltages represented by Vmid for CM2 and DM2 were -141 and -146 mV, respectively. The reversal potentials observed in ND96 for CM2 and DM2 were -21 and -22 mV, respectively. Compared with intact DM2, DM2 variants with mutation in the YxxxK motif, namely Y72A and K76A DM2, showed lower Vmid values while showing no change in reversal potential.

CONCLUSION

The M2 protein from newly isolated influenza D virus showed ion channel activity similar to that of CM2. The gating voltage was shown to be affected by the YxxxK motif and by the hydrophobicity and bulkiness of the carboxyl end of the molecule.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

A Comprehensive Analysis of Anion-Quadrupole Interactions in Protein Structures

The edgewise interactions of anions with phenylalanine (Phe) aromatic rings in proteins, known as anion-quadrupole interactions, have been well studied. However, the anion-quadrupole interactions of the tyrosine (Tyr) and tryptophan (Trp) rings have been less well studied, probably because these have been considered weaker than interactions of anions hydrogen bonded to Trp/Tyr side chains. Distinguishing such hydrogen bonding interactions, we comprehensively surveyed the edgewise interactions of certain anions (aspartate, glutamate, and phosphate) with Trp, Tyr, and Phe rings in high-resolution, nonredundant protein single chains and interfaces (protein-protein, DNA/RNA-protein, and membrane-protein). Trp/Tyr anion-quadrupole interactions are common, with Trp showing the highest propensity and average interaction energy for this type of interaction. The energy of an anion-quadrupole interaction (-15.0 to 0.0 kcal/mol, based on quantum mechanical calculations) depends not only on the interaction geometry but also on the ring atom. The phosphate anions at DNA/RNA-protein interfaces interact with aromatic residues with energies comparable to that of aspartate/glutamate anion-quadrupole interactions. At DNA-protein interfaces, the frequency of aromatic ring participation in anion-quadrupole interactions is comparable to that of positive charge participation in salt bridges, suggesting an underappreciated role for anion-quadrupole interactions at DNA-protein (or membrane-protein) interfaces. Although less frequent than salt bridges in single-chain proteins, we observed highly conserved anion-quadrupole interactions in the structures of remote homologues, and evolutionary covariance-based residue contact score predictions suggest that conserved anion-quadrupole interacting pairs, like salt bridges, contribute to polypeptide folding, stability, and recognition.

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).

Cation-π interactions in protein-ligand binding: theory and data-mining reveal different roles for lysine and arginine

The interactions of neutral aromatic ligands with cationic arginine, histidine and lysine amino acid residues have been studied with ab initio calculations, symmetry adapted perturbation theory (SAPT), and a systematic meta-analysis of X-ray structures.

We have studied the cation–π interactions of neutral aromatic ligands with the cationic amino acid residues arginine, histidine and lysine using ab initio calculations, symmetry adapted perturbation theory (SAPT), and a systematic meta-analysis of all available Protein Data Bank (PDB) X-ray structures. Quantum chemical potential energy surfaces (PES) for these interactions were obtained at the DLPNO-CCSD(T) level of theory and compared against the empirical distribution of 2012 unique protein–ligand cation–π interactions found in X-ray crystal structures. We created a workflow to extract these structures from the PDB, filtering by interaction type and residue pK_a. The gas phase cation–π interaction of lysine is the strongest by more than 10 kcal mol^–1, but the empirical distribution of 582 X-ray structures lies away from the minimum on the interaction PES. In contrast, 1381 structures involving arginine match the underlying calculated PES with good agreement. SAPT analysis revealed that underlying differences in the balance of electrostatic and dispersion contributions are responsible for this behavior in the context of the protein environment. The lysine–arene interaction, dominated by electrostatics, is greatly weakened by a surrounding dielectric medium and causes it to become essentially negligible in strength and without a well-defined equilibrium separation. The arginine–arene interaction involves a near equal mix of dispersion and electrostatic attraction, which is weakened to a much smaller degree by the surrounding medium. Our results account for the paucity of cation–π interactions involving lysine, even though this is a more common residue than arginine. Aromatic ligands are most likely to interact with cationic arginine residues as this interaction is stronger than for lysine in higher polarity surroundings.

Structure of the Guanidine III Riboswitch

Riboswitches are structural elements found in mRNA molecules that couple small-molecule binding to regulation of gene expression, usually by controlling transcription or translation. We have determined high-resolution crystal structures of the ykkC guanidine III riboswitch from Thermobifida fusca. The riboswitch forms a classic H-type pseudoknot that includes a triple helix that is continuous with a central core of conserved nucleotides. These form a left-handed helical ramp of inter-nucleotide interactions, generating the guanidinium cation binding site. The ligand is hydrogen bonded to the Hoogsteen edges of two guanine bases. The binding pocket has a side opening that can accommodate a small side chain, shown by structures with bound methylguanidine, aminoguanidine, ethylguanidine, and agmatine. Comparison of the new structure with those of the guanidine I and II riboswitches reveals that evolution generated three different structural solutions for guanidine binding and subsequent gene regulation, although with some common elements.

Structure-based prediction of transcription factor binding specificity using an integrative energy function

Transcription factors (TFs) regulate gene expression through binding to specific target DNA sites. Accurate annotation of transcription factor binding sites (TFBSs) at genome scale represents an essential step toward our understanding of gene regulation networks. In this article, we present a structure-based method for computational prediction of TFBSs using a novel, integrative energy (IE) function. The new energy function combines a multibody (MB) knowledge-based potential and two atomic energy terms (hydrogen bond and π interaction) that might not be accurately captured by the knowledge-based potential owing to the mean force nature and low count problem. We applied the new energy function to the TFBS prediction using a non-redundant dataset that consists of TFs from 12 different families. Our results show that the new IE function improves the prediction accuracy over the knowledge-based, statistical potentials, especially for homeodomain TFs, the second largest TF family in mammals.

CONTACT

jguo4@uncc.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

An efficient algorithm for improving structure-based prediction of transcription factor binding sites

Background

Gene expression is regulated by transcription factors binding to specific target DNA sites. Understanding how and where transcription factors bind at genome scale represents an essential step toward our understanding of gene regulation networks. Previously we developed a structure-based method for prediction of transcription factor binding sites using an integrative energy function that combines a knowledge-based multibody potential and two atomic energy terms. While the method performs well, it is not computationally efficient due to the exponential increase in the number of binding sequences to be evaluated for longer binding sites. In this paper, we present an efficient pentamer algorithm by splitting DNA binding sequences into overlapping fragments along with a simplified integrative energy function for transcription factor binding site prediction.

Results

A DNA binding sequence is split into overlapping pentamers (5 base pairs) for calculating transcription factor-pentamer interaction energy. To combine the results from overlapping pentamer scores, we developed two methods, Kmer-Sum and PWM (Position Weight Matrix) stacking, for full-length binding motif prediction. Our results show that both Kmer-Sum and PWM stacking in the new pentamer approach along with a simplified integrative energy function improved transcription factor binding site prediction accuracy and dramatically reduced computation time, especially for longer binding sites.

Conclusion

Our new fragment-based pentamer algorithm and simplified energy function improve both efficiency and accuracy. To our knowledge, this is the first fragment-based method for structure-based transcription factor binding sites prediction.

Electronic supplementary material

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

The Structure of the Guanidine-II Riboswitch

The guanidine-II (mini-ykkC) riboswitch is the smallest of the guanidine-responsive riboswitches, comprising two stem loops of similar sequence. We have solved high-resolution crystal structures of both stem loops for the riboswitch from Gloeobacter violaceus. The stem loops have a strong propensity to dimerize by intimate loop-loop interaction. The dimerization creates specific binding pockets for two guanidine molecules, explaining their cooperative binding. Within the binding pockets the ligands are hydrogen bonded to a guanine at O6 and N7, and to successive backbone phosphates. Additionally they are each stacked upon a guanine nucleobase. One side of the pocket has an opening to the solvent, slightly lowering the specificity of ligand binding, and structures with bound methylguanidine, aminoguanidine, and agmatine show how this is possible.

Statistical analysis of structural determinants for protein-DNA-binding specificity

DNA-binding proteins play critical roles in biological processes including gene expression, DNA packaging and DNA repair. They bind to DNA target sequences with different degrees of binding specificity, ranging from highly specific (HS) to nonspecific (NS). Alterations of DNA-binding specificity, due to either genetic variation or somatic mutations, can lead to various diseases. In this study, a comparative analysis of protein-DNA complex structures was carried out to investigate the structural features that contribute to binding specificity. Protein-DNA complexes were grouped into three general classes based on degrees of binding specificity: HS, multispecific (MS), and NS. Our results show a clear trend of structural features among the three classes, including amino acid binding propensities, simple and complex hydrogen bonds, major/minor groove and base contacts, and DNA shape. We found that aspartate is enriched in HS DNA binding proteins and predominately binds to a cytosine through a single hydrogen bond or two consecutive cytosines through bidentate hydrogen bonds. Aromatic residues, histidine and tyrosine, are highly enriched in the HS and MS groups and may contribute to specific binding through different mechanisms. To further investigate the role of protein flexibility in specific protein-DNA recognition, we analyzed the conformational changes between the bound and unbound states of DNA-binding proteins and structural variations. The results indicate that HS and MS DNA-binding domains have larger conformational changes upon DNA-binding and larger degree of flexibility in both bound and unbound states. Proteins 2016; 84:1147-1161. © 2016 Wiley Periodicals, Inc.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

Molecular modeling and docking of small molecule inhibitors against NEK2

Aberrant expression of NEK2 (NIMA-related kinase 2) is indicated in a wide variety of human cancers. NEK2 is highly correlated to multi drug resistance by activating drug efflux activity. Identification of new small molecule inhibitors targeted against NEK2 therefore, facilitates to increase drug sensitivity of cancer cells, by stabilizing drug influx and minimizes the dose of therapeutic drug. Our work investigates to screen for optimal small molecule inhibitors against NEK2. In this study, we used a computational approach by modeling NEK2 protein using I-TASSER (Iterative Threading ASSEmbly Refinement) software. The modeled structure was subjected to protein preparation wizard; to add hydrogens and to optimize the protonation states of His, Gln and Asn residues. Active site of the modeled protein was identified using SiteMap tool of Schrodinger package. We further carried out docking studies by means of Glide, with various ligands downloaded from EDULISS database. Based on glide score, potential ligands were screened and their interaction with NEK2 was identified. The best hits were further screened for Lipinski's rule for drug-likeliness, bioactivity scoring and ADME properties. Thus, we report two (didemethylchlorpromazine and 2-[5-fluoro-1Hindol- 3-yl] propan-1-amine) compounds that have successfully satisfied all in silico parameters, necessitating further in vitro and in vivo studies.

A New Extrudable Form of Hypromellose: AFFINISOL™ HPMC HME

Hypromellose is a hydrophilic polymer widely used in immediate- and modified-release oral pharmaceutical dosage forms. However, currently available grades of hypromellose are difficult, if not impossible, to process by hot melt extrusion (HME) because of their high glass transition temperature, high melt viscosity, and low degradation temperature. To overcome these challenges, a modified grade of hypromellose, AFFINISOL™ HPMC HME, was recently introduced. It has a significantly lower glass transition temperature and melt viscosity as compared to other available grades of hypromellose. The objective of this paper is to assess the extrudability and performance of AFFINISOL™ HPMC HME (100LV and 4M) as compared to other widely used polymers in HME, including HPMC 2910 100cP (the currently available hypromellose), Soluplus®, Kollidon® VA 64, and EUDRAGIT® E PO. Formulations containing polymer and carbamazepine (CBZ) were extruded on a co-rotating 16-mm twin-screw extruder, and the effect of temperature, screw speed, and feed rate was investigated. The performance of the solid dispersions was evaluated based on Flory-Huggins modeling and characterized by differential scanning calorimetry (DSC), X-ray powder diffraction (XRD), Raman spectroscopy, Fourier-transform infrared (FTIR) spectroscopy, and dissolution. All formulations extruded well except for HPMC 2910 100cP, which resulted in over-torqueing the extruder (machine overloading because the motor cannot provide efficient energy to rotate the shaft). Among the HME extrudates, only the EUDRAGIT® E PO formulation was crystalline as confirmed by DSC, XRD, and Raman, which agreed with predictions from Flory-Huggins modeling. Dissolution testing was conducted under both sink and non-sink conditions. Sink dissolution testing in neutral media revealed that amorphous CBZ in the HME extrudates completely dissolved within 15 min, which was much more rapid than the time for complete dissolution of bulk CBZ (60 min) and EUDRAGIT® E PO solid dispersion (more than 6 h). Non-sink dissolution in acidic media testing revealed that only CBZ contained in the AFFINISOL™ HPMC HME, and EUDRAGIT® E PO solid dispersions rapidly supersaturated after 15 min, reaching a twofold drug concentration compared to the CBZ equilibrium solubility. In summary, AFFINISOL™ HPMC HME 100LV and AFFINISOL™ HPMC HME 4M are useful in the pharmaceutical HME process to increase wetting and dissolution properties of poorly water-soluble drugs like CBZ.

Cation–π interactions in CREBBP bromodomain inhibition: an electrostatic model for small-molecule binding affinity and selectivity

A new model is presented to explain and predict binding affinity of aromatic and heteroaromatic ligands for the CREBBP bromodomain based on cation–π interaction strength.

Histidine-Aromatic Interactions in Proteins and Protein-Ligand Complexes: Quantum Chemical Study of X-ray and Model Structures

His-aromatic complexes, with the His located above the aromatic plane, are stabilized by π-π, δ(+)-π and/or cation-π interactions according to whether the His is neutral or protonated and the partners are in stacked or T-shape conformations. Here we attempt to probe the relative strength of these interactions as a function of the geometry and protonation state, in gas phase, in water and protein-like environments (acetone, THF and CCl4), by means of quantum chemistry calculations performed up to second order of the Møller-Plesset pertubation theory. Two sets of conformations are considered for that purpose. The first set contains 89 interactions between His and Phe, Tyr, Trp, or Ade, observed in X-ray structures of proteins and protein-ligand complexes. The second set contains model structures obtained by moving an imidazolium/imidazole moiety above a benzene ring or an adenine moiety. We found that the protonated complexes are much more stable than the neutral ones in gas phase. This higher stability is due to the electrostatic contributions, the electron correlation contributions being equally important in the two forms. Thus, π-π and δ(+)-π interactions present essentially favorable electron correlation energy terms, whereas cation-π interactions feature in addition favorable electrostatic energies. The protonated complexes remain more stable than the neutral ones in protein-like environments, but the difference is drastically reduced. Furthermore, the T-shape conformation is undoubtedly more favorable than the stacked one in gas phase. This advantage decreases in the solvents, and the stacked conformation becomes even slightly more favorable in water. The frequent occurrence of His-aromatic interactions in catalytic sites, at protein-DNA or protein-ligand interfaces and in 3D domain swapping proteins emphasize their importance in biological processes.

Modeling Aromatic Liquids: Toluene, Phenol, and Pyridine

Aromatic groups are now acknowledged to play an important role in many systems of interest. However, existing molecular mechanics methods provide a poor representation of these groups. In a previous paper, we have shown that the molecular mechanics treatment of benzene can be improved by the incorporation of an explicit representation of the aromatic π electrons. Here, we develop this concept further, developing charge-separation models for toluene, phenol, and pyridine. Monte Carlo simulations are used to parametrize the models, via the reproduction of experimental thermodynamic data, and our models are shown to outperform an existing atom-centered model. The models are then used to make predictions about the structures of the liquids at the molecular level and are tested further through their application to the modeling of gas-phase dimers and cation-π interactions.