Amino/aromatic interactions in proteins: is the evidence stacked against hydrogen bonding?

Aromatic and arginine content drives multiphasic condensation of protein-RNA mixtures

Multiphasic architectures are found ubiquitously in biomolecular condensates and are thought to have important implications for the organization of multiple chemical reactions within the same compartment. Many of these multiphasic condensates contain RNA in addition to proteins. Here, we investigate the importance of different interactions in multiphasic condensates comprising two different proteins and RNA using computer simulations with a residue-resolution coarse-grained model of proteins and RNA. We find that in multilayered condensates containing RNA in both phases, protein-RNA interactions dominate, with aromatic residues and arginine forming the key stabilizing interactions. The total aromatic and arginine content of the two proteins must be appreciably different for distinct phases to form, and we show that this difference increases as the system is driven toward greater multiphasicity. Using the trends observed in the different interaction energies of this system, we demonstrate that we can also construct multilayered condensates with RNA preferentially concentrated in one phase. The "rules" identified can thus enable the design of synthetic multiphasic condensates to facilitate further study of their organization and function.

Risk mitigation strategy and mechanism analysis of neonicotinoid pesticides on earthworms

Neonicotinoid insecticides (NNIs) are a new class of widely used insecticides with certain risks to non-target organisms, like earthworms. The gray correlation method was used to calculate the comprehensive risk effect value of acute toxicity (LC50) and bioaccumulation (logKow) of NNIs on earthworms. A comprehensive effects three-dimensional quantitative structure-activity relationship (3D-QSAR) model was constructed, using NNIs molecular structures and the comprehensive effect value as the independent and dependent variables, respectively. One of the representatives guadipyr (GUA) was selected as the template molecule for the molecular design and modification. A total of 63 NNIs alternatives were designed with a reduced comprehensive value higher than 10%, and as high as 42%. After screening, 15 NNIs alternatives were screened with decreased acute toxicity to earthworms, bioaccumulation effects and improved functional property. The calculated primary acute risk quotient of earthworms shows that the designed NNIs alternatives have lower earthworm risks (reduction of 70.48-99.99%). Results also found that the electronic, geometric and topological parameters of NNIs are the key descriptors that affect NNIs alternatives' toxicity. The number of hydrophobic interaction amino acid residues in NNIs molecules also contributes to the acute toxicity and the bioaccumulation of NNIs alternatives on earthworms. This study aims to design and screen functionally improved and environmentally friendly NNIs alternatives that have low risk to earthworms and provide theoretical methods and new ideas for the risk control and development of pesticides represented by NNIs.

Systematic enhancement of protein crystallization efficiency by bulk lysine-to-arginine (KR) substitution

Structural genomics consortia established that protein crystallization is the primary obstacle to structure determination using x-ray crystallography. We previously demonstrated that crystallization propensity is systematically related to primary sequence, and we subsequently performed computational analyses showing that arginine is the most overrepresented amino acid in crystal-packing interfaces in the Protein Data Bank. Given the similar physicochemical characteristics of arginine and lysine, we hypothesized that multiple lysine-to-arginine (KR) substitutions should improve crystallization. To test this hypothesis, we developed software that ranks lysine sites in a target protein based on the redundancy-corrected KR substitution frequency in homologs. This software can be run interactively on the worldwide web at https://www.pxengineering.org/. We demonstrate that three unrelated single-domain proteins can tolerate 5-11 KR substitutions with at most minor destabilization, and, for two of these three proteins, the construct with the largest number of KR substitutions exhibits significantly enhanced crystallization propensity. This approach rapidly produced a 1.9 Å crystal structure of a human protein domain refractory to crystallization with its native sequence. Structures from Bulk KR-substituted domains show the engineered arginine residues frequently make hydrogen-bonds across crystal-packing interfaces. We thus demonstrate that Bulk KR substitution represents a rational and efficient method for probabilistic engineering of protein surface properties to improve crystallization.

Applying the bioisosterism strategy to obtain lead compounds against SARS-CoV-2 cysteine proteases: An in-silico approach

SARS-CoV-2 cysteine proteases are essential nonstructural proteins due to their role in the formation of the virus multiple enzyme replication-transcription complex. As a result, those functional proteins are extremely relevant targets in the development of a new drug candidate to fight COVID-19. Based on this fact and guided by the bioisosterism strategy, the present work has selected 126 out of 1050 ligands from DrugBank website. Subsequently, 831 chemical analogs containing bioisosteres, some of which became structurally simplified, were created using the MB-Isoster software, and molecular docking simulations were performed using AutoDock Vina. Finally, a study of physicochemical properties, along with pharmacokinetic profiles, was carried out through SwissADME and ADMETlab 2.0 platforms. The promising results obtained with the molecules encoded as DB00549_BI_005, DB04868_BI_003, DB11984_BI_002, DB12364_BI_006 and DB12805_BI_004 must be confirmed by molecular dynamics studies, followed by in vitro and in vivo empirical tests that ratify the advocated in-silico results.

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.

CIPDB: A biological structure databank for studying cation and π interactions

As major forces for modulating protein folding and molecular recognition, cation and π interactions are extensively identified in protein structures. They are even more competitive than hydrogen bonds in molecular recognition, thus, are vital in numerous biological processes. In this review, we introduce the methods for the identification and quantification of cation and π interactions, provide insights into the characteristics of cation and π interactions in the natural state, and reveal their biological function together with our developed database (Cation and π Interaction in Protein Data Bank; CIPDB; http://chemyang.ccnu.edu.cn/ccb/database/CIPDB). This review lays the foundation for the in-depth study of cation and π interactions and will guide the use of molecular design for drug discovery.

The biological role of charge distribution in linear antimicrobial peptides

INTRODUCTION

Antimicrobial peptides (AMP) have received particular attention due to their capacity to kill bacteria. Although much is known about them, peptides are currently being further researched. A large number of AMPs have been discovered, but only a few have been approved for topical use, due to their promiscuity and other challenges, which need to be overcome.

AREAS COVERED

AMPs are diverse in structure. Consequently, they have varied action mechanisms when targeting microorganisms or eukaryotic cells. Herein, the authors focus on linear peptides, particularly those that are alpha-helical structured, and examine how their charge distribution and hydrophobic amino acids could modulate their biological activity.

EXPERT OPINION

The world currently needs urgent solutions to the infective problems caused by resistant pathogens. In order to start the race for antimicrobial development from the charge distribution viewpoint, bioinformatic tools will be necessary. Currently, there is no software available that allows to discriminate charge distribution in AMPs and predicts the biological effects of this event. Furthermore, there is no software available that predicts the side-chain length of residues and its role in biological functions. More specialized software is necessary.

Piezoresistive Fibers with Large Working Factors for Strain Sensing Applications

Piezoresistive fibers with large working factors remain of great interest for strain sensing applications involving large strains, yet difficult to achieve. Here, we produced strain-sensitive fibers with large working factors by dip-coating nanocomposite piezoresistive inks on surface-modified polyether block amide (PEBA) fibers. Surface modification of neat PEBA fibers was carried out with polydopamine (PDA) while nanocomposite conductive inks consisted of styrene-ethylene-butylene-styrene (SEBS) elastomer and carbon black (CB). As such, the deposition of piezoresistive coatings was enabled through nonconventional hydrogen-bonding interactions. The resultant fibers demonstrated well-defined piezoresistive linear relationships, which increased with CB filler loading in SEBS. In addition, gauge factors decreased with increasing CB mass fractions from ∼15 to ∼7. Furthermore, we used the fatigue theory to predict the endurance limit (C_e) of our fibers toward resistance signal stability. Such a piezoresistive performance allowed us to explore the application of our fibers as strain sensors for monitoring the movement of finger joints.

Through-Space Stabilization of an Imidazolium Cation by Aromatic Rings

Imidazole-based compounds are widely found in natural products, synthetic molecules, and biomolecules. Noncovalent interactions between the imidazole ring and other functional groups play an important role in determining the function of diverse molecules. However, there is a limited understanding of the underlying noncovalent interactions between imidazoles and aromatic systems. In this work, we report physical-organic chemistry studies on 2-(2,6-diarylphenyl)-1H-imidazoles and their protonated forms to investigate the noncovalent interactions between the central imidazole ring and two flanking aromatic rings possessing substituents at the para/meta position. Hammett analysis revealed that pK_a values and proton affinities correlate well with Hammett σ values of para-substituents at the flanking rings. Additional quantitative Kohn-Sham molecular orbital and energy decomposition analyses reveal that through-space π-π interactions and NH-π interactions contribute to the intramolecular stabilization of the imidazolium cation. The results are important because they clearly demonstrate that the imidazolium cation forms energetically favorable noncovalent interactions with aromatic rings via the through-space effect, a knowledge that can be used in rational drug and catalyst design.

3D Interaction Homology: Hydropathic Analyses of the "π-Cation" and "π-π" Interaction Motifs in Phenylalanine, Tyrosine, and Tryptophan Residues

Three-dimensional (3D) maps of the hydropathic environments of protein amino acid residues are information-rich descriptors of preferred conformations, interaction types and energetics, and solvent accessibility. The interactions made by each residue are the primary factor for rotamer selection and the secondary, tertiary, and even quaternary protein structure. Our evolving basis set of environmental data for each residue type can be used to understand the protein structure. This work focuses on the aromatic residues phenylalanine, tyrosine, and tryptophan and their structural roles. We calculated and analyzed side chain-to-environment 3D maps for over 70,000 residues of these three types that reveal, with respect to hydrophobic and polar interactions, the environment around each. After binning with backbone ϕ/ψ and side chain χ₁, we clustered each bin by 3D similarities between map-map pairs. For each of the three residue types, four bins were examined in detail: one in the β-pleat, two in the right-hand α-helix, and one in the left-hand α-helix regions of the Ramachandran plot. For high degrees of side chain burial, encapsulation of the side chain by hydrophobic interactions is ubiquitous. The more solvent-exposed side chains are more likely to be involved in polar interactions with their environments. Evidence for π-π interactions was observed in about half of the residues surveyed [phenylalanine (PHE): 53.3%, tyrosine (TYR): 34.1%, and tryptophan (TRP): 55.7%], but on an energy basis, this contributed to only ∼4% of the total. Evidence for π-cation interactions was observed in 14.1% of PHE, 8.3% of TYR, and 26.8% of TRP residues, but on an energy basis, this contributed to only ∼1%. This recognition of even these subtle interactions in the 3D hydropathic environment maps is key support for our interaction homology paradigm of protein structure elucidation and possibly prediction.

Structure of Phycobilisomes

Phycobilisomes (PBSs) are extremely large chromophore-protein complexes on the stromal side of the thylakoid membrane in cyanobacteria and red algae. The main function of PBSs is light harvesting, and they serve as antennas and transfer the absorbed energy to the reaction centers of two photosynthetic systems (photosystems I and II). PBSs are composed of phycobiliproteins and linker proteins. How phycobiliproteins and linkers are organized in PBSs and how light energy is efficiently harvested and transferred in PBSs are the fundamental questions in the study of photosynthesis. In this review, the structures of the red algae Griffithsia pacifica and Porphyridium purpureum are discussed in detail, along with the functions of linker proteins in phycobiliprotein assembly and in fine-tuning the energy state of chromophores.

Cation-π Interactions and Their Contribution to Mussel Underwater Adhesion Studied Using a Surface Forces Apparatus: A Mini-Review

Mussel underwater adhesion is a model phenomenon important for the understanding of broader biological adhesion and the development of biomimetic wet adhesives. The catechol moiety of 3,4-dihydroxyphenyl-l-alanine (DOPA) is known to be actively involved in the mechanism of mussel underwater adhesion; however, other underwater adhesion mechanisms are also crucial. The surface forces apparatus (SFA) has often been used to explore the contributions of other mechanisms to mussel underwater adhesion; e.g., recent SFA-based nanomechanical studies have revealed that cation-π interactions, one of the strongest intermolecular interactions in water, are the pivotal interactions of adhesive proteins involved in underwater mussel adhesion. This mini-review surveys recent research on cation-π interactions and their contributions to strong mussel underwater adhesion, shedding light on some biological processes and facilitating the development of biomedical adhesives.

DFT investigation of role of N - H⋯O and N - H⋯π interactions in the stabilization of the hydrogen bonded complexes of anisole with aromatic amines

Theoretical investigations have been performed on hydrogen (H-) bonded complexes of two aromatic amines with anisole to investigate the effect of the methyl substituent on N – H⋯O and N – –H⋯π interactions. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses were done to elucidate the nature of H- bonding. In 1:1 complexes, the total interaction energy of N-methylaniline complex is higher than that of aniline complex. The existence of bond critical point between N–H of amine and oxygen of anisole confirms weak hydrogen bonding. The energy decomposition analysis showed the role of CT in stabilizing complexes.

A homozygous variant in RRM2B is associated with severe metabolic acidosis and early neonatal death

RRM2B encodes the crucial p53-inducible ribonucleotide reductase small subunit 2 homolog (p53R2), which is required for DNA synthesis throughout the cell cycle. Mutations in this gene have been associated with a lethal mitochondrial depletion syndrome. Here we present the case of an infant with a novel homozygous p.Asn221Ser mutation in RRM2B who developed hypotonia, failure to thrive, sensorineural hearing loss, and severe metabolic lactic acidosis, ultimately progressing to death at 3 months of age. Through molecular modeling using the X-ray crystal structure of p53R2, we demonstrate that this mutation likely causes disruption of a highly conserved helix region of the protein by altering intramolecular interactions. This report expands our knowledge of potential pathogenic RRM2B mutations as well as our understanding of the molecular function of p53R2 and its role in the pathogenesis of mitochondrial DNA depletion.

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.

Stacking Interactions of Heterocyclic Drug Fragments with Protein Amide Backbones

Stacking interactions can be important enthalpic contributors to drug binding. Among the less well-studied stacking interactions are those occurring between an arene and the π-face of an amide group. Given the ubiquity of heterocycles in drugs, combined with the abundance of amides in the protein backbone, optimizing these noncovalent interactions can provide a potential route to enhanced drug binding. Previously, Diederich et al. (ChemMedChem 2013, 8, 397-404) studied stacked dimers of a model amide with a set of 18 heterocycles, showing that computed interaction energies correlate with the dipole moments of the heterocycles and providing guidelines for the optimization of these interactions. We considered stacked dimers of the same model amide with a larger set of 28 heterocycles common in pharmaceuticals, by using more robust ab initio methods. While the overall trends in these new data corroborate many of the results of Diederich et al., these data provide a more refined view of the nature of amide stacking interactions. We present a robust scoring function for amide stacking interaction energies based on the molecular dipole moment and strength of the electric field above the arene.

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

Total Synthesis and Conformational Study of Callyaerin A: Anti-Tubercular Cyclic Peptide Bearing a Rare Rigidifying (Z)-2,3- Diaminoacrylamide Moiety

The first synthesis of the anti-TB cyclic peptide callyaerin A (1), containing a rare (Z)-2,3-diaminoacrylamide bridging motif, is reported. Fmoc-formylglycine-diethylacetal was used as a masked equivalent of formylglycine in the synthesis of the linear precursor to 1. Intramolecular cyclization between the formylglycine residue and the N-terminal amine in the linear peptide precursor afforded the macrocyclic natural product 1. Synthetic 1 possessed potent anti-TB activity (MIC₁₀₀ =32 μm) while its all-amide congener was inactive. Variable-temperature NMR studies of both the natural product and its all-amide analogue revealed the extraordinary rigidity imposed by this diaminoacrylamide unit on peptide conformation. The work reported herein pinpoints the intrinsic role that the (Z)-2,3-diaminoacrylamide moiety confers on peptide bioactivity.

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.

The structures and properties of anionic tryptophan complexes

IRMPD spectroscopy and electronic structure calculations are employed to identify π–π interactions in ionic tryptophan clusters.

A systematic analysis of atomic protein-ligand interactions in the PDB

As the protein databank (PDB) recently passed the cap of 123 456 structures, it stands more than ever as an important resource not only to analyze structural features of specific biological systems, but also to study the prevalence of structural patterns observed in a large body of unrelated structures, that may reflect rules governing protein folding or molecular recognition. Here, we compiled a list of 11 016 unique structures of small-molecule ligands bound to proteins - 6444 of which have experimental binding affinity - representing 750 873 protein-ligand atomic interactions, and analyzed the frequency, geometry and impact of each interaction type. We find that hydrophobic interactions are generally enriched in high-efficiency ligands, but polar interactions are over-represented in fragment inhibitors. While most observations extracted from the PDB will be familiar to seasoned medicinal chemists, less expected findings, such as the high number of C-H···O hydrogen bonds or the relatively frequent amide-π stacking between the backbone amide of proteins and aromatic rings of ligands, uncover underused ligand design strategies.

A Paradigm for CH Bond Cleavage: Structural and Functional Aspects of Transition State Stabilization by Mandelate Racemase

Mandelate racemase (MR) from Pseudomonas putida catalyzes the Mg²⁺-dependent, 1,1-proton transfer reaction that racemizes (R)- and (S)-mandelate. MR shares a partial reaction (i.e., the metal ion-assisted, Brønsted base-catalyzed proton abstraction of the α-proton of carboxylic acid substrates) and structural features ((β/α)₇β-barrel and N-terminal α + β capping domains) with a vast group of homologous, yet functionally diverse, enzymes in the enolase superfamily. Mechanistic and structural studies have developed this enzyme into a paradigm for understanding how enzymes such as those of the enolase superfamily overcome kinetic and thermodynamic barriers to catalyze the abstraction of an α-proton from a carbon acid substrate with a relatively high pK_a value. Structural studies on MR bound to intermediate/transition state analogues have delineated those structural features that MR uses to stabilize transition states and enhance reaction rates of proton abstraction. Kinetic, site-directed mutagenesis, and structural studies have also revealed that the phenyl ring of the substrate migrates through the hydrophobic cavity within the active site during catalysis and that the Brønsted acid-base catalysts (Lys 166 and His 297) may be utilized as binding determinants for inhibitor recognition. In addition, structural studies on the adduct formed from the irreversible inhibition of MR by 3-hydroxypyruvate revealed that MR can form and deprotonate a Schiff-base with 3-hydroxypyruvate to yield an enol(ate)-aldehyde adduct, suggesting a possible evolutionary link between MR and the Schiff-base forming aldolases. As the archetype of the enolase superfamily, mechanistic and structural studies on MR will continue to enhance our understanding of enzyme catalysis and furnish insights into the evolution of enzyme function.

Role of 6-Mercaptopurine in the potential therapeutic targets DNA base pairs and G-quadruplex DNA: insights from quantum chemical and molecular dynamics simulations

The theoretical studies on DNA with the anticancer drug 6-Mercaptopurine (6-MP) are investigated using theoretical methods to shed light on drug designing. Among the DNA base pairs considered, 6-MP is stacked with GC with the highest interaction energy of -46.19 kcal/mol. Structural parameters revealed that structure of the DNA base pairs is deviated from the planarity of the equilibrium position due to the formation of hydrogen bonds and stacking interactions with 6-MP. These deviations are verified through the systematic comparison between X-H bond contraction and elongation and the associated blue shift and red shift values by both NBO analysis and vibrational analysis. Bent's rule is verified for the C-H bond contraction in the 6-MP interacted base pairs. The AIM results disclose that the higher values of electron density (ρ) and Laplacian of electron density (∇²ρ) indicate the increased overlap between the orbitals that represent the strong interaction and positive values of the total electron density show the closed-shell interaction. The relative sensitivity of the chemical shift values for the DNA base pairs with 6-MP is investigated to confirm the hydrogen bond strength. Molecular dynamics simulation studies of G-quadruplex DNA d(TGGGGT)₄ with 6-MP revealed that the incorporation of 6-MP appears to cause local distortions and destabilize the G-quadruplex DNA.

Interplay among Electrostatic, Dispersion, and Steric Interactions: Spectroscopy and Quantum Chemical Calculations of π-Hydrogen Bonded Complexes

π-Hydrogen bonding interactions are ubiquitous in both materials and biology. Despite their relatively weak nature, great progress has been made in their investigation by experimental and theoretical methods, but this becomes significantly more complicated when secondary intermolecular interactions are present. In this study, the effect of successive methyl substitution on the supramolecular structure and interaction energy of indole⋅⋅⋅methylated benzene (ind⋅⋅⋅n-mb, n=1-6) complexes is probed through a combination of supersonic jet experiments and benchmark-quality quantum chemical calculations. It is demonstrated that additional secondary interactions introduce a subtle interplay among electrostatic and dispersion forces, as well as steric repulsion, which fine-tunes the overall structural motif. Resonant two-photon ionization and IR-UV double-resonance spectroscopy techniques are used to probe jet-cooled ind⋅⋅⋅n-mb (n=2, 3, 6) complexes, with redshifting of the N-H IR stretching frequency showing that increasing the degree of methyl substitution increases the strength of the primary N-H⋅⋅⋅π interaction. Ab initio harmonic frequency and binding energy calculations confirm this trend for all six complexes. Electronic spectra of the three dimers are broad and structureless, with quantum chemical calculations revealing that this is likely to be due to multiple tilted conformations of each dimer possessing similar stabilization energies.

Antimicrobial and anti-inflammatory activities of three chensinin-1 peptides containing mutation of glycine and histidine residues

The natural peptide chensinin-1 doesnot exhibit its desired biological properties. In this study, the mutant MC1-1 was designed by replacing Gly in the chensinin-1 sequence with Trp. Mutants MC1-2 and MC1-3 were designed based on the MC1-1 sequence to investigate the specific role of His residues. The mutated peptides presented α-helicity in a membrane-mimetic environment and exhibited broad-spectrum antimicrobial activities; in contrast to Trp residues, His residues were dispensable for interacting with the cell membrane. The interactions between the mutant peptides and lipopolysaccharide (LPS) facilitated the ingestion of peptides by Gram-negative bacteria. The binding affinities of the peptides were similar, at approximately 10 μM, but ΔH for MC1-2 was −7.3 kcal.mol⁻¹, which was 6-9 folds higher than those of MC1-1 and MC1-3, probably due to the conformational changes. All mutant peptides demonstrated the ability to inhibit LPS-induced tumour-necrosis factor-α (TNF-α) and interleukin-6 (IL-6) release from murine RAW264.7 cells. In addition, the representative peptide MC1-1showed better inhibition of serum TNF-α and IL-6 levels compared to polymyxin B (PMB), a potent binder and neutralizer of LPS as positive control in LPS-challenged mice model. These data suggest that the mutant peptides could be promising molecules for development as chensinin-based therapeutic agents against sepsis.

Elucidating a chemical defense mechanism of Antarctic sponges: A computational study

In 2000, a novel secondary metabolite (erebusinone, Ereb) was isolated from the Antarctic sea sponge, Isodictya erinacea. The bioactivity of Ereb was investigated, and it was found to inhibit molting when fed to the arthropod species Orchomene plebs. Xanthurenic acid (XA) is a known endogenous molt regulator present in arthropods. Experimental studies have confirmed that XA inhibits molting by binding to either (or both) of two P450 enzymes (CYP315a1 or CYP314a1) that are responsible for the final two hydroxylations in the production of the molt-inducing hormone, 20-hydroxyecdysone (20E). The lack of crystal structures and biochemical assays for CYP315a1 or CYP314a1, has prevented further experimental exploration of XA and Ereb's molt inhibition mechanisms. Herein, a wide array of computational techniques - homology modeling, molecular dynamics simulations, binding site bioinformatics, flexible receptor-flexible ligand docking, and molecular mechanics-generalized Born surface area calculations - have been employed to elucidate the structure-function relationships between the aforementioned P450s and the two described small molecule inhibitors (Ereb and XA). Results indicate that Ereb likely targets CYP315a1 by interacting with a network of aromatic residues in the binding site, while XA may inhibit both CYP315a1 and CYP314a1 because of its aromatic, as well as charged nature.

Removal of antibiotics from water in the coexistence of suspended particles and natural organic matters using amino-acid-modified-chitosan flocculants: A combined experimental and theoretical study

Contamination of trace antibiotics is widely found in surface water sources. This work delineates removal of trace antibiotics (norfloxacin (NOR), sulfadiazine (SDZ) or tylosin (TYL)) from synthetic surface water by flocculation, in the coexistence of inorganic suspended particles (kaolin) and natural organic matter (humic acid, HA). To avoid extra pollution caused by petrochemical products-based modification reagents, environmental-friendly amino-acid-modified-chitosan flocculants, Ctrp and Ctyr, with different functional aromatic-rings structures were employed. Jar tests at various pHs exhibited that, Ctyr, owning phenol groups as electron donors, was favored for elimination of cationic NOR (∼50% removal; optimal pH: 6; optimal dosage: 4mg/L) and TYL (∼60% removal; optimal pH: 7; optimal dosage: 7.5mg/L), due to π-π electron donator-acceptor (EDA) effect and unconventional H-bonds. Differently, Ctrp with indole groups as electron acceptor had better removal rate (∼50%) of SDZ anions (electron donator). According to correlation analysis, the coexisted kaolin and HA played positive roles in antibiotics' removal. Detailed pairwise interactions in molecular level among different components were clarified by spectral analysis and theoretical calculations (density functional theory), which are important for both the structural design of new flocculants aiming at targeted contaminants and understanding the environmental behaviors of antibiotics in water.

Essential dynamics analysis captures the concerted motion of the integrin-binding site in jerdostatin, an RTS disintegrin

Disintegrins, small molecular weight proteins contained in the venom of vipers and rattlesnakes, are high-affinity and selectivity integrin antagonists. Disintegrins inhibitory epitope mainly consists in a tripeptide sequence localized in a mobile loop protruding from the protein core. RTS and/or KTS tripeptide characterizes the most recently discovered group of disintegrins that selectively block α1β1 integrin receptor. A NMR study dedicated to structure and dynamics properties of jerdostatin, an RTS disintegrin, demonstrated that the substitution of the native RTS with KTS motif impaired flexibility and inhibitory activity of the molecule. Here we add atomic details to the experimental profiles of jerdostatin and its R24K mutant by analyzing the dynamics behavior of the molecules through computational methods. For jerdostatin wild type, molecular dynamics simulations and essential dynamics analyses showed that Y31 residue acts as hinge element in the concerted motions involving the active loop and the C-terminal tail. R24 side chain ability to engage both cation-π and H-bond interactions with Y31 residue was found crucial for that breathing mechanism. Less significant loop-tail concerted motions were observed for the R24K mutant. The description at atomic resolution of jerdostatin dynamics is useful for decoding the influence of specific residues on disintegrin functional properties.

Activation of Benzyl Aryl Carbonates: The Role of Cation-π Interactions

Benzyl aryl carbonates can react with a nucleophile to yield an activated electrophile and an aryloxide anion. Previously, we had utilized this in the synthesis of α-nitro esters from nitroalkanes. To further understand the process of activation of these carbonates by nucleophiles, we have performed kinetic studies on the hydrolysis of carbonates using nucleophiles. Rate constants for the hydrolysis were obtained under pseudo-first-order conditions with DABCO as the nucleophile. A comparison of rate constant for hydrolysis of isobutyl phenyl carbonate with benzyl phenyl carbonate shows that the presence of benzyl group results in a 16-fold acceleration of hydrolysis rate. This indicates that the transition state for activation of carbonate is stabilized by cation-π interactions. A comparison of the rate constant for various aromatic rings indicates that electron-donating substituents on the benzyl groups accelerate the rate of hydrolysis. Studies were also carried out with DMAP as nucleophile and the results are presented. Our studies show that stable carbonates can be activated using nucleophiles. Activated acyl groups generated from acid anhydrides have been used in several enantioselective reactions. Our studies show that carbonates can be stable alternatives to acid anhydrides.

The n → π* interaction: a rapidly emerging non-covalent interaction

This perspective describes the current status of a recently discovered non-covalent interaction named as the n → π* interaction, which is very weak and counterintuitive in nature. In this review, we have provided a brief overview of the widespread presence of this interaction in biomacromolecules, small biomolecules and materials, as well as the physical nature of this interaction explored using various experimental and theoretical techniques. It has been found that this interaction is equally important to other non-covalent interactions for the stability and specific structures of biomolecules and materials. An in-depth understanding of this interaction can help in designing more efficient functional materials as well as drugs. The review also provides a future outlook in terms of exploring the detailed functional role of this interaction in biological processes and its direct spectroscopic evidence, which other commonly known non-covalent interactions (conventional hydrogen bonding, π-hydrogen bonding, π-stacking, etc.) have.

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.

UV Action Spectroscopy of Gas-Phase Peptide Radicals

UV photodissociation (UVPD) action spectroscopy is reported to provide a sensitive tool for the detection of radical sites in gas-phase peptide ions. UVPD action spectra of peptide cation radicals of the z-type generated by electron-transfer dissociation point to the presence of multiple structures formed as a result of spontaneous isomerizations by hydrogen atom migration. N-terminal Cα radicals are identified as the dominant components, but the content of isomers differing in the radical defect position in the backbone or side chain depends on the nature of the aromatic residue with phenylalanine being more prone to isomerization than tryptophan. These results illustrate that spontaneous hydrogen atom migrations can occur in peptide cation-radicals upon electron-transfer dissociation.

The SPL7013 dendrimer destabilizes the HIV-1 gp120-CD4 complex

The poly (l-lysine)-based SPL7013 dendrimer with naphthalene disulphonate surface groups blocks the entry of HIV-1 into target cells and is in clinical trials for development as a topical microbicide. Its mechanism of action against R5 HIV-1, the HIV-1 variant implicated in transmission across individuals, remains poorly understood. Using docking and fully atomistic MD simulations, we find that SPL7013 binds tightly to R5 gp120 in the gp120-CD4 complex but weakly to gp120 alone. Further, the binding, although to multiple regions of gp120, does not occlude the CD4 binding site on gp120, suggesting that SPL7013 does not prevent the binding of R5 gp120 to CD4. Using MD simulations to compute binding energies of several docked structures, we find that SPL7013 binding to gp120 significantly weakens the gp120-CD4 complex. Finally, we use steered molecular dynamics (SMD) to study the kinetics of the dissociation of the gp120-CD4 complex in the absence of the dendrimer and with the dendrimer bound in each of the several stable configurations to gp120. We find that SPL7013 significantly lowers the force required to rupture the gp120-CD4 complex and accelerates its dissociation. Taken together, our findings suggest that SPL7013 compromises the stability of the R5 gp120-CD4 complex, potentially preventing the accrual of the requisite number of gp120-CD4 complexes across the virus-cell interface, thereby blocking virus entry.

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.

Contribution of cation-π interactions to the stability of Sm/LSm oligomeric assemblies

In this work, we have analyzed the influence of cation-π interactions to the stability of Sm/LSm assemblies and their environmental preferences. The number of interactions formed by arginine is higher than lysine in the cationic group, while histidine is comparatively higher than phenylalanine and tyrosine in the π group. Arg-Tyr interactions are predominant among the various pairs analyzed. The furcation level of multiple cation-π interactions is much higher than that of single cation-π interactions in Sm/LSm interfaces. We have found hot spot residues forming cation-π interactions, and hot spot composition is similar for all aromatic residues. The Arg-Phe pair has the strongest interaction energy of -8.81 kcal mol(-1) among all the possible pairs of amino acids. The extent of burial of the residue side-chain correlates with the ΔΔG of binding for residues in the core and also for hot spot residues cation-π bonded across the interface. Secondary structure of the cation-π residues shows that Arg and Lys preferred to be in strand. Among the π residues, His prefers to be in helix, Phe prefers to be in turn, and Tyr prefers to be in strand. Stabilization centers for these proteins showed that all the five residues found in cation-π interactions are important in locating one or more of such centers. More than 50 % of the cation-π interacting residues are highly conserved. It is likely that the cation-π interactions contribute significantly to the overall stability of Sm/LSm proteins.

M2 proton channel: toward a model of a primitive proton pump

Transmembrane proton transfer was essential to early cellular systems in order to transduce energy for metabolic functions. The reliable, efficient and controlled generation of proton gradients became possible only with the emergence of active proton pumps. On the basis of features shared by most modern proton pumps we identify the essential mechanistic steps in active proton transport. Further, we discuss the mechanism of action of a small, transmembrane M2 proton channel from influenza A virus as a model for proton transport in protocells. The M2 channel is a 94-residue long, α-helical tetramer that is activated at low pH and exhibits high selectivity and directionality. A shorter construct, built of transmembrane fragments that are only 24 amino acids in length, exhibits very similar proton transport properties. Molecular dynamics simulations on the microsecond time-scale carried out for the M2 channel provided atomic level details on the activation of the channel in response to protonation of the histidine residue, His37. The pathway of proton conduction is mediated by His37, which accepts and donates protons at different interconverting conformation states when pH is lower than 6.5. The Val27 and Trp41 gates and the salt bridge between Asp44 and Arg45 further enhance the directionality of proton transport. It is argued that the architecture and the mechanism of action similar to that found in the M2 channel might have been the perfect starting point for evolution towards the earliest proton pumps, indicating that active proton transport could have readily emerged from simple, passive proton channels.

Aromatic Amino Acids-Guanidinium Complexes through Cation-π Interactions

Continuing with our interest in the guanidinium group and the different interactions than can establish, we have carried out a theoretical study of the complexes formed by this cation and the aromatic amino acids (phenylalanine, histidine, tryptophan and tyrosine) using DFT methods and PCM-water solvation. Both hydrogen bonds and cation-π interactions have been found upon complexation. These interactions have been characterized by means of the analysis of the molecular electron density using the Atoms-in-Molecules approach as well as the orbital interactions using the Natural Bond Orbital methodology. Finally, the effect that the cation-π and hydrogen bond interactions exert on the aromaticity of the corresponding amino acids has been evaluated by calculating the theoretical NICS values, finding that the aromatic character was not heavily modified upon complexation.

Competition between n → π(Ar)* and conventional hydrogen bonding (N-H···N) interactions: an ab initio study of the complexes of 7-azaindole and fluorosubstituted pyridines

In this work, we have investigated a subtle competition between a very weak n → π(Ar)* interaction and a very strong hydrogen bond (N-H···N) interaction present in the complexes of 7-azaindole with a series of 2,6-substituted fluoropyridines and observed how the weak interaction modulates the overall structural motif of these complexes in the presence of the strong interaction. We have studied the structures and binding energies of these complexes using MP2 as well as dispersion-corrected DFT calculations. It has been found that the strength of the N-H···N interaction in these complexes decreases with increasing fluorination in the fluoropyridine ring while the proximity between the nitrogen atom in 7-azaindole and the aromatic ring of fluoropyridine increases through n → π(Ar)* interaction. Comparison of the binding energy values as well as structural parameters of these complexes calculated at the B3LYP level with those obtained at the MP2, M05-2X, and B97-D levels of theory clearly indicates that the dispersion effect is mostly responsible for this attractive n → π(Ar)* interaction. This conclusion is also supported by localized molecular orbital-energy decomposition analysis (LMO-EDA). The current investigation is the first theoretical study on the n → π(Ar)* interaction in the presence of a conventional strong hydrogen bonding interaction in the molecular system. Thus the present study has great significance for understanding the structures of the biomolecules as well as materials, as these interactions are very often present there simultaneously.