Assessment of chemical-crosslink-assisted protein structure modeling in CASP13

With the advance of experimental procedures obtaining chemical crosslinking information is becoming a fast and routine practice. Information on crosslinks can greatly enhance the accuracy of protein structure modeling. Here, we review the current state of the art in modeling protein structures with the assistance of experimentally determined chemical crosslinks within the framework of the 13th meeting of Critical Assessment of Structure Prediction approaches. This largest-to-date blind assessment reveals benefits of using data assistance in difficult to model protein structure prediction cases. However, in a broader context, it also suggests that with the unprecedented advance in accuracy to predict contacts in recent years, experimental crosslinks will be useful only if their specificity and accuracy further improved and they are better integrated into computational workflows.     

A phylogenetic study of krill (Crustacea: Euphausiacea) reveals new taxa and co-evolution of morphological characters

The first comprehensive phylogenetic study of Euphausiacea (all 86 valid species) is presented. It is based on four molecular markers and 168 morphological characters (including 58 characters of the petasma). Phylogenetic analyses support the monophyly and robustness of the families Bentheuphausidae and Euphausiidae and reveal three major clades for which we erect three new subfamilies: Thysanopodinae, Euphausiinae and Nematoscelinae. All genus-level clades are statistically supported (except Thysanopoda in molecular analyses), deeply nested within the subfamily-level clades, and encompass 14 new species groups. Copulatory structures have a major impact on tree topology in the morphological analysis, the removal of which resulted in only half the number of supported clades and genera. We revealed three groups of morphological characters, which are probably coupled with the same biological role and thus interlinked evolutionarily: (i) antennular peduncle and petasma (copulation); (ii) eyes and anterior thoracopods (feeding); and (iii) shape of carapace and pleon (defence). We analysed the evolutionary pathways of the clades into main oceanic biotopes and compared them with morphological adaptations most likely to be coupled with this process.     

Latitudinal‐diversity gradients can be shaped by biotic processes: new insights from an eco‐evolutionary model

The processes involved in shaping latitudinal‐diversity gradients (LDGs) have been a longstanding source of debate and research. Climatic, historical and evolutionary factors have all been shown to contribute to the formation of LDGs. However, meta‐analyses have shown that different clades have LDG slopes that may vary in more than one order of magnitude. Such large variation cannot be explained solely by climatic or historical factors (e.g. difference in surface area between temperate and tropical zones) given that all clades within a geographic region are subject to the same conditions. Therefore, biotic processes intrinsic to each taxonomic group could be relevant in explaining rate differences in diversity decline across latitudinal gradients among groups. In this study, we developed a model simulating multiple competing species subjected (or not) to a demographic Allee effect. We simulated the range expansion of these species across an environmental gradient to show how these two overlooked factors (competition and Allee effects) are capable of modulating LDGs. Allee effects resulted in a steeper LDG given a higher probability of local extinction and lower colonization capacity compared to species without Allee effects. Likewise, stronger competition also led to a steeper decline in species diversity compared to scenarios with weaker species antagonistic interactions. This pattern occurred mostly due to the strength of priority effects, wherein scenarios with strong competition, species that dispersed earlier in the landscape were able to secure many patches whereas late‐arriving species were progressively precluded from expanding their ranges. Overall, our results suggest that the effect of biotic processes in shaping macroecological patterns could be more important than it is currently appreciated.     

Exploring protein-protein interactions using the site-identification by ligand competitive saturation methodology

Protein docking methods are powerful computational tools to study protein-protein interactions (PPI). While a significant number of docking algorithms have been developed, they are usually based on rigid protein models or with limited considerations of protein flexibility and the desolvation effect is rarely considered in docking energy functions, which may lower the accuracy of the predictions. To address these issues, we introduce a PPI energy function based on the site-identification by ligand competitive saturation (SILCS) framework and utilize the fast Fourier transform (FFT) correlation approach. The free energy content of the SILCS FragMaps represent an alternative to traditional energy grids and they can be efficiently utilized to guide FFT-based protein docking. Application of the approach to eight diverse test cases, including seven from Protein Docking Benchmark 5.0, showed the PPI prediction using SILCS approach (SILCS-PPI) to be competitive with several commonly used protein docking methods indicating that the method has the ability to both qualitatively and quantitatively inform the prediction of PPI. Results show the utility of the SILCS-PPI docking approach for determination of probability distributions of PPI interactions over the surface of both partner proteins, allowing for identification of alternate binding poses. Such binding poses are confirmed by experimental crystal contacts in our test cases. While more computationally demanding than available PPI docking technologies, we anticipate that the SILCS-PPI docking approach will offer an alternative methodology for improved evaluation of PPIs that could be used in a variety of fields from systems biology to excipient design for biologics-based drugs.     

Calorimetric study of the order-disorder conformational transition in succinoglycan

Thermally induced order-disorder conformational transition in succinoglycan was studied using the method of high-sensitivity differential scanning microcalorimetry within the range of polysaccharide concentrations from 0.1 to 3.5 mg mL⁻¹ at NaCl concentrations 0, 0.01, and 0.1M. The positions and shapes of the excess heat capacity curves depended substantially on both the NaCl and polysaccharide concentrations. At low polysaccharide concentrations in salt-free solution the experimental curves were closely approximated by the two-state model suggesting the transition mechanism to be of the single helix-coil type. With increasing polysaccharide and/or NaCl concentration, the experimental curves changed significantly in symmetry, which indicated a changing transition mechanism. At high polysaccharide concentrations or in the presence of the salt, the order-disorder transition of succinoglycan was shown to include two stages: the cooperative dissociation of the helix dimer and subsequent two-state melting of the helix monomer. The dependence of thermodynamic parameters for the dissociation and melting of helix structures in succinoglycan on NaCl and polysaccharide concentrations was obtained by fitting the experimental excess heat capacity curves. The cooperativity parameter σ for the single helix-coil transition as well as the average length of the helix segment of succinoglycan were calculated. Some features of succinoglycan ordering in solution are discussed. © 1996 John Wiley & Sons, Inc.     

Maturity and fertility of rhesus monkey oocytes collected at different intervals after an ovulatory stimulus (human chorionic gonadotropin) in in vitro fertilization cycles

In rhesus monkeys undergoing ovarian stimulation for in vitro fertilization (IVF), a midcycle injection of human chorionic gonadotropin (hCG) substitutes for the LH surge and induces preovulatory oocyte maturation. The time interval between injection and oocyte collection, ideally, allows for the completion of oocyte maturation without ovulation, which would reduce the number of oocytes available for harvest. To evaluate the influence of this time interval on oocyte parameters following hCG administration, we conducted a series of gonadotropin treatment protocols in 51 animals in which the interval from hCG administration to follicular aspiration was systematically varied from 27 to 36 hr. Follicle number and size, evaluated prior to hCG administration by sonography, did not vary significantly or consistently with preovulatory maturation time. Oocytes were harvested by laparotomy or laparoscopy, and scored for maturity before insemination. The percentage of mature, metaphase II (MII) oocytes at recovery increased significantly with increasing preovulatory time and was inversely proportional to that of metaphase I (MI) oocytes. However, oocyte yield tended toward a progressive decrease with increasing preovulatory maturation times from a high of 27 oocytes at 27 hr to a low of 17 oocytes/animal at the 36 hr time interval. Fertilization levels declined significantly from a high of 50% at 27 hr to a low of 30% at 36 hr. Thus, although higher percentages of mature oocytes were recovered at the longer time intervals, optimal oocyte/embryo harvests were realized after the shorter time intervals (27 and 32 hr) and are most compatible with the goal of achieving high yields of fertile oocytes and embryos following gonadotropin stimulation in rhesus monkeys. © 1996 Wiley-Liss, Inc.     

Asymmetric and transient properties of reciprocal activity of antagonists during the paw-shake response in the cat

Mutually inhibitory populations of neurons, half-center oscillators (HCOs), are commonly involved in the dynamics of the central pattern generators (CPGs) driving various rhythmic movements. Previously, we developed a multifunctional, multistable symmetric HCO model which produced slow locomotor-like and fast paw-shake-like activity patterns. Here, we describe asymmetric features of paw-shake responses in a symmetric HCO model and test these predictions experimentally. We considered bursting properties of the two model half-centers during transient paw-shake-like responses to short perturbations during locomotor-like activity. We found that when a current pulse was applied during the spiking phase of one half-center, let’s call it #1, the consecutive burst durations (BDs) of that half-center increased throughout the paw-shake response, while BDs of the other half-center, let’s call it #2, only changed slightly. In contrast, the consecutive interburst intervals (IBIs) of half-center #1 changed little, while IBIs of half-center #2 increased. We demonstrated that this asymmetry between the half-centers depends on the phase of the locomotor-like rhythm at which the perturbation was applied. We suggest that the fast transient response reflects functional asymmetries of slow processes that underly the locomotor-like pattern; e.g., asymmetric levels of inactivation across the two half-centers for a slowly inactivating inward current. We compared model results with those of in-vivo paw-shake responses evoked in locomoting cats and found similar asymmetries. Electromyographic (EMG) BDs of anterior hindlimb muscles with flexor-related activity increased in consecutive paw-shake cycles, while BD of posterior muscles with extensor-related activity did not change, and vice versa for IBIs of anterior flexors and posterior extensors. We conclude that EMG activity patterns during paw-shaking are consistent with the proposed mechanism producing transient paw-shake-like bursting patterns found in our multistable HCO model. We suggest that the described asymmetry of paw-shaking responses could implicate a multifunctional CPG controlling both locomotion and paw-shaking.     

Alterations of the glutathione redox cycle status in fumonisin B1-treated pig kidney cells

Fumonisin B₁ (FB₁) causes equine leukoencephalomalacia, porcine pulmonary edema, and liver tumors and chronic nephritis in rats. To investigate mechanisms by which FB₁ induces toxicity, effects of FB₁ on cellular glutathione (GSH) redox status and GSH depletion on FB₁ toxicity in pig kidney (LLC-PK₁) cells were studied. Treatment of LLC-PK₁ cells with 50 μM FB₁ for 24 hours significantly decreased cellular GSH contents from 56 ± 3.2 to 42.7 ± 4.4 nmol/mg protein (p < 0.05) and increased the activities of glutathione reductase (GR) from 25.7 ± 2.4 to 35.7 ± 5.0 μmol NADPH/mg protein (p < 0.05). The activities of glutathione peroxidase (GSHpx), catalase, and Cu,Zn-superoxide dismutase (SOD) were not changed by this treatment. Treatment of LLC-PK₁ cells for 12 hours with 0.1. mM buthionine sulfoximine (BSO), a selective inhibitor of the enzyme γ-glutamylcysteine synthetase that catalyzes the rate-limiting reaction in de novo GSH synthesis, decreased cellular GSH levels to about 20% of that found in the control cells. The cells pretreated with 0.1 mM BSO for 12 hours were significantly sensitized to the FB₁ cytotoxicity as determined by a long-term survival assay (p < 0.05). The results demonstrate that FB₁ changes GSH redox cycle status in LLC-PK₁ cells, and GSH may play a role in cytoprotection against FB₁ toxicity. © 1997 John Wiley & Sons, Inc.     

Endogenously produced catecholamines improve the regulatory function of TLR9-activated B cells

The sympathetic nervous system (SNS) contributes to immune balance by promoting anti-inflammatory B cells. However, whether B cells possess a self-regulating mechanism by which they modulate regulatory B cell (Breg) function is not well understood. In this study, we investigated the ability of B cells to synthesize their own catecholamines upon stimulation with different B cell activators and found that expression of the enzyme tyrosine hydroxylase (TH), required to generate catecholamines, is up-regulated by Toll-like receptor (TLR)9. This TLR9-dependent expression of TH correlated with up-regulation of adrenergic receptors (ADRs), enhanced interleukin (IL)-10 production, and overexpression of the co-inhibitory ligands programmed death ligand 1 (PD-L1) and Fas ligand (FasL). Moreover, concomitant stimulation of ß1-3-ADRs together with a B cell receptor (BCR)/TLR9 stimulus clearly enhances the anti-inflammatory potential of Bregs to suppress CD4 T cells, a crucial population in the pathogenesis of autoimmune diseases, like rheumatoid arthritis (RA). Furthermore, TH up-regulation was also demonstrated in B cells during the course of collagen-induced arthritis (CIA), a mouse model for the investigation of RA. In conclusion, our data show that B cells possess an autonomous mechanism to modulate their regulatory function in an autocrine and/or paracrine manner. These findings help to better understand the function of B cells in the regulation of autoimmune diseases and the interplay of SNS.

The sympathetic nervous system produces neurotransmitters such as catecholamines which contribute to immune balance by promoting anti-inflammatory B cells. This study shows that mouse B cells can themselves synthesize, sense, and transport catecholamines, which in turn modulate regulatory B cell function in an autocrine and/or paracrine manner to suppress T cell proliferation.