Performance of criteria for selecting evolutionary models in phylogenetics: a comprehensive study based on simulated datasets

Background  

Explicit evolutionary models are required in maximum-likelihood and Bayesian inference, the two methods that are overwhelmingly used in phylogenetic studies of DNA sequence data. Appropriate selection of nucleotide substitution models is important because the use of incorrect models can mislead phylogenetic inference. To better understand the performance of different model-selection criteria, we used 33,600 simulated data sets to analyse the accuracy, precision, dissimilarity, and biases of the hierarchical likelihood-ratio test, Akaike information criterion, Bayesian information criterion, and decision theory.     

Indirect readout of DNA sequence at the primary-kink site in the CAP-DNA complex: recognition of pyrimidine-purine and purine-purine steps

The catabolite activator protein (CAP) bends DNA in the CAP-DNA complex, typically introducing a sharp DNA kink, with a roll angle of approximately 40 degrees and a twist angle of approximately 20 degrees, between positions 6 and 7 of the DNA half-site, 5'-A1A2A3T4G5T6G7A8T9C10T11 -3' ("primary kink"). In previous work, we showed that CAP recognizes the nucleotide immediately 5' to the primary-kink site, T6, through an "indirect-readout" mechanism involving sequence effects on energetics of primary-kink formation. Here, to understand further this example of indirect readout, we have determined crystal structures of CAP-DNA complexes containing each possible nucleotide at position 6. The structures show that CAP can introduce a DNA kink at the primary-kink site with any nucleotide at position 6. The DNA kink is sharp with the consensus pyrimidine-purine step T6G7 and the non-consensus pyrimidine-purine step C6G7 (roll angles of approximately 42 degrees, twist angles of approximately 16 degrees ), but is much less sharp with the non-consensus purine-purine steps A6G7 and G6G7 (roll angles of approximately 20 degrees, twist angles of approximately 17 degrees). We infer that CAP discriminates between consensus and non-consensus pyrimidine-purine steps at positions 6-7 solely based on differences in the energetics of DNA deformation, but that CAP discriminates between the consensus pyrimidine-purine step and non-consensus purine-purine steps at positions 6-7 both based on differences in the energetics of DNA deformation and based on qualitative differences in DNA deformation. The structures further show that CAP can achieve a similar, approximately 46 degrees per DNA half-site, overall DNA bend through a sharp DNA kink, a less sharp DNA kink, or a smooth DNA bend. Analysis of these and other crystal structures of CAP-DNA complexes indicates that there is a large, approximately 28 degrees per DNA half-site, out-of-plane component of CAP-induced DNA bending in structures not constrained by end-to-end DNA lattice interactions and that lattice contacts involving CAP tend to involve residues in or near biologically functional surfaces.     

Engagement of ubiquitination and de-ubiquitination at rostral ventrolateral medulla in experimental brain death

Background

Whereas brain death is a vitally important clinical phenomenon, our contemporary understanding on its underlying cellular mechanisms remains elusive. This study evaluated whether the ubiquitin-proteasome system (UPS) in the rostral ventrolateral medulla (RVLM), a neural substrate that our laboratory identified previously to be intimately related to brain death, is engaged in this fatal process.

Methods

We performed proteomics, Western Blot, real-time PCR, ELISA and pharmacological experiments in conjunction with a clinically relevant experimental endotoxemia model of brain death based on intravenous administration of Escherichia coli lipopolysaccharide in adult male Sprague–Dawley rats.

Results

Proteomics, Western blot and enzyme activity analyses demonstrated that polyubiquitination was preserved and de-ubiquitination by ubiquitin C-terminal hydrolase isozyme-L1 (UCH-L1) was sustained, alongside increased monoubiquitin availability or proteasome activity in RVLM over the course of experimental endotoxemia. However, real-time PCR revealed no significant alteration in proteasome subunit alpha type-1, ubiquitin or UCH-L1 at mRNA level. Functionally, whereas microinjection into the bilateral RVLM of proteasome inhibitors (lactacystin or proteasome inhibitor II) potentiated survival, an inhibitor of ubiquitin-recycling (ubiquitin aldehyde) or an UCH-L1 inhibitor exacerbated mortality.

Conclusions

We proposed previously that the progression towards brain death entails a tug-of-war between pro-death and pro-life programs in RVLM. It is conceivable that ubiquitination or de-ubiquitination in RVLM participate in brain death by regulating the degradation of the proteins involved in those programs.     

Use of "loss-of-contact" substitutions to identify residues involved in an amino acid-base pair contact: effect of substitution of Gln18 of lac repressor by Gly, Ser, and Leu

A procedure to identify which base pair of lac operator (lacO) a suspected contacting amino acid of Lac repressor (LacR) interacts with is presented. The procedure is to eliminate the ability of the amino acid under study to contact DNA, and then to determine at which base pair--if any--specificity is eliminated. To implement this procedure, four sets of Escherichia coli K-12 strains have been constructed. These strains permit: (i) the substitution of a selected amino acid of LacR by, respectively, Gly, Ser, Leu, or Gln, and (ii) the analysis of the specificity of the resulting substituted LacR with respect to base pairs 5, 6, 7, 8, 9, and 10 of lacO. This procedure has been applied to Gln18 of LacR. The preliminary data indicate that LacR (Gln18----Gly) is unable to distinguish between the O+ base pair G:C and the Oc base pair T:A at position 7 of lacO (KDOc/KDO+ = 0.93). In contrast, LacR(Gln18----Gly) discriminates O+ from Oc by a factor of 13 to 23 at each other position. The same qualitative pattern of results was obtained with LacR(Gln18----Ser) and LacR (Gln18----Leu). Therefore, I propose that Gln18 contacts base pair 7 of lacO. This proposal is consistent with the contact predicted in Ebright, R. in Protein Structure, Folding, and Design. D. Oxender ed., Alan R. Liss, New York (1985), in press.     

Indirect readout of DNA sequence at the primary-kink site in the CAP-DNA complex: alteration of DNA binding specificity through alteration of DNA kinking.

The catabolite activator protein (CAP) sharply bends DNA in the CAP-DNA complex, introducing a DNA kink, with a roll angle of approximately 40 degrees and a twist angle of approximately 20 degrees, between positions 6 and 7 of the DNA half-site, 5'-A(1)A(2)A(3)T(4)G(5)T(6)G(7)A(8)T(9)C(10)T(11)-3' ("primary kink"). CAP recognizes the base-pair immediately 5' to the primary-kink site, T:A(6), through an "indirect-readout" mechanism involving sequence effects on the energetics of primary-kink formation. CAP recognizes the base-pair immediately 3' to the primary-kink site, G:C(7), through a "direct-readout" mechanism involving formation of a hydrogen bond between Glu181 of CAP and G:C(7). Here, we report that substitution of the carboxylate side-chain of Glu181 of CAP by the one-methylene-group-shorter carboxylate side-chain of Asp changes DNA binding specificity at position 6 of the DNA half site, changing specificity for T:A(6) to specificity for C:G(6), and we report a crystallographic analysis defining the structural basis of the change in specificity. The Glu181-->Asp substitution eliminates the primary kink and thus eliminates indirect-readout-based specificity for T:A(6). The Glu181-->Asp substitution does not eliminate hydrogen-bond formation with G:C(7), and thus does not eliminate direct-readout-based specificity for G:C(7).     

Upstream promoter sequences and alphaCTD mediate stable DNA wrapping within the RNA polymerase-promoter open complex

We show that the extent of stable DNA wrapping by Escherichia coli RNA polymerase (RNAP) in the RNAP-promoter open complex depends on the sequence of the promoter and, in particular, on the sequence of the upstream region of the promoter. We further show that the extent of stable DNA wrapping depends on the presence of the RNAP alpha-subunit carboxy-terminal domain and on the presence and length of the RNAP alpha-subunit interdomain linker. Our results indicate that the extensive stable DNA wrapping observed previously in the RNAP-promoter open complex at the lambda P(R) promoter is not a general feature of RNAP-promoter open complexes.     

HEN1 recognizes 21-24 nt small RNA duplexes and deposits a methyl group onto the 2' OH of the 3' terminal nucleotide

microRNAs (miRNAs) and small interfering RNAs (siRNAs) in plants bear a methyl group on the ribose of the 3′ terminal nucleotide. We showed previously that the methylation of miRNAs and siRNAs requires the protein HEN1 in vivo and that purified HEN1 protein methylates miRNA/miRNA* duplexes in vitro. In this study, we show that HEN1 methylates both miRNA/miRNA* and siRNA/siRNA* duplexes in vitro with a preference for 21–24 nt RNA duplexes with 2 nt overhangs. We also demonstrate that HEN1 deposits the methyl group on to the 2′ OH of the 3′ terminal nucleotide. Among various modifications that can occur on the ribose of the terminal nucleotide, such as 2′-deoxy, 3′-deoxy, 2′-O-methyl and 3′-O-methyl, only 2′-O-methyl on a small RNA inhibits the activity of yeast poly(A) polymerase (PAP). These findings indicate that HEN1 specifically methylates miRNAs and siRNAs and implicate the importance of the 2′-O-methyl group in the biology of RNA silencing.