Sterically controlled recognition of macromolecular sequence information by molecular tweezers

Sequence-specific binding is demonstrated between pyrene-based tweezer molecules and soluble, high molar mass copolyimides. The binding involves complementary pi - pi stacking interactions, polymer chain-folding, and hydrogen bonding and is extremely sensitive to the steric environment around the pyromellitimide binding-site. A detailed picture of the intermolecular interactions involved has been obtained through single-crystal X-ray studies of tweezer complexes with model diimides. Ring-current magnetic shielding of polyimide protons by the pyrene '' arms '' of the tweezer molecule induces large complexation shifts of the corresponding H-1 NMR resonances, enabling specific triplet sequences to be identified by their complexation shifts. Extended comonomer sequences (triplets of triplets in which the monomer residues differ only by the presence or absence of a methyl group) can be '' read '' by a mechanism which involves multiple binding of tweezer molecules to adjacent diimide residues within the copolymer chain. The adjacent-binding model for sequence recognition has been validated by two conceptually different sets of tweezer binding experiments. One approach compares sequence-recognition events for copolyimides having either restricted or unrestricted triple-triplet sequences, and the other makes use of copolymers containing both strongly binding and completely nonbinding diimide residues. In all cases the nature and relative proportions of triple-triplet sequences predicted by the adjacent-binding model are fully consistent with the observed H-1 NMR data.

Effect of headgroup size, charge, and solvent structure on polymer−micelle interactions, studied by molecular dynamics simulations

We performed atomistic molecular dynamics simulations of anionic and cationic micelles in the presence of poly(ethylene oxide) (PEO) to understand why nonionic water-soluble polymers such as PEO interact strongly with anionic micelles but only weakly with cationic micelles. Our micelles include sodium n-dodecyl sulfate (SDS), n-dodecyl trimethylammonium chloride (DTAC), n-dodecyl ammonium chloride (DAC), and micelles in which we artificially reverse the sign of partial charges in SDS and DTAC. We observe that the polymer interacts hydrophobically with anionic SDS but only weakly with cationic DTAC and DAC, in agreement with experiment. However, the polymer also interacts with the artificial anionic DTAC but fails to interact hydrophobically with the artificial cationic SDS, illustrating that large headgroup size does not explain the weak polymer interaction with cationic micelles. In addition, we observe through simulation that this preference for interaction with anionic micelles still exists in a dipolar "dumbbell" solvent, indicating that water structure and hydrogen bonding alone cannot explain this preferential interaction. Our simulations suggest that direct electrostatic interactions between the micelle and polymer explain the preference for interaction with anionic micelles, even though the polymer overall carries no net charge. This is possible given the asymmetric distribution of negative charges on smaller atoms and positive charges oil larger units in the polymer chain.

Infrared linear dichroism spectroscopy on amyloid fibrils aligned by molecular combing

We report the use of molecular combing as an alignment method to obtain macroscopically oriented amyloid fibrils on planar surfaces. The aligned fibrils are studied by polarized infrared spectroscopy. This gives structural information that cannot be definitively obtained from standard infrared experiments on isotropic samples, for example, confirmation of the characteristic cross-beta amyloid core structure, the side-chain orientation from specific amino acids, and the arrangement of the strands within the fibrils, as we demonstrate here. We employed amyloid fibrils from hen egg white lysozyme (HEWL) and from a model octapeptide. Our results demonstrate molecular combing as a straightforward method to align amyloid fibrils, producing highly anisotropic infrared linear dichroism (IRLD) spectra.

Sulphate transport by rat ileum. Effect of molybdate and other anions

Kinetic constants for SO42− transport by upper and lower rat ileum in vitro have been determined by computer fitting of rate vs concentration data obtained using the everted sac technique. MoO42− inhibition of this transport is competitive, and kinetic constants for the inhibition were similarly determined. Transport is also inhibited by the anions WO42−, S2O32− and SeO42−, in the order . These anions have no effect on the transport of l-valine. Low SO42− transport rates were observed in sacs from animals fed a high-molybdenum diet. The significance of the results with respect to the problem of molybdate toxicity in animals is discussed, and related to the known protective effect of SO42−.

Estrogen and thyroid hormone receptor interactions: physiological flexibility by molecular specificity

The influence of thyroid hormone on estrogen actions has been demonstrated both in vivo and in vitro. In transient transfection assays, the effects of liganded thyroid hormone receptors (TR) on transcriptional facilitation by estrogens bound to estrogen receptors (ER) display specificity according to the following: 1) ER isoform, 2) TR isoform, 3) the promoter through which transcriptional facilitation occurs, and 4) cell type. Some of these molecular phenomena may be related to thyroid hormone signaling of seasonal limitations upon reproduction. The various combinations of these molecular interactions provide multiple and flexible opportunities for relations between two major hormonal systems important for neuroendocrine feedbacks and reproductive behaviors.

The SV2A ligand levetiracetam inhibits presynaptic Ca2+ channels via an intracellular pathway

Levetiracetam (LEV) is a prominent antiepileptic drug (AED) which binds to neuronal synaptic vesicle glycoprotein 2A (SV2A) protein and has reported effects on ion channels, but retains a poorly-defined mechanism of action. Here, we investigate inhibition of voltage-dependent Ca2+ (CaV) channels as a potential mechanism by which LEV imparts effects on neuronal activity. We used electrophysiological methods to investigate the effects of LEV on cholinergic synaptic transmission and CaV channel activity in superior cervical ganglion neurons (SCGNs). In parallel, we investigated effects of the LEV ‘inactive’ R-enantiomer, UCB L060. Thus, LEV, but not UCB L060 (each 100 μM), inhibited synaptic transmission between SCGNs in long-term culture in a time-dependent manner, significantly reducing excitatory postsynaptic potentials (EPSP) following ≥30 min application. In isolated SCGNs, LEV pretreatment (≥1 h), but not acute (5 min) application, significantly inhibited whole-cell IBa amplitude. In current clamp recordings, LEV reduced the amplitude of the afterhyperpolarizing potential (AHP) in a Ca2+-dependent manner, but also increased action potential (AP) latency in a Ca2+-independent manner, suggesting further mechanisms associated with reduced excitability. Intracellular LEV application (4-5 min) caused a rapid inhibition of IBa amplitude to an extent comparable to that seen following extracellular LEV pretreatment ( ≥ 1 h). Neither pretreatment nor intracellular application of UCB L060 produced any inhibitory effects on IBa amplitude. These results identify a stereospecific intracellular pathway by which LEV inhibits presynaptic CaV channels; resultant reductions in neuronal excitability are proposed to contribute to the anticonvulsant effects of LEV.

Microbe-associated molecular pattern (MAMP) signatures, synergy, size and charge: influences on perception or mobility and host defence responses

Triggering of defences by microbes has mainly been investigated using single elicitors or microbe-associated molecular patterns (MAMPs), but MAMPs are released in planta as complex mixtures together with endogenous oligogalacturonan (OGA) elicitor. We investigated the early responses in Arabidopsis of calcium influx and oxidative burst induced by non-saturating concentrations of bacterial MAMPs, used singly and in combination: flagellin peptide (flg22), elongation factor peptide (elf18), peptidoglycan (PGN) and component muropeptides, lipo-oligosaccharide (LOS) and core oligosaccharides. This revealed that some MAMPs have additive (e.g. flg22 with elf18) and even synergistic (flg22 and LOS) effects, whereas others mutually interfere (flg22 with OGA). OGA suppression of flg22-induced defences was not a result of the interference with the binding of flg22 to its receptor flagellin-sensitive 2 (FLS2). MAMPs induce different calcium influx signatures, but these are concentration dependent and unlikely to explain the differential induction of defence genes [pathogenesis-related gene 1 (PR1), plant defensin gene 1.2 (PDF1.2) and phenylalanine ammonia lyase gene 1 (PAL1)] by flg22, elf18 and OGA. The peptide MAMPs are potent elicitors at subnanomolar levels, whereas PGN and LOS at high concentrations induce low and late host responses. This difference might be a result of the restricted access by plant cell walls of MAMPs to their putative cellular receptors. flg22 is restricted by ionic effects, yet rapidly permeates a cell wall matrix, whereas LOS, which forms supramolecular aggregates, is severely constrained, presumably by molecular sieving. Thus, MAMPs can interact with each other, whether directly or indirectly, and with the host wall matrix. These phenomena, which have not been considered in detail previously, are likely to influence the speed, magnitude, versatility and composition of plant defences.

Ribonucleocapsid formation of severe acute respiratory syndrome coronavirus through molecular action of the N-terminal domain of N protein

Conserved among all coronaviruses are four structural proteins: the matrix (M), small envelope (E), and spike (S) proteins that are embedded in the viral membrane and the nucleocapsid phosphoprotein (N), which exists in a ribonucleoprotein complex in the lumen. The N-terminal domain of coronaviral N proteins (N-NTD) provides a scaffold for RNA binding, while the C-terminal domain (N-CTD) mainly acts as oligomerization modules during assembly. The C terminus of the N protein anchors it to the viral membrane by associating with M protein. We characterized the structures of N-NTD from severe acute respiratory syndrome coronavirus (SARS-CoV) in two crystal forms, at 1.17 A (monoclinic) and at 1.85 A (cubic), respectively, resolved by molecular replacement using the homologous avian infectious bronchitis virus (IBV) structure. Flexible loops in the solution structure of SARS-CoV N-NTD are now shown to be well ordered around the beta-sheet core. The functionally important positively charged beta-hairpin protrudes out of the core, is oriented similarly to that in the IBV N-NTD, and is involved in crystal packing in the monoclinic form. In the cubic form, the monomers form trimeric units that stack in a helical array. Comparison of crystal packing of SARS-CoV and IBV N-NTDs suggests a common mode of RNA recognition, but they probably associate differently in vivo during the formation of the ribonucleoprotein complex. Electrostatic potential distribution on the surface of homology models of related coronaviral N-NTDs suggests that they use different modes of both RNA recognition and oligomeric assembly, perhaps explaining why their nucleocapsids have different morphologies.

Selectivity of bis(calix[4]diquinone) ionophores towards metal ions in solvent dimethylsulfoxide: a molecular mechanics and molecular dynamics study

Molecular modelling studies have been carried out on two bis(calix[4]diqu(inone) ionophores, each created from two (calix[4]diquinone)arenes bridged at their bottom rims via alkyl chains (CH2)(n), 1: n = 3, 2; n = 4, in order to understand the reported selectivity of these ligands towards different sized metal ions such as Na+, K+, Rb+, and Cs+ in dmso solution. Conformational. analyses have been carried out which show that in the lowest energy conformations of the two macrocycles, the individual calix[4]diquinones exhibit a combination of partial cone, 1,3-alternate and cone conformations. The interactions of these alkali metals with the macrocycles have been studied in the gas phase and in a periodic box of solvent dmso by molecular mechanics and molecular dynamics calculations. Molecular mechanics calculations have been carried out on the mode of entry of the ions into the macrocycles and suggest that this is likely to occur from the side of the central cavity, rather than through the main axis of the calix[4]diquinones. There are energy barriers of ca. 19 kcal mol(-1) for this entry path in the gas phase, but in solution no energy barrier is found. Molecular dynamics simulations show that in both 1 and 2, though particularly in the latter macrocycle, one or two solvent molecules are bonded to the metal throughout the course of the simulation, often to the exclusion, of one or more of the ether oxygen atoms. By contrast the carbonyl oxygen atoms remain bonded to the metal atoms throughout with bond lengths that remain significantly less than those to the ether oxygen atoms. Free energy perturbation studies have been carried out in dmso and indicate that for 1, the selectivity follows the order Rb+ approximate to K+ > Cs+ >> Na+, which is partially in agreement with the experimental results. The energy differences are small and indeed the ratio between stability constants found for Cs+ and K+ complexes is only 0.60, showing that 1 has only a slight preference for K+. For the larger receptor 2, which is better suited to metal complexation, the binding affinity follows the pattern Cs+ >> Rb+ >> K+ >> Na+, with energy differences of 5.75, 2.61, 2.78 kcal mol(-1) which is perfectly consistent with experimental results.

Bistability of the Atlantic overturning circulation in a global climate model and links to ocean freshwater transport

The possibility of a rapid collapse in the strength of the Atlantic meridional overturning circulation (AMOC), with associated impacts on climate, has long been recognized. The suggested basis for this risk is the existence of two stable regimes of the AMOC (‘on’ and ‘off’), and such bistable behaviour has been identified in a range of simplified climate models. However, up to now, no state-of-the-art atmosphere-ocean coupled global climate model (AOGCM) has exhibited such behaviour, leading to the interpretation that the AMOC is more stable than simpler models indicate. Here we demonstrate AMOC bistability in the response to freshwater perturbations in the FAMOUS AOGCM - the most complex AOGCM to exhibit such behaviour to date. The results also support recent suggestions that the direction of the net freshwater transport at the southern boundary of the Atlantic by the AMOC may be a useful physical indicator of the existence of bistability. We also present new estimates for this net freshwater transport by the AMOC from a range of ocean reanalyses which suggest that the Atlantic AMOC is currently in a bistable regime, although with large uncertainties. More accurate observational constraints, and an improved physical understanding of this quantity, could help narrow uncertainty in the future evolution of the AMOC and to assess the risk of a rapid AMOC collapse.

The influence of eddy parameterizations on the transport of the Antarctic Circumpolar Current in coupled climate models

The transport of the Antarctic Circumpolar Current (ACC) varies strongly across the coupled GCMs (general circulation models) used for the IPCC AR4. This note shows that a large fraction of this across-model variance can be explained by relating it to the parameterization of eddy-induced transports. In the majority of models this parameterization is based on the study by Gent and McWilliams (1990). The main parameter is the quasi-Stokes diffusivity kappa (often referred to less accurately as ’’thickness diffusion’’). The ACC transport and the meridional density gradient both correlate strongly with kappa across those models where kappa is a prescribed constant. In contrast, there is no correlation with the isopycnal diffusivity jiso across the models. The sensitivity of the ACC transport to kappa is larger than to the zonal wind stress maximum. Experiments with the fast GCM FAMOUS show that changing kappa directly affects the ACC transport by changing the density structure throughout the water column. Our results suggest that this limits the role of the wind stress magnitude in setting the ACC transport in FAMOUS. The sensitivities of the ACC and the meridional density gradient are very similar across the AR4 GCMs (for those models where kappa is a prescribed constant) and among the FAMOUS experiments. The strong sensitivity of the ACC transport to kappa needs careful assessment in climate models.

The vertical structure of ocean heat transport

One of the most important contributions the ocean makes to Earth's climate is through its poleward heat transport: about 1.5 PW or more than 30% of that accomplished by the ocean-atmosphere system (Trenberth and Caron, 2001). Recently, concern has arisen over whether global warming could affect this heat transport (Watson et al., 2001), for example, reducing high latitude convection and triggering a collapse of the deep overturning circulation (Rahmstorf, 1995). While the consequences of abrupt changes in oceanic circulation should be of concern, we argue that the attention devoted to deep circulations is disproportionate to their role in heat transport. For this purpose, we introduce a heat function which identifies the contribution to the heat transport by different components of the oceanic circulation. A new view of the ocean emerges in which a shallow surface intensified circulation dominates the poleward heat transport.

Notable differences between oxidized diruthenium complexes bridged by four isomeric diethynyl benzodithiophene ligands

Four new diruthenium complexes [{(η5-C5Me5)Ru(dppe)}2(μ-CuC–L–CuC)] featuring different bridging isomeric diethynyl benzodithiophenes viz. L = benzo[1,2-b;4,5-b’]dithiophene (complex 1), benzo[2,1-b;4,5b’]dithiophene (complex 2), benzo[1,2-b;3,4-b’]dithiophene (complex 3) and benzo[1,2-b;4,3-b’]-dithiophene (complex 4), were synthesized and characterized by molecular spectroscopic and crystallographicmethods. The subtle changes in the molecular structure introduced by the diethynyl benzodithiophene isomers have a notable impact on the stability of the oxidized complexes and their absorption characteristics in the visible-NIR and IR spectral domains. Electronic properties of stable oxidized complexes[1]n+ and [4]n+ (n = 1, 2) were investigated by cyclic voltammetry, UV-vis-NIR and IR spectroelectrochemistry as well as DFT and TDDFT calculations. The results document the largely bridgelocalized character of the oxidation of parents 1 and 4. Cations [2]+ and [3]+ are too unstable at ambient temperature to afford their unambiguous characterization. UV-vis-NIR absorption spectral data combined with TDDFT calculations (BLYP35) reveal that the broad electronic absorption of [1]+ and [4]+ in the NIR region has a mixed intraligand π–π* and MLCT character, with similar contribution from their spin-delocalized trans and cis conformers. A spin-localized (mixed-valence) rotamer was only observed for [1]+ at ambient temperature as a minor component on the time scale of IR spectroscopy.

An overview of Earth's global electric circuit and atmospheric conductivity

The Earth’s global atmospheric electric circuit depends on the upper and lower atmospheric boundaries formed by the ionosphere and the planetary surface. Thunderstorms and electrified rain clouds drive a DC current (∼1 kA) around the circuit, with the current carried by molecular cluster ions; lightning phenomena drive the AC global circuit. The Earth’s near-surface conductivity ranges from 10−7 S m−1 (for poorly conducting rocks) to 10−2 S m−1 (for clay or wet limestone), with a mean value of 3.2 S m−1 for the ocean. Air conductivity inside a thundercloud, and in fair weather regions, depends on location (especially geomagnetic latitude), aerosol pollution and height, and varies from ∼10−14 S m−1 just above the surface to 10−7 S m−1 in the ionosphere at ∼80 km altitude. Ionospheric conductivity is a tensor quantity due to the geomagnetic field, and is determined by parameters such as electron density and electron–neutral particle collision frequency. In the current source regions, point discharge (coronal) currents play an important role below electrified clouds; the solar wind-magnetosphere dynamo and the unipolar dynamo due to the terrestrial rotating dipole moment also apply atmospheric potential differences. Detailed measurements made near the Earth’s surface show that Ohm’s law relates the vertical electric field and current density to air conductivity. Stratospheric balloon measurements launched from Antarctica confirm that the downward current density is ∼1 pA m−2 under fair weather conditions. Fortuitously, a Solar Energetic Particle (SEP) event arrived at Earth during one such balloon flight, changing the observed atmospheric conductivity and electric fields markedly. Recent modelling considers lightning discharge effects on the ionosphere’s electric potential (∼+250 kV with respect to the Earth’s surface) and hence on the fair weather potential gradient (typically ∼130 V m−1 close to the Earth’s surface. We conclude that cloud-to-ground (CG) lightning discharges make only a small contribution to the ionospheric potential, and that sprites (namely, upward lightning above energetic thunderstorms) only affect the global circuit in a miniscule way. We also investigate the effects of mesoscale convective systems on the global circuit.

On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy

In this paper, the available potential energy (APE) framework of Winters et al. (J. Fluid Mech., vol. 289, 1995, p. 115) is extended to the fully compressible Navier– Stokes equations, with the aims of clarifying (i) the nature of the energy conversions taking place in turbulent thermally stratified fluids; and (ii) the role of surface buoyancy fluxes in the Munk & Wunsch (Deep-Sea Res., vol. 45, 1998, p. 1977) constraint on the mechanical energy sources of stirring required to maintain diapycnal mixing in the oceans. The new framework reveals that the observed turbulent rate of increase in the background gravitational potential energy GPEr , commonly thought to occur at the expense of the diffusively dissipated APE, actually occurs at the expense of internal energy, as in the laminar case. The APE dissipated by molecular diffusion, on the other hand, is found to be converted into internal energy (IE), similar to the viscously dissipated kinetic energy KE. Turbulent stirring, therefore, does not introduce a new APE/GPEr mechanical-to-mechanical energy conversion, but simply enhances the existing IE/GPEr conversion rate, in addition to enhancing the viscous dissipation and the entropy production rates. This, in turn, implies that molecular diffusion contributes to the dissipation of the available mechanical energy ME =APE +KE, along with viscous dissipation. This result has important implications for the interpretation of the concepts of mixing efficiency γmixing and flux Richardson number Rf , for which new physically based definitions are proposed and contrasted with previous definitions. The new framework allows for a more rigorous and general re-derivation from the first principles of Munk & Wunsch (1998, hereafter MW98)’s constraint, also valid for a non-Boussinesq ocean: G(KE) ≈ 1 − ξ Rf ξ Rf Wr, forcing = 1 + (1 − ξ )γmixing ξ γmixing Wr, forcing , where G(KE) is the work rate done by the mechanical forcing, Wr, forcing is the rate of loss of GPEr due to high-latitude cooling and ξ is a nonlinearity parameter such that ξ =1 for a linear equation of state (as considered by MW98), but ξ <1 otherwise. The most important result is that G(APE), the work rate done by the surface buoyancy fluxes, must be numerically as large as Wr, forcing and, therefore, as important as the mechanical forcing in stirring and driving the oceans. As a consequence, the overall mixing efficiency of the oceans is likely to be larger than the value γmixing =0.2 presently used, thereby possibly eliminating the apparent shortfall in mechanical stirring energy that results from using γmixing =0.2 in the above formula.