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

Comparative Structural and Dynamics Study of Free and gRNA-bound FnCas9 and SpCas9 Proteins

•

gRNA binding caused Cas9 to open up, exposing hydrophobic residues.

•

Concerted domain movement was observed in both SpCas9 and FnCas9.

•

gRNA binding led to decrease in helicity and increase in structural transitions.

•

High binding affinity of FnCas9 with gRNA is mostly due to electrostatic interaction.

•

Arginine-helix has shown a pivotal role in the binding of gRNA in FnCas9 protein.

The introduction of CRISPR/Cas9 based gene editing has greatly accelerated therapeutic genome editing. However, the off-target DNA cleavage by CRISPR/Cas9 protein hampers its clinical translation, hindering its widespread use as a programmable genome editing tool. Although Cas9 variants with better mismatch discrimination have been developed, they have significantly lower rates of on-target DNA cleavage. Here, we have compared the dynamics of a more specific naturally occurring Cas9 from Francisella novicida (FnCas9) to the most widely used, SpCas9 protein. Long-scale atomistic MD simulation of free and gRNA bound forms of both the Cas9 proteins was performed, and their domain rearrangements and binding affinity with gRNA were compared to decipher the possible reason behind the enhanced specificity of FnCas9 protein. The greater binding affinity with gRNA, high domain electrostatics, and more volatility of FnCas9 than SpCas9 may explain its increased specificity and lower tolerance for mismatches.

Predicted Hotspot Residues Involved in Allosteric Signal Transmission in Pro-Apoptotic Peptide-Mcl1 Complexes

Mcl1 is a primary member of the Bcl–2 family—anti–apoptotic proteins (AAP)—that is overexpressed in several cancer pathologies. The apoptotic regulation is mediated through the binding of pro-apoptotic peptides (PAPs) (e.g., Bak and Bid) at the canonical hydrophobic binding groove (CBG) of Mcl1. Although all PAPs form amphipathic α-helices, their amino acid sequences vary to different degree. This sequence variation exhibits a central role in the binding partner selectivity towards different AAPs. Thus, constructing a novel peptide or small organic molecule with the ability to mimic the natural regulatory process of PAP is essential to inhibit various AAPs. Previously reported experimental binding free energies (BFEs) were utilized in the current investigation aimed to understand the mechanistic basis of different PAPs targeted to mMcl1. Molecular dynamics (MD) simulations used to estimate BFEs between mMcl1—PAP complexes using Molecular Mechanics-Generalized Born Solvent Accessible (MMGBSA) approach with multiple parameters. Predicted BFE values showed an excellent agreement with the experiment (R² = 0.92). The van–der Waals (ΔG_vdw) and electrostatic (ΔG_ele) energy terms found to be the main energy components that drive heterodimerization of mMcl1—PAP complexes. Finally, the dynamic network analysis predicted the allosteric signal transmission pathway involves more favorable energy contributing residues. In total, the results obtained from the current investigation may provide valuable insights for the synthesis of a novel peptide or small organic inhibitor targeting Mcl1.

Probing the Behaviour of Cas1-Cas2 upon Protospacer Binding in CRISPR-Cas Systems using Molecular Dynamics Simulations

Adaptation in CRISPR-Cas systems enables the generation of an immunological memory to defend against invading viruses. This process is driven by foreign DNA spacer (termed protospacer) selection and integration mediated by Cas1-Cas2 protein. Recently, different states of Cas1-Cas2, in its free form and in complex with protospacer DNAs, were solved by X-ray crystallography. In this paper, molecular dynamics (MD) simulations are employed to study crystal structures of one free and two protospacer-bound Cas1-Cas2 complexes. The simulated results indicate that the protospacer binding markedly increases the system stability, in particular when the protospacer containing the PAM-complementary sequence. The hydrogen bond and binding free energy calculations explain that PAM recognition introduces more specific interactions to increase the cleavage activity of Cas1. By using principal component analysis (PCA) and intramolecular angle calculation, this study observes two dominant slow motions associated with the binding of Ca1-Cas2 to the protospacer and potential target DNAs respectively. The comparison of DNA structural deformation further implies a cooperative conformational change of Cas1-Cas2 and protospacer for the target DNA capture. We propose that this cooperativity is the intrinsic requirement of the CRISPR integration complex formation. This study provides some new insights into the understanding of CRISPR-Cas adaptation.

Insight into tetrodotoxin blockade and resistance mechanisms of Nav 1.2 sodium channel by theoretical approaches

Na_v 1.2, a member of voltage-gated sodium channels (Na_v s) that are responsible for the generation and propagation of action potentials along the cell membrane, and play a vital role in the process of information transmission within the nervous system and muscle contraction, is preferentially expressed in the central nervous system. As a potent and selective blocker of Na_v s, tetrodotoxin (TTX) has been extensively studied in biological and chemical sciences, whereas the detailed mechanism by which it blocks nine Na_v 1 channel subtypes remain elusive. Despite the high structural similarity, the TTX metabolite 4,9-anhydro-TTX is 161 times less effective toward the mammalian Na_v 1.2, which puzzled us to ask a question why such a subtle structural variation results in the largely binding affinity difference. In the current work, an integrated computational strategy, including homology modeling, induced fit docking, explicit-solvent MD simulations, and free energy calculations, was employed to investigate the binding mechanism and conformational determinants of TTX analogs. Based on the computational results, the H-bond interactions between C4-OH and C9-OH of TTX and the outer ring carboxylates of the selectivity-filter residues, and the cation-π interaction between the primary amine of guanidinium of TTX and Phe385 determine the difference of their binding affinities. Moreover, the computationally simulations were carried out for the D384N and E945K mutants of hNa_v 1.2-TTX, and the rank of the predicted binding free energies is in accordance with the experimental data. These observations provide a valuable model to design potent and selective neurotoxins of Na_v 1.2 and shed light on the blocking mechanism of TTX to sodium channels.

Molecular Simulation Studies on the Binding Selectivity of Type-I Inhibitors in the Complexes with ROS1 versus ALK

ROS1 and ALK are promising targets of anticancer drugs for non-small-cell lung cancer. Since they have 49% amide acid sequence homology in the kinases domain and 77% identity at the ATP binding area, some ALK inhibitors also showed some significant responses for ROS1 in the clinical trial, such as the type-I binding inhibitor crizotinib and PF-06463922. As a newly therapeutic target, the selective ROS1 inhibitor is relatively rare. Moreover, the molecular basis for the selectivity of ROS1 versus ALK still remains unclear. In order to disclose the binding preference toward ROS1 over ALK and to aid the design of selective ROS1 inhibitors, the specific interactions and difference of conformational changes in the dual and selective ROS1/ALK inhibitors systems were investigated by molecular dynamics (MD) simulation and principle component analysis (PCA) in our work. Afterward, binding free energies (MM/GBSA) and binding free energies decomposition analysis indicated that the dominating effect of Van der Waals interaction drives the specific binding process of the type-I inhibitor, and residues of the P-loop and the DFG motif would play an important role in selectivity. On the basis of the modeling results, the new designed compound 14c was verified as a selective ROS1 inhibitor versus ALK, and SMU-B was a dual ROS1/ALK inhibitor by the kinase inhibitory study. These results are expected to facilitate the discovery and rational design of novel and specific ROS1 inhibitors.

Exploring the Mutational Robustness of Nucleic Acids by Searching Genotype Neighborhoods in Sequence Space

To assess the mutational robustness of nucleic acids, many genome- and protein-level studies have been performed, where nucleic acids are treated as genetic information carriers and transferrers. However, the molecular mechanisms through which mutations alter the structural, dynamic, and functional properties of nucleic acids are poorly understood. Here we performed a SELEX in silico study to investigate the fitness distribution of the l-Arm-binding aptamer genotype neighborhoods. Two novel functional genotype neighborhoods were isolated and experimentally verified to have comparable fitness as the wild-type. The experimental aptamer fitness landscape suggests the mutational robustness is strongly influenced by the local base environment and ligand-binding mode, whereas bases distant from the binding pocket provide potential evolutionary pathways to approach the global fitness maximum. Our work provides an example of successful application of SELEX in silico to optimize an aptamer and demonstrates the strong sensitivity of mutational robustness to the site of genetic variation.

Metal binding mediated conformational change of XPA protein:a potential cytotoxic mechanism of nickel in the nucleotide excision repair

Nucleotide excision repair (NER) is a pivotal life process for repairing DNA nucleotide mismatch caused by chemicals, metal ions, radiation, and other factors. As the initiation step of NER, the xeroderma pigmentosum complementation group A protein (XPA) recognizes damaged DNA molecules, and recruits the replication protein A (RPA), another important player in the NER process. The stability of the Zn(2+)-chelated Zn-finger domain of XPA center core portion (i.e., XPA98-210) is the foundation of its biological functionality, while the displacement of the Zn(2+) by toxic metal ions (such as Ni(2+), a known human carcinogen and allergen) may impair the effectiveness of NER and hence elevate the chance of carcinogenesis. In this study, we first calculated the force field parameters for the bonded model in the metal center of the XPA98-210 system, showing that the calculated results, including charges, bonds, angles etc., are congruent with previously reported results measured by spectrometry experiments and quantum chemistry computation. Then, comparative molecular dynamics simulations using these parameters revealed the changes in the conformation and motion mode of XPA98-210 Zn-finger after the substitution of Zn(2+) by Ni(2+). The results showed that Ni(2+) dramatically disrupted the relative positions of the four Cys residues in the Zn-finger structure, forcing them to collapse from a tetrahedron into an almost planar structure. Finally, we acquired the binding mode of XPA98-210 with its ligands RPA70N and DNA based on molecular docking and structural alignment. We found that XPA98-210's Zn-finger domain primarily binds to a V-shaped cleft in RPA70N, while the cationic band in its C-terminal subdomain participates in the recognition of damaged DNA. In addition, this article sheds light on the multi-component interaction pattern among XPA, DNA, and other NER-related proteins (i.e., RPA70N, RPA70A, RPA70B, RPA70C, RPA32, and RPA14) based on previously reported structural biology information. Thus, we derived a putative cytotoxic mechanism associated with the nickel ion, where the Ni(2+) disrupts the conformation of the XPA Zn-finger, directly weakening its interaction with RPA70N, and thus lowering the effectiveness of the NER process. In sum, this work not only provides a theoretical insight into the multi-protein interactions involved in the NER process and potential cytotoxic mechanism associated with Ni(2+) binding in XPA, but may also facilitate rational anti-cancer drug design based on the NER mechanism.

Unraveling the conformational determinants of LARP7 and 7SK small nuclear RNA by theoretical approaches

LARP7, a member of the La-related proteins (LARPs), shares a conserved La module comprising the La-motif (LAM) and the RNA-recognition motif (RRM1), binding exclusively to the non-coding RNA 7SK.

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.