Positively charged residues in DNA-binding domains of structural proteins follow sequence-specific positions of DNA phosphate groups

Exploring the molecular aspect and updating evolutionary approaches to the DNA polymerase enzymes for biotechnological needs: A comprehensive review

DNA polymerases are essential enzymes that play a key role in living organisms, as they participate in the synthesis and maintenance of the DNA molecule. The intrinsic properties of these enzymes have been widely observed and studied to understand their functions, activities, and behavior, which has allowed their natural power in DNA synthesis to be exploited in modern biotechnology, to the point of making them true pillars of the field. In this context, the laboratory evolution of these enzymes, either by directed evolution or rational design, has led to the generation of a wide range of new DNA polymerases with novel properties, suitable for a variety of biotechnological needs. In this review, we examine DNA polymerases at the molecular level, their biotechnological use, and their evolutionary methods in relation to the novel properties sought, providing a chronological selection of evolved DNA polymerases cited in the literature that we consider to be of great interest. To our knowledge, this work is the first to bring together the molecular, functional and evolutionary aspects of the DNA polymerase enzyme. We believe it will be of great interest to researchers whose aim is to produce new lines of evolved DNA polymerases.

A patient-based murine model recapitulates human STAT3 gain-of-function syndrome

STAT3 gain-of-function (GOF) variants results in a heterogeneous clinical syndrome characterized by early onset immunodeficiency, multi-organ autoimmunity, and lymphoproliferation. While 191 documented cases with STAT3 GOF variants have been reported, the impact of individual variants on immune regulation and the broad clinical spectrum remains unclear. We developed a Stat3p.L387R mouse model, mirroring a variant identified in a family exhibiting common STAT3 GOF symptoms, and rare phenotypes including pulmonary hypertension and retinal vasculitis. In vitro experiments revealed increased STAT3 phosphorylation, nuclear migration, and DNA binding of the variant. Our Stat3p.L387R model displayed similar traits from previous Stat3GOF strains, such as splenomegaly and lymphadenopathy. Notably, Stat3p.L387R/+ mice exhibited heightened embryonic lethality compared to prior Stat3GOF/+ models and ocular abnormalities were observed. This research underscores the variant-specific pathology in Stat3p.L387R/+ mice, highlighting the ability to recapitulate human STAT3 GOF syndrome in patient-specific transgenic murine models. Additionally, such models could facilitate tailored treatment development.

Equilibria in the aqueous system of cobalt(II) based on 2-picolinehydroxamic acid and N-(2-hydroxybenzyl)phenylalanine and its ability to inhibit the propagation of cancer cells

Mixed-ligand complexes of cobalt(II) with two bioligands, viz. 2-picolinehydroxamic acid and the reduced Schiff base N-(2-hydroxybenzyl)phenylalanine, were studied in aqueous solution by potentiometry and UV-Vis spectroscopic analysis. The coordination mode of the complexes and their stability were determined and compared to their parent species. Stacking interactions between the rings present in the ligands influence the stability of the complexes. Also, UV-Vis spectroscopy revealed that the stacking interactions affected the intercalation of DNA and mixed-ligand complexes. The in vitro anticancer activity of the free ligand 2-picolinehydroxamic acid and the complexes was tested against cervical and gastric human adenocarcinoma epithelial cell lines. At concentrations of 0.06 and 0.11 mM, the mixed-ligand structures showed the ability to reduce gastric cancer cells with no inhibitory effect on mouse fibroblasts. The cytotoxic effect was accompanied by damage to the cell nuclei, which may confirm that the complexes demonstrate effective binding to DNA. No determination of minimal inhibitory and bactericidal/fungicidal concentrations against the test organisms was possible at higher complex concentrations due to precipitation.

Adaptive laboratory evolution of Salmonella enterica in acid stress

Introduction

Adaptive laboratory evolution (ALE) studies play a crucial role in understanding the adaptation and evolution of different bacterial species. In this study, we have investigated the adaptation and evolution of Salmonella enterica serovar Enteritidis to acetic acid using ALE.

Materials and methods

Acetic acid concentrations below the minimum inhibitory concentration (sub-MIC) were used. Four evolutionary lineages (EL), namely, EL1, EL2, EL3, and EL4, of S. Enteritidis were developed, each demonstrating varying levels of resistance to acetic acid.

Results

The acetic acid MIC of EL1 remained constant at 27 mM throughout 70 days, while the MIC of EL2, EL3, and EL4 increased throughout the 70 days. EL4 was adapted to the highest concentration of acetic acid (30 mM) and demonstrated the highest increase in its MIC against acetic acid throughout the study, reaching an MIC of 35 mM on day 70. The growth rates of the evolved lineages increased over time and were dependent on the concentration of acetic acid used during the evolutionary process. EL4 had the greatest increase in growth rate, reaching 0.33 (h-1) after 70 days in the presence of 30 mM acetic acid as compared to EL1, which had a growth rate of 0.2 (h-1) after 70 days with no exposure to acetic acid. Long-term exposure to acetic acid led to an increased MIC of human antibiotics such as ciprofloxacin and meropenem against the S. enterica evolutionary lineages. The MIC of ciprofloxacin for EL1 stayed constant at 0.016 throughout the 70 days while that of EL4 increased to 0.047. Bacterial whole genome sequencing revealed single-nucleotide polymorphisms in the ELs in various genes known to be involved in S. enterica virulence, pathogenesis, and stress response including phoP, phoQ, and fhuA. We also observed genome deletions in some of the ELs as compared to the wild-type S. Enteritidis which may have contributed to the bacterial acid adaptation.

Discussion

This study highlights the potential for bacterial adaptation and evolution under environmental stress and underscores the importance of understanding the development of cross resistance to antibiotics in S. enterica populations. This study serves to enhance our understanding of the pathogenicity and survival strategies of S. enterica under acetic acid stress.

Protein DEK and DTA Aptamers: Insight Into the Interaction Mechanisms and the Computational Aptamer Design

By blocking the DEK protein, DEK-targeted aptamers (DTAs) can reduce the formation of neutrophil extracellular traps (NETs) to reveal a strong anti-inflammatory efficacy in rheumatoid arthritis. However, the poor stability of DTA has greatly limited its clinical application. Thus, in order to design an aptamer with better stability, DTA was modified by methoxy groups (DTA_OMe) and then the exact DEK–DTA interaction mechanisms were explored through theoretical calculations. The corresponding 2′-OCH3-modified nucleotide force field was established and the molecular dynamics (MD) simulations were performed. It was proved that the 2′-OCH3-modification could definitely enhance the stability of DTA on the premise of comparative affinity. Furthermore, the electrostatic interaction contributed the most to the binding of DEK–DTA, which was the primary interaction to maintain stability, in addition to the non-specific interactions between positively-charged residues (e.g., Lys and Arg) of DEK and the negatively-charged phosphate backbone of aptamers. The H-bond network analysis reminded that eight bases could be mutated to probably enhance the affinity of DTA_OMe. Therein, replacing the 29th base from cytosine to thymine of DTA_OMe was theoretically confirmed to be with the best affinity and even better stability. These research studies imply to be a promising new aptamer design strategy for the treatment of inflammatory arthritis.

Multifunctional amphiphilic peptide dendrimer as nonviral gene vectors for effective cancer therapy via combined gene/photodynamic therapies

Gene therapy holds great promise for treatment of gene-associated diseases. However, safe and successful clinical application urgently requires further advancement of constructing efficient delivery systems. Herein, three amphiphilic peptide dendrimers (TTC-L-KRR/KKK/KHH), containing the natural amino acid residues (lysine K, arginine R, and histidine H) and AIE-based photosensitizer (tetraphenylethenethiophene modified cyanoacrylate, TTC) modified with alkyl chain (L), have been designed and prepared for improving therapeutic potency via the combination of gene therapy (GT) and photodynamic therapy (PDT). All three compounds possessed typical aggregation-induced emission (AIE) characteristics and ultralow critical micelle concentrations (CMCs). The liposomes consisting of amphiphilic peptide dendrimers and dioleoylphosphatidylethanolamine (DOPE) can effectively bind DNA into nanoparticles with appropriate sizes, regular morphology and good biocompatibility. Among them, liposomes TTC-L-KKK/DOPE exhibited the highest transfection efficiency up to 5.7-fold as compared with Lipo2000 in HeLa cells. Meanwhile, rapid endocytosis, successful endo/lysosomal escape, gene release and rapid nuclear delivery of DNA revealed the superiority of liposomes TTC-L-KKK/DOPE during gene delivery process. More importantly, efficient reactive oxygen species (ROS) generation by TTC-L-KKK/DOPE led to effective PDT, thus improving therapeutic potency via combining with p53 mediated-gene therapy. Our work brought novel insight and direction for the construction of bio-safe and bio-imaging liposome as the multifunctional nonviral gene vectors for the effective combined gene/photodynamic therapies.

Adsorption of lysozyme into a charged confining pore

Several applications arise from the confinement of proteins on surfaces because their stability and biological activity are enhanced. It is also known that the way in which a protein adsorbs on the surface is important for its biological function since its active sites should not be obstructed. In this study, the adsorption properties of hen egg-white lysozyme, HEWL, into a negatively charged silica pore is examined by employing a coarse-grained model and constant-pH Monte Carlo simulations. The role of electrostatic interactions is taken into account via including the Debye-Hückel potentials into the Cα structure-based model. We evaluate the effects of pH, salt concentration, and pore radius on the protein preferential orientation and spatial distribution of its residues regarding the pore surface. By mapping the residues that stay closer to the pore surface, we find that the increase of pH leads to orientational changes of the adsorbed protein when the solution pH gets closer to the HEWL isoelectric point. Under these conditions, the pK_a shift of these important residues caused by the adsorption into the charged confining surface results in a HEWL charge distribution that stabilizes the adsorption in the observed protein orientation. We compare our observations to the results of the pK_a shift for HEWL available in the literature and to some experimental data.

Induced DNA bending by unique dimerization of HigA antitoxin

The bacterial toxin-antitoxin (TA) system regulates cell growth under various environmental stresses. Mycobacterium tuberculosis, the causative pathogen of tuberculosis (TB), has three HigBA type II TA systems with reverse gene organization, consisting of the toxin protein HigB and labile antitoxin protein HigA. Most type II TA modules are transcriptionally autoregulated by the antitoxin itself. In this report, we first present the crystal structure of the M. tuberculosis HigA3 antitoxin (MtHigA3) and MtHigA3 bound to its operator DNA complex. We also investigated the interaction between MtHigA3 and DNA using NMR spectroscopy. The MtHigA3 antitoxin structure is a homodimer that contains a structurally well conserved DNA-binding domain at the N-terminus and a dimerization domain at the C-terminus. Upon comparing the HigA homologue structures, a distinct difference was found in the C-terminal region that possesses the β-lid, and diverse orientations of two helix-turn-helix (HTH) motifs from HigA homologue dimers were observed. The structure of MtHigA3 bound to DNA reveals that the promoter DNA is bound to two HTH motifs of the MtHigA3 dimer presenting 46.5° bending, and the distance between the two HTH motifs of each MtHigA3 monomer was increased in MtHigA3 bound to DNA. The β-lid, which is found only in the tertiary structure of MtHigA3 among the HigA homologues, causes the formation of a tight dimerization network and leads to a unique arrangement for dimer formation that is related to the curvature of the bound DNA. This work could contribute to the understanding of the HigBA system of M. tuberculosis at the atomic level and may contribute to the development of new antibiotics for TB treatment.

Changes in gene expression predictably shift and switch genetic interactions

Non-additive interactions between mutations occur extensively and also change across conditions, making genetic prediction a difficult challenge. To better understand the plasticity of genetic interactions (epistasis), we combine mutations in a single protein performing a single function (a transcriptional repressor inhibiting a target gene). Even in this minimal system, genetic interactions switch from positive (suppressive) to negative (enhancing) as the expression of the gene changes. These seemingly complicated changes can be predicted using a mathematical model that propagates the effects of mutations on protein folding to the cellular phenotype. More generally, changes in gene expression should be expected to alter the effects of mutations and how they interact whenever the relationship between expression and a phenotype is nonlinear, which is the case for most genes. These results have important implications for understanding genotype-phenotype maps and illustrate how changes in genetic interactions can often—but not always—be predicted by hierarchical mechanistic models.

Non-additive genetic interactions are plastic and can complicate genetic prediction. Here, using deep mutagenesis of the lambda repressor, Li et al. reveal that changes in gene expression can alter the strength and direction of genetic interactions between mutations in many genes and develop mathematical models for predicting them.

Tunable order-disorder continuum in protein-DNA interactions

DNA-binding protein domains (DBDs) sample diverse conformations in equilibrium facilitating the search and recognition of specific sites on DNA over millions of energetically degenerate competing sites. We hypothesize that DBDs have co-evolved to sense and exploit the strong electric potential from the array of negatively charged phosphate groups on DNA. We test our hypothesis by employing the intrinsically disordered DBD of cytidine repressor (CytR) as a model system. CytR displays a graded increase in structure, stability and folding rate on increasing the osmolarity of the solution that mimics the non-specific screening by DNA phosphates. Electrostatic calculations and an Ising-like statistical mechanical model predict that CytR exhibits features of an electric potential sensor modulating its dimensions and landscape in a unique distance-dependent manner, while DNA plays the role of a non-specific macromolecular chaperone. Accordingly, CytR binds its natural half-site faster than the diffusion-controlled limit and even random DNA conforming to an electrostatic-steering binding mechanism. Our work unravels for the first time the synergistic features of a natural electrostatic potential sensor, a novel binding mechanism driven by electrostatic frustration and disorder, and the role of DNA in promoting distance-dependent protein structural transitions critical for switching between specific and non-specific DNA-binding modes.

Protein Interaction with Charged Macromolecules: From Model Polymers to Unfolded Proteins and Post-Translational Modifications

Interaction of proteins with charged macromolecules is involved in many processes in cells. Firstly, there are many naturally occurred charged polymers such as DNA and RNA, polyphosphates, sulfated glycosaminoglycans, etc., as well as pronouncedly charged proteins such as histones or actin. Electrostatic interactions are also important for "generic" proteins, which are not generally considered as polyanions or polycations. Finally, protein behavior can be altered due to post-translational modifications such as phosphorylation, sulfation, and glycation, which change a local charge of the protein region. Herein we review molecular modeling for the investigation of such interactions, from model polyanions and polycations to unfolded proteins. We will show that electrostatic interactions are ubiquitous, and molecular dynamics simulations provide an outstanding opportunity to look inside binding and reveal the contribution of electrostatic interactions. Since a molecular dynamics simulation is only a model, we will comprehensively consider its relationship with the experimental data.

A systematic analysis of nucleosome core particle and nucleosome-nucleosome stacking structure

Chromatin condensation is driven by the energetically favourable interaction between nucleosome core particles (NCPs). The close NCP-NCP contact, stacking, is a primary structural element of all condensed states of chromatin in vitro and in vivo. However, the molecular structure of stacked nucleosomes as well as the nature of the interactions involved in its formation have not yet been systematically studied. Here we undertake an investigation of both the structural and physico-chemical features of NCP structure and the NCP-NCP stacking. We introduce an “NCP-centred” set of parameters (NCP-NCP distance, shift, rise, tilt, and others) that allows numerical characterisation of the mutual positions of the NCPs in the stacking and in any other structures formed by the NCP. NCP stacking in more than 140 published NCP crystal structures were analysed. In addition, coarse grained (CG) MD simulations modelling NCP condensation was carried out. The CG model takes into account details of the nucleosome structure and adequately describes the long range electrostatic forces as well as excluded volume effects acting in chromatin. The CG simulations showed good agreement with experimental data and revealed the importance of the H2A and H4 N-terminal tail bridging and screening as well as tail-tail correlations in the stacked nucleosomes.

Acceleration of bursty multiprotein target search kinetics on DNA by colocalisation

Proteins are capable of locating specific targets on DNA by employing a facilitated diffusion process with intermittent 1D and 3D search steps. We here uncover the implications of colocalisation of protein production and DNA binding sites via computer simulations.

Modeling competitive substitution in a polyelectrolyte complex

We have simulated the invasion of a polyelectrolyte complex made of a polycation chain and a polyanion chain, by another longer polyanion chain, using the coarse-grained united atom model for the chains and the Langevin dynamics methodology. Our simulations reveal many intricate details of the substitution reaction in terms of conformational changes of the chains and competition between the invading chain and the chain being displaced for the common complementary chain. We show that the invading chain is required to be sufficiently longer than the chain being displaced for effecting the substitution. Yet, having the invading chain to be longer than a certain threshold value does not reduce the substitution time much further. While most of the simulations were carried out in salt-free conditions, we show that presence of salt facilitates the substitution reaction and reduces the substitution time. Analysis of our data shows that the dominant driving force for the substitution process involving polyelectrolytes lies in the release of counterions during the substitution.

Theory on the mechanism of site-specific DNA-protein interactions in the presence of traps

The speed of site-specific binding of transcription factor (TFs) proteins with genomic DNA seems to be strongly retarded by the randomly occurring sequence traps. Traps are those DNA sequences sharing significant similarity with the original specific binding sites (SBSs). It is an intriguing question how the naturally occurring TFs and their SBSs are designed to manage the retarding effects of such randomly occurring traps. We develop a simple random walk model on the site-specific binding of TFs with genomic DNA in the presence of sequence traps. Our dynamical model predicts that (a) the retarding effects of traps will be minimum when the traps are arranged around the SBS such that there is a negative correlation between the binding strength of TFs with traps and the distance of traps from the SBS and (b) the retarding effects of sequence traps can be appeased by the condensed conformational state of DNA. Our computational analysis results on the distribution of sequence traps around the putative binding sites of various TFs in mouse and human genome clearly agree well the theoretical predictions. We propose that the distribution of traps can be used as an additional metric to efficiently identify the SBSs of TFs on genomic DNA.

Protein search for multiple targets on DNA

Protein-DNA interactions are crucial for all biological processes. One of the most important fundamental aspects of these interactions is the process of protein searching and recognizing specific binding sites on DNA. A large number of experimental and theoretical investigations have been devoted to uncovering the molecular description of these phenomena, but many aspects of the mechanisms of protein search for the targets on DNA remain not well understood. One of the most intriguing problems is the role of multiple targets in protein search dynamics. Using a recently developed theoretical framework we analyze this question in detail. Our method is based on a discrete-state stochastic approach that takes into account most relevant physical-chemical processes and leads to fully analytical description of all dynamic properties. Specifically, systems with two and three targets have been explicitly investigated. It is found that multiple targets in most cases accelerate the search in comparison with a single target situation. However, the acceleration is not always proportional to the number of targets. Surprisingly, there are even situations when it takes longer to find one of the multiple targets in comparison with the single target. It depends on the spatial position of the targets, distances between them, average scanning lengths of protein molecules on DNA, and the total DNA lengths. Physical-chemical explanations of observed results are presented. Our predictions are compared with experimental observations as well as with results from a continuum theory for the protein search. Extensive Monte Carlo computer simulations fully support our theoretical calculations.

Adsorption of highly charged Gaussian polyelectrolytes onto oppositely charged surfaces

In many biological processes highly charged biopolymers are adsorbed onto oppositely charged surfaces of macroions and membranes. They form strongly correlated structures close to the surface which cannot be explained by the conventional Poisson-Boltzmann theory. In this work strong coupling theory is used to study the adsorption of highly charged Gaussian polyelectrolytes. Two cases of adsorptions are considered, when the Gaussian polyelectrolytes are confined (a) by one charged wall, and (b) between two charged walls. The effects of salt and the geometry of the polymers on their adsorption-depletion transitions in the strong coupling regime are discussed.

Arginine-phosphate salt bridges between histones and DNA: intermolecular actuators that control nucleosome architecture

Structural bioinformatics and van der Waals density functional theory are combined to investigate the mechanochemical impact of a major class of histone-DNA interactions, namely, the formation of salt bridges between arginine residues in histones and phosphate groups on the DNA backbone. Principal component analysis reveals that the configurational fluctuations of the sugar-phosphate backbone display sequence-specific directionality and variability, and clustering of nucleosome crystal structures identifies two major salt-bridge configurations: a monodentate form in which the arginine end-group guanidinium only forms one hydrogen bond with the phosphate, and a bidentate form in which it forms two. Density functional theory calculations highlight that the combination of sequence, denticity, and salt-bridge positioning enables the histones to apply a tunable mechanochemical stress to the DNA via precise and specific activation of backbone deformations. The results suggest that selection for specific placements of van der Waals contacts, with high-precision control of the spatial distribution of intermolecular forces, may serve as an underlying evolutionary design principle for the structure and function of nucleosomes, a conjecture that is corroborated by previous experimental studies.

Which way up? Recognition of homologous DNA segments in parallel and antiparallel alignments

Homologous gene shuffling between DNA molecules promotes genetic diversity and is an important pathway for DNA repair. For this to occur, homologous genes need to find and recognize each other. However, despite its central role in homologous recombination, the mechanism of homology recognition has remained an unsolved puzzle of molecular biology. While specific proteins are known to play a role at later stages of recombination, an initial coarse grained recognition step has, however, been proposed. This relies on the sequence dependence of the DNA structural parameters, such as twist and rise, mediated by intermolecular interactions, in particular, electrostatic ones. In this proposed mechanism, sequences that have the same base pair text, or are homologous, have lower interaction energy than those sequences with uncorrelated base pair texts. The difference between the two energies is termed the "recognition energy." Here, we probe how the recognition energy changes when one DNA fragment slides past another, and consider, for the first time, homologous sequences in antiparallel alignment. This dependence on sliding is termed the "recognition well." We find there is a recognition well for anti-parallel, homologous DNA tracts, but only a very shallow one, so that their interaction will differ little from the interaction between two nonhomologous tracts. This fact may be utilized in single molecule experiments specially targeted to test the theory. As well as this, we test previous theoretical approximations in calculating the recognition well for parallel molecules against MC simulations and consider more rigorously the optimization of the orientations of the fragments about their long axes upon calculating these recognition energies. The more rigorous treatment affects the recognition energy a little, when the molecules are considered rigid. When torsional flexibility of the DNA molecules is introduced, we find excellent agreement between the analytical approximation and simulations.

Prediction of protein folding rates from simplified secondary structure alphabet

Protein folding is a very complicated and highly cooperative dynamic process. However, the folding kinetics is likely to depend more on a few key structural features. Here we find that secondary structures can determine folding rates of only large, multi-state folding proteins and fails to predict those for small, two-state proteins. The importance of secondary structures for protein folding is ordered as: extended β strand > α helix > bend > turn > undefined secondary structure>310 helix > isolated β strand > π helix. Only the first three secondary structures, extended β strand, α helix and bend, can achieve a good correlation with folding rates. This suggests that the rate-limiting step of protein folding would depend upon the formation of regular secondary structures and the buckling of chain. The reduced secondary structure alphabet provides a simplified description for the machine learning applications in protein design.

Inferring coarse-grain histone-DNA interaction potentials from high-resolution structures of the nucleosome

The histone-DNA interaction in the nucleosome is a fundamental mechanism of genomic compaction and regulation, which remains largely unknown despite increasing structural knowledge of the complex. In this paper, we propose a framework for the extraction of a nanoscale histone-DNA force-field from a collection of high-resolution structures, which may be adapted to a larger class of protein-DNA complexes. We applied the procedure to a large crystallographic database extended by snapshots from molecular dynamics simulations. The comparison of the structural models first shows that, at histone-DNA contact sites, the DNA base-pairs are shifted outwards locally, consistent with locally repulsive forces exerted by the histones. The second step shows that the various force profiles of the structures under analysis derive locally from a unique, sequence-independent, quadratic repulsive force-field, while the sequence preferences are entirely due to internal DNA mechanics. We have thus obtained the first knowledge-derived nanoscale interaction potential for histone-DNA in the nucleosome. The conformations obtained by relaxation of nucleosomal DNA with high-affinity sequences in this potential accurately reproduce the experimental values of binding preferences. Finally we address the more generic binding mechanisms relevant to the 80% genomic sequences incorporated in nucleosomes, by computing the conformation of nucleosomal DNA with sequence-averaged properties. This conformation differs from those found in crystals, and the analysis suggests that repulsive histone forces are related to local stretch tension in nucleosomal DNA, mostly between adjacent contact points. This tension could play a role in the stability of the complex.

Ubiquitous human 'master' origins of replication are encoded in the DNA sequence via a local enrichment in nucleosome excluding energy barriers

As the elementary building block of eukaryotic chromatin, the nucleosome is at the heart of the compromise between the necessity of compacting DNA in the cell nucleus and the required accessibility to regulatory proteins. The recent availability of genome-wide experimental maps of nucleosome positions for many different organisms and cell types has provided an unprecedented opportunity to elucidate to what extent the DNA sequence conditions the primary structure of chromatin and in turn participates in the chromatin-mediated regulation of nuclear functions, such as gene expression and DNA replication. In this study, we use in vivo and in vitro genome-wide nucleosome occupancy data together with the set of nucleosome-free regions (NFRs) predicted by a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix, to investigate the role of intrinsic nucleosome occupancy in the regulation of the replication spatio-temporal programme in human. We focus our analysis on the so-called replication U/N-domains that were shown to cover about half of the human genome in the germline (skew-N domains) as well as in embryonic stem cells, somatic and HeLa cells (mean replication timing U-domains). The 'master' origins of replication (MaOris) that border these megabase-sized U/N-domains were found to be specified by a few hundred kb wide regions that are hyper-sensitive to DNase I cleavage, hypomethylated, and enriched in epigenetic marks involved in transcription regulation, the hallmarks of localized open chromatin structures. Here we show that replication U/N-domain borders that are conserved in all considered cell lines have an environment highly enriched in nucleosome-excluding-energy barriers, suggesting that these ubiquitous MaOris have been selected during evolution. In contrast, MaOris that are cell-type-specific are mainly regulated epigenetically and are no longer favoured by a local abundance of intrinsic NFRs encoded in the DNA sequence. At the smaller few hundred bp scale of gene promoters, CpG-rich promoters of housekeeping genes found nearby ubiquitous MaOris as well as CpG-poor promoters of tissue-specific genes found nearby cell-type-specific MaOris, both correspond to in vivo NFRs that are not coded as nucleosome-excluding-energy barriers. Whereas the former promoters are likely to correspond to high occupancy transcription factor binding regions, the latter are an illustration that gene regulation in human is typically cell-type-specific.

Role of indirect readout mechanism in TATA box binding protein-DNA interaction

Gene expression generally initiates from recognition of TATA-box binding protein (TBP) to the minor groove of DNA of TATA box sequence where the DNA structure is significantly different from B-DNA. We have carried out molecular dynamics simulation studies of TBP-DNA system to understand how the DNA structure alters for efficient binding. We observed rigid nature of the protein while the DNA of TATA box sequence has an inherent flexibility in terms of bending and minor groove widening. The bending analysis of the free DNA and the TBP bound DNA systems indicate presence of some similar structures. Principal coordinate ordination analysis also indicates some structural features of the protein bound and free DNA are similar. Thus we suggest that the DNA of TATA box sequence regularly oscillates between several alternate structures and the one suitable for TBP binding is induced further by the protein for proper complex formation.

Electrostatic effect of H1-histone protein binding on nucleosome repeat length

Within a simple biophysical model we describe the effect of electrostatic binding of H1 histone proteins on the nucleosome repeat length in chromatin. The length of wrapped DNA optimizes its binding energy to the histone core and the elastic energy penalty of DNA wrapping. The magnitude of the effect predicted from our model is in agreement with the systematic experimental data on the linear variation of nucleosome repeat lengths with H1/nucleosome ratio (Woodcock C L et al 2006 Chromos. Res. 14 17-25). We compare our model to the data for different cell types and organisms, with a widely varying ratio of bound H1 histones per nucleosome. We underline the importance of this non-specific histone-DNA charge-balance mechanism in regulating the positioning of nucleosomes and the degree of compaction of chromatin fibers in eukaryotic cells.

Critical adsorption of polyelectrolytes onto charged Janus nanospheres

The conditions of critical polyelectrolyte adsorption onto spherical charged Janus nano-particles are exploited by Monte-Carlo computer simulations and theoretically.

Structure-driven homology pairing of chromatin fibers: the role of electrostatics and protein-induced bridging

Chromatin domains formed in vivo are characterized by different types of 3D organization of interconnected nucleosomes and architectural proteins. Here, we quantitatively test a hypothesis that the similarities in the structure of chromatin fibers (which we call "structural homology") can affect their mutual electrostatic and protein-mediated bridging interactions. For example, highly repetitive DNA sequences in heterochromatic regions can position nucleosomes so that preferred inter-nucleosomal distances are preserved on the surfaces of neighboring fibers. On the contrary, the segments of chromatin fiber formed on unrelated DNA sequences have different geometrical parameters and lack structural complementarity pivotal for stable association and cohesion. Furthermore, specific functional elements such as insulator regions, transcription start and termination sites, and replication origins are characterized by strong nucleosome ordering that might induce structure-driven iterations of chromatin fibers. We propose that shape-specific protein-bridging interactions facilitate long-range pairing of chromatin fragments, while for closely-juxtaposed fibers electrostatic forces can in addition yield fine-tuned structure-specific recognition and pairing. These pairing effects can account for some features observed for mitotic and inter-phase chromatins.

Speed-selectivity paradox in the protein search for targets on DNA: is it real or not?

Protein search for targets on DNA starts all major biological processes. Although significant experimental and theoretical efforts have been devoted to investigation of these phenomena, mechanisms of protein-DNA interactions during the search remain not fully understood. One of the most surprising observations is known as a speed-selectivity paradox. It suggests that experimentally observed fast findings of targets require smooth protein-DNA binding potentials, while the stability of the specific protein-DNA complex imposes a large energy gap which should significantly slow down the protein molecule. We developed a discrete-state stochastic approach that allowed us to investigate explicitly target search phenomena and to analyze the speed-selectivity paradox. A general dynamic phase diagram for different search regimes is constructed. The effect of the target position on search dynamics is investigated. Using experimentally observed parameters, it is found that slow protein diffusion on DNA does not lead to an increase in the search times. Thus, our theory resolves the speed-selectivity paradox by arguing that it does not exist. It is just an artifact of using approximate continuum theoretical models for analyzing protein search in the region of the parameter space beyond the range of validity of these models. In addition, the presented method, for the first time, provides an explanation for fast target search at the level of single protein molecules. Our theoretical predictions agree with all available experimental observations, and extensive Monte Carlo computer simulations are performed to support analytical calculations.

Regulation of the H4 tail binding and folding landscapes via Lys-16 acetylation

Intrinsically disordered proteins (IDP) are a broad class of proteins with relatively flat energy landscapes showing a high level of functional promiscuity, which are frequently regulated through posttranslational covalent modifications. Histone tails, which are the terminal segments of the histone proteins, are prominent IDPs that are implicated in a variety of signaling processes, which control chromatin organization and dynamics. Although a large body of work has been done on elucidating the roles of posttranslational modifications in functional regulation of IDPs, molecular mechanisms behind the observed behaviors are not fully understood. Using extensive atomistic molecular dynamics simulations, we found in this work that H4 tail mono-acetylation at LYS-16, which is a key covalent modification, induces a significant reorganization of the tail's conformational landscape, inducing partial ordering and enhancing the propensity for alpha-helical segments. Furthermore, our calculations of the potentials of mean force between the H4 tail and a DNA fragment indicate that contrary to the expectations based on simple electrostatic reasoning, the Lys-16 mono-acetylated H4 tail binds to DNA stronger than the unacetylated protein. Based on these results, we propose a molecular mechanism for the way Lys-16 acetylation might lead to experimentally observed disruption of compact chromatin fibers.

Polyelectrolyte adsorption onto oppositely charged interfaces: image-charge repulsion and surface curvature

We analyze theoretically the influence of low-dielectric boundaries on the adsorption of flexible polyelectrolytes onto planar and spherical oppositely charged surfaces in electrolyte solutions. We rationalize to what extent polymer chains are depleted from adsorbing interfaces by repulsive image forces. We employ the WKB (Wentzel-Kramers-Brillouin) quantum mechanical method for the Green function of the Edwards equation to determine the adsorption equilibrium. Scaling relations are determined for the critical adsorption strength required to initiate polymer adsorption onto these low-dielectric supports. Image-force repulsion shifts the equilibrium toward the desorbed state, demanding larger surface charge densities and polyelectrolyte linear charge densities for the adsorption to take place. The effect is particularly pronounced for a planar interface in a low-salt regime, where a dramatic change in the scaling behavior for the adsorption-desorption transition is predicted. For the adsorbed state, polymers with higher charge densities are displaced further from the interface by image-charge repulsions. We discuss relevant experimental evidence and argue about possible biological applications of the results.

How to accelerate protein search on DNA: location and dissociation

One of the most important features of biological systems that controls their functioning is the ability of protein molecules to find and recognize quickly specific target sites on DNA. Although these phenomena have been studied extensively, detailed mechanisms of protein-DNA interactions during the search are still not well understood. Experiments suggest that proteins typically find their targets fast by combining three-dimensional and one-dimensional motions, and most of the searching time proteins are non-specifically bound to DNA. However these observations are surprising since proteins diffuse very slowly on DNA, and it seems that the observed fast search cannot be achieved under these conditions for single proteins. Here we propose two simple mechanisms that might explain some of these controversial observations. Using first-passage time analysis, it is shown explicitly that the search can be accelerated by changing the location of the target and by effectively irreversible dissociations of proteins. Our theoretical predictions are supported by Monte Carlo computer simulations.

The physics of protein-DNA interaction networks in the control of gene expression

Protein-DNA interaction networks play a central role in many fundamental cellular processes. In gene regulation, physical interactions and reactions among the molecular components together with the physical properties of DNA control how genes are turned on and off. A key player in all these processes is the inherent flexibility of DNA, which provides an avenue for long-range interactions between distal DNA elements through DNA looping. Such versatility enables multiple interactions and results in additional complexity that is remarkably difficult to address with traditional approaches. This topical review considers recent advances in statistical physics methods to study the assembly of protein-DNA complexes with loops, their effects in the control of gene expression, and their explicit application to the prototypical lac operon genetic system of the E. coli bacterium. In the last decade, it has been shown that the underlying physical properties of DNA looping can actively control transcriptional noise, cell-to-cell variability, and other properties of gene regulation, including the balance between robustness and sensitivity of the induction process. These physical properties are largely dependent on the free energy of DNA looping, which accounts for DNA bending and twisting effects. These new physical methods have also been used in reverse to uncover the actual in vivo free energy of looping double-stranded DNA in living cells, which was not possible with existing experimental techniques. The results obtained for DNA looping by the lac repressor inside the E. coli bacterium showed a more malleable DNA than expected as a result of the interplay of the simultaneous presence of two distinct conformations of looped DNA.

The interplay of mutations and electronic properties in disease-related genes

Electronic properties of DNA are believed to play a crucial role in many phenomena in living organisms, for example the location of DNA lesions by base excision repair (BER) glycosylases and the regulation of tumor-suppressor genes such as p53 by detection of oxidative damage. However, the reproducible measurement and modelling of charge migration through DNA molecules at the nanometer scale remains a challenging and controversial subject even after more than a decade of intense efforts. Here we show, by analysing 162 disease-related genes from a variety of medical databases with a total of almost 20,000 observed pathogenic mutations, a significant difference in the electronic properties of the population of observed mutations compared to the set of all possible mutations. Our results have implications for the role of the electronic properties of DNA in cellular processes, and hint at the possibility of prediction, early diagnosis and detection of mutation hotspots.

Towards a molecular view of transcriptional control

The accumulation of experimental data over recent years has fueled theoretical work on how transcription factors (TFs) search for and recognise their DNA target sites, how they interact with one another, or with other DNA-binding proteins, and how they cope with the compaction of DNA within bacterial nucleoids or within eukaryotic chromatin. Many models have been built to study the kinetic, thermodynamic and mechanistic aspects of these questions. In some cases they have resulted in a relatively clear consensus view, but a number of questions remain controversial. We present an overview of recent work, with an emphasis on models that provide, or can inspire, a better understanding of transcriptional control at a detailed molecular level.

Critical polyelectrolyte adsorption under confinement: planar slit, cylindrical pore, and spherical cavity

We explore the properties of adsorption of flexible polyelectrolyte chains in confined spaces between the oppositely charged surfaces in three basic geometries. A method of approximate uniformly valid solutions for the Green function equation for the eigenfunctions of polymer density distributions is developed to rationalize the critical adsorption conditions. The same approach was implemented in our recent study for the "inverse" problem of polyelectrolyte adsorption onto a planar surface, and on the outer surface of rod-like and spherical obstacles. For the three adsorption geometries investigated, the theory yields simple scaling relations for the minimal surface charge density that triggers the chain adsorption, as a function of the Debye screening length and surface curvature. The encapsulation of polyelectrolytes is governed by interplay of the electrostatic attraction energy toward the adsorbing surface and entropic repulsion of the chain squeezed into a thin slit or small cavities. Under the conditions of surface-mediated confinement, substantially larger polymer linear charge densities are required to adsorb a polyelectrolyte inside a charged spherical cavity, relative to a cylindrical pore and to a planar slit (at the same interfacial surface charge density). Possible biological implications are discussed briefly in the end.

Deconstructing honeybee vitellogenin: novel 40 kDa fragment assigned to its N terminus

Vitellogenin, an egg-yolk protein precursor common to oviparous animals, is found abundantly in honeybee workers - a caste of helpers that do not usually lay eggs. Instead, honeybee vitellogenin (180 kDa) participates in processes other than reproduction: it influences hormone signaling, food-related behavior, immunity, stress resistance and longevity. The molecular basis of these functions is largely unknown. Here, we establish and compare the molecular properties of vitellogenin from honeybee hemolymph (blood) and abdominal fat body, two compartments that are linked to vitellogenin functions. Our results reveal a novel 40 kDa vitellogenin fragment in abdominal fat body tissue, the main site for vitellogenin synthesis and storage. Using MALDI-TOF combined with MS/MS mass-spectroscopy, we assign the 40 kDa fragment to the N terminus of vitellogenin, whereas a previously observed 150 kDa fragment corresponded to the remainder of the protein. We show that both protein units are N glycosylated and phosphorylated. Focusing on the novel 40 kDa fragment, we present a homology model based on the structure of lamprey lipovitellin that includes a conserved β-barrel-like shape, with a lipophilic cavity in the interior and two insect-specific loops that have not been described before. Our data indicate that the honeybee fat body vitellogenin experiences cleavage unlike hemolymph vitellogenin, a pattern that can suggest a tissue-specific role. Our experiments advance the molecular understanding of vitellogenin, of which the multiple physiological and behavioral effects in honeybees are well established.