Allostery through protein-induced DNA bubbles

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

Influence of a Single Deuterium Substitution for Protium on the Frequency Generation of Different-Size Bubbles in IFNA17

The influence of a single 2H/1H replacement on the frequency generation of different-size bubbles in the human interferon alpha-17 gene (IFNA17) under various energies was studied by a developed algorithm and mathematical modeling without simplifications or averaging. This new approach showed the efficacy of researching DNA bubbles and open states both when all hydrogen bonds in nitrogenous base pairs are protium and after an 2H-substitution. After a single deuterium substitution under specific energies, it was demonstrated that the non-coding region of IFNA17 had a more significant regulatory role in bubble generation in the whole gene than the promoter had. It was revealed that a single deuterium substitution for protium has an influence on the frequency generation of DNA bubbles, which also depends on their size and is always higher for the smaller bubbles under the largest number of the studied energies. Wherein, compared to the natural condition under the same critical value of energy, the bigger raises of the bubble frequency occurrence (maximums) were found for 11-30 base pair (bp) bubbles (higher by 319%), 2-4 bp bubbles (higher by 300%), and 31 bp and over ones (higher by 220%); whereas the most significant reductions of the indicators (minimums) were observed for 11-30 bp bubbles (lower by 43%) and bubbles size over 30 bp (lower by 82%). In this study, we also analyzed the impact of several circumstances on the AT/GC ratio in the formation of DNA bubbles, both under natural conditions and after a single hydrogen isotope exchange. Moreover, based on the obtained data, substantial positive and inverse correlations were revealed between the AT/GC ratio and some factors (energy values, size of DNA bubbles). So, this modeling and variant of the modified algorithm, adapted for researching DNA bubbles, can be useful to study the regulation of replication and transcription in the genes under different isotopic substitutions in the nucleobases.

Disease causing mutation (P178L) in mitochondrial transcription factor A results in impaired mitochondrial transcription initiation

Mitochondrial transcription factor A (TFAM) is essential for the maintenance, expression, and packaging of mitochondrial DNA (mtDNA). Recently, a pathogenic homozygous variant in TFAM (P178L) has been associated with a severe mtDNA depletion syndrome leading to neonatal liver failure and early death. We have performed a biochemical characterization of the TFAM variant P178L in order to understand the molecular basis for the pathogenicity of this mutation. We observe no effects on DNA binding, and compaction of DNA is only mildly affected by the P178L amino acid change. Instead, the mutation severely impairs mtDNA transcription initiation at the mitochondrial heavy and light strand promoters. Molecular modeling suggests that the P178L mutation affects promoter sequence recognition and the interaction between TFAM and the tether helix of POLRMT, thus explaining transcription initiation deficiency.

Bubble lifetimes in DNA gene promoters and their mutations affecting transcription

Relative lifetimes of inherent double stranded DNA openings with lengths up to ten base pairs are presented for different gene promoters and corresponding mutants that either increase or decrease transcriptional activity in the framework of the Peyrard-Bishop-Dauxois model. Extensive microcanonical simulations are used with energies corresponding to physiological temperature. The bubble lifetime profiles along the DNA sequences demonstrate a significant reduction of the average lifetime at the mutation sites when the mutated promoter decreases transcription, while a corresponding enhancement of the bubble lifetime is observed in the case of mutations leading to increased transcription. The relative difference in bubble lifetimes between the mutated and wild type promoters at the position of mutation varies from 20% to more than 30% as the bubble length decreases.

Interactions of Mitochondrial Transcription Factor A with DNA Damage: Mechanistic Insights and Functional Implications

Mitochondria have a plethora of functions in eukaryotic cells, including cell signaling, programmed cell death, protein cofactor synthesis, and various aspects of metabolism. The organelles carry their own genomic DNA, which encodes transfer and ribosomal RNAs and crucial protein subunits in the oxidative phosphorylation system. Mitochondria are vital for cellular and organismal functions, and alterations of mitochondrial DNA (mtDNA) have been linked to mitochondrial disorders and common human diseases. As such, how the cell maintains the integrity of the mitochondrial genome is an important area of study. Interactions of mitochondrial proteins with mtDNA damage are critically important for repairing, regulating, and signaling mtDNA damage. Mitochondrial transcription factor A (TFAM) is a key player in mtDNA transcription, packaging, and maintenance. Due to the extensive contact of TFAM with mtDNA, it is likely to encounter many types of mtDNA damage and secondary structures. This review summarizes recent research on the interaction of human TFAM with different forms of non-canonical DNA structures and discusses the implications on mtDNA repair and packaging.

Ideas and methods of nonlinear mathematics and theoretical physics in DNA science: the McLaughlin-Scott equation and its application to study the DNA open state dynamics

The review is devoted to a new and rapidly developing area related to the application of ideas and methods of nonlinear mathematics and theoretical physics to study the internal dynamics of DNA and, in particular, the behavior of the open states of DNA. There are two main competing approaches to this research. The first approach is based on the molecular dynamics method, which takes into account the motions of all structural elements of the DNA molecule and all interactions between them. The second approach is based on prior selection of the main (dominant) motions and their mathematical description using a small number of model equations. This review describes the results of the study of the open states of DNA performed within the framework of the second approach using the McLaughlin-Scott equation. We present the results obtained both in the case of homogeneous sequences: poly (A), poly (T), poly (G) and poly (C), and in the inhomogeneous case when the McLaughlin-Scott equation has been used for studying the dynamics of open states activated in the promoters A1, A2 and A3 of the bacteriophage T7 genome, in the genes IFNA17, ADRB2, NOS1 and IL-5, in the pBR322 and pTTQ18 plasmids. Particular attention is paid to the results concerning the effect of various external fields on the behavior of open states. In the concluding part of the review, new possibilities and prospects for the development of the considered approach and especially of the McLaughlin-Scott equation are discussed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12551-021-00801-0.

Distributions of bubble lifetimes and bubble lengths in DNA

We investigate the distribution of bubble lifetimes and bubble lengths in DNA at physiological temperature, by performing extensive molecular dynamics simulations with the Peyrard-Bishop-Dauxois (PBD) model, as well as an extended version (ePBD) having a sequence-dependent stacking interaction, emphasizing the effect of the sequences' guanine-cytosine (GC)/adenine-thymine (AT) content on these distributions. For both models we find that base pair-dependent (GC vs AT) thresholds for considering complementary nucleotides to be separated are able to reproduce the observed dependence of the melting temperature on the GC content of the DNA sequence. Using these thresholds for base pair openings, we obtain bubble lifetime distributions for bubbles of lengths up to ten base pairs as the GC content of the sequences is varied, which are accurately fitted with stretched exponential functions. We find that for both models the average bubble lifetime decreases with increasing either the bubble length or the GC content. In addition, the obtained bubble length distributions are also fitted by appropriate stretched exponential functions and our results show that short bubbles have similar likelihoods for any GC content, but longer ones are substantially more likely to occur in AT-rich sequences. We also show that the ePBD model permits more, longer-lived, bubbles than the PBD system.

Organization of DNA in Mammalian Mitochondria

As with all organisms that must organize and condense their DNA to fit within the limited volume of a cell or a nucleus, mammalian mitochondrial DNA (mtDNA) is packaged into nucleoprotein structures called nucleoids. In this study, we first introduce the general modes of DNA compaction, especially the role of the nucleoid-associated proteins (NAPs) that structure the bacterial chromosome. We then present the mitochondrial nucleoid and the main factors responsible for packaging of mtDNA: ARS- (autonomously replicating sequence-) binding factor 2 protein (Abf2p) in yeast and mitochondrial transcription factor A (TFAM) in mammals. We summarize the single-molecule manipulation experiments on mtDNA compaction and visualization of mitochondrial nucleoids that have led to our current knowledge on mtDNA compaction. Lastly, we discuss the possible regulatory role of DNA packaging by TFAM in DNA transactions such as mtDNA replication and transcription.

Heterogeneity and chaos in the Peyrard-Bishop-Dauxois DNA model

We discuss the effect of heterogeneity on the chaotic properties of the Peyrard-Bishop-Dauxois nonlinear model of DNA. Results are presented for the maximum Lyapunov exponent and the deviation vector distribution. Different compositions of adenine-thymine (AT) and guanine-cytosine (GC) base pairs are examined for various energies up to the melting point of the corresponding sequence. We also consider the effect of the alternation index, which measures the heterogeneity of the DNA chain through the number of alternations between different types (AT or GC) of base pairs, on the chaotic behavior of the system. Biological gene promoter sequences have been also investigated, showing no distinct behavior of the maximum Lyapunov exponent.

Acetylation and phosphorylation of human TFAM regulate TFAM-DNA interactions via contrasting mechanisms

Mitochondrial transcription factor A (TFAM) is essential for the maintenance, expression and transmission of mitochondrial DNA (mtDNA). However, mechanisms for the post-translational regulation of TFAM are poorly understood. Here, we show that TFAM is lysine acetylated within its high-mobility-group box 1, a domain that can also be serine phosphorylated. Using bulk and single-molecule methods, we demonstrate that site-specific phosphoserine and acetyl-lysine mimics of human TFAM regulate its interaction with non-specific DNA through distinct kinetic pathways. We show that higher protein concentrations of both TFAM mimics are required to compact DNA to a similar extent as the wild-type. Compaction is thought to be crucial for regulating mtDNA segregation and expression. Moreover, we reveal that the reduced DNA binding affinity of the acetyl-lysine mimic arises from a lower on-rate, whereas the phosphoserine mimic displays both a decreased on-rate and an increased off-rate. Strikingly, the increased off-rate of the phosphoserine mimic is coupled to a significantly faster diffusion of TFAM on DNA. These findings indicate that acetylation and phosphorylation of TFAM can fine-tune TFAM-DNA binding affinity, to permit the discrete regulation of mtDNA dynamics. Furthermore, our results suggest that phosphorylation could additionally regulate transcription by altering the ability of TFAM to locate promoter sites.

Maintenance and Expression of Mammalian Mitochondrial DNA

Mammalian mitochondrial DNA (mtDNA) encodes 13 proteins that are essential for the function of the oxidative phosphorylation system, which is composed of four respiratory-chain complexes and adenosine triphosphate (ATP) synthase. Remarkably, the maintenance and expression of mtDNA depend on the mitochondrial import of hundreds of nuclear-encoded proteins that control genome maintenance, replication, transcription, RNA maturation, and mitochondrial translation. The importance of this complex regulatory system is underscored by the identification of numerous mutations of nuclear genes that impair mtDNA maintenance and expression at different levels, causing human mitochondrial diseases with pleiotropic clinical manifestations. The basic scientific understanding of the mechanisms controlling mtDNA function has progressed considerably during the past few years, thanks to advances in biochemistry, genetics, and structural biology. The challenges for the future will be to understand how mtDNA maintenance and expression are regulated and to what extent direct intramitochondrial cross talk between different processes, such as transcription and translation, is important.

Adaptability of protein structures to enable functional interactions and evolutionary implications

Several studies in recent years have drawn attention to the ability of proteins to adapt to intermolecular interactions by conformational changes along structure-encoded collective modes of motions. These so-called soft modes, primarily driven by entropic effects, facilitate, if not enable, functional interactions. They represent excursions on the conformational space along principal low-ascent directions/paths away from the original free energy minimum, and they are accessible to the protein even before protein-protein/ligand interactions. An emerging concept from these studies is the evolution of structures or modular domains to favor such modes of motion that will be recruited or integrated for enabling functional interactions. Structural dynamics, including the allosteric switches in conformation that are often stabilized upon formation of complexes and multimeric assemblies, emerge as key properties that are evolutionarily maintained to accomplish biological activities, consistent with the paradigm sequence→structure→dynamics→function where 'dynamics' bridges structure and function.