Molecular Dynamics Simulations Shed Light on the Enthalpic and Entropic Driving Forces That Govern the Sequence Specific Recognition between Netropsin and DNA

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.

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.

Dynamic hydrogen bonding and DNA flexibility in minor groove binders: molecular dynamics simulation of the polyamide f-ImPyIm bound to the Mlu1 (MCB) sequence 5'-ACGCGT-3' in 2:1 motif

Molecular dynamics simulations of the DNA 10-mer 5'-CCACGCGTGG-3' alone and complexed with the formamido-imidazole-pyrrole-imidazole (f-ImPyIm) polyamide minor groove binder in a 2:1 fashion were conducted for 50 ns using the pbsc0 parameters within the AMBER 12 software package. The change in DNA structure upon binding of f-ImPyIm was evaluated via minor groove width and depth, base pair parameters of Slide, Twist, Roll, Stretch, Stagger, Opening, Propeller, and x-displacement, dihedral angle distributions of ζ, ε, α, and γ determined using the Curves+ software program, and hydrogen bond formation. The dynamic hydrogen bonding between the f-ImPyIm and its cognate DNA sequence was compared to the static image used to predict sequence recognition by polyamide minor groove binders. Many of the predicted hydrogen bonds were present in less than 50% of the simulation; however, persistent hydrogen bonds between G5/15 and the formamido group of f-ImPyIm were observed. It was determined that the DNA is wider in the Complex than without the polyamide binder; however, there is flexibility in this particular sequence, even in the presence of the f-ImPyIm as evidenced by the range of minor groove widths the DNA exhibits and the dynamics of the hydrogen bonding that binds the two f-ImPyIm ions to the minor groove. The Complex consisting of the DNA and the 2 f-ImPyIm binders shows slight fraying of the 5' end of the 10-mer at the end of the simulation, but the portion of the oligomer responsible for recognition and binding is stable throughout the simulation. Several structural changes in the Complex indicate that minor groove binders may have a more active role in inhibiting transcription than just preventing binding of important transcription factors.

Mechanical Model of DNA Allostery

The importance of allosteric effects in DNA is becoming increasingly appreciated, but the underlying mechanisms remain poorly understood. In this work, we propose a general modeling framework to study DNA allostery. We describe DNA in a coarse-grained manner by intra-base pair and base pair step coordinates, complemented by groove widths. Quadratic deformation energy is assumed, yielding linear relations between the constraints and their effect. Model parameters are inferred from standard unrestrained, explicit-solvent molecular dynamics simulations of naked DNA. We applied the approach to study minor groove binding of diamidines and pyrrole-imidazole polyamides. The predicted DNA bending is in quantitative agreement with experiment and suggests that diamidine binding to the alternating TA sequence brings the DNA closer to the A-tract conformation, with potentially important functional consequences. The approach can be readily applied to other allosteric effects in DNA and generalized to model allostery in various molecular systems.

Truncated variants of the GCN4 transcription activator protein bind DNA with dramatically different dynamical motifs

The yeast protein GCN4 is a transcriptional activator in the basic leucine zipper (bZip) family, whose distinguishing feature is the "chopstick-like" homodimer of alpha helices formed at the DNA-binding interface. While experiments have shown that truncated versions of the protein retain biologically relevant DNA-binding affinity, we present the results of a computational study revealing that these variants show a wide variety of dynamical modes in their interaction with the target DNA sequence. We have performed all-atom molecular dynamics simulations of the full-length GCN4 protein as well as three truncated variants; our data indicate that the truncated mutants show dramatically different correlation patterns. We conclude that although the truncated mutants still retain DNA-binding ability, the bZip interface present in the full-length protein provides important stability for the protein-DNA complex.

Bi and tri-substituted phenyl rings containing bisbenzimidazoles bind differentially with DNA duplexes: a biophysical and molecular simulation study

Recently synthesis of programmable DNA ligands which can regulate transcription factors have increased the interest of researchers on the functional ability of DNA interacting compounds. A series of DNA interacting compounds are being designed which can differentiate between GC and AT rich DNA. In this study, we have studied the specificity of a few novel bisbenzimidazoles having different bi/tri-substituted phenyl rings, with DNA duplexes using spectroscopic methods. This study entails an integrative approach where we combine biophysical methods and molecular dynamics simulation studies to establish suitable scaffolds to target A/T DNA. We have designed a few analogues of Hoechst 33342 viz.; dimethoxy (DMA), trimethoxy (TMA), dichloro (DCA) and difluoro (DFA) functionalities and performed molecular docking of newly designed analogues with biologically relevant AT and GC rich DNA sequences. The docking studies, along with molecular dynamics (MD) simulations of d(ATATATATATATATAT)2, d(GA4T4C)2, d(GT4A4C)2 and GC rich sequence: d(GCGCGCGCGCGCGCGC)2 complexed with DMA, TMA and DFA, showed that these molecules have higher binding affinity towards AT rich DNA. None of these compounds exhibited an affinity to GC rich DNA rather we observed that these compounds destabilize GC rich DNA. The binding was characterized by strong stabilization of the polynucleotides against thermal strand separation in thermal melting experiments. New insights into the molecules binding to DNA have emerged from these studies. All the DNA binding ligands stabilized d(GA4T4C)2 and d(GT4A4C)2 more out of the five oligomers used for the study, suggesting that these ligands bind 'A4T4' and 'T4A4' strongly as compared to 'ATAT' base pairs.

The allosteric modulation of thyroxine-binding globulin affinity is entropy driven

BACKGROUND

Thyroxine-binding globulin (TBG) is a non-inhibitory member of the serpin family of proteins whose main structural element is the reactive center loop (RCL), that, upon cleavage by proteases, is inserted into the protein core adopting a β-strand conformation (stressed to relaxed transition, S-to-R). After S-to-R transition thyroxine (T4) affinity decreases. However, crystallographic studies in the presence or absence of the hormone in different states are unable to show significant differences in the structure and interactions of the binding site. Experimental results also suggest the existence of several S states (differing in the number of inserted RCL residues), associated with a differential affinity.

METHODS

To shed light into the molecular basis that regulates T4 affinity according to the degree of RCL insertion in TBG, we performed extended molecular dynamics simulations combined with several thermodynamic analysis of the T4 binding to TBG in three different S states, and in the R state.

RESULTS

Our results show that, despite T4 binding in the protein by similar interactions in all states, a good correlation between the degree of RCL insertion and the binding affinity, driven by a change in TBG conformational entropy, was observed.

CONCLUSION

TBG allosteric regulation is entropy driven. The presence of multiple S states may allow more efficient T4 release due to protease activity.

GENERAL SIGNIFICANCE

The presented results are clear examples of how computer simulation methods can reveal the thermodynamic basis of allosteric effects, and provide a general framework for understanding serpin allosteric affinity regulation.

Melting thermodynamics of reversible DNA/ligand complexes at interfaces

A variety of solution methods exist for analysis of interactions between small molecule ligands and nucleic acids; however, accomplishing this task economically at the scale of hundreds to thousands of sequences remains challenging. Surface assays offer a prospective solution through array-based multiplexing, capable of mapping out the full sequence context of a DNA/ligand interaction in a single experiment. However, relative to solution assays, accurate quantification of DNA/ligand interactions in a surface format must contend with limited understanding of molecular activities and interactions at a solid-liquid interface. We report a surface adaptation of a solution method in which shifts in duplex stability, induced by ligand binding and quantified from melting transitions, are used for thermodynamic analysis of DNA/ligand interactions. The results are benchmarked against solution calorimetric data. Equilibrium operation is confirmed through superposition of denaturation/hybridization transitions triggered by heating and cooling. The antibiotic compound netropsin, which undergoes electrostatic and sequence-specific minor groove interactions with DNA, is used as a prototypical small molecule. DNA/netropsin interactions are investigated as a function of ionic strength and drug concentration through electrochemical tracing of surface melt transitions. Comparison with solution values finds excellent agreement in free energy, though reliable separation into enthalpic and entropic contributions proves more difficult. The results establish key guidelines for analysis of DNA-ligand interactions via reversible melting denaturation at surfaces.

Conformational fluctuations of a protein-DNA complex and the structure and ordering of water around it

Protein-DNA binding is an important process responsible for the regulation of genetic activities in living organisms. The most crucial issue in this problem is how the protein recognizes the DNA and identifies its target base sequences. Water molecules present around the protein and DNA are also expected to play an important role in mediating the recognition process and controlling the structure of the complex. We have performed atomistic molecular dynamics simulations of an aqueous solution of the protein-DNA complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. The conformational fluctuations of the protein and DNA and the microscopic structure and ordering of water around them in the complex have been explored. In agreement with experimental studies, the calculations reveal conformational immobilization of the terminal segments of the protein on complexation. Importantly, it is discovered that both structural adaptations of the protein and DNA, and the subsequent correlation between them to bind, contribute to the net entropy loss associated with the complex formation. Further, it is found that water molecules around the DNA are more structured with significantly higher density and ordering than that around the protein in the complex.

Dynamic properties of water around a protein-DNA complex from molecular dynamics simulations

Formation of protein-DNA complex is an important step in regulation of genes in living organisms. One important issue in this problem is the role played by water in mediating the protein-DNA interactions. In this work, we have carried out atomistic molecular dynamics simulations to explore the heterogeneous dynamics of water molecules present in different regions around a complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. It is demonstrated that such heterogeneous water motions around the complex are correlated with the relaxation time scales of hydrogen bonds formed by those water molecules with the protein and DNA. The calculations reveal the existence of a fraction of extraordinarily restricted water molecules forming a highly rigid thin layer in between the binding motifs of the protein and DNA. It is further proved that higher rigidity of water layers around the complex originates from more frequent reformations of broken water-water hydrogen bonds. Importantly, it is found that the formation of the complex affects the transverse and longitudinal degrees of freedom of surrounding water molecules in a nonuniform manner.