Probing the interactions of the solvated electron with DNA by molecular dynamics simulations: II. bromodeoxyuridine-thymidine mismatched DNA

The interaction of solvated electrons with DNA results in various types of DNA lesions. The in vitro and in vivo sensitisation of DNA to -induced damage is achieved by incorporation of the electron-affinity radiosensitiser bromodeoxyuridine (BUdR) in place of thymidine. However, in DNA duplexes containing single-stranded regions (bulged BUdR-DNA), the type of lesion is different and the efficiency of damage is enhanced. In particular, DNA interstrand crosslinks (ICL) form at high efficiency in bulged DNA but are not detectable in completely duplex DNA. Knowledge about the processes and interactions leading to these differences is obscure. Previously, we addressed the problem by applying molecular modelling and molecular dynamics (MD) simulations to a system of normal (BUdR·A)-DNA and a hydrated electron, where the excess electron was modelled as a localised eˉ(H₂O)₆ anionic cluster. The goal of the present study was to apply the same MD simulation to a wobble system, containing a pyrimidine–pyrimidine mismatched base pair, BUdR·T. The results show an overall dynamic pattern similar to that of the motion around normal DNA. However, the number of configuration states when was particularly close to DNA is different. Moreover, in the (BUdR·T)-wobble DNA system, the electron frequently approaches the brominated strand, including BUdR, which was not observed with the normal (BUdR·A)-DNA. The structure and exchange of water at the sites of immobilisation near DNA were also characterised. The structural dynamics of the wobble DNA is prone to more extensive perturbations, including frequent formation of cross-strand (cs) interatomic contacts. The structural deviations correlated with approaching DNA from the major groove side, with sodium ions trapped deep in the minor groove. Altogether, the obtained results confirm and/or throw light on dynamic-structure determinants possibly responsible for the enhanced radiation damage of wobble DNA. Figure The structure of the tightly bound single water-layer between the DNA and the electron (Site-8, five H₂O molecules, bold capped sticks); the rest of the “second” shell waters (lines, in atom type colour) surround the ˉ(H₂O)₆ cluster (yellow, space fill). Orange dashed lines H-bonds; only one of the five molecules from the single H₂O layer mediates a single-step H-bond bridge with N7(A8); the other four present a network of two(three)-step H-bond bridges between DNA/ partner atoms     

Topological research on diamagnetic susceptibilities of organic compounds

A novel molecular connectivity index, , based on the adjacency matrix of molecular graphs and novel atomic valence connectivities, , for predicting the molar diamagnetic susceptibilities of organic compounds is proposed. The is defined as: , where and E_i are the atomic valence connectivity and the valence orbital energy of atom i, respectively. A good QSPR model for molar diamagnetic susceptibilities can be constructed from and using multivariate linear regression (MLR). The correlation coefficient r, standard error, and average absolute deviation of the MLR model are 0.9918, 5.56 cgs, and 4.26 cgs, respectively, for the 721 organic compounds tested (training set). Cross-validation using the leave-one-out method demonstrates that the MLR model is highly reliable statistically. Using the MLR model, the average absolute deviations of the predicted values of molar diamagnetic susceptibility of another 360 organic compounds (test set) is 4.34 cgs. The results show that the current method is more effective than literature methods for estimating the molar diamagnetic susceptibility of an organic compound. The MLR method thus provides an acceptable model for the prediction of molar diamagnetic susceptibilities of organic compounds. Figure Plot of calculated vs experimental values of molar diamagnetic susceptibilities using the multivariate linear regression (MLR) model (Eq. 8)     

Effects of common soil anions and pH on the uptake and accumulation of perchlorate in lettuce

A mechanistic understanding of perchlorate () entry into plants is important for establishing the human health risk associated with consumption of contaminated produce and for assessing the effectiveness of phytoremediation. To determine whether common soil anions affect uptake and accumulation in higher plants, a series of competition experiments using lettuce (Lactuca sativa L.) were conducted between (50 nM) and (4–12 mM), (1–10 mM), or Cl⁻ (5–15 mM) in hydroponic solution. The effects of (0–5 mM) and pH (5.5–7.5) on uptake were also examined. Increasing in solution significantly reduced the amount of taken up by green leaf, butter head, and crisphead lettuces. Sulfate and Cl⁻ had no significant effects on uptake in lettuce over the concentrations tested. Increasing pH significantly reduced the amount of taken up by crisphead and green leaf lettuces, whereas increasing significantly reduced uptake in butter head lettuce. The inhibition by across all lettuce genotypes suggests that may share an ion carrier with , and the decrease in uptake with increasing pH or provides macroscopic evidence for cotransport across the plasma membrane.     

Theoretical study of AlH(2+)n (n=1-7) dications

The structures and stability of 1–7 dications were calculated at the ab initio MP2/aug-cc-pVTZ level of theory. The dications AlH²⁺ 1 and 2 were characterized to be unstable thermodynamically. However, these and the stable dications, 3–7 have considerable kinetic barriers for deprotonation. Each of the structures 3–7 contains one or more two-electron three-center (2e–3c) bonds. Aluminum atoms of these dications carry most of the positive charges, as indicated by NBO charge calculations.Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Long-term changes in nutrient concentrations of the Changjiang River and principal tributaries

We present long-term nutrient data on the Changjiang River (Yangtze River) at six hydrological stations and eight principal tributaries during the period 1958–1985. Three patterns of temporal changes were observed in nitrate and nitrite : minimal variations in the upper catchment area, rapid increases in the middle watershed towards the end of the 1970s, and a gradual increase in the lower drainage basin. Prior to the 1970s, the level of throughout the Changjiang River system remained fairly constant. In the 1980s, however, this changed, with the lowest values in the upper Changjiang changing rapidly to the highest in the middle reaches and then declining slowly but steadily in the lower courses. Compared to and ammonium and soluble reactive phosphorus (SRP) showed smaller increases or no long-term variations, while dissolved silica (DSi) concentration generally decreased at most stations. These three patterns of and changes in the Changjiang River system were reflective of the difference in chemical fertilizer use and landscape features (e.g., slope, soil type and water body area) of the drainage basins of the primary tributaries. The decreases in DSi were most likely attributed to a reduction in suspended sediment loading due to dam constructions and increasing diatom consumption. The increase in and with a reduction in DSi concentrations in the Changjiang River could have significant effects on the stoichiometric balance of nutrients delivered to the East China Sea and the ecosystem in this dynamic region.     

Anaerobic Ammonium Oxidation (Anammox) in Chesapeake Bay Sediments

Anaerobic ammonium oxidation (anammox) has recently been recognized as a pathway for the removal of fixed N from aquatic ecosystems. However, the quantitative significance of anammox in estuarine sediments is variable, and measurements have been limited to a few estuaries. We measured anammox and conventional denitrification activities in sediments along salinity gradients in the Chesapeake Bay and two of its sub-estuaries, the Choptank River and Patuxent River. Homogenized sediments were incubated with ^14/15N amendments of , , and to determine relative activities of anammox and denitrification. The percent of N₂ production due to anammox (ra%) ranged from 0 to 22% in the Chesapeake system, with the highest ra% in the freshwater portion of the main stem of upper Chesapeake Bay, where water column concentrations are consistently high. Intermediate levels of relative anammox (10%) were detected at locations corresponding to tidal freshwater and mesohaline locations in the Choptank River, whereas anammox was not detected in the tidal freshwater location in the Patuxent River. Anammox activity was also not detected in the seaward end of Chesapeake Bay, where water column concentrations are consistently low. The ra% did not correlate with accumulation rate in anoxic sediment incubations, but ra% was related to water column concentrations and salinity. Anammox bacterial communities were also examined by amplifying DNA extracted from the upper Chesapeake Bay sediment with polymerase chain reaction (PCR) primers that are specific for 16S rRNA genes of anammox organisms. A total of 35 anammox-like sequences were detected, and phylogenetic analysis grouped the sequences in two distinct clusters belonging to the Candidatus “Scalindua” genus.     

Regional Assessment of N Saturation using Foliar and Root $\varvec {\delta}^{\bf 15}{\bf N}$\varvec {\delta}^{\bf 15}{\bf N}

N saturation induced by atmospheric N deposition can have serious consequences for forest health in many regions. In order to evaluate whether foliar may be a robust, regional-scale measure of the onset of N saturation in forest ecosystems, we assembled a large dataset on atmospheric N deposition, foliar and root and N concentration, soil C:N, mineralization and nitrification. The dataset included sites in northeastern North America, Colorado, Alaska, southern Chile and Europe. Local drivers of N cycling (net nitrification and mineralization, and forest floor and soil C:N) were more closely coupled with foliar than the regional driver of N deposition. Foliar increased non-linearly with nitrification:mineralization ratio and decreased with forest floor C:N. Foliar was more strongly related to nitrification rates than was foliar N concentration, but concentration was more strongly correlated with N deposition. Root was more tightly coupled to forest floor properties than was foliar . We observed a pattern of decreasing foliar values across the following species: American beech>yellow birch>sugar maple. Other factors that affected foliar included species composition and climate. Relationships between foliar and soil variables were stronger when analyzed on a species by species basis than when many species were lumped. European sites showed distinct patterns of lower foliar , due to the importance of ammonium deposition in this region. Our results suggest that examining values of foliage may improve understanding of how forests respond to the cascading effects of N deposition.     

Determination of the chemical composition of ulvan, a cell wall polysaccharide from Ulva spp. (Ulvales, Chlorophyta) by FT-IR and chemometrics

The water-soluble cell wall polysaccharides from green seaweeds of Ulva spp. (Ulvales, Chlorophyta), referred to as ulvan, demonstrate composition- and structure-related functional properties. Mid-infrared spectroscopy combined with chemometric techniques was investigated as a means to rapidly predict the chemical composition of ulvan extracts. A calibration was realized with 41 ulvan extracts from two Ulva species. The variables studied included the constituent sugars (rhamnose, xylose, glucose, galactose, glucuronic acid, iduronic acid), protein, and sulfate contents. The correlation between Fourier transform infrared and chemical data was developed using partial least squares (PLS) regression with full cross-validation (leave one out). The coefficients of determination in cross-validation () and the standard error in cross-validation were determined for each variable. The PLS model validation resulted in a coefficient of determination () and a standard error in prediction. Good predictions were obtained for rhamnose ( = 0.9244), xylose ( = 0.8758), glucuronic acid ( = 0.9415), and sulfate ( = 0.9218), which are the main ulvan constituents. However, minor components such as proteins, glucose, galactose, and iduronic acid were not correctly predicted. This study showed that mid-infrared spectroscopy combined with PLS regression is a reliable and fast method for the quantification of the main chemical constituents of ulvan extracts.     

Factors affecting denitrification in agricultural headwater streams in Northeast Ohio,USA

As a result of increased anthropogenic nitrogen (N) loading in surface waters of agricultural watersheds, there is enhanced interest to understand and quantify N removal mechanisms. Denitrification, an important N removal mechanism in aquatic systems, may contribute to reducing N pollution in agricultural headwater streams. However, the key factors controlling this process in lotic systems remain unclear. The objective of our study was to examine the factors regulating rates of denitrification in the sediments of agricultural headwater streams in the mid-western USA. Denitrification rates were variable among streams and treatments (<0.1–28.0 μg N g AFDM⁻¹ h⁻¹) and on average, were higher than those reported for similar headwater streams. Carbon quantity and quality, and pH had no effect on denitrification, while temperature and nitrate ( ) concentrations had a positive effect on rates of denitrification. Specifically, controlled denitrification following Michaelis-Menten kinetics. We calculated a value of k_m (1.0 mg -N L^-1) that was comparable to other studies in aquatic sediments but was well below the median in-stream concentrations (5.2–17.4 mg -N L⁻¹) observed at the study sites. Despite high rates of denitrification, this removal mechanism is most likely saturated in the agricultural headwater streams we examined, suggesting that these systems are not effective at removing in-stream N. Handling editor: D. Ryder     

Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of various benzodiazepine analogues of γ-secretase inhibitors

A 3D QSAR analysis has been performed on a series of 67 benzodiazepine analogues reported as γ-secretase inhibitors using molecular field analysis (MFA), with G/PLS to predict steric and electrostatic molecular field interaction for the activity. The MFA study was carried out using a training set of 54 compounds. The predictive ability of model developed was assessed using a test set of 13 compounds ( as high as 0.729). The analyzed MFA model has demonstrated a good fit, having r² value of 0.858 and cross validated coefficient, value as 0.790. The analysis of the best MFA model provided insight into possible modification of the molecules for better activity.   Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A bias correction for estimates of effective population size based on linkage disequilibrium at unlinked gene loci*

Analysis of linkage disequilibrium (=mean squared correlation of allele frequencies at different gene loci) provides a means of estimating effective population size (N _e) from a single sample, but this method has seen much less use than the temporal method (which requires at least two samples). It is shown that for realistic numbers of loci and alleles, the linkage disequilibrium method can provide precision comparable to that of the temporal method. However, computer simulations show that estimates of N _e based on for unlinked, diallelic gene loci are sharply biased downwards ( in some cases) if sample size (S) is less than true N _e. The bias is shown to arise from inaccuracies in published formula for when S and/or N _e are small. Empirically derived modifications to for two mating systems (random mating and lifetime monogamy) effectively eliminate the bias (residual bias in % in most cases). The modified method also performs well in estimating N _e in non-ideal populations with skewed sex ratio or non-random variance in reproductive success. Recent population declines are not likely to seriously affect , but if N has recently increased from a bottleneck can be biased downwards for a few generations. These results should facilitate application of the disequilibrium method for estimating contemporary N _e in natural populations. However, a comprehensive assessment of performance of with highly polymorphic markers such as microsatellites is needed.The US Governmentȁ9s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

Regional Assessment of N Saturation using Foliar and Root \varvec {\delta}^{\bf 15}{\bf N}

N saturation induced by atmospheric N deposition can have serious consequences for forest health in many regions. In order to evaluate whether foliar may be a robust, regional-scale measure of the onset of N saturation in forest ecosystems, we assembled a large dataset on atmospheric N deposition, foliar and root and N concentration, soil C:N, mineralization and nitrification. The dataset included sites in northeastern North America, Colorado, Alaska, southern Chile and Europe. Local drivers of N cycling (net nitrification and mineralization, and forest floor and soil C:N) were more closely coupled with foliar than the regional driver of N deposition. Foliar increased non-linearly with nitrification:mineralization ratio and decreased with forest floor C:N. Foliar was more strongly related to nitrification rates than was foliar N concentration, but concentration was more strongly correlated with N deposition. Root was more tightly coupled to forest floor properties than was foliar . We observed a pattern of decreasing foliar values across the following species: American beech>yellow birch>sugar maple. Other factors that affected foliar included species composition and climate. Relationships between foliar and soil variables were stronger when analyzed on a species by species basis than when many species were lumped. European sites showed distinct patterns of lower foliar , due to the importance of ammonium deposition in this region. Our results suggest that examining values of foliage may improve understanding of how forests respond to the cascading effects of N deposition.     

Quantifying Lipari–Szabo modelfree parameters from 13CO NMR relaxation experiments

It is proposed to obtain effective Lipari–Szabo order parameters and local correlation times for relaxation vectors of protein ¹³CO nuclei by carrying out a ¹³CO-R₁ auto relaxation experiment, a transverse CSA/dipolar cross correlation and a transverse ¹³CO CSA/¹³CO–¹⁵N CSA/dipolar cross correlation experiment. Given the global rotational correlation time from ¹⁵N relaxation experiments, a new program COMFORD (CO-Modelfree Fitting Of Relaxation Data) is presented to fit the ¹³CO data to an effective order parameter , an effective local correlation time and the orientation of the CSA tensor with respect to the molecular frame. It is shown that the effective is least sensitive to rotational fluctuations about an imaginary axis and most sensitive to rotational fluctuations about an imaginary axis parallel to the NH bond direction. As such, the information is fully complementary to the ¹⁵N relaxation order parameter, which is least sensitive to fluctuations about the NH axis and most sensitive to fluctuations about the axis. The new paradigm is applied on data of Ca²⁺ saturated Calmodulin, and on available literature data for Ubiquitin. Our data indicate that the order parameters rapport on slower, and sometimes different, motions than the ¹⁵N relaxation order parameters. The CO local correlation times correlate well with the calmodulin’s secondary structure. Electronic Supplementary Material Supplementary material is available to authorized users in the online version of this article at .     

Use of isotopic signatures to assess the food web in a tropical shallow marine ecosystem of Southeastern Brazil

A dual isotope approach was used to assess the relative importance of terrestrial vegetation detritus and other primary producers in the trophic web of Flamengo Sound (Ubatuba, SP), SE Brazil, surrounded by the Atlantic Rain Forest. Primary producers showed distinct C signatures and the observed values suggest that little terrestrial or bulk sediment organic matter enter the food web of the sound. Suspended particulate organic matter (POM, supports the bulk of the consumers, with some contribution by macroalgae . Consumers C values ranged from −17.4 to . At least three trophic levels were detectable in the food web. The N value of POM was , while that of sediment and detritus was . The N values of suspension feeding benthic invertebrates were 8.2–, deposit feeders 8.3–, and carnivores 10.7–. Values for fishes were for detritivore, 11.4– for benthic feeders, 12.4– for zooplanktivores, and for piscivores/benthic invertebrate feeders. Squid mean value was . There is a reasonable agreement between feeding habits information from the literature and N values from this study. In the sound, the first and second trophic steps seem to be about 1– higher than those of similar organisms studied in temperate waters and this may reflect an input of allochtonous anthropogenic nitrogen enriched in ¹⁵N from human activities.     

Numerical treatment of two-center overlap integrals

Among the two-center integrals occurring in the molecular context, the two-center overlap integrals are numerous and difficult to evaluate to a level of high accuracy. The analytical and numerical difficulties arise mainly from the presence of the spherical Bessel integrals in the analytic expressions of these molecular integrals. Different approaches have been used to develop efficient algorithms for the numerical evaluation of the molecular integrals under consideration. These approaches are based on quadrature rules, Levin’s u transform, or the epsilon-algorithm of Wynn. In the present work, we use the nonlinear transformation of Sidi. This transformation is shown to be highly efficient in improving the convergence of highly oscillatory integrals, and it has been applied to molecular multicenter integrals, namely three-center attraction, hybrid, two-, three-, and four-center two-electron Coulomb and exchange integrals over B functions and over Slater-type functions. It is also been shown that when evaluating these molecular multicenter integrals the transformation is more efficient compared with the methods cited above. It is now proven that the integrand occurring in the analytic expression of the two-center overlap integrals satisfies all the conditions required to apply the transformation. A highly accurate algorithm based on this transformation is now developed. Special cases are presented and discussed for a better optimization of the algorithm. The numerical results section illustrates clearly the high efficiency of our algorithm.     

The mechanism and kinetics of decomposition of 5-aminotetrazole

The pathway and ab initio direct kinetics of the decomposition 5-aminotetrazole (5-ATZ) to HN₃ and NH₂CN was investigated. Reactant, products and transition state were optimized with MP2 and B3LYP methods using 6–311G** and aug-cc-pVDZ basis sets. The intrinsic reaction coordinate curve of the reaction was calculated using the MP2 method with 6–311G** basis set. The energies were refined using CCSD(T)/6–311G**. Rate constants were evaluated by conventional transition-state theory (CVT) and canonical variational transition-state theory (TST), with tunneling effect over 300 to 2,500 K. The results indicated that the tunneling effect and the variational effect are small for the calculated rate constants. The fitted three-parameter expression calculated using the CVT and TST methods are and , respectively. Figure The mechanism of the decomposition process of 5-ATZ to HN₃ and NH₂CN     

High groundwater nitrate concentrations inhibit eutrophication of sulphate-rich freshwater wetlands

During the last 60 years, pollution of the groundwater with has greatly increased in many parts of Europe, as a consequence of excessive use of manure and synthetic fertilisers. Monitoring of groundwater-fed wetlands indicated that sediments with high concentrations had the lowest Fe and concentrations in the pore water. A comparison of two restored open water fens, differing in supply via the groundwater, indicated that the redox potential and the sulphate ( ) reduction rate were lower when the groundwater contained not only but also high concentrations. The lower reduction rates in the -rich open water fen were associated with lower concentrations and the presence of plant species characteristic of clear water. In contrast, the higher reduction rates in the -poor open water fen were associated with very high concentrations and massive development of plant species characteristic of eutrophic environments. Investigations at -rich seepage sites in black alder carrs, showed that high concentrations in the pore water caused chlorosis in the alder carr vegetation, due to lower availability of Fe in the pore water and less Fe uptake by the plants. Experimental desiccation of sediments proved that the -rich seepage sites contained no oxidisable FeS _x , contrary to -poor locations, which became acidified and mobilised extremely high amounts of due to FeS _x oxidation. A laboratory experiment showed that addition to sediments led to reduced releases of Fe, and S^2–, very likely due to the oxidation of reduced Fe and S compounds. Overall, the results confirmed that is an energetically more favourable electron acceptor in anaerobic sediments than Fe and , and that high loads function as a redox buffer, preventing reduction of Fe and . Limited reduction prevents S^2– -mediated mobilisation of from Fe- complexes. At a higher redox potential, reduced Fe, including FeS _x , was oxidised, increasing the content of Fe(III) capable of binding . This prevented increased availability and the concomitant massive development of plant species characteristic of eutrophic environments.     

Inheritance of foliar stable carbon isotope discrimination and third-year height in Pinus taeda clones on contrasting sites in Florida and Georgia

Quantifying foliar stable carbon isotope discrimination (Δ) is a powerful approach for understanding genetic variation in gas exchange traits in large populations. The genetic architecture of Δ and third-year height is described for more than 1,000 clones of Pinus taeda tested on two contrasting sites. for Δ was 0.14 (±0.03), 0.20 (±0.07), and 0.09 (±0.04) at Florida, Georgia, and across sites, respectively. for stable carbon isotope discrimination ranged from 0.25 (±0.03) at the Florida site to 0.33 (±0.03) at the Georgia site, while the across-site estimate of was 0.19 (±0.02). For third-year height, ranged from 0.13 (±0.05) at the Georgia site to 0.20 (±0.06) at the Florida site with an across-site estimate of 0.09 (±0.05). Broad-sense heritability estimates for third-year height were 0.23 (±0.03), 0.28 (±0.03), and 0.13 (±0.02) at the Florida site, Georgia site, and across sites, respectively. Type B total genetic correlation for Δ was 0.70 ± 0.06, indicating that clonal rankings were relatively stable across sites, while for third-year height, rankings of clones were more unstable across the two trials . Third-year height and Δ were negatively correlated at the parental , full-sib family , and clonal levels, suggesting that genetic variation for Δ in P. taeda may be a result of differences in photosynthetic capacity. We conclude that Δ may be a useful selection trait to improve water-use efficiency and for guiding deployment decisions in P. taeda.     

Model-consistent estimation of the basic reproduction number from the incidence of an emerging infection

We investigate the merit of deriving an estimate of the basic reproduction number early in an outbreak of an (emerging) infection from estimates of the incidence and generation interval only. We compare such estimates of with estimates incorporating additional model assumptions, and determine the circumstances under which the different estimates are consistent. We show that one has to be careful when using observed exponential growth rates to derive an estimate of , and we quantify the discrepancies that arise.