Protein-DNA hydrophobic recognition in the minor groove is facilitated by sugar switching

Traversing DNA-Protein Interactions Between Mesophilic and Thermophilic Bacteria: Implications from Their Cold Shock Response

Cold shock proteins (CSPs) are small, acidic proteins which contain a conserved nucleic acid-binding domain. These perform mRNA translation acting as "RNA chaperones" when triggered by low temperatures initiating their cold shock response. CSP- RNA interactions have been predominantly studied. Our focus will be CSP-DNA interaction examination, to analyse the diverse interaction patterns such as electrostatic, hydrogen and hydrophobic bonding in both thermophilic and mesophilic bacteria. The differences in the molecular mechanism of these contrasting bacterial proteins are studied. Computational techniques such as modelling, energy refinement, simulation and docking were operated to obtain data for comparative analysis. The thermostability factors which stabilise a thermophilic bacterium and their effect on their molecular regulation is investigated. Conformational deviation, atomic residual fluctuations, binding affinity, Electrostatic energy and Solvent Accessibility energy were determined during stimulation along with their conformational study. The study revealed that mesophilic bacteria E. coli CSP have higher binding affinity to DNA than thermophilic G. stearothermophilus. This was further evident by low conformation deviation and atomic fluctuations during simulation.

A Computational Approach for Designing and Validating Small Interfering RNA against SARS-CoV-2 Variants

AIMS

The aim of this study is to develop a novel antiviral strategy capable of efficiently targeting a broad set of SARS-CoV-2 variants.

BACKGROUND

Since the first emergence of SARS-CoV-2, it has rapidly transformed into a global pandemic, posing an unprecedented threat to public health. SARS-CoV-2 is prone to mutation and continues to evolve, leading to the emergence of new variants capable of escaping immune protection achieved due to previous SARS-CoV-2 infections or by vaccination.

OBJECTIVE

RNA interference (RNAi) is a remarkable biological mechanism that can induce gene silencing by targeting complementary mRNA and inhibiting its translation.

METHODS

In this study, using the computational approach, we predicted the most efficient siRNA capable of inhibiting SARS-CoV-2 variants of concern (VoCs).

RESULTS

The presented siRNA was characterized and evaluated for its thermodynamic properties, offsite-target hits, and in silico validation by molecular docking and molecular dynamics simulations (MD) with Human AGO2 protein.

CONCLUSION

The study contributes to the possibility of designing and developing an effective response strategy against existing variants of concerns and preventing further.

Assessment of A- to B- DNA Transitions Utilizing the Drude Polarizable Force Field

In addition to the well-characterized B-form of DNA, duplex DNA can adopt various conformations, such as A or Z-DNA. Though less common, these structures can be induced biologically through protein or ligand interactions or experimentally with niche environmental conditions, such as high salt concentrations or in mixed water-ethanol. Reproducing these alternate structures through molecular dynamics simulations in recent years has been quite challenging with the currently available force fields, simulation techniques, and time scales. In this study, the Drude polarizable force field is tested for its ability to facilitate transitions between A-DNA and B-DNA or maintain A-DNA. Though transitions away from B-DNA were observed in high concentrations of ethanol, the resulting structures had hybrid properties taken from both B-DNA and A-DNA structures. This was also true for A-DNA in ethanol, which lost some of the A-DNA properties that it was expected to maintain. When B-DNA was tested in high salt environments, the resulting B-DNA structures showed no distinguishable differences with the increasing salt concentrations tested. These results with the Drude FF and recent results with additive force fields suggest that at present the current additive and polarizable force fields do not facilitate a complete transition between B- to A-DNA conformations under the conditions simulated. At present, the Drude FF favors A-B DNA hybrid structures when simulated in nonphysiological conditions.

The molecular grammar of protein disorder guiding genome-binding locations

Intrinsically disordered regions (IDRs) direct transcription factors (TFs) towards selected genomic occurrences of their binding motif, as exemplified by budding yeast's Msn2. However, the sequence basis of IDR-directed TF binding selectivity remains unknown. To reveal this sequence grammar, we analyze the genomic localizations of >100 designed IDR mutants, each carrying up to 122 mutations within this 567-AA region. Our data points at multivalent interactions, carried by hydrophobic-mostly aliphatic-residues dispersed within a disordered environment and independent of linear sequence motifs, as the key determinants of Msn2 genomic localization. The implications of our results for the mechanistic basis of IDR-based TF binding preferences are discussed.

An In-silico Approach to Design and Validate siRNA against Monkeypox Virus

INTRODUCTION

The monkeypox virus has emerged as an uncommon zoonotic infection. The recent outbreak of MPXV in Europe and abroad in 2022 presented a major threat to individuals at risk. At present, no specific MPXV vaccinations or medications are available.

METHODS

In this study, we predicted the most effective siRNA against the conserved region of the MPXV and validated the activity by performing molecular docking studies.

RESULTS

Ultimately, the most efficient siRNA molecule was shortlisted against the envelope protein gene (B6R) based on its toxicity, effectivity, thermodynamic stability, molecular interaction, and molecular dynamics simulations (MD) with the Human Argonaute 2 protein.

CONCLUSION

Thus, the strategy may offer a platform for the development of potential antiviral RNA therapeutics that target MPXV at the genomic level.

MA'AT Analysis of Aldofuranosyl Rings: Unbiased Modeling of Conformational Equilibria and Dynamics in Solution

MA'AT analysis has been applied to methyl β-d-ribofuranoside (3) and methyl 2-deoxy-β-d-erythro-pentofuranoside (4) to demonstrate the ability of this new experimental method to determine multi-state conformational equilibria in solution. Density functional theory (DFT) was used to obtain parameterized equations for >20 NMR spin-coupling constants sensitive to furanose ring conformation in 3 and 4, and these equations were used in conjunction with experimental spin-couplings to produce unbiased MA'AT models of ring pseudorotation. These models describe two-state north-south conformational exchange consistent with results obtained from traditional treatments of more limited sets of NMR spin-couplings (e.g., PSEUROT). While PSEUROT, MA'AT, and aqueous molecular dynamics models yielded similar two-state models, MA'AT analysis gives more reliable results since significantly more experimental observables are employed compared to PSEUROT, and no assumptions are needed to render the fitting tractable. MA'AT models indicate a roughly equal distribution of north and south ring conformers of 4 in aqueous (²H₂O) solution compared to ∼80% north forms for 3. Librational motion about the mean pseudorotation phase angles P of the preferred north and south conformers of 3 in solution is more constrained than that for 4. The greater rigidity of the β-ribo ring may be caused by synergistic stereoelectronic effects and/or noncovalent (e.g., hydrogen-bonding) interactions in solution that preferentially stabilize north forms of 3. MA'AT analysis of oligonucleotides and other furanose ring-containing biomolecules promises to improve current experimental models of sugar ring behavior in solution and help reveal context effects on ring conformation in more complex biologically important systems.

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.

Variability of the DNA Backbone Geometry in DNA-Protein Complexes: Experimental Data Analysis

We have analyzed and compared the available experimental data (PDB) on the backbone geometry of the DNA in solution (NMR), in crystals (X-rays), and in complexes with proteins (X-rays and cryo-electron microscopy). The deoxyribose (pseudorotational angle τ₀) and ε/ζ (BI-BII transition in phosphates) flexibilities are practically the same in the four samples. The α/γ mobility is minimal in crystalline DNA: on the histograms, there is one canonical and one noncanonical t/t peak. The α/γ mobility increases in DNA solutions (three more noncanonical peaks) and is maximal in DNA-protein complexes (another additional peak). On a large amount of data, we have confirmed that the three main degrees of freedom of the sugar-phosphate backbone are "orthogonal": changes in any of the angles τ₀, (ζ-ε), and (γ-α) occur, as a rule, at a constant (usually canonical) value of any other. In the DNA-protein complexes, none of the geometrical parameters commonly used to distinguish the A and B forms of DNA, except for Zp and its simpler analog Zp', show an unambiguous correlation with τ₀. Proteins, binding to DNA, in 59% of cases change the local shape of the helix up to the characteristic of the A-form without switching the deoxyribose conformation from south to north. However, we have found simple local characteristics of one nucleotide that correlate with the angles τ₀ and (ζ-ε). These are the angles C3'C1'N* and C4'C3'P(2), respectively. They are orthogonal in DNA-protein complexes exactly as the pair τ₀ and (ζ-ε). Most characteristics of DNA in complexes with proteins are the same in X-ray and in cryo-EM data, except for the histogram for the angle τ₀. We offer a possible explanation for this difference. We also discuss the artifacts on the ε/ζ histogram for DNA in solutions caused by the currently used NMR refinement protocols.

The importance of hydrophobic interactions in the structure of transcription systems

Hydrophobic forces play a crucial role in both the stability of B DNA and its interactions with proteins. In the present study, we postulate that the hydrophobic effect is an essential component in establishing specificity in the interaction transcription factor proteins with their consensus DNA sequence partners. The PDB coordinates of more than 50 transcription systems have been used to analyze the hydrophobic attraction of proteins towards their DNA consensus. This analysis includes computing the hydrophobic energy of the interacting molecules by means of their hydrophobic moments. Hydrophobic moments have successfully been used in previous studies involving self-assembly protein systems. In the present case, in spite of some variability, we found specificity in transcription factors when interacting with their respective consensus DNA sequences. By applying our model of biological membrane pattern for hydrophobic interactions, we postulate that hydrophobic forces constitute the necessary intermediate interaction between the unspecific electrostatic attraction for DNA phosphate groups and the very short-range interaction promoting hydrogen bonds. We conclude that hydrophobic interactions serve as the intermediate force guiding transcriptions factors towards the proper hydrogen bonds to their DNAs.

Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA

Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.

Kinetics of the Conformational Transformation between B- and A-Forms in the Drew-Dickerson Dodecamer

Some DNA sequences in crystals and in complexes with proteins can exist in the forms intermediate between the B- and A-DNA. Based on this, it was implied that the B-to-A transition for any DNA molecule should go through these intermediate forms also in kinetics. More precisely, the helix parameter Slide has to change first, and the molecule should take the E-form. After that, the Roll parameter changes. In the present work, we simulated the kinetics of the B-A transition in the Drew-Dickerson dodecamer, a known B-philic DNA oligomer. We used the "sugar" coarse-grained model that reproduces ribose flexibility, preserves sequence specificity, employs implicit water and explicit ions, and offers the possibility to vary friction. As the control parameter of the transition, we chose the volume available for a counterion and considered the change from a large to a small volume. In the described system, the B-to-A conformational transformation proved to correspond to a first-order phase transition. The molecule behaves like a small cluster in the region of such a transition, jumping between the A- and B-forms in a wide range of available volumes. The viscosity of the solvent does not affect the midpoint of the transition but only the overall mobility of the system. All helix parameters change synchronously on average, we have not observed the sequence "Slide first, Roll later" in kinetics, and the E-DNA is not a necessary step for the transition between the B- and A-forms in the studied system. So, the existence of the intermediate DNA forms requires specific conditions, shifting the common balance of interactions: certain nucleotide sequence in specific solution or/and the interaction with some protein.

Aptamer-Based Affinity Chromatography for Protein Extraction and Purification

Aptamers are oligonucleotide molecules able to recognize very specifically proteins. Among the possible applications, aptamers have been used for affinity chromatography with effective results and advantages over most advanced protein separation technologies. This chapter first discusses the context of the affinity chromatography with aptamer ligands. With the adaptation of SELEX, the chemical modifications of aptamers to comply with the covalent coupling and the separation process are then extensively presented. A focus is then made about the most important applications for protein separation with real-life examples and the comparison with immunoaffinity chromatography. In spite of well-advanced demonstrations and the extraordinary potential developments, a significant optimization work is still due to deserve large-scale applications with all necessary validations. Graphical Abstract Aptamer-protein complexes by X-ray crystallography.

Hydrophobic Amino Acids as Universal Elements of Protein-Induced DNA Structure Deformation

Interaction with the DNA minor groove is a significant contributor to specific sequence recognition in selected families of DNA-binding proteins. Based on a statistical analysis of 3D structures of protein–DNA complexes, we propose that distortion of the DNA minor groove resulting from interactions with hydrophobic amino acid residues is a universal element of protein–DNA recognition. We provide evidence to support this by associating each DNA minor groove-binding amino acid residue with the local dimensions of the DNA double helix using a novel algorithm. The widened DNA minor grooves are associated with high GC content. However, some AT-rich sequences contacted by hydrophobic amino acids (e.g., phenylalanine) display extreme values of minor groove width as well. For a number of hydrophobic amino acids, distinct secondary structure preferences could be identified for residues interacting with the widened DNA minor groove. These results hold even after discarding the most populous families of minor groove-binding proteins.

Interaction of ASOs with PC4 Is Highly Influenced by the Cellular Environment and ASO Chemistry

The activity of PS-ASOs is strongly influenced by association with both inter- and intracellular proteins. The sequence, chemical nature, and structure of the ASO can have profound influences on the interaction of PS-ASOs with specific proteins. A more thorough understanding of how these pharmacological agents interact with various proteins and how chemical modifications, sequence, and structure influence interactions with proteins is needed to inform future ASO design efforts. To better understand the chemistry of PS-ASO interactions, we have focused on human positive cofactor 4 (PC4). Although several studies have investigated the in vitro binding properties of PC4 with endogenous nucleic acids, little is known about the chemistry of interaction of PS-ASOs with this protein. Here we examine in detail the impact of ASO backbone chemistry, 2'-modifications, and buffer environment on the binding affinity of PC4. In addition, using site-directed mutagenesis, we identify those amino acids that are specifically required for ASO binding interactions, and by substitution of abasic nucleotides we identify the positions on the ASO that most strongly influence affinity for PC4. Finally, to confirm that the interactions observed in vitro are biologically relevant, we use a recently developed complementation reporter system to evaluate the kinetics and subcellular localization of the interaction of ASO and PC4 in live cells.

Kinetic and subcellular analysis of PS-ASO/protein interactions with P54nrb and RNase H1

The rapid RNase H1-dependent mislocalization of heterodimer proteins P54nrb and PSF to nucleoli is an early event in the pathway that explains the effects of most toxic phosphorothioate ASOs (PS-ASOs). Using a recently developed NanoLuciferace (NLuc)-based structural complementation reporter system which allows us to observe ASO/protein interactions in real time in live cells, we have determined that safe and toxic PS-ASOs associate with these proteins with kinetics and impact on subcellular localization that differ. Toxic PS-ASOs interact in a complex that includes RNase H1, P54nrb and PSF; but RNase H1/P54nrb complexes were observed in only the cells treated with toxic, but not safe PS-ASOs. In addition, experiments performed in vitro suggest that RNA is also a required component of the complex. The protein–protein interaction between P54nrb and RNase H1 requires the spacer region of RNAse H1, while the P54nrb core domains are required for association with RNase H1. In addition, we have determined that PS-ASOs bind P54nrb via RRM1 and RRM2, while they bind RNase H1 primarily via the hybrid binding domain, however catalytic domain interactions also contribute to overall affinity. These ASO–protein interactions are highly influenced by the chemistry of the PS-ASO binding environment, however little correlation between affinity for specific proteins and PS-ASO toxicity was observed.

AMOEBA Polarizable Atomic Multipole Force Field for Nucleic Acids

The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, 3 J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.

A Multidimensional B-Spline Correction for Accurate Modeling Sugar Puckering in QM/MM Simulations

The computational efficiency of approximate quantum mechanical methods allows their use for the construction of multidimensional reaction free energy profiles. It has recently been demonstrated that quantum models based on the neglect of diatomic differential overlap (NNDO) approximation have difficulty modeling deoxyribose and ribose sugar ring puckers and thus limit their predictive value in the study of RNA and DNA systems. A method has been introduced in our previous work to improve the description of the sugar puckering conformational landscape that uses a multidimensional B-spline correction map (BMAP correction) for systems involving intrinsically coupled torsion angles. This method greatly improved the adiabatic potential energy surface profiles of DNA and RNA sugar rings relative to high-level ab initio methods even for highly problematic NDDO-based models. In the present work, a BMAP correction is developed, implemented, and tested in molecular dynamics simulations using the AM1/d-PhoT semiempirical Hamiltonian for biological phosphoryl transfer reactions. Results are presented for gas-phase adiabatic potential energy surfaces of RNA transesterification model reactions and condensed-phase QM/MM free energy surfaces for nonenzymatic and RNase A-catalyzed transesterification reactions. The results show that the BMAP correction is stable, efficient, and leads to improvement in both the potential energy and free energy profiles for the reactions studied, as compared with ab initio and experimental reference data. Exploration of the effect of the size of the quantum mechanical region indicates the best agreement with experimental reaction barriers occurs when the full CpA dinucleotide substrate is treated quantum mechanically with the sugar pucker correction.

Development of a Quantitative BRET Affinity Assay for Nucleic Acid-Protein Interactions

Protein-nucleic acid interactions play a crucial role in the regulation of diverse biological processes. Elucidating the roles that protein-nucleic acid complexes play in the regulation of transcription, translation, DNA replication, repair and recombination, and RNA processing continues to be a crucial aspect of understanding of cell biology and the mechanisms of disease. In addition, proteins have been demonstrated to interact with antisense oligonucleotide therapeutics in a sequence and chemistry dependent manner, influencing ASO potency and distribution in cells and in vivo. While many assays have been developed to measure protein-nucleic acid interactions, many suffer from lack of throughput and sensitivity, or challenges with protein purification and scalability. In this report we present a new BRET assay for the analysis of DNA-protein interactions which makes use of an extremely bright luciferase as a tag for the binding protein, along with a long-wavelength fluorophore conjugated to the nucleic acid. The resulting assay is high throughput, sensitive, does not require protein purification, and even allows for quantitative characterization of these interactions within the biologically relevant context of whole cells.

Photochromic switching of the DNA helicity induced by azobenzene derivatives

The photochromic properties of azobenzene, involving conformational changes occurring upon interaction with light, provide an excellent tool to establish new ways of selective regulation applied to biosystems. We report here on the binding of two water-soluble 4-(phenylazo)benzoic acid derivatives (Azo-2N and Azo-3N) with double stranded DNA and demonstrate that the photoisomerization of Azo-3N leads to changes in DNA structure. In particular, we show that stabilization and destabilization of the B-DNA secondary structure can be photochemically induced in situ by light. This photo-triggered process is fully reversible and could be an alternative pathway to control a broad range of biological processes. Moreover, we found that the bicationic Azo-3N exhibited a higher DNA-binding constant than the monocationic Azo-2N pointing out that the number of positive charges along the photosensitive polyamines chain plays a pivotal role in stabilizing the photochrome-DNA complex.

High-resolution biophysical analysis of the dynamics of nucleosome formation

We describe a biophysical approach that enables changes in the structure of DNA to be followed during nucleosome formation in in vitro reconstitution with either the canonical "Widom" sequence or a judiciously mutated sequence. The rapid non-perturbing photochemical analysis presented here provides 'snapshots' of the DNA configuration at any given moment in time during nucleosome formation under a very broad range of reaction conditions. Changes in DNA photochemical reactivity upon protein binding are interpreted as being mainly induced by alterations in individual base pair roll angles. The results strengthen the importance of the role of an initial (H3/H4)2 histone tetramer-DNA interaction and highlight the modulation of this early event by the DNA sequence. (H3/H4)2 binding precedes and dictates subsequent H2A/H2B-DNA interactions, which are less affected by the DNA sequence, leading to the final octameric nucleosome. Overall, our results provide a novel, exciting way to investigate those biophysical properties of DNA that constitute a crucial component in nucleosome formation and stabilization.

Novel nucleosomal particles containing core histones and linker DNA but no histone H1

Eukaryotic chromosomal DNA is assembled into regularly spaced nucleosomes, which play a central role in gene regulation by determining accessibility of control regions. The nucleosome contains ∼147 bp of DNA wrapped ∼1.7 times around a central core histone octamer. The linker histone, H1, binds both to the nucleosome, sealing the DNA coils, and to the linker DNA between nucleosomes, directing chromatin folding. Micrococcal nuclease (MNase) digests the linker to yield the chromatosome, containing H1 and ∼160 bp, and then converts it to a core particle, containing ∼147 bp and no H1. Sequencing of nucleosomal DNA obtained after MNase digestion (MNase-seq) generates genome-wide nucleosome maps that are important for understanding gene regulation. We present an improved MNase-seq method involving simultaneous digestion with exonuclease III, which removes linker DNA. Remarkably, we discovered two novel intermediate particles containing 154 or 161 bp, corresponding to 7 bp protruding from one or both sides of the nucleosome core. These particles are detected in yeast lacking H1 and in H1-depleted mouse chromatin. They can be reconstituted in vitro using purified core histones and DNA. We propose that these ‘proto-chromatosomes’ are fundamental chromatin subunits, which include the H1 binding site and influence nucleosome spacing independently of H1.

The 3D structure of Kaposi sarcoma herpesvirus LANA C-terminal domain bound to DNA

Kaposi sarcoma herpesvirus (KSHV) persists as a latent nuclear episome in dividing host cells. This episome is tethered to host chromatin to ensure proper segregation during mitosis. For duplication of the latent genome, the cellular replication machinery is recruited. Both of these functions rely on the constitutively expressed latency-associated nuclear antigen (LANA) of the virus. Here, we report the crystal structure of the KSHV LANA DNA-binding domain (DBD) in complex with its high-affinity viral target DNA, LANA binding site 1 (LBS1), at 2.9 Å resolution. In contrast to homologous proteins such as Epstein-Barr virus nuclear antigen 1 (EBNA-1) of the related γ-herpesvirus Epstein-Barr virus, specific DNA recognition by LANA is highly asymmetric. In addition to solving the crystal structure, we found that apart from the two known LANA binding sites, LBS1 and LBS2, LANA also binds to a novel site, denoted LBS3. All three sites are located in a region of the KSHV terminal repeat subunit previously recognized as a minimal replicator. Moreover, we show that the LANA DBD can coat DNA of arbitrary sequence by virtue of a characteristic lysine patch, which is absent in EBNA-1 of the Epstein-Barr virus. Likely, these higher-order assemblies involve the self-association of LANA into supermolecular spirals. One such spiral assembly was solved as a crystal structure of 3.7 Å resolution in the absence of DNA. On the basis of our data, we propose a model for the controlled nucleation of higher-order LANA oligomers that might contribute to the characteristic subnuclear KSHV microdomains ("LANA speckles"), a hallmark of KSHV latency.

Detection of an en masse and reversible B- to A-DNA conformational transition in prokaryotes in response to desiccation

The role that DNA conformation plays in the biochemistry of cells has been the subject of intensive research since DNA polymorphism was discovered. B-DNA has long been considered the native form of DNA in cells although alternative conformations of DNA are thought to occur transiently and along short tracts. Here, we report the first direct observation of a fully reversible en masse conformational transition between B- and A-DNA within live bacterial cells using Fourier transform infrared (FTIR) spectroscopy. This biospectroscopic technique allows for non-invasive and reagent-free examination of the holistic biochemistry of samples. For this reason, we have been able to observe the previously unknown conformational transition in all four species of bacteria investigated. Detection of this transition is evidence of a previously unexplored biological significance for A-DNA and highlights the need for new research into the role that A-DNA plays as a cellular defence mechanism and in stabilizing the DNA conformation. Such studies are pivotal in understanding the role of A-DNA in the evolutionary pathway of nucleic acids. Furthermore, this discovery demonstrates the exquisite capabilities of FTIR spectroscopy and opens the door for further investigations of cell biochemistry with this under-used technique.

Bioinformatic analysis of the protein/DNA interface

To investigate the principles driving recognition between proteins and DNA, we analyzed more than thousand crystal structures of protein/DNA complexes. We classified protein and DNA conformations by structural alphabets, protein blocks [de Brevern, Etchebest and Hazout (2000) (Bayesian probabilistic approach for predicting backbone structures in terms of protein blocks. Prots. Struct. Funct. Genet., 41:271-287)] and dinucleotide conformers [Svozil, Kalina, Omelka and Schneider (2008) (DNA conformations and their sequence preferences. Nucleic Acids Res., 36:3690-3706)], respectively. Assembling the mutually interacting protein blocks and dinucleotide conformers into 'interaction matrices' revealed their correlations and conformer preferences at the interface relative to their occurrence outside the interface. The analyzed data demonstrated important differences between complexes of various types of proteins such as transcription factors and nucleases, distinct interaction patterns for the DNA minor groove relative to the major groove and phosphate and importance of water-mediated contacts. Water molecules mediate proportionally the largest number of contacts in the minor groove and form the largest proportion of contacts in complexes of transcription factors. The generally known induction of A-DNA forms by complexation was more accurately attributed to A-like and intermediate A/B conformers rare in naked DNA molecules.

Sequence-specific transitions of the torsion angle gamma change the polar-hydrophobic profile of the DNA grooves: implication for indirect protein-DNA recognition

Variations of the shape and polarity of the DNA grooves caused by changes of the DNA conformation play an important role in the DNA readout. Despite the fact that non-canonical trans and gauche- conformations of the DNA backbone angle γ (O5'-C5'-C4'-C3') are frequently found in the DNA crystal structures, their possible role in the DNA recognition has not been studied systematically. In order to fill in this gap, we analyze the available high-resolution crystal structures of the naked and complexed DNA. The analysis shows that the non-canonical γ angle conformations are present both in the naked and bound DNA, more often in the bound vs. naked DNA, and in the nucleotides with the A-like vs. the B-like sugar pucker. The alternative angle γ torsions are more frequently observed in the purines with the A-like sugar pucker and in the pyrimidines with the B-like sugar conformation. The minor groove of the nucleotides with non-canonical γ angle conformation is more polar, while the major groove is more hydrophobic than in the nucleotides with the classical γ torsions due to variations in exposure of the polar and hydrophobic groups of the DNA backbone. The propensity of the nucleotides with different γ angle conformations to participate in the protein-nucleic acid contacts in the minor and major grooves is connected with their sugar pucker and sequence-specific. Our findings imply that the angle γ transitions contribute to the process of the protein-DNA recognition due to modification of the polar/hydrophobic profile of the DNA grooves.

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.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

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.

Codons support the maintenance of intrinsic DNA polymer flexibility over evolutionary timescales

Despite our long familiarity with how the genetic code specifies the amino acid sequence, we still know little about why it is organized in the way that it is. Contrary to the view that the organization of the genetic code is a "frozen accident" of evolution, recent studies have demonstrated that it is highly nonrandom, with implications for both codon assignment and usage. We hypothesize that this inherent nonrandomness may facilitate the coexistence of both sequence and structural information in DNA. Here, we take advantage of a simple metric of intrinsic DNA flexibility to analyze mutational effects on the four phosphate linkages present in any given codon. Application of a simple evolutionary neutral model of substitution to random sequences, translated with alternative genetic codes, reveals that the standard code is highly optimized to favor synonymous substitutions that maximize DNA polymer flexibility, potentially counteracting neutral evolutionary drift toward stiffer DNA caused by spontaneous deamination. Comparison to existing mutational patterns in yeast also demonstrates evidence of strong selective constraint on DNA flexibility, especially at so-called "silent" sites. We also report a fundamental relationship between DNA flexibility, codon usage bias, and several important evolutionary descriptors of comparative genomics (e.g., base composition, transition/transversion ratio, and nonsynonymous vs. synonymous substitution rate). Recent advances in structural genomics have emphasized the role of the DNA polymer's flexibility in both gene function and whole genome folding, thereby implicating possible reasons for codons to facilitate the multiplexing of both genetic and structural information within the same molecular context.

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.

Monitoring the reversible B to A-like transition of DNA in eukaryotic cells using Fourier transform infrared spectroscopy

The ability to detect DNA conformation in eukaryotic cells is of paramount importance in understanding how some cells retain functionality in response to environmental stress. It is anticipated that the B to A transition might play a role in resistance to DNA damage such as heat, desiccation and toxic damage. To this end, conformational detail about the molecular structure of DNA has been derived primarily from in vitro experiments on extracted or synthetic DNA. Here, we report that a B- to A-like DNA conformational change can occur in the nuclei of intact cells in response to dehydration. This transition is reversible upon rehydration in air-dried cells. By systematically monitoring the dehydration and rehydration of single and double-stranded DNA, RNA, extracted nuclei and three types of eukaryotic cells including chicken erythrocytes, mammalian lymphocytes and cancerous rodent fibroblasts using Fourier transform infrared (FTIR) spectroscopy, we unequivocally assign the important DNA conformation marker bands within these cells. We also demonstrate that by applying FTIR spectroscopy to hydrated samples, the DNA bands become sharper and more intense. This is anticipated to provide a methodology enabling differentiation of cancerous from non-cancerous cells based on the increased DNA content inherent to dysplastic and neoplastic tissue.

Structural insights into the cTAR DNA recognition by the HIV-1 nucleocapsid protein: role of sugar deoxyriboses in the binding polarity of NC

An essential step of the reverse transcription of the HIV-1 genome is the first strand transfer that requires the annealing of the TAR RNA hairpin to the cTAR DNA hairpin. HIV-1 nucleocapsid protein (NC) plays a crucial role by facilitating annealing of the complementary hairpins. Using nuclear magnetic resonance and gel retardation assays, we investigated the interaction between NC and the top half of the cTAR DNA (mini-cTAR). We show that NC(11-55) binds the TGG sequence in the lower stem that is destabilized by the adjacent internal loop. The 5′ thymine interacts with residues of the N-terminal zinc knuckle and the 3′ guanine is inserted in the hydrophobic plateau of the C-terminal zinc knuckle. The TGG sequence is preferred relative to the apical and internal loops containing unpaired guanines. Investigation of the DNA–protein contacts shows the major role of hydrophobic interactions involving nucleobases and deoxyribose sugars. A similar network of hydrophobic contacts is observed in the published NC:DNA complexes, whereas NC contacts ribose differently in NC:RNA complexes. We propose that the binding polarity of NC is related to these contacts that could be responsible for the preferential binding to single-stranded nucleic acids.

Understanding the sequence-dependence of DNA groove dimensions: implications for DNA interactions

Background

The B-DNA major and minor groove dimensions are crucial for DNA-protein interactions. It has long been thought that the groove dimensions depend on the DNA sequence, however this relationship has remained elusive. Here, our aim is to elucidate how the DNA sequence intrinsically shapes the grooves.

Methodology/Principal Findings

The present study is based on the analysis of datasets of free and protein-bound DNA crystal structures, and from a compilation of NMR ³¹P chemical shifts measured on free DNA in solution on a broad range of representative sequences. The ³¹P chemical shifts can be interpreted in terms of the BI↔BII backbone conformations and dynamics. The grooves width and depth of free and protein-bound DNA are found to be clearly related to the BI/BII backbone conformational states. The DNA propensity to undergo BI↔BII backbone transitions is highly sequence-dependent and can be quantified at the dinucleotide level. This dual relationship, between DNA sequence and backbone behavior on one hand, and backbone behavior and groove dimensions on the other hand, allows to decipher the link between DNA sequence and groove dimensions. It also firmly establishes that proteins take advantage of the intrinsic DNA groove properties.

Conclusions/Significance

The study provides a general framework explaining how the DNA sequence shapes the groove dimensions in free and protein-bound DNA, with far-reaching implications for DNA-protein indirect readout in both specific and non specific interactions.

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

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

Origins of specificity in protein-DNA recognition

Specific interactions between proteins and DNA are fundamental to many biological processes. In this review, we provide a revised view of protein-DNA interactions that emphasizes the importance of the three-dimensional structures of both macromolecules. We divide protein-DNA interactions into two categories: those when the protein recognizes the unique chemical signatures of the DNA bases (base readout) and those when the protein recognizes a sequence-dependent DNA shape (shape readout). We further divide base readout into those interactions that occur in the major groove from those that occur in the minor groove. Analogously, the readout of the DNA shape is subdivided into global shape recognition (for example, when the DNA helix exhibits an overall bend) and local shape recognition (for example, when a base pair step is kinked or a region of the minor groove is narrow). Based on the >1500 structures of protein-DNA complexes now available in the Protein Data Bank, we argue that individual DNA-binding proteins combine multiple readout mechanisms to achieve DNA-binding specificity. Specificity that distinguishes between families frequently involves base readout in the major groove, whereas shape readout is often exploited for higher resolution specificity, to distinguish between members within the same DNA-binding protein family.

Kinetics of bacterial fluorescence staining with 3,3'-diethylthiacyanine

For more than a century, colorimetric and fluorescence staining have been the foundation of a broad range of key bioanalytical techniques. The dynamics of such staining processes, however, still remains largely unexplored. We investigated the kinetics of fluorescence staining of two gram-negative and two gram-positive species with 3,3'-diethylthiacyanine (THIA) iodide. An increase in the THIA fluorescence quantum yield, induced by the bacterial dye uptake, was the principal reason for the observed emission enhancement. The fluorescence quantum yield of THIA depended on the media viscosity and not on the media polarity, which suggested that the microenvironment of the dye molecules taken up by the cells was restrictive. The kinetics of fluorescence staining did not manifest a statistically significant dependence neither on the dye concentration, nor on the cell count. In the presence of surfactant additives, however, the fluorescence-enhancement kinetic patterns manifested species specificity with statistically significant discernibility.

Prediction of functional class of proteins and peptides irrespective of sequence homology by support vector machines

Various computational methods have been used for the prediction of protein and peptide function based on their sequences. A particular challenge is to derive functional properties from sequences that show low or no homology to proteins of known function. Recently, a machine learning method, support vector machines (SVM), have been explored for predicting functional class of proteins and peptides from amino acid sequence derived properties independent of sequence similarity, which have shown promising potential for a wide spectrum of protein and peptide classes including some of the low- and non-homologous proteins. This method can thus be explored as a potential tool to complement alignment-based, clustering-based, and structure-based methods for predicting protein function. This article reviews the strategies, current progresses, and underlying difficulties in using SVM for predicting the functional class of proteins. The relevant software and web-servers are described. The reported prediction performances in the application of these methods are also presented.

Cavities in protein-DNA and protein-RNA interfaces

An analysis of cavities present in protein–DNA and protein–RNA complexes is presented. In terms of the number of cavities and their total volume, the interfaces formed in these complexes are akin to those in transient protein–protein heterocomplexes. With homodimeric proteins protein–DNA interfaces may contain cavities involving both the protein subunits and DNA, and these are more than twice as large as cavities involving a single protein subunit and DNA. A parameter, cavity index, measuring the degree of surface complementarity, indicates that the packing of atoms in protein–protein/DNA/RNA is very similar, but it is about two times less efficient in the permanent interfaces formed between subunits in homodimers. As within the tertiary structure and protein–protein interfaces, protein–DNA interfaces have a higher inclination to be lined by β-sheet residues; from the DNA side, base atoms, in particular those in minor grooves, have a higher tendency to be located in cavities. The larger cavities tend to be less spherical and solvated. A small fraction of water molecules are found to mediate hydrogen-bond interactions with both the components, suggesting their primary role is to fill in the void left due to the local non-complementary nature of the surface patches.

Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs

BACKGROUND

An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.

RESULTS

Most of the dinucleotide steps in free DNA display high flexibility, assuming different conformations in a sequence-dependent fashion. With the exception of the AA/TT and GA/TC steps, which are 'A-phobic', and the GG/CC step, which is 'A-philic', the dinucleotide steps show no preference for A or B forms of DNA. Protein-bound DNA adopts the B-conformation most often. However, in certain cases, protein-binding causes the DNA backbone to take up energetically unfavourable conformations. At the gross structural level, several protein-bound DNA duplexes are observed to assume a curved conformation in the absence of any large distortions, indicating that a series of normal structural parameters at the dinucleotide and trinucleotide level, similar to the ones in free B-DNA, can give rise to curvature at the overall level.

CONCLUSION

The results illustrate that the free DNA molecule, even in the crystalline state, samples a large amount of conformational space, encompassing both the A and the B-forms, in the absence of any large ligands. A-form as well as some non-A, non-B, distorted geometries are observed for a small number of dinucleotide steps in DNA structures bound to the proteins belonging to a few specific families. However, for most of the bound DNA structures, across a wide variety of protein families, the average step parameters for various dinucleotide sequences as well as backbone torsion angles are observed to be quite close to the free 'B-like' DNA oligomer values, highlighting the flexibility and biological significance of this structural form.

Dissection, residue conservation, and structural classification of protein-DNA interfaces

The basic DNA-binding modules of 128 protein-DNA interfaces have been analyzed. Although these are less planar, like the protein-protein interfaces, the protein-DNA interfaces can also be dissected into core regions in which all the fully-buried atoms are located, and rim regions having atoms with residual accessibilities. The sequence entropy of the core residues is smaller than those in the rim, indicating that the former are better conserved and possibly contribute more towards the binding free energy, as has been implicated in protein-protein interactions. On the protein side, 1014 A(2) of the surface is buried of which 63% belong to the core. There are some differences in the propensities of residues to occur in the core and the rim. In the DNA strands, the nucleotide(s) containing fully-buried atoms in all three components usually occupy central positions of the binding region. A new classification scheme for the interfaces has been introduced based on the composition of secondary structural elements of residues and the results compared with the conventional classification of DNA-binding proteins, as well as the protein class of the molecule. It appears that a common framework may be developed to understand both protein-protein and protein-DNA interactions.