Friend or foe: inter‐specific interactions and conflicts of interest within the family

1. Interactions between species can vary from mutually beneficial to evolutionarily neutral to antagonistic, even when the same two species are involved. Similarly, social interactions between members of the same species can lie on a spectrum from conflict to cooperation. 2. The aim of the present study was to investigate whether variation in the two types of social behaviour are interconnected. Is the fitness of the various classes of social partner within species (such as parent and offspring, or male and female) differently affected by interactions with a second species? Moreover, can inter‐specific interactions influence the outcome of social interactions within species? 3. The present experiments focus on the interactions between the burying beetle Nicrophorus vespilloides Herbst and the phoretic mite Poecilochirus carabi G. Canestrini & R. Canestrini. The approach was to measure the fitness of burying beetle mothers, fathers, and offspring after reproduction, which took place either in the presence or absence of mites. 4. We found that male, female, and larval burying beetles derive contrasting fitness costs and benefits from their interactions with the mite, despite sharing a common family environment. From the mite's perspective, its relationship with the burying beetle can, therefore, be simultaneously antagonistic, neutral, and possibly even mutualistic, depending on the particular family member involved. We also found that mites can potentially change the outcome of evolutionary conflicts within the family. 5. We conclude that inter‐specific interactions can explain some of the variation in social interactions seen within species. It is further suggested that intra‐specific interactions might contribute to variation in the outcome of interactions between species.     

Do Insectivorous Birds use Volatile Organic Compounds from Plants as Olfactory Foraging Cues? Three Experimental Tests

Some insectivorous birds orient towards insect‐defoliated trees even when they do not see the foliar damage or the herbivores. There are, however, only a few studies that have examined the mechanisms behind this foraging behaviour. Previous studies suggest that birds can use olfactory foraging cues (e.g. volatile organic compounds (VOCs) emitted by defoliated plants), indirect visual cues or a combination of the two sensory cues. VOCs from insect‐defoliated plants are known to attract natural enemies of herbivores, and researchers have hypothesized that VOCs could also act as olfactory foraging cues for birds. We conducted three experiments across a range of spatial scales to test this hypothesis. In each experiment, birds were presented with olfactory cues and their behavioural responses or foraging outcomes were observed. In the first experiment, two different VOC blends, designed to simulate the volatile emissions of mountain birch (Betula pubescens ssp. czerepanovii) after defoliation by autumnal moth (Epirrita autumnata) larvae, were used in behavioural experiments in aviaries with pied flycatchers (Ficedula hypoleuca). The second experiment was a field‐based trial of bird foraging efficiency; the same VOC blends were applied to mountain birches, silver birches (B. pendula) and European white birches (B. pubescens) with plasticine larvae attached to the trees to serve as artificial prey for birds and provide a means to monitor predation rate. In the third experiment, the attractiveness of silver birch saplings defoliated by autumnal moth larvae versus intact controls was tested with great tits (Parus major) and blue tits (Cyanistes caeruleus) in an aviary. Birds did not orient towards either artificial or real trees with VOC supplements or towards herbivore‐damaged saplings when these saplings and undamaged alternatives were hidden from view. These findings do not support the hypothesis that olfactory foraging cues are necessary in the attraction of birds to herbivore‐damaged trees.     

Predators marked with chemical cues from one prey have increased attack success on another prey species

1. To reduce the risk of being eaten by predators, prey alter their morphology or behaviour. This response can be tuned to the current danger if chemical or other cues associated with predators inform the prey about the risks involved. 2. It is well known that various prey species discriminate between chemical cues from predators that fed on conspecific prey and those that fed on heterospecific prey, and react stronger to the first. It is therefore expected that generalist predators are more successful in capturing a given prey species when they are contaminated with chemical cues from another prey species instead of cues from the same prey species. 3. Here, a generalist predatory mite was studied that feeds on thrips larvae as well as on whitefly eggs and crawlers. Mites were marked with cues (i.e. body fluids) of one of these two prey species and were subsequently offered thrips larva. 4. Predators marked with thrips cues killed significantly fewer thrips than predators marked with whitefly cues, even though the predator's tendency to attack was the same. In addition, more thrips larvae sought refuge in the presence of a predatory mite marked with thrips cues instead of whitefly cues. 5. This suggests that generalist predators may experience improved attack success when switching prey species.     

Loss of intramolecular electrostatic interactions and limited conformational ensemble may promote self‐association of cis–tau peptide

Self‐association of proteins can be triggered by a change in the distribution of the conformational ensemble. Posttranslational modification, such as phosphorylation, can induce a shift in the ensemble of conformations. In the brain of Alzheimer's disease patients, the formation of intra‐cellular neurofibrillary tangles deposition is a result of self‐aggregation of hyper‐phosphorylated tau protein. Biochemical and NMR studies suggest that the cis peptidyl prolyl conformation of a phosphorylated threonine‐proline motif in the tau protein renders tau more prone to aggregation than the trans isomer. However, little is known about the role of peptidyl prolyl cis/trans isomerization in tau aggregation. Here, we show that intra‐molecular electrostatic interactions are better formed in the trans isomer. We explore the conformational landscape of the tau segment containing the phosphorylated‐Thr²³¹‐Pro²³² motif using accelerated molecular dynamics and show that intra‐molecular electrostatic interactions are coupled to the isomeric state of the peptidyl prolyl bond. Our results suggest that the loss of intra‐molecular interactions and the more restricted conformational ensemble of the cis isomer could favor self‐aggregation. The results are consistent with experiments, providing valuable complementary atomistic insights and a hypothetical model for isomer specific aggregation of the tau protein. Proteins 2015; 83:436–444. © 2014 Wiley Periodicals, Inc.     

Conformational dynamics of nonsynonymous variants at protein interfaces reveals disease association

Recent studies have shown that the protein interface sites between individual monomeric units in biological assemblies are enriched in disease‐associated non‐synonymous single nucleotide variants (nsSNVs). To elucidate the mechanistic underpinning of this observation, we investigated the conformational dynamic properties of protein interface sites through a site‐specific structural dynamic flexibility metric (dfi) for 333 multimeric protein assemblies. dfi measures the dynamic resilience of a single residue to perturbations that occurred in the rest of the protein structure and identifies sites contributing the most to functionally critical dynamics. Analysis of dfi profiles of over a thousand positions harboring variation revealed that amino acid residues at interfaces have lower average dfi (31%) than those present at non‐interfaces (50%), which means that protein interfaces have less dynamic flexibility. Interestingly, interface sites with disease‐associated nsSNVs have significantly lower average dfi (23%) as compared to those of neutral nsSNVs (42%), which directly relates structural dynamics to functional importance. We found that less conserved interface positions show much lower dfi for disease nsSNVs as compared to neutral nsSNVs. In this case, dfi is better as compared to the accessible surface area metric, which is based on the static protein structure. Overall, our proteome‐wide conformational dynamic analysis indicates that certain interface sites play a critical role in functionally related dynamics (i.e., those with low dfi values), therefore mutations at those sites are more likely to be associated with disease. Proteins 2015; 83:428–435. © 2014 Wiley Periodicals, Inc.     

Quantitative delineation of how breathing motions open ligand migration channels in myoglobin and its mutants

Ligand migration and binding are central to the biological functions of many proteins such as myoglobin (Mb) and it is widely thought that protein breathing motions open up ligand channels dynamically. However, how a protein exerts its control over the opening and closing of these channels through its intrinsic dynamics is not fully understood. Specifically, a quantitative delineation of the breathing motions that are needed to open ligand channels is lacking. In this work, we present and apply a novel normal mode‐based method to quantitatively delineate what and how breathing motions open ligand migration channels in Mb and its mutants. The motivation behind this work springs from the observation that normal mode motions are closely linked to the breathing motions that are thought to open ligand migration channels. In addition, the method provides a direct and detailed depiction of the motions of each and every residue that lines a channel and can identify key residues that play a dominating role in regulating the channel. The all‐atom model and the full force‐field employed in the method provide a realistic energetics on the work cost required to open a channel, and as a result, the method can be used to efficiently study the effects of mutations on ligand migration channels and on ligand entry rates. Our results on Mb and its mutants are in excellent agreement with MD simulation results and experimentally determined ligand entry rates. Proteins 2015; 83:757–770. © 2015 Wiley Periodicals, Inc.     

Inferring the microscopic surface energy of protein–protein interfaces from mutation data

Mutations at protein–protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical values, yielding mean energies of ?48 cal mol^?1 Å^?2 for interactions between hydrophobic surfaces, ?51 to ?80 cal mol^?1 Å^?2 for surfaces of complementary charge, and 66–83 cal mol^?1 Å^?2 for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at ?24 cal mol^?1 Å^?2. Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid‐body interactions (r = 0.66). When used to rerank docking poses, it can place near‐native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50–95% of the cohesive energy. The model is available open‐source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/ . Proteins 2015; 83:640–650. © 2015 Wiley Periodicals, Inc.     

Survey of phosphorylation near drug binding sites in the Protein Data Bank (PDB) and their effects

While it is currently estimated that 40 to 50% of eukaryotic proteins are phosphorylated, little is known about the frequency and local effects of phosphorylation near pharmaceutical inhibitor binding sites. In this study, we investigated how frequently phosphorylation may affect the binding of drug inhibitors to target proteins. We examined the 453 non‐redundant structures of soluble mammalian drug target proteins bound to inhibitors currently available in the Protein Data Bank (PDB). We cross‐referenced these structures with phosphorylation data available from the PhosphoSitePlus database. Three hundred twenty‐two of 453 (71%) of drug targets have evidence of phosphorylation that has been validated by multiple methods or labs. For 132 of 453 (29%) of those, the phosphorylation site is within 12 Å of the small molecule‐binding site, where it would likely alter small molecule binding affinity. We propose a framework for distinguishing between drug‐phosphorylation site interactions that are likely to alter the efficacy of drugs versus those that are not. In addition we highlight examples of well‐established drug targets, such as estrogen receptor alpha, for which phosphorylation may affect drug affinity and clinical efficacy. Our data suggest that phosphorylation may affect drug binding and efficacy for a significant fraction of drug target proteins. Proteins 2015; 83:25–36. © 2014 Wiley Periodicals, Inc.     

3d interaction homology: The structurally known rotamers of tyrosine derive from a surprisingly limited set of information‐rich hydropathic interaction environments described by maps

Sidechain rotamer libraries are obtained through exhaustive statistical analysis of existing crystallographic structures of proteins and have been applied in multiple aspects of structural biology, for example, crystallography of relatively low‐resolution structures, in homology model building and in biomolecular NMR. Little is known, however, about the driving forces that lead to the preference or suitability of one rotamer over another. Construction of 3D hydropathic interaction maps for nearly 30,000 tyrosines reveals the environment around each, in terms of hydrophobic (π–π stacking, etc.) and polar (hydrogen bonding, etc.) interactions. After partitioning the tyrosines into backbone‐dependent (?, ψ) bins, a map similarity metric based on the correlation coefficient was applied to each map‐map pair to build matrices suitable for clustering with k‐means. The first bin (?200° ≤ ? < –155°; ?205° ≤ ψ < –160°), representing 631 tyrosines, reduced to 14 unique hydropathic environments, with most diversity arising from favorable hydrophobic interactions with many different residue partner types. Polar interactions for tyrosine include surprisingly ubiquitous hydrogen bonding with the phenolic OH and a handful of unique environments surrounding the tyrosine backbone. The memberships of all but one of the 14 environments are dominated (>50%) by a single χ₁/χ₂ rotamer. The last environment has weak or no interactions with the tyrosine ring and its χ₁/χ₂ rotamer is indeterminate, which is consistent with it being composed of mostly surface residues. Each tyrosine residue attempts to fulfill its hydropathic valence and thus, structural water molecules are seen in a variety of roles throughout protein structure. Proteins 2015; 83:1118–1136. © 2015 Wiley Periodicals, Inc.     

Effects of initial settings on computational protein–ligand docking accuracies for several docking programs

In this study, the influences of initial settings, i.e. initial conformations, configurations and docking parameters, on docking results were investigated. The conformations used in the study were generated by the CAMDAS program. After the conformational search calculations, five structures were selected from the conformer groups according to their conformation energies and root mean square deviations against crystal structures; for example, the lowest energy conformer, as well as the closest and farthest conformers to the crystal structure, was retrieved. Several docking parameter settings were used (default, high speed, generating 50 poses). In this study, docking calculations were conducted using the GOLD, eHiTS, AutoDock, AutoDock vina, FRED and DOCK programs. The success rates of GOLD, eHiTS and FRED were better than those of AutoDock, AutoDock vina and DOCK. The docking results using the farthest conformations were worse than those obtained using other conformations, indicating that some conformation search for the ligand molecule should be performed before the docking calculations.