Non-classical helix-stabilizing interactions: C-H...O H-bonding between Phe and Glu side chains in alpha-helical peptides

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.

Comparative Study of Three Dyes' Adsorption onto Activated Carbon from Chenopodium quinoa Willd and Quillaja saponaria

The present study shows porous activated carbon obtained from Chenopodium quinoa Willd and Quillaja saponaria and their use as potential adsorbents to remove three types of dyes from aqueous solutions. The adsorption results were compared with commercial charcoal to check their efficiency. All porous carbon materials were activated using carbon dioxide and steam and fully characterized. Moreover, the steam-activated samples exhibited a high total pore volume with a BET surface area of around 800 m2 g−1. Batch adsorption experiments showed that commercial charcoal is the charcoal that offered the best adsorption efficiency for tartrazine and sunset yellow FCF. However, in the case of crystal violet, all activated carbons obtained from Chenopodium quinoa Willd and Quillaja saponaria showed the best captures, outperforming commercial charcoal. Molecular dockings of the dyes on the commercial charcoal surface were performed using AutoDock Vina. The kinetic results of the three isotherm’s models for the present data follow the order: Langmuir~Freundlich > Temkin.

Secondary Forces in Protein Folding

A complete inventory of the forces governing protein folding is critical for productive protein modeling, including structure prediction and de novo design, as well as understanding protein misfolding diseases of clinical significance. The dominant contributors to protein folding include the hydrophobic effect and conventional hydrogen bonding, along with Coulombic and van der Waals interactions. Over the past few decades, important additional contributors have been identified, including C-H···O hydrogen bonding, n→π* interactions, C5 hydrogen bonding, chalcogen bonding, and interactions involving aromatic rings (cation-π, X-H···π, π-π, anion-π, and sulfur-arene). These secondary contributions fall into two general classes: (1) weak but abundant interactions of the protein main chain and (2) strong but less frequent interactions involving protein side chains. Though interactions with high individual energies play important roles in specifying nonlocal molecular contacts and ligand binding, we estimate that weak but abundant interactions are likely to make greater overall contributions to protein folding, particularly at the level of secondary structure. Further research is likely to illuminate additional roles of these noncanonical interactions and could also reveal contributions yet unknown.

Highly Dynamic Anion-Quadrupole Networks in Proteins

The dynamics of anion-quadrupole (or anion-π) interactions formed between negatively charged (Asp/Glu) and aromatic (Phe) side chains are for the first time computationally characterized in RmlC (Protein Data Bank entry 1EP0 ), a homodimeric epimerase. Empirical force field-based molecular dynamics simulations predict anion-quadrupole pairs and triplets (anion-anion-π and anion-π-π) are formed by the protein during the simulated trajectory, which suggests that the anion-quadrupole interactions may provide a significant contribution to the overall stability of the protein, with an average of -1.6 kcal/mol per pair. Some anion-π interactions are predicted to form during the trajectory, extending the number of anion-quadrupole interactions beyond those predicted from crystal structure analysis. At the same time, some anion-π pairs observed in the crystal structure exhibit marginal stability. Overall, most anion-π interactions alternate between an "on" state, with significantly stabilizing energies, and an "off" state, with marginal or null stabilizing energies. The way proteins possibly compensate for transient loss of anion-quadrupole interactions is characterized in the RmlC aspartate 84-phenylalanine 112 anion-quadrupole pair observed in the crystal structure. A double-mutant cycle analysis of the thermal stability suggests a possible loss of anion-π interactions compensated by variations of hydration of the residues and formation of compensating electrostatic interactions. These results suggest that near-planar anion-quadrupole pairs can exist, sometimes transiently, which may play a role in maintaining the structural stability and function of the protein, in an otherwise very dynamic interplay of a nonbonded interaction network as well as solvent effects.

First steps in growth of a polypeptide toward β-sheet structure

The full conformational energy surface is examined for a molecule in which a dipeptide is attached to the same spacer group as another peptide chain, so as to model the seminal steps of β-sheet formation. This surface is compared with the geometrical preferences of the isolated dipeptide to extract the perturbations induced by interactions with the second peptide strand. These interpeptide interactions remove any tendency of the dipeptide to form a C5 ring structure, one of its two normally stable geometries. A C7 structure, the preferred conformation of the isolated dipeptide, remains as the global minimum in the full molecule. However, the stability of this structure is highly dependent upon interpeptide H-bonds with the second chain. The latter forces include not only the usual NH···O interaction, but also a pair of CH···O H-bonds. The secondary minimum is also of C7 type and likewise depends in part upon CH···O H-bonds for its stability. The latter interactions also play a part in the tertiary minimum. A two-strand β-sheet structure is not yet in evidence for this small model system, requiring additional peptide units to be added to each chain.

Carbon-oxygen hydrogen bonding in biological structure and function

Carbon-oxygen (CH···O) hydrogen bonding represents an unusual category of molecular interactions first documented in biological structures over 4 decades ago. Although CH···O hydrogen bonding has remained generally underappreciated in the biochemical literature, studies over the last 15 years have begun to yield direct evidence of these interactions in biological systems. In this minireview, we provide a historical context of biological CH···O hydrogen bonding and summarize some major advancements from experimental studies over the past several years that have elucidated the importance, prevalence, and functions of these interactions. In particular, we examine the impact of CH···O bonds on protein and nucleic acid structure, molecular recognition, and enzyme catalysis and conclude by exploring overarching themes and unresolved questions regarding unconventional interactions in biomolecular structure.

JC virus agnoprotein enhances large T antigen binding to the origin of viral DNA replication: evidence for its involvement in viral DNA replication

Agnoprotein is required for the successful completion of the JC virus (JCV) life cycle and was previously shown to interact with JCV large T-antigen (LT-Ag). Here, we further characterized agnoprotein's involvement in viral DNA replication. Agnoprotein enhances the DNA binding activity of LT-Ag to the viral origin (Ori) without directly interacting with DNA. The predicted amphipathic α-helix of agnoprotein plays a major role in this enhancement. All three phenylalanine (Phe) residues of agnoprotein localize to this α-helix and Phe residues in general are known to play critical roles in protein-protein interaction, protein folding and stability. The functional relevance of all Phe residues was investigated by mutagenesis. When all were mutated to alanine (Ala), the mutant virus (F31AF35AF39A) replicated significantly less efficiently than each individual Phe mutant virus alone, indicating the importance of Phe residues for agnoprotein function. Collectively, these studies indicate a close involvement of agnoprotein in viral DNA replication.

A survey of aspartate-phenylalanine and glutamate-phenylalanine interactions in the protein data bank: searching for anion-π pairs

Protein structures are stabilized using noncovalent interactions. In addition to the traditional noncovalent interactions, newer types of interactions are thought to be present in proteins. One such interaction, an anion-π pair, in which the positively charged edge of an aromatic ring interacts with an anion, forming a favorable anion-quadrupole interaction, has been previously proposed [Jackson, M. R., et al. (2007) J. Phys. Chem. B111, 8242-8249]. To study the role of anion-π interactions in stabilizing protein structure, we analyzed pairwise interactions between phenylalanine (Phe) and the anionic amino acids, aspartate (Asp) and glutamate (Glu). Particular emphasis was focused on identification of Phe-Asp or -Glu pairs separated by less than 7 Å in the high-resolution, nonredundant Protein Data Bank. Simplifying Phe to benzene and Asp or Glu to formate molecules facilitated in silico analysis of the pairs. Kitaura-Morokuma energy calculations were performed on roughly 19000 benzene-formate pairs and the resulting energies analyzed as a function of distance and angle. Edgewise interactions typically produced strongly stabilizing interaction energies (-2 to -7.3 kcal/mol), while interactions involving the ring face resulted in weakly stabilizing to repulsive interaction energies. The strongest, most stabilizing interactions were identified as preferentially occurring in buried residues. Anion-π pairs are found throughout protein structures, in helices as well as β strands. Numerous pairs also had nearby cation-π interactions as well as potential π-π stacking. While more than 1000 structures did not contain an anion-π pair, the 3134 remaining structures contained approximately 2.6 anion-π pairs per protein, suggesting it is a reasonably common motif that could contribute to the overall structural stability of a protein.

Amino acid pair- and triplet-wise groupings in the interior of α-helical segments in proteins

A statistical approach has been applied to analyse primary structure patterns at inner positions of α-helices in proteins. A systematic survey was carried out in a recent sample of non-redundant proteins selected from the Protein Data Bank, which were used to analyse α-helix structures for amino acid pairing patterns. Only residues more than three positions apart from both termini of the α-helix were considered as inner. Amino acid pairings i, i+k (k=1, 2, 3, 4, 5), were analysed and the corresponding 20×20 matrices of relative global propensities were constructed. An analysis of (i, i+4, i+8) and (i, i+3, i+4) triplet patterns was also performed. These analysis yielded information on a series of amino acid patterns (pairings and triplets) showing either high or low preference for α-helical motifs and suggested a novel approach to protein alphabet reduction. In addition, it has been shown that the individual amino acid propensities are not enough to define the statistical distribution of these patterns. Global pair propensities also depend on the type of pattern, its composition and orientation in the protein sequence. The data presented should prove useful to obtain and refine useful predictive rules which can further the development and fine-tuning of protein structure prediction algorithms and tools.

Noncovalent keystone interactions controlling biomembrane structure

There is a biomedical need to develop molecular recognition systems that selectively target the interfaces of protein and lipid aggregates in biomembranes. This is an extremely challenging problem in supramolecular chemistry because the biological membrane is a complex dynamic assembly of multifarious molecular components with local inhomogeneity. Two simplifying concepts are presented as a framework for basing molecular design strategies. The first generalization is that association of two binding partners in a biomembrane will be dominated by one type of non-covalent interaction which is referred to as the keystone interaction. Structural mutations in membrane proteins that alter the strength of this keystone interaction will likely have a major effect on biological activity and often will be associated with disease. The second generalization is to view the structure of a cell membrane as three spatial regions, that is, the polar membrane surface, the midpolar interfacial region and the non-polar membrane interior. Each region has a distinct dielectric, and the dominating keystone interaction between binding partners will be different. At the highly polar membrane surface, the keystone interactions between charged binding partners are ion-ion and ion-dipole interactions; whereas, ion-dipole and ionic hydrogen bonding are very influential at the mid-polar interfacial region. In the non-polar membrane interior, van der Waals forces and neutral hydrogen bonding are the keystone interactions that often drive molecular association. Selected examples of lipid and transmembrane protein association systems are described to illustrate how the association thermodynamics and kinetics are dominated by these keystone noncovalent interactions.

Effect of solvent upon CH...O hydrogen bonds with implications for protein folding

The series of CH...O bonds formed between CF(n)H(4-n) (n = 0-3) and water are studied by quantum calculations under vacuum and in various solvents, including aqueous environment. The results are compared with the OH...O bond of the water dimer in the same solvents. Increasing polarity of the solvent leads in all cases to a lessening of the H-bond interaction energy, in a uniform fashion such that the CH...O bonds all remain weaker than OH...O in any solvent. These H-bond weakenings are coupled to a shortening of the inter-subunit separation. The contraction of the covalent CH bond to the bridging proton is reduced as the solvent becomes more polar, and the blue shift of its stretching vibration is likewise diminished. A process is considered that simulates protein folding by starting from a pair of noninteracting subunits in aqueous solvent and then goes to a H-bonded pair within the confines of a protein environment. This process is found to be energetically more favorable for some of the CH...O H-bonds than for the nominally stronger conventional OH...O H-bond. This finding suggests that CH...O bonds can make important energetic contributions to protein folding, on par with those made by traditional H-bonds.

Aromatic interactions in peptides: Impact on structure and function

Aromatic interactions, including pi-pi, cation-pi, aryl-sulfur, and carbohydrate-pi interactions, have been shown to be prevalent in proteins through protein structure analysis, suggesting that they are important contributors to protein structure. However, the magnitude and significance of aromatic interactions is not defined by such studies. Investigation of aromatic interactions in the context of structured peptides has complemented studies of protein structure and has provided a wealth of information regarding the role of aromatic interactions in protein structure and function. Recent advances in this area are reviewed.

Association of a model transmembrane peptide containing gly in a heptad sequence motif

A peptide containing glycine at a and d positions of a heptad motif was synthesized to investigate the possibility that membrane-soluble peptides with a Gly-based, left-handed helical packing motif would associate. Based on analytical ultracentrifugation in C14-betaine detergent micelles, the peptide did associate in a monomer-dimer equilibrium, although the association constant was significantly less than that reported for the right-handed dimer of the glycophorin A transmembrane peptide in similar detergents. Fluorescence resonance energy transfer (FRET) experiments conducted on peptides labeled at their N-termini with either tetramethylrhodamine (TMR) or 7-nitrobenz-2-oxa-1,3-diazole (NBD) also indicated association. However, analysis of the FRET data using the usual assumption of complete quenching for NBD-TMR pairs in the dimer could not be quantitatively reconciled with the analytical ultracentrifugation-measured dimerization constant. This led us to develop a general treatment for the association of helices to either parallel or antiparallel structures of any aggregation state. Applying this treatment to the FRET data, constraining the dimerization constant to be within experimental uncertainty of that measured by analytical ultracentrifugation, we found the data could be well described by a monomer-dimer equilibrium with only partial quenching of the dimer, suggesting that the helices are most probably antiparallel. These results also suggest that a left-handed Gly heptad repeat motif can drive membrane helix association, but the affinity is likely to be less strong than the previously reported right-handed motif described for glycophorin A.

Cation-pi interactions studied in a model coiled-coil peptide

Cation-pi interactions between aromatic amino acids and the positively charged residues lysine and arginine have been proposed to play an important role in stabilizing protein structure. We have used a peptide that adopts a coiled coil structure as a model system to evaluate the energetic contribution of cation-pi interactions to protein folding. Peptides were designed in which phenylalanine, tyrosine, and tryptophan were placed at a solvent-exposed position of the helix, one turn removed from an arginine residue that could provide a favorable cation-pi interaction. Only the arginine-phenylalanine pairing provided significant stabilization of the peptide structure and it appears that hydrophobic packing, rather than the cation-pi effect, is more likely to be responsible for the stability of this peptide. We conclude that any stabilizing effect of cation-pi interactions in these peptides is much smaller than that predicted from computational studies.