Selecting high-dimensional mixed graphical models using minimal AIC or BIC forests

Background  

Chow and Liu showed that the maximum likelihood tree for multivariate discrete distributions may be found using a maximum weight spanning tree algorithm, for example Kruskal's algorithm. The efficiency of the algorithm makes it tractable for high-dimensional problems.     

Nuclear magnetic resonance studies of isolated structural domains of yeast phosphoglycerate kinase

The structural integrity and substrate binding properties of the two genetically engineered domains of yeast phosphoglycerate kinase were investigated using one- and two-dimensional nuclear magnetic resonance techniques. Both domains were found to fold with regions of native-like structure, with the N-domain showing greater conformational flexibility than the C-domain. The 'basic patch' region of the N-domain is, however, clearly perturbed by removal of the C-domain. This is most likely due to the absence of stabilizing interactions between the C-terminal peptide (including alpha-helices XIII and XIV) and the N-domain. The C-domain is able to bind nucleotide with an affinity only three times less than that of the native protein.     

Flexibility and folding of phosphoglycerate kinase

Flexibility and folding of phosphoglycerate kinase, a two-domain monomeric enzyme, have been studied using a wide variety of methods including theoretical approaches. Mutants of yeast phosphoglycerate kinase have been prepared in order to introduce cysteinyl residues as local probes throughout the molecule without perturbating significantly the structural or the functional properties of the enzyme. The apparent reactivity of a unique cysteine in each mutant has been used to study the flexibility of PGK. The regions of larger mobility have been found around residue 183 on segment beta F in the N-domain and residue 376 on helix XII in the C-domain. These regions are also parts of the molecule which unfold first. Ligand binding induces conformational motions in the molecule, especially in the regions located in the cleft. Moreover, the results obtained by introducing a fluorescent probe covalently linked to a cysteine are in agreement with the helix scissor motion of helices 7 and 14 assumed by Blake to direct the hinge bending motion of the domains during the catalytic cycle. The folding process of both horse muscle and yeast phosphoglycerate kinases involves intermediates. These intermediates are more stable in the horse muscle than in the yeast enzyme. In both enzymes, domains behave as structural modules capable of folding and stabilizing independently, but in the horse muscle enzyme the C-domain is more stable and refolds prior to the N-domain, contrary to that which has been observed in the yeast enzyme. A direct demonstration of the independence of domains in yeast phosphoglycerate kinase has been provided following the obtention of separated domains by site-directed mutagenesis. These domains have a native-like structure and refold spontaneously after denaturation by guanidine hydrochloride.     

Inactivation of 3 alpha-hydroxysteroid dehydrogenase by superoxide radicals. Modification of histidine and cysteine residues causes the conformational change.

3 alpha-Hydroxysteroid dehydrogenase (EC 1.1.1.50) from Pseudomonas testosterone was inactivated by superoxide radicals generated by the aerobic xanthine oxidase reaction. Superoxide dismutase, NAD+, bovine serum albumin and histidine and cysteine as free amino acids partially protected the enzyme from inactivation. NADH-binding properties were determined by fluorescence spectroscopy, and no variation was found between native enzyme and the unmodified fraction of the partly inactivated one. The fluorescence emission maximum for the completely inactivated enzyme was shifted 10 nm to a longer wavelength when compared with the native one, and it seems possible that the modification of histidine and cysteine residues by superoxide radicals causes the conformational change of the enzyme and the consequent loss of catalytic activity.     

Microbial transformation of ethylpyridines

Pyridine and its derivatives have been found as pollutants in the environment. Although alkylpyridines constitute the largest class of pyridines contaminating the environment, little information is available concerning the fate and transformation of these compounds. In this investigation ethylpyridines have been used as model compounds for investigating the biodegradability of alkylpyridines. A mixed culture of ethylpyridine-degrading microorganisms was obtained from a soil that had been exposed to a variety of pyridine derivatives for several decades. The enrichment culture was able to degrade 2-, 3-, and 4-ethylpyridine (100 mg/L) at 28° C and pH 7 within two weeks under aerobic conditions. The degradation rate was greatest for 2-ethylpyridine and least for 3-ethylpyridine. Transformation of ethylpyridines was dependent on substrate concentration, pH, and incubation temperature. Studies on the metabolic pathway of 4-ethylpyridine revealed two products; these chemicals were identified by MS and NMR analyses as 4-ethyl-2(1H)-pyridone and 4-ethyl-2-piperidone. 6-Ethyl-2(1H)-pyridone was determined to be a product of 2-ethylpyridine degradation. These results indicate that the transformation mechanism of ethylpyridines involves hydroxylation and reduction of the aromatic ring before ring cleavage.     

Unfolding-refolding of the domains in yeast phosphoglycerate kinase: comparison with the isolated engineered domains

The role of domains as folding units was investigated with a two-domain protein, yeast phosphoglycerate kinase. Each of the domains was produced independently by site-directed mutagenesis. It has been previously demonstrated by several criteria that these domains are able to fold in vivo into a quasi-native structure [Minard et al. (1989a) Protein Eng. 3, 55-60; Fairbrother et al. (1989) Protein Eng. 3, 5-11]. In the present study, the reversibility of the unfolding-refolding process induced by guanidine hydrochloride was investigated for the intact protein and the isolated domains. The transitions were followed by circular dichroism for both domains and the intact protein and by the variations in enzyme activity for the intact protein. Tryptophan residues were used as intrinsic conformational probes of the C-domain. An extrinsic fluorescent probe, N-[[(iodoacetyl)amino]ethyl]-8-naphthylamine-1-sulfonic acid (IAEDANS), was bound to the unique cysteinyl residue Cys97 to observe the conformational events in the N-domain. The unfolding-refolding transitions of each domain in the intact protein and in the isolated domains prepared by site-directed mutagenesis were compared. It was shown that the two domains are able to refold in a fully reversible process. A hyperfluorescent intermediate was detected during the folding of both the isolated C-domain and the intact yeast phosphoglycerate kinase. The stability of each isolated domain was found to be similar, the free energy of unfolding being approximately half that of the intact molecule.     

Alpha repeat proteins (αRep) as expression and crystallization helpers

We have previously described a highly diverse library of artificial repeat proteins based on thermostable HEAT-like repeats, named αRep. αReps binding specifically to proteins difficult to crystallize have been selected and in several examples, they made possible the crystallization of these proteins. To further simplify the production and crystallization experiments we have explored the production of chimeric proteins corresponding to covalent association between the targets and their specific binders strengthened by a linker. Although chimeric proteins with expression partners are classically used to enhance expression, these fusions cannot usually be used for crystallization. With specific expression partners like a cognate αRep this is no longer true, and chimeric proteins can be expressed purified and crystallized. αRep selection by phage display suppose that at least a small amount of the target protein should be produced to be used as a bait for selection and this might, in some cases, be difficult. We have therefore transferred the αRep library in a new construction adapted to selection by protein complementation assay (PCA). This new procedure allows to select specific binders by direct interaction with the target in the cytoplasm of the bacteria and consequently does not require preliminary purification of target protein. αRep binders selected by PCA or by phage display can be used to enhance expression, stability, solubility and crystallogenesis of proteins that are otherwise difficult to express, purify and/or crystallize.     

Murine immunization by cesium-137 irradiation attenuated Schistosoma mansoni cercariae

Cesium-137, becoming a more readily available ionizing gamma radiation source for laboratory use, was shown to effectively attenuate Schistosoma mansoni cercariae for vaccine production. In parallel comparison studies with the murine model, cesium-137 attenuated cercariae consistently afforded better (P greater than 0.05) protection than did the cobalt-60 prepared vaccine. Dose-response data indicated that the optimal total irradiation with cesium-137 was between 45 and 50 Krad.     

Enlarging a training set for genomic selection by imputation of un-genotyped animals in populations of varying genetic architecture

Background

The most common application of imputation is to infer genotypes of a high-density panel of markers on animals that are genotyped for a low-density panel. However, the increase in accuracy of genomic predictions resulting from an increase in the number of markers tends to reach a plateau beyond a certain density. Another application of imputation is to increase the size of the training set with un-genotyped animals. This strategy can be particularly successful when a set of closely related individuals are genotyped.

Methods

Imputation on completely un-genotyped dams was performed using known genotypes from the sire of each dam, one offspring and the offspring’s sire. Two methods were applied based on either allele or haplotype frequencies to infer genotypes at ambiguous loci. Results of these methods and of two available software packages were compared. Quality of imputation under different population structures was assessed. The impact of using imputed dams to enlarge training sets on the accuracy of genomic predictions was evaluated for different populations, heritabilities and sizes of training sets.

Results

Imputation accuracy ranged from 0.52 to 0.93 depending on the population structure and the method used. The method that used allele frequencies performed better than the method based on haplotype frequencies. Accuracy of imputation was higher for populations with higher levels of linkage disequilibrium and with larger proportions of markers with more extreme allele frequencies. Inclusion of imputed dams in the training set increased the accuracy of genomic predictions. Gains in accuracy ranged from close to zero to 37.14%, depending on the simulated scenario. Generally, the larger the accuracy already obtained with the genotyped training set, the lower the increase in accuracy achieved by adding imputed dams.

Conclusions

Whenever a reference population resembling the family configuration considered here is available, imputation can be used to achieve an extra increase in accuracy of genomic predictions by enlarging the training set with completely un-genotyped dams. This strategy was shown to be particularly useful for populations with lower levels of linkage disequilibrium, for genomic selection on traits with low heritability, and for species or breeds for which the size of the reference population is limited.