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

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 ...

Beijing Computational Science Research Center

Shenzhen JL computational science and applied research institute

Rotary sequential hydrolysis of the metabolic machine F-ATPase is a prominent manifestation of high coordination among multiple chemical sites in ring-shaped molecular machines, and it is also functionally essential for F to tightly couple chemical reactions and central γ-shaft rotation. High-speed AFM experiments have identified that sequential hydrolysis is maintained in the F stator ring even in the absence of the γ-rotor. To explore the origins of intrinsic sequential performance, we computationally investigated essential inter-subunit couplings on the hexameric ring of mitochondrial and bacterial F. We first reproduced in stochastic Monte Carlo simulations the experimentally determined sequential hydrolysis schemes by kinetically imposing inter-subunit couplings and following subsequent tri-site ATP hydrolysis cycles on the F ring. We found that the key couplings to support the sequential hydrolysis are those that accelerate neighbor-site ADP and Pi release upon a certain ATP binding or hydrolysis reaction. The kinetically identified couplings were then examined in atomistic molecular dynamics simulations at a coarse-grained level to reveal the underlying structural mechanisms. To do that, we enforced targeted conformational changes of ATP binding or hydrolysis to one chemical site on the F ring and monitored the ensuing conformational responses of the neighboring sites using structure-based simulations. Notably, we found asymmetrical neighbor-site opening that facilitates ADP release upon enforced ATP binding. We also captured a complete charge-hopping process of the Pi release subsequent to enforced ATP hydrolysis in the neighbor site, confirming recent single-molecule analyses with regard to the role of ATP hydrolysis in F. Our studies therefore elucidate both the coordinated chemical kinetics and structural dynamics mechanisms underpinning the sequential operation of the F ring.

Deciphering Intrinsic Inter-subunit Couplings that Lead to Sequential Hydrolysis of F-ATPase Ring.

Green fluorescent protein (GFP) is a widely used biomarker that demands systematical rational approaches to its structure function redesign. In this work, we mainly utilized atomistic molecular dynamics simulations to inspect and visualize internal fluctuation and coordination around chromophore inside GFP, from water to nonpolar octane solvent. We found that GFP not only maintains its β-barrel structure well into the octane, but also sustains internal residue and water coordination to position the chromophore stably while suppress dihedral fluctuations of the chromophore, so that functional robustness of GFP is achieved. Our accompanied fluorescence microscope measurements accordingly confirmed the GFP functioning into the octane. Furthermore, we identified that crucial water sites inside GFP along with permeable pores on the β-barrel of the protein are largely preserved from the water to the octane solvent, which allows sufficiently fast exchanges of internal water with the bulk or with the water layer kept on the surface of the protein. By additionally pulling GFP from bulk water to octane, we suggest that the GFP function can be well maintained into the nonpolar solvent as long as, first, the protein does not denature in the nonpolar solvent nor across the polar-nonpolar solvent interface; second, a minimal set of water molecules are in accompany with the protein; third, the nonpolar solvent molecules may need to be large enough to be nonpermeable via the water pores on the β-barrel.

Inspecting fluctuation and coordination around chromophore inside green fluorescent protein from water to nonpolar solvent.

Oncogenic protein Myc serves as a transcription factor to control cell metabolisms. Myc dimerizes via leucine zipper with its associated partner protein Max to form a heterodimer structure, which then binds target DNA sequences to regulate gene transcription. The regulation depends on Myc-Max binding to DNA and searching for target sequences via diffusional motions along DNA. Here, we conduct structure-based molecular dynamics (MD) simulations to investigate the diffusion dynamics of the Myc-Max heterodimer along DNA. We found that the heterodimer protein slides on the DNA in a rotation-uncoupled manner in coarse-grained simulations, as its two helical DNA binding basic regions (BRs) alternate between open and closed conformations via inchworm stepping motions. In such motions, the two BRs of the heterodimer step across the DNA strand one by one, with step sizes reaching about half of a DNA helical pitch length. Atomic MD simulations of the Myc-Max heterodimer in complex with DNA have also been conducted. Hydrogen bond interactions are revealed between the two BRs and two complementary DNA strands, respectively. In the non-specific DNA binding, the BR from Myc shows an onset of stepping on one association DNA strand and starts detaching from the other strand. Overall, our simulation studies suggest that the inchworm stepping motions of the Myc-Max heterodimer can be achieved during the protein diffusion along DNA.

Inchworm stepping of Myc-Max heterodimer protein diffusion along DNA.

Transcription factor (TF) target search on genome is highly essential for gene expression and regulation. High-resolution determination of TF diffusion along DNA remains technically challenging. Here, we constructed a TF model system using the plant WRKY domain protein in complex with DNA from crystallography and demonstrated microsecond diffusion dynamics of WRKY on DNA by employing all-atom molecular-dynamics (MD) simulations. Notably, we found that WRKY preferentially binds to one strand of DNA with significant energetic bias compared with the other, or nonpreferred strand. The preferential DNA-strand binding becomes most prominent in the static process, from nonspecific to specific DNA binding, but less distinct during diffusive movements of the domain protein on the DNA. Remarkably, without employing acceleration forces or bias, we captured a complete one-base-pair stepping cycle of the protein tracking along major groove of DNA with a homogeneous poly-adenosine sequence, as individual hydrogen bonds break and reform at the protein-DNA binding interface. Further DNA-groove tracking motions of the protein forward or backward, with occasional sliding as well as strand crossing to minor groove of DNA, were also captured. The processive diffusion of WRKY along DNA has been further sampled via coarse-grained MD simulations. The study thus provides structural dynamics details on diffusion of a small TF domain protein, suggests how the protein approaches a specific recognition site on DNA, and supports further high-precision experimental detection. The stochastic movements revealed in the TF diffusion also provide general clues about how other protein walkers step and slide along DNA.

Revealing atomic-scale molecular diffusion of a plant-transcription factor WRKY domain protein along DNA.

Cas1 and Cas2 are highly conserved proteins across clustered-regularly-interspaced-short-palindromic-repeat-Cas systems and play a significant role in protospacer acquisition. Based on crystal structure of twofold symmetric Cas1-Cas2 in complex with dual-forked protospacer DNA (psDNA), we conducted all-atom molecular dynamics simulations to study the psDNA binding, recognition, and response to cleavage on the protospacer-adjacent-motif complementary sequence, or PAMc, of Cas1-Cas2. In the simulation, we noticed that two active sites of Cas1 and Cas1' bind asymmetrically to two identical PAMc on the psDNA captured from the crystal structure. For the modified psDNA containing only one PAMc, as that to be recognized by Cas1-Cas2 in general, our simulations show that the non-PAMc association site of Cas1-Cas2 remains destabilized until after the stably bound PAMc being cleaved at the corresponding association site. Thus, long-range correlation appears to exist upon the PAMc cleavage between the two active sites (∼10 nm apart) on Cas1-Cas2, which can be allosterically mediated by psDNA and Cas2 and Cas2' in bridging. To substantiate such findings, we conducted repeated runs and further simulated Cas1-Cas2 in complex with synthesized psDNA sequences psL and psH, which have been measured with low and high frequency in acquisition, respectively. Notably, such intersite correlation becomes even more pronounced for the Cas1-Cas2 in complex with psH but remains low for the Cas1-Cas2 in complex with psL. Hence, our studies demonstrate that PAMc recognition and cleavage at one active site of Cas1-Cas2 may allosterically regulate non-PAMc association or even cleavage at the other site, and such regulation can be mediated by noncatalytic Cas2 and DNA protospacer to possibly support the ensued psDNA acquisition.

Allosteric regulation in CRISPR/Cas1-Cas2 protospacer acquisition mediated by DNA and Cas2.

In this work, we computationally investigated how a viral RNA polymerase (RNAP) from bacteriophage T7 evolves into RNAP variants under lab-directed evolution to switch recognition from T7 promoter to T3 promoter in transcription initiation. We first constructed a closed initiation complex for the wild-type T7 RNAP and then for six mutant RNAPs discovered from phage-assisted continuous evolution experiments. All-atom molecular dynamics simulations up to 1 μs each were conducted on these RNAPs in a complex with the T7 and T3 promoters. Our simulations show notably that protein-DNA electrostatic interactions or stabilities at the RNAP-DNA promoter interface well dictate the promoter recognition preference of the RNAP and variants. Key residues and structural elements that contribute significantly to switching the promoter recognition were identified. Followed by a first point mutation N748D on the specificity loop to slightly disengage the RNAP from the promoter to hinder the original recognition, we found an auxiliary helix (206-225) that takes over switching the promoter recognition upon further mutations (E222K and E207K) by forming additional charge interactions with the promoter DNA and reorientating differently on the T7 and T3 promoters. Further mutations on the AT-rich loop and the specificity loop can fully switch the RNAP-promoter recognition to the T3 promoter. Overall, our studies reveal energetics and structural dynamics details along an exemplary directed evolutionary path of the phage RNAP variants for a rewired promoter recognition function. The findings demonstrate underlying physical mechanisms and are expected to assist knowledge and data learning or rational redesign of the protein enzyme structure function.

Liqiang Dai

Beijing Computational Science Research Center

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