Structure model of a complex between the factor for inversion stimulation (FIS) and DNA: Modeling protein-DNA complexes with dyad symmetry and known protein structures

A method is presented to predict overall conformations of protein-DNA complexes on the basis of the known three-dimensional structures of the proteins. The method is restricted to proteins with a common twofold symmetry axis, which show only minor conformational changes upon binding to DNA. The method uses a numerical finite difference solution of the linearized Poisson-Boltzmann equation and subsequent energy minimization cycles. Structural parameters—the rotation angle of the DNA relative to the protein around the common symmetry axis, the protein-DNA distance, and intermolecular hydrogen-bonding contacts—are presented for two test cases, DNA bound to CAP (catabolite gene activator protein) and to the Cro-repressor of bacteriophage 434. The DNA curvature in the starting model of the docking procedure was chosen as a smoothed approximation of the conformation found in the X-ray structures of these complexes. The method is further used to predict the unknown structure of the complex between the factor for inversion stimulation (FIS) and DNA, which is bent upon binding to FIS. In contrast to the test cases, the unknown curvature of the starting model is derived from a calibration of electrostatic precalculations for different proteins according to crystallographically observed DNA bending. The results of the modeling are in good accordance with the experimentally observed overall structure of protein-DNA complexes for the two test cases; for FIS, they correspond to several of the experimentally proposed protein-DNA contacts. © 1996 Wiley-Liss, Inc.     

Molecular dynamics simulations of papilloma virus E2 DNA sequences: dynamical models for oligonucleotide structures in solution

The specificity of papilloma virus E2 protein-DNA binding depends critically upon the sequence of a region of the DNA not in direct contact with the protein, and represents one of the simplest known examples of indirect readout. A detailed characterization of this system in solution is important to the further investigation hypothesis of a structural code for DNA recognition by regulatory proteins. In the crystalline state, the E2 DNA oligonucleotide sequence, d(ACCGAATTCGGT), exhibits three different structural forms. We report herein studies of the structure of E2 DNA in solution based on a series of molecular dynamics (MD) simulations including counterions and water, utilizing both the canonical and various crystallographic structures as initial points of departure. All MDs converged on a single dynamical structure of d(ACCGAATTCGGT) in solution. The predicted structure is in close accord with two of the three crystal structures, and indicates that a significant kink in the double helix at the central ApT step in the other crystal molecule may be a packing effect. The dynamical fine structure was analyzed on the basis of helicoidal parameters. The calculated curvature in the sequence was found to originate primarily from YPR steps in the regions flanking the central AATT tract. In order to study the role of structural adaptation of the DNA in the binding process, a subsequent simulation on the 16-mer cognate sequence d(CAACCGAATTCGGTTG) was initiated from the crystallographic coordinates of the bound DNA in the crystal structure of the protein DNA complex. MD simulations starting with the protein-bound form relaxed rapidly back to the dynamical structure predicted from the previous simulations on the uncomplexed DNA. The MD results show that the bound form E2 DNA is a dynamically unstable structure in the absence of protein, and arises as a consequence of both structural changes intrinsic to the sequence and induced by the interaction with protein.     

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.     

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.     

Incorporation of an EDTA-metal complex at a rationally selected site within a protein: application to EDTA-iron DNA affinity cleaving with catabolite gene activator protein (CAP) and Cro.

We have developed a simple procedure to incorporate an EDTA-metal complex at a rationally selected site within a full-length protein. Our procedure has two steps: In step 1, we use site-directed mutagenesis to introduce a unique solvent-accessible cysteine residue at the site of interest. In step 2, we derivatize the resulting protein with S-(2-pyridylthio)cysteaminyl-EDTA-metal, a novel aromatic disulfide derivative of EDTA-metal. We have used this procedure to incorporate an EDTA-iron complex at amino acid 2 of the helix-turn-helix motif of each of two helix-turn-helix motif sequence-specific DNA binding proteins, catabolite gene activator protein (CAP) and Cro, and we have analyzed EDTA-iron-mediated DNA affinity cleavage by the resulting protein derivatives. The CAP derivative cleaves DNA at base pair 2 of the DNA half-site in the protein-DNA complex, and the Cro derivative cleaves DNA at base pairs -3 to 5 of the DNA half-site in the protein-DNA complex. We infer that amino acid 2 of the helix-turn-helix motif of CAP is close to base pair 2 of the DNA half-site in the CAP-DNA complex in solution and that amino acid 2 of the helix-turn-helix motif of Cro is close to base pairs -3 to 5 of the DNA half-site in the Cro-DNA complex in solution.(ABSTRACT TRUNCATED AT 250 WORDS)     

Charge neutralization and DNA bending by the Escherichia coli catabolite activator protein

We are interested in the role of asymmetric phosphate neutralization in DNA bending induced by proteins. We describe an experimental estimate of the actual electrostatic contribution of asymmetric phosphate neutralization to the bending of DNA by the Escherichia coli catabolite activator protein (CAP), a prototypical DNA-bending protein. Following assignment of putative electrostatic interactions between CAP and DNA phosphates based on X-ray crystal structures, appropriate phosphates in the CAP half-site DNA were chemically neutralized by methylphosphonate substitution. DNA shape was then evaluated using a semi-synthetic DNA electrophoretic phasing assay. Our results confirm that the unmodified CAP DNA half-site sequence is intrinsically curved by 26° in the direction enhanced in the complex with protein. In the absence of protein, neutralization of five appropriate phosphates increases DNA curvature to 32° (~23% increase), in the predicted direction. Shifting the placement of the neutralized phosphates changes the DNA shape, suggesting that sequence-directed DNA curvature can be modified by the asymmetry of phosphate neutralization. We suggest that asymmetric phosphate neutralization contributes favorably to DNA bending by CAP, but cannot account for the full DNA deformation.     

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

Molecular dynamics simulations of DNA and a protein-DNA complex including solvent

The results of a recent nanosecond (ns) molecular dynamics (MD) simulation of the d(CGCGAATTCGCG) double helix in water and a 100 ps MD study of the repressor-operator complex are described. The DNA simulations are analyzed in terms of the structural dynamics, fluctuations in the groove width and bending of the helical axis. The results indicate that the ns dynamical trajectory progresses through a series of three substates of B form DNA, with lifetimes of the order of hundreds of picoseconds (ps). An incipient dynamical equilibrium is evident. A comparison of the calculated axis bending with that observed in corresponding crystal structure data is presented. Simulation of the DNA in complex with the protein and that of the free DNA in solution, starting from the crystal conformation, reveal the dynamical changes that occur on complex formation.     

Dynamic Allostery of the Catabolite Activator Protein Revealed by Interatomic Forces

The Catabolite Activator Protein (CAP) is a showcase example for entropic allostery. For full activation and DNA binding, the homodimeric protein requires the binding of two cyclic AMP (cAMP) molecules in an anti-cooperative manner, the source of which appears to be largely of entropic nature according to previous experimental studies. We here study at atomic detail the allosteric regulation of CAP with Molecular dynamics (MD) simulations. We recover the experimentally observed entropic penalty for the second cAMP binding event with our recently developed force covariance entropy estimator and reveal allosteric communication pathways with Force Distribution Analyses (FDA). Our observations show that CAP binding results in characteristic changes in the interaction pathways connecting the two cAMP allosteric binding sites with each other, as well as with the DNA binding domains. We identified crucial relays in the mostly symmetric allosteric activation network, and suggest point mutants to test this mechanism. Our study suggests inter-residue forces, as opposed to coordinates, as a highly sensitive measure for structural adaptations that, even though minute, can very effectively propagate allosteric signals.     

Functional domain motions in proteins on the ~1-100 ns timescale: comparison of neutron spin-echo spectroscopy of phosphoglycerate kinase with molecular-dynamics simulation

Protein function often requires large-scale domain motion. An exciting new development in the experimental characterization of domain motions in proteins is the application of neutron spin-echo spectroscopy (NSE). NSE directly probes coherent (i.e., pair correlated) scattering on the ~1-100 ns timescale. Here, we report on all-atom molecular-dynamics (MD) simulation of a protein, phosphoglycerate kinase, from which we calculate small-angle neutron scattering (SANS) and NSE scattering properties. The simulation-derived and experimental-solution SANS results are in excellent agreement. The contributions of translational and rotational whole-molecule diffusion to the simulation-derived NSE and potential problems in their estimation are examined. Principal component analysis identifies types of domain motion that dominate the internal motion's contribution to the NSE signal, with the largest being classic hinge bending. The associated free-energy profiles are quasiharmonic and the frictional properties correspond to highly overdamped motion. The amplitudes of the motions derived by MD are smaller than those derived from the experimental analysis, and possible reasons for this difference are discussed. The MD results confirm that a significant component of the NSE arises from internal dynamics. They also demonstrate that the combination of NSE with MD is potentially useful for determining the forms, potentials of mean force, and time dependence of functional domain motions in proteins.     

Characterization of the DNA- and dNTP-binding activities of the human cytomegalovirus DNA polymerase catalytic subunit UL54

The catalytic subunit of the human cytomegalovirus DNA polymerase is critical for the replication of the virus. In the present study, we report the expression and purification of a recombinant catalytic subunit of the human cytomegalovirus DNA polymerase expressed in bacteria which retains polymerase activity. As a first step towards elucidating the nature of the interaction between the enzyme, DNA and dNTPs, we have utilized endogenous tryptophan fluorescence to evaluate the binding of ligands to the enzyme. Using this technique, we demonstrate that the minimal DNA-binding site of the enzyme is 6 nt. We also report the first detailed study of the binding kinetics and thermodynamic parameters involved in the interaction between the enzyme, DNA and dNTPs. Our thermodynamic analyses indicate that the initial formation of the enzyme-DNA binary complex is driven by a favourable entropy change, but is also clearly associated with an unfavourable enthalpic contribution. In contrast, the interaction of dNTPs to the binary complex was shown to depend on a completely different mode of binding that is dominated by a favourable enthalpy change and associated with an unfavourable entropy change. In order to provide additional insights into the structural modifications that occur during catalysis, we correlated the effect of DNA and dNTP binding on protein structure using CD. Our results indicate that the enzyme undergoes a first conformational change upon the formation of the protein-DNA binary complex, which is followed by a second structural modification upon dNTP binding. The present study provides a better understanding of the molecular basis of DNA and dNTP recognition by the catalytic subunit of the human cytomegalovirus DNA polymerase.     

A new method for evaluating the specificity of indirect readout in protein-DNA recognition

Proteins recognize a specific DNA sequence not only through direct contact (direct readout) with base pairs but also through sequence-dependent conformation and/or flexibility of DNA (indirect readout). However, it is difficult to assess the contribution of indirect readout to the sequence specificity. What is needed is a straightforward method for quantifying its contributions to specificity. Using Bayesian statistics, we derived the probability of a particular sequence for a given DNA structure from the trajectories of molecular dynamics (MD) simulations of DNAs containing all possible tetramer sequences. Then, we quantified the specificity of indirect readout based on the information entropy associated with the probability. We tested this method with known structures of protein-DNA complexes. This method enabled us to correctly predict those regions where experiments suggested the involvement of indirect readout. The results also indicated new regions where the indirect readout mechanism makes major contributions to the recognition. The present method can be used to estimate the contribution of indirect readout without approximations to the distributions in the conformational ensembles of DNA, and would serve as a powerful tool to study the mechanism of protein-DNA recognition.     

Psoralen photofootprinting of protein-binding sites on DNA

Using a BAL31 exonuclease assay to determine the sites of 4,5',8-trimethylpsoralen photocrosslinking in DNA we have shown that 5'-TA sites which are accessible to psoralen DNA interstrand photocrosslinking in naked DNA become inaccessible when protein, in casu, lambda-repressor E. coli or RNA polymerase are bound at their recognition DNA sequences (OR1 operator or deo1 promoter, respectively). These results show that psoralens can be used as photofootprinting reagents to study specific protein-DNA interactions.     

Examining the contribution of a dA+dT element to the conformation of Escherichia coli integration host factor-DNA complexes.

DNA binding proteins that induce structural changes in DNA are common in both prokaryotes and eukaryotes. Integration host factor (IHF) is a multi-functional DNA binding and bending protein of Escherichia coli that can mediate protein-protein and protein-DNA interactions by bending DNA. Previously we have shown that the presence of a dA+dT element 5'-proximal to an IHF consensus sequence can affect the binding of IHF to a particular site. In this study the contribution of various sequence elements to the formation of IHF-DNA complexes was examined. We show that IHF bends DNA more when it binds to a site containing a dA+dT element upstream of its core consensus element than to a site lacking a dA+dT element. We demonstrate that IHF can be specifically crosslinked to DNA with binding sites either containing or lacking this dA+dT element. These results indicate the importance of flanking DNA and a dA+dT element in the binding and bending of a site by IHF.     

Characterization of NMR relaxation-active motions of a partially folded A-state analogue of ubiquitin

The dominant dynamics of a partially folded A-state analogue of ubiquitin that give rise to NMR 15N spin relaxation have been investigated using molecular dynamics (MD) computer simulations and reorientational quasiharmonic analysis. Starting from the X-ray structure of native ubiquitin with a protonation state corresponding to a low pH, the A-state analogue was generated by a MD simulation of a total length of 33 ns in a 60%/40% methanol/water mixture using a variable temperature scheme to control and speed up the structural transformation. The N-terminal half of the A-state analogue consists of loosely coupled native-like secondary structural elements, while the C-terminal half is mostly irregular in structure. Analysis of dipolar N-H backbone correlation functions reveals reorientational amplitudes and time-scale distributions that are comparable to those observed experimentally. Thus, the trajectory provides a realistic picture of a partially folded protein that can be used for gaining a better understanding of the various types of reorientational motions that are manifested in spin-relaxation parameters of partially folded systems. For this purpose, a reorientational quasiharmonic reorientational analysis was performed on the final 5 ns of the trajectory of the A-state analogue, and for comparison on a 5 ns trajectory of native ubiquitin. The largest amplitude reorientational modes show a markedly distinct behavior for the two states. While for native ubiquitin, such motions have a more local character involving loops and the C-terminal end of the polypeptide chain, the A-state analogue shows highly collective motions in the nanosecond time-scale range corresponding to larger-scale movements between different segments. Changes in reorientational backbone entropy between the A-state analogue and the native state of ubiquitin, which were computed from the reorientational quasiharmonic analyses, are found to depend significantly on motional correlation effects.     

Interfacial water as a "hydration fingerprint" in the noncognate complex of BamHI

The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a "hydration fingerprint" of a given DNA sequence.     

Axis Curvature and Ligand Induced Bending in the CAP-DNA Oligomers

The origin of DNA axis curvature in complexes of the catabolite activator protein with DNA is studied using multiple molecular dynamics (MD) simulations of the free and protein-bound forms of the DNA. The results are compared to available solution and crystal structure data. The MD simulations reproduce the experimentally observed bend in DNA and indicate that ∼40% of the bending observed in the complex is intrinsic to the DNA sequence, whereas ∼60% is induced on protein binding. The MD provides a model for the dynamical structure of the DNA free in solution and for ligand-induced bending.     

Refinement of docked protein-ligand and protein-DNA structures using low frequency normal mode amplitude optimization

Prediction of structural changes resulting from complex formation, both in ligands and receptors, is an important and unsolved problem in structural biology. In this work, we use all-atom normal modes calculated with the Elastic Network Model as a basis set to model structural flexibility during formation of macromolecular complexes and refine the non-bonded intermolecular energy between the two partners (protein-ligand or protein-DNA) along 5-10 of the lowest frequency normal mode directions. The method handles motions unrelated to the docking transparently by first applying the modes that improve non-bonded energy most and optionally restraining amplitudes; in addition, the method can correct small errors in the ligand position when the first six rigid-body modes are switched on. For a test set of six protein receptors that show an open-to-close transition when binding small ligands, our refinement scheme reduces the protein coordinate cRMS by 0.3-3.2 A. For two test cases of DNA structures interacting with proteins, the program correctly refines the docked B-DNA starting form into the expected bent DNA, reducing the DNA cRMS from 8.4 to 4.8 A and from 8.7 to 5.4 A, respectively. A public web server implementation of the refinement method is available at http://lorentz.immstr.pasteur.fr.