Cooperative binding of DNA and CBFbeta to the Runt domain of the CBFalpha studied via MD simulations

Proteogenomic analysis reveals cytoplasmic sequestration of RUNX1 by the acute myeloid leukemia-initiating CBFB::MYH11 oncofusion protein

Several canonical translocations produce oncofusion genes that can initiate acute myeloid leukemia (AML). Although each translocation is associated with unique features, the mechanisms responsible remain unclear. While proteins interacting with each oncofusion are known to be relevant for how they act, these interactions have not yet been systematically defined. To address this issue in an unbiased fashion, we fused a promiscuous biotin ligase (TurboID) in-frame with 3 favorable-risk AML oncofusion cDNAs (PML::RARA, RUNX1::RUNX1T1, and CBFB::MYH11) and identified their interacting proteins in primary murine hematopoietic cells. The PML::RARA- and RUNX1::RUNX1T1-TurboID fusion proteins labeled common and unique nuclear repressor complexes, implying their nuclear localization. However, CBFB::MYH11-TurboID-interacting proteins were largely cytoplasmic, probably because of an interaction of the MYH11 domain with several cytoplasmic myosin-related proteins. Using a variety of methods, we showed that the CBFB domain of CBFB::MYH11 sequesters RUNX1 in cytoplasmic aggregates; these findings were confirmed in primary human AML cells. Paradoxically, CBFB::MYH11 expression was associated with increased RUNX1/2 expression, suggesting the presence of a sensor for reduced functional RUNX1 protein, and a feedback loop that may attempt to compensate by increasing RUNX1/2 transcription. These findings may have broad implications for AML pathogenesis.

Super-sensitive bifunctional nanoprobe: Self-assembly of peptide-driven nanoparticles demonstrating tumor fluorescence imaging and therapy

The development of nanomedicine has recently achieved several breakthroughs in the field of cancer treatment; however, biocompatibility and targeted penetration of these nanomaterials remain as limitations, which lead to serious side effects and significantly narrow the scope of their application. The self-assembly of intermediate filaments with arginine-glycine-aspartate (RGD) peptide (RGD-IFP) was triggered by the hydrophobic cationic molecule 7-amino actinomycin D (7-AAD) to synthesize a bifunctional nanoparticle that could serve as a fluorescent imaging probe to visualize tumor treatment. The designed RGD-IFP peptide possessed the ability to encapsulate 7-AAD molecules through the formation of hydrogen bonds and hydrophobic interactions by a one-step method. This fluorescent nanoprobe with RGD peptide could be targeted for delivery into tumor cells and released in acidic environments such as endosomes/lysosomes, ultimately inducing cytotoxicity by arresting tumor cell cycling with inserted DNA. It is noteworthy that the RGD-IFP/7-AAD nanoprobe tail-vein injection approach demonstrated not only high tumor-targeted imaging potential, but also potent antitumor therapeutic effects in vivo. The proposed strategy may be used in peptide-driven bifunctional nanoparticles for precise imaging and cancer therapy.

Inclusion of methoxy groups inverts the thermodynamic stabilities of DNA-RNA hybrid duplexes: A molecular dynamics simulation study

Modified nucleic acids have found profound applications in nucleic acid based technologies such as antisense and antiviral therapies. Previous studies on chemically modified nucleic acids have suggested that modifications incorporated in furanose sugar especially at 2'-position attribute special properties to nucleic acids when compared to other modifications. 2'-O-methyl modification to deoxyribose sugars of DNA-RNA hybrids is one such modification that increases nucleic acid stability and has become an attractive class of compounds for potential antisense applications. It has been reported that modification of DNA strands with 2'-O-methyl group reverses the thermodynamic stability of DNA-RNA hybrid duplexes. Molecular dynamics simulations have been performed on two hybrid duplexes (DR and RD) which differ from each other and 2'-O-methyl modified counterparts to investigate the effect of 2'-O-methyl modification on their duplex stability. The results obtained suggest that the modification drives the conformations of both the hybrid duplexes towards A-RNA like conformation. The modified hybrid duplexes exhibit significantly contrasting dynamics and hydration patterns compared to respective parent duplexes. In line with the experimental results, the relative binding free energies suggest that the introduced modifications stabilize the less stable DR hybrid, but destabilize the more stable RD duplex. Binding free energy calculations suggest that the increased hydrophobicity is primarily responsible for the reversal of thermodynamic stability of hybrid duplexes. Free energy component analysis further provides insights into the stability of modified duplexes.

A wrench in the works of human acetylcholinesterase: soman induced conformational changes revealed by molecular dynamics simulations

Irreversible inactivation of human acetylcholinesterase (hAChE) by organophosphorous pesticides (OPs) and chemical weapon agents (CWA) has severe morbidity and mortality consequences. We present data from quantum mechanics/molecular mechanics (QM/MM) and 80 classical molecular dynamics (MD) simulations of the apo and soman-adducted forms of hAChE to investigate the effects on the dynamics and protein structure when the catalytic Serine 203 is phosphonylated. We find that the soman phosphonylation of the active site Ser203 follows a water assisted addition-elimination mechanism with the elimination of the fluoride ion being the highest energy barrier at 6.5 kcal/mole. We observe soman-dependent changes in backbone and sidechain motions compared to the apo form of the protein. These alterations restrict the soman-adducted hAChE to a structural state that is primed for the soman adduct to be cleaved and removed from the active site. The altered motions and resulting structures provide alternative pathways into and out of the hAChE active site. In the soman-adducted protein both side and back door pathways are viable for soman adduct access. Correlation analysis of the apo and soman adducted MD trajectories shows that the correlation of gorge entrance and back door motion is disrupted when hAChE is adducted. This supports the hypothesis that substrate and product can use two different pathways as entry and exit sites in the apo form of the protein. These alternative pathways have important implications for the rational design of medical countermeasures.

Understanding the recognition mechanisms of Zα domain of human editing enzyme ADAR1 (hZα(ADAR1)) and various Z-DNAs from molecular dynamics simulation

The Z-DNA-binding domain of human double-stranded RNA adenosine deaminase I (hZαADAR1) can specifically recognize the left-handed Z-DNA which preferentially occurs at alternating purine-pyrimidine repeats, especially the CG-repeats. The interactions of hZαADAR1 and Z-DNAs in different sequence contexts can affect many important biological functions including gene regulation and chromatin remodeling. Therefore it is of great necessity to fully understand their recognition mechanisms. However, most existing studies are aimed at the standard CG-repeat Z-DNA rather than the non-CG-repeats, and whether the molecular basis of hZαADAR1 binding to various Z-DNAs are identical or not is still unclear on the atomic level. Here, based on the recently determined crystal structures of three representative non-CG-repeat Z-DNAs (d(CACGTG)2, d(CGTACG)2 and d(CGGCCG)2) in complex with hZαADAR1, 40 ns molecular dynamics simulation together with binding free energy calculation were performed for each system. For comparison, the standard CG-repeat Z-DNA (d(CGCGCG)2) complexed with hZαADAR1 was also simulated. The consistent results demonstrate that nonpolar interaction is the driving force during the protein-DNA binding process, and that polar interaction mainly from helix α3 also provides important contributions. Five common hot-spot residues were identified, namely Lys169, Lys170, Asn173, Arg174 and Tyr177. Hydrogen bond analysis coupled with surface charge distribution further reveal the interfacial information between hZαADAR1 and Z-DNA in detail. All of the analysis illustrate that four complexes share the common key features and the similar binding modes irrespective of Z-DNA sequences, suggesting that Z-DNA recognition by hZαADAR1 is conformation-specific rather than sequence-specific. Additionally, by analyzing the conformational changes of hZαADAR1, we found that the binding of Z-DNA could effectively stabilize hZαADAR1 protein. Our study can provide some valuable information for better understanding the binding mechanism between hZαADAR1 or even other Z-DNA-binding protein and Z-DNA.

RUNX3 is inactivated by promoter hypermethylation in malignant transformation of ovarian endometriosis

The aim of the present study was to investigate the role of epigenetic inactivation of the runt-related transcription factor 3 gene (RUNX3) in the malignant transformation of ovarian endometriosis. Samples obtained by microdissection and scraping included 30 malignant ovarian endometriotic cyst tissues and 30 corresponding eutopic endometrium tissues from the endometriosis-associated ovarian carcinoma (EAOC) group, 19 benign ovarian endometriotic cyst tissues and 22 corresponding eutopic endometrium tissues from the endometriosis (EM) group and 22 normal eutopic endometrium tissues from the control endometrium (CE) group. RUNX3 methylation status was determined by methylation-specific PCR and bisulfite sequencing, while levels of RUNX3 and ERα protein expression were evaluated using immunohistochemistry. The percentage of RUNX3 methylation and negative RUNX3 protein expression in the malignant ovarian endometriotic cysts from the EAOC group was significantly higher than that in the benign ovarian endometriotic cysts from the EM group. The percentage of RUNX3 methylation and negative RUNX3 protein expression in the eutopic endometrium from the EAOC group was significantly higher than that in the EM and CE groups. An inverse correlation between positive RUNX3 protein expression and methylation was observed and a positive correlation was shown between RUNX3 methylation and ERα protein expression. In the malignant ovarian endometriotic cysts from the EAOC group, there was no significant correlation between methylation frequency of the RUNX3 gene and histological type. However, the percentage of RUNX3 gene methylation was significantly higher in the tissue samples from patients with surgical stage IC EAOC than the percentage in patients with stage IA and IB disease. These results suggest that RUNX3 inactivation by promoter hypermethylation plays a role in the progression of malignant transformation of ovarian EM and is closely related to estrogen metabolism. Negative protein expression and abnormal RUNX3 methylation in the eutopic endometrium could be used as diagnostic markers in patients with ovarian EM who may be at an increased risk of developing EAOC.

Structures, dynamics, and stabilities of fully modified locked nucleic acid (β-D-LNA and α-L-LNA) duplexes in comparison to pure DNA and RNA duplexes

Locked nucleic acid (LNA) is a chemical modification which introduces a -O-CH2- linkage in the furanose sugar of nucleic acids and blocks its conformation in a particular state. Two types of modifications, namely, 2'-O,4'-C-methylene-β-D-ribofuranose (β-D-LNA) and 2'-O,4'-C-methylene-α-L-ribofuranose (α-L-LNA), have been shown to yield RNA and DNA duplex-like structures, respectively. LNA modifications lead to increased melting temperatures of DNA and RNA duplexes, and have been suggested as potential therapeutic agents in antisense therapy. In this study, molecular dynamics (MD) simulations were performed on fully modified LNA duplexes and pure DNA and RNA duplexes sharing a similar sequence to investigate their structure, stabilities, and solvation properties. Both LNA duplexes undergo unwinding of the helical structure compared to the pure DNA and RNA duplexes. Though the α-LNA substituent has been proposed to mimic deoxyribose sugar in its conformational properties, the fully modified duplex was found to exhibit unique structural and dynamic properties with respect to the other three nucleic acid structures. Free energy calculations accurately capture the enhanced stabilization of the LNA duplex structures compared to DNA and RNA molecules as observed in experiments. π-stacking interaction between bases from complementary strands is shown to be one of the contributors to enhanced stabilization upon LNA substitution. A combination of two factors, namely, nature of the -O-CH2- linkage in the LNAs vs their absence in the pure duplexes and similar conformations of the sugar rings in DNA and α-LNA vs the other two, is suggested to contribute to the stark differences among the four duplexes studied here in terms of their structural, dynamic, and energetic properties.

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480 K) and MM-PBSA at 300 and 360 K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.

Molecular simulations on the thermal stabilization of DNA by hyperthermophilic chromatin protein Sac7d, and associated conformational transitions

Sac7d belongs to a family of chromosomal proteins, which are crucial for thermal stabilization of DNA at higher growth temperatures. It is capable of binding DNA nonspecifically, and is responsible for the increase in the melting temperature of DNA in the bound form up to 85 °C. Molecular dynamics (MD) simulations were performed at different temperatures on two protein-DNA complexes of Sac7d. Various structural and energetic parameters were calculated to examine the DNA stability and to investigate the conformational changes in DNA and the protein-DNA interactions. Room temperature simulations indicated very good agreement with the experimental structures. The protein structure is nearly unchanged at both 300 and 360 K, and only up to five base pairs of the DNA are stabilized by Sac7d at 360 K. However, the MD simulations on DNA alone systems show that they lose their helical structures at 360 K further supporting the role of Sac7d in stabilizing the oligomers. At higher temperatures (420 and 480 K), DNA undergoes denaturation in the presence and the absence of the protein. The DNA molecules were found to undergo B- to A-form transitions consistent with experimental studies, and the extent of these transitions are examined in detail. The extent of sampling B- and A-form regions was found to show temperature and sequence dependence. Multiple MD simulations yielded similar results validating the proposed model. Interaction energy calculations corresponding to protein-DNA binding indicates major contribution due to DNA backbone, explaining the nonspecific interactions of Sac7d.

Cooperativity and specificity of Cys2His2 zinc finger protein-DNA interactions: a molecular dynamics simulation study

Cys(2)His(2) zinc finger proteins are one of the most frequently observed DNA-binding motifs in eukaryotes. They have been widely used as a framework for designing new DNA-binding proteins. In this work, the binding affinity and conformational change of the Zif268-DNA complex were successfully reproduced with MD simulations and MM-PBSA analysis. The following new discoveries on the zinc finger protein-DNA interactions were obtained by careful energy decomposition analysis. First, a dramatic increase in the binding affinity was observed when the third zinc finger is added, indicating a cooperative nature. This cooperativity is shown to be a consequence of the small but distinctive conformational change of DNA, which enables a tight fit of the protein into the major groove of DNA. Second, specificity of the amino acid-nucleotide recognitions observed in the crystal structure is explained as originating from the ability of specific side chains and bases to take the optimal geometries for favorable interactions between polar groups. The success of the current approach implies that similar methods could be further applied to the study of protein-DNA interactions involving longer polyfingers or different linkers between fingers to provide insights for design of novel zinc finger proteins.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

Relationship between tumor suppressor gene RUNX3 expression and cell proliferation and apoptosis in colonic cancer cell line Lovo

AIM: To investigate the expression and methylation status of tumor suppressor gene RUNX3 in human colon cancer cell line Lovo and explore the effects of 5-aza-2'-deoxycytidine (5-Aza-CdR) on the proliferation and apoptosis of Lovo cells and the expression of RUNX3 gene.

METHODS: Human colon cancer cell line Lovo was treated with 5-Aza-CdR, a specific methyltransferase inhibitor, at the concentrations of 0.4, 4 and 40 μmol/L for 3 d, and then cultured in RPMI 1640 medium for 5 d. The activation of Lovo cells was respectively observed by Tetrazolium salt colorimetric (MTT) assay before and after 5-Aza-CdR treatment. The change in expression of RUNX3 mRNA was observed by semi-quantitative reverse transcription-polymerase chain reaction (RT-PCR). The apoptosis was analyzed by flow cytometry. The methylation status of gene promoter was determined by methylation-specific PCR (MSP).

RESULTS: Lovo cells treated with 5-Aza-CdR (0.4, 4, 40 μmol/L) displayed a slowed growth rate in different degrees in contrast with those in the control group and their growth rates decreased accordingly with the increase of 5-Aza-CdR concentration. There were significant increases in RUNX3 mRNA expression (0.46 ± 0.06, 0.71 ± 0.06, 0.84 ± 0.07 vs 0, P < 0.01) and apoptotic rates of Lovo cells (10.95% ± 2.09%, 17.61% ± 1.51%, 26.60% ± 1.89% vs 2.92% ± 0.93%, P < 0.01) after 5-Aza-CdR treatment in comparison with those in the control group. The level of RUNX3 mRNA expression and the apoptotic rates of Lovo cells were increased in correlation with 5-Aza-CdR concentration (F = 168.4, F = 145.7, P < 0.01). Methylation of RUNX3 promoter region was confirmed in Lovo cells of control group and detected partly in 5-Aza-CdR-treated group.

CONCLUSION: 5-Aza-CdR is able to reverse the methylation status of RUNX3 promoter region. The re-expression of RUNX3 gene can inhibit Lovo cell growth and partly induce Lovo cell apoptosis.

Molecular dynamics simulations of nucleic acid-protein complexes

Molecular dynamics simulation studies of protein-nucleic acid complexes are more complicated than studies of either component alone-the force field has to be properly balanced, the systems tend to become very large, and a careful treatment of solvent and of electrostatic interactions is necessary. Recent investigations into several protein-DNA and protein-RNA systems have shown the feasibility of the simulation approach, yielding results of biological interest not readily accessible to experimental methods.

RUN-CBFβ Interaction inC. elegans: Computational Prediction and Experimental Verification

The Runt domain proteins are eukaryotic transcription factors that regulate major developmental pathways. All members of this family contain a highly-conserved sequence-specific DNA binding domain: the Runt domain (RD). Structural and biochemical studies have shown that the Runt domain undergoes a conformational transition upon binding to DNA and that this process is regulated by an unrelated partner protein CBFbeta that enhances the DNA binding affinity of RD. Most of the reported studies on the Runt domain transcription factors were performed on proteins from mammals and Drosophila whereas very little has been known about the C. elegans RD protein, RUN, which provides the simplest model system for understanding the function of this class of transcription factors. We performed computational studies on RD domains from various species including C. elegans, Drosophila, and human, using the atom-atom contact surface area scoring method. The scoring analysis indicates that the DNA binding regulation of the C. elegans RD protein (CeRD) occurs via its interaction with a CBFbeta-like partner, as found for the human proteins, whereas a different mode of regulation may occur in the Drosophila system. Sequence, secondary structure and fold analyses of a putative CBFbeta protein identified in the C. elegans genome, CeCBFbeta, sharing a 22% identity with the human protein, predict a similar structure of this protein to that of the human CBFbeta protein. We produced the C. elegans proteins CeRD and CeCBFbeta in bacteria and confirmed their physical interaction as well as cross interactions with the corresponding human proteins. We also confirmed the structural similarity of CBFbeta and CeCBFbeta by circular dichroism analysis. The combined results suggest that a similar mechanism of regulation operates for the human and the C. elegans RD proteins despite the low sequence identity between their CBFbeta proteins and the evolutionary distance between the two systems.