Triathlon for energy functions: Who is the winner for design of protein–protein interactions?

Computational prediction of stabilizing mutations into monomeric proteins has become an almost ordinary task. Yet, computational stabilization of protein–protein complexes remains a challenge. Design of protein–protein interactions (PPIs) is impeded by the absence of an energy function that could reliably reproduce all favorable interactions between the binding partners. In this work, we present three energy functions: one function that was trained on monomeric proteins, while the other two were optimized by different techniques to predict side-chain conformations in a dataset of PPIs. The performances of these energy functions are evaluated in three different tasks related to design of PPIs: predicting side-chain conformations in PPIs, recovering native binding-interface sequences, and predicting changes in free energy of binding due to mutations. Our findings show that both functions optimized on side-chain repacking in PPIs are more suitable for PPI design compared to the function trained on monomeric proteins. Yet, no function performs best at all three tasks. Comparison of the three energy functions and their performances revealed that (1) burial of polar atoms should not be penalized significantly in PPI design as in single-protein design and (2) contribution of electrostatic interactions should be increased several-fold when switching from single-protein to PPI design. In addition, the use of a softer van der Waals potential is beneficial in cases when backbone flexibility is important. All things considered, we define an energy function that captures most of the nuances of the binding energetics and hence, should be used in future for design of PPIs.     

Are protein–protein interfaces more conserved in sequence than the rest of the protein surface?

Protein interfaces are thought to be distinguishable from the rest of the protein surface by their greater degree of residue conservation. We test the validity of this approach on an expanded set of 64 protein-protein interfaces using conservation scores derived from two multiple sequence alignment types, one of close homologs/orthologs and one of diverse homologs/paralogs. Overall, we find that the interface is slightly more conserved than the rest of the protein surface when using either alignment type, with alignments of diverse homologs showing marginally better discrimination. However, using a novel surface-patch definition, we find that the interface is rarely significantly more conserved than other surface patches when using either alignment type. When an interface is among the most conserved surface patches, it tends to be part of an enzyme active site. The most conserved surface patch overlaps with 39% (+/- 28%) and 36% (+/- 28%) of the actual interface for diverse and close homologs, respectively. Contrary to results obtained from smaller data sets, this work indicates that residue conservation is rarely sufficient for complete and accurate prediction of protein interfaces. Finally, we find that obligate interfaces differ from transient interfaces in that the former have significantly fewer alignment gaps at the interface than the rest of the protein surface, as well as having buried interface residues that are more conserved than partially buried interface residues.     

Why is an energy metabolic defect the common outcome in BMFS?

Inherited bone marrow failure syndromes (BMFS) are rare, distressing, inherited blood disorders of children. Although the genetic origin of these pathologies involves genes with different functions, all are associated with progressive haematopoietic impairment and an excessive risk of malignancies. Defects in energy metabolism induce oxidative stress, impaired energy production and an unbalanced ratio between ATP and AMP. This assumes an important role in self-renewal and differentiation in haematopoietic stem cells (HSC) and can play an important role in bone marrow failure. Defects in energetic/respiratory metabolism, in particular in FA and SDS cells, have been described recently and seem to be a pertinent argument in the discussion of the haematopoietic defect in BMFS, as an alternative to the hypotheses already established on this subject, which may shed new light on the evolution of these diseases.     

Is the bovine lysosomal phospholipase B‐like protein an amidase?

The main function of lysosomal proteins is to degrade cellular macromolecules. We purified a novel lysosomal protein to homogeneity from bovine kidneys. By gene annotation, this protein is defined as a bovine phospholipase B‐like protein 1 (bPLBD1) and, to better understand its biological function, we solved its structure at 1.9 Å resolution. We showed that bPLBD1 has uniform noncomplex‐type N‐glycosylation and that it localized to the lysosome. The first step in lysosomal protein transport, the initiation of mannose‐6‐phosphorylation by a N‐acetylglucosamine‐1‐phosphotransferase, requires recognition of at least two distinct lysines on the protein surface. We identified candidate lysines by analyzing the structural and sequentially conserved N‐glycosylation sites and lysines in bPLBD1 and in the homologous mouse PLBD2. Our model suggests that N408 is the primarily phosphorylated glycan, and K358 a key residue for N‐acetylglucosamine‐1‐phosphotransferase recognition. Two other lysines, K334 and K342, provide the required second site for N‐acetylglucosamine‐1‐phosphotransferase recognition. bPLBD1 is an N‐terminal nucleophile (Ntn) hydrolase. By comparison with other Ntn‐hydrolases, we conclude that the acyl moiety of PLBD1 substrate must be small to fit the putative binding pocket, whereas the space for the rest of the substrate is a large open cleft. Finally, as all the known substrates of Ntn‐hydrolases have amide bonds, we suggest that bPLBD1 may be an amidase or peptidase instead of lipase, explaining the difficulty in finding a good substrate for any members of the PLBD family. Proteins 2014; 82:300–311. © 2013 Wiley Periodicals, Inc.     

“How pendulum‐like are siamangs? energy exchange during brachiation”

Hylobatidae (gibbons and siamangs) are known for their brachiation skills. The comparison of brachiation with a pendulum is made several times in the literature, and the costs and benefits of being pendulum-like are well described. However, the amount of energy exchange during brachiation of gibbons has rarely been determined. In this study, the amount of energy recovery (ER) during brachiation is assessed for three siamangs in a seminatural environment. The animals were recorded by four cameras while voluntarily brachiating on three different setups. The effects of locomotion speed, brachiation type, and setup on ER as well as on the external mechanical work during brachiation are determined. It is hypothesized that the amount of ER decreases with an increasing setup complexity while the external mechanical work increases. Additionally, we expect that support arm kinematics will be adjusted according to spatial complexity in order to maintain high recovery percentages. Our results show that ER is mainly determined by brachiation speed. Regardless of type of brachiation or setup, brachiation is done with a lower ER when brachiating faster. Within our limited range of setup variation, the expected effect of increasing complexity is not found. Although there is significant variation in support arm joint angles, no clear relation with speed, brachiation type, or setup is observed.     

H‐bond network optimization in protein–protein complexes: Are all‐atom force field scores enough?

Structural prediction of protein-protein complexes given the structures of the two interacting compounds in their unbound state is a key problem in biophysics. In addition to the problem of sampling of near-native orientations, one of the modeling main difficulties is to discriminate true from false positives. Here, we present a hierarchical protocol for docking refinement able to discriminate near native poses from a group of docking candidates. The main idea is to combine an efficient sampling of the full system hydrogen bond network and side chains, together with an all-atom force field and a surface generalized born implicit solvent. We tested our method on a set of twenty two complexes containing a near-native solution within the top 100 docking poses, obtaining a near native solution as the top pose in 70% of the cases. We show that all atom force fields optimized H-bond networks do improve significantly state of the art scoring functions.     

How common is common? Rapidly assessing population size and structure of the frog Mantidactylus betsileanus at a site in east‐central Madagascar

Population‐level data are urgently needed for amphibians in light of the ongoing amphibian extinction crisis. Studies focused on population dynamics are not only important for rare species but also for common species which shape ecosystems to a greater degree than those that are rare. Some of the greatest global amphibian species diversity is found in Madagascar, yet there are few studies on the ecology of frog species on the island. We carried out a mark‐recapture study on the widespread frog Mantidactylus betsileanus (Mantellidae: Mantellinae: Mantidactylus) at two adjacent rainforest sites in east‐central Madagascar to assess its population size and structure. To do so, we validated and implemented an individual identification protocol using photographs of the ventral patterns of frogs and identified individuals with photographic‐matching software. Using this rapid, non‐invasive survey method, we were able to estimate a density of 26 and 28 frogs per 100 m² at each of the two sites sampled. Our results show the rainforests near the village of Andasibe, Madagascar support remarkably high amphibian abundance, helping illustrate the significant ecological role of frogs in this ecosystem. Further, individual frog markings allowed us to develop more precise estimates than traditional survey methods. This study provides a blueprint to augment existing population studies or develop new monitoring programs in Madagascar and beyond.     

Dual energy landscape: the functional state of the β-barrel outer membrane protein G molds its unfolding energy landscape

We applied dynamic single-molecule force spectroscopy to quantify the parameters (free energy of activation and distance of the transition state from the folded state) characterizing the energy barriers in the unfolding energy landscape of the outer membrane protein G (OmpG) from Escherichia coli. The pH-dependent functional switching of OmpG directs the protein along different regions on the unfolding energy landscape. The two functional states of OmpG take the same unfolding pathway during the sequential unfolding of β-hairpins I-IV. After the initial unfolding events, the unfolding pathways diverge. In the open state, the unfolding of β-hairpin V in one step precedes the unfolding of β-hairpin VI. In the closed state, β-hairpin V and β-strand S11 with a part of extracellular loop L6 unfold cooperatively, and subsequently β-strand S12 unfolds with the remaining loop L6. These two unfolding pathways in the open and closed states join again in the last unfolding step of β-hairpin VII. Also, the conformational change from the open to the closed state witnesses a rigidified extracellular gating loop L6. Thus, a change in the conformational state of OmpG not only bifurcates its unfolding pathways but also tunes its mechanical properties for optimum function.     

Do electrostatic interactions destabilize protein–nucleic acid binding?

The negatively charged phosphates of nucleic acids are often paired with positively charged residues upon binding proteins. It was thus counter-intuitive when previous Poisson-Boltzmann (PB) calculations gave positive energies from electrostatic interactions, meaning that they destabilize protein-nucleic acid binding. Our own PB calculations on protein-protein binding have shown that the sign and the magnitude of the electrostatic component are sensitive to the specification of the dielectric boundary in PB calculations. A popular choice for the boundary between the solute low dielectric and the solvent high dielectric is the molecular surface; an alternative is the van der Waals (vdW) surface. In line with results for protein-protein binding, in this article, we found that PB calculations with the molecular surface gave positive electrostatic interaction energies for two protein-RNA complexes, but the signs are reversed when the vdW surface was used. Therefore, whether destabilizing or stabilizing effects are predicted depends on the choice of the dielectric boundary. The two calculation protocols, however, yielded similar salt effects on the binding affinity. Effects of charge mutations differentiated the two calculation protocols; PB calculations with the vdW surface had smaller deviations overall from experimental data.     

Extending the stress‐gradient hypothesis – is competition among animals less common in harsh environments?

The role of positive interactions has become widely accepted as a mechanism shaping community dynamics. Most empirical evidence comes from plant communities and sessile marine organisms. However, evidence for the relative role of positive interactions in organizing terrestrial animal communities is more limited, and a general framework that includes positive interactions among animals is lacking. The ‘stress gradient hypothesis’ (SGH) developed by plant ecologists predicts that the balance between positive and negative interactions will vary along gradients of biotic and abiotic stress, with positive interactions being more important in stressful environments. Paralleling the SGH, stress gradients for terrestrial herbivores could be equated to inverse primary productivity gradients, so we would expect positive interactions to prevail in more stressful, low productivity environments. However, this contradicts the typical view of terrestrial animal ecology that low primary productivity systems will foster intense competition for resources among consumers. Here we use alpine herbivores as a case study to test one of the predictions of the SGH in animal communities, namely the prevalence of positive interactions in low productivity environments. We identify potential mechanisms of facilitation and review the limited number of examples of interspecific interactions among alpine herbivores to assess the role of positive and negative interactions in structuring their communities. A meta‐analysis showed no clear trend in the strength and direction of interactions among alpine herbivores. Although studies were biased towards reporting significant negative inter actions, we found no evidence of competition dominating in harsh environments. Thus, our results only partially support the SGH, but directly challenge the dominant view among animal ecologists. Clearly, a sound theoretical framework is needed to include competition, positive and neutral interactions as potential mechanisms determining the structure of animal communities under differing environmental conditions, and the stress‐gradient hypothesis can provide a solid starting point.     

How thick is the layer of thermal volume surrounding the protein?

Investigation on the volume properties of protein hydration layers is reported. Presented results are based on combination of Monte Carlo modeling and available experimental data. Six globular proteins with known data are chosen for analysis. Analyzing the model and the experimental results we found that water molecules bound to proteins by hydrogen bond are preferentially located at the places with local depressions on the protein surface. Consequently, the hydration level is not strictly proportional to the area of charged and polar surfaces, but also depends on the shape of the molecular surface. The thickness of the thermal volume layer as calculated in the framework of the scaled particle theory is 0.6-0.65 A for chosen proteins. The obtained value is significantly lower than that presented for proteins in earlier papers (where proportionality between the hydration level and the area of charged and polar surfaces was assumed), but is close to the value published for small solute molecules. Discussion including the influence of protein size and the thermal motion of the surface is presented.     

How does climate warming affect plant‐pollinator interactions?

Climate warming affects the phenology, local abundance and large-scale distribution of plants and pollinators. Despite this, there is still limited knowledge of how elevated temperatures affect plant-pollinator mutualisms and how changed availability of mutualistic partners influences the persistence of interacting species. Here we review the evidence of climate warming effects on plants and pollinators and discuss how their interactions may be affected by increased temperatures. The onset of flowering in plants and first appearance dates of pollinators in several cases appear to advance linearly in response to recent temperature increases. Phenological responses to climate warming may therefore occur at parallel magnitudes in plants and pollinators, although considerable variation in responses across species should be expected. Despite the overall similarities in responses, a few studies have shown that climate warming may generate temporal mismatches among the mutualistic partners. Mismatches in pollination interactions are still rarely explored and their demographic consequences are largely unknown. Studies on multi-species plant-pollinator assemblages indicate that the overall structure of pollination networks probably are robust against perturbations caused by climate warming. We suggest potential ways of studying warming-caused mismatches and their consequences for plant-pollinator interactions, and highlight the strengths and limitations of such approaches.     

When is translocation required for the population recovery of old‐growth epiphytes in a reforested landscape?

It is often assumed that species recolonization follows from the restoration of key habitat structure. Thus, forest restoration focuses on the recovery of trees into deforested landscapes, so that a multitude of associated organisms can achieve “colonization credit” and recolonize from remnant source populations into restored habitat. This opportunity for recolonization exists because species vulnerable to habitat loss may experience an “extinction debt,” during which their remnant populations decline only slowly to equilibrium with a deforested landscape. These persistent but declining populations become propagule sources for recolonization. To test limits to “colonization credit,” this study focused on old‐growth dependent lichen epiphytes, using a simulation to identify a hypothetical threshold at which: (1) the number of remnant populations, and (2) their population sizes, are too low to achieve recolonization and population recovery, despite efforts placed into forest restoration. The results show that for a landscape scenario relevant to the industrialized temperate zone, with less than 5% of old‐growth forest remaining, and ambitions for restoration to circa 10–15% forest cover, there is a failure to achieve population recovery over long timescales (i.e. within 600 years), making translocation a necessary option. This delay represents a “colonization deficit” that may be a common feature in ecological restoration more generally.     

How important are Toll‐like receptors for antimicrobial responses?

The innate immune system is the primary line of defence against invading pathogenic microbes. Toll-like receptors (TLRs) are a family of membrane receptors which play a pivotal role in sensing a wide range of invading pathogens including bacteria, fungi and viruses. TLR-deficient mice have provided us with immense knowledge on the functioning of individual TLRs. Dysregulation of TLR signalling is linked with a number of disease conditions. Disease models have helped show that targeting components of TLR signalling cascades could lead to novel therapies in the treatment of infectious diseases. In this review we focus on the evidence provided to date to explain just how important TLRs are in host defence against microbial pathogens.     

How essential is the ‘essential’ active-site lysine in dihydrodipicolinate synthase?

Dihydrodipicolinate synthase (DHDPS, E.C. 4.2.1.52), a validated antibiotic target, catalyses the first committed step in the lysine biosynthetic pathway: the condensation reaction between (S)-aspartate β-semialdehyde [(S)-ASA] and pyruvate via the formation of a Schiff base intermediate between pyruvate and the absolutely conserved active-site lysine. Escherichia coli DHDPS mutants K161A and K161R of the active-site lysine were characterised for the first time. Unexpectedly, the mutant enzymes were still catalytically active, albeit with a significant decrease in activity. The k_cat values for DHDPS-K161A and DHDPS-K161R were 0.06 ± 0.02 s⁻¹ and 0.16 ± 0.06 s⁻¹ respectively, compared to 45 ± 3 s⁻¹ for the wild-type enzyme. Remarkably, the K_M values for pyruvate increased by only 3-fold for DHDPS-K161A and DHDPS-K161R (0.45 ± 0.04 mM and 0.57 ± 0.06 mM, compared to 0.15 ± 0.01 mM for the wild-type DHDPS), while the K_M values for (S)-ASA remained the same for DHDPS-K161R (0.12 ± 0.01 mM) and increased by only 2-fold for DHDPS-K161A (0.23 ± 0.02 mM) and the K_i for lysine was unchanged. The X-ray crystal structures of DHDPS-K161A and DHDPS-K161R were solved at resolutions of 2.0 and 2.1 Å respectively and showed no changes in their secondary or tertiary structures when compared to the wild-type structure. The crystal structure of DHDPS-K161A with pyruvate bound at the active site was solved at a resolution of 2.3 Å and revealed a defined binding pocket for pyruvate that is thus not dependent upon lysine 161. Taken together with ITC and NMR data, it is concluded that although lysine 161 is important in the wild-type DHDPS-catalysed reaction, it is not absolutely essential for catalysis.     

How is the balance between protein synthesis and degradation achieved?

Unlike most substances that cells manufacture, proteins are not produced and broken down by a common series of chemical reactions, but by completely different (independent and disconnected) mechanisms that possess no intrinsic means of making the rates of the two processes equal and attaining steady state concentrations. Balance between them is achieved extrinsically and is often imagined today to be the result of the actions of chemical feedback agents. But however instantiated, chemical feedback or any similar mechanism can only rectify induced imbalances in a system previously balanced by other means. Those "other means" necessarily involve reversible mass action or equilibrium-based interactions between native and altered forms of protein molecules somewhere in time and space between their synthesis and degradation.     

How consistent is RAD‐seq divergence with DNA‐barcode based clustering in insects?

Promoted by the barcoding approach, mitochondrial DNA is more than ever used as a molecular marker to identify species boundaries. Yet, it has been repeatedly argued that it may be poorly suited for this purpose, especially in insects where mitochondria are often associated with invasive intracellular bacteria that may promote their introgression. Here, we inform this debate by assessing how divergent nuclear genomes can be when mitochondrial barcodes indicate very high proximity. To this end, we obtained RAD‐seq data from 92 barcode‐based species‐like units (operational taxonomic units [OTUs]) spanning four insect orders. In 100% of the cases, the observed median nuclear divergence was lower than 2%, a value that was recently estimated as one below which nuclear gene flow is not uncommon. These results suggest that although mitochondria may occasionally leak between species, this process is rare enough in insects to make DNA barcoding a reliable tool for clustering specimens into species‐like units.