Structural Mechanism behind Distinct Efficiency of Oct4/Sox2 Proteins in Differentially Spaced DNA Complexes

Implementing accelerated dynamics to unravel the effects of high-fidelity Cas9 mutants on target DNA and guide RNA hybrid stability

The clustered regularly interspersed short palindromic repeats (CRISPR) and its associated nuclease (Cas9) offers a unique and easily reprogrammable system for editing eukaryotic genomes. Cas9 is guided to the target by an RNA strand, and precise edits are created by introducing double-stranded breaks. However, nuclease activity of Cas9 is also triggered at other sites other than the target sit, which is a major limitation for various applications. Cas9 variants have been designed to improve the efficacy of the tool by introducing certain mutations. However, the on-target activity of such Cas9 variants is often seen as compromised. Hence, understanding the sub-molecular differences in the variants is essential to elucidate the factors that contribute to efficiency. The study reveals distortions in the PAM-distal regions of the nucleic hybrids as well as changes in the interactions between the Cas9 variants and RNA-DNA hybrid, contributing to the explanation for differences in on-target activity.Communicated by Ramaswamy H. Sarma.

In silico identification of natural product inhibitors against Octamer-binding transcription factor 4 (Oct4) to impede the mechanism of glioma stem cells

Octamer-binding transcription factor 4 (Oct4) is a core regulator in the retention of stemness, invasive, and self-renewal properties in glioma initiating cells (GSCs) and its overexpression inhibits the differentiation of glioma cells promoting tumor cell proliferation. The Pit-Oct-Unc (POU) domain comprising POU-specific domain (POUS) and POU-type homeodomain (POUHD) subdomains is the most critical part of the Oct4 for the generation of induced pluripotent stem cells from somatic cells that lead to tumor initiation, invasion, posttreatment relapse, and therapeutic resistance. Therefore, the present investigation hunts for natural product inhibitors (NPIs) against the POUHD domain of Oct4 by employing receptor-based virtual screening (RBVS) followed by binding free energy calculation and molecular dynamics simulation (MDS). RBVS provided 13 compounds with acceptable ranges of pharmacokinetic properties and good docking scores having key interactions with the POUHD domain. More Specifically, conformational and interaction stability analysis of 13 compounds through MDS unveiled two compounds ZINC02145000 and ZINC32124203 which stabilized the backbone of protein even in the presence of linker and POUS domain. Additionally, ZINC02145000 and ZINC32124203 exhibited stable and strong interactions with key residues W277, R242, and R234 of the POUHD domain even in dynamic conditions. Interestingly, ZINC02145000 and ZINC32124203 established communication not only with the POUHD domain but also with the POUS domain indicating their incredible potency toward thwarting the function of Oct4. ZINC02145000 and ZINC32124203 also reduced the flexibility and escalated the correlations between the amino acid residues of Oct4 evidenced by PCA and DCCM analysis. Finally, our examination proposed two NPIs that can impede the Oct4 function and may help to improve overall survival, diminish tumor relapse, and achieve a cure not only in deadly disease GBM but also in other cancers with minimal side effects.

Computational Insights into the Deleterious Impacts of Missense Variants on N-Acetyl-d-glucosamine Kinase Structure and Function

An enzyme of the mammalian amino-sugar metabolism pathway, N-acetylglucosamine kinase (NAGK), that synthesizes N-acetylglucosamine (GlcNAc)-6-phosphate, is reported to promote dynein functions during mitosis, axonal and dendritic growth, cell migration, and selective autophagy, which all are unrelated to its enzyme activity. As non-enzymatic structural functions can be altered by genetic variation, we made an effort in this study aimed at deciphering the pathological effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in NAGK gene. An integrated computational approach, including molecular dynamics (MD) simulation and protein-protein docking simulation, was used to identify the damaging nsSNPs and their detailed structural and functional consequences. The analysis revealed the four most damaging variants (G11R, G32R, G120E, and A156D), which are highly conserved and functional, positioned in both small (G11R and G32R) and large (G120E and A156D) domains of NAGK. G11R is located in the ATP binding region, while variants present in the large domain (G120E and A156D) were found to induce substantial alterations in the structural organizations of both domains, including the ATP and substrate binding sites. Furthermore, all variants were found to reduce binding energy between NAGK and dynein subunit DYNLRB1, as revealed by protein-protein docking and MM-GBSA binding energy calculation supporting their deleteriousness on non-canonical function. We hope these findings will direct future studies to gain more insight into the role of these variants in the loss of NAGK function and their role in neurodevelopmental disorders.

Molecular Mechanism of the Intercalation of the SOX-4 Protein into DNA Inducing Bends and Kinks

DNA-protein interactions regulate several biophysical functions, yet the mechanism of only a few is investigated in molecular detail. An important example is the intercalation of transcription factor proteins into DNA that produce bent and kinked DNA. Here, we have studied the molecular mechanism of the intercalation of a transcription factor SOX4 into DNA with a goal to understand the sequence of molecular events that precede the bending and kinking of the DNA. Our long well-tempered metadynamics and molecular dynamics (MD) simulations show that the protein primarily binds to the backbone of DNA and rotates around it to form an intercalative native state. We show that although there are multiple pathways for intercalation, the deintercalation pathway matches with the most probable intercalation pathway. In both cases, bending and kinking happen simultaneously, driven by the onset of the intercalation of the amino acid.

Prevention of SARS-CoV-2 Proliferation with a Novel and Potent Main Protease Inhibitor by Docking, ADMET, MM-PBSA, and Molecular Dynamics Simulation

COVID-19 is the last disease caused by SARS-CoV-2 associated with a severe immune response and lung damage. The main protease (Mpro) has a vital role in SARS-CoV-2 proliferation. Moreover, humans lack homologous Mpro, which makes the Mpro a suitable drug target for the development of SARS-CoV-2 drugs. The purchasable L5000 library (Selleckchem Inc) includes 99,040 compounds that were used for virtual screening. After molecular docking and ADME studies, we selected a compound (WAY-604395) with a potent binding affinity to the Mpro active site and acceptable ADME properties compared to the reference drug (nelfinavir). Molecular dynamics (MD) simulation outcomes have proved that the Mpro-WAY604395 complex possesses a considerable value of flexibility, stability, compactness and binding energy. Our Molecular Mechanics Poisson–Boltzmann Surface Area (MM-PBSA) calculation demonstrates that WAY-604395 is more potent ([Formula: see text]272.19[Formula: see text]kcal mol[Formula: see text]) in comparison with nelfinavir ([Formula: see text]173.39[Formula: see text]kcal[Formula: see text]mol[Formula: see text]) against SARS-CoV-2 Mpro. In conclusion, we suggest that WAY-604395 has the potential for the treatment of SARS-CoV-2 by inhibition of the Mpro.

Whole proteome screening and identification of potential epitopes of SARS-CoV-2 for vaccine design-an immunoinformatic, molecular docking and molecular dynamics simulation accelerated robust strategy

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the most cryptic pandemic outbreak of the 21st century, has gripped more than 1.8 million people to death and infected almost eighty six million. As it is a new variant of SARS, there is no approved drug or vaccine available against this virus. This study aims to predict some promising cytotoxic T lymphocyte epitopes in the SARS-CoV-2 proteome utilizing immunoinformatic approaches. Firstly, we identified 21 epitopes from 7 different proteins of SARS-CoV-2 inducing immune response and checked for allergenicity and conservancy. Based on these factors, we selected the top three epitopes, namely KAYNVTQAF, ATSRTLSYY, and LTALRLCAY showing functional interactions with the maximum number of MHC alleles and no allergenicity. Secondly, the 3D model of selected epitopes and HLA-A*29:02 were built and Molecular Docking simulation was performed. Most interestingly, the best two epitopes predicted by docking are part of two different structural proteins of SARS-CoV-2, namely Membrane Glycoprotein (ATSRTLSYY) and Nucleocapsid Phosphoprotein (KAYNVTQAF), which are generally target of choice for vaccine designing. Upon Molecular Docking, interactions between selected epitopes and HLA-A*29:02 were further validated by 50 ns Molecular Dynamics (MD) simulation. Analysis of RMSD, Rg, SASA, number of hydrogen bonds, RMSF, MM-PBSA, PCA, and DCCM from MD suggested that ATSRTLSYY is the most stable and promising epitope than KAYNVTQAF epitope. Moreover, we also identified B-cell epitopes for each of the antigenic proteins of SARS CoV-2. Findings of our work will be a good resource for wet lab experiments and will lessen the timeline for vaccine construction.Communicated by Ramaswamy H. Sarma.

The SARS-Cov-2 Proliferation Blocked by a Novel and Potent Main Protease Inhibitor via Computer-aided Drug Design

The recent prevalence of novel "coronavirus disease 2019" has expanded quickly globally, causing a universal pandemic. Herein, an effort was constructed to design a potent drug to inhibit the main protease of SARS-Cov-2 (3CLp) by means of structure-based drug design. A large library of the compounds was used for virtual screening. After molecular docking and ADME studies, we selected a compound with a better binding affinity to the 3CLp active site and acceptable ADME properties compared to the selected positive control drug. Molecular dynamic (MD) simulation (200 ns) and Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) were used for further analysis. MD simulation outcomes have proved that the 3CLp-ZINC31157475 complex possesses a considerable value of dynamic properties such as flexibility, stability, compactness, and binding energy. Our MM-PBSA computation illustrates that ZINC31157475 is more potent (-88.03 kcal mol-1) than nelfinavir (-19.54 kcal mol-1) against COVID-19 3CLp. Further, we have determined that the main residues of the 3CLp interact with ligands from per-residue binding energy. In conclusion, we suggest that ZINC31157475 can potentially treat COVID-19 by inhibition of the 3CLp. However, in-vitro and in-vivo study is essential for approval of this suggestion.

Exploring the effect of aplidin on low molecular weight protein tyrosine phosphatase by molecular docking and molecular dynamic simulation study

The low molecular weight protein tyrosine phosphatase (LMW-PTP) could regulate many signaling pathways, and it had drawn attention as a potential target for cancer. As previous report has indicated that the aplidin could inhibit the LMW-PTP, and thus, the relevant cancer caused by the abnormal regulation of the LMW-PTP could be remission. However, the molecular mechanism of inhibition of the LMW-PTP by the aplidin had not been fully understood. In this study, various computational approaches, namely molecular docking, MDs and post-dynamic analyses were utilized to explore the effect of the aplidin on the LMW-PTP. The results suggested that the intramolecular interactions of the residues in the two sides of the active site (Ser43-Ala55 and Pro121-Asn134) and the P-loop region (Leu13-Ser19) in the LMW-PTP was disturbed owing to the aplidin, meanwhile, the π-π interaction between Tyr131 and Tyr132 might be broken. The Asn15 might be the key residue to break the residues interactions. In a word, this study may provide more information for understanding the effect of inhibition of the aplidin on the LMW-PTP.

Design, synthesis, biological evaluation and molecular dynamics studies of 4-thiazolinone derivatives as protein tyrosine phosphatase 1B (PTP1B) inhibitors

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of insulin signaling pathway, and more and more studies have shown that it is a potential target for the treatment of type 2 diabetes mellitus (T2DM). In this study, 17 new 4-thiazolinone derivatives were designed and synthesized as novel PTP1B inhibitors, and ADMET prediction confirmed that these compounds were to be drug-like. In vitro enzyme activity experiments were performed on these compounds, and it was found that a plurality of compounds had good inhibitory activity and high selectivity against PTP1B protein. Among them, compound 7p exhibited the best inhibitory activity with an IC₅₀ of 0.92 μM. The binding mode of compound 7p and PTP1B protein was explored, revealing the reason for its high efficiency. In addition, molecular dynamics simulations for the PTP1B^WT and PTP1B^comp#7p systems revealed the effects of compound 7p on PTP1B protein at the molecular level. In summary, the study reported for the first time that 4-thiazolinone derivatives as a novel PTP1B inhibitor had good inhibitory activity and selectivity for the treatment of T2DM, providing more options for the development of PTP1B inhibitors. AbbreviationsBBBblood-brain barrierCDC25Bcell division cycle 25 homolog BCYP2D6Cytochrome P450 2D6 bindingDCCMdynamic cross-correlation mapDSDiscovery StudioH bondhydrogen bondHIAhuman intestinal absorptionLARleukocyte antigen-related phosphataseMDmolecular dynamicsMEG-2maternal-effect germ-cell defective 2MM-PBSAmolecular mechanics Poisson Boltzmann surface area)PCAprincipal component analysisPDBProtein Data BankpNPPp-nitrophenyl phosphatePPBplasma protein bindingPTP1Bprotein tyrosine phosphotase 1BRMSDroot mean square deviationRMSFroot mean square fluctuationSHP-1src homologous phosphatase-1SHP-2src homologous phosphatase-2SPCsingle-point chargeTCPTPT cell protein tyrosine phosphataseT2DMType 2 diabetes mellitusVDWvan der WaalsCommunicated by Ramaswamy H. Sarma.

Exploring the effect of inhibitor AKB-9778 on VE-PTP by molecular docking and molecular dynamics simulation

Diabetic macular edema, also known as diabetic eye disease, is mainly caused by the overexpression of vascular endothelial protein tyrosine phosphatase (VE-PTP) at hypoxia/ischemic. AKB-9778 is a known VE-PTP inhibitor that can effectively interact with the active site of VE-PTP to inhibit the activity of VE-PTP. However, the binding pattern of VE-PTP with AKB-9778 and the dynamic implications of AKB-9778 on VE-PTP system at the molecular level are poorly understood. Through molecular docking, it was found that the AKB-9778 was docked well in the binding pocket of VE-PTP by the interactions of hydrogen bond and Van der Waals. Furthermore, after molecular dynamic simulations on VE-PTP system and VE-PTP ^AKB-9778 system, a series of postdynamic analyses found that the flexibility and conformation of the active site undergone an obvious transition after VE-PTP binding with AKB-9778. Moreover, by constructing the RIN, it was found that the different interactions in the active site were the detailed reasons for the conformational differences between these two systems. Thus, the finding here might provide a deeper understanding of AKB-9778 as VE-PTP Inhibitor.

Structural mechanism of DNA-mediated Nanog–Sox2 cooperative interaction

The efficiency of stem cell transcriptional regulation always depends on the cooperative association and expression of transcription factors (TFs). Among these, Oct4, Sox2, and Nanog play major roles. Their cooperativity is facilitated via direct protein–protein interactions or DNA-mediated interactions, yet the mechanism is not clear. Most biochemical studies have examined Oct4/Sox2 cooperativity, whereas few studies have evaluated how Nanog competes in the connection between these TFs. In this study, using computational models and molecular dynamics simulations, we built a framework representing the DNA-mediated cooperative interaction between Nanog and Sox2 and analyzed the plausible interaction factors experienced by Nanog because of Sox2, its cooperative binding partner. Comparison of a wild-type and mutant Nanog/Sox2 model with the Nanog crystal structure revealed the regulatory structural mechanism between Nanog/Sox2–DNA-mediated cooperative bindings. Along with the transactivation domains interaction, the DNA-mediated allosteric interactions are also necessary for Nanog cooperative binding. DNA-mediated Nanog–Sox2 cooperativity influences the protein conformational changes and a stronger interaction profile was observed for Nanog-Mut (L103E) in comparison with the Nanog-WT complex.

The efficiency of stem cell transcriptional regulation always depends on the cooperative association and expression of transcription factors (TFs).

Exploring the effect of N308D mutation on protein tyrosine phosphatase-2 cause gain-of-function activity by a molecular dynamics study

One of the most common protein tyrosine phosphatase-2 (SHP2) mutations in Noonan syndrome is the N308D mutation, and it increases the activity of the protein. However, the molecular basis of the activation of N308D mutation on SHP2 conformations is poorly understood. Here, molecular dynamic simulations were performed on SHP2 and SHP2-N308D to explore the effect of N308D mutation on SHP2 cause gain of function activity, respectively. The principal component analysis, dynamic cross-correlation map, secondary structure analysis, residue interaction networks, and solvent accessible surface area analysis suggested that the N308D mutation distorted the residues interactions network between the allosteric site (residue Gly244-Gly246) and C-SH2 domain, including the hydrogen bond formation and the binding energy. Meanwhile, the activity of catalytic site (residue Gly503-Val505) located in the Q-loop in mutant increased due to this region's high fluctuations. Therefore, the substrate had more chances to access to the catalytic activity site of the precision time protocol domain of SHP2-N308D, which was easy to be exposed. In addition, we had speculated that the Lys244 located in the allosteric site was the key residue which lead to the protein conformation changes. Consequently, overall calculations presented in this study ultimately provide a useful understanding of the increased activity of SHP2 caused by the N308D mutation.

Evaluation of Sox2 binding affinities for distinct DNA patterns using steered molecular dynamics simulation

Transcription factors (TFs) are gene expression regulators that bind to DNA in a sequence‐specific manner and determine the functional characteristics of the gene. It is worthwhile to study the unique characteristics of such specific TF‐binding pattern in DNA. Sox2 recognizes a 6‐ to 7‐base pair consensus DNA sequence; the central four bases of the binding site are highly conserved, whereas the two to three flanking bases are variable. Here, we attempted to analyze the binding affinity and specificity of the Sox2 protein for distinct DNA sequence patterns via steered molecular dynamics, in which a pulling force is employed to dissociate Sox2 from Sox2–DNA during simulation to study the behavior of a complex under nonequilibrium conditions. The simulation results revealed that the first two stacking bases of the binding pattern have an exclusive impact on the binding affinity, with the corresponding mutant complexes showing greater binding and longer dissociation time than the experimental complexes do. In contrast, mutation of the conserved bases tends to reduce the affinity, and mutation of the complete conserved region disrupts the binding. It might pave the way to identify the most likely binding pattern recognized by Sox2 based on the affinity of each configuration. The α2‐helix of Sox2 was found to be the key player in the Sox2–DNA association. The characterization of Sox2's binding patterns for the target genes in the genome helps in understanding of its regulatory functions.

Structural and conformational insights into SOX2/OCT4-bound enhancer DNA: a computational perspective

The roles of SOX2 and OCT4 are critical in stem cell maintenance either in the context of iPSCs generation or cancer stem cell growth; therefore, it is imperative to study their cooperative binding and SOX2/OCT4-induced DNA conformational switching.