This work has been supported by NSFC Grant #11775016 and #11635002. JY has been supported by the CMCF of UCI via NSF DMS 1763272 and the Simons Foundation grant #594598 and start-up fund from UCI. LTD has been supported by Natural Science Foundation of Shanghai #20ZR1425400 and #21JC1403100. We also acknowledge the computational support from the Beijing Computational Science Research Center (CSRC).

The authors have no conflict of interests.

This work addresses how to conduct structure-based computational simulation and samplings to reveal a transcription factor or TF protein moving along DNA, not only at atomic detail of stepping, but also in the processive diffusion, which is essential for the facilitated diffusion of TF in the DNA target search. To do that, the Markov state model or MSM of a small TF domain protein WRKY stepping for 1-bp along homogeneous poly-A DNA was first constructed, so that an ensemble of protein conformations on the DNA along with ...

Transcription factors (TF) search for the target DNA to bind and regulate gene transcription and related activities1. Aside from the three-dimensional (3D) diffusion, the facilitated diffusion of TF has been suggested to be essential for target DNA search, in which the proteins can also slide or hop along one-dimensional (1D) DNA, or jump with intersegmental transfer on the DNA2,3,4,5,6,7.
In a recent study, we have conducted tens of microseconds (&#181;s) all-atom equilibrium molecular dynamics (MD) simulations on a plant TF - the WRKY domain protein on the DNA8. A complete 1-bp stepping of WRKY on poly-A DNA within microseconds has been captured. The movements of the protein along the DNA groove and hydrogen bonds (HBs) breaking-reforming dynamics have been observed. While such a trajectory represents one sampled path, an overall protein stepping landscape is still lack of. Here, we show how to expand computational samplings around the initially captured protein stepping path with the constructed Markov state model (MSM), which have been implemented widely for simulating a variety of biomolecular systems involving substantial conformational changes and time-scale separation9,10,11,12,13,14,15,16,17,18,19. The purpose is to reveal the conformational ensemble and meta-stable states of the TF protein diffusion along DNA for one cyclic step.
While the above MD simulation reveals atomic resolution of the protein movements for 1 bp on the DNA, the structural dynamics of long-time processive diffusion of the TF along DNA at the same high-resolution is hardly accessible. Conducting coarse-grained (CG) MD simulations at residue level is however technically approachable. The CG simulation time scale can be effectively extended to tens or hundreds of times longer than the atomic simulations20,21,22,23,24,25,26,27,28,29. Here, we show the CG simulations conducted by implementing the CafeMol software developed by Takada lab30.
In current protocol, we present the atomic simulations of the WRKY domain protein along poly-A DNA and the MSM construction first, which focus on sampling the protein stepping motions for only 1 bp along DNA. Then we present the CG modeling and simulations of the same protein-DNA system, which extend the computational sampling to the protein processive diffusion over tens of bps along DNA.
Here, we use GROMACS31,32,33 software to conduct MD simulations and MSMbuilder34 to construct the MSM for sampled conformational snapshots, as well as to use VMD35 to visualize the biomolecules. The protocol requires that the user to be able to install and implement the software above. The installation and implementation of the CafeMol30 software is then necessary for conducting the CG MD simulations. Further analyses of the trajectories and visualization are also conducted in VMD.

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>CafeMol</td>
			<td>Kyoto University</td>
			<td></td>
			<td>coarse-grained (CG) simulations</td>
		</tr>
		<tr>
			<td>GROMACS</td>
			<td>University of Groningen Royal Institute of Technology Uppsala University</td>
			<td></td>
			<td>molecular dynamics simulations software</td>
		</tr>
		<tr>
			<td>Matlab</td>
			<td>MathWorks</td>
			<td></td>
			<td>Numerical calculation software</td>
		</tr>
		<tr>
			<td>MSMbuilder</td>
			<td>Stanford University</td>
			<td></td>
			<td>build MSM</td>
		</tr>
		<tr>
			<td>VMD</td>
			<td>UNIVERSITY OF ILLINOIS AT URBANA-CHAMPAIGN</td>
			<td></td>
			<td>molecular visualization program</td>
		</tr>
	</tbody>
</table>

structure-based simulation and sampling of transcription factor protein movements along dna from atomic-scale stepping to coarse-grained diffusion

One-dimensional (1-D) sliding of transcription factor (TF) protein along DNA is essential for facilitated diffusion of the TF to locate target DNA site for genetic regulation. Detecting base-pair (bp) resolution of the TF sliding or stepping on the DNA is still experimentally challenging. We have recently performed all-atom molecular dynamics (MD) simulations capturing spontaneous 1-bp stepping of a small WRKY domain TF protein along DNA. Based on the 10 &#181;s WRKY stepping path obtained from such simulations, the protocol here shows how to conduct more extensive conformational samplings of the TF-DNA systems, by constructing the Markov state model (MSM) for the 1-bp protein stepping, with various numbers of micro- and macro-states tested for the MSM construction. In order to examine processive 1-D diffusional search of the TF protein along DNA with structural basis, the protocol further shows how to conduct coarse-grained (CG) MD simulations to sample long-time scale dynamics of the system. Such CG modeling and simulations are particularly useful to reveal the protein-DNA electrostatic impacts on the processive diffusional motions of the TF protein above tens of microseconds, in comparison with sub-microseconds to microseconds protein stepping motions revealed from the all-atom simulations.

Rotation-coupled sliding or 1 bp stepping of WRKY from the MSM construction 
All protein conformations on the DNA are mapped to the longitudinal movement X and rotation angle of the protein COM along DNA&#160;(see Figure 3A). The linear coupling of these two degrees indicates rotation-coupled stepping of the WRKY domain protein on the DNA. The conformations can be further clustered into 3 macrostates (S1, S2, and S3) in the MSM. The forward stepping of WRKY then follows...

Watch this Scientific Journal Video about Structure-Based Simulation and Sampling of Transcription Factor Protein Movements along DNA from Atomic-Scale Stepping to Coarse-Grained Diffusion at JoVE.com

Structure-Based Simulation and Sampling of Transcription Factor Protein Movements along DNA from Atomic-Scale Stepping to Coarse-Grained Diffusion

The goal of this protocol is to reveal structural dynamics of one-dimensional diffusion of protein along DNA, using a plant transcription factor WRKY domain protein as an exemplary system. To do this, both atomistic and coarse-grained molecular dynamics simulations along with extensive computational samplings have been implemented.

One-dimensional (1-D) sliding of transcription factor (TF) protein along DNA is essential for facilitated diffusion of the TF to locate target DNA site ...

structure-based-simulation-sampling-transcription-factor-protein

Research

JoVE Journal

Biology

3.0K Views.  Beijing Computational Science Research Center. The goal of this protocol is to reveal structural dynamics of one-dimensional diffusion of protein along DNA, using a plant transcription factor WRKY domain protein as an exemplary system. To do this, both atomistic and coarse-grained molecular dynamics simulations along with extensive computational samplings have been implemented.

Video: Structure-Based Simulation and Sampling of Transcription Factor Protein Movements along DNA from Atomic-Scale Stepping to Coarse-Grained Diffusion

Isolation of Genomic DNA from Mouse Tails

Preparing T Cell Growth Factor from Rat Splenocytes

Maintenance of antigen-specific T cell lines or clones in culture requires rounds of antigen-induced activation separated by phases of cell expansion 1,2. Addition of interleukin 2 to the culture media during the expansion phase is necessary to prevent cell death and sufficient to maintain short-term T cell lines but has been shown to increase Th1 polarization 3. Replacement of interleukin 2 by T cell growth factor (TCGF) which contains a mix of cytokines is more effective than interleukin 2 in maintaining long-term T cell lines in vitro 3. Moreover, TCGF can easily be prepared in large amounts in the laboratory and is much cheaper than recombinant interleukin 2. 
Here, we show how to prepare TCGF from rat splenocyte culture supernatants. For this procedure, we harvest spleens from naive Lewis rats euthanized for thymus and blood collection. We prepare single-cell suspensions from the spleens, lyze the red blood cells by osmotic shock, and seed the splenocytes in culture medium. The cells are stimulated with concanavalin A, a mitogen that non-selectively activates all rat T lymphocytes, inducing the production of cytokines. The culture supernantant is collected 48 hours later andexcess concanavalin A is bound to alpha methyl mannoside to prevent it from activating T cell lines to which TCGF will be added. The TCGF is then sterile-filtered, aliquoted, and stored at -20&deg;C.

We describe the preparation of T cell growth factor used for the in vitro expansion of antigen-specific rat T lymphocyte lines.

Maintenance of antigen-specific T cell lines or clones in culture requires rounds of antigen-induced activation separated by phases of cell expansion 1,2. ...

Application of Stopped-flow Kinetics Methods to Investigate the Mechanism of Action of a DNA Repair Protein

Transient kinetic analysis is indispensable for understanding the workings of biological macromolecules, since this approach yields mechanistic information including active site concentrations and intrinsic rate constants that govern macromolecular function. In case of enzymes, for example, transient or pre-steady state measurements identify and characterize individual events in the reaction pathway, whereas steady state measurements only yield overall catalytic efficiency and specificity. Individual events such as protein-protein or protein-ligand interactions and rate-limiting conformational changes often occur in the millisecond timescale, and can be measured directly by stopped-flow and chemical-quench flow methods. Given an optical signal such as fluorescence, stopped-flow serves as a powerful and accessible tool for monitoring reaction progress from substrate binding to product release and catalytic turnover1,2.
Here, we report application of stopped-flow kinetics to probe the mechanism of action of Msh2-Msh6, a eukaryotic DNA repair protein that recognizes base-pair mismatches and insertion/deletion loops in DNA and signals mismatch repair (MMR)3-5. In doing so, Msh2-Msh6 increases the accuracy of DNA replication by three orders of magnitude (error frequency decreases from ~10-6 to10-9 bases), and thus helps preserve genomic integrity. Not surprisingly, defective human Msh2-Msh6 function is associated with hereditary non-polyposis colon cancer and other sporadic cancers6-8. In order to understand the mechanism of action of this critical DNA metabolic protein, we are probing the dynamics of Msh2-Msh6 interaction with mismatched DNA as well as the ATPase activity that fuels its actions in MMR. DNA binding is measured by rapidly mixing Msh2-Msh6 with DNA containing a 2-aminopurine (2-Ap) fluorophore adjacent to a G:T mismatch and monitoring the resulting increase in 2-aminopurine fluorescence in real time. DNA dissociation is measured by mixing pre-formed Msh2-Msh6 G:T(2-Ap) mismatch complex with unlabeled trap DNA and monitoring decrease in fluorescence over time9. Pre-steady state ATPase kinetics are measured by the change in fluorescence of 7-diethylamino-3-((((2-maleimidyl)ethyl)amino)carbonyl) coumarin)-labeled Phosphate Binding Protein (MDCC-PBP) on binding phosphate (Pi) released by Msh2-Msh6 following ATP hydrolysis9,10.
The data reveal rapid binding of Msh2-Msh6 to a G:T mismatch and formation of a long-lived Msh2-Msh6 G:T complex, which in turn results in suppression of ATP hydrolysis and stabilization of the protein in an ATP-bound form. The reaction kinetics provide clear support for the hypothesis that ATP-bound Msh2-Msh6 signals DNA repair on binding a mismatched base pair in the double helix.
F. Noah Biro and Jie Zhai contributed to this paper equally.

Msh2-Msh6 is responsible for initiating repair of replication errors in DNA. Here we present a transient kinetics approach towards understanding how this critical protein works. The report illustrates stopped-flow experiments for measuring the coupled DNA binding and ATPase kinetics underlying Msh2-Msh6 mechanism of action in DNA repair.

Transient kinetic analysis is indispensable for understanding the workings of biological macromolecules, since this approach yields mechanistic ...

A Noninvasive Hair Sampling Technique to Obtain High Quality DNA from Elusive Small Mammals

Noninvasive genetic sampling approaches are becoming increasingly important to study wildlife populations. A number of studies have reported using noninvasive sampling techniques to investigate population genetics and demography of wild populations1. This approach has proven to be especially useful when dealing with rare or elusive species2. While a number of these methods have been developed to sample hair, feces and other biological material from carnivores and medium-sized mammals, they have largely remained untested in elusive small mammals. In this video, we present a novel, inexpensive and noninvasive hair snare targeted at an elusive small mammal, the American pika (Ochotona princeps). We describe the general set-up of the hair snare, which consists of strips of packing tape arranged in a web-like fashion and placed along travelling routes in the pikas&#x2019; habitat. We illustrate the efficiency of the snare at collecting a large quantity of hair that can then be collected and brought back to the lab. We then demonstrate the use of the DNA IQ system (Promega) to isolate DNA and showcase the utility of this method to amplify commonly used molecular markers including nuclear microsatellites, amplified fragment length polymorphisms (AFLPs), mitochondrial sequences (800bp) as well as a molecular sexing marker. Overall, we demonstrate the utility of this novel noninvasive hair snare as a sampling technique for wildlife population biologists. We anticipate that this approach will be applicable to a variety of small mammals, opening up areas of investigation within natural populations, while minimizing impact to study organisms.

We present a noninvasive sampling approach to efficiently collect hair samples from elusive small mammals, as shown for the American pika. We demonstrate the utility of this method by extracting DNA from sampled hair and amplifying several types of molecular markers commonly used in studies of wildlife ecology and conservation.

Noninvasive genetic sampling approaches are becoming increasingly important to study wildlife populations. A number of studies have reported using ...

A Practical and Novel Method to Extract Genomic DNA from Blood Collection Kits for Plasma Protein Preservation

Laboratory tests can be done on the cellular or fluid portions of the blood. The use of different blood collection tubes determines the portion of the blood that can be analyzed (whole blood, plasma or serum). Laboratories involved in studying the genetic basis of human disorders rely on anticoagulated whole blood collected in EDTA-containing vacutainer as the source of DNA for genetic / genomic analysis. Because most clinical laboratories perform biochemical, serologic and viral testing as a first step in phenotypic outcome investigation, anticoagulated blood is also collected in heparin-containing tube (plasma tube). Therefore when DNA and plasma are needed for simultaneous and parallel analyses of both genomic and proteomic data, it is customary to collect blood in both EDTA and heparin tubes. If blood could be collected in a single tube and serve as a source for both plasma and DNA, that method would be considered an advancement to existing methods. The use of the compacted blood after plasma extraction represents an alternative source for genomic DNA, thus minimizing the amount of blood samples processed and reducing the number of samples required from each patient. This would ultimately save time and resources.
The BD P100 blood collection system for plasma protein preservation were created as an improved method over previous plasma or serum collection tubes1, to stabilize the protein content of blood, enabling better protein biomarker discovery and proteomics experimentation from human blood. The BD P100 tubes contain 15.8 ml of spray-dried K2EDTA and a lyophilized proprietary broad spectrum cocktail of protease inhibitors to prevent coagulation and stabilize the plasma proteins. They also include a mechanical separator, which provides a physical barrier between plasma and cell pellets after centrifugation. Few methods have been devised to extract DNA from clotted blood samples collected in old plasma tubes2-4. Challenges from these methods were mainly associated with the type of separator inside the tubes (gel separator) and included difficulty in recovering the clotted blood, the inconvenience of fragmenting or dispersing the clot, and obstruction of the clot extraction by the separation gel.
We present the first method that extracts and purifies genomic DNA from blood drawn in the new BD P100 tubes. We compare the quality of the DNA sample from P100 tubes to that from EDTA tubes. Our approach is simple and efficient. It involves four major steps as follows: 1) the use of a plasma BD P100 (BD Diagnostics, Sparks, MD, USA) tube with mechanical separator for blood collection, 2) the removal of the mechanical separator using a combination of sucrose and a sterile paperclip metallic hook, 3) the separation of the buffy coat layer containing the white cells and 4) the isolation of the genomic DNA from the buffy coat using a regular commercial DNA extraction kit or a similar standard protocol.

We are describing a new method of isolating genomic DNA from whole blood collected for plasma/serology. After plasma collection, the compacted blood is usually discarded. Our novel method represents a significant improvement over existing methods and makes DNA and plasma available from a single collection, without requesting additional blood.

Laboratory tests can be done on the cellular or  fluid portions of the blood. The use of different blood collection tubes  determines the portion of the ...

Sampling Blood from the Lateral Tail Vein of the Rat

Blood samples are commonly obtained in many experimental contexts to measure targets of interest, including hormones, immune factors, growth factors, proteins, and glucose, yet the composition of the blood is dynamically regulated and easily perturbed. One factor that can change the blood composition is the stress response triggered by the sampling procedure, which can contribute to variability in the measures of interest. Here we describe a procedure for blood sampling from the lateral tail vein in the rat. This procedure offers significant advantages over other more commonly used techniques. It permits rapid sampling with minimal pain or invasiveness, without anesthesia or analgesia. Additionally, it can be used to obtain large volume samples (upwards of 1 ml in some rats), and it may be used repeatedly across experimental days. By minimizing the stress response and pain resulting from blood sampling, measures can more accurately reflect the true basal state of the animal, with minimal influence from the sampling procedure itself.

Blood samples are useful for assessing biomarkers of physiological states or disease in vivo. Here we describe the methodology to sample blood from the lateral tail vein in the rat. This method provides rapid samples with minimal pain and invasiveness.

Blood samples are commonly obtained in many experimental contexts to measure targets of interest, including hormones, immune factors, growth factors, ...

Amplification, Next-generation Sequencing, and Genomic DNA Mapping of Retroviral Integration Sites

Retroviruses exhibit signature integration preferences on both the local and global scales. Here, we present a detailed protocol for (1) generation of diverse libraries of retroviral integration sites using ligation-mediated PCR (LM-PCR) amplification and next-generation sequencing (NGS), (2) mapping the genomic location of each virus-host junction using BEDTools, and (3) analyzing the data for statistical relevance. Genomic DNA extracted from infected cells is fragmented by digestion with restriction enzymes or by sonication. After suitable DNA end-repair, double-stranded linkers are ligated onto the DNA ends, and semi-nested PCR is conducted using primers complementary to both the long terminal repeat (LTR) end of the virus and the ligated linker DNA. The PCR primers carry sequences required for DNA clustering during NGS, negating the requirement for separate adapter ligation. Quality control (QC) is conducted to assess DNA fragment size distribution and adapter DNA incorporation prior to NGS. Sequence output files are filtered for LTR-containing reads, and the sequences defining the LTR and the linker are cropped away. Trimmed host cell sequences are mapped to a reference genome using BLAT and are filtered for minimally 97% identity to a unique point in the reference genome. Unique integration sites are scrutinized for adjacent nucleotide (nt) sequence and distribution relative to various genomic features. Using this protocol, integration site libraries of high complexity can be constructed from genomic DNA in three days. The entire protocol that encompasses exogenous viral infection of susceptible tissue culture cells to integration site analysis can therefore be conducted in approximately one to two weeks. Recent applications of this technology pertain to longitudinal analysis of integration sites from HIV-infected patients.

We describe a protocol for amplifying retroviral integration sites from the genomic DNA of infected cells, sequencing the amplified virus-host junctions, and then mapping these sequences to a reference genome. We also describe techniques to quantify the distribution of integration sites relative to various genomic annotations using BEDTools.

Retroviruses exhibit signature integration preferences on both the local and global scales. Here, we present a detailed protocol for (1) generation of ...

Visualization of Surface-tethered Large DNA Molecules with a Fluorescent Protein DNA Binding Peptide

Large DNA molecules tethered on the functionalized glass surface have been utilized in polymer physics and biochemistry particularly for investigating interactions between DNA and its binding proteins. Here, we report a method that uses fluorescent microscopy for visualizing large DNA molecules tethered on the surface. First, glass coverslips are biotinylated and passivated by coating with biotinylated polyethylene glycol, which specifically binds biotinylated DNA via avidin protein linkers and significantly reduces undesirable binding from non-specific interactions of proteins or DNA molecules on the surface. Second, the DNA molecules are biotinylated by two different methods depending on their terminals. The blunt ended DNA is tagged with biotinylated dUTP at its 3' hydroxyl terminus, by terminal transferase, while the sticky ended DNA is hybridized with biotinylated complimentary oligonucleotides by DNA ligase. Finally, a microfluidic shear flow makes single DNA molecules stretch to their full contour lengths after being stained with fluorescent protein-DNA binding peptide (FP-DBP).

We present an approach for visualizing fluorescent protein DNA binding peptide (FP-DBP)-stained large DNA molecules tethered on the polyethylene glycol (PEG) and avidin-coated glass surface and stretched with microfluidic shear flows.

Large DNA molecules tethered on the functionalized glass surface have been utilized in polymer physics and biochemistry particularly for investigating ...

Advanced Confocal Microscopy Techniques to Study Protein-protein Interactions and Kinetics at DNA Lesions

Local microirradiation with lasers represents a useful tool for studies of DNA-repair-related processes in live cells. Here, we describe a methodological approach to analyzing protein kinetics at DNA lesions over time or protein-protein interactions on locally microirradiated chromatin. We also show how to recognize individual phases of the cell cycle using the Fucci cellular system to study cell-cycle-dependent protein kinetics at DNA lesions. A methodological description of the use of two UV lasers (355 nm and 405 nm) to induce different types of DNA damage is also presented. Only the cells microirradiated by the 405-nm diode laser proceeded through mitosis normally and were devoid of cyclobutane pyrimidine dimers (CPDs). We also show how microirradiated cells can be fixed at a given time point to perform immunodetection of the endogenous proteins of interest. For the DNA repair studies, we additionally describe the use of biophysical methods including FRAP (Fluorescence Recovery After Photobleaching) and FLIM (Fluorescence Lifetime Imaging Microscopy) in cells with spontaneously occurring DNA damage foci. We also show an application of FLIM-FRET (Fluorescence Resonance Energy Transfer) in experimental studies of protein-protein interactions.

Laser microirradiation is a useful tool for studies of DNA repair in living cells. A methodological approach for the use of UVA lasers to induce various DNA lesions is shown. We have optimized a method for local microirradiation that maintains the normal cell cycle; thus, irradiated cells proceed through mitosis.

Local microirradiation with lasers represents a useful tool for studies of DNA-repair-related processes in live cells. Here, we describe a methodological ...

Phloem Sap Sampling from Brassica napus for 3D-PAGE of Protein and Ribonucleoprotein Complexes

Sampling the phloem of higher plants is often laborious and significantly dependent on the plant species. However, proteome studies under denaturing conditions could be achieved in different plant species. Native protein:protein and protein:nucleic acid complexes from phloem samples have as yet scarcely been analyzed, although they might play important roles in maintenance of this specialized compartment or in long-distance signaling. Large molecular assemblies can be isolated using a blue native gel electrophoresis (BN-PAGE). Their protein components can be separated by a subsequent sodium dodecyl sulfate PAGE (SDS-PAGE). However, proteins with similar molecular weights co-migrate, what can hinder protein identification by mass spectrometry. Combining BN-PAGE with two different denaturing gel electrophoresis steps, namely Tris-Tricine-urea and SDS-PAGE, enables the additional separation of proteins according to their hydrophilicity/hydrophobicity and thus increases resolution and the success of protein identification. It even allows distinguishing proteins that only differ in their posttranslational modifications. In addition, blue native northern blotting can be applied to identify the RNA components in macromolecular complexes. We show that our protocol is suitable to unravel the protein and RNA components of native protein:protein and ribonucleoprotein (RNP) complexes occurring in phloem samples. Combining a blue native PAGE with two different denaturing PAGE steps can help to separate different kinds of large protein complexes, and also enables an increased identification rate of their components by mass spectrometry. Furthermore, the protocol is robust enough to simultaneously detect potentially bound nucleic acids within single protein complexes.

Phloem Sap Sampling from Brassica napus for 3D-PAGE of Protein and Ribonucleoprotein Complexes

Here we present a protocol to analyze the protein composition of large native protein:protein and protein:nucleic acid complexes from oilseed rape (B. napus) phloem exudate using a 3D polyacrylamide gel electrophoresis (PAGE) approach combining blue native (BN) with two denaturing PAGEs followed by mass spectrometric identification.

Sampling the phloem of higher plants is often laborious and significantly dependent on the plant species. However, proteome studies under denaturing ...