Interactions Among Fuel Management,Species Composition,Bark Beetles,and Climate Change and the Potential Effects on Forests of the Lake Tahoe Basin

Climate-driven increases in wildfires, drought conditions, and insect outbreaks are critical threats to forest carbon stores. In particular, bark beetles are important disturbance agents although their long-term interactions with future climate change are poorly understood. Droughts and the associated moisture deficit contribute to the onset of bark beetle outbreaks although outbreak extent and severity is dependent upon the density of host trees, wildfire, and forest management. Our objective was to estimate the effects of climate change and bark beetle outbreaks on ecosystem carbon dynamics over the next century in a western US forest. Specifically, we hypothesized that (a) bark beetle outbreaks under climate change would reduce net ecosystem carbon balance (NECB) and increase uncertainty and (b) these effects could be ameliorated by fuels management. We also examined the specific tree species dynamics—competition and release—that determined NECB response to bark beetle outbreaks. Our study area was the Lake Tahoe Basin (LTB), CA and NV, USA, an area of diverse forest types encompassing steep elevation and climatic gradients and representative of mixed-conifer forests throughout the western United States. We simulated climate change, bark beetles, wildfire, and fuels management using a landscape-scale stochastic model of disturbance and succession. We simulated the period 2010–2100 using downscaled climate projections. Recurring droughts generated conditions conducive to large-scale outbreaks; the resulting large and sustained outbreaks significantly increased the probability of LTB forests becoming C sources over decadal time scales, with slower-than-anticipated landscape-scale recovery. Tree species composition was substantially altered with a reduction in functional redundancy and productivity. Results indicate heightened uncertainty due to the synergistic influences of climate change and interacting disturbances. Our results further indicate that current fuel management practices will not be effective at reducing landscape-scale outbreak mortality. Our results provide critical insights into the interaction of drivers (bark beetles, wildfire, fuel management) that increase the risk of C loss and shifting community composition if bark beetle outbreaks become more frequent.     

Molecular dynamics simulations of "loop closing" in the enzyme triose phosphate isomerase

We present molecular dynamics simulations on the active site region of dimeric triose phosphate isomerase (TIM) using the co-ordinates of native chicken muscle TIM as a starting point and performing simulations with no substrate, with dihydroxyacetone phosphate (DHAP), the natural substrate, and with dihydroxyacetone sulfate (DHAS), a substrate analog. Whereas most of the protein moves less than 1 A during the simulation, some residues in the active site loop move more than 8 A during the 10.5 picoseconds of dynamics for each of the simulations. Most interestingly, the nature of the loop motion depends on the substrate, with the largest motion found in the presence of DHAP, and only in the presence of DHAP does the loop move to "close off" the active site pocket. The final structure found for the DHAP-chicken TIM complex is qualitatively similar to that described by Alber et al. for DHAP-yeast TIM. Simulations on the monomeric protein gives insight into why the molecule is active only as a dimer.     

Effects of nucleotide bromination on the stabilities of Z-RNA and Z-DNA: a molecular mechanics/thermodynamic perturbation study

The structures of ZI- and ZII-form RNA and DNA oligonucleotides were energy minimized in vacuum using the AMBER molecular mechanics force field. Alternating C-G sequences were studied containing either unmodified nucleotides, 8-bromoguanosine in place of all guanosine residues, 5-bromocytidine in place of all cytidine residues, or all modified residues. Some molecules were also energy minimized in the presence of H2O and cations. Free energy perturbation calculations were done in which G8 and C5 hydrogen atoms in one or two residues of Z-form RNAs and DNAs were replaced in a stepwise manner by bromines. Bromination had little effect on the structures of the energy-minimized molecules. Both the minimized molecular energies and the results of the perturbation calculations indicate that bromination of guanosine at C8 will stabilize the Z forms of RNA and DNA relative to the nonbrominated Z form, while bromination of cytidine at C5 stabilizes Z-DNA and destabilizes Z-RNA. These results are in agreement with experimental data. The destabilizing effect of br5C in Z-RNAs is apparently due to an unfavorable interaction between the negatively charged C5 bromine atom and the guanosine hydroxyl group. The vacuum-minimized energies of the ZII-form oligonucleotides are lower than those of the corresponding ZI-form molecules for both RNA and DNA. Previous x-ray diffraction, nmr, and molecular mechanics studies indicate that hydration effects may favor the ZI conformation over the ZII form in DNA. Molecular mechanics calculations show that the ZII-ZI energy differences for the RNAs are greater than three times those obtained for the DNAs. This is due to structurally reinforcing hydrogen-bonding interactions involving the hydroxyl groups in the ZII form, especially between the guanosine hydroxyl hydrogen atom and the 3'-adjacent phosphate oxygen. In addition, the cytidine hydroxyl oxygen forms a hydrogen bond with the 5'-adjacent guanosine amino group in the ZII-form molecule. Both of these interactions are less likely in the ZI-form molecule: the former due to the orientation of the GpC phosphate away from the guanosine ribose in the ZI form, and the latter apparently due to competitive hydrogen bonding of the cytidine 2'-hydroxyl hydrogen with the cytosine carbonyl oxygen in the ZI form. The hydrogen-bonding interaction between the cytidine hydroxyl oxygen and the 5'-adjacent guanosine amino group in Z-RNA twists the amino group out of the plane of the base. This may be responsible for differences in the CD and Raman spectra of Z-RNA and Z-DNA.     

Conformations of deoxydinucleoside phosphates containing pyrimidine-pyrimidine photodimers

Conformational analysis of deoxydinucleoside monophosphates with the sequences TpT and CpC have been carried out with the incorporation of both cyclobutane type pyrimidine dimers and 6-4 photoadducts using the methods of molecular mechanics energy minimization. The effect of flexibility with respect to sugar geometries and glycosidic torsions have been studied and the relative energies of a large variety of structures have been compared. The salient features obtained from these calculations have been compared with the crystallographic and spectroscopic data on pyrimidine dimer incorporated deoxydinucleoside monophosphates. Effects of "inserting" the energetically favourable conformations of such structures into B-DNA helices have been discussed in terms of the distortions in helical structures.     

Clustering of hypervariable minisatellites in the proterminal regions of human autosomes

Six of the human minisatellites detected by DNA fingerprint probes have been localized by in situ hybridization to human metaphase chromosomes. These hypervariable loci are not dispersed at random in the human genome, but show preferential, though not exclusive, localization to terminal G-bands of human autosomes. Two of the proterminal minisatellites are very closely linked to other variable loci. Sequence analysis of one of these additional minisatellites suggests that the two linked minisatellites arose by independent amplification of different repeat units. The proterminal regions of human autosomes may therefore be rich in minisatellites, analogous to the pseudoautosomal terminal pairing region of human sex chromosomes that is similarly abundant in hypervariable minisatellites.     

A molecular mechanical study of netropsin-DNA interactions

We present molecular mechanical calculations on the complexes of netropsin with dA₆·dT₆, d(TATATA)₂, d(CGCGCG)₂, and d(CGCGAATTCGCG)₂. The complexes were model built using computer graphics and then completely energy refined. Our calculations are consistent with the observed AT preference for netropsin and suggest that mixed sugar pucker geometries should be more stable than uniform in netropsin complexes with poly[d(A-T)]·poly[d(A-T)] and poly(dA)·poly(dt). The netropsin·d(TATATA) and netropsin·dA₆·dT₆ complexes are significantly different in structure, leading to a possible reason why the observed thermodynamics of netropsin-association with poly[d(A-T)]·poly[d(A-T)] and with poly(dA)·poly(dT) are so different. We also model built and energy refined a structure of netropsin-d(CGCGAATTCGCG)₂ using as a guide the nmr data of Patel [(1982) Proc. Natl. Acad. Sci. USA, 79 , 6424–6428] and found a three-dimensional structure qualitatively consistent with the NOE enhancements observed by him. After our calculations were completed, we learned of an x-ray structure of a netropsin:d(CGCGAATTCGCG)₂ complex, and we compared the structure found in our calculation with the x-ray structure.     

Computer modeling of actinomycin D interactions with double-helical DNA

We have performed molecular mechanical calculations on intercalation complexes of actinomycin D with a series of base-paired hexanucleoside pentaphosphates; d(GCGCGC)2, d(GCCGGC)2, d(GCATGC)2, d(GCTAGC)2 and d(ATGCAT)2. Our results are in good agreement with previous experimental work on sequence selectivity. The results provide a rationalization for the strong preference of actinomycin D to intercalate on the 3' side of guanine residues, consistent with previously proposed models. Finally, the computed structures for d(ATGCAT)2-actinomycin D complexes have been compared with two-dimensional nuclear magnetic resonance nuclear Overhauser effect experimental results. To our knowledge, this is the first extensive comparison of molecular mechanical model structures for a drug-DNA complex with experimental solution phase data. We find generally good agreement between our computational models and the experimental solution phase structures.