Contribution of phenylalanine side chain intercalation to the TATA-box binding protein–DNA interaction: molecular dynamics and dispersion-corrected density functional theory studies

Substance P and Neurokinin-1 Receptor System in Thyroid Cancer: Potential Targets for New Molecular Therapies

In recent years, numerous approaches have been developed to comprehend the molecular alterations underlying thyroid cancer (TC) oncogenesis and explore novel therapeutic strategies for TC. It is now well established that the neurokinin-1 receptor (NK-1R) is overexpressed in cancer cells and that NK-1R is essential for the viability of cancer cells. The binding of substance P (SP) to NK-1R in neoplastic cells plays a pivotal role in cancer progression by promoting neoplastic cell growth, protecting tumor cells from apoptosis, triggering invasion and metastasis through the enhanced migration of cancer cells, and stimulating endothelial cell proliferation for tumor angiogenesis. Remarkably, all types of human TC (papillary, follicular, medullary, anaplastic), as well as metastatic lesions, exhibit the overexpression of SP and NK-1R compared to the normal thyroid gland. TC cells synthesize and release SP, which exerts its multiple functions through autocrine, paracrine, intracrine, and neuroendocrine processes, including the regulation of tumor burden. Consequently, the secretion of SP from TC results in increased SP levels in plasma, which are significantly higher in TC patients compared to controls. Additionally, NK-1R antagonists have demonstrated a dose-dependent antitumor action. They impair cancer cell proliferation on one side and induce apoptosis of tumor cells on the other side. Furthermore, it has been demonstrated that NK-1R antagonists inhibit neoplastic cell migration, thereby impairing both invasiveness and metastatic abilities, as well as angiogenesis. Given the consistent overexpression of NK-1R in all types of TC, targeting this receptor represents a promising therapeutic approach for TC. Therefore, NK-1R antagonists, such as the drug aprepitant, may represent novel drugs for TC treatment.

C-B-A Test of DNA Force Fields

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

DNA mismatches reveal conformational penalties in protein-DNA recognition

Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.

Photon Dissipation as the Origin of Information Encoding in RNA and DNA

Ultraviolet light incident on organic material can initiate its spontaneous dissipative structuring into chromophores which can catalyze their own replication. This may have been the case for one of the most ancient of all chromophores dissipating the Archean UVC photon flux, the nucleic acids. Oligos of nucleic acids with affinity to particular amino acids which foment UVC photon dissipation would most efficiently catalyze their own reproduction and thus would have been selected through non-equilibrium thermodynamic imperatives which favor dissipation. Indeed, we show here that those amino acids with characteristics most relevant to fomenting UVC photon dissipation are precisely those with greatest chemical affinity to their codons or anticodons. This could provide a thermodynamic basis for the specificity in the amino acid-nucleic acid interaction and an explanation for the accumulation of information in nucleic acids since this information is relevant to the optimization of dissipation of the externally imposed thermodynamic potentials. The accumulation of information in this manner provides a link between evolution and entropy production.

Heterogeneous Dynamical Environment at the Interface of a Protein-DNA Complex

Binding between protein and DNA is an essential process to regulate different biological activities. Two puzzling questions in protein-DNA recognition are (i) how the protein's binding domain identifies the DNA sequence in an aqueous solution and (ii) how the formation of the complex alters the dynamical environment around it. In this work, we present results obtained from molecular dynamics simulations of the N-terminal α-helical domain of the λ-repressor protein (in dimeric form) bound to the corresponding operator DNA. Effects of formation of the complex in modifying the microscopic dynamics of water as well as the kinetics of hydrogen bonds at the interface have been explored. Locally heterogeneous restricted water motions at the complex interface have been observed, the extent of restriction being more significant around the directly bound residues of the protein and the DNA. In particular, the calculation revealed the existence of significantly constrained motionally restricted water layer that can form either bridges around the directly bound residues of the protein and DNA or are engaged in forming water-mediated contacts between a fraction of the unbound residues. More importantly, it is observed that the restricted water motion around the complex is correlated with the hydrogen bond relaxation time scale at the interface. It is further demonstrated that the kinetics of water-water hydrogen bonds involving the bridged water are influenced more due to complex formation.

Flexibility of the Binding Regions of a Protein-DNA Complex and the Structure and Ordering of Interfacial Water

Noncovalent interactions between protein and DNA are important to comprehend different biological activities in living organisms. One important issue is how the protein identifies the target DNA and the influence of the resulting protein-DNA complex on the hydration environment around it. In this study, we have carried out atomistic molecular dynamics simulations of the protein-DNA complex formed by the dimeric form of the α-helical N-terminal domain of the λ-repressor protein with its operator DNA. Local heterogeneous flexibilities of the residues of the protein and the DNA components that are involved in binding and the microscopic structure and ordering of water around those have been investigated in detail. The calculations revealed concurrent existence of highly ordered as well as disordered water molecules at the interface. It is found that a fraction of doubly coordinated water molecules exhibit high degree of ordering at the interface, while the randomly oriented ones are coordinated with three water molecules. The effect has been found to be more around the protein and DNA residues that are in contact in the complexed state. We believe that such highly ordered two-coordinated water molecules are likely to act as an adhesive to facilitate the formation of a protein-DNA complex and maintain its structural stability.

Dynamics changes of CRISPR-Cas9 systems induced by high fidelity mutations

CRISPR-Cas9, a powerful genome editing tool, has widely been applied in biological fields.

Specific DNA sequences allosterically enhance protein-protein interaction in a transcription factor through modulation of protein dynamics: implications for specificity of gene regulation

Most genes are regulated by multiple transcription factors, often assembling into multi-protein complexes in the gene regulatory region. Understanding of the molecular origin of specificity of gene regulatory complex formation in the context of the whole genome is currently inadequate. A phage transcription factor λ-CI forms repressive multi-protein complexes by binding to multiple binding sites in the genome to regulate the lifecycle of the phage. The protein-protein interaction between two DNA-bound λ-CI molecules is stronger when they are bound to the correct pair of binding sites, suggesting allosteric transmission of recognition of correct DNA sequences to the protein-protein interaction interface. Exploration of conformation and dynamics by time-resolved fluorescence anisotropy decay and molecular dynamics suggests a change in protein dynamics to be a crucial factor in mediating allostery. A lattice-based model suggests that DNA-sequence induced allosteric effects could be crucial underlying factors in differentially stabilizing the correct site-specific gene regulatory complexes. We conclude that transcription factors have evolved multiple mechanisms to augment the specificity of DNA-protein interactions in order to achieve an extraordinarily high degree of spatial and temporal specificities of gene regulatory complexes, and DNA-sequence induced allostery plays an important role in the formation of sequence-specific gene regulatory complexes.

Microscopic understanding of the conformational features of a protein–DNA complex

Protein–DNA interactions play crucial roles in different stages of genetic activities, such as replication of genome, initiation of transcription,etc.

Investigating the effect of charged amino acids on DNA conformation in EcoRI-DNA complex: a molecular dynamics simulation study

Sequence-specific binding of proteins to DNA is essential for almost all the cellular processes like transcription, translation, replication, etc. One among the various mechanisms that has been identified so far that contributes to the specificity in protein-DNA interaction is the DNA conformational change. Electrostatic neutralization of the phosphate groups by the positively charged amino acid residues in proteins is thought to bring about such conformational changes in DNA. Here, we employ molecular dynamics simulations to examine the effect of charge on amino acids Lys113, Arg145, and Asp91 which are attached to the scissile phosphate on the conformation of DNA in EcoRI-DNA complex. The results indicate that the charge of these amino acids is essential for maintaining the local conformation of DNA in the EcoRI-bound form. Interestingly, we observe that the positively charged amino acids Lys113 and Arg145 have a long-range influence on the DNA conformation, whereas the negatively charged amino acid Asp91 has only a localized effect on the DNA conformation. The charge on the amino acids also alters the collective dynamics of EcoRI. Collectively, the results shed light on the diversity of the effect of charges on DNA conformation as well as on protein dynamics.

A real-time study of the interaction of TBP with a TATA box-containing duplex identical to an ancestral or minor allele of human gene LEP or TPI

It is known that only a single-nucleotide substitution (SNP: a single nucleotide polymorphism) in the sequence of a TATA box can influence the affinity of the interaction of TBP with the TATA box and contribute to the pathogenesis of complex hereditary human diseases and sometimes may be a cause of monogenic diseases (for instance, β-thalassemia). In the present work, we studied the interaction of human TBP with a double-stranded oligodeoxyribonucleotide (ODN) 15 or 26 bp long identical to a TATA box of promoters of a real-life human gene, TPI or LEP, and labeled with fluorophores TAMRA and FAM. To analyze the interaction of TBP with a TATA box of an ancestral or minor allele (SNP in the TATA box) in real time, we used the stopped-flow method with detection of a Förster resonance energy transfer (FRET) signal. The nature of the resulting kinetic curves reflecting changes in the FRET signal (and therefore of DNA conformation during the interaction with TBP) pointed to a multistage mechanism of the formation of the TBP complex with the TATA-containing ODN. The results showed that with the increasing concentration and length of the ODN, heterogeneity of conformational changes (taking place during the first second of the interaction with TBP) in DNA also increases. In contrast to the initial nonspecific interaction, the subsequent phases strictly depend on TBP concentration: at the TBP:ODN ratio of 10:1, the velocity of change of the FRET signal increases approximately 100-fold.

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.