In-solution hybrid capture of bisulfite-converted DNA for targeted bisulfite sequencing of 174 ADME genes

DNA methylation is one of the most important epigenetic alterations involved in the control of gene expression. Bisulfite sequencing of genomic DNA is currently the only method to study DNA methylation patterns at single-nucleotide resolution. Hence, next-generation sequencing of bisulfite-converted DNA is the method of choice to investigate DNA methylation profiles at the genome-wide scale. Nevertheless, whole genome sequencing for analysis of human methylomes is expensive, and a method for targeted gene analysis would provide a good alternative in many cases where the primary interest is restricted to a set of genes.Here, we report the successful use of a custom Agilent SureSelect Target Enrichment system for the hybrid capture of bisulfite-converted DNA. We prepared bisulfite-converted next-generation sequencing libraries, which are enriched for the coding and regulatory regions of 174 ADME genes (i.e. genes involved in the metabolism and distribution of drugs). Sequencing of these libraries on Illumina’s HiSeq2000 revealed that the method allows a reliable quantification of methylation levels of CpG sites in the selected genes, and validation of the method using pyrosequencing and the Illumina 450K methylation BeadChips revealed good concordance.     

Using stable MutS dimers and tetramers to quantitatively analyze DNA mismatch recognition and sliding clamp formation

The process of DNA mismatch repair is initiated when MutS recognizes mismatched DNA bases and starts the repair cascade. The Escherichia coli MutS protein exists in an equilibrium between dimers and tetramers, which has compromised biophysical analysis. To uncouple these states, we have generated stable dimers and tetramers, respectively. These proteins allowed kinetic analysis of DNA recognition and structural analysis of the full-length protein by X-ray crystallography and small angle X-ray scattering. Our structural data reveal that the tetramerization domains are flexible with respect to the body of the protein, resulting in mostly extended structures. Tetrameric MutS has a slow dissociation from DNA, which can be due to occasional bending over and binding DNA in its two binding sites. In contrast, the dimer dissociation is faster, primarily dependent on a combination of the type of mismatch and the flanking sequence. In the presence of ATP, we could distinguish two kinetic groups: DNA sequences where MutS forms sliding clamps and those where sliding clamps are not formed efficiently. Interestingly, this inability to undergo a conformational change rather than mismatch affinity is correlated with mismatch repair.     

Allowance for Spatial Dispersion of Dielectric Permittivity in Polyelectrolyte Model of DNA

Abstract

In order to allow for real dielectric properties of a solvent in calculating of electrostatic characteristics of strongly charged polyions such as DNA in salt solution we consider a simple model of linear dielectric response of a medium. The interactions between charged particles are treated in the framework of self-consistent-field approximation. The basic characteristic of the problem, electrostatic potential, can be found from the solution of non-linear integro-differential equation. Specifically we consider so-called quasimacroscopic model where dielectric response of a medium depends only on the distance from the polyion. Application of the approach for calculating of the B-to-Z free energy qualitatively retains the main conclusion obtained previously within the model with fixed dielectric constant: non-monotonous behavior of the free energy difference as a function of ionic strength. At the same time, essential sensitivity of the results to specific values of dielectric parameters is observed.     

Computational biomechanics of a lumbar motion segment in pure and combined shear loads

Anterior shear has been implicated as a risk factor in spinal injuries. A 3D nonlinear poroelastic finite element model study of a lumbar motion segment L4-L5 was performed to predict the temporal shear response under various single and combined shear loads. Effects of nucleotomy and facetectomy as well as changes in the posture and facet gap distance were analyzed as well.     

Early stage of cytomegalovirus infection suppresses host microRNA expression regulation in human fibroblasts

Responses to human cytomegalovirus (HCMV) infection are largely individual and cell type specific. We investigated molecular profiles in 2 primary cell cultures of human fibroblasts, which are highly or marginally sensitive to HCMV infection, respectively. We screened expression of genes and microRNAs (miRs) at the early (3 hours) stage of infection. To assess molecular pathway activation profiles, we applied bioinformatic algorithms OncoFinder and MiRImpact. In both cell types, pathway regulation properties at mRNA and miR levels were markedly different. Surprisingly, in the infected highly sensitive cells, we observed a “freeze” of miR expression profiles compared to uninfected controls. Our results evidence that in the sensitive cells, HCMV blocks intracellular regulation of microRNA expression already at the earliest stage of infection. These data suggest somewhat new functions for HCMV products and demonstrate dependence of miR expression arrest on the host-encoded factors.     

Scale-dependence in geographic variation in a freshwater gastropod across the Palearctic

ABSTRACT

The reality of spatial clinal variation in morphological traits of freshwater pulmonate snails (Gastropoda: Pulmonata) has repeatedly been questioned or totally discounted. There is a lack of sound statistical evidence in the articles hitherto published on this subject supporting these claims. Here, by means of different analytical methods (analysis of spatial autocorrelation, linear regression analysis, canonical correlation analysis and others), we demonstrate that shell variation in the dwarf pond snail, Galba truncatula, is patterned in space throughout the northern and central Palearctic, with latitudinally-oriented clines in body size and in some shell proportions. Shell size in G. truncatula decreases with latitude and temperature, representing a special case of converse Bergmann cline. However, the temperature itself is hardly the main driver of shell size variation. It is argued that the shorter growing seasons at high latitudes may represent a better explanation for the observed trend. Shell proportions in the dwarf pond snails vary weakly at the macrogeographic scale, being spatially patterned at lower (mesogeographic) scales around 1200–1500?km. In general, spatial variation in G. truncatula shell size is decoupled from variation in shell shape, demonstrating clear scale-dependence similar to that found in different species of terrestrial (non-aquatic) pulmonate snails.     

Human MICAL1: Activation by the small GTPase Rab8 and small‐angle X‐ray scattering studies on the oligomerization state of MICAL1 and its complex with Rab8

Human MICAL1 is a member of a recently discovered family of multidomain proteins that couple a FAD‐containing monooxygenase‐like domain to typical protein interaction domains. Growing evidence implicates the NADPH oxidase reaction catalyzed by the flavoprotein domain in generation of hydrogen peroxide as a second messenger in an increasing number of cell types and as a specific modulator of actin filaments stability. Several proteins of the Rab families of small GTPases are emerging as regulators of MICAL activity by binding to its C‐terminal helical domain presumably shifting the equilibrium from the free – auto‐inhibited – conformation to the active one. We here extend the characterization of the MICAL1–Rab8 interaction and show that indeed Rab8, in the active GTP‐bound state, stabilizes the active MICAL1 conformation causing a specific four‐fold increase of k_cat of the NADPH oxidase reaction. Kinetic data and small‐angle X‐ray scattering (SAXS) measurements support the formation of a 1:1 complex between full‐length MICAL1 and Rab8 with an apparent dissociation constant of approximately 8 μM. This finding supports the hypothesis that Rab8 is a physiological regulator of MICAL1 activity and shows how the protein region preceding the C‐terminal Rab‐binding domain may mask one of the Rab‐binding sites detected with the isolated C‐terminal fragment. SAXS‐based modeling allowed us to propose the first model of the free full‐length MICAL1, which is consistent with an auto‐inhibited conformation in which the C‐terminal region prevents catalysis by interfering with the conformational changes that are predicted to occur during the catalytic cycle.     

Local polarity and hydrogen bonding inside the Sec14p phospholipid-binding cavity: high-field multi-frequency electron paramagnetic resonance studies

Sec14p promotes the energy-independent transfer of either phosphatidylinositol (PtdIns) or phosphatidylcholine (PtdCho) between lipid bilayers in vitro and represents the major PtdIns/PtdCho transfer protein in the budding yeast Saccharomyces cerevisiae. Herein, we employ multi-frequency high-field electron paramagnetic resonance (EPR) to analyze the electrostatic and hydrogen-bonding microenvironments for series of doxyl-labeled PtdCho molecules bound by Sec14p in a soluble protein-PtdCho complex. A structurally similar compound, 5-doxyl stearic acid dissolved in a series of solvents, was used for experimental calibration. The experiments yielded two-component rigid limit 130- and 220-GHz EPR spectra with excellent resolution in the gx region. Those components were assigned to hydrogen-bonded and nonhydrogen-bonded nitroxide species. Partially resolved 130-GHz EPR spectra from n-doxyl-PtdCho bound to Sec14p were analyzed using this two-component model and allowed quantification of two parameters. First, the fraction of hydrogen-bonded nitroxide species for each n-doxyl-PtdCho was calculated. Second, the proticity profile along the phospholipid-binding cavity of Sec14p was characterized. The data suggest the polarity gradient inside the Sec14p cavity is a significant contributor to the driving molecular forces for extracting a phospholipid from the bilayer. Finally, the enhanced g-factor resolution of EPR at 130 and 220 GHz provides researchers with a spectroscopic tool to deconvolute two major contributions to the x-component of the nitroxide g-matrix: hydrogen-bond formation and local electrostatic effects.