My research is mainly on physical modeling and simulation of biomolecular machinery, particularly those that regulate genetic process. I have recently studied a viral DNA packaging motor, DNA/RNA helicases, polymerases, etc.. I'm interested in understanding how these microscopic machines react to thermal fluctuations and environmental noises, how they achieve regular movements and high-fidelity functions, and how they emerge, evolve through and impact on evolution. I'm also interested in bringing biomolecular redesign from in silico to in vitro and vivo. Supervisors: George Oster and Carlos Bustamante
ABSTRACTCas1 and Cas2 are highly conserved proteins across CRISPR-Cas systems and play a signific... more ABSTRACTCas1 and Cas2 are highly conserved proteins across CRISPR-Cas systems and play a significant role in protospacer acquisition. Here we study the protospacer (or ps) DNA binding, recognition, and response to cleavage on the protospacer-adjacent-motif complementary sequence or PAMc by Cas1-Cas2, implementing all-atom molecular dynamics simulations. First, we noticed that two active sites of Cas1&1’ bind asymmetrically to two identical PAMc in the simulation. For psDNA containing only one PAMc to be recognized, it is then found that the non-PAMc association site remains destabilized until after the bound PAMc being cleaved. Thus, correlation appears to exist between the two active sites, which can be allosterically mediated by psDNA and Cas2&2’ in bridging. To substantiate such findings, we further simulated Cas1-Cas2 in complex with synthesized psDNA sequences psL and psH, which have been measured with low and high efficiency in acquisition, respectively. Notably, such inter-si...
Oncogenic protein Myc serves as a transcription factor to control cell metabolisms. Myc dimerizes... more 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 by 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 up about half of a DNA helical pitch length. Atomic MD simulations of the Myc-Max heterodimer in complex with DNA have ...
Oncogenic protein Myc serves as a transcription factor to control cell metabolisms. Myc dimerizes... more 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 by 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 up about half of a DNA helical pitch length. Atomic MD simulations of the Myc-Max heterodimer in complex with DNA have ...
Transcription factor (TF) target search on genome is highly essential for gene expression and reg... more 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 of the plant WRKY domain protein in complex with DNA from crystallography and demonstrated microsecond diffusion dynamics of WRKY on the DNA employing all-atom molecular dynamics (MD) simulations. Notably, we found that WRKY preferentially binds to the Crick strand of DNA with significantly stronger energetic association than to the Watson strand. The preferential binding becomes highly prominent from non-specific to specific DNA binding, but less distinct from static binding to diffusive movements of WRKY on the DNA. Remarkably, without employing acceleration forces or bias, we captured a complete one-base pair (bp) stepping cycle of WRKY tracking along major groove of DNA with homogenous (AT)n sequence, as individual protein-DNA contacts break and refo...
Polymerases select nucleotides according to a template before incorporating them for chemical syn... more Polymerases select nucleotides according to a template before incorporating them for chemical synthesis during gene replication or transcription. Efficient selection to achieve sufficiently high fidelity and speed is essential for polymerase function. Due to multiple kinetic steps detected in a polymerase elongation cycle, there exist multiple selection checkpoints to allow different strategies of fidelity control. In our current work, we examined step-by-step selections in an elongation cycle that have conformational transition rates tuned one at a time, with a controlled differentiation free energy between the right and wrong nucleotides at each checkpoint. The elongation is sustained at non-equilibrium steady state with constant free energy input and heat dissipation. It is found that a selection checkpoint in the later stage of a reaction path has less capability for error reduction. Hence, early selection is essential to achieve an efficient fidelity control. In particular, for...
Rotary sequential hydrolysis of the metabolic machine F1-ATPase is a prominent manifestation of h... more Rotary sequential hydrolysis of the metabolic machine F1-ATPase is a prominent manifestation of high coordination among multiple chemical sites in ring-shaped molecular machines, and it is also functionally essential for F1 to tightly couple chemical reactions and central γ-shaft rotation. High-speed AFM experiments have identified that sequential hydrolysis is maintained in the F1 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 F1. 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 F1 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 ...
Here, we studied the complete process of a viral T7 RNA polymerase (RNAP) translocation on DNA du... more Here, we studied the complete process of a viral T7 RNA polymerase (RNAP) translocation on DNA during transcription elongation by implementing extensive all-atom molecular dynamics (MD) simulations to construct a Markov state model (MSM). Our studies show that translocation proceeds in a Brownian motion, and the RNAP thermally transits among multiple metastable states. We observed non-synchronized backbone movements of the nucleic acid (NA) chains with the RNA translocation accomplished first, while the template DNA lagged. Notably, both the O-helix and Y-helix on the fingers domain play key roles in facilitating NA translocation through the helix opening. The helix opening allows a key residue Tyr639 to become inserted into the active site, which pushes the RNA-DNA hybrid forward. Another key residue, Phe644, coordinates the downstream template DNA motions by stacking and un-stacking with a transition nucleotide (TN) and its adjacent nucleotide. Moreover, the O-helix opening at pre...
Pyrophosphate ion (PPi) release during transcription elongation is a signature step in each nucle... more Pyrophosphate ion (PPi) release during transcription elongation is a signature step in each nucleotide addition cycle. The kinetics and energetics of the process as well as how it proceeds with substantial conformational changes of the polymerase complex determine the mechano-chemical coupling mechanism of the transcription elongation. Here we investigated detailed dynamics of the PPi release process in a single-subunit RNA polymerase (RNAP) from bacteriophage T7, implementing all-atom molecular dynamics (MD) simulations. We obtained a jump-from-cavity kinetic model of the PPi release utilizing extensive nanosecond MD simulations. We found that the PPi release in T7 RNAP is initiated by the PPi dissociation from two catalytic aspartic acids, followed by a comparatively slow jump-from-cavity activation process. Combining with a number of microsecond long MD simulations, we also found that the activation process is hindered by charged residue associations as well as by local steric an...
Nucleotide selection is essential for fidelity control in gene replication and transcription. Rec... more Nucleotide selection is essential for fidelity control in gene replication and transcription. Recent work on T7 RNA polymerase suggested that a small posttranslocation free energy bias stabilizes Tyr(639) in the active site to aid nucleotide selection. However, it was not clear exactly how Tyr(639) assists the selection. Here we report a molecular-dynamics simulation study revealing atomistic detail of this critical selectivity. The study shows first that Tyr(639) blocks the active site at posttranslocation by marginally stacking to the end basepair of the DNA-RNA hybrid. The study then demonstrates that at the nucleotide preinsertion state, a cognate RNA nucleotide does not affect the local Tyr(639) stabilization, whereas a noncognate nucleotide substantially stabilizes Tyr(639) so that Tyr(639) keeps blocking the active site. As a result, further nucleotide insertion into the active site, which requires moving Tyr(639) out of the site, would be hindered for the noncognate nucleoti...
The RNA polymerase II elongation is central in eukaryotic transcription. Although multiple interm... more The RNA polymerase II elongation is central in eukaryotic transcription. Although multiple intermediates of the elongation complex have been identified, the dynamical mechanisms remain elusive or controversial. Here we build a structure-based kinetic model of a full elongation cycle of polymerase II, taking into account transition rates and conformational changes characterized from both single molecule experimental studies and computational simulations at atomistic scale. Our model suggests a force-dependent slow transition detected in the single molecule experiments corresponds to an essential conformational change of a trigger loop (TL) opening prior to the polymerase translocation. The analyses on mutant study of E1103G and on potential sequence effects of the translocation substantiate this proposal. Our model also investigates another slow transition detected in the transcription elongation cycle which is independent of mechanical force. If this force-independent slow transitio...
The pentameric ATPase motor gp16 packages double-stranded DNA into the bacteriophage phi29 virus ... more The pentameric ATPase motor gp16 packages double-stranded DNA into the bacteriophage phi29 virus capsid. On the basis of the results of single-molecule experimental studies, we propose a push and roll mechanism to explain how the packaging motor translocates the DNA in bursts of four 2.5 bp power strokes, while rotating the DNA. In this mechanism, each power stroke accompanies P(i) release after ATP hydrolysis. Since the high-resolution structure of the gp16 motor is not available, we borrowed characterized features from the P4 RNA packaging motor in bacteriophage phi12. For each power stroke, a lumenal lever from a single subunit is electrostatically steered to the DNA backbone. The lever then pushes sterically, orthogonal to the backbone axis, such that the right-handed DNA helix is translocated and rotated in a left-handed direction. The electrostatic association allows tight coupling between the lever and the DNA and prevents DNA from slipping back. The lever affinity for DNA de...
The pentameric ATPase motor gp16 packages double-stranded DNA into the bacteriophage phi29 virus ... more The pentameric ATPase motor gp16 packages double-stranded DNA into the bacteriophage phi29 virus capsid. On the basis of the results of single-molecule experimental studies, we propose a push and roll mechanism to explain how the packaging motor translocates the DNA in bursts of four 2.5 bp power strokes, while rotating the DNA. In this mechanism, each power stroke accompanies P(i) release after ATP hydrolysis. Since the high-resolution structure of the gp16 motor is not available, we borrowed characterized features from the P4 RNA packaging motor in bacteriophage phi12. For each power stroke, a lumenal lever from a single subunit is electrostatically steered to the DNA backbone. The lever then pushes sterically, orthogonal to the backbone axis, such that the right-handed DNA helix is translocated and rotated in a left-handed direction. The electrostatic association allows tight coupling between the lever and the DNA and prevents DNA from slipping back. The lever affinity for DNA de...
ABSTRACTCas1 and Cas2 are highly conserved proteins across CRISPR-Cas systems and play a signific... more ABSTRACTCas1 and Cas2 are highly conserved proteins across CRISPR-Cas systems and play a significant role in protospacer acquisition. Here we study the protospacer (or ps) DNA binding, recognition, and response to cleavage on the protospacer-adjacent-motif complementary sequence or PAMc by Cas1-Cas2, implementing all-atom molecular dynamics simulations. First, we noticed that two active sites of Cas1&1’ bind asymmetrically to two identical PAMc in the simulation. For psDNA containing only one PAMc to be recognized, it is then found that the non-PAMc association site remains destabilized until after the bound PAMc being cleaved. Thus, correlation appears to exist between the two active sites, which can be allosterically mediated by psDNA and Cas2&2’ in bridging. To substantiate such findings, we further simulated Cas1-Cas2 in complex with synthesized psDNA sequences psL and psH, which have been measured with low and high efficiency in acquisition, respectively. Notably, such inter-si...
Oncogenic protein Myc serves as a transcription factor to control cell metabolisms. Myc dimerizes... more 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 by 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 up about half of a DNA helical pitch length. Atomic MD simulations of the Myc-Max heterodimer in complex with DNA have ...
Oncogenic protein Myc serves as a transcription factor to control cell metabolisms. Myc dimerizes... more 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 by 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 up about half of a DNA helical pitch length. Atomic MD simulations of the Myc-Max heterodimer in complex with DNA have ...
Transcription factor (TF) target search on genome is highly essential for gene expression and reg... more 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 of the plant WRKY domain protein in complex with DNA from crystallography and demonstrated microsecond diffusion dynamics of WRKY on the DNA employing all-atom molecular dynamics (MD) simulations. Notably, we found that WRKY preferentially binds to the Crick strand of DNA with significantly stronger energetic association than to the Watson strand. The preferential binding becomes highly prominent from non-specific to specific DNA binding, but less distinct from static binding to diffusive movements of WRKY on the DNA. Remarkably, without employing acceleration forces or bias, we captured a complete one-base pair (bp) stepping cycle of WRKY tracking along major groove of DNA with homogenous (AT)n sequence, as individual protein-DNA contacts break and refo...
Polymerases select nucleotides according to a template before incorporating them for chemical syn... more Polymerases select nucleotides according to a template before incorporating them for chemical synthesis during gene replication or transcription. Efficient selection to achieve sufficiently high fidelity and speed is essential for polymerase function. Due to multiple kinetic steps detected in a polymerase elongation cycle, there exist multiple selection checkpoints to allow different strategies of fidelity control. In our current work, we examined step-by-step selections in an elongation cycle that have conformational transition rates tuned one at a time, with a controlled differentiation free energy between the right and wrong nucleotides at each checkpoint. The elongation is sustained at non-equilibrium steady state with constant free energy input and heat dissipation. It is found that a selection checkpoint in the later stage of a reaction path has less capability for error reduction. Hence, early selection is essential to achieve an efficient fidelity control. In particular, for...
Rotary sequential hydrolysis of the metabolic machine F1-ATPase is a prominent manifestation of h... more Rotary sequential hydrolysis of the metabolic machine F1-ATPase is a prominent manifestation of high coordination among multiple chemical sites in ring-shaped molecular machines, and it is also functionally essential for F1 to tightly couple chemical reactions and central γ-shaft rotation. High-speed AFM experiments have identified that sequential hydrolysis is maintained in the F1 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 F1. 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 F1 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 ...
Here, we studied the complete process of a viral T7 RNA polymerase (RNAP) translocation on DNA du... more Here, we studied the complete process of a viral T7 RNA polymerase (RNAP) translocation on DNA during transcription elongation by implementing extensive all-atom molecular dynamics (MD) simulations to construct a Markov state model (MSM). Our studies show that translocation proceeds in a Brownian motion, and the RNAP thermally transits among multiple metastable states. We observed non-synchronized backbone movements of the nucleic acid (NA) chains with the RNA translocation accomplished first, while the template DNA lagged. Notably, both the O-helix and Y-helix on the fingers domain play key roles in facilitating NA translocation through the helix opening. The helix opening allows a key residue Tyr639 to become inserted into the active site, which pushes the RNA-DNA hybrid forward. Another key residue, Phe644, coordinates the downstream template DNA motions by stacking and un-stacking with a transition nucleotide (TN) and its adjacent nucleotide. Moreover, the O-helix opening at pre...
Pyrophosphate ion (PPi) release during transcription elongation is a signature step in each nucle... more Pyrophosphate ion (PPi) release during transcription elongation is a signature step in each nucleotide addition cycle. The kinetics and energetics of the process as well as how it proceeds with substantial conformational changes of the polymerase complex determine the mechano-chemical coupling mechanism of the transcription elongation. Here we investigated detailed dynamics of the PPi release process in a single-subunit RNA polymerase (RNAP) from bacteriophage T7, implementing all-atom molecular dynamics (MD) simulations. We obtained a jump-from-cavity kinetic model of the PPi release utilizing extensive nanosecond MD simulations. We found that the PPi release in T7 RNAP is initiated by the PPi dissociation from two catalytic aspartic acids, followed by a comparatively slow jump-from-cavity activation process. Combining with a number of microsecond long MD simulations, we also found that the activation process is hindered by charged residue associations as well as by local steric an...
Nucleotide selection is essential for fidelity control in gene replication and transcription. Rec... more Nucleotide selection is essential for fidelity control in gene replication and transcription. Recent work on T7 RNA polymerase suggested that a small posttranslocation free energy bias stabilizes Tyr(639) in the active site to aid nucleotide selection. However, it was not clear exactly how Tyr(639) assists the selection. Here we report a molecular-dynamics simulation study revealing atomistic detail of this critical selectivity. The study shows first that Tyr(639) blocks the active site at posttranslocation by marginally stacking to the end basepair of the DNA-RNA hybrid. The study then demonstrates that at the nucleotide preinsertion state, a cognate RNA nucleotide does not affect the local Tyr(639) stabilization, whereas a noncognate nucleotide substantially stabilizes Tyr(639) so that Tyr(639) keeps blocking the active site. As a result, further nucleotide insertion into the active site, which requires moving Tyr(639) out of the site, would be hindered for the noncognate nucleoti...
The RNA polymerase II elongation is central in eukaryotic transcription. Although multiple interm... more The RNA polymerase II elongation is central in eukaryotic transcription. Although multiple intermediates of the elongation complex have been identified, the dynamical mechanisms remain elusive or controversial. Here we build a structure-based kinetic model of a full elongation cycle of polymerase II, taking into account transition rates and conformational changes characterized from both single molecule experimental studies and computational simulations at atomistic scale. Our model suggests a force-dependent slow transition detected in the single molecule experiments corresponds to an essential conformational change of a trigger loop (TL) opening prior to the polymerase translocation. The analyses on mutant study of E1103G and on potential sequence effects of the translocation substantiate this proposal. Our model also investigates another slow transition detected in the transcription elongation cycle which is independent of mechanical force. If this force-independent slow transitio...
The pentameric ATPase motor gp16 packages double-stranded DNA into the bacteriophage phi29 virus ... more The pentameric ATPase motor gp16 packages double-stranded DNA into the bacteriophage phi29 virus capsid. On the basis of the results of single-molecule experimental studies, we propose a push and roll mechanism to explain how the packaging motor translocates the DNA in bursts of four 2.5 bp power strokes, while rotating the DNA. In this mechanism, each power stroke accompanies P(i) release after ATP hydrolysis. Since the high-resolution structure of the gp16 motor is not available, we borrowed characterized features from the P4 RNA packaging motor in bacteriophage phi12. For each power stroke, a lumenal lever from a single subunit is electrostatically steered to the DNA backbone. The lever then pushes sterically, orthogonal to the backbone axis, such that the right-handed DNA helix is translocated and rotated in a left-handed direction. The electrostatic association allows tight coupling between the lever and the DNA and prevents DNA from slipping back. The lever affinity for DNA de...
The pentameric ATPase motor gp16 packages double-stranded DNA into the bacteriophage phi29 virus ... more The pentameric ATPase motor gp16 packages double-stranded DNA into the bacteriophage phi29 virus capsid. On the basis of the results of single-molecule experimental studies, we propose a push and roll mechanism to explain how the packaging motor translocates the DNA in bursts of four 2.5 bp power strokes, while rotating the DNA. In this mechanism, each power stroke accompanies P(i) release after ATP hydrolysis. Since the high-resolution structure of the gp16 motor is not available, we borrowed characterized features from the P4 RNA packaging motor in bacteriophage phi12. For each power stroke, a lumenal lever from a single subunit is electrostatically steered to the DNA backbone. The lever then pushes sterically, orthogonal to the backbone axis, such that the right-handed DNA helix is translocated and rotated in a left-handed direction. The electrostatic association allows tight coupling between the lever and the DNA and prevents DNA from slipping back. The lever affinity for DNA de...
Uploads
Papers by Jin Yu