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

How DNA affects the hyperthermophilic protein Ape10b2 for oligomerization: an investigation using multiple short molecular dynamics simulations

Alba2 is a hyperthermophilic DNA-binding protein, and DNA plays a crucial role in the Alba2 oligomerization process. It is a pity that there is limited research in terms of how DNA affects the conformational change of Alba2 in oligomerization. Herein, we complement the crystal structure of the Ape10b2 (belongs to Alba2)-dsDNA complex (PDB ID: 3U6Y) and employ multiple short molecular dynamics (MSMD) simulations to illuminate the influence of DNA on Ape10b2 at four temperatures (300, 343, 363, and 373 K). Our results indicate that DNA could cause the conformational changes of two important regions (loop1 and loop5), which may be beneficial for protein oligomerization. The results of hydrogen bond analysis show that the increasing number of hydrogen bonds between two monomers of Ape10b2 may also be a favorable factor for oligomerization. In addition, Ape10b2 can stabilize DNA by electrostatic interactions with an increase in temperature, and five residues (Arg40, Arg42, Asn43, Asn45, and Arg46) play a stabilizing role during protein binding to DNA. Our findings could help in understanding the favorable factors leading to protein oligomerization, which contributes to enzyme engineering research from an industrial perspective.

Probing the role of intercalating protein sidechains for kink formation in DNA

Protein binding can induce DNA kinks, which are for example important to enhance the specificity of the interaction and to facilitate the assembly of multi protein complexes. The respective proteins frequently exhibit amino acid sidechains that intercalate between the DNA base steps at the site of the kink. However, on a molecular level there is only little information available about the role of individual sidechains for kink formation. To unravel structural principles of protein-induced DNA kinking we have performed molecular dynamics (MD) simulations of five complexes that varied in their architecture, function, and identity of intercalated residues. Simulations were performed for the DNA complexes of wildtype proteins (Sac7d, Sox-4, CcpA, TFAM, TBP) and for mutants, in which the intercalating residues were individually or combined replaced by alanine. The work revealed that for systems with multiple intercalated residues, not all of them are necessarily required for kink formation. In some complexes (Sox-4, TBP), one of the residues proved to be essential for kink formation, whereas the second residue has only a very small effect on the magnitude of the kink. In other systems (e.g. Sac7d) each of the intercalated residues proved to be individually capable of conferring a strong kink suggesting a partially redundant role of the intercalating residues. Mutation of the key residues responsible for kinking either resulted in stable complexes with reduced kink angles or caused conformational instability as evidenced by a shift of the kink to an adjacent base step. Thus, MD simulations can help to identify the role of individual inserted residues for kinking, which is not readily apparent from an inspection of the static structures. This information might be helpful for understanding protein-DNA interactions in more detail and for designing proteins with altered DNA binding properties in the future.

Temperature effect on poly(dA).poly(dT): molecular dynamics simulation studies of polymeric and oligomeric constructs

Understanding unwinding and melting of double helical DNA is very important to characterize role of DNA in replication, transcription, translation etc. Sequence dependent melting thermodynamics is used extensively for detecting promoter regions but melting studies are generally done for short oligonucleotides. This study reports several molecular dynamics (MD) simulations of homopolymeric poly(dA).poly(dT) as regular oligonucleotide fragments as well as its corresponding polymeric constructs with water and charge-neutralizing counterions at different temperatures ranging from 300 to 400 K. We have eliminated the end-effect or terminal peeling propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The dynamic properties such as basepairing and stacking geometry, groove width, backbone conformational parameters, bending, distribution of counter ions and number of hydrogen bonds of oligomeric and polymeric constructs of poly(dA).poly(dT) have been analyzed. The oligomer shows terminal fraying or peeling effect at temperatures above 340 K. The polymer shows partial melting at elevated temperatures although complete denaturations of basepairs do not take place. The analysis of cross strand hydrogen bonds shows that the number of N-H···O hydrogen bonds increases with increase in temperature while C-H···O hydrogen bond frequencies decrease with temperature. Restructuring of counterions in the minor groove with temperature appear as initiation of melting in duplex structures.

DNA–RNA hybrid duplexes with decreasing pyrimidine content in the DNA strand provide structural snapshots for the A- to B-form conformational transition of nucleic acids

A gradual increase in the deoxypyrimidine content in DNA–RNA hybrids leads to B- to A-form nucleic acid transition. Possible factors that govern nuclease activity on hybrid duplexes are presented.

Conformational changes of enzymes and DNA in molecular dynamics: influenced by pH, temperature, and ligand

Protein conformation, which has been a research hotspot for human diseases, is an important factor of protein properties. Recently, a series of approaches have been utilized to investigate the conformational changes under different conditions. Some of them have gained promising achievements, but it is still deficient in the detail researches at the atomic level. In this chapter, a series of computational examples of protein conformational changes under different pH environment, temperature, and ligand binding are described. We further show some useful methods, such as constant pH molecular dynamics simulations, molecular docking, and molecular mechanics Poisson-Boltzmann surface area/generalized Born surface area calculations. In comparison with the experimental results, the methods mentioned above are reasonable to detect and predict the interaction between residue and residue, residue and DNA, and residue and ligand. Additionally, some crucial interactions that cause protein conformational changes are discovered and discussed in this chapter. In summary, our work can give penetrating information to understand the pH-, temperature-, and ligand-induced conformational change mechanisms.

Crenarchaeal chromatin proteins Cren7 and Sul7 compact DNA by inducing rigid bends

Archaeal chromatin proteins share molecular and functional similarities with both bacterial and eukaryotic chromatin proteins. These proteins play an important role in functionally organizing the genomic DNA into a compact nucleoid. Cren7 and Sul7 are two crenarchaeal nucleoid-associated proteins, which are structurally homologous, but not conserved at the sequence level. Co-crystal structures have shown that these two proteins induce a sharp bend on binding to DNA. In this study, we have investigated the architectural properties of these proteins using atomic force microscopy, molecular dynamics simulations and magnetic tweezers. We demonstrate that Cren7 and Sul7 both compact DNA molecules to a similar extent. Using a theoretical model, we quantify the number of individual proteins bound to the DNA as a function of protein concentration and show that forces up to 3.5 pN do not affect this binding. Moreover, we investigate the flexibility of the bending angle induced by Cren7 and Sul7 and show that the protein–DNA complexes differ in flexibility from analogous bacterial and eukaryotic DNA-bending proteins.

Influence of hyperthermophilic protein Cren7 on the stability and conformation of DNA: insights from molecular dynamics simulation and free energy analysis

Cren7, a novel chromatin protein highly conserved among crenarchaea, plays an important role in genome packaging and gene regulation. However, the detail dynamical structural characteristic of the Cren7-DNA complex and the detail study of the DNA in the complex have not been done. Focused on two specific Cren7-DNA complexes (PDB codes 3LWH and 3LWI ), we applied molecular dynamics (MD) simulations at four different temperatures (300, 350, 400, and 450 K) and the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculation at 300 and 350 K to examine the role of Cren7 protein in enhancing the stability of DNA duplexes via protein-DNA interactions, and to study the structural transition in DNA. The simulation results indicate that Cren7 stabilizes DNA duplex in a certain temperature range in the binary complex compared with the unbound DNA molecules. At the same time, DNA molecules were found to undergo B-like to A-like form transitions with increased temperature. The results of statistical analyses of the H-bond and hydrophobic contacts show that some residues have significant influence on the structure of DNA molecules. Our work can give important information to understand the interactions of proteins with nucleic acids and other ligands.

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.