Specific DNA recognition by the Antp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.     

The Influence of Spike Rate and Stimulus Duration on Noradrenergic Neurons

We model spiking neurons in locus coeruleus (LC), a brain nucleus involved in modulating cognitive performance, and compare with recent experimental data. Extracellular recordings from LC of monkeys performing target detection and selective attention tasks show varying responses dependent on stimuli and performance accuracy. From membrane voltage and ion channel equations, we derive a phase oscillator model for LC neurons. Average spiking probabilities of a pool of cells over many trials are then computed via a probability density formulation. These show that: (1) Post-stimulus response is elevated in populations with lower spike rates; (2) Responses decay exponentially due to noise and variable pre-stimulus spike rates; and (3) Shorter stimuli preferentially cause depressed post-activation spiking. These results allow us to propose mechanisms for the different LC responses observed across behavioral and task conditions, and to make explicit the role of baseline firing rates and the duration of task-related inputs in determining LC response.     

Vegetation of the Lago de Sanabria area (NW Iberia) since the end of the Pleistocene: a palaeoecological reconstruction on the basis of two new pollen sequences

Various pollen sequences from lacustrine deposits close to Lago de Sanabria (NW Iberia) have for several decades been a key source of information for palaeoenvironmental reconstructions of SW Europe, though their interpretation has been the subject of some controversy. Here we present two new pollen sequences obtained from this area, and a new palaeoenvironmental reconstruction of the region. The available pollen data reach back to before 18,000 b.p., a period of very harsh climate with seasonal (non continuous) sedimentation and a landscape characterised by herbaceous formations dominated by Gramineae and Artemisia, and scrub formations dominated by Ericaceae and Cistaceae. Subsequently sedimentation became continuous, and various regional forest expansions are apparent. At a local level, the first forest expansion began about 12,000 b.p., when Betula pollen reached 70% followed by peaks in Pinus sylvestris-type (>80%) and Quercus robur-type (40%). The Younger Dryas saw a retreat of woodland formations in the area around the lake, with broadleaved deciduous woodland (largely oak) retreating at mid and low altitudes, but with pine woodland persisting in more sheltered sites. The climatic improvement in the Early Holocene favoured re-expansion of woodland, dominated by Pinus sylvestris-type at higher and Quercus robur and Q. pyrenaica at lower altitudes, until anthropogenic deforestation commenced around 4,000 b.p. The disappearance of natural pine woodlands in this region is probably largely attributable to human interference.     

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

Dealing with Uncertainty in Spatially Explicit Population Models

It has been argued that spatially explicit population models (SEPMs) cannot provide reliable guidance for conservation biology because of the difficulty of obtaining direct estimates for their demographic and dispersal parameters and because of error propagation. We argue that appropriate model calibration procedures can access additional sources of information, compensating the lack of direct parameter estimates. Our objective is to show how model calibration using population-level data can facilitate the construction of SEPMs that produce reliable predictions for conservation even when direct parameter estimates are inadequate. We constructed a spatially explicit and individual-based population model for the dynamics of brown bears (Ursus arctos) after a reintroduction program in Austria. To calibrate the model we developed a procedure that compared the simulated population dynamics with distinct features of the known population dynamics (=patterns). This procedure detected model parameterizations that did not reproduce the known dynamics. Global sensitivity analysis of the uncalibrated model revealed high uncertainty in most model predictions due to large parameter uncertainties (coefficients of variation CV 0.8). However, the calibrated model yielded predictions with considerably reduced uncertainty (CV 0.2). A pattern or a combination of various patterns that embed information on the entire model dynamics can reduce the uncertainty in model predictions, and the application of different patterns with high information content yields the same model predictions. In contrast, a pattern that does not embed information on the entire population dynamics (e.g., bear observations taken from sub-areas of the study area) does not reduce uncertainty in model predictions. Because population-level data for defining (multiple) patterns are often available, our approach could be applied widely.     

Chaotic pairwise competition

We investigate a kind of competition possible in a system of at least three populations competing for the same limited resource. As a model we use generalised Volterra equations in which the growth rates and competition coefficients of populations depend on the number of members of all populations. Because of the nonconstant values of the last quantities the system could be repelled from the state of cyclic pairwise competition described by May and Leonard (SIAM J. Appl. Math. 29 (1975) 243.). We investigate the competition in a chaotic regime of evolution of the number of members of populations. We show that the nonconstant competition coefficients can lead to a regularisation of the time intervals of domination of each population and the non-constant growth rates can lead to decreasing length of the time intervals of domination as well as to chaotisation of the occurrence of these intervals. A quantity characterising the time intervals between the successive maxima of the number of the populations individuals is discussed. By means of the wavelet transform modulus maxima method we calculate the tau(q)-spectrum and the H?lder exponent for the time series of this quantity. The results of the theory are illustrated by an example of competition among the three main political parties in Bulgaria and we discuss qualitative aspects of the dynamics of change of preferences of voters.     

Long-term behaviour of 15N in an alpine grassland ecosystem

Nitrogen dynamics in semi-natural environments is crucial for the development and ecological stability of these systems. The present paper shows the results of the reinvestigation of a ¹⁵N-tracer experiment, which was established in the Grossglockner massif in Austria at 2300ma.s.l. in 1974/1975. We show that large quantities of nitrogen introduced by a single pulse labelling (amounting to approximately 1.7% of the nitrogen in the system) into an alpine grassland remain in the soil–plant system, with only 55% being lost during 27–28 years. In the first 10cm of the four investigated soil profiles 40% of ¹⁵N was recovered, being mainly bound in organic forms. A simple site specific model was established on the basis of the results considering a biological, residual and labile N-pool, the latter being the source for N-losses. By the model a long mean residence time close to 100 years was derived for the remaining ¹⁵N.     

Nitrogen processing in the hyporheic zone of a pastoral stream

The distribution of nitrogen-transforming processes, and factors controlling their rates, were determined within the hyporheic zone of a lowland stream draining agricultural land. In the field, physicochemical parameters were measured along a 10m-long hyporheic flow line between downwelling and upwelling zones. Sediment cores were retrieved from the stream bed surface, and from 20, 40 and 60cm deep in each zone, and in the laboratory, water from the corresponding depth was percolated through each core at the natural flow rate. Concentrations of nitrogen species and oxygen were measured before and after flow through each core. Denitrification was measured using a ¹⁵N-nitrate tracer. Shallow and downwelling zone samples were clearly distinct from deeper and upwelling zone samples in terms of physicochemical conditions, microbial processes and factors controlling nitrogen processing. Denitrification was highest in surface and downwelling zone cores, despite high oxygen levels, probably due to high pore-water nitrate concentrations in these cores and isolation of the denitrifying bacteria from oxygen in the bulk water by the hyporheic biofilms. Denitrification was limited by oxygen inhibition in the downwelling group, and by nitrate availability in the upwelling group. Strong evidence indicated that dissimilatory nitrate reduction to ammonium, occurred in almost all cores, and outcompeted denitrification for nitrate. In contrast, nitrification was undetectable in all but two cores, probably because of intense competition for oxygen. Field patterns and lab experiments indicated that the hyporheic zone at this moderately N-rich site is a strong sink for nitrate, fitting current theories that predict where hyporheic zones are nitrate sinks or nitrate sources.     

Significance of low-frequency local fluctuation motions in the transmembrane B and C α-helices of bacteriorhodopsin,to facilitate efficient proton uptake from the cytoplasmic surface,as revealed by site-directed solid-state <Superscript>13</Superscript>C NMR

¹³C NMR spectra of [1-¹³C]Val- or -Pro-labeled bacteriorhodopsin (bR) and its single or double mutants, including D85N, were recorded at various pH values to reveal conformation and dynamics changes in the transmembrane -helices, in relation to proton release and uptake between bR and the M-like state caused by modified charged states at Asp85 and the Schiff base (SB). It was found that the D85N mutant acquired local fluctuation motion with a frequency of 10⁴ Hz in the transmembrane B -helix, concomitant with deprotonation of SB in the M-like state at pH 10, as manifested from a suppressed ¹³C NMR signal of the [1-¹³C]-labeled Val49 residue. Nevertheless, local dynamics at Pro50 neighboring with Val49 turned out to be unchanged, irrespective of the charged state of SB as viewed from the ¹³C NMR of [1-¹³C]-labeled Pro50. This means that the transmembrane B -helix is able to acquire the fluctuation motion with a frequency of 10⁴ Hz beyond the kink at Pro50 in the cytoplasmic side. Concomitantly, fluctuation motion at the C helix with frequency in the order of 10⁴ Hz was found to be prominent, due to deprotonation of SB at pH 10, as viewed from the ¹³C NMR signal of Pro91. Accordingly, we have proposed here a novel mechanism as to proton uptake and transport based on a dynamic aspect that a transient environmental change from a hydrophobic to hydrophilic nature at Asp96 and SB is responsible for the reduced pK_a value which makes proton uptake efficient, as a result of acquisition of the fluctuation motion at the cytoplasmic side of the transmembrane B and C -helices in the M-like state. Further, it is demonstrated that the presence of a van der Waals contact of Val49 with Lys216 at the SB is essential to trigger this sort of dynamic change, as revealed from the ¹³C NMR data of the D85N/V49A mutant.     

Anti-actin antibodies generated against profilin:actin distinguish between non-filamentous and filamentous actin, and label cultured cells in a dotted pattern

Actin polymerization is a prominent feature of migrating cells, where it powers the protrusion of the leading edge. Many studies have characterized the well-ordered and dynamic arrangement of filamentous actin in this submembraneous space. However, less is known about the organization of unpolymerized actin. Previously, we reported on the use of covalently coupled profilin:actin to study actin dynamics and presented evidence that profilin-bound actin is a major source of actin for filament growth. To locate profilin:actin in the cell we have now used this non-dissociable complex for antibody generation, and obtained monospecific anti-actin and anti-profilin antibodies from two separate immunizations. Fluorescence microscopy revealed drastic differences in the staining pattern generated by the anti-actin antibody preparations. With one, distinct puncta appeared at the actin-rich leading edge and sometimes aligned with microtubules in the interior of the lamella, while the other displayed typical actin filament staining. Labelling experiments in vitro demonstrated failure of the first antibody to recognize filamentous actin and none of the two bound microtubules. The two anti-profilin antibodies purified in parallel generated a punctated pattern similar to that seen with the first anti-actin antibody. All antibody preparations labelled the nuclei.