Monitoring protein-protein interactions between the mammalian integral membrane transporters and PDZ-interacting partners using a modified split-ubiquitin membrane yeast two-hybrid system

PDZ-binding motifs are found in the C-terminal tails of numerous integral membrane proteins where they mediate specific protein-protein interactions by binding to PDZ-containing proteins. Conventional yeast two-hybrid screens have been used to probe protein-protein interactions of these soluble C termini. However, to date no in vivo technology has been available to study interactions between the full-length integral membrane proteins and their cognate PDZ-interacting partners. We previously developed a split-ubiquitin membrane yeast two-hybrid (MYTH) system to test interactions between such integral membrane proteins by using a transcriptional output based on cleavage of a transcription factor from the C terminus of membrane-inserted baits. Here we modified MYTH to permit detection of C-terminal PDZ domain interactions by redirecting the transcription factor moiety from the C to the N terminus of a given integral membrane protein thus liberating their native C termini. We successfully applied this "MYTH 2.0" system to five different mammalian full-length renal transporters and identified novel PDZ domain-containing partners of the phosphate (NaPi-IIa) and sulfate (NaS1) transporters that would have otherwise not been detectable. Furthermore this assay was applied to locate the PDZ-binding domain on the NaS1 protein. We showed that the PDZ-binding domain for PDZK1 on NaS1 is upstream of its C terminus, whereas the two interacting proteins, NHERF-1 and NHERF-2, bind at a location closer to the N terminus of NaS1. Moreover NHERF-1 and NHERF-2 increased functional sulfate uptake in Xenopus oocytes when co-expressed with NaS1. Finally we used MYTH 2.0 to demonstrate that the NaPi-IIa transporter homodimerizes via protein-protein interactions within the lipid bilayer. In summary, our study establishes the MYTH 2.0 system as a novel tool for interactive proteomics studies of membrane protein complexes.

Positive and negative mutant selection in the human histone hairpin-binding protein using the yeast three-hybrid system.

We have used the yeast three-hybrid system in a positive selection for mutants of the human histone hairpin-binding protein (HBP) capable of interacting with non-canonical hairpins and in a negative selection for loss-of-binding mutants. Interestingly, all mutations from the positive selection are located in the N- and C-terminal regions flanking a minimal RNA-binding domain (RBD) previously defined between amino acids 126 and 198. Further, in vitro binding studies demonstrate that the RBD, which shows no obvious similarity to other RNA-binding motifs, has a relaxed sequence specificity compared to full-length HBP, allowing it to bind to mutant hairpin RNAs not normally found in histone genes. These findings indicate that the sequences flanking the RBD are important for restricting binding to the highly conserved histone hairpin structure. Among the loss-of-binding mutations, about half are nonsense mutations distributed throughout the N-terminal part and the RBD whereas the other half are missense mutations restricted to the RBD. Whereas the nonsense mutations permit a more precise definition of the C-terminal border of the RBD, the missense mutations identify critical residues for RNA binding within the RBD.

A real time simulation model of glucose-insulin metabolism for type 1 diabetes patients

In this paper, a simulation model of glucose-insulin metabolism for Type 1 diabetes patients is presented. The proposed system is based on the combination of Compartmental Models (CMs) and artificial Neural Networks (NNs). This model aims at the development of an accurate system, in order to assist Type 1 diabetes patients to handle their blood glucose profile and recognize dangerous metabolic states. Data from a Type 1 diabetes patient, stored in a database, have been used as input to the hybrid system. The data contain information about measured blood glucose levels, insulin intake, and description of food intake, along with the corresponding time. The data are passed to three separate CMs, which produce estimations about (i) the effect of Short Acting (SA) insulin intake on blood insulin concentration, (ii) the effect of Intermediate Acting (IA) insulin intake on blood insulin concentration, and (iii) the effect of carbohydrate intake on blood glucose absorption from the gut. The outputs of the three CMs are passed to a Recurrent NN (RNN) in order to predict subsequent blood glucose levels. The RNN is trained with the Real Time Recurrent Learning (RTRL) algorithm. The resulted blood glucose predictions are promising for the use of the proposed model for blood glucose level estimation for Type 1 diabetes patients.

GABARAP deficiency modulates expression of NaPi-IIa in renal brush-border membranes

Renal reabsorption of inorganic phosphate (P(i)) is mainly mediated by the Na(+)-dependent P(i)-cotransporter NaPi-IIa that is expressed in the brush-border membrane (BBM) of renal proximal tubules. Regulation and apical expression of NaPi-IIa are known to depend on a network of interacting proteins. Most of the interacting partners identified so far associate with the COOH-terminal PDZ-binding motif (TRL) of NaPi-IIa. In this study GABA(A) receptor-associated protein (GABARAP) was identified as a novel interacting partner of NaPi-IIa applying a membrane yeast-two-hybrid system (MYTH 2.0) to screen a mouse kidney library with the TRL-truncated cotransporter as bait. GABARAP mRNA and protein are present in renal tubules, and the interaction of NaPi-IIa and GABARAP was confirmed by using glutathione S-transferase pulldowns from BBM and coimmunoprecipitations from transfected HEK293 cells. Amino acids 36-68 of GABARAP were identified as the determinant for the described interaction. The in vivo effects of this interaction were studied in a murine model. GABARAP(-/-) mice have reduced urinary excretion of P(i), higher Na(+)-dependent (32)P(i) uptake in BBM vesicles, and increased expression of NaPi-IIa in renal BBM compared with GABARAP(+/+) mice. The expression of Na(+)/H(+) exchanger regulatory factor (NHERF)1, an important scaffold for the apical expression of NaPi-IIa, is also increased in GABARAP(-/-) mice. The absence of GABARAP does not interfere with the regulation of the cotransporter by either parathyroid hormone or acute changes of dietary P(i) content.

A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

Paleoclimate data-model comparison and the role of climate forcings over the past 1500 years

The past 1500 years provide a valuable opportunity to study the response of the climate system to external forcings. However, the integration of paleoclimate proxies with climate modeling is critical to improving the understanding of climate dynamics. In this paper, a climate system model and proxy records are therefore used to study the role of natural and anthropogenic forcings in driving the global climate. The inverse and forward approaches to paleoclimate data–model comparison are applied, and sources of uncertainty are identified and discussed. In the first of two case studies, the climate model simulations are compared with multiproxy temperature reconstructions. Robust solar and volcanic signals are detected in Southern Hemisphere temperatures, with a possible volcanic signal detected in the Northern Hemisphere. The anthropogenic signal dominates during the industrial period. It is also found that seasonal and geographical biases may cause multiproxy reconstructions to overestimate the magnitude of the long-term preindustrial cooling trend. In the second case study, the model simulations are compared with a coral δ18O record from the central Pacific Ocean. It is found that greenhouse gases, solar irradiance, and volcanic eruptions all influence the mean state of the central Pacific, but there is no evidence that natural or anthropogenic forcings have any systematic impact on El Niño–Southern Oscillation. The proxy climate relationship is found to change over time, challenging the assumption of stationarity that underlies the interpretation of paleoclimate proxies. These case studies demonstrate the value of paleoclimate data–model comparison but also highlight the limitations of current techniques and demonstrate the need to develop alternative approaches.

Influence of the Central American Seaway and Drake Passage on ocean circulation and neodymium isotopes: A model study

The sensitivity of the neodymium isotopic composition (ϵNd) to tectonic rearrangements of seaways is investigated using an Earth System Model of Intermediate Complexity. The shoaling and closure of the Central American Seaway (CAS) is simulated, as well as the opening and deepening of Drake Passage (DP). Multiple series of equilibrium simulations with various intermediate depths are performed for both seaways, providing insight into ϵNd and circulation responses to progressive throughflow evolutions. Furthermore, the sensitivity of these responses to the Atlantic Meridional Overturning Circulation (AMOC) and the neodymium boundary source is examined. Modeled ϵNd changes are compared to sediment core and ferromanganese (Fe-Mn) crust data. The model results indicate that the North Atlantic ϵNd response to the CAS shoaling is highly dependent on the AMOC state, i.e., on the AMOC strength before the shoaling to shallow depths (preclosure). Three scenarios based on different AMOC forcings are discussed, of which the model-data agreement favors a shallow preclosure (Miocene) AMOC (∼6 Sv). The DP opening causes a rather complex circulation response, resulting in an initial South Atlantic ϵNd decrease preceding a larger increase. This feature may be specific to our model setup, which induces a vigorous CAS throughflow that is strongly anticorrelated to the DP throughflow. In freshwater experiments following the DP deepening, ODP Site 1090 is mainly influenced by AMOC and DP throughflow changes, while ODP Site 689 is more strongly influenced by Southern Ocean Meridional Overturning Circulation and CAS throughflow changes. The boundary source uncertainty is largest for shallow seaways and at shallow sites.

Earth system commitments due to delayed mitigation

As long as global CO₂ emissions continue to increase annually, long-term committed Earth system changes grow much faster than current observations. A novel metric linking this future growth to policy decisions today is the mitigation delay sensitivity (MDS), but MDS estimates for Earth system variables other than peak temperature (ΔT max) are missing. Using an Earth System Model of Intermediate Complexity, we show that the current emission increase rate causes a ΔT max increase roughly 3–7.5 times as fast as observed warming, and a millenial steric sea level rise (SSLR) 7–25 times as fast as observed SSLR, depending on the achievable rate of emission reductions after the peak of emissions. These ranges are only slightly affected by the uncertainty range in equilibrium climate sensitivity, which is included in the above values. The extent of ocean acidification at the end of the century is also strongly dependent on the starting time and rate of emission reductions. The preservable surface ocean area with sufficient aragonite supersaturation for coral reef growth is diminished globally at an MDS of roughly 25%–80% per decade. A near-complete loss of this area becomes unavoidable if mitigation is delayed for a few years to decades. Also with respect to aragonite, 12%–18% of the Southern Ocean surface become undersaturated per decade, if emission reductions are delayed beyond 2015–2040. We conclude that the consequences of delaying global emission reductions are much better captured if the MDS of relevant Earth system variables is communicated in addition to current trends and total projected future changes.

Past and future evolution of Abies alba forests in Europe - comparison of a dynamic vegetation model with palaeo data and observations

Information on how species distributions and ecosystem services are impacted by anthropogenic climate change is important for adaptation planning. Palaeo data suggest that Abies alba formed forests under significantly warmer-than-present conditions in Europe and might be a native substitute for widespread drought-sensitive temperate and boreal tree species such as beech (Fagus sylvatica) and spruce (Picea abies) under future global warming conditions. Here, we combine pollen and macrofossil data, modern observations, and results from transient simulations with the LPX-Bern dynamic global vegetation model to assess past and future distributions of A. alba in Europe. LPX-Bern is forced with climate anomalies from a run over the past 21 000 years with the Community Earth System Model, modern climatology, and with 21st-century multimodel ensemble results for the high-emission RCP8.5 and the stringent mitigation RCP2.6 pathway. The simulated distribution for present climate encompasses the modern range of A. alba, with the model exceeding the present distribution in north-western and southern Europe. Mid-Holocene pollen data and model results agree for southern Europe, suggesting that at present, human impacts suppress the distribution in southern Europe. Pollen and model results both show range expansion starting during the Bølling–Allerød warm period, interrupted by the Younger Dryas cold, and resuming during the Holocene. The distribution of A. alba expands to the north-east in all future scenarios, whereas the potential (currently unrealized) range would be substantially reduced in southern Europe under RCP8.5. A. alba maintains its current range in central Europe despite competition by other thermophilous tree species. Our combined palaeoecological and model evidence suggest that A. alba may ensure important ecosystem services including stand and slope stability, infrastructure protection, and carbon sequestration under significantly warmer-than-present conditions in central Europe.

The stress response protein Gls24 is induced by copper and interacts with the CopZ copper chaperone of Enterococcus hirae

Intracellular copper routing in Enterococcus hirae is accomplished by the CopZ copper chaperone. Under copper stress, CopZ donates Cu(+) to the CopY repressor, thereby releasing its bound zinc and abolishing repressor-DNA interaction. This in turn induces the expression of the cop operon, which encodes CopY and CopZ, in addition to two copper ATPases, CopA and CopB. To gain further insight into the function of CopZ, the yeast two-hybrid system was used to screen for proteins interacting with the copper chaperone. This led to the identification of Gls24, a member of a family of stress response proteins. Gls24 is part of an operon containing eight genes. The operon was induced by a range of stress conditions, but most notably by copper. Gls24 was overexpressed and purified, and was shown by surface plasmon resonance analysis to also interact with CopZ in vitro. Circular dichroism measurements revealed that Gls24 is partially unstructured. The current findings establish a novel link between Gls24 and copper homeostasis.

Interaction of the receptor FGFRL1 with the negative regulator Spred1

FGFRL1 is a member of the fibroblast growth factor receptor family. It plays an essential role during branching morphogenesis of the metanephric kidneys, as mice with a targeted deletion of the Fgfrl1 gene show severe kidney dysplasia. Here we used the yeast two-hybrid system to demonstrate that FGFRL1 binds with its C-terminal, histidine-rich domain to Spred1 and to other proteins of the Sprouty/Spred family. Members of this family are known to act as negative regulators of the Ras/Raf/Erk signaling pathway. Truncation experiments further showed that FGFRL1 interacts with the SPR domain of Spred1, a domain that is shared by all members of the Sprouty/Spred family. The interaction could be verified by coprecipitation of the interaction partners from solution and by codistribution at the cell membrane of COS1 and HEK293 cells. Interestingly, Spred1 increased the retention time of FGFRL1 at the plasma membrane where the receptor might interact with ligands. FGFRL1 and members of the Sprouty/Spred family belong to the FGF synexpression group, which also includes FGF3, FGF8, Sef and Isthmin. It is conceivable that FGFRL1, Sef and some Sprouty/Spred proteins work in concert to control growth factor signaling during branching morphogenesis of the kidneys and other organs.

Variations of the Atlantic meridional overturning circulation in control and transient simulations of the last millennium

The variability of the Atlantic meridional overturing circulation (AMOC) strength is investigated in control experiments and in transient simulations of up to the last millennium using the low-resolution Community Climate System Model version 3. In the transient simulations the AMOC exhibits enhanced low-frequency variability that is mainly caused by infrequent transitions between two semi-stable circulation states which amount to a 10 percent change of the maximum overturning. One transition is also found in a control experiment, but the time-varying external forcing significantly increases the probability of the occurrence of such events though not having a direct, linear impact on the AMOC. The transition from a high to a low AMOC state starts with a reduction of the convection in the Labrador and Irminger Seas and goes along with a changed barotropic circulation of both gyres in the North Atlantic and a gradual strengthening of the convection in the Greenland-Iceland-Norwegian (GIN) Seas. In contrast, the transition from a weak to a strong overturning is induced by decreased mixing in the GIN Seas. As a consequence of the transition, regional sea surface temperature (SST) anomalies are found in the midlatitude North Atlantic and in the convection regions with an amplitude of up to 3 K. The atmospheric response to the SST forcing associated with the transition indicates a significant impact on the Scandinavian surface air temperature (SAT) in the order of 1 K. Thus, the changes of the ocean circulation make a major contribution to the Scandinavian SAT variability in the last millennium.

Zyxin interacts with the SH3 domains of the cytoskeletal proteins LIM-nebulette and Lasp-1

Zyxin is a versatile component of focal adhesions in eukaryotic cells. Here we describe a novel binding partner of zyxin, which we have named LIM-nebulette. LIM-nebulette is an alternative splice variant of the sarcomeric protein nebulette, which, in contrast to nebulette, is expressed in non-muscle cells. It displays a modular structure with an N-terminal LIM domain, three nebulin-like repeats, and a C-terminal SH3 domain and shows high similarity to another cytoskeletal protein, Lasp-1 (LIM and SH3 protein-1). Co-precipitation studies and results obtained with the two-hybrid system demonstrate that LIM-nebulette and Lasp-1 interact specifically with zyxin. Moreover, the SH3 domain from LIM-nebulette is both necessary and sufficient for zyxin binding. The SH3 domains from Lasp-1 and nebulin can also interact with zyxin, but the SH3 domains from more distantly related proteins such as vinexin and sorting nexin 9 do not. On the other hand, the binding site in zyxin is situated at the extreme N terminus as shown by site-directed mutagenesis. LIM-nebulette and Lasp-1 use the same linear binding motif. This motif shows some similarity to a class II binding site but does not contain the classical PXXP sequence. LIM-nebulette reveals a subcellular distribution at focal adhesions similar to Lasp-1. Thus, LIM-nebulette, Lasp-1, and zyxin may play an important role in the organization of focal adhesions.

Comparison of climate field reconstruction techniques: Application to Europe

This paper presents a comparison of principal component (PC) regression and regularized expectation maximization (RegEM) to reconstruct European summer and winter surface air temperature over the past millennium. Reconstruction is performed within a surrogate climate using the National Center for Atmospheric Research (NCAR) Climate System Model (CSM) 1.4 and the climate model ECHO-G 4, assuming different white and red noise scenarios to define the distortion of pseudoproxy series. We show how sensitivity tests lead to valuable “a priori” information that provides a basis for improving real world proxy reconstructions. Our results emphasize the need to carefully test and evaluate reconstruction techniques with respect to the temporal resolution and the spatial scale they are applied to. Furthermore, we demonstrate that uncertainties inherent to the predictand and predictor data have to be more rigorously taken into account. The comparison of the two statistical techniques, in the specific experimental setting presented here, indicates that more skilful results are achieved with RegEM as low frequency variability is better preserved. We further detect seasonal differences in reconstruction skill for the continental scale, as e.g. the target temperature average is more adequately reconstructed for summer than for winter. For the specific predictor network given in this paper, both techniques underestimate the target temperature variations to an increasing extent as more noise is added to the signal, albeit RegEM less than with PC regression. We conclude that climate field reconstruction techniques can be improved and need to be further optimized in future applications.