Team:TEC-Costa Rica/Project/Modelling Simulations/Structural Modelling

Modelling & Simulations
Structural Modelling

On the structural modelling area, we worked on the visualization of different Streptococcus pyogenes Cas9 Structures on varying conformations (Sternberg, LaFrance, Kaplan, & Doudna, 2015). We selected the PDB files 4cmp (Cas9), 4zt0 (Cas9+sgRNA) and 4un3 (Cas9+sgRNA+tDNA) as described on the article and visualized the structures using the molecular modelling package Chimera, by the University of California in San Francisco Resource for Biocomputing, Visualization, and Informatics.


View of the PDB file 4un3 as seen when fetched from the database

View of the PDB file 4un3 as seen when fetched from the database

The first thing we did, was to color the domains of the Cas9 protein, using as reference the information in the Uniprot accession Q99ZW2. The domains were colored in all the three conformations of the Cas9.

The domains were colored as:

  • Yellow: HNH.

  • Blue: RuvC.

  • Green: PAM Interacting.

  • Gray: REC lobe.

  • Red: sgRNA.

  • Beige: DNA.

In order to be able to align the models, we had to have the same number of chains in each model, so we removed the sgRNA and the DNA from 4zt0 and 4un3, and were left with these models:

First conformation of Cas9 protein, no sgRNA or DNA attached.

First conformation of Cas9 protein, no sgRNA or DNA attached.

Second conformation of Cas9, after the binding of the sgRNA.

Second conformation of Cas9, after the binding of the sgRNA.

Third conformation of Cas9 protein, after the binding of the sgRNA and the target DNA.

Third conformation of Cas9 protein, after the binding of the sgRNA and the target DNA.


After this was done, the three models were aligned in their PAM Interacting and RuvC domains using the “MatchMaker” structure comparison tool. These domains remain relatively stable through all the conformational change, being the HNH and the REC lobes the one which present a bigger change. In order to see more clearly the behaviour of the protein during its conformational change, the tool “Morph Conformations” was used to create an animated morph of the changing protein. This animation shows that the REC lobe suffers an aggressive reconformation when it interacts with the sgRNA, and is brought together with the HNH domain. This suggests that the HNH and REC zones of the Cas9 protein are good candidates to engineer the protein in order to be able to detect the recognition of the target DNA through conformational change.


The sites for the engineering of the protein were selected based on the article by Oakes, et. al., (2016). The article showed the sites throughout the protein that can tolerate insertions without affecting the functionality of the protein. We made a screening of the hotspots in the REC lobe and the HNH domain and highlighted them in Chimera. Then, using the conformation we made, and the Distance tool, we selected the pairs of sites that were far in the first conformation and came together as the sgRNA and the target binded to the dCas9.


We found two pairs of sites that conformed to the specifications: L390 & E802, and N588 & N888. The first pair of sites has a big distance in the first conformation: 55Å (Figure 8). After that sgRNA and the target bind to the Cas9, this distance reduces to 17Å (Figure 9).


Figure 8. L390 and E802 on the first conformation.

Figure 8. L390 and E802 on the first conformation.

Figure 9. L390 and E802 after the conformational change.

Figure 9. L390 and E802 after the conformational change.


The second pair, N588 and N888 also comes together as an effect of the conformational change, but the pair comes close with only adding the sgRNA, so the change is not specific to the target. When the Cas is not bound, the distance is of 95Å (Figure 10), then drops to 16Å (Figure 11) after the binding of the sgRNA and lastly increases to 17Å after the binding of the target (Figure 12).


Figure 10. First conformation of N588 and N888.

Figure 10. First conformation of N588 and N888.

Figure 11. Second conformation.

Figure 11. Second conformation.

Figure 12. Last conformation of N588 and N888.

Figure 12. Last conformation of N588 and N888.



These two pairs of sites are the ones we used for the design of our modular dCas9. We are using both pairs as a way to prove that conformational change can be used as a method of detection

Mathematical Modelling

Mathematical Modelling

Our mathematical model was based on Boolean logic gates, because we want a signal response only if all our inputs are present, and an analogical output is not our main focus...


Read More