An Overview of Computational Tools of Nucleic Acid Binding Site Prediction for Site-specific Proteins and Nucleases.

Understanding the interaction mechanism of proteins and nucleic acids is one of the most fundamental problems for genome editing with engineered nucleases. Due to some limitations of experimental investigations, computational methods have played an important role in obtaining the knowledge of protein-nucleic acid interaction. Over the past few years, dozens of computational tools have been used for identification of nucleic acid binding site for site-specific proteins and design of site-specific nucleases because of their significant advantages in genome editing. Here, we review existing widely-used computational tools for target prediction of site-specific proteins as well as off-target prediction of site-specific nucleases. This article provides a list of on-line prediction tools according to their features followed by the description of computational methods used by these tools, which range from various sequence mapping algorithms (like Bowtie, FetchGWI and BLAST) to different machine learning methods (such as Support Vector Machine, hidden Markov models, Random Forest, elastic network and deep neural networks). We also make suggestions on the further development in improving the accuracy of prediction methods. This survey will provide a reference guide for computational biologists working in the field of genome editing.

Potential Role of the Last Half Repeat in TAL Effectors Revealed by a Molecular Simulation Study.

TAL effectors (TALEs) contain a modular DNA-binding domain that is composed of tandem repeats. In all naturally occurring TALEs, the end of tandem repeats is invariantly a truncated half repeat. To investigate the potential role of the last half repeat in TALEs, we performed comparative molecular dynamics simulations for the crystal structure of DNA-bound TALE AvrBs3 lacking the last half repeat and its modeled structure having the last half repeat. The structural stability analysis indicates that the modeled system is more stable than the nonmodeled system. Based on the principle component analysis, it is found that the AvrBs3 increases its structural compactness in the presence of the last half repeat. The comparison of DNA groove parameters of the two systems implies that the last half repeat also causes the change of DNA major groove binding efficiency. The following calculation of hydrogen bond reveals that, by stabilizing the phosphate binding with DNA at the C-terminus, the last half repeat helps to adopt a compact conformation at the protein-DNA interface. It further mediates more contacts between TAL repeats and DNA nucleotide bases. Finally, we suggest that the last half repeat is required for the high-efficient recognition of DNA by TALE.

Molecular Dynamics Simulations of Ternary Complexes: Comparisons of LEAFY Protein Binding to Different DNA Motifs.

LEAFY (LFY) is a plant-specific transcription factor, with a variety of roles in different species. LFY contains a conserved DNA-binding domain (DBD) that determines its DNA-binding specificity. Recently, the structures of the dimeric LFY-DBD bound to different DNA motifs were successively solved by X-ray crystallography. In this article, molecular dynamics (MD) simulations are employed to study two crystal structures of DNA-bound LFY protein from angiosperms and the moss Physcomitrella patens, respectively. The comparison of stabilities of the two systems is consistent with the experimental data of binding affinities. The calculation of hydrogen bonds showed that position 312 in LFY determines the difference of DNA-binding specificity. By using principal component analysis (PCA) and free energy landscape (FEL) methods, the open-close conformational change of the dimerization interface was found to be important for the system stability. At the dimerization interface, the protein-protein interaction has multiple influences on the cooperative DNA binding of LFY. The following analysis of DNA structural parameters further revealed that the protein-protein interaction contributes varying roles according to the specific DNA-binding efficiency. We propose that the protein-protein interaction serves a dual function as a connector between LFY monomers and a regulator of DNA-binding specificity. It will improve the robustness and adaptivity of the LFY-DNA ternary structure. This study provides some new insights into the understanding of the dynamics and interaction mechanism of dimeric LFY-DBD bound to DNA at the atomic level.

Molecular dynamics simulations of DNA-free and DNA-bound TAL effectors.

TAL (transcriptional activator-like) effectors (TALEs) are DNA-binding proteins, containing a modular central domain that recognizes specific DNA sequences. Recently, the crystallographic studies of TALEs revealed the structure of DNA-recognition domain. In this article, molecular dynamics (MD) simulations are employed to study two crystal structures of an 11.5-repeat TALE, in the presence and absence of DNA, respectively. The simulated results indicate that the specific binding of RVDs (repeat-variable diresidues) with DNA leads to the markedly reduced fluctuations of tandem repeats, especially at the two ends. In the DNA-bound TALE system, the base-specific interaction is formed mainly by the residue at position 13 within a TAL repeat. Tandem repeats with weak RVDs are unfavorable for the TALE-DNA binding. These observations are consistent with experimental studies. By using principal component analysis (PCA), the dominant motions are open-close movements between the two ends of the superhelical structure in both DNA-free and DNA-bound TALE systems. The open-close movements are found to be critical for the recognition and binding of TALE-DNA based on the analysis of free energy landscape (FEL). The conformational analysis of DNA indicates that the 5' end of DNA target sequence has more remarkable structural deformability than the other sites. Meanwhile, the conformational change of DNA is likely associated with the specific interaction of TALE-DNA. We further suggest that the arrangement of N-terminal repeats with strong RVDs may help in the design of efficient TALEs. This study provides some new insights into the understanding of the TALE-DNA recognition mechanism.

Molecular dynamics simulations of wild type and mutants of human complement receptor 2 complexed with C3d.

The interaction between human complement receptor type 2 (CR2) and antigen-bound C3d can bridge the innate and adaptive immune systems. The recently determined structure of the CR2(SCR1-2):C3d complex has revealed the expected binding interface of CR2-C3d. In this article, wild type (WT) and three mutants of the new structure are studied by molecular dynamics (MD) simulations. The differently decreased structural stabilities of the mutants relative to WT are shown to be consistent with the experimental data, which can be explained by the different hydrogen bond patterns at the interfaces. It is also found that two clusters of residues (D36/E37/E39 and E160/D163/E166) in the acidic pocket of C3d are important for CR2-C3d interactions, which is in good agreement with previous mutagenesis study. In addition, functional dynamics and the conformational change of CR2 are explored by using domain cross-correlation map (DCCM), principal component analysis (PCA), and free energy landscape (FEL) methods. The conformational change mainly corresponds to the opening of a V-shaped structure of CR2, which is consistent with the previously reported high interdomain flexibility of CR2. We further suppose that the opening of a V-shaped structure of CR2 may favor the binding stability of CR2(SCR1-2):C3d. This study would provide some new insights into the understanding of the CR2-C3d interaction mechanism.

Network models reveal stability and structural rearrangement of signal recognition particle.

The signal recognition particle (SRP) and its receptors (SR) mediate the cotranslational targeting of the membrane and secretory proteins in all cells. In Escherichia coli, SRP is composed of the Ffh protein and the 4.5S SRP RNA. Ffh is a multidomain protein comprising a methionine-rich (M) domain, a helical N domain, and a Ras-like guanine triphosphatase (GTPase) (G) domain. The N and G domains are commonly referred to as one structural unit, the NG domain. In this article, the complex structure of SRP and SR is investigated with the Gaussian network model (GNM) and anisotropic network model (ANM). GNM provides the information of structure stability. It is found that the intermolecular interactions between SRP and SR can obviously decrease the fluctuation of NG domains. Nevertheless, the large structural rearrangement will take place during the cotranslational protein targeting cycle. Hence, the moving directions of fluctuation regions are further ascertained by using cross-correlation analysis and the ANM. The NG domain of Ffh undergoes a clockwise rotation around the GM linker and the M domain of Ffh shows an opposite direction to the NG domain. These functional movements will facilitate the SRP structure to transform into the free form and the sequence-bound form. These simple coarse-grained analyses can be used as a general and quick method for the mechanism studies of protein assembly and supramolecular systems.

Computational model for protein unfolding simulation.

The protein folding problem is one of the fundamental and important questions in molecular biology. However, the all-atom molecular dynamics studies of protein folding and unfolding are still computationally expensive and severely limited by the time scale of simulation. In this paper, a simple and fast protein unfolding method is proposed based on the conformational stability analyses and structure modeling. In this method, two structure-based conditions are considered to identify the unstable regions of proteins during the unfolding processes. The protein unfolding trajectories are mimicked through iterative structure modeling according to conformational stability analyses. Two proteins, chymotrypsin inhibitor 2 (CI2) and α -spectrin SH3 domain (SH3) were simulated by this method. Their unfolding pathways are consistent with the previous molecular dynamics simulations. Furthermore, the transition states of the two proteins were identified in unfolding processes and the theoretical Φ values of these transition states showed significant correlations with the experimental data (the correlation coefficients are >0.8). The results indicate that this method is effective in studying protein unfolding. Moreover, we analyzed and discussed the influence of parameters on the unfolding simulation. This simple coarse-grained model may provide a general and fast approach for the mechanism studies of protein folding.

Stability and folding behavior analysis of zinc-finger using simple models.

Zinc-fingers play crucial roles in regulating gene expression and mediating protein-protein interactions. In this article, two different proteins (Sp1f2 and FSD-1) are investigated using the Gaussian network model and anisotropy elastic network model. By using these simple coarse-grained methods, we analyze the structural stabilization and establish the unfolding pathway of the two different proteins, in good agreement with related experimental and molecular dynamics simulation data. From the analysis, it is also found that the folding process of the zinc-finger motif is predominated by several factors. Both the zinc ion and C-terminal loop affect the folding pathway of the zinc-finger motif. Knowledge about the stability and folding behavior of zinc-fingers may help in understanding the folding mechanisms of the zinc-finger motif and in designing new zinc-fingers. Meanwhile, these simple coarse-grained analyses can be used as a general and quick method for mechanistic studies of metalloproteins.

Substrate recognition and transport behavior analyses of amino acid antiporter with coarse-grained models.

The L-arginine (Arg)/agmatine (Agm) antiporter AdiC is a vital transport protein of the arginine-dependent extreme acid resistance system of enteric bacteria. Recently, both substrate-free and Arg-bound structures of AdiC were determined by X-ray crystallography. In this article, the two different proteins were investigated with three simple models. Gaussian network model provided the information of conformational changes. It is found that Arg binding induces structural rearrangement in the extracellular domain, and transmembrane helix 6 (TM6) has the most pronounced trend of conformational changes. The moving directions of fluctuation regions were further ascertained by using anisotropy elastic network model and cross-correlation analysis. Interestingly, the two substrate-binding sites hypothesis of AdiC was confirmed directly by molecular docking. Furthermore, the binding preferences of these two sites were explained from the aspects of electrostatic complementarity and geometric matching. These simple coarse-grained analyses can be used as a general and quick method for the mechanism studies of transport proteins.