Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.

Short-Term Plasticity in a Computational Model of the Tail-Withdrawal Circuit in Aplysia.

The tail-withdrawal circuit of Aplysia provides a useful model system for investigating synaptic dynamics. Sensory neurons within the circuit manifest several forms of synaptic plasticity. Here, we developed a model of the circuit and investigated the ways in which depression (DEP) and potentiation (POT) contributed to information processing. DEP limited the amount of motor neuron activity that could be elicited by the monosynaptic pathway alone. POT within the monosynaptic pathway did not compensate for DEP. There was, however, a synergistic interaction between POT and the polysynaptic pathway. This synergism extended the dynamic range of the network, and the interplay between DEP and POT made the circuit responded preferentially to long-duration, low-frequency inputs.

Eluding antibiotic resistance: capitalizing on antimicrobial peptides interaction with the lipid bilayer

It is widely accepted that the emergence of drug-resistant pathogens is the result of the overuse and misuse of antibiotics. Infectious Disease Society of America, Center for Disease Control and World Health Organization continue to view, with concern, the lack of antibiotics in development, especially those against Gram-negative bacteria. Antimicrobial peptides (AMPs) have been proposed as an alternative to antibiotics due to their selective activity against microbes and minor ability to induce resistance. For example, the Food and Drug Administration approved Daptomycin (DAP) in 2003 for treatment of severe skin infections caused by susceptible Gram-positive organisms. Currently, there are 12 to 15 examples of modified natural and synthetic AMPs in clinical development. But most of these agents are against Gram-positive bacteria. Therefore, there is unmet medical need for antimicrobials used to treat infections caused by Gram-negative bacteria. In this study, we show that a pro-apoptotic peptide predominantly used in cancer therapy, (KLAKLAK)2, is an effective antimicrobial against Gram-negative laboratory strains and clinical isolates. Despite the therapeutic promise, AMPs development is hindered by their susceptibility to proteolysis. Here, we demonstrate that an all-D enantiomer of (KLAKLAK)2, resistant to proteolysis, retains its activity against Gram-negative pathogens. In addition, we have elucidated the specific site and mechanism of action of D(KLAKLAK)2 through a repertoire of whole-cell and membrane-model assays. Although it is considered that development of resistance does not represent an obstacle for AMPs clinical development, strains with decreased susceptibility to these compounds have been reported. Staphylococci resistance to DAP was observed soon after its approval for use and has been linked to alterations of the cell wall (CW) and cellular membrane (CM) properties. Immediately following staphylococcal resistance, Enterococci resistance to DAP was seen, yet the mechanism of resistance in enterococci remains unknown. Our findings demonstrate that, similar to S. aureus, development of DAP-resistance in a vancomycin-resistant E. faecalis isolate is associated with alterations of the CW and properties of the CM. However, the genes linked to these changes in enterococci appear to be different from those described in S. aureus.

A convolution model for energy transport in a therapeutic fast neutron beam

A three-dimensional model has been proposed that uses Monte Carlo and fast Fourier transform convolution techniques to calculate the dose distribution from a fast neutron beam. This method transports scattered neutrons and photons in the forward, lateral, and backward directions and protons, electrons, and positrons in the forward and lateral directions by convolving energy spread kernels with initial interaction available energy distributions. The primary neutron and photon spectrums have been derived from narrow beam attenuation measurements. The positions and strengths of the effective primary neutron, scattered neutron, and photon sources have been derived from dual ion chamber measurements. The size of the effective primary neutron source has been measured using a copper activation technique. Heterogeneous tissue calculations require a weighted sum of two convolutions for each component since the kernels must be invariant for FFT convolution. Comparisons between calculations and measurements were performed for several water and heterogeneous phantom geometries. ^

The interaction of P160 BCR with BCR-ABL gene products unique to Philadelphia chromosome leukemias

The BCR gene is involved in the pathogenesis of Philadelphia chromosome-positive (Ph$\sp1$) leukemias. Typically, the 5$\sp\prime$ portion of BCR on chromosome 22 becomes fused to a 5$\sp\prime$ truncated ABL gene from chromosome 9 resulting in a chimeric BCR-ABL gene. To investigate the role of the BCR gene product, a number of BCR peptide sequences were used to generate anti-BCR antibodies for detection of BCR and BCR-ABL proteins. Since both BCR and ABL proteins have kinase activity, the anti-BCR antibodies were tested for their ability to immunoprecipitate BCR and BCR-ABL proteins from cellular lysates by use of an immunokinase assay. Antisera directed towards the C-terminal portions of P160 BCR, sequences not present in BCR-ABL proteins, were capable of co-immunoprecipitating P210 BCR-ABL from the Ph$\sp1$- positive cell line K562. Re-immunoprecipitation studies following complete denaturation showed that C-terminal BCR antisera specifically recognized P160 BCR but not P210 BCR-ABL. These and other results indicated the presence of a P160 BCR/P210 BCR-ABL protein complex in K562 cells. Experiments performed with Ph$\sp1$-positive ALL cells and uncultured Ph$\sp1$-positive patient white blood cells established the general presence of BCR/BCR-ABL protein complexes in BCR-ABL expressing cells. However, two cell lines derived from Ph$\sp1$-positive patients lacked P160 BCR/P210 BCR-ABL complexes. Lysates from one of these cell lines mixed with lysates from a cell line that expresses only P160 BCR failed to generate BCR/BCR-ABL protein complexes in vitro indicating that P160 BCR and P210 BCR-ABL do not simply oligomerize.^ Two-dimensional tryptic maps were performed on both BCR and BCR-ABL proteins labeled in vitro with $\sp{32}$P. These maps indicate that the autophosphorylation sites in BCR-ABL proteins are primarily located within BCR exon 1 sequences in both P210 and P185 BCR-ABL, and that P160 BCR is phosphorylated in trans in similar sites by the activated ABL kinase of both BCR-ABL proteins. These results provide strong evidence that P160 BCR serves as a target for the BCR-ABL oncoprotein.^ K562 cells, induced to terminally differentiate with the tumor promoter TPA, show a loss of P210 BCR-ABL kinase activity 12-18 hours after addition of TPA. This loss coincides with the loss of activity in P160 BCR/P210 BCR-ABL complexes but not with the loss of the P210 BCR-ABL, suggesting the existence of an inactive form of P210 BCR-ABL. However, a degraded BCR-ABL protein served as the kinase active form preferentially sequestered within the remaining BCR/BCR-ABL protein complex.^ The results described in this thesis form the basis for a model for BCR-ABL induced leukemias which is presented and discussed. ^

Protein-protein interaction between phototaxis receptor sensory Rhodopsin I and its transducer HtrI

The molecular complex containing the seven transmembrane helix photoreceptor S&barbelow;ensory R&barbelow;hodopsin I&barbelow; (SRI) and transducer protein HtrI (H&barbelow;alobacterial Transducer for SRI&barbelow;) mediates color-sensitive phototaxis responses in the archaeon Halobacterium salinarum. Orange light causes an attractant response by a one-photon reaction and white light (orange + UV light) a repellent response by a two-photon reaction. Three aspects of SRI-HtrI structure/function and the signal transduction pathway were explored. First, the coupling of HtrI to the photoactive site of SRI was analyzed by mutagenesis and kinetic spectroscopy. Second, SRI-HtrI mutations and suppressors were selected and characterized to elucidate the color-sensing mechanism. Third, the signal relay through the transducer-bound histidine kinase was analyzed using an in vitro reconstitution system with known and newly identified taxis components. ^ Twenty-one mutations on HtrI were introduced by site-directed mutagenesis. Several replacements of charged residues perturbed the photochemical kinetics of SRI which led to the finding of a cluster of residues at the membrane/cytoplasm interface in HtrI electrostatically coupled to the photoactive site of SRI. We found by laser-flash kinetic spectroscopy that the transducer and these residues have specific effects on the light-induced proton transfer between the retinal chromophore and the protein. ^ One of the mutations showed an unusual mutant phenotype we called “inverted” signaling, in which the cell produces a repellent response to normally attractant light. Therefore, this mutant (E56Q of HtrI) had lost the color-discrimination by the SRI-HtrI complex. We used suppressor analysis to better understand the phenotype. Certain suppressors resulted in return of attractant responses to orange light but with inversion of the normally repellent response to white light to an attractant response. To explain this and other results, we formulated the Conformational Shuttling model in which the HtrI-SRI complex is poised in a metastable equilibrium of two conformations shifted in opposite directions by orange and white light. We tested this model by behavioral analysis (computerized cell tracking and motion study) of double mutants of inverting and suppressing mutations and the results confirmed the equilibrium-shift explanation. ^ We developed an in vitro system for measuring the effect of purified transducer on the histidine-kinase CheAH that controls the flagellar motor switch. The rate of kinase autophosphorylation was stimulated >2 fold in the reconstitution of the complete signal transduction system from purified components from H. salinarum. The in vitro assay also showed that the kinase activity was reduced in the absence and in the presence of high levels of linker protein CheWH. (Abstract shortened by UMI.) ^

A cardiovascular mathematical model of graded head-up tilt

A lumped parameter model of the cardiovascular system has been developed and optimized using experimental data obtained from 13 healthy subjects during graded head-up tilt (HUT) from the supine position to [Formula: see text]. The model includes descriptions of the left and right heart, direct ventricular interaction through the septum and pericardium, the systemic and pulmonary circulations, nonlinear pressure volume relationship of the lower body compartment, arterial and cardiopulmonary baroreceptors, as well as autoregulatory mechanisms. A number of important features, including the separate effects of arterial and cardiopulmonary baroreflexes, and autoregulation in the lower body, as well as diastolic ventricular interaction through the pericardium have been included and tested for their significance. Furthermore, the individual effect of parameter associated with heart failure, including LV and RV contractility, baseline systemic vascular resistance, pulmonary vascular resistance, total blood volume, LV diastolic stiffness and reflex gain on HUT response have also been investigated. Our fitted model compares favorably with our experimental measurements and published literature at a range of tilt angles, in terms of both global and regional hemodynamic variables. Compared to the normal condition, a simulated congestive heart failure condition produced a blunted response to HUT with regards to the percentage changes in cardiac output, stroke volume, end diastolic volume and effector response (i.e., heart contractility, venous unstressed volume, systemic vascular resistance and heart rate) with progressive tilting.

Biokinetics of nanoparticles and susceptibility to particulate exposure in a murine model of cystic fibrosis.

BACKGROUND Persons with cystic fibrosis (CF) are at-risk for health effects from ambient air pollution but little is known about the interaction of nanoparticles (NP) with CF lungs. Here we study the distribution of inhaled NP in a murine CF model and aim to reveal mechanisms contributing to adverse effects of inhaled particles in susceptible populations. METHODS Chloride channel defective CftrTgH (neoim) Hgu mice were used to analyze lung function, lung distribution and whole body biokinetics of inhaled NP, and inflammatory responses after intratracheal administration of NP. Distribution of 20-nm titanium dioxide NP in lungs was assessed on ultrathin sections immediately and 24 h after a one-hour NP inhalation. NP biokinetics was deduced from total and regional lung deposition and from whole body translocation of inhaled 30-nm iridium NP within 24 h after aerosol inhalation. Inflammatory responses were assessed within 7 days after carbon NP instillation. RESULTS Cftr mutant females had moderately reduced lung compliance and slightly increased airway resistance compared to wild type mice. We found no genotype dependent differences in total, regional and head deposition or in secondary-organ translocation of inhaled iridium NP. Titanium dioxide inhalation resulted in higher NP uptake by alveolar epithelial cells in Cftr mutants. Instillation of carbon NP induced a comparable acute and transient inflammatory response in both genotypes. The twofold increase of bronchoalveolar lavage (BAL) neutrophils in Cftr mutant compared to wild type mice at day 3 but not at days 1 and 7, indicated an impaired capacity in inflammation resolution in Cftr mutants. Concomitant to the delayed decline of neutrophils, BAL granulocyte-colony stimulating factor was augmented in Cftr mutant mice. Anti-inflammatory 15-hydroxyeicosatetraenoic acid was generally significantly lower in BAL of Cftr mutant than in wild type mice. CONCLUSIONS Despite lacking alterations in lung deposition and biokinetics of inhaled NP, and absence of significant differences in lung function, higher uptake of NP by alveolar epithelial cells and prolonged, acute inflammatory responses to NP exposure indicate a moderately increased susceptibility of lungs to adverse effects of inhaled NP in Cftr mutant mice and provides potential mechanisms for the increased susceptibility of CF patients to air pollution.

Development of a high- versus low-pathogenicity model of the free-living amoeba Naegleria fowleri

Species in the genus Naegleria are free-living amoebae of the soil and warm fresh water. Although around 30 species have been recognized, Naegleria fowleri is the only one that causes primary amoebic meningoencephalitis (PAM) in humans. PAM is an acute and fast progressing disease affecting the central nervous system. Most of the patients die within 1-2 weeks of exposure to the infectious water source. The fact that N. fowleri causes such fast progressing and highly lethal infections has opened many questions regarding the relevant pathogenicity factors of the amoeba. In order to investigate the pathogenesis of N. fowleri under defined experimental conditions, we developed a novel high- versus low-pathogenicity model for this pathogen. We showed that the composition of the axenic growth media influenced growth behaviour and morphology, as well as in vitro cytotoxicity and in vivo pathogenicity of N. fowleri. Trophozoites maintained in Nelson's medium were highly pathogenic for mice, demonstrated rapid in vitro proliferation, characteristic expression of surface membrane vesicles and a small cell diameter, and killed target mouse fibroblasts by both contact-dependent and -independent destruction. In contrast, N. fowleri cultured in PYNFH medium exhibited a low pathogenicity, slower growth, increased cell size and contact-dependent target cell destruction. However, cultivation of the amoeba in PYNFH medium supplemented with liver hydrolysate (LH) resulted in trophozoites that were highly pathogenic in mice, and demonstrated an intermediate proliferation rate in vitro, diminished cell diameter and contact-dependent target cell destruction. Thus, in this model, the presence of LH resulted in increased proliferation of trophozoites in vitro and enhanced pathogenicity of N. fowleri in mice. However, neither in vitro cytotoxicity mechanisms nor the presence of membrane vesicles on the surface correlated with the pathologic potential of the amoeba. This indicated that the pathogenicity of N. fowleri remains a complex interaction between as-yet-unidentified cellular mechanisms.

Air-Sea Interactions During the Passage of a Winter Storm Over the Gulf Stream: A Three-Dimensional Coupled Atmosphere-Ocean Model Study

A three-dimensional, regional coupled atmosphere-ocean model with full physics is developed to study air-sea interactions during winter storms off the U. S. east coast. Because of the scarcity of open ocean observations, models such as this offer valuable opportunities to investigate how oceanic forcing drives atmospheric circulation and vice versa. The study presented here considers conditions of strong atmospheric forcing (high wind speeds) and strong oceanic forcing (significant sea surface temperature (SST) gradients). A simulated atmospheric cyclone evolves in a manner consistent with Eta reanalysis, and the simulated air-sea heat and momentum exchanges strongly affect the circulations in both the atmosphere and the ocean. For the simulated cyclone of 19-20 January 1998, maximum ocean-to-atmosphere heat fluxes first appear over the Gulf Stream in the South Atlantic Bight, and this results in rapid deepening of the cyclone off the Carolina coast. As the cyclone moves eastward, the heat flux maximum shifts into the region near Cape Hatteras and later northeast of Hatteras, where it enhances the wind locally. The oceanic response to the atmospheric forcing is closely related to the wind direction. Southerly and southwesterly winds tend to strengthen surface currents in the Gulf Stream, whereas northeasterly winds weaken the surface currents in the Gulf Stream and generate southwestward flows on the shelf. The oceanic feedback to the atmosphere moderates the cyclone strength. Compared with a simulation in which the oceanic model always passes the initial SST to the atmospheric model, the coupled simulation in which the oceanic model passes the evolving SST to the atmospheric model produces higher ocean-to-atmosphere heat flux near Gulf Stream meander troughs. This is due to wind-driven lateral shifts of the stream, which in turn enhance the local northeasterly winds. Away from the Gulf Stream the coupled simulation produces surface winds that are 5 similar to 10% weaker. Differences in the surface ocean currents between these two experiments are significant on the shelf and in the open ocean.

Pentafluorophenyl phenyl interaction in biphenyl DNA

We prepared and investigated oligonucleotide duplexes of the sequence d(GATGAC(X)(n)GCTAG)d(CTAGC(Y)(n)GTCATC), in which X and Y designate biphenyl- (bph) and pentafluorobiphenyl- ((5F)bph) C-nucleotides, respectively, and n varies from 0-4. These hydrophobic base substitutes are expected to adopt a zipperlike, interstrand stacking motif, in which not only bph/bph or (5F)bph/(5F)bph homo pairs, but also (5F)bph/bph mixed pairs can be formed. By performing UV-melting curve analysis we found that incorporation of a single (5F)bph/(5F)bph pair leads to a duplex that is essentially as stable as the unmodified duplex (n=0), and 2.4 K more stable than the duplex with the nonfluorinated bph/bph pair. The T(m) of the mixed bph/(5F)bph pair was in between the T(m) values of the respective homo pairs. Additional, unnatural aromatic pairs increased the T(m) by +3.0-4.4 K/couple, irrespective of the nature of the aromatic residue. A thermodynamic analysis using isothermal titration calorimetry (ITC) of a series of duplexes with n=3 revealed lower (less negative) duplex formation enthalpies (DeltaH) in the (5F)bph/(5F)bph case than in the bph/bph case, and confirmed the higher thermodynamic stability (DeltaG) of the fluorinated duplex, suggesting it to be of entropic origin. Our data are compatible with a model in which the stacking of (5F)bph versus bph is dominated by dehydration of the aromatic units upon duplex formation. They do not support a model in which van der Waals dispersive forces (induced dipoles) or electrostatic (quadrupole) interactions play a dominant role

Polychlorinated Biphenyls in Glaciers. 2. Model Results of Deposition and Incorporation Processes

In previous work, Alpine glaciers have been identified as a secondary source of persistent organic pollutants (POPs). However, detailed understanding of the processes organic chemicals undergo in a glacial system was missing. Here, we present results from a chemical fate model describing deposition and incorporation of polychlorinated biphenyls (PCBs) into an Alpine glacier (Fiescherhorn, Switzerland) and an Arctic glacier (Lomonosovfonna, Norway). To understand PCB fate and dynamics, we investigate the interaction of deposition, sorption to ice and particles in the atmosphere and within the glacier, revolatilization, diffusion and degradation, and discuss the effects of these processes on the fate of individual PCB congeners. The model is able to reproduce measured absolute concentrations in the two glaciers for most PCB congeners. While the model generally predicts concentration profiles peaking in the 1970s, in the measurements, this behavior can only be seen for higher-chlorinated PCB congeners on Fiescherhorn glacier. We suspect seasonal melt processes are disturbing the concentration profiles of the lower-chlorinated PCB congeners. While a lower-chlorinated PCB congener is mainly deposited by dry deposition and almost completely revolatilized after deposition, a higher-chlorinated PCB congener is predominantly transferred to the glacier surface by wet deposition and then is incorporated into the glacier ice. The incorporated amounts of PCBs are higher on the Alpine glacier than on the Arctic glacier due to the higher precipitation rate and aerosol particle concentration on the former. Future studies should include the effects of seasonal melt processes, calculate the quantities of PCBs incorporated into the entire glacier surface, and estimate the quantity of chemicals released from glaciers to determine the importance of glaciers as a secondary source of organic chemicals to remote aquatic ecosystems.

Mechanisms of Ligand–Protein Interaction in Sec-14-like Transporters Investigated by Computer Simulations

We review our recent work on protein-ligand interactions in vitamin transporters of the Sec-14-like protein. Our studies focused on the cellular-retinaldehyde binding protein (CRALBP) and the alpha-tocopherol transfer protein (alpha-TTP). CRALBP is responsible for mobilisation and photo-protection of short-chain cis-retinoids in the dim-light visual cycle or rod photoreceptors. alpha-TTP is a key protein responsible for selection and retention of RRR-alpha-tocopherol, the most active isoform of vitamin E in superior animals. Our simulation studies evidence how subtle chemical variations in the substrate can lead to significant distortion in the structure of the complex, and how these changes can either lead to new protein function, or be used to model engineered protein variants with tailored properties. Finally, we show how integration of computational and experimental results can contribute in synergy to the understanding of fundamental processes at the biomolecular scale.

Loop formulation of the supersymmetric nonlinear O(N) sigma model

We derive the fermion loop formulation for the supersymmetric nonlinear O(N) sigma model by performing a hopping expansion using Wilson fermions. In this formulation the fermionic contribution to the partition function becomes a sum over all possible closed non-oriented fermion loop configurations. The interaction between the bosonic and fermionic degrees of freedom is encoded in the constraints arising from the supersymmetry and induces flavour changing fermion loops. For N ≥ 3 this leads to fermion loops which are no longer self-avoiding and hence to a potential sign problem. Since we use Wilson fermions the bare mass needs to be tuned to the chiral point. For N = 2 we determine the critical point and present boson and fermion masses in the critical regime.

Delay effects in the response of low-grade gliomas to radiotherapy: a mathematical model and its therapeutical implications

Low-grade gliomas (LGGs) are a group of primary brain tumours usually encountered in young patient populations. These tumours represent a difficult challenge because many patients survive a decade or more and may be at a higher risk for treatment-related complications. Specifically, radiation therapy is known to have a relevant effect on survival but in many cases it can be deferred to avoid side effects while maintaining its beneficial effect. However, a subset of LGGs manifests more aggressive clinical behaviour and requires earlier intervention. Moreover, the effectiveness of radiotherapy depends on the tumour characteristics. Recently Pallud et al. (2012. Neuro-Oncology, 14: , 1-10) studied patients with LGGs treated with radiation therapy as a first-line therapy and obtained the counterintuitive result that tumours with a fast response to the therapy had a worse prognosis than those responding late. In this paper, we construct a mathematical model describing the basic facts of glioma progression and response to radiotherapy. The model provides also an explanation to the observations of Pallud et al. Using the model, we propose radiation fractionation schemes that might be therapeutically useful by helping to evaluate tumour malignancy while at the same time reducing the toxicity associated to the treatment.