Mapping protein-DNA interactions using ChIP-sequencing

Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.

Metal-Nucleotide Interactions: Crystal Structures of Alkali (Li+, Na+, K+) and Alkaline Earth (Ca2+, Mg2+) Metal Complexes of Adenosine 2'-Monophosphate

Crystal structures of lithium, sodium, potassium, calcium and magnesium salts of adenosine 2'-monophosphate (2'-AMP) have been obtained at atomic resolution by X-ray crystallographic methods. 2'-AMP.Li belongs to the monoclinic space group P21 with a = 7.472(3)Å, b = 26.853(6) Å, c = 9.184(1)Å, b = 113.36(1)Å and Z= 4. 2'-AMP.Na and 2'-AMP.K crystallize in the trigonal space groups P31 and P3121 with a = 8.762(1)Å, c = 34.630(5)Å, Z= 6 and a = 8.931(4), Åc = 34.852(9)Å and Z= 6 respectively while 2'-AMP.Ca and 2'-AMP.Mg belong to space groups P6522 and P21 with cell parameters a = 9.487(2), c = 74.622(13), Z = 12 and a = 4.973(1), b = 10.023(2), c = 16.506(2), beta = 91.1(0) and Z = 2 respectively. All the structures were solved by direct methods and refined by full matrix least-squares to final R factors of 0.033, 0.028, 0.075, 0.069 and 0.030 for 2'-AMP.Li, 2'-AMP.Na, 2'- AMP.K, 2'-AMP.Ca and 2'-AMP.Mg, respectively. The neutral adenine bases in all the structures are in syn conformation stabilized by the O5'-N3 intramolecular hydrogen bond as in free acid and ammonium complex reported earlier. In striking contrast, the adenine base is in the anti geometry (cCN = -156.4(2)°) in 2'-AMP.Mg. Ribose moieties adopt C2'-endo puckering in 2'-AMP.Li and 2'-AMP.Ca, C2'-endo-C3'-exo twist puckering in 2'-AMP.Na and 2'-AMP.K and a C3'-endo-C2'-exo twist puckering in 2'-AMP.Mg structure. The conformation about the exocyclic C4'-C5' bond is the commonly observed gauche-gauche (g+) in all the structures except the gauche- trans (g-) conformation observed in 2'-AMP.Mg structure. Lithium ions coordinate with water, ribose and phosphate oxygens at distances 1.88 to 1.99Å. Na+ ions and K+ ions interact with phosphate and ribose oxygens directly and with N7 indirectly through a water oxygen. A distinct feature of 2'-AMP.Na and 2'-AMP.K structures is the involvement of ribose O4' in metal coordination. The calcium ion situated on a two-fold axis coordinates directly with three oxygens OW1, OW2 and O2 and their symmetry mates at distances 2.18 to 2.42Å forming an octahedron. A classic example of an exception to the existence of the O5'-N3 intramolecular hydorgen bond is the 2'-AMP.Mg strucure. Magnesium ion forms an octahedral coordination with three water and three phosphate oxygens at distances ranging from 2.02 to 2.11Å. A noteworthy feature of its coordination is the indirect link with N3 through OW3 oxygen resulting in macrochelation between the base and the phosphate group. Greater affnity of metal clays towards 5' compared to 2' and 3' nucleotides (J. Lawless, E. Edelson, and L. Manring, Am. Chem. Soc. Northwest Region Meeting, Seattle. 1978) due to macrochelation infered from solution studies (S. S. Massoud, H. Sigel, Eur. J. Biochem. 179, 451-458 (1989)) and interligand hydrogen bonding induced by metals postulated from metal-nucleotide structures in solid state (V. Swaminathan and M. Sundaralingam, CRC. Crit. Rev. Biochem. 6, 245-336 (1979)) are borne out by our structures also. The stacking patterns of adenine bases of both 2'-AMP.Na and 2'-AMP.K structures resemble the 2'-AMP.NH4 structure reported in the previous article. 2'-AMP.Li, 2'-AMP.Ca and 2'-AMP.Mg structures display base-ribose O4' stacking. An overview of interaction of monovalent and divalent cations with 2' and 5'-nucleotides has been presented.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

Chemical studies of viral entry mechanisms: I. Hydrophobic protein-lipid interactions during sendai virus membrane fusion. II. Kinetics of bacteriophage λ DNA injection.

Viruses possess very specific methods of targeting and entering cells. These methods would be extremely useful if they could also be applied to drug delivery, but little is known about the molecular mechanisms of the viral entry process. In order to gain further insight into mechanisms of viral entry, chemical and spectroscopic studies in two systems were conducted, examining hydrophobic protein-lipid interactions during Sendai virus membrane fusion, and the kinetics of bacteriophage λ DNA injection.

Sendai virus glycoprotein interactions with target membranes during the early stages of fusion were examined using time-resolved hydrophobic photoaffinity labeling with the lipid-soluble carbene generator3-(trifluoromethyl)-3-(m-^(125 )I] iodophenyl)diazirine (TID). The probe was incorporated in target membranes prior to virus addition and photolysis. During Sendai virus fusion with liposomes composed of cardiolipin (CL) or phosphatidylserine (PS), the viral fusion (F) protein is preferentially labeled at early time points, supporting the hypothesis that hydrophobic interaction of the fusion peptide at the N-terminus of the F_1 subunit with the target membrane is an initiating event in fusion. Correlation of the hydrophobic interactions with independently monitored fusion kinetics further supports this conclusion. Separation of proteins after labeling shows that the F_1 subunit, containing the putative hydrophobic fusion sequence, is exclusively labeled, and that the F_2 subunit does not participate in fusion. Labeling shows temperature and pH dependence consistent with a need for protein conformational mobility and fusion at neutral pH. Higher amounts of labeling during fusion with CL vesicles than during virus-PS vesicle fusion reflects membrane packing regulation of peptide insertion into target membranes. Labeling of the viral hemagglutinin/neuraminidase (HN) at low pH indicates that HN-mediated fusion is triggered by hydrophobic interactions, after titration of acidic amino acids. HN labeling under nonfusogenic conditions reveals that viral binding may involve hydrophobic as well as electrostatic interactions. Controls for diffusional labeling exclude a major contribution from this source. Labeling during reconstituted Sendai virus envelope-liposome fusion shows that functional reconstitution involves protein retention of the ability to undergo hydrophobic interactions.

Examination of Sendai virus fusion with erythrocyte membranes indicates that hydrophobic interactions also trigger fusion between biological membranes, and that HN binding may involve hydrophobic interactions as well. Labeling of the erythrocyte membranes revealed close membrane association of spectrin, which may play a role in regulating membrane fusion. The data show that hydrophobic fusion protein interaction with both artificial and biological membranes is a triggering event in fusion. Correlation of these results with earlier studies of membrane hydration and fusion kinetics provides a more detailed view of the mechanism of fusion.

The kinetics of DNA injection by bacteriophage λ. into liposomes bearing reconstituted receptors were measured using fluorescence spectroscopy. LamB, the bacteriophage receptor, was extracted from bacteria and reconstituted into liposomes by detergent removal dialysis. The DNA binding fluorophore ethidium bromide was encapsulated in the liposomes during dialysis. Enhanced fluorescence of ethidium bromide upon binding to injected DNA was monitored, and showed that injection is a rapid, one-step process. The bimolecular rate law, determined by the method of initial rates, revealed that injection occurs several times faster than indicated by earlier studies employing indirect assays.

It is hoped that these studies will increase the understanding of the mechanisms of virus entry into cells, and to facilitate the development of virus-mimetic drug delivery strategies.

Computational Inference of Genome-Wide Protein-DNA Interactions Using High-Throughput Genomic Data

Transcriptional regulation has been studied intensively in recent decades. One important aspect of this regulation is the interaction between regulatory proteins, such as transcription factors (TF) and nucleosomes, and the genome. Different high-throughput techniques have been invented to map these interactions genome-wide, including ChIP-based methods (ChIP-chip, ChIP-seq, etc.), nuclease digestion methods (DNase-seq, MNase-seq, etc.), and others. However, a single experimental technique often only provides partial and noisy information about the whole picture of protein-DNA interactions. Therefore, the overarching goal of this dissertation is to provide computational developments for jointly modeling different experimental datasets to achieve a holistic inference on the protein-DNA interaction landscape.

We first present a computational framework that can incorporate the protein binding information in MNase-seq data into a thermodynamic model of protein-DNA interaction. We use a correlation-based objective function to model the MNase-seq data and a Markov chain Monte Carlo method to maximize the function. Our results show that the inferred protein-DNA interaction landscape is concordant with the MNase-seq data and provides a mechanistic explanation for the experimentally collected MNase-seq fragments. Our framework is flexible and can easily incorporate other data sources. To demonstrate this flexibility, we use prior distributions to integrate experimentally measured protein concentrations.

We also study the ability of DNase-seq data to position nucleosomes. Traditionally, DNase-seq has only been widely used to identify DNase hypersensitive sites, which tend to be open chromatin regulatory regions devoid of nucleosomes. We reveal for the first time that DNase-seq datasets also contain substantial information about nucleosome translational positioning, and that existing DNase-seq data can be used to infer nucleosome positions with high accuracy. We develop a Bayes-factor-based nucleosome scoring method to position nucleosomes using DNase-seq data. Our approach utilizes several effective strategies to extract nucleosome positioning signals from the noisy DNase-seq data, including jointly modeling data points across the nucleosome body and explicitly modeling the quadratic and oscillatory DNase I digestion pattern on nucleosomes. We show that our DNase-seq-based nucleosome map is highly consistent with previous high-resolution maps. We also show that the oscillatory DNase I digestion pattern is useful in revealing the nucleosome rotational context around TF binding sites.

Finally, we present a state-space model (SSM) for jointly modeling different kinds of genomic data to provide an accurate view of the protein-DNA interaction landscape. We also provide an efficient expectation-maximization algorithm to learn model parameters from data. We first show in simulation studies that the SSM can effectively recover underlying true protein binding configurations. We then apply the SSM to model real genomic data (both DNase-seq and MNase-seq data). Through incrementally increasing the types of genomic data in the SSM, we show that different data types can contribute complementary information for the inference of protein binding landscape and that the most accurate inference comes from modeling all available datasets.

This dissertation provides a foundation for future research by taking a step toward the genome-wide inference of protein-DNA interaction landscape through data integration.

Probing protein−tannin interactions by isothermal titration microcalorimetry

Isothermal titration microcalorimetry (ITC) has been applied to investigate protein−tannin interactions. Two hydrolyzable tannins were studied, namely myrabolan and tara tannins, for their interaction with bovine serum albumin (BSA), a model globular protein, and gelatin, a model proline-rich random coil protein. Calorimetry data indicate that protein−tannin interaction mechanisms are dependent upon the nature of the protein involved. Tannins apparently interact nonspecifically with the globular BSA, leading to binding saturation at estimated tannin/BSA molar ratios of 48:1 for tara- and 178:1 for myrabolan tannins. Tannins bind to the random coil protein gelatin by a two-stage mechanism. The energetics of the first stage show evidence for cooperative binding of tannins to the protein, while the second stage indicates gradual saturation of binding sites as observed for interaction with BSA. The structure and flexibility of the tannins themselves alters the stoichiometry of the interaction, but does not appear to have any significant affect on the overall binding mechanism observed. This study demonstrates the potential of ITC for providing an insight into the nature of protein−tannin interactions.

Probing protein-tannin interactions by isothermal titration microcalorimetry

Isothermal titration microcalorimetry (ITC) has been applied to investigate protein-tannin interactions. Two hydrolyzable tannins were studied, namely myrabolan and tara tannins, for their interaction with bovine serum albumin (BSA), a model globular protein, and gelatin, a model proline-rich random coil protein. Calorimetry data indicate that protein-tannin interaction mechanisms are dependent upon the nature of the protein involved. Tannins apparently interact nonspecifically with the globular BSA, leading to binding saturation at estimated tannin/BSA molar ratios of 48:1 for tara- and 178:1 for myrabolan tannins. Tannins bind to the random coil protein gelatin by a two-stage mechanism. The energetics of the first stage show evidence for cooperative binding of tannins to the protein, while the second stage indicates gradual saturation of binding sites as observed for interaction with BSA. The structure and flexibility of the tannins themselves alters the stoichiometry of the interaction, but does not appear to have any significant affect on the overall binding mechanism observed. This study demonstrates the potential of ITC for providing an insight into the nature of protein-tannin interactions.

Characterization of protein-polyphenol interactions

Dietary polyphenols have received attention for their biologically significant functions as antioxidants, anticarcinogens or antimutagens, which have led to their recognition as potential nutraceuticals. Polyphenols also characteristically possess a significant binding affinity for proteins, which can lead to the formation of soluble and insoluble protein-polyphenol complexes. Questions remain concerning whether and to what extent the protein-polyphenol interaction influences functionality. For example, is the formation of protein-polyphenol complexes an obstacle to the nutritional bioavailability of either species? This article discusses the development of suitable methodologies to investigate the physicochemical basis of protein-polyphenol interactions and the influence of structure-activity relationships on binding affinities. (C) 2003 Elsevier Ltd. All rights reserved.